GROMACS version: 2021.4
GROMACS modification: Yes/No
Here post your question
Hello,
Before proceeding to the production run, I was performing membrane equilibration. Since I expected the deprotonated lipid to migrate toward the membrane core, I planned to monitor box size changes as an indicator of system stabilization. I’ve previously had success using this method for equilibration.
Below is the .mdp
file I used:
title = KALP15-DPPC Production MD
; Run parameters
integrator = md
nsteps = 50000000
dt = 0.002
; Output control
nstxout = 0
nstvout = 0
nstxtcout = 30000
nstenergy = 30000
nstlog = 30000
; Bond parameters
continuation = yes
constraint_algorithm = lincs
constraints = h-bonds
lincs_iter = 1
lincs_order = 4
; Neighborsearching
ns_type = grid
nstlist = 5
rlist = 1.2
rcoulomb = 1.2
rvdw = 1.2
; Electrostatics
coulombtype = PME
pme_order = 4
fourierspacing = 0.16
; Temperature coupling
tcoupl = v-rescale
tc-grps = SOL_ION MEMB
tau_t = 0.5 0.5
ref_t = 300 300
; Pressure coupling
pcoupl = parrinello-rahman
pcoupltype = semiisotropic
tau_p = 12.0
ref_p = 1.0 1.0
compressibility = 3.0e-5 3.0e-5
; PBC and dispersion correction
pbc = xyz
DispCorr = EnerPres
; Velocity generation
gen_vel = no
; Center of mass motion removal
nstcomm = 1
comm-mode = Linear
comm-grps = SOL_ION MEMB
refcoord_scaling = com
At 300 K, the membrane appeared as follows (attached image).
At 300 K, the membrane appeared to be stable and there was no noticeable change in box size or lipid distribution. However, to accelerate equilibration, I increased the temperature to 373 K. Upon doing so, I observed a drastic shrinkage of the box size and elongation of the membrane in the Z-direction, as shown in the image below.
I’m wondering whether this is an expected effect of the elevated temperature or if there might be something incorrectly set in my .mdp
file. If this issue is due to the temperature increase, would it be more appropriate to simply extend the equilibration time at 300 K instead?
Thank you in advance for your insights!