water_box288.lmp

#LAMMPS Input file generated by Avogadro
# Title

# Intialization
read_data      water_box.lmpdat
units          real
dimension      3
boundary       p p p
atom_style     full

# Atom Definition
replicate      1 1 1

#The SPC/E water potential
pair_style      lj/cut/coul/long 9.8 9.8
kspace_style    pppm 1.0e-5
pair_coeff      1 1 0.15535 3.5533
pair_coeff      * 2 0.00000 0.0000
bond_style      harmonic
angle_style     harmonic
dihedral_style  none
improper_style  none
bond_coeff      1 100.00   1.000
angle_coeff     1 100.00 109.47
special_bonds   lj/coul 0.0 0.0 0.5

# Minimize 

minimize 1.e-4 1.e-6 1000 10000

# Constraints (SPC/E)
fix             RigidOHBonds all shake 0.0001 20 0 b 1 a 1

# Settings
velocity       all create 288 4928459 rot yes mom yes dist gaussian
fix            ensemble all npt temp 288 288 200

timestep       2.0

# Output
dump           traj all dcd 1000 water_box288.dcd
thermo_style   one
thermo         50

# Run the simulation
run            100000