Correlated-Participating-Orbitals Pair-Density Functional Method and Application to Multiplet Energy Splittings of Main-Group Divalent Radicals

Universal Framework for Multiconfigurational DFT

Strong correlation remains a significant challenge for DFT with no satisfying solutions found yet within the standard Kohn-Sham framework. Instead, for decades, a number of different approaches have been suggested to combine the accuracy of multiconfigurational methods with the efficiency of DFT. In this article, we demonstrate that many of these methods are or would be significantly improved by being reformulated as variants of multiconfigurational pair-density functional theory (MC-PDFT). This work presents the first implementation of these methods within the recently proposed variational formulation of MC-PDFT. It also provides for the first time a systematic comparison of their accuracy across representative examples of strongly correlated systems. By analyzing their accuracy and formal properties, we provide design guidelines to inform the development of future functionals.

Reliable Diradical Characterization via Precise Singlet-Triplet Gap Calculations: Application to Thiele, Chichibabin, and Müller Analogous Diradicals

Accurately calculating the diradical character (y0) of molecular systems remains a significant challenge due to the scarcity of experimental data and the inherent multireference nature of the electronic structure. In this study, various quantum mechanical approaches, including broken symmetry density functional theory (BS-DFT), spin-flip time-dependent density functional theory (SF-TDDFT), mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT), complete active space self-consistent field (CASSCF), complete active space second-order perturbation theory (CASPT2), and multiconfigurational pair-density functional theory (MCPDFT), are employed to compute the singlet-triplet energy gaps (EST) and y0 values in Thiele, Chichibabin, and Müller analogous diradicals. By systematically comparing the results from these computational methods, we identify optimally tuned long-range corrected functional CAM-B3LYP in the BS-DFT framework as a most efficient method for accurately and affordably predicting both EST and y0 values. Additionally, our results demonstrate that (i) MRSF-TDDFT performs much better than SF-TDDFT; (ii) the MCPDFT method is robust in determining EST with minimal dependence on the choice of active space. These findings provide insight into the electronic structure and diradical character of the investigated molecules and highlight effective computational strategies for future studies in this domain. Thus, this work not only advances our understanding of diradical systems but also offers practical guidelines for their computational investigation.

A hybrid meta on-top functional for multiconfiguration pair-density functional theory

Multiconfiguration pair-density functional theory (MC-PDFT) was proposed a decade ago, but it is still in the early stage of density functional development. MC-PDFT uses functionals that are called on-top functionals; they depend on the density and the on-top pair density. Most MC-PDFT calculations to date have been unoptimized translations of generalized gradient approximations (GGAs) of Kohn-Sham density functional theory (KS-DFT). A hybrid MC-PDFT has also been developed, in which one includes a fraction of the complete active space self-consistent-field wave function energy in the total energy. Meta-GGA functionals, which use kinetic-energy densities in addition to GGA ingredients, have shown higher accuracy than GGAs in KS-DFT, yet the translation of meta-GGAs has not been previously proposed for MC-PDFT. In this paper, we propose a way to include kinetic energy density in a hybrid on-top functional for MC-PDFT, and we optimize the parameters of the resulting functional by training with a database developed as part of the present work that contains a wide variety of systems with diverse characters. The resulting hybrid meta functional is called the MC23 functional. We find that MC23 has improved performance as compared to KS-DFT functionals for both strongly and weakly correlated systems. We recommend MC23 for future MC-PDFT calculations.

Spin-symmetry-enforced solution of the many-body Schrödinger equation with a deep neural network

The integration of deep neural networks with the variational Monte Carlo (VMC) method has marked a substantial advancement in solving the Schrödinger equation. In this work we enforce spin symmetry in the neural-network-based VMC calculation using a modified optimization target. Our method is designed to solve for the ground state and multiple excited states with target spin symmetry at a low computational cost. It predicts accurate energies while maintaining the correct symmetry in strongly correlated systems, even in cases in which different spin states are nearly degenerate. Our approach also excels at spin-gap calculations, including the singlet-triplet gap in biradical systems, which is of high interest in photochemistry. Overall, this work establishes a robust framework for efficiently calculating various quantum states with specific spin symmetry in correlated systems.

SCINE-Software for chemical interaction networks

The software for chemical interaction networks (SCINE) project aims at pushing the frontier of quantum chemical calculations on molecular structures to a new level. While calculations on individual structures as well as on simple relations between them have become routine in chemistry, new developments have pushed the frontier in the field to high-throughput calculations. Chemical relations may be created by a search for specific molecular properties in a molecular design attempt, or they can be defined by a set of elementary reaction steps that form a chemical reaction network. The software modules of SCINE have been designed to facilitate such studies. The features of the modules are (i) general applicability of the applied methodologies ranging from electronic structure (no restriction to specific elements of the periodic table) to microkinetic modeling (with little restrictions on molecularity), full modularity so that SCINE modules can also be applied as stand-alone programs or be exchanged for external software packages that fulfill a similar purpose (to increase options for computational campaigns and to provide alternatives in case of tasks that are hard or impossible to accomplish with certain programs), (ii) high stability and autonomous operations so that control and steering by an operator are as easy as possible, and (iii) easy embedding into complex heterogeneous environments for molecular structures taken individually or in the context of a reaction network. A graphical user interface unites all modules and ensures interoperability. All components of the software have been made available as open source and free of charge.

Variational Pair-Density Functional Theory: Dealing with Strong Correlation at the Protein Scale

