Multilevel Approach for Direct VSCF/VCI MULTIMODE Calculations with Applications to Large “Zundel” Cations

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.

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.

Approaching the "Zundel" Limit: Tuning the Vibrational Coupling in N2H+Ng, Ng = {He, Ne, Ar, Kr, Xe, and Rn}

The diazenylium ion (N₂H⁺) is a ubiquitous ion in dense molecular clouds. This ion is often used as a dense gas tracer in outer space. Most of the previous works on diazenylium ion have focused on the shared-proton stretch band, ν_H⁺. In this work, we have performed reduced-dimensional calculations to investigate the vibrational structure of N₂H⁺Ng, Ng = {He, Ne, Ar, Kr, Xe, and Rn}. We demonstrate a few interesting things about this system. First, the vibrational coupling in N₂H⁺ can be tuned to switch on interesting anharmonic effects such as Fermi resonance or combination bands by tagging it with different noble gases. Second, a comparison of the vibrational spectrum from N₂H⁺He to N₂H⁺Rn shows that the ν_H⁺ can be swept from an "Eigen-like" to a "Zundel-like" limiting case. Anharmonic calculations were performed using a multilevel approach, which utilized the MP2 and CCSD(T) levels of theories. Binding energies for the elimination of Ng in N₂H⁺Ng are also reported.

Anharmonic IR spectra of solvated ammonium and aminium ions: resemblance between water and bisulfate solvations

In this work, we analyze the vibrational spectra of ammonium, methylammonium, and dimethylammonium ions solvated by either water molecules or bisulfate anions using anharmonic vibrational algorithms. Rich and complicated spectral features in the 2700-3200 cm^-1 region of the experimental spectra of these clusters are attributed to originate from strong Fermi resonance between hydrogen-bonded NH stretching fundamentals and NH bending overtones. Additional weaker bands around 2500-2600 cm^-1 in solvated aminium ions are assigned to the combination tones involving the CH-NH (methyl-amino) rocking modes. Furthermore, the qualitative resemblance in band positions and spectral patterns between two-water-solvated and two-bisulfate-solvated cations suggest a common vibrational coupling scheme beneath the two seemingly different micro-solvation environments.

Spectral Signatures of Protonated Noble Gas Clusters of Ne, Ar, Kr, and Xe: From Monomers to Trimers

The structures and spectral features of protonated noble gas clusters are examined using a first principles approach. Protonated noble gas monomers (NgH⁺) and dimers (NgH⁺Ng) have a linear structure, while the protonated noble gas trimers (Ng₃H⁺) can have a T-shaped or linear structure. Successive binding energies for these complexes are calculated at the CCSD(T)/CBS level of theory. Anharmonic simulations for the dimers and trimers unveil interesting spectral features. The symmetric NgH⁺Ng are charactized by a set of progression bands, which involves one quantum of the asymmetric Ng-H⁺ stretch with multiple quanta of the symmetric Ng-H⁺ stretch. Such a spectral signature is very robust and is predicted to be observed in both T-shaped and linear isomers of Ng₃H⁺. Meanwhile, for selected asymmetric NgH⁺Ng’, a Fermi resonance interaction involving the first overtone of the proton bend with the proton stretch is predicted to occur in ArH⁺Kr and XeH⁺Kr.

Quantum vibration perturbation approach with polyatomic probe in simulating infrared spectra

The quantitative prediction of vibrational spectra of chromophore molecules in solution is challenging and numerous methods have been developed. In this work, we present a quantum vibration perturbation (QVP) approach, which is a procedure that combines molecular quantum vibration and molecular dynamics with perturbation theory. In this framework, an initial Newtonian molecular dynamics simulation is performed, followed by a substitution process to embed molecular quantum vibrational wave functions into the trajectory. The instantaneous vibrational frequency shift at each time step is calculated using the Rayleigh-Schrödinger perturbation theory, where the perturbation operator is the difference in the vibrational potential between the reference chromophore and the perturbed chromophore in the environment. Semi-classical statistical mechanics is employed to obtain the spectral lineshape function. We validated our method using HCOOH·nH₂O (n = 1-2) clusters and HCOOH aqueous solution as examples. The QVP method can be employed for rapid prediction of the vibrational spectrum of a specific mode in solution.

