Second-order Vibrational Perturbation Theory Research Articles

Use of deperturbed and generalized second-order vibrational perturbation theories to study vibrational spectra intensity and the role of cubic force constant and total anharmonic matrix contributions to anharmonicity: [formula omitted]-thujone and S-camphor

Herein, we employ deperturbed second-order vibrational perturbation theory (DVPT2) and generalized VPT2 (GVPT2) methodologies to predict and analyze two distinct spectral regions in the infrared, VCD, and Raman spectra of α-thujone and S-camphor. The first spectral region spans from 3200 to 2900 cm−1, while the second region extends from 1800 to 500 cm−1. We aim to introduce and validate a program for the simultaneous analysis of vibrational spectra intensity in medium-sized molecules using GVPT2 and DVPT2. Furthermore, employing a graphical presentation (heat map) to illustrate the important contribution of cubic force constant and χ-matrix to anharmonicity.

Infrared Fingerprint and Unimolecular Decay Dynamics of the Hydroperoxyalkyl Intermediate (•QOOH) in Cyclopentane Oxidation.

A transient carbon-centered hydroperoxyalkyl intermediate (•QOOH) in the oxidation of cyclopentane is identified by IR action spectroscopy with time-resolved unimolecular decay to hydroxyl (OH) radical products that are detected by UV laser-induced fluorescence. Two nearly degenerate •QOOH isomers, β- and γ-QOOH, are generated by H atom abstraction of the cyclopentyl hydroperoxide precursor. Fundamental and first overtone OH stretch transitions and combination bands of •QOOH are observed and compared with anharmonic frequencies computed by second-order vibrational perturbation theory. An OH stretch transition is also observed for a conformer arising from torsion about a low-energy CCOO barrier. Definitive identification of the β-QOOH isomer relies on its significantly lower transition state (TS) barrier to OH products, which results in rapid unimolecular decay and near unity branching to OH products. A benchmarking approach is utilized to compute high-accuracy stationary point energies, most importantly TS barriers, for cyclopentane oxidation (C5H9O2), building on higher level reference calculations for ethane oxidation (C2H5O2). The experimental OH product appearance rates are compared with computed statistical microcanonical rates using RRKM theory, including heavy-atom tunneling, thereby validating the computed TS barrier. The results are extended to thermal unimolecular decay rate constants at temperatures and pressures relevant to cyclopentane combustion via master-equation modeling. The various torsional and ring puckering states of the wells and transition states are explicitly considered in these calculations.

Vibrational Features of Oxyamines: A Comparative Study of N,N-Diethylhydroxylamine and N,N-Diethylacetyloxyamine.

The ability of oxyamines to undergo homolytic cleavage of the O-C bond, leading to the formation of stable radicals, is widely used in polymerization processes and to prevent oxidative stress in materials. We present a mid and near-infrared spectroscopy study on two model compounds, N,N-diethylhydroxyloxyamine (C4H11NO) and its acetyl derivative N,N-diethylacetyloxyamine (C6H13NO2) in the liquid phase. The analysis of the spectra is based on a complete exploration of the conformational space, coupled to harmonic and anharmonic calculations performed using the generalized second-order vibrational perturbation theory (GVPT2) formalism at the B3LYP-D3(BJ)/Def2-TZVP level of calculation and potential energy distribution analysis. In the most stable species out of 25, the three amine chains present an all-anti arrangement, with the carbonyl oxygen atom pointing towards the nitrogen lone pair. The simulated spectra are in overall good agreement with the experimental ones, and suitable for the assignment of the main observed bands. Furthermore, similarities and divergences between the two molecules are discussed.

Toward the identification of cyano-astroCOMs via vibrational features: benzonitrile as a test case.

