计算化学公社

 找回密码 Forget password
 注册 Register
Views: 1236|回复 Reply: 0
打印 Print 上一主题 Last thread 下一主题 Next thread

[ORCA] 使用ORCA做DLPNO-NEVPT2计算后得到Singlet和Triplet能量高低翻转

[复制链接 Copy URL]

87

帖子

1

威望

3919

eV
积分
4026

Level 6 (一方通行)

本帖最后由 阿甘 于 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没收敛到所需的态,导致后面错了


本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-19 19:56 , Processed in 0.145813 second(s), 20 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list