计算化学公社

 找回密码 Forget password
 注册 Register
楼主 Author: sobereva
打印 Print 上一主题 Last thread 下一主题 Next thread

[综合交流] 分子动力学初始结构构建程序Packmol的使用

  [复制链接 Copy URL]

5万

帖子

99

威望

5万

eV
积分
112349

管理员

公社社长

16#
 楼主 Author| 发表于 Post on 2020-4-1 03:20:19 | 只看该作者 Only view this author
xzb 发表于 2020-3-31 11:31
请教老师,怎么才能实现在packmol中插入两个不同的分子,并控制它们之间的距离使它们靠近一些呢?

没有对分子间距离进行约束的设定,没法实现。如果只是希望出现在相近的空间区域倒是可以
北京科音自然科学研究中心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!

11

帖子

0

威望

45

eV
积分
56

Level 2 能力者

17#
发表于 Post on 2020-4-1 09:52:38 | 只看该作者 Only view this author
sobereva 发表于 2020-4-1 03:20
没有对分子间距离进行约束的设定,没法实现。如果只是希望出现在相近的空间区域倒是可以

好的谢谢老师

23

帖子

0

威望

111

eV
积分
134

Level 2 能力者

18#
发表于 Post on 2020-4-17 12:33:21 | 只看该作者 Only view this author
请问packmol建立石墨烯模型 好建立吗

5万

帖子

99

威望

5万

eV
积分
112349

管理员

公社社长

19#
 楼主 Author| 发表于 Post on 2020-4-19 00:40:09 | 只看该作者 Only view this author
lxc 发表于 2020-4-17 12:33
请问packmol建立石墨烯模型 好建立吗

packmol根本不是用来构建单个分子的
石墨烯用VMD自带的Extensions - Modeling - Nanotube builder就能建立
北京科音自然科学研究中心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!

7

帖子

0

威望

99

eV
积分
106

Level 2 能力者

20#
发表于 Post on 2020-6-27 12:55:20 | 只看该作者 Only view this author
本帖最后由 yinzi 于 2020-6-27 12:59 编辑

请问老师,packmol建立混合溶剂体系(多肽在水和甲醇)后,再用gromacs模拟时,npt出现问题,请问这是盒子设置问题出现了错误吗?
Fatal error:
The X-size of the box (5.297948) times the triclinic skew factor (1.000000) is
smaller than the number of DD cells (4) times the smallest allowed cell size
(1.325000)

5万

帖子

99

威望

5万

eV
积分
112349

管理员

公社社长

21#
 楼主 Author| 发表于 Post on 2020-6-27 18:00:15 | 只看该作者 Only view this author
yinzi 发表于 2020-6-27 12:55
请问老师,packmol建立混合溶剂体系(多肽在水和甲醇)后,再用gromacs模拟时,npt出现问题,请问这是盒子 ...

mdrun先结合-ntmpi 1跑试试
北京科音自然科学研究中心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!

63

帖子

0

威望

201

eV
积分
264

Level 3 能力者

22#
发表于 Post on 2020-7-8 16:39:59 | 只看该作者 Only view this author
新手求问packmo可以用于NAMD吗,利用packmol建模后怎么生成psf文件呢。

5万

帖子

99

威望

5万

eV
积分
112349

管理员

公社社长

23#
 楼主 Author| 发表于 Post on 2020-7-8 17:01:53 | 只看该作者 Only view this author
孙宝刚 发表于 2020-7-8 16:39
新手求问packmo可以用于NAMD吗,利用packmol建模后怎么生成psf文件呢。

我不用NAMD
北京科音自然科学研究中心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!

7

帖子

0

威望

99

eV
积分
106

Level 2 能力者

24#
发表于 Post on 2020-7-20 22:45:42 | 只看该作者 Only view this author
老师,我模拟多肽在乙腈和水的均匀混合体系中的md,但是在跑完1ns的npt之后,会发现原本均匀的溶剂,水分子会抱团。请问这是什么原因。inp文件如下:
#
# CsA inside water and methanol
#

