MC-PDFT can calculate singlet–triplet splittings of organic diradicals

Multireference Correlated Oscillator Strengths from Adiabatic Connection Approaches Based on Extended Random Phase Approximation

We show that accurate oscillator strengths can be obtained from adiabatic connection (AC) approaches based on the extended random phase approximation (ERPA) combined with multireference (complete active space, CAS) wave functions. The oscillator strengths calculated using the perturbation-corrected ERPA transition density matrices, proposed in this work, and the excitation energies calculated with recently introduced AC correlation energy methods, AC0 and AC0D, compete with accuracy in the perturbational CASPT2 approach and require less computational effort. AC0 and AC0D methods scale more favorably with the number of active orbitals than multiconfigurational perturbation approaches like CASPT2 and NEVPT2 thanks to their dependence on reduced density matrices up to the order of 2. Importantly, the newly developed approach for computing correlated transition dipole moments does not entail any additional costs, as all intermediate quantities become available when AC0 energies are being computed. We also test the performance of the recently proposed AC method corrected for the negative-transition contributions to the correlation energy, AC0D, for triplet excitation energies. Similarly, as for the singlet excitations, the correction improves the performance of the AC0 method, particularly for the low-lying excited states.

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.

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.

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.

Multiconfiguration Pair-Density Functional Theory for Chromium(IV) Molecular Qubits

Pseudotetrahedral organometallic complexes containing chromium(IV) and aryl ligands have been experimentally identified as promising molecular qubit candidates. Here we present a computational protocol based on multiconfiguration pair-density functional theory for computing singlet-triplet gaps and zero-field splitting (ZFS) parameters in Cr(IV) aryl complexes. We find that two multireference methods, multistate complete active space second-order perturbation theory (MS-CASPT2) and hybrid multistate pair-density functional theory (HMS-PDFT), perform better than Kohn-Sham density functional theory for singlet-triplet gaps. Despite the very small values of the ZFS parameters, both multireference methods performed qualitatively well. MS-CASPT2 and HMS-PDFT performed particularly well for predicting the trend in the ratio of the rhombic and axial ZFS parameters, |E/D|. We have also investigated the dependence and sensitivity of the calculated ZFS parameters on the active space and the molecular geometry. The methodologies outlined here can guide future prediction of ZFS parameters in molecular qubit candidates.

Difficulty of the evaluation of the barrier height of an open-shell transition state between closed shell minima: The case of small C4n rings

C_4n cyclacenes exhibit strong bond-alternation in their equilibrium geometry. In the two equivalent geometries, the system keeps an essentially closed-shell character. The two energy minima are separated by a transition state suppressing the bond-alternation, where the wave function is strongly diradical. This paper discusses the physical factors involved in this energy difference and possible evaluations of the barrier height. The barrier given as the energy difference between the restricted density functional theory (DFT)/B3LYP for the equilibrium and the broken symmetry DFT/B3LYP of the transition state is either negative or small, in contradiction with the most reliable Wave Function Theory calculations. The minimal (two electrons in two molecular orbitals) Complete Active Space self-consistent field (CASSCF) overestimates the barrier, and the subsequent second-order perturbation cancels it. Due to the collective character of the spin-polarization effect, it is necessary to perform a full π CASSCF + second-order perturbation to reach a reasonable value of the barrier, but this type of treatment cannot be applied to large molecules. DFT procedures treating on an equal foot the closed-shell and open-shell geometries have been explored, such as Mixed-Reference Spin-Flip Time-dependent-DFT and a new spin-decontamination proposal, namely, DFT-dressed configuration interaction, but the results still depend on the density functional. M06-2X without or with spin-decontamination gives the best agreement with the accurate wave function results.

Efficient Adiabatic Connection Approach for Strongly Correlated Systems: Application to Singlet-Triplet Gaps of Biradicals

Strong electron correlation can be captured with multireference wave function methods, but an accurate description of the electronic structure requires accounting for the dynamic correlation, which they miss. In this work, a new approach for the correlation energy based on the adiabatic connection (AC) is proposed. The AC_n method accounts for terms up to order n in the coupling constant, and it is size-consistent and free from instabilities. It employs the multireference random phase approximation and the Cholesky decomposition technique, leading to a computational cost growing with the fifth power of the system size. Because of the dependence on only one- and two-electron reduced density matrices, AC_n is more efficient than existing ab initio multireference dynamic correlation methods. AC_n affords excellent results for singlet-triplet gaps of challenging organic biradicals. The development presented in this work opens new perspectives for accurate calculations of systems with dozens of strongly correlated electrons.

Diffusion quantum Monte Carlo method on diradicals using single- and multi-determinant-Jastrow trial wavefunctions and different orbitals

