Infrared cavity ringdown spectroscopy of acid–water clusters: HCl–H2O, DCl–D2O, and DCl–(D2O)2

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.

Femtochemistry of bimolecular reactions from weakly bound complexes: computational study of the H + H'OD → H'OH + D or HOD + H' exchange reactions

A full-dimensional wavepacket propagation describing the bimolecular exchange reactions H + H'OD → H'OH + D or HOD + H' initiated by photolysis of HCl in the hydrogen-bound complex (HCl)⋯(HOD) is reported. The dynamics of this reaction is carried out with the MCTDH method on an ab initio potential energy surface (PES) of H₃O and the initial state is derived from the ground state wavefunction of the complex obtained by relaxation on its own electronic ground state ab initio PES. The description of the system makes use of polyspherical coordinates parametrizing a set of Radau and Jacobi vectors. The calculated energy- and time-resolved reaction probabilities show, owing to the large collision energies at play stemming from the (almost full) photolysis of HCl, that the repulsion between oxygen in the H'OD molecule and the incoming hydrogen atom is the main feature of the collision and leads to non-reactive scattering. No abstraction reaction products are observed. However, both exchange processes are still observable, with a preference in O-H' bond dissociation over that of O-D. The selectivity is reversed upon vibrational pre-excitation of the O-D stretching mode in the H'OD molecule. It is shown that, after the collision, the hydrogen atom of HCl does most likely not encounter the almost stationary chlorine atom again but we also consider the limit case where the H atom is forced to collide multiple times against H'OD as a result of being pushed back by the Cl atom.

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.

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.

Influence of Spatial Inhomogeneity of Detector Temporal Responses on the Spectral Fidelity in Continuous Wave Cavity Ringdown Spectroscopy

Due to their advantages of having a wide bandwidth, low cost, and being easy to obtain, traditional photodetectors (PDs) are being widely applied in measurements of transient signals. The spatial inhomogeneity of such PD temporal responses was measured directly to account for the PD spatial effect of decay rate due to poor alignment in continuous wave cavity ringdown spectroscopy (CW-CRDS) experiments. Based on the measurements of three PDs (i.e., model 1611 (Newport), model 1811 (Newport), and model PDA10CF-EC (Thorlabs)), all the temporal responses followed a tendency of declining first and then rising, and steady platforms existed for the last two PDs. Moreover, as we expected, the closer the PD center was, the faster the response. On the other hand, the initial shut-off amplitude generally reached a larger value for a faster temporal response. As a result, the spatial effect can strongly influence the spectral line shape and value, which will introduce more errors into the precise measurements of spectral parameters using the CRDS technique if this effect is not considered. The defined effective detection area (EDA) of the PDs, which was close to the active area given by manufacturers, was the key parameter that should be paid more attention by researchers. Therefore, the PD should be aligned perfectly to make sure that the EDA covers the laser spot completely.

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.

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.

Anharmonicity of Vibrational Modes in Hydrogen Chloride-Water Mixtures

A thorough analysis of molecular vibrations in the binary system hydrogen chloride/water is presented considering a set of small mixed and pure clusters. In addition to the conventional normal-mode analysis based on the diagonalization of the Hessian, anharmonic frequencies were obtained from the perturbative VPT2 and PT2-VSCF method using hybrid density functional theory. For all normal modes, potential energy curves were modeled by displacing the atoms from the minimum geometry along the normal mode vectors. Three model potentials, a harmonic potential, a Morse potential, and a fourth order polynomial, were applied to fit these curves. From these data, it was possible not only to characterize distinct vibrations as mainly harmonic, anharmonic, or involving higher order terms but also to extract force constants, k, and anharmonicity constants, x_e. By investigating all different types of intramolecular vibrations including covalent stretching or bending vibrations and intermolecular vibrations such as librations, we could demonstrate that while vibrational frequencies can be obtained applying scaling factors to harmonic results, useful anharmonicity constants cannot be predicted in such a way and the usage of more elaborate vibrational methods is necessary. For each particular type of molecular vibration, we could however determine a relationship between the wavenumber or wavenumber shift and the anharmonicity constant, which allows us to estimate mode dependent anharmonicity constants for larger clusters in the future.

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.

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.

Aggregation-induced chemical reactions: acid dissociation in growing water clusters

