Structure of Hydronium (H3O+)/Chloride (Cl−) Contact Ion Pairs in Aqueous Hydrochloric Acid Solution: A Zundel-like Local Configuration

Fast Structural Dynamics in Concentrated HCl Solutions: From Proton Hopping to the Bulk Viscosity

Concentrated acid solutions, particularly HCl, have been studied extensively to examine the proton hopping and infrared spectral signatures of hydronium ions. Much less attention has been given to the structural dynamics of concentrated HCl solutions. Here, we apply optical heterodyne detected-optical Kerr effect (OHD-OKE) measurements to examine HCl concentration-dependent dynamics from moderate (0.8 m) to very high (15.5 m) concentrations and compare the results to the dynamics of NaCl solutions, as Na+ is similar in size to the hydronium cation. Both HCl and NaCl OHD-OKE signals decay as triexponentials at all concentrations, in contrast to pure water, which decays as a biexponential. Two remarkable features of the HCl dynamics are the following: (1) the bulk viscosity is linearly related to the slowest decay constant, t3, and (2) the concentration-dependent proton hopping times, determined by ab initio MD simulations and 2D IR chemical exchange experiments, both obtained from the literature, fall on the same line as the slowest structural dynamics relaxation time, t3, within experimental error. The structural dynamics of hydronium/chloride/water clusters, with relaxation times t3, are responsible for the concentration dependence of microscopic property of proton hopping and the macroscopic bulk viscosity. The slowest time constant (t3), which does not have a counterpart in pure water, is 3 ps at 0.8 m and increases by a factor of ∼2 by 15.5 m. The two fastest HCl decay constants, t1 and t2, are similar to those of pure water and increase mildly with the concentration.

Solvent Dynamics of Aqueous Halides before and after Photoionization

Electron transfer reactions can be strongly influenced by solvent dynamics. We study the photoionization of halides in water as a model system for such reactions. There are no internal nuclear degrees of freedom in the solute, allowing the dynamics of the solvent to be uniquely identified. We simulate the equilibrium solvent dynamics for Cl^-, Br^-, I^-, and their respective neutral atoms in water, comparing quantum mechanical/molecular mechanical (QM/MM) and classical molecular dynamics (MD) methods. On the basis of the obtained configurations, we calculate the extended X-ray absorption fine structure (EXAFS) spectra rigorously based on the MD snapshots and compare them in detail with other theoretical and experimental results available in the literature. We find our EXAFS spectra based on QM/MM MD simulations in good agreement with their experimental counterparts for the ions. Classical MD simulations for the ions lead to EXAFS spectra that agree equally well with the experiment when it comes to the oscillatory period of the signal, even though they differ from the QM/MM radial distribution functions extracted from the MD. The amplitude is, however, considerably overestimated. This suggests that to judge the reliability of theoretical simulation methods or to elucidate fine details of the atomistic dynamics of the solvent based on EXAFS spectra, the amplitude as well as the oscillatory period need to be considered. If simulations fail qualitatively, as does the classical MD for the aqueous neutral halogen atoms, the resulting EXAFS will also be strongly affected in both oscillatory period and amplitude. The good reliability of QM/MM-based EXAFS simulations, together with clear qualitative differences in the EXAFS spectra found between halides and their atomic counterparts, suggests that a combined theory and experimental EXAFS approach is suitable for elucidating the nonequilibrium solvent dynamics in the photoionization of halides and possibly also for electron transfer reactions in more complex systems.

X-Ray absorption spectroscopy on airborne aerosols

