| 5 | 1/1 | 返回列表 |
| 查看: 8791 | 回復(fù): 8 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||||
spphins510新蟲 (小有名氣)
|
[交流]
分享與交流CASTEP 計算中截斷能和K值的選擇經(jīng)驗 已有6人參與
|
|||
|
一直在小木蟲索取,這次也想分享一下~雖然很基礎(chǔ)的內(nèi)容,但感覺對初學(xué)者,尤其是基礎(chǔ)比較差的(我自己)還是會有幫助的。這些內(nèi)容來源于軟件說明書、小木蟲的討論等資源。而且如果有不對的地方,還請路過的大家能夠指點一下~謝謝啦。 截斷能 Energy cutoff:控制平面波基組的大小,平衡準確性和計算速度。在做吸附計算的時候,要分別計算氣體分子,潔凈表面和吸附后表面的,為了消除最大誤差,該值不能夠采取默認值,而應(yīng)該統(tǒng)一。 一般400eV對于計算表面科學(xué)足夠了。 K點 k-points: 確定布里淵區(qū)域采樣的準確性;k點的質(zhì)量對于金屬系統(tǒng)尤其重要,在金屬系統(tǒng)中,電子結(jié)構(gòu)的快速變化可能發(fā)生在跨越費米能級的能帶上。 經(jīng)驗:實空間越大,K空間就越小。金屬體系需要精細的網(wǎng)格,絕緣體使用很少的k 點通常就可以。小單胞需要精細格點,大單胞很可能不需要。因此:單位晶胞內(nèi)原子數(shù)很多(比如40-60 個)的絕緣體,可能僅需要一個(移動后的)k 點。另一方面,面心立方的鋁可能需要上萬個k點以獲得好的DOS。對于孤立原子或分子的超晶胞,僅需要在Gamma點計算。對于表面(層面)的超晶胞計算,僅需要(垂直于表面)z 方向上有1 個k 點。 取值依據(jù):(1)a×ka ≈ b×kb ≈ c×kc 或者說 1/a : 1/b : 1/c ≈ ka : kb : kc (2)ka~15 Å 為絕緣體 ka~20 Å 為半導(dǎo)體 ka~25 Å 為simple metals ka~30 Å 為d band metals 截斷能和k值的收斂性測試: (1)方法:首先測試Ecut,把K點設(shè)置成足夠高的值,以步長50eV逐步增加Ecut,查看總能的變化,當(dāng)變化小于0.001 eV/atom(dE/dEcut/atom)的時候就已經(jīng)足夠了。然后固定Ecut為測試好的值,逐步增加K值,以同樣的收斂標準測試就行了。 (2)對于固體表面計算,和在切面之前對晶體做收斂性分析相比,在切面后做收斂性測試更合理。 (3)但要注意改變k值能量可能不是平穩(wěn)的收斂。 |
MS學(xué)習(xí)積累 |
木蟲 (文壇精英)

新蟲 (小有名氣)
木蟲 (正式寫手)

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

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 304求調(diào)劑(085602一志愿985) +4 | 化工人999 2026-03-09 | 4/200 |
|
|---|---|---|---|---|
|
[考研]
|
簡木ChuFront 2026-03-09 | 4/200 |
|
|
[考研] 材料與化工 一志愿山大 321分 求調(diào)劑 +4 | 每天散步 2026-03-09 | 5/250 |
|
|
[考研] 0817一志愿蘇州大學(xué)280復(fù)試調(diào)劑 +4 | kk扛 2026-03-07 | 4/200 |
|
|
[考研] 數(shù)二英二309分請求調(diào)劑 +3 | dtdxzxx 2026-03-09 | 4/200 |
|
|
[考研] 334求調(diào)劑 +5 | Trying] 2026-03-07 | 5/250 |
|
|
[考研] 296求調(diào)劑 +4 | Xinyu Wu311 2026-03-09 | 4/200 |
|
|
[考研] 083000環(huán)境科學(xué)與工程調(diào)劑 +5 | 加油呀fxy 2026-03-07 | 6/300 |
|
|
[考研] 269求調(diào)劑 +3 | 朔朔話 2026-03-08 | 4/200 |
|
|
[考研] 化學(xué)工程求調(diào)劑 +12 | 化工人999 2026-03-04 | 12/600 |
|
|
[考研] 材料化工求調(diào)劑 +7 | 博斯特525 2026-03-06 | 8/400 |
|
|
[考研] 0307化學(xué)求調(diào)劑 +6 | 0ok0k 2026-03-07 | 6/300 |
|
|
[考研] 材料調(diào)劑 +10 | ounce. 2026-03-04 | 12/600 |
|
|
[考研] 求調(diào)劑推薦 +4 | 微辣不吃 2026-03-06 | 4/200 |
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +4 | 邱gl 2026-03-06 | 4/200 |
|
|
[考研] 085600材料調(diào)劑 總分330 +6 | 池池丶 2026-03-03 | 6/300 |
|
|
[考研] 085600,一志愿鄭州大學(xué),280分求調(diào)劑 +7 | Wuqi725 2026-03-05 | 7/350 |
|
|
[考研] 求材料調(diào)劑 +4 | berdmond 2026-03-05 | 4/200 |
|
|
[考研] 295求調(diào)劑 +4 | 小賽不吃香菜 2026-03-04 | 4/200 |
|
|
[考研] 085601 材料305分求助 +4 | 泡泡郵件 2026-03-03 | 6/300 |
|