Nanometer-Scale Correlations in Aqueous Salt Solutions

Detecting underscreening and generalized Kirkwood transitions in aqueous electrolytes

We establish the connection between the measured small angle x-ray scattering signal and the charge-charge correlations underlying Kirkwood transitions (KTs) in 1:1, 2:1, and 3:1 aqueous electrolytes. These measurements allow us to obtain underscreening lengths for bulk electrolytes independently verified by theory and simulations. Furthermore, we generalize the concept of KTs beyond those theoretically predicted for 1:1 electrolytes, which involves the inverse screening length, a0, and the inverse periodicity length, Q0. Above the KTs, we find a universal scaling of a0∝c-ζ/3 and Q0 ∝ c1/3 for the studied electrolyte solutions, where ζ is the ionic strength factor.

Resolving intermediates during the growth of aluminum deuteroxide (Hydroxide) polymorphs in high chemical potential solutions

Aluminum hydroxide polymorphs are of widespread importance yet their kinetics of nucleation and growth remain beyond the reach of current models. Here we attempt to unveil the reaction processes underlying the polymorphs formation at high chemical potential. We examine their formation in-situ from supersaturated alkaline sodium aluminate solutions using deuteration and time-resolved neutron pair distribution function analyses, which indicate the formation of individual Al(OD)3 layers as an intermediate particle phase. These layers ultimately stack to form gibbsite- or bayerite-like layered heterostructures. Ex-situ characterization of the recovered precipitates using 27Al magic angle spinning nuclear magnetic resonance spectroscopy, Raman, X-ray diffraction, and scanning electron microscopy, suggests the presence of additional intermediate states, an amorphous compound bearing both tetrahededrally- and penta-coordinated Al3+. These observations reveal the complex pathways to form Al(OD)3 monolayers via either transient pentacoordinate species or amorphous-to-ordered transitions. The subsequent crystallization of admixed gibbsite/bayerite is followed by an Al(OD)3 monolayer attachment process.

Dynamic Ion Correlations and Ion-Pair Lifetimes in Aqueous Alkali Metal Chloride Electrolytes

Electrolytes are central to many technological applications, as well as life itself. The behavior and properties of electrolytes are often described in terms of ion pairs, whereby ions associate as either contact ion pairs (in which ions are "touching") solvent-separated ion pairs (in which ions' solvent shells overlap) or solvent-solvent-separated ion pairs (in which ions' solvent shells are distinct). However, this paradigm is generally restricted to statistically averaged descriptions of solution structure and ignores temporal behavior. Here we elucidate the time-resolved dynamics of these ion-ion interactions in aqueous metal chloride electrolytes using the partial van Hove correlation function, based on polarizable molecular dynamics simulations. Our results show that the existence and persistence of ion pairs in aqueous metal chloride electrolytes should not be assumed a priori, but in fact are ion specific features of the solution with lifetimes on subpicosecond time scales.

Dynamic Second Harmonic Imaging of Proton Translocation Through Water Needles in Lipid Membranes

Proton translocation through lipid membranes is a fundamental process in the field of biology. Several theoretical models have been developed and presented over the years to explain the phenomenon, yet the exact mechanism is still not well understood. Here, we show that proton translocation is directly related to membrane potential fluctuations. Using high-throughput wide-field second harmonic (SH) microscopy, we report apparently universal transmembrane potential fluctuations in lipid membrane systems. Molecular simulations and free energy calculations suggest that H+ permeation proceeds predominantly across a thin, membrane-spanning water needle and that the transient transmembrane potential drives H+ ions across the water needle. This mechanism differs from the transport of other cations that require completely open pores for transport and follows naturally from the well-known Grotthuss mechanism for proton transport in bulk water. Furthermore, SH imaging and conductivity measurements reveal that the rate of proton transport depends on the structure of the hydrophobic core of bilayer membranes.

Reexpansion of charged nanoparticle assemblies in concentrated electrolytes

