Dear All,
I am using lammps-20Dec07 version on an Intel cluster.
I am working on a molecular liquid (fully flexible + LJ + Coulomb). I find the “etot” to be constant of motion (1 part in 10^4) for the following two runs.
- NVT with run_style verlet (dt=0.5fs)
- NVE with “run_style respa” with the following attributes (dt=1.0fs)
run_style respa 2 2 bond 1 angle 1 dihedral 1 improper 1 pair 2 kspace 2
However, under NVT conditions with respa turned on (same respa break-up as above, and dt=1.0fs), etot is stable (no drift), but with a much larger RMS deviation (around 3 parts in 100). I provide below the input file and the top portion of the parm file for this run.
I would appreciate your inputs.
thanks,
bala
==============INPUT File==========
units real
dimension 3
atom_style full
read_restart run01.rest
include parm.bmimpf6_256
neighbor 1.0 bin
neigh_modify delay 0 every 20 check no
timestep 1.0
run_style respa 2 2 bond 1 angle 1 dihedral 1 improper 1 pair 2 kspace 2
velocity all create 300.0 23482341 mom yes rot yes dist gaussian
fix 1 all nvt 300.0 300.0 500.0
fix_modify 1 energy yes
thermo_style custom step vol temp etotal pe evdwl ecoul elong ebond eangle edihed eimp f_1
thermo 10
thermo_modify flush yes
dump 1 all custom 20000 dump.bmimpf6_256 tag mol type xu yu zu
restart 5000 rA.rest rB.rest
run 20000
write_restart r.rest