| 查看: 4140 | 回復(fù): 32 | ||||
物理毀一生新蟲 (初入文壇)
|
[交流]
關(guān)于密度泛函計(jì)算的疑惑(初學(xué)) 已有5人參與
|
|||
|
一、原來求解薛定諤方程的波函數(shù)(本證函數(shù))和能量(本征值)有3N個(gè)變量,KS證明通過密度(僅3個(gè)變量)也能求出總能量(不過是基態(tài)的),所以現(xiàn)在要對(duì)含有電子密度的能量泛函理進(jìn)行求解。 二、給出K-S方程,能量可有密度表示,E[ro] = Ts[ro] + Exc[ro] + EH[ro] + Eext[ro],其中交換關(guān)聯(lián)項(xiàng)可采用LDA,GGA等多種近似,而其他幾項(xiàng)都有清晰的表達(dá)式。 三、利用泛函分析,利用SCF求解使能量最低的密度函數(shù)。(其中數(shù)學(xué)方法暫時(shí)不討論,因?yàn)槲铱床欢?br /> 四、利用得到的基態(tài)能量和電子密度計(jì)算體系的光學(xué)、電學(xué)等性質(zhì)。 我的理解是這樣的,那么,問題來了,挖掘機(jī)技術(shù)哪家強(qiáng)...... Q1:在計(jì)算過程中,我只能求得基態(tài)能量和電子密度,那么體系的波函數(shù)我能夠求得嗎?如果不知道波函數(shù),怎么可能了解其它的物理性質(zhì)呢?只有波函數(shù)才能反映體系的狀態(tài)啊! Q2:K-S方程中動(dòng)能項(xiàng)Ts[ro]的表達(dá)式需要各電子的波函數(shù),而不僅僅是電子密度,求解過程中能夠得到嗎? Q3:LDA的物理圖像是什么,能用LDA近似需要滿足什么樣的條件? Q4:在計(jì)算過程中(使用軟件時(shí)),自洽求解時(shí)方法是不是唯一的?是不是我們只用選擇交換關(guān)聯(lián)參數(shù)和邊界限制? (任意一個(gè)問題都可,求大神指導(dǎo)) |
光催化 |

專家顧問 (著名寫手)
Ab Initio Amateur
![]() |
專家經(jīng)驗(yàn): +224 |
|
Q1:KS的框架雖然是“密度”泛函,但是求解過程還是借助了“(贗)波函數(shù)”來進(jìn)行計(jì)算。但是這個(gè)波函數(shù)在Exc非精確的情況下是非物理的,但目前沒有一個(gè)Exc是真正意義上的精確。 Q2:如上面所說,自洽過程借助的是波函數(shù),計(jì)算密度則通過<ψ|ψ>進(jìn)行,都是基于波函數(shù)進(jìn)行的。 Q3:LDA的導(dǎo)出基于均勻電子氣,一般比較適合電子密度變化較為緩慢的體系。 Q4:這個(gè)沒明白你指的方法是?自洽求解可以采用不同的展開基,比如LCAO或者是平面波等,求解的方式從數(shù)學(xué)的角度稍有差異,廣義來說沒什么區(qū)別。 |

鐵桿木蟲 (著名寫手)

鐵桿木蟲 (著名寫手)
|
LZ對(duì)密度泛函和K-S方程的關(guān)系理解錯(cuò)了。 密度泛函是H-K的定理,K-S方程是“在密度泛函理論框架下求解問題”的方法。 如樓上已經(jīng)說了,KS方法借用了一個(gè)非相互作用系統(tǒng)的波函數(shù),用來表達(dá)動(dòng)能。這個(gè)波函數(shù)和真實(shí)體系的波函數(shù)沒啥必然聯(lián)系。但是,這個(gè)波函數(shù)仍然可以用來計(jì)算各種物理性質(zhì)。 Q4:自洽是個(gè)數(shù)學(xué)問題,有多種可行的數(shù)學(xué)方法,一般不需要操心。程序一般會(huì)給選項(xiàng)來選方法。 |
新蟲 (初入文壇)

新蟲 (初入文壇)

專家顧問 (著名寫手)
Ab Initio Amateur
![]() |
專家經(jīng)驗(yàn): +224 |

鐵桿木蟲 (著名寫手)
鐵桿木蟲 (著名寫手)

鐵桿木蟲 (著名寫手)
鐵桿木蟲 (著名寫手)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 326求調(diào)劑 +5 | 諾貝爾化學(xué)獎(jiǎng)覬?/a> 2026-03-15 | 8/400 |
|
|---|---|---|---|---|
|
[基金申請(qǐng)] 山東省面上項(xiàng)目限額評(píng)審 +3 | 石瑞0426 2026-03-19 | 3/150 |
|
|
[考研] 材料求調(diào)劑 +4 | @taotao 2026-03-21 | 4/200 |
|
|
[考研] 材料學(xué)碩333求調(diào)劑 +3 | 北道巷 2026-03-18 | 3/150 |
|
|
[基金申請(qǐng)]
學(xué)校已經(jīng)提交到NSFC,還能修改嗎?
40+4
|
babangida 2026-03-19 | 9/450 |
|
|
[考研] 316求調(diào)劑 +6 | 梁茜雯 2026-03-19 | 6/300 |
|
|
[考研] 化學(xué)求調(diào)劑 +4 | 臨澤境llllll 2026-03-17 | 5/250 |
|
|
[考研] 南京大學(xué)化學(xué)376求調(diào)劑 +3 | hisfailed 2026-03-19 | 6/300 |
|
|
[考研] 廣西大學(xué)家禽遺傳育種課題組2026年碩士招生(接收計(jì)算機(jī)專業(yè)調(diào)劑) +3 | 123阿標(biāo) 2026-03-17 | 3/150 |
|
|
[考研] 281求調(diào)劑(0805) +14 | 煙汐憶海 2026-03-16 | 25/1250 |
|
|
[論文投稿]
申請(qǐng)回稿延期一個(gè)月,編輯同意了。但系統(tǒng)上的時(shí)間沒變,給編輯又寫郵件了,沒回復(fù)
10+3
|
wangf9518 2026-03-17 | 4/200 |
|
|
[考研] 288求調(diào)劑,一志愿華南理工大學(xué)071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|
[考研] 材料專碩306英一數(shù)二 +10 | z1z2z3879 2026-03-16 | 13/650 |
|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個(gè)二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
[考研]
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
[考研] 機(jī)械專碩325,尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
[考研] 0856求調(diào)劑 +3 | 劉夢(mèng)微 2026-03-15 | 3/150 |
|