Nonadiabatic vibrational dynamics in the HCO2−⋅H2O complex

Can classical mechanics sense conical intersection?

Conical intersection (CI) leads to fast electronic energy transfer. However, Hamm and Stock [Phys. Rev. Lett. 109, 173201 (2012)] showed the existence of a vibrational CI and its role in vibrational energy relaxation. In this paper, we further investigate the vibrational energy relaxation using an isolated model Hamiltonian system of four vibrational modes with two distinctively different timescales (two fast modes and two slow modes). We show that the excitation of the slow modes plays a crucial role in the energy relaxation mechanism. We also analyze the system from a mixed quantum-classical (surface hopping method) and a completely classical point of view. Notably, surface hopping and even classical simulations also capture fast energy relaxation, which is a signature of CI's existence.

Two-Dimensional Infrared Spectroscopy of Isolated Molecular Ions

Two-dimensional infrared (2D IR) spectroscopy of mass-selected, cryogenically cooled molecular ions is presented. Nonlinear response pathways, encoded in the time-domain photodissociation action response of weakly bound N2 messenger tags, were isolated using pulse shaping techniques following excitation with four collinear ultrafast IR pulses. 2D IR spectra of Re(CO)3(CH3CN)3+ ions capture off-diagonal cross-peak bleach signals between the asymmetric and symmetric carbonyl stretching transitions. These cross peaks display intensity variations as a function of pump-probe delay time due to coherent coupling between the vibrational modes. Well-resolved 2D IR features in the congested fingerprint region of protonated caffeine (C8H10N4O2H+) are also reported. Importantly, intense cross-peak signals were observed at 3 ps waiting time, indicating that tag-loss dynamics are not competing with the measured nonlinear signals. These demonstrations pave the way for more precise studies of molecular interactions and dynamics that are not easily obtainable with current condensed-phase methodologies.

Vibrational Dynamics of the Intramolecular H-Bond in Acetylacetone Investigated with Transient and 2D IR Spectroscopy

Acetylacetone (AcAc) has proven to be a fruitful but highly challenging model system for the experimental and computational interrogation of strong intramolecular hydrogen bonds. Key questions remain, however, regarding the identity of the minimum-energy structure of AcAc and the dynamics of intramolecular proton transfer. Here, we investigate the OH/OD stretch and bend regions of the enol tautomer of AcAc and its deuterated isotopologue with transient absorption and 2D IR spectroscopy. The OH bend region reveals a single dominant diagonal transition near 1625 cm^-1 with intense cross peaks to lower-frequency modes, demonstrating highly mixed fingerprint transitions that contain OH bend character. The anharmonic coupling of the OH bend results in a highly elongated OH bend excited-state absorption transition that indicates a large manifold of OH bend overtone/combination bands in the OH stretch region that leads to strong bend-stretch Fermi resonance interactions. The OH and OD stretch regions consist of broad ground-state bleach signals, but there is no clear evidence of ω₂₁ excited-state absorptions due to rapid population relaxation arising from strong intramolecular coupling to bending, fingerprint, and low-frequency H-bond modes. Orientational relaxation dynamics persist for timescales longer than the vibrational lifetimes, with polarization anisotropy components decaying within approximately 2 and 10 periods of the O-O oscillation for the OH and OD stretch, respectively. The significant isotopic dependence of the orientational dynamics is discussed in the context of intramolecular mode coupling, diffusional processes, and contributions from proton/deuteron transfer dynamics.

Water Network Shape-Dependence of Local Interactions with the Microhydrated -NO2- and -CO2- Anionic Head Groups by Cold Ion Vibrational Spectroscopy

