计算化学公社

 找回密码 Forget password
 注册 Register
楼主 Author: sobereva
打印 Print 上一主题 Last thread 下一主题 Next thread

[综合交流] 用NWChem做SODFT在DFT计算中考虑旋轨耦合效应

[复制链接 Copy URL]

729

帖子

21

威望

5233

eV
积分
6382

Level 6 (一方通行)

16#
发表于 Post on 2024-8-18 12:01:29 | 只看该作者 Only view this author
本帖最后由 beefly 于 2024-8-18 12:04 编辑
冰释之川 发表于 2017-4-18 13:01
问题解决了,是Linear dependencies的锅:

Unfortunately, the spin-orbit code does not correctly h ...

nwchem判断基组线性依赖的默认阈值是1.0e-5。由于sodft不支持消除冗余积函数的算法,这么高的阈值显然不适用。可以用手册里没有讲的一个命令降低阈值。根据出错信息里的“Smallest S eigenvalue”,假如此值在1.0e-8以上,可以用

set lindep:tol 1.0e-8

但是如果出现更低的特征值,降低阈值的做法就不保险了,可能会带来其他问题,比如scf不收敛。

125

帖子

0

威望

468

eV
积分
593

Level 4 (黑子)

17#
发表于 Post on 2025-3-12 22:11:58 | 只看该作者 Only view this author
sob老师您好!我想用NWChem考虑锕系元素旋轨耦合,请教一下Stuttgart赝势的官网上蓝色的这一部分(图1)是旋轨势吗?手动加旋轨势是不是加这一部分?

QQ20250312-220831.png (77.46 KB, 下载次数 Times of downloads: 5)

图1

图1

6万

帖子

99

威望

5万

eV
积分
120148

管理员

公社社长

18#
 楼主 Author| 发表于 Post on 2025-3-13 04:33:32 | 只看该作者 Only view this author
Xx_DerEnTh 发表于 2025-3-12 22:11
sob老师您好!我想用NWChem考虑锕系元素旋轨耦合,请教一下Stuttgart赝势的官网上蓝色的这一部分(图1)是 ...

这是赝势而非旋轨势
北京科音自然科学研究中心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

125

帖子

0

威望

468

eV
积分
593

Level 4 (黑子)

19#
发表于 Post on 2025-3-13 09:54:30 | 只看该作者 Only view this author
sobereva 发表于 2025-3-13 04:33
这是赝势而非旋轨势

感谢sob老师回复,请问这种没给旋轨势的情况是不是就不能用NWChem分析了?

6万

帖子

99

威望

5万

eV
积分
120148

管理员

公社社长

20#
 楼主 Author| 发表于 Post on 2025-3-14 03:51:14 | 只看该作者 Only view this author
Xx_DerEnTh 发表于 2025-3-13 09:54
感谢sob老师回复,请问这种没给旋轨势的情况是不是就不能用NWChem分析了?

不能用。对于基态DFT计算,此时需要做严格的二分量相对论计算体现旋轨耦合效应
北京科音自然科学研究中心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

125

帖子

0

威望

468

eV
积分
593

Level 4 (黑子)

21#
发表于 Post on 2025-3-14 09:02:40 | 只看该作者 Only view this author
sobereva 发表于 2025-3-14 03:51
不能用。对于基态DFT计算,此时需要做严格的二分量相对论计算体现旋轨耦合效应

好的,谢谢sob老师!

125

帖子

0

威望

468

eV
积分
593

Level 4 (黑子)

22#
发表于 Post on 2025-3-28 15:49:42 | 只看该作者 Only view this author
sobereva 发表于 2025-3-14 03:51
不能用。对于基态DFT计算,此时需要做严格的二分量相对论计算体现旋轨耦合效应
sob老师您好!对于这篇博文算TlCl键解离能的例子,我如果用Dirac的X2C二分量哈密顿考虑耦合,单独算Cl的时候为了保证能量可相加减,是不是同样也得用X2C?如果是这样的话,对像Cl原子这样通常不需要考虑相对论的原子,用X2C这样处理相对论的方法会不会造成什么负面影响?

6万

帖子

99

威望

5万

eV
积分
120148

管理员

公社社长

23#
 楼主 Author| 发表于 Post on 2025-3-28 16:30:13 | 只看该作者 Only view this author
Xx_DerEnTh 发表于 2025-3-28 15:49
sob老师您好!对于这篇博文算TlCl键解离能的例子,我如果用Dirac的X2C二分量哈密顿考虑耦合,单独算Cl的时 ...



不需要特意考虑不代表不能考虑,考虑了当然从精度来说有益,代价就是多花时间而已
北京科音自然科学研究中心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

125

帖子

0

威望

468

eV
积分
593

Level 4 (黑子)

24#
发表于 Post on 2025-3-28 17:05:35 | 只看该作者 Only view this author
sobereva 发表于 2025-3-28 16:30


不需要特意考虑不代表不能考虑,考虑了当然从精度来说有益,代价就是多花时间而已

好的,谢谢sob老师!

125

帖子

0

威望

468

eV
积分
593

Level 4 (黑子)

25#
发表于 Post on 2025-3-29 21:18:43 | 只看该作者 Only view this author
sobereva 发表于 2025-3-28 16:30


不需要特意考虑不代表不能考虑,考虑了当然从精度来说有益,代价就是多花时间而已

sob老师您好!想再请教一下Dirac输出的电子能量具体是在哪里读?是下面这个吗?
Sum of all contributions to the energy
   Total energy                             :    -31332.519464472796
感觉数值好大

6万

帖子

99

威望

5万

eV
积分
120148

管理员

公社社长

26#
 楼主 Author| 发表于 Post on 2025-3-29 23:09:19 | 只看该作者 Only view this author
Xx_DerEnTh 发表于 2025-3-29 21:18
sob老师您好!想再请教一下Dirac输出的电子能量具体是在哪里读?是下面这个吗?
Sum of all contributio ...

内核电子能量很低,数量级自然很大
北京科音自然科学研究中心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

125

帖子

0

威望

468

eV
积分
593

Level 4 (黑子)

27#
发表于 Post on 2025-3-30 18:49:38 | 只看该作者 Only view this author
sobereva 发表于 2025-3-29 23:09
内核电子能量很低,数量级自然很大

sob老师,电子能量是在这里(Total energy)读对吗?

6万

帖子

99

威望

5万

eV
积分
120148

管理员

公社社长

28#
 楼主 Author| 发表于 Post on 2025-3-30 19:36:11 | 只看该作者 Only view this author
Xx_DerEnTh 发表于 2025-3-30 18:49
sob老师,电子能量是在这里(Total energy)读对吗?

我忘了,如果只有这一处看起来是电子能量,那就只可能是这个
北京科音自然科学研究中心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

125

帖子

0

威望

468

eV
积分
593

Level 4 (黑子)

29#
发表于 Post on 2025-3-30 20:11:06 | 只看该作者 Only view this author
sobereva 发表于 2025-3-30 19:36
我忘了,如果只有这一处看起来是电子能量,那就只可能是这个

好的,谢谢sob老师

本版积分规则 Credits rule

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

GMT+8, 2025-8-16 02:12 , Processed in 0.170626 second(s), 22 queries , Gzip On.

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