计算化学公社

标题: 使用ORCA进行NEBTS计算不收敛 [打印本页]

---

作者
Author: Iris_1    时间: 2023-6-13 16:46
标题: 使用ORCA进行NEBTS计算不收敛


使用ORCA进行2-fluoro-2,2-dinitroethylamine分子NO2裂解的过渡态搜索，反应物和产物结构都优化到了最低能量状态。NEBTS任务采用64核算了4天CI一直不收敛，不知道问题出在哪里。发现out文件里CI迭代的过程中从Image4变成了Image6又变成了Image3，是什么情况呢。另外，不收敛是否和Image数有关呢，看网上说Image数需要被CPU核数整除。附件为反应物和产物优化过的结构。感谢指点！
NEB-TS.inp:
! B3LYP D4 def2-TZVP def2/J RIJCOSX NEB-TS FREQ verytightSCF noautostart miniprint nopop
%neb
  Product "prod.xyz"
  NImages 6
end

*xyz 0 1
C   2.19862367435602      0.42665204669430      0.36639359440675
C   0.95828244111602      0.90190052283573      1.09291183104736
N   0.63275739981328     -0.02184113959044      2.28949375040034
O   1.50285903012036     -0.06209816986893      3.13587356137438
O  -0.41783219572831     -0.60531338370234      2.27214240951224
N   1.11652506510750      2.31020680143357      1.71748244741355
O   0.11702818629361      2.77312433803615      2.21579143714780
O   2.21272210764266      2.82194357502574      1.64439751211242
F  -0.11133945651520      0.94896846387698      0.29806444460181
N   1.98895272236662     -0.94864192661775     -0.04106268473325
H   3.03367168489204      0.48623100882471      1.05928549815136
H   2.38330339133929      1.14250017613629     -0.44287830421236
H   2.86228621379426     -1.44125906796708     -0.15616839055244
H   1.44983973540186     -1.02686324511692     -0.89258710666995
*

---

作者
Author: wzkchem5    时间: 2023-6-13 17:04
首先，产物是双自由基，所以过渡态也很可能有双自由基特征。除非你已经发现过渡态附近波函数是稳定的，否则应该寻找自旋破缺波函数，不然即使收敛了结果也不能用

---

作者
Author: Iris_1    时间: 2023-6-13 17:23

wzkchem5 发表于 2023-6-13 17:04
首先，产物是双自由基，所以过渡态也很可能有双自由基特征。除非你已经发现过渡态附近波函数是稳定的，否则 ...

感谢指点，请问如何寻找自旋破缺波函数呢

---

作者
Author: wzkchem5    时间: 2023-6-13 19:16

Iris_1 发表于 2023-6-13 10:23
感谢指点，请问如何寻找自旋破缺波函数呢

参见http://sobereva.com/82

---

欢迎光临 计算化学公社 (http://bbs.keinsci.com/)

Powered by Discuz! X3.3