[gmx-users] Meaningless pressure values: NVT and NVE ensembles
Francesco Mambretti
francesco.mambretti at studenti.unimi.it
Wed Nov 5 14:14:04 CET 2014
Sure! here it is. Am I missing anything?
Thanks
; RUN CONTROL PARAMETERS
integrator = md
; Start time and timestep in ps
tinit = 0
dt = 0.001
nsteps = 5000000
; mode for center of mass motion removal
comm-mode = Linear
; number of steps for center of mass motion removal
nstcomm = 1
; OUTPUT CONTROL OPTIONS
; Output frequency for coords (x), velocities (v) and forces (f)
nstxout = 0
nstvout = 0
nstfout = 0
; Output frequency and precision for xtc file
nstxtcout = 10000
xtc-precision = 1000
; Output frequency for energies to log file and energy file
nstlog = 1000
nstcalcenergy = 1000
nstenergy = 1000
;energygrps
; NEIGHBORSEARCHING PARAMETERS
; nblist update frequency
nstlist = 1
; ns algorithm (simple or grid)
ns-type = grid
; Periodic boundary conditions: xyz (default), no (vacuum)
; or full (infinite systems only)
pbc = xyz
;cutoff-scheme
cutoff-scheme = Verlet
verlet-buffer-tolerance = 0.0001
; OPTIONS FOR ELECTROSTATICS AND VDW
; Method for doing electrostatics
coulombtype = Cut-off
rcoulomb-switch = 0
coulomb-modifier = potential-shift
; cutoff 2.5*sigma = 2.5*0.34nm = 0.85nm 2.2*sigma =0.75
rcoulomb = 0.75
; Dielectric constant (DC) for cut-off or DC of reaction field
epsilon-r = 1
; Method for doing Van der Waals
vdw-type = Cut-off
vdw-modifier = potential-shift
; cut-off lengths
rvdw-switch = 0
rvdw = 0.75
; Apply long range dispersion corrections for Energy and Pressure
DispCorr = EnerPres
; OPTIONS FOR WEAK COUPLING ALGORITHMS
; Temperature coupling
tcoupl = No
; Pressure coupling
Pcoupl = No
; GENERATE VELOCITIES FOR STARTUP RUN
gen-vel = yes
gen-temp = 96
gen-seed = 2223
More information about the gromacs.org_gmx-users
mailing list