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

Accurate structures and rotational constants of bicyclic monoterpenes at DFT cost by means of the bond-corrected Pisa composite scheme (BPCS)

The structural, conformational, and spectroscopic properties in the gas phase of 20 bicyclic monoterpenes and monoterpenoids have been analyzed by a new accurate, reduced-cost computational strategy. In detail, the revDSD-PBEP86 double-hybrid functional in conjunction with the D3BJ empirical dispersion corrections and a suitable triple-zeta basis set provides accurate geometrical parameters, whence equilibrium rotational constants, which are further improved by proper account of core-valence correlation. Average deviations within 0.1% between computed and experimental rotational constants are reached when taking into account the vibrational corrections obtained by the B3LYP functional in conjunction with a double-zeta basis set in the framework of second-order vibrational perturbation theory. In addition to their intrinsic interest, the studied terpenes further extend the panel of systems for which the proposed strategy has provided accurate results at density functional theory cost. Therefore, a very accurate yet robust and user-friendly tool is now available for systematic investigations of the role of stereo-electronic effects on the properties of large systems of current technological and/or biological interest by experimentally oriented researchers.

Picking up Good Vibrations through Quartic Force Fields and Vibrational Perturbation Theory

Quartic force fields (QFFs) define sparse potential energy surfaces (compared to semiglobal surfaces) that are the cheapest and easiest means of computing anharmonic vibrational frequencies, especially when utilized with second-order vibrational perturbation theory (VPT2). However, flat and shallow potential surfaces are exceedingly difficult for QFFs to treat through a combination of numerical noise in the often numerically computed derivatives and in competing energy factors in the composite energies often utilized to provide high-level spectroscopic predictions. While some of these issues can be alleviated with analytic derivatives, hybrid QFFs, and intelligent choices in coordinate systems, the best practice is for predicting good molecular vibrations via QFFs is to understand what they cannot do, and this manuscript documents such cases where QFFs may fail.

Toward Accurate Quantum Chemical Methods for Molecules of Increasing Dimension: The New Family of Pisa Composite Schemes

The new versions of the Pisa composite scheme introduced in the present paper are based on the careful selection of different quantum chemical models for energies, geometries, and vibrational frequencies, with the aim of maximizing the accuracy of the overall description while retaining a reasonable cost for all the steps. In particular, the computation of accurate electronic energies has been further improved introducing more reliable complete basis set extrapolations and estimation of core-valence correlation, together with improved basis sets for third-row atoms. Furthermore, the reduced-cost frozen natural orbital (FNO) model has been introduced and validated for large molecules. Accurate molecular structures can be obtained avoiding complete basis set extrapolation and evaluating core-valence correlation at the MP2 level. Unfortunately, analytical gradients are not available for the FNO version of the model. Therefore, for large molecules, an accurate reduced-cost alternative is offered by evaluation of valence contributions with a double-hybrid functional in conjunction with the same MP2 contribution for core-valence correlation or by means of a one-parameter approximation. The same double-hybrid functional and basis set are employed to evaluate zero-point energies and partition functions. After the validation of the new models for small systems, a panel of molecular bricks of life has been used to analyze their performances for problems of current fundamental or technological interest. The fully black-box implementation of the computational workflow paves the way toward the accurate yet not prohibitively expensive study of medium- to large-sized molecules also by experimentally oriented researchers.

Rovibrational analysis of AlCO3, OAlO2, and HOAlO2 for possible atmospheric detection

The lack of observational data for the AlO molecule in the mesosphere/lower thermosphere may be due to ablated aluminum reacting quickly to form other species. Previously proposed reaction pathways show that aluminum could be ablated in the atmosphere from meteoritic activity, but there currently exist very limited spectroscopic data on the intermediates in these reactions, limiting the possible detection of said molecules. As such, rovibrational spectroscopic data are computed herein using quartic force field methodology at four different levels of theory for the neutral intermediates AlCO3, OAlO2, and HOAlO2. Each molecule exhibits multiple vibrational modes with large vibrational transition intensities. For instance, the C-O stretch (ν1) in AlCO3 has a harmonic intensity of 536 km mol-1, the Al-O stretch (ν2) in OAlO2 has an intensity of 678 km mol-1, and the out-of-plane torsion (ν9) in HOAlO2 has an intensity of 158 km mol-1. All three molecules have exceptionally large dipole moments of 6.27, 4.21, and 5.04 D, respectively. These properties indicate that all three molecules are good candidates for potential atmospheric observation utilizing vibrational and/or rotational spectroscopic techniques.

Assigning the CH stretch overtone spectrum of benzene and naphthalene with extension to anthracene and tetracene using 2- and 3-quanta anharmonic quantum chemical computations

The CH stretch overtone region (5750-6300 cm-1) of benzene and naphthalene is assigned herein using anharmonic quantum chemical computations, and the trend of how this extends to larger polycyclic aromatic hydrocarbons (PAHs) is established. The assignment of all experimental bands to specific vibrational states is performed for the first time. Resonance polyads and the inclusion of 3-quanta vibrational states are both needed to compute accurate vibrational frequencies with the proper density-of-states to match the experimental band shape. Hundreds of 3-quanta states produce the observed band structure in naphthalene, anthracene, and tetracene, and this number is expected to increase drastically for larger PAHs. The width and shape of the main peak are consistent from naphthalene to anthracene, necessitating further exploration of this trend to confirm whether it is representative of all PAHs in the CH stretch overtone region. Understanding observations of PAH sources in the 1-3 μm region from the NIRSpec instrument aboard JWST requires new computational data, and this study provides a benchmark and foundation for their computation.

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.

Infrared Spectrum of the Pyrene Anion in the CH Stretching Region

In this work, we report the infrared spectrum of the pyrene anion, measured using messenger tagging with up to three Ar atoms. We assign the spectrum using density functional theory and vibrational perturbation theory. We discuss our results in the context of computed and experimental spectra from the literature as well as recent observations from astronomical sources, addressing the question of whether polycyclic aromatic hydrocarbon anions could contribute to the strong infrared emission bands at 3.29 μm from carbon-rich regions of space.

Infrared spectroscopy of the syn-methyl-substituted Criegee intermediate: A combined experimental and theoretical study

