| 5 | 1/1 | 返回列表 |
| 查看: 26189 | 回復(fù): 160 | ||
| 【獎(jiǎng)勵(lì)】 本帖被評(píng)價(jià)126次,作者identation增加金幣 111.05 個(gè) | ||
| 當(dāng)前只顯示滿(mǎn)足指定條件的回帖,點(diǎn)擊這里查看本話(huà)題的所有回帖 | ||
identation金蟲(chóng) (著名寫(xiě)手)
|
[資源]
【原創(chuàng)】第一原理計(jì)算的一些心得
|
|
|
把自己關(guān)于第一原理的一些理解寫(xiě)出: (1)第一性原理其實(shí)是包括基于密度泛函的從頭算和基于Hartree-Fock自洽計(jì)算的從頭算,前者以電子密度作為基本變量(霍亨伯格-科洪定理),通過(guò)求解Kohn-Sham方程,迭代自洽得到體系的基態(tài)電子密度,然后求體系的基態(tài)性質(zhì);后者則通過(guò)自洽求解Hartree-Fock方程,獲得體系的波函數(shù),求基態(tài)性質(zhì); 評(píng)述:K-S方程的計(jì)算水平達(dá)到了H-F水平,同時(shí)還考慮了電子間的交換關(guān)聯(lián)作用。 (2)關(guān)于DFT中密度泛函的Functional,其實(shí)是交換關(guān)聯(lián)泛函 包括LDA,GGA,雜化泛函等等 一般LDA為局域密度近似,在空間某點(diǎn)用均勻電子氣密度作為交換關(guān)聯(lián)泛函的唯一變量,多數(shù)為參數(shù)化的CA-PZ方案; GGA為廣義梯度近似,不僅將電子密度作為交換關(guān)聯(lián)泛函的變量,也考慮了密度的梯度為變量,包括PBE,PW,RPBE等方案,BLYP泛函也屬于GGA; 此外還有一些雜化泛函,B3LYP等。 (3)關(guān)于贗勢(shì) 在處理計(jì)算體系中原子的電子態(tài)時(shí),有兩種方法,一種是考慮所有電子,叫做全電子法,比如WIEN2K中的FLAPW方法(線(xiàn)性綴加平面波);此外還有一種方法是只考慮價(jià)電子,而把芯電子和原子核構(gòu)成離子實(shí)放在一起考慮,即贗勢(shì)法,一般贗勢(shì)法是選取一個(gè)截?cái)喟霃,截(cái)喟霃揭詢(xún)?nèi),波函數(shù)變化較平滑,和真實(shí)的不同,截?cái)喟霃揭酝鈩t和真實(shí)情況相同,而且贗勢(shì)法得到的能量本征值和全電子法應(yīng)該相同。 贗勢(shì)包括模守恒和超軟,模守恒較硬,一般需要較大的截?cái)嗄,超軟?shì)則可以用較小的截?cái)嗄芗纯伞A硗,模守恒?shì)的散射特性和全電子相同,因此一般紅外,拉曼等光譜的計(jì)算需要用模守恒勢(shì)。 贗勢(shì)的測(cè)試標(biāo)準(zhǔn)應(yīng)是贗勢(shì)與全電子法計(jì)算結(jié)果的匹配度,而不是贗勢(shì)與實(shí)驗(yàn)結(jié)果的匹配度,因?yàn)楹蛯?shí)驗(yàn)結(jié)果的匹配可能是偶然的。 (4)關(guān)于收斂測(cè)試 (a)Ecut,也就是截?cái)嗄,一般情況下,總能相對(duì)于不同Ecut做計(jì)算,當(dāng)Ecut增大時(shí)總能變化不明顯了即可;然而,在需要考慮體系應(yīng)力時(shí),還需對(duì)應(yīng)力進(jìn)行收斂測(cè)試,而且應(yīng)力相對(duì)于Ecut的收斂要比總能更為苛刻,也就是某個(gè)截?cái)嗄芟驴偰芤呀?jīng)收斂了,但應(yīng)力未必收斂。 (b)K-point,即K網(wǎng)格,一般金屬需要較大的K網(wǎng)格,采用超晶胞時(shí)可以選用相對(duì)較小的K網(wǎng)格,但實(shí)際上還是要經(jīng)過(guò)測(cè)試。 (5)關(guān)于磁性 一般何時(shí)考慮自旋呢?舉例子,例如BaTiO3中,Ba、Ti和O分別為+2,+4和-2價(jià),離子全部為各個(gè)軌道滿(mǎn)殼層的結(jié)構(gòu),就不必考慮自旋了;對(duì)于BaMnO3中,由于Mn+3價(jià)時(shí)d軌道還有電子,但未滿(mǎn),因此需考慮Mn的自旋,至于Ba和O則不必考慮。其實(shí)設(shè)定自旋就是給定一個(gè)原子磁矩的初始值,只在剛開(kāi)始計(jì)算時(shí)作為初始值使用,具體的可參照磁性物理。 (6)關(guān)于幾何優(yōu)化 包括很多種了,比如晶格常數(shù)和原子位置同時(shí)優(yōu)化,只優(yōu)化原子位置,只優(yōu)化晶格常數(shù),還有晶格常數(shù)和原子位置分開(kāi)優(yōu)化等等。 在PRL一篇文章中見(jiàn)到過(guò)只優(yōu)化原子位置,晶格常數(shù)用實(shí)驗(yàn)值的例子(PRL 100, 186402 (2008));也見(jiàn)到過(guò)晶格常數(shù)先優(yōu)化,之后固定晶格常數(shù)優(yōu)化原子位置的情況;更多的情況則是Full geometry optimization。 一般情況下,也有不優(yōu)化幾何結(jié)構(gòu)直接計(jì)算電子結(jié)構(gòu)的,但是對(duì)于缺陷形成能的計(jì)算則往往要優(yōu)化。 (7)關(guān)于軟件 軟件大致分為基于平面波的軟件,如CASTEP、PWSCF和ABINIT等等,計(jì)算量大概和體系原子數(shù)目的三次方相關(guān);還有基于原子軌道線(xiàn)性組合的軟件(LCAO),比如openmx,siesta,dmol等,計(jì)算量和體系原子數(shù)目相關(guān),一般可模擬較多原子數(shù)目的體系。 參考文獻(xiàn): Designing meaningful density functional theory calculations in materials science—a primer |
木蟲(chóng) (著名寫(xiě)手)
木蟲(chóng) (著名寫(xiě)手)
至尊木蟲(chóng) (著名寫(xiě)手)
木蟲(chóng) (正式寫(xiě)手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 297求調(diào)劑 +5 | 戲精丹丹丹 2026-03-17 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 收復(fù)試調(diào)劑生 +3 | 雨后秋荷 2026-03-18 | 3/150 |
|
|
[考研] 0703化學(xué)求調(diào)劑 總分331 +3 | ZY-05 2026-03-13 | 3/150 |
|
|
[考研] 0703化學(xué)336分求調(diào)劑 +6 | zbzihdhd 2026-03-15 | 7/350 |
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
[考研] 328求調(diào)劑,英語(yǔ)六級(jí)551,有科研經(jīng)歷 +3 | 生物工程調(diào)劑 2026-03-16 | 8/400 |
|
|
[考研] 301求調(diào)劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|
[考研] 275求調(diào)劑 +4 | 太陽(yáng)花天天開(kāi)心 2026-03-16 | 4/200 |
|
|
[考研] 285化工學(xué)碩求調(diào)劑(081700) +9 | 柴郡貓_ 2026-03-12 | 9/450 |
|
|
[考研] 085600材料與化工 求調(diào)劑 +13 | enenenhui 2026-03-13 | 14/700 |
|
|
[考研]
|
笨笨兔子 2026-03-12 | 3/150 |
|
|
[考研] 中科大材料專(zhuān)碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
[考研] 297一志愿上交085600求調(diào)劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|
|
[考研] 255求調(diào)劑 +3 | 李嘉慧, 2026-03-12 | 4/200 |
|
|
[考研] 材料080500調(diào)劑求收留 +3 | 一顆meteor 2026-03-13 | 3/150 |
|
|
[考研] 266求調(diào)劑 +4 | 學(xué)員97LZgn 2026-03-13 | 4/200 |
|
|
[基金申請(qǐng)]
有必要更換申報(bào)口嗎
20+3
|
fannyamoy 2026-03-11 | 3/150 |
|
|
[考研] 一志愿西南交大,材料專(zhuān)碩317求調(diào)劑 +5 | lx8568 2026-03-11 | 5/250 |
|
|
[考研] 315求調(diào)劑 +9 | 小羊小羊_ 2026-03-11 | 10/500 |
|
|
[考研] 328化工專(zhuān)碩求調(diào)劑 +4 | 。,。,。,。i 2026-03-12 | 4/200 |
|