“第10届量子化学波函数分析与Multiwfn程序培训班将于5月4-8日于北京举办,这是一次性完整、系统学习波函数分析的各种理论知识和全面掌握强大的Multiwfn波函数分析程序使用的最不可错过的机会!请点击此链接查看详情和报名方式,欢迎参加!

“第18届北京科音分子动力学与GROMACS培训班” 将于5月23-26日于北京举办。这是一次性全面、系统学习分子动力学模拟知识和最流行的分子动力学程序GROMACS的关键机会!报名正在进行中,请点击此链接查看详情,欢迎参加!

计算化学公社

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

[GROMACS] 求助:用Gromacs做多肽穿过POPC膜时,出现伞形势能不断增加问题该怎么解决

[复制链接 Copy URL]

4

帖子

0

威望

93

eV
积分
97

Level 2 能力者

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 coisiniiiiiii 于 2026-3-28 08:43 编辑

请教各位老师。我在计算多肽穿膜过程的所需要的能量来衡量穿膜的难易程度。先进行SMD拉伸,拉伸轨迹中提取出120个文件进行采样。发现得到的伞形势是不断增加的,与图1的伞形势能不相符。
想请问各位老师,这个问题应该如何解决?感谢各位老师!



image-2.png (122.75 KB, 下载次数 Times of downloads: 1)

图1:伞形势能

图1:伞形势能

202603280828329165..png (15.67 KB, 下载次数 Times of downloads: 1)

图2:伞形势能

图2:伞形势能

202603280827485460..png (23.62 KB, 下载次数 Times of downloads: 1)

图3:采样窗口

图3:采样窗口

输入文件.zip

1.83 MB, 下载次数 Times of downloads: 0

输入文件

pull_2.bat

1.19 KB, 下载次数 Times of downloads: 1

计算用到的脚本文件

331

帖子

0

威望

2072

eV
积分
2403

Level 5 (御坂)

2#
发表于 Post on 2026-3-28 09:16:09 | 只看该作者 Only view this author
应该先检查一下smd的轨迹是否合理

4

帖子

0

威望

93

eV
积分
97

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2026-3-28 10:19:10 | 只看该作者 Only view this author
pal 发表于 2026-3-28 09:16
应该先检查一下smd的轨迹是否合理

您好,SMD轨迹如下。感谢

SMD.mp4

5.28 MB, 下载次数 Times of downloads: 4

107

帖子

2

威望

1878

eV
积分
2025

Level 5 (御坂)

4#
发表于 Post on 2026-3-28 11:58:13 | 只看该作者 Only view this author
1. 我这图牵引的乙醇分子, 这个分子是非常小的, 并且本身就可以被动扩散进去细胞膜. 他的跨膜PMF曲线当然和你的多肽的跨膜PMF曲线不相似.
2. 针对多肽这种的分子, 你应该和其他做多肽跨膜的论文进行比较.
3. 强烈建议用colvars做增强采样, 然后用WHAM程序计算PMF!
分子模拟玩家

4

帖子

0

威望

93

eV
积分
97

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2026-3-28 12:36:28 | 只看该作者 Only view this author
Acee 发表于 2026-3-28 11:58
1. 我这图牵引的乙醇分子, 这个分子是非常小的, 并且本身就可以被动扩散进去细胞膜. 他的跨膜PMF曲线当然和 ...

感谢解答,我现在去试试

本版积分规则 Credits rule

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

GMT+8, 2026-4-13 23:32 , Processed in 0.174754 second(s), 23 queries , Gzip On.

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