I will try to calculate N2/CO2 selectivity for CO2 capture in MOF. I would like to compute gcmc in lammps for CO2 capture. Do you know how to make a ratio of molecule? For instance, I will give a CO2 to N2 molar ratio of 15 to 85. Should I change molar chemical potential because of this ratio?
If I change pressure, should I change molar chemical potential as well?
If you want to set the number of molecules, then just create the number of molecules you want. Your chemical potentials and pressure are not going to be controled though.
Thanks a lot for answer. I prefer to give a certain pressure, like 0.1 atm, and temperature (298.15 K). The system can control the number of molecules depending on a ratio (15:85). What is the command for it?
Are you aware of a Monte Carlo method that allows enforcing a strict molecule ratio? I’m not. In typical GCMC, you fix the chemical potentials (or partial pressures), and the composition results from that. LAMMPS won’t let you fix both the molecule numbers and the chemical potential at the same time, that wouldn’t make thermodynamic sense.
It’s not very different to have two species; you just need to use two GCMC commands. I think you don’t see many examples of such GCMC applications because there are faster codes than LAMMPS for GCMC, and people tend to perform these types of calculations using other software like GOMC or RASPA.
We will look at the Working Capacity, Isosteric Adsorption Heat (QST) and Selectivity (S) quantities to be obtained by simulating real operating conditions for N2, CO2 and H2O gas mixtures using MOF5. Can I make this in Lammps?
Actually, RASPA is too slow because it runs serially with one core. So I do not prefer to use RASPA. But I do not know GOMC.