Electrostatic forces in solutions are highly relevant to a variety of fields, ranging from electrochemical energy storage to biology. However, their manifestation in concentrated electrolytes is not fully understood, as exemplified by counterintuitive observations of colloidal stability and long-ranged repulsions in molten salts. Highly charged biomolecules, such as DNA, respond sensitively to ions in dilute solutions. Here, we use non-base-pairing DNA-coated nanoparticles (DNA-NP) to analyze electrostatic interactions in concentrated salt solutions. Despite their negative charge, these conjugates form colloidal crystals in solutions of sufficient divalent cation concentration. We utilize small-angle X-ray scattering (SAXS) to study such DNA-NP assemblies across the full accessible concentration ranges of aqueous CaCl2, MgCl2, and SrCl2 solutions. SAXS shows that the crystallinity and phases of the assembled structures vary with cation type. For all tested salts, the aggregates contract with added ions at low salinities and then begin expanding above a cation-dependent threshold salt concentration. Wide-angle X-ray scattering (WAXS) reveals enhanced positional correlations between ions in the solution at high salt concentrations. Complementary molecular dynamics simulations show that these ion-ion interactions reduce the favorability of dense ion configurations within the DNA brushes below that of the bulk solution. Measurements in solutions with lowered permittivity demonstrate a simultaneous increase in ion coupling and decrease in the concentration at which aggregate expansion begins, thus confirming the connection between these phenomena. Our work demonstrates that interactions between charged objects continue to evolve considerably into the high-concentration regime, where classical theories project electrostatics to be of negligible consequence.

Ionization Difference between Weak and Strong Electrolytes as Perturbed by Conductivity Spectra Analysis

While the static structure of aqueous electrolytes has been studied for decades, their dynamic microscopic structure remains unresolved yet critical in many areas. We report a comparative study of dc and ac (1 Hz to 20 GHz) conductivity data of weak and strong electrolytes, highlighting previously missing differences and similarities. Based on these results, we introduce into consideration the intrinsic short-lived ions of water, namely, excess protons (H₃O⁺) and proton holes (OH^-). We show that the model accounting for the neutralization of these ions by the species of electrolyte explains the conductivity of aqueous solutions in the concentration range 10^-7-10 M. Based on independent experimental data, we hypothesize that the aggregation of the species in weak electrolytes may determine the main difference between the conductivity of weak and strong electrolytes. Our results push forward the understanding of the dynamic structure of aqueous electrolyte solutions and are important to nanofluidic, biological, and electrochemical systems.

Ion-Induced Transient Potential Fluctuations Facilitate Pore Formation and Cation Transport through Lipid Membranes

Unassisted ion transport through lipid membranes plays a crucial role in many cell functions without which life would not be possible, yet the precise mechanism behind the process remains unknown due to its molecular complexity. Here, we demonstrate a direct link between membrane potential fluctuations and divalent ion transport. High-throughput wide-field non-resonant second harmonic (SH) microscopy of membrane water shows that membrane potential fluctuations are universally found in lipid bilayer systems. Molecular dynamics simulations reveal that such variations in membrane potential reduce the free energy cost of transient pore formation and increase the ion flux across an open pore. These transient pores can act as conduits for ion transport, which we SH image for a series of divalent cations (Cu2+, Ca2+, Ba2+, Mg2+) passing through giant unilamellar vesicle (GUV) membranes. Combining the experimental and computational results, we show that permeation through pores formed via an ion-induced electrostatic field is a viable mechanism for unassisted ion transport.

Coexistence of Transient Liquid Droplets and Amorphous Solid Particles in Nonclassical Crystallization of Cerium Oxalate

