24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 2923  |  回复: 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

---------

运行结果

最小能量路径搜索软件MEPsearcher发布
Graph1.jpg
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangz00

禁虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
本帖内容被屏蔽

2楼2016-02-19 17:03:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chenxin1992

银虫 (正式写手)

引用回帖:
2楼: Originally posted by wangz00 at 2016-02-19 17:03:15
你好,请问你的图片的PES是解析的吧?有没有将这个程序用到真实的多维体系?

不是解析的,是离散的。
多维的可以,不过意义不大,因为多维没法作图。

可以做跟其他软件的接口,找到反应路径上结构,不过这个就是另一个问题了。
目前问题是优化算法写的不好,迭代收敛慢,调用gaussian,vasp之类的很耗时。正在改善
3楼2016-02-19 19:04:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangz00

禁虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
本帖内容被屏蔽

4楼2016-02-19 19:47:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chenxin1992

银虫 (正式写手)

引用回帖:
4楼: Originally posted by wangz00 at 2016-02-19 19:47:53
有时间会试一下,多谢分享...

我的是测试版本,只测试了几个势能面。
如果有使用不通的地方,可以随时联系我,我好改进他
5楼2016-02-19 19:49:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

virtualzx

木虫 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
和一般string method的主要差别在哪里呢

发自小木虫IOS客户端
6楼2016-02-20 14:04:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chenxin1992

银虫 (正式写手)

引用回帖:
6楼: Originally posted by virtualzx at 2016-02-20 14:04:19
和一般string method的主要差别在哪里呢

没有区别的。我的是string方法最简单的实现,最速下降加上差值重分布。

发自小木虫Android客户端
7楼2016-02-21 17:02:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

virtualzx

木虫 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
一定要点阵输入吗?高维度体系不大现实啊

发自小木虫IOS客户端
8楼2016-02-22 01:05:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chenxin1992

银虫 (正式写手)

引用回帖:
8楼: Originally posted by virtualzx at 2016-02-22 01:05:48
一定要点阵输入吗?高维度体系不大现实啊

是的,高纬度也没有什么好办法。解析的表达式改改程序就可以

发自小木虫Android客户端
9楼2016-02-22 15:59:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hkhtlxf

新虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
内容已删除
10楼2018-05-08 18:23:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 chenxin1992 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085600材料与化工,求调剂 +14 won_qii 2026-04-07 14/700 2026-04-13 22:21 by pies112
[考研] 22专硕求调剂 +8 haoyun上岸 2026-04-11 10/500 2026-04-13 22:15 by zhq0425
[考研] 300分求调剂 (085501机械专硕,本科扬大) +8 xu@841019 2026-04-11 9/450 2026-04-13 21:36 by xu@841019
[考研] 一志愿沪9,326求生物学调剂 +9 刘墨墨 2026-04-13 9/450 2026-04-13 18:59 by lbsjt
[考研] 0854调剂 +10 长弓傲 2026-04-12 13/650 2026-04-13 17:32 by 长弓傲
[考研] 一志愿2110,化学学硕310分,本科重点双非求调剂 +20 努力奋斗112 2026-04-08 20/1000 2026-04-13 14:24 by 张zhihao
[考研] 化工学硕 285求调剂 +29 Wisjxn 2026-04-07 29/1450 2026-04-13 13:47 by floriea
[考研] 211本科材料化工求调剂 +17 YHLAH 2026-04-11 20/1000 2026-04-13 08:46 by Sammy2
[考研] 求调剂,一志愿大连理工大学354分 +5 雨声余生 2026-04-11 6/300 2026-04-11 16:12 by 雨声余生
[考研] 085600材料与化工329分求调剂 +16 叶zilin 2026-04-10 16/800 2026-04-11 11:04 by may_新宇
[考研] 298求调剂 +5 残荷新柳 2026-04-07 5/250 2026-04-11 11:02 by 紫曦紫棋
[考研] 085500求调剂材料 +10 易11122 2026-04-09 10/500 2026-04-11 10:39 by maddjdld
[考研] 281求调剂 +11 觉得好的吧 2026-04-10 11/550 2026-04-11 09:35 by 逆水乘风
[考研] 22408调剂求助 +7 毂12 2026-04-09 9/450 2026-04-11 09:23 by 哦哦123
[考研] 材料专硕调剂 +16 哈哈哈吼吼吼哈 2026-04-07 17/850 2026-04-09 21:16 by wutongshun
[考研] 283电子信息求调剂 +4 三石WL 2026-04-08 4/200 2026-04-09 10:21 by wp06
[考研] 考研求调剂 +4 雯??? 2026-04-08 4/200 2026-04-08 21:44 by 土木硕士招生
[考研] 388求调剂 +6 四川王涛 2026-04-07 8/400 2026-04-08 00:17 by JourneyLucky
[考研] 307求调剂 +3 Youth@@ 2026-04-07 3/150 2026-04-07 22:00 by hemengdong
[考研] 318求调剂 +5 李青山山山 2026-04-07 5/250 2026-04-07 18:24 by 蓝云思雨
信息提示
请填处理意见