High resolution spectroscopy of HCl–water clusters: IR bands of undissociated and dissociated clusters revisited

Electric nuclear quadrupole coupling reveals dissociation of HCl with a few water molecules

Investigating the dissociation of acids in the presence of a limited number of water molecules is crucial for understanding various elementary chemical processes. In our study, focusing on HCl(H2O)n clusters (where HCl is hydrogen chloride and H2O is water) formed in a cold and isolated jet expansion, we used the nuclear quadrupole coupling tensor obtained through rotational spectroscopy to decipher the nature of the hydrogen-chlorine (H-Cl) chemical bond in a microaqueous environment. For n = 1 to 4, the H-Cl bond is covalent. At n = 5 and 7, the contact ion pair of H3O+Cl- is spontaneously formed within the hydrogen bond networks of book and cube acid-water clusters, respectively.

Why do dipole moments of HCl-water clusters fail to determine acid dissociation?

This paper quantitatively examines why dipole moments of HCl(H2O)n=1-8 cannot serve as the dissociation criterion for acid molecules using the Hirshfeld-I approach. Also, we propose the possible experimental parameter 〈P(HCl)〉, whose statistical average enables the assessment of acid dissociation in mixed clusters. Furthermore, our calculations reveal that a minimum of four water molecules are necessary to dissociate an HCl molecule.

Mechanism of ionic dissociation of HCl in the smallest water clusters

The dissociation of strong acids into water is a fundamental process in chemistry and biology. Determining the minimum number of water molecules that can result in an ionic dissociation of hydrochloric acid (HCl → H+ + Cl-) remains a challenging subject. In this study, the reactions of H2O with HCl(H2O)n-1 (HCl-H2O cluster), i.e., HCl(H2O)n-1 + H2O (n = 3-7), were investigated by using the direct ab initio molecular dynamics (AIMD) method. Direct AIMD calculations were performed to set the collision energy of H2O to zero for all trajectories. For n = 3, no reaction occurred. In contrast, HCl dissociated to H+ + Cl- at n = 4, forming a contact ion pair (cIP) and solvent-separated ion pair (ssIP) as products. The reactions were expressed as HCl(H2O)3 + H2O → H3O+(H2O)2Cl- (ssIP), and HCl(H2O)3 + H2O → H3O+(Cl-)(H2O)2 (cIP). The ion pair (IP) products were dependent on the collision site of H2O relative to HCl(H2O)3. For n = 5-7, both IPs were formed through the reaction between H2O and HCl(H2O)n-1 (n = 5-7). The reaction between HCl and (H2O)4 (HCl + (H2O)4 → HCl(H2O)4) was non-reactive in IP formation. The reaction mechanism was discussed based on the theoretical results.

Noncovalently bound molecular complexes beyond diatom–diatom systems: full-dimensional, fully coupled quantum calculations of rovibrational states

The methodological advances made in recent years have significantly extended the range and dimensionality of noncovalently bound molecular complexes for which full-dimensional quantum calculations of their rovibrational states are feasible.

Dissociation of HCl in water nanoclusters: an energy decomposition analysis perspective

As known, small HCl-water nanoclusters display a particular dissociation behaviour, whereby at least four water molecules are required for the ionic dissociation of HCl. In this work, we examine how intermolecular interactions promote the ionic dissociation of such nanoclusters. To this end, a set of 45 HCl-water nanoclusters with up to four water molecules is introduced. Energy decomposition analysis based on absolutely localized molecular orbitals (ALMO-EDA) is employed in order to study the importance of frozen interaction, dispersion, polarization, and charge-transfer for the dissociation. The vertical ALMO-EDA scheme is applied to HCl-water clusters along a proton-transfer coordinate varying the amount of spectator water molecules. The corresponding ALMO-EDA results show a clear preference for the dissociated cluster only in the presence of four water molecules. Our analysis of adiabatic ALMO-EDA results reveals a push-pull mechanism for the destabilization of the HCl bond based on the synergy between forward and backward charge-transfer.

Infrared spectroscopy of the protonated HCl dimer and trimer

