Index of multi-determinantal and multi-reference character in coupled-cluster theory

Fulminic acid: a quasibent spectacle

Fulminic acid (HCNO) played a critical role in the early development of organic chemistry, and chemists have sought to discern the structure and characteristics of this molecule and its isomers for over 200 years. The mercurial nature of the extremely flat H-C-N bending potential of fulminic acid, with a nearly vanishing harmonic vibrational frequency at linearity, remains enigmatic and refractory to electronic structure theory, as dramatic variation with both orbital basis set and electron correlation method is witnessed. To solve this problem using rigorous electronic wavefunction methods, we have employed focal point analyses (FPA) to ascertain the ab initio limit of optimized linear and bent geometries, corresponding vibrational frequencies, and the HCN + O(3P) → HCNO reaction energy. Electron correlation treatments as extensive as CCSDT(Q), CCSDTQ(P), and even CCSDTQP(H) were employed, and complete basis set (CBS) extrapolations were performed using the cc-pCVXZ (X = 4-6) basis sets. Core electron correlation, scalar relativistic effects (MVD1), and diagonal Born-Oppenheimer corrections (DBOC) were all included and found to contribute significantly in determining whether vibrationless HCNO is linear or bent. At the all-electron (AE) CCSD(T)/CBS level, HCNO is a linear molecule with ω5(H-C-N bend) = 120 cm-1. However, composite AE-CCSDT(Q)/CBS computations give an imaginary frequency (51i cm-1) at the linear optimized geometry, leading to an equilibrium structure with an H-C-N angle of 173.9°. Finally, at the AE-CCSDTQ(P)/CBS level, HCNO is once again linear with ω5 = 45 cm-1, and inclusion of both MVD1 and DBOC effects yields ω5 = 32 cm-1. A host of other topics has also been investigated for fulminic acid, including the dependence of re and ωi predictions on a variety of CBS extrapolation formulas, the question of multireference character, the N-O bond energy and enthalpy of formation, and issues that give rise to the quasibent appellation.

Factorized Quadruples and a Predictor of Higher-Level Correlation in Thermochemistry

Coupled cluster theory has had a momentous impact on the ab initio prediction of molecular properties, and remains a staple ingratiate in high-accuracy thermochemical model chemistries. However, these methods require inclusion of at least some connected quadruple excitations, which generally scale at best as O ( N 9 ) with the number of basis functions. It is very difficult to predict, a priori, the effect correlation of past CCSD(T) on a given reaction energy. The purpose of this work is to examine cost-effective quadruple corrections based on the factorization theorem of the many-body perturbation theory that may address these challenges. We show that the O ( N 7 ) factorized CCSD(TQf) method introduces minimal error to predicted correlation and reaction energies as compared to the O ( N 9 ) CCSD(TQ). Further, we examine the performance of Goodson's continued fraction method in the estimation of CCSDT(Q)Λ contributions to reaction energies as well as a "new" method related to %TAE[(T)] that we refer to as a scaled perturbation estimator. We find that the scaled perturbation estimator based upon CCSD(TQf)/cc-pVDZ is capable of predicting CCSDT(Q)Λ/cc-pVDZ contributions to reaction energies with an average error of 0.07 kcal mol-1 and an L2D of 0.52 kcal mol-1 when applied to a test-suite of nearly 3000 reactions. This offers a means by which to reliably "ballpark" how important post-CCSD(T) contributions are to reaction energies while incurring no more than CCSD(T) formal cost and a little mental math.

Perspective on Coupled-cluster Theory. The evolution toward simplicity in quantum chemistry

Coupled-cluster theory has revolutionized quantum chemistry. It has provided the framework to effectively solve the problem of electron correlation, the main focus of the field for over 60 years. This has enabled ab initio quantum chemistry to provide predictive quality results for most quantities of interest that are obtainable from first-principle calculations. The best that one can do in a basis is the 'full CI,' the exact solution of the non-relativistic Schrödinger equation or, if need be, the relativistic Dirac equation. With due regard to converging the basis set and adequate consideration of higher clusters and relativity in a calculation, virtually predictive results can be obtained. But in addition to its numerical performance, coupled-cluster theory also offers a conceptually new, many-body foundation for the theory that should be appreciated by all practitioners. The latter is emphasized in this perspective, leading to the 'evolution toward simplicity' in the title. The ultimate theory will benefit from the several features that are uniquely exact in coupled-cluster theory and its equation-of-motion (EOM-CC) extensions.

