Electron correlation: The many-body problem at the heart of chemistry

Coupled cluster-inspired geminal wavefunctions

Electron pairs have an illustrious history in chemistry, from powerful concepts to understanding structural stability and reactive changes to the promise of serving as building blocks of quantitative descriptions of the electronic structure of complex molecules and materials. However, traditionally, two-electron wavefunctions (geminals) have not enjoyed the popularity and widespread use of the more standard single-particle methods. This has changed recently, with a renewed interest in the development of geminal wavefunctions as an alternative to describing strongly correlated phenomena. Hence, there is a need to find geminal methods that are accurate, computationally tractable, and do not demand significant input from the user (particularly via cumbersome and often ill-behaved orbital optimization steps). Here, we propose new families of geminal wavefunctions inspired by the pair coupled cluster doubles ansatz. We present a new hierarchy of two-electron wavefunctions that extends the one-reference orbital idea to other geminals. Moreover, we show how to incorporate single-like excitations in this framework without leaving the quasiparticle picture. We explore the role of imposing seniority restrictions on these wavefunctions and benchmark these new methods on model strongly correlated systems.

Exploration of the Two-Electron Excitation Space with Data-Driven Coupled Cluster

Computational cost limits the applicability of post-Hartree-Fock methods such as coupled-cluster on larger molecular systems. The data-driven coupled-cluster (DDCC) method applies machine learning to predict the coupled-cluster two-electron amplitudes (t2) using data from second-order perturbation theory (MP2). One major limitation of the DDCC models is the size of training sets that increases exponentially with the system size. Effective sampling of the amplitude space can resolve this issue. Five different amplitude selection techniques that reduce the amount of data used for training were evaluated, an approach that also prevents model overfitting and increases the portability of data-driven coupled-cluster singles and doubles to more complex molecules or larger basis sets. In combination with a localized orbital formalism to predict the CCSD t2 amplitudes, we have achieved a 10-fold error reduction for energy calculations.

Developing correlation-consistent numeric atom-centered orbital basis sets for krypton: Applications in RPA-based correlated calculations

Localized atomic orbitals are the preferred basis set choice for large-scale explicit correlated calculations, and high-quality hierarchical correlation-consistent basis sets are a prerequisite for correlated methods to deliver numerically reliable results. At present, numeric atom-centered orbital (NAO) basis sets with valence correlation consistency (VCC), designated as NAO-VCC-nZ, are only available for light elements from hydrogen (H) to argon (Ar) [Zhang et al., New J. Phys. 15, 123033 (2013)]. In this work, we extend this series by developing NAO-VCC-nZ basis sets for krypton (Kr), a prototypical element in the fourth row of the periodic table. We demonstrate that NAO-VCC-nZ basis sets facilitate the convergence of electronic total-energy calculations using the Random Phase Approximation (RPA), which can be used together with a two-point extrapolation scheme to approach the complete basis set limit. Notably, the Basis Set Superposition Error (BSSE) associated with the newly generated NAO basis sets is minimal, making them suitable for applications where BSSE correction is either cumbersome or impractical to do. After confirming the reliability of NAO basis sets for Kr, we proceed to calculate the Helmholtz free energy for Kr crystal at the theoretical level of RPA plus renormalized single excitation correction. From this, we derive the pressure-volume (P-V) diagram, which shows excellent agreement with the latest experimental data. Our work demonstrates the capability of correlation-consistent NAO basis sets for heavy elements, paving the way toward numerically reliable correlated calculations for bulk materials.

Corrigendum: Coupled cluster theory on modern heterogeneous supercomputers

[This corrects the article DOI: 10.3389/fchem.2023.1154526.].

Coupled cluster theory on modern heterogeneous supercomputers

This study examines the computational challenges in elucidating intricate chemical systems, particularly through ab-initio methodologies. This work highlights the Divide-Expand-Consolidate (DEC) approach for coupled cluster (CC) theory-a linear-scaling, massively parallel framework-as a viable solution. Detailed scrutiny of the DEC framework reveals its extensive applicability for large chemical systems, yet it also acknowledges inherent limitations. To mitigate these constraints, the cluster perturbation theory is presented as an effective remedy. Attention is then directed towards the CPS (D-3) model, explicitly derived from a CC singles parent and a doubles auxiliary excitation space, for computing excitation energies. The reviewed new algorithms for the CPS (D-3) method efficiently capitalize on multiple nodes and graphical processing units, expediting heavy tensor contractions. As a result, CPS (D-3) emerges as a scalable, rapid, and precise solution for computing molecular properties in large molecular systems, marking it an efficient contender to conventional CC models.

