diaolanxinyu 发表于 2015-3-10 20:28
谢谢sobereva老师，我试试。

想请问一下楼主，我想做一个分子内的BOMD，想计算一个原子进攻另一个原子，这种模拟应该怎么写输入文件啊

谢谢sobereva老师，我试试。

1 建议用M062X/TZVP算算试试，如果计算量有富余，也可以结合更贵的def2-TZVP，热力学数据计算结果应该比MP2/6-31+G*好，计算量也更小。MP2虽然可以，但6-31+G*相比之下显得略穷酸了。

2 两种情况的动力学步进算法不一样，gradientonly下是固定每一步时间，BOMD默认情况下是固定每一步步长。非要用时间来讨论的话，要么gradientonly，要么用ADMP。

sobereva 发表于 2015-3-9 12:31
1 用nosymm避免自动调整朝向

2 可以再加大入射速度看看是否反应。也可以先做个IRC研究摸一下情况。

首先感谢sobereva老师！
1.加了nosymm后解决了自动朝向的问题，而坐标原点会变，这应该是Gaussian默认的，包括在计算完成后也会将原点放置到整个分子的中心位置。

2. 如sobereva所料，我加大了入射速度后反应顺利进行。

3. 我把所有的iop(5/13)都去掉了，加了int=ultrafine或SCF(conver=6)后将SCF收敛的问题解决了。

对于您建议的B3LYP/6-311G**我也试过了，有时也会存在SCF不收敛的问题，但是可以解决。
我不准备用这个方法和基组的原因是，此法在计算静态的热力学数据如delta H，文献中是放热，为负值，而此法计算的delta H为正值。

还有几个问题请教一下：
问题1：对于我这个体系用高斯的BOMD计算动力学，用mp2/6-31+g(d) 的计算精度够精确么？
          我选择此法的依据就是计算静态的热力学数据（delta H和delta E）跟实验值比较后较接近。
        不知道我这样考虑是否充分，是否还需要跟文献中用过的一些更高水平的方法做一下比较呢？

问题二：关于BOMD选项中的StepSize，当StepSize中加GradientOnly后步长变为0.025飞秒，不加为0.25 amu^1/2*Bohr,后者像是长度单位，请问后者0.25 amu^1/2*Bohr较前者有什么优势吗？
        一般在动力学计算的讨论中我们一般说在某时刻（如2.0ps）的一些变化，即以时间为步长；当步长为0.25 amu^1/2*Bohr的时候我突然感觉分析的时候不知该怎么表达了？

1 用nosymm避免自动调整朝向

2 可以再加大入射速度看看是否反应。也可以先做个IRC研究摸一下情况。


不要用iop(5/13=1)，这极为危险。把弥散函数去掉，自然就好收敛了。对于当前体系，弥散函数不是必要的，对于B3LYP，用6-311G**就足够精确了。

114步没有正常结束，停止是因为SCF没收敛，输出文件中可见
>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(UB3LYP) =  -115.693996015     A.U. after   65 cycles
            NFock= 64  Conv=0.80D-04     -V/T= 2.0035

问题三：och4-s-bomd3.out中改变方法后为什么只跑了114步就正常结束了，我设的是500步