求助：ORCA做AIMD为什么未达到最大步数就正常退出

djjj148

学习sob老师的帖子"使用ORCA做从头算动力学(AIMD)的简单例子"用一个14原子的体系练习了一下，inp文件是用Multiwfn生成的，机器56核，192G内存。设置run 2000结果跑了不到100步就显示"正常退出"。inp和out文件的节选如下：
inp：
! B97-3c noautostart miniprint nopop
%maxcore  3000
%pal nprocs   14 end
%md
#restart ifexists  # Continue MD by reading [basename].mdrestart if it exists. In this case "initvel" should be commented
#minimize  # Do minimization prior to MD simulation
timestep 1.0_fs  # This stepsize is very safe
initvel 298.15_K  # Assign velocity to atoms according to temperature
thermostat berendsen 298.15_K timecon 30.0_fs  # Target temperature and coupling time constant
dump position stride 1 format xyz filename "pos.xyz"  # Dump position every "stride" steps
#dump force stride 1 format xyz filename "force.xyz"  # Dump force every "stride" steps
#dump velocity stride 1 format xyz filename "vel.xyz"  # Dump velocity every "stride" steps
#dump gbw stride 20 filename "wfn"  # Dump wavefunction to "wfn[step].gbw" files every "stride" steps
constraint add center 0..13  # Fix center of mass at initial position
run 2000  # Number of MD steps
end
* xyz   0   1


       86      86.0      8.9      1.1     323.47    0.019975     -344.477523     -344.457548     -344.470350     -0.54
       87      87.0      8.5      1.1     333.56    0.020598     -344.478116     -344.457517     -344.470276      0.01
       88      88.0      8.7      1.2     317.56    0.019610     -344.476953     -344.457343     -344.470094      1.39

out：
Timings for individual modules:

Sum of individual times         ...      689.173 sec (=  11.486 min)
Molecular Dynamics run          ...      689.173 sec (=  11.486 min) 100.0 %
                             ****ORCA TERMINATED NORMALLY****
TOTAL RUN TIME: 0 days 0 hours 11 minutes 29 seconds 882 msec



对应的文件已上传，麻烦大家帮看下为什么没跑到2000步就正常退出呢？

djjj148

sobereva 发表于 2021-1-3 22:28
AIMD计算过程中对每一步都会把SCF迭代情况、计算受力的情况的信息输出到相应文件里，打开各个相关文件看末 ...

谢谢sob老师，查看了scf.log，的确是scf出问题了。

djjj148

wzkchem5 发表于 2021-1-3 22:08
看看是不是最后一步的SCF或梯度计算出问题了

谢谢老师，查看了scf.log，的确是scf出问题了。

sobereva

AIMD计算过程中对每一步都会把SCF迭代情况、计算受力的情况的信息输出到相应文件里，打开各个相关文件看末尾的信息判断。如果没有异常，那就是ORCA的bug了，可以去ORCA官方论坛问开发者

wzkchem5

看看是不是最后一步的SCF或梯度计算出问题了