On the use of quartic force fields in variational calculations

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.

Potential energy surfaces and dynamic properties via ab initio composite and density functional approaches

Vibrational spectroscopy enables critical insight into the structural and dynamic properties of molecules. Presently, the majority of theoretical approaches to spectroscopy employ wavefunction-based ab initio or density functional methods that rely on the harmonic approximation. This approximation breaks down for large molecules with strongly anharmonic bonds or for molecules with large internuclear separations. An alternative to these methods involves generating molecular anharmonic potential energy surfaces (potentials) and using them to extrapolate the vibrational frequencies. This study examines the efficacy of density functional theory (DFT) and the correlation consistent Composite Approach (ccCA) in generating anharmonic frequencies from potentials of small main group molecules. Vibrational self-consistent field Theory (VSCF) and post-VSCF methods were used to calculate the fundamental frequencies of these molecules from their potentials. Functional choice, basis set selection, and mode-coupling are also examined as factors in influencing accuracy. The absolute deviations for the calculated frequencies using potentials at the ccCA level of theory were lower than the potentials at the DFT level. With DFT resulting in bending modes that are better described than those of ccCA, a multilevel DFT:ccCA approach where DFT potentials are used for single vibrational mode potentials and ccCA is used for vibrational mode-mode couplings can be utilized for larger polyatomic systems. The frequencies obtained with this multilevel approach using VCIPSI-PT2 were closer to experimental frequencies than the scaled harmonic frequencies, indicating the success of utilizing post-VSCF methods to generate more accurate representations of computed infrared spectra.

Efficient vibrationally correlated calculations using n-mode expansion-based kinetic energy operators

Due to its efficiency and flexibility, the n-mode expansion is a frequently used tool for representing molecular potential energy surfaces in quantum chemical simulations. In this work, we investigate the performance of n-mode expansion-based models of kinetic energy operators in general polyspherical coordinate systems. In particular, we assess the operators with respect to accuracy in vibrationally correlated calculations and their effect on potential energy surface construction with the adaptive density guided approach. Our results show that the n-mode expansion-based operator variants are reliable and systematically improvable approximations of the full kinetic energy operator. Moreover, we introduce a workflow to generate the n-mode expanded kinetic energy operators on-the-fly within the adaptive density guided approach. This scheme can be applied in studies of species and coordinate systems, for which an analytical form of the kinetic energy operator is not available.

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.

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.

Computing vibrational spectra using a new collocation method with a pruned basis and more points than basis functions: Avoiding quadrature

We present a new collocation method for computing the vibrational spectrum of a polyatomic molecule. Some form of quadrature or collocation is necessary when the potential energy surface does not have a simple form that simplifies the calculation of the potential matrix elements required to do a variational calculation. With quadrature, better accuracy is obtained by using more points than basis functions. To achieve the same advantage with collocation, we introduce a collocation method with more points than basis functions. Critically important, the method can be used with a large basis because it is incorporated into an iterative eigensolver. Previous collocation methods with more points than functions were incompatible with iterative eigensolvers. We test the new ideas by computing energy levels of molecules with as many as six atoms. We use pruned bases but expect the new method to be advantageous whenever one uses a basis for which it is not possible to find an accurate quadrature with about as many points as there are basis functions. For our test molecules, accurate energy levels are obtained even using non-optimal, simple, equally spaced points.

