Multi-state pair-density 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.

Accurate quantum chemical methods for the prediction of spin-state energy gaps for strongly correlated systems are computationally expensive and scale poorly with the size of the system. This makes calculations for many experimentally interesting molecules impractical even with abundant computational resources. In previous work, we have shown that the localized active space (LAS) self-consistent field (SCF) method is an efficient way to obtain multi-configuration SCF wave functions of comparable quality to the corresponding complete active space (CAS) ones. To obtain quantitative results, a post-SCF method is needed to estimate the complete correlation energy. One such method is multiconfiguration pair-density functional theory (PDFT), which calculates the energy based on the density and on-top pair density obtained from a multiconfiguration wave function. In this work we introduce localized-active-space pair-density functional theory, which uses a LAS wave function for subsequent PDFT calculations. The method is tested for computing spin-state energy gaps in conjugated organic molecules and bimetallic compounds and is shown to give results within 0.05 eV of the corresponding CAS-PDFT results at a significantly lower cost.

We analyzed static and dynamic electron correlation by decomposing the total electronic energy of calculations by restricted Hartree-Fock theory, complete active-space self-consistent field (CASSCF) theory, and multireference configuration interaction (MRCI). We used three different schemes to break down the relative energy contributions to the potential energy curves for the dissociation of H2, F2, and N2. The first decomposition scheme involves the classical and nonclassical components of the energy. The second and third recognize the part of the energy that is not expressible in terms of the one-body reduced density matrix; this is called the connected energy. The unconnected component is further decomposed into a part calculable from the density and the part calculable from the density coherence. The first decomposition scheme shows that the sum of the one-electron energy and the classical two-electron energy contains a negligible portion of the static correlation. This quantity has a relatively small variance between the three levels, especially for CASSCF and MRCI. This provides an explanation of why multiconfiguration pair-density functional theory and multiconfiguration density-coherence functional theory are able to improve the CASSCF energy. The latter two decompositions show that the connected energy contains a significant portion of static correlation. The energy representable by either the density or the density coherence is significantly different at the three levels. Mixing the density and density coherence between different methods may lead to a systematic error in the bond dissociation energy and the equilibrium bond distance, indicating that the density energy component and the density coherence energy component both include a significant amount of both static and dynamic correlation. These wave function decompositions can be useful for developing new functionals for density functional theory, density-coherence functional theory, density matrix functional theory, and pair-density functional theory and for guiding expectations for these theories.

Spin-projected unrestricted Hartree-Fock (SUHF) theory is a valuable method that effectively addresses static correlation. To further enhance its accuracy, it is important to augment it with dynamic correlation. Based on SUHF theory, we propose a pair-density functional theory, namely, SU-PDFT, which formally follows the concept of multiconfiguration pair-density functional theory (MC-PDFT). SU-PDFT shows comparable accuracy to MC-PDFT but avoids the issue of the exponential scaling of the multiconfiguration (MC) approach. The hybrid version of SU-PDFT with one or two parameters is also developed, displaying further improved accuracy.

Accurate quantum chemical methods for the prediction of spin-state energy gaps for strongly correlated systems are computationally expensive and scale poorly with the size of the system. This makes calculations for many experimentally interesting molecules impractical even with abundant computational resources. Previous work has shown that the localized active space (LAS) self-consistent field (SCF) method can be an efficient way to obtain multiconfiguration SCF wave functions of comparable quality to the corresponding complete active space (CAS) ones. To obtain quantitative results, a post-SCF method is needed to estimate the complete correlation energy. One such method is multiconfiguration pair-density functional theory (PDFT), which calculates the energy based on the density and on-top pair density obtained from a multiconfiguration wave function. In this work, we introduce localized-active-space PDFT, which uses a LAS wave function for subsequent PDFT calculations. The method is tested by computing spin-state energies and gaps in conjugated organic molecules and a bimetallic compound and comparing to the corresponding CAS-PDFT values.

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.

