| 5 | 1/1 | 返回列表 |
| 查看: 1894 | 回復(fù): 11 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
[求助]
VASP計算流程疑惑求助 已有1人參與
|
|||
|
蟲子們好,近期計算晶體材料的時候遇到一些問題,特向蟲字們求助。 目標(biāo)是求高溫高壓下礦物晶體的結(jié)構(gòu)和彈性常數(shù)。目前已用MS第一性原理計算了0K、高壓下的礦物結(jié)構(gòu)。 接下來就要用VASP5.2的MD模式計算高溫下的結(jié)構(gòu)和彈性性質(zhì)了。 具體流程我是這樣想的: 1,那0K下的結(jié)構(gòu)作為原始結(jié)構(gòu),使用MD模式的INCAR文件,計算高溫下的結(jié)構(gòu)。 2,拿上一步的輸出結(jié)構(gòu),使用求彈性模量的INCAR模板,計算對應(yīng)的彈性常數(shù)。 不知道這個流程是不是合理的? 還有就是我在castep下的“ENCUT和K點的設(shè)置”需要和MD計算相同嗎? MD計算晶體必須使用超胞嗎? 最后,本人不才,求INCAR模板。 |
送紅花一朵 |
謝謝你的答案,不過我發(fā)現(xiàn),vasp5.2現(xiàn)在支持直接計算彈性模量了,用的原理是你說的那種嗎? 還有就是,我提的那個流程能用嗎? 發(fā)自小木蟲Android客戶端 |
新蟲 (小有名氣)
| 據(jù)我了解,準(zhǔn)確計算高溫下的彈性常數(shù)很困難,一般會簡化處理,采用所謂準(zhǔn)靜態(tài)的方法,先用德拜模型求出給定結(jié)構(gòu)的體積與溫度的關(guān)系,例如,知道300K下的體積是多少,然后以此體積去寫相應(yīng)的POSCAR, 然后計算彈性常數(shù)。參看,A first-principles approach to finite temperature elastic constants,J. Phys.: Condens. Matter 22 (2010) 225404 (8pp)。其他的不知,見諒。 |
新蟲 (小有名氣)
| 你說的可能是用MD的方法計算有限溫度下的彈性常數(shù),例如,Temperature dependence of TiN elastic constants from ab initio molecular dynamics simulations,PHYSICAL REVIEW B 87, 094114 (2013);這方面抱歉不了解。VASP似乎不能直接計算有限溫度下的彈性常數(shù),而是需要跟據(jù)我前面說的,分別確定每個溫度下的POSCAR,然后采用應(yīng)變應(yīng)力法計算。至于你說的利用MD先確定給定溫度下的結(jié)構(gòu),然后設(shè)置0K時一樣的INCAR文件,方法與我說的類似,只是確定結(jié)構(gòu)的方法不一樣。這方面不了解,你可以查查有無相關(guān)的論文。 |
|
請問我要用VASP進(jìn)行彈性常數(shù)的計算,我直接設(shè)置IBRION=6和ISIF=7就可以了還是說要先進(jìn)行IBRION=2的普通優(yōu)化然后將CONTCAR作為下一個的POSCAR? 發(fā)自小木蟲Android客戶端 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|