Electrolytes induce long-range orientational order and free energy changes in the H-bond network of bulk water

Kinetic trapping of nanoparticles by solvent-induced interactions

Theoretical analysis based on mean field theory indicates that solvent-induced interactions (i.e. structural forces due to the rearrangement of wetting solvent molecules) not considered in DLVO theory can induce the kinetic trapping of nanoparticles at finite nanoscale separations from a well-wetted surface, under a range of ubiquitous physicochemical conditions for inorganic nanoparticles of common materials (e.g., metal oxides) in water or simple molecular solvents. This work proposes a simple analytical model that is applicable to arbitrary materials and simple solvents to determine the conditions for direct particle-surface contact or kinetic trapping at finite separations, by using experimentally measurable properties (e.g., Hamaker constants, interfacial free energies, and nanoparticle size) as input parameters. Analytical predictions of the proposed model are verified by molecular dynamics simulations and numerical solution of the Smoluchowski diffusion equation.

Multiscale Modeling of Aqueous Electric Double Layers

From the stability of colloidal suspensions to the charging of electrodes, electric double layers play a pivotal role in aqueous systems. The interactions between interfaces, water molecules, ions and other solutes making up the electrical double layer span length scales from Ångströms to micrometers and are notoriously complex. Therefore, explaining experimental observations in terms of the double layer's molecular structure has been a long-standing challenge in physical chemistry, yet recent advances in simulations techniques and computational power have led to tremendous progress. In particular, the past decades have seen the development of a multiscale theoretical framework based on the combination of quantum density functional theory, force-field based simulations and continuum theory. In this Review, we discuss these theoretical developments and make quantitative comparisons to experimental results from, among other techniques, sum-frequency generation, atomic-force microscopy, and electrokinetics. Starting from the vapor/water interface, we treat a range of qualitatively different types of surfaces, varying from soft to solid, from hydrophilic to hydrophobic, and from charged to uncharged.

Free-energy decomposition of salt effects on the solubilities of small molecules and the role of excluded-volume effects

The roles of cations and anions are different in the perturbation on solvation, and thus, the analyses of the separated contributions from cations and anions are useful to establish molecular pictures of ion-specific effects. In this work, we investigate the effects of cations, anions, and water separately in the solvation of n-alcohols and n-alkanes by free-energy decomposition. By utilising energy-representation theory of solvation, we address the contributions arising from the direct solute-solvent interactions and the excluded-volume effects. It is found that the change in solvation of n-alcohols and n-alkanes upon addition of salt depends primarily on the anion species. The direct interaction between the anion and solute is in agreement with the Setschenow coefficient in terms of the ranking of salting-in and salting-out for n-alkanes, which corresponds to the extent of accumulation of the anion on the solute surface. For each of the n-alcohols and n-alkanes examined, the excluded-volume component in the Setschenow coefficient is well correlated to the (total) Setschenow coefficient when the salt effects are concerned. The ranking of the excluded-volume component in the variation of the salt species is parallel to the water contribution, which is correlated further to the change in the water density upon the addition of the salt.

Agglomeration Drives the Reversed Fractionation of Aqueous Carbonate and Bicarbonate at the Air-Water Interface

In the course of our investigations of the adsorption of ions to the air-water interface, we previously reported the surprising result that doubly charged carbonate anions exhibit a stronger surface affinity than singly charged bicarbonate anions. In contrast to monovalent, weakly hydrated anions, which generally show enhanced concentrations in the interfacial region, multivalent (and strongly hydrated) anions are expected to show a much weaker surface propensity. In the present work, we use resonantly enhanced deep-UV second-harmonic generation spectroscopy to measure the Gibbs free energy of adsorption of both carbonate (CO32-) and bicarbonate (HCO3-) anions to the air-water interface. Contrasting the predictions of classical electrostatic theory and in support of our previous findings from X-ray photoelectron spectroscopy, we find that carbonate anions do indeed exhibit much stronger surface affinity than do the bicarbonate anions. Extensive computer simulations reveal that strong ion pairing of CO32- with the Na+ countercation in the interfacial region results in the formation of near-neutral agglomerate clusters, consistent with a theory of interfacial ion adsorption based on hydration free energy and capillary waves. Simulated X-ray photoelectron spectra predict a 1 eV shift in the carbonate spectra compared to that of bicarbonate, further confirming our experiments. These findings not only advance our fundamental understanding of ion adsorption chemistry but also impact important practical processes such as ocean acidification, sea-spray aerosol chemistry, and mammalian respiration physiology.

Hydrophobic dye solubilization via hybrid imogolite nanotubes probed using second harmonic scattering

This article explores the organization and interactions of Disperse Orange 3 (DO3) hydrophobic dye molecules within hybrid organic-inorganic imogolite nanotubes. In pure water, the DO3 dye molecules self assemble into large insoluble 2D nanosheets whose structure is also explored by molecular dynamics simulations. The dye molecules are however efficiently solubilized in the presence of hybrid imogolite nanotubes. The filling of the internal hydrophobic cavity of the nanotubes is quantified. The organization of the molecules inside the nanotube is probed using the polarization resolved second harmonic scattering (SHS) technique coupled with simulation. At the highest loading, the dyes fill the nanotube with their principal axis parallel to the nanotube walls showing a strong SHS signal due to this encapsulation.

Self-Driven Photo-Polarized Water Molecule-Triggered Graphene-Based Photodetector

Flowing water can be used as an energy source for generators, providing a major part of the energy for daily life. However, water is rarely used for information or electronic devices. Herein, we present the feasibility of a polarized liquid-triggered photodetector in which polarized water is sandwiched between graphene and a semiconductor. Due to the polarization and depolarization processes of water molecules driven by photogenerated carriers, a photo-sensitive current can be repeatedly produced, resulting in a high-performance photodetector. The response wavelength of the photodetector can be fine-tuned as a result of the free choice of semiconductors as there is no requirement of lattice match between graphene and the semiconductors. Under zero voltage bias, the responsivity and specific detectivity of Gr/NaCl (0.5 M)W/N-GaN reach values of 130.7 mA/W and 2.3 × 10⁹ Jones under 350 nm illumination, respectively. Meanwhile, using a polar liquid photodetector can successfully read the photoplethysmography signals to produce accurate oxygen blood saturation and heart rate. Compared with the commercial pulse oximetry sensor, the average errors of oxygen saturation and heart rate in the designed photoplethysmography sensor are ~1.9% and ~2.1%, respectively. This study reveals that water can be used as a high-performance photodetector in informative industries.

