计算化学公社

标题: 利用GAMESS计算旋轨耦合常数问题 [打印本页]

作者
Author: xiaobogaga 时间: 2018-5-10 18:56
标题: 利用GAMESS计算旋轨耦合常数问题
想请教大家一下，怎样用GAMESS计算旋轨耦合常数（SOC），有没有正在做类似研究的同学，希望能给一个输入文件看一看，参考一下，都有哪些关键词来控制，怎么分析输出文件？谢谢大家！

作者
Author: 西柚西柚啊 时间: 2018-5-10 19:54
你好，可以用dalton计算旋轨耦合的。

作者
Author: 我本是个娃娃 时间: 2018-5-10 20:36
免费软件考虑一下http://sobereva.com/411

作者
Author: xiaobogaga 时间: 2018-5-10 20:47

我本是个娃娃发表于 2018-5-10 20:36
免费软件考虑一下http://sobereva.com/411

我看过这个了，也安装编译了PySOC，从帖子来看，对于金属参与的大体系来说，这个方法很好，可是我的体系很小，一般只有5-6个原子，而且重原子也不多。如果用TDDFT+PySOC来算的话，怕是误差太大。所以打算用GAMESS 或者Molpro的 CAS来计算，得到高精度结果。请问你有了解吗？

作者
Author: 我本是个娃娃 时间: 2018-5-10 22:50

xiaobogaga 发表于 2018-5-10 20:47
我看过这个了，也安装编译了PySOC，从帖子来看，对于金属参与的大体系来说，这个方法很好，可是我的体系 ...

不甚了解，推荐你到时候参加一下量化高级班，我也会到时候仔细学习

作者
Author: sobereva 时间: 2018-5-11 08:31
不如用免费的ORCA，手册里就有例子，用着也省事

作者
Author: xiaobogaga 时间: 2018-5-11 08:56

西柚西柚啊发表于 2018-5-10 19:54
你好，可以用dalton计算旋轨耦合的。

服务器上还没有安装编译dalton软件，请问您有dalton计算SOC的输入文件吗？参考一下，谢谢！

作者
Author: xiaobogaga 时间: 2018-5-11 08:56

sobereva 发表于 2018-5-11 08:31
不如用免费的ORCA，手册里就有例子，用着也省事

谢谢sob老师，我这就去了解一下！

作者
Author: sobereva 时间: 2018-5-11 09:19

xiaobogaga 发表于 2018-5-11 08:56
服务器上还没有安装编译dalton软件，请问您有dalton计算SOC的输入文件吗？参考一下，谢谢！

用Dalton算SOC的人大多都是在TDDFT级别算，怎么操作在本版dalton分类里我有过充分回复
CASSCF/多参考级别没必要用Dalton

作者
Author: xiaobogaga 时间: 2018-5-11 09:24

sobereva 发表于 2018-5-11 09:19
用Dalton算SOC的人大多都是在TDDFT级别算，怎么操作在本版dalton分类里我有过充分回复
CASSCF/多参考级 ...

谢谢sob老师，那我准备用ORCA计算SOC了，在看它的手册，还没找到算例

作者
Author: beefly 时间: 2018-5-11 17:03
本帖最后由 beefly 于 2018-5-11 17:08 编辑

如果是一个闭壳层基态/参考态 + n个单电子激发态这样的体系，可以用DFT/TDDFT做。然后在ADF，Dalton，或者Gaussian+PySOC里面选一个。

一般情况：mcscf+soc或mrci+soc。gamess虽然能做，但是输入很罗嗦，不如molpro，molcas，orca简单方便

作者
Author: xiaobogaga 时间: 2018-5-11 17:07

beefly 发表于 2018-5-11 17:03
如果是一个闭壳层基态/参考态 + n个单电子激发态这样的体系，可以用TDDFT甚至EOM-CC这类方法去做。然后在AD ...

我是准备计算单-三重态MECP处的旋轨耦合常数。用sobMECP计算得到了MECP，接下来计算SOC的时候，不知道怎么算了。。。一开始想用GAMESS，但是手册实在看不懂，然后准备用ORCA或者Molpro,上午安装好了ORCA，实验室买的也有MolPro，现在又卡在输入文件上了。。。。

作者
Author: rayinyin 时间: 2018-5-11 17:14
ORCA的CASSCF或MRCI都能做SOC，他们自己说是在QDPT下作的(准简并微扰理论)，其实就是考虑SOC的CI。
对应的MOLCAS里的RASSI也能做。
MOLPRO肯定也没问题，只是我只用过ORCA和MOLCAS

作者
Author: xiaobogaga 时间: 2018-5-11 17:19

rayinyin 发表于 2018-5-11 17:14
ORCA的CASSCF或MRCI都能做SOC，他们自己说是在QDPT下作的(准简并微扰理论)，其实就是考虑SOC的CI。
对应的 ...

您还有ORCA的输入文件吗？

作者
Author: paramecium86 时间: 2018-5-12 01:16
我用GAMESS在算镧系金属化合物的旋轨耦合。基本思路是用ZAPT/ROHF 做OPT和生成分子Ground state的MO信息。以这个为基础计算State-Average MCSCF, 然后再进一步计算SO-MCQDPT。附件里放了CeO 分子的 SA-MCSCF和SO-MCQDPT的input文件.

作者
Author: xiaobogaga 时间: 2018-5-15 08:47

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

非常感谢！

作者
Author: xiaobogaga 时间: 2018-5-16 15:56

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的输入文件可否看一下！万分感谢！！

作者
Author: paramecium86 时间: 2018-5-17 03:45

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信息。

作者
Author: xiaobogaga 时间: 2018-5-17 14:54

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

非常感谢！！

作者
Author: qazwer 时间: 2020-2-24 23:32

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

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

作者
Author: qazwer 时间: 2020-2-24 23:36

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

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

作者
Author: 猫头鹰 时间: 2024-5-24 15:33
我的体系是自旋三重态基态，80个原子，想计算它的三重态激发态和单重态激发态之间的旋轨耦合，请问用哪个软件计算好呢？

作者
Author: wzkchem5 时间: 2024-5-24 15:40

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

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

作者
Author: 猫头鹰 时间: 2024-5-24 16:00

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

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

作者
Author: wzkchem5 时间: 2024-5-24 17:51

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

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

作者
Author: 猫头鹰 时间: 2024-5-31 15:07

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

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

作者
Author: 猫头鹰 时间: 2024-5-31 22:13

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

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

(, 下载次数 Times of downloads: 29)

作者
Author: wzkchem5 时间: 2024-6-4 16:06

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

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

作者
Author: 猫头鹰 时间: 2024-6-7 17:21

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

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)