计算化学公社

 找回密码 Forget password
 注册 Register

新手求助:真空环境下分子动力学,npt下退火出现错误问题

查看数: 427 | 评论数: 14 | 收藏 Add to favorites 0
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2025-3-24 16:19

正文摘要:

各位老师好,我在进行真空环境下分子动力学模拟时,在npt情况下进行退火,跑完gmx mdrun后的log文件出现下面的错误,这种该怎么解决,新手小白希望各位大佬能给予指导,万分感谢

回复 Reply

科研小白0126 发表于 Post on 2025-3-25 17:37:39
snljty2 发表于 2025-3-25 17:09
尽量接近你想模拟的实验条件。(以下只是举个例子,不代表就是最合理的做法)如果比如是溶液加工得到的, ...

好的明白了谢谢您
科研小白0126 发表于 Post on 2025-3-25 17:37:18
CrysW555 发表于 2025-3-25 16:25
我觉得你可以试试用packmol的时候把盒子尺寸弄大一点,让初始结构松散一点,然后不管是NVT还是NPT,分子会 ...

好的谢谢
snljty2 发表于 Post on 2025-3-25 17:09:50
科研小白0126 发表于 2025-3-25 15:52
您好我是新手,不是很熟悉,那我这种具体应该怎么操作

尽量接近你想模拟的实验条件。(以下只是举个例子,不代表就是最合理的做法)如果比如是溶液加工得到的,可以考虑加一堆溶剂分子,跑单分子的NPT得到一段轨迹,必要的时候可以用周期性退火,找一系列轨迹中势能低点处的结构,当做考虑了无序的初始结构,再依次插入一个足够大的盒子里,用溶剂填充,MD跑平衡后再用各种方式跑溶剂逐渐消失的非平衡动力学来得到一种比较符合实际情况的构象。这类分子尤其如果烷基链很长,直接摆一堆一样的上来就跑NPT,肯定和实验条件得到的形貌差距非常大的。
CrysW555 发表于 Post on 2025-3-25 16:25:46
我觉得你可以试试用packmol的时候把盒子尺寸弄大一点,让初始结构松散一点,然后不管是NVT还是NPT,分子会自己聚集,然后你再退火,我觉得NVT会好一点。
科研小白0126 发表于 Post on 2025-3-25 16:17:11
海棠依旧cdut 发表于 2025-3-25 16:11
其实别人已经回答你了,如果不组内之前的"文章"这样做的话,建议你不要这样做。目的是告诉你应该从类似的 ...

好的
海棠依旧cdut 发表于 Post on 2025-3-25 16:11:16
科研小白0126 发表于 2025-3-25 15:52
您好我是新手,不是很熟悉,那我这种具体应该怎么操作

其实别人已经回答你了,如果不组内之前的"文章"这样做的话,建议你不要这样做。目的是告诉你应该从类似的文章中去找被人的处理方案。以别人的模版进行你自己的创新
科研小白0126 发表于 Post on 2025-3-25 15:52:17
snljty2 发表于 2025-3-25 15:40
你是说,这个是用packmol直接把一个单分子,按照完全相同的构象,直接插了50个到盒子里?如果不是组里文 ...

您好我是新手,不是很熟悉,那我这种具体应该怎么操作
snljty2 发表于 Post on 2025-3-25 15:40:46
科研小白0126 发表于 2025-3-25 15:15
您好这是我的gro文件能帮我看看吗,谢谢您

你是说,这个是用packmol直接把一个单分子,按照完全相同的构象,直接插了50个到盒子里?如果不是组里文章之前就这么做的,我不推荐这么干。
科研小白0126 发表于 Post on 2025-3-25 15:15:30
FrancisCho 发表于 2025-3-25 08:51
试试去掉h-bonds的约束?

好的我去试试
科研小白0126 发表于 Post on 2025-3-25 15:15:05
snljty2 发表于 2025-3-25 12:04
如果你想模拟的是一个明显非致密的体系,当然不能用NPT系综,控压算法会把真空区压没了的。建议上传一个g ...

您好这是我的gro文件能帮我看看吗,谢谢您

7_box.gro

489.79 KB, 下载次数 Times of downloads: 3

用packmol制作盒子的gro

emminim.gro

489.69 KB, 下载次数 Times of downloads: 0

能量最小化的gro

snljty2 发表于 Post on 2025-3-25 12:04:26
科研小白0126 发表于 2025-3-25 09:42
我没加溶剂这算真空吗

如果你想模拟的是一个明显非致密的体系,当然不能用NPT系综,控压算法会把真空区压没了的。建议上传一个gro方便大家检查。
科研小白0126 发表于 Post on 2025-3-25 09:42:50

我没加溶剂这算真空吗
student0618 发表于 Post on 2025-3-25 09:09:23
真空用npt?
FrancisCho 发表于 Post on 2025-3-25 08:51:23
试试去掉h-bonds的约束?

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

GMT+8, 2025-8-15 00:36 , Processed in 0.163176 second(s), 25 queries , Gzip On.

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