利用GAMESS计算旋轨耦合常数问题

xiaobogaga

paramecium86 发表于 2018-5-12 01:16
我用GAMESS在算 镧系金属化合物的旋轨耦合。基本思路是 用ZAPT/ROHF 做OPT和 生成分子Ground state的MO信息 ...

非常感谢！

xiaobogaga

paramecium86 发表于 2018-5-12 01:16
我用GAMESS在算 镧系金属化合物的旋轨耦合。基本思路是 用ZAPT/ROHF 做OPT和 生成分子Ground state的MO信息 ...

您好，在您输入文件的后一部分“----- RESULTS FROM SUCCESSFUL ROHF     GEOMETRY SEARCH -----
----- COORDS, ORBS, GRADIENT, AND APPROX. HESSIAN -----” 这些应该是来自于上一步OPT产生的文件吧。请问是OPT结束后的哪一个文件呢？OPT的输入文件可否看一下！万分感谢！！

paramecium86

xiaobogaga 发表于 2018-5-16 15:56
您好，在您输入文件的后一部分“----- RESULTS FROM SUCCESSFUL ROHF     GEOMETRY SEARCH -----
-----  ...

是的， state average mcscf 的轨道初猜来自于 ROHF 计算。我第一步是ROHF/ZAPT 级别的optimization. 然后单独用得到的 geometry 再重新单独算一次 modified virtual orbitals （MVOs）的轨道信息 作为MCSCF的初猜轨道信息。 所有轨道信息 是保留在后缀是 .dat 的文件里。 通常是存在GAMESS的scratch 文件夹里。 打开.dat 文件后 找到 $vec 开头的部分就是要用的MOs信息。如果是OPT运算， 每一步OPT cycle 就会生成一个$VEC， 一般就不止一个$vec。 只要找到最后一个就行。

MCSCF用MVOs，作为初猜
SO-MCQDPT用的轨道信息 是MCSCF生成的。 而且一定要用对应 .dat 文件里 Optimized MCSCF orbitals 而不要用 MCSCF生成的 Natural orbitals信息。

xiaobogaga

paramecium86 发表于 2018-5-17 03:45
是的， state average mcscf 的轨道初猜来自于 ROHF 计算。我第一步是ROHF/ZAPT 级别的optimization. 然 ...

非常感谢！！

qazwer

beefly 发表于 2018-5-11 17:03
如果是一个闭壳层基态/参考态 + n个单电子激发态这样的体系，可以用DFT/TDDFT做。然后在ADF，Dalton，或者G ...

您好，我的这篇帖子里说的情况感觉和您说的一般情况很类似
http://bbs.keinsci.com/thread-16387-1-1.html，在这种情况下，用什么软件，方法比较好呢？

qazwer

beefly 发表于 2018-5-11 17:03
如果是一个闭壳层基态/参考态 + n个单电子激发态这样的体系，可以用DFT/TDDFT做。然后在ADF，Dalton，或者G ...

您好，我的这篇帖子，http://bbs.keinsci.com/thread-16387-1-1.html，和您说的一般情况很像。
这种情况用什么软件和方法比较好呢？

猫头鹰

我的体系是自旋三重态基态，80个原子，想计算它的三重态激发态和单重态激发态之间的旋轨耦合，请问用哪个软件计算好呢？

wzkchem5

猫头鹰 发表于 2024-5-24 08:33
我的体系是自旋三重态基态，80个原子，想计算它的三重态激发态和单重态激发态之间的旋轨耦合，请问用哪个软 ...

如果这个体系不需要很大的活性空间就可以准确描述，那么可以做NEVPT2等多参考计算，这个很多程序可以做，ORCA、BDF、Molpro等。
如果多参考态方法算不动，而且这个体系不存在适合做参考态的单重态，那么可以用ORCA的DFT/ROCIS方法算。
如果多参考态方法算不动，但是这个体系存在适合做参考态的单重态，那么除了DFT/ROCIS以外，也可以用TDDFT算SOC，例如用ORCA、BDF、Dalton等

猫头鹰

wzkchem5 发表于 2024-5-24 15:40
如果这个体系不需要很大的活性空间就可以准确描述，那么可以做NEVPT2等多参考计算，这个很多程序可以做， ...

“适合做参考态的单重态”，这个指的是基态吗？我能不能这么理解：因为我的体系基态是自旋三重态，所以不存在适合做参考态的单重态，可以用orca的DFT/ROCIS算

wzkchem5

猫头鹰 发表于 2024-5-24 09:00
“适合做参考态的单重态”，这个指的是基态吗？我能不能这么理解：因为我的体系基态是自旋三重态，所以不 ...

不是，参考态不一定是基态。选择参考态的依据是用单参考态方法可以算得比较准，且基本没有自旋污染。
仔细阅读https://doi.org/10.1063/1.2047553了解参考态和基态的区别，以及什么时候需要选择一个不是基态的态作为参考态

猫头鹰

wzkchem5 发表于 2024-5-24 17:51
不是，参考态不一定是基态。选择参考态的依据是用单参考态方法可以算得比较准，且基本没有自旋污染。
仔 ...

好的，非常感谢老师的指点

猫头鹰

wzkchem5 发表于 2024-5-24 15:40
如果这个体系不需要很大的活性空间就可以准确描述，那么可以做NEVPT2等多参考计算，这个很多程序可以做， ...

老师您好，我用ROCIS计算SOC时出现了一个问题：SCF收敛之后就报错了，输出文件和报错信息如图所示，这个是内存不够还是ROCIS的关键词写的有问题呢？【基本信息：80个原子的分子团簇，基态是自旋三重态，想算第一三重态激发态和单重态之间的SOC】
输入文件：

报错信息：

wzkchem5

猫头鹰 发表于 2024-5-31 15:13
老师您好，我用ROCIS计算SOC时出现了一个问题：SCF收敛之后就报错了，输出文件和报错信息如图所示，这个 ...

maxcore太多了，maxcore（以MB为单位）乘以nprocs必须小于你机器的可用物理内存

猫头鹰

wzkchem5 发表于 2024-6-4 16:06
maxcore太多了，maxcore（以MB为单位）乘以nprocs必须小于你机器的可用物理内存

谢谢老师，我再调试一下。