Exploring the Rich Energy Landscape of Sulfate–Water Clusters SO42– (H2O)n=3–7: An Electronic Structure Approach

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.

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.

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.

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.

Cryogenic ion trap vibrational spectroscopy of the microhydrated sulfate dianions SO42-(H2O)3-8

Infrared photodissociation spectra of the D2-tagged microhydrated sulfate dianions with three to eight water molecules are presented over a broad spectral range that covers the OH stretching and H2O bending modes of the solvent molecules at higher energies, the sulfate stretching modes of the solute at intermediate energies and the intermolecular solute librational modes at the lowest energies. A low ion temperature combined with messenger-tagging ensures well-resolved vibrational spectra that allow for structure assignments based on a comparison to harmonic and anharmonic IR spectra from density functional theory (DFT) calculations. DFT ab initio molecular dynamics simulations are required to disentangle the broad and complex spectral signatures of microhydrated sulfate dianions in the OH stretching region and to identify systematic trends in the correlation of the strength and evolution of the solute-solvent and solvent-solvent interactions with cluster size. The onset for the formation of the second solvation shell is observed for n = 8.

Matrix Isolation Spectroscopy: Aqueous p-Toluenesulfonic Acid Solvation

Interaction between p-toluenesulfonic acid (pTSA) and water is studied at -20 °C in a CCl₄ matrix. In CCl₄ water exists as monomers with restricted rotational motion about its symmetry axis. Additionally, CCl₄ is transparent in the hydrogen-bonded region; CCl₄ thus constitutes an excellent ambient thermal energy matrix isolation medium for diagnosing interactions with water. Introducing pTSA-nH₂O gives rise to two narrow resonances at 3642 cm^-1 and at 2835 cm^-1 plus a broad 3000-3550 cm^-1 absorption. In addition, negative monomer symmetric and asymmetric stretch features relative to nominally dry CCl₄ indicate that fewer water monomers exist in the cooled (-20 °C) acid solution than in room-temperature anhydrous CCl₄. The negative peaks along with the broad absorption band indicate that water monomers are incorporated into clusters. The 3642 cm^-1 resonance is assigned to the OH-π interaction with a cluster containing many water molecules per acid molecule. The 2835 cm^-1 resonance is assigned to the (S-)O-H stretch of pTSA-dihydrate. The coexistence of these two species provides insights into interactions in this acid-water CCl₄ system.

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.

Glyoxal and Methylglyoxal Setschenow Salting Constants in Sulfate, Nitrate, and Chloride Solutions: Measurements and Gibbs Energies

Knowledge about Setschenow salting constants, KS, the exponential dependence of Henry's Law coefficients on salt concentration, is of particular importance to predict secondary organic aerosol (SOA) formation from soluble species in atmospheric waters with high salt concentrations, such as aerosols. We have measured KS of glyoxal and methylglyoxal for the atmospherically relevant salts (NH4)2SO4, NH4NO3, NaNO3, and NaCl and find that glyoxal consistently "salts-in" (KS of -0.16, -0.06, -0.065, -0.1 molality(-1), respectively) while methylglyoxal "salts-out" (KS of +0.16, +0.075, +0.02, +0.06 molality(-1)). We show that KS values for different salts are additive and present an equation for use in atmospheric models. Additionally, we have performed a series of quantum chemical calculations to determine the interactions between glyoxal/methylglyoxal monohydrate with Cl(-), NO3(-), SO4(2-), Na(+), and NH4(+) and find Gibbs free energies of water displacement of -10.9, -22.0, -22.9, 2.09, and 1.2 kJ/mol for glyoxal monohydrate and -3.1, -10.3, -7.91, 6.11, and 1.6 kJ/mol for methylglyoxal monohydrate with uncertainties of 8 kJ/mol. The quantum chemical calculations support that SO4(2-), NO3(-), and Cl(-) modify partitioning, while cations do not. Other factors such as ion charge or partitioning volume effects likely need to be considered to fully explain salting effects.

Hydration of the sulfate dianion in cold nanodroplets: SO42−(H2O)12and SO42−(H2O)13

The sulfate ion is found to lie at the cluster surface in both cases.

Separate electronic attenuation allowing a spin-component-scaled second-order Møller-Plesset theory to be effective for both thermochemistry and noncovalent interactions

Spin-component-scaled (SCS) second-order Møller-Plesset perturbation theory (MP2) improves the treatment of thermochemistry and noncovalent interactions relative to MP2, although the optimal scaling coefficients are quite different for thermochemistry versus noncovalent interactions. This work reconciles these two different scaling regimes for SCS-MP2 by using two different length scales for electronic attenuation of the two spin components. The attenuation parameters and scaling coefficients are optimized in the aug-cc-pVTZ (aTZ) basis using the S66 database of intermolecular interactions and the W4-11 database of thermochemistry. Transferability tests are performed for atomization energies and barrier heights, as well as on further test sets for inter- and intramolecular interactions. SCS dual-attenuated MP2 in the aTZ basis, SCS-MP2(2terfc, aTZ), performs similarly to SCS-MP2/aTZ for thermochemistry while frequently outperforming MP2 at the complete basis set limit (CBS) for nonbonded interactions.