Reconciling spectroscopy with dynamics in global potential energy surfaces: The case of the astrophysically relevant SiC2.

SiC2 is a fascinating molecule due to its unusual bonding and astrophysical importance. In this work, we report the first global potential energy surface (PES) for ground-state SiC2 using the combined-hyperbolic-inverse-power-representation method and accurate ab initio energies. The calibration grid data are obtained via a general dual-level protocol developed afresh herein that entails both coupled-cluster and multi-reference configuration interaction energies jointly extrapolated to the complete basis set limit. Such an approach is specially devised to recover much of the spectroscopy from the PES, while still permitting a proper fragmentation of the system to allow for reaction dynamics studies. Besides describing accurately the valence strongly bound region that includes both the cyclic global minimum and isomerization barriers, the final analytic PES form is shown to properly reproduce dissociation energies, diatomic potentials, and long-range interactions at all asymptotic channels, in addition to naturally reflect the correct permutational symmetry of the potential. Bound vibrational state calculations have been carried out, unveiling an excellent match of the available experimental data on c-SiC2(A11). To further exploit the global nature of the PES, exploratory quasi-classical trajectory calculations for the endothermic C2 + Si → SiC + C reaction are also performed, yielding thermalized rate coefficients for temperatures up to 5000 K. The results hint for the prominence of this reaction in the innermost layers of the circumstellar envelopes around carbon-rich stars, hence conceivably playing therein a key contribution to the gas-phase formation of SiC, and eventually, solid SiC dust.

High-level ab initio quartic force fields and spectroscopic characterization of C2N.

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 C2N-, is largely hindered by a lack of accurate spectroscopic signatures. In this work, we provide such data for both ground -CCN-(3Σ-) and low-lying c-CNC-(1A1) 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-(3Σ-), 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 C2N-. To this end and using the theoretically-predicted spectroscopic constants, the rotational spectra of both -CCN-(3Σ-) and c-CNC-(1A1) 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-(1A1) at about 15.3 kcal mol-1 above the ground-state -CCN-(3Σ-) species.

A global CHIPR potential energy surface for ground-state C3H and exploratory dynamics studies of reaction C2 + CH → C3 + H.

