A systematic study of molecular vibrational anharmonicity and vibration-rotation interaction by self-consistent-field higher-derivative methods. Linear polyatomic molecules

MARVEL analysis of high-resolution rovibrational spectra of 16O12C18O

Empirical rovibrational energy levels are presented for the third most abundant, asymmetric carbon dioxide isotopologue, 16O12C18O, based on a compiled dataset of experimental rovibrational transitions collected from the literature. The 52 literature sources utilized provide 19,438 measured lines with unique assignments in the wavenumber range of 2-12,676 cm-1. The MARVEL (Measured Active Rotational-Vibrational Energy Levels) protocol, which is built upon the theory of spectroscopic networks, validates the great majority of these transitions and outputs 8786 empirical rovibrational energy levels with an uncertainty estimation based on the experimental uncertainties of the transitions. Issues found in the literature data, such as misassignment of quantum numbers, typographical errors, and misidentifications, are fixed before including them in the final MARVEL dataset and analysis. Comparison of the empirical energy-level data of this study with those in the line lists CDSD-2019 and Ames-2021 shows good overall agreement, significantly better for CDSD-2019; some issues raised by these comparisons are discussed.

Emerging Conformational-Analysis Protocols from the RTCONF55-16K Reaction Thermochemistry Conformational Benchmark Set

RTCONF55-16K is a new, reactive conformational data set based on cost-efficient methods to assess different conformational analysis protocols. Our reference calculations underpinned the accuracy of the CENSO (Grimme et al. J. Phys. Chem. A, 2021, 125, 4039) procedure and resulted in alternative recipes with different cost-accuracy compromises. Our general-purpose and economical protocols (CENSO-light and zero, respectively) were found to be 10-30 times faster than the original algorithm, adding only 0.4-0.7 kcal/mol absolute error to the relative free energy estimates.

Accurate Vibrational and Ro-Vibrational Contributions to the Properties of Large Molecules by a New Engine Employing Curvilinear Internal Coordinates and Vibrational Perturbation Theory to Second Order

The unbiased comparison between theory and experiment requires approaches more sophisticated than the basic harmonic-oscillator rigid-rotor model, for taking into account vibrational averaging effects and ro-vibrational couplings in molecules of increasing size. Second-order vibrational perturbation theory based on curvilinear internal coordinates (ICs) offers a remarkable compromise between accuracy and computational cost, thanks to the reduction of mode-mode couplings with respect to their counterparts based on Cartesian coordinates. Therefore, we have developed, implemented, and validated a general engine employing ICs, which allows the accurate evaluation of vibrational averages and ro-vibrational couplings for molecules containing up to about 50 atoms beyond the harmonic approximation. After validation of the new tool for relatively small molecules, the effectiveness of ICs has been demonstrated for some flexible and/or quite large molecular bricks of life.

MARVEL analysis of high-resolution rovibrational spectra of 13 C 16 O 2

A set of empirical rovibrational energy levels, obtained through the MARVEL (measured active rotational-vibrational energy levels) procedure, is presented for the  13 C  16 O  2 isotopologue of carbon dioxide. This procedure begins with the collection and analysis of experimental rovibrational transitions from the literature, allowing for a comprehensive review of the literature on the high-resolution spectroscopy of  13 C  16 O  2 , which is also presented. A total of 60 sources out of more than 750 checked provided 14,101 uniquely measured and assigned rovibrational transitions in the wavenumber range of 579-13,735 cm  - 1 . This is followed by a weighted least-squares refinement yielding the energy levels of the states involved in the measured transitions. Altogether 6318 empirical rovibrational energies have been determined for  13 C  16 O  2 . Finally, estimates have been given for the uncertainties of the empirical energies, based on the experimental uncertainties of the transitions. The detailed analysis of the lines and the spectroscopic network built from them, as well as the uncertainty estimates, all serve to pinpoint possible errors in the experimental data, such as typos, misassignment of quantum numbers, and misidentifications. Errors found in the literature data were corrected before including them in the final MARVEL dataset and analysis.

Unusual Dynamics and Vibrational Fingerprints of van der Waals Dimers Formed by Linear Molecules and Rare-Gas Atoms