Here we demonstrate a method for performing X-ray absorption spectroscopy (XAS) on airborne aerosols. XAS provides unique insight into elemental composition, chemical and phase state, local coordination and electronic structure of both crystalline and amorphous matter. The aerosol is generated from different salt solutions using a commercial atomizer and dried using a diffusion drier. Embedded in a carrier gas, the aerosol is guided into the experimental chamber for XAS analysis. Typical particle sizes range from some 10 to a few 100 nm. Inside the chamber the aerosol bearing gas is then confined into a region of about 1-2 cm³ in size, by a pure flow of helium, generating a stable free-flowing stream of aerosol. It is hit by a monochromatic X-ray beam, and the emitted fluorescent light is used for spectroscopic analysis. Using an aerosol generated from CaCl₂, KCl, and (NH₄)₂SO₄ salt solutions, we demonstrate the functionality of the system in studying environmentally relevant systems. In addition, we show that the detection limits are sufficient to also observe subtle spectroscopic signatures in XAS spectra with integration times of about 1-2 hours using a bright undulator beamline. This novel setup opens new research opportunities for studying the nucleation of new phases in multicomponent aerosol systems in situ, and for investigating (photo-) chemical reactions on airborne matter, as relevant to both atmospheric science and also for general chemical application.

Water Defect Stabilizes the Bi3+ Lone-Pair Electronic State Leading to an Unusual Aqueous Hydration Structure

The aqueous hydration structure of the Bi³⁺ ion is probed using a combination of extended X-ray absorption fine structure (EXAFS) spectroscopy and density functional theory (DFT) simulations of ion-water clusters and condensed-phase solutions. Anomalous features in the EXAFS spectra are found to be associated with a highly asymmetric first-solvent water shell. The aqueous chemistry and structure of the Bi³⁺ ion are dramatically controlled by the water stabilization of a lone-pair electronic state involving the mixed 6s and 6p orbitals. This leads to a distinct multimodal distribution of water molecules in the first shell that are separated by about 0.2 Å. The lone-pair structure is stabilized by a collective response of multiple waters that are localized near the lone-pair anti-bonding site. The findings indicate that the lone-pair stereochemistry of aqueous Bi³⁺ ions plays a major role in the binding of water and ligands in aqueous solutions.

Spectral signatures of excess-proton waiting and transfer-path dynamics in aqueous hydrochloric acid solutions

The theoretical basis for linking spectral signatures of hydrated excess protons with microscopic proton-transfer mechanisms has so far relied on normal-mode analysis. We introduce trajectory-decomposition techniques to analyze the excess-proton dynamics in ab initio molecular-dynamics simulations of aqueous hydrochloric-acid solutions beyond the normal-mode scenario. We show that the actual proton transfer between two water molecules involves for relatively large water-water separations crossing of a free-energy barrier and thus is not a normal mode, rather it is characterized by two non-vibrational time scales: Firstly, the broadly distributed waiting time for transfer to occur with a mean value of 200-300 fs, which leads to a broad and weak shoulder in the absorption spectrum around 100 cm-1, consistent with our experimental THz spectra. Secondly, the mean duration of a transfer event of about 14 fs, which produces a rather well-defined spectral contribution around 1200 cm-1 and agrees in location and width with previous experimental mid-infrared spectra.

Dynamics of Ion Pairing in Dilute Aqueous HCl Solutions by Spectroscopic Measurements of Hydroxyl Radical Conversion into Dichloride Radical Anions

The rate of formation of dichloride anions (Cl₂^•–) in dilute aqueous solutions of HCl (2–100 mmol·kg^–1) was measured by the technique of pulse radiolysis over the temperature range of 288–373 K. The obtained Arrhenius dependence shows a concentration averaged activation energy of 7.3 ± 1.8 kJ·mol^–1, being half of that expected from the mechanism assuming the ^•OHCl^– intermediate and supporting the ionic equilibrium-based mechanism, i.e., the formation of Cl₂^•– in the reaction of ^•OH with a hydronium–chloride (Cl^–·H₃O⁺) contact ion pair. Assuming diffusion-controlled encounter of the hydronium and chloride ions and including the effect of the ionic atmosphere, we showed that the reciprocal of τ, the lifetime of (Cl^–·H₃O⁺), follows an Arrhenius dependence with an activation energy of 23 ± 4 kJ·mol^–1, independent of the acid concentration. This result indicates that the contact pair is stabilized by hydrogen bonding interaction of the solvent molecules. We also found that at a fixed temperature, τ is noticeably increased in less-concentrated solutions (m_HCl < 0.01 m). Since this concentration effect is particularly pronounced at near ambient temperatures, the increasing pair lifetime may result from the solvent cage effect enhanced by the presence of large supramolecular structures (patches) formed by continuously connected four-bonded water molecules.

