|
|
各位老师，我计算锂离子穿过单层COF膜的PMF，反应坐标选取的是锂离子(示意图中的红色粒子)到虚拟点(示意图中的黄色粒子)之间的距离。在每一个窗口中锂离子距离虚拟点的距离通过shell脚本自动替换mdp文件中pull-coord1-init后的值实现，力常数为1000 kJ/mol/nm2，每个窗口的间距是0.1 nm。整个模拟过程分为以下两个步骤：首先是平衡相(对应eq.mdp文件)过程，然后是产生相(对应prod.mdp文件)过程。PMF数据来自产生相过程，利用gmx中wham命令计算。计算的命令为gmx wham -it tpr.dat -if pullf.dat -bins 200 -bsres -nBootstrap 100 -unit kJ。按照我的理解，PMF曲线应该是左右对称的，但是实际计算出来的PMF曲线却不是这样的。不知道哪儿出了问题，麻烦各位老师帮我看看哪儿出了问题。以下是相关的文件。谢谢各位老师。
|
-
示意图.png
(215.25 KB, 下载次数 Times of downloads: 8)
示意图
-
histo-xvg.png
(24.19 KB, 下载次数 Times of downloads: 7)
histo.xvg图像
-
profile-xvg.png
(15.11 KB, 下载次数 Times of downloads: 7)
profile.xvg图像
-
-
eq.mdp
1.64 KB, 下载次数 Times of downloads: 18
eq过程的mdp
-
-
prod.mdp
1.45 KB, 下载次数 Times of downloads: 22
prod过程的mdp
|
发表于 Post on 2024-1-2 16:06:17
|
只看该作者 Only view this author
|只看大图 View big picture
|倒序浏览 Reverse view
|阅读模式 Reading model
收藏 Add to favorites
楼主 Author
，应该只在xy方向上施加位置限制(因为离子是沿着z轴穿越COF)。当我把z轴上的位置限制关掉后，离子在靠近COF附近PMF有一个峰，但是曲线的波动性比较大。