Constrained-pairing mean-field theory. II. Exact treatment of dissociations to nondegenerate orbitals

Particle-Breaking Unrestricted Hartree-Fock Theory for Open Molecular Systems

We recently introduced the particle-breaking restricted Hartree-Fock (PBRHF) model, a mean-field approach to address the fractional charging of molecules when they interact with an electronic environment. In this paper, we present an extension of the model referred to as particle-breaking unrestricted Hartree-Fock (PBUHF). The unrestricted formulation contains odd-electron states necessary for a realistic description of fractional charging. Within the PBUHF parametrization, we use two-body operators as they yield convenient operator transformations. However, two-body operators can change only the particle number by two. Therefore, we include noninteracting zero-energy bath orbitals to generate a linear combination of even and odd electron states. Depending on whether the occupied or virtual orbitals of a molecule interact with the environment, the average number of electrons is either decreased or increased. Without interaction, PBUHF reduces to the unrestricted Hartree-Fock wave function.

Particle-Breaking Hartree-Fock Theory for Open Molecular Systems

In this work we present the particle-breaking Hartree-Fock (PBHF) model which is a mean-field approach to open molecular systems. The interaction of a system with the environment is parametrized through a particle-breaking term in the molecular Hamiltonian. The PBHF wave function is constructed through an exponential unitary transformation of a Slater determinant with a given number of electrons. We consider only the closed-shell formalism. The parametrization results in a linear combination of Slater determinants with different numbers of electrons, i.e., the PBHF wave function is not an eigenfunction of the number operator. As a result, the density matrix may have fractional occupations which are, because of the unitary parametrization, always between 0.0 and 2.0. The occupations are optimized simultaneously with the orbitals, using the trust-region optimization procedure. In the limit of a particle-conserving Hamiltonian, the PBHF optimization will converge to a standard Hartree-Fock wave function. We show that the average number of electrons may be decreased or increased depending on whether the particle-breaking term affects occupied or virtual orbitals.

Coupled Cluster Benchmark of New DFT and Local Correlation Methods: Mechanisms of Hydroarylation and Oxidative Coupling Catalyzed by Ru(II, III) Chloride Carbonyls

We have evaluated a set of accurate canonical CCSD(T) energies for stationary points on the potential energy surface for Ru(II, III) chloride carbonyl catalysis of two competing reactions between benzene and methyl acrylate (MA), namely, hydroarylation and oxidative coupling. We have then applied this set to evaluate the performance of localized orbital coupled-cluster methods and several new and common density functionals. We find that (a) DLPNO-CCSD(T) with TightPNO cutoffs is an acceptable substitute for full canonical CCSD(T) calculations on this system; (b) for the closed-shell systems where it could be applied, LNO-CCSD(T) with tight convergence criteria is very close to the canonical results; (c) the recent ωB97X-V and ωB97M-V functionals exhibit superior performance to commonly used DFT functionals in both closed- and open-shell calculations; (d) the revDSD-PBEP86 revision of the DSD-PBEP86 double hybrid represents an improvement over the original, even though transition metals were not involved in its parametrization; and (e) DSD-SCAN and DOD-SCAN show comparable efficiency. Most tested (meta)-GGA and hybrid density functionals perform better for open-shell than for closed-shell complexes; this is not the case for the double hybrids considered.

Data-Driven Approaches Can Overcome the Cost-Accuracy Trade-Off in Multireference Diagnostics

High-throughput computational screening typically employs methods (i.e., density functional theory or DFT) that can fail to describe challenging molecules, such as those with strongly correlated electronic structure. In such cases, multireference (MR) correlated wavefunction theory (WFT) would be the appropriate choice but remains more challenging to carry out and automate than single-reference (SR) WFT or DFT. Numerous diagnostics have been proposed for identifying when MR character is likely to have an effect on the predictive power of SR calculations, but conflicting conclusions about diagnostic performance have been reached on small data sets. We compute 15 MR diagnostics, ranging from affordable DFT-based to more costly MR-WFT-based diagnostics, on a set of 3165 equilibrium and distorted small organic molecules containing up to six heavy atoms. Conflicting MR character assignments and low pairwise linear correlations among diagnostics are also observed over this set. We evaluate the ability of existing diagnostics to predict the percent recovery of the correlation energy, %E_corr. None of the DFT-based diagnostics are nearly as predictive of %E_corr as the best WFT-based diagnostics. To overcome the limitation of this cost-accuracy trade-off, we develop machine learning (ML, i.e., kernel ridge regression) models to predict WFT-based diagnostics from a combination of DFT-based diagnostics and a new, size-independent 3D geometric representation. The ML-predicted diagnostics correlate as well with MR effects as their computed (i.e., with WFT) values, significantly improving over the DFT-based diagnostics on which the models were trained. These ML models thus provide a promising approach to improve upon DFT-based diagnostic accuracy while remaining suitably low cost for high-throughput screening.

