francesca vitalini wrote:
Tryed the double precision without any success... Again the same type of error message. Might need to try other minimization algorithm. Any suggestions on that?


What is the Fmax you achieve, and what are your intents for this process? Minimizing very low is only strictly necessary for certain types of very sensitive calculations. For normal MD, convergence below 100-1000 kJ/mol-nm is typical.

Other algorithms are listed in the manual.  They may or may not be necessary.

-Justin

2012/2/9 Justin A. Lemkul <jalem...@vt.edu <mailto:jalem...@vt.edu>>



    francesca vitalini wrote:

        Hi all!
        I'm trying to run some energy minimization with gromacs 3.3.1 on
        a system that I have transformed from cg to aa. After the
        reverse transformation I solvate the system and I try to
        minimize just the proteins (my system is a dimer) keeping fixed
        the water. Then the idea would be to minimize the water keeping
        the protein in harmonic distance constraints as I just want to
        obtain an aa system which describes the configuration obtained
        previously in the cg representation.
        So, I first use genrestr to create the .itp file for fixing the
        position of the water, then I run 10 ps of EM and what I get is
        the following error message:

        Stepsize too small, or no change in energy.
        Converged to machine precision,
        but not to the requested precision Fmax < 10

        Double precision normally gives you higher accuracy.
        You might need to increase your constraint accuracy, or turn
        off constraints alltogether (set constraints = none in mdp file)

         From what I've read in previous posts the problem might be
        either in the mdp file or in the topology, but I don't
        understand what have I made wrong. Here is my mdp file


    The other possibility is that there is no problem at all.  Some
    systems cannot be minimized below certain thresholds, based on the
    input geometry.  You may also want to try a different minimizer or
    double precision.

    
http://www.gromacs.org/__Documentation/Errors#Stepsize___too_small.2c_or_no_change_in___energy._Converged_to_machine___precision.2c_but_not_to_the___requested_precision
    
<http://www.gromacs.org/Documentation/Errors#Stepsize_too_small.2c_or_no_change_in_energy._Converged_to_machine_precision.2c_but_not_to_the_requested_precision>

    -Justin

         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          = -DEPOSRES
        ;define           = -DFLEXIBLE
        ; Parameters describing what to do, when to stop and what to save
        integrator      = steep         ; Algorithm (steep = steepest
        descent minimization)
        emtol           = 10.0          ; Stop minimization when the
        maximum force < 1.0 kJ/mol
        dt                       = 0.005
        ;nsteps                   = 2000
        nsteps          = 5000          ; 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
        nstlist         = 1             ; Frequency to update the
        neighbor list
        ns_type         = grid          ; Method to determine neighbor
        list (simple, grid)
        coulombtype     = Reaction-Field ; Treatment of long range
        electrostatic interactions
        epsilon_rf      = 78
        rcoulomb        = 1.4           ; long range electrostatic cut-off
        rvdw            = 1.4           ; long range Van der Waals cut-off
        constraints     = none          ; Bond types to replace by
        constraints
        pbc             = xyz           ; Periodic Boundary Conditions
        (yes/no)

        table-extension = 1.2

        and this is my topology

        ;
        ;       File 'dynamin_dimer_fg.top' was generated
        ;       By user: vitalini (213883)
        ;       On host: goat
        ;       At date: Mon Jan 30 11:02:21 2012
        ;
        ;       This is your topology file
        ;       "The Poodle Chews It" (F. Zappa)
        ;
        ; Include forcefield parameters
        #include
        
