An experimental determination of the heat of formation of C2 and the CH bond dissociation energy in C2H

Toward Compact Selected Configuration Interaction Wave Functions with Quantum Monte Carlo─A Case Study of C2

The ¹Σ_g⁺ ground state of C₂ is investigated using truncated CIPSI-Jastrow CSF wave functions with Hartree-Fock orbitals within the framework of variational and diffusion quantum Monte Carlo. The truncation is performed based on the absolute value of the CI coefficients, and the Jastrow, molecular orbitals, and CI parameters are either partially or fully reoptimized with respect to the variational energy. Excellent absolute as well as bond dissociation energies are obtained at DMC level with very compact, fully optimized wave functions. By studying the expansions in more detail, we observe a change in the CI picture when reoptimizing the antisymmetric part of the CIPSI-Jastrow wave functions. Furthermore, we demonstrate that a decrease in the VMC energy as well as an improvement of the nodal surface quality can be achieved─with the same expansion size─if the CSFs are selected in the presence of a Jastrow correlation function, laying the foundation for a Jastrow selected CI scheme with quantum Monte Carlo.

The Spectroscopy of C2: A Cosmic Beacon

ConspectusDicarbon, the molecule formed from two carbon atoms, is among the most abundant molecules in the universe. Said by some to exhibit a quadruple bond, it is bound by more than 6 eV and supports a large number of valence electronic states. It thus has a rich spectroscopy, with 19 one-photon band systems, four of which were discovered by the author and co-workers. Its spectrum was among the first to be described: Wollaston reported the emission spectra from blue flames in 1802.C₂ is observed in a variety of astronomical objects, including stars, circumstellar shells, nebulae, comets and the interstellar medium. It is responsible for the green color of cometary comae but is not observed in the comet tail. It can be observed in absorption and emission by optical spectroscopy in the infrared, visible, and ultraviolet regions of the spectrum, and because it has no electric-dipole-allowed vibrational or rotational transitions, its spectral signature is a sensitive probe of the local environment.Before the work described in this Account, models of C₂ photophysics included the thitherto-unobserved c³Σ_u⁺ state and parametrized the strength of spin-forbidden intercombination transitions. Furthermore, they did not account for photodissociation of C₂, even though it was identified in the 1930s as a key process. Inspired by the observation of C₂ in the Red Rectangle nebula, the author was motivated to instill rigor into C₂ models and embarked on a spectroscopic and computational journey that has lasted 15 years.We were the first to identify the c³Σ_u⁺ state through the d³Π_g-c³Σ_u⁺ transitions, which were to become known as the "Duck" system. This minor partner to the well-known Swan bands is a key part of astrophysical C₂ models and can now be included with rigor. We identified the e³Π_g-c³Σ_u⁺ system, and the c³Σ_u⁺ state is now well-studied. Meanwhile others described the singlet-triplet and triplet-quintet interactions in exquisite detail, allowing rigorous modeling of the a-X and c-X intercombination transitions.The final piece of the C₂ puzzle would be understanding how long it survives before being broken into carbon atom fragments. Though predicted by Herzberg, predissociation in the e³Π_g state had never been observed. To find it would require the complicated ultraviolet spectroscopy of C₂ to be disentangled. In so doing, we identified the 4³Π_g and 3³Π_g states of C₂, thus uncovering two new band systems. The 4³Π_g state allowed the first accurate determination of the ionization energy of C₂. With these new band systems secure, we extracted new levels of the D¹Σ_u⁺ state (Mulliken bands) and the e³Π_g state (Fox-Herzberg bands) from our spectra. Upon climbing the energy ladder in the e³Π_g state to v = 12, we finally identified the route to predissociation of C₂ via non-adiabatic coupling to the d³Π_g state. This observation provided the first laboratory evidence for why C₂ is observed in the coma of a comet but not the tail.

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 C₃H 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 C₂H(²A') and C₃(¹A',³A'), 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 C₂ + CH → C₃ + 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 C₃(³A') 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 C₃(ã ³Π_u), while assisting with IR band positions for C₃(³A₂^') that may serve as a guide for its laboratory and astronomical detection.

Can density cumulant functional theory describe static correlation effects?

We evaluate the performance of density cumulant functional theory (DCT) for capturing static correlation effects. In particular, we examine systems with significant multideterminant character of the electronic wave function, such as the beryllium dimer, diatomic carbon, m-benzyne, 2,6-pyridyne, twisted ethylene, as well as the barrier for double-bond migration in cyclobutadiene. We compute molecular properties of these systems using the ODC-12 and DC-12 variants of DCT and compare these results to multireference configuration interaction and multireference coupled-cluster theories, as well as single-reference coupled-cluster theory with single, double (CCSD), and perturbative triple excitations [CCSD(T)]. For all systems the DCT methods show intermediate performance between that of CCSD and CCSD(T), with significant improvement over the former method. In particular, for the beryllium dimer, m-benzyne, and 2,6-pyridyne, the ODC-12 method along with CCSD(T) correctly predict the global minimum structures, while CCSD predictions fail qualitatively, underestimating the multireference effects. Our results suggest that the DC-12 and ODC-12 methods are capable of describing emerging static correlation effects but should be used cautiously when highly accurate results are required. Conveniently, the appearance of multireference effects in DCT can be diagnosed by analyzing the DCT natural orbital occupations, which are readily available at the end of the energy computation.

Insights into the Perplexing Nature of the Bonding in C2 from Generalized Valence Bond Calculations

Diatomic carbon, C2, has been variously described as having a double, triple, or quadruple bond. In this article, we report full generalized valence bond (GVB) calculations on C2. The GVB wave function-more accurate than the Hartree-Fock wave function and easier to interpret than traditional multiconfiguration wave functions-is well-suited for characterizing the bonding in C2. The GVB calculations show that the electronic wave function of C2 is not well described by a product of singlet-coupled, shared electron pairs (perfect pairing), which is the theoretical basis for covalent chemical bonds. Rather, C2 is best described as having a traditional covalent σ bond with the electrons in the remaining orbitals of the two carbon atoms antiferromagnetically coupled. However, even this description is incomplete as the perfect pairing spin function also makes a significant contribution to the full GVB wave function. The complicated structure of the wave function of C2 is the source of the uncertainty about the nature of the bonding in this molecule.

Understanding the stability, bonding nature and chemical reactivity of 3d-substituted heterofullerenes C58TM (TM = Sc–Zn) from DFT studies

A systematic DFT study for the stability, bonding nature and reactivity of 3d-substituted heterofullerenes C₅₈TM (TM = Sc–Zn).

Probing the limits of accuracy in electronic structure calculations: Is theory capable of results uniformly better than “chemical accuracy”?

Current limitations in electronic structure methods are discussed from the perspective of their potential to contribute to inherent uncertainties in predictions of molecular properties, with an emphasis on atomization energies (or heats of formation). The practical difficulties arising from attempts to achieve high accuracy are illustrated via two case studies: the carbon dimer (C2) and the hydroperoxyl radical (HO2). While the HO2 wave function is dominated by a single configuration, the carbon dimer involves considerable multiconfigurational character. In addition to these two molecules, statistical results will be presented for a much larger sample of molecules drawn from the Computational Results Database. The goal of this analysis will be to determine if a combination of coupled cluster theory with large 1-particle basis sets and careful incorporation of several computationally expensive smaller corrections can yield uniform agreement with experiment to better than "chemical accuracy" (+/-1 kcal/mol). In the case of HO2, the best current theoretical estimate of the zero-point-inclusive, spin-orbit corrected atomization energy (SigmaD0=166.0+/-0.3 kcal/mol) and the most recent Active Thermochemical Table (ATcT) value (165.97+/-0.06 kcal/mol) are in excellent agreement. For C2 the agreement is only slightly poorer, with theory (D0=143.7+/-0.3 kcal/mol) almost encompassing the most recent ATcT value (144.03+/-0.13 kcal/mol). For a larger collection of 68 molecules, a mean absolute deviation of 0.3 kcal/mol was found. The same high level of theory that produces good agreement for atomization energies also appears capable of predicting bond lengths to an accuracy of +/-0.001 A.

CH and C2 Measurements Imply a Radical Pool within a Pool in Acetylene Flames

Measured CH and C2 profiles show a striking resemblance as a function of time in a series of seven well-characterized fuel-rich (phi=1.2-2.0) non-sooting acetylene flames. This implied commonality and interrelationship are unexpected as these radicals have dissimilar chemical kinetic natures. As a result, a rigorous examination was undertaken of the behavior of each of the hydrocarbon species known to be present, C, CH, CH2, CH3, CH4, CHO, CHOH, CH2O, CH2OH, CH3O, CH3OH, C2, C2H, C2H2, CHCO, CH2CO, and C2O. This emphasized the main region where CH and C2 are observed (50-600 micros) and reduced the kinetic reactions to only those that operate efficiently and are dominant. It was immediately apparent that this region of the flame reflects the nature of a hydrogen flame heavily doped with CO and CO2 and containing traces of hydrocarbons. The radical species, H, OH, O, along with H2, H2O, and O2, form an important controlling radical pool that is in partial equilibrium, and the concentrations of each of the hydrocarbon radicals are minor to this, playing secondary roles. As a result, the dominant fast reactions are those between the hydrocarbons and the basic hydrogen/oxygen radicals. Hydrocarbon-hydrocarbon reactions are unimportant here at these equivalence ratios. CH and C2 are formed and destroyed on a sub-microsecond time scale so that their flame profiles are the reflection of a complex kinetically dynamic system. This is found to be the case for all of the hydrocarbon species examined. As might be expected, these rapidly form steady-state distributions. However, with the exceptions of C, CHO, CHOH, and CH2O, which are irreversibly being oxidized, the others all form an interconnected hydrocarbon pool that is under the control of the larger hydrogen radical pool. The hydrocarbon pool can rapidly adjust, and the CH and C2 decay together as the pool is drained. This is either by continuing oxidation in less rich mixtures, or in richer flames where this is negligible by the onset of hydrocarbon-hydrocarbon reactions. The implications of such a hydrocarbon pool are significant. It introduces a buffering effect on their distribution and provides the indirect connection between CH and C2. Moreover, because they are members of this radical pool, flame studies alone cannot answer questions concerning their specific importance in combustion other than their contributing role to this pool. The presence of such a pool modifies the exactness that is needed for kinetic mechanisms, and knowledge of every species in the system no longer is necessary. Furthermore, as rate constants become refined, it will allow for the calculation of the relative concentrations of the hydrocarbon species and facilitate reduced kinetic mechanisms. It provides an explanation for previous isotopically labeled experiments and illustrates the difficulty of exactly identifying in flames the role of individual species. It resolves the fact that differing kinetic models can show similar levels of accuracy and has implications for sensitivity analyses. It finally unveils the mechanism of the flame ionization detector and has implications for the differing interpretations of diamond formation mechanisms.

Sigma- and pi-bond strengths in main group 3-5 compounds

The sigma- and pi-bond strengths for the molecules BH2NH2, BH2PH2, AlH2NH2, and AlH2PH2 have been calculated by using ab initio molecular electronic structure theory at the CCSD(T)/CBS level. The adiabatic pi-bond energy is defined as the rotation barrier between the equilibrium ground-state configuration and the C(s)symmetry transition state for torsion about the A-X bond. We also report intrinsic pi-bond energies corresponding to the adiabatic rotation barrier corrected for the inversion barrier at N or P. The adiabatic sigma-bond energy is defined as the dissociation energy of AH2XH2 to AH2 + XH2 in their ground states minus the adiabatic pi-bond energy. The adiabatic sigma-bond strengths for the molecules BH2NH2, BH2PH2, AlH2NH2, and AlH2PH2 are 109.8, 98.8, 77.6, and 68.3 kcal/mol, respectively, and the corresponding adiabatic pi-bond strengths are 29.9, 10.5, 9.2, and 2.7 kcal/mol, respectively.

Correlation energy extrapolation by intrinsic scaling. IV. Accurate binding energies of the homonuclear diatomic molecules carbon, nitrogen, oxygen, and fluorine

The method of extrapolation by intrinsic scaling, recently introduced to obtain correlation energies, is generalized to multiconfigurational reference functions and used to calculate the binding energies of the diatomic molecules C2, N2, O2, and F2. First, accurate approximations to the full configuration interaction energies of the individual molecules and their constituent atoms are determined, employing Dunning's correlation consistent double-, triple- and quadruple zeta basis sets. Then, these energies are extrapolated to their full basis set limits. Chemical accuracy is attained for the binding energies of all molecules.