sobereva 发表于 2024-7-2 19:28
SCF收敛了，但不一定确保收敛到了稳定波函数。虚频肯定得消掉
如果你没那么多ORCA的使用经验就都用Gauss ...

感谢回复！
使用orca是因为体系原子数多200多个用orca优化更快
我是想得到色散作用和静电作用的比值以及色散在总相互作用的占比
谢谢sob老师指导！

zmeeyhnn 发表于 2024-7-2 09:53
谢谢sob老师！
结构是在B97-3c下优化好的，out文件中最后的SCF CONVERGED是下面这样的，请问这样是不是 ...

SCF收敛了，但不一定确保收敛到了稳定波函数。虚频肯定得消掉
如果你没那么多ORCA的使用经验就都用Gaussian就完了

“色散相互作用进行定量分析” 语义不明
二者的差异在此文里明确说了
使用sobEDA和sobEDAw方法做非常准确、快速、方便、普适的能量分解分析
http://sobereva.com/685（http://bbs.keinsci.com/thread-39446-1-1.html）
sobEDAw只适合用于弱相互作用

sobereva 发表于 2024-7-1 17:49
要比就在相同级别下比

如果结构已经经过合理优化，相互作用能为那么大值明显不正常。检查SCF收敛情况。 ...

谢谢sob老师！
结构是在B97-3c下优化好的，out文件中最后的SCF CONVERGED是下面这样的，请问这样是不是不正常的收敛？这个结构刚刚计算出的频率也是有两个虚频。。是不是要在verytightSCF再次优化？
想对色散相互作用进行定量分析是不是应该用sobEDAw而不是sobEDA？
out文件已上传

---

--------------------------------------------------------------------------------------------
   Iter.        energy            ||Error||_2        Shift      TRadius  Mac/Mic        Rej.
--------------------------------------------------------------------------------------------
WARNING: 3 diagonal Hessian elements are negative!
    0      -6104.664219525172     2.655487e-03                     0.400 (TRAH MAcro)   No
WARNING : negative HOMO - LUMO gap : -0.007899
    0     dE    -3.574273e-07     2.856620e-03    -7.1486e-07      0.000 (TRAH MIcro)
    0     dE    -5.109495e-07     3.665924e-04    -1.0219e-06      0.000 (TRAH MIcro)
    0     dE    -5.204952e-07     7.667045e-05    -1.0410e-06      0.000 (TRAH MIcro)
WARNING: 3 diagonal Hessian elements are negative!
    1      -6104.664220031228     7.712922e-05                           (NR   MAcro)
WARNING : negative HOMO - LUMO gap : -0.007897
    1     dE    -3.623618e-10     3.579171e-05                           (NR   MIcro)
    1     dE    -4.254563e-10     2.060147e-05                           (NR   MIcro)
    1     dE    -4.381144e-10     1.044774e-05                           (NR   MIcro)
    1     dE    -4.420205e-10     3.876396e-06                           (NR   MIcro)
    1     dE    -4.427767e-10     1.511761e-06                           (NR   MIcro)
    1     dE    -4.429818e-10     8.936825e-07                           (NR   MIcro)
WARNING: 3 diagonal Hessian elements are negative!
    2      -6104.664220031331     1.414020e-06                           (NR   MAcro)
WARNING : negative HOMO - LUMO gap : -0.007897

               *****************************************************
               *                     SUCCESS                       *
               *           SCF CONVERGED AFTER  32 CYCLES          *
               *****************************************************

Total Energy       :        -6104.66422003 Eh         -166116.35862 eV
  Last Energy change         ...   -1.0277e-10  Tolerance :   1.0000e-08
  Last Orbital Gradient      ...    1.4140e-06  Tolerance :   1.0000e-05
  Last Orbital Rotation      ...    1.7387e-05

要比就在相同级别下比

如果结构已经经过合理优化，相互作用能为那么大值明显不正常。检查SCF收敛情况。是否、怎么考虑BSSE是相对次要的事

sobEDA牵扯到Cs时用def2-TZVP没问题。甭用LanL2DZ那么烂的基组