Transition states for Diels-Alder reactions are strongly correlated, as evidenced by high-to-very-high M diagnostics, and therefore they require treatment by multireference methods. Multiconfiguration pair-density functional theory (MC-PDFT) combines a multiconfiguration wave function with a functional of the electron density and the on-top pair density to calculate the electronic energy for strongly correlated systems at a much lower cost than wave function methods that do not employ density functionals. Here we apply MC-PDFT to the Diels-Alder cycloaddition reaction of 1,3-butadiene with ethylene, where two kinds of reaction paths have been widely studied: concerted synchronous paths and diradical stepwise paths. The lowest-energy reaction path is now known to be a concerted synchronous one, and a method's ability to predict this is an important test. By comparison to the best available theoretical results in the literature, we test the accuracy of MC-PDFT with several choices of on-top functional for geometries and enthalpies of stable structures along both paths and for the transition state geometries. We also calculate the Arrhenius activation energies for both paths and compare these to experiment. We also compare to Kohn-Sham density functional theory (KS-DFT) with selected exchange-correlation functionals. CAS-PDFT gives consistently good energies and geometries for both the concerted and stepwise mechanisms, but none of the KS-DFT functionals gives accurate activation energies for both. The stepwise transition state is very strongly correlated, and MC-PDFT can treat it, but KS-DFT (which involves a single-configuration treatment) has larger errors. The results confirm that using a multiconfigurational reference function for strongly correlated transition states can significantly improve the reliability and that MC-PDFT can provide good accuracy at a much lower computational cost than competing multireference methods.

The use of flavins and flavoproteins in photocatalytic, sensing, and biotechnological applications has led to a growing interest in computationally modeling the excited-state electronic structure and photophysics of flavin. However, there is limited consensus regarding which computational methods are appropriate for modeling flavin's photophysics. We compare the energies of low-lying excited states of flavin computed with time-dependent density functional theory (TD-DFT), equation-of-motion coupled cluster (EOM-EE-CCSD), scaled opposite-spin configuration interaction [SOS-CIS(D)], multiconfiguration pair-density functional theory (MC-PDFT), and several multireference perturbation theory (MR-PT2) methods. In the first part, we focus on excitation energies of the first singlet excited state (S1) of five different redox and protonation states of flavin, with the goal of finding a suitable active space for MR-PT2 calculations. In the second part, we construct two sets of one-dimensional potential energy surfaces connecting the S0 and S1 equilibrium geometries (S0-S1 path) and the S1 (π,π*) and S2 (n,π*) equilibrium geometries (S1-S2 path). The first path therefore follows a Franck-Condon active mode of flavin while the second path maps crossings points between low-lying singlet and triplet states in flavin. We discuss the similarities and differences in the TD-DFT, EOM-EE-CCSD, SOS-CIS(D), MC-PDFT and MR-PT2 energy profiles along these paths. We find that (TD-)DFT methods are suitable for applications such as simulating the spectra of flavins but are inconsistent with several other methods when used for some geometry optimizations and when describing the energetics of dark (n,π*) states. MR-PT2 methods show promise for the simulation of flavin's low-lying excited states, but the selection of orbitals for the active space and the number of roots used for state averaging must be done carefully to avoid artifacts. Some properties, such as the intersystem crossing geometry and energy between the S1 (π,π*) and T2 (n,π*) states, may require additional benchmarking before they can be determined quantitatively.

In the present work we propose two distinct state-averaging (SA)-based methodologies for the calculation of excited states, in conjunction with the long-range complete active space self-consistent field (CASSCF) short-range density functional theory (DFT) approach (CAS-srDFT). The state-specific density ansatz, termed SA-CAS-srDFT, initially determines the variational parameters of an approximate srDFT functional that operates with state-averaged densities. Subsequent to convergence, the CAS-srDFT energies of each state are computed from the state-specific one- and two-body densities. The second approach is termed configuration interaction (CI)-srDFT, for which a first-order correction is added to the approximate SA density CAS-srDFT functionals. Unlike the state-specific density approach SA-CAS-srDFT, diagonalisation of the first-order corrected effective Hamiltonian CI matrix yields orthonormal CI solutions for every state. In both approaches, the total one-body and on-top pair density (OTPD) was employed for the final energy evaluation. It was observed that the CI-srDFT approach gives physically correct potential curves for ethylene, in contrast to SA-CAS-srDFT. Moreover, the CI-srDFT approach demonstrates a reduced dependence of excitation energies on the number of states in the average when compared to the SA-CAS-srDFT method. The accuracy of the various CAS-srDFT methods was investigated for 139 singlet excitation energies of 28 typical organic chromophores. The two excited-state approaches in conjunction with multiconfiguration pair-density functional theory (MC-PDFT) were also employed in the benchmark study for comparing the accuracy with CAS-srDFT. It was found that CI-srDFT methods are more accurate than their SA counterparts, giving a mean absolute error of just 0.17 eV when using the sr-ctPBE functional. The accuracy of the new SA-based CAS-srDFT methods was observed to be impressive for organic molecules; however, this was not found to be transferable when investigating excited states of transition-metal complexes. In fact, none of the CASSCF-DFT excited-state methods introduced in this study, and also MC-PDFT, were found to provide a consistent improvement of CASSCF excitation energies.

