Demystifying the Diffuse Vibrational Spectrum of Aqueous Protons Through Cold Cluster Spectroscopy

Correlated rotational alignment spectroscopy: a new tool for high-resolution spectroscopy and the analysis of heterogeneous samples

Correlated rotational alignment spectroscopy correlates observables of ultrafast gas-phase spectroscopy with high-resolution, broad-band rotational Raman spectra. This article reviews the measurement principle of CRASY, existing implementations for mass-correlated measurements, and the potential for future developments. New spectroscopic capabilities are discussed in detail: signals for individual sample components can be separated even in highly heterogeneous samples. Isotopologue rotational spectra can be observed at natural isotope abundance. Fragmentation channels are readily assigned in molecular and cluster mass spectra. And finally, rotational Raman spectra can be measured with sub-MHz resolution, an improvement of several orders-of-magnitude as compared to preceding experiments.

Observation of the Infrared-Induced Structural Change in the Microscopic Hydrogen Bond Network of Phenol-Methanol Cluster Cations in a Cold-Ion Trap

To gain insight into microscopic hydrogen bond networks, we measured ultraviolet photodissociation (UVPD) spectra of the phenol-methanol 1:3 cluster cation, [PhOH(MeOH)3]+ trapped in a variable temperature ion trap. At low temperatures, an isomer with a ring-type hydrogen bond structure dominates, whereas at higher temperatures the chain-type isomers become dominant due to the flexibility of their hydrogen bond structures. We also found a clear temperature dependence of the spectral features, such as band position and width. In addition to the above measurement, we observed the infrared (IR) induced isomerization of [PhOH(MeOH)3]+ to study the dynamical aspects of hydrogen bond networks. We succeeded in observing IR-induced isomerization from the ring to chain forms of [PhOH(MeOH)3]+ at low temperature. The isomerization was clearly identified as a change in the UVPD spectra. The time evolution of the UVPD spectra after IR excitation indicated that the IR-induced isomerization occurs within a nanosecond. The chain-type isomers produced by the IR-induced isomerization are then converted back to the ring-type form by collisions with cold He buffer gas in the trap. This backward isomerization proceeds with a time constant of 100 μs under our experimental conditions. In this study, we evaluated the temperatures of the chain isomers during the backward isomerization on the basis of the spectral features.

The Near-Sightedness of Many-Body Interactions in Anharmonic Vibrational Couplings

Couplings between vibrational motions are driven by electronic interactions, and these couplings carry special significance in vibrational energy transfer, multidimensional spectroscopy experiments, and simulations of vibrational spectra. In this investigation, the many-body contributions to these couplings are analyzed computationally in the context of clathrate-like alkali metal cation hydrates, including Cs+(H2O)20, Rb+(H2O)20, and K+(H2O)20, using both analytic and quantum-chemistry potential energy surfaces. Although the harmonic spectra and one-dimensional anharmonic spectra depend strongly on these many-body interactions, the mode-pair couplings were, perhaps surprisingly, found to be dominated by one-body effects, even in cases of couplings to low-frequency modes that involved the motion of multiple water molecules. The origin of this effect was traced mainly to geometric distortion within water monomers and cancellation of many-body effects in differential couplings, and the effect was also shown to be agnostic to the identity of the ion. These outcomes provide new understanding of vibrational couplings and suggest the possibility of improved computational methods for the simulation of infrared and Raman spectra.

Switching of Protonation Sites in Hydrated Nicotine via a Grotthuss Mechanism

The switching of the protonation sites in hydrated nicotine, probed by experimental infrared (IR) spectroscopy and theoretical ab initio calculations, is facilitated via a Grotthuss instead of a bimolecular proton transfer (vehicle) mechanism at the experimental temperature (T = 130 K) as unambiguously confirmed by experiments with deuterated water. In contrast, the bimolecular vehicle mechanism is preferred at higher temperatures (T = 300 K) as determined by theory. The Grotthuss mechanism for the concerted proton transfer results in the production of nicotine's bioactive and addictive pyrrolidine-protonated (Pyrro-H+) protomer with just 5 water molecules. Theoretical analysis suggests that the concerted proton transfer occurs via hydrogen-bonded bridges consisting of a 3 water molecule "core" that connects the pyridine protonated (Pyri-H+) with the pyrrolidine-protonated (Pyrro-H+) protomers. Additional water molecules attached as acceptors to the hydrogen-bonded "core" bridge result in lowering the reaction barrier of the concerted proton transfer down to less than 6 kcal/mol, which is consistent with the experimental observations.

Pitfalls in the n-mode representation of vibrational potentials

