<div dir="ltr"><div class="gmail_default" style="font-family:georgia,serif">Hello Dr. Meng,<br><br></div><div class="gmail_default" style="font-family:georgia,serif">Thanks for your explanations on this. Yes, of course Chimera really isn't even intended to handle the sort of molecules of interest. <br><br>The transition metal complexes have organic ligands with asymmetric centers attached. The idea was to compare the ligand steric interactions while considering the point chirality about the central tetrahedral copper(I) center. Not so much molecular orbitals or electronics.<br><br></div><div class="gmail_default" style="font-family:georgia,serif">As an organic chemist, I'm more accustomed to seeing energies having positive values, i.e. stability increases as energy drops. Since Chimera optimizes according to energies becoming increasingly negative, then it looks like you're using the opposite scale. For example, for -20 vs -14 kJ/mol , -20 is more stable than -14.<br><br></div><div class="gmail_default" style="font-family:georgia,serif">Thanks again and have a nice evening.<br><br></div><div class="gmail_default" style="font-family:georgia,serif">Fenton<br></div></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><span style="font-family:georgia,serif"><font size="2"><i>Dr. Fenton
Heirtzler</i><br>-------------------------------<br>Phone: +1-301-312-5145</font></span></div><div><span style="font-family:georgia,serif"><a href="http://www.linkedin.com/in/fentonhorgchem/" target="_blank">LinkedIn</a></span></div></div></div></div></div></div></div></div></div></div></div></div>
<br><div class="gmail_quote">On Mon, Jun 5, 2017 at 4:25 PM, Elaine Meng <span dir="ltr"><<a href="mailto:meng@cgl.ucsf.edu" target="_blank">meng@cgl.ucsf.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hello Dr. Fenton Heirtzler,<br>
I agree it can be difficult to find specifics amongst a lot of text.  I don’t know if this is what you meant, but from the “Minimize Structure” page:<br>
<<a href="http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/minimize/minimize.html" rel="noreferrer" target="_blank">http://www.rbvi.ucsf.edu/<wbr>chimera/docs/<wbr>ContributedSoftware/minimize/<wbr>minimize.html</a>><br>
<br>
"Energies (kJ/mol) are reported in the Reply Log. **Step numbers reported by MMTK are 2 greater than the actual numbers of minimization steps performed. The additional “steps” are not minimization steps but operations required to obtain gradient values and updated coordinates.**”<br>
<br>
As I understand it there is no energy decomposition available directly from MMTK to Chimera, and thus only the total energy is reported.  Total energies may be positive or negative depending on whether the positive interactions outweigh the negative interactions and conformational strain.  As there is no decomposition, there is no way to tell how much of the total is from the strain vs. nonbonded interactions, or electrostatic vs. VDW interactions.<br>
<br>
That said, my understanding is that simple molecular mechanics with atoms as nonpolarizable point charges is unlikely to provide an accurate representation of transition metal complexes, and Antechamber parametrization was not developed to handle such cases.  This reply from several years ago somewhat addresses this issue, although the URL within it is no longer valid:<br>
<<a href="http://www.cgl.ucsf.edu/pipermail/chimera-users/2007-June/001665.html" rel="noreferrer" target="_blank">http://www.cgl.ucsf.edu/<wbr>pipermail/chimera-users/2007-<wbr>June/001665.html</a>><br>
<br>
In short, you may need to use some other program to accurately model these situations.  I hope this clarifies the situation,<br>
Elaine<br>
----------<br>
Elaine C. Meng, Ph.D.<br>
UCSF Chimera(X) team<br>
Department of Pharmaceutical Chemistry<br>
University of California, San Francisco<br>
<div class="HOEnZb"><div class="h5"><br>
<br>
<br>
> On Jun 5, 2017, at 11:04 AM, Fenton Heirtzler <<a href="mailto:heirtzlerchempro@gmail.com">heirtzlerchempro@gmail.com</a>> wrote:<br>
><br>
> Hello,<br>
> I'm using Chimera to optimize the structures of several small-molecule transition metal complexes. The goal is to compare the steric energies of different diastereomers.<br>
><br>
> I​'ve searched the local documentation on version 1.11.2 (Macintosh) for the term "Minimize structure" but found nothing which actually described the meaning,number format,  reference points and units which are produced by the "Minimize Structure" command and appear in the "Reply log".<br>
><br>
> For example, the potential energy values which are returned can be either < 0 or > 0 (!!)<br>
><br>
> If a minimization is occurring, then the force field will describe those parameters - but the preceding values make no sense.<br>
><br>
> Some advice and instruction would be appreciated.<br>
><br>
> Thanks<br>
> Dr. Fenton Heirtzler<br>
<br>
</div></div></blockquote></div><br></div>