Dynamics of ions in a water drop using the AMOEBA polarizable force field

Variation of Surface Propensity of Halides with Droplet Size and Temperature: The Planar Interface Limit

The radial number density profiles of halide and alkali ions in aqueous clusters with equimolar radius ≲1.4 nm, which correspond to ≲255 H2O molecules, have been extensively studied by computations. However, the surface abundance of Cl-, Br-, and I- relative to the bulk interior in these smaller clusters may not be representative of the larger systems. Indeed, here we show that the larger the cluster is, the lower the relative surface abundance of chaotropic halides is. In droplets with an equimolar radius of ≈2.45 nm, which corresponds to ≈2000 H2O molecules, the polarizable halides show a clear number density maximum in the droplet's bulk-like interior. A similar pattern is observed in simulations of the aqueous planar interface with halide salts at room temperature. At elevated temperature the surface propensity of Cl- decreases gradually, while that of I- is partially preserved. The change in the chaotropic halide location at higher temperatures than the room temperature may considerably affect photochemical reactivity in atmospheric aerosols, vapor-liquid nucleation and growth mechanisms, and salt crystallization via solvent evaporation. We argue that the commonly used approach of nullifying parameters in a force field in order to find the factors that determine the ion location does not provide transferable insight into other force fields.

Modeling Gd3+ Complexes for Molecular Dynamics Simulations: Toward a Rational Optimization of MRI Contrast Agents

The correct parametrization of lanthanide complexes is of the utmost importance for their characterization using computational tools such as molecular dynamics simulations. This allows the optimization of their properties for a wide range of applications, including medical imaging. Here we present a systematic study to establish the best strategies for the correct parametrization of lanthanide complexes using [Gd(DOTA)]^- as a reference, which is used as a contrast agent in MRI. We chose the bonded model to parametrize the lanthanide complexes, which is especially important when considering the study of the complex as a whole (e.g., for the study of the dynamics of its interaction with proteins or membranes). We followed two strategies: a so-called heuristic approach employing strategies already published by other authors and another based on the more recent MCPB.py tool. Adjustment of the Lennard-Jones parameters of the metal was required. The final topologies obtained with both strategies were able to reproduce the experimental ion to oxygen distance, vibrational frequencies, and other structural properties. We report a new strategy to adjust the Lennard-Jones parameters of the metal ion in order to capture dynamic properties such as the residence time of the capping water (τ_m). For the first time, the correct assessment of the τ_m value for Gd-based complexes was possible by recording the dissociative events over up to 10 μs all-atom simulations. The MCPB.py tool allowed the accurate parametrization of [Gd(DOTA)]^- in a simpler procedure, and in this case, the dynamics of the water molecules in the outer hydration sphere was also characterized. This sphere was divided into the first hydration layer, an intermediate region, and an outer hydration layer, with a residence time of 18, 10 and 19 ps, respectively, independent of the nonbonded parameters chosen for Gd³⁺. The Lennard-Jones parameters of Gd³⁺ obtained here for [Gd(DOTA)]^- may be used with similarly structured gadolinium MRI contrast agents. This allows the use of molecular dynamics simulations to characterize and optimize the contrast agent properties. The characterization of their interaction with membranes and proteins will permit the design of new targeted contrast agents with improved pharmacokinetics.

Manifolds of low energy structures for a magic number of hydrated sulfate: SO42-(H2O)24

Low energy structures of SO₄^2-(H₂O)₂₄ have been obtained using a combination of classical molecular dynamics simulations and refinement of structures and energies by quantum chemical calculations. Extensive exploration of the potential energy surface led to a number of low-energy structures, confirmed by accurate calibration calculations. An overall analysis of this large set was made after devising appropriate structural descriptors such as the numbers of cycles and their combinations. Low energy structures bear common motifs, the most prominent being fused cycles involving alternatively four and six water molecules. The latter adopt specific conformations which ensure the appropriate surface curvature to form a closed cage without dangling O-H bonds and at the same time provide 12-coordination of the sulfate ion. A prominent feature to take into account is isomerism via inversion of hydrogen bond orientations along cycles. This generates large families of ca. 100 isomers for this cluster size, spanning energy windows of 10-30 kJ mol^-1. This relatively ignored isomerism must be taken into account to identify reliably the lowest energy minima. The overall picture is that the magic number cluster SO₄^2-(H₂O)₂₄ does not correspond to formation of a single, remarkable structure, but rather to a manifold of structural families with similar stabilities. Extensive calculations on isomerization mechanisms within a family indicate that large barriers are associated to direct inversion of hydrogen bond networks. Possible implications of these results for magic number clusters of other anions are discussed.

