计算化学公社

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

[GROMACS] 求助:关于限制部分分子间距离的设置

[复制链接 Copy URL]

116

帖子

0

威望

1091

eV
积分
1207

Level 4 (黑子)

各位老师好,我之前看到了一篇文章是限制蛋白质中某两个特定原子的相对距离进行MD模拟,想重现这种方法,因此想请教各位老师gromacs有没有办法限制我选择的原子的距离进行模拟呢?(constraint限制,距离改变时受到额外的力惩罚)

69

帖子

1

威望

982

eV
积分
1071

Level 4 (黑子)

2#
发表于 Post on 2021-10-14 12:27:50 | 只看该作者 Only view this author
相对距离的话,不知道能不能当作键来处理,比如在[bond]项里面写一个较大的距离,然后用简谐力来控制距离的变化。你可以拿一个小体系试试这个方法,看体系会不会崩。
信息不灭,熵永增

6万

帖子

99

威望

6万

eV
积分
125131

管理员

公社社长

3#
发表于 Post on 2021-10-14 13:02:43 | 只看该作者 Only view this author
你搞混了constraint和restraint,看此文
解析gromacs的restraint、constraint和freeze
http://sobereva.com/10
你要求的是restraint。constraint是直接通过约束算法保持原子间距离精确不变

添加[distance_restraint]字段,或者用functype=6或10的bond项可以实现。
北京科音自然科学研究中心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

116

帖子

0

威望

1091

eV
积分
1207

Level 4 (黑子)

4#
 楼主 Author| 发表于 Post on 2021-10-14 21:07:50 | 只看该作者 Only view this author
EdwardLimit 发表于 2021-10-14 12:27
相对距离的话,不知道能不能当作键来处理,比如在项里面写一个较大的距离,然后用简谐力来控制距离的变化。 ...

谢谢老师

116

帖子

0

威望

1091

eV
积分
1207

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2021-10-14 21:09:03 | 只看该作者 Only view this author
sobereva 发表于 2021-10-14 13:02
你搞混了constraint和restraint,看此文
解析gromacs的restraint、constraint和freeze
http://sobereva.c ...

明白了,谢谢sob老师指点

226

帖子

0

威望

878

eV
积分
1104

Level 4 (黑子)

6#
发表于 Post on 2021-10-14 21:25:27 | 只看该作者 Only view this author
sobereva 发表于 2021-10-14 13:02
你搞混了constraint和restraint,看此文
解析gromacs的restraint、constraint和freeze
http://sobereva.c ...

我看手册中[distance_restraint]用的函数为1,且有一个fac项(力场数的系数),请教老师只要一个fac但是力场数在[[distance_restraint]]并没有设置,实在哪里设置或有什么规则吗。

functype=6或10使只能[intermolecular_interactions]中设置吗。还是说在[distance_restraint]也可以用

6万

帖子

99

威望

6万

eV
积分
125131

管理员

公社社长

7#
发表于 Post on 2021-10-15 05:14:39 | 只看该作者 Only view this author
sun666 发表于 2021-10-14 21:25
我看手册中[distance_restraint]用的函数为1,且有一个fac项(力场数的系数),请教老师只要一个fac但是 ...

对于用[distance_restraint],力常数通过mdp里的disre-fc设置,而[ distance_restraint ]里面每两个原子之间还可以设fac(multiplication factor),它被乘到disre-fc上作为这一对原子的实际用的力常数

functype=6、10显然是在[moleculetype]或[intermolecular_interactions]中的[bonds]里写
北京科音自然科学研究中心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

226

帖子

0

威望

878

eV
积分
1104

Level 4 (黑子)

8#
发表于 Post on 2021-10-15 13:55:37 | 只看该作者 Only view this author
sobereva 发表于 2021-10-15 05:14
对于用[distance_restraint],力常数通过mdp里的disre-fc设置,而[ distance_restraint ]里面每两个原子 ...

对于[distance_restraint]使用的为函数1形式(需要定义的参数为b0以及Kb),但是手册给的例子虽然用的函数1 形式,但是要设置 low   up1   up2这三个数值,我怎么觉得这个是函数10平底势需要设置的值呢。
到底在[distance_restraint]中应设置的是平衡键长还是平地势需要的low   up1   up这三个数值呢

6万

帖子

99

威望

6万

eV
积分
125131

管理员

公社社长

9#
发表于 Post on 2021-10-16 03:19:52 | 只看该作者 Only view this author
sun666 发表于 2021-10-15 13:55
对于[distance_restraint]使用的为函数1形式(需要定义的参数为b0以及Kb),但是手册给的例子虽然用的函数 ...

如果不需要平底势,用functype=6的[bonds]项就完了
北京科音自然科学研究中心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

226

帖子

0

威望

878

eV
积分
1104

Level 4 (黑子)

10#
发表于 Post on 2021-10-16 11:58:16 | 只看该作者 Only view this author
sobereva 发表于 2021-10-16 03:19
如果不需要平底势,用functype=6的项就完了

在[intermolecular_interaction]中的functype=6的[bonds]项,键的力常数使用陈老师的forcefit算。是不是这样与amber的MCPB达到的效果完全一致了。

6万

帖子

99

威望

6万

eV
积分
125131

管理员

公社社长

11#
发表于 Post on 2021-10-16 18:32:53 | 只看该作者 Only view this author
sun666 发表于 2021-10-16 11:58
在中的functype=6的项,键的力常数使用陈老师的forcefit算。是不是这样与amber的MCPB达到的效果完全一致 ...

MCPB我了解不多,不好说
光是一个bond项未必能完全达到目的(可以先跑跑试试),可能还得添加angle项
北京科音自然科学研究中心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

114

帖子

0

威望

1034

eV
积分
1148

Level 4 (黑子)

12#
发表于 Post on 2022-2-5 08:49:45 | 只看该作者 Only view this author
sobereva 发表于 2021-10-14 13:02
你搞混了constraint和restraint,看此文
解析gromacs的restraint、constraint和freeze
http://sobereva.c ...

老师,您在这篇博文里提到
“[ constraints ]
1 8 2 0.153    //原子1 原子2 type 约束的距离(0.153nm)。
之后1 2原子距离会固定保持在0.153nm,即便1 2原子在[ bonds ]中定义是成键的也是如此。
注意,.mdp里如果写constraint=none,只是说批量约束hbonds、all-bonds等没有了,在拓扑文件中自定义的constraint约束项仍然生效。”
我想问一下,此处的数字1和2代表原子序号1和2的话,那么8是指的什么呢?

6万

帖子

99

威望

6万

eV
积分
125131

管理员

公社社长

13#
发表于 Post on 2022-2-5 21:39:19 | 只看该作者 Only view this author
ljc050512 发表于 2022-2-5 08:49
老师,您在这篇博文里提到
“[ constraints ]
1 8 2 0.153    //原子1 原子2 type 约束的距离(0.153nm) ...

1和8是原子序号
2是constraint的类型,决定是否影响键连关系,手册里写了


北京科音自然科学研究中心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

114

帖子

0

威望

1034

eV
积分
1148

Level 4 (黑子)

14#
发表于 Post on 2022-2-5 23:51:05 | 只看该作者 Only view this author
sobereva 发表于 2022-2-5 21:39
1和8是原子序号
2是constraint的类型,决定是否影响键连关系,手册里写了

谢谢老师的指导。

本版积分规则 Credits rule

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

GMT+8, 2026-2-20 07:25 , Processed in 0.173090 second(s), 22 queries , Gzip On.

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