<HTML dir=ltr><HEAD>
<META http-equiv=Content-Type content="text/html; charset=unicode">
<META content="MSHTML 6.00.2800.1649" name=GENERATOR></HEAD>
<BODY>
<DIV><FONT face=Arial color=#000000 size=2>Dear Gromacs-users,</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT>&nbsp;</DIV>
<DIV><FONT face=Arial size=2>I would like to add a new bond type to the ffbonded.itp of the Gromos53a6 force field. For that,&nbsp;I tried to obtain the force constant for my new bond by performing vibrational frequency calculations with Gaussian. With that, I obtained the force constant in mDyne/A.</FONT></DIV>
<DIV><FONT face=Arial size=2>My problem is now: How can&nbsp;I convert this&nbsp;value to&nbsp;kJ/(mol*nm(4)) as required for the force field? I already tried, but&nbsp;I only end up at kJ/(mol*nm(2)).&nbsp;Also arround the web&nbsp;I could only find a factor for conversion of mDyne/A to&nbsp;kJ/(mol*nm(2)).&nbsp;Is there any trick or hint you can give me about that?</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT>&nbsp;</DIV>
<DIV><FONT face=Arial size=2>Thanks for your help,</FONT></DIV>
<DIV><FONT face=Arial size=2>Felix</FONT></DIV>
<DIV><FONT face=Arial size=2>&nbsp;&nbsp;</FONT></DIV></BODY></HTML>