I did not try running your script, but I am guessing that the small volume:
region zone1 block 0 10 0 10 0 10 units box
and the small GCMC target pressure
pressure 0.88 fugacity_coeff 0.99
result in a very small average particle count. If you want to see the time-average pressure reported by LAMMPS converge to your target pressure, you will need to do one of the following:
-
Increase target pressure
-
Increase box volume
-
Increase duration of your simulation
I suggest that pick the pressure and volume so that you get a few hundred argon atoms on average, and then invoke the fix many times i.e. nsteps/N should be large. Note that invoking the fix 1000/100=10 times with large values of X and M is not a good idea, because you only measure the particle count 10 times, even though you performed 10^7 exchange moves.
Aidan