| 查看: 716 | 回复: 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
» 猜你喜欢
今天放榜吗？
已经有12人回复
-
filecode，4个jtjc了
已经有9人回复
-
今天基金会出结果吗？20260819
已经有8人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有115人回复
-
2027广东省杰青
已经有10人回复
-
朋友圈看到的
已经有8人回复
-
时间戳变了，能看出什么问题？
已经有14人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有30人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有7人回复
-
什么时候开奖？
已经有6人回复
» 本主题相关商家推荐: (我也要在这里推广)
