[gmx-users] Error in NVT equilibration

[Justin Lemkul]

On 9/18/18 5:33 AM, Gonzalez Fernandez, Cristina wrote:
> Hi,
>
>
> As I am using the GROMOS froce field, I constrained all bonds as suggested. However, I obtained the same error.
>
>
> I copy the .mdp file in case there are errors and you can help me to find them:
>
>
> title                   = Lipid NVT equilibration
> define                  = -DPOSRES  ; position restrain the Lipid
> ; Run parameters
> integrator              = md        ; leap-frog integrator
> nsteps                  = 6667      ; 2 * 50000 = 100 ps
> dt                      = 0.015     ; 2 fs
> ; Output control
> nstxout                 = 500       ; save coordinates every 1.0 ps
> nstvout                 = 500       ; save velocities every 1.0 ps
> nstenergy               = 500       ; save energies every 1.0 ps
> nstlog                  = 500       ; update log file every 1.0 ps
> ; Bond parameters
> continuation            = no        ; first dynamics run
> constraint_algorithm    = lincs     ; holonomic constraints
> constraints             = all-bonds   ; bonds involving H are constrained
> ;shake_tol              = 0.0001
> lincs_iter            = 1            ; accuracy of LINCS
> lincs_order            = 4            ; also related to accuracy
> ; Nonbonded settings
> cutoff-scheme           = Verlet    ; Buffered neighbor searching
> ns_type                 = grid      ; search neighboring grid cells
> nstlist                 = 10        ; 20 fs, largely irrelevant with Verlet
> rlist                   = 1.4
> rcoulomb                = 1.2       ; short-range electrostatic cutoff (in nm)
> rvdw                    = 1.2       ; short-range van der Waals cutoff (in nm)
> DispCorr                = EnerPres  ; account for cut-off vdW scheme
> vdwtype                 = Cut-off
> vdw_modifier            = Force-switch

You're using nonbonded settings for the CHARMM force field, so you need 
to correct these settings to use the proper GROMOS values.

Beyond that, follow 
http://www.gromacs.org/Documentation/Terminology/Blowing_Up#Diagnosing_an_Unstable_System

-Justin

-- 
==================================================

Justin A. Lemkul, Ph.D.
Assistant Professor
Virginia Tech Department of Biochemistry

303 Engel Hall
340 West Campus Dr.
Blacksburg, VA 24061

jalemkul at vt.edu | (540) 231-3129
http://www.thelemkullab.com

==================================================