计算化学公社

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

[新手求助] 求助:电子轨道占据数

[复制链接 Copy URL]

112

帖子

0

威望

471

eV
积分
583

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
请教一下各位老师:(1)何为电子轨道占据数,其数值都是整数吗?(2)电子轨道占据数与自然轨道占据数是同一个概念吗?(3)已知后HF计算下用pop=NO可以输出自然轨道占据数,请问在DFT计算下用什么关键词输出自然轨道占据数呢?
路漫漫其修远兮,吾将上下而求索

4104

帖子

4

威望

8869

eV
积分
13053

Level 6 (一方通行)

MOKIT开发者

2#
发表于 Post on 2025-5-14 16:44:54 | 只看该作者 Only view this author
本帖最后由 zjxitcc 于 2025-5-14 16:47 编辑

没有“电子轨道占据数”,只有“自然轨道占据数”,简称 轨道占据数。虽然你说“电子轨道占据数”也能听明白是自然轨道占据数,但不太常用。

“DFT计算”比较笼统,没说明是用限制性闭壳层、限制性开壳层,还是非限制性开壳层方法做计算,因此需要先分为三类:
(1)RHF/RKS方法有无数套等价的自然轨道,轨道占据数一定是整数2,2,...,0,0,其中含nocc个2,nvir个0。nocc是双占据轨道数目,nvir是空轨道数目。既然占据数是确定的值,轨道系数不是唯一的,故而无需输出、输出了也没意义。
(2)ROHF/ROKS方法有无数套等价的自然轨道,轨道占据数一定是整数2,2,...,1,1,...,0,0,其中含nocc个2,nopen个1,nvir个0。nocc是双占据轨道数目,nopen是单占据轨道数目,nvir是空轨道数目。同样无需输出、输出了也没意义。
(3)UHF/UKS空间自然轨道,轨道占据数为小数,在0.0~2.0之间。如何输出/打印,请见下文中的“UHF自然轨道”

看我写的这两篇即可,一篇原理,一篇实操
从密度矩阵产生自然轨道-理论篇》,《从密度矩阵产生自然轨道_实战篇(上)
如果都看不懂,则需要先学习基本量化知识,比如阅读Szabo的《Modern Quantum Chemistry》一书前3章。

自动做多参考态计算的程序MOKIT

112

帖子

0

威望

471

eV
积分
583

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2025-5-14 16:54:15 | 只看该作者 Only view this author
zjxitcc 发表于 2025-5-14 16:44
没有“电子轨道占据数”,只有“自然轨道占据数”,简称 轨道占据数。虽然你说“电子轨道占据数”也能听明 ...

嗯嗯,好的,谢谢邹老师,回头我确实得好好学习一下《Modern Quantum Chemistry》这本书。
路漫漫其修远兮,吾将上下而求索

6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

4#
发表于 Post on 2025-5-14 18:56:07 | 只看该作者 Only view this author
DFT得说清楚什么泛函

非双杂化泛函下,DFT的轨道占据数非整数就这些情况:
(1)自旋极化计算(UKS),然后用比如以下方法产生的UNO轨道
在Multiwfn中基于fch产生自然轨道的方法与激发态波函数、自旋自然轨道分析实例
http://sobereva.com/403http://bbs.keinsci.com/thread-9001-1-1.html
(2)有限温度DFT计算,如CP2K、ORCA都可以做,直接会输出非整数占据数
(3)有的程序如ORCA支持占据数在简并轨道上平摊(%scf FracOcc true end),比如这样以UKS方式计算基态碳原子,可以得到每个beta p轨道占据数是1/3
以上(1)属于自然轨道。(2)和(3)虽然也是非整数占据数但不属于自然轨道。普通DFT计算产生的占据数为整数的MO完全扯不上自然轨道

如果是双杂化泛函,可以要求产生相应级别的密度矩阵并得到自然轨道,占据数也为小数,类似于MP2的情况。


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

112

帖子

0

威望

471

eV
积分
583

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2025-5-14 20:00:36 | 只看该作者 Only view this author
sobereva 发表于 2025-5-14 18:56
DFT得说清楚什么泛函

非双杂化泛函下,DFT的轨道占据数非整数就这些情况:

好的,谢谢老师!
路漫漫其修远兮,吾将上下而求索

本版积分规则 Credits rule

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

GMT+8, 2025-8-14 04:41 , Processed in 0.158633 second(s), 20 queries , Gzip On.

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