Liquid Water: When Hyperpolarizability Fluctuations Boost and Reshape the Second Harmonic Scattering Intensities

Second harmonic scattering (SHS) is a method of choice to investigate the molecular structure of liquids. While a clear interpretation of SHS intensity exists for diluted solutions of dyes, the scattering due to solvents remains difficult to interpret quantitatively. Here, we report a quantum mechanics/molecular mechanics (QM/MM) approach to model the polarization-resolved SHS intensity of liquid water, quantifying different contributions to the signal. We point out that the molecular hyperpolarizability fluctuations and correlations cannot be neglected. The intermolecular orientational and hyperpolarizability correlations up to the third solvation layer strongly increase the scattering intensities and modulate the polarization-resolved oscillation that is predicted here by QM/MM without fitting parameters. Our approach can be generalized to other pure liquids to provide a quantitative interpretation of SHS intensities in terms of short-range molecular ordering.

High throughput wide field second harmonic imaging of giant unilamellar vesicles

Cell-sized giant unilamellar vesicles (GUVs) are an ideal tool for understanding lipid membrane structure and properties. Label-free spatiotemporal images of their membrane potential and structure would greatly aid the quantitative understanding of membrane properties. In principle, second harmonic imaging is a great tool to do so, but the low degree of spatial anisotropy that arises from a single membrane limits its application. Here, we advance the use of wide-field high throughput SH imaging by SH imaging with the use of ultrashort laser pulses. We achieve a throughput improvement of 78% of the maximum theoretical value and demonstrate subsecond image acquisition times. We show how the interfacial water intensity can be converted into a quantitative membrane potential map. Finally, for GUV imaging, we compare this type of nonresonant SH imaging to resonant SH imaging and two photon imaging using fluorophores.

Decreased Water Mobility Contributes To Increased α-Synuclein Aggregation

The solvation shell is essential for the folding and function of proteins, but how it contributes to protein misfolding and aggregation has still to be elucidated. We show that the mobility of solvation shell H2 O molecules influences the aggregation rate of the amyloid protein α-synuclein (αSyn), a protein associated with Parkinson's disease. When the mobility of H2 O within the solvation shell is reduced by the presence of NaCl, αSyn aggregation rate increases. Conversely, in the presence CsI the mobility of the solvation shell is increased and αSyn aggregation is reduced. Changing the solvent from H2 O to D2 O leads to increased aggregation rates, indicating a solvent driven effect. We show the increased aggregation rate is not directly due to a change in the structural conformations of αSyn, it is also influenced by a reduction in both the H2 O mobility and αSyn mobility. We propose that reduced mobility of αSyn contributes to increased aggregation by promoting intermolecular interactions.

Decreased Water Mobility Contributes To Increased α-Synuclein Aggregation

The solvation shell is essential for the folding and function of proteins, but how it contributes to protein misfolding and aggregation has still to be elucidated. We show that the mobility of solvation shell H2O molecules influences the aggregation rate of the amyloid protein α-synuclein (αSyn), a protein associated with Parkinson's disease. When the mobility of H2O within the solvation shell is reduced by the presence of NaCl, αSyn aggregation rate increases. Conversely, in the presence CsI the mobility of the solvation shell is increased and αSyn aggregation is reduced. Changing the solvent from H2O to D2O leads to increased aggregation rates, indicating a solvent driven effect. We show the increased aggregation rate is not directly due to a change in the structural conformations of αSyn, it is also influenced by a reduction in both the H2O mobility and αSyn mobility. We propose that reduced mobility of αSyn contributes to increased aggregation by promoting intermolecular interactions.

Accurate Modeling of Bromide and Iodide Hydration with Data-Driven Many-Body Potentials

Ion-water interactions play a central role in determining the properties of aqueous systems in a wide range of environments. However, a quantitative understanding of how the hydration properties of ions evolve from small aqueous clusters to bulk solutions and interfaces remains elusive. Here, we introduce the second generation of data-driven many-body energy (MB-nrg) potential energy functions (PEFs) representing bromide-water and iodide-water interactions. The MB-nrg PEFs use permutationally invariant polynomials to reproduce two-body and three-body energies calculated at the coupled cluster level of theory, and implicitly represent all higher-body energies using classical many-body polarization. A systematic analysis of the hydration structure of small Br^-(H₂O)_n and I^-(H₂O)_n clusters demonstrates that the MB-nrg PEFs predict interaction energies in quantitative agreement with "gold standard" coupled cluster reference values. Importantly, when used in molecular dynamics simulations carried out in the isothermal-isobaric ensemble for single bromide and iodide ions in liquid water, the MB-nrg PEFs predict extended X-ray absorption fine structure (EXAFS) spectra that accurately reproduce the experimental spectra, which thus allows for characterizing the hydration structure of the two ions with a high level of confidence.

Three-Dimensional Confinement of Water: H2O Exhibits Long-Range (>50 nm) Structure while D2O Does Not

Water is the liquid of life thanks to its three-dimensional adaptive hydrogen (H)-bond network. Confinement of this network may lead to dramatic structural changes influencing chemical and physical transformations. Although confinement effects occur on a <1 nm length scale, the upper length scale limit is unknown. Here, we investigate 3D-confinement over lengths scales ranging from 58-140 nm. By confining water in zwitterionic liposomes of different sizes and measuring the change in H-bond network conformation using second harmonic scattering (SHS), we determined long-range confinement effects in light and heavy water. D₂O displays no detectable 3D-confinement effects <58 nm (<3 × 10⁶ D₂O molecules). H₂O is distinctly different. The vesicle enclosed inner H-bond network has a different conformation compared to the outside network and the SHS response scales with the volume of the confining space. H₂O displays confinement effects over distances >100 nm (>2 × 10⁷ H₂O molecules).

