计算化学公社

 找回密码 Forget password
 注册 Register

Ne元素计算ELF结果和角度无关的问题

查看数: 623 | 评论数: 7 | 收藏 Add to favorites 0
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2023-4-14 11:51

正文摘要:

各位老师好!我在复现ELF的那篇原始文献(A simple measure of electron localization in atomic and molecular systems)中的图,用和原始文献相同的Roothaan-HF计算结果输入计算ELF,我用python自己写程序画了Ne波 ...

回复 Reply

Rurutiaxy 发表于 Post on 2023-5-10 17:11:38

感谢老师的回复!按照您的方法确实可以证明Ne的动能密度不随角度变化
978142355 发表于 Post on 2023-5-1 00:13:40
具体如下:

202305010013141769..png (35.24 KB, 下载次数 Times of downloads: 21)

202305010013141769..png
wzkchem5 发表于 Post on 2023-4-29 16:06:29
Rurutiaxy 发表于 2023-4-27 08:05
接上一楼,因为s轨道本身就是角度无关的,那个球谐函数是个常数,所以计算ELF的话,只要看p轨道的动能密度 ...

你是把球谐函数从极坐标转到直角坐标,再求的nabla吧?
此外,你用的是球谐函数的实数形式还是复数形式?如果是复数形式,有没有注意到psi的梯度的模平方不是(x^2+y^2+z^2)的形式,而是(|x|^2 + |y|^2 + |z|^2)的形式?
Rurutiaxy 发表于 Post on 2023-4-27 15:05:24
接上一楼,因为s轨道本身就是角度无关的,那个球谐函数是个常数,所以计算ELF的话,只要看p轨道的动能密度和单粒子密度就好,单粒子密度被证明是角度无关的,那么ELF的第二项,就是对单粒子密度求梯度模方的那一项肯定是角度无关的,角度相关性就只剩下动能密度一项影响了,我也查了不少书,问了一些老师同学,他们也都觉得如果是我这么写的这种形式,动能密度应该是和角度有关的。所以我就是在这一点上和论文对不上。
granvia 发表于 Post on 2023-4-14 21:04:45
本帖最后由 granvia 于 2023-4-14 21:09 编辑

从物理意义上说,全充满壳层,各电子轨道角动量完全抵消,总轨道角动量为零,从而体系是球对称的(空间旋转不变性)。从数学上也可以严格证明,也就是Uns”old定理
sobereva 发表于 Post on 2023-4-14 15:01:50
只考虑一个或两个p时和角度有关,同时考虑三个p和角度无关。如同x^2+y^2+z^2=1对应的是个半径为1的完全对称的球面,角度无关

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

GMT+8, 2024-11-27 17:55 , Processed in 0.259976 second(s), 32 queries , Gzip On.

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