An IR-vacuum ultraviolet (VUV) ion-dip spectroscopy method is utilized to examine the IR spectrum of acetaldehyde oxide (CH3CHOO) in the overtone CH stretch (2νCH) spectral region. IR activation creates a depletion of the ground state population that reduces the VUV photoionization signal on the parent mass channel. IR activation of the more stable and populated syn-CH3CHOO conformer results in rapid unimolecular decay to OH + vinoxy products and makes the most significant contribution to the observed spectrum. The resultant IR-VUV ion-dip spectrum of CH3CHOO is similar to that obtained previously for syn-CH3CHOO using IR action spectroscopy with UV laser-induced fluorescence detection of OH products. The prominent IR features at 5984 and 6081 cm-1 are also observed using UV + VUV photoionization of OH products. Complementary theoretical calculations utilizing a general implementation of second-order vibrational perturbation theory provide new insights on the vibrational transitions that give rise to the experimental spectrum in the overtone CH stretch region. The introduction of physically motivated small shifts of the harmonic frequencies yields remarkably improved agreement between experiment and theory in the overtone CH stretch region. The prominent features are assigned as highly mixed states with contributions from two quanta of CH stretch and/or a combination of CH stretch with an overtone in mode 4. The generality of this approach is demonstrated by applying it to three different levels of electronic structure theory/basis sets, all of which provide spectra that are virtually indistinguishable despite showing large deviations prior to introducing the shifts to the harmonic frequencies.

Exploring the Tl 2 H 2 potential energy surface: A comparative analysis with group 13 systems and experiment

Thallium chemistry is experiencing unprecedented importance. Therefore, it is valuable to characterize some of the simplest thallium compounds. Stationary points along the singlet and triplet Tl  2 H  2 potential energy surface have been characterized. Stationary point geometries were optimized with the CCSD(T)/aug-cc-pwCVQZ-PP method. Harmonic vibrational frequencies were computed at the same level of theory while anharmonic vibrational frequencies were computed at the CCSD(T)/aug-cc-pwCVTZ-PP level of theory. Final energetics were obtained with the CCSDT(Q) method. Basis sets up to augmented quintuple-zeta cardinality (aug-cc-pwCV5Z-PP) were employed to obtain energetics in order to extrapolate to the complete basis set limits using the focal point approach. Zero-point vibrational energy corrections were appended to the extrapolated energies in order to determine relative energies at 0 K. It was found that the planar dibridged isomer lies lowest in energy while the linear structure lies highest in energy. The results were compared to other group 13 M  2 H  2 (M = B, Al, Ga, In, and Tl) theoretical studies and some interesting variations are found. With respect to experiment, incompatibilities exist.

Rotational spectroscopy of methyl tert-butyl ether with a new Ka band chirped-pulse Fourier transform microwave spectrometer

Chirped-pulse Fourier transform microwave (CP-FTMW) spectroscopy is a powerful tool for performing broadband gas-phase rotational spectroscopy, and its applications include discovery of new molecules, complex mixture analysis, and exploration of fundamental molecular physics. Here we report the development of a new Ka band (26.5-40 GHz) CP-FTMW spectrometer that is equipped with a pulsed supersonic expansion source and a heated reservoir for low-volatility samples. The spectrometer is built around a 150 W traveling wave tube amplifier and has an instantaneous bandwidth that covers the entire Ka band spectral range. To test the performance of the spectrometer, the rotational spectrum of methyl tert-butyl ether (MTBE), a former gasoline additive and environmental pollutant, has been measured for the first time in this spectral range. Over 1000 spectroscopic transitions have been measured and assigned to the vibrational ground state and a newly-identified torsionally excited state; all transitions were fit using the XIAM program to a root-mean-square deviation of 22 kHz. The spectrum displays internal rotation splitting, nominally forbidden transitions, and an intriguing axis-switching effect between the ground and torsionally excited state that is a consequence of MTBE's extreme near-prolate nature. Finally, the sensitivity of the spectrometer enabled detection of all singly-substituted 13C and 18O isotopologues in natural abundance. This set of isotopic spectra allowed for a partial r0 structure involving the heavy atoms to be derived, resolving a structural discrepancy in the literature between previous microwave and electron diffraction measurements.

Harmonic and anharmonic vibrational computations for biomolecular building blocks: Benchmarking DFT and basis sets by theoretical and experimental IR spectrum of glycine conformers

Advanced vibrational spectroscopic experiments have reached a level of sophistication that can only be matched by numerical simulations in order to provide an unequivocal analysis, a crucial step to understand the structure-function relationship of biomolecules. While density functional theory (DFT) has become the standard method when targeting medium-size or larger systems, the problem of its reliability and accuracy are well-known and have been abundantly documented. To establish a reliable computational protocol, especially when accuracy is critical, a tailored benchmark is usually required. This is generally done over a short list of known candidates, with the basis set often fixed a priori. In this work, we present a systematic study of the performance of DFT-based hybrid and double-hybrid functionals in the prediction of vibrational energies and infrared intensities at the harmonic level and beyond, considering anharmonic effects through vibrational perturbation theory at the second order. The study is performed for the six-lowest energy glycine conformers, utilizing available "state-of-the-art" accurate theoretical and experimental data as reference. Focusing on the most intense fundamental vibrations in the mid-infrared range of glycine conformers, the role of the basis sets is also investigated considering the balance between computational cost and accuracy. Targeting larger systems, a broad range of hybrid schemes with different computational costs is also tested.

Scaling-up VPT2: A feasible route to include anharmonic correction on large molecules

Vibrational analysis plays a crucial role in the investigation of molecular systems. Though methodologies like second-order vibrational perturbation theory (VPT2) have paved the way to more accurate simulations, the computational cost remains a difficult barrier to overcome when the molecular size increases. Building upon recent advances in the identification of resonances, we propose an approach making anharmonic simulations possible for large-size systems, typically unreachable by standard means. This relies on the fact that, often, only portions of the whole spectra are of actual interest. Therefore, the anharmonic corrections can be included selectively on subsets of normal modes directly related to the regions of interest. Starting from the VPT2 equations, we evaluate rigorously and systematically the impact of the truncated anharmonic treatment onto simulations. The limit and feasibility of the reduced-dimensionality approach are detailed, starting on a smaller model system. The methodology is then challenged on the IR absorption and vibrational circular dichroism spectra of an organometallic complex in three different spectral ranges.

Direct Probe of Vibrational Fingerprint and Combination Band Coupling