We report the structural evolutions of water networks and solvatochromic response of the CH₃NO₂^- radical anion in the OH and CH stretching regions by analysis of the vibrational spectra displayed by cryogenically cooled CH₃NO₂^-·(H₂O)_n=1-6 clusters. The OH stretching bands evolve with a surprisingly large discontinuity at n = 6, which features the emergence of an intense, strongly red-shifted band along with a weaker feature that appears in the region assigned to a free OH fundamental. Very similar behavior is displayed by the perdeuterated carboxylate clusters, RCO₂^-·(H₂O)_n=5-7 (R = CD₃CD₂), indicating that this behavior is a general feature in the microhydration of the triatomic anionic domain and not associated with CH oscillators. Electronic structure calculations trace this behavior to the formation of a "book" isomer of the water hexamer that adopts a configuration in which one of the water molecules resides in an acceptor-acceptor-donor (AAD) (A = acceptor, D = donor) H-bonding site. Excitation of the bound OH in the AAD site explores the local network topology best suited to stabilize an incipient -XO₂H-OH^-(H₂O)₂ intracluster proton-transfer reaction. These systems thus provide particularly clear examples where the network shape controls the potential energy landscape that governs water network-mediated, intracluster proton transfer. The CH stretching bands of the CH₃NO₂^-·(H₂O)_n=1-6 clusters also exhibit strong solvatochromic shifts, but in this case, they smoothly blue-shift with increasing hydration with no discontinuity at n = 6. This behavior is analyzed in the context of the solute-ion polarizability response and partial charge transfer to the water networks.

Exploring the Origins of Spectral Signatures of Strong Hydrogen Bonding in Protonated Water Clusters

The effects of anharmonicity on the spectral features of strong ionic hydrogen bonds are explored through reduced dimensional studies of the couplings between the hydrogen bonding OH and the donor-acceptor OO stretching vibrations in protonated water clusters with 2-4 water molecules. Specifically, this study focuses on how the anharmonicities and couplings in these ions are reflected in the vibrational spectra by exploring the intensities of the transitions to states with excitation in both the OH and the OO stretching vibrations and changes in the frequency of the OO stretching vibration when the OH stretching vibration is excited. These questions are addressed through the application of several approximate treatments that are based on an adiabatic separation of the high-frequency OH and low-frequency OO stretching vibrations as well as low-order expansions of the potential and dipole surfaces. While an adiabatic approximation captures most of the trends found in the spectra and from an analysis of the two-dimensional model, a vibrational Franck-Condon approach fails to capture the intensities of these transitions. Of the terms in the expansion of the dipole moment function, those that are proportional to Δr_OH and Δr_OH² are found to provide the largest contributions to the calculated intensities of the transitions involving excitation of both the OH and the OO stretches. This leads to the conclusion that the intensities of these transitions encode information about the frequency and anharmonicity of the OH stretching vibration and how they are affected by changes in the OO distance. The anharmonicity of the potential also leads to changes in the OO stretching frequency with excitation of the OH stretching vibration. The direction of this change in frequency encodes additional information about the strength of the ionic hydrogen bond.

Two-Dimensional Adiabatic Model for Calculating Progressions Resulting from Stretch-Rock Coupling in Vibrational Spectra of Anion-Water Complexes

Several anion-water dimers feature a distinct progression in the OH stretch region of their vibrational spectra. This progression arises from strong anharmonic couplings between the OH stretch and low-frequency intermolecular modes. In this work, we introduce a two-dimensional adiabatic model accounting explicitly for the water and anion rock degrees of freedom and use it to calculate the vibrational spectra of HCO₂^-·(H₂O) and NO₃^-·(H₂O). The spectra calculated by using this model are in excellent agreement with experiment, in terms of both peak spacings and lengths of the progressions, and represent a substantial improvement over earlier models.

Infrared spectroscopy of CO3 •-(H2O)1,2 and CO4 •-(H2O)1,2

Hydrated molecular anions are present in the atmosphere. Revealing the structure of the microsolvation is key to understanding their chemical properties. The infrared spectra of CO₃ ^•-(H₂O)_1,2 and CO₄ ^•-(H₂O)_1,2 were measured via infrared multiple photon dissociation spectroscopy in both warm and cold environments. Redshifted from the free O-H stretch frequency, broad, structured spectra were observed in the O-H stretching region for all cluster ions, which provide information on the interaction of the hydrogen atoms with the central ion. In the C-O stretching region, the spectra exhibit clear maxima, but dissociation of CO₃ ^•-(H₂O)_1,2 was surprisingly inefficient. While CO₃ ^•-(H₂O)_1,2 and CO₄ ^•-(H₂O) dissociate via loss of water, CO₂ loss is the dominant dissociation channel for CO₄ ^•-(H₂O)₂. The experimental spectra are compared to calculated spectra within the harmonic approximation and from analysis of molecular dynamics simulations. The simulations support the hypothesis that many isomers contribute to the observed spectrum at finite temperatures. The highly fluxional nature of the clusters is the main reason for the spectral broadening, while water-water hydrogen bonding seems to play a minor role in the doubly hydrated species.

