Out of Range Atoms of Water Molecules in a Box

Dear all LAMMPS user,

I have been trying to simulate water molecules in a box with dimension of 30x30x30 A. I am using packmol to distribute the location of each water molecule then I have used minimize command to update the coordinate of each atom before I run the simulation.
The problem I met are:

  1. Error with computation of PPPM, it is said that the atom was out of range that lammps cannot compute the PPPM. I have been trying to fix this issue by doing some suggestion from the manual (like changing the neighbor_modify), but this issue still has not fixed. I tried to change the pair_style into lj/cut/tip4p/cut and delete the kspace. It worked but the simulation ended without any comment (capture.jpg).
  2. Warning of inconsistent image flags, I have no idea about this issue and the error part in my script.

also, I’d like to ask whether there is any other website or book that can be used to undertanding more about lammps.

Here is my input script:

dimension 3
boundary p p p
units real

atom_style full
bond_style harmonic
angle_style harmonic

pair_style lj/cut/tip4p/long 1 2 1 1 0.15 10.0
kspace_style pppm/tip4p 1.0e-5

Data-Water-Only.txt (161 KB)

Input-Water-Box.txt (1.02 KB)

Capture.JPG

My guess on the out-or-range error is that you have bad dynamics,

i.e. the water is moving too far too fast.

For image flags: from doc/Section_error.html

{Inconsistent image flags}

The image flags for a pair on bonded atoms appear to be inconsistent.
Inconsistent means that when the coordinates of the two atoms are
unwrapped using the image flags, the two atoms are far apart.
Specifically they are further apart than half a periodic box length.
Or they are more than a box length apart in a non-periodic dimension.
This is usually due to the initial data file not having correct image
flags for the 2 atoms in a bond that straddles a periodic boundary.
They should be different by 1 in that case. This is a warning because
inconsistent image flags will not cause problems for dynamics or most
LAMMPS simulations. However they can cause problems when such atoms
are used with the fix rigid or replicate commands.

also, I’d like to ask whether there is any other website or book that can be used to undertanding more >about lammps.

The manual is 1000+ pages, there is a doc/Developer.pdf, and the entire

web site. That’s all there is that I’m aware of.

Steve