Near-Infrared Spectroscopy and Anharmonic Theory of Protonated Water Clusters: Higher Elevations in the Hydrogen Bonding Landscape

Hydrogen bond network structures of protonated dimethylamine clusters H+(DMA)n (n = 3-7)

Infrared spectroscopy of protonated dimethylamine clusters, H+(DMA)n, (n = 3-7), and their Ar-tagged clusters was performed in the NH and CH stretching vibrational region to explore their hydrogen bond network structures. A stable isomer search and vibrational spectral simulations of the clusters were also carried out to support the interpretations of the observed spectra. Weakly hydrogen-bonded NH stretching vibrational bands, which are characteristic of cyclic structures of small-sized protonated clusters, are observed in the spectra of the Ar-tagged clusters of n ≥ 5, while only linear chain type structures are suggested for the Ar-tagged clusters of n = 3-4 and the bare clusters of all the sizes. These results demonstrate that the size and temperature dependence of the hydrogen bond network structures of the protonated dimethylamine clusters is analogous to that of protonated monohydric alcohol clusters.

Near-infrared spectroscopy of H3O+⋯Xn (X = Ar, N2, and CO, n = 1-3)

Near-infrared (NIR) spectra of H3O+⋯Xn (X = Ar, N2, and CO, n = 1-3) in the first overtone region of OH-stretching vibrations (4800-7000 cm-1) were measured. Not only OH-stretching overtones but also several combination bands are major features in this region, and assignments of these observed bands are not obvious at a glance. High-precision anharmonic vibrational simulations based on the discrete variable representation approach were performed. The simulated spectra show good agreement with the observed ones and provide firm assignments of the observed bands, except in the case of X = CO, in which higher order vibrational mode couplings seem significant. This agreement demonstrates that the present system can be a benchmark for high precision anharmonic vibrational computations of NIR spectra. Band broadening in the observed spectra becomes remarkable with an increase of the interaction with the solvent molecule (X). The origin of the band broadening is explored by rare gas tagging experiments and anharmonic vibrational simulations of hot bands.

Calculating Vibrational Excited State Absorptions with Excited State Constrained Minimized Energy Surfaces

The modeling and interpretation of vibrational spectra are crucial for studying reaction dynamics using vibrational spectroscopy. Most prior theoretical developments focused on describing fundamental vibrational transitions while fewer developments focused on vibrational excited state absorptions. In this study, we present a new method that uses excited state constrained minimized energy surfaces (CMESs) to describe vibrational excited state absorptions. The excited state CMESs are obtained similarly to the previous ground state CMES development in our group but with additional wave function orthogonality constraints. Using a series of model systems, including the harmonic oscillator, Morse potential, double-well potential, quartic potential, and two-dimensional anharmonic potential, we demonstrate that this new procedure provides good estimations of the transition frequencies for vibrational excited state absorptions. These results are significantly better than those obtained from harmonic approximations using conventional potential energy surfaces, demonstrating the promise of excited state CMES-based methods for calculating vibrational excited state absorptions in real systems.

Fragment imaging in the infrared photodissociation of the Ar-tagged protonated water clusters H3O+-Ar and H+(H2O)2-Ar

Infrared photodissociation of protonated water clusters with an Ar atom, namely H₃O⁺-Ar and H⁺(H₂O)₂-Ar, was investigated by an imaging technique for mass-selected ions, to reveal the intra- and intermolecular vibrational dynamics. The presented system has the advantage of achieving fragment ion images with the cluster size- and mode-selective photoexcitation of each OH stretching vibration. Translational energy distributions of photofragments were obtained from the images upon the excitation of the bound (ν_b) and free (ν_f) OH stretching vibrations. The energy fractions in the translational motion were compared between ν_b^I and ν_f^I in H₃O⁺-Ar or between ν_b^II and ν_f^II in H⁺(H₂O)₂-Ar, where the labels "I" and "II" represent H₃O⁺-Ar and H⁺(H₂O)₂-Ar, respectively. In H₃O⁺-Ar, the ν_f^I excitation exhibited a smaller translational energy than ν_b^I. This result can be explained by the higher vibrational energy of ν_f^I, which enabled it to produce bending (ν₄) excited H₃O⁺ fragments that should be favored according to the energy-gap model. In contrast to H₃O⁺-Ar, the ν_b^II excitation of an Ar-tagged H₂O subunit and the ν_f^II excitation of an untagged H₂O subunit resulted in very similar translational energy distributions in H⁺(H₂O)₂-Ar. The similar energy fractions independent of the excited H₂O subunits suggested that the ν_b^II and ν_f^II excited states relaxed into a common intermediate state, in which the vibrational energy was delocalized within the H₂O-H⁺-H₂O moiety. However, the translational energy distributions for H⁺(H₂O)₂-Ar did not agree with a statistical dissociation model, which implied another aspect of the process, that is, Ar dissociation via incomplete energy randomization in the whole H⁺(H₂O)₂-Ar cluster.