All-Purpose Measure of Electron Correlation for Multireference Diagnostics

We present an analytical relationship between two natural orbital occupancy-based indices, I N D ¯ and INDmax, and two established electron correlation metrics: the leading term of a configuration interaction expansion, c0, and the D2 diagnostic. Numerical validation revealed that I N D ¯ and INDmax can effectively substitute for c0 and D2, respectively. These indices offer three distinct advantages: (i) they are universally applicable across all electronic structure methods, (ii) their interpretation is more intuitive, and (iii) they can be readily incorporated into the development of hybrid electronic structure methods. Additionally, we draw a distinction between correlation measures and correlation diagnostics, establishing MP2 and CCSD numerical thresholds for INDmax, which are to be used as a multireference diagnostic. Our findings further demonstrate that establishing thresholds for other electronic structure methods can be easily accomplished using small data sets.

S-Diagnostic─An a Posteriori Error Assessment for Single-Reference Coupled-Cluster Methods

We propose a novel a posteriori error assessment for the single-reference coupled-cluster (SRCC) method called the S-diagnostic. We provide a derivation of the S-diagnostic that is rooted in the mathematical analysis of different SRCC variants. We numerically scrutinized the S-diagnostic, testing its performance for (1) geometry optimizations, (2) electronic correlation simulations of systems with varying numerical difficulty, and (3) the square-planar copper complexes [CuCl4]2-, [Cu(NH3)4]2+, and [Cu(H2O)4]2+. Throughout the numerical investigations, the S-diagnostic is compared to other SRCC diagnostic procedures, that is, the T1, D1, max T2, and D2 diagnostics as well as different indices of multideterminantal and multireference character in coupled-cluster theory. Our numerical investigations show that the S-diagnostic outperforms the T1, D1, max T2 and D2 diagnostics and is comparable to the indices of multideterminantal and multireference character in coupled-cluster theory in their individual fields of applicability. The experiments investigating the performance of the S-diagnostic for geometry optimizations using SRCC reveal that the S-diagnostic correlates well with different error measures at a high level of statistical relevance. The experiments investigating the performance of the S-diagnostic for electronic correlation simulations show that the S-diagnostic correctly predicts strong multireference regimes. The S-diagnostic, moreover, correctly detects the successful SRCC computations for [CuCl4]2-, [Cu(NH3)4]2+, and [Cu(H2O)4]2+, which have been known to be misdiagnosed by T1 and D1 diagnostics in the past. This shows that the S-diagnostic is a promising candidate for an a posteriori diagnostic for SRCC calculations.

Excited states with pair coupled cluster doubles tailored coupled cluster theory

It is known that some non-dynamic effects of electron correlation can be included in coupled cluster theory using a tailoring technique that separates the effects of non-dynamic and dynamic correlations. Recently, the simple pCCD (pair coupled cluster doubles) wavefunction was shown to provide good results for some non-dynamic correlation problems, such as bond-breaking, in a spin-adapted way with no active space selection. In this paper, we report a study of excited states using "tailored coupled cluster singles and doubles," to attempt to use pCCD as a kernel for more complete coupled-cluster singles and doubles (CCSD) results for excited states. Several excited states are explored from those primarily due to single excitations to those dominated by doubly excited states and from singlet-triplet splittings for some diradical states. For the first two situations, tailored pCCD-TCCSD offers no improvement over equation of motion-CCSD. However, when we explore the singlet-triplet gap of diradical molecules that are manifestly multi-reference, a pCCD kernel provides improved results, particularly with generalized valence bond orbitals.

Are Accelerated and Enhanced Wave Function Methods Accurate to Compute Static Linear and Nonlinear Optical Properties?

