Neutron diffraction studies of electrolytes in null water: a direct determination of the first hydration zone of ions

High-Dimensional Parameter Search Method to Determine Force Field Mixing Terms in Molecular Simulations

Molecular dynamics (MD) force fields for lipids and ions are typically developed independently of one another. In simulations consisting of both lipids and ions, lipid-ion interaction energies are estimated using a predefined set of mixing rules for Lennard-Jones (LJ) interactions. This, however, does not guarantee their reliability. In fact, compared to the quantum mechanical reference data, Lorentz-Berthelot mixing rules substantially underestimate the binding energies of Na⁺ ions with small-molecule analogues of lipid headgroups, yielding errors on the order of 80 and 130 kJ/mol, respectively, for methyl acetate and diethyl phosphate. Previously, errors associated with mixing force fields have been reduced using approaches such as "NB-fix" in which LJ interactions are computed using explicit cross terms rather than those from mixing rules. Building on this idea, we derive explicit lipid-ion cross terms that also may implicitly include many-body cooperativity effects. Additionally, to account for the interdependency between cross terms, we optimize all cross terms simultaneously by performing high-dimensional searches using our ParOpt software. The cross terms we obtain reduce the errors due to mixing rules to below 10 kJ/mol. MD simulation of the lipid bilayer conducted using these optimized cross terms resolves the structural discrepancies between our previous simulations and small-angle X-ray and neutron scattering experiments. These results demonstrate that simulations of lipid bilayers with ions that are accurate up to structural data from scattering experiments can be performed without explicit polarization terms. However, it is worth noting that such NB-fix cross terms are not based on any physical principle; a polarizable lipid model would be more realistic and is still desired. Our approach is generic and can be applied to improve the accuracies of simulations employing mixed force fields.

Structural Analysis of Molecular Materials Using the Pair Distribution Function

This is a review of atomic pair distribution function (PDF) analysis as applied to the study of molecular materials. The PDF method is a powerful approach to study short- and intermediate-range order in materials on the nanoscale. It may be obtained from total scattering measurements using X-rays, neutrons, or electrons, and it provides structural details when defects, disorder, or structural ambiguities obscure their elucidation directly in reciprocal space. While its uses in the study of inorganic crystals, glasses, and nanomaterials have been recently highlighted, significant progress has also been made in its application to molecular materials such as carbons, pharmaceuticals, polymers, liquids, coordination compounds, composites, and more. Here, an overview of applications toward a wide variety of molecular compounds (organic and inorganic) and systems with molecular components is presented. We then present pedagogical descriptions and tips for further implementation. Successful utilization of the method requires an interdisciplinary consolidation of material preparation, high quality scattering experimentation, data processing, model formulation, and attentive scrutiny of the results. It is hoped that this article will provide a useful reference to practitioners for PDF applications in a wide realm of molecular sciences, and help new practitioners to get started with this technique.

Ionic radii of hydrated sodium cation from QTAIM

The sodium cation is ubiquitous in aqueous chemistry and biological systems. Yet, in spite of numerous studies, the (average) distance between the sodium cation and its water ligands, and the corresponding ionic radii, are still controversial. Recent experimental values in solution are notably smaller than those from previous X-ray studies and ab initio molecular dynamics. Here we adopt a "bottom-up" approach of obtaining these distances from quantum chemistry calculations [full MP2 with the 6-31++G(d,p) and cc-pVTZ basis-sets] of gas-phase Na⁺(H₂O)_n clusters, as a function of the sodium coordination number (CN = 2-6). The bulk limit is obtained by the polarizable continuum model, which acts to increase the interatomic distances at small CN, but has a diminishing effect as the CN increases. This extends the CN dependence of the sodium-water distances from crystal structures (CN = 4-12) to lower CN values, revealing a switch between two power laws, having a small exponent at small CNs and a larger one at large CNs. We utilize Bader's theory of atoms in molecules to bisect the Na⁺-O distances into Na⁺ and water radii. Contrary to common wisdom, the water radius is not constant, decreasing even more than that of Na⁺ as the CN decreases. We also find that the electron density at the bond critical point increases exponentially as the sodium radius decreases.

Evidence for Reduced Hydrogen-Bond Cooperativity in Ionic Solvation Shells from Isotope-Dependent Dielectric Relaxation

We find that the reduction in dielectric response (depolarization) of water caused by solvated ions is different for H_{2}O and D_{2}O. This isotope dependence allows us to reliably determine the kinetic contribution to the depolarization, which is found to be significantly smaller than predicted by existing theory. The discrepancy can be explained from a reduced hydrogen-bond cooperativity in the solvation shell: we obtain quantitative agreement between theory and experiment by reducing the Kirkwood correlation factor of the solvating water from 2.7 (the bulk value) to ∼1.6 for NaCl and ∼1 (corresponding to completely uncorrelated motion of water molecules) for CsCl.