The James Webb Space Telescope (JWST) opened a new era for the identification of molecular systems in the interstellar medium (ISM) by vibrational features. One group of molecules of increasing interest is cyano-derivatives of aromatic organic molecules, which have already been identified in the ISM on the basis of the analysis of rotational signatures, and so, are plausible candidates for the detection by the JWST. Benzonitrile considered in this work represents a suitable example for the validation of a computational strategy, which can be further applied for different, larger, and not-yet observed molecules. For this purpose, anharmonic simulations of infrared (IR) spectra have been compared with recent FTIR experimental studies. The anharmonic computations using the generalized second-order vibrational perturbation theory (GVPT2) in conjunction with a hybrid force field combining the harmonic part of revDSD-PBEP86-D3/jun-cc-pVTZ with anharmonic corrections from B3LYP-D3/SNSD show very good agreement with those in the experiment, with a mean error of for all fundamental transitions overall and only for the stretching fundamental at 4.49 . The inclusion of overtones up to three-quanta transitions also allowed the prediction of spectra in the near-infrared region, which shows distinct features due to overtones at the 2.26 and 1.52 . The remarkable accuracy of the GVPT2 results opens a pathway for the reliable prediction of spectra for a broader range of cyano-astroCOMs.

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.

Infrared signature of the hydroperoxyalkyl intermediate (·QOOH) in cyclohexane oxidation: An isomer-resolved spectroscopic study.

Infrared (IR) action spectroscopy is utilized to characterize carbon-centered hydroperoxy-cyclohexyl radicals (·QOOH) transiently formed in cyclohexane oxidation. The oxidation pathway leads to three nearly degenerate ·QOOH isomers, β-, γ-, and δ-QOOH, which are generated in the laboratory by H-atom abstraction from the corresponding ring sites of the cyclohexyl hydroperoxide (CHHP) precursor. The IR spectral features of jet-cooled and stabilized ·QOOH radicals are observed from 3590 to 7010cm-1 (∼10-20kcal mol-1) at energies in the vicinity of the transition state (TS) barrier leading to OH radicals that are detected by ultraviolet laser-induced fluorescence. The experimental approach affords selective detection of β-QOOH, arising from its significantly lower TS barrier to OH products compared to γ and δ isomers, which results in rapid unimolecular decay and near unity branching to OH products. The observed IR spectrum of β-QOOH includes fundamental and overtone OH stretch transitions, overtone CH stretch transitions, and combination bands involving OH or CH stretch with lower frequency modes. The assignment of β-QOOH spectral features is guided by anharmonic frequencies and intensities computed using second-order vibrational perturbation theory. The overtone OH stretch (2νOH) of β-QOOH is shifted only a few wavenumbers from that observed for the CHHP precursor, yet they are readily distinguished by their prompt vs slow dissociation rates to OH products.

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.

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 6081cm-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.

Probing Ion-Receptor Interactions in Halide Complexes of Octamethyl Calix[4]Pyrrole.

Ion receptors are molecular hosts that bind ionic guests, often with great selectivity. The interplay of solvation and ion binding in anion host-guest complexes in solution governs the binding efficiency and selectivity of such ion receptors. To gain molecular-level insight into the intrinsic binding properties of octamethyl calix[4]pyrrole (omC4P) host molecules with halide guest ions, we performed cryogenic ion vibrational spectroscopy (CIVS) of omC4P in complexes with fluoride, chloride, and bromide ions. We interpret the spectra using density functional theory, describing the infrared spectra of these complexes with both harmonic and anharmonic second-order vibrational perturbation theory (VPT2) calculations. The NH stretching modes of the pyrrole moieties serve as sensitive probes of the ion binding properties, as their frequencies encode the ion-receptor interactions. While scaled harmonic spectra reproduce the experimental NH stretching modes of the chloride and bromide complexes in broad strokes, the high proton affinity of fluoride introduces strong anharmonic effects. As a result, the spectrum of F-·omC4P is not even qualitatively captured by harmonic calculations, but it is recovered very well by VPT2 calculations. In addition, the VPT2 calculations recover the intricate coupling of the NH stretching modes with overtones and combination bands of CH stretching and NH bending modes and with low-frequency vibrations of the omC4P macrocycle, which are apparent for all of the halide ion complexes investigated here. A comparison of the CIVS spectra with infrared spectra of solutions of the same ion-receptor complexes in d3-acetonitrile and d6-acetone shows how ion solvation changes the ion-receptor interactions for the different halide ions.

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 states of HCF, HCCl, HSiF, HNO, and HPO. However, some states have exceptional accuracy in the fundamentals, most notably for ν2 of 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 :CCH2, a molecule of potential astrochemical interest, in order to aid in its potential future experimental rovibronic characterization.