Key components of organic-based electro-optic devices are challenging to design or optimize because they exhibit nonlinear optical responses, which are difficult to model or rationalize. Computational chemistry furnishes the tools to investigate extensive collections of molecules in the quest for target compounds. Among the electronic structure methods that provide static nonlinear optical properties (SNLOPs), density functional approximations (DFAs) are often preferred because of their low cost/accuracy ratio. However, the accuracy of the SNLOPs critically depends on the amount of exact exchange and electron correlation included in the DFA, precluding the reliable calculation of many molecular systems. In this scenario, wave function methods such as MP2, CCSD, and CCSD(T) constitute a reliable alternative to compute SNLOPs. Unfortunately, the computational cost of these methods significantly restricts the size of molecules to study, a limitation that hampers the identification of molecules with significant nonlinear optical responses. This paper analyzes various flavors and alternatives to MP2, CCSD, and CCSD(T) methods that either drastically reduce the computational cost or improve their performance but were scarcely and unsystematically employed to compute SNLOPs. In particular, we have tested RI-MP2, RIJK-MP2, RIJCOSX-MP2 (with GridX2 and GridX4 setups), LMP2, SCS-MP2, SOS-MP2, DLPNO-MP2, LNO-CCSD, LNO-CCSD(T), DLPNO-CCSD, DLPNO-CCSD(T0), and DLPNO-CCSD(T1). Our results indicate that all these methods can be safely employed to calculate the dipole moment and the polarizability with average relative errors below 5% with respect to CCSD(T). On the other hand, the calculation of higher-order properties represents a challenge for LNO and DLPNO methods, which present severe numerical instabilities in computing the single-point field-dependent energies. RI-MP2, RIJK-MP2, or RIJCOSX-MP2 are cost-effective methods to compute first and second hyperpolarizabilities with a marginal average error with respect to canonical MP2 (up to 5% for β and up to 11% for γ). More accurate hyperpolarizabilities can be obtained with DLPNO-CCSD(T1); however, this method cannot be employed to obtain reliable second hyperpolarizabilities. These results open the way to obtain accurate nonlinear optical properties at a computational cost that can compete with current DFAs.

Dark Doubly Excited States with Modified Coupled Cluster Models: A Reliable Compromise between Cost and Accuracy?

To treat doubly excited states, the treatment of triple excitations is considered necessary in the framework of equation-of-motion coupled cluster (EOM-CC) methods. We investigate models without explicit triples and seek quantitative measure for the performance of EOM based on CC with singles and doubles (CCSD) or modified CCSD (Distinguishable Cluster Approximation) approaches for states with predominant double excitation character. We also test the efficacy of including triples in perturbative manner through EOM-CCSD(T) and in an iterative way through EOM-CCSDT-3 method. Extended similarity transformed EOM-CCSD(EXT-STEOM-CCSD) method is also tested and provides superior quality results at comparatively low cost. We use the QUEST2 benchmark set of double excitations proposed by Loos et al. [ J. Chem. Theory Comput.2019, 15, 1939] to investigate the performance of methods such as EOM-CCSD, EOM-DCSD, EXT-STEOM-CCSD, ΔCCSD, and ΔDCSD. We also test a tailored CC approach, ΔpairCCD-TCCSD.

Sensitivity of coupled cluster electronic properties on the reference determinant: Can Kohn-Sham orbitals be more beneficial than Hartree-Fock orbitals?

Coupled cluster calculations are traditionally performed over Hartree-Fock reference orbitals (HF-CC methodology). However, in the literature it has been repeatedly raised whether the use of a Kohn-Sham reference (KS-CC methodology) might result in improved performance relative to HF-CC. In the present study, we re-examine the relation of HF-CC and KS-CC methods by comparing the results of widely applied truncated CC calculations (CCSD, CCSD(T), CCSDT) to the limit of full configuration interaction (FCI), which serves as an undebatable reference point of accuracy. Based on a series of CC calculations on diatoms and transition metal complexes, we demonstrate that no systematic improvement of coupled cluster electronic energies, densities and chemical reaction energies is expected when changing from HF to a KS reference. Nevertheless, fortuitous error cancellations might occasionally result in illusory improvement compared to HF-CC. Altogether, the application of KS-CC is not advantageous over HF-CC, but it is also not unreasonable as the choice of reference has negligible influence on the results at sufficiently high CC levels. In addition, KS-CC can be a particularly useful alternative if difficulties are encountered in HF or HF-CC convergence. It is also notable that KS-CC results are found to be practically independent of the chosen density functional, which implies that almost any KS-CC method can be used in place of HF-CC.