Performance of Localized Coupled Cluster Methods in a Moderately Strong Correlation Regime: Hückel-Möbius Interconversions in Expanded Porphyrins

Localized orbital coupled cluster theory has recently emerged as a nonempirical alternative to DFT for large systems. Intuitively, one might expect such methods to perform less well for highly delocalized systems. In the present work, we apply both canonical CCSD(T) approximations and a variety of localized approximations to a set of flexible expanded porphyrins—macrocycles that can switch between Hückel, figure-eight, and Möbius topologies under external stimuli. Both minima and isomerization transition states are considered. We find that Möbius(-like) structures have much stronger static correlation character than the remaining structures, and that this causes significant errors in DLPNO-CCSD(T) and even DLPNO-CCSD(T₁) approaches, unless TightPNO cutoffs are employed. If sub-kcal mol^–1 accuracy with respect to canonical relative energies is required even for Möbius-type systems (or other systems plagued by strong static correlation), then Nagy and Kallay’s LNO-CCSD(T) method with “tight” settings is the suitable localized approach. We propose the present POLYPYR21 data set as a benchmark for localized orbital methods or, more broadly, for the ability of lower-level methods to handle energetics with strongly varying degrees of static correlation.

Extended separated-pair approximation for transition metal potential energy curves

Developing a computational method that is both affordable and accurate for transition-metal chemistry is a major challenge. The bond dissociation energies and the potential energy curves are two important targets for theoretical prediction. Here, we investigate the performance of multiconfiguration pair-density functional theory (MC-PDFT) based on wave functions calculated by the complete-active-space (CAS) and generalized active space (GAS) self-consistent-field (SCF) methods for three transition-metal diatomics (TiC, TiSi, and WCl) for which accurate bond energies are available from recent experiments. We compare the results to those obtained by CAS second-order perturbation theory (CASPT2) and Kohn-Sham DFT (KS-DFT). We use six systematic methods to choose the active spaces: (1) we put the bonding orbitals, antibonding orbitals, and singly occupied nonbonding orbitals into the active space in the first method; (2) we also put s and p valence orbitals into the active space; we tried two levels of correlated participating orbitals (CPO) active spaces: (3) nominal CPO (nom-CPO) and (4) extended CPO (ext-CPO); and we used (5) the separated-pair (SP) approximation and (6) a new method presented here called extended separate pairs (ESP) approximation to divide the nom-CPO active space into subspaces. Schemes 1-4 are carried out within the CAS framework, and schemes 5 and 6 are carried out in the GAS framework to eliminate deadwood configurations. For TiC and TiSi, we used all six kinds of active spaces. For WCl, we used three active spaces (nom-CPO, SP, and ESP). We found that MC-PDFT performs better than both CASPT2 and KS-DFT. We also found that the SP (for TiSi) and ESP (for TiC and WCl) approximations are particularly appealing because they make the potential curves smoother and significantly decrease the computational cost of CASSCF calculations. Furthermore, ESP-PDFT can be as accurate as CAS-PDFT.

Geminal perturbation theory based on the unrestricted Hartree-Fock wavefunction

A perturbative correction exploiting natural orbitals and the pair function structure of the unrestricted Hartree-Fock (UHF) wavefunction is devised. The method offers a simple framework for describing multireference systems where static correlation is captured by UHF. The UHF wavefunction is built of two-electron fragments (geminals), involving both singlet and triplet (m_s = 0) parts. At order zero of the perturbative treatment, configuration interaction coefficients of UHF geminals are relaxed. The zero order Hamiltonian is of the Dyall-type, including explicit two-electron interaction within geminals and leading to a formal 6th power scaling. Adopting an effective one-electron zero order Hamiltonian term for the subset of virtual orbitals reduces scaling of the correction step to 4th power. Formal properties of the proposed schemes are discussed. Energetic data and natural occupation numbers of illustrative test systems are used to assess the new approach. The cases where the wavefunction becomes essentially spin pure at the level of reference show good performance. Spin contamination remaining at order zero is found to undermine the perturbative correction.

Do CCSD and approximate CCSD-F12 variants converge to the same basis set limits? The case of atomization energies

