本帖最后由 gzu401 于 2024-11-25 10:56 编辑
各位同行们，正常的TD-DFT计算rISC (或ISC)已按照Orca手册成功，现在想计算基于STEOM-DLPNO-CCSD的rISC，不知道步骤和输入文件怎么编写，遂来求助。
手册中写道："Theinput is similar to those discussed above. Here ESD(ISC) should be used on themain input to indicate an InterSystem Crossing calculation, and the Hessians should be provided byISCISHESSIAN and ISCFSHESSIAN for the initial and final states, respectively. Please note that the geometryused in the input file should correspond to that of the FINAL state, specified through the ISCFSHESSIAN flag. The relevantmatrix elements can be computed using any method available in ORCA.As a simple example, one could compute theexcited singlet and ground triplet geometries and Hessians for anthracene usingTD-DFT. Then, compute the spin-orbit coupling (SOC) matrix elements for aspecific triplet spin-sublevel using the same method (see the details below),potentially employing methods like CASSCF, MRCI,STEOM-CCSD"
即：先得到GS，ES和三重态的结构和hessians，然后基于GS结构做STEOM-DLPNO-CCSD下的单线态和三重态的能量以用于DELE。
但是，手册中没有给出具体的输入例子，就有点无从下手了，于是产生了如下问题和想法，麻烦各位老师帮忙建议。
1. “GS，ES和三重态的结构和hessians”可以用TD-DFT方法得到，后直接使用？
2. 基于上述得到的GS结构，做STEOM-DLPNO-CCSD计算得到S1和T1的能量，是如下结果不呢？
STEOM-CCSDRESULTS (SINGLETS)
IROOT= 1: 0.109392 au 2.977 eV 24008.9 cm**-1
STEOM-CCSD RESULTS (TRIPLETS)
IROOT= 1: 0.089216 au 2.428 eV 19580.7 cm**-1
所以DELE for rISC be : 19580.7 - 24008.9 = -4,428.2 cm-1?
3. 最后，最重要的步骤计算rISC:
a)是用手册中提供的例子作为输入文件（TD-DFT），将上述TD-DFT得到的S1-T1结构和hessian， STEOM-DLPNO-CCSD得到的能量，写入如下输入文件中来计算？
%TDDFT
NROOTS 5
SROOT 1
TROOT 1
TROOTSSL 0 1 -1 #T1的三个亚能级
Dosoc
%ESD
ISCISHESS "t1.hess"
ISCFSHESS "s1.hess"
DELE xxx # 用STEOM-DLPNO-CCSD计算得到的
（应该是3个job，此处省略了）
b)还是直接用TD-DFT得到的S1-T1结构和hessian，按照STEOM-DLPNO-CCSD的输入文件计算？（这里不会写，尝试着写了如下输入文件，不知道对还是错）
!STEOM-DLPNO-CCSDDEF2-TZVP DEF2-TZVP/C DEF2/J TightSCF TightPNO ESD(ISC)
%maxcore 5000
%scf
MaxIter 7000
end
%mdci NRoots 3 #如何控制T1的3个亚能级呢？（直接按照手册中的TD-DFT方式写3个job？）
DoRootWise true
OThresh 0.001
VThresh 0.001
TCutPNOSingles 1e-11
NDAV 400
DoStoreSTEOM true
DoSimpleDens false
AddL2Term True
DTol 1e-5
DoTriplet true
STEOMSOC true
end
%ESD
ISCISHESS "t1-opt.hess"
ISCFSHESS "s1-opt.hess"
APPROXADEN TRUE # 这个 keyword会直接去自动获取DELE不呢？
USEJ TRUE
DOHT false
TEMP 298.15
%pal nprocs 24
end
* xyzfile s0-opt.xyz 0 1
最后，麻烦各位行家了，祝大家生活愉快，科研顺利。
| 
发表于 Post on 
收藏 Add to favorites
楼主 Author