Yes, that&#39;s exactly the article I had in mind!<br><br>There are also people studying stuff like dissociative hydration of, e.g., silica, they have some water models that can dissociate when needed, as well...<br><br><a href="http://pubs.acs.org/cgi-bin/abstract.cgi/nalefd/2003/3/i10/abs/nl034480e.html">http://pubs.acs.org/cgi-bin/abstract.cgi/nalefd/2003/3/i10/abs/nl034480e.html</a><br>
<a href="http://www.rsc.org/Publishing/Journals/DT/article.asp?doi=b516258k">http://www.rsc.org/Publishing/Journals/DT/article.asp?doi=b516258k</a><br><br><div class="gmail_quote">2008/4/28 David van der Spoel &lt;<a href="mailto:spoel@xray.bmc.uu.se">spoel@xray.bmc.uu.se</a>&gt;:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Vasilii Artyukhov wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<br>
2008/4/28 Mark Abraham &lt;<a href="mailto:Mark.Abraham@anu.edu.au" target="_blank">Mark.Abraham@anu.edu.au</a> &lt;mailto:<a href="mailto:Mark.Abraham@anu.edu.au" target="_blank">Mark.Abraham@anu.edu.au</a>&gt;&gt;:<div>
<div></div><div class="Wj3C7c"><br>
<br>
 &nbsp; &nbsp; Trying to model hydronium is a pretty flawed thing to do, since the<br>
 &nbsp; &nbsp;hydrogens in real water are labile - so the molecule you&#39;d label as<br>
 &nbsp; &nbsp;&quot;hydronium&quot; moves much faster than any atoms do. <br>
<br>
I actually seem to remember reading something about some dissociative water models that could exchange protons between molecules, so that one could really have something like a pH in the simulation... doesn&#39;t really account for the quantum stuff like tunneling, though :)<br>

</div></div></blockquote>
<br>
<br>
@Article{ Mahadevan2007a,<br>
 &nbsp;author = &nbsp; &nbsp; &nbsp; &quot;T. S. Mahadevan and S. H. Garofalini&quot;,<br>
 &nbsp;title = &nbsp; &nbsp; &nbsp; &nbsp;&quot;Dissociative Water Potential for Molecular Dynamics<br>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;Simulations&quot;,<br>
 &nbsp;journal = &nbsp; &nbsp; &nbsp;&quot;J. Phys. Chem. B&quot;,<br>
 &nbsp;year = &nbsp; &nbsp; &nbsp; &nbsp; 2007,<br>
 &nbsp;url =<br>
<br>
&quot;<a href="http://pubs.acs.org/cgi-bin/abstract.cgi/jpcbfk/asap/abs/jp072530o.html" target="_blank">http://pubs.acs.org/cgi-bin/abstract.cgi/jpcbfk/asap/abs/jp072530o.html</a>&quot;,<br>
 &nbsp;doi = &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;&quot;10.1021/jp072530o&quot;,<br>
 &nbsp;abstract = &nbsp; &nbsp; &quot;A new interatomic potential for dissociative water<br>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;was developed for use in molecular dynamics<br>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;simulations. The simulations use a multibody<br>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;potential, with both pair and three-body terms, and<br>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;the Wolf summation method for the long-range Coulomb<br>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;interactions. A major feature in the potential is<br>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;the change in the short-range O-H repulsive<br>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;interaction as a function of temperature and/or<br>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;pressure in order to reproduce the<br>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;density-temperature curve between 273 K and 373 at 1<br>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;atm, as well as high-pressure data at various<br>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;temperatures. Using only the change in this one<br>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;parameter, the simulations also reproduce<br>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;room-temperature properties of water, such as the<br>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;structure, cohesive energy, diffusion constant, and<br>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;vibrational spectrum, as well as the liquid-vapor<br>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;coexistence curve. Although the water molecules<br>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;could dissociate, no dissociation is observed at<br>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;room temperature. However, behavior of the hydronium<br>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;ion was studied by introduction of an extra H+ into<br>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;a cluster of water molecules. Both Eigen and Zundel<br>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;configurations, as well as more complex<br>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;configurations, are observed in the migration of the<br>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;hydronium.&quot;<br>
}<br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<br>
<br>
------------------------------------------------------------------------<div class="Ih2E3d"><br>
<br>
_______________________________________________<br>
gmx-users mailing list &nbsp; &nbsp;<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<a href="http://www.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://www.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/search</a> before posting!<br>
Please don&#39;t post (un)subscribe requests to the list. Use the www interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>.<br>
Can&#39;t post? Read <a href="http://www.gromacs.org/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a><br>
</div></blockquote><font color="#888888">
<br>
<br>
-- <br>
David van der Spoel, Ph.D.<br>
Molec. Biophys. group, Dept. of Cell &amp; Molec. Biol., Uppsala University.<br>
Box 596, 75124 Uppsala, Sweden. Phone: &nbsp;+46184714205. Fax: +4618511755.<br>
<a href="mailto:spoel@xray.bmc.uu.se" target="_blank">spoel@xray.bmc.uu.se</a> &nbsp; &nbsp;<a href="mailto:spoel@gromacs.org" target="_blank">spoel@gromacs.org</a> &nbsp; <a href="http://folding.bmc.uu.se" target="_blank">http://folding.bmc.uu.se</a></font><div>
<div></div><div class="Wj3C7c"><br>
_______________________________________________<br>
gmx-users mailing list &nbsp; &nbsp;<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<a href="http://www.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://www.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/search</a> before posting!<br>
Please don&#39;t post (un)subscribe requests to the list. Use the www interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>.<br>
Can&#39;t post? Read <a href="http://www.gromacs.org/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a><br>
</div></div></blockquote></div><br>