我想用Gaussian做个bomd,各位老师，同学能帮我看看吗？

蒲昊贤

各位老师同学好，我做吸附氢分子，想看它在300k能不能释放，77k稳不稳定
我的gif文件开头是
%nprocshared=16
%mem=2000MB
#p wb97xd/6-31++g(2d,2p) nosymm
bomd=(update=5,maxpoints=10000,stepsize=10000,rtemp=300,readstop)


我的问题是
Maxpoints=10000 每个轨迹的最大步数

Stepsize=10000  是 10000*0.0001=1fs  相当于 步长为1fs

总的就是10000fs =10 ps

RTemp=300这个就是温度300k

ReadStop 读入另外的停止标准（我这个后面没有跟什么内容，是不是跑到最大步数就停止了呀？或者是我直接不加）

update=5 BOMD默认在每一个点进行二阶求导，使用Update关键词会使程序在计算新的Hessian之前，在n个梯度点上更新hessian。GradientOnly定义只做梯度计算，并且总是更新hessian（不计算完全的二阶微分）。ReCalcFC是该选项的同义词,我不太明白这个

NTraj=N 计算N个轨迹  我这个没有是不是不加也没有关系呢？

算出log文件是不是在log中  例如Predicted information for step  19415 这个地方去找每步的能量

sobereva

除非你要做准经典动力学，否则建议用ORCA做MD，明显易用，看
使用ORCA做从头算动力学(AIMD)的简单例子
http://sobereva.com/576（http://bbs.keinsci.com/thread-20800-1-1.html）

PS：Gaussian、ORCA、xtb做动力学的方法在北京科音高级量子化学培训班（http://www.keinsci.com/KAQC）里有极为全面系统、掰开了揉碎了的讲授，如果学一遍在AIMD方面就精通了，并且你问的这些问题也就一次性全都明白了