Charge Gradients around Dendritic Voids Cause Nanoscale Inhomogeneities in Liquid Water

Water is the matrix of life and serves as a solvent for numerous physical and chemical processes. The origins of the nature of inhomogeneities that exist in liquid water and the time scales over which they occur remains an open question. Here, we report femtosecond elastic second harmonic scattering (fs-ESHS) of liquid water in comparison to an isotropic liquid (CCl₄) and show that water is indeed a nonuniform liquid. The coherent fs-ESHS intensity was interpreted, using molecular dynamics simulations, as arising from charge density fluctuations with enhanced nanoscale polarizabilities around transient voids having an average lifetime of 300 fs. Although voids were also present in CCl₄, they were not characterized by hydrogen bond defects and did not show strong polarizability fluctuations, leading to fs-ESHS of an isotropic liquid. The voids increased in number at higher temperatures above room temperature, in agreement with the fs-ESHS results.

Physical insight into the entropy-driven ion association

The ion association is widely believed to be dominated by the favorable entropy change arising from the release of water molecules from ion hydration shells. However, no direct thermodynamic evidence exists to validate the reliability and suitability of this view. Herein, we employ complicated free energy calculations to rigorously split the free energy including its entropic and enthalpic components into the water-induced contributions and ion-ion interaction terms for several ion pairs from monatomic to polyatomic ions, spanning the size range from small kosmotropes to large chaotropes (Na⁺ , Cs⁺ , Ca²⁺ , F^- , I^- , CO₃ ^2- , and HPO₄ ^2- ). Our results successfully reveal that though ion associations are indeed determined by a delicate balance between the favorable entropy variation and the repulsive enthalpy change, the entropy gain dominated by the solvent occurs only for the monatomic ion pairing. The water-induced entropic contribution significantly goes against the ion pairing between polyatomic anion and cation, which is, alternatively, dominated by the favorable entropy from the ion-ion interaction term, due to the configurational arrangement of polyatomic anions involved in ion association. The structural and dynamic analysis demonstrates that the entropy penalty from the water phase is primarily ascribed to the enhanced stability of water molecules around the cation imposed by the incoming anion. Our study successfully provides a fundamental understanding of water-mediated ion associations and highlights disparate lengthscale dependencies of the dehydration thermodynamics on the specific types of ions.

The Role of Surface Chemistry in the Orientational Behavior of Water at an Interface

Molecular dynamics studies have demonstrated that molecular water at an interface, with either a gas or a solid, displays anisotropic orientational behavior in contrast to its bulk counterpart. This effect has been recently implicated in the like-charge attraction problem for colloidal particles in solution. Here, negatively charged particles in solution display a long-ranged attraction where continuum electrostatic theory predicts monotonically repulsive interactions, particularly in solutions with monovalent salt ions at low ionic strength. Anisotropic orientational behavior of solvent molecules at an interface gives rise to an excess interfacial electrical potential which we suggest generates an additional solvation contribution to the total free energy that is traditionally overlooked in continuum descriptions of interparticle interactions in solution. In the present investigation we perform molecular dynamics simulation based calculations of the interfacial potential using realistic surface models representing various chemistries as well as different solvents. Similar to previous work that focused on simple model surfaces constructed by using oxygen atoms, we find that solvents at more realistic model surfaces exhibit substantial anisotropic orientational behavior. We explore the dependence of the interfacial solvation potential on surface properties such as surface group chemistry and group density at silica and carboxylated polystyrene interfaces. For water, we note surprisingly good agreement between results obtained for a simple O-atom wall and more complex surface models, suggesting a general qualitative consistency of the interfacial solvation effect for surfaces in contact with water. In contrast, for an aprotic solvent such as DMSO, surface chemistry appears to exert a stronger influence on the sign and magnitude of the interfacial solvation potential. The study carries broad implications for molecular-scale interactions and may find relevance in explaining a range of phenomena in soft-matter physics and cell biology.

Time-resolved terahertz-Raman spectroscopy reveals that cations and anions distinctly modify intermolecular interactions of water

The solvation of ions changes the physical, chemical and thermodynamic properties of water, and the microscopic origin of this behaviour is believed to be ion-induced perturbation of water's hydrogen-bonding network. Here we provide microscopic insights into this process by monitoring the dissipation of energy in salt solutions using time-resolved terahertz-Raman spectroscopy. We resonantly drive the low-frequency rotational dynamics of water molecules using intense terahertz pulses and probe the Raman response of their intermolecular translational motions. We find that the intermolecular rotational-to-translational energy transfer is enhanced by highly charged cations and is drastically reduced by highly charged anions, scaling with the ion surface charge density and ion concentration. Our molecular dynamics simulations reveal that the water-water hydrogen-bond strength between the first and second solvation shells of cations increases, while it decreases around anions. The opposite effects of cations and anions on the intermolecular interactions of water resemble the effects of ions on the stabilization and denaturation of proteins.

Connection between water's dynamical and structural properties: Insights from ab initio simulations

SignificanceFirst-principles calculations, which explicitly account for the electronic structure of matter, can shed light on the molecular structure and dynamics of water in its supercooled state. In this work, we use density functional theory, which relies on a functional to describe electronic exchange and correlations, to evaluate which functional best describes the temperature evolution of bulk water transport coefficients. We also assess the validity of the Stokes-Einstein relation for all the functionals in the temperature range studied, and explore the link between structure and dynamics. Based on these results, we show how transport coefficients can be computed from structural descriptors, which require shorter simulation times to converge, and we point toward strategies to develop better functionals.

Water does not dance as ions sing: A new approach in elucidation of ion-invariant water fluctuations