The accurate description of reaction barrier heights is challenging for quantum mechanical methods due to the need for a balanced treatment of dynamic and static correlation energies because their importance varies during the course of a chemical reaction. While some regions of potential energy surfaces are well-described by a single-reference wave function or by Kohn-Sham density functional theory, in other cases a multireference treatment is needed. For systems with many active electrons, most accurate multireference methods have prohibitive computational scalings with system size. Multiconfiguration pair-density functional theory, MC-PDFT, is a more affordable multireference approach that computes the total electron correlation energy in a single step by using the multiconfiguration kinetic energy, density, and on-top pair density and an on-top density functional. In this work, we apply MC-PDFT to a benchmark database (DBH24/18) of 24 diverse reaction barrier heights. We explore the role of active space and basis set selection on the performance of MC-PDFT. We find that MC-PDFT is able to calculate reaction barrier heights with a similar accuracy to complete active space second order perturbation theory, CASPT2, but at a lower computational cost, and we find that MC-PDFT is less dependent on basis set selection than CASPT2.

The accurate description of ground- and excited-state potential energy surfaces poses a challenge for many electronic structure methods, especially in regions where strong electronic state interaction occurs. Here we introduce a new methodology, state-interaction pair-density functional theory (SI-PDFT), to target molecular systems exhibiting strong interaction of electronic states. SI-PDFT is an extension of multiconfiguration pair-density functional theory in which a set of N electronic states is generated through the diagonalization of an N × N effective Hamiltonian. We demonstrate the accuracy of the method by performing calculations on the ionic-neutral avoided crossing in lithium fluoride and the 1ππ-1πσ* avoided crossing in the H-O bond photodissociation in phenol. We show that SI-PDFT can be a useful tool in the study of photochemistry and nonadiabatic dynamics.

Quantum Chemical Study of Coupled Molecular Photoswitches

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.

In this article, the theoretical photoionization spectroscopy of the aziridine (C2H5N) molecule is investigated. To start with, we have optimized the geometry of this molecule at the neutral electronic ground state at the density functional theory/augmented correlation-consistent polarized valence triple zeta level of theory using the G09 program. The electronic structure calculations were restricted to the first six low-lying electronic states in order to account for the experimental photoelectron spectrum of the C2H5N molecule. The first six low-lying electronic states (X̃2A', Ã2A', B̃2A″, C̃2A″, D̃2A', and Ẽ2A') of the potential energy surfaces (PESs) are calculated by both equation of motion-ionization potential-coupled cluster singles and doubles and multi-configuration quasi-degenerate perturbation theory abinitio quantum chemistry methods along the dimensionless normal displacement coordinates in which multiple conical intersections were established among the considered electronic states. A (6 × 6) model vibronic Hamiltonian is constructed on a diabatic electronic basis, using the symmetry selection rules and Taylor series expansion. The Cs symmetry point group of the aziridine molecule leads to electronic states symmetry of either A' or A″, and these states are close in energy, due to which the same symmetry electronic states avoid each other. To get a smooth diabatic PES, a fourfold diabatization scheme is used, which is implemented in the General Atomic and Molecular Electronic Structure Systems suite of programs. All the parameters used in the diabatic vibronic coupling model Hamiltonian are calculated in terms of the normal modes of vibrational coordinates. Finally, the vibronic model Hamiltonian constructed for the coupled six electronic states is used to solve both time-independent and time-dependent Schrödinger equations using the multi-configuration time-dependent Hartree program module to obtain the dynamical observables. The theoretical vibronic band structure is found to be in good accord with the available experimental results.

Accurate rate constants are crucial for understanding and optimizing catalytic reactions mediated by enzymes, metalloproteins, and heterogeneous catalysts. These systems frequently present a dual computational challenge. Multiconfigurational reaction sites require multireference techniques for the accurate treatment of the electronic structure, and high activation barriers prevent efficient sampling of unbiased reactive transitions. In this work, we combine multiconfiguration pair-density functional theory (MC-PDFT) as an accurate and efficient multireference electronic structure method with on-the-fly probability-enhanced sampling flooding (OPESf) as an enhanced sampling method capable of accelerating reactive transitions. We demonstrate the approach on the Diels-Alder [4+2] cycloaddition between cis-butadiene and ethene as a reaction characterized by a large activation barrier and multireference character. MC-PDFT-OPESf provides reaction rates in agreement with experiments at a fraction of the computational cost required by conventional unbiased ab initio calculations. We propose MC-PDFT-OPESf as an efficient approach for computing kinetics in strongly correlated molecular systems.