计算化学公社

标题: 求助MOlCAS 自旋轨道耦合(SOC）的计算 [打印本页]
作者
Author: 函数与激情    时间: 2022-6-11 20:31
标题: 求助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
复制代码




作者
Author: beefly    时间: 2022-6-12 07:55
本帖最后由 beefly 于 2022-6-12 07:58 编辑

参照这篇文章第二页左侧理论部分的写法
Journal of Quantitative Spectroscopy & Radiative Transfer 288 (2022) 108247
不过你这个STO-3G基组，说明计算毫无意义
作者
Author: scf    时间: 2022-6-12 08:04
beefly 发表于 2022-6-12 07:55
参照这篇文章第二页左侧理论部分的写法
Journal of Quantitative Spectroscopy & Radiative Transfer 288  ...

问下,这个体系MRCI+Q算得动吗?会比CASPT2更准吗?
作者
Author: 函数与激情    时间: 2022-6-12 08:28
beefly 发表于 2022-6-12 07:55
参照这篇文章第二页左侧理论部分的写法
Journal of Quantitative Spectroscopy & Radiative Transfer 288  ...

已看到，十分感谢您的帮助！
作者
Author: beefly    时间: 2022-6-12 18:21
本帖最后由 beefly 于 2022-6-12 18:22 编辑
scf 发表于 2022-6-12 08:04
问下,这个体系MRCI+Q算得动吗?会比CASPT2更准吗?

取决于MRCI的算法，体系大小，以及活性轨道个数。效率最高的Molpro一般能处理10个以内的CAS轨道，30个以内的关联电子和300个以下的基函数，如果用了RI技术能算的更大
作者
Author: scf    时间: 2022-6-23 11:48
Xian CI的MRCI据说不用RDM，https://www.tandfonline.com/doi/abs/10.1080/00268976.2018.1441464
相比molpro效率如何？



欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3