In this work, the diffusion quantum Monte Carlo (DMC) method is employed to calculate the energies of singlet and triplet states for a series of organic diradicals and diatomic diradicals with π2 configuration. Single-determinant-Jastrow (SDJ) trial wavefunctions for triplet states, two-determinant-Jastrow (2DJ) trial wavefunctions for the singlet states, and multi-determinant-Jastrow (MDJ) trial wavefunctions are employed in DMC calculations using restricted open-shell B3LYP (ROB3LYP) orbitals, complete-active-space self-consistent field (CASSCF) orbitals, state-average CASSCF orbitals, or frozen-CASSCF orbitals. Our results show that DMC energies using either SDJ/2DJ or MDJ with ROB3LYP orbitals are close to or lower than those with the other orbitals for organic diradicals, while they are not very sensitive to the employed orbitals for diatomic diradicals. Furthermore, using MDJ can reduce DMC energies to some extent for most of the investigated organic diradicals and some diatomic diradicals. The importance of MDJ on DMC energies can be estimated based on the percentage of main determinants in the CASCI wavefunction. On the other hand, singlet–triplet gaps can be calculated reasonably with DMC using MDJ with a mean absolute error of less than 2 kcal/mol with all these orbitals. CASCI wavefunctions using density functional theory orbitals are preferred in constructing MDJ trial wavefunctions in practical DMC calculations since it is easier to obtain such wavefunctions than CASSCF methods.

Machine-Learned Energy Functionals for Multiconfigurational Wave Functions

We introduce multiconfiguration data-driven functional methods (MC-DDFMs), a group of methods which aim to correct the total or classical energy of a qualitatively accurate multiconfigurational wave function using a machine-learned functional of some featurization of the wave function such as its density, on-top density, or both. On a data set of carbene singlet-triplet energy splittings, we show that MC-DDFMs are able to achieve near-benchmark performance on systems not used for training with a robust degree of active-space independence. Beyond demonstrating that the density and on-top density hold the information necessary to correct the singlet-triplet energy splittings of multiconfigurational wave functions, this approach shows great promise for the development of functionals for MC-PDFT because corrections to the classical energy appear to be more transferable to types of molecules not included in the training data than corrections to total energies such as those yielded by CASSCF or NEVPT2.

Role of Triplet States in the Photodynamics of Aniline

The dynamics of excited heteroaromatic molecules is a key to understanding the photoprotective properties of many biologically relevant chromophores that dissipate their excitation energy nonreactively and thereby prevent the detrimental effects of ultraviolet radiation. Despite their structural variability, most substituted aromatic compounds share a common feature of a repulsive ¹πσ* potential energy surface. This surface can lead to photoproducts, and it can also facilitate the population transfer back to the ground electronic state by means of a ¹πσ*/S₀ conical intersection. Here, we explore a hidden relaxation route involving the triplet electronic state of aniline, which has recently been discovered by means of time-selected photofragment translational spectroscopy [J. Chem. Phys. 2019, 151, 141101]. By using the recently available analytical gradients for multiconfiguration pair-density functional theory, it is now possible to locate the minimum-energy crossing points between states of different spin and therefore compute the intersystem crossing rates with a multireference method, rather than with the less reliable single-reference methods. Using such calculations, we demonstrate that the population loss of aniline in the T₁(³ππ*) state is dominated by C₆H₅NH₂ → C₆H₅NH· + H· dissociation, and we explain the long nonradiative lifetimes of the T₁(³ππ*) state at the excitation wavelengths of 294-264 nm.

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.

Spin-Flip Pair-Density Functional Theory: A Practical Approach To Treat Static and Dynamical Correlations in Large Molecules

We present a practical approach to treat static and dynamical correlation accurately in large multiconfigurational systems. The static correlation is taken into account by using the spin-flip approach, which is well-known for capturing static correlation accurately at low-computational expense. Unlike previous approaches to add dynamical correlation to spin-flip models which use perturbation theory or coupled-cluster theory, we explore the ability to use the on-top pair-density functional theory approaches recently developed by Gagliardi and co-workers (J. Comput. Theor. Chem., 2014, 10, 3669). External relaxations are performed in the spin-flip calculations through a restricted active space framework for which a truncation scheme for the orbitals used in the external excitation is presented. The performance of the approach is demonstrated by computing energy gaps between ground and excited states for diradicals, triradicals, and linear polyacene chains ranging from naphthalene to dodecacene. Accurate results are obtained using the new approach for these challenging open-shell molecular systems.

Analytic gradients for multiconfiguration pair-density functional theory with density fitting: Development and application to geometry optimization in the ground and excited states

