求助：lammps控压问题，导致主体相的水密度一直达不到1

查看数: 1744 | 评论数: 8 | 收藏 Add to favorites 0

关灯 | 提示：支持键盘翻页<-左右->

zcw

发布时间: 2023-7-14 16:31

求助：各位老师，使用lammps进行固体板在水溶液中的NPT模拟时，发现计算的水密度在主体相部分一直不为1，从log文件也发现最开始压强也是突然下降导致盒子变大。想要知道我控压的方式是否有问题，有没有更好的控压方 ...

回复 Reply

zcw 发表于 Post on 2025-6-12 09:33:51

liangxf 发表于 2025-6-8 10:11
请问后续解决了吗？我在计算类似问题的时候，在最小化是周期性边界条件，控压是沿z方向，计算结果总是出现 ...

我后面更换使用的gromacs软件，计算出来的密度在0.997左右。相同的体系，相同的控温孔压方式。

liangxf 发表于 Post on 2025-6-8 10:11:44

请问后续解决了吗？我在计算类似问题的时候，在最小化是周期性边界条件，控压是沿z方向，计算结果总是出现盒子沿z方向扩张导致质量密度总是达不到1

Lacrimosa 发表于 Post on 2023-7-15 13:15:22

zcw 发表于 2023-7-14 23:32
你好，我看了文章SPC水的密度大概计算在0.975g/cm3，可我计算的只有0.93g/cm3左右。因此想知道是不是控压 ...

建个纯水体系以同样的控压参数跑个模拟算下密度就知道了，如果纯水体系密度和文献一致，界面体系不一致，那就做个水分子沿z方向的密度分布，看看是不是界面造成的影响。

zcw 发表于 Post on 2023-7-14 23:32:14

牧生发表于 2023-7-14 21:15
4℃的时候，水的密度最大，为1 g/cm3

其次，看看这里 http://bbs.keinsci.com/thread-25479-1-1.html

你好，我看了文章SPC水的密度大概计算在0.975g/cm3，可我计算的只有0.93g/cm3左右。因此想知道是不是控压有问题，有没有更佳的控压方式。

zcw 发表于 Post on 2023-7-14 22:01:23

牧生发表于 2023-7-14 21:15
4℃的时候，水的密度最大，为1 g/cm3

其次，看看这里 http://bbs.keinsci.com/thread-25479-1-1.html

好的，感谢您的答复，我去拜读一下文章

牧生 发表于 Post on 2023-7-14 21:15:50

4℃的时候，水的密度最大，为1 g/cm3

其次，看看这里 http://bbs.keinsci.com/thread-25479-1-1.html

Lacrimosa 发表于 Post on 2023-7-14 18:47:29

信息不足不太好判断，你用的单位制是什么，水模型是什么，以及实际计算出的密度是多少，这些信息得提供一下。另外，最好放一张你体系的结构图

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register