Simulations of anharmonic vibrational motion rely on computationally expedient representations of the governing potential energy surface. The n-mode representation (n-MR)-effectively a many-body expansion in the space of molecular vibrations-is a general and efficient approach that is often used for this purpose in vibrational self-consistent field (VSCF) calculations and correlated analogues thereof. In the present analysis, a lack of convergence in many VSCF calculations is shown to originate from negative and unbound potentials at truncated orders of the n-MR expansion. For cases of strong anharmonic coupling between modes, the n-MR can both dip below the true global minimum of the potential surface and lead to effective single-mode potentials in VSCF that do not correspond to bound vibrational problems, even for bound total potentials. The present analysis serves mainly as a pathology report of this issue. Furthermore, this insight into the origin of VSCF non-convergence provides a simple, albeit ad hoc, route to correct the problem by "painting in" the full representation of groups of modes that exhibit these negative potentials at little additional computational cost. Somewhat surprisingly, this approach also reasonably approximates the results of the next-higher n-MR order and identifies groups of modes with particularly strong coupling. The method is shown to identify and correct problematic triples of modes-and restore SCF convergence-in two-mode representations of challenging test systems, including the water dimer and trimer, as well as protonated tropine.

Accelerating Anharmonic Spectroscopy Simulations via Local-Mode, Multilevel Methods

Ab initio computer simulations of anharmonic vibrational spectra provide nuanced insight into the vibrational behavior of molecules and complexes. The computational bottleneck in such simulations, particularly for ab initio potentials, is often the generation of mode-coupling potentials. Focusing specifically on two-mode couplings in this analysis, the combination of a local-mode representation and multilevel methods is demonstrated to be particularly symbiotic. In this approach, a low-level quantum chemistry method is employed to predict the pairwise couplings that should be included at the target level of theory in vibrational self-consistent field (and similar) calculations. Pairs that are excluded by this approach are "recycled" at the low level of theory. Furthermore, because this low-level pre-screening will eventually become the computational bottleneck for sufficiently large chemical systems, distance-based truncation is applied to these low-level predictions without substantive loss of accuracy. This combination is demonstrated to yield sub-wavenumber fidelity with reference vibrational transitions when including only a small fraction of target-level couplings; the overhead of predicting these couplings, particularly when employing distance-based, local-mode cutoffs, is a trivial added cost. This combined approach is assessed on a series of test cases, including ethylene, hexatriene, and the alanine dipeptide. Vibrational self-consistent field (VSCF) spectra were obtained with an RI-MP2/cc-pVTZ potential for the dipeptide, at approximately a 5-fold reduction in computational cost. Considerable optimism for increased accelerations for larger systems and higher-order couplings is also justified, based on this investigation.

Magic Numbers and Stabilities of Photoionized Water Clusters: Computational and Experimental Characterization of the Nanosolvated Hydronium Ion

The stability and distributions of small water clusters generated in a supersonic beam expansion are interrogated by tunable vacuum ultraviolet (VUV) radiation generated at a synchrotron. Time-of-flight mass spectrometry reveals enhanced population of various protonated water clusters (H⁺(H₂O)_n) based upon ionization energy and photoionization distance from source, suggesting there are "magic" numbers below the traditional n = 21 that predominates in the literature. These intensity distributions suggest that VUV threshold photoionization (11.0-11.5 eV) of neutral water clusters close to the nozzle exit leads to a different nonequilibrium state compared to a skimmed molecular beam. This results in the appearance of a new magic number at 14. Metadynamics conformer searches coupled with modern density functional calculations are used to identify the global minimum energy structures of protonated water clusters between n = 2 and 21, as well as the manifold of low-lying metastable minima. New lowest energy structures are reported for the cases of n = 5, 6, 11, 12, 16, and 18, and special stability is identified by several measures. These theoretical results are in agreement with the experiments performed in this work in that n = 14 is shown to exhibit additional stability, based on the computed second-order stabilization energy relative to most cluster sizes, though not to the extent of the well-known n = 21 cluster. Other cluster sizes that show some additional energetic stability are n = 7, 9, 12, 17, and 19. To gain insight into the balance between ion-water and water-water interactions as a function of the cluster size, an analysis of the effective two-body interactions (which sum exactly to the total interaction energy) was performed. This analysis reveals a crossover as a function of cluster size between a water-hydronium-dominated regime for small clusters and a water-water-dominated regime for larger clusters around n = 17.

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.

The primary gas phase hydration shell of hydroxide

