计算化学公社

标题: 求助Gaussian软件做锕系配合物全电子相对论的计算的设置问题 [打印本页]

作者
Author: wangzh 时间: 2025-6-19 10:34
标题: 求助Gaussian软件做锕系配合物全电子相对论的计算的设置问题
老师们好，以Th(H2O2)4+为例，我想用Gaussian软件进行单点能计算，用PBE1PBE/SARC-DKH2/DKH-def2-TZVP int=dkh2全电子相对论，输入文件已经上传。我有些问题不太懂，想请教各位老师
1.我的输入文件所使用的方法是否正确或者合理，如果有错误或者不合理的地方请老师指点?
2.我的单点能计算想要考虑标量相对论和SOC,请问选择所使用的方法是否已经考虑了SOC?
3.Gaussian软件用全电子二分量相对论计算时，我是否可以认为int=dkh2和用给二分量的全电子基组考虑了标量相对论和SOC，而int=dkh2和用标量的全电子基组考虑的是标量相对论
4.这篇“使用Gaussian做镧系金属配合物的量子化学计算”文章中说，“此处DKHSO代表做DKH4标量相对论计算。x2c-TZVPall基组定义是从BSE上拷来的。”我是否可以认为int=DKHSO和x2c-TZVPall-2c是四分量的全电子相对论计算。考虑了标量相对论和SOC。

作者
Author: wzkchem5 时间: 2025-6-19 12:29
建议用sf-X2C。DKH早就过时了

作者
Author: wxsxjj 时间: 2025-6-19 14:10
1. sf-x2c更好
2. 没有soc，只有标量相对论
3. 二分量重收缩的基组和考虑soc没直接关系，你不用二分量重收缩的基组也能考虑soc。使用更合适基组只是改善计算结果
4. 都是二分量，何来的四分量
ps. 你对相对论量化的理解有些问题。软件的话改用ORCA可能更好一些

作者
Author: sobereva 时间: 2025-6-19 15:04
1 合理。虽然X2C更好，但DKH2的精度也能够满足一般需要

2、3 你当前输入文件的写法做的是spin-free的DKH2，只能考虑标量相对论效应。若要Gaussian做全电子二分量相对论计算同时考虑标量和旋轨耦合效应，需要int=DKHSO结合GKS，并建议结合x2c的-2c版本基组计算

4 int=DKHSO结合标准的KS-DFT，只考虑了DKH4的标量部分

另外，你写了guess=alter但没写轨道序号，毫无意义。
引用博文时应当给出具体链接，给你加上了。

PS：北京科音高级量子化学培训班（http://www.keinsci.com/KAQC）里专门有一节“相对论量子化学”，100多页幻灯片详细讲解了相对论量子化学的背景知识和计算应用，系统学过一遍后对于相对论计算的各种细节、程序用法就全明白了。

如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示，求助帖必须在帖子标题明确体现出此帖内容是求助或提问，避免有任何歧义和含糊性，仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题 “Gaussian软件中锕系配合物全电子相对论的计算” 改了，以后务必注意！

作者
Author: wangzh 时间: 2025-6-20 09:46
本帖最后由 wangzh 于 2025-6-20 11:11 编辑

wzkchem5 发表于 2025-6-19 12:29
建议用sf-X2C。DKH早就过时了

好的，请问老师如果用sf-X2C的话，用Gaussian软件，对钍有推荐的基组吗？因为与X2C配套的x2c-系列基组只到Rn。

作者
Author: wzkchem5 时间: 2025-6-20 12:35

wangzh 发表于 2025-6-20 09:46
好的，请问老师如果用sf-X2C的话，用Gaussian软件，对钍有推荐的基组吗？因为与X2C配套的x2c ...

应当改用ORCA或BDF。Gaussian目前还不支持sf-X2C

作者
Author: wangzh 时间: 2025-6-20 15:10

wzkchem5 发表于 2025-6-20 12:35
应当改用ORCA或BDF。Gaussian目前还不支持sf-X2C

好的，感谢老师

作者
Author: wangzh 时间: 2025-6-20 22:42

sobereva 发表于 2025-6-19 15:04
1 合理。虽然X2C更好，但DKH2的精度也能够满足一般需要

2、3 你当前输入文件的写法做的是spin-free的DKH ...

感谢老师您的指点和修改，我以后会注意的，我这次先报名了七月份的中级量子化学班，等高级班开了我会报名参加。
老师您这里说的同时考虑标量和旋轨耦合效应，需要int=DKHSO结合GKS，请问老师GKS应该怎么考虑，或者说这里具体输入形式应该是怎样的？感谢老师

作者
Author: sobereva 时间: 2025-6-21 04:17

wangzh 发表于 2025-6-20 22:42
感谢老师您的指点和修改，我以后会注意的，我这次先报名了七月份的中级量子化学班，等高级班开了我会报名 ...

效仿北京科音高级量子化学培训班（http://www.keinsci.com/KAQC）这页幻灯片里的例子的写法
(, 下载次数 Times of downloads: 28)

作者
Author: wangzh 时间: 2025-6-23 09:08

sobereva 发表于 2025-6-21 04:17
效仿北京科音高级量子化学培训班（http://www.keinsci.com/KAQC）这页幻灯片里的例子的写法

非常感谢老师

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