24小时热门版块排行榜    

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

jiangning198511

木虫 (正式写手)

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

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

[ Last edited by cadick on 2009-12-12 at 00:15 ]
回复此楼
已阅   回复此楼   关注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的回帖
查看全部 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的回帖

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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 工科材料085601 279求调剂 +6 困于星晨 2026-03-17 8/400 2026-03-20 11:24 by kkcoco25
[考研] 317求调剂 +4 申子申申 2026-03-19 9/450 2026-03-20 11:08 by 申子申申
[考研] 一志愿吉林大学材料学硕321求调剂 +11 Ymlll 2026-03-18 14/700 2026-03-20 10:46 by EBSD
[考研] 08工学调剂 +4 用户573181 2026-03-20 4/200 2026-03-20 10:32 by 朗月清风1
[考研] 271材料工程求调剂 +7 .6lL 2026-03-18 7/350 2026-03-20 09:10 by xingguangj
[考研] 329求调剂 +6 想上学吖吖 2026-03-19 6/300 2026-03-20 09:00 by 每天只摆一小会
[考研] 一志愿苏州大学材料求调剂,总分315(英一) +3 sbdksD 2026-03-19 3/150 2026-03-19 23:21 by fmesaito
[考研] 307求调剂 +9 冷笙123 2026-03-17 9/450 2026-03-19 22:44 by 学员8dgXkO
[考研] 288求调剂 +15 于海海海海 2026-03-19 15/750 2026-03-19 22:41 by 学员8dgXkO
[考研] 294求调剂材料与化工专硕 +14 陌の森林 2026-03-18 14/700 2026-03-19 22:38 by 学员8dgXkO
[考博] 申博26年 +3 八6八68 2026-03-19 3/150 2026-03-19 19:43 by nxgogo
[考研] 085600材料与化工调剂 324分 +10 llllkkkhh 2026-03-18 12/600 2026-03-19 14:33 by llllkkkhh
[考研] 一志愿华中科技大学,080502,354分求调剂 +4 守候夕阳CF 2026-03-18 4/200 2026-03-18 22:16 by li123456789.
[考研] 311求调剂 +4 冬十三 2026-03-18 4/200 2026-03-18 21:47 by 尽舜尧1
[考研] 070300化学319求调剂 +6 锦鲤0909 2026-03-17 6/300 2026-03-18 13:22 by Iveryant
[考研] 312求调剂 +8 陌宸希 2026-03-16 9/450 2026-03-18 12:39 by Linda Hu
[基金申请] 被我言中:新模板不强调格式了,假专家开始管格式了 +4 beefly 2026-03-14 4/200 2026-03-17 22:04 by 黄鸟于飞Chao
[考研] 268求调剂 +8 一定有学上- 2026-03-14 9/450 2026-03-17 17:47 by laoshidan
[基金申请] 今年的国基金是打分制吗? 50+3 zhanghaozhu 2026-03-14 3/150 2026-03-16 17:07 by 北京莱茵润色
[考研] 304求调剂 +4 ahbd 2026-03-14 4/200 2026-03-16 16:48 by 我的船我的海
信息提示
请填处理意见