SAPT2+(3)dMP2/aug-cc-pVTZPSI4下计算二重态报错求助

younggood

各位老师们，我用“金”级别的SAPT2+(3)δMP2/aug-cc-pVTZ做EDA，其中有一个片段是0，2，我也在输入文件里加入了set reference uhf，但是有如下报错： Only SAPT0 supports a reference different from "reference rhf". !
!                                                                  !
!------------------------------------------------------------------!

    Psi4 stopped on: Thursday, 07 October 2021 07:53PM
    Psi4 wall time for execution: 0:00:00.07

*** Psi4 encountered an error. Buy a developer more coffee!
*** Resources and help at github.com/psi4/psi4.

然后我就尝试了用sSAPT0/jun-cc-pVDZ，但是能量分解的结果显示如下：

Electrostatics                -45.70299450 [mEh]     -28.67906203 [kcal/mol]    -119.99319552 [kJ/mol]
      Elst10,r                    -45.70299450 [mEh]     -28.67906203 [kcal/mol]    -119.99319552 [kJ/mol]

    Exchange                      128.96926342 [mEh]      80.92943462 [kcal/mol]     338.60875445 [kJ/mol]
      Exch10                      128.96926342 [mEh]      80.92943462 [kcal/mol]     338.60875445 [kJ/mol]
      Exch10(S^2)                 119.69789009 [mEh]      75.11156003 [kcal/mol]     314.26676715 [kJ/mol]

    Induction                     -48.59268163 [mEh]     -30.49236808 [kcal/mol]    -127.58006804 [kJ/mol]
      Ind20,r                     -90.24619225 [mEh]     -56.63034061 [kcal/mol]    -236.94134510 [kJ/mol]
      Exch-Ind20,r                 73.03402884 [mEh]      45.82954501 [kcal/mol]     191.75081631 [kJ/mol]
      delta HF,r (2)              -31.38051823 [mEh]     -19.69157248 [kcal/mol]     -82.38953925 [kJ/mol]

    Dispersion                    -17.82670007 [mEh]     -11.18642318 [kcal/mol]     -46.80399458 [kJ/mol]
      Disp20                      -23.66413869 [mEh]     -14.84947122 [kcal/mol]     -62.13018757 [kJ/mol]
      Exch-Disp20                   5.83743862 [mEh]       3.66304804 [kcal/mol]      15.32619299 [kJ/mol]

  Total HF                         34.67358729 [mEh]      21.75800451 [kcal/mol]      91.03549089 [kJ/mol]
  Total SAPT0                      16.84688722 [mEh]      10.57158134 [kcal/mol]      44.23149631 [kJ/mol]

这个相互作用能总和是10.57158134 kcal/mol，所以很疑惑，相互作用能总和不应该是负值吗，是不是我用的这个计算级别太低，导致结果不太好呢？另外，双重态片段是只能用sSAPT0/jun-cc-pVDZ级别吗？谢谢老师们？

sobereva

只有SAPT0才支持开壳层体系，因此没法用SAPT2+(3)dMP2

下文介绍的我的文章中，只有C18和O2做的SAPT不是SAPT2+(3)dMP2而是SAPT0，也是因为这个原因
全面探究18碳环独特的分子间相互作用与pi-pi堆积特征
http://sobereva.com/572（http://bbs.keinsci.com/thread-19660-1-1.html）

yflchx

检查一下二聚体结构、自洽场收敛。

如果画一条相互作用能量曲线，随结构的变化相互作用有排斥有吸引。

sSAPT0/jun-cc-pVDZ表现较好是因为误差相抵，说不定也会出例外。

younggood

sobereva 发表于 2021-10-7 22:17
只有SAPT0才支持开壳层体系，因此没法用SAPT2+(3)dMP2

下文介绍的我的文章中，只有C18和O2做的SAPT不是S ...

谢谢sob老师，另外，我用sSAPT0/jun-cc-pVDZ算出来的相互作用能总和是正值，请问这个是合理的吗？非常感谢！

sobereva

younggood 发表于 2021-10-8 10:03
谢谢sob老师，另外，我用sSAPT0/jun-cc-pVDZ算出来的相互作用能总和是正值，请问这个是合理的吗？非常感 ...

和高精度方法算的相互作用能对比，如果符号都不同显然不合理

younggood

sobereva 发表于 2021-10-8 10:34
和高精度方法算的相互作用能对比，如果符号都不同显然不合理

老师，因为刚接触EDA分析，还是不太理解，如果我这个体系用PSI4做EDA不能用高精度的SAPT2+(3)δMP2/aug-cc-pVTZ,那么怎么将这个正值的结果与高精度的方法做对比呢？谢谢老师！

喵星大佬

高精度方法指用比如CCSD(T)/CBS之类的级别片段相减得到的相互作用能

喵星大佬

当然了用LMO-EDA或者之类的方法显然也是可以的分解的