Vibrational fingerprints and combination bands are a direct measure of couplings that control molecular properties. However, most combination bands possess small transition dipoles. Here we use multiple, ultrafast coherent infrared pulses to resolve vibrational coupling between CH3CN fingerprint modes at 918 and 1039 cm-1 and combination bands in the 2750-6100 cm-1 region via doubly vibrationally enhanced (DOVE) coherent multidimensional spectroscopy (CMDS). This approach provides a direct probe of vibrational coupling between fingerprint modes and near-infrared combination bands of large and small transition dipoles in a molecular system over a large frequency range.

Unbiased Comparison between Theoretical and Experimental Molecular Structures and Properties: Toward an Accurate Reduced-Cost Evaluation of Vibrational Contributions

The tremendous development of hardware and software is constantly increasing the role of quantum chemical (QC) computations in the assignment and interpretation of experimental results. However, an unbiased comparison between theory and experiment requires the proper account of vibrational averaging effects. In particular, high-resolution spectra in the gas phase are now available for molecules containing up to about 50 atoms, which are too large for a brute-force approach with the available QC methods of sufficient accuracy. In the present paper, we introduce hybrid approaches, which allow the accurate evaluation of vibrational averaging effects for molecules of this size beyond the harmonic approximation, with special attention being devoted to rotational constants. After the validation of new tools for relatively small molecules, the β-estradiol hormone and a prototypical molecular motor have been considered to witness the feasibility of accurate computations for large molecules.

Accurate Structures and Rotational Constants of Steroid Hormones at DFT Cost: Androsterone, Testosterone, Estrone, β-Estradiol, and Estriol

A comprehensive analysis of the structural, conformational, and spectroscopic properties in the gas phase has been performed for five prototypical steroid hormones, namely, androsterone, testosterone, estrone, β-estradiol, and estriol. The revDSD-PBEP86 double-hybrid functional in conjunction with the D3BJ empirical dispersion and a suitable triple-ζ basis set provides accurate conformational energies and equilibrium molecular structures, with the latter being further improved by proper account of core-valence correlation. Average deviations within 0.1% between computed and experimental ground state rotational constants are reached when adding to those equilibrium values vibrational corrections obtained at the cost of standard harmonic frequencies thanks to the use of a new computational tool. Together with the intrinsic interest of the studied hormones, the accuracy of the results obtained at DFT cost for molecules containing about 50 atoms paves the way toward the accurate investigations of other flexible bricks of life.

The infrared absorption spectrum of phenylacetylene and its deuterated isotopologue in the mid- to far-IR

Anharmonicity strongly influences the absorption and emission spectra of polycyclic aromatic hydrocarbon (PAH) molecules. Here, IR-UV ion-dip spectroscopy experiments together with detailed anharmonic computations reveal the presence of fundamental, overtone, as well as 2- and 3-quanta combination band transitions in the far- and mid-infrared absorption spectra of phenylacetylene and its singly deuterated isotopologue. Strong absorption features in the 400-900 cm-1 range originate from CH(D) in-plane and out-of-plane wags and bends, as well as bending motions including the C≡C and CH bonds of the acetylene substituent and the aromatic ring. For phenylacetylene, every absorption feature is assigned either directly or indirectly to a single or multiple vibrational mode(s). The measured spectrum is dense, broad, and structureless in many regions but well characterized by computations. Upon deuteration, large isotopic shifts are observed. At frequencies above 1500 cm-1 for d1-phenylacetylene, a one-to-one match is seen when comparing computations and experiments with all features assigned to combination bands and overtones. The C≡C stretch observed in phenylacetylene is not observed in d1-phenylacetylene due to a computed 40-fold drop in intensity. Overall, a careful treatment of anharmonicity that includes 2- and 3-quanta modes is found to be crucial to understand the rich details of the infrared spectrum of phenylacetylene. Based on these results, it can be expected that such an all-inclusive anharmonic treatment will also be key for unraveling the infrared spectra of PAHs in general.

From the Automated Calculation of Potential Energy Surfaces to Accurate Infrared Spectra

Advances in the development of quantum chemical methods and progress in multicore architectures in computer science made the simulation of infrared spectra of isolated molecules competitive with respect to established experimental methods. Although it is mainly the multidimensional potential energy surface that controls the accuracy of these calculations, the subsequent vibrational structure calculations need to be carefully converged in order to yield accurate results. As both aspects need to be considered in a balanced way, we focus on approaches for molecules of up to 12-15 atoms with respect to both parts, which have been automated to some extent so that they can be employed in routine applications. Alternatives to machine learning will be discussed, which appear to be attractive, as long as local regions of the potential energy surface are sufficient. The automatization of these methods is still in its infancy, and the generalization to molecules with large amplitude motions or molecular clusters is far from trivial, but many systems relevant for astrophysical studies are already in reach.

Performance of EOM-CCSD(T)(a)*-Based Quartic Force Fields in Computing Fundamental, Anharmonic Vibrational Frequencies of Molecular Electronically Excited States with Application to the Ã1A″ State of :CCH2 (Vinylidene)

Highly accurate anharmonic vibrational frequencies of electronically excited states are not as easily computed as their ground electronic state counterparts, but recently developed approximate triple excited state methods may be changing that. One emerging excited state method is equation of motion coupled cluster theory at the singles and doubles level with perturbative triples computed via the (a)* formalism, EOMEE-CCSD(T)(a)*. One of the most employed means for the ready computation of vibrational anharmonic frequencies for ground electronic states is second-order vibrational perturbation theory (VPT2), a theory based on quartic force fields (QFFs),fourth-order Taylor series expansions of the potential portion of the internuclear Watson Hamiltonian. The combination of these two is herein benchmarked for its performance for use as a means of computing rovibrational spectra of electronically excited states. Specifically, the EOMEE-CCSD(T)(a)* approach employing a complete basis set extrapolation along with core electron inclusion and relativity (the so-called "CcCR" approach) defining the QFF produces anharmonic fundamental vibrational frequencies within 2.83%, on the average, of reported gas-phase experimentally assigned values for the test set including the A ~ 1 A ″ states of HCF, HCCl, HSiF, HNO, and HPO. However, some states have exceptional accuracy in the fundamentals, most notably for ν2 of A ~ 1 A ″ HCCl in which the CcCR QFF value is within 1.8 cm-1 at 927.9 cm-1 (or 0.2%) of the experiment. Additionally, this approach produces rotational constants to, on the absolute average, within 0.41% of available experimental data, showcasing notable accuracy in the computation of rovibronic spectral data. Furthermore, utilizing a hybrid approach composed of harmonic CcCR force constants along with a set of simple EOMEE-CCSD(T)(a)*/aug-cc-pVQZ QFF cubic and quartic force constants is faster than using pure CcCR and better represents those modes that suffer from numerical instability in the anharmonic portion of the QFF, implying that this so-called "CcCR + QZ" QFF approach may be the best for future applications. Finally, complete, rovibrational spectral data are provided for A ~ 1 A 2 :CCH2, a molecule of potential astrochemical interest, in order to aid in its potential future experimental rovibronic characterization.

