24小时热门版块排行榜    

查看: 2884  |  回复: 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的回帖

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的回帖
查看全部 12 个回答

wangz00

禁虫 (正式写手)


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

2楼2016-02-19 17:03:15
已阅   回复此楼   关注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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿山东大学药学学硕求调剂 +3 开开心心没烦恼 2026-03-23 4/200 2026-03-24 00:06 by 开开心心没烦恼
[考研] 341求调剂(一志愿湖南大学070300) +5 番茄头--- 2026-03-22 6/300 2026-03-23 23:45 by Txy@872106
[考研] 材料292调剂 +4 橘颂思美人 2026-03-23 4/200 2026-03-23 23:16 by peike
[考研] 材料专硕找调剂 +3 哈哈哈吼吼吼哈 2026-03-23 3/150 2026-03-23 23:13 by peike
[考研] 一志愿国科过程所081700,274求调剂 +3 三水研0水立方 2026-03-23 3/150 2026-03-23 23:11 by MajorWen
[考研] 一志愿陕师大生物学071000,298分,求调剂 +3 SYA! 2026-03-23 3/150 2026-03-23 19:09 by macy2011
[考研] 350求调剂 +6 weudhdk 2026-03-19 6/300 2026-03-23 15:47 by tangyuan0840221
[考研] 070300,一志愿北航320求调剂 +3 Jerry0216 2026-03-22 5/250 2026-03-23 09:16 by 。。堂堂
[考研] 317求调剂 +12 申子申申 2026-03-19 18/900 2026-03-22 22:23 by luoyongfeng
[考研] 求调剂 +4 要好好无聊 2026-03-21 4/200 2026-03-21 18:57 by 学员8dgXkO
[考研] 296求调剂 +4 www_q 2026-03-20 4/200 2026-03-21 17:26 by 学员8dgXkO
[考研] 268求调剂 +9 简单点0 2026-03-17 9/450 2026-03-21 15:37 by lature00
[考研] 332求调剂 +3 凤凰院丁真 2026-03-20 3/150 2026-03-21 10:27 by luoyongfeng
[考研] 085700资源与环境308求调剂 +12 墨墨漠 2026-03-18 13/650 2026-03-21 01:42 by JourneyLucky
[考研] 一志愿西南交大,求调剂 +5 材化逐梦人 2026-03-18 5/250 2026-03-21 00:26 by JourneyLucky
[考研] 296求调剂 +6 www_q 2026-03-18 10/500 2026-03-20 23:56 by JourneyLucky
[考研] 一志愿武汉理工材料工程专硕调剂 +9 Doleres 2026-03-19 9/450 2026-03-20 22:36 by JourneyLucky
[考研] 261求B区调剂,科研经历丰富 +3 牛奶很忙 2026-03-20 4/200 2026-03-20 19:34 by JourneyLucky
[考研] 086500 325 求调剂 +3 领带小熊 2026-03-19 3/150 2026-03-20 18:38 by 尽舜尧1
[考研] 【同济软件】软件(085405)考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
信息提示
请填处理意见