Measuring Electron Correlation: The Impact of Symmetry and Orbital Transformations

In this perspective, the various measures of electron correlation used in wave function theory, density functional theory and quantum information theory are briefly reviewed. We then focus on a more traditional metric based on dominant weights in the full configuration solution and discuss its behavior with respect to the choice of the N-electron and the one-electron basis. The impact of symmetry is discussed, and we emphasize that the distinction among determinants, configuration state functions and configurations as reference functions is useful because the latter incorporate spin-coupling into the reference and should thus reduce the complexity of the wave function expansion. The corresponding notions of single determinant, single spin-coupling and single configuration wave functions are discussed and the effect of orbital rotations on the multireference character is reviewed by analyzing a simple model system. In molecular systems, the extent of correlation effects should be limited by finite system size and in most cases the appropriate choices of one-electron and N-electron bases should be able to incorporate these into a low-complexity reference function, often a single configurational one.

Calculation of the ELF in the excited state with single-determinant methods

Since its first definition, back in 1990, the electron localization function (ELF) has settled as one of the most commonly employed techniques to characterize the nature of the chemical bond in real space. Although most of the work using the ELF has focused on the study of ground-state chemical reactivity, a growing interest has blossomed to apply these techniques to the nearly unexplored realm of excited states and photochemistry. Since accurate excited electronic states usually require to account appropriately for electron correlation, the standard single-determinant ELF formulation cannot be blindly applied to them, and it is necessary to turn to correlated ELF descriptions based on the two-particle density matrix (2-PDM). The latter requires costly wavefunction approaches, unaffordable for most of the systems of current photochemical interest. Here, we compare the exact, 2-PDM-based ELF results with those of approximate 2-PDM reconstructions taken from reduced density matrix functional theory. Our approach is put to the test in a wide variety of representative scenarios, such as those provided by the lowest-lying excited electronic states of simple diatomic and polyatomic molecules. Altogether, our results suggest that even approximate 2-PDMs are able to accurately reproduce, on a general basis, the topological and statistical features of the ELF scalar field, paving the way toward the application of cost-effective methodologies, such as time-dependent-Hartree-Fock or time-dependent density functional theory, in the accurate description of the chemical bonding in excited states of photochemical relevance.

Machine learning, artificial intelligence, and chemistry: how smart algorithms are reshaping simulation and the laboratory

Machine learning and artificial intelligence are increasingly gaining in prominence through image analysis, language processing, and automation, to name a few applications. Machine learning is also making profound changes in chemistry. From revisiting decades-old analytical techniques for the purpose of creating better calibration curves, to assisting and accelerating traditional in silico simulations, to automating entire scientific workflows, to being used as an approach to deduce underlying physics of unexplained chemical phenomena, machine learning and artificial intelligence are reshaping chemistry, accelerating scientific discovery, and yielding new insights. This review provides an overview of machine learning and artificial intelligence from a chemist’s perspective and focuses on a number of examples of the use of these approaches in computational chemistry and in the laboratory.

Probing the molecular structure of aqueous triiodide via X-ray photoelectron spectroscopy and correlated electron phenomena

Liquid-microjet-based X-ray photoelectron spectroscopy was applied to aqueous triiodide solutions, I₃^-_(aq.), to investigate the anion's valence- and core-level electronic structure, ionization dynamics, associated electron-correlation effects, and nuclear geometric structure. The roles of multi-active-electron (shake-up) ionization processes - with noted sensitivity to the solute geometric structure - were investigated through I₃^-_(aq.) solution valence, I 4d, and I 3d core-level measurements. The experimental spectra were interpreted with the aid of simulated photoelectron spectra, built upon multi-reference ab initio electronic structure calculations associated with different I₃^-_(aq.) molecular geometries. A comparison of the single-to-multi-active-electron ionization signal ratios extracted from the experimental and theoretical core-level photoemission spectra suggests that the ground state of the solute adopts a near-linear average geometry in aqueous solutions. This contrasts with the interpretation of time-resolved X-ray solution scattering studies, but is found to be fully consistent with the rest of the solution-phase I₃^-_(aq.) literature. Comparing the results of low- and high-photon-energy photoemission measurements, we further suggest that the aqueous anion adopts a more asymmetric geometry at the aqueous-solution-gas interface than in the aqueous bulk.