Structure and Vibrational Spectra of ArnH+ (n = 2-3)

The structure and vibrational spectra of protonated Ar clusters Ar_nH⁺ (n = 2-3) are studied using potential energy surfaces at the CCSD(T)/aug-cc-pVTZ level and basis set. Ar binding energies, as well as position isomerism in Ar₃H⁺, were investigated. In our previous work, the spectra of Ar₂H⁺ reveal a strong progression of combination bands, which involves the asymmetric Ar-H⁺ stretch with multiple quanta of the symmetric Ar-H⁺ stretch. In this work, insights on the origin of such progression were examined using an adiabatic model. In addition, contributions from mechanical and electrical anharmonicity on the progressions' intensities were also examined. Comparison of the calculated spectrum for the bare and Ar-tagged ions reveals that the reduction of the symmetry group, from D_∞h to either C_∞v or C_2v, results in a richer vibrational structure in the 500-1700 cm^-1 region. When compared with previously reported action spectra (D. C. McDonald III, D. T. Mauney, D. Leicht, J. H. Marks, J. A. Tan, J.-L. Kuo, and M. A. Duncan, J. Chem. Phys., 2016, 145, 231,101), it appears that the position isomers, because of the binding of the weakly bound Ar messenger, are needed to account for the additional bands in the infrared photodissociation spectrum for Ar₃H⁺. These findings demonstrate the active role of the messenger atom in relaxing some of the selection rules for the bare ion's vibrational transitions - resulting in an augmentation of the bands in the action spectrum.

Spectroscopic Signatures of Mode-Dependent Tunnel Splitting in the Iodide-Water Binary Complex

The gas-phase vibrational spectrum of the isolated iodide-water cluster ion (I^-·H₂O), first reported in 1996, presents one of the most difficult, long-standing spectroscopic puzzles involving ion microhydration. Although the spectra of the smaller halides are well described in the context of an asymmetrical ground-state structure in which only one OH group is hydrogen-bonded to the ion, the I^-·H₂O spectrum displays multiplet structures with partially resolved rotational patterns that are additionally influenced by quantum nuclear spin statistics. In this study, this complex behavior is unraveled with a combination of experimental methods, including ion preparation in a temperature-controlled ion trap and spectral simplification through applications of tag-free, two-color IR-IR double-resonance spectroscopy. Analysis of the double-resonance spectra reveals a vibrational ground-state tunneling splitting of about 20 cm^-1, which is on the same order as the spacing between the peaks that comprise the multiplet structure. These findings are further supported by the results obtained from a fully coupled, six-dimensional calculation of the vibrational spectrum. The underlying level structure can then be understood as a consequence of experimentally measurable, vibrational mode-dependent tunneling splittings (which, in the case of the ground vibrational state, is comparable to the rotational energy spacing between levels with K_a = 0 and 1), as well as Fermi resonance interactions. The latter include the hydrogen-bonded OH stretches and combination bands that involve the HOH bend overtones and soft-mode excitations of frustrated translation and rotation displacements of the water molecule relative to the ion. These anharmonic couplings yield closely spaced bands that are activated in the IR by borrowing intensity from the OH stretch fundamentals.

Capturing intrinsic site-dependent spectral signatures and lifetimes of isolated OH oscillators in extended water networks

The extremely broad infrared spectrum of water in the OH stretching region is a manifestation of how profoundly a water molecule is distorted when embedded in its extended hydrogen-bonding network. Many effects contribute to this breadth in solution at room temperature, which raises the question as to what the spectrum of a single OH oscillator would be in the absence of thermal fluctuations and coupling to nearby OH groups. We report the intrinsic spectral responses of isolated OH oscillators embedded in two cold (~20 K), hydrogen-bonded water cages adopted by the Cs⁺·(HDO)(D₂O)₁₉ and D₃O⁺·(HDO)(D₂O)₁₉ clusters. Most OH oscillators yield single, isolated features that occur with linewidths that increase approximately linearly with their redshifts. Oscillators near 3,400 cm^-1, however, occur with a second feature, which indicates that OH stretch excitation of these molecules drives low-frequency, phonon-type motions of the cage. The excited state lifetimes inferred from the broadening are considered in the context of fluctuations in the local electric fields that are available even at low temperature.

