Search

Anion hydration is complicated by H-bond donation between neighboring water molecules in addition to H-bond donation to the anion. This situation can lead to competing structures for chemically simple clusters like (H$_2$O)$_n$Cl$^-$ and to anharmonic vibrational motions. Quasi-chemical theory builds from electronic structure treatment of isolated ion-water clusters, partitions the hydration free energy into inner-shell and outer-shell contributions, and provides a general statistical mechanical framework to study complications of anion hydration. The present study exploits dynamics calculations on isolated (H$_2$O)$_n$Cl$^-$ clusters to account for anharmonicity, utilizing ADMP (atom-centered basis sets and density-matrix propagation) tools. Comparing singly hydrated F$^-$ and Cl$^-$ clusters, classic OH-bond donation to the anion occurs for F$^-$, while Cl$^-$ clusters exhibit more flexible but dipole-dominated interactions between ligand and ion. The predicted Cl$^-$ -- F$^-$ hydration free energy difference agrees well with experiment, a significant theoretical step for addressing issues like Hofmeister ranking and selectivity in ion channels., Comment: 17 pages, including Supporting Information, 10 figures

Accurate predictions of the hydration free energy for anions typically have been more challenging than for cations. Hydrogen bond donation to the anion in hydrated clusters such as $\mathrm{F(H_2O)}_n{}^-$ can lead to delicate structures. Consequently, the energy landscape contains many local minima, even for small clusters, and these minima present a challenge for computational optimization. Utilization of cluster experimental results for the free energies of gas-phase clusters shows that, even though anharmonic effects are interesting, they need not be troublesome magnitudes for careful applications of quasi-chemical theory to ion hydration. Energy-optimized cluster structures for anions can leave the central ion highly exposed and application of implicit solvation models to these structures can incur more serious errors than for metal cations. Utilizing cluster structures sampled from ab initio molecular dynamics simulations substantially fixes those issues., Comment: 6 pages, 5 figures

Laying a basis for molecularly specific theory for the mobilities of ions in solutions of practical interest, we report a broad survey of velocity autocorrelation functions (VACFs) of Li$^+$ and PF$_6{}^-$ ions in water, ethylene carbonate, propylene carbonate, and acetonitrile solutions. We extract the memory function, $\gamma(t)$, which characterizes the random forces governing the mobilities of ions. We provide comparisons, controlling for electrolyte concentration and ion-pairing, for van~der~Waals attractive interactions and solvent molecular characteristics. For the heavier ion (PF$_6{}^-$), velocity relaxations are all similar: negative tail relaxations for the VACF and a clear second relaxation for $\gamma\left(t\right)$, observed previously also for other molecular ions and with \emph{n}-pentanol as solvent. For the light Li$^+$ ion, short time-scale oscillatory behavior masks simple, longer time-scale relaxation of $\gamma\left(t\right)$. But the corresponding analysis of the \emph{solventberg} Li$^+\left(\mathrm{H}_2\mathrm{O}\right)_4$ does conform to the standard picture set by all the PF$_6{}^-$ results., Comment: 10 pages, 8 figures, minor revisions to improve presentation. Now in press \emph{J. Chem. Phys.}

Extensive all-atom molecular dynamics calculations on the water-squalane interface for nine different loadings with sorbitan monooleate (SPAN80), at $T=300$K, are analyzed for the surface tension equation of state, desorption free energy profiles as they depend on loading, and to evaluate escape times for absorbed SPAN80 into the bulk phases. These results suggest that loading only weakly affects accommodation of a SPAN80 molecule by this squalane-water interface. Specifically, the surface tension equation of state is simple through the range of high tension to high loading studied, and the desorption free energy profiles are weakly dependent on loading here. The perpendicular motion of the centroid of the SPAN80 head-group ring is well-described by a diffusional model near the minimum of the desorption free energy profile. Lateral diffusional motion is weakly dependent on loading. Escape times evaluated on the basis of a diffusional model and the desorption free energies are $7\times 10^{-2}$ s (into the squalane) and $3\times 10^2$ h (into the water). The latter value is consistent with irreversible absorption observed by related experimental work., Comment: 6 pages, 8 figures

The role of solute attractive forces on hydrophobic interactions is studied by coordinated development of theory and simulation results for Ar atoms in water. We present a concise derivation of the local molecular field (LMF) theory for the effects of solute attractive forces on hydrophobic interactions, a derivation that clarifies the close relation of LMF theory to the EXP approximation applied to this problem long ago. The simulation results show that change from purely repulsive atomic solute interactions to include realistic attractive interactions \emph{diminishes} the strength of hydrophobic bonds. For the Ar-Ar rdfs considered pointwise, the numerical results for the effects of solute attractive forces on hydrophobic interactions are of opposite sign and larger in magnitude than predicted by LMF theory. That comparison is discussed from the point of view of quasi-chemical theory, and it is suggested that the first reason for this difference is the incomplete evaluation within LMF theory of the hydration energy of the Ar pair. With a recent suggestion for the system-size extrapolation of the required correlation function integrals, the Ar-Ar rdfs permit evaluation of osmotic second virial coefficients $B_2$. Those $B_2$ also show that incorporation of attractive interactions leads to more positive (repulsive) values. With attractive interactions in play, $B_2$ can change from positive to negative values with increasing temperatures. This is consistent with the historical work of Watanabe, \emph{et al.,} that $B_2 \approx 0$ for intermediate cases. In all cases here, $B_2$ becomes more attractive with increasing temperature., Comment: 9 pages, 7 figures

