计算化学公社

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

[GROMACS] 油水两相界面,拉伸分子动力学模拟和采样相关问题

[复制链接 Copy URL]

97

帖子

0

威望

634

eV
积分
731

Level 4 (黑子)

模拟体系为油水两相界面,先经过NPT之后使得油相中活性分子吸附在界面,然后把界面吸附的一个活性分子拉下来,参数设置如下:

拉伸完后通过可视化发现被拉伸的分子几乎没有移动,反复在界面附近波动,但是通过对该分子位置坐标随着时间的变化分析发现分子确实沿着z轴发生了移动,请问能否说明该分子已经被拉动?还有其他方法可以验证吗?

另外,拉伸结束后进行了伞状采样,采样思路是对拉伸过程生成一系列构象然后对每个构象计算拉伸组和参考组之间的距离,然后根据这一组距离进行筛选,每隔一段距离采样一次,最后得到的采样直方图和能量如下(绘制成一幅图了,能量是图里的红色粗线,对应右边纵坐标):


请问这个采样思路是否合理?有些部分并没有采样,但是有些地方又被多次采样为什么?有什么推荐的采样方法吗?



395

帖子

8

威望

3912

eV
积分
4467

Level 6 (一方通行)

石墨

2#
发表于 Post on 2025-3-21 17:38:18 | 只看该作者 Only view this author
是不是油相和小分子一起向z方向移动了,z绝对坐标完美移动,小分子和油相绝对坐标只是小范围波动。
也就是周期性和相对位置的问题,例如是三轴都是周期性,跑的时候整个体系的绝对坐标如图2下一样都在向z方向移动,只不过因为vmd等软件显示的时候以第0帧为基准,显示的是相对运动,如图2上的情况。
先解决这个问题,然后第二步采样的时候每个采样窗口调整k,令它在每个窗口都固定好,做到比较均匀的采样。
自在飞花轻似梦,无边丝雨细如愁。

全自动反应动力学(ReaxFF、AIMD、NEP等)后处理工具网页版:http://cc-portal.xyz/reax_tools

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

GMT+8, 2025-8-16 06:45 , Processed in 0.152146 second(s), 24 queries , Gzip On.

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