I am trying to calculate the radial distribution function (rdf) for
TIP3P water molecules: atom info is in the file wa.atoms. I want
to equilibrate with an npt calculation, and then calculate the
rdf with nvt. Based on the manual, I assume that my "pair_style"
must include lj/cut and coul/cut, even though I would prefer
lj/cut/coul/long. (I am taking short cuts in this input file
while I debug the rdf error.)
I am getting an "ERROR: Illegal fix style" in the LAMMPS output file
when I use the following LAMMPS input file (hydrogen-1, oxygen-2):
dimension 3
boundary p p p
units real
atom_style full
pair_style lj/cut/coul/cut 5.0 10.0
bond_style harmonic
angle_style harmonic
read_data wa.atoms
fix 1 all shake 0.001 10 100 b 1 a 1
fix 2 all npt 300.0 300.0 100.0 xyz 1.0 1.0 1000.0 drag 0.2
velocity water create 300.0 527245 rot yes dist gaussian
timestep 2.0
run 1000
unfix 2
fix 3 all nvt 300.0 300.0 100.0
fix 4 all rdf 20 wa.rdf 100 2 2
run 2000
I know that I could do the rdf calculation using the dump file, but the
rdf is for a qualitative check, and I don't want to expend the effort for
a post-processor rdf calculation, if I don't have to.
Ken.