Detailed structural, dynamical, and vibrational analyses have been performed for systems composed of linear triatomic molecules solvated by a single rare-gas atom, He, Ne, or Ar. Among the chromophores of these van der Waals (vdW) dimers, there are four neutral molecules (CO2, CS2, N2O, and OCS) and six molecular cations (HHe2+, HNe2+, HAr2+, HHeNe+, HHeAr+, and HNeAr+), both of apolar and polar nature. Following the exploration of bonding preferences, high-level four-dimensional (4D) potential energy surfaces (PESs) have been developed for 24 vdW dimers, keeping the two intramonomer bond lengths fixed. For these 24 complexes, over 1500 bound vibrational states have been obtained via quasi-variational nuclear-motion computations, employing exact kinetic-energy operators together with the accurate 4D PESs and their 2D/3D cuts. The reduced-dimensional (2D to 4D) dimer models have been compared with full-dimensional (6D) ones in the cases of the neutral CO2·Ar and charged HHe2+·He dimers, corroborating the high accuracy of the 2D to 4D vibrational energies. The reduced-dimensional models suggest that (a) while the equilibrium structures are T-shaped and planar, the effective ground-state structures are nonplanar, (b) certain bound states belong to collinear molecular structures, even when they are not minima, (c) the vdW vibrations are heavily mixed and many states have amplitudes corresponding to both the T-shaped and collinear structures, (d) there are a few dimers, for which even some of the vdW fundamentals lie above the first dissociation limit, and (e) the vdW vibrations are almost fully decoupled from the intramonomer bending motion.

Ab Initio Rovibrational Spectroscopy of the Acetylide Anion

In this work the rovibrational spectrum of the acetylide anion HCC- is investigated using high-level electronic structure methods and variational rovibrational calculations. Using a composite approach the potential energy surface and dipole surface is constructed from explicitly correlated coupled-cluster accounting for corrections due to core-valence correlation, scalar relativistic effects and higher-order excitation effects. Previous approaches for approximating the latter are critically evaluated. Employing the composite potential, accurate spectroscopic parameters determined from variational calculations are presented. In comparison to the few available reference data the present results show excellent agreement with ground state rotational constants within 0.005% of the experimental value. Intensities determined from the variational calculations suggest the bending fundamental transition ν2 around 510 cm-1 to be the best target for detection. The rather weak CD stretching fundamental ν1 in deuterated isotopologues show a second-order resonance with the (0,20,1) state and the consequences are discussed in some detail. The spectroscopic parameters and band intensities provided for a number of vibrational bands in isotopologues of the acetylide anion should facilitate future spectroscopic investigations.

Structure, energetics, and spectroscopy of the chromophores of HHe+n, H2He+n, and He+n clusters and their deuterated isotopologues

The linear molecular ions H₂He⁺, HHe+2, and He+3 are the central units (chromophores) of certain He-solvated complexes of the H₂He+n, HHe+n, and He+n families, respectively. These are complexes which do exist, according to mass-spectrometry studies, up to very high n values. Apparently, for some of the H₂He+n and He+n complexes, the linear symmetric tetratomic H₂He+2 and the diatomic He+2 cations, respectively, may also be the central units. In this study, definitive structures, relative energies, zero-point vibrational energies, and (an)harmonic vibrational fundamentals, and, in some cases, overtones and combination bands, are established mostly for the triatomic chromophores. The study is also extended to the deuterated isotopologues D₂He⁺, DHe+2, and D₂He+2. To facilitate and improve the electronic-structure computations performed, new atom-centered, fixed-exponent, Gaussian-type basis sets called MAX, with X = T(3), Q(4), P(5), and H(6), are designed for the H and He atoms. The focal-point-analysis (FPA) technique is employed to determine definitive relative energies with tight uncertainties for reactions involving the molecular ions. The FPA results determined include the 0 K proton and deuteron affinities of the ⁴He atom, 14 875(9) cm^-1 [177.95(11) kJ mol^-1] and 15 229(8) cm^-1 [182.18(10) kJ mol^-1], respectively, the dissociation energies of the He+2 → He⁺ + He, HHe+2 → HHe⁺ + He, and He+3 → He+2 + He reactions, 19 099(13) cm^-1 [228.48(16) kJ mol^-1], 3948(7) cm^-1 [47.23(8) kJ mol^-1], and 1401(12) cm^-1 [16.76(14) kJ mol^-1], respectively, the dissociation energy of the DHe+2 → DHe⁺ + He reaction, 4033(6) cm^-1 [48.25(7) kJ mol^-1], the isomerization energy between the two linear isomers of the [H, He, He]⁺ system, 3828(40) cm^-1 [45.79(48) kJ mol^-1], and the dissociation energies of the H₂He⁺ → H+2 + He and the H₂He+2 → H₂He⁺ + He reactions, 1789(4) cm^-1 [21.40(5) kJ mol^-1] and 435(6) cm^-1 [5.20(7) kJ mol^-1], respectively. The FPA estimates of the first dissociation energy of D₂He⁺ and D₂He+2 are 1986(4) cm^-1 [23.76(5) kJ mol^-1] and 474(5) cm^-1 [5.67(6) kJ mol^-1], respectively. Determining the vibrational fundamentals of the triatomic chromophores with second-order vibrational perturbation theory (VPT2) and vibrational configuration interaction (VCI) techniques, both built around the Eckart-Watson Hamiltonian, proved unusually challenging. For the species studied, VPT2 has difficulties yielding dependable results, in some cases even for the fundamentals of the H-containing molecular cations, while carefully executed VCI computations yield considerably improved spectroscopic results. In a few cases unusually large anharmonic corrections to the fundamentals, on the order of 15% of the harmonic value, have been observed.