Benchmarking GNOF against FCI in challenging systems in one, two, and three dimensions

This work assesses the reliability of the recently proposed [M. Piris, Phys. Rev. Lett. 127, 233001 (2021)] global natural orbital functional (GNOF) in the treatment of the strong electron correlation regime. First, we use an H10 benchmark set of four hydrogen model systems of different dimensionalities and distinctive electronic structures: a 1D chain, a 2D ring, a 2D sheet, and a 3D close-packed pyramid. Second, we study two paradigmatic models for strongly correlated Mott insulators, namely, a 1D H50 chain and a 4 × 4 × 4 3D H cube. We show that GNOF, without hybridization to other electronic structure methods and free of tuned parameters, succeeds in treating weak and strong correlation in a more balanced way than the functionals that have preceded it.

Quantum Chemical Approaches to the Calculation of NMR Parameters: From Fundamentals to Recent Advances

Quantum chemical methods for the calculation of indirect NMR spin–spin coupling constants and chemical shifts are always in progress. They never stay the same due to permanently developing computational facilities, which open new perspectives and create new challenges every now and then. This review starts from the fundamentals of the nonrelativistic and relativistic theory of nuclear magnetic resonance parameters, and gradually moves towards the discussion of the most popular common and newly developed methodologies for quantum chemical modeling of NMR spectra.

Electronic Structure and Electron Delocalization in Bare and Dressed Boron Pentamer Clusters

The electronic structures of the lowest energy spin-states of the cationic, neutral and anionic bare boron pentamer clusters have been investigated by means of high level multiconfigurational type calculations, in view of the large static and dynamical electron correlation effects for these species. We found that resembles a singlet spin-state perfect pentagon, which bears no intra-annular chemical bonding interactions, as shown by our analysis of the electron delocalization carried out in terms of the normalized Giambiagi ring-current index, and the total and adjacent atom-pair delocalization indices. However, its lowest-energy triplet and quintet spin-state isomers have C_2v symmetry, with large intra-annular chemical bonding interactions. This geometrical feature extends to both the neutral and the anionic species. Namely, the lowest-energy isomers of boron pentamer neutral and anionic clusters have peripheral and intra-annular sizable bonding interactions reflected in the delocalization of both π- and σ-type valence natural orbitals over the whole molecular plane, which impart large structural stability. In accordance to our calculations, the lowest energy triplet spin-state isomer of the anionic boron pentamer cluster has C₂ symmetry, and consequently, it should show optical activity. Finally, we have studied the change of the geometrical structure of the boron pentamer clusters from planar to compact three-dimensional structures caused by the bonding of ligands to the boron atoms. Our explicit all-electron calculations have been rationalized in terms of the shell-closure of the delocalized valence orbitals of the clusters as predicted by the jellium model extended to nonspherical confinement potentials, circumscribing the role of the ligand to modulate the total number of valence electrons assigned to the core cluster.

Measuring correlated electron motion in atoms with the momentum-balance density

Three new measures of relative electron motion are introduced: equimomentum, antimomentum, and momentum-balance. The equimomentum is the probability that two electrons have the exact same momentum, whereas the antimomentum is the probability that their momenta are the exact opposite. Momentum-balance (MB) is the difference between the equimomentum and antimomentum and, therefore, indicates if equal or opposite momentum is more probable in a system of electrons. The equimomentum, antimomentum, and MB densities are also introduced, which are the local contribution to each quantity. The MB and MB density of the extrapolated-full configuration interaction wave functions of atoms of the first three rows of the periodic table are analyzed, with a particular focus on contrasting the correlated motion of electrons with opposite-spin and parallel-spin. Coulomb correlation between opposite-spin electrons leads to a higher probability of equimomentum, whereas Fermi correlation between parallel-spin electrons leads to a higher probability of antimomentum. The local contribution to MB, given an electron is present, is a minimum at the nucleus and generally increases as the distance from the nucleus increases. There are also interesting similarities between the effects of Fermi correlation and Coulomb correlation (of opposite-spin electrons) on MB.

