On 10/1/12 12:43 AM, Ravi Kumar Venkatraman wrote:
Dear All,
To find out gas phase torsional energetics I have restraints
dihedral and calculated the energies. But the Dihedral restraints are not
working properly, its changes +/- 2-5 degree from the specified value. Is
there any other way we can run torsional energetics are I have chenge the
modifications in the *.top and *.mdp files.
Restraints are simply biasing potentials; they do not fix any parameter at a
given value. The dihedral fluctuation is entirely normal. If you want to
reduce the fluctuations, increase the value of dihre_fc in the .mdp file.
-Justin
I made following modifications in my *.top file and *.mdp file to make
dihedral restraints.
#include "dab.itp"
[ dihedral_restraints ]
; ai aj ak al type label phi dphi kfac power
26 23 18 6 1 1 0 0 1 2
;dihedral restraints
dihre = yes
dihre_fc = 100
Thank you
*
With Regards,
Ravi Kumar Venkatraman,
IPC Dept., IISc,
Bangalore, INDIA.
+91-9686933963.*
--
========================================
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