A full-dimensional global potential-energy surface (PES) is first reported for ground-state doublet C3H using the combined-hyperbolic-inverse-power-representation (CHIPR) method and accurate ab initio energies extrapolated to the complete basis set limit. The PES is based on a many-body expansion-type development where the two-body and three-body energy terms are from our previously reported analytic potentials for C2H(2A') and C3(1A',3A'), while the effective four-body one is calibrated using an extension of the CHIPR formalism for tetratomics. The final form is shown to accurately reproduce all known stationary structures of the PES, some of which are unreported thus far, and their interconversion pathways. Moreover, it warrants by built-in construction the appropriate permutational symmetry and describes in a physically reasonable manner all long-range features and the correct asymptotic behavior at dissociation. Exploratory quasi-classical trajectory calculations for the reaction C2 + CH → C3 + H are also performed, yielding thermalized rate coefficients for temperatures up to 4000 K.

Accurate CHIPR Potential Energy Surface for the Lowest Triplet State of C3.

We report the first global ab initio-based potential energy surface (PES) for ground-state triplet C3(3A') based on accurate energies extrapolated to the complete basis set (CBS) limit, and using the combined-hyperbolic-inverse-power-representation method for the analytical modeling. By relying on a cost-effective CBS(D,T) protocol, we ensure that the final form reproduces all topographical features of the PES, including its cyclic-linear isomerization barrier, with CBS(5,6)-quality. To partially account for the incompleteness of the N-electron basis and other minor effects, the available accurate experimental data on the relevant diatomics were used to obtain direct-fit curves that replace the theoretical ones in the many-body expansion. Besides describing properly long-range interactions at all asymptotic channels and permutational symmetry by built-in construction, the PES reported here reproduces the proper exothermicities at dissociation regions as well as the spectroscopy of the diatomic fragments. Bound vibrational state calculations in both linear and cyclic isomers have also been carried out, unveiling a good match of the available data on C3(ã 3Πu), while assisting with IR band positions for C3(3A2') that may serve as a guide for its laboratory and astronomical detection.

Difficulties and Virtues in Assessing the Potential Energy Surfaces of Carbon Clusters via DMBE Theory: Stationary Points of Cκ (κ = 2-10) at the Focal Point.

The previously reported potential energy surfaces (PESs) of ground-state singlet C3 and triplet C4 are here utilized as input for the construction of approximate cluster expansions for larger Cκ (κ = 5-10) species. Relying primarily on the double many-body expansion (DMBE) approach, global potentials are obtained by summing the total interaction energies of all atomic subclusters up to four-body terms. The notable capability of the final forms in predicting good estimates of the linear global minima and their thermochemical/structural properties may provide important insights into the structure-determining nature of the (2 + 3 + 4) terms. The main difficulties and virtues in assessing Cκ's via DMBE theory are analyzed and new prospects given for the construction of global reliable PESs for the target species.

Multiple conical intersections in small linear parameter Jahn-Teller systems: the DMBE potential energy surface of ground-state C3 revisited.

A new single-sheeted DMBE potential energy surface for ground-state C3 is reported. The novel analytical form accurately describes the three symmetry-equivalent C2v disjoint seams, in addition to the symmetry-required D3h one, over the entire configuration space. The present formalism warrants by built-in construction the confluence of the above crossings, and the rotation-in-plane of the C2v seams when the perimeter of the molecule fluctuates. Up to 1050 ab initio energies have been employed in the calibration procedure, of which 421 map the loci of intersection. The calculated energies have been scaled to account for the incompleteness of the basis set and truncation of the MRCI expansion, and fitted analytically with chemical accuracy. The novel form is shown to accurately mimic the region defined by the 4 conical intersections, while exhibiting similar attributes to the previously reported one [J. Chem. Phys., 2015, 143, 074302] at the regions of configuration space away from the crossing seams. Despite being mainly addressed to C3, the present approach should be applicable to adiabatic PESs of any X3 system experiencing similar topological attributes, in particular the small-linear-parameter Jahn-Teller molecules.

Cn (n=2-4): current status.

The major aspects of the C2, C3 and C4 elemental carbon clusters are surveyed. For C2, a brief analysis of its current status is presented. Regarding C3, the most recent results obtained in our group are reviewed with emphasis on modelling its potential energy surface which is particularly complicated due to the presence of multiple conical intersections. As for C4, the most stable isomeric forms of both triplet and singlet spin states and their possible interconversion pathways are examined afresh by means of accurate ab initio calculations. The main strategies for modelling the ground triplet C4 potential are also discussed. Starting from a truncated cluster expansion and a previously reported DMBE form for C3, an approximate four-body term is calibrated from the ab initio energies. The final six-dimensional global DMBE form so obtained reproduces all known topographical aspects while providing an accurate description of the C4 linear-rhombic isomerization pathway. It is therefore commended for both spectroscopic and reaction dynamics studies.This article is part of the theme issue 'Modern theoretical chemistry'.

The Jahn-Teller plus pseudo-Jahn-Teller vibronic problem in the C3 radical and its topological implications.

The combined Jahn-Teller plus pseudo-Jahn-Teller [(E'+A1')⊗e'] problem is discussed for the tricarbon radical (C3) by means of ab initio calculations at the multireference configuration interaction level of theory. For the (1)E' electronic state arising from a e'(2) valence configuration, three additional symmetry-equivalent C2v seams are found to lie in close proximity to the D3h symmetry-required seam over the entire range of the breathing coordinate here considered. As the perimeter of the molecule increases, the C2v disjoint seams approach the D3h one almost linearly and ultimately coalesce with it at Q1 = 5.005 a0, thence forming an intersection node or confluence. By further increasing the size of the molecular triangle, the C2v seams get rotated by ±π in the g-h plane. A three-state vibronic Hamiltonian is also proposed to model locally the title system and shown to accurately mimic the calculated data over the region close to the minimum energy crossing point. No net geometric phase effect is observed when the associated electronic wave functions are adiabatically transported along closed paths encircling the four singularity points. For all paths enclosing the intersection node, the sign reversal criterion is shown to be not fulfilled, even for infinitesimal loops. The results so obtained are expected to be valid for other ring systems experiencing similar topological attributes.

Accurate ab initio-based double many-body expansion potential energy surface for the adiabatic ground-state of the C3 radical including combined Jahn-Teller plus pseudo-Jahn-Teller interactions.

A fully ab initio-based potential energy surface is first reported for the ground electronic state of the C3 radical using the double many-body expansion (DMBE) method. The DMBE form so obtained mimics the full set of energies calculated at the multireference configuration interaction level of theory with chemical accuracy. To account for the incompleteness of the one- and N-electron bases, the calculated external correlation energies have been scaled prior to the fitting procedure via DMBE-scaled external correlation method. Furthermore, the novel potential energy surface reproduces accurately dissociation energies, diatomic potentials, long-range interactions at all asymptotic channels, and the correct topological behavior at the region of 4 conical intersections with the partner state of the same symmetry near equilateral triangular geometries due to combined Jahn-Teller (E' ⊗ e') plus pseudo-Jahn-Teller [(E'+A1')⊗e'] interactions. Rovibrational calculations have also been performed, unveiling a good match of the vibrational spectrum of C3 for 53 calculated levels. The present DMBE form is, therefore, commended for both spectroscopic and reaction dynamics studies, some also performed in the present work.