Aqueous solutions of salts composed from monovalent ions are explored using temperature-dependent FT-IR spectroscopy in transmission. Water combination band, being extremely sensitive to the network of hydrogen bonds due to the contribution of water librations (ρ_LH₂O), is analyzed in uni- and multivariate fashion. Univariate analysis of the combination band maximum (ν_max) reveals that perturbation of water hydrogen bond network by ions is primary driven by electrostatic interactions between water and ions. Using multivariate curve resolution with alternating least squares and evolving factor analysis this band is separated into two components that represent low- and high-density water. The observed asymmetry in their behavior is interpreted in terms of fluctuations of a hydrogen bond network of two water components. The significance of the found phenomenon is unambiguously confirmed by performing analogous analysis in the spectral range that contains partial signature of water linear bending (δHOH) and is free from ρ_LH₂O, in which the asymmetry is absent. Additionally, we show that this phenomenon, namely ion-invariant behavior of water fluctuations, persists even in the regime of high ionic strengths. Although ions indeed participate in shaping of water hydrogen bond network, this straightforward approach shows that its temperature-dependent fluctuations are ion-independent.

Fluorescence of KCl Aqueous Solution: A Possible Spectroscopic Signature of Nucleation

Ion pairing in water solutions alters both the water hydrogen-bond network and ion solvation, modifying the dynamics and properties of electrolyte water solutions. Here, we report an anomalous intrinsic fluorescence of KCl aqueous solution at room temperature and show that its intensity increases with the salt concentration. From the ab initio density functional theory (DFT) and time-dependent DFT modeling, we propose that the fluorescence emission could originate from the stiffening of the hydrogen bond network in the hydration shell of solvated ion-pairs that suppresses the fast nonradiative decay and allows the slower radiative channel to become a possible decay pathway. Because computations suggest that the fluorophores are the local ion-water structures present in the prenucleation phase, this band could be the signature of the incoming salt precipitation.

Dissolving salt is not equivalent to applying a pressure on water

Salt water is ubiquitous, playing crucial roles in geological and physiological processes. Despite centuries of investigations, whether or not water’s structure is drastically changed by dissolved ions is still debated. Based on density functional theory, we employ machine learning based molecular dynamics to model sodium chloride, potassium chloride, and sodium bromide solutions at different concentrations. The resulting reciprocal-space structure factors agree quantitatively with neutron diffraction data. Here we provide clear evidence that the ions in salt water do not distort the structure of water in the same way as neat water responds to elevated pressure. Rather, the computed structural changes are restricted to the ionic first solvation shells intruding into the hydrogen bond network, beyond which the oxygen radial-distribution function does not undergo major change relative to neat water. Our findings suggest that the widely cited pressure-like effect on the solvent in Hofmeister series ionic solutions should be carefully revisited.

By advanced machine learning techniques, first-principles simulations find that dissolving salt in water does not change water structure drastically. It is contrary to the notion of “pressure effect” which has been widely applied over past 25 years.

Driving Forces of Cationic Dye Adsorption, Confinement, and Long-Range Correlation in Zeolitic Materials

Zeolitic materials are commonly used to capture emergent contaminants in water or complex aqueous effluents. The efficiency of this adsorption depends strongly on the guest-host interactions and on the surrounding environment with possible coadsorption of the solvent. Only a few experimental techniques are available to probe in situ the sequestration processes at the solid/liquid interface. We propose in the present work to combine the second harmonic scattering technique with isothermal titration calorimetry in order to investigate the adsorption and the confinement of a hemicyanine dye adsorbed inside faujasite materials. The methodology described here permits the quantification of the correlations between the confined dyes in the material and thus gives local information about the organization at the nanometer scale. Various impacts, such as the effect of the solvent type and the silicon to aluminum ratio of the zeolitic adsorbent, are quantitatively estimated and discussed. This work highlights that the most correlated system matches the higher adsorption capacity associated with the lower entropic contribution.

Nonmonotonic Pair Potentials in the Interaction of Like-Charged Objects in Solution

We consider the long-standing like-charge attraction problem, wherein under certain conditions, similarly charged spheres suspended in aqueous electrolyte have been observed to display a minimum in their interaction potential, contrary to the intuitively expected monotonically varying repulsion. Recently, we described an interfacial mechanism invoking the molecular nature of the solvent that explains this anomalous experimental observation. In our model for the interaction of negatively charged particles in water, the minimum in the pair potential results from the superposition of competing contributions to the total free energy. One of these contributions is the canonical repulsive electrostatic term, whereas the other is a solvation-induced attractive contribution. We find that whereas both contributions grow approximately exponentially with decreasing interparticle separation, the occurrence of a stable, long-ranged minimum in the pair potential arises from differences in the precise interparticle separation dependence of the two terms. Specifically, the interfacial solvation term exhibits a more gradual decay with distance than the electrostatic repulsion, permitting the attractive contribution to dominate the interaction at large distances. Importantly, these disparities become evident in quantities calculated from exact numerical solutions to the governing nonlinear Poisson-Boltzmann (PB) equation for the spatial electrical potential distribution in the system. In marked contrast, we find that the linearized PB equation, applicable in the regime of low surface electrical potentials, does not support nonmonotonic trends in the total interaction free energy within the present model. Our results point to the importance of exact descriptions of electrostatic interactions in real systems that most often do not subscribe to particular mathematical limits where analytical approximations may provide a sufficiently accurate description of the problem.

Multistep Micellization of Standard Surfactants Evidenced by Second Harmonic Scattering

Processes involving in solution a reduced number of molecules are difficult to identify and characterize. Here, we show that micellization of standard surfactants, namely sodium dodecyl sulfate and trimethyl tetradecyl ammonium bromide, two nonefficient compounds for quadratic nonlinear optics, can be investigated by second harmonic scattering (SHS). In particular, the formation of aggregates at concentrations smaller than the critical micellar concentration is evidenced through a nonmonotonic behavior of the SHS intensity as a function of the surfactant concentration. A simple model based on chemical equilibria between monomers and micelles is proposed to account for the experimental observations. Signature of long-range molecular orientation correlation is revealed by polarization resolved experiments and is discussed regarding micellization and charge-induced effects.

