简单量子化学问题答疑专帖

NHX

sobereva 发表于 2026-4-9 01:53
涉及势能面二分的时候过渡态会连接两个极小点
参考北京科音中级量子化学培训班（http://www.keinsci.com ...

谢谢sob老师

LorneMeng

sobereva 发表于 2026-4-9 01:55
优化完的位置和IRC端点不完全相同，有时候可能差得不少，这么估计不合理，除非优化过程中结构确实变化甚 ...

谢谢sob老师的解答，我尝试在ORCA6.1中进行NEB-TS计算，思路如下：

首先通过IRC计算，并对两端结构进行几何优化获得反应复合物和产物复合物结构信息；在ORCA进行NEB-TS计算，输入文件如下(目前在个人电脑上尝试，因而使用B3LYP/6-31G*)，计算完成后通过.final.interp文件信息，获得反应复合物和产物复合物相对过渡态位置。

老师，我这样写输入文件，以及这样处理是否合理。

! B3LYP D3 6-31G* NEB-TS TightSCF noautostart miniprint nopop
%maxcore     2000
%pal nprocs  20 end
%neb
  neb_end_xyzfile "ProCom.xyz"
  NImages 15
end
*xyzfile 0 2 ReaCom.xyz


.final.interp文件内容
Iteration: -1
Images: Distance  (Bohr), Energy (Eh)
0.0000   0.00000000      0.00000000
0.0541   0.34967968     -0.00105932
0.1173   0.75876072     -0.00098704
0.2125   1.37463705     -0.00040066
0.3099   2.00485470      0.00071167
0.3792   2.45306505      0.00240670
0.4250   2.74968353      0.00443544
0.4557   2.94791946      0.00642049
0.4893   3.16543194      0.00824847
0.5096   3.29650346      0.00698975
0.5359   3.46692619      0.00191014
0.6110   3.95283265     -0.00893399
0.7446   4.81709105     -0.01350540
0.8314   5.37823299     -0.01412350
0.8938   5.78178369     -0.01413301
0.9485   6.13567259     -0.01397512
1.0000   6.46909210     -0.01281765

sobereva

LorneMeng 发表于 2026-4-9 21:06
谢谢sob老师的解答，我尝试在ORCA6.1中进行NEB-TS计算，思路如下：

首先通过IRC计算，并对两端结构进 ...

能量不是先单调上升后单调下降，不完全正确

LorneMeng

sobereva 发表于 2026-4-10 06:05
能量不是先单调上升后单调下降，不完全正确

感谢sob老师的解答，上次在ORCA6.1中进行NEB-TS计算后，为了快速验证计算思路，ORCA中NEB-TS计算级别与IRC和过渡态计算并不一致，使得计算结果异常。我更新了计算级别，计算完成后.final.interp文件信息显示在过渡态两侧能量单调变化。根据.final.interp文件内容，从而确定反应复合物距离过渡态(s=0)为3.186 Bohr，产物复合物距离过渡态(s=0)为3.402 Bohr

我想咨询下，sob老师，我这样写输入文件，以及处理思路是否合理

orca输入文件
! M062X D3ZERO 6-311++g(d,p) NEB-TS TightSCF noautostart miniprint nopop
%maxcore     2000
%pal nprocs  32 end
%neb
  neb_end_xyzfile "ProCom.xyz"
  NImages 15
end
*xyzfile 0 2 ReaCom.xyz


.final.interp文件内容
Iteration: -1
Images: Distance  (Bohr), Energy (Eh)
0.0000   0.00000000      0.00000000
0.0587   0.38657730      0.00002461
0.1375   0.90583615      0.00025783
0.2418   1.59288762      0.00099287
0.3177   2.09284943      0.00261966
0.3773   2.48553502      0.00494737
0.4185   2.75750497      0.00778956
0.4521   2.97834248      0.01134712
0.4836   3.18628638      0.01339817
0.5134   3.38275719      0.01012776
0.5369   3.53729566      0.00303164
0.5961   3.92703752     -0.01192154
0.7438   4.90056804     -0.02170310
0.8390   5.52774534     -0.02319339
0.9038   5.95438446     -0.02349828
0.9554   6.29448426     -0.02358712
1.0000   6.58838221     -0.02359610

sobereva

LorneMeng 发表于 2026-4-10 09:25
感谢sob老师的解答，上次在ORCA6.1中进行NEB-TS计算后，为了快速验证计算思路，ORCA中NEB-TS计算级别与IR ...

没明显问题