请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 5967|回复: 13

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

[复制链接]

1万

帖子

25

威望

1万

eV
积分
34206

管理员

公社社长

发表于 2014-12-14 22:54:03 | 显示全部楼层 |阅读模式
使用Gromacs模拟碳纳米管的一个简单例子

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


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

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

评分

参与人数 2eV +10 收起 理由
hlmkh + 5 <font><font>با &
ter20 + 5 谢谢

查看全部评分

北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

290

帖子

0

威望

2860

eV
积分
3150

Level 5 (御坂)

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

4

帖子

0

威望

66

eV
积分
70

Level 2 能力者

发表于 2015-3-5 09:18:01 | 显示全部楼层
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
非常感谢!!!

1万

帖子

25

威望

1万

eV
积分
34206

管理员

公社社长

 楼主| 发表于 2015-3-5 18:50:11 | 显示全部楼层
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程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

4

帖子

0

威望

66

eV
积分
70

Level 2 能力者

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

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

7

帖子

0

威望

30

eV
积分
37

Level 2 能力者

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

1万

帖子

25

威望

1万

eV
积分
34206

管理员

公社社长

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

单独构造一个纳米管的itp文件。然后用一般方法构建蛋白质的top文件,把纳米管的include进去并把纳米管的结构和蛋白质结构合并。然后加水。
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

32

帖子

0

威望

2087

eV
积分
2119

Level 5 (御坂)

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

7

帖子

0

威望

30

eV
积分
37

Level 2 能力者

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

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

1万

帖子

25

威望

1万

eV
积分
34206

管理员

公社社长

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

按此贴的方法用x2top就得到了
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

7

帖子

0

威望

30

eV
积分
37

Level 2 能力者

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

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

290

帖子

0

威望

2860

eV
积分
3150

Level 5 (御坂)

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

评分

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

查看全部评分

7

帖子

0

威望

30

eV
积分
37

Level 2 能力者

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

非常感谢

6

帖子

0

威望

47

eV
积分
53

Level 2 能力者

发表于 2017-10-29 20:35:56 | 显示全部楼层
看来这个软件用途广,学习
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2018-9-26 13:44 , Processed in 0.219785 second(s), 27 queries .

快速回复 返回顶部 返回列表