Quantum chemistry meets high-resolution spectroscopy for characterizing the molecular bricks of life in the gas-phase

Computation of accurate geometrical structures and spectroscopic properties of large flexible molecules in the gas-phase is tackled at an affordable cost using a general exploration/exploitation strategy. The most distinctive feature of the approach is the careful selection of different quantum chemical models for energies, geometries and vibrational frequencies with the aim of maximizing the accuracy of the overall description while retaining a reasonable cost for all the steps. In particular, a composite wave-function method is used for energies, whereas a double-hybrid functional (with the addition of core-valence correlation) is employed for geometries and harmonic frequencies and a cheaper hybrid functional for anharmonic contributions. A thorough benchmark based on a wide range of prototypical molecular bricks of life shows that the proposed strategy is close to the accuracy of state-of-the-art composite wave-function methods, and is applicable to much larger systems. A freely available web-utility post-processes the geometries optimized by standard electronic structure codes paving the way toward the accurate yet not prohibitively expensive study of medium- to large-sized molecules by experimentally-oriented researchers.

Correlations between the Structures and Spectra of Protonated Water Clusters

Badger's rule-like correlations between OH stretching frequencies and intensities and the OH bond length are used to develop a spectral mapping procedure for studies of pure and protonated water clusters. This approach utilizes the vibrationally averaged OH bond lengths, which were obtained from diffusion Monte Carlo simulations that were performed using the general potential developed by Yu and Bowman. Good agreement is achieved between the spectra obtained using this approach and previously reported spectra for H+(H2O)n clusters, with n = 3, 4, and 5, as well as their perdeuterated analogues. The analysis of the spectra obtained by this spectral mapping approach supports previous work that assigned the spectrum of H+(H2O)6 to a mixture of Eigen and Zundel-like structures. Analysis of the calculated spectra also suggests a reassignment of the frequency of one of the transitions that involves the OH stretching vibration of the OH bonds in the hydronium core in the Eigen-like structure of H+(H2O)6 from 1917 cm-1 to roughly 2100 cm-1. For D+(D2O)6, comparison of the measured spectrum to those obtained by using the spectral mapping approach suggests that the carrier of the measured spectrum is one or more of the isomers of D+(D2O)6 that contain a four-membered ring and two flanking water molecules. While there are several candidate structures, the two flanking water molecules most likely form a chain that is bound to the hydronium core.

First-Principles Anharmonic Infrared and Raman Vibrational Spectra of Materials: Fermi Resonance in Dry Ice

We introduce a computational tool for the quantum-mechanical simulation of anharmonic infrared and Raman vibrational spectra of materials. The approach, implemented in the CRYSTAL software, stems from Taylor's expansion of the potential energy surface (PES) on the basis of normal modes up to cubic and quartic terms. The PES can be sampled with four different numerical schemes at the level of density functional theory (DFT), with local, generalized-gradient, and hybrid density functional approximations. Anharmonic states are obtained by solving Shrödinger's nuclear equation with either the vibrational self-consistent field (VSCF) or vibrational configuration interaction (VCI) methods. Nuclear quantum effects (NQEs) are thus fully accounted for. Infrared intensities are computed numerically through a Berry phase approach or analytically through a coupled-perturbed (CP) approach. Raman intensities are computed analytically via the CP approach. A variety of anharmonic features of vibrational spectra of materials can be simulated, including band shifts, combination bands, overtones, resonances (first-order Fermi, second-order Darling-Dennison), and hot bands. We showcase the effectiveness of the approach on the description of a first-order Fermi resonance (FR) in CO2 dry ice: a challenging test-case given that the FR occurs in the Raman spectrum, requires NQEs, and involves two- and three-mode couplings. Fundamental mechanistic differences with respect to the well-known FR in molecular CO2 are addressed. This application represents the first quantum-mechanical, periodic description of FR in dry ice.

DFT + F12 QFFs for Cost-Effective Rovibrational Spectral Data Predictions of Ground and Excited Electronic States

The quest for faster computation of anharmonic vibrational frequencies of both ground and excited electronic states has led to combining coupled cluster theory harmonic force constants with density functional theory cubic and quartic force constants for defining a quartic force field (QFF) utilized in conjunction with vibrational perturbation theory at second order (VPT2). This work shows that explicitly correlated coupled cluster theory at the singles, doubles, and perturbative triples levels [CCSD(T)-F12] provides accurate anharmonic vibrational frequencies and rotational constants when conjoined with any of B3LYP, CAM-B3LYP, BHandHLYP, PBE0, and ωB97XD for roughly one-quarter of the computational time of the CCSD(T)-F12 QFF alone for our test set. As the number of atoms in the molecule increases, however, the anharmonic terms become a greater portion of the QFF, and the cost comparison improves with HOCO+ and formic acid, requiring less than 15 and 10% of the time, respectively. In electronically excited states, PBE0 produces more consistently accurate results. Additionally, as the size of the molecule and, in turn, QFF increase, the cost savings for utilizing such a hybrid approach for both ground- and excited-state computations grows. As such, these methods are promising for predicting accurate rovibrational spectral properties for electronically excited states. In cases where well-behaved potentials for a small selection of targeted excited states are needed, such an approach should reduce the computational cost compared to that of methods requiring semiglobal potential surfaces or variational treatments of the rovibronic Hamiltonian. Such applications include spectral characterization of comets, exoplanets, or any situation in which gas phase molecules are being excited by UV-vis radiation.

Influence of fourth-order vibrational corrections on semi-experimental (reSE) structures of linear molecules

