Can you tell me whether the pair style for tip4p water model works well?
I’ve tested the pair style tip4p water model of lammps. It seems that due to the pppm method, the tip4p model is very difficult to use. The program always complains “ERROR: Out of range atoms - cannot compute PPPM”. So at first I guessed the time step was too large. But it only works when the time step is extremely small (about 0.01 fs).
Below is my input file. The file input.lammps is in the attachment of this mail.
input.lammps (60.9 KB)