Theoretical Rovibrational Spectroscopy of Magnesium Tricarbide-Multireference Character Thwarts a Full Analysis of All Isomers

Magnesium tricarbide isomers are studied herein with coupled cluster theory and multireference configuration interaction to support their possible detection in astrochemical environments such as the circumstellar envelope surrounding the star IRC +10216 or in terrestrial laboratories. Magnesium-bearing species may abound in the interstellar medium (ISM), but only eight (MgNC, MgCN, HMgNC, MgC₂H, MgC₃N, MgC₄H, MgC₅N, and MgC₆H) have been directly identified thus far. Several possible isomers for the related MgC₃ system are explored in their singlet and triplet spin multiplicities. Overall, this work offers quantum chemical insight of rovibrational spectroscopic data for MgC₃ using quartic force fields (QFFs) based on the CCSD(T) and CCSD(T)-F12 levels of theory at the complete basis set (CBS) limit. Additional corrections with small basis set CCSDT(Q) and scalar relativistic effects are also included in the analysis. Salient multireference character is found in the singlet diamond electronic state, which makes a definitive assignment of the ground state challenging. Nevertheless, coupled cluster-based composite energies and multireference configuration interaction both predict that the ¹A₁ diamond isomer is 1.6-2.2 kcal mol^-1 lower in energy than the ³A₁ diamond isomer. Furthermore, highly accurate binding energies of various isomers MgC₃ are provided for comparison to photodetachment experiments. Dipole moments along with harmonic infrared intensities will guide efforts for astronomical and spectroscopic characterization.

Examining fundamental and excitation gaps at the thermodynamic limit: A combined (QTP) DFT and coupled cluster study on trans-polyacetylene and polyacene

Interest in ab initio property prediction of π-conjugated polymers for technological applications places significant demand on “cost-effective” and conceptual computational methods, particularly effective, one-particle theories. This is particularly relevant in the case of Kohn–Sham Density Functional Theory (KS-DFT) and its new competitors that arise from correlated orbital theory, the latter defining the QTP family of DFT functionals. This study presents large, ab initio equation of motion-coupled cluster calculations using the massively parallel ACESIII to target the fundamental bandgap of two prototypical organic polymers, trans-polyacetylene (tPA) and polyacene (Ac), and provides an assessment of the new quantum theory project (QTP) functionals for this problem. Further results focusing on the 1 A g (1 A g), 1 B u (1 B2 u), and 3 B u (3 B2 u) excited states of tPA (Ac) are also presented. By performing calculations on oligomers of increasing size, extrapolations to the thermodynamic limit for the fundamental and all excitation gaps, as well as estimations of the exciton binding energy, are made. Thermodynamic-limit results for a combination of “optimal” and model geometries are presented. Calculated results for excitations that are adequately described using a single-particle model illustrate the benefits of requiring a KS-DFT functional to satisfy the Bartlett ionization potential theorem.

Formation and Evolution of H2C3O+• Radical Cations: A Computational and Matrix Isolation Study

The family of isomeric H₂C₃O^+• radical cations is of great interest for physical organic chemistry and chemistry occurring in extraterrestrial media. In this work, we have experimentally examined a unique synthetic route to the generation of H₂C₃O^+• from the C₂H₂···CO intermolecular complex and also considered the relative stability and monomolecular transformations of the H₂C₃O^+• isomers through high-level ab initio calculations. The structures, energetics, harmonic frequencies, hyperfine coupling constants, and isomerization pathways for several of the most important H₂C₃O^+• isomers were calculated at the UCCSD(T) level of theory. The complementary FTIR and EPR studies in argon matrices at 5 K have demonstrated that the ionized C₂H₂···CO complex transforms into the E-HCCHCO^+• isomer, and this latter species is supposed to be the key intermediate in further chemical transformations, providing a remarkable piece of evidence for kinetic control in low-temperature chemistry. Photolysis of this species at λ = 410-465 nm results in its transformation to the thermodynamically most stable H₂CCCO^+• isomer. Possible implications of the results and potentiality of the proposed synthetic strategy to the preparation of highly reactive organic radical cations are discussed.

