计算化学公社

 找回密码 Forget password
 注册 Register
Views: 77616|回复 Reply: 59

[GROMACS] 使用Gromacs模拟碳纳米管的一个简单例子

  [复制链接 Copy URL]

6万

帖子

99

威望

6万

eV
积分
130245

管理员

公社社长

发表于 Post on 2014-12-14 22:54:03 | 显示全部楼层 Show all |阅读模式 Reading model
注:本文的做法如今已经完全过时,如今强烈建议用Sobtop(http://sobereva.com/soft/Sobtop)建立纳米球/管/板的拓扑文件,比x2top灵活得多,参考网页里的例子。另外,北京科音分子动力学与GROMACS培训班(http://www.keinsci.com/workshop/KGMX_content.html)里有模拟纳米管、纳米板的完整讲解和详细的例子。


使用Gromacs模拟碳纳米管的一个简单例子
A simple example of using Gromacs to simulate carbon nanotubes

文/Sobereva @北京科音   2014-Dec-15


很早以前写过一个帖子《amber与gromacs读入碳纳米管的方法》(http://sobereva.com/39)介绍了怎么用gromacs处理碳纳米管,不过貌似还是有很多人不会。这里就提供一个在真空中和在水中模拟碳纳米管的简单且完整的例子。这个小教程本来是给一个国际友人写的,所以都是英文。本文用的是Gromacs 4.6.5,文中的CNT10指的是一个很普通的碳纳米管,cnt10.pdb是其结构文件,这可以用Nanotube Modeler很容易地产生。

文中涉及到的文件在此:
md.mdp (898 Bytes, 下载次数 Times of downloads: 417)

评分 Rate

参与人数
Participants 5
eV +23 收起 理由
Reason
一条君 + 5 好物!
heqiandai + 4 谢谢
Cladoss + 4
hlmkh + 5 <font><font>با &
ter20 + 5 谢谢

查看全部评分 View all ratings

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

294

帖子

0

威望

3541

eV
积分
3835

Level 5 (御坂)

发表于 Post on 2014-12-15 07:33:34 | 显示全部楼层 Show all
强烈顶贴!Sob老师就是强!

4

帖子

0

威望

66

eV
积分
70

Level 2 能力者

发表于 Post on 2015-3-5 09:18:01 | 显示全部楼层 Show all
Sob老师,问一个比较傻的问题:这里如果用
g_x2top -f cnt10.pdb -o CNT10.top -ff select -nopbc -name CNT10 -kb 400000 -kt 600 -kp 150
atomname2type.n2t可以让距离0.142nm的碳看成是CR1的碳, 但是怎么得到键,角的信息呢?
因为我在用的时候老是显示: can not find forcefield for atom C with 0 bonds
非常感谢!!!

6万

帖子

99

威望

6万

eV
积分
130245

管理员

公社社长

 楼主 Author| 发表于 Post on 2015-3-5 18:50:11 | 显示全部楼层 Show all
shanshan 发表于 2015-3-5 09:18
Sob老师,问一个比较傻的问题:这里如果用
g_x2top -f cnt10.pdb -o CNT10.top -ff select -nopbc -name C ...

你先严格、完整地重复一下帖子中的例子,需要的文件都提供了。如果不行,换成文中的gmx版本(有的版本在判断时可能有bug)
只要能根据n2t正确判断出原子,在top里bond、angle、dihedral项就都会有了
北京科音自然科学研究中心(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

4

帖子

0

威望

66

eV
积分
70

Level 2 能力者

发表于 Post on 2015-3-6 03:24:41 | 显示全部楼层 Show all
sobereva 发表于 2015-3-5 18:50
你先严格、完整地重复一下帖子中的例子,需要的文件都提供了。如果不行,换成文中的gmx版本(有的版本在 ...

非常感谢~~~试了无数遍终于发现是版本问题。

7

帖子

0

威望

32

eV
积分
39

Level 2 能力者

发表于 Post on 2015-4-9 10:04:38 | 显示全部楼层 Show all
你好,我想请问一下,如果初始的pdb是一个蛋白质和纳米管的复合物,然后再加溶剂进行动力学,那么生成top的过程是怎么样的?

6万

帖子

99

威望

6万

eV
积分
130245

管理员

公社社长

 楼主 Author| 发表于 Post on 2015-4-9 13:57:10 | 显示全部楼层 Show all
lingogo 发表于 2015-4-9 10:04
你好,我想请问一下,如果初始的pdb是一个蛋白质和纳米管的复合物,然后再加溶剂进行动力学,那么生成top的 ...

单独构造一个纳米管的itp文件。然后用一般方法构建蛋白质的top文件,把纳米管的include进去并把纳米管的结构和蛋白质结构合并。然后加水。
北京科音自然科学研究中心(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

82

帖子

0

威望

7313

eV
积分
7395

Level 6 (一方通行)

发表于 Post on 2015-4-9 18:00:16 | 显示全部楼层 Show all
谢谢分享!

7

帖子

0

威望

32

eV
积分
39

Level 2 能力者

发表于 Post on 2015-4-14 10:34:42 | 显示全部楼层 Show all
sobereva 发表于 2015-4-9 13:57
单独构造一个纳米管的itp文件。然后用一般方法构建蛋白质的top文件,把纳米管的include进去并把纳米管的 ...

sob老师,再请问一下,碳纳米管的itp文件如何构造呢?

6万

帖子

99

威望

6万

eV
积分
130245

管理员

公社社长

 楼主 Author| 发表于 Post on 2015-4-18 05:16:28 | 显示全部楼层 Show all
lingogo 发表于 2015-4-14 10:34
sob老师,再请问一下,碳纳米管的itp文件如何构造呢?

按此贴的方法用x2top就得到了
北京科音自然科学研究中心(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

7

帖子

0

威望

32

eV
积分
39

Level 2 能力者

发表于 Post on 2015-4-18 10:14:11 | 显示全部楼层 Show all
sobereva 发表于 2015-4-18 05:16
按此贴的方法用x2top就得到了

sob老师,不好意思我比较笨,我用了这个方法,最后只得到了CNT10.top,itp文件需要用什么命令生成或者转化吗?

294

帖子

0

威望

3541

eV
积分
3835

Level 5 (御坂)

发表于 Post on 2015-4-18 14:25:30 | 显示全部楼层 Show all
你把CNT10.top改名为CNT10.itp,然后自己写一个top文件将这个itp文件include进去就可以了!建议看看手册第5章的内容

评分 Rate

参与人数
Participants 1
eV +1 收起 理由
Reason
sobereva + 1

查看全部评分 View all ratings

7

帖子

0

威望

32

eV
积分
39

Level 2 能力者

发表于 Post on 2015-4-19 08:55:23 | 显示全部楼层 Show all
ruanyang 发表于 2015-4-18 14:25
你把CNT10.top改名为CNT10.itp,然后自己写一个top文件将这个itp文件include进去就可以了!建议看看手册第5 ...

非常感谢

6

帖子

0

威望

98

eV
积分
104

Level 2 能力者

发表于 Post on 2017-10-29 20:35:56 | 显示全部楼层 Show all
看来这个软件用途广,学习

15

帖子

0

威望

111

eV
积分
126

Level 2 能力者

发表于 Post on 2019-6-14 23:16:32 | 显示全部楼层 Show all
本帖最后由 wh.shen 于 2019-6-14 23:38 编辑

sob老师,按照您的方法在gromos54a7.ff力场下可以做碳纳米管的模拟。我要做的是在amber03力场下丝蛋白加入纳米管,这个时候我在主top文件中include CNT.itp(CNT.itp是在gromos54a7.ff力场下生成的),但是在生成运行文件时出现截图中错误。
当我尝试在amber03力场下生成CNT.itp时,(同样将CNT.itpinclude到主top中)错误变为找不到CR1类型原子。
请教老师!@sobereva
2019-06-14 23-11-27屏幕截图.png

本版积分规则 Credits rule

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

GMT+8, 2026-10-15 21:41 , Processed in 0.809169 second(s), 32 queries , Gzip On.

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