Towards a Holomorphic Density Functional Theory

Self-consistent-field (SCF) approximations formulated using Hartree-Fock (HF) or Kohn-Sham density-functional theory (KS-DFT) have the potential to yield multiple solutions. However, the formal relationship between multiple solutions identified using HF or KS-DFT remains generally unknown. We investigate the connection between multiple SCF solutions for HF or KS-DFT by introducing a parameterized functional that scales between the two representations. Using the hydrogen molecule and a model of electron transfer, we continuously map multiple solutions from the HF potential to a KS-DFT description. We discover that multiple solutions can coalesce and vanish as the functional changes, forming a direct analogy with the disappearance of real HF solutions along a change in molecular structure. To overcome this disappearance of solutions, we develop a complex-analytic extension of DFT-the "holomorphic DFT" approach-that allows every SCF stationary state to be analytically continued across all molecular structures and exchange-correlation functionals.

Spin-state dependence of exchange-correlation holes

Applications of density-functional theory (DFT) in computational chemistry rely on an approximate exchange-correlation (xc) functional. However, existing approximations can fail dramatically for open-shell molecules, in particular for transition-metal complexes or radicals. Most importantly, predicting energy differences between different spin-states with approximate exchange-correlation functionals remains extremely challenging. Formally, it is known that the exact xc functional should be spin-state dependent, but none of the available approximations feature such an explicit spin-state dependence [C. R. Jacob and M. Reiher, Int. J. Quantum Chem., 2012, 112, 3661-3684]. Thus, to find novel approximations for the xc functional for open-shell systems, the development of spin-state dependent xc functionals appears to be a promising avenue. Here, we set out to shed light on the spin-state dependence of the xc functional by investigating the underlying xc holes, which we extract from configuration interaction calculations for model systems. We analyze the similarities and differences between the xc holes of the lowest-energy singlet and triplet states of the dihydrogen molecule, the helium atom, and the lithium dimer. To shed further light on the spin-state dependence of these xc holes we also discuss exact conditions that can be derived from the spin structure of the reduced two-electron density matrix. Altogether, our results suggest several possible routes towards the construction of explicitly spin-state dependent approximations for the xc functional.

RPA(D) and HRPA(D): Two new models for calculations of NMR indirect nuclear spin-spin coupling constants

In this article, the RPA(D) and HRPA(D) models for the calculation of linear response functions are presented. The performance of the new RPA(D) and HRPA(D) models is compared to the performance of the established RPA, HRPA, and SOPPA models in calculations of indirect nuclear spin-spin coupling constants using the CCSD model as a reference. The doubles correction offers a significant improvement on both the RPA and HRPA models; however, the improvement is more dramatic in the case of the RPA model. For all coupling types investigated in this study, the results obtained using the HRPA(D) model are comparable in accuracy to those given by the SOPPA model, while requiring between 30% and 90% of the calculation time needed for SOPPA. The RPA(D) model, while of slightly lower accuracy compared to the CCSD model than HRPA(D), offered calculation times of only approximately 25% of those required for SOPPA for all the investigated molecules. © 2018 Wiley Periodicals, Inc.

Data-Driven Acceleration of the Coupled-Cluster Singles and Doubles Iterative Solver

Solving the coupled-cluster (CC) equations is a cost-prohibitive process that exhibits poor scaling with system size. These equations are solved by determining the set of amplitudes (t) that minimize the system energy with respect to the coupled-cluster equations at the selected level of truncation. Here, a novel approach to predict the converged coupled-cluster singles and doubles (CCSD) amplitudes, thus the coupled-cluster wave function, is explored by using machine learning and electronic structure properties inherent to the MP2 level. Features are collected from quantum chemical data, such as orbital energies, one-electron Hamiltonian, Coulomb, and exchange terms. The data-driven CCSD (DDCCSD) is not an alchemical method because the actual iterative coupled-cluster equations are solved. However, accurate energetics can also be obtained by bypassing solving the CC equations entirely. Our preliminary data show that it is possible to achieve remarkable speedups in solving the CCSD equations, especially when the correct physics are encoded and used for training of machine learning models.

Orbitals and the Interpretation of Photoelectron Spectroscopy and (e,2e) Ionization Experiments