Multiconfigurational pair-density functional theory (MC-PDFT) offers a promising solution to the challenges faced by traditional density functional theory (DFT) in addressing molecular systems containing transition metals, open-shells, or strong correlations in general. By utilizing both the density and on-top pair-density, MC-PDFT can make use of a more flexible multiconfigurational wave function to capture the necessary static correlation, while the pair-density functional also includes the effect of dynamic correlation. So far, MC-PDFT has been used after a multiconfigurational self-consistent field (MCSCF) step, using the orbitals and configuration interaction coefficients from the converged MCSCF wave function to compute PDFT energies and properties. Here, instead, we propose to perform a direct optimization of the wave function using the pair-density functionals, resulting in a variational formulation of MC-PDFT. We derive the expressions for the wave function gradient and illustrate their similarity to standard MCSCF equations. Furthermore, we illustrate the accuracy on a set of singlet-triplet gaps as well as dissociation curves. Our findings highlight one of MC-PDFT's standout features: a reduced dependency on the active space size compared to conventional multiconfigurational wave function methodologies. Additionally, we show that the computational cost of MC-PDFT is potentially lower than MCSCF and often on-par with standard Kohn-Sham DFT, which is demonstrated by performing a MC-PDFT calculation of the entire ferredoxin protein with 1447 atoms and nearly 12 000 basis functions.

Molecules in Environments: Toward Systematic Quantum Embedding of Electrons and Drude Oscillators

We develop a quantum embedding method that enables accurate and efficient treatment of interactions between molecules and an environment, while explicitly including many-body correlations. The molecule is composed of classical nuclei and quantum electrons, whereas the environment is modeled via charged quantum harmonic oscillators. We construct a general Hamiltonian and introduce a variational Ansatz for the correlated ground state of the fully interacting molecule-environment system. This wave function is optimized via the variational Monte Carlo method and the ground state energy is subsequently estimated through the diffusion Monte Carlo method. The proposed scheme allows an explicit many-body treatment of electrostatic, polarization, and dispersion interactions between the molecule and the environment. We study solvation energies and excitation energies of benzene derivatives, obtaining excellent agreement with explicit ab initio calculations and experiments.

Variational Active Space Selection with Multiconfiguration Pair-Density Functional Theory

The selection of an adequate set of active orbitals for modeling strongly correlated electronic states is difficult to automate because it is highly dependent on the states and molecule of interest. Although many approaches have shown some success, no single approach has worked well in all cases. In light of this, we present the "discrete variational selection" (DVS) approach to active space selection, in which one generates multiple trial wave functions from a diverse set of systematically constructed active spaces and then selects between these wave functions variationally. We apply this DVS approach to 207 vertical excitations of small-to-medium-sized organic and inorganic molecules (with 3 to 18 atoms) in the QUESTDB database by (i) constructing various sets of active space orbitals through diagonalization of parametrized operators and (ii) choosing the result with the lowest average energy among the states of interest. This approach proves ineffective when variationally selecting between wave functions using the density matrix renormalization group (DMRG) or complete active space self-consistent field (CASSCF) energy but is able to provide good results when variationally selecting between wave functions using the energy of the translated PBE (tPBE) functional from multiconfiguration pair-density functional theory (MC-PDFT). Applying this DVS-tPBE approach to selection among state-averaged DMRG wave functions, we obtain a mean unsigned error of only 0.17 eV using hybrid MC-PDFT. This result matches that of our previous benchmark without the need to filter out poor active spaces and with no further orbital optimization following active space selection of the SA-DMRG wave functions. Furthermore, we find that DVS-tPBE is able to robustly and effectively select between the new SA-DMRG wave functions and our previous SA-CASSCF results.

Multiconfigurational Pair-Density Functional Theory Is More Complex than You May Think

Multiconfigurational pair-density functional theory (MC-PDFT) is a promising way to describe both strong and dynamic correlations in an inexpensive way. The functionals in MC-PDFT are often "translated" from standard spin density functionals. However, these translated functionals can in principle lead to "translated spin densities" with a nonzero imaginary component. Current developments so far neglect this imaginary part by simply setting it to zero. In this work, we show how this imaginary component is actually needed to reproduce the correct physical behavior in a range of cases, especially low-spin open shells. We showcase the resulting formalism on both local density approximation and generalized gradient approximation functionals and illustrate the numerical behavior by benchmarking a number of singlet-triplet splittings (ST gaps) of organic diradicals and low-lying excited states of some common organic molecules. The results demonstrate that this scheme improves existing translated functionals and gives more accurate results, even with minimal active spaces.

