zgdong 发表于 2023-7-10 00:18 您好，请问您解决了吗？我现在也想做这个 |
| 您好，有相关WHAM的资料吗，能否分享一份。我想求解气体分子从气相到液相过程的自由能变化，不知道怎么入手，谢谢啦 |
k64_cc 发表于 2019-3-1 13:59 嗯 或许我应该试一下Gromacs |
k64_cc 发表于 2019-3-1 13:59 好的好的 我再试试，谢谢您！ |
wasngsimin 发表于 2019-3-1 13:48 我没用过Lammps，不知道……最简单的就是把所有加bias的操作都转移到PLUMED上，PLUMED操作方便，接口也好用。 fix com那个也是，Lammps很大概率能做，但我真的是不会lammps…… 以及GMX比Lammps计算效率高啊，为啥不用？ |
本帖最后由 wasngsimin 于 2019-3-1 13:56 编辑 k64_cc 发表于 2019-2-28 15:54 yes, you are right. You are a kind and smart person who capture the key point. I do want to calculate the absorption interaction between the rightside soild and the small methane molecule in the water solution. I have run a simulation in Lammps using fix recenter to fix the center of the rightside solid this Tuesday. I defined the distance beteween the center of the solid and the center of the small methane molecule as the collective coodinate(CV). The code is still ruunning so I donnot known whether the outcomes is reasonable now. I looked through the Lammps manual carefully but to be honestly I did not found the command fix COM you refered here. And I want to know how to count CV? Could you please elabrate on it? Maybe it is a stupid question.Thanks for your suggestions! |
|
本帖最后由 k64_cc 于 2019-2-28 15:58 编辑 WHAM和MetaD和ABF都差不多，想拿放心数据的话需要的data数量也没啥太大差别——eABF采样做的比较好，MetaD系列就别想太多了，Umbrella Sampling有的毛病它都有。 |
|
本帖最后由 k64_cc 于 2019-2-28 15:58 编辑 你这哪行啊…… 1.你CV选的没物理意义。我猜你是想改变小分子和固液界面的距离，那CV应该是固体质心和小分子的x distance。体系整体移动的话fix COM可破，涉及到你thermostat不同，fix COM对体系的影响也不一样。建议用Langevin. 2.多拿轨迹就平衡了…具体就是做到即使加数据error bar也不减小。 3.Count CV, draw histogram. 以及PMF……全称是potential of mean force, 是mean force的积分——平均力的积分，你体会一下…… |
fhh2626 发表于 2019-2-25 17:18 谢谢您的回复！我是算整个体系对小分子的作用 |
wasngsimin 发表于 2019-2-25 12:04 我可能是理解错了，你不是算小分子的作用，而是算水中两坨分子的作用？ 总之约束大的那一坨的质心，让小分子动就行啦，不要让反应坐标的组成部分跨越pbc 都是理论严谨的方法（MtD还没有数学上证明过），没有精度一说，只能说收敛速度的快慢 PMF是自由能随反应坐标的变化，是热力学量 |
fhh2626 发表于 2019-2-25 11:19 另外我还有个疑问：PMF应该是瞬时的性质吧？体现了把分子放到水溶液里它会往哪边移动，那MD模拟那么长的时间及采样过程，构型会发生变化，这样说来计算的PMF会随时间变化吧？即后面计算的PMF还是初始构型的吗？ |
|
0、你用的应该是distanceZ作为反应坐标吧 1、建议约束水的质心，并且把PBC的盒子放大一些，不要让水跨过PBC边界，要不然colvars计算质心会出错 2和3、不建议用伞状采样，在colvars里面直接用eABF或者WT-MtD简单多了 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2025-8-17 03:03 , Processed in 0.173151 second(s), 25 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites