求助MOlCAS 自旋轨道耦合(SOC）的计算

函数与激情

在最近投稿的一篇文章中，我使用Molcas手册中的模板（输入见下文）计算来MS-CASPT2级别的SOC,审稿人提问了一个问题，如下：
It iswritten that the SOC was calculated at the MSCASPT2 level using Openmolcas.Does this mean that the matrix elements of the SOC matrix are computed at theMSCASPT2 level? Or is the SOC calculated with SACASSCF and only the diagonalelements of the SOC matrix are replaced with MSCASPT2 energy and diagonalized?Please describe clearly.


我不太了解程序到底是怎么算的，之前粗浅的认识是第一种算法，MS-CASPT2直接计算得到SOC。如有错误，请论坛有熟悉Molcas程序的前辈不吝赐教，十分感谢！

1. *CASPT2/RASSI calculation on singlet-triplet gaps in acrolein
   
2. *File: CASPT2.S-T_gap.acrolein
   
3. *
   
4. &GATEWAY
   
5. Title= Acrolein molecule
   
6. coord = acrolein.xyz; basis = STO-3G; group= c1
   
7. &SEWARD
   
8. AMFI
   
9. &SCF
   
10. &RASSCF
    
11. Spin= 1; Nactel= 6 0 0; Inactive= 12; Ras2= 5
    
12. CiRoot= 1 1 1
    
13. &CASPT2
    
14. Frozen= 4
    
15. MultiState= 1 1
    
16. >>COPY $Project.JobMix JOB001
    
17. &RASSCF
    
18. LumOrb
    
19. Spin= 3; Nactel= 6 0 0; Inactive= 12; Ras2= 5
    
20. CiRoot= 5 5 1
    
21. &CASPT2
    
22. Frozen= 4
    
23. MultiState= 5 1 2 3 4 5
    
24. >>COPY $Project.JobMix JOB002
    
25. &RASSI
    
26. Nr of JobIph= 2 1 5; 1; 1 2 3 4 5
    
27. Spin
    
28. EJob

复制代码

beefly

参照这篇文章第二页左侧理论部分的写法
Journal of Quantitative Spectroscopy & Radiative Transfer 288 (2022) 108247
不过你这个STO-3G基组，说明计算毫无意义

scf

beefly 发表于 2022-6-12 07:55
参照这篇文章第二页左侧理论部分的写法
Journal of Quantitative Spectroscopy & Radiative Transfer 288  ...

问下,这个体系MRCI+Q算得动吗?会比CASPT2更准吗?

函数与激情

beefly 发表于 2022-6-12 07:55
参照这篇文章第二页左侧理论部分的写法
Journal of Quantitative Spectroscopy & Radiative Transfer 288  ...

已看到，十分感谢您的帮助！

beefly

scf 发表于 2022-6-12 08:04
问下,这个体系MRCI+Q算得动吗?会比CASPT2更准吗?

取决于MRCI的算法，体系大小，以及活性轨道个数。效率最高的Molpro一般能处理10个以内的CAS轨道，30个以内的关联电子和300个以下的基函数，如果用了RI技术能算的更大

scf

Xian CI的MRCI据说不用RDM，https://www.tandfonline.com/doi/abs/10.1080/00268976.2018.1441464
相比molpro效率如何？