计算化学公社

标题: THF盒子平衡相模拟问题 [打印本页]

作者
Author: deltazheng 时间: 2018-12-2 12:56
标题: THF盒子平衡相模拟问题
我尝试用gromacs模拟一个纯THF盒子，在做平衡相模拟时遇到一个问题。如果使用eq.mdp文件可以正常进行前处理和模拟，但是如果将eq.mdp中compressibility从4.6e-5改为9.7e-5（THF的等温压缩系数）则会出现如下错误提示：

$ gmx grompp -f eq.mdp -c em.gro -p THFbox.top -o eq.tpr
...
$ gmx mdrun -v -deffnm eq
...
Step 11 Warning: pressure scaling more than 1%, mu: 0.987387 0.987387 0.987387

Step 21 Warning: pressure scaling more than 1%, mu: 0.986063 0.986063 0.986063
...

请问应该如何解决？

作者
Author: sobereva 时间: 2018-12-2 13:17
其实改这个意义不大，因为各种液体的可压缩系数差的并不多，而且这个参数主要影响控压过程的盒子波动情况，对最终感兴趣的问题影响不大。
虽然有warning，但最终模拟的结果没问题，那也无妨。可压缩系数大了，控压造成的盒子体积波动就会变大，就更容易出现这个提示。你的任务前期是升温，等温度稳定后，估计这种warning也就基本没了。

作者
Author: deltazheng 时间: 2018-12-2 13:25

sobereva 发表于 2018-12-2 13:17
其实改这个意义不大，因为各种液体的可压缩系数差的并不多，而且这个参数主要影响控压过程的盒子波动情况， ...

感谢Sob老师的解答~

作者
Author: deltazheng 时间: 2018-12-2 13:39

sobereva 发表于 2018-12-2 13:17
其实改这个意义不大，因为各种液体的可压缩系数差的并不多，而且这个参数主要影响控压过程的盒子波动情况， ...

的确，平衡相模拟9000步以后就没有这种提示了。

作者
Author: 牧生 时间: 2021-7-28 21:16
近几日也遇到了这个问题，改成如下参数，就好了。

tau_t = 2

tau_p = 4

欢迎光临计算化学公社 (http://bbs.keinsci.com/)