Dear Lammps developers and users,
I’m going to model salt (CaSO4) precipitation in aqueous solution using LAMMPS, stable version 12 Aug 2018.
To do so, firstly a large system consisting of 66000 water molecules dissolving 10 CaSO4 was generated. Then, to model super-saturation, fix evaporate was employed to reduce the number of water molecules, that follows, and consequently increase the concentration of dissolved CaSO4:
fix 3 owater evaporate 100 9 surface 38277 molecules yes
here, owater is group_ID of the water’s oxygen atoms and “surface” is the whole simulation box. In this way, 3 water molecules would be deleted at every 100-time steps.
Due to the restriction imposed by mailing list on the size of attachments, I cannot send associating input and data file, but, it will be OK to send anyone if want.
Unfortunately, the simulation goes wrong! the number of removed atoms is more than the desired one, i.e., 9. Also, the deleting manner is apparently erratic, deleted number of molecules varies at every time step.
I tried that scenario for a smaller system with 6000 waters and 25 CaSO4 and surprisingly the system works well, deleting atoms as invoked by preceding fix evaporate.
I suppose this issue is a bug with the large system.
Any hint is welcome.
Many thanks in advance