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

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

木蟲(chóng) (文壇精英)
|
截?cái)嗄,改為custom(定值),多設(shè)幾個(gè),做出總能量隨截?cái)嗄荜P(guān)系曲線,斜率趨于0即可 發(fā)自小木蟲(chóng)Android客戶(hù)端 |

木蟲(chóng) (文壇精英)

木蟲(chóng) (文壇精英)
|
還有,研究表面,最后一個(gè)數(shù)字改成1,比如,如果晶體用k值999得到能量最小,表面就改為991 發(fā)自小木蟲(chóng)Android客戶(hù)端 |

鐵蟲(chóng) (初入文壇)

| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 歡迎211本科同學(xué),過(guò)A區(qū)國(guó)家線,A區(qū)非偏遠(yuǎn)一本,交叉學(xué)科課題組 +10 | lisimayy 2026-03-04 | 14/700 |
|
|---|---|---|---|---|
|
[考研] 085600材料與化工(高分子)290分求調(diào)劑 +8 | wengyujian 2026-03-04 | 8/400 |
|
|
[考研] 調(diào)劑材料學(xué)碩 +6 | 詞凝Y 2026-03-02 | 6/300 |
|
|
[考研] 一志愿中石油,化學(xué)工程,301求調(diào)劑 +3 | 豚饅枳枳 2026-03-04 | 3/150 |
|
|
[考研] 求調(diào)劑 +4 | 熬夜的貓頭鷹 2026-03-02 | 4/200 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +3 | 將未有期 2026-03-04 | 3/150 |
|
|
[考博] 26申博 求博導(dǎo) +3 | 愛(ài)讀書(shū)的小帥 2026-02-28 | 5/250 |
|
|
[考研] 298求調(diào)劑 +3 | 人間唯你是清歡 2026-03-03 | 4/200 |
|
|
[考研] 299求調(diào)劑 +5 | kkcoco25 2026-03-02 | 9/450 |
|
|
[考研] 化學(xué)0703求調(diào)劑 學(xué)碩 理/工科均可 總分279 +3 | 1一11 2026-03-03 | 5/250 |
|
|
[考研]
材料工程專(zhuān)碩283求調(diào)劑
5+8
|
,? 2026-03-02 | 10/500 |
|
|
[考研] 化工335求調(diào)劑 +5 | 摸摸貓貓頭 2026-03-02 | 5/250 |
|
|
[考研] 【2026 碩士調(diào)劑】課題組 招收調(diào)劑生 +3 | 考研版棒棒 2026-03-02 | 5/250 |
|
|
[考研] 11408,學(xué)碩276求調(diào)劑 +3 | 崔wj 2026-03-02 | 5/250 |
|
|
[考研] 一志愿東北大學(xué)材料專(zhuān)碩328,求調(diào)劑 +3 | shs1083 2026-03-02 | 3/150 |
|
|
[考研] 0856化工專(zhuān)碩求調(diào)劑 +15 | 董boxing 2026-03-01 | 15/750 |
|
|
[基金申請(qǐng)] 此成果不能導(dǎo)入原因:元數(shù)據(jù)必填信息不完整,可 進(jìn)行補(bǔ)充。 +4 | Kittylucky 2026-03-02 | 5/250 |
|
|
[考研] 275求調(diào)劑 +3 | L-xin? 2026-03-01 | 6/300 |
|
|
[考博] 26申博 +4 | 想申博! 2026-02-26 | 6/300 |
|
|
[考研] 307求調(diào)劑 +4 | 73372112 2026-02-28 | 6/300 |
|