Semi-experimental structures (reSE) are derived from experimental ground state rotational constants combined with theoretical vibrational corrections. They permit a meaningful comparison with equilibrium structures based on high-level ab initio calculations. Typically, the vibrational corrections are evaluated with second-order vibrational perturbation theory (VPT2). The amount of error introduced by this approximation is generally thought to be small; however, it has not been thoroughly quantified. Herein, we assess the accuracy of theoretical vibrational corrections by extending the treatment to fourth order (VPT4) for a series of small linear molecules. Typical corrections to bond distances are on the order of 10-5 Å. Larger corrections, nearly 0.0002 Å, are obtained to the bond lengths of NCCN and CNCN. A borderline case is CCCO, which will likely require variational computations for a satisfactory answer. Treatment of vibrational effects beyond VPT2 will thus be important when one wishes to know bond distances confidently to four decimal places (10-4 Å). Certain molecules with shallow bending potentials, e.g., HOC+, are not amenable to a VPT2 description and are not improved by VPT4.

Comprehensive quantum chemical analysis of the (ro)vibrational spectrum of thiirane and its deuterated isotopologue

The (ro)vibrational spectra of thiirane, c-C2H4S, and its fully deuterated isotopologue, c-C2D4S, have been studied by means of vibrational configuration interaction theory, VCI, its incremental variant, iVCI, and subsequent variational rovibrational calculations, RVCI, which rely on multidimensional potential energy surfaces of coupled-cluster quality including up to four-mode coupling terms. Accurate geometrical parameters, fundamental vibrational transitions and first overtones, rovibrational spectra and rotational spectroscopic constants have been determined from these calculations and were compared with experimental results whenever available. A number of tentative misassignments in the vibrational spectra could be resolved and most results for the deuterated thiirane are high-level predictions, which may guide experiments to come. Besides this, a new implementation of infrared intensities within the iVCI framework has been tested for the transitions of the title compounds and are compared with results obtained from standard VCI calculations.

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.

Formic acid-methanol complexation vs. esterification: elusive pre-reactive species identified by vibrational spectroscopy

Vibrational spectra of the mixed dimer and the two mixed trimers of methanol and formic acid as well as some of their isotopologues are presented. Out of the eight expected OH stretching fundamentals of these three pre-reactive hydrogen-bonded complexes, the three modes mainly involving an acid OH group bound to the alcohol appear to be missing in the jet-cooled spectra despite a combination of infrared and Raman probing. A possibility of spectral overlap is discussed in the mixed dimer case, but largely discarded. The missing modes correspond to (fractional) concerted elongation of all engaged OH bonds, promoting synchronous degenerate proton transfer between the molecules. One other trimer mode is very tentatively attributed to a broad spectral feature, whereas all OH bonds contacting carbonyl groups can be unambiguously identified by four relatively narrow infrared absorptions. The spectral features are confirmed by vibrational perturbation theory and deviate in a subtle but systematic way from scaled harmonic predictions which were previously validated for the formic acid complex with a more acidic alcohol. Despite being exothermic and exergonic, ester formation can only be detected in the rarefied gas expansions after extended pre-mixture of the gases, which somewhat contrasts the recent microwave spectroscopic evidence of in situ ester formation and in particular the lack of pre-reactive complex signals.

Cryogenic Ion Vibrational Predissociation (CIVP) Spectroscopy of Aryl Cobinamides in the Gas Phase: How Good Are the Calculations for Vitamin B12 Derivatives?

Aryl corrins represent a novel class of designed B12 derivatives with biological properties of "antivitamins B12". In our previous study, we experimentally determined bond strength in a series of aryl-corrins by the threshold collision-induced dissociation experiments (T-CID) and compared the measured bond dissociation energies (BDEs) with those calculated with density functional theory (DFT). We found that the BDEs are modulated by the side chains around the periphery of the corrin unit. Given that aryl cobinamides have many side chains that increase their conformational space and that the question of a specific structure, measured in the gas phase, was important for further evaluation of our T-CID experiment, we proceeded to analyze structural properties of aryl cobinamides using cryogenic ion vibrational predissociation (CIVP) spectroscopy, static DFT, and Born-Oppenheimer molecular dynamic (BOMD) simulations. We found that none of the examined DFT models could reproduce the CIVP spectra convincingly; both "static" DFT calculations and "dynamic" BOMD simulations provide a surprisingly poor representation of the vibrational spectra, specifically of the number, position, and intensity of bands assigned to hydrogen-bonded versus non-hydrogen-bonded NH and OH moieties. We conclude that, for a flexible molecule with ca. 150 atoms, more accurate approaches are needed before definitive conclusions about computed properties, specifically the structure of the ground-state conformer, may be made.

Approximating large-basis coupled-cluster theory vibrational frequencies using focal-point approximations

The focal-point approximation can be used to estimate a high-accuracy, slow quantum chemistry computation by combining several lower-accuracy, faster computations. We examine the performance of focal-point methods by combining second-order Møller-Plesset perturbation theory (MP2) with coupled-cluster singles, doubles, and perturbative triples [CCSD(T)] for the calculation of harmonic frequencies and that of fundamental frequencies using second-order vibrational perturbation theory (VPT2). In contrast to standard CCSD(T), the focal-point CCSD(T) method approaches the complete basis set (CBS) limit with only triple-ζ basis sets for the coupled-cluster portion of the computation. The predicted harmonic and fundamental frequencies were compared with the experimental values for a set of 20 molecules containing up to six atoms. The focal-point method combining CCSD(T)/aug-cc-pV(T + d)Z with CBS-extrapolated MP2 has mean absolute errors vs experiment of only 7.3 cm-1 for the fundamental frequencies, which are essentially the same as the mean absolute error for CCSD(T) extrapolated to the CBS limit using the aug-cc-pV(Q + d)Z and aug-cc-pV(5 + d)Z basis sets. However, for H2O, the focal-point procedure requires only 3% of the computation time as the extrapolated CCSD(T) result, and the cost savings will grow for larger molecules.

Mid-infrared Laser Spectroscopy of Jet-Cooled Formic Acid Trimer: Mode-Dependent Line Broadening in the C-O Stretching Region