Unusual Influence of Fluorinated Anions on the Stretching Vibrations of Liquid Water

Infrared (IR) spectroscopy is a commonly used and invaluable tool in the studies of solvation phenomena in aqueous solutions. Concurrently, ab initio molecular dynamics (AIMD) simulations deliver the solvation shell picture at a molecular detail level and allow for a consistent decomposition of the theoretical IR spectrum into underlying spatial correlations. Here, we demonstrate how the novel spectral decomposition techniques can extract important information from the computed IR spectra of aqueous solutions of BF₄^- and PF₆^-, interesting weakly coordinating anions that have been known for a long time to alter the IR spectrum of water in an unusual manner. The distance-dependent spectra of both ions are analyzed using the spectral similarity method that provides a quantitative picture of both the spectrum of the solute-affected solvent and the number of solvent molecules thus altered. We find, in accordance with previous experiments, a considerable blue shift of the ν_OH stretching band of liquid water by 264 cm^-1 for BF₄^- and 306 cm^-1 for PF₆^-, with the affected numbers being 3.7 and 4.2, respectively. Considering also the additional information on solute-solvent dipolar couplings delivered by radially and spatially resolved IR spectra, the computational IR spectroscopy based on AIMD simulations is shown to be a viable predictive tool with strong interpretative power.

The conformation of a ribose derivative in aqueous solution: a neutron-scattering and molecular dynamics study

The structure of aqueous solutions of methyl β-D-ribofuranoside was investigated by coupling molecular dynamics (MD) simulations and neutron scattering measurements with isotopic substitution. Using a sample of the sugar isotopically-labeled at a single unique position, neutron scattering structure factors and radial distribution functions can be compared with MD simulations constrained to different conformations to determine which conformer best fits the experimental results. Three different simulations were performed with the methyl ether group of the sugar unconstrained and constrained in each of its staggered orientations. The results of the unconstrained simulation showed that the methyl ester group occupied predominantly the 300° position, which is in agreement with the diffraction experimental results. This result suggests that the molecular mechanics force field used in the simulation adequately describes the conformation of the 1-methyl ether group in the methyl β-D-ribofuranoside.

The Radial Distribution Functions of Water as Derived from Radiation Total Scattering Experiments: Is There Anything We Can Say for Sure?

The present paper reviews the investigation of ambient water structure and focusses in particular on the determination of the radial distribution functions of water from total experimental radiation scattering experiments. A novel method for removing the inelastic scattering from neutron data is introduced, and the effect of Compton scattering on X-ray data is discussed. In addition the extent to which quantum effects can be discerned between heavy and light water is analysed against these more recent data. It is concluded that, with the help of modern data analysis and computer simulation tools to interrogate the scattering data, a considerable degree of consistency can be obtained between recent and past scattering experiments on water. That consistency also gives a realistic estimate of the likely uncertainties in the extracted radial distribution functions, as well as offering a benchmark against which future experiments can be judged.

UNDERSTANDING STRUCTURAL EFFECTS OF MULTIPOLE MOMENTS ON AQUEOUS SOLVATION OF IONS USING THE SOFT-STICKY DIPOLE-QUADRUPOLE-OCTUPOLE WATER MODEL

The effects of water multipole moments on the aqueous solvation of ions were determined in Monte Carlo simulations using soft-sticky dipole-quadrupole-octupole (SSDQO) water. Water molecules formed linear hydrogen bonds to Cl(-) using the new SSDQO1 parameters, similar to multi-site models. However, the dipole vector was tilted rather than parallel to the oxygen-Na(+) internuclear vector as in most multi-site model, while experiment and ab initio molecular dynamics simulations generally indicate a range of values between tilted and parallel. By varying the multipoles in SSDQO, the octupole was found to determine the orientation around Na(+). Moreover, analysis of the multipoles of more conventional models is predictive of their performance as solvents.

Ion selectivity from local configurations of ligands in solutions and ion channels

Probabilities of numbers of ligands proximal to an ion lead to simple, general formulae for the free energy of ion selectivity between different media. That free energy does not depend on the definition of an inner shell for ligand-counting, but other quantities of mechanistic interest do. If analysis is restricted to a specific coordination number, then two distinct probabilities are required to obtain the free energy in addition. The normalizations of those distributions produce partition function formulae for the free energy. Quasi-chemical theory introduces concepts of chemical equilibrium, then seeks the probability that is simplest to estimate, that of the most probable coordination number. Quasi-chemical theory establishes the utility of distributions of ligand-number, and sharpens our understanding of quasi-chemical calculations based on electronic structure methods. This development identifies contributions with clear physical interpretations, and shows that evaluation of those contributions can establish a mechanistic understanding of the selectivity in ion channels.