Electron momentum spectroscopy, scanning tunneling microscopy, and photoelectron spectroscopy provide unique information about electronic structure, but their interpretation has been controversial. This essay discusses a framework for interpretation. Although this interpretation is not new, we believe it is important to present this framework in light of recent publications. The key point is that these experiments provide information about how the electron distribution changes upon ionization, not how electrons behave in the pre-ionized state. Therefore, these experiments do not lead to a "selection of the correct orbitals" in chemistry and do not overturn the well-known conclusion that both delocalized molecular orbitals and localized molecular orbitals are useful for interpreting chemical structure and dynamics. The two types of orbitals can produce identical total molecular electron densities and therefore molecular properties. Different types of orbitals are useful for different purposes.

Computational Approach to Molecular Catalysis by 3d Transition Metals: Challenges and Opportunities

Computational chemistry provides a versatile toolbox for studying mechanistic details of catalytic reactions and holds promise to deliver practical strategies to enable the rational in silico catalyst design. The versatile reactivity and nontrivial electronic structure effects, common for systems based on 3d transition metals, introduce additional complexity that may represent a particular challenge to the standard computational strategies. In this review, we discuss the challenges and capabilities of modern electronic structure methods for studying the reaction mechanisms promoted by 3d transition metal molecular catalysts. Particular focus will be placed on the ways of addressing the multiconfigurational problem in electronic structure calculations and the role of expert bias in the practical utilization of the available methods. The development of density functionals designed to address transition metals is also discussed. Special emphasis is placed on the methods that account for solvation effects and the multicomponent nature of practical catalytic systems. This is followed by an overview of recent computational studies addressing the mechanistic complexity of catalytic processes by molecular catalysts based on 3d metals. Cases that involve noninnocent ligands, multicomponent reaction systems, metal-ligand and metal-metal cooperativity, as well as modeling complex catalytic systems such as metal-organic frameworks are presented. Conventionally, computational studies on catalytic mechanisms are heavily dependent on the chemical intuition and expert input of the researcher. Recent developments in advanced automated methods for reaction path analysis hold promise for eliminating such human-bias from computational catalysis studies. A brief overview of these approaches is presented in the final section of the review. The paper is closed with general concluding remarks.

Towards a formal definition of static and dynamic electronic correlations

Some of the most spectacular failures of density-functional and Hartree-Fock theories are related to an incorrect description of the so-called static electron correlation. Motivated by recent progress in the N-representability problem of the one-body density matrix for pure states, we propose a method to quantify the static contribution to the electronic correlation. By studying several molecular systems we show that our proposal correlates well with our intuition of static and dynamic electron correlation. Our results bring out the paramount importance of the occupancy of the highest occupied natural spin-orbital in such quantification.

Separation of dynamic and nondynamic correlation

The account of electron correlation and its efficient separation into dynamic and nondynamic parts plays a key role in the development of computational methods. In this paper we suggest a physically-sound matrix formulation to split electron correlation into dynamic and nondynamic parts using the two-particle cumulant matrix and a measure of the deviation from idempotency of the first-order density matrix. These matrices are applied to a two-electron model, giving rise to a simplified electron correlation index that (i) depends only on natural orbitals and their occupancies, (ii) can be straightforwardly decomposed into orbital contributions and (iii) splits into dynamic and nondynamic correlation parts that (iv) admit a local version. These expressions are shown to account for dynamic and nondynamic correlation in a variety of systems containing different electron correlation regimes, thus providing the first separation of dynamic and nondynamic correlation using solely natural orbital occupancies.

On the importance of initial conditions for excited-state dynamics

The vast majority of ab initio excited-state simulations are performed within semiclassical, trajectory-based approaches. Apart from the underlying electronic-structure theory, the reliability of the simulations is controlled by a selection of initial conditions for the classical trajectories. We discuss appropriate choices of initial conditions for simulations of different experimental arrangements: dynamics initiated by continuum-wave (CW) laser fields or triggered by ultrashort laser pulses.

The role of topology in organic molecules: origin and comparison of the radical character in linear and cyclic oligoacenes and related oligomers

We discuss the nature of electron-correlation effects in carbon nanorings and nanobelts by a combined approach based on FT-DFT and RAS-SF methods.

Resolving the Notorious Case of Conical Intersections for Coupled Cluster Dynamics