"/home/cocktail/vitalini/__gromacs_special/share/gromacs/__top/ffG53a6m.itp"

        ; Include chain topologies
        #include "dynamin_dimer_fg_A.itp"
        #include "dynamin_dimer_fg_B.itp"

        ; Include water topology
        #include
        "/home/cocktail/vitalini/__gromacs_special/share/gromacs/__top/spc.itp"

        #ifdef POSRES_WATER
        ; Position restraint for each water oxygen
        [ position_restraints ]
        ;  i funct       fcx        fcy        fcz
          1    1       1000       1000       1000
        #endif


        ;Include restraints on water
        #ifdef POSRES_WAT
        #include "posre_wat.itp"
        #endif

        ; Include generic topology for ions
        #include
        "/home/cocktail/vitalini/__gromacs_special/share/gromacs/__top/ions.itp"

        [ system ]
        ; Name
        Protein in water

        [ molecules ]
        ; Compound        #mols
        Protein_A           1
        Protein_B           1
        ;SOL              390842
        SOL             390820
        NA+                22

        while this is the posre_wat.itp file I generated with genrestr

        ; position restraints for non-Protein of Protein in water

        [ position_restraints ]
        ;  i funct       fcx        fcy        fcz
        15231    1       1000       1000       1000
        15232    1       1000       1000       1000
        15233    1       1000       1000       1000
        15234    1       1000       1000       1000
        15235    1       1000       1000       1000
        15236    1       1000       1000       1000
        15237    1       1000       1000       1000
        15238    1       1000       1000       1000
        15239    1       1000       1000       1000
        15240    1       1000       1000       1000
        15241    1       1000       1000       1000
        15242    1       1000       1000       1000
        15243    1       1000       1000       1000
        15244    1       1000       1000       1000
        15245    1       1000       1000       1000
        15246    1       1000       1000       1000
        15247    1       1000       1000       1000
        15248    1       1000       1000       1000
        15249    1       1000       1000       1000
        15250    1       1000       1000       1000
        15251    1       1000       1000       1000
        ..........................
        187712    1       1000       1000       1000

        Any help?
        Thanks Francesca


-- Francesca Vitalini

        PhD student at Computational Molecular Biology Group,
        Department of Mathematics and Informatics, FU-Berlin
        Arnimallee 6 14195 Berlin

        vital...@zedat.fu-berlin.de <mailto:vital...@zedat.fu-berlin.de>
        <mailto:vital...@zedat.fu-__berlin.de
        <mailto:vital...@zedat.fu-berlin.de>>
        francesca.vitalini@fu-berlin.__de
        <mailto:francesca.vital...@fu-berlin.de>
        <mailto:francesca.vitalini@fu-__berlin.de
        <mailto:francesca.vital...@fu-berlin.de>>

        +49 3083875776 <tel:%2B49%203083875776>
        +49 3083875412 <tel:%2B49%203083875412>


-- ==============================__==========

    Justin A. Lemkul
    Ph.D. Candidate
    ICTAS Doctoral Scholar
    MILES-IGERT Trainee
    Department of Biochemistry
    Virginia Tech
    Blacksburg, VA
    jalemkul[at]vt.edu <http://vt.edu> | (540) 231-9080
    <tel:%28540%29%20231-9080>
    http://www.bevanlab.biochem.__vt.edu/Pages/Personal/justin
    <http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin>

    ==============================__==========
-- gmx-users mailing list gmx-users@gromacs.org
    <mailto:gmx-users@gromacs.org>
    http://lists.gromacs.org/__mailman/listinfo/gmx-users
    <http://lists.gromacs.org/mailman/listinfo/gmx-users>
    Please search the archive at
    http://www.gromacs.org/__Support/Mailing_Lists/Search
    <http://www.gromacs.org/Support/Mailing_Lists/Search> before posting!
    Please don't post (un)subscribe requests to the list. Use the www
    interface or send it to gmx-users-requ...@gromacs.org
    <mailto:gmx-users-requ...@gromacs.org>.
    Can't post? Read http://www.gromacs.org/__Support/Mailing_Lists
    <http://www.gromacs.org/Support/Mailing_Lists>




--
Francesca Vitalini

PhD student at Computational Molecular Biology Group,
Department of Mathematics and Informatics, FU-Berlin
Arnimallee 6 14195 Berlin

vital...@zedat.fu-berlin.de <mailto:vital...@zedat.fu-berlin.de>
francesca.vital...@fu-berlin.de <mailto:francesca.vital...@fu-berlin.de>

+49 3083875776
+49 3083875412


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

Justin A. Lemkul
Ph.D. Candidate
ICTAS Doctoral Scholar
MILES-IGERT Trainee
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin

========================================
--
gmx-users mailing list    gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

Reply via email to