Thermodynamic Modeling of Salting Effects in Solvent Extraction of Cobalt(II) from Chloride Media by the Basic Extractant Methyltrioctylammonium Chloride

The design and optimization of solvent extraction processes for metal separations are challenging tasks due to the large number of adjustable parameters. A quantitative predictive solvent extraction model could help to determine the optimal parameters for solvent extraction flow sheets, but such predictive models are not available yet. The main difficulties for such models are the large deviations from ideal thermodynamic behavior in both the aqueous and organic phases due to high solute concentrations. We constructed a molecular thermodynamic model for the extraction of CoCl₂ from different chloride salts by 0.2 mol L^-1 trioctylmethylammonium chloride in toluene using the OLI mixed-solvent electrolyte (OLI-MSE) framework. This was accomplished by analyzing the water and hydrochloric acid content of the organic phase, measuring the water activity of the system, and using metal complex speciation and solvent extraction data. The full extractant concentration range cannot be modeled by the OLI-MSE framework as this framework lacks a description for reversed micelle formation. Nevertheless, salting effects and the behavior of hydrochloric acid can be accurately described with the presented extraction model, without determining specific Co(II)-salt cation interaction parameters. The resulting model shows that the salting effects originate from indirect salt cation-solvent interactions that influence the availability of water in the aqueous and organic phases.

Cation Effect of Chloride Salting Agents on Transition Metal Ion Hydration and Solvent Extraction by the Basic Extractant Methyltrioctylammonium Chloride

The addition of a nonextractable salt has an important influence on the solvent extraction of metal ions, but the underlying principles are not completely understood yet. However, relating solute hydration mechanisms to solvent extraction equilibria is key to understanding the mechanism of solvent extraction of metal ions as a whole. We have studied the speciation of Co(II), Zn(II), and Cu(II) in aqueous solutions containing different chloride salts to understand their extraction to the basic extractant methyltrioctylammonium chloride (TOMAC). This includes the first speciation profile of Zn(II) in chloride media with the three Zn(II) species [Zn(H₂O)₆]²⁺, [ZnCl₃H₂O]⁻, and [ZnCl₄]^2–. The observed differences in extraction efficiency for a given transition metal ion can be explained by transition metal ion hydration due to ion–solvent interactions, rather than by ion–solute interactions or by differences in speciation. Chloride salting agents bearing a cation with a larger hydration Gibbs free energy reduce the free water content more, resulting in a lower hydration for the transition metal ion. This destabilizes the transition metal chloro complex in the aqueous phase and increases the extraction efficiency. Salting agents with di- and trivalent cations reduce the transition metal chloro complex hydration less than expected, resulting in a lower extraction efficiency. The cations of these salting agents have a very large hydration Gibbs free energy, but the overall hydration of these salts is reduced due to significant salt ion pair formation. The general order of salting-out strength for the extraction of metal ions from chloride salt solutions is Cs⁺ < Rb⁺ < NH₄⁺ ≈ K⁺ < Al³⁺ ≈ Mg²⁺ ≈ Ca²⁺ ≈ Na⁺ < Li⁺. These findings can help in predicting the optimal conditions for metal separation by solvent extraction and also contribute to a broader understanding of the effects of dissolved salts on solutes.

Addition of a nonextractable salt influences the stability and solvent extraction efficiency of metal complexes. Cations of different chloride salts reduce the solution free water content as a function of their increasing hydration energy and decreasing tendency for ion pair formation with chloride anions. These ion−solvent interactions reduce the hydration of metal complexes, increasing their distribution ratios. These effects influence aqueous transition metal complexes more than direct ion−solute interactions and changes in complex speciation.