Development of a universal method for vibrational analysis of the terminal alkyne C≡C stretch.

The terminal alkyne C≡C stretch has a large Raman scattering cross section in the "silent" region for biomolecules. This has led to many Raman tag and probe studies using this moiety to study biomolecular systems. A computational investigation of these systems is vital to aid in the interpretation of these results. In this work, we develop a method for computing terminal alkyne vibrational frequencies and isotropic transition polarizabilities that can easily and accurately be applied to any terminal alkyne molecule. We apply the discrete variable representation method to a localized version of the C≡C stretch normal mode. The errors of (1) vibrational localization to the terminal alkyne moiety, (2) anharmonic normal mode isolation, and (3) discretization of the Born-Oppenheimer potential energy surface are quantified and found to be generally small and cancel each other. This results in a method with low error compared to other anharmonic vibrational methods like second-order vibrational perturbation theory and to experiments. Several density functionals are tested using the method, and TPSS-D3, an inexpensive nonempirical density functional with dispersion corrections, is found to perform surprisingly well. Diffuse basis functions are found to be important for the accuracy of computed frequencies. Finally, the computation of vibrational properties like isotropic transition polarizabilities and the universality of the localized normal mode for terminal alkynes are demonstrated.

Exploring Alternate Methods for the Calculation of High-Level Vibrational Corrections of NMR Spin-Spin Coupling Constants.

Traditional nuclear magnetic resonance (NMR) calculations typically treat systems with a Born-Oppenheimer-derived electronic wave function that is solved for a fixed nuclear geometry. One can numerically account for this neglected nuclear motion by averaging over property values for all nuclear geometries with a vibrational wave function and adding this expectation value as a correction to an equilibrium geometry property value. Presented are benchmark coupled-cluster singles and doubles (CCSD) vibrational corrections to spin-spin coupling constants (SSCCs) computed at the level of vibrational second-order perturbation theory (VPT2) using the vibrational averaging driver of the CFOUR program. As CCSD calculations of vibrational corrections are very costly, cheaper electronic structure methods are explored via a newly developed Python vibrational averaging program within the Dalton Project. Namely, results obtained with the second-order polarization propagator approximation (SOPPA) and density functional theory (DFT) with the B3LYP and PBE0 exchange-correlation functionals are compared to the benchmark CCSD//CCSD(T) and experimental values. CCSD//CCSD(T) corrections are also combined with literature CC3 equilibrium geometry values to form the highest-order vibrationally corrected values available (i.e., CC3//CCSD(T) + CCSD//CCSD(T)). CCSD//CCSD(T) statistics showed favorable statistics in comparison to experimental values, albeit at an unfavorably high computational cost. A cheaper CCSD//CCSD(T) + B3LYP method showed quite similar mean absolute deviation (MAD) values as CCSD//CCSD(T), concluding that CCSD//CCSD(T) + B3LYP is optimal in terms of cost and accuracy. With reference to experimental values, a vibrational correction was not worth the cost for all of the other methods tested. Finally, deviation statistics showed that CC3//CCSD(T) + CCSD//CCSD(T) vibrational-corrected equilibrium values deteriorated in comparison to CCSD//CCSD(T) attributed to the use of a smaller basis set or lack of solvation effects for the CC3 equilibrium calculations.

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.

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 abinitio 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.

Fundamental vibrational frequencies of pnictogen (Pn) containing linear triatomic anions: OCPn - and SCPn - where Pn = N, P, As and Sb