The motion of electrons and nuclei in photochemical events often involves conical intersections, or degeneracies between electronic states. They serve as funnels in nuclear relaxation processes where the electrons and nuclei couple nonadiabatically. Accurate ab initio quantum chemical models are essential for interpreting experimental measurements of such phenomena. In this Letter, we resolve a long-standing problem in coupled cluster theory, presenting the first formulation of the theory that correctly describes conical intersections between excited electronic states of the same symmetry. This new development demonstrates that the highly accurate coupled cluster theory can be applied to describe dynamics on excited electronic states involving conical intersections.

Revisiting the thermochemistry of chlorine fluorides

In this work, accurate calculations of standard enthalpies of formation of chlorine fluorides (ClF_n, n = 1-7; Cl₂ F and Cl₃ F₂ ) were performed through the isodesmic reactions scheme. It is argued that, for many chlorine fluorides, the gold standard method of quantum chemistry (CCSD(T)) is not capable to predict enthalpy values nearing chemical accuracy if atomization scheme is used. This is underpinned by a thorough analysis of total atomization energy results and the inspection of multireference features of these compounds. Other thermodynamic quantities were also calculated at different temperatures. To complement the energetic description, elimination curves were studied through density functional theory as a computationally affordable alternative to highly correlated wave function-based methods. © 2017 Wiley Periodicals, Inc.

From Electrons and Nuclei to Assemblies of Correlated Atoms

What constitutes an atom? The answer may depend on whether you ask a physicist or a chemist. But the complementary views lead to productive overlap in atomic theory. Limiting cases, under high pressures and temperatures, have had a profound impact in a wider context.

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.

Separated-pair approximation and separated-pair pair-density functional theory

Multi-configuration pair-density functional theory (MC-PDFT) has proved to be a powerful way to combine the capabilities of multi-configuration self-consistent-field theory to represent the an electronic wave function with a highly efficient way to include dynamic correlation energy by density functional theory. All applications reported previously involved complete active space self-consistent-field (CASSCF) theory for the reference wave function. For treating large systems efficiently, it is necessary to ask whether good accuracy is retained when using less complete configuration interaction spaces. To answer this question, we present here calculations employing MC-PDFT with the separated pair (SP) approximation, which is a special case (defined in this article) of generalized active space self-consistent-field (GASSCF) theory in which no more than two orbitals are included in any GAS subspace and in which inter-subspace excitations are excluded. This special case of MC-PDFT will be called SP-PDFT. In SP-PDFT, the electronic kinetic energy and the classical Coulomb energy, the electronic density and its gradient, and the on-top pair density and its gradient are obtained from an SP approximation wave function; the electronic energy is then calculated from the first two of these quantities and an on-top density functional of the last four. The accuracy of the SP-PDFT method for predicting the structural properties and bond dissociation energies of twelve diatomic molecules and two triatomic molecules is compared to the SP approximation itself and to CASSCF, MC-PDFT based on CASSCF, CASSCF followed by second order perturbation theory (CASPT2), and Kohn-Sham density functional theory with the PBE exchange-correlation potential. We show that SP-PDFT reproduces the accuracy of MC-PDFT based on the corresponding CASSCF wave function for predicting C-H bond dissociation energies, the reaction barriers of pericyclic reactions and the properties of open-shell singlet systems, all at only a small fraction of the computational cost.

Extrapolation to the Gold-Standard in Quantum Chemistry: Computationally Efficient and Accurate CCSD(T) Energies for Large Molecules Using an Automated Thermochemical Hierarchy

The CCSD(T) method is known as the gold-standard in quantum chemistry and has been the method of choice in quantum chemistry for over 20 years to obtain accurate bond energies and molecular properties. Its computational cost formally scales as the seventh power of the size of the system and can be prohibitive for large molecules. As part of our efforts to reduce the computational cost of the CCSD(T) method yet retain its accuracy, we present a simple, efficient, and user-friendly protocol to extrapolate to CCSD(T) energies in conjunction with MP2 energies. The method is based on the automated error-canceling thermochemical hierarchy previously developed by us called the Connectivity-Based Hierarchy (CBH). For a test set containing 30 diverse nonaromatic organic molecules and biomonomers, we obtain highly accurate extrapolated CCSD(T) energies (with a mean absolute error of only 0.2-0.3 kcal/mol with different basis-set). Additionally, the work also features the successful extrapolation to CCSD energies using a similar protocol.