Crystallization from solution often occurs via "nonclassical" routes; that is, it involves transient, non-crystalline states like reactant-rich liquid droplets and amorphous particles. However, in mineral crystals, the well-defined thermodynamic character of liquid droplets and whether they convert─or not─into amorphous phases have remained unassessed. Here, by combining cryo-transmission electron microscopy and X-ray scattering down to a 250 ms reaction time, we unveil that crystallization of cerium oxalate involves a metastable chemical equilibrium between transient liquid droplets and solid amorphous particles: contrary to the usual expectation, reactant-rich droplets do not evolve into amorphous solids. Instead, at concentrations above 2.5 to 10 mmol L^-1, both amorphous and reactant-rich liquid phases coexist for several tens of seconds and their molar fractions remain constant and follow the lever rule in a multicomponent phase diagram. Such a metastable chemical equilibrium between solid and liquid precursors has been so far overlooked in multistep nucleation theories and highlights the interest of rationalizing phase transformations using multicomponent phase diagrams not only when designing and recycling rare earths materials but also more generally when describing nonclassical crystallization.

Electric fields and potentials in condensed phases

The electric fields and potentials inside and at the interface of matter are relevant to many branches of physics, chemistry, and biology. Accurate quantification of these fields and/or potentials is essential to control and exploit chemical and physical transformations. Before we understand the response of matter to external fields, it is first important to understand the intrinsic interior and interfacial fields and potentials, both classically and quantum mechanically, as well as how they are probed experimentally. Here we compare and contrast, beginning with the hydrogen atom in vacuum and ending with concentrated aqueous NaCl electrolyte, both classical and quantum mechanical electric potentials and fields. We make contact with experimental vibrational Stark, electrochemical, X-ray, and electron spectroscopic probes of these potentials and fields, outline relevant conceptual difficulties, and underscore the advantage of electron holography as a basis to better understand electrostatics in matter.

Toward a First-Principles Framework for Predicting Collective Properties of Electrolytes

Given the universal importance of electrolyte solutions, it is natural to expect that we have a nearly complete understanding of the fundamental properties of these solutions (e.g., the chemical potential) and that we can therefore explain, predict, and control the phenomena occurring in them. In fact, reality falls short of these expectations. But, recent advances in the simulation and modeling of electrolyte solutions indicate that it should soon be possible to make progress toward these goals. In this Account, we will discuss the use of first-principles interaction potentials based in quantum mechanics (QM) to enhance our understanding of electrolyte solutions. Specifically, we will focus on the use of quantum density functional theory (DFT) combined with molecular dynamics simulation (DFT-MD) as the foundation for our approach. The overarching concept is to understand and accurately reproduce the balance between local or short-ranged (SR) structural details and long-range (LR) correlations, allowing the prediction of the thermodynamics of both single ions in solution as well as the collective interactions characterized by activity/osmotic coefficients. In doing so, relevant collective motions and driving forces characterized by chemical potentials can be determined.In this Account, we will make the case that understanding electrolyte solutions requires a faithful QM representation of the SR nature of the ion-ion, ion-water, and water-water interactions. However, the number of molecules that is required for collective behavior makes the direct application of high-level QM methods that contain the best SR physics untenable, making methods that balance accuracy and efficiency a practical goal. Alternatives such as continuum solvent models (CSMs) and empirically based classical molecular dynamics have been extensively employed to resolve this problem but without yet overcoming the fundamental issue of SR accuracy. We will demonstrate that accurately describing the SR interaction is imperative for predicting both intrinsic properties, namely, at infinite dilution, and collective properties of electrolyte solutions.DFT has played an important role in our understanding of condensed phase systems, e.g., bulk liquid water, the air-water interface, ions in bulk, and at the air-water interface. This approach holds huge promise to provide benchmark calculations of electrolyte solution properties that will allow for the development and improvement of more efficient methods, as well as an enhanced understanding of fundamental phenomena. However, the standard protocol using the generalized gradient approximation with van der Waals (vdW) correction requires improvement in order to achieve a high level of quantitative accuracy. Simply simulating with higher level DFT functionals may not be the best route considering the significant computational cost. Alternative methods of incorporating information from higher levels of QM should be explored; e.g., using force matching techniques on small clusters, where high level benchmark calculations are possible, to develop ideal correction terms to the DFT functional is a promising possibility. We argue that DFT with statistical mechanics is becoming an increasingly useful framework enabling the prediction of collective electrolyte properties.