| 查看: 2920 | 回复: 11 | ||||
| 当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖 | ||||
chenxin1992银虫 (正式写手)
|
[交流]
最小能量路径搜索软件MEPsearcher发布 已有4人参与
|
|||
|
简介 MEPsearcher基于string计算方法,用于搜索二维全局离散势能面上的能量最小路径Mminimun Energy Path(MEP),通常MEP可以完整揭示反应机理。过渡态、活化能等信息可以从MEP中获取。能量最低路径的搜索可以利用于蛋白质折叠、化学反应、相变等过程的机理研究。 下载地址: https://github.com/chenxin199261/MEPSearcher 具体信息请参考: http://chen-xin.tk/blog/index.php/archives/44/ ------ 使用方法 1. 准备文件 你只需要准备两个输入文件,势能面文件PES.data,初猜生成文件guess.data。 势能面文件写成三列X Y Energy即可。可参考test文件夹中PES.data格式 准备初始猜测需要对势能面有一定的直观感觉,给出起始点、终止点以及初猜珠子(beads)的个数。 2.参数设定 需要进入MEPsearcher.py中对参数进行设定,只有三个参数需要测定,分别是最大迭代次数,收敛限以及优化步长。 最大迭代次数,Max_iter默认1000步,可以根据具体需要调大或者调小。 收敛限,用与上次优化结果的距离差加和表示。Tol默认1e-5,beads过多,可以适当调小、调大这个值。 优化步长值,h,是最速下降优化方法一个必需参数。这个值的设定麻烦些。过小导致优化缓慢,过大导致直接优化到势能面错误的区域。需要使用者反复调试。希望在下一个版本中,能启用其它不需要步长的优化方法。 3.运行 赋予MEPsearch.py执行权限 chmod +x MEPsearcher.py 直接运行 ./MEPsearcher.py 计算过程每隔50步迭代,会打印在屏幕上。 运行结束后,产生out.data文件。 这个文件存储了MEP的最终结果。分三列表示,分别为: X Y Energy --------- 运行结果 Graph1.jpg |
» 收录本帖的淘帖专辑推荐
多孔材料 |
» 猜你喜欢
淮北师范大学 化学与化工科学学院 电催化课题组还有1名硕士研究生招生名额
已经有0人回复
化学相关硕士研究生招收调剂-湖北师范大学 (快来捡漏,最新调剂信息,过国家线即可)
已经有1人回复
物理化学论文润色/翻译怎么收费?
已经有289人回复
河南理工大学化学化工学院招收学硕调剂生
已经有0人回复
黑龙江科技大学化工 环境工程调剂
已经有22人回复
图书求助:Fluorescent Probes for Bioactive Species: Design
已经有0人回复
莆田学院环境与生物工程学院资源与环境专硕(环境工程方向) 开放调剂招生
已经有12人回复
浙江海洋大学石油化工与环境学院招收工科调剂生(学硕、专硕均可,线上面试)
已经有0人回复
浙江海洋大学石油化工学院接收工科调剂(学硕、专硕均可,近期线上面试)
已经有8人回复
浙江海洋大学石油与天然气工程专业接收工科调剂(13日早上截止,线上面试)
已经有0人回复
考研调剂~湖南工程学院材料与化工学院
已经有0人回复
chenxin1992
银虫 (正式写手)
- 应助: 18 (小学生)
- 金币: 1030.1
- 散金: 126
- 红花: 11
- 帖子: 337
- 在线: 103.1小时
- 虫号: 1813050
- 注册: 2012-05-12
- 性别: GG
- 专业: 无机合成和制备化学
3楼2016-02-19 19:04:11
★
小木虫: 金币+0.5, 给个红包,谢谢回帖
小木虫: 金币+0.5, 给个红包,谢谢回帖
|
本帖内容被屏蔽 |
2楼2016-02-19 17:03:15
chenxin1992
银虫 (正式写手)
- 应助: 18 (小学生)
- 金币: 1030.1
- 散金: 126
- 红花: 11
- 帖子: 337
- 在线: 103.1小时
- 虫号: 1813050
- 注册: 2012-05-12
- 性别: GG
- 专业: 无机合成和制备化学
5楼2016-02-19 19:49:14
virtualzx
木虫 (著名写手)
- QC强帖: 2
- 应助: 263 (大学生)
- 金币: 7161.3
- 红花: 54
- 帖子: 1605
- 在线: 317.6小时
- 虫号: 2069080
- 注册: 2012-10-18
- 性别: GG
- 专业: 理论和计算化学
6楼2016-02-20 14:04:19














回复此楼