Hi,
I want to raise the temperature of a Polystyrebe melt systen that has run for 1000ns in 500K to 800K.Below is the gro.mdp I used to raise temperature with Temperature ramp.
I would like to know T ramp method is good to raise the T in a system or is there any other reommended ways to do that in gromacs?
In below script I used 800K in ref_t and gen_t and alsoannealing T as 500 800 based on my understanding.
Please advice.
; Time step in ps
dt = 0.001
; Number of steps to run
nsteps = 1000000000
; 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 = 800
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 = 800
;Use when ramping temperature
annealing = single;ramps a single time to the final T and remains constant at final T
annealing-npoints = 2 ;number of reference points (match the number of temps, ie 300K to 550K would be 2, one for each temp)
annealing-time = 0 1000 ;ramping time in ps - (1000 is 1 ns)
annealing-temp = 500 800 ;start and stop temp
; 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 = 1