Galván I, Alavi A, Aleotti F, Aquilante F, Autschbach J, Avagliano D, Baiardi A, Bao JJ, Battaglia S, Birnoschi L, Blanco-González A, Bokarev SI, Broer R, Cacciari R, Calio PB, Carlson RK, Carvalho Couto R, Cerdán L, Chibotaru LF, Chilton NF, Church JR, Conti I, Coriani S, Cuéllar-Zuquin J, Daoud RE, Dattani N, Decleva P, de Graaf C, Delcey M, De Vico L, Dobrautz W, Dong SS, Feng R, Ferré N, Filatov(Gulak) M, Gagliardi L, Garavelli M, González L, Guan Y, Guo M, Hennefarth MR, Hermes MR, Hoyer CE, Huix-Rotllant M, Jaiswal VK, Kaiser A, Kaliakin DS, Khamesian M, King DS, Kochetov V, Krośnicki M, Kumaar AA, Larsson ED, Lehtola S, Lepetit MB, Lischka H, López Ríos P, Lundberg M, Ma D, Mai S, Marquetand P, Merritt ICD, Montorsi F, Mörchen M, Nenov A, Nguyen VHA, Nishimoto Y, Oakley MS, Olivucci M, Oppel M, Padula D, Pandharkar R, Phung QM, Plasser F, Raggi G, Rebolini E, Reiher M, Rivalta I, Roca-Sanjuán D, Romig T, Safari AA, Sánchez-Mansilla A, Sand AM, Schapiro I, Scott TR, Segarra-Martí J, Segatta F, Sergentu DC, Sharma P, Shepard R, Shu Y, Staab JK, Straatsma TP, Sørensen LK, Tenorio BNC, Truhlar DG, Ungur L, Vacher M, Veryazov V, Voß TA, Weser O, Wu D, Yang X, Yarkony D, Zhou C, Zobel JP, Lindh R. The OpenMolcas Web: A Community-Driven Approach to Advancing Computational Chemistry

The developments of the open-source OpenMolcas chemistry software environment since spring 2020 are described, with a focus on novel functionalities accessible in the stable branch of the package or via interfaces with other packages. These developments span a wide range of topics in computational chemistry and are presented in thematic sections: electronic structure theory, electronic spectroscopy simulations, analytic gradients and molecular structure optimizations, ab initio molecular dynamics, and other new features. This report offers an overview of the chemical phenomena and processes OpenMolcas can address, while showing that OpenMolcas is an attractive platform for state-of-the-art atomistic computer simulations.

Comprehensive Theoretical Study on Four Typical Intramolecular Hydrogen Shift Reactions of Peroxy Radicals: Multireference Character, Recommended Model Chemistry, and Kinetics

Intramolecular hydrogen shift reactions in peroxy radicals (RO₂^• → ^•QOOH) play key roles in the low-temperature combustion and in the atmospheric chemistry. In the present study, we found that a mild-to-moderate multireference character of a potential energy surface (PES) is widely present in four typical hydrogen shift reactions of peroxy radicals (RO₂^•, R = ethyl, vinyl, formyl methyl, and acetyl) by a systematic assessment based on the T₁ diagnostic, %TAE diagnostic, M diagnostic, and contribution of the dominant configuration of the reference CASSCF wavefunction (C₀²). To assess the effects of these inherent multireference characters on electronic structure calculations, we compared the PESs of the four reactions calculated by the multireference method CASPT2 in the complete basis set (CBS) limit, single-reference method CCSD(T)-F12, and single-reference-based composite method WMS. The results showed that ignoring the multireference character will introduce a mean unsigned deviation (MUD) of 0.46-1.72 kcal/mol from CASPT2/CBS results by using the CCSD(T)-F12 method or a MUD of 0.49-1.37 kcal/mol by WMS for three RO₂^• reactions (R = vinyl, formyl methyl, and acetyl) with a stronger multireference character. Further tests by single-reference Kohn-Sham (KS) density functional theory methods showed even larger deviations. Therefore, we specifically developed a new hybrid meta-generalized gradient approximation (GGA) functional M06-HS for the four typical H-shift reactions of peroxy radicals based on the WMS results for the ethyl peroxy radical reaction and on the CASPT2/CBS results for the others. The M06-HS method has an averaged MUD of 0.34 kcal/mol over five tested basis sets against the benchmark PESs, performing best in the tested 38 KS functionals. Last, in a temperature range of 200-3000 K, with the new functional, we calculated the high-pressure-limit rate coefficients of these H-shift reactions by the multi-structural variational transition-state theory with the small-curvature tunneling approximation (MS-CVT/SCT) and the thermochemical properties of all of the involved key radicals by the multi-structural torsional (MS-T) anharmonicity approximation method.

Automated Active Space Selection with Dipole Moments

Multireference calculations can provide accurate information of systems with strong correlation, which have increasing importance in the development of new molecules and materials. However, selecting a suitable active space for multireference calculations is nontrivial, and the selection of an unsuitable active space can sometimes lead to results that are not physically meaningful. Active space selection often requires significant human input, and the selection that leads to reasonable results often goes beyond chemical intuition. In this work, we have developed and evaluated two protocols for automated selection of the active space for multireference calculations based on a simple physical observable, the dipole moment, for molecules with nonzero ground-state dipole moments. One protocol is based on the ground-state dipole moment, and the other is based on the excited-state dipole moments. To evaluate the protocols, we constructed a dataset of 1275 active spaces from 25 molecules, each with 51 active space sizes considered, and have mapped out the relationship between the active space, dipole moments, and vertical excitation energies. We have demonstrated that, within this dataset, our protocols allow one to choose among a number of accessible active spaces one that is likely to give reasonable vertical excitation energies, especially for the first three excitations, with no parameters manually decided by the user. We show that, with large active spaces removed from consideration, the accuracy is similar and the time-to-solution can be reduced by more than 10 fold. We also show that the protocols can be applied to potential energy surface scans and determining the spin states of transition metal oxides.

Chemical Bonding in Isoelectronic NdO2 and SmO22