The dependence on volume fraction $\varphi$ of the Flory-Huggins $\chi_{\mathrm{wp}}\left(\varphi\right)$ describing the free energy of mixing of polymers in water is obtained by exploiting the connection of $\chi_{\mathrm{wp}}\left(\varphi\right)$ to the chemical potential of the water, for which quasi-chemical theory is satisfactory. We test this theoretical approach with simulation data for aqueous solutions of capped PEO oligomers. For CH$_3$(CH$_2$-O-CH$_2$)$_m$CH$_3$ ($m$=11), $\chi_{\mathrm{wp}}\left(\varphi\right)$ depends strongly on $\varphi$, consistent with experiment. These results identify coexisting water-rich and water-poor solutions at $T$ = 300K and $p$ = 1atm. Direct observation of the coexistence of these two solutions on simulation time scales supports that prediction for the system studied. This approach directly provides the osmotic pressures. The osmotic second virial coefficient for these chains is positive, reflecting repulsive interactions between the chains in the water, a good solvent for these chains., Comment: 4 pages, 4 figures

Molecular dynamics simulations of water with both multi-Kr and single Kr atomic solutes are carried out to implement quasi-chemical theory evaluation of the hydration free energy of Kr(aq). This approach obtains free energy differences reflecting Kr-Kr interactions at higher concentrations. Those differences are negative changes in hydration free energies with increasing concentrations at constant pressure. The changes are due to a slight reduction of packing contributions in the higher concentration case. The observed Kr-Kr distributions, analyzed with the extrapolation procedure of Kr\"{u}ger, \emph{et al.}, yield a modestly attractive osmotic second virial coefficient, $B_2\approx -60~\mathrm{cm}^3$/mol. The thermodynamic analysis interconnecting these two approaches shows that they are closely consistent with each other, providing support for both., Comment: 6 pages, 7 figures. Revision follows the extrapolation procedure of Refs. 33 and 34 which works nicely. The thermodynamic results are now clearly consistent. The $k \rightarrow 0$ extrapolation of the Fourier transform was not was satisfactory

Parallel-tempering MD results for a CH$_3$(CH$_2$-O-CH$_2$)$_m$CH$_3$ chain in water are exploited as a data-base for analysis of collective structural characteristics of the PEO globule with a goal of defining models permitting statistical thermodynamic analysis of dispersants of Corexit type. The chain structure factor, relevant to neutron scattering from a deuterated chain in neutral water, is considered specifically. The traditional continuum-Gaussian structure factor is inconsistent with the simple $k \rightarrow \infty$ behavior, but we consider a discrete-Gaussian model that does achieve that consistency. Shifting-and-scaling the discrete-Gaussian model helps to identify the low-$k$ to high-$k$ transition near $k \approx 2\pi/0.6 \mathrm{nm}$ when an empirically matched number of Gaussian links is about one-third of the total number of effective-atom sites. This short distance-scale boundary of 0.6 nm is directly verified with the $r$-space distributions, and this distance is thus identified with a natural size for coarsened monomers. The probability distribution of $R_g{}^2$ is compared with the classic predictions for both Gaussian model and freely-jointed chains. $\left\langle R_g{}^2(j)\right\rangle$, the contribution of the $j$-th chain segment to $\left\langle R_g{}^2\right\rangle$, depends on contour index about as expected for Gaussian chains despite significant quantitative discrepancies which express the swelling of these chains in water. Monomers central to the chain contour occupy the center of the chain globule. The density profiles of chain segments relative to their center of mass can show distinctive density structuring for smaller chains due close proximity of central elements to the globule center. But that density structuring washes-out for longer chains where many chain elements additively contribute to the density profiles., Comment: 5 pages, 6 figures

This work tests the statistical mechanical theory of hydrophobic interactions, isolates consequences of excluded volume interactions, and obtains B2 for those purposes. Cavity methods that are particularly appropriate for study of hydrophobic interactions between atomic-size hard spheres in liquid water are developed and applied to test aspects of the Pratt-Chandler (PC) theory that have not been tested. Contact hydrophobic interactions between Ar-size hard-spheres in water are significantly more attractive than predicted by the PC theory. The corresponding results for the osmotic second virial coefficient are attractive (B2 <0), and more attractive with increasing temperature (Delta B2/Delta T < 0) in the temperature range 300K < T < 360K. This information has not been available previously, but is essential for development of the molecular-scale statistical mechanical theory of hydrophobic interactions, particularly for better definition of the role of attractive intermolecular interactions associated with the solutes., Comment: five pages, four figures