|
|
[求助]
espresso中CPMD溫度漂移的非常嚴(yán)重如何解決,求高人指點(diǎn)。
我目前主要做的是BN材料的脫氫機(jī)理問題,相通過espresso中的cp模塊解決其釋氫機(jī)理問題,事實(shí)上是想模仿這篇文章“First-Principles Study on the Mechanisms for H2 Formation in Ammonia Borane at Ambient and High Pressure”(J. Phys. Chem. C 2012, 116, 2146–2152)解決自己體系中的問題。
按照cp.x的tutorial做了如下的步驟 (在本體系中,固定住晶體中大部分非關(guān)鍵的原子,只想觀察局部的化學(xué)反應(yīng)):1)我們的晶體結(jié)構(gòu)是在castep中建模和優(yōu)化的,之后導(dǎo)入espresso中; 2)通過damp方法優(yōu)化了波函數(shù), 確定各原子力在10^-4左右; 3)對(duì)原子做randomization,范圍是0.01au;4)再次測(cè)試了波函數(shù)(damp);5)做MD的初始化,v=0, 沒有加溫度等任何限制; 6)在5的基礎(chǔ)上,去掉速度項(xiàng);7)在6的基礎(chǔ)上,加上溫度項(xiàng),用nose加溫項(xiàng),考慮到實(shí)際反應(yīng)是373K左右,我們?cè)O(shè)定的溫度為373K,fnosep設(shè)為10,溫度緩慢升至373K左右就開始飄,最高飄到1000K。
附件中為各項(xiàng)能量結(jié)果,random為randomiztion的測(cè)試,AMgB-wav4是random后的波函數(shù)測(cè)試,md1,md2,md3是上述5,6,7的結(jié)果;附上最后一步md的輸入文件。謝謝,求指點(diǎn)。 |
» 本帖附件資源列表
-
歡迎監(jiān)督和反饋:小木蟲僅提供交流平臺(tái),不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布,如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問題,其責(zé)任在于用戶本人,如對(duì)版權(quán)有異議,請(qǐng)聯(lián)系郵箱:xiaomuchong@tal.com
- 附件 1 : AMgB-md1.evp
2014-06-05 22:06:48, 113.28 K
- 附件 2 : AMgB-md2.evp
2014-06-05 22:06:51, 113.28 K
- 附件 3 : AMgB-md3.evp
2014-06-05 22:06:53, 113.28 K
- 附件 4 : AMgB-md3.in
2014-06-05 22:06:55, 11.93 K
- 附件 5 : AMgB-wav4.evp
2014-06-05 22:06:57, 17.9 K
» 猜你喜歡
|