Dear Christopher,<br>Considering the amount of human time vs. computer time, and if the molecules are smallish as you said, why don&#39;t you go for a QM method? Semi-empirical methods can treat prettly large molecules and you will not need to parameterize anything. <br>
Of course, this can&#39;t be done with MMTK. Before choorsing another software, make sure it has analytical second derivatives of the energy. Otherwise, frequencies are rather slow to calculate. If you do need to work with MMTK afterwards, you can always write a parser that reads the QM outputs and generates ScientificPython vectors. I think this is still less work than a thorough parameterization.<br>
Cheers,<br>Ramon<br><br><br><div class="gmail_quote">2011/7/7 Konrad Hinsen <span dir="ltr">&lt;<a href="mailto:research@khinsen.fastmail.net">research@khinsen.fastmail.net</a>&gt;</span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div class="im">On 6 Jul 2011, at 15:13, Christopher Drost wrote:<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
(1) The first question I have regards sp-hybridization. My molecule<br>
has both C-C and C-N triple bonds, which can be something like 50-100%<br>
stiffer than double bonds, and prefer a different bending geometry. So<br>
I have to specify the sp-hybridization somehow. I&#39;ve found<br>
documentation for specifying these in Allinger&#39;s MM models, but I&#39;m<br>
still uncertain how to specify a triple bond for the Amber99 force<br>
field as implemented in MMTK. Can anyone help with how this might be<br>
specified? I have trouble believing that Amber99 would be ~15 years<br>
behind the MM models.<br>
</blockquote>
<br></div>
Amber doesn&#39;t have the notion of a double or triple bond. The molecule topology specifies just a bond, without any parameters. The interactions are then specified in terms of atom types. To get a triple bond, you thus have to introduce new atom types for the atoms on either side of the bond, and then define an appropriate force constant for a bond between the atoms of the corresponding types. This would typically be done in the form of a &quot;modification file&quot; (in Amber terminology), i.e. an extension to the list of standard Amber parameters. MMTK uses standard Amber modification files, so you can follow the Amber documentation for the details:<br>

<br>
        <a href="http://ambermd.org/#ff" target="_blank">http://ambermd.org/#ff</a><div class="im"><br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
with an initial &quot;configuration&quot;. You might want to fix this to be a<br>
more verbose error message.<br>
</blockquote>
<br></div>
Indeed. There isn&#39;t much I can do about the fundamental problem that energy minimization requires an initial configuration, but I can add a check for this condition.<div class="im"><br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
(3) One limitation: my molecules have a net dipole moment. (Probably<br>
the benzene C-H bonds do as well?) I don&#39;t yet know how to calculate<br>
the proper &#39;amber_charge&#39; values in such a case.<br>
</blockquote>
<br></div>
There are no special precautions. If you know the dipole moment of benzene, you could use it as a constraint in the charge fitting procedure. Otherwise you just fit the charges and accept whatever dipole moment they result in.<div class="im">
<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
There is some documentation about the ChargeFit class but I can&#39;t really divine how<br>
it works, exactly. So, does anyone have an example for calculating<br>
these amber_charge parameters, something like the way that Amber&#39;s<br>
energy minimization algorithm can be used to calculate their<br>
positions?<br>
</blockquote>
<br></div>
Charge fitting first of all requires values for the electrostatic potential around the molecule, which you would typically get from quantum chemistry calculations. Do you have those? If yes, have a look at Examples/Miscellaneous/charge_<u></u>fit.py for a simple example of how to use the charge fit module.<br>

<br>
Of course it might be worth checking first if someone else has already computed Amber parameters for benzene. Or for the OPLS force field, which is also supported by MMTK and often a better choice for working with non-biological molecules. OPLS uses the same functional form as Amber, but has different atom types and different parameters.<br>

<br>
Konrad.<br>
--<br>
------------------------------<u></u>------------------------------<u></u>---------<br><font color="#888888">
Konrad Hinsen<br>
Centre de Biophysique Moléculaire, CNRS Orléans<br>
Synchrotron Soleil - Division Expériences<br>
Saint Aubin - BP 48<br>
91192 Gif sur Yvette Cedex, France<br>
Tel. +33-1 69 35 97 15<br>
E-Mail: research at  khinsen dot fastmail dot net<br>
------------------------------<u></u>------------------------------<u></u>---------</font><div><div></div><div class="h5"><br>
<br>
<br>
______________________________<u></u>_________________<br>
mmtk maillist  -  <a href="mailto:mmtk@starship.python.net" target="_blank">mmtk@starship.python.net</a><br>
<a href="http://starship.python.net/mailman/listinfo/mmtk" target="_blank">http://starship.python.net/<u></u>mailman/listinfo/mmtk</a><br>
</div></div></blockquote></div><br>