| 查看: 740 | 回复: 0 | ||
[求助]
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的输入文件。谢谢，求指点。 |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : AMgB-md1.evp
- 附件 2 : AMgB-md2.evp
- 附件 3 : AMgB-md3.evp
- 附件 4 : AMgB-md3.in
- 附件 5 : AMgB-wav4.evp
2014-06-05 22:06:48, 113.28 K
2014-06-05 22:06:51, 113.28 K
2014-06-05 22:06:53, 113.28 K
2014-06-05 22:06:55, 11.93 K
2014-06-05 22:06:57, 17.9 K
» 猜你喜欢
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有8人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有8人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有6人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
» 本主题相关商家推荐: (我也要在这里推广)
