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

A Merger of the Spin-Flip ORMAS Approach and the MC-PDFT Method

The SF-ORMAS-PDFT (spin-flip occupation restricted multiple active space-pair density functional theory) approach combines the SF-ORMAS-CI method with the MC-PDFT method to treat both static and dynamic correlation in multiconfigurational systems. The static correlation description is generated via the spin-flip approach, which uses a high-spin single reference determinant to treat excited states with multiconfigurational characters. The on-top pair density functional theory uses a translation scheme applied to GGA density functionals. The SF-ORMAS-PDFT scheme has also been combined with virtual valence orbitals (VVO), a well-defined subspace of the virtual molecular orbitals, giving rise to significant speedups relative to the use of the full virtual space. The accuracy of the SF-ORMAS-PDFT method is tested by calculating 65 vertical excitation energies of 12 small- and medium-sized organic molecules. The SF-ORMAS-PDFT vertical excitation energies calculated with VVOs are comparable to those calculated with the full virtual space. The SF-ORMAS-PDFT/6-31G(d) level of theory predicts the rotational barrier of ethylene to be 65.5 and 65.9 kcal/mol, with full virtual space and VVOs, respectively. These predicted barrier heights compare well with the experimental value of 65 kcal/mol.

Accurate Electronic and Optical Properties of Organic Doublet Radicals Using Machine Learned Range-Separated Functionals

Luminescent organic semiconducting doublet-spin radicals are unique and emergent optical materials because their fluorescent quantum yields (Φfl) are not compromised by the spin-flipping intersystem crossing (ISC) into a dark high-spin state. The multiconfigurational nature of these radicals challenges their electronic structure calculations in the framework of single-reference density functional theory (DFT) and introduces room for method improvement. In the present study, we extended our earlier development of ML-ωPBE [J. Phys. Chem. Lett., 2021, 12, 9516-9524], a range-separated hybrid (RSH) exchange-correlation (XC) functional constructed using the stacked ensemble machine learning (SEML) algorithm, from closed-shell organic semiconducting molecules to doublet-spin organic semiconducting radicals. We assessed its performance for a new test set of 64 doublet-spin radicals from five categories while placing all previously compiled 3926 closed-shell molecules in the new training set. Interestingly, ML-ωPBE agrees with the nonempirical OT-ωPBE functional regarding the prediction of the molecule-dependent range-separation parameter (ω), with a small mean absolute error (MAE) of 0.0197 a0-1, but saves the computational cost by 2.46 orders of magnitude. This result demonstrates an outstanding domain adaptation capacity of ML-ωPBE for diverse organic semiconducting species. To further assess the predictive power of ML-ωPBE in experimental observables, we also applied it to evaluate absorption and fluorescence energies (Eabs and Efl) using linear-response time-dependent DFT (TDDFT), and we compared its behavior with nine popular XC functionals. For most radicals, ML-ωPBE reproduces experimental measurements of Eabs and Efl with small MAEs of 0.299 and 0.254 eV, only marginally different from those of OT-ωPBE. Our work illustrates a successful extension of the SEML framework from closed-shell molecules to doublet-spin radicals and will open the venue for calculating optical properties for organic semiconductors using single-reference TDDFT.

Chemically accurate singlet-triplet gaps of organic chromophores and linear acenes by the random phase approximation and σ-functionals

Predicting the energy differences between different spin-states is challenging for many widely used ab initio electronic structure methods. We here assess the ability of the direct random phase approximation (dRPA), dRPA plus two different screened second-order exchange (SOX) corrections, and σ-functionals to predict adiabatic singlet-triplet gaps. With mean absolute deviations of below 0.1 eV to experimental reference values, independent of the Kohn-Sham starting point, dRPA and σ-functionals accurately predict singlet-triplet gaps of 18 organic chromophores. The addition of SOX corrections to dRPA considerably worsens agreement with experiment, adding to the mounting evidence that dRPA+SOX methods are not generally applicable beyond-RPA methods. Also for a series of linear acene chains with up to ten fused rings, dRPA, and σ-functionals are in excellent agreement with coupled-cluster single double triple reference data. In agreement with advanced multi-reference methods, dRPA@PBE and σ-functional@PBE predict a singlet ground state for all chain lengths, while dRPA@PBE0 and σ-functional@PBE0 predict a triplet ground state for longer acenes. Our work shows dRPA and σ-functionals to be reliable methods for calculating singlet-triplet gaps in aromatic molecules.

Efficient Implementation of Block-Correlated Coupled Cluster Theory Based on the Generalized Valence Bond Reference for Strongly Correlated Systems

An optimized implementation of block-correlated coupled cluster theory based on the generalized valence bond wave function (GVB-BCCC) for the singlet ground state of strongly correlated systems is presented. The GVB-BCCC method with two-pair correlation (GVB-BCCC2b) or up to three-pair correlation (GVB-BCCC3b) will be the focus of this work. Three major techniques have been adopted to dramatically accelerate GVB-BCCC2b and GVB-BCCC3b calculations. First, the GVB-BCCC2b and GVB-BCCC3b codes are noticeably optimized by removing redundant calculations. Second, independent amplitudes are identified by constraining excited configurations to be pure singlet states and only independent amplitudes need to be solved. Third, an incremental updating scheme for the amplitudes in solving the GVB-BCCC equations is adopted. With these techniques, accurate GVB-BCCC3b calculations are now accessible for systems with relatively large active spaces (50 electrons in 50 orbitals) and GVB-BCCC2b calculations are affordable for systems with much larger active spaces. We have applied GVB-BCCC methods to investigate three typical kinds of systems: polyacenes, pentaprismane, and [Cu₂O₂]²⁺ isomers. For polyacenes, we demonstrate that GVB-BCCC3b can capture more than 94% of the total correlation energy even for 12-acene with 50 π electrons. For the potential energy curve of simultaneously stretching 15 C-C bonds in pentaprismane, our calculations show that the GVB-BCCC3b results are quite close to the results from the density matrix renormalization group (DMRG) over the whole range. For two dinuclear copper oxide isomers, their relative energy predicted by GVB-BCCC3b is also in good accord with the DMRG result. All calculations show that the inclusion of three-pair correlation in GVB-BCCC is critical for accurate descriptions of strongly correlated systems.

Perspective on Simplified Quantum Chemistry Methods for Excited States and Response Properties

We review recent developments in the framework of simplified quantum chemistry for excited state and optical response properties (sTD-DFT) and present future challenges for new method developments to improve accuracy and extend the range of application. In recent years, the scope of sTD-DFT was extended to molecular response calculations of the polarizability, optical rotation, first hyperpolarizability, two-photon absorption (2PA), and excited-state absorption for large systems with hundreds to thousands of atoms. The recently introduced spin-flip simplified time-dependent density functional theory (SF-sTD-DFT) variant enables an ultrafast treatment for diradicals and related strongly correlated systems. A few drawbacks were also identified, specifically for the computation of 2PA cross sections. We propose solutions to this problem and how to generally improve the accuracy of simplified schemes. New possible simplified schemes are also introduced for strongly correlated systems, e.g., with a second-order perturbative correlation correction. Interpretation tools that can extract chemical structure-property relationships from excited state or response calculations are also discussed. In particular, the recently introduced method-agnostic RespA approach based on natural response orbitals (NROs) as the key concept is employed.