Neodymium dioxide (NdO₂) and samarium dioxide cation (SmO₂²⁺) are isoelectronic molecules. Here we used calculations of the spin-orbit-free wave functions to study and compare their geometries, spin states, and bonding. We used Kohn-Sham density functional theory with the B97-1 exchange-correlation functional to optimize the geometries and found that the two molecules have different ground spin states and structures. NdO₂ favors a linear ONdO triplet structure, and SmO₂²⁺ favors a linear SmOO²⁺ quintet structure. We then used state-averaged complete-active-space self-consistent-field (SA-CASSCF) calculations to investigate the bonding characteristics of NdO₂ and SmO₂²⁺ in various geometries. We found that in NdOO, one electron is transferred from Nd to O, while in SmO₂²⁺, there is no electron transfer between Sm and O. The SA-CASSCF calculation also shows that ONdO has a stronger bonding orbital between a 4f orbital of Nd and a p_z orbital of oxygen atoms. We compared three multireference methods, namely, extended multistate complete active space second-order perturbation theory (XMS-CASPT2), extended multistate pair-density functional theory (XMS-PDFT), and compressed multistate pair-density functional theory (CMS-PDFT), for calculating the spin-orbit-free energies of various isomers of both molecules. We found that although XMS-PDFT and CMS-PDFT are at the same cost level as SA-CASSCF, they give results with the same accuracy as given for the much more demanding XMS-CASPT2 calculation. Between the two multistate PDFT methods, CMS-PDFT is better at giving good degeneracies for states that should be degenerate.

6,6'-Biindeno[1,2-b]anthracene: An Open-Shell Biaryl with High Diradical Character

We report in situ generation of a 6,6'-biindeno[1,2-b]anthracene (BIA) derivative as an open-shell biaryl with high diradical character, which could be identified by mass spectrometry, NMR spectroscopy, single-crystal X-ray analysis, UV-vis-NIR absorption spectroscopy, and electron paramagnetic resonance (EPR) spectroscopy. Theoretical calculations by various methods and variable-temperature EPR analyses were performed to tackle the elusive ground state of BIA diradical, suggesting a singlet ground state with a nearly degenerate triplet state. These results provide insight into the design of unique open-shell biaryls.

Large-Scale Benchmarking of Multireference Vertical-Excitation Calculations via Automated Active-Space Selection

We have calculated state-averaged complete-active-space self-consistent-field (SA-CASSCF), multiconfiguration pair-density functional theory (MC-PDFT), hybrid MC-PDFT (HMC-PDFT), and n-electron valence state second-order perturbation theory (NEVPT2) excitation energies with the approximate pair coefficient (APC) automated active-space selection scheme for the QUESTDB benchmark database of 542 vertical excitation energies. We eliminated poor active spaces (20-40% of calculations) by applying a threshold to the SA-CASSCF absolute error. With the remaining calculations, we find that NEVPT2 performance is significantly impacted by the size of the basis set the wave functions are converged in, regardless of the quality of their description, which is a problem absent in MC-PDFT. Additionally, we find that HMC-PDFT is a significant improvement over MC-PDFT with the translated PBE (tPBE) density functional and that it performs about as well as NEVPT2 and second-order coupled cluster on a set of 373 excitations in the QUESTDB database. We optimized the percentage of SA-CASSCF energy to include in HMC-PDFT when using the tPBE on-top functional, and we find the 25% value used in tPBE0 to be optimal. This work is by far the largest benchmarking of MC-PDFT and HMC-PDFT to date, and the data produced in this work are useful as a validation of HMC-PDFT and of the APC active-space selection scheme. We have made all the wave functions produced in this work (orbitals and CI vectors) available to the public and encourage the community to utilize this data as a tool in the development of further multireference model chemistries.

Dual-Level Training of Gaussian Processes with Physically Inspired Priors for Geometry Optimizations

Gaussian process (GP) regression has been recently developed as an effective method in molecular geometry optimization. The prior mean function is one of the crucial parts of the GP. We design and validate two types of physically inspired prior mean functions: force-field-based priors and posterior-type priors. In this work, we implement a dual-level training (DLT) optimizer for the posterior-type priors. The DLT optimizers can be considered as a class of optimization algorithms that belong to the delta-machine learning paradigm but with several major differences compared to the previously proposed algorithms in the same paradigm. In the first level of the DLT, we incorporate the classical mechanical descriptions of the equilibrium geometries into the prior function, which enhances the performance of the GP optimizer as compared to the one using a constant (or zero) prior. In the second level, we utilize the surrogate potential energy surfaces (PESs), which incorporate the physics learned in the first-level training, as the prior function to refine the model performance further. We find that the force-field-based priors and posterior-type priors reduce the overall optimization steps by a factor of 2-3 when compared to the limited-memory Broyden-Fletcher-Goldfarb-Shanno (L-BFGS) optimizer as well as the constant-prior GP optimizer proposed in previous works. We also demonstrate the potential of recovering the real PESs with GP with a force-field prior. This work shows the importance of including domain knowledge as an ingredient in the GP, which offers a potentially robust learning model for molecular geometry optimization and for exploring molecular PESs.

Elucidating the molecular orbital dependence of the total electronic energy in multireference problems

