计算化学公社

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

[GROMACS] 求助:有关氢的力场修改

[复制链接 Copy URL]

76

帖子

0

威望

413

eV
积分
489

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
请教各位前辈老师,在论文中看到的氢的力场参数是这样的
H  +0.4932   0.0000   0.000
HCM  -0.9864  3.038  epsilon  0.2852  sigma

其中HCM是质心
但是我是在ATB网站下载的H的pdb和itp文件,里面只有两个氢原子。


请教各位我是不是需要自己改pdb文件和itp文件,加一个HCM质心并且对应上文献里的参数那?还是只需要把已有的两个氢原子的电荷分别改为+0.4932和-0.4932以及文献里的epsilon和sigma那?

1376

帖子

0

威望

3984

eV
积分
5360

Level 6 (一方通行)

2#
发表于 Post on 2022-6-25 21:07:59 | 只看该作者 Only view this author
oplsaa和ATB给出的gromos的形式是不一样的,一般来说,不通用。
又菜又爱玩

5万

帖子

99

威望

5万

eV
积分
112349

管理员

公社社长

3#
发表于 Post on 2022-6-25 23:56:51 | 只看该作者 Only view this author
发帖之后把自己的帖子好好看看,别留着那么多没解析的标签,看着极乱。重新编辑帖子把多余的标签删掉
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

76

帖子

0

威望

413

eV
积分
489

Level 3 能力者

4#
 楼主 Author| 发表于 Post on 2022-6-26 10:38:28 | 只看该作者 Only view this author
牧生 发表于 2022-6-25 21:07
oplsaa和ATB给出的gromos的形式是不一样的,一般来说,不通用。

格式可以自己改,就是只用他的电荷和lj参数是否有问题那?

76

帖子

0

威望

413

eV
积分
489

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2022-6-26 10:41:14 | 只看该作者 Only view this author
sobereva 发表于 2022-6-25 23:56
发帖之后把自己的帖子好好看看,别留着那么多没解析的标签,看着极乱。重新编辑帖子把多余的标签删掉

好的Sob老师,已经改好了,可以请教一下是否可以直接修改itp文件那?然后lj算好HCM和其他原子的直接加在nonbond params里那?

5万

帖子

99

威望

5万

eV
积分
112349

管理员

公社社长

6#
发表于 Post on 2022-6-27 04:54:38 | 只看该作者 Only view this author
我没看过你说的论文,应该是用三个点描述每个H2,每个氢原子核都有一个点,H2中心再加一个点。
自行创建个H2的[moleculetype],按照上述情况自己写拓扑信息。[atoms]里每个氢的电荷都是+0.4932,H2中心的HCM的电荷是-0.9864。给这两类粒子分别定义两种原子类型,范德华参数在[atomtypes]里设。对应氢核的粒子的范德华参数都为0,HCM的写成文献里给的sigma和epsilon。

并且必须令结构文件里每个H2也是靠三个点来描述,从而和[moleculetype]里的[atoms]里的粒子数和顺序相对应
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

76

帖子

0

威望

413

eV
积分
489

Level 3 能力者

7#
 楼主 Author| 发表于 Post on 2022-6-27 12:44:47 | 只看该作者 Only view this author
sobereva 发表于 2022-6-27 04:54
我没看过你说的论文,应该是用三个点描述每个H2,每个氢原子核都有一个点,H2中心再加一个点。
自行创建个 ...

懂了,谢谢sob老师,想问一下这样的三点氢分子的pdb文件是需要自己写吗?如果自己写的话是不是HCM的坐标位置在两个氢原子的之间吗?键连接是不是只需要写两个氢原子的就可以那?

5万

帖子

99

威望

5万

eV
积分
112349

管理员

公社社长

8#
发表于 Post on 2022-6-29 01:49:17 | 只看该作者 Only view this author
sj12138 发表于 2022-6-27 12:44
懂了,谢谢sob老师,想问一下这样的三点氢分子的pdb文件是需要自己写吗?如果自己写的话是不是HCM的坐标 ...


你那篇文章可能是把HCM作为虚拟点对待的。看gmx手册了解虚拟点是什么,以及有哪些规则根据实际原子位置产生虚拟点的位置
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

76

帖子

0

威望

413

eV
积分
489

Level 3 能力者

9#
 楼主 Author| 发表于 Post on 2022-6-29 15:48:16 | 只看该作者 Only view this author
sobereva 发表于 2022-6-29 01:49

你那篇文章可能是把HCM作为虚拟点对待的。看gmx手册了解虚拟点是什么,以及有哪些规则根据实际原子位 ...

好的谢谢sob老师

76

帖子

0

威望

413

eV
积分
489

Level 3 能力者

10#
 楼主 Author| 发表于 Post on 2022-7-2 12:52:58 | 只看该作者 Only view this author
sobereva 发表于 2022-6-29 01:49

你那篇文章可能是把HCM作为虚拟点对待的。看gmx手册了解虚拟点是什么,以及有哪些规则根据实际原子位 ...

老师您好,我想问一下,为什么gromos54a7力场无法加载tip4p的水模型呢?

tip4p.itp

1.12 KB, 下载次数 Times of downloads: 2

5万

帖子

99

威望

5万

eV
积分
112349

管理员

公社社长

11#
发表于 Post on 2022-7-3 03:47:17 | 只看该作者 Only view this author
sj12138 发表于 2022-7-2 12:52
老师您好,我想问一下,为什么gromos54a7力场无法加载tip4p的水模型呢?

二者根本毫无联系

“加载”语义不明
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

76

帖子

0

威望

413

eV
积分
489

Level 3 能力者

12#
 楼主 Author| 发表于 Post on 2022-7-3 10:19:45 | 只看该作者 Only view this author
本帖最后由 sj12138 于 2022-7-3 10:29 编辑
sobereva 发表于 2022-7-3 03:47
二者根本毫无联系

“加载”语义不明

就是top文件include进去之后top报错意思说这个tip4p是不能用的,然后pdb2gmx的交互里还没有tip4p的选项

5万

帖子

99

威望

5万

eV
积分
112349

管理员

公社社长

13#
发表于 Post on 2022-7-4 01:16:24 | 只看该作者 Only view this author
sj12138 发表于 2022-7-3 10:19
就是top文件include进去之后top报错意思说这个tip4p是不能用的,然后pdb2gmx的交互里还没有tip4p的选项

直接贴报错,别自己翻译

没有特殊情况,别用GROMOS+TIP4P这种莫名其妙的组合,TIP4P本来就很烂,而且GROMOS又不是结合TIP4P参数化的力场
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

22

帖子

0

威望

225

eV
积分
247

Level 3 能力者

14#
发表于 Post on 2024-9-4 16:18:26 | 只看该作者 Only view this author
sj12138 发表于 2022-7-3 10:19
就是top文件include进去之后top报错意思说这个tip4p是不能用的,然后pdb2gmx的交互里还没有tip4p的选项

你好,能分享下这篇关于氢力场的原文么

本版积分规则 Credits rule

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

GMT+8, 2024-11-23 05:06 , Processed in 0.178488 second(s), 24 queries , Gzip On.

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