The protonated HCl dimer and trimer complexes were prepared by pulsed discharges in supersonic expansions of helium or argon doped with HCl and hydrogen. The ions were mass selected in a reflectron time-of-flight spectrometer and investigated with photodissociation spectroscopy in the IR and near-IR regions. Anharmonic vibrational frequencies were computed with VPT2 at the MP2/cc-pVTZ level of theory. The Cl-H stretching fundamentals and overtones were measured in addition to stretch-torsion combinations. VPT2 theory at this level confirms the proton-bound structure of the dimer complex and provides a reasonably good description of the anharmonic vibrations in this system. The trimer has a HCl-HClH⁺-ClH structure in which a central chloronium ion is solvated by two HCl molecules via hydrogen bonding. VPT2 reproduces anharmonic frequencies for this system, including several combinations involving core ion Cl-H stretches, but fails to describe the relative band intensities.

DCl-H2O, HCl-D2O, and DCl-D2O Dimers: Inter- and Intramolecular Vibrational States and Frequency Shifts from Fully Coupled Quantum Calculations on a Full-Dimensional Neural Network Potential Energy Surface

We report full-dimensional and fully coupled quantum calculations of the inter- and intramolecular vibrational states of three isotopologues of the hydrogen chloride-water dimer: DCl-H₂O (DH), HCl-D₂O (HD), and DCl-D₂O (DD). The present study extends our recent theoretical investigation of the nine-dimensional (9D) vibrational level structure of the HCl-H₂O (HH) dimer [Liu, Y.; Li, J.; Felker, P. M.; Bačić, Z. Phys. Chem. Chem. Phys. 2021, 23, 7101-7114]. It employs the same accurate 9D permutation invariant polynomial-neural network potential energy surface and the highly efficient bound-state methodology. The objective of this work is to elucidate the isotopologue variations of a range of bound-state properties of the hydrogen chloride-water dimer and compare them to those of the HH dimer. In order to achieve this, for the isotopologues considered, the rigorous 9D quantum calculations performed encompass all intramolecular vibrational fundamentals, and their frequency shifts relative to the isolated monomer values, together with the low-lying intermolecular vibrational states in each of the intramolecular vibrational manifolds of interest. Moreover, for the ground state of each isotopologue, several informative vibrationally averaged intermolecular geometric properties of the dimer are computed, as well as the three rotational constants. The energies of the intermolecular inversion and rock modes, which mainly involve the motions of the water moiety, differ greatly for H₂O and D₂O, but are much less sensitive to whether the hydrogen chloride isotopologue is HCl or DCl. On the other hand, the excitation of the HCl/DCl stretch changes significantly the energies of the water inversion and rock modes. The DCl stretch frequency shift computed in 9D for the DD dimer, -114.91 cm^-1, agrees extremely well with the corresponding experimental value of -115.20 cm^-1 measured by Saykally and co-workers.

Unveiling Zwitterionization of Glycine in the Microhydration Limit

Charge separation under solvation stress conditions is a fundamental process that comes in many forms in doped water clusters. Yet, the mechanism of intramolecular charge separation, where constraints due to the molecular structure might be intricately tied to restricted solvation structures, remains largely unexplored. Microhydrated amino acids are such paradigmatic molecules. Ab initio simulations are carried out at 300 K in the frameworks of metadynamics sampling and thermodynamic integration to map the thermal mechanisms of zwitterionization using Gly(H₂O) _n with n = 4 and 10. In both cases, a similar water-mediated proton transfer chain mechanism is observed; yet, detailed analyses of thermodynamics and kinetics demonstrate that the charge-separated zwitterion is the preferred species only for n = 10 mainly due to kinetic stabilization. Structural analyses disclose that bifurcated H-bonded water bridges, connecting the cationic and anionic sites in the fluctuating microhydration network at room temperature, are enhanced in the transition-state ensemble exclusively for n = 10 and become overwhelmingly abundant in the stable zwitterion. The findings offer potential insights into charge separation under solvation stress conditions beyond the present example.

HCl–H2O dimer: an accurate full-dimensional potential energy surface and fully coupled quantum calculations of intra- and intermolecular vibrational states and frequency shifts

The present work reports a new full-dimensional potential energy surface (PES) of the HCl–H₂O dimer, and the first fully coupled 9D quantum calculations of the intra- and intermolecular vibrational states of the complex, utilizing this PES.

Generating Excess Protons in Microsolvated Acid Clusters under Ambient Conditions: An Issue of Configurational Entropy versus Internal Energy