The accurate resolution of the chemical properties of strongly correlated systems, such as biradicals, requires the use of electronic structure theories that account for both multi-reference and dynamic correlation effects. A variety of methods exist that aim to resolve the dynamic correlation in multi-reference problems, commonly relying on an exponentially scaling complete-active-space self-consistent-field (CASSCF) calculation to generate reference molecular orbitals (MOs). However, while CASSCF orbitals provide the optimal solution for a selected set of correlated (active) orbitals, their suitability in the quest for the resolution of the total correlation energy has not been thoroughly investigated. Recent research has shown the ability of Kohn-Shan density functional theory to provide improved orbitals for coupled cluster (CC) and Møller-Plesset perturbation theory (MP) calculations. Here, we extend the search for optimal and more cost effective MOs to post-configuration-interaction [post-(CI)] methods, surveying the ability of the MOs obtained with various density functional theory (DFT) functionals, as well as Hartree-Fock and CC and MP calculations to accurately capture the total electronic correlation energy. Applying the anti-Hermitian contracted Schrödinger equation to the dissociation of N₂, the calculation of biradical singlet-triplet gaps, and the transition states of bicylobutane isomerization, we demonstrate that DFT provides a cost-effective alternative to CASSCF in providing reference orbitals for post-CI dynamic correlation calculations.

Electronic structure of strongly correlated systems: recent developments in multiconfiguration pair-density functional theory and multiconfiguration nonclassical-energy functional theory

Strong electron correlation plays an important role in transition-metal and heavy-metal chemistry, magnetic molecules, bond breaking, biradicals, excited states, and many functional materials, but it provides a significant challenge for modern electronic structure theory. The treatment of strongly correlated systems usually requires a multireference method to adequately describe spin densities and near-degeneracy correlation. However, quantitative computation of dynamic correlation with multireference wave functions is often difficult or impractical. Multiconfiguration pair-density functional theory (MC-PDFT) provides a way to blend multiconfiguration wave function theory and density functional theory to quantitatively treat both near-degeneracy correlation and dynamic correlation in strongly correlated systems; it is more affordable than multireference perturbation theory, multireference configuration interaction, or multireference coupled cluster theory and more accurate for many properties than Kohn–Sham density functional theory. This perspective article provides a brief introduction to strongly correlated systems and previously reviewed progress on MC-PDFT followed by a discussion of several recent developments and applications of MC-PDFT and related methods, including localized-active-space MC-PDFT, generalized active-space MC-PDFT, density-matrix-renormalization-group MC-PDFT, hybrid MC-PDFT, multistate MC-PDFT, spin–orbit coupling, analytic gradients, and dipole moments. We also review the more recently introduced multiconfiguration nonclassical-energy functional theory (MC-NEFT), which is like MC-PDFT but allows for other ingredients in the nonclassical-energy functional. We discuss two new kinds of MC-NEFT methods, namely multiconfiguration density coherence functional theory and machine-learned functionals.

This feature article overviews recent work on active spaces, matrix product reference states, treatment of quasidegeneracy, hybrid theory, density-coherence functionals, machine-learned functionals, spin–orbit coupling, gradients, and dipole moments.

Multiconfiguration Pair-Density Functional Theory for Transition Metal Silicide Bond Dissociation Energies, Bond Lengths, and State Orderings

Transition metal silicides are promising materials for improved electronic devices, and this motivates achieving a better understanding of transition metal bonds to silicon. Here we model the ground and excited state bond dissociations of VSi, NbSi, and TaSi using a complete active space (CAS) wave function and a separated-pair (SP) wave function combined with two post-self-consistent field techniques: complete active space with perturbation theory at second order and multiconfiguration pair-density functional theory. The SP approximation is a multiconfiguration self-consistent field method with a selection of configurations based on generalized valence bond theory without the perfect pairing approximation. For both CAS and SP, the active-space composition corresponds to the nominal correlated-participating-orbital scheme. The ground state and low-lying excited states are explored to predict the state ordering for each molecule, and potential energy curves are calculated for the ground state to compare to experiment. The experimental bond dissociation energies of the three diatomic molecules are predicted with eight on-top pair-density functionals with a typical error of 0.2 eV for a CAS wave function and a typical error of 0.3 eV for the SP approximation. We also provide a survey of the accuracy achieved by the SP and extended separated-pair approximations for a broader set of 25 transition metal–ligand bond dissociation energies.

A Ranked-Orbital Approach to Select Active Spaces for High-Throughput Multireference Computation

The past decade has seen a great increase in the application of high-throughput computation to a variety of important problems in chemistry. However, one area which has been resistant to the high-throughput approach is multireference wave function methods, in large part due to the technicalities of setting up these calculations and in particular the not always intuitive challenge of active space selection. As we look toward a future of applying high-throughput computation to all areas of chemistry, it is important to prepare these methods for large-scale automation. Here, we propose a ranked-orbital approach to select active spaces with the goal of standardizing multireference methods for high-throughput computation. This method allows for the meaningful comparison of different active space selection schemes and orbital localizations, and we demonstrate the utility of this approach across 1120 multireference calculations for the excitation energies of small molecules. Our results reveal that it is helpful to distinguish the method used to generate orbitals from the method of ranking orbitals in terms of importance for the active space. Additionally, we propose our own orbital ranking scheme that estimates the importance of an orbital for the active space through a pair-interaction framework from orbital energies and features of the Hartree-Fock exchange matrix. We call this new scheme the "approximate pair coefficient" (APC) method and we show that it performs quite well for the test systems presented.

Multiconfiguration Pair-Density Functional Theory

Kohn-Sham density functional theory with the available exchange-correlation functionals is less accurate for strongly correlated systems, which require a multiconfigurational description as a zero-order function, than for weakly correlated systems, and available functionals of the spin densities do not accurately predict energies for many strongly correlated systems when one uses multiconfigurational wave functions with spin symmetry. Furthermore, adding a correlation functional to a multiconfigurational reference energy can lead to double counting of electron correlation. Multiconfiguration pair-density functional theory (MC-PDFT) overcomes both obstacles, the second by calculating the quantum mechanical part of the electronic energy entirely by a functional, and the first by using a functional of the total density and the on-top pair density rather than the spin densities. This allows one to calculate the energy of strongly correlated systems efficiently with a pair-density functional and a suitable multiconfigurational reference function. This article reviews MC-PDFT and related background information.