Building on recent progress in the vibrational spectroscopy of the formic acid trimer, we present the first high-resolution measurements of the jet-cooled laser absorption spectrum of (HCOOH)3. The spectra of the lowest- and highest-frequency C-O stretching fundamentals are analyzed whereas the third band is not observed, complicated by monomer and dimer absorptions at 1219 cm-1 (8.2 μm). Vibration-rotation parameters are obtained for the band at 1172.31512(68) cm-1 whereas the C-O stretch at 1246.33(5) cm-1 exhibits a significantly larger breadth, allowing only resolution of the coarse PQR structure. Vibrational predissociation can be ruled out, and intramolecular vibrational redistribution mechanisms are discussed, particularly coupling to the concerted proton exchange within the cyclic dimer subunit. Ultimately, the question remains open. The prospects of high-resolution measurements of other trimer bands or isotope substitution experiments, which might assist in revealing the mode-specificity of the underlying broadening mechanisms, are discussed.

Accuracy of Different Electronic Basis Set Families for Anharmonic Molecular Vibrations: A Comprehensive Benchmark Study

In this work, the accuracy and convergence of different electronic basis set families for the computation of anharmonic molecular vibrational spectroscopic calculations are benchmarked. A series of 39 different basis sets from different families following their hierarchy are assessed on VSCF and VSCF-PT2 algorithms with commonly used MP2 and DFT based B3LYP-D potentials for a set of molecular systems. Such an effort has been validated in a previous work ( J. Phys. Chem. A 2020, 124, 9203-9221) with split-valence basis sets for fundamentals and intensities. Here, fundamental transitions, vibrationally excited states, and intensities are compared with the experimental data to estimate the accuracy for a series of Jensen, Dunning, Calendar, Karlsruhe, and Sapporo basis set families. The convergence of basis sets are also compared with the large ANO basis set. Comprehensive statistical error analysis in terms of accuracy and precision was carried out to assess the performance of each basis set. It is observed that the improvement for the calculated harmonic and anharmonic values from the smaller basis sets to the medium (i.e., triple-ξ) is considerable. Beyond this, from medium to large basis sets, the convergence is slow and mostly posits nearly converged values. Basis sets with and without diffuse functions offer characteristically different accuracies and convergence patterns. Finally, recommendations are given on the choice of basis set chosen as black-box which can balance between accuracy and computational time, estimation of the errors, and their selections especially for large molecules.

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.

Assessment of Density Functional Theory Methods for the Structural Prediction of Transition and Post-Transition Metal-Nucleic Acid Complexes

Understanding the structure of metal-nucleic acid systems is important for many applications such as the design of new pharmaceuticals, metal detection platforms, and nanomaterials. Herein, we explore the ability of 20 density functional theory (DFT) functionals to reproduce the crystal structure geometry of transition and post-transition metal-nucleic acid complexes identified in the Protein Data Bank and Cambridge Structural Database. The environmental extremes of the gas phase and implicit water were considered, and analysis focused on the global and inner coordination geometry, including the coordination distances. Although gas-phase calculations were unable to describe the structure of 12 out of the 53 complexes in our test set regardless of the DFT functional considered, accounting for the broader environment through implicit solvation or constraining the model truncation points to crystallographic coordinates generally afforded agreement with the experimental structure, suggesting that functional performance for these systems is likely due to the models rather than the methods. For the remaining 41 complexes, our results show that the reliability of functionals depends on the metal identity, with the magnitude of error varying across the periodic table. Furthermore, minimal changes in the geometries of these metal-nucleic acid complexes occur upon use of the Stuttgart-Dresden effective core potential and/or inclusion of an implicit water environment. The overall top three performing functionals are ωB97X-V, ωB97X-D3(BJ), and MN15, which reliably describe the structure of a broad range of metal-nucleic acid systems. Other suitable functionals include MN15-L, which is a cheaper alternative to MN15, and PBEh-3c, which is commonly used in QM/MM calculations of biomolecules. In fact, these five methods were the only functionals tested to reproduce the coordination sphere of Cu²⁺-containing complexes. For metal-nucleic acid systems that do not contain Cu²⁺, ωB97X and ωB97X-D are also suitable choices. These top-performing methods can be utilized in future investigations of diverse metal-nucleic acid complexes of relevance to biology and material science.

Isotope Effects in the Zundel-Eigen Isomerization of H+(H2O)6

The isomerization pathway between the energetically low-lying Zundel and Eigen isomers of the protonated water hexamer was investigated using high-level ab initio calculations including a treatment of zero-point corrections. On the basis of these calculations, the Zundel-Eigen isomerization was found to proceed through a stable intermediate isomer, which consists of a four-membered ring with two single acceptor water molecules. The inclusion of vibrational zero-point energy is shown to be important for accurately establishing the relative energies of the three relevant isomers involved in the Zundel-Eigen isomerization. Diffusion Monte Carlo calculations including anharmonic vibrational effects show that all three isomers of H⁺(H₂O)₆ and D⁺(D₂O)₆ have well-defined structures. The energetic ordering of the three isomers changes upon deuteration. The implications of these results for the vibrational spectra of H⁺(H₂O)₆ and D⁺(D₂O)₆ are also discussed.

Spectroscopy of Radicals, Clusters, and Transition States Using Slow Electron Velocity-Map Imaging of Cryogenically Cooled Anions

Slow electron velocity-map imaging of cryogenically cooled anions (cryo-SEVI) is a high-resolution variant of anion photoelectron spectroscopy that has been applied with considerable success over the years to the study of radicals, size-selected clusters, and transition states for unimolecular and bimolecular reactions. Cryo-SEVI retains the versatility of conventional anion photoelectron spectroscopy while offering sub-meV resolution, thereby enabling the resolution of vibrational structure in the photoelectron spectra of complex anions. This Feature Article describes recent experiments in our laboratory using cryo-SEVI, including a new research direction in which anions are vibrationally pre-excited with an infrared laser pulse prior to photodetachment.

pbqff: Push-Button Quartic Force Fields

pbqff is an open-source program for fully automating the production of quartic force fields (QFFs) and their corresponding anharmonic spectroscopic data. Rather than being a monolithic piece of code, it consists of several key modules including a generic interface to quantum chemistry codes and, notably, queuing systems; a molecular point group symmetry library; an internal-to-Cartesian coordinate conversion module; a module for the ordinary least-squares fitting of potential energy surfaces; and an improved second-order rotational and vibrational perturbation theory package for asymmetric and symmetric tops that handles type-1 and -2 Fermi resonances, Fermi resonance polyads, and Coriolis resonances. All of these pieces are written in Rust, a modern, safe, and performant programming language that has much to offer for scientific programming. This work introduces pbqff and its surrounding ecosystem, in addition to reporting new anharmonic vibrational data for c-(C)C₃H₂ and describing how the components of pbqff can be leveraged in other projects.

