[gmx-users] Related to energy minimization
Bratin Kumar Das
177cy500.bratin at nitk.edu.in
Mon Sep 10 11:53:37 CEST 2018
I think you should remove define=-DFLEXIBLE
and define=-DPOSRES in the .mdp file. Give 1000 force constant. Don't
give large force costant.
On Mon, Sep 10, 2018, 3:14 PM ISHRAT JAHAN <jishrat17 at gmail.com> wrote:
> Dear all,
> I am trying to do energy minimization of two docked protein. At the energy
> minimization step during MD simulation process, the em.gro file is found to
> be broken due to which equilibration of protein does not occur and core
> dumped error has been found.
> Please help me in this regard. I am attaching my em.mdp file and snapshot
> of em.gro file.
> ; Lines starting with ';' ARE COMMENTS
> ; Everything following ';' is also comment
>
> title = Energy Minimization ; Title of run
>
> ; The following line tell the program the standard locations where to find
> certain files
> cpp = /lib/cpp ; Preprocessor
>
> ; Define can be used to control processes
> define = -DFLEXIBLE
>
> ; Parameters describing what to do, when to stop and what to save
> integrator = steep ; Algorithm (steep = steepest descent
> minimization)
> emtol = 1000.0 ; Stop minimization when the
> maximum force < 1.0 kJ/mol
> nsteps = 2000 ; Maximum number of (minimization) steps to
> perform
> nstenergy = 1 ; Write energies to disk every nstenergy
> steps
> energygrps = system ; Which energy group(s) to write to disk
>
> ; Parameters describing how to find the neighbors of each atom and how to
> calculate the interactions
> ns_type = grid ; Method to determine neighbor list (simple, grid)
> coulombtype = PME ; Longrange electrostatics (Ewald)
> rvdw = 1.0
> rlist = 1.0
> rcoulomb = 1.0
> fourierspacing = 0.12
> pme_order = 4
> ewald_rtol = 1e-5
> constraints = none ; Bond types to replace by constraints
> pbc = xyz ; Periodic Boundary Conditions (yes/no)
> [image: Screenshot_from_2018-09-10_14_59_11.jpg]
>
> --
> Ishrat Jahan
> Research Scholar
> Department Of Chemistry
> A.M.U Aligarh
> --
> Gromacs Users mailing list
>
> * Please search the archive at
> http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
> posting!
>
> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>
> * For (un)subscribe requests visit
> https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or
> send a mail to gmx-users-request at gromacs.org.
More information about the gromacs.org_gmx-users
mailing list