miaof 发表于 2023-11-25 21:40
**** 作者被禁止或删除 内容自动屏蔽 ****

lammps自带的例子in.gcmc.lj中用variable             muex equal ${mu}-${temp}*ln(density*${lambda}+${nugget})直接计算出了过剩的化学势，想问一下，这里算出的过剩化学势和用fix widom计算的过剩化学势有联系吗？

dreamer123 发表于 2023-10-18 21:20
**** 作者被禁止或删除 内容自动屏蔽 ****

miaof 发表于 2023-11-25 21:40
Lammps中fix的widom命令可以算过剩的化学势在加上理想化学势就行了

我还有几个问题。
1、GCMC命令插入和删除粒子，并根据metropolis准则决定是否接受新构型，为什么还需要化学势呢？
2、计算化学势需要将粒子插入吸附的位置，并根据能量差异得到过化学势，再加上理想气体化学势。但是gcmc接受了一个构型之后，体系已经发生了变化，初始指定的化学势还合适吗？如果根据计算实时发生变化，是不是过于麻烦？
3、LAMMPS提供了一种忽略化学势的方法，即指定理想气藏压力的方法，并且可以指定逸度系数。纯水溶液也可以使用吗？在MS的sorption板块中，也只是输入逸度或者loading即可以进行吸附。

刘小胖 发表于 2023-10-24 13:22
楼主问一下甲烷的化学势怎么计算的 我在另一个软件上算出来化学势 转化成lammps的对不上

Lammps中fix的widom命令可以算过剩的化学势在加上理想化学势就行了

刘小胖 发表于 2023-10-24 13:22
楼主问一下甲烷的化学势怎么计算的 我在另一个软件上算出来化学势 转化成lammps的对不上

我的方法不一定对的，而且也很笨，就是不断调整化学势，直到平衡时甲烷的密度在目标密度附近波动

楼主问一下甲烷的化学势怎么计算的 我在另一个软件上算出来化学势 转化成lammps的对不上

已经解决了，fix gcmc中的group ID要和动态组的parent group相同，才能统计到gcmc插入的分子。