计算化学公社

标题: SAPT2+(3)dMP2/aug-cc-pVTZPSI4下计算二重态报错求助 [打印本页]

作者
Author: younggood 时间: 2021-10-7 20:19
标题: SAPT2+(3)dMP2/aug-cc-pVTZPSI4下计算二重态报错求助
各位老师们，我用“金”级别的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级别吗？谢谢老师们？

作者
Author: sobereva 时间: 2021-10-7 22:17
只有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）

作者
Author: yflchx 时间: 2021-10-8 10:00
检查一下二聚体结构、自洽场收敛。

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

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

作者
Author: younggood 时间: 2021-10-8 10:03

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

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

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

作者
Author: sobereva 时间: 2021-10-8 10:34

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

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

作者
Author: younggood 时间: 2021-10-8 15:06

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

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

作者
Author: 喵星大佬 时间: 2021-10-8 15:12
高精度方法指用比如CCSD(T)/CBS之类的级别片段相减得到的相互作用能

作者
Author: 喵星大佬 时间: 2021-10-8 15:15
当然了用LMO-EDA或者之类的方法显然也是可以的分解的

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