Accurate singlet-triplet gaps in biradicals via the spin averaged anti-Hermitian contracted Schrödinger equation

The accurate description of biradical systems, and in particular the resolution of their singlet-triplet gaps, has long posed a major challenge to the development of electronic structure theories. Biradicaloid singlet ground states are often marked by strong correlation and, hence, may not be accurately treated by mainstream, single-reference methods such as density functional theory or coupled cluster theory. The anti-Hermitian contracted Schrödinger equation (ACSE), whose fundamental quantity is the two-electron reduced density matrix rather than the N-electron wave function, has previously been shown to account for both dynamic and strong correlations when seeded with a strongly correlated guess from a complete active space (CAS) calculation. Here, we develop a spin-averaged implementation of the ACSE, allowing it to treat higher multiplicity states from the CAS input without additional state preparation. We apply the spin-averaged ACSE to calculate the singlet-triplet gaps in a set of small main group biradicaloids, as well as the organic four-electron biradicals trimethylenemethane and cyclobutadiene, and naphthalene, benchmarking the results against other state-of-the-art methods reported in the literature.

Toward a Resolution of the Static Correlation Problem in Density Functional Theory from Semidefinite Programming

Kohn-Sham density functional theory (DFT) has long struggled with the accurate description of strongly correlated and open shell systems, and improvements have been minor even in the newest hybrid functionals. In this Letter we treat the static correlation in DFT when frontier orbitals are degenerate by the means of using a semidefinite programming (SDP) approach to minimize the system energy as a function of the N-representable, non-idempotent 1-electron reduced density matrix. While showing greatly improved singlet-triplet gaps for local density approximation and generalized gradient approximation (GGA) functionals, the SDP procedure reveals flaws in modern meta and hybrid GGA functionals, which show no major improvements when provided with an accurate electron density.

In silico prediction of annihilators for triplet-triplet annihilation upconversion via auxiliary-field quantum Monte Carlo

The energy of the lowest-lying triplet state (T1) relative to the ground and first-excited singlet states (S0, S1) plays a critical role in optical multiexcitonic processes of organic chromophores. Focusing on triplet-triplet annihilation (TTA) upconversion, the S0 to T1 energy gap, known as the triplet energy, is difficult to measure experimentally for most molecules of interest. Ab initio predictions can provide a useful alternative, however low-scaling electronic structure methods such as the Kohn-Sham and time-dependent variants of Density Functional Theory (DFT) rely heavily on the fraction of exact exchange chosen for a given functional, and tend to be unreliable when strong electronic correlation is present. Here, we use auxiliary-field quantum Monte Carlo (AFQMC), a scalable electronic structure method capable of accurately describing even strongly correlated molecules, to predict the triplet energies for a series of candidate annihilators for TTA upconversion, including 9,10 substituted anthracenes and substituted benzothiadiazole (BTD) and benzoselenodiazole (BSeD) compounds. We compare our results to predictions from a number of commonly used DFT functionals, as well as DLPNO-CCSD(T0), a localized approximation to coupled cluster with singles, doubles, and perturbative triples. Together with S1 estimates from absorption/emission spectra, which are well-reproduced by TD-DFT calculations employing the range-corrected hybrid functional CAM-B3LYP, we provide predictions regarding the thermodynamic feasibility of upconversion by requiring (a) the measured T1 of the sensitizer exceeds that of the calculated T1 of the candidate annihilator, and (b) twice the T1 of the annihilator exceeds its S1 energetic value. We demonstrate a successful example of in silico discovery of a novel annihilator, phenyl-substituted BTD, and present experimental validation via low temperature phosphorescence and the presence of upconverted blue light emission when coupled to a platinum octaethylporphyrin (PtOEP) sensitizer. The BTD framework thus represents a new class of annihilators for TTA upconversion. Its chemical functionalization, guided by the computational tools utilized herein, provides a promising route towards high energy (violet to near-UV) emission.

Predicting Bond Dissociation Energies and Bond Lengths of Coordinatively Unsaturated Vanadium-Ligand Bonds

Understanding the electronic structure of coordinatively unsaturated transition-metal compounds and predicting their physical properties are of great importance for catalyst design. Bond dissociation energy D_e and bond length r_e are two of the fundamental quantities for which good predictions are important for a successful design strategy. In the present work, recent experimentally measured bond energies and bond lengths of VX diatomic molecules (X = C, N, S) are used as a gauge to consider the utility of a number of electronic structure methods. Single-reference methods are one focus because of their efficiency and utility in practical calculations, and multireference configuration interaction (MRCISD) methods and a composite coupled cluster (CCC) method are a second focus because of their potential high accuracy. The comparison is especially challenging because of the large multireference M diagnostics of these molecules, in the range 0.15-0.19. For the single-reference methods, Kohn-Sham density functional theory (KS-DFT) has been tested with a variety of approximate exchange-correlation functionals. Of these, MOHLYP provides the bond dissociation energies in best agreement with experiments, and BLYP provides the bond lengths that are in best agreement with experiments; but by requiring good performance for both the D_e and r_e of the vanadium compounds, MOHLYP, MN12-L, MGGA_MS1, MGGA_MS0, O3LYP, and M06-L are the most highly recommended functionals. The CCC calculations include up to connected pentuple excitations for the valence electrons and up to connected quadruple excitations for the core-valence terms; this results in highly accurate dissociation energies and good bond lengths. Averaged over the three molecules, the mean unsigned deviation of CCC bond energies from experimental ones is only 0.4 kcal/mol, demonstrating excellent convergence of theory and experiments.