Accurate Structures and Spectroscopic Parameters of Phenylalanine and Tyrosine in the Gas Phase: A Joint Venture of DFT and Composite Wave-Function Methods

A general strategy for the accurate computation of conformational and spectroscopic properties of flexible molecules in the gas phase is applied to two representative proteinogenic amino acids with aromatic side chains, namely, phenylalanine and tyrosine. The main features of all the most stable conformers predicted by this computational strategy closely match those of the species detected in microwave and infrared experiments. Together with their intrinsic interest, the accuracy of the results obtained with reasonable computer times paves the route for accurate investigations of other flexible bricks of life.

Bringing Machine-Learning Enhanced Quantum Chemistry and Microwave Spectroscopy to Conformational Landscape Exploration: the Paradigmatic Case of 4-Fluoro-Threonine

A combined experimental and theoretical study has been carried out on 4-fluoro-threonine, the only naturally occurring fluorinated amino acid. Fluorination of the methyl group significantly increases the conformational complexity with respect to the parent amino acid threonine. The conformational landscape has been characterized in great detail, with special attention given to the inter-conversion pathways between different conformers. This led to the identification of 13 stable low-energy minima. The equilibrium population of so many conformers produces a very complicated and congested rotational spectrum that could be assigned through a strategy that combines several levels of quantum chemical calculations with the principles of machine learning. Twelve conformers out of 13 could be experimentally characterized. The results obtained from the analysis of the intra-molecular interactions can be exploited to accurately model fluorine-substitution effects in biomolecules.

Rotamers of Methanediol: Composite Ab Initio Predictions of Structures, Frequencies, and Rovibrational Constants

Geminal diols are known to be important intermediates in atmospheric ozonolysis and the aerosol cycle. Recently, the simplest member of this class, methanediol, was interrogated in the gas phase with infrared spectroscopy. To aid in future spectroscopic investigations of methanediol, including in the interstellar medium, we report fundamental frequencies and rovibrational constants for the two rotamers of this molecule using ab initio composite methods along with vibrational perturbation theory. Sensitivity of the predictions to the level of theory and the treatment of anharmonic resonances are carefully assessed. The OH stretching harmonic frequencies of both rotamers are particularly sensitive to the level of theory. The CH stretches of the C_s rotamer are sensitive to the treatment of anharmonic resonances with VPT2-based effective Hamiltonian models. Equilibrium bond distances and harmonic frequencies are converged conservatively to within 0.0005 Å and 3 cm^-1, respectively. The effect of tunneling on the rotational constants is investigated with a 2D variational calculation, based on a relaxed hydroxyl torsional potential energy surface. Tunneling is found to be negligible in the lower energy C₂ rotamer but should modify the rotational constants of the C_s rotamer on the order of MHz, giving rise to rotational line splittings of the same order. The rovibrational constants of the C_s rotamer are dominated by hydroxyl torsional effects, and here we see evidence for the breakdown of vibrational perturbation theory.

Benchmark Structures and Conformational Landscapes of Amino Acids in the Gas Phase: A Joint Venture of Machine Learning, Quantum Chemistry, and Rotational Spectroscopy

The accurate characterization of prototypical bricks of life can strongly benefit from the integration of high resolution spectroscopy and quantum mechanical computations. We have selected a number of representative amino acids (glycine, alanine, serine, cysteine, threonine, aspartic acid and asparagine) to validate a new computational setup rooted in quantum-chemical computations of increasing accuracy guided by machine learning tools. Together with low-lying energy minima, the barriers ruling their interconversion are evaluated in order to unravel possible fast relaxation paths. Vibrational and thermal effects are also included in order to estimate relative free energies at the temperature of interest in the experiment. The spectroscopic parameters of all the most stable conformers predicted by this computational strategy, which do not have low-energy relaxation paths available, closely match those of the species detected in microwave experiments. Together with their intrinsic interest, these accurate results represent ideal benchmarks for more approximate methods.

Phosphine Reactivity and Its Implications for Pyrolysis Experiments and Astrochemistry

Despite the importance of phosphorus-bearing molecules for life and their abundance outside Earth, the chemistry of those compounds still is poorly described. The present study investigates phosphine (PH₃) decomposition and formation pathways. The reactions studied include phosphine thermal dissociation, conversion into PO(²Π), PN(¹Σ₊), and reactions in the presence of H₂O⁺. The thermodynamic and rate coefficients of all reactions are calculated in the range of 50-2000 K considering the CCSD(T)/6-311G(3df,3pd)//ωB97xD/6-311G(3df,3pd) electronic structure data. The rate coefficients were calculated by RRKM and semiclassical transition-state theory (SCTST). According to our results, PH₃ is stable to thermal decomposition at T < 100 K and can be formed promptly by a reaction network involving PH(³Σ_-), PO(²Π), and PN(¹Σ₊). In the presence of radiation or ions, PH₃ is readily decomposed. For this reason, it should be mainly associated with dust grains or icy mantles to be observed under the physical conditions prevailing in the interstellar medium (ISM). The intersystem crossing associated with the dissociation of the isomers PON, NPO, and PNO is accessed by multireference methods, and its importance for the gas-phase PH₃ formation/destruction is discussed. Also, the implications of the present outcomes on phosphorus astrochemistry are highlighted.

Vibrationally excited states of 1H- and 2H-1,2,3-triazole isotopologues analyzed by millimeter-wave and high-resolution infrared spectroscopy with approximate state-specific quartic distortion constants

