计算化学公社

 找回密码 Forget password
 注册 Register
Views: 349|回复 Reply: 9
打印 Print 上一主题 Last thread 下一主题 Next thread

[Gaussian/gview] 求助Gaussian软件做锕系配合物全电子相对论的计算的设置问题

[复制链接 Copy URL]

62

帖子

0

威望

651

eV
积分
713

Level 4 (黑子)

老师们好,以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。

th(h20)2.gjf

6.2 KB, 下载次数 Times of downloads: 4

1万

帖子

0

威望

8991

eV
积分
20739

Level 6 (一方通行)

2#
发表于 Post on 2025-6-19 12:29:24 | 只看该作者 Only view this author
建议用sf-X2C。DKH早就过时了
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?hl=zh-CN&user=XW6C6eQAAAAJ&view_op=list_works&sortby=pubdate
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1034/1702.htm
本团队长期招收研究生,有意者可私信联系

98

帖子

1

威望

1932

eV
积分
2050

Level 5 (御坂)

3#
发表于 Post on 2025-6-19 14:10:36 | 只看该作者 Only view this author
1. sf-x2c更好
2. 没有soc,只有标量相对论
3. 二分量重收缩的基组和考虑soc没直接关系,你不用二分量重收缩的基组也能考虑soc。使用更合适基组只是改善计算结果
4. 都是二分量,何来的四分量
ps. 你对相对论量化的理解有些问题。软件的话改用ORCA可能更好一些

6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

4#
发表于 Post on 2025-6-19 15:04:11 | 只看该作者 Only view this author
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软件中锕系配合物全电子相对论的计算” 改了,以后务必注意!
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

62

帖子

0

威望

651

eV
积分
713

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2025-6-20 09:46:35 | 只看该作者 Only view this author
本帖最后由 wangzh 于 2025-6-20 11:11 编辑
wzkchem5 发表于 2025-6-19 12:29
建议用sf-X2C。DKH早就过时了

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

1万

帖子

0

威望

8991

eV
积分
20739

Level 6 (一方通行)

6#
发表于 Post on 2025-6-20 12:35:12 | 只看该作者 Only view this author
wangzh 发表于 2025-6-20 09:46
好的,请问老师如果用sf-X2C的话,用Gaussian软件,对钍有推荐的基组吗?因为与X2C配套的x2c ...

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

评分 Rate

参与人数
Participants 1
eV +1 收起 理由
Reason
zjxitcc + 1 PySCF亦可

查看全部评分 View all ratings

Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?hl=zh-CN&user=XW6C6eQAAAAJ&view_op=list_works&sortby=pubdate
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1034/1702.htm
本团队长期招收研究生,有意者可私信联系

62

帖子

0

威望

651

eV
积分
713

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2025-6-20 15:10:05 | 只看该作者 Only view this author
wzkchem5 发表于 2025-6-20 12:35
应当改用ORCA或BDF。Gaussian目前还不支持sf-X2C

好的,感谢老师

62

帖子

0

威望

651

eV
积分
713

Level 4 (黑子)

8#
 楼主 Author| 发表于 Post on 2025-6-20 22:42:25 | 只看该作者 Only view this author
sobereva 发表于 2025-6-19 15:04
1 合理。虽然X2C更好,但DKH2的精度也能够满足一般需要

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

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

6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

9#
发表于 Post on 2025-6-21 04:17:04 | 只看该作者 Only view this author
wangzh 发表于 2025-6-20 22:42
感谢老师您的指点和修改,我以后会注意的,我这次先报名了七月份的中级量子化学班,等高级班开了我会报名 ...

效仿北京科音高级量子化学培训班(http://www.keinsci.com/KAQC)这页幻灯片里的例子的写法

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

62

帖子

0

威望

651

eV
积分
713

Level 4 (黑子)

10#
 楼主 Author| 发表于 Post on 2025-6-23 09:08:56 | 只看该作者 Only view this author
sobereva 发表于 2025-6-21 04:17
效仿北京科音高级量子化学培训班(http://www.keinsci.com/KAQC)这页幻灯片里的例子的写法

非常感谢老师

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-14 13:19 , Processed in 0.168000 second(s), 24 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list