| 9 | 1/1 | 返回列表 |
| 查看: 8792 | 回復(fù): 8 | ||||
spphins510新蟲 (小有名氣)
|
[交流]
分享與交流CASTEP 計(jì)算中截?cái)嗄芎蚄值的選擇經(jīng)驗(yàn) 已有6人參與
|
|
一直在小木蟲索取,這次也想分享一下~雖然很基礎(chǔ)的內(nèi)容,但感覺對(duì)初學(xué)者,尤其是基礎(chǔ)比較差的(我自己)還是會(huì)有幫助的。這些內(nèi)容來源于軟件說明書、小木蟲的討論等資源。而且如果有不對(duì)的地方,還請(qǐng)路過的大家能夠指點(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)。另一方面,面心立方的鋁可能需要上萬個(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 或者說 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í)積累 |
新蟲 (小有名氣)
木蟲 (正式寫手)

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

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

木蟲 (文壇精英)

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

鐵蟲 (初入文壇)

| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0857環(huán)境調(diào)劑 +3 | 熠熠_11 2026-03-10 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 一志愿東北大學(xué)080500,英一數(shù)二343分求調(diào)劑 +8 | Seele_v5 2026-03-04 | 8/400 |
|
|
[考研] 一志愿:武漢理工,材料工程,英二數(shù)二 總分314 +3 | 2202020125 2026-03-10 | 4/200 |
|
|
[基金申請(qǐng)]
PDF在線壓縮軟件
10+3
|
haxia 2026-03-08 | 4/200 |
|
|
[考研] 材料與化工 一志愿山大 321分 求調(diào)劑 +6 | 每天散步 2026-03-09 | 7/350 |
|
|
[考研] 0703化學(xué)調(diào)劑 +7 | 小奶蒼 2026-03-07 | 7/350 |
|
|
[考研] 294 英二數(shù)二物化 求調(diào)劑 +6 | 米飯團(tuán)不好吃 2026-03-09 | 6/300 |
|
|
[考研] 材料工程085601調(diào)劑求老師收留 +5 | 強(qiáng)木木木 2026-03-07 | 5/250 |
|
|
[考研] 296求調(diào)劑 +4 | Xinyu Wu311 2026-03-09 | 4/200 |
|
|
[考研] 中科大材料299求調(diào)劑 +10 | DAIjiayo 2026-03-05 | 16/800 |
|
|
[考研] 求0856調(diào)劑 +5 | 小力氣珂珂 2026-03-08 | 5/250 |
|
|
[教師之家] 交大前校長(zhǎng)王樹國:現(xiàn)在最先進(jìn)的科技并不在大學(xué)實(shí)驗(yàn)室,而是在企業(yè)研究院 +4 | zju2000 2026-03-08 | 6/300 |
|
|
[考研] 322分 085600求調(diào)劑,有互聯(lián)網(wǎng)+國金及主持省級(jí)大創(chuàng)經(jīng)歷 +4 | 熊境喆 2026-03-05 | 4/200 |
|
|
[考研] 誠求調(diào)劑,323分 有專利/科研/實(shí)習(xí)經(jīng)歷 +6 | 橙子cyx 2026-03-06 | 8/400 |
|
|
[考研] 276求調(diào)劑 本科天津工業(yè)大學(xué) 一志愿北京化工大學(xué)0855 +4 | 路lyh123 2026-03-06 | 6/300 |
|
|
[考研] 復(fù)試調(diào)劑 +7 | 呼呼?~+123456 2026-03-05 | 10/500 |
|
|
[考研] 085602高分子方向求調(diào)劑 +7 | tlgudy 2026-03-04 | 7/350 |
|
|
[考研] 0856材料專碩274能調(diào)劑去哪里? +3 | 22735 2026-03-04 | 4/200 |
|
|
[考研] 293求調(diào)劑 +3 | 是樂渝哇 2026-03-04 | 3/150 |
|
|
[考研] 085601 材料305分求助 +4 | 泡泡郵件 2026-03-03 | 6/300 |
|