The fundamental vibrational frequencies of the isolated OCAs − and SCAs − ions, as well as their heavier antimony analogues, are reported for the first time. A thorough analysis of basis set convergence of the bond lengths and harmonic vibrational frequencies, as well as the examination of anharmonic corrections, is conducted using the MP2 and CCSD(T) ab initio methods with robust basis sets for the OC P n − and SC P n − anions moving down the pnictogen series from Pn = N to Sb. Second-order vibrational perturbation theory (VPT2) and vibrational configuration interaction (VCI) theory are used to compute the fundamental vibrational frequencies of each triatomic anion approaching the MP2 and CCSD(T) complete basis set (CBS) limits. While fundamental vibrational frequencies have been reported for OCN − and SCN − in the gas phase, only experimental vibrational frequencies of the salts for the heavier pnictogen analogues (P and As) have been reported. Experimental Raman vibrational frequencies for OCAs − and SCAs − salts were found to be coincidentally in good agreement with the CCSD(T) VCI frequencies near the CBS limit, differing by at most 13 cm − 1 from the former and 9 cm − 1 from the latter. These modest differences are likely due to the presence of counter-ions and other environmental effects in the solid state, and the overall agreement between the salts and gas-phase frequencies is quite similar to that observed for the OCP − and SCP − ions.

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.

Accuracy Meets Feasibility for the Structures and Rotational Constants of the Molecular Bricks of Life: A Joint Venture of DFT and Wave-Function Methods.

A fully unsupervised computational protocol is proposed with the aim of obtaining reliable structural properties for molecular bricks of life in the gas phase. The results of the new composite scheme approach spectroscopic accuracy at a moderate cost without any empirical parameter in addition to those of the underlying electronic structure method. The whole workflow is fully automated and provides optimized geometries and equilibrium rotational constants. Direct comparison with experimental ground state rotational constants can be performed thanks to the effective computation of vibrational corrections in the framework of second-order vibrational perturbation theory. The results for all the nucleic acid bases and several flexible molecules of biological or medicinal interest show that the accuracy of the new tool is close to that delivered by state-of-the-art composite wave function methods for small semirigid molecules.

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)C3H2 and describing how the components of pbqff can be leveraged in other projects.

Gas-phase formation and spectroscopic characterization of the disubstituted cyclopropenylidenes c-C3(C2H)2, c-C3(CN)2, and c-C3(C2H)(CN)

Aims. The detection of c-C3HC2H and possible future detection of c-C3HCN provide new molecules for reaction chemistry in the dense interstellar medium (ISM) where R-C2H and R-CN species are prevalent. Determination of chemically viable c-C3HC2H and c-C3HCN derivatives and their prominent spectral features can accelerate potential astrophysical detection of this chemical family. This work characterizes three such derivatives: c-C3(C2H)2, c-C3(CN)2, and c-C3(C2H)(CN). Methods. Interstellar reaction pathways of small carbonaceous species are well replicated through quantum chemical means. Highly accurate cc-pVXZ-F12/CCSD(T)-F12 (X = D,T) calculations generate the energetics of chemical formation pathways as well as the basis for quartic force field and second-order vibrational perturbation theory rovibrational analysis of the vibrational frequencies and rotational constants of the molecules under study. Results. The formation of c-C3(C2H)2 is as thermodynamically and, likely, as stepwise favorable as the formation of c-C3HC2H, rendering its detectability to be mostly dependent on the concentrations of the reactants. Both c-C3(C2H)2 and c-C3(C2H)(CN) will be detectable through radioastronomical observation with large dipole moments of 2.84 D and 4.26 D, respectively, while c-C3(CN)2 has an exceedingly small and likely unobservable dipole moment of 0.08 D. The most intense frequency for c-C3(C2H)2 is v2 at 3316.9 cm–1 (3.01 μm), with an intensity of 140 km mol–1. The mixed-substituent molecule c-C3(C2H)(CN) has one frequency with a large intensity, v1, at 3321.0 cm–1 (3.01 μm), with an intensity of 82 km mol–1. The molecule c-C3(CN)2 lacks intense vibrational frequencies within the range that current instrumentation can readily observe. Conclusions. Both c-C3(C2H)2 and c-C3(C2H)(CN) are viable candidates for astrophysical observation, with favorable reaction profiles and spectral data produced herein, but c-C3(CN)2 will not be directly observable through any currently available remote sensing means, even if it forms in large abundances.