Single-Molecule Frequency Fingerprint for Ion Interaction Networks in a Confined Nanopore

The transport of molecules and ions through biological nanopores is governed by interaction networks among restricted ions, transported molecules, and residue moieties at pore inner walls. However, identification of such weak ion fluctuations from only few tens of ions inside nanopore is hard to achieve owing to electrochemical measurement limitations. Here, we developed an advanced frequency method to achieve qualitative and spectral analysis of ion interaction networks inside a nanopore. The peak frequency f_m reveals the dissociation rate between nanopore and ions; the peak amplitude a_m depicts the amount of combined ions with the nanopore after interaction equilibrium. A mathematical model for single-molecule frequency fingerprint achieved the prediction of interaction characteristics of mutant nanopores. This single-molecule frequency fingerprint is important for classification, characterization, and prediction of synergetic interaction networks inside nanoconfinement.

Probing Interfacial Effects on Ionization Energies: The Surprising Banality of Anion-Water Hydrogen Bonding at the Air/Water Interface

Liquid microjet photoelectron spectroscopy is an increasingly common technique to measure vertical ionization energies (VIEs) of aqueous solutes, but the interpretation of these experiments is subject to questions regarding sensitivity to bulk versus interfacial solvation environments. We have computed aqueous-phase VIEs for a set of inorganic anions, using a combination of molecular dynamics simulations and electronic structure calculations, with results that are in excellent agreement with experiment regardless of whether the simulation data are restricted to ions at the air/water interface or to those in bulk aqueous solution. Although the computed VIEs are sensitive to ion-water hydrogen bonding, we find that the short-range solvation structure is sufficiently similar in both environments that it proves impossible to discriminate between the two on the basis of the VIE, a conclusion that has important implications for the interpretation of liquid-phase photoelectron spectroscopy. More generally, analysis of the simulation data suggests that the surface activity of soft anions is largely a second or third solvation shell effect, arising from disruption of water-water hydrogen bonds and not from significant changes in first-shell anion-water hydrogen bonding.

Quadrupole-mediated dielectric response and the charge-asymmetric solvation of ions in water

Treating water as a linearly responding dielectric continuum on molecular length scales allows very simple estimates of the solvation structure and thermodynamics for charged and polar solutes. While this approach can successfully account for basic length and energy scales of ion solvation, computer simulations indicate not only its quantitative inaccuracies but also its inability to capture some basic and important aspects of microscopic polarization response. Here, we consider one such shortcoming, a failure to distinguish the solvation thermodynamics of cations from that of otherwise-identical anions, and we pursue a simple, physically inspired modification of the dielectric continuum model to address it. The adaptation is motivated by analyzing the orientational response of an isolated water molecule whose dipole is rigidly constrained. Its free energy suggests a Hamiltonian for dipole fluctuations that accounts implicitly for the influence of higher-order multipole moments while respecting constraints of molecular geometry. We propose a field theory with the suggested form, whose nonlinear response breaks the charge symmetry of ion solvation. An approximate variational solution of this theory, with a single adjustable parameter, yields solvation free energies that agree closely with simulation results over a considerable range of solute size and charge.

Isomer-Specific Vibrational Spectroscopy of Microhydrated Lithium Dichloride Anions: Spectral Fingerprint of Solvent-Shared Ion Pairs

The vibrational spectroscopy of lithium dichloride anions microhydrated with one to three water molecules, [LiCl2 (H2 O)1-3 ]- , is studied in the OH stretching region (3800-2800 cm-1 ) using isomer-specific IR/IR double-resonance population labelling experiments. The spectroscopic fingerprints of individual isomers can only be unambiguously assigned after anharmonic effects are considered, but then yield molecular level insight into the onset of salt dissolution in these gas phase model systems. Based on the extent of the observed frequency shifts ΔνOH of the hydrogen-bonded OH stretching oscillators solvent-shared ion pair motifs (<3200 cm-1 ) can be distinguished from intact-core structures (>3200 cm-1 ). The characteristic fingerprint of a water molecule trapped directly in-between two ions of opposite charge provides an alternative route to evaluate the extent of ion pairing in aqueous electrolyte solutions.

The many-body expansion for aqueous systems revisited: III. Hofmeister ion-water interactions

We report a Many Body Energy (MBE) analysis of aqueous ionic clusters containing anions and cations at the two opposite ends of the Hofmeister series, viz. the kosmotropes Ca²⁺ and SO₄^2- and the chaotropes NH₄⁺ and ClO₄^-, with 9 water molecules to quantify how these ions alter the interaction between the water molecules in their immediate surroundings. We specifically aim at quantifying how various ions (depending on their position in the Hofmeister series) affect the interaction between the surrounding water molecules and probe whether there is a qualitatively different behavior between kosmotropic vs. chaotropic ions. The current results when compared to the ones reported earlier for water clusters [J. P. Heindel and S. S. Xantheas, J. Chem. Theor. Comput., 2020, 16, 6843-6855] as well as for alkali metal and halide ion aqueous clusters of the same size [J. P. Heindel and S. S. Xantheas, J. Chem. Theor. Comput., 2021, 17, 2200-2216], which lie in the middle of the Hofmeister series, offer a complete account of the effect an ion across the Hofmeister series from "kosmotropes" to "chaotropes" has on the interaction between the neighboring water molecules. Through this analysis, noteworthy differences between the MBE of kosmotropes and chaotropes were identified. The MBE of kosmotropes is dominated by ion-water interactions that extend beyond the 4-body term, the rank at which the MBE of pure water converges. The percentage contribution of the 2-B term to the total cluster binding energy is noticeably larger. The disruption of the hydrogen bonded network due to the dominant ion-water interactions results in weak, unfavorable water-water interactions. The MBE for chaotropes, on the other hand, was found to converge more quickly as it more closely resembles that of pure water clusters. Chaotropes exhibit weaker overall binding energies and weaker ion-water interactions in favor of water-water interactions, somewhat recovering the pattern of the 2-4 body terms exemplified by pure water clusters. A remarkable anti-correlation between the 2-B ion-water (I-W) and water-water (W-W) interactions as well as between the 3-B (I-W-W) and (I-W) interactions was found for both kosmotropic and chaotropic ions. This anti-correlation is linear for both monatomic anions and monatomic cations, suggesting the existence of underlying physical mechanisms that were previously unexplored. The consideration of two different structural arrangements (ion inside and outside of a water cluster) suggests that fully solvated (ion inside) chaotropes disrupt the hydrogen bonding network in a similar manner to partially solvated (ion outside) kosmotropes and offers useful insights into the modeling requirements of bulk vs. interfacial ion solvation. It is noteworthy that the 2-B contribution to the total Basis Set Superposition Error (BSSE) correction for both kosmotropic and chaotropic ions follows the universal erf profile vs. intermolecular distance previously reported for pure water, halide ion-water and alkali metal ion-water clusters. When scaled for the corresponding dimer energies and distances, a single profile fits the current results together with all previously reported ones for pure water and halide water clusters. This finding lends further support to schemes for accurately estimating the 2-B BSSE correction in condensed environments.

