<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;"><BR>
<DIV>Hi Mark and everyone!</DIV>
<DIV> </DIV>
<DIV>This has been a previous topic though and thanks Mark for the reply.</DIV>
<DIV>What I need to do is to calculate the potential energy (energy minimzation) of a molecule at every fixed pairs of dihedrals (phi, psi). I did the suggestions on <A href="http://www.gromacs.org/Documentation/How-tos/Dihedral_Restraints"><U><FONT color=#0000ff>http://www.gromacs.org/Documentation/How-tos/Dihedral_Restraints</U></FONT></A>. In fact what follows is the entry I placed in my topol.top </DIV>
<DIV> </DIV>
<DIV>[ dihedral_restraints ]</DIV><SPAN lang="">
<DIV>; ai aj ak al type label phi dphi kfac power</DIV></SPAN>
<DIV>5 1 40 38 1 1 30 0 1 2 </DIV>
<DIV>12 40 38 34 1 1 10 0 1 2 </DIV>
<DIV> </DIV>
<DIV>I also added the following in my mdp file </DIV>
<DIV>dihre = yes</DIV>
<DIV>dihre_fc = 10000</DIV>
<DIV> </DIV>
<DIV>I removed disre_tau and nstdihreout as they are no longer needed in the recent gromacs versions. Unfortunately, after checking the dihedral angles, they appear to be far from the set value of (phi=30, psi=10). I understand that we are restraining the angles and not constraining/fixing them. Is there a way we can have the resulting phi and psi as close to the required values? I get values 15 degrees away from the original settings. Your help is greatly appreciated. </DIV>
<DIV>One more question. This restraining will not affect the final potential energy right?</DIV>
<DIV>Bernard</DIV><BR>--- On <B>Mon, 6/27/11, Mark Abraham <I><Mark.Abraham@anu.edu.au></I></B> wrote:<BR>
<BLOCKQUOTE style="BORDER-LEFT: rgb(16,16,255) 2px solid; PADDING-LEFT: 5px; MARGIN-LEFT: 5px"><BR>From: Mark Abraham <Mark.Abraham@anu.edu.au><BR>Subject: Re: [gmx-users] placing dihedral constraints<BR>To: "Discussion list for GROMACS users" <gmx-users@gromacs.org><BR>Date: Monday, June 27, 2011, 5:38 PM<BR><BR>
<DIV id=yiv350602181>On 06/27/2011 07:31 PM, Mr Bernard Ramos wrote:
<BLOCKQUOTE type="cite">
<TABLE border=0 cellSpacing=0 cellPadding=0>
<TBODY>
<TR>
<TD vAlign=top>Hi!
<DIV><BR></DIV>
<DIV>I would like to freeze two dihedrals/torsion angles but allow the rest (i.e. bonds, angles) to relax during energy minimization and MD. How do I go over applying constraints? Thanks.</DIV></TD></TR></TBODY></TABLE></BLOCKQUOTE><BR>See <A href="http://www.gromacs.org/Documentation/Terminology/Constraints_and_Restraints" rel=nofollow target=_blank>http://www.gromacs.org/Documentation/Terminology/Constraints_and_Restraints</A> to sort out some nomenclature. You probably want a dihedral restraint with a strong force constant (though other approaches are possible). See relevant parts of chapters 4 and 5. Work by analogy from the [dihedrals] and [position_restraints] sections in your existing [moleculetype].<BR><BR>Mark<BR></DIV><BR>-----Inline Attachment Follows-----<BR><BR>
<DIV class=plainMail>-- <BR>gmx-users mailing list <A href="http://us.mc1615.mail.yahoo.com/mc/compose?to=gmx-users@gromacs.org" ymailto="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</A><BR><A href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target=_blank>http://lists.gromacs.org/mailman/listinfo/gmx-users</A><BR>Please search the archive at <A href="http://www.gromacs.org/Support/Mailing_Lists/Search" target=_blank>http://www.gromacs.org/Support/Mailing_Lists/Search</A> before posting!<BR>Please don't post (un)subscribe requests to the list. Use the <BR>www interface or send it to <A href="http://us.mc1615.mail.yahoo.com/mc/compose?to=gmx-users-request@gromacs.org" ymailto="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</A>.<BR>Can't post? Read <A href="http://www.gromacs.org/Support/Mailing_Lists" target=_blank>http://www.gromacs.org/Support/Mailing_Lists</A></DIV></BLOCKQUOTE></td></tr></table>