Acid dissociation, and thus liberation of excess protons in small water droplets, impacts on diverse fields such as interstellar, atmospheric or environmental chemistry. At cryogenic temperatures below 1 K, it is now well established that as few as four water molecules suffice to dissociate the generic strong acid HCl, yet temperature-driven recombination sets in simply upon heating that cluster. Here, the fundamental question is posed of how many more water molecules are required to stabilize a hydrated excess proton at room temperature. Ab initio path integral simulations disclose that not five, but six water molecules are needed at 300 K to allow for HCl dissociation independently from nuclear quantum effects. In order to provide the molecular underpinnings of these observations, the classical and quantum free energy profiles were decomposed along the dissociation coordinate in terms of the corresponding internal energy and entropy profiles. What decides in the end about acid dissociation, and thus ion pair formation, in a specific microsolvated water cluster at room temperature is found to be a fierce competition between classical configurational entropy and internal energy, where the former stabilizes the undissociated state whereas the latter favors dissociation. It is expected that these are generic findings with broad implications on acid-base chemistry depending on temperature in small water assemblies.

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.

Probing the Partial Activation of Water by Open-Shell Interactions, Cl(H2O)1-4

The partial chemical activation of water by reactive radicals was examined computationally for small clusters of chlorine and water, Cl^•(H₂O)_n=1-4. Using an automated isomer-search procedure, dozens of unique, stable structures were computed. Among the resulting structural classes were intact, hydrated-chlorine isomers, as well as hydrogen-abstracted (HCl)(OH)(H₂O)_n-1 configurations. The latter showed increased stability as the degree of hydration increased, until n = 4, where a new class of structures was discovered with a chloride ion bound to an oxidized water network. The electronic structure of these three structural classes was investigated, and spectral signatures of this hydration-based evolution were connected to these electronic properties. An ancillary outcome of this detailed computational analysis, including coupled-cluster benchmarks, was the calibration of cost-effective quantum chemistry methods for future studies of these radical-water complexes.

Acid solvation versus dissociation at "stardust conditions": Reaction sequence matters

Chemical reactions at ultralow temperatures are of fundamental importance to primordial molecular evolution as it occurs on icy mantles of dust nanoparticles or on ultracold water clusters in dense interstellar clouds. As we show, studying reactions in a stepwise manner in ultracold helium nanodroplets by mass-selective infrared (IR) spectroscopy provides an avenue to mimic these "stardust conditions" in the laboratory. In our joint experimental/theoretical study, in which we successively add H₂O molecules to HCl, we disclose a unique IR fingerprint at 1337 cm^-1 that heralds hydronium (H₃O⁺) formation and, thus, acid dissociation generating solvated protons. In stark contrast, no reaction is observed when reversing the sequence by allowing HCl to interact with preformed small embryonic ice-like clusters. Our ab initio simulations demonstrate that not only reaction stoichiometry but also the reaction sequence needs to be explicitly considered to rationalize ultracold chemistry.

Size-Dependent Formation of an Ion Pair in HSO4-(H2O) n: A Molecular Model for Probing the Microsolvation of Acid Dissociation

With a p K_a of 2.0, HSO₄^- is not a strong acid as its dissociation percentage is only ∼10% in a solution of 1 M. However, our ab initio molecular dynamics and density functional theory calculations show that acid dissociation in the hydrate clusters, HSO₄^-(H₂O) _n, is quite facile at moderate sizes. It starts at around n = 12 and is completed by n = 16 when the energetics becomes very favorable. The dissociation explains the significant broadening at n = 16 in the S═O stretching region of the previously reported infrared photodissociation spectra for HSO₄^-(H₂O) _n as the solvation shell is tightened around the sulfate dianion and the proton. More importantly, HSO₄^-(H₂O) _n should provide an ideal model to probe the molecular details involved in an acid dissociation process.

Hydration motifs of ammonium bisulfate clusters of relevance to atmospheric new particle formation

We have analyzed the binding motifs of water bound to a prototypical cluster containing three ammonium cations and two bisulfate anions using mass-selective vibrational spectroscopy and quantum chemical calculations.

Acid Dissociation in HCl-Water Clusters is Temperature Dependent and Cannot be Detected Based on Dipole Moments