In this work, we present the spectral analysis of 1H- and 2H-1,2,3-triazole vibrationally excited states alongside provisional and practical computational predictions of the excited-state quartic centrifugal distortion constants. The low-energy fundamental vibrational states of 1H-1,2,3-triazole and five of its deuteriated isotopologues ([1-²H]-, [4-²H]-, [5-²H]-, [4,5-²H]-, and [1,4,5-²H]-1H-1,2,3-triazole), as well as those of 2H-1,2,3-triazole and five of its deuteriated isotopologues ([2-²H]-, [4-²H]-, [2,4-²H]-, [4,5-²H]-, and [2,4,5-²H]-2H-1,2,3-triazole), are studied using millimeter-wave spectroscopy in the 130-375 GHz frequency region. The normal and [2-²H]-isotopologues of 2H-1,2,3-triazole are also analyzed using high-resolution infrared spectroscopy, determining the precise energies of three of their low-energy fundamental states. The resulting spectroscopic constants for each of the vibrationally excited states are reported for the first time. Coupled-cluster vibration-rotation interaction constants are compared with each of their experimentally determined values, often showing agreement within 500 kHz. Newly available coupled-cluster predictions of the excited-state quartic centrifugal distortion constants based on fourth-order vibrational perturbation theory are benchmarked using a large number of the 1,2,3-triazole tautomer isotopologues and vibrationally excited states studied.

Toward Accurate yet Effective Computations of Rotational Spectroscopy Parameters for Biomolecule Building Blocks

The interplay of high-resolution rotational spectroscopy and quantum-chemical computations plays an invaluable role in the investigation of biomolecule building blocks in the gas phase. However, quantum-chemical methods suffer from unfavorable scaling with the dimension of the system under consideration. While a complete characterization of flexible systems requires an elaborate multi-step strategy, in this work, we demonstrate that the accuracy obtained by quantum-chemical composite approaches in the prediction of rotational spectroscopy parameters can be approached by a model based on density functional theory. Glycine and serine are employed to demonstrate that, despite its limited cost, such a model is able to predict rotational constants with an accuracy of 0.3% or better, thus paving the way toward the accurate characterization of larger flexible building blocks of biomolecules.

On the vibrations of formic acid predicted from first principles

In this article, we review recent first principles, anharmonic studies on the molecular vibrations of gaseous formic acid in its monomer form. Transitions identified as fundamentals for both cis- and trans form reported in these studies are collected and supported by results from high-resolution experiments. Attention is given to the effect of coordinate coupling on the convergence of the computed vibrational states.

Toward the Observation of the Tin and Lead Analogs of Formaldehyde

Heavy aldehyde and ketone analogues, R₂X═O (X = Si, Ge, Sn, or Pb), differ from their R₂C═O counterparts due to their greater tendency to oligeramize as the X═O bond polarity increases as one goes down the periodic table. To date, H₂Sn═O and H₂Pb═O have eluded experimental detection. Herein we present the most rigorous theoretical study to date on these structures, providing CCSD(T)/pwCVTZ fundamental frequencies computed on CCSD(T)/CBS optimized structures for the H₂X═O (X = Sn, Pb) potential energy surface. The focal point approach is employed to produce the CCSDTQ/CBS relative energies. For the Sn and Pb structures, the carbene-like cis-HXOH was the global minima, with the trans species being less than 0.6 and 1.1 kcal mol^-1 above the cis structures, respectively. The formaldehyde-like H₂X═O structure is in an energy well of at least 34.8 and 25.4 kcal mol^-1 for Sn and Pb, respectively. Our results provide guidance for future work that may detect H₂Sn═O or H₂Pb═O for the first time.

Exploring Expansions of the Potential and Dipole Surfaces Used for Vibrational Perturbation Theory

A scheme for evaluating expansions of the potential and dipole moment surfaces for vibrational perturbation theory is described. The approach is based on numerical differentiation of the Hessian in the coordinates of interest. It is shown that performing these calculations in internal coordinates generates expansions that are transferable among isotopologues of the molecule of interest. Additionally, re-expressing the expansion of the potential in terms of functions of the internal coordinates, for example, cosines of angles or exponential functions of the bond length displacements, provides expansions that can be used for higher-order perturbation theory calculations. The approach is explored and the results are discussed for water, HOD, ammonia, isomers of HNO₃, and halogenated methane.

Helium Droplet Infrared Spectroscopy of the Butyl Radicals

Butyl radicals ( n-, s-, i-,} and tert-butyl) are formed from the pyrolysis of stable precursors (1-pentyl nitrite, 2-methyl-1-butyl nitrite, isopentyl nitrite, and azo- tert-butane, respectively). The radicals are doped into a beam of liquid helium droplets and probed with infrared action spectroscopy from 2700-3125 cm-1, allowing for a low temperature measurement of the CH stretching region. The presence of anharmonic resonance polyads in the 2800-3000 cm-1 region complicates its interpretation. To facilitate spectral assignment, the anharmonic resonances are modeled with two model Hamiltonian approaches that explicitly couple CH stretch fundamentals to HCH bend overtones and combinations: a VPT2+K normal mode model based on CCSD(T) quartic force fields and a semi-empirical local mode model. Both of these computational methods provide generally good agreement with the experimental spectra.

Rotational and Vibrational Spectra of the Pyridyl Radicals: A Coupled-Cluster Study

Pyridyl is a prototypical nitrogen-containing aromatic radical that may be a key intermediate in the formation of nitrogen-containing aromatic molecules under astrophysical conditions. On meteorites, a variety of complex molecules with nitrogen-containing rings have been detected with nonterrestrial isotopic abundances, and larger nitrogen-containing polycyclic aromatic hydrocarbons (PANHs) have been proposed to be responsible for certain unidentified infrared emission bands in the interstellar medium. In this work, the three isomers of pyridyl (2-, 3-, and 4-pyridyl) have been investigated with coupled cluster methods. For each species, structures were optimized at the CCSD(T)/cc-pwCVTZ level of theory and force fields were calculated at the CCSD(T)/ANO0 level of theory. Second-order vibrational perturbation theory (VPT2) was used to derive anharmonic vibrational frequencies and vibrationally corrected rotational constants, and resonances among vibrational states below 3500 cm^-1 were treated variationally with the VPT2+K method. The results yield a complete set of spectroscopic parameters needed to simulate the pure rotational spectrum of each isomer, including electron-spin, spin-spin, and nuclear hyperfine interactions, and the calculated hyperfine parameters agree well with the limited available data from electron paramagnetic resonance spectroscopy. For the handful of experimentally measured vibrational frequencies determined from photoelectron spectroscopy and matrix isolation spectroscopy, the typical agreement is comparable to experimental uncertainty. The predicted parameters for rotational spectroscopy reported here can guide new experimental investigations into the yet-unobserved rotational spectra of these radicals.