Low Density Interior in Supercooled Aqueous Nanodroplets Expels Ions to the Subsurface

The interaction between water and ions within droplets plays a key role in the chemical reactivity of atmospheric and man-made aerosols. Here we report direct computational evidence that in supercooled aqueous nanodroplets a lower density core of tetrahedrally coordinated water expels the cosmotropic ions to the denser and more disordered subsurface. In contrast, at room temperature, depending on the nature of the ion, the radial distribution in the droplet core is nearly uniform or elevated toward the center. We analyze the spatial distribution of a single ion in terms of a reference electrostatic model. The energy of the system in the analytical model is expressed as the sum of the electrostatic and surface energy of a deformable droplet. The model predicts that the ion is subject to a harmonic potential centered at the droplet's center of mass. We name this effect "electrostatic confinement". The model's predictions are consistent with the simulation findings for a single ion at room temperature but not at supercooling. We anticipate this study to be the starting point for investigating the structure of supercooled (electro)sprayed droplets that are used to preserve the conformations of macromolecules originating from the bulk solution.

Two-Dimensional Infrared Spectroscopy of Aqueous Solutions of Metal Nitrates: Slowdown of Spectral Diffusion in the Presence of Divalent Cations

The hydrogen-bonded network of water can be affected both structurally and dynamically by the presence of ions. In the present study, we have considered three aqueous solutions of metal nitrates to investigate the effects of divalent cations (Mg²⁺ and Ca²⁺), compared to that of monovalent Na⁺ ions, on hydrogen-bond fluctuations and vibrational spectral diffusion through calculations of linear and two-dimensional infrared spectra of these solutions at room temperature. We have employed the methods of molecular dynamics simulations using effective polarizable models of ions combined with quantum mechanical calculations of transition variables and statistical mechanical calculations of spectral response functions of vibrational spectroscopy. Divalent cations are found to have much stronger and longer-ranged effects on the structure and dynamics of the hydrogen-bonded network than that induced by the monovalent sodium ions. The blue shifts in the calculated linear spectra are found to follow the Hofmeister trend for the cations. The 2D-IR spectral lineshape and intensity corresponding to three-pulse echo peak shift (3PEPS) experiments are calculated. The timescales of these nonlinear spectral responses and also frequency-time correlations show significant slowing down of spectral diffusion for solutions containing divalent Mg²⁺ and Ca²⁺ ions compared to the corresponding dynamics of the solution containing Na⁺ ions. Unlike NaNO₃ solution, the relaxation of frequency and dipole orientational fluctuations of anion-bound water in Mg(NO₃)₂ and Ca(NO₃)₂ solutions are found to be somewhat slower than bulk water, which can be attributed to the presence of divalent cations whose effects go beyond their first solvation shells. This is also seen in the dynamics of bulk water in these solutions which is found to be notably slower for the solutions containing divalent cations than that in the NaNO₃ solution. Unlike Mg²⁺ and Ca²⁺ ions, no specific cationic effect is observed for the Na⁺ ions.

Assessing cluster models of solvation for the description of vibrational circular dichroism spectra: synergy between static and dynamic approaches

Solvation effects are essential for defining the shape of vibrational circular dichroism (VCD) spectra.

Strategy for Modeling the Infrared Spectra of Ion-Containing Water Drops

Hydrated ions are ubiquitous in environmental and biological media. Understanding the perturbation exerted by an ion on the water hydrogen bond network is possible in the nanodrop regime by recording vibrational spectra in the O-H bond stretching region. This has been achieved experimentally in recent years by forming gaseous ions containing tens to hundreds of water molecules and recording their infrared photodissociation spectra. In this paper, we demonstrate the capabilities of a modeling strategy based on an extension of the AMOEBA polarizable force field to implement water atomic charge fluctuations along with those of intramolecular structure along the dynamics. This supplementary flexibility of nonbonded interactions improves the description of the hydrogen bond network and, therefore, the spectroscopic response. Finite temperature IR spectra are obtained from molecular dynamics simulations by computing the Fourier transform of the dipole moment autocorrelation function. Simulations of 1-2 ns are required for extensive sampling in order to reproduce the experimental spectra. Furthermore, bands are assigned with the driven molecular dynamics approach. This method package is shown to compare successfully with experimental spectra for 11 ions in water drops containing 36-100 water molecules. In particular, band frequency shifts of the free O-H stretching modes at the cluster surface are well reproduced as a function of both ion charge and drop size.