The dissociation of acids in aqueous environments at low temperatures in the presence of a limited amount of water is underlying a wealth of processes from atmospheric to interstellar science. For the paradigmatic case of HCl(H_{2}O)_{n} clusters, our extensive ab initio path integral simulations quantify in terms of free energy differences and barriers that n=4 water molecules are indeed required to dissociate HCl at low temperatures. Increasing the temperature, however, reverses the process and thus counteracts dissociation by fluctuation-driven recombination. The size of the electric dipole moment is shown to not correlate with the acid being in its dissociated or molecular state, thus rendering its measurement as a function of n unable to detect the dissociation transition.

Concerted transfer of multiple protons in acid-water clusters: [(HCl)(H2O)]2 and [(HF)(H2O)]4

Molecular dynamics simulations using directly ab initio potentials are carried out for the ionically bonded clusters [(Cl^-)(H₃O⁺)]₂ and [(F^-)(H₃O⁺)]₄ to explore their transitions to the hydrogen-bonded [(HCl)(H₂O)]₂ and [(HF)(H₂O)]₄ structures during the first picosecond of simulation. Both the ionic and the H-bonded structures that are formed are highly symmetric. It is found that proton transfers are concerted in all trajectories for [(Cl^-)(H₃O⁺)]₂. For [(F^-)(H₃O⁺)]₄, the fully concerted mechanism is dominant but partially concerted transfers of two or three protons at the same time also occur. The concerted mechanism also holds for the reverse process of ionization of neutral acid molecules. It is suggested that the high symmetry of the ionic and the H-bonded structures plays a role in the preference for concerted transfers. Possible implications of the results for proton transfers in other systems are discussed.

Perturbation Approach for Computing Infrared Spectra of the Local Mode of Probe Molecules

Linear and two-dimensional infrared (IR) spectroscopy of site-specific probe molecules provides an opportunity to gain a molecular-level understanding of the local hydrogen-bonding network, conformational dynamics, and long-range electrostatic interactions in condensed-phase and biological systems. A challenge in computation is to determine the time-dependent vibrational frequencies that incorporate explicitly both nuclear quantum effects of vibrational motions and an electronic structural representation of the potential energy surface. In this paper, a nuclear quantum vibrational perturbation (QVP) method is described for efficiently determining the instantaneous vibrational frequency of a chromophore in molecular dynamics simulations. Computational efficiency is achieved through the use of (a) discrete variable representation of the vibrational wave functions, (b) a perturbation theory to evaluate the vibrational energy shifts due to solvent dynamic fluctuations, and (c) a combined QM/MM potential for the systems. It was found that first-order perturbation is sufficiently accurate, enabling time-dependent vibrational frequencies to be obtained on the fly in molecular dynamics. The QVP method is illustrated in the mode-specific linear and 2D-IR spectra of the H-Cl stretching frequency in the HCl-water clusters and the carbonyl stretching vibration of acetone in aqueous solution. To further reduce computational cost, a hybrid strategy was proposed, and it was found that the computed vibrational spectral peak position and line shape are in agreement with experimental results. In addition, it was found that anharmonicity is significant in the H-Cl stretching mode, and hydrogen-bonding interactions further enhance anharmonic effects. The present QVP method complements other computational approaches, including path integral-based molecular dynamics, and represents a major improvement over the electrostatics-based spectroscopic mapping procedures.

Insights into acid dissociation of HCl and HBr with internal electric fields

A critical electric field exerted by the solvent on the ionizable group leads to acid dissociation.

Vibrational Predissociation of the HCl-(H2O)3 Tetramer

The vibrational predissociation of the HCl-(H₂O)₃ tetramer, the largest HCl-(H₂O)_n cluster for which HCl is not predicted to be ionized, is reported. This work focuses on the predissociation pathway giving rise to H₂O + HCl-(H₂O)₂ following IR laser excitation of the H-bonded OH stretch fundamental. H₂O fragments are monitored state selectively by 2 + 1 resonance-enhanced multiphoton ionization (REMPI) combined with time-of-flight mass spectrometry (TOF-MS). Velocity map images of H₂O in selected rotational levels are used to determine translational energy distributions from which the internal energy distributions in the pair-correlated cofragments are derived. From the maximum translational energy release, the bond dissociation energy, D₀ = 2400 ± 100 cm^-1, is determined for the investigated channel. The energy distributions in the fragments are broad, encompassing the entire range of allowed states. The importance of cooperative (nonpairwise) interactions is discussed.