Variational Multistate Density Functional Theory for a Balanced Treatment of Static and Dynamic Correlations

Here, an approach to variational multistate density functional theory (vMSDFT) is explored. In this approach, the Kohn-Sham orbitals as well as configuration coefficients were simultaneously optimized, thus yielding a full variational minimum. Furthermore, this work also proposes two important improvements on the MSDFT framework. First, a "point-to-point correction" is used to correct the static correlation present in the DFT framework. Therefore, double counting of static correlation in vMSDFT is mitigated. Second, a general form to construct the transition density functional in the vMSDFT framework is proposed, which allows for the properties of vMSDFT wave functions to be standardized to the complete active space self-consistent field properties. The utility of vMSDFT is illustrated on molecular systems of interest including bond breaking, diradicals, excited states, and conical intersections. The numerical results suggest that the accuracy of vMSDFT is in close agreement with the high-level multireference methods.

Analytic gradients for state-averaged multiconfiguration pair-density functional theory

Analytic gradients are important for efficient calculations of stationary points on potential energy surfaces, for interpreting spectroscopic observations, and for efficient direct dynamics simulations. For excited electronic states, as are involved in UV-Vis spectroscopy and photochemistry, analytic gradients are readily available and often affordable for calculations using a state-averaged complete active space self-consistent-field (SA-CASSCF) wave function. However, in most cases, a post-SA-CASSCF step is necessary for quantitative accuracy, and such calculations are often too expensive if carried out by perturbation theory or configuration interaction. In this work, we present the analytic gradients for multiconfiguration pair-density functional theory based on SA-CASSCF wave functions, which is a more affordable alternative. A test set of molecules has been studied with this method, and the stationary geometries and energetics are compared to values in the literature as obtained by other methods. Excited-state geometries computed with state-averaged pair-density functional theory have similar accuracy to those from complete active space perturbation theory at the second-order.

Spin Splitting Energy of Transition Metals: A New, More Affordable Wave Function Benchmark Method and Its Use to Test Density Functional Theory

Accurately predicting the spin splitting energy of chemical species is important for understanding their reactivity and magnetic properties, but it is very challenging, especially for molecules containing transition metals. One impediment to progress is the scarcity of accurate benchmark data. Here we report a set of calculations designed to yield reliable benchmarks for simple transition-metal complexes that can be used to test density functional methods that are affordable for large systems of more practical interest. Various wave function methods are tested against experiment for Fe²⁺, Fe³⁺, and Co³⁺, including CASSCF, CASPT2, CASPT3, MRCISD, MRCISD+Q, ACPF, AQCC, CCSD(T), and CASPT2/CCSD(T) and also a new method called CASPT2.5, which is performed by taking the average of the CASPT2 and CASPT3 energies. We find that MRCISD+Q, ACPF, and AQCC require smaller active spaces for good accuracy than are required by CASPT2 and CASPT3, and this aspect may be important for calculations on larger molecules; here we find that CASPT2.5 extrapolated to a complete basis set is the most suitable method-in terms of computational cost and in terms of accuracy on monatomic systems-and therefore we chose this method for molecular benchmarks. Then Kohn-Sham density functional calculations with 60 exchange-correlation functionals are tested for FeF₂, FeCl₂, and CoF₂. We find that MN15-L, M06-SX, and revM06 have very good agreement with CASPT2.5 benchmarks in terms of both the spin splitting energy and the optimized geometry for each spin state. In addition, we recommend def2-TZVP as the most suitable basis set to perform density functional calculations for molecular spin splitting energies; extra polarization functions in the basis set do not help to increase the accuracy of the spin splitting energy in KS calculations.

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.

Weak Interactions in Alkaline Earth Metal Dimers by Pair-Density Functional Theory

Alkaline earth dimers have small bond energies (less than 5 kcal/mol) that provide a difficult challenge for electronic structure calculations. They are especially challenging for Kohn-Sham density functional theory (KS-DFT) using generalized gradient approximations (GGAs) as the exchange-correlation density functional because GGAs often do not provide accurate results for weak interactions. Here we treat alkaline earth dimers from six different rows of the periodic table. We show that the dominant correlating configurations are the same in all six dimers. We also show that multiconfiguration pair-density functional theory (MC-PDFT) using a fully translated GGA as the on-top density functional not only performs much better than KS-DFT with GGAs in predicting equilibrium distances and dissociation energies but also performs better than the more computationally demanding complete active space second-order perturbation theory (CASPT2) with large basis sets and performs even better than CASPT2 with smaller basis sets.

How well can density functional theory and pair-density functional theory predict the correct atomic charges for dissociation and accurate dissociation energetics of ionic bonds?

The accuracy of density functional theory (DFT) is often judged by predicted dissociation energies, but one should also consider charge densities as illustrated here for dissociation of heteronuclear diatomic molecules, including ionic bonds for which local density functionals yield erroneous results. Some hybrid density functionals with 100% exact exchange in Kohn-Sham DFT and the local functionals in multiconfiguration pair-density functional theory give relatively acurate dissociation energies for NaCl, and they correctly yield uncharged dissociated atoms.