Protein phase separation and its role in chromatin organization and diseases

In the physical sciences, solid, liquid, and gas are the most familiar phase states, whose essence is their existence reflecting the different spatial distribution of molecular components. The biological molecules in the living cell also have differences in spatial distribution. The molecules organized in the form of membrane-bound organelles are well recognized. However, the biomolecules organized in membraneless compartments called biomolecular condensates remain elusive. The liquid-liquid phase separation (LLPS), as a new emerging scientific breakthrough, describes the biomolecules assembled in special distribution and appeared as membraneless condensates in the form of a new "phase" compared with the surrounding liquid milieu. LLPS provides an important theoretical basis for explaining the composition of biological molecules and related biological reactions. Mounting evidence has emerged recently that phase-separated condensates participate in various biological activities. This article reviews the occurrence of LLPS and underlying regulatory mechanisms for understanding how multivalent molecules drive phase transitions to form the biomolecular condensates. And, it also summarizes recent major progress in elucidating the roles of LLPS in chromatin organization and provides clues for the development of new innovative therapeutic strategies for related diseases.

The Many-Body Expansion for Aqueous Systems Revisited: II. Alkali Metal and Halide Ion-Water Interactions

We present a detailed study of the many-body expansion (MBE) for alkali metal and halide ion-water interactions and quantify the effect of these ions on the strength of the surrounding aqueous hydrogen bonding environment. Building on our previous work on neutral water clusters [J. P. Heindel and S. S. Xantheas, J. Chem. Theor. Comput. 16 (11), 6843-6855 (2020)], we carry out the MBE for the alkali metal and halide ion-water clusters, Z^+/-(H₂O)₉, where Z = Li⁺, Na⁺, K⁺, Rb⁺, Cs⁺, F^-, Cl^-, Br^-, and I^- and compare them with the results for a pure water cluster with the same number of "bodies", viz., (H₂O)₁₀. The 2-B ion-water (I-W) interaction accounts for a larger percentage of the total cluster binding energy compared to a pure water cluster of the same size, with the total 3-B term being smaller and of opposite sign (repulsive), whereas higher order terms are essentially negligible. The same oscillating behavior around zero for the MBE terms higher than the 5-B with a basis set that was reported for water clusters is also observed for the ion-water clusters considered here, with the basis set superposition error (BSSE) corrections amending this as in the water cluster case. A remarkable, linear anticorrelation between the total 2-B (I-W), the total 2-B (W-W), and also the 3-B (W-W-W) interactions is found, quantifying the effect of the different ions in disrupting and altering (weakening) the neighboring hydrogen bonded water network: stronger (I-W) interactions result in weaker (W-W) interactions. Additional linear correlations across the two series of alkali metals and halide ions were found between the 3-B (I-W-W) and the 2-B (I-W) as well as between the 3-B (I-W-W) and the 3-B (W-W-W) interactions, suggesting the existence of previously unrealized underlying physics governing these 2-B intermolecular and 3-B collective interactions. Our results further suggest a universal behavior of the two different families of ions (alkali metals and halides) for both the correlations of the various components of the total binding energies and the estimate of the 2-B BSSE correction, which is reported to follow a common profile for ion-water and water-water interactions when cast in terms of reduced distances and energies of the respective dimers. We expect the current results that quantify the interplay between ion-water and water-water interactions in aqueous clusters to impact the development of classical, ab initio-based force fields for monatomic ion solvation, whereas the insights into the nature of the BSSE to be critical in future ab initio-based, many-body molecular dynamics studies.

Convergence behaviour of solvation shells in simulated liquids

A convenient way to analyse solvent structure around a solute is to use solvation shells, whereby solvent position around the solute is discretised by the size of a solvent molecule, leading to multiple shells around the solute. The two main ways to define multiple shells around a solute are either directly with respect to the solute, called solute-centric, or locally for both solute and solvent molecules alike. It might be assumed that both methods lead to solvation shells with similar properties. However, our analysis suggests otherwise. Solvation shells are analysed in a series of simulations of five pure liquids of differing polarity. Shells are defined locally working outwards from each molecule treated as a reference molecule using two methods: the cutoff at the first minimum in the radial distribution function and the parameter-free Relative Angular Distance method (RAD). The molecular properties studied are potential energy, coordination number and coordination radius. Rather than converging to bulk values, as might be expected for pure solvents, properties are found to deviate as a function of shell index. This behaviour occurs because molecules with larger coordination numbers and radius have more neighbours, which make them more likely to be connected to the reference molecule via fewer shells. The effect is amplified for RAD because of its more variable coordination radii and for water with its more open structure and stronger interactions. These findings indicate that locally defined shells should not be thought of as directly comparable to solute-centric shells or to distance. As well as showing how box size and cutoff affect the non-convergence, to restore convergence we propose a hybrid method by defining a new set of shells with boundaries at the uppermost distance of each locally derived shell.

