计算化学公社

标题: CP2K指定占据的问题 [打印本页]

作者
Author: lujun 时间: 2023-9-16 16:10
标题: CP2K指定占据的问题
各位CP2K专家，
  大家好！我最近看了卢天老师讲义，了解到cp2k可以采用dft+u及指定占据方法来处理一些磁性体系。
我现在打算用这个方法计算一下NiO体系。向大家请教两个简单的问题：
（1）cp2k可以指定不同的d轨道占据，同时可以打印出占据数。这个占据数是否就是原子轨道的占据呢？
  （2）我看占据数有时候会超过1.0占据，为1.02这样的，其分析方法是mulliken。按理说不应该超过1.0的占据。这个是为什么呢？
  我的量子化学基础比较薄弱，问的不太专业，还请各位专家见谅！
  谢谢大家！

作者
Author: sobereva 时间: 2023-9-16 16:15
1 我不知道你说的打印是哪里、具体怎么打印的
Multiwfn基于CP2K的molden文件做Mulliken布居分析倒是可以给你最终状态的各个原子轨道的占据数信息

2 Mulliken分析自身的理论层面的问题，得到的原子轨道对分子轨道的贡献可能超过100%，下文有充分的讨论
《分子轨道成分的计算》（http://sioc-journal.cn/Jwk_hxxb/CN/abstract/abstract340458.shtml）

作者
Author: lujun 时间: 2023-9-17 14:02
本帖最后由 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原子轨道上的占据。这是不是有问题呢？

作者
Author: sobereva 时间: 2023-9-18 10:30

lujun 发表于 2023-9-17 14:02
谢谢卢天老师！
我看cp2k讲义的"DFT+U的计算"章节的45页上，列出了打印的UO2的f shell上的电子占据信息 ...

根据基函数与原子轨道的对应描述关系，可以从前者的占据数判断后者的占据数

作者
Author: lujun 时间: 2023-11-4 22:11

sobereva 发表于 2023-9-18 10:30
根据基函数与原子轨道的对应描述关系，可以从前者的占据数判断后者的占据数

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

作者
Author: sobereva 时间: 2023-11-5 11:56

lujun 发表于 2023-11-4 22:11
谢谢卢天老师！MOLOPT的基组的基函数和原子轨道是怎么对应的哇？MOLOPT的基函数是不是就是原子轨道？

不是

对应关系按下文判断
利用布居分析判断基函数与原子轨道的对应关系
http://sobereva.com/418（http://bbs.keinsci.com/thread-10094-1-1.html）

作者
Author: lujun 时间: 2023-11-5 12:55

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轨道
不知道我的理解是否正确，还请卢天老师批评指正！

作者
Author: sobereva 时间: 2023-11-5 13:02

lujun 发表于 2023-11-5 12:55
谢谢卢天老师的指导！我按照这个帖子做了一下：
我测试的是Sc原子，采用DZVP的MOLOPT基组，使用Lowdin方 ...

正确

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