计算化学公社

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

[GROMACS] 有机分子模拟通过分子间作用自组装时力场的选取

[复制链接 Copy URL]

24

帖子

0

威望

647

eV
积分
671

Level 4 (黑子)

之前用gromacs模拟的体系比较固定,在力场选取的时候都是参考文献。最近做了一些尝试性的工作,在力场选取遇到了问题。模拟的是有机分子进行共组装的一个体系,主要是看生成团簇的形貌和进行氢键、π-π之前的弱作用分析。选取哪种力场比较合适,原子包括CHONS五种。gromos54a7 和amber99sb-ildn哪个更合适?原子电荷用resp

1560

帖子

0

威望

4995

eV
积分
6555

Level 6 (一方通行)

2#
发表于 Post on 2021-9-8 20:38:10 | 只看该作者 Only view this author
amber是全原子力场,gromos为联合原子力场。。所以原则上,对于CH3,CH2基团居多的体系,amber力场精度高于gromos力场。

但对于氢键,π-π的弱作用分析,感觉这两种力场区别不是很大。
又菜又爱玩

24

帖子

0

威望

647

eV
积分
671

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2021-9-8 20:39:36 | 只看该作者 Only view this author
牧生 发表于 2021-9-8 20:38
amber是全原子力场,gromos为联合原子力场。。所以原则上,amber力场精度高于gromos力场。

对于氢键,π ...

嗯,谢谢您,对于氢键,π-π的弱作用分析,比较适合的力场都有什么呢?

1560

帖子

0

威望

4995

eV
积分
6555

Level 6 (一方通行)

4#
发表于 Post on 2021-9-8 20:50:16 | 只看该作者 Only view this author
zwn 发表于 2021-9-8 20:39
嗯,谢谢您,对于氢键,π-π的弱作用分析,比较适合的力场都有什么呢?

针对低浓度有机小分子水溶液,个人觉得在力场的精度方面,amber≈opls力场,其次gromos。

愚见,不一定正确。
又菜又爱玩

226

帖子

0

威望

878

eV
积分
1104

Level 4 (黑子)

5#
发表于 Post on 2021-9-8 23:23:03 | 只看该作者 Only view this author
牧生 发表于 2021-9-8 20:50
针对低浓度有机小分子水溶液,个人觉得在力场的精度方面,amber≈opls力场,其次gromos。

愚见,不一 ...

请问前辈gaff力场如何

1171

帖子

7

威望

6858

eV
积分
8169

Level 6 (一方通行)

6#
发表于 Post on 2021-9-8 23:56:02 | 只看该作者 Only view this author
你确定你说的是有机分子?你说的这两种力场都不能用于有机分子吧...

有机分子要用OPLS,GAFF,CGENFF这几种力场

6万

帖子

99

威望

6万

eV
积分
125128

管理员

公社社长

7#
发表于 Post on 2021-9-9 01:19:18 | 只看该作者 Only view this author
amber99sb-ildn面向生物分子体系,根本就不适合小分子。小分子有专用的GAFF

老老实实用acpype产生GAFF(实质上是GAFF2)的就完了,方便省事
北京科音自然科学研究中心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

24

帖子

0

威望

647

eV
积分
671

Level 4 (黑子)

8#
 楼主 Author| 发表于 Post on 2021-9-9 08:45:30 | 只看该作者 Only view this author
fhh2626 发表于 2021-9-8 23:56
你确定你说的是有机分子?你说的这两种力场都不能用于有机分子吧...

有机分子要用OPLS,GAFF,CGENFF这 ...

谢谢

24

帖子

0

威望

647

eV
积分
671

Level 4 (黑子)

9#
 楼主 Author| 发表于 Post on 2021-9-9 08:47:09 | 只看该作者 Only view this author
sobereva 发表于 2021-9-9 01:19
amber99sb-ildn面向生物分子体系,根本就不适合小分子。小分子有专用的GAFF

老老实实用acpype产生GAFF( ...

谢谢sob老师

349

帖子

0

威望

905

eV
积分
1258

Level 4 (黑子)

10#
发表于 Post on 2022-10-24 22:39:43 | 只看该作者 Only view this author
sobereva 发表于 2021-9-9 01:19
amber99sb-ildn面向生物分子体系,根本就不适合小分子。小分子有专用的GAFF

老老实实用acpype产生GAFF( ...

老师,您好,我想模拟有机分子在甲醇溶液中的自组装过程。然后有机分子的参数通过ambertools+acpype获得,那么甲醇溶剂的参数能否也通过这个方式获得?在进行packmol搭建初始模型时,甲醇添加的数目或者位置是否有要求,还是我可以作为单独一项按照下面方式进行添加呢?辛苦老师了
tolerance 2.0
add_box_sides 1.2
output mix.pdb
structure TB.pdb
number 100
end structure
structure methanol.pdb
  number 1000
end structure

6万

帖子

99

威望

6万

eV
积分
125128

管理员

公社社长

11#
发表于 Post on 2022-10-24 23:51:34 | 只看该作者 Only view this author
xxzj 发表于 2022-10-24 22:39
老师,您好,我想模拟有机分子在甲醇溶液中的自组装过程。然后有机分子的参数通过ambertools+acpype获得 ...

ambertools+acpype已经过时了,用sobtop(http://sobereva.com/soft/Sobtop)方便灵活普适得多。显然也适用于甲醇这么普通的分子

让packmol随机放置甲醇,填满一个不是很小的盒子就完了。
北京科音自然科学研究中心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

349

帖子

0

威望

905

eV
积分
1258

Level 4 (黑子)

12#
发表于 Post on 2022-10-25 18:21:30 | 只看该作者 Only view this author
sobereva 发表于 2022-10-24 23:51
ambertools+acpype已经过时了,用sobtop(http://sobereva.com/soft/Sobtop)方便灵活普适得多。显然也适 ...

老师,还需要再请教您三个问题,辛苦您
(1)水分子的参数文件是否也可以通过sobtop进行获取呢?
(2)进行有机小分子在溶剂中自组装过程是否和正常的gromacs模拟一样呢(packmol搭建模型随时放置溶剂分子和有机小分子,我目前小分子放了80个,甲醇让其填满是否就可以了)
(3)进行有机分子在溶剂中的自组装模拟,最后md模拟之后是否还需要其它操作去除盒子中的溶剂分子呢?还是不用去除

82

帖子

2

威望

793

eV
积分
915

Level 4 (黑子)

分子模拟晶戈

13#
发表于 Post on 2022-10-26 11:29:21 | 只看该作者 Only view this author
xxzj 发表于 2022-10-25 18:21
老师,还需要再请教您三个问题,辛苦您
(1)水分子的参数文件是否也可以通过sobtop进行获取呢?
(2) ...

1.不用 用amber力场目录下自带的tip3p就行
2.一样,甚至都不需要packmol,直接用insert-molecules即可,自组装是一个自驱动的过程,尽量让初始的体系没有人为的“组装”倾向。
3.没必要,甚至有副作用。

349

帖子

0

威望

905

eV
积分
1258

Level 4 (黑子)

14#
发表于 Post on 2022-10-26 11:33:51 | 只看该作者 Only view this author
对抗路达摩 发表于 2022-10-26 11:29
1.不用 用amber力场目录下自带的tip3p就行
2.一样,甚至都不需要packmol,直接用insert-molecules即可, ...

好嘞,谢谢您,辛苦啦

本版积分规则 Credits rule

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

GMT+8, 2026-2-20 05:08 , Processed in 0.249967 second(s), 20 queries , Gzip On.

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