; NVT equilibration (field-free)
integrator              = md
dt                      = 0.002
nsteps                  = 500000       ; set time as needed

constraints             = h-bonds
constraint-algorithm    = lincs
lincs_iter              = 1
lincs_order             = 4

cutoff-scheme           = Verlet
nstlist                 = 20
ns_type                 = grid
rlist                   = 1.2
rcoulomb                = 1.2
rvdw                    = 1.2
DispCorr                = EnerPres

coulombtype             = PME
pme_order               = 4
fourierspacing          = 0.12

; Temperature coupling (proper equil)
tcoupl                  = Berendsen
tc-grps                 = System
tau_t                   = 1.0
ref_t                   = 298

; No pressure coupling in NVT
pcoupl                  = no

; COM motion removal
nstcomm                 = 100
comm-mode               = linear
comm-grps               = System

; Output (light)
nstenergy               = 1000
nstlog                  = 1000
nstxout                 = 0
nstvout                 = 0
nstfout                 = 0
nstxout-compressed      = 5000         ; 10 ps during equil
compressed-x-precision  = 1000

pbc                     = xyz