On 1/29/13 4:43 AM, SANTU BISWAS wrote:
dear users,
I am performing an energy minimization of
polypeptide(formed by 10-residues of alanine, lysine,leucine,glutamic
acid ) in vacuum box by using
Steepest Descent(initially) and then Conjugate Gradient methods in
gromacs_4.5.5_doubleprecision
using charmm-27 force field for all system.I saw that for lysine and
alanine potential energy
value is positive but for leucine and glutamic acid potential energy
value is nagative for both energy-
minimisation methods.
I do not understand why this is happening.
How are you extracting this information? You don't even have energygrps set in
the .mdp file, so I'm not sure how you've got energies per amino acid, which are
relatively useless numbers anyway. In vacuo, you probably have a lot of
unsatisfied interactions that would not occur in the condensed phase. Positive
potential energies can thus arise. As long as the EM converged to an acceptable
value of Fmax, that's what you should be most concerned with.
-Justin
Input file for energy-minimisation is
;title =
cpp = /lib/cpp
;include =-I../top/
define = -DFLEXIBLE
; Input file
integrator = cg
nsteps = 10000
pbc = no
ns_type = simple
nstlist = 0
rlist = 0
rcoulomb = 0
rvdw = 0
; Energy minimizing stuff
emtol = 0.001
emstep = 0.1
nstcgsteep = 10
---
santu
--
========================================
Justin A. Lemkul, Ph.D.
Research Scientist
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