Hyaluronan orders water molecules in its nanoscale extended hydration shells

Hyaluronan (HA) is an anionic, highly hydrated bio-polyelectrolyte found in the extracellular environment, like the synovial fluid between joints. We explore the extended hydration shell structure of HA in water using femtosecond elastic second-harmonic scattering (fs-ESHS). HA enhances orientational water-water correlations. Angle-resolved fs-ESHS measurements and nonlinear optical modeling show that HA behaves like a flexible chain surrounded by extended shells of orientationally correlated water. We describe several ways to determine the concentration-dependent size and shape of a polyelectrolyte in water, using the amount of water oriented by the polyelectrolyte charges as a contrast agent. The spatial extent of the hydration shell is determined via temperature-dependent measurements and can reach up to 475 nm, corresponding to a length of 1600 water molecules. A strong isotope effect, stemming from nuclear quantum effects, is observed when light water (H2O) is replaced by heavy water (D2O), amounting to a factor of 4.3 in the scattered SH intensity.

NaCl aggregation in water at elevated temperatures and pressures: Comparison of classical force fields

The properties of water vary dramatically with temperature and density. This can be exploited to control its effectiveness as a solvent. Thus, supercritical water is of keen interest as solvent in many extraction processes. The low solubility of salts in lower density supercritical water has even been suggested as a means of desalination. The high temperatures and pressures required to reach supercritical conditions can present experimental challenges during collection of required physical property and phase equilibria data, especially in salt-containing systems. Molecular simulations have the potential to be a valuable tool for examining the behavior of solvated ions at these high temperatures and pressures. However, the accuracy of classical force fields under these conditions is unclear. We have, therefore, undertaken a parametric study of NaCl in water, comparing several salt and water models at 200 bar-600 bar and 450 K-750 K for a range of salt concentrations. We report a comparison of structural properties including ion aggregation, hydrogen bonding, density, and static dielectric constants. All of the force fields qualitatively reproduce the trends in the liquid phase density. An increase in ion aggregation with decreasing density holds true for all of the force fields. The propensity to aggregate is primarily determined by the salt force field rather than the water force field. This coincides with a decrease in the water static dielectric constant and reduced charge screening. While a decrease in the static dielectric constant with increasing NaCl concentration is consistent across all model combinations, the salt force fields that exhibit more ionic aggregation yield a slightly smaller dielectric decrement.

PySHS: Python Open Source Software for Second Harmonic Scattering

The PySHS package is a new python open source software tool which simulates the second harmonic scattering (SHS) of different kinds of colloidal nano-objects in various experimental configurations. This package is able to compute polarizations resolved at a fixed scattered angle or angular distribution for different polarization configurations. This article presents the model implemented in the PySHS software and gives some computational examples. A comparison between computational results and experimental data concerning molecular dye intercalated inside liposomes membrane is presented to illustrate the possibilities with PySHS.

Multi-scale approach for the prediction of atomic scale properties

Electronic nearsightedness is one of the fundamental principles that governs the behavior of condensed matter and supports its description in terms of local entities such as chemical bonds. Locality also underlies the tremendous success of machine-learning schemes that predict quantum mechanical observables - such as the cohesive energy, the electron density, or a variety of response properties - as a sum of atom-centred contributions, based on a short-range representation of atomic environments. One of the main shortcomings of these approaches is their inability to capture physical effects ranging from electrostatic interactions to quantum delocalization, which have a long-range nature. Here we show how to build a multi-scale scheme that combines in the same framework local and non-local information, overcoming such limitations. We show that the simplest version of such features can be put in formal correspondence with a multipole expansion of permanent electrostatics. The data-driven nature of the model construction, however, makes this simple form suitable to tackle also different types of delocalized and collective effects. We present several examples that range from molecular physics to surface science and biophysics, demonstrating the ability of this multi-scale approach to model interactions driven by electrostatics, polarization and dispersion, as well as the cooperative behavior of dielectric response functions.

Long-Range Orientational Organization of Dipolar and Steric Liquids

Long-range orientational correlations in liquids have received recent renewed interest, in particular for the neat water case. These long-range orientational correlations, exceeding several tens of nanometers, originate from the presence of the strong permanent water dipolar moment. However, the exact dependence with the dipolar moment and the role of other local forces like steric hindrance has never been addressed. In this work, we experimentally measure long-range correlations for a set of liquids differing by their molecular weight and dipolar moment, in order to reveal the origin of their long-range organization. Hence, we show that the dipolar moment of a solvent molecule is not the unique feature determining the orientational correlation. Steric hindrance significantly helps to structure the liquids as well. In order to quantify these long-range correlations, we also derive theoretically the polarization resolved second harmonic scattering intensity as a function of the rotational invariants describing the dipolar and octupolar interaction.

Hydrogen-Bond Structure and Low-Frequency Dynamics of Electrolyte Solutions: Hydration Numbers from ab Initio Water Reorientation Dynamics and Dielectric Relaxation Spectroscopy

We present an atomistic simulation scheme for the determination of the hydration number (h) of aqueous electrolyte solutions based on the calculation of the water dipole reorientation dynamics. In this methodology, the time evolution of an aqueous electrolyte solution generated from ab initio molecular dynamics simulations is used to compute the reorientation time of different water subpopulations. The value of h is determined by considering whether the reorientation time of the water subpopulations is retarded with respect to bulk-like behavior. The application of this computational protocol to magnesium chloride (MgCl2 ) solutions at different concentrations (0.6-2.8 mol kg-1 ) gives h values in excellent agreement with experimental hydration numbers obtained using GHz-to-THz dielectric relaxation spectroscopy. This methodology is attractive because it is based on a well-defined criterion for the definition of hydration number and provides a link with the molecular-level processes responsible for affecting bulk solution behavior. Analysis of the ab initio molecular dynamics trajectories using radial distribution functions, hydrogen bonding statistics, vibrational density of states, water-water hydrogen bonding lifetimes, and water dipole reorientation reveals that MgCl2 has a considerable influence on the hydrogen bond network compared with bulk water. These effects have been assigned to the specific strong Mg-water interaction rather than the Cl-water interaction.