While the title question is a clear "yes" from purely theoretical arguments, the case is less clear for practical calculations with finite (one-particle) basis sets. To shed further light on this issue, the convergence to the basis set limit of CCSD (coupled cluster theory with all single and double excitations) and of different approximate implementations of CCSD-F12 (explicitly correlated CCSD) has been investigated in detail for the W4-17 thermochemical benchmark. Near the CBS ([1-particle] complete basis set) limit, CCSD and CCSD(F12^*) agree to within their respective uncertainties (about ±0.04 kcal/mol) due to residual basis set incompleteness error, but a nontrivial difference remains between CCSD-F12b and CCSD(F12^*), which is roughly proportional to the degree of static correlation. The observed basis set convergence behavior results from the superposition of a rapidly converging, attractive, CCSD[F12]-CCSD-F12b difference (consisting mostly of third-order terms) and a more slowly converging, repulsive, fourth-order difference between CCSD(F12^*) and CCSD[F12]. For accurate thermochemistry, we recommend CCSD(F12^*) over CCSD-F12b if at all possible. There are some indications that the nZaPa family of basis sets exhibits somewhat smoother convergence than the correlation consistent family.

Multiconfiguration Pair-Density Functional Theory: A New Way To Treat Strongly Correlated Systems

The electronic energy of a system provides the Born-Oppenheimer potential energy for internuclear motion and thus determines molecular structure and spectra, bond energies, conformational energies, reaction barrier heights, and vibrational frequencies. The development of more efficient and more accurate ways to calculate the electronic energy of systems with inherently multiconfigurational electronic structure is essential for many applications, including transition metal and actinide chemistry, systems with partially broken bonds, many transition states, and most electronically excited states. Inherently multiconfigurational systems are called strongly correlated systems or multireference systems, where the latter name refers to the need for using more than one ("multiple") configuration state function to provide a good zero-order reference wave function. This Account describes multiconfiguration pair-density functional theory (MC-PDFT), which was developed as a way to combine the advantages of wave function theory (WFT) and density functional theory (DFT) to provide a better treatment of strongly correlated systems. First we review background material: the widely used Kohn-Sham DFT (which uses only a single Slater determinant as reference wave function), multiconfiguration WFT methods that treat inherently multiconfigurational systems based on an active space, and previous attempts to combine multiconfiguration WFT with DFT. Then we review the formulation of MC-PDFT. It is a generalization of Kohn-Sham DFT in that the electron kinetic energy and classical electrostatic energy are calculated from a reference wave function, while the rest of the energy is obtained from a density functional. However, there are two main differences with respent to Kohn-Sham DFT: (i) The reference wave function is multiconfigurational rather than being a single Slater determinant. (ii) The density functional is a function of the total density and the on-top pair density rather than being a function of the spin-up and spin-down densities. In work carried out so far, the multiconfigurational wave function is a multiconfiguration self-consistent-field wave function. The new formulation has the advantage that the reference wave function has the correct spatial and spin symmetry and can describe bond dissociation (of both single and multiple bonds) and electronic excitations in a formally and physically correct way. We then review the formulation of density functionals in terms of the on-top pair density. Finally we review successful applications of the theory to bond energies and bond dissociation potential energy curves of main-group and transition metal bonds, to barrier heights (including pericyclic reactions), to proton affinities, to the hydrogen bond energy of water dimer, to ground- and excited-state charge transfer, to valence and Rydberg excitations of molecules, and to singlet-triplet splittings of radicals. We find that that MC-PDFT can give accurate results not only with complete-active-space multiconfiguration wave functions but also with generalized-active-space multiconfiguration wave functions, which are practical for larger numbers of active electrons and active orbitals than are complete-active-space wave functions. The separated-pair approximation, which is a special case of generalized active space self-consistent-field theory, is especially promising. MC-PDFT, because it requires much less computer time and storage than pure WFT methods, has the potential to open larger and more complex strongly correlated systems to accurate simulation.

Electron pair density in the lowest 1Σ(u)(+) and 1Σ(g)(+) states of H2

We demonstrate and advocate the use of observable quantities derived from the two-electron reduced density matrix - pair densities, conditional densities, and exchange-correlation holes--as signatures of the type of electron correlation in a chemical bond. The prototype cases of the lowest (1)Σ(u)(+) and (1)Σ(g)(+) states of H(2), which exhibit large variation in types of bonding, ranging from strongly ionic to covalent, are discussed. Both the excited (1)Σ(g)(+) and (1)Σ(u)(+) states have been interpreted as essentially consisting of (natural) orbital configurations with an inner electron in a contracted 1sσ(g) orbital and an outer electron in a diffuse (united atom type, Rydberg) orbital. We show that nevertheless totally different correlation behavior is encountered in various states when comparing them at a common internuclear distance. Also when following one state along the internuclear distance coordinate, strong variation in correlation behavior is observed, as expected. Switches between ionic to covalent character of a state occur till very large distances (40 bohrs for states approaching the 1s3[script-l] asymptotic limit, and 282 bohrs for states approaching the 1s4[script-l] limit).