计算化学公社

标题: 求助:ORCA做AIMD为什么未达到最大步数就正常退出 [打印本页]

作者
Author:
djjj148    时间: 2021-1-3 20:54
标题: 求助:ORCA做AIMD为什么未达到最大步数就正常退出
学习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步就正常退出呢?



作者
Author:
wzkchem5    时间: 2021-1-3 22:08
看看是不是最后一步的SCF或梯度计算出问题了
作者
Author:
sobereva    时间: 2021-1-3 22:28
AIMD计算过程中对每一步都会把SCF迭代情况、计算受力的情况的信息输出到相应文件里,打开各个相关文件看末尾的信息判断。如果没有异常,那就是ORCA的bug了,可以去ORCA官方论坛问开发者
作者
Author:
djjj148    时间: 2021-1-4 10:05
wzkchem5 发表于 2021-1-3 22:08
看看是不是最后一步的SCF或梯度计算出问题了

谢谢老师,查看了scf.log,的确是scf出问题了。
作者
Author:
djjj148    时间: 2021-1-4 10:06
sobereva 发表于 2021-1-3 22:28
AIMD计算过程中对每一步都会把SCF迭代情况、计算受力的情况的信息输出到相应文件里,打开各个相关文件看末 ...

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




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3