CCSDT(Q)/CBS thermochemistry for the D5h → D10h isomerization in the C10 carbon cluster: Getting the right answer for the right reason

First-Principles Path Integral Monte Carlo Studies of the Pseudo Jahn-Teller Effect in the Aromatic Cyclo[10]carbon

There has been renewed interest in carbon nanoscale structures. Experimental measurements at 4.7 K and subsequent first-principles-based vibrational diffusion Monte Carlo simulations at 0 K recently showed that the aromatic cyclo[10]carbon prefers a D5h pentagon-like structure to a regular D10h decagon. This symmetry breaking is due to the second-order Jahn-Teller effect (JTE) and has been amply described in the literature for the cumulenic cyclo[4m + 2]carbon clusters. Yet temperature dependence of the JTE in cyclo[4m + 2]carbon clusters in general and the cyclo[10]carbon in particular has not been studied systematically. In this work, we employ path integral Monte Carlo simulations on a first-principles-derived permutationally invariant potential energy surface (PES) to examine the JTE in cyclo[10]carbon as a function of temperature. The PES was trained on a set of τHCTH/cc-pVQZ energies sampled up to ∼7.7 eV above the D5h global minimum and locally adjusted to a high-level benchmark (reported by others) of the 812 cm-1 electronic energy difference between the D5h global minimum and the D10h transition state. The calculations show a strong JTE at lower temperatures with a dominant D5h composition at 100 K and a gradually diminishing JTE at higher temperatures with a washed-out pentagonal structure above 300 K.

A Low-Order Permutationally Invariant Polynomial Approach to Learning Potential Energy Surfaces Using the Bond-Order Charge-Density Matrix: Application to Cn Clusters for n = 3-10, 20

A representation for learning potential energy surfaces (PESs) in terms of permutationally invariant polynomials (PIPs) using the Hartree-Fock expression for electronic energy is proposed. Our approach is based on the one-electron core Hamiltonian weighted by the configuration-dependent elements of the bond-order charge density matrix (CDM). While the previously reported model used an s-function Gaussian basis for the CDM, the present formulation is expanded with p-functions, which are crucial for describing chemical bonding. Detailed results are demonstrated on linear and cyclic Cn clusters (n = 3-10) trained on extensive B3LYP/aug-cc-pVTZ data. The described method facilitates PES learning by reducing the root mean squared error (RMSE) by a factor of 5 relative to the s-function formulation and by a factor of 20 relative to the conventional PIP approach. This is equivalent to using CDM and an sp basis with a PIP of order M to achieve the same RMSE as with the conventional method with a PIP of order M + 2. Implications for large-scale problems are discussed using the case of the PES of the C20 fullerene in full permutational symmetry.

Quantum Monte Carlo Simulations of the Vibrational Wavefunction of the Aromatic Cyclo[10]carbon Using a Full Dimensional Permutationally Invariant Potential Energy Surface

New experimental measurements [Sun et al., Nature 2023, 623, 972] of the cyclic C10 reveal a cumulenic pentagon-like D5h structure at ∼5 K. However, the long-standing presumption that a large zero-point vibrational energy combined with an extremely flat D5h ↔ D10h ↔ D5h isomerization pathway washes out the pentagonal D5h structure and yields a symmetric D10h decagon remains at odds with the experiment. We resolve this issue with our fitting approach based on a bond-order charge-density matrix expressed in permutationally invariant polynomials. We train the model on τHCTH/cc-pVQZ data morphed to reproduce a relativistic all-electron CCSDT(Q)/CBS D5h-D10h potential energy barrier (benchmarked previously by others). Large scale diffusion Monte Carlo simulations in full dimensionality show that the vibrational ground state of C10 has compositional character of more than 96% D5h, fully reflecting the experimental imaging data. Quantum mechanical variational calculations in 1-D further suggest persistence of the D5h symmetry structure at higher temperatures.

Endocircular Li Carbon Rings

By employing accurate state-of-the-art many-electron quantum-chemistry methods, we establish that monocyclic carbon rings can accommodate Li guest atoms. The low-lying electronic states of these endocircular systems are analyzed and found to include both charge-separated states where the guest Li atom appears as a cation and the ring as an anion and encircled-electron states where Li and the ring are neutral. The electron binding energies of the encircled-electron states increase drastically at their highly symmetric equilibrium geometries with increasing size of the ring, and in Li@C24 , this state becomes the ground state. Li is very weakly bound vertical to the rings in the low-lying encircled-electron states, hinting to van-der-Waals binding. Applcations are mentioned.

Si2C5H2 isomers - search algorithms versus chemical intuition

The pros and cons of using search algorithms alone in identifying new geometries have been discussed by using the Si₂C₅H₂ elemental composition as an example. Within 30 kcal mol^-1 at the CCSD(T)/def2-TZVP//PBE0/def2-TZVP level of theory, the coalescence kick and cuckoo methods postulate merely four isomers (1, 3, 6, and 7) for Si₂C₅H₂ (O. Yañez et. al., Chem. Commun., 2017, 53, 12112). On the contrary, chemical intuition yields fourteen (2, 4, 5, and 8-18) new isomers within the same energy range at the B3LYP/6-311++G(2d,2p) level of theory. Based on the relative energies of the first eleven isomers of Si₂C₅H₂ (1, C_2v, 0.00; 2, C_s, 21.39; 3, C_s, 21.95; 4, C_s, 22.76; 5, C_s, 24.74; 6, C_s, 25.34; 7, C_s, 25.64; 8, C_s, 25.79; 9, C_s, 27.20; 10, C_2v, 28.59; and 11, C_2v, 29.16 kcal mol^-1) calculated at the CCSD(T)/cc-pVTZ level of theory, it is evident that the search algorithms had missed at least seven isomers in the same energy range. The relative energy gaps of isomers 12-18 fall in the range of 30-40 kcal mol^-1 at the latter level of theory. Consequentially, this scenario triggers a speculation going forward with search algorithms alone in the search of all new isomers. While one cannot underestimate the power of these algorithms, the role of chemical intuition may not be completely neglected. Retrospectively, the fourteen new isomers found by chemical intuition may help in writing better search algorithms. All eighteen isomers - including the most stable isomer with a planar tetracoordinate carbon atom 1- remain elusive in the laboratory to date. Thus, structural and spectroscopic parameters have been presented here, which may possibly aid the future experimental studies.