计算化学公社

 找回密码 Forget password
 注册 Register
Views: 11314|回复 Reply: 7
打印 Print 上一主题 Last thread 下一主题 Next thread

[GROMACS] 如何在模拟过程中根据体系自由能实时改变系统温度?

[复制链接 Copy URL]

19

帖子

0

威望

197

eV
积分
216

Level 3 能力者

大概是每隔几纳秒就要改变热库的温度,此处温度不是模拟退火里的可以提前设定的温度,而是要根据模拟完成后体系自由能来计算,请问有没有什么好方法?

6万

帖子

99

威望

6万

eV
积分
125148

管理员

公社社长

2#
发表于 Post on 2021-4-8 03:22:45 | 只看该作者 Only view this author
只能自己写个shell脚本,反复产生输入文件、运行任务、根据提取的数据自动修改mdp
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

19

帖子

0

威望

197

eV
积分
216

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2021-4-8 15:47:05 | 只看该作者 Only view this author
sobereva 发表于 2021-4-8 03:22
只能自己写个shell脚本,反复产生输入文件、运行任务、根据提取的数据自动修改mdp

好吧,谢谢sob老师,
这个方法太慢了,今天和导师讨论了一下,认为这个方法太慢了,决定换成LAMMPS,用里边现成的GCE包

6万

帖子

99

威望

6万

eV
积分
125148

管理员

公社社长

4#
发表于 Post on 2021-4-10 06:14:07 | 只看该作者 Only view this author
怎么就太慢了?
gromacs计算速度远远快于lammps,自己写个脚本自动做这个事(也就超不过20行),每次自动重新产生tpr文件然后继续开始计算也就耗几秒钟的时间,不比换成lammps快得多得多?假设一天跑100ns,重新跑的次数也就几十次,也就多花一分钟额外的耗时
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

19

帖子

0

威望

197

eV
积分
216

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2021-4-10 13:54:29 | 只看该作者 Only view this author
sobereva 发表于 2021-4-10 06:14
怎么就太慢了?
gromacs计算速度远远快于lammps,自己写个脚本自动做这个事(也就超不过20行),每次自动 ...

sob老师,是这样的,导师修改了我的想法,是每隔几十步甚至几步(步长2fs)就要修改系统温度,lammps上实现了每步改变温度的包,但是gromacs还没有

313

帖子

2

威望

3900

eV
积分
4253

Level 6 (一方通行)

6#
发表于 Post on 2021-4-10 14:47:17 | 只看该作者 Only view this author
fs级别的改变温度不是shell脚本能够干的事,也不现实。除非自己有能力改gmx源码,否则另辟蹊径

19

帖子

0

威望

197

eV
积分
216

Level 3 能力者

7#
 楼主 Author| 发表于 Post on 2021-4-11 23:21:29 | 只看该作者 Only view this author
lyj714 发表于 2021-4-10 14:47
fs级别的改变温度不是shell脚本能够干的事,也不现实。除非自己有能力改gmx源码,否则另辟蹊径

嗯嗯,我改用lammps了,谢谢您

71

帖子

0

威望

1148

eV
积分
1219

Level 4 (黑子)

8#
发表于 Post on 2021-4-12 15:56:46 | 只看该作者 Only view this author
gromacs退火模拟可以控温,多长时间改变温度

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-22 23:51 , Processed in 0.171528 second(s), 20 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list