Nonadiabatic effects in electronic and nuclear dynamics

Due to their very nature, ultrafast phenomena are often accompanied by the occurrence of nonadiabatic effects. From a theoretical perspective, the treatment of nonadiabatic processes makes it necessary to go beyond the (quasi) static picture provided by the time-independent Schrödinger equation within the Born-Oppenheimer approximation and to find ways to tackle instead the full time-dependent electronic and nuclear quantum problem. In this review, we give an overview of different nonadiabatic processes that manifest themselves in electronic and nuclear dynamics ranging from the nonadiabatic phenomena taking place during tunnel ionization of atoms in strong laser fields to the radiationless relaxation through conical intersections and the nonadiabatic coupling of vibrational modes and discuss the computational approaches that have been developed to describe such phenomena. These methods range from the full solution of the combined nuclear-electronic quantum problem to a hierarchy of semiclassical approaches and even purely classical frameworks. The power of these simulation tools is illustrated by representative applications and the direct confrontation with experimental measurements performed in the National Centre of Competence for Molecular Ultrafast Science and Technology.

An extended E⊗e Jahn-Teller Hamiltonian for large-amplitude motion: Application to vibrational conical intersections in CH3SH and CH3OH

An extended E⊗e Jahn-Teller Hamiltonian is presented for the case where the (slow) nuclear motion extends far from the symmetry point and may be described approximately as motion on a sphere. Rather than the traditional power series expansion in the displacement from the C_3v symmetry point, an expansion in the spherical harmonics is employed. Application is made to the vibrational Jahn-Teller effect in CH₃XH, with X = S, O, where the equilibrium CXH angles are 83° and 72°, respectively. In addition to the symmetry-required conical intersection (CI) at the C_3v symmetry point, ab initio calculations reveal sets of six symmetry-allowed vibrational CIs in each molecule. The CIs for each molecule are arranged differently in the large-amplitude space, and that difference is reflected in the infrared spectra. The CIs in CH₃SH are found in both eclipsed and staggered geometries, whereas those for CH₃OH are found only in the eclipsed geometry near the torsional saddle point. This difference between the two molecules is reflected in the respective high-resolution spectra in the CH stretch fundamental region.

Hidden role of intermolecular proton transfer in the anomalously diffuse vibrational spectrum of a trapped hydronium ion

We report the vibrational spectra of the hydronium and methyl-ammonium ions captured in the C_3v binding pocket of the 18-crown-6 ether ionophore. Although the NH stretching bands of the CH₃NH₃⁺ ion are consistent with harmonic expectations, the OH stretching bands of H₃O⁺ are surprisingly broad, appearing as a diffuse background absorption with little intensity modulation over 800 cm^-1 with an onset ∼400 cm^-1 below the harmonic prediction. This structure persists even when only a single OH group is present in the HD₂O⁺ isotopologue, while the OD stretching region displays a regular progression involving a soft mode at about 85 cm^-1 These results are rationalized in a vibrationally adiabatic (VA) model in which the motion of the H₃O⁺ ion in the crown pocket is strongly coupled with its OH stretches. In this picture, H₃O⁺ resides in the center of the crown in the vibrational zero-point level, while the minima in the VA potentials associated with the excited OH vibrational states are shifted away from the symmetrical configuration displayed by the ground state. Infrared excitation between these strongly H/D isotope-dependent VA potentials then accounts for most of the broadening in the OH stretching manifold. Specifically, low-frequency motions involving concerted motions of the crown scaffold and the H₃O⁺ ion are driven by a Franck-Condon-like mechanism. In essence, vibrational spectroscopy of these systems can be viewed from the perspective of photochemical interconversion between transient, isomeric forms of the complexes corresponding to the initial stage of intermolecular proton transfer.