I posted a similar question on fix gcmc as I’m too trying to adsorb water into a porous substrate and having issues with it. I’m using 31 March 2017 version of lammps. The current version has examples on gcmc but that’s for a single system of CO2 where the bond types and angle types are similar in the host system and the new molecules. For such a system I have no trouble running gcmc and it works fine for just water molecules in my case (as Axel mentioned with preexisting water in the box). However, when I try to adsorb water in a different system (for example Ca-silicate) then lammps gives ERROR: Molecule auto special bond generation overflow. For the host system I have atom types 1-4 with no bonds or angles. For the incoming SPC water molecules I have bonds and angles with types 5 and 6 for O and H of water respectively. I kept the provision for these extra types in the substrate data file. But I’m still seeing this error. Is there something else that need to be mentioned here? I am copying my inputfile with the water molecule and the data file of the host system that I’m reading with read_data. Please let me know what I’m doing wrong about this. I have very little experience with fix gcmc. I tried looking up for similar questions in the lammps user archive but couldn’t find a solution for my problem. So now I definitely need help on this from the lammps users who have experience with fix gcmc. Any example for adsorbing water or some other molecule into a different substrate would be also fine. Thanks in advance.
LAMMPS INPUT DECK
units metal
atom_style full
bond_style harmonic
dihedral_style none
improper_style none
angle_style harmonic
pair_style lj/cut 12.00
read_data data.anhydrous
molecule watermol mol.water
group water type 5 6
mass 1 28.065 #Si
mass 2 16.00 #O
mass 3 40.078 #Ca1
mass 4 40.078 #Ca2
mass 5 16.00 #Ow
mass 6 1.00 #Hw
set type 1 charge 1.72
set type 2 charge -1.1608219
set type 3 charge 1.43
set type 4 charge 1.70
set type 5 charge -0.82
set type 6 charge 0.41
pair_coeff * * 0.000000 0.100000 0.0000000
pair_coeff 1 2 0.2433e-4 3.6716
pair_coeff 2 2 0.5392e-1 3.0687
pair_coeff 2 3 0.3773e-4 4.898
pair_coeff 2 4 0.6331e-4 4.898
pair_coeff 1 5 0.2299e-4 3.6298
pair_coeff 3 5 0.3802e-4 4.898
pair_coeff 4 5 0.2617e-4 5.0168
pair_coeff 5 5 6.74E-03 3.5532 #from clayff
pair_coeff 2 5 0.2278e-3 4.7557
#H2O SPC/Fw
bond_coeff 1 45.93 1.012
angle_coeff 1 3.29136 113.24
kspace_style none
neighbor 1.0 bin
neigh_modify delay 4 every 2 check yes
timestep 0.001
dump id all atom 1000 dump.${name}
thermo 10
thermo_style custom step pe ke temp press density vol atoms
variable tfac equal 5.0/3.0
fix g1 all gcmc 10 100 100 0 27868 300 0.1 0.01 mol watermol tfac_insert ${tfac} group water #rigid m1 # maxangle 180 full_energy
fix r3 all nvt temp 300.0 300.0 1.0
run 1000
DATA FILE (data.anhydrous)
LAMMPS data file written using python script
70 atoms
6 atom types
0 bonds
1 bond types
0 angles
1 angle types
-0.849377 28.998874 xlo xhi
-0.018042 29.630416 ylo yhi
-1.064126 29.490964 zlo zhi
Masses
1 26.980000
2 16.000000
3 40.078000
4 40.078000
5 16.000000
6 1.000000
Atoms
1 0 1 0.000000 11.944669 5.947721 18.383305
2 0 2 0.000000 13.540550 6.356191 19.344793
3 0 2 0.000000 12.390029 5.412430 16.601149
4 0 2 0.000000 10.828668 7.498327 18.311030
5 0 2 0.000000 11.015977 4.531274 19.260118
6 0 1 0.000000 26.455698 5.114050 9.789356
7 0 2 0.000000 24.884309 4.843987 8.742483
8 0 2 0.000000 26.518044 3.819125 11.196568
9 0 2 0.000000 26.420799 6.866727 10.552225
10 0 2 0.000000 27.998874 4.927879 8.684124
11 0 1 0.000000 10.155061 15.784010 6.449462
12 0 2 0.000000 9.578487 17.072908 7.731821
13 0 2 0.000000 8.630115 14.982098 5.617330
14 0 2 0.000000 11.177258 14.426079 7.324791
…
MOLECULE FILE
A single water molecule file
3 atoms
2 bonds
1 angles
Masses
1 1.000000
2 1.000000
3 16.00000
Coords
1 0.000000 0.000000 0.000000
2 0 0.000000 1.633045 0.000000
3 0.577415 0.816522 0.000000
Types
1 6
2 6
3 5
Charges
1 0.41
2 0.41
3 -0.82
Bonds
1 1 1 3
2 1 2 3
Angles
1 1 1 3 2