Density fitting reduces the computational cost of both energy and gradient calculations by avoiding the computation and manipulation of four-index electron repulsion integrals. With this algorithm, one can efficiently optimize the geometries of large systems with an accurate multireference treatment. Here, we present the derivation of multiconfiguration pair-density functional theory for energies and analytic gradients with density fitting. Six systems are studied, and the results are compared to those obtained with no approximation to the electron repulsion integrals and to the results obtained by complete active space second-order perturbation theory. With the new approach, there is an increase in the speed of computation with a negligible loss in accuracy. Smaller grid sizes have also been used to reduce the computational cost of multiconfiguration pair-density functional theory with little effect on the optimized geometries and gradient values.

Unconventional singlet fission materials

Singlet fission (SF) is a photophysical downconversion pathway, in which a singlet excitation transforms into two triplet excited states. As such, it constitutes an exciton multiplication generation process, which is currently at the focal point for future integration into solar energy conversion devices. Beyond this, various other exciting applications were proposed, including quantum cryptography or organic light emitting diodes. Also, the mechanistic understanding evolved rapidly during the last year. Unfortunately, the number of suitable SF-chromophores is still limited. This is per se problematic, considering the wide range of envisaged applicability. With that in mind, we emphasize uncommon SF-scaffolds and outline requirements as well as strategies to expand the chromophore pool of SF-materials.

Spin-state dependence of exchange-correlation holes

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

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.

Transition Metal Spin-State Energetics by MC-PDFT with High Local Exchange

The energetics of the spin states of transition metal complexes have been explored with a variety of electronic structure methods, but the calculations require a compromise between accuracy and affordability. In this work, the spin splittings of several iron complexes are studied with multiconfiguration pair-density functional theory (MC-PDFT). The results are compared to previously published results obtained by complete active space second-order perturbation theory (CASPT2) and CASPT2 with coupled-cluster semicore correlation (CASPT2/CC). In contrast to CASPT2's systematic overstabilization of high-spin states with respect to the CASPT2/CC reference, MC-PDFT with the tPBE on-top functional understabilizes high-spin states. This systematic understabilization is largely corrected by revising the exchange and correlation contributions to the on-top functional using the high local-exchange approximation (tPBE-HLE). Moreover, tPBE-HLE correctly predicts the spin of the ground state in most cases, while CASPT2 incorrectly predicts high-spin ground states in all cases. This is encouraging for practical work because tPBE and tPBE-HLE are faster than CASPT2 by a factor of 50 even in a moderately sized example.

Capturing the Dynamic Correlation for Arbitrary Spin-Symmetry CASSCF Reference with Adiabatic Connection Approaches: Insights into the Electronic Structure of the Tetramethyleneethane Diradical

The recently proposed approach to multireference dynamic correlation energy based on the adiabatic connection (AC) is extended to an arbitrary spin symmetry of the reference state. We show that both the spin-free AC approach and its computationally inexpensive approximation, AC0, when combined with a complete active space wave function, constitute viable alternatives to the perturbation-based and density-functional-based multiconfiguration methods. In particular, the AC0 approach, thanks to its favorable scaling with the system size and the size of the active space, allows for treating larger systems than its perturbation-based counterparts while maintaining comparable accuracy. We show the method's robustness on illustrative chemical systems, including the elusive tetramethyleneethane (TME) diradical, potential energy surfaces of which present a challenge to most computational approaches. For the latter system, AC0 outperforms other methods, staying in close agreement with the full configuration interaction quantum Monte Carlo benchmark. A careful analysis of the contributions to the correlation energy of TME's lowest singlet and triplet states reveals the subtle interplay of the dynamic and static correlation as the key to understanding the shape of the diradical's potential energy surfaces.

Structural diradical character

A reliable first-principles description of singlet diradical character is essential for predicting nonlinear optical and magnetic properties of molecules. As diradical and closed-shell electronic structures differ in their distribution of single, double, triple, and aromatic bonds, modeling electronic diradical character requires accurate bond-length patterns, in addition to accurate absolute bond lengths. We therefore introduce structural diradical character, which we suggest as an additional measure for comparing first-principles calculations with experimental data. We employ this measure to identify suitable exchange-correlation functionals for predicting the bond length patterns and electronic diradical character of a biscobaltocene with the potential for photoswitchable nonlinear optical activity. Out of four popular approximate exchange-correlation functionals with different exact-exchange admixtures (BP86, TPSS, B3LYP, TPSSh), the two hybrid functionals TPSSh and B3LYP perform best for diradical bond length patterns, with TPSSh being best for the organometallic validation systems and B3LYP for the organic ones (for which the D3 dispersion correction was included). Still, none of the functionals is suitable for correctly describing relative bond lengths across the range of molecules studied, so that none can be recommended for predictive studies of (potential) diradicals without reservation. © 2018 Wiley Periodicals, Inc.

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.