Character of the OH Bend-Stretch Combination Band in the Vibrational Spectra of the "Magic" Number H3O+(H2O)20 and D3O+(D2O)20 Cluster Ions

The fundamental transitions that contribute to the diffuse OH stretching spectrum of water are known to increase in width and intensity with increasing red shift from the free OH frequency. In contrast, the profile of the higher-energy combination band involving the OH stretching and the intramolecular HOH bending modes displays a qualitatively different spectral shape with a much faster falloff on the lower-energy side. We elucidate the molecular origin of this difference by analyzing the shapes of the combination bands in the IR spectra of cryogenically cooled H₃O⁺(H₂O)₂₀ and D₃O⁺(D₂O)₂₀ clusters. The difference in the shapes of the bands is traced to differences in the dependence of their transition dipole matrix elements on the hydrogen-bonding environment. The fact that individual transitions across the combination band envelope have similar intensities makes it a useful way to determine the participation of various sites in extended H-bonding networks.

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.

Using Diffusion Monte Carlo Wave Functions to Analyze the Vibrational Spectra of H7O3+ and H9O4. J Phys Chem A 2021;125:7185-7197. [PMID: 34433268 DOI: 10.1021/acs.jpca.1c05025]

An approach for evaluating spectra from ground state probability amplitudes (GSPA) obtained from diffusion Monte Carlo (DMC) simulations is extended to improve the description of excited state energies and allow for coupling among vibrational excited states. This approach is applied to studies of the protonated water trimer and tetramer, and their deuterated analogs. These ions provide models for solvated hydronium, and analysis of these spectra provides insights into spectral signatures of proton transfer in aqueous environments. In this approach, we obtain a separable set of internal coordinates from the DMC ground state probability amplitude. A basis is then developed from products of the DMC ground state wave function and low-order polynomials in these internal coordinates. This approach provides a compact basis in which the Hamiltonian and dipole moment matrix are evaluated and used to obtain the spectrum. The resulting spectra are in good agreement with experiment and in many cases provide comparable agreement to the results obtained using much larger basis sets. In addition, the compact basis allows for interpretation of the spectral features and how they evolve with cluster size and deuteration.

Weighted-Graph-Theoretic Methods for Many-Body Corrections within ONIOM: Smooth AIMD and the Role of High-Order Many-Body Terms

We present a weighted-graph-theoretic approach to adaptively compute contributions from many-body approximations for smooth and accurate post-Hartree-Fock (pHF) ab initio molecular dynamics (AIMD) of highly fluxional chemical systems. This approach is ONIOM-like, where the full system is treated at a computationally feasible quality of treatment (density functional theory (DFT) for the size of systems considered in this publication), which is then improved through a perturbative correction that captures local many-body interactions up to a certain order within a higher level of theory (post-Hartree-Fock in this publication) described through graph-theoretic techniques. Due to the fluxional and dynamical nature of the systems studied here, these graphical representations evolve during dynamics. As a result, energetic "hops" appear as the graphical representation deforms with the evolution of the chemical and physical properties of the system. In this paper, we introduce dynamically weighted, linear combinations of graphs, where the transition between graphical representations is smoothly achieved by considering a range of neighboring graphical representations at a given instant during dynamics. We compare these trajectories with those obtained from a set of trajectories where the range of local many-body interactions considered is increased, sometimes to the maximum available limit, which yields conservative trajectories as the order of interactions is increased. The weighted-graph approach presents improved dynamics trajectories while only using lower-order many-body interaction terms. The methods are compared by computing dynamical properties through time-correlation functions and structural distribution functions. In all cases, the weighted-graph approach provides accurate results at a lower cost.

How to VPT2: Accurate and Intuitive Simulations of CH Stretching Infrared Spectra Using VPT2+K with Large Effective Hamiltonian Resonance Treatments

