24小时热门版块排行榜

返回列表

Trouton

新虫 (小有名气)

应助: 2 (幼儿园)
金币: 259.4
散金: 12
红花: 1
帖子: 54
在线: 59.5小时
虫号: 1207174
注册: 2011-02-20
性别: GG
专业: 理论和计算化学

[求助] espresso中CPMD温度漂移的非常严重如何解决，求高人指点。

我目前主要做的是BN材料的脱氢机理问题，相通过espresso中的cp模块解决其释氢机理问题，事实上是想模仿这篇文章“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做了如下的步骤（在本体系中，固定住晶体中大部分非关键的原子，只想观察局部的化学反应）：1）我们的晶体结构是在castep中建模和优化的，之后导入espresso中； 2）通过damp方法优化了波函数，确定各原子力在10^-4左右； 3）对原子做randomization，范围是0.01au；4）再次测试了波函数（damp）；5）做MD的初始化，v=0，没有加温度等任何限制； 6）在5的基础上，去掉速度项；7）在6的基础上，加上温度项，用nose加温项，考虑到实际反应是373K左右，我们设定的温度为373K，fnosep设为10，温度缓慢升至373K左右就开始飘，最高飘到1000K。
附件中为各项能量结果，random为randomiztion的测试，AMgB-wav4是random后的波函数测试，md1,md2,md3是上述5,6,7的结果；附上最后一步md的输入文件。谢谢，求指点。

回复此楼