Ion Pairing in HCl-Water Clusters: From Electronic Structure Investigations to Multiconfigurational Force-Field Development

In the bulk, condensed-phase HCl exists as a dissociated Cl^- ion and a proton that is delocalized over solvating water molecules. However, in the gas phase, HCl is covalent, and even on the introduction of hydrating water molecules, the HCl covalent state dominates small clusters and is relevant at larger clusters including 21 water molecules. Electronic structure calculations (at the MP2 level) and ab initio metadynamics simulations (at the DFT level) have been carried out on HCl-(H₂O)_n clusters with n = 2-22 to investigate distinct solvation environments in clusters from covalent HCl structure, to contact ion pairs and solvent-separated ion pairs. The data were further used to train and validate a multiconfigurational force-field for HCl-water clusters that incorporates covalent HCl states into the MS-EVB3.2 formalism. Additionally, the many-body interaction of the Cl^- ion with water and the excess proton was modeled by the introduction of two geometric three-body terms that incorporates the dominant many-body interaction in an efficient noniterative manner.

C70 Fullerene Cage as a Novel Catalyst for Efficient Proton Transfer Reactions between Small Molecules: A Theoretical study

When acids are supplied with an excess electron (or placed in an Ar or the more polarizable N₂ matrix) in the presence of species such as NH₃, the formation of ion-pairs is a likely outcome. Using density functional theory and first-principles calculations, however, we show that, without supplying an external electron or an electric field, or introducing photo-excitation and -ionization, a single molecule of HCl or HBr in the presence of a single molecule of water inside a C₇₀ fullerene cage is susceptible to cleavage of the σ-bond of the Brønsted-Lowry acid into X⁻ and H⁺ ions, with concomitant transfer of the proton along the reaction coordinate. This leads to the formation of an X⁻···⁺HOH₂ (X = Cl, Br) conjugate acid-base ion-pair, similar to the structure in water of a Zundel ion. This process is unlikely to occur in other fullerene derivatives in the presence of H₂O without significantly affecting the geometry of the carbon cage, suggesting that the interior of C₇₀ is an ideal catalytic platform for proton transfer reactions and the design of related novel materials. By contrast, when a single molecule of HF is reacted with a single molecule of H₂O inside the C₇₀ cage, partial proton transfers from HF to H₂O is an immediate consequence, as recently observed experimentally. The geometrical, energetic, electron density, orbital, optoelectronic and vibrational characteristics supporting these observations are presented. In contrast with the views that have been advanced in several recent studies, we show that the encaged species experiences significant non-covalent interaction with the interior of the cage. We also show that the inability of current experiments to detect many infrared active vibrational bands of the endo species in these systems is likely to be a consequence of the substantial electrostatic screening effect of the cage.

Effects of H+ and OH− on H2O as probed by the 1-propanol probing methodology: differential thermodynamic approach

Two-dimensional characterization map of H⁺ and OH⁻, together with other ions.

Ab initio and density functional theory (DFT) studies on triflic acid with water and protonated water clusters

The structure, stability and infrared spectral signatures of triflic acid (TA) with water clusters (W_n) and protonated water clusters (TAH⁺W_n, n = 1 - 6) were computed using DFT and MP2 methods. Our calculations show that a minimum of three water molecules are necessary to stabilize the dissociated zwitterionic form of TA. It can be seen from the results that there is no significant movement of protons in smaller (n = 1 and 2) and linear (n = 1 - 6) types of water clusters. Further, the geometries of TAW_n clusters first form a neutral pair (NP) to contact ion pair (CIP), then form a solvent separated ion pair (SSIP) in a water hexamer. These findings reveal that proton transfer may take place through NP to CIP and then CIP to SSIP. The calculated binding energies (BEs) of ion pair clusters is always higher than that of NP clusters (i.e., more stable than the NP). Existing excess proton linear chain clusters transfer a proton to adjacent water molecules via a Grotthuss mechanism, whereas the same isomers in the branched motifs do not conduct protons. Examination of geometrical parameters and infrared frequencies reveals hydronium ion (H₃O⁺ also called Eigen cation) formation in both TAW_n and protonated TAW_n clusters. The stability of Eigen water clusters is three times higher than that of other non-Eigen water clusters. Our study shows clearly that formation of ion pairs in TAW_n and TAH⁺W_n clusters greatly favors proton transfer to neighboring water molecules and also enhances the stability of these complexes.