High-level ab initio quartic force fields and spectroscopic characterization of C2N. Phys Chem Chem Phys 2021;23:26227-26240. [PMID: 34787132 DOI: 10.1039/d1cp03505c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

While it is now well established that large carbon chain species and radiative electron attachment (REA) are key ingredients triggering interstellar anion chemistry, the role played by smaller molecular anions, for which REA appears to be an unlikely formation pathway, is as yet elusive. Advancing this research undoubtedly requires the knowledge (and modeling) of their astronomical abundances which, for the case of C₂N^-, is largely hindered by a lack of accurate spectroscopic signatures. In this work, we provide such data for both ground -CCN^-(³Σ^-) and low-lying c-CNC^-(¹A₁) isomers and their singly-substituted isotopologues by means of state-of-the-art rovibrational quantum chemical techniques. Their quartic force fields are herein calibrated using a high-level composite energy scheme that accounts for extrapolations to both one-particle and (approximate) -particle basis set limits, in addition to relativistic effects, with the final forms being subsequently subject to nuclear motion calculations. Besides standard spectroscopic attributes, the full set of computed properties includes fine and hyperfine interaction constants and can be readily introduced as guesses in conventional experimental data reduction analyses through effective Hamiltonians. On the basis of benchmark calculations performed anew for a minimal test set of prototypical triatomics and limited (low-resolution) experimental data for -CCN^-(³Σ^-), the target accuracies are determined to be better than 0.1% of experiment for rotational constants and 0.3% for vibrational fundamentals. Apart from laboratory investigations, the results here presented are expected to also prompt future astronomical surveys on C₂N^-. To this end and using the theoretically-predicted spectroscopic constants, the rotational spectra of both -CCN^-(³Σ^-) and c-CNC^-(¹A₁) are derived and their likely detectability in the interstellar medium is further explored in connection with working frequency ranges of powerful astronomical facilities. Our best theoretical estimate places c-CNC^-(¹A₁) at about 15.3 kcal mol^-1 above the ground-state -CCN^-(³Σ^-) species.

Quantum Chemistry Calculations for Metabolomics

A primary goal of metabolomics studies is to fully characterize the small-molecule composition of complex biological and environmental samples. However, despite advances in analytical technologies over the past two decades, the majority of small molecules in complex samples are not readily identifiable due to the immense structural and chemical diversity present within the metabolome. Current gold-standard identification methods rely on reference libraries built using authentic chemical materials ("standards"), which are not available for most molecules. Computational quantum chemistry methods, which can be used to calculate chemical properties that are then measured by analytical platforms, offer an alternative route for building reference libraries, i.e., in silico libraries for "standards-free" identification. In this review, we cover the major roadblocks currently facing metabolomics and discuss applications where quantum chemistry calculations offer a solution. Several successful examples for nuclear magnetic resonance spectroscopy, ion mobility spectrometry, infrared spectroscopy, and mass spectrometry methods are reviewed. Finally, we consider current best practices, sources of error, and provide an outlook for quantum chemistry calculations in metabolomics studies. We expect this review will inspire researchers in the field of small-molecule identification to accelerate adoption of in silico methods for generation of reference libraries and to add quantum chemistry calculations as another tool at their disposal to characterize complex samples.

A rectangular collocation multi-configuration time-dependent Hartree (MCTDH) approach with time-independent points for calculations on general potential energy surfaces

We introduce a collocation-based multi-configuration time-dependent Hartree (MCTDH) method that uses more collocation points than basis functions. We call it the rectangular collocation MCTDH (RC-MCTDH) method. It does not require that the potential be a sum of products. RC-MCTDH has the important advantage that it makes it simple to use time-independent collocation points. When using time-independent points, it is necessary to evaluate the potential energy function only once and not repeatedly during an MCTDH calculation. It is inexpensive and straightforward to use RC-MCTDH with combined modes. Using more collocation points than basis functions enables one to reduce errors in energy levels without increasing the size of the single-particle function basis. On the contrary, whenever a discrete variable representation is used, the only way to reduce the quadrature error is to increase the basis size, which then also reduces the basis-set error. We demonstrate that with RC-MCTDH and time-independent points, it is possible to calculate accurate eigenenergies of CH₃ and CH₄.

What It Takes to Compute Highly Accurate Rovibrational Line Lists for Use in Astrochemistry

ConspectusWe review the Best Theory + Reliable High-Resolution Experiment (BTRHE) strategy for obtaining highly accurate molecular rovibrational line lists with InfraRed (IR) intensities. The need for highly accurate molecular rovibrational line lists is twofold: (a) assignment of the many rovibrational lines for common stable molecules especially those that exhibit a large amplitude motion, such as NH₃, or have a high density of states such as SO₂; and (b) characterization of the atmospheres of exoplanets, which will be one of the main areas of research in astronomy in the coming decades. The first motivation arises due to the need to eliminate lines due to common molecules in an astronomical observation in order to identify lines from new molecules, while the second motivation arises due to the need to obtain accurate molecular opacities in order to characterize the atmosphere of an exoplanet. The BTRHE strategy first consists of using high-quality ab initio quantum-chemical methods to obtain a global potential energy surface (PES) and dipole moment surface (DMS) that contains the proper physics. The global PES is then refined using a subset of the reliable high-resolution experimental data. The refined PES then gives energy-level predictions to an accuracy similar to the reproduction accuracy of the experimental data used in the refinement step in the interpolation region (i.e., within the range of the experimental data used in the refinement step). The accuracy of the energy levels will slowly degrade as they are extrapolated to spectral regions beyond the high-resolution experimental data used in the refinement step. However, because the degradation is slow, the predicted energy levels can be used to assign new high-resolution experiments, and the data from these can then be used in a subsequent refinement step. In this way, the global PES eventually can yield highly accurate energy levels for all desired spectral regions including to very high energies and high J values. We show that IR intensities computed with the BTRHE rovibrational wave functions and the DMS can be very accurate provided one has minimized the fitting error of the DMS and tested the completeness of the DMS. Some examples of our work on NH₃, CO₂, and SO₂ are given to highlight the usefulness of the BTRHE strategy and to provide ideas on how to further improve its predictive power in the future. In particular, it is shown how successive refinement steps, once new high-resolution data are available, can lead to PESs that yield highly accurate transition energies to larger spectral regions. The importance of including nonadiabatic corrections to reduce the J-dependence of errors for H-containing molecules is shown with work on NH₃. Another very important aspect of the BTRHE approach is the consistency across isotopologues, which allows for highly accurate line lists for any isotopologue once one is obtained for the main isotopologue (which has more high-resolution data available for refinement).

Analysis of the Proton Transfer Bands in the Infrared Spectra of Linear N2H+···OC and N2D+···OC Complexes Using Electric Field-Driven Classical Trajectories

In this work, we describe ab initio calculations and assignment of infrared (IR) spectra of hydrogen-bonded ion-molecular complexes that involve a fluxional proton: the linear N₂H⁺···OC and N₂D⁺···OC complexes. Given the challenges of describing fluxional proton dynamics and especially its IR activity, we use electric field-driven classical trajectories, i.e., the driven molecular dynamics (DMD) method that was developed by us in recent years and for similar applications, in conjunction with high-level electronic structure theory. Namely, we present a modified and a numerically efficient implementation of DMD specifically for direct (or "on the fly") calculations, which we carry out at the MP2-F12/AVDZ level of theory for the potential energy surface (PES) and MP2/AVDZ for the dipole moment surfaces (DMSs). Detailed analysis of the PES, DMS, and the time-dependence of the first derivative of the DMS, referred to as the driving force, for the highly fluxional vibrations involving H⁺/D⁺ revealed that the strongly non-harmonic PES and non-linear DMS yield remarkably complex vibrational spectra. Interestingly, the classical trajectories reveal a doublet in the proton transfer part of the spectrum with the two peaks at 1800 and 1980 cm^-1. We find that their shared intensity is due to a Fermi-like resonance interaction, within the classical limit, of the H⁺ parallel stretch fundamental and an H⁺ perpendicular bending overtone. This doublet is also observed in the deuterated species at 1360 and 1460 cm^-1.

Highly Accurate Quartic Force Field and Rovibrational Spectroscopic Constants for the Azirinyl Cation (c-C2NH2+) and Its Isomers

The azirinyl cation is an aromatic cyclic molecule that is isoelectronic with cyclopropenylidene, c-C₃H₂, and c-C₃H₃⁺. Cyclopropenylidene has been shown to be ubiquitous, existing in many different astrophysical environments. Given the similar chemistry between C and N, and the relative abundances between C and N in astrophysical environments, it is expected that there should be aromatic ringed molecules that incorporate N in the ring, but as yet, no such molecule has been identified. To address this issue, the present study uses high levels of electronic structure theory to compute a highly accurate quartic force field (QFF) for the azirinyl cation and its two lowest lying isomers, the cyanomethyl and isocyanomethyl cations. The theoretical approach uses the singles and doubles coupled-cluster method that includes a perturbative correction for connected triple excitations, CCSD(T), together with extrapolation to the one-particle basis set limit and corrections for scalar relativity and core-correlation. The QFF is then used in a second-order vibrational perturbation theory analysis (VPT2) to compute the fundamental vibrational frequencies and rovibrational spectroscopic constants for all three C₂NH₂⁺ isomers. The reliability of the VPT2 vibrational frequencies is tested by comparison to vibrational configuration interaction (VCI) calculations, and excellent agreement is found between the two approaches. Fundamental vibrational frequencies and rovibrational spectroscopic constants for all singly substituted ¹³C, ¹⁵N, and D isotopologues are also reported. It is expected that the highly accurate spectroscopic data reported herein will be useful in the identification of these cations in high-resolution experimental or astronomical observations.

Spectroscopic investigation of [Al,N,C,O] refractory molecules

As of yet, unexamined aluminum bearing molecules may help elucidate aluminum chemistry and associated refractory atom reactions in the interstellar medium. The flexibility of modern quantum chemistry in the construction and analysis of novel molecules makes it perfectly suited to analyze molecules of astrochemical significance. In this paper, high level ab initio electronic structure calculations using the coupled cluster CCSD(T) and explicitly correlated coupled cluster CCSD(T)-F12 methods with large basis sets extrapolated to the complete basis set limit have been performed on the various [Al,N,C,O] isomers. The anharmonic rotational and vibrational spectroscopic parameters for all isomers are produced with these same levels of theory via quartic force fields and vibrational perturbation theory in order to aid in their potential laboratory or even astrophysical identification. The most stable isomer is determined here to be the aluminum isocyanate radical with linear equilibrium geometry AlNCO (X¹Σ⁺). The NCO antisymmetric stretch of AlNCO has an intensity of 1500 km/mol, which should greatly aid in its infrared detection in the region around 2305 cm^-1. Additionally, the AlOCN isomer is relatively low lying, possesses a 5.12 D dipole moment, and has a notable kinetic stability, making it a viable candidate for astronomical observation. All isomers are characterized by small frequencies, which indicates that these are floppy molecules. Isomers with a terminal aluminum atom are especially floppy, with bending modes less than 100 cm^-1.

Anharmonic vibrational spectroscopy of polycyclic aromatic hydrocarbons (PAHs)

While powerful techniques exist to accurately account for anharmonicity in vibrational molecular spectroscopy, they are computationally very expensive and cannot be routinely employed for large species and/or at non-zero vibrational temperatures. Motivated by the study of Polycyclic Aromatic Hydrocarbon (PAH) emission in space, we developed a new code, which takes into account all modes and can describe all infrared transitions including bands becoming active due to resonances as well as overtone, combination, and difference bands. In this article, we describe the methodology that was implemented and discuss how the main difficulties were overcome, so as to keep the problem tractable. Benchmarking with high-level calculations was performed on a small molecule. We carried out specific convergence tests on two prototypical PAHs, pyrene (C₁₆H₁₀) and coronene (C₂₄H₁₂), aiming at optimising tunable parameters to achieve both acceptable accuracy and computational costs for this class of molecules. We then report the results obtained at 0 K for pyrene and coronene, comparing the calculated spectra with available experimental data. The theoretical band positions were found to be significantly improved compared to harmonic density functional theory calculations. The band intensities are in reasonable agreement with experiments, the main limitation being the accuracy of the underlying calculations of the quartic force field. This is a first step toward calculating moderately high-temperature spectra of PAHs and other similarly rigid molecules using Monte Carlo sampling.

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.

Symmetry breaking and spectral considerations of the surprisingly floppy c-C3H radical and the related dipole-bound excited state of c-C3H. J Chem Phys 2017;146:224303. [PMID: 29166048 DOI: 10.1063/1.4985095]

The C₃H radical is believed to be prevalent throughout the interstellar medium and may be involved in the formation of polycyclic aromatic hydrocarbons. C₃H exists as both a linear and a cyclic isomer. The C_2v cyclopropenylidenyl radical isomer was detected in the dark molecular cloud TMC-1, and the linear propenylidenyl radical isomer has been observed in various dark molecular clouds. Even though the c-C₃H radical has been classified rotationally, the vibrational frequencies of this seemingly important interstellar molecule have never been directly observed. Established, highly accurate quartic force field methodologies are employed here to compute useful geometrical data, spectroscopic constants, and vibrational frequencies. The computed rotational constants are consistent with the experimental results. Consequently, the three a₁ (ν₁, ν₂, and ν₃) and one b₁ (ν₆) anharmonic vibrational frequencies at 3117.7 cm^-1, 1564.3 cm^-1, 1198.5 cm^-1, and 826.7 cm^-1, respectively, are reliable predictions for these, as of yet unseen, observables. Unfortunately, the two b₂ fundamentals (ν₄ and ν₅) cannot be treated adequately in the current approach due to a flat and possible double-well potential described in detail herein. The dipole-bound excited state of the anion suffers from the same issues and may not even be bound. However, the trusted fundamental vibrational frequencies described for the neutral radical should not be affected by this deformity and are the first robustly produced for c-C₃H. The insights gained here will also be applicable to other structures containing three-membered bare and exposed carbon rings that are surprisingly floppy in nature.

Vibrational frequencies and spectroscopic constants of three, stable noble gas molecules: NeCCH+, ArCCH+, and ArCN+

The search for possible, natural, noble gas molecules has led to quantum chemical, spectroscopic analysis of NeCCH⁺, ArCCH⁺, and ArCN⁺.

The performance of low-cost commercial cloud computing as an alternative in computational chemistry

The growth of commercial cloud computing (CCC) as a viable means of computational infrastructure is largely unexplored for the purposes of quantum chemistry. In this work, the PSI4 suite of computational chemistry programs is installed on five different types of Amazon World Services CCC platforms. The performance for a set of electronically excited state single-point energies is compared between these CCC platforms and typical, "in-house" physical machines. Further considerations are made for the number of cores or virtual CPUs (vCPUs, for the CCC platforms), but no considerations are made for full parallelization of the program (even though parallelization of the BLAS library is implemented), complete high-performance computing cluster utilization, or steal time. Even with this most pessimistic view of the computations, CCC resources are shown to be more cost effective for significant numbers of typical quantum chemistry computations. Large numbers of large computations are still best utilized by more traditional means, but smaller-scale research may be more effectively undertaken through CCC services.

Quartic force fields for excited electronic states: rovibronic reference data for the 1 (2)A' and 1 (2)A″ states of the isoformyl radical, HOC

Quartic force fields (QFFs) have been shown to be an effective, accurate, and relatively compact means of computing rovibrational spectroscopic data for numerous molecules with numerous applications. However, excited states have been nearly excluded from the this approach since most accurate QFFs are based on the "gold standard" coupled cluster singles, doubles, and perturbative triples [CCSD(T)] method which is not readily extended to excited states. In this work, rovibronic spectroscopic data is provided for the isoformyl radical, a molecule of significance in combustion and astrochemistry, both through the traditional means of variational access to excited states with CCSD(T) and in the novel extension of QFFs routinely to treat electronically excited states through the standard coupled cluster excited state approach, equation of motion (EOM) CCSD. It is shown here that the new EOM-based QFF provides structural parameters and rotational constants that are quite close to those from a related CCSD(T)-based QFF for the 1 (2)A(″) excited state of HOC. The anharmonic vibrational frequency percent differences between the two QFFs are less than 0.4% for the O-H stretch, less than 1.9% for the C-O stretch, and around 3.0% for the bend. Even so, the pure excited state EOM-QFF anharmonic frequencies are still very good abinitio representations that may be applied to systems where electronically excited states are not variationally accessible. Additionally, rovibrational spectroscopic data is provided for the 1 (2)A(') ground state of HOC and for both the ground and excited state of DOC.

Anharmonic rovibrational calculations of singlet cyclic C4 using a new ab initio potential and a quartic force field

We report a CCSD(T)/cc-pCV5Z quartic force field (QFF) and a semi-global CCSD(T)-F12b/aug-cc-pVTZ potential energy surface (PES) for singlet, cyclic C4. Vibrational fundamentals, combinations, and overtones are obtained using vibrational second-order perturbation theory (VPT2) and the vibrational configuration-interaction (VCI) approach. Agreement is within 10 cm(-1) between the VCI calculated fundamentals on the QFF and PES using the MULTIMODE (MM) program, and VPT2 and VCI results agree for the fundamentals. The agreement between VPT2-QFF and MM-QFF results is also good for the C4 combinations and overtones. The J = 1 and J = 2 rovibrational energies are reported from both VCI (MM) on the PES and VPT2 on the QFF calculations. The spectroscopic constants of (12)C4 and two C2v-symmetry, single (13)C-substituted isotopologues are presented, which may help identification of cyclic C4 in future experimental analyses or astronomical observations.

Vibrational frequencies and spectroscopic constants for 1 3A' HNC and 1 3A' HOC+ from high-accuracy quartic force fields

The spectroscopic constants and vibrational frequencies for the 1 (3)A' states of HNC, DNC, HOC(+), and DOC(+) are computed and discussed in this work. The reliable CcCR quartic force field based on high-level coupled cluster ab initio quantum chemical computations is exclusively utilized to provide the anharmonic potential. Then, second-order vibrational perturbation theory and vibrational configuration interaction methods are employed to treat the nuclear Schrödinger equation. Second-order perturbation theory is also employed to provide spectroscopic data for all molecules examined. The relationship between these molecules and the corresponding 1 (3)A' HCN and HCO(+) isomers is further developed here. These data are applicable to laboratory studies involving formation of HNC and HOC(+) as well as astronomical observations of chemically active astrophysical environments.