计算化学公社

标题: 在lammps中如何使用gcmc模拟获取气体密度 [打印本页]

作者
Author: hhhhhhhhh 时间: 2024-12-29 12:45
标题: 在lammps中如何使用gcmc模拟获取气体密度
本帖最后由 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是否是导致密度无法一直的主要原因
(, 下载次数 Times of downloads: 7)

附件包含输入文件，CO2的分子模型，以及两种不同设置的fix gcmc模拟的log文件
(, 下载次数 Times of downloads: 27)
(, 下载次数 Times of downloads: 11)

(, 下载次数 Times of downloads: 8)
(, 下载次数 Times of downloads: 8)

作者
Author: 杨星元 时间: 2025-3-6 20:49
兄弟，我用lammps算145.66atm,304.2K的二氧化碳的密度与NIST也是差了0.17g/cm^3左右，我用的是NPT系综弛豫的，EPM2模型，你知道这是为啥吗

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