GROMACS version:
GROMACS modification: Yes/No
Here post your question
Hi,I have run lammps diblock PS-PEP simulation in dt=0.5fs.The system has compressed and run in NPT and then in NVE in lammps.
When I convert to gromacs in NVE I can not run it with dt=0.002=2fs with constraints and lincs it gave below error:
Fatal error:
3 particles communicated to PME rank 82 are more than 2/3 times the cut-off
out of the domain decomposition cell of their charge group in dimension y.
This usually means that your system is not well equilibrated.
In gromacs I could run dt=0.0005 =0.5 fs without constraints in NVE for few nanoseconds.
I want to know is it ok to continue the system in dt=0.0005 =0.5 fs without constraints.Will that affect my simulation?how can I convert the simulation to dt=0.002=2fs?
Below is the gro.mdp file with dt=0.oo2 that gives above error:
; Time step in ps
dt = 0.002
; Number of steps to run
nsteps = 100000000
; Remove center of mass motion every 100 steps
comm-mode = Linear
nstcomm = 100
; Write positions and velocities to the trr file every 5000 steps
; Output frequency for coords (x), velocities (v) and forces (f)
nstxout = 500000
nstvout = 500000
; Write to the log and energy file every 5000 steps
nstlog = 500000
nstenergy = 500000
; Calculate energies every 100 steps
nstcalcenergy = 100
; Write positions to three decimal places to xtc file every 5000 steps
nstxout-compressed = 500000
compressed-x-precision = 1000
; Use the Verlet cutoff scheme with grid neighbor searching; you should basically always use this when you can (the only time you wouldn’t would be if you were using tabulated potentials, buckingham interactions, or energy group exclusions; see manual for more detail)
cutoff-scheme = Verlet
ns_type = grid
; Periodic in all three dimensions
pbc = xyz
; Minimum cutoffs in nm for neighbor list, electrostatics, and LJ interactions. mdrun will tune rcoulomb for improved performance.
rlist = 1.2
rcoulomb = 1.2
rvdw = 1.2
; Use the PME algorithm (an alternative to PPPM) (this is very likely the electrostatics algorithm you want to use)
coulombtype = PME
; Smoothly shift vdw interactions to zero at the cutoff
vdw-type = Cut-off
vdw-modifier = Potential-shift-Verlet
; GENERATE VELOCITIES FOR STARTUP RUN =
gen-vel = yes
gen-temp = 500
gen-seed = 173529
; Use the v-rescale thermostat with a time constant of 0.1 ps and temperature of 600. V-rescale correctly samples the canonical distribution. Nose-hoover would be another good choice.
tcoupl = V-rescale
tc-grps = System
tau_t = 0.1
ref_t = 500
; Don’t use pressure coupling. For NPT this should be set. For equilibration, Berendsen is best since it is more stable. For production runs, Parrinello-Rahman is best since it samples the correct ensemble.
pcoupl = no
pcoupltype = Isotropic
tau-p = 1
refcoord-scaling = No
; Use the LINCS algorithm to turn all bonds with hydrogen atoms into constraints. This allows a larger timestep. (can also use shake if needed, but lincs is faster)
constraint_algorithm = lincs
constraints = H-bonds
; Do not generate initial velocities; rather this run is a continuation from a previous run. Alternately, set gen_vel to yes and give it a temperature to sample the maxwell distribution at, and set continuation to no.
;gen_vel = no
;continuation = yes
; Dielectric constant (DC) for cut-off or DC of reaction field =
epsilon-r = 30
~
~