This article primarily discusses the utility of vibrational perturbation theory for the prediction of X-H stretching vibrations with particular focus on the specific variant, second-order vibrational perturbation theory with resonances (VPT2+K). It is written as a tutorial, reprinting most important formulas and providing numerous simple examples. It discusses the philosophy and practical considerations behind vibrational simulations with VPT2+K, including but not limited to computational method selection, cost-saving approximations, approaches to evaluating intensity, resonance identification, and effective Hamiltonian structure. Particular attention is given to resonance treatments, beginning with simple Fermi dyads and gradually progressing to arbitrarily large polyads that describe both Fermi and Darling-Dennison resonances. VPT2+K combined with large effective Hamiltonians is shown to be a reliable framework for modeling the complicated CH stretching spectra of alkenes. An error is also corrected in the published analytic formula for the VPT2 transition moment between the vibrational ground state and triply excited states.

Infrared photodissociation spectroscopy and anharmonic vibrational study of the HO4 + molecular ion

Molecular cations of HO₄ ⁺ and DO₄ ⁺ are produced in a supersonic expansion. They are mass-selected, and infrared photodissociation spectra of these species are measured with the aid of argon-tagging. Although previous theoretical studies have modeled these systems as proton-bound dimers of molecular oxygen, infrared spectra have free OH stretching bands, suggesting other isomeric structures. As a consequence, we undertook extensive computational studies. Our conformer search used a composite method based on an economical combination of single- and multi-reference theories. Several conformers were located on the quintet, triplet, and singlet surfaces, spanning in energy of only a few thousand wavenumbers. Most of the singlet and triplet conformers have pronounced multiconfigurational character. Previously unidentified covalent-like structures (H-O-O-O-O) on the singlet and triplet surfaces likely represent the global minima. In our experiments, HO₄ ⁺ is formed in a relatively hot environment, and similar experiments have been shown capable of producing multiple conformers in low-lying electronic states. None of the predicted HO₄ ⁺ isomers can be ruled out a priori based on energetic arguments. We interpret our argon-tagged spectra with Second-Order Vibrational Perturbation Theory with Resonances (VPT2+K). The presence of one or more covalent-like isomers is the only reasonable explanation for the spectral features observed.

Bend-to-Break: Curvilinear Proton Transfer in Phenol-Ammonia Clusters

The electric field experienced by the OH group of phenol embedded in the cluster of ammonia molecules depends on the relative orientation of the ammonia molecules, and a critical field of 236 MV cm^-1 is essential for the transfer of a proton from phenol to the surrounding ammonia cluster. However, exceptions to this rule were observed, which indicates that the projection of the solvent electric field over the O-H bond is not a definite descriptor of the proton transfer reaction. Therefore, a critical electric field is necessary, but it is not a sufficient condition for the proton abstraction. This, in combination with an adequate solvation of the acceptor ammonia molecule in a triple donor motif that energetically favors the proton transfer process, constitutes necessary and sufficient conditions for the spontaneous proton abstraction. The proton transfer process in phenol-(ammonia)_n clusters is statistically favored to occur away from the plane of the phenyl ring and follows a curvilinear path which includes the O-H bond elongation and out-of-plane movement of the proton. Colloquially, this proton transfer can be referred to as a "bend-to-break" process.

Microsolvation in V+(H2O)n Clusters Studied with Selected-Ion Infrared Spectroscopy

Gas-phase ion-molecule clusters of the form V⁺(H₂O)_n (n = 1-30) are produced by laser vaporization in a supersonic expansion. These ions are analyzed and mass-selected with a time-of-flight mass spectrometer and investigated with infrared laser photodissociation spectroscopy. The small clusters (n ≤ 7) are studied with argon tagging, while the larger clusters are studied via the elimination of water molecules. The vibrational spectra for the small clusters include only free O-H stretching vibrations, while larger clusters exhibit redshifted hydrogen bonding vibrations. The spectral patterns reveal that the coordination around V⁺ ions is completed with four water molecules. A symmetric square-planar structure forms for the n = 4 ion, and this becomes the core ion in larger structures. Clusters up to n = 8 have mostly two-dimensional structures, but hydrogen bonding networks evolve to three-dimensional structures in larger clusters. The free O-H vibration of acceptor-acceptor-donor (AAD)-coordinated surface molecules converges to a frequency near that of bulk water by the cluster size of n = 30. However, the splitting of this vibration for AAD- versus AD-coordinated molecules is still different compared to other singly charged or doubly charged cation-water clusters. This indicates that cation identity and charge-site location in the cluster can produce discernable spectral differences for clusters in this size range.

Spectral signatures of proton delocalization in H+(H2O)n=1−4 ions

Vibrational couplings in protonated water clusters are described by harmonic analysis, vibrational perturbation theory (VPT2) and diffusion Monte Carlo (DMC) approaches.