tolerance 2.0
filetype pdb
output CsA-5-1086-1044-cubic-5-center.pdb

structure CsA-em.pdb
  number 1
  center
  fixed 25. 25. 25. 0. 0. 0.
end structure

structure jz2.pdb
  number 1086
  inside box 0. 0. 0. 50. 50. 50.
end structure

structure water.pdb
  number 1044
  inside box 0. 0. 0. 50. 50. 50.
end structure

5万

帖子

99

威望

5万

eV
积分
112349

管理员

公社社长

25#
 楼主 Author| 发表于 Post on 2020-7-21 00:31:07 | 只看该作者 Only view this author
yinzi 发表于 2020-7-20 22:45
老师,我模拟多肽在乙腈和水的均匀混合体系中的md,但是在跑完1ns的npt之后,会发现原本均匀的溶剂,水分子 ...

这只体现初始建模情况,而完全没有提供动力学过程的信息

乙腈和水的亲和性也不算特别高,出现局部相分离是正常的
北京科音自然科学研究中心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!

7

帖子

0

威望

99

eV
积分
106

Level 2 能力者

26#
发表于 Post on 2020-7-21 09:37:34 | 只看该作者 Only view this author
本帖最后由 yinzi 于 2020-7-22 15:06 编辑
sobereva 发表于 2020-7-21 00:31
这只体现初始建模情况,而完全没有提供动力学过程的信息

乙腈和水的亲和性也不算特别高,出现局部相分 ...

谢谢老师回复!动力学模拟的参数文件如下,在进行30ns的md模拟之后,发现水分子会渐渐成为柱状穿插在乙腈中,如图(图片不太能明确表示我所描述的状态)。模拟过程中的压力平均值在2,如下图所示,最终的主链骨架rmsd值如下图所示。这样的模拟结果是否值得信任。

npt.mdp

2.23 KB, 下载次数 Times of downloads: 11

5-md.mdp

2.34 KB, 下载次数 Times of downloads: 7

5万

帖子

99

威望

5万

eV
积分
112349

管理员

公社社长

27#
 楼主 Author| 发表于 Post on 2020-7-21 22:05:03 | 只看该作者 Only view this author
yinzi 发表于 2020-7-21 09:37
谢谢老师回复!动力学模拟的参数文件如下,在进行30ns的md模拟之后,发现水分子会渐渐成为柱状穿插在乙腈 ...

你应当贴考虑PBC的结构图,当前图什么也看不出来。
按照你的描述,就是形成了两相,这没什么异常
蛋白质结构合理与否一看结构便知。如果你要发表文章的话,MD最好再跑跑,RMSD还没有充分收敛
北京科音自然科学研究中心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!

7

帖子

0

威望

99

eV
积分
106

Level 2 能力者

28#
发表于 Post on 2020-7-22 15:06:17 | 只看该作者 Only view this author
sobereva 发表于 2020-7-21 22:05
你应当贴考虑PBC的结构图,当前图什么也看不出来。
按照你的描述,就是形成了两相,这没什么异常
蛋白 ...

谢谢老师的回答,我再磨一磨。

74

帖子

0

威望

247

eV
积分
321

Level 3 能力者

29#
发表于 Post on 2020-7-28 10:51:20 | 只看该作者 Only view this author
老师请问一下win版本和linux版本功能有差异吗

5万

帖子

99

威望

5万

eV
积分
112349

管理员

公社社长

30#
 楼主 Author| 发表于 Post on 2020-7-29 03:17:58 | 只看该作者 Only view this author
WB1040720479 发表于 2020-7-28 10:51
老师请问一下win版本和linux版本功能有差异吗

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

本版积分规则 Credits rule

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

GMT+8, 2024-11-23 00:49 , Processed in 0.185619 second(s), 23 queries , Gzip On.

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