在lammps中如何使用gcmc模拟获取气体密度

hhhhhhhhh · 发表于 Post on 2024-12-29 12:45:25

本帖最后由 hhhhhhhhh 于 2024-12-29 12:45 编辑

各位老师好，我按照lammps中自带的gcmc教程案例尝试模拟333.15k，20MPa下CO2的密度，但是始终无法得到与NIST数据库相匹配的密度值，想请各位老师帮助我看看是什么原因？
备注：1. 333.15K，20MPa下，CO2的密度是0.72368g/cm3，逸度系数是0.45718（数据来源refprop）
2. 图是文献中压力-化学势曲线，按照图上的信息可以得到此条件下的化学势是7.88kJ/mol（1.8834kcal/mol）
3. gcmc的原教程文件是mc/in.gcmc.co2.lmp文件，修改了盒子大小和晶格大小，以及化学势和fix gcmc命令
4. fix gcmc尝试两种控制a.设置化学势，b.设置压强以及逸度系数（此时化学势设置失效），两者均无法得到想要的密度。上传的附件是设置压强和逸度系数的，使用化学势时，注释了这两项参数
5. 文献使用gcmc+md的过程，我认为md主要是用于弛豫体系，且目前我对使用gcmc过程的正确性无法保证，因此省略了md的过程，不知道md是否是导致密度无法一直的主要原因

附件包含输入文件，CO2的分子模型，以及两种不同设置的fix gcmc模拟的log文件

gcmc.in (3 KB, 下载次数 Times of downloads: 27)

CO2.txt (436 Bytes, 下载次数 Times of downloads: 11)

pressure-log.lammps (177.48 KB, 下载次数 Times of downloads: 8)

mu-log.lammps (75.34 KB, 下载次数 Times of downloads: 8)

杨星元 · 发表于 Post on 2025-3-6 20:49:06

兄弟，我用lammps算145.66atm,304.2K的二氧化碳的密度与NIST也是差了0.17g/cm^3左右，我用的是NPT系综弛豫的，EPM2模型，你知道这是为啥吗

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

[Lammps] 在lammps中如何使用gcmc模拟获取气体密度