| 7 | 1/1 | 返回列表 |
| 查看: 2425 | 回復(fù): 6 | ||||
gaven0728新蟲 (小有名氣)
|
[求助]
【求助】castep計(jì)算彈性常數(shù)和熱力學(xué)性質(zhì)問(wèn)題 已有1人參與
|
我參考網(wǎng)上教程用castep計(jì)算了BN的彈性常數(shù)和Ge的熱力學(xué)性質(zhì),發(fā)現(xiàn)原胞和晶胞(單胞)的結(jié)果有差距,不明白怎么回事。(1)計(jì)算BN的彈性常數(shù)時(shí)用了兩種方法,第一種是先把數(shù)據(jù)庫(kù)導(dǎo)入的BN晶胞轉(zhuǎn)為原胞,然后結(jié)構(gòu)優(yōu)化,利用結(jié)構(gòu)優(yōu)化好之后的原胞直接進(jìn)行彈性常數(shù)計(jì)算。第二種方法是利用原胞結(jié)構(gòu)優(yōu)化,優(yōu)化好之后轉(zhuǎn)為單胞,然后進(jìn)行彈性常數(shù)計(jì)算。計(jì)算結(jié)果如下,請(qǐng)問(wèn)各位原胞和單胞的計(jì)算結(jié)果差距大的原因?以及到底應(yīng)該用原胞算還是用單胞算晶格常數(shù)?? 單胞: id i j Cij (GPa) 1 1 1 713.55526 +/- 42.164 4 4 4 407.07452 +/- 5.991 7 1 2 128.73305 +/- 19.964 原胞: id i j Cij (GPa) 1 1 1 690.59274 +/- 26.342 4 4 4 381.91248 +/- 12.475 7 1 2 108.47395 +/- 10.319 (2)同樣用Ge計(jì)算了熱力學(xué)性質(zhì),參數(shù)一樣的情況下,原胞和 單胞的德拜溫度一樣,但熵、含、比熱容等單胞是原胞的4倍(因?yàn)镚e的單胞原子數(shù)為8,原胞的原子數(shù)為2,相差四倍)。請(qǐng)問(wèn)我到底是用原胞還是單胞的結(jié)果?? 總之,到底用原胞還是單胞來(lái)算晶格常數(shù)和熱力學(xué)性質(zhì) ????? |
新蟲 (小有名氣)
|
非常感謝你的回答。是不是算彈性常數(shù)時(shí)原胞和單胞一樣,算熱力學(xué)性質(zhì)時(shí)原胞用,然后結(jié)果分析時(shí)統(tǒng)一成單胞的? 發(fā)自小木蟲IOS客戶端 |
專家顧問(wèn) (著名寫手)
Ab Initio Amateur
![]() |
專家經(jīng)驗(yàn): +224 |
|
強(qiáng)度量原則上單胞、原胞、超胞結(jié)果一樣,廣延量則正比原子數(shù),相當(dāng)于使用對(duì)稱性多或少的問(wèn)題,但偶爾會(huì)因?yàn)閿?shù)值誤差造成些許差異(收斂性測(cè)試沒(méi)做好和收斂標(biāo)準(zhǔn)不夠高等因素)。彈性常數(shù)應(yīng)該是強(qiáng)度量,所以不應(yīng)有差異,造成差異的原因比較可能的原因我覺(jué)得還是收斂性測(cè)試沒(méi)做好,尤其是k點(diǎn)使用會(huì)受到晶格大小影響。另外,算熱力學(xué)量最好檢查下虛頻,并且聲子計(jì)算應(yīng)當(dāng)與構(gòu)型優(yōu)化使用相同的計(jì)算參數(shù)。 |

新蟲 (小有名氣)
|
卡卡老師,感謝您的回答。我剛接觸第一性原理計(jì)算,還有兩個(gè)問(wèn)題想請(qǐng)教下。 第一個(gè)是利用castep算phonons時(shí),有些晶體結(jié)構(gòu)的某個(gè)位置可能被多個(gè)原子占據(jù),這個(gè)時(shí)候計(jì)算phonons時(shí)會(huì)提示disordered crystal,對(duì)于這種情況如何計(jì)算phonons? 另外一個(gè)問(wèn)題是,在計(jì)算phonos提示以下錯(cuò)誤,請(qǐng)問(wèn)如何改參數(shù)? Warning in parameters_validate: current value of ELEC_EIGENVALUE_TOL = 0.100000E-08eV is too large to achieve desired level of convergence of response properties. This may cause convergence failures and/or inaccuracy of results of PHONON calculations - recommend you use a smaller value, e.g. ELEC_EIGENVALUE_TOL = 0.416667E-10eV |
|
| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600 材料與化工 295 求調(diào)劑 +9 | dream…… 2026-03-10 | 10/500 |
|
|---|---|---|---|---|
|
[基金申請(qǐng)]
NSFC申報(bào)書里申請(qǐng)人簡(jiǎn)歷中代表性論著還需要在申報(bào)書最后的附件里面再上傳一遍嗎
20+4
|
NSFC2026我來(lái)了 2026-03-10 | 10/500 |
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-11 | 4/200 |
|
|
[考研] 一志愿中科大材料與化工。353分能調(diào)劑到哪? +5 | MayUxw1 2026-03-08 | 7/350 |
|
|
[考研] 0857環(huán)境調(diào)劑 +5 | 熠熠_11 2026-03-10 | 5/250 |
|
|
[考研] 279求調(diào)劑 +3 | 莫xiao 2026-03-10 | 4/200 |
|
|
[考研] 298求調(diào)劑 +3 | Vv呀! 2026-03-10 | 3/150 |
|
|
[考研] 085602化工求調(diào)劑 +7 | 董boxing 2026-03-10 | 7/350 |
|
|
[考研] 一志愿天津大學(xué),英一數(shù)二305分求調(diào)劑,四六級(jí)已過(guò) +7 | 小小番的茄 2026-03-09 | 7/350 |
|
|
[考研] 08工科 +5 | li李樂(lè)成 2026-03-06 | 5/250 |
|
|
[考研] 物理學(xué)求調(diào)劑 +4 | sx宋 2026-03-05 | 4/200 |
|
|
[考研] 求0856調(diào)劑 +5 | 小力氣珂珂 2026-03-08 | 5/250 |
|
|
[考研] 334求調(diào)劑 +8 | Trying] 2026-03-06 | 8/400 |
|
|
[考研] 誠(chéng)求調(diào)劑,323分 有專利/科研/實(shí)習(xí)經(jīng)歷 +6 | 橙子cyx 2026-03-06 | 8/400 |
|
|
[考研] 304求調(diào)劑 +4 | 52hz~~ 2026-03-05 | 5/250 |
|
|
[考研] 276求調(diào)劑 本科天津工業(yè)大學(xué) 一志愿北京化工大學(xué)0855 +4 | 路lyh123 2026-03-06 | 6/300 |
|
|
[考研] 求調(diào)劑推薦 +4 | 微辣不吃 2026-03-06 | 4/200 |
|
|
[考研] 哈爾濱理工大學(xué)2026年研究生調(diào)劑,材料科學(xué)與化學(xué)工程學(xué)院研究生調(diào)劑 +3 | xinliu866 2026-03-06 | 3/150 |
|
|
[考研] 316求調(diào)劑 +3 | 林小星發(fā)大財(cái) 2026-03-05 | 5/250 |
|
|
[考博] 2026申博自薦 六級(jí)440電催化方向 +4 | 櫻落成影花成雙 2026-03-05 | 4/200 |
|