DFT Functionals for Modeling of Polyethylene Chains Cross-Linked by Metal Atoms. DLPNO-CCSD(T) Benchmark Calculations

Density functional theory (DFT) functionals for calculations of binding energies (BEs) of the polyethylene (PE) chains cross-linked by selected metal atoms (M) are benchmarked against DLPNO-CCSD(T) and DLPNO-CCSD(T1) data. PEX-M-PEX complexes as compared with plain parallel PEX···PEX chains with X = 3-9 carbon atoms are model species characterized by a cooperative effect of covalent C-M-C bonds and interchain dispersion interactions. The accuracy of DLPNO-CC methods was assessed by a comparison of BEs with the canonical CCSD(T) results for small PE3-M-PE3 complexes. Functionals for PEX···PEX and closed-shell PEX-M-PEX complexes (M = Be, Mg, Zn) were benchmarked against DLPNO-CCSD(T) BEs; open-shell complexes (M = Li, Ag, Au) were benchmarked against the DLPNO-CCSD(T1) method with iterative triples. Three dispersion corrections were combined with 25 DFT functionals for calculations of BEs with respect to PEX-M and PEX fragments employing def2-TZVPP and def2-QZVPP basis sets. Accuracy to within 5% for the closed-shell PEX-M-PEX complexes was achieved with five functionals. Less accurate are functionals for the open-shell PEX-M-PEX complexes; only two functionals deviate by less than 15% from DLPNO-CCSD(T1). Particularly problematic were PEX-Li-PEX complexes. A reasonable overall performance across all complexes in terms of the mean absolute percentage error is found for the range-separated hybrid functionals ωB97X-D3 and CAM-B3LYP/D3(BJ)-ABC.

High-Level ab Initio Predictions for the Ionization Energies, Bond Dissociation Energies, and Heats of Formation of Vanadium Methylidene, Vanadium Methyl Species, and Their Cations (VCH2/VCH2+, VCH3/VCH3+)

The ionization energies of VCH₂ and VCH₃, the various 0 K bond dissociation energies (D₀s) in their neutrals and cations, and their respective heats of formation at 0 and 298 K are computed by the single-reference, wave function-based CCSDTQ/CBS procedure. The core of the composite method is the approximation to the complete basis set (CBS) limit at the coupled cluster (CC) level which includes up to full quadruple excitations. The zero-point vibrational energy, core-valence correlation, spin-orbit coupling, and scalar relativistic effects have their contributions incorporated in an additive manner. For the species in the current study, this protocol requires geometry optimizations and harmonic frequency calculations practically no higher than the CCSD(T)/aug-cc-pwCVTZ and CCSD(T)/aug-cc-pVTZ levels, respectively. The present calculations successfully predict D₀(V⁺-CH₃) = 2.126 eV and D₀(V⁺-CH₂) = 3.298 eV in remarkable agreement with the data recently measured by a spin-orbit state selected V⁺ + CH₄ collision experiment (Phys. Chem. Chem. Phys. 2021, 23, 273-286). The good accord encourages the use of CCSDTQ/CBS protocol in thermochemical predictions of various feasible product channels identified in methane activation by transition metal species.

On the Use of Normalized Metrics for Density Sensitivity Analysis in DFT

In the past years, there has been a discussion about how the errors in density functional theory might be related to errors in the self-consistent densities obtained from different density functional approximations. This, in turn, brings up the discussion about the different ways in which we can measure such errors and develop metrics that assess the sensitivity of calculated energies to changes in the density. It is important to realize that there cannot be a unique metric in order to look at this density sensitivity, simultaneously needing size-extensive and size-intensive metrics. In this study, we report two metrics that are widely applicable to any density functional approximation. We also show how they can be used to classify different chemical systems of interest with respect to their sensitivity to small variations in the density.