计算化学公社

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

[CP2K] CP2K指定占据的问题

[复制链接 Copy URL]

72

帖子

0

威望

297

eV
积分
369

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
各位CP2K专家,
  大家好!我最近看了卢天老师讲义,了解到cp2k可以采用dft+u及指定占据方法来处理一些磁性体系。
我现在打算用这个方法计算一下NiO体系。向大家请教两个简单的问题:
(1)cp2k可以指定不同的d轨道占据,同时可以打印出占据数。这个占据数是否就是原子轨道的占据呢?
  (2)我看占据数有时候会超过1.0占据,为1.02这样的,其分析方法是mulliken。按理说不应该超过1.0的占据。这个是为什么呢?
  我的量子化学基础比较薄弱,问的不太专业,还请各位专家见谅!
  谢谢大家!

6万

帖子

99

威望

6万

eV
积分
125122

管理员

公社社长

8#
发表于 Post on 2023-11-5 13:02:05 | 只看该作者 Only view this author
lujun 发表于 2023-11-5 12:55
谢谢卢天老师的指导!我按照这个帖子做了一下:
我测试的是Sc原子,采用DZVP的MOLOPT基组,使用Lowdin方 ...

正确
北京科音自然科学研究中心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

72

帖子

0

威望

297

eV
积分
369

Level 3 能力者

7#
 楼主 Author| 发表于 Post on 2023-11-5 12:55:38 | 只看该作者 Only view this author
sobereva 发表于 2023-11-5 11:56
不是

对应关系按下文判断

谢谢卢天老师的指导!我按照这个帖子做了一下:
我测试的是Sc原子,采用DZVP的MOLOPT基组,使用Lowdin方法分析了基函数上的占据数情况:
1S:2.00
2S:1.84
3S:0.00
1P:5.99
2P:0.00
1D:1.147
2D:0.003
Sc原子赝化后的电子组态是3s2 3p6 4s2 3d1,根据基函数布居分析结果以及Sc原子的电子组态是不是可以得出以下结论:
1S对应的是3s原子轨道,2S+3S对应的是4s轨道,1P+2P对应的是3p轨道,1D+2D对应的是3d轨道
不知道我的理解是否正确,还请卢天老师批评指正!

6万

帖子

99

威望

6万

eV
积分
125122

管理员

公社社长

6#
发表于 Post on 2023-11-5 11:56:43 | 只看该作者 Only view this author
lujun 发表于 2023-11-4 22:11
谢谢卢天老师!MOLOPT的基组的基函数和原子轨道是怎么对应的哇?MOLOPT的基函数是不是就是原子轨道?

不是

对应关系按下文判断
利用布居分析判断基函数与原子轨道的对应关系
http://sobereva.com/418http://bbs.keinsci.com/thread-10094-1-1.html
北京科音自然科学研究中心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

72

帖子

0

威望

297

eV
积分
369

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2023-11-4 22:11:14 | 只看该作者 Only view this author
sobereva 发表于 2023-9-18 10:30
根据基函数与原子轨道的对应描述关系,可以从前者的占据数判断后者的占据数

谢谢卢天老师!MOLOPT的基组的基函数和原子轨道是怎么对应的哇?MOLOPT的基函数是不是就是原子轨道?

6万

帖子

99

威望

6万

eV
积分
125122

管理员

公社社长

4#
发表于 Post on 2023-9-18 10:30:32 | 只看该作者 Only view this author
lujun 发表于 2023-9-17 14:02
谢谢卢天老师!
我看cp2k讲义的"DFT+U的计算"章节的45页上,列出了打印的UO2的f shell上的电子占据信息 ...

根据基函数与原子轨道的对应描述关系,可以从前者的占据数判断后者的占据数
北京科音自然科学研究中心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

72

帖子

0

威望

297

eV
积分
369

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2023-9-17 14:02:04 | 只看该作者 Only view this author
本帖最后由 lujun 于 2023-9-17 14:06 编辑
sobereva 发表于 2023-9-16 16:15
1 我不知道你说的打印是哪里、具体怎么打印的
Multiwfn基于CP2K的molden文件做Mulliken布居分析倒是可以给 ...

谢谢卢天老师!
我看cp2k讲义的"DFT+U的计算"章节的45页上,列出了打印的UO2的f shell上的电子占据信息,使用的ENFORCE_OCCUPATION方法指定占据的。
我的困惑是实际上程序打印的其实只是基函数的占据信息,而不是原子轨道的占据的。但是讲义里面给出的参考文献Phys. Scr. 90, 094014(2015)说的是5f原子轨道上的占据。这是不是有问题呢?

6万

帖子

99

威望

6万

eV
积分
125122

管理员

公社社长

2#
发表于 Post on 2023-9-16 16:15:33 | 只看该作者 Only view this author
1 我不知道你说的打印是哪里、具体怎么打印的
Multiwfn基于CP2K的molden文件做Mulliken布居分析倒是可以给你最终状态的各个原子轨道的占据数信息

2 Mulliken分析自身的理论层面的问题,得到的原子轨道对分子轨道的贡献可能超过100%,下文有充分的讨论
《分子轨道成分的计算》(http://sioc-journal.cn/Jwk_hxxb/CN/abstract/abstract340458.shtml
北京科音自然科学研究中心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

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

GMT+8, 2026-2-17 21:50 , Processed in 0.175858 second(s), 21 queries , Gzip On.

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