Microhydration of LiOH: Insight from electronic decays of core-ionized states

We compute and compare the autoionization spectra of a core-ionized LiOH molecule both in its isolated and microhydrated states. Stepwise microhydration of LiOH leads to gradual elongation of the Li-OH bond length and finally to molecular dissociation. The accompanying changes in the local environment of the OH(-) and Li(+) counterions are reflected in the computed O 1s and Li 1s spectra. The role of solvent water molecules and the counterion in the spectral shape formation is assessed. Electronic decays of the microhydrated LiOH are found to be mostly intermolecular since the majority of the populated final states have at least one outer-valence vacancy outside the initially core-ionized ion, mainly on a neighboring water molecule. The charge delocalization occurs through the intermolecular Coulombic and electron transfer mediated decays. Both mechanisms are highly efficient that is partly attributed to hybridization of molecular orbitals. The computed spectral shapes are sensitive to the counterion separation as well as to the number and arrangement of solvent molecules. These sensitivities can be used for studying the local hydration structure of solvated ions in aqueous solutions.

How Many Water Molecules Does it Take to Dissociate HCl?

The potential energy surfaces of the HCl(H2O)n (n is the number of water molecules) clusters are systematically explored using density functional theory and high-level ab initio computations. On the basis of electronic energies, the number of water molecules needed for HCl dissociation is four as reported by some experimental groups. However, this number is five owing to the inclusion of entropic factors. Wiberg bond indices are calculated and analyzed, and the results provide a quadratic correlation and classification of clusters according to the nondissociated, partially dissociated, and fully dissociated character of the H-Cl bond. Our computations show that if temperature is not controlled during the experiment, the values obtained for the dipole moment (or for any measurable property) are susceptible to change, providing a different picture of the number of water molecules needed for HCl dissociation in a nanoscopic droplet.

Electric dipole moments of nanosolvated acid molecules in water clusters

The electric dipole moments of (H2O)nDCl (n=3-9) clusters have been measured by the beam-deflection method. Reflecting the (dynamical) charge distribution within the system, the dipole moment contributes information about the microscopic structure of nanoscale solvation. The addition of a DCl molecule to a water cluster results in a strongly enhanced susceptibility. There is evidence for a noticeable rise in the dipole moment occurring at n≈5-6. This size is consistent with predictions for the onset of ionic dissociation. Additionally, a molecular-dynamics model suggests that even with a nominally bound impurity an enhanced dipole moment can arise due to the thermal and zero-point motion of the proton and the water molecules. The experimental measurements and the calculations draw attention to the importance of fluctuations in defining the polarity of water-based nanoclusters and generally to the essential role played by motional effects in determining the response of fluxional nanoscale systems under realistic conditions.

Isolating the spectral signature of H3O+ in the smallest droplet of dissociated HCl acid

The smallest droplet of HCl acid, H₃O⁺(H₂O)₃Cl⁻, and its isolated H₃O⁺ infrared signature.

Effects of Zero-Point Delocalization on the Vibrational Frequencies of Mixed HCl and Water Clusters

We demonstrate the significant effect that large-amplitude zero-point vibrational motion can have on the high-frequency fundamental vibrations of molecular clusters, specifically small (HCl)n-(H2O)m clusters. Calculations were conducted on a many-body potential, constructed from a mix of new and previously reported semiempirical and high-level ab initio potentials. Diffusion Monte Carlo simulations were performed to determine ground-state wave functions. Visualization of these wave functions indicates that the clusters exhibit delocalized ground states spanning multiple stationary point geometries. The ground states are best characterized by planar ring configurations, despite the clusters taking nonplanar configurations at their global minima. Vibrational calculations were performed at the global minima and the Diffusion Monte Carlo predicted configurations and also using an approach that spans multiple stationary points along a rectilinear normal-mode reaction path. Significantly better agreement was observed between the calculated vibrational frequencies and experimental peak positions when the delocalized ground state was accounted for.

Alcohol dimers – how much diagonal OH anharmonicity?

The massive infrared attenuation and increased anharmonicity of OH stretching overtones in alcohol dimers is experimentally quantified.