<html> <head> <title>MCCCS Towhee (Dick and Ritchie 1994)</title> </head> <body bgcolor="#FFFFFF" text="#000000"> <table width="675" border="0" cellspacing="0" cellpadding="0" height="590"> <tr> <td colspan="2" height="29"> <div align="center"> <font size="5"> <b><font face="Arial, Helvetica, sans-serif"><a name="top"></a>MCCCS Towhee (Dick and Ritchie 1994)</font></b> </font> </div> </td> </tr> <A href="http://sourceforge.net"> <IMG src="http://sourceforge.net/sflogo.php?group_id=87039&type=5" width="210" height="62" border="0" alt="SourceForge.net Logo" align="right"/> </A> <tr> <td width="18" height="371" valign="top"> <p> </p> <p> </p> </td> <td width="697" valign="top"> <b>Overview</b> <ul> This section covers an extension to the <a href="weiner1986.html">Weiner <i>et al.</i> 1986</a> force field. Unlike most force fields, Dick and Ritchie 1994 is only useful in combination with the <a href="weiner1986.html">Weiner <i>et al.</i> 1986</a> force field. All of the additional Towhee atom types for the Dick and Ritchie 1994 extension are listed , along with a short description of their meanings. For more information about the Amber family of force fields see the <a href="http://amber.scripps.edu">Amber web site</a>. Note that Dick and Ritchie 1994 is a 12-6 plus 12-10 H-bond force field and can only be combined with other '12-6 plus 12-10 H-bond' force fields. Any discrepencies (especially typos) from the published force field values are the sole responsibility of Marcus G. Martin, and I welcome feedback on how this implementation compares with other programs. <p> </p> </ul> <b>References for Dick and Ritchie 1994</b> <ul> <ul> <li><a href="../references.html#dick_ritchie_1994">Dick and Ritchie 1994.</a></li> </ul> </ul> <b>Dick and Ritchie 1994 in Towhee</b> <ul> The official force field name for Dick and Ritchie 1996 in Towhee is 'Weiner1986'. This indicates that Dick and Ritchie 1994 is not a stand alone force field and is only useful in combination with the <a href="weiner1986.html">Weiner <i>et al.</i> 1986</a> force field. In order to use the Dick and Ritchie 1994 extensions you need to include the towhee_ff_Dick1994 file from the ForceFields directory. Here I list all of the atom names for use in the towhee_input file, along with a brief description taken from <a href="../references.html#dick_ritchie_1994">Dick and Ritchie 1994</a>. These atom types are designed to be mixed freely with the atom types described in the <a href="weiner1986.html">Weiner <i>et al.</i> 1986</a> documentation. Please note that the capitalization and spacing pattern is important and must be followed exactly as listed here. <ul> <dt><font color="red">Nitrogen</font> (atomic)</dt> <li><b>'NN'</b> : sp<sup>2</sup> nitrogen in a nitro group.</li> </ul> </ul> <b>Coulombic Interactions</b> <ul> Point charges for pentaerythritol tetranitrate (PETN) are provided in Table III of the <a href="../references.html#dick_ritchie_1994">Dick and Ritchie 1994</a> paper. However, the charges listed there add up to -0.003 for the PETN molecule so there must be some sort of typo in that table as the molecule should be net neutral. They also list different charges for the two symmetric oxygens on the nitro group, which might make sense for crystals, but does not make sense for fluids. Here I list their suggested charges for PETN and also the modifications I am currently using. <ul> <li>Atom: Dick and Ritchie 1996 charges [Martin unpublished charges]</li> <li>CT (central carbon): -0.479 [-0.476]</li> <li>CT (branch carbon): -0.030 [-0.030]</li> <li>OS (ether-like oxygen): -0.476 [-0.476]</li> <li>NN (nitro nitrogen): 0.771 [0.771]</li> <li>O (nitro oxygens): -0.346 and -0.388 [-0.367]</li> </ul> </ul> <b>Improper torsions</b> <ul> Improper torsions are not automatically generated by the Towhee code as the rules for determining where they are applied are not always straight-forward. Dick and Ritchie do not specify any additional improper torsion parameters so I do not believe any apply to the PETN molecule. For more information about improper torsions on other molecules please consult the <a href="weiner1986.html">Weiner <i>et al.</i> 1986</a> web manual. </ul> <a href="../towhee_capabilities.html">Return to the Towhee Capabilities web page</a> <p> </p> </td> </tr> </table> <hr width="715" align="left"> <i><font size="2">Send comments to:</font></i> <font size="2"> <a href="mailto:marcus_martin@users.sourceforge.net">Marcus G. Martin</a><br> <i>Last updated:</i> <!-- #BeginDate format:Am1 -->September 07, 2005<!-- #EndDate --> </font> <br> </body> </html>