24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2935  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 304求调剂 +7 castLight 2026-04-16 7/350 2026-04-17 20:05 by 关一盏灯cd
[考研] 295分求调剂 +5 ?要上岸? 2026-04-17 5/250 2026-04-17 16:51 by fenglj492
[考研] 322求调剂 +6 tekuzu 2026-04-17 6/300 2026-04-17 13:48 by Espannnnnol
[考研] 收到复试调剂但是去不了 +6 小蜗牛* 2026-04-16 6/300 2026-04-17 10:05 by 涵竹刘
[考研] 297,工科调剂? +4 河南农业大学-能 2026-04-14 4/200 2026-04-16 22:52 by wulijun2012
[考研] 一志愿沪9,生物学326求调剂 +9 刘墨墨 2026-04-15 9/450 2026-04-16 17:14 by 崔崔崔cccc
[考研] 291求调剂 +11 关忆北. 2026-04-14 11/550 2026-04-16 15:18 by jiahl2024
[考博] 申博自荐 +3 Linxia林夏 2026-04-13 3/150 2026-04-16 12:55 by 墨荷之露
[基金申请] RY:中国产出的科学垃圾论文,绝对数量和比例都世界第一 +7 zju2000 2026-04-14 18/900 2026-04-16 11:36 by 欢乐颂叶蓁
[考研] 279学硕食品专业求调剂院校 20+7 孤独的狼爱吃羊 2026-04-12 29/1450 2026-04-16 09:00 by screening
[考研] 求调剂推荐 +8 小聂爱学习 2026-04-14 8/400 2026-04-16 07:22 by 学员JpLReM
[考研] 289 分105500药学专硕求调剂(找B区学校) +4 白云123456789 2026-04-13 4/200 2026-04-16 00:18 by 粉沁若尘
[考研] 求调剂 +11 小聂爱学习 2026-04-11 15/750 2026-04-15 21:57 by noqvsozv
[考研] 一志愿A区211,22408 321求调剂 +6 随心所欲☆ 2026-04-15 7/350 2026-04-15 21:45 by lbsjt
[考研] 085404 22408 309分求调剂 +9 lzmk 2026-04-14 10/500 2026-04-15 20:02 by 学员JpLReM
[考研] 材料工程281还有调剂机会吗 +43 xaw. 2026-04-11 44/2200 2026-04-15 12:46 by 西北望—风沙
[考研] 366求调剂 +11 不知名的小卅 2026-04-11 11/550 2026-04-14 15:50 by zs92450
[考研] 考研调剂 +13 长弓傲 2026-04-13 14/700 2026-04-14 14:44 by zs92450
[考研] 085408光电信息工程专硕355一志愿长春光机所调剂 +6 王ymaa 2026-04-13 13/650 2026-04-14 11:33 by 王ymaa
[考研] 调剂结束 +6 floriea 2026-04-12 8/400 2026-04-12 18:13 by zhouxiaoyu
信息提示
请填处理意见