Dear LAMMPS-USERS,
I tried to simulate a small model (lx=ly=lz=10) consists of 2 SPC/E water model inside two carbon atoms walls. When I run the simulation it works fine with SHAKE algorithm.
But, when I want to run a big simulation model (lx=ly=30, lz=20) with 800 SPC/E water I come across to the error as
fix ensemble_set all nvt temp 300.0 300.0 0.05
fix spce_model water shake 0.0001 20 0 b 1 a 1#
ERROR: Unknown fix style (…/modify.cpp:762)
Please share your comments.
Thanks in advance
My input
***************** Initialization *********************
units real
dimension 3
boundary p p p
atom_style full
read_data atm_data.dat
group oxygen type 1
group hydrogen type 2
group carbon type 3
group water type 1 2
set group oxygen charge -0.8476
set group hydrogen charge 0.4238
set group carbon charge 0.0000
#pair_style dpd/tstat 300.0 300.0 10.0 34387
#comm_style brick
#comm_modify vel yes
#pair_style hybrid lj/cut/coul/long 10.0 10.0 tersoff
pair_style lj/cut/coul/long 10.0 10.0
pair_coeff 1 1 0.155300 3.16900 #lj/cut/coul/long
pair_coeff 1 2 0.00000 0.00000
pair_coeff 1 3 0.093627 3.1900
pair_coeff 2 2 0.000000 0.0000
pair_coeff 2 3 0.000000 0.0000
pair_coeff 3 3 0.000000 0.0000
#pair_coeff * * tersoff SiC_1994.tersoff NULL NULL C
bond_style harmonic
bond_coeff 1 529.581 1.0
angle_style harmonic
angle_coeff 1 37.95 109.47
fix force_zero carbon setforce 0.0 0.0 0.0
fix vel_zero carbon move linear 0.0 0.0 0.0 units lattice
#velocity oxygen create 300.0 34387 rot yes dist gaussian # for Oxygen
#velocity hydrogen create 300.0 34387 rot yes dist gaussian # for Hydrogen
#velocity water create 300.0 34387 rot yes dist gaussian # for water
kspace_style ewald 1.0e-6
*************** Setting ******************************
neighbor 2.0 bin
neigh_modify delay 0 every 10 check yes
reset_timestep 0
thermo 100
thermo_style custom step elapsed dt time atoms
thermo_modify norm no flush yes
fix ensemble_set all nvt temp 300.0 300.0 0.05
fix spce_model water shake 0.0001 20 0 b 1 a 1#
#compute 1 all pe/atom
#minimize 1.0e-4 1.0e-6 100 1000
#min_style cg
dump trj all cfg 100 wat_lammpstrj.*.cfg id type xs ys zs
dump_modify trj element O H C
#write_data waterinfo.data
#write_restart waterinfo.restart
timestep 0.2
run 1000