Combining Wave Function Methods with Density Functional Theory for Excited States

We review state-of-the-art electronic structure methods based both on wave function theory (WFT) and density functional theory (DFT). Strengths and limitations of both the wave function and density functional based approaches are discussed, and modern attempts to combine these two methods are presented. The challenges in modeling excited-state chemistry using both single-reference and multireference methods are described. Topics covered include background, combining density functional theory with single-configuration wave function theory, generalized Kohn-Sham (KS) theory, global hybrids, range-separated hybrids, local hybrids, using KS orbitals in many-body theory (including calculations of the self-energy and the GW approximation), Bethe-Salpeter equation, algorithms to accelerate GW calculations, combining DFT with multiconfigurational WFT, orbital-dependent correlation functionals based on multiconfigurational WFT, building multiconfigurational wave functions from KS configurations, adding correlation functionals to multiconfiguration self-consistent-field (MCSCF) energies, combining DFT with configuration-interaction singles by means of time-dependent DFT, using range separation to combine DFT with MCSCF, embedding multiconfigurational WFT in DFT, and multiconfiguration pair-density functional theory.

Self-Interaction Error in Density Functional Theory: An Appraisal

Self-interaction error (SIE) is considered to be one of the major sources of error in most approximate exchange-correlation functionals for Kohn-Sham density-functional theory (KS-DFT), and it is large with all local exchange-correlation functionals and with some hybrid functionals. In this work, we consider systems conventionally considered to be dominated by SIE. For these systems, we demonstrate that by using multiconfiguration pair-density functional theory (MC-PDFT), the error of a translated local density-functional approximation is significantly reduced (by a factor of 3) when using an MCSCF density and on-top density, as compared to using KS-DFT with the parent functional; the error in MC-PDFT with local on-top functionals is even lower than the error in some popular KS-DFT hybrid functionals. Density-functional theory, either in MC-PDFT form with local on-top functionals or in KS-DFT form with some functionals having 50% or more nonlocal exchange, has smaller errors for SIE-prone systems than does CASSCF, which has no SIE.

Multiconfiguration Pair-Density Functional Theory and Complete Active Space Second Order Perturbation Theory. Bond Dissociation Energies of FeC, NiC, FeS, NiS, FeSe, and NiSe

We investigate the performance of multiconfiguration pair-density functional theory (MC-PDFT) and complete active space second-order perturbation theory for computing the bond dissociation energies of the diatomic molecules FeC, NiC, FeS, NiS, FeSe, and NiSe, for which accurate experimental data have become recently available [Matthew, D. J.; Tieu, E.; Morse, M. D. J. Chem. Phys. 2017, 146, 144310-144320]. We use three correlated participating orbital (CPO) schemes (nominal, moderate, and extended) to define the active spaces, and we consider both the complete active space (CAS) and the separated-pair (SP) schemes to specify the configurations included for a given active space. We found that the moderate SP-PDFT scheme with the tPBE on-top density functional has the smallest mean unsigned error (MUE) of the methods considered. This level of theory provides a balanced treatment of the static and dynamic correlation energies for the studied systems. This is encouraging because the method is low in cost even for much more complicated systems.

Multiconfiguration Pair-Density Functional Theory Is Free From Delocalization Error

Delocalization error has been singled out by Yang and co-workers as the dominant error in Kohn-Sham density functional theory (KS-DFT) with conventional approximate functionals. In this Letter, by computing the vertical first ionization energy for well separated He clusters, we show that multiconfiguration pair-density functional theory (MC-PDFT) is free from delocalization error. To put MC-PDFT in perspective, we also compare it with some Kohn-Sham density functionals, including both traditional and modern functionals. Whereas large delocalization errors are almost universal in KS-DFT (the only exception being the very recent corrected functionals of Yang and co-workers), delocalization error is removed by MC-PDFT, which bodes well for its future as a step forward from KS-DFT.

Spin-Multiplet Components and Energy Splittings by Multistate Density Functional Theory

Kohn-Sham density functional theory has been tremendously successful in chemistry and physics. Yet, it is unable to describe the energy degeneracy of spin-multiplet components with any approximate functional. This work features two contributions. (1) We present a multistate density functional theory (MSDFT) to represent spin-multiplet components and to determine multiplet energies. MSDFT is a hybrid approach, taking advantage of both wave function theory and density functional theory. Thus, the wave functions, electron densities and energy density-functionals for ground and excited states and for different components are treated on the same footing. The method is illustrated on valence excitations of atoms and molecules. (2) Importantly, a key result is that for cases in which the high-spin components can be determined separately by Kohn-Sham density functional theory, the transition density functional in MSDFT (which describes electronic coupling) can be defined rigorously. The numerical results may be explored to design and optimize transition density functionals for configuration coupling in multiconfigurational DFT.

Comparative Study of Nonhybrid Density Functional Approximations for the Prediction of 3d Transition Metal Thermochemistry

The utility of several nonhybrid density functional approximations (DFAs) is considered for the prediction of gas phase enthalpies of formation for a large set of 3d transition metal-containing molecules. Nonhybrid DFAs can model thermochemical values for 3d transition metal-containing molecules with accuracy comparable to that of hybrid functionals. The GAM-generalized gradient approximation (GGA); the TPSS, M06-L, and MN15-L meta-GGAs; and the Rung 3.5 PBE+ΠLDA(s) DFAs all give root-mean-square deviations below that of the widely used B3LYP hybrid. Modern nonhybrid DFAs continue to show utility for transition metal thermochemistry.