| 查看: 679 | 回复: 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
» 猜你喜欢
26/27申博自荐
已经有10人回复
-
东北林业大学材料科学与工程学院“一流”A+学科国家级人才团队课题组招收2026级博士生
已经有3人回复
-
医学类期刊求推荐
已经有5人回复
-
生活琐事由它去
已经有4人回复
-
提交了我也来说说感想
已经有12人回复
-
青B发送上会通知了吗
已经有9人回复
-
西安交大新媒学院副院长用撤稿论文结题
已经有6人回复
-
论文撤稿了
已经有8人回复
-
化学专业申博
已经有4人回复
-
某211大学教师把个人教师官方主页改成:我跑了我跑了我跑了!官宣跑路！
已经有5人回复
» 本主题相关商家推荐: (我也要在这里推广)