The number of water molecules in hydroxide's primary hydration shell has been long debated to be three from the interpretation of experimental data and four from theoretical studies. Here, we provide direct evidence for the presence of a fourth water molecule in hydroxide's primary hydration shell from a combined study based on high-resolution cryogenic experimental photoelectron spectroscopy and high-level quantum chemical computations. Well-defined spectra of OH^-(H₂O)_n clusters (n = 2 to 5) yield accurate electron binding energies, which are, in turn, used as key signatures of the underlying molecular conformations. Although the smaller OH^-(H₂O)₃ and OH^-(H₂O)₄ clusters adopt close-lying conformations with similar electron binding energies that are hard to distinguish, the OH^-(H₂O)₅ cluster clearly has a predominant conformation with a four-coordinated hydroxide binding motif, a finding that unambiguously determines the gas phase coordination number of hydroxide to be four.

State-resolved infrared spectrum of the protonated water dimer: revisiting the characteristic proton transfer doublet peak

The infrared (IR) spectra of protonated water clusters encode precise information on the dynamics and structure of the hydrated proton. However, the strong anharmonic coupling and quantum effects of these elusive species remain puzzling up to the present day. Here, we report unequivocal evidence that the interplay between the proton transfer and the water wagging motions in the protonated water dimer (Zundel ion) giving rise to the characteristic doublet peak is both more complex and more sensitive to subtle energetic changes than previously thought. In particular, hitherto overlooked low-intensity satellite peaks in the experimental spectrum are now unveiled and mechanistically assigned. Our findings rely on the comparison of IR spectra obtained using two highly accurate potential energy surfaces in conjunction with highly accurate state-resolved quantum simulations. We demonstrate that these high-accuracy simulations are important for providing definite assignments of the complex IR signals of fluxional molecules.

Atomistic Simulations for Reactions and Vibrational Spectroscopy in the Era of Machine Learning─Quo Vadis?

Atomistic simulations using accurate energy functions can provide molecular-level insight into functional motions of molecules in the gas and in the condensed phase. This Perspective delineates the present status of the field from the efforts of others and some of our own work and discusses open questions and future prospects. The combination of physics-based long-range representations using multipolar charge distributions and kernel representations for the bonded interactions is shown to provide realistic models for the exploration of the infrared spectroscopy of molecules in solution. For reactions, empirical models connecting dedicated energy functions for the reactant and product states allow statistically meaningful sampling of conformational space whereas machine-learned energy functions are superior in accuracy. The future combination of physics-based models with machine-learning techniques and integration into all-purpose molecular simulation software provides a unique opportunity to bring such dynamics simulations closer to reality.

Stabilizing the Exotic Carbonic Acid by Bisulfate Ion

Carbonic acid is an important species in a variety of fields and has long been regarded to be non-existing in isolated state, as it is thermodynamically favorable to decompose into water and carbon dioxide. In this work, we systematically studied a novel ionic complex [H₂CO₃·HSO₄]^- using density functional theory calculations, molecular dynamics simulations, and topological analysis to investigate if the exotic H₂CO₃ molecule could be stabilized by bisulfate ion, which is a ubiquitous ion in various environments. We found that bisulfate ion could efficiently stabilize all the three conformers of H₂CO₃ and reduce the energy differences of isomers with H₂CO₃ in three different conformations compared to the isolated H₂CO₃ molecule. Calculated isomerization pathways and ab initio molecular dynamics simulations suggest that all the optimized isomers of the complex have good thermal stability and could exist at finite temperatures. We also explored the hydrogen bonding properties in this interesting complex and simulated their harmonic infrared spectra to aid future infrared spectroscopic experiments. This work could be potentially important to understand the fate of carbonic acid in certain complex environments, such as in environments where both sulfuric acid (or rather bisulfate ion) and carbonic acid (or rather carbonic dioxide and water) exist.

Coupled Local-Mode Approach for the Calculation of Vibrational Spectra: Application to Protonated Water Clusters

Spectroscopic studies of protonated water clusters (PWCs) have yielded enormous insights into the fundamental nature of the hydrated proton. Here, we introduce a new coupled local-mode (CLM) approach to calculate PWC OH stretch vibrational spectra. The CLM method combines a sampling of representative configurations from density functional theory (DFT)-based ab initio molecular dynamics (AIMD) simulations with DFT calculations of local-mode vibrational frequencies and couplings. Calculations of inhomogeneous OH stretch vibrational spectra for H⁺(H₂O)₄ and H⁺(H₂O)₂₁ agree well with experiment and higher-level calculations, and decompositions of the calculated spectra in terms of the coupled modes aids in the interpretation of the spectra. This observation is consistent with the idea that capturing anharmonicity and coupling is as important to accuracy as the underlying level of electronic structure theory. The CLM calculations can easily discern the configuration that dominates the experimental measurement for H⁺(H₂O)₅, which can adopt several low-energy conformations.