<html>
<body>
I would be extra cautious before using any of the early OPLS ion
parameters which are derived<br>
from really early Hartree-Fock calculations.&nbsp; Note that in standard
OPLS usage many of the aforementioned<br>
parameters are superceded by OPLS types 406-410 which are aqvist-derived
(and therefore are almost<br>
identical to the aqvist-derived AMBER ions except for the differing
combination rule).&nbsp; <br><br>
If you check out the 1984 OPLS paper you will also notice that these
parameters are meant <i>only</i> to be<br>
used with TIP4P (to get the right free energy of solvation), while aqvist
is meant to be used with TIP3P/SPC.<br><br>
I've been doing similar literature digging,&nbsp; and I suspect the
reason there was early consensus on Cl- is simply<br>
because there was lots of scattering data available at that time of
aqueous chloride salts but none on fluoride salts:<br>
the location of the first ion-water oxygen peak in the RDF seems to be
the primary sanity check that a reasonable size<br>
was chosen. Without a size constraint, you can imagine getting the
&quot;right&quot;&nbsp; free energy of hydration with many <br>
degenerate pairs of sigma and epsilon.<br><br>
Lastly, I'll mention that the extra-deep OPLS F- potential leads to
really screwy behavior - if you run a 100 mm box of NaF<br>
in tip3p water (298K 1Atm) you'll find almost all of them form really,
really long lived contact ion pairs (half-life &gt; 10ns).<br>
This is a full two orders of magnitude longer than the the transient and
relatively rare&nbsp; NaCl pairs (aqvist) you<br>
would find at the same concentration . . .<br><br>
<br>
Cheers!<br><br>
Alan Chen<br><br>
Graduate Student<br>
Pappu Lab<br>
Washington University in St. Louis<br><br>
<br><br>
<br>
<blockquote type=cite class=cite cite=""><pre>Erik and all,

&gt; No, I don't think so, or at least it's not our typo :-)
&gt;
&gt; It is quite common that ions end up with strange radii when you
&gt; parametrize for free energy of solvation, and in addition those
&gt; parameters are somewhat old (1984). Here's the extract from Bill
&gt; Jorgensens file that we translated them automatically from:
&gt;
&gt;
&gt; 400 09 F-&nbsp;&nbsp; -1.0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
2.73295&nbsp;&nbsp; 0.72000&nbsp;&nbsp; F-&nbsp; JACS 106, 903 (1984)
&gt; 401 17 ??&nbsp;&nbsp; -1.0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
4.41724&nbsp;&nbsp; 0.11779&nbsp;&nbsp; Cl- JACS 106, 903 (1984)
&gt; 402 35 ??&nbsp;&nbsp; -1.0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
4.62376&nbsp;&nbsp; 0.09000&nbsp;&nbsp; Br- JACS 107, 7793(1985)
&gt; 404 03 Li&nbsp;&nbsp;&nbsp; 1.0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
1.25992&nbsp;&nbsp; 6.25000&nbsp;&nbsp; Li+ JACS 106, 903 (1984)
&gt; 405 11 Na&nbsp;&nbsp;&nbsp; 1.0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
1.89744&nbsp;&nbsp; 1.60714&nbsp;&nbsp; Na+ JACS 106, 903 (1984)

OK, that is weird. If I plot the Cl potential energies from AMBER and
OPLS, they are almost identical. But if I plot F from OPLS versus F
from AMBER, the well is roughly, well, 10x as deep in OPLS as AMBER...

I guess you're off the hook in terms of it being your typo... But I'm
still not convinced it's right. Ugh. My understanding is that the LJ
parameters from AMBER and OPLS have roughly similar origins, so it
doesn't make a lot of sense to me to just say that AMBER is wrong --
especially when Cl and the other ion parameters are pretty consistent.

I suppose I will have to dig into the literature. Anyway, thanks for
the response.

David</pre><font face="Courier New, Courier"></font></blockquote></body>
</html>