bobchi999 发表于 2026-1-29 11:41 抱歉，没用过MS，我说的是Gromacs |
bobchi999 发表于 2026-1-29 11:49 我没用过ms干过类似的事情，无法解答 |
牧生 发表于 2026-1-29 11:20 您好，看来MS里只能冻结+NVT了是吧？ 那由于cell尺寸不能变，分子团聚会形成真空层，这大概就是我遇到的情况。 那有什么方法，在MS里可以让水分子的运动更合理呢？ |
reid 发表于 2026-1-28 21:42 您好， 我刚试了一下，如果按您说的给所有Fe原子添加限制，在MS里就需要测量每一个Fe原子和相邻三个Fe原子的距离，然后添加限制势。这个操作太繁琐了，只能写Script来实现。 因此我想问一下，您在MS里做过类似操作吗？ |
|
实际经验，铁不能用限制，因为仍然会变形。还是得用冻结才能保持不变。 从晶体制造者里面弄出来的铁胞体非常好用，用了冻结，在gmx跑起来，且结合npt毫无问题。 从m$里面弄出来的铁胞，如果放在gmx里面用，无论如何nvt，npt，冻结，也都有问题。 |
reid 发表于 2026-1-28 21:42 好，我试一下 |
sobereva 发表于 2026-1-28 15:29 谢谢老师！ 那我这个cell里水分子的运动情况是什么原因导致的呢？ |
| 给每个Fe原子一个1500~2000的限制势，3个方向都加上，然后用NPT |
|
谈谈做分子动力学模拟什么时候适合用NPT和NVT系综 http://sobereva.com/761（http://bbs.keinsci.com/thread-57835-1-1.html） 上文专门说了 非要用冻结就用NVT，要用NPT就把冻结改成位置限制势 |
bobchi999 发表于 2026-1-28 08:52 有时候冻结会使NPT不稳定，因为NPT涉及到坐标的放缩，是否冻结看研究目的吧 我不用MS，对于MS中模型的构建没法帮到你 |
本帖最后由 bobchi999 于 2026-1-28 08:55 编辑 KazusaT 发表于 2026-1-27 22:18 我并没有主动地在水和分子间加真空层，只是按照通常做法把分子和水分别设成一个layer添加进cell，这应该是MS里很常见的操作吧？ 如果真的有真空层存在，那肯定是我设置cell时搞出了问题，现在求各路大神帮忙出出主意。 之前逛论坛，卢老师在回复一个帖子时好像提到过，冻结原子后最好不要用NPT。也可能是我记错了。 |
|
为什么初始设置要在水和分子之间加真空层而不是直接用水分子填满整个盒子？ 冻结铁原子做NPT没有什么问题，如果Fe只是充当一个基底，结构变化不重要的话，冻结+半各向同性控压的话XY方向应该是不动的（其实对于这种周期性基底来说XY方向本来也不会变，冻结只是避免在Z方向上的移动吧） |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-2-16 20:23 , Processed in 0.181278 second(s), 25 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites