| 查看: 2450 | 回复: 5 | |||
[交流]
PWmat计算NEB的详细教程
|
|||
|
PWmat计算NEB的过程，这里以Si64B1体系为例说明： 1.对于NEB计算，首先建立两个文件夹分别命名initial（初态的意思），final（末态的意思），每个文件夹放入PWnat计算的必备文件（etot.input，UPF赝势文件，atom.config），initail和final文件的atom.config对应未优化的初末态，优化后得到稳定的结构文件final.config。要保证二者优化时除结构文件以为的其他文件一致。etot.input的设置如下 http://muchongimg.xmcimg.com/data/bcs/2017/0727/w365h6383842_1501145275_672.jpg 下图中的（a），（b）分别是初末态优化后的结构 http://muchongimg.xmcimg.com/data/bcs/2017/0727/bw304h6383842_1501145303_920.jpg 2.NEB计算，在当前目录下放入UPF赝势文件，将初态优化后的结构final.config拷贝到当前目录并命名为atom.config，将末态优化后的结构文件final.config也拷贝到当前目录并命名为atom_final.config，并在etot.input设置参数“NEB_DETAIL=IMTH，NSTEP，FORCE_TOL，NIMAGE，AK，TYPE_STRING，E0，EN，ITYPE_AT2，ATOM2.CONFIG“。其中第一个参数表示原子弛豫方法；对于NEB计算推荐用3；NSTEP表示最大的优化步数；FORCE_TOL表示力的收敛标准；NIMAGE表示插入的点数；AK表示弹性力的大小；TYPE_STRING表示弹性力的类型，PWmat提供了两种类型的弹性力供选择（TYPE_SPRING=1,2），而且弹性力的大小可以通过参数AK来调节。对于复杂的体系，我们首先选择TYPE_SPRING为2，并且使用了较大的ak，该设置有利于计算的收敛，但是计算结果准确性不高。为了得到准确的势垒值，需要把TYPE_SPRING设置为1或者调小ak，该体系通过将TYPE_SPRING设置为1，并将TYPE_SPRING为2得到的收敛结构作为再次NEB计算的初始结构。E0和EN分别是初末态的能量；ATOM2.CONFIG表示末态优化后的结构文件。etot.input的设置如下 http://muchongimg.xmcimg.com/data/bcs/2017/0727/bw333h6383842_1501146837_486.jpg 3.数据处理 计算结束后，得到NEB.BARRIER文件，其记录了每步优化后，各个结构的能量，具体如下图所示，其中第二列表示能量，第三列表示两个点之间的距离，第四列表示角度 http://muchongimg.xmcimg.com/data/bcs/2017/0727/w385h6383842_1501147550_409.jpg 可以利用上述数据画出扩散曲线，如下图所示 http://muchongimg.xmcimg.com/data/bcs/2017/0727/bw142h6383842_1501147640_340.jpg  initail.jpg  tot.jpg  nebetot.jpg  BARRIER.jpg  d.jpg |
回复此楼» 猜你喜欢
小木虫没落了，除了祈祷帖子，几乎看不到有价值的帖子
已经有15人回复
-
面上的水很深
已经有18人回复
-
国社科通讯评审
已经有5人回复
-
5个%2F，非固定段22个字母
已经有9人回复
-
生命口C13面上会评时间
已经有11人回复
-
内心匮乏
已经有16人回复
-
生命口面上会评结束了吗大家有知道结果吗？现在还没有结果正常不？
已经有9人回复
-
投英文期刊如何查重
已经有6人回复
-
不要再数国自然申请书的 filecode 的分隔符个数了
已经有21人回复
-
伯胺上甲基求助
已经有7人回复
» 抢金币啦！回帖就可以得到:
-
工材电工至今没任何消息，散金币以求心安吧~
+2/1136
-
北京工业大学材料科学与工程学院李建荣教授课题组诚聘电池相关领域应届毕业生和博士后
+1/578
-
坐标山东，诚征女友
+1/411
-
苏州B+—>C轮 胶带/胶黏剂公司 双休 食宿免 空调洁净室聘胶带研发
+3/249
-
重磅分享｜全链路AI材料研发闭环平台PULSE即将发布，重构新材料研发范式
+2/126
-
年薪40W+项目论文绩效，广州医科大学招收骨科、生物材料相关背景博士后2名
+1/81
-
哈工大医康学院材料模拟计算方向人才（全职正高、副高和博士后人员）招聘
+1/78
-
北京理工大学杰青沈国震团队招聘柔性智能电子博士后
+1/33
-
祈福，求顺利！
+1/33
-
西安交通大学补亚忠课题组招聘博士生
+1/28
-
希望面上有好消息:D:D:D
+2/20
-
中国药科大学韩凯课题组博士后及博士招聘启事（药剂/生物医用材料/免疫治疗等方向）
+2/8
-
人生失去与收获
+1/8
-
双一流高校招收博士（现代技术管理/环境技术经济管理）及硕士研究生（长期有效）
+1/6
-
诚邀生物医用材料领域学者加入Biofunctional Materials青年编委（Scopus收录）
+1/5
-
双一流高校招收博士（现代技术管理/环境技术经济管理）及硕士研究生（长期有效）
+1/5
-
东北大学理学院力学李健课题组招收2027级硕士推免研究生
+1/2
-
福州大学吴衔誉教授课题组 诚招2027届博士生
+1/2
-
上海交通大学医学院附属第一人民医院陈思思课题组招聘
+1/1
-
北理工集成电路杰青团队 | 诚招科助理
+1/1
3楼2020-05-30 11:12:00
4楼2020-06-21 10:46:10
5楼2020-07-30 10:43:48
6楼2020-07-30 10:43:57
简单回复
dsctg2楼
2017-07-27 17:43
回复
HZ_Tian(金币+1): 谢谢参与