Structural and vibrational characterization of HCO+ and Rg-HCO+, Rg = {He, Ne, Ar, Kr, and Xe}

The structures of the formyl ion (HCO⁺) and its rare gas tagged counterparts (Rg-HCO⁺, Rg = He, Ne, Ar, Kr, and Xe) were studied at the coupled-cluster singles, doubles, and perturbative triples [CCSD(T)]/aug-cc-pVTZ level of theory and basis set. A linear structure for these tagged complexes was predicted. The Rg binding energies for Rg-HCO⁺ are also examined at the CCSD(T) level. It was found that the binding interaction increases from He-HCO⁺ to Xe-HCO⁺. A multilevel potential energy surface built at the CCSD(T) and second-order Møller-Plesset perturbation levels of theory were used to study these species' vibrational spectra. By changing the Rg in the first-solvation shell for HCO⁺, the Fermi resonance interaction between the first H⁺ bend overtone and the asymmetric and symmetric H-C-O stretches can be modulated. This Fermi resonance modulation is demonstrated by examining a series of rare gas solvated HCO⁺.

Fermi resonance switching in KrH+Rg and XeH+Rg (Rg = Ne, Ar, Kr, and Xe)

Matrix isolation experiments have been successfully employed to extensively study the infrared spectrum of several proton-bound rare gas complexes. Most of these studies have focused on the spectral signature for the H⁺ stretch (ν₃) and its combination bands with the intermolecular stretch coordinate (ν₁). However, little attention has been paid to the Fermi resonance interaction between the H⁺ stretch (ν₃) and H⁺ bend overtone (2ν₂) in the asymmetric proton-bound rare gas dimers, RgH⁺Rg'. In this work, we have investigated this interaction on KrH⁺Rg and XeH⁺Rg with Rg = (Ne, Ar, Kr, and Xe). A multilevel potential energy surface (PES) was used to simulate the vibrational structure of these complexes. This PES is a dual-level comprising of second-order Møller-Plesset perturbation theory and coupled-cluster singles doubles with perturbative triples [CCSD(T)] levels of ab initio theories. We found that when both the combination bands (nν₁ + ν₃) and bend overtone 2ν₂ compete to borrow intensity from the ν₃ band, the latter wins over the former, which then results in the suppression of the nν₁ + ν₃ bands. The current simulations offer new assignments for the ArH⁺Xe and KrH⁺Xe spectra. Complete basis set (CBS) binding energies for these complexes were also calculated at the CCSD(T)/CBS level.

A theoretical study on the infrared signatures of proton-bound rare gas dimers (Rg-H+-Rg), Rg = {Ne, Ar, Kr, and Xe}

The infrared spectrum of proton-bound rare gas dimers has been extensively studied via matrix isolation spectroscopy. However, little attention has been paid on their spectrum in the gas phase. Most of the Rg₂H⁺ has not been detected outside the matrix environment. Recently, Ar_nH⁺ (n = 3-7) has been first detected in the gas-phase [D. C. McDonald et al., J. Chem. Phys. 145, 231101 (2016)]. In that work, anharmonic theory can reproduce the observed vibrational structure. In this paper, we extend the existing theory to examine the vibrational signatures of Rg₂H⁺, Rg = {Ne, Ar, Kr, and Xe}. The successive binding of Rg to H⁺ was investigated through the calculation of stepwise formation energies. It was found that this binding is anti-cooperative. High-level full-dimensional potential energy surfaces at the CCSD(T)/aug-cc-pVQZ//MP2/aug-cc-pVQZ were constructed and used in the anharmonic calculation via discrete variable representation. We found that the potential coupling between the symmetric and asymmetric Rg-H⁺ stretch (ν₁ and ν₃ respectively) causes a series of bright n₁ν₁ + ν₃ progressions. From Ne₂H⁺ to Xe₂H⁺, an enhancement of intensities for these bands was observed.