Understanding chemical reactivity at ultracold conditions, thus enabling molecular syntheses via interstellar and atmospheric processes, is a key issue in cryochemistry. In particular, acid dissociation and proton transfer reactions are ubiquitous in aqueous microsolvation environments. Here, the full dissociation of a HCl molecule upon stepwise solvation by a small number of water molecules at low temperatures, as relevant to helium nanodroplet isolation (HENDI) spectroscopy, is analyzed in mechanistic detail. It is found that upon successive aggregation of HCl with H(2)O molecules, a series of cyclic heteromolecular structures, up to and including HCl(H(2)O)(3), are initially obtained before a precursor state for dissociation, HCl(H(2)O)(3)···H(2)O, is observed upon addition of a fourth water molecule. The latter partially aggregated structure can be viewed as an "activated species", which readily leads to dissociation of HCl and to the formation of a solvent-shared ion pair, H(3)O(+)(H(2)O)(3)Cl(-). Overall, the process is mostly downhill in potential energy, and, in addition, small remaining barriers are overcome by using kinetic energy released as a result of forming hydrogen bonds due to aggregation. The associated barrier is not ruled by thermal equilibrium but is generated by athermal non-equilibrium dynamics. These "aggregation-induced chemical reactions" are expected to be of broad relevance to chemistry at ultralow temperature much beyond HENDI spectroscopy.

Comparison of IR and UV cavity ring-down spectroscopy detection of transient intermediates: pyrolysis of methyl azide to form methyleneimine

Mid-infrared cavity ring-down spectroscopy (CRDS) has been employed in this study to examine the hydride stretching region of methyl azide and its pyrolysis product methyleneimine. The absorption spectrum of methyl azide over 2835-3085 cm(-1) was recorded, and the integrated absorption cross section was determined. The pyrolysis of methyl azide and subsequent production of methyleneimine was observed at various wavenumbers. Using IR CRDS, we were able to observe vibrational transitions of methyleneimine without interference from the methyl azide precursor. Our previous UV CRDS study showed that electronic transitions of methyleneimine overlapped with those of methyl azide. IR CRDS should thus be useful for the detection of polyatomic transient intermediates without interference from precursors.

Fourier transform infrared spectroscopy and ab initio theory of acid–hydrogen sulfide clusters: H2S–HCl, D2S–DCl and H2S–(HCl)2

The rotationally resolved Fourier transform infrared (FTIR) spectrum of the nu(s) HCl and DCl stretching bands for the hydrogen bonded complex H2S-HCl and its isotopomer D2S-DCl have been observed in a supersonic jet at 0.02 cm(-1) resolution. In the same experimental conditions, two additional bands observed without rotational structure in the HCl range of the dimer have been assigned to the cyclic trimer H2S-(HCl)(2). The multidimensional coupling picture involving the donor stretch mode nu(s) and low frequency intermolecular modes already evidenced in several medium strength hydrogen bonded complexes is beautifully confirmed by the observation of completely separated hot band progressions in the 198 K cell spectrum of both dimers. Based on our anharmonic adiabatic approach for the treatment of the coupled vibrations, absolute vibrational frequencies, diagonal and off-diagonal anharmonicities as well as rovibrational coupling constants obtained from analyses of several 2-D subspaces at MP2 and CCSD(T) level are in excellent agreement with spectroscopic results. In the case of small light complexes, the combination of elevated rotational constants and a negligible contribution of intramolecular vibrational redistribution (IVR) improve the reliability of predissociation lifetime measurements, estimated to 180 ps for H2S-HCl and above 200 ps for D2S-DCl.

Vibrational dynamics of the hydrogen bond in H2S–HF: Fourier-transform-infrared spectra and ab initio theory

The rotationally resolved infrared spectrum of the hydrogen bonded complex H(2)S-HF and of its isotopomer D(2)S-DF in the HF/DF stretching range have been observed in a supersonic jet Fourier-transform infrared (FTIR) experiment and indicate a predissociation lifetime of 130 ps for H(2)S-HF. Complementary spectra taken at a temperature of 190 K in a cell without resolved rotational structure indicate the presence of strong anharmonic couplings between low frequency intermolecular modes and the HF donor stretch mode previously observed in other complexes with heavier acceptor molecules without rotational fine structure. The anharmonic analysis of the hot band progressions and of the rotational data confirm the coupling mechanism. The coupling constants and the absolute frequency of the hydrogen bonded stretch mode are in excellent agreement with theoretical predictions based on adiabatic variational calculations on potential surfaces computed at MP2 and CCSD(T) level. Complementary calculations with a perturbational approach further confirm the coupling model.