Solvation structure of the halides from x-ray absorption spectroscopy

Three-dimensional models for the aqueous solvation structures of chloride, bromide, and iodide are reported. K-edge extended X-ray absorption fine structure (EXAFS) and Minuit X-ray absorption near edge (MXAN) analyses found well-defined single shell solvation spheres for bromide and iodide. However, dissolved chloride proved structurally distinct, with two solvation shells needed to explain its strikingly different X-ray absorption near edge structure (XANES) spectrum. Final solvation models were as follows: iodide, 8 water molecules at 3.60 ± 0.13 Å and bromide, 8 water molecules at 3.40 ± 0.14 Å, while chloride solvation included 7 water molecules at 3.15 ± 0.10 Å, and a second shell of 7 water molecules at 4.14 ± 0.30 Å. Each of the three derived solvation shells is approximately uniformly disposed about the halides, with no global asymmetry. Time-dependent density functional theory calculations simulating the chloride XANES spectra following from alternative solvation spheres revealed surprising sensitivity of the electronic state to 6-, 7-, or 8-coordination, implying a strongly bounded phase space for the correct structure during an MXAN fit. MXAN analysis further showed that the asymmetric solvation predicted from molecular dynamics simulations using halide polarization can play no significant part in bulk solvation. Classical molecular dynamics used to explore chloride solvation found a 7-water solvation shell at 3.12 (-0.04/+0.3) Å, supporting the experimental result. These experiments provide the first fully three-dimensional structures presenting to atomic resolution the aqueous solvation spheres of the larger halide ions.

Protons and Hydroxide Ions in Aqueous Systems

Understanding the structure and dynamics of water's constituent ions, proton and hydroxide, has been a subject of numerous experimental and theoretical studies over the last century. Besides their obvious importance in acid-base chemistry, these ions play an important role in numerous applications ranging from enzyme catalysis to environmental chemistry. Despite a long history of research, many fundamental issues regarding their properties continue to be an active area of research. Here, we provide a review of the experimental and theoretical advances made in the last several decades in understanding the structure, dynamics, and transport of the proton and hydroxide ions in different aqueous environments, ranging from water clusters to the bulk liquid and its interfaces with hydrophobic surfaces. The propensity of these ions to accumulate at hydrophobic surfaces has been a subject of intense debate, and we highlight the open issues and challenges in this area. Biological applications reviewed include proton transport along the hydration layer of various membranes and through channel proteins, problems that are at the core of cellular bioenergetics.

A THz/FTIR fingerprint of the solvated proton: evidence for Eigen structure and Zundel dynamics

Zundel (orange), Eigen (red) and hydration water (light blue) contributions to the THz/FIR extinction of the solvated proton.

A variable ultra-short-pathlength solution cell for XAFS transmission spectroscopy of light elements

An X-ray absorption fine-structure spectroscopy (XAFS) cell that is suitable for solution-phase studies of the light elements in the series from Na(+) and Ca(2+) is described. This cell has an ultra-short pathlength that can be remotely adjusted using a miniature stepper-motor drive and thereby readily provides transmission pathlengths in the range from submicrometer to several hundred micrometers. The flexibility to vary the pathlength enables acquisition of high-quality XAFS spectra and also allows one to check for potential distortions in the spectra from thickness effects. The primary components are mostly commercially available optical parts. The performance of this device is demonstrated at the Cl K-edge (2.8 keV) for several different aqueous Cl(-) solutions.