Hi dear lammps users:
I want to simulate a nylon6 polymer with md and reach a suitable equilibrium state with temperature changes, when I run my system with the following input script:
units real
neigh_modify delay 2 every 1
atom_style full
bond_style harmonic
angle_style harmonic
pair_style soft 10.0
pair_modify mix arithmetic
kspace_style none
read_data Polyny6-15.data
################ lattice #######################
lattice sc 1
#################################################
log log.Soft
thermo 100
thermo_style custom step temp press etotal vol
fix Soft all nve 298.15 298.15
fix ShakeE1 water shake 1e-6 1000000 0 m 1.008 a 154
#################### dump ######################
dump dumpSoft all custom 1000 atoms.Soft.xyz id mol type q x y z
dump_modify dumpSoft sort id
restart 500000 restartSoft.*
run 1000000
I faced with this error:
ERROR: Could not find fix group ID water (…/modify.cpp:822)
Last command: fix ShakeE1 water shake 1e-6 1000000 0 m 1.008 a 154
I read lammps manual bun I can’t find my problem. I need any suggestion/advice about this errors and about something obvious mistake in my input script.
Thanks very much
Moein