24小时热门版块排行榜    

查看: 1654  |  回复: 5
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

jiangning198511

木虫 (正式写手)

[交流] 【求助】密度泛函理论中密度函数的定义问题

在看关于密度泛函理论书籍的时候,不同的书对密度函数的定义不太一样,其中主要有以下两种,希望高手解答以下两者的联系,最后解释一下

[ Last edited by cadick on 2009-12-12 at 00:15 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiangning198511

木虫 (正式写手)


fegg7502(金币+1,VIP+0):鼓励交流2 12-13 17:10
谢谢,请问能不能在详细的推导一下
为什么
<\Psi|\sum_{i}\delta(r-r_i)|\Psi>
=\int{dr}\int{dr'}<\Psi|r>

=(1)
4楼2009-09-14 13:48:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

zxyan

至尊木虫 (文坛精英)

★ ★ ★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
mingdong(金币+1,VIP+0):感谢交流! 9-12 00:27
fegg7502(金币+1,VIP+0):thank you very much! 9-12 00:30
jiangning198511(金币+5,VIP+0):解答 9-14 13:48
一种是积分形式,一种是矩阵形式,不过两种都是对电子波函数的求解。
2楼2009-09-11 18:10:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yongleli

木虫 (正式写手)

★ ★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
mingdong(金币+1,VIP+0):感谢交流! 9-12 18:44
jiangning198511(金币+5,VIP+0):评分 9-14 13:48
两个公式是等价的,有什么不一样?

  <\Psi|\sum_{i}\delta(r-r_i)|\Psi>
=\int{dr}\int{dr'}<\Psi|r>

=(1)

其中,
<\Psi|r>=\Psi(r)^{\star}
=\Psi(r')

完全是根据定义推导。
3楼2009-09-12 14:45:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

c111999

铜虫 (小有名气)

★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
fegg7502(金币+5,VIP+0):thank you very much! 12-13 17:10
引用回帖:
Originally posted by jiangning198511 at 2009-9-10 21:44:
在看关于密度泛函理论书籍的时候,不同的书对密度函数的定义不太一样,其中主要有以下两种,希望高手解答以下两者的联系,最后解释一下
[ Last edited by cadick on 2009-12-12 at 00:15 ]

自认为不是高手,请谨慎使用下列信息。

(1): 单电子 自旋轨道 函数定义式 积分形式;
(1)=N* ρ(ri), ri=(x,y,z,s)i
对(1)继续积分则得到整个空间所有电子数,=N.
也可以认为,如果有波函数ψ(r1,...rN)能够描述整个系统的话,它能够被分解为多个单电子函数。即,对于ψψ*进行除了i 电子的积分,将得到 i 电子的密度 ρ(ri),。这一招是迫不得已,因为谁也不知道 ψ,所以只是定义式。
(2) 意义同(1),矩阵形式
< > = 对于ψ所有空间积分,
∫_r1∫_r2…∫∣ψ(r1,r2...rN)∣^2 dr1dr2....rN

与(1)表达不同; (1)只积分到r2, r1部分还未进行。但是由于使用了迪拉克单点函数,则使(2)展开为:(其中第一项,i=1)

∫_r1*[∫_r2…∫∣ψ(r1,r2...rN)∣^2 dr2....rN]*delta(r-r1)*dr1
= ρ(r1)
= ∫_r2…∫_rN∣ψ(r1,r2...rN)∣^2 dr2....rN]* = (1)式/N

迪拉克单点函数(∫f(x)delta(x-x0)dx=f(x0))的使用,实际上将ψ 又化为了N个单电子部分。所有空间r将会被影射/分配到r1~rN的空间中(可以认为是一坨电子云被投影到(x,y,z,s)1,。。。,(x,y,z,s)_N的空间去了)。如果把求和拿出<>来,则电子密度展开为:ρ(r)= ρ(r1)+ρ(r2)...+ρ(rN),(二阶项都被忽略了 ρ(ri,rj),原子间重叠电荷没有涉及 )
矩阵式的下一步是用ψ的展开式(如LCAO),将电子密度变为矩阵元(素),分类(奇偶性)逐个表达(存于内存中,不必重复计算……),出于实际操作的需要。

请结合:【讨论】密度泛函理论中“密度”究竟指什么
http://muchong.com/bbs/viewthread.php?tid=1642769&fpage=1

[ Last edited by c111999 on 2009-12-14 at 05:24 ]
5楼2009-12-12 12:45:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +11 医学老男孩 2026-08-13 29/1450 2026-08-14 00:04 by 大快人心啦
[基金申请] 关于Filecode分析方法 +9 majunge000 2026-08-10 12/600 2026-08-13 23:20 by iwuli
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +3 启萌科技 2026-08-12 11/550 2026-08-13 22:57 by 启萌科技
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +22 Lanmanbaby 2026-08-09 36/1800 2026-08-12 22:50 by sdfapple719
[基金申请] 好奇怪的filecode +5 布布和一二 2026-08-08 6/300 2026-08-12 16:11 by 云上清扬
[论文投稿] 职称评审,求友友推荐见刊最快的期刊 +4 工厂打螺丝 2026-08-08 4/200 2026-08-12 09:18 by hansi2025
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +7 布布和一二 2026-08-10 10/500 2026-08-12 09:12 by 布布和一二
[基金申请] 是这周出结果还是下周出结果? +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[硕博家园] 读博的好处 +3 lnee 2026-08-11 3/150 2026-08-11 18:10 by 希望我好好的
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
信息提示
请填处理意见