计算化学公社

标题: 使用ORCA做DLPNO-NEVPT2计算后得到Singlet和Triplet能量高低翻转 [打印本页]

作者
Author: 阿甘 时间: 2023-3-9 10:59
标题: 使用ORCA做DLPNO-NEVPT2计算后得到Singlet和Triplet能量高低翻转
本帖最后由阿甘于 2023-3-9 11:39 编辑

各位老师好，学生最近在使用ORCA 5.0.3做DLPNO-NEVPT2计算有机物（体系共90原子，全是CHNO原子）的exchange coupling时发现，在CASSCF时三重态是基态，而加了DLPNO-NEVPT2之后，单重态却是基态了，不知如何理解这个结果。是不是因为活性空间没取够的原因呢？烦请各位老师帮学生看看，附上CASSCF和DLPNO-NEVPT2的输出结果和输入行：

输出结果: CASSCF:
CAS-SCF STATES FOR BLOCK  1 MULT= 3 NROOTS= 1
---------------------------------------------

ROOT 0:  E= -2760.1609496896 Eh
   1.00000 [    0]: 11

---------------------------------------------
CAS-SCF STATES FOR BLOCK  2 MULT= 1 NROOTS= 1
---------------------------------------------

ROOT 0:  E= -2760.1606329749 Eh
   0.56896 [    0]: 20
   0.43104 [    2]: 02

-----------------------------
SA-CASSCF TRANSITION ENERGIES
------------------------------

LOWEST ROOT (ROOT 0 ,MULT 3) =  -2760.160949690 Eh -75107.798 eV

STATE ROOT MULT  DE/a.u. DE/eV DE/cm**-1
1: 0 1 0.000317    0.009    69.5

DLPNO-NEVPT2:
-----------------------------
NEVPT2 TRANSITION ENERGIES
------------------------------

LOWEST ROOT (ROOT 0, MULT 1) =  -2769.290202144 Eh -75356.217 eV

STATE ROOT MULT  DE/a.u.    DE/eV DE/cm**-1
1: 0 3 0.000246    0.007    54.0

输入行：
CASSCF:
%casscf nel 2
      norb 2
      nroots 1,1
      mult 3,1
      OrbStep SuperCI_PT
      SwitchStep DIIS
      SwitchConv 0.003
      SwitchIter 30
      shiftup 2
      shiftdn 2
      end
* xyz 0 3
DLPNO-NEVPTS:
%casscf nel 2
      norb 2
      mult 3,1
      OrbStep SuperCI_PT
      SwitchStep DIIS
      SwitchConv 0.003
      SwitchIter 30
      shiftup 2
      shiftdn 2
      TrafoStep RI
      PTMethod DLPNO_NEVPT2
      PTSettings
      TCutPNO 1e-8
      MaxIter 20
      MaxDIIS 7
      D4Tpre 1e-14
      D4Step efficient
      end
      end
* xyz 0 3

刚好像发现了问题所在了，是MULT=1的CASSCF没收敛到所需的态，导致后面错了

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