Rovibrational and dynamical properties of the hydrogen bonded complex (CH2)2S-HF: A combined free jet, cell, and neon matrix-Fourier transform infrared study

Fourier transform infrared spectra of the nu(s) (HF stretching) band of the (CH(2))(2)S-HF complex have been recorded at 0.1-0.5 cm(-1) resolution in a cooled cell, in a supersonic jet expansion seeded with argon and in a neon matrix at 4.5 K. The combination of controlled temperature effects over a range of 40-250 K and a sophisticated band contour simulation program allows the separation of homogeneous and inhomogeneous contributions and reveals significant anharmonic couplings between intramolecular and intermolecular vibrational modes similar to our previous work on (CH(2))(2)S-DF. The sign of the coupling constants is consistent with the expected strengthening of the hydrogen bond upon vibrational excitation of HF which also explains the observed small variations of the geometrical parameters in the excited state. The analysis of sum and difference combination bands involving nu(s) provides accurate values of intermolecular harmonic frequencies and anharmonicities and a good estimate of the dissociation energy of the complex. Frequencies and coupling parameters derived from gas phase spectra compare well with results from neon matrix experiments. The effective linewidth provides a lower bound for the predissociation lifetime of 10 ps. The comparison between effective linewidths and vibrational densities of states for (CH(2))(2)S-HF and -DF complexes highlights the important role of intramolecular vibrational redistribution in the vibrational dynamics of medium strength hydrogen bonds.

Vibrational dynamics of medium strength hydrogen bonds: Fourier transform infrared spectra and band contour analysis of the DF stretching region of (CH[sub 2])[sub 2]S–DF

Fourier transform infrared spectra of the nu(s) band of the (CH2)(2)S-DF complex have been recorded at 0.1-0.5 cm(-1) resolution in a cooled cell and in a supersonic jet expansion seeded with argon. A sufficient density of (CH(2))(2)S-DF heterodimers is produced by a double injection nozzle device, which limits the possibility of reaction between thiirane and DF before the expansion. The observation of partially resolved PQR branch structures at cell temperatures as high as 252 K indicates relatively small effective line widths, which allow a detailed analysis of the underlying vibrational couplings and of the structural properties of the complex. The analysis of cell and free jet spectra in the temperature range 50-250 K is performed with a software package for the simulation and fitting of multiple hot band progressions in asymmetric rotors. The analysis reveals that the three low frequency hydrogen-bond modes are strongly coupled to the DF stretch with anharmonic coupling constants, which indicates a strengthening of the hydrogen bond upon vibrational excitation of DF. Rovibrational parameters and a reliable upper bound for the homogeneous line width have been extracted.

Experimental isotherms of HCl on H[sub 2]O ice under stratospheric conditions: Connections between bulk and interfacial thermodynamics

The adsorption of HCl on the surface of H(2)O ice has been measured at temperatures and pressures relevant to the upper troposphere and lower stratosphere. The measured HCl surface coverage is found to be at least 100 times lower than currently assumed in models of chlorine catalyzed ozone destruction in cold regions of the upper atmosphere. Measurements were conducted in a closed system by simultaneous application of surface spectroscopy and gas phase mass spectrometry to fully characterize vapor/solid equilibrium. Surface adsorption is clearly distinguished from bulk liquid or solid phases. From 180 to 200 K, submonolayer adsorption of HCl is well described by a Bragg-Williams modified Langmuir model which includes the dissociation of HCl into H(+) and Cl(-) ions. Furthermore, adsorption is consistent with two distinct states on the ice substrate, one in which the ions only weakly adsorb on separate sites, and another where the ions adsorb as an H(+)-Cl(-) pair on a single site with adsorption energy comparable to the bulk trihydrate. The number of substrate H(2)O molecules per adsorption site is also consistent with the stoichiometry of bulk hydrates under these conditions. The ionic states exist in equilibrium, and the total adsorption energy is a function of the relative population of both states. These observations and model provide a quantitative connection between the thermodynamics of the bulk and interfacial phases of HCl/H(2)O, and represent a consistent physicochemical model of the equilibrium system.