计算化学公社

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

[GROMACS] 求助:模拟界面吸附过程跑NVT分子吸附不下来

[复制链接 Copy URL]

8

帖子

0

威望

637

eV
积分
645

Level 4 (黑子)

本帖最后由 Black-Mamba2018 于 2019-5-4 16:59 编辑

各位老师好!本人在模拟界面吸附过程,跑NVT(mdp文件如图所示)500ps、2ns分子始终吸附不下来(如图所示),只有极少数分子下去了,甚至穿到底部去了,目标分子A和B距离底层分子的高度为2nm,且底层分子全都固定了,请问各位老师,怎样才能使目标分子A和B全部吸附上去呢?

mdp.png (22.69 KB, 下载次数 Times of downloads: 27)

mdp文件

mdp文件

md2.png (152.21 KB, 下载次数 Times of downloads: 37)

初始结构

初始结构

md2-500ps.png (98.36 KB, 下载次数 Times of downloads: 36)

跑完500ps后的结构

跑完500ps后的结构

md2-2ns.png (99.94 KB, 下载次数 Times of downloads: 38)

跑完2ns后的结构

跑完2ns后的结构

6万

帖子

99

威望

5万

eV
积分
120206

管理员

公社社长

2#
发表于 Post on 2019-5-4 17:08:52 | 只看该作者 Only view this author
rvdw再设大点,初始结构间距摆小点,控温温度设低点再试
北京科音自然科学研究中心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

8

帖子

0

威望

637

eV
积分
645

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2019-5-12 09:58:55 | 只看该作者 Only view this author
sobereva 发表于 2019-5-4 17:08
rvdw再设大点,初始结构间距摆小点,控温温度设低点再试

感谢sob老师的指点,我按照您给的建议,继续尝试,但跑完仍有分子跑到底部和顶部去,我在VMD中打开Z方向的镜像看了下,结果发现有成键,请问老师这该怎么处理呢?是添加periodic-molecules=yes(已在尝试) ?对于这句话"One host or emitter molecule was deposited every 100 ps within a 2 nm height from the top of the relaxed substrate or previously deposited
molecules."分子相当于以0.02 nm/ps沉积到基底去,这是怎么做到的呢?我参照咱们培训班的教材中的一个例子(金原子撞击金板),给所有要沉积分子(原子)添加速度矢量(0.0000  0.0000 -0.0200),跑完还是没有全部沉积到基底表面,学生学艺不精,还请老师指点一二。

6万

帖子

99

威望

5万

eV
积分
120206

管理员

公社社长

4#
发表于 Post on 2019-5-12 10:19:14 | 只看该作者 Only view this author
你的图我看不清楚
如果跑到底下的是跑到了下方分子的中间去,说明下方分子排列不够致密
如果是跑到了下方分子低端,那应当是由于PBC条件所致(具体通过轨迹考察)

你的体系里没有两头都跨越盒子的分子,用periodic-molecules=yes毫无意义。培训班例子里那种诸如无限延展的石墨烯才需要用这个

gro里设速度就行。但是如果你用了控温,实际速度就不会和你预设的一致。
北京科音自然科学研究中心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

8

帖子

0

威望

637

eV
积分
645

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2019-5-12 11:02:54 | 只看该作者 Only view this author
sobereva 发表于 2019-5-12 10:19
你的图我看不清楚
如果跑到底下的是跑到了下方分子的中间去,说明下方分子排列不够致密
如果是跑到了下方 ...

感谢老师的悉心指点

本版积分规则 Credits rule

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

GMT+8, 2025-8-18 10:07 , Processed in 1.187882 second(s), 23 queries , Gzip On.

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