Nanoscale kinetics of amorphous calcium carbonate precipitation in H2O and D2O

Calcium carbonate (CaCO₃) is one of the most well-studied and abundant natural materials on Earth. Crystallisation of CaCO₃ is often observed to proceed via an amorphous calcium carbonate (ACC) phase, as a precursor to more stable crystalline polymorphs such as vaterite and calcite. Despite its importance, the kinetics of ACC formation have proved difficult to study, in part due to rapid precipitation at moderate supersaturations, and the instability of ACC with respect to all other polymorphs. However, ACC can be stabilised under confinement conditions, such as those provided by a nanopipette. This paper demonstrates electrochemical mixing of a Ca²⁺ salt (CaCl₂) and a HCO₃^- salt (NaHCO₃) in a nanopipette to repeatedly and reversibly precipitate nanoparticles of ACC under confined conditions, as confirmed by scanning transmission electron microscopy (STEM). Measuring the current as a function of applied potential across the end of the nanopipette and time provides millisecond-resolved measurements of the induction time for ACC precipitation. We demonstrate that under conditions of electrochemical mixing, ACC precipitation is extremely fast, and highly pH sensitive with an apparent third order dependence on CO₃^2- concentration. Furthermore, the rate is very similar for the equivalent CO₃^2- concentrations in D₂O, suggesting that neither ion dehydration nor HCO₃^- deprotonation represent significant energetic barriers to the formation of ACC. Finite element method simulations of the electrochemical mixing process enable the supersaturation to be estimated for all conditions and accurately predict the location of precipitation.

pH-Dependent Water Clusters in Photoacid Solution: Real-Time Observation by ToF-SIMS at a Submicropore Confined Liquid-Vacuum Interface

Water clusters are ubiquitously formed in aqueous solutions by hydrogen bonding, which is quite sensitive to various environment factors such as temperature, pressure, electrolytes, and pH. Investigation of how the environment has impact on water structure is important for further understanding of the nature of water and the interactions between water and solutes. In this work, pH-dependent water structure changes were studied by monitoring the changes for the size distribution of protonated water clusters by in-situ liquid ToF-SIMS. In combination with a light illumination system, in-situ liquid ToF-SIMS was used to real-time measure the changes of a light-activated organic photoacid under different light illumination conditions. Thus, the proton transfer and pH-mediated water cluster changes were analyzed in real-time. It was found that higher concentration of free protons could lead to a strengthened local hydrogen bonding network as well as relatively larger protonated water clusters in both organic acid and inorganic acid. Besides, the accumulation of protons at the liquid-vacuum interface under light illumination was observed owing to the affinity of organic molecules to the low-pressure gas phase. The application of in-situ liquid ToF-SIMS analysis in combination with in-situ light illumination system opened up an avenue to real-time investigate light-activated reactions. Besides, the results regarding water structure changes in acidic solutions showed important insights in related atmospheric and physiochemical processes.

Slowing Down of the Molecular Reorientation of Water in Concentrated Alkaline Solutions

It is generally accepted that the hydroxide ion (OH^-) is a strong hydrogen bond acceptor and that its anomalously high diffusion constant in water results from a Grotthuss-like structural diffusion mechanism. However, the spatial extent over which OH^- ions influence the dynamics of the hydrogen-bond network of water remained largely unclear. Here, we measure the ultrafast dynamics of OH groups of HDO molecules interacting with the deuterated hydroxide ion OD^-. For solutions with OD^- concentrations up to 4 M, we find that HDO molecules that are not directly interacting with the ions have a reorientation time constant of ∼2.7 ps, similar to that of pure liquid water. When the concentration of OD^- ions is increased, the reorientation time constant increases, indicating a strong slowing down of the structural dynamics of the solution.

Dielectric response with short-ranged electrostatics

The dielectric nature of polar liquids underpins much of their ability to act as useful solvents, but its description is complicated by the long-ranged nature of dipolar interactions. This is particularly pronounced under the periodic boundary conditions commonly used in molecular simulations. In this article, the dielectric properties of a water model whose intermolecular electrostatic interactions are entirely short-ranged are investigated. This is done within the framework of local molecular-field theory (LMFT), which provides a well-controlled mean-field treatment of long-ranged electrostatics. This short-ranged model gives a remarkably good performance on a number of counts, and its apparent shortcomings are readily accounted for. These results not only lend support to LMFT as an approach for understanding solvation behavior, but also are relevant to those developing interaction potentials based on local descriptions of liquid structure.

Effect of Divalent Cations on the Structure and Mechanics of Vimentin Intermediate Filaments

Divalent cations behave as effective cross-linkers of intermediate filaments (IFs) such as vimentin IF (VIF). These interactions have been mostly attributed to their multivalency. However, ion-protein interactions often depend on the ion species, and these effects have not been widely studied in IFs. Here, we investigate the effects of two biologically important divalent cations, Zn²⁺ and Ca²⁺, on VIF network structure and mechanics in vitro. We find that the network structure is unperturbed at micromolar Zn²⁺ concentrations, but strong bundle formation is observed at a concentration of 100 μM. Microrheological measurements show that network stiffness increases with cation concentration. However, bundling of filaments softens the network. This trend also holds for VIF networks formed in the presence of Ca²⁺, but remarkably, a concentration of Ca²⁺ that is two orders higher is needed to achieve the same effect as with Zn²⁺, which suggests the importance of salt-protein interactions as described by the Hofmeister effect. Furthermore, we find evidence of competitive binding between the two divalent ion species. Hence, specific interactions between VIFs and divalent cations are likely to be an important mechanism by which cells can control their cytoplasmic mechanics.

Aqueous solvation of the chloride ion revisited with density functional theory: impact of correlation and exchange approximations

Aqueous chloride is simulated using PBE-D3, PBE0-D3, and SCAN to investigate the impact of exchange and correlation approximations; we find the exact exchange fraction strongly impacts the energetics and polarizability of solvated chloride.