Potential Catalytic Role of Small Heterocycles in Interstellar H2 Formation: A Laboratory Astrochemistry Study on Furan and Its Hydrogenated Forms

It is now well-accepted in astrochemistry that the formation of interstellar H₂ is taking place on the surface of interstellar grains. It has also been suggested a long time ago that polyaromatic hydrocarbons (PAHs) can catalyze this process by subsequent H atom addition and H abstraction reactions. Recent quantum chemical computations suggested that small heterocycles can be better catalysts than PAHs. In this study, the reaction of H atoms with furan, 2,3- and 2,5-dihydrofurans, and tetrahydrofuran were studied in solid para-H₂ at 3.1 K. The reactions were followed by Fourier transform infrared (FTIR) spectroscopy. By the analysis of spectra, 2-hydrofuran-3-yl, 3-hydrofuran-2-yl, 2,3,4-trihydrofuran-5-yl, and 2,3,5-trihydrofuran-4-yl radicals were identified among the products. The experiments revealed that all the possible H atom addition and H abstraction cycles connecting furan and tetrahydrofuran proceed effectively in both directions at a low temperature. This indicates the possible important role of small heterocycles in interstellar H₂ formation. Furthermore, it also indicates that, in the case of H atom excess, a quasi-equilibrium exists between the c-C₄H_xO (x = 4-8) species, and the ratios of these species in an astrophysical object are determined by the rate of the different H atom addition and H abstraction reaction steps.

Vibrational Dephasing along the Reaction Coordinate of an Electron Transfer Reaction

The role of molecular vibration in photoinduced electron transfer (ET) reactions has been extensively debated in recent years. In this study, we investigated vibrational wavepacket dynamics in a model ET system consisting of an organic dye molecule as an electron acceptor dissolved in various electron donating solvents. By using broad band pump-probe (BBPP) spectroscopy with visible laser pulses of sub-10 fs duration, coherent vibrational wavepackets of naphthacene dye with frequencies spanning 170-1600 cm-1 were observed in the time domain. The coherence properties of 11 vibrational modes were analyzed by an inverse Fourier filtering procedure, and we discovered that the dephasing times of some vibrational coherences are reduced with increasing ET rates. Density functional theory calculations indicated that the corresponding vibrational modes have a large Huang-Rhys factor between the reactant and the product states, supporting the hypothesis that the loss of phase coherence along certain vibrational modes elucidates that those vibrations are coupled to the reaction coordinate of an ET reaction.

Ultrafast "end-on"-to-"side-on" binding-mode isomerization of an iron-carbon dioxide complex

Carbon dioxide (CO₂) binding by transition metals is a captivating phenomenon with a tremendous impact in environmental science and technology, most notably, for establishing circular economies based on greenhouse gas emissions. The molecular and electronic structures of coordination compounds containing CO₂ can be studied in great detail using photochemical precursors bearing the photolabile oxalato-ligand. Here, we study the photoinduced elementary dynamics of the ferric complex, [Fe^III(cyclam)(C₂O₄)]⁺, in dimethyl sulfoxide solution using femtosecond mid-infrared spectroscopy following oxalate-to-iron charge transfer excitation with 266 nm pulses. The pump-probe response in the ν₃-region of carbon dioxide gives unequivocal evidence that a CO₂-molecule is detached from the metal within only 500 fs and with a primary quantum yield of 38%. Simultaneously, a primary ferrous product is formed that carries a carbon dioxide radical anion ligand absorbing at 1649 cm^-1, which is linked to the metal in a bent-O-"end-on" fashion. This primary η_O,bent¹-product is formed with substantial excess vibrational energy, which relaxes on a time scale of several picoseconds. Prior to full thermalization, however, a fraction of the ferrous primary product can structurally isomerize at a rate of 1/(3.5 ps) to a secondary η_CO²-product absorbing at 1727 cm^-1, which features a bent carbon dioxide ligand that is linked to the metal in a "side-on" fashion. The η_O,bent¹-to-η_CO² isomerization requires an intersystem crossing from the sextet to the quartet state, which rationalizes a partial trapping of the system in the metastable bent-O-"end-on" geometry. Finally, a fraction (62%) of the initially photoexcited complexes can return without structural changes to the parent's electronic ground state, but dressed with excess kinetic energy, which relaxes again on a time scale of several picoseconds.

Exploring the anharmonic vibrational structure of carbon dioxide trimers

Our previously developed mbCO2 potential [O. Sode and J. N. Cherry, J. Comput. Chem. 38, 2763 (2017)] is used to describe the vibrational structure of the intermolecular motions of the CO₂ trimers: barrel-shaped and cyclic trimers. Anharmonic corrections are accounted for using the vibrational self-consistent field theory, vibrational second-order Møller-Plesset perturbation (VMP2) theory, and vibrational configuration interaction (VCI) methods and compared with experimental observations. For the cyclic structure, we revise the assignments of two previously observed experimental peaks based on our VCI and VMP2 results. We note that the experimental band observed near 13 cm^-1 is the out-of-phase out-of-plane degenerate motion with E″ symmetry, while the peak observed at 18 cm^-1 likely corresponds to the symmetric out-of-plane torsion A″ vibration. Since the VCI treatment of the vibrational motions accounts for vibrational mixing and delocalization, overtones and combination bands were also observed and quantified in the intermolecular regions of the two trimer isomers.

Cryogenic Vibrationally Resolved Photoelectron Spectroscopy of OH-(H2O): Confirmation of Multidimensional Franck-Condon Simulation Results for the Transition State of the OH + H2O Reaction

We present a transition state spectroscopic study of the OH + H₂O reaction using the experimental technique of cryogenic negative ion photoelectron spectroscopy (NIPES). The recorded NIPE spectrum at 193 nm exhibits multiple vibrational progressions that include excitations to the shared H atom antisymmetric stretching mode with an interval of 0.32 eV as well as other progressions, mainly involving the H bending and O···O symmetric stretching modes. The vertical detachment energy (VDE) was measured at 3.53 eV, whereas an upper limit for the adiabatic detachment energy (ADE) was estimated at 2.90 eV. These values are in excellent agreement with the theoretically computed values of 3.51 and 2.87 eV, respectively, obtained at the CCSD(T)/aug-cc-pV5Z level of theory. The recorded NIPE spectrum is in very good agreement when compared to the one recently reported from four-dimensional Franck-Condon simulations, in which a similar spectral profile was predicted. Besides observing the ground state, we identified a charge-transfer excited state in the form of [OH^-(H₂O)⁺] with a relative energy of 1.39 eV, well matching the previous prediction of 1.36 eV.

Energetics and kinetics of various cyano radical hydrogen abstractions

The cyano radical (CN) is an abundant, open-shell molecule found in a variety of environments, including the atmosphere, the interstellar medium and combustion processes. In these environments, it often reacts with small, closed-shell molecules via hydrogen abstraction. Both carbon and nitrogen atoms of the cyano radical are reactive sites, however the carbon is more reactive with reaction barrier heights generally between 2-15 kcal mol^-1 lower than those of the analogous nitrogen. The CN + HX → HCN/HNC + X, with X = H, CH₃, NH₂, OH, F, SiH₃, PH₂, SH, Cl, C₂H, CN reactions have been studied at a high-level of theory, including CCSD(T)-F12a. Furthermore, kinetics were obtained over the 100-1000 K temperature range, showing excellent agreement with those rate constants that have been determined experimentally.

The structure of ScC2 (X̃2A1): A combined Fourier transform microwave/millimeter-wave spectroscopic and computational study

Pure rotational spectra of Sc¹³C₂ (X̃²A₁) and Sc¹²C¹³C (X̃²A') have been measured using Fourier transform microwave/millimeter-wave methods. These molecules were synthesized in a DC discharge from the reaction of scandium vapor, produced via laser ablation, with ¹³CH₄ or ¹³CH₄/¹²CH₄, diluted in argon. The N_Ka,Kc = 1_0,1 → 0_0,0, 2_0,2 → 1_0,1, 3_0,3 → 2_0,2, and 4_0,4 → 3_0,3 transitions in the frequency range of 14 GHz-61 GHz were observed for both species, each exhibiting hyperfine splittings due to the nuclear spins of ¹³C (I = 1/2) and/or Sc (I = 7/2). These data have been analyzed with an asymmetric top Hamiltonian, and rotational, spin-rotation, and hyperfine parameters have been determined for Sc¹³C₂ and Sc¹²C¹³C. In addition, a quartic force field was calculated for ScC₂ and its isotopologues using a highly accurate coupled cluster-based composite method, incorporating complete basis set extrapolation, scalar relativistic corrections, outer core and inner core electron correlation, and higher-order valence correlation effects. The agreement between experimental and computed rotational constants, including the effective constant (B + C), is ∼0.5% for all three isotopologues. This remarkable agreement suggests promise in predicting rotational spectra of new transition metal-carbon bearing molecules. In combination with previous work on Sc¹²C₂, an accurate structure for ScC₂ has been established using combined experimental (B, C) and theoretical (A) rotational constants. The radical is cyclic (or T-shaped) with r_(Sc-C) = 2.048(2) Å, r_(C-C) = 1.272(2) Å, and ∠_(C-Sc-C) = 36.2(1)°. The experimental and theoretical results also suggest that ScC₂ contains a C₂ ^- moiety and is largely ionic.

Understanding the anharmonic vibrational structure of the carbon dioxide dimer

Understanding the vibrational structure of the CO₂ system is important to confirm the potential energy surface and interactions in such van der Waals complexes. In this work, we use our previously developed mbCO2 potential function to explore the vibrational structure of the CO₂ monomer and dimer. The potential function has been trained to reproduce the potential energies at the CCSD(T)-F12b/aug-cc-pVTZ level of electronic structure theory. The harmonic approximation, as well as anharmonic corrections using vibrational structure theories such as vibrational self-consistent field, vibrational second-order Møller-Plesset perturbation, and vibrational configuration interaction (VCI), is applied to address the vibrational motions. We compare the vibrational results using the mbCO2 potential function with traditional electronic structure theory results and to experimental frequencies. The anharmonic results for the monomer most closely match the experimental data to within 3 cm^-1, including the Fermi dyad frequencies. The intermolecular and intramolecular dimer frequencies were treated separately and show good agreement with the most recent theoretical and experimental results from the literature. The VCI treatment of the dimer vibrational motions accounts for vibrational mixing and delocalization, such that we observe the dimer Fermi resonance phenomena, both in the intramolecular and intermolecular regions.

Vibrational circular dichroism under the quantum magnifying glass: from the electronic flow to the spectroscopic observable

A chemically intuitive method to analyse and interpret vibrational circular dichroism spectra based on the vibrational transition current density.

Fundamental Vibrational Analyses of the HCN Monomer, Dimer and Associated Isotopologues

In this work we provide high level ab initio treatments of the structures, vibrational frequencies, and electronic energies of the HCN monomer and dimer systems along with several isotopologues. The plethora of information related to this system within the literature is summarized and serves as a basis for comparison with the results of this paper. The geometry of the dimer and monomer are reported at the all electroncoupled-cluster singles, doubles, and perturbative triples level of theory [AE-CCSD(T)] with the correlation consistent quadruple-zeta quality basis sets with extra core functions (cc-pCVQZ) from Dunning. The theoretical geometries and electronic structures are further analyzed through the use of the Natural Bond Orbital (NBO) method and Natural Resonance Theory (NRT). At the AE-CCSD(T)/cc-pCVQZ level of theory, the full cubic with semi-diagonal quartic force field for nine dimer and four monomer isotopologues (the parent isotopologue along with 15 N, 13 C, and D derivatives) were obtained to treat the anharmonicity of the vibrations via second order vibrational perturbation theory (VPT2). Lastly, the enthalpy change associated with the formation of the dimer from two monomer units was determined using the focal point analysis. Computations including coupled-cluster through perturbative quadruples as well as basis sets up to six-zeta quality, including core functions (cc-pCVXZ, X=D,T,Q,5,6) were used to extrapolate to the AE-CCSDT(Q)/CBS energy associated with this hydrogen-bond forming process. After appending anharmonic zero-point vibrational, relativistic, and diagonal Born-Oppenheimer corrections, we report a value of -3.93 kcal mol-1 for the enthalpy of formation. To our knowledge, each set of results (geometries, vibrational frequencies, and energetics) reported in this study represents the highest-level and most reliable theoretical predictions reported for this system.

Femtosecond infrared spectroscopy reveals the primary events of the ferrioxalate actinometer

Ultrafast vibrational spectroscopy unravels the primary processes of Hatchard and Parker's ferrioxalate actinometer – a widely used tool in the photochemistry community.

Towards a quantum chemical protocol for the prediction of rovibrational spectroscopic data for transition metal molecules: Exploration of CuCN, CuOH, and CuCCH

High accuracy electronic structure computations for small transition metal-containing molecules have been a long term challenge. Due to coupling between electronic and nuclear wave functions, even experimental/theoretical identification of the ground electronic state requires tremendous efforts. Quartic force fields (QFFs) are effective ab initio tools for obtaining reliable anharmonic spectroscopic properties. However, the method that employs complete basis set limit extrapolation ("C"), consideration of core electron correlation ("cC"), and inclusion of scalar relativity ("R") to produce the energy points on the QFF, the composite CcCR methodology, has not yet been utilized to study inorganic spectroscopy. This work takes the CcCR methodology and adapts it to test whether such an approach is conducive for the closed-shell, copper-containing molecules CuCN, CuOH, and CuCCH. Gas phase rovibrational data are provided for all three species in their ground electronic states. Equilibrium geometries and many higher-order rovibrational properties show good agreement with earlier studies. However, there are notable differences, especially in computation of fundamental vibrational frequencies. Even with further additive corrections for the inner core electron correlation and coupled cluster with full single, double, and triple substitutions (CCSDT), the differences are still larger than expected indicating that more work should follow for predicting rovibrational properties of transition metal molecules.

Effective Inclusion of Mechanical and Electrical Anharmonicity in Excited Electronic States: VPT2-TDDFT Route

We present a reliable and cost-effective procedure for the inclusion of anharmonic effects in excited-state energies and spectroscopic intensities by means of second-order vibrational perturbation theory. This development is made possible thanks to a recent efficient implementation of excited-state analytic Hessians and properties within the time-dependent density functional theory framework. As illustrated in this work, by taking advantage of such algorithmic developments, it is possible to perform calculations of excited-state infrared spectra of medium-large isolated molecular systems, with anharmonicity effects included in both the energy and property surfaces. We also explore the use of this procedure for the inclusion of anharmonic effects in the simulation of vibronic bandshapes of electronic spectra and compare the results with previous, more approximate models.

Characterizing a nonclassical carbene with coupled cluster methods: cyclobutylidene

Carbenes represent a special class of reactive compounds that possess a lone pair of electrons on a carbon atom. Among the myriad examples of carbenes in the literature, cyclobutylidene stands out as a unique nonclassical compound that includes transannular interaction between opposing C1 and C3 carbon atoms within a four-membered ring. On its lowest potential energy surface (X[combining tilde](1)A'), cyclobutylidene quickly rearranges, following three reaction paths: (i) 1,2-H migration; (ii) 1,2-C migration; and, (iii) 1,3-H migration. Herein, this reactivity is examined with high-level coupled-cluster methods [up to CCSDT(Q)]. At this level of theory, combined with extrapolation techniques to obtain energies at the complete basis set (CBS) limit, the long-standing disparity between theoretical and experimental results is resolved. Specifically, cyclobutylidene is predicted to prefer 1,2-C migration rather than 1,2-H migration. Rate constants for the three reaction paths are obtained from canonical variational transition state theory (CVT) and yield reasonable agreement with existing experimental results. Further characterization of cyclobutylidene is also reported: the singlet-triplet gap (ΔES-T) is found to be -9.3 kcal mol(-1) at the CCSDT(Q)/CBS level of theory, and anharmonic vibrational frequencies are determined with second-order vibrational perturbation theory (VPT2).

High-level theoretical rovibrational spectroscopy beyond fc-CCSD(T): The C3 molecule

An accurate local (near-equilibrium) potential energy surface (PES) is reported for the C3 molecule in its electronic ground state (X̃(1)Σg (+)). Special care has been taken in the convergence of the potential relative to high-order correlation effects, core-valence correlation, basis set size, and scalar relativity. Based on the aforementioned PES, several rovibrational states of all (12)C and (13)C substituted isotopologues have been investigated, and spectroscopic parameters based on term energies up to J = 30 have been calculated. Available experimental vibrational term energies are reproduced to better than 1 cm(-1) and rotational constants show relative errors of not more than 0.01%. The equilibrium bond length has been determined in a mixed experimental/theoretical approach to be 1.294 07(10) Å in excellent agreement with the ab initio composite value of 1.293 97 Å. Theoretical band intensities based on a newly developed electric dipole moment function also suggest that the infrared active (1, 1(1), 0)←(0, 0(0), 0) combination band might be observable by high-resolution spectroscopy.

Zero-Cost Estimation of Zero-Point Energies

An additive, linear, atom-type-based (ATB) scheme is developed allowing no-cost estimation of zero-point vibrational energies (ZPVE) of neutral, closed-shell molecules in their ground electronic states. The atom types employed correspond to those defined within the MM2 molecular mechanics force field approach. The reference training set of 156 molecules cover chained and branched alkanes, alkenes, cycloalkanes and cycloalkenes, alkynes, alcohols, aldehydes, carboxylic acids, amines, amides, ethers, esters, ketones, benzene derivatives, heterocycles, nucleobases, all the natural amino acids, some dipeptides and sugars, as well as further simple molecules and ones containing several structural units, including several vitamins. A weighted linear least-squares fit of atom-type-based ZPVE increments results in recommended values for the following atoms, with the number of atom types defined in parentheses: H(8), D(1), B(1), C(6), N(7), O(3), F(1), Si(1), P(2), S(3), and Cl(1). The average accuracy of the ATB ZPVEs is considerably better than 1 kcal mol(-1), that is, better than chemical accuracy. The proposed ATB scheme could be extended to many more atoms and atom types, following a careful validation procedure; deviation from the MM2 atom types seems to be necessary, especially for third-row elements.

CC/DFT Route toward Accurate Structures and Spectroscopic Features for Observed and Elusive Conformers of Flexible Molecules: Pyruvic Acid as a Case Study

The structures and relative stabilities as well as the rotational and vibrational spectra of the three low-energy conformers of pyruvic acid (PA) have been characterized using a state-of-the-art quantum-mechanical approach designed for flexible molecules. By making use of the available experimental rotational constants for several isotopologues of the most stable PA conformer, Tc-PA, the semiexperimental equilibrium structure has been derived. The latter provides a reference for the pure theoretical determination of the equilibrium geometries for all conformers, thus confirming for these structures an accuracy of 0.001 Å and 0.1 deg for bond lengths and angles, respectively. Highly accurate relative energies of all conformers (Tc-, Tt-, and Ct-PA) and of the transition states connecting them are provided along with the thermodynamic properties at low and high temperatures, thus leading to conformational enthalpies accurate to 1 kJ mol(-1). Concerning microwave spectroscopy, rotational constants accurate to about 20 MHz are provided for the Tt- and Ct-PA conformers, together with the computed centrifugal-distortion constants and dipole moments required to simulate their rotational spectra. For Ct-PA, vibrational frequencies in the mid-infrared region accurate to 10 cm(-1) are reported along with theoretical estimates for the transitions in the near-infrared range, and the corresponding infrared spectrum including fundamental transitions, overtones, and combination bands has been simulated. In addition to the new data described above, theoretical results for the Tc- and Tt-PA conformers are compared with all available experimental data to further confirm the accuracy of the hybrid coupled-cluster/density functional theory (CC/DFT) protocol applied in the present study. Finally, we discuss in detail the accuracy of computational models fully based on double-hybrid DFT functionals (mainly at the B2PLYP/aug-cc-pVTZ level) that avoid the use of very expensive CC calculations.