Restricted Open-shell Kohn-Sham Research Articles

Benchmarking Charge-Transfer Excited States in TADF Emitters: ΔDFT Outperforms TD-DFT for Emission Energies.

Charge-transfer (CT) excited states are crucial to organic light-emitting diodes (OLEDs), particularly to those based on thermally activated delayed fluorescence (TADF). However, accurately modeling CT states remains challenging, even with modern implementations of (time-dependent) density functional theory [(TD-)DFT], especially in a dielectric environment. To identify shortcomings and improve the methodology, we previously established the STGABS27 benchmark set with highly accurate experimental references for the adiabatic energy gap between the lowest singlet and triplet excited states (ΔEST). Here, we diversify this set to the STGABS27-EMS benchmark by including experimental emission energies (Eem) and use this new set to (re)-evaluate various DFT-based approaches. Surprisingly, these tests demonstrate that a state-specific (un)restricted open-shell Kohn-Sham (U/ROKS) DFT coupled with a polarizable continuum model for perturbative state-specific nonequilibrium solvation (ptSS-PCM) provides exceptional accuracy for predicting Eem over a wide range of density functionals. In contrast, the main workhorse of the field, Tamm-Dancoff-approximated TD-DFT (TDA-DFT) paired with the same ptSS-PCM, is distinctly less accurate and strongly functional-dependent. More importantly, while TDA-DFT requires the choice of two very different density functionals for good performance on either ΔEST or Eem, the time-independent U/ROKS/PCM approaches deliver excellent accuracy for both quantities with a wide variety of functionals.

Predicting the X-ray Absorption Spectrum of Ozone with Single Configuration State Functions.

X-ray absorption spectra (XAS) of biradicaloid species are often thought to represent a challenge to theoretical methods. This has led to the testing of recently developed multireference techniques on the XAS of ozone, but reproduction of the experimental spectral profile has proven difficult. We utilize a minimal model consisting of a single configuration state function (CSF) per excited state to model core-level excitations of ozone, with the orbitals of each CSF optimized using the restricted open-shell Kohn-Sham (ROKS) method. This protocol leads to semiquantitative agreement with experimental XAS. In fact, we find that low-lying core-hole excited states in biradicaloids can be approximated with individual CSFs, despite the presence of multireference character in the ground state. We also report that the 1s → π* and 1s → σ* transitions have quite distinct widths for O3. This reveals the importance of sampling over a representative range of geometries from the vibrational ground state for properly assessing the accuracy of electronic structure methods against experiments instead of the popular procedure of uniformly broadening stick spectra at the equilibrium geometry.

Exciton Nature of Plasma Phase Transition in Warm Dense Fluid Hydrogen: ROKS Simulation.

The transition of warm dense fluid hydrogen from an insulator to a conducting state at pressures of about 20-400 GPa and temperatures of 500-5000 K has been the subject of active scientific research over the past few decades. However, various experimental and theoretical methods do not provide consistent results. In this work, we have applied the restricted open-shell Kohn-Sham (ROKS) method for first principles molecular dynamics of dense hydrogen after thermal excitation to the first singlet excited state. The Wannier localization method has allowed us to analyze the exciton dynamics in this system. The model shows that a key mechanism of the transition is associated with the dissociation of electron-hole pairs, which allows explaining several stages of the transition of fluid H2 from molecular state to plasma. This mechanism is able to give a quantitative description of several experimental results as well as to resolve the discrepancies between experimental studies.

Computing x-ray absorption spectra from linear-response particles atop optimized holes.

State specific orbital optimized density functional theory (OO-DFT) methods, such as restricted open-shell Kohn-Sham (ROKS), can attain semiquantitative accuracy for predicting x-ray absorption spectra of closed-shell molecules. OO-DFT methods, however, require that each state be individually optimized. In this Communication, we present an approach to generate an approximate core-excited state density for use with the ROKS energy ansatz, which is capable of giving reasonable accuracy without requiring state-specific optimization. This is achieved by fully optimizing the core-hole through the core-ionized state, followed by the use of electron-addition configuration interaction singles to obtain the particle level. This hybrid approach can be viewed as a DFT generalization of the static-exchange (STEX) method and can attain ∼0.6eV rms error for the K-edges of C-F through the use of local functionals, such as PBE and OLYP. This ROKS(STEX) approach can also be used to identify important transitions for full OO ROKS treatment and can thus help reduce the computational cost of obtaining OO-DFT quality spectra. ROKS(STEX), therefore, appears to be a useful technique for the efficient prediction of x-ray absorption spectra.

A theoretical perspective of the ultrafast transient absorption dynamics of CsPbBr3.

Transient absorption spectra (TAS) of lead halide perovskites can provide important insights into the nature of the photoexcited state dynamics of this prototypical class of materials. Here, we perform ground and excited state molecular dynamics (MD) simulations within a restricted open shell Kohn‐Sham (ROKS) approach in order to interpret the characteristic features of the TAS of CsPbBr3. Our results reveal that properties such as the finite temperature band gap, the Stokes shift, and therefore, also the TAS are strongly size‐dependent. Our TAS simulations show an early positive red‐shifted feature on the fs scale that can be explained by geometric relaxation in the excited state. As excited‐state processes can crucially affect the electronic properties of this class of photoactive materials, our observations are an important ingredient for further optimization of lead halide based optoelectronic devices.

PCM-ROKS for the Description of Charge-Transfer States in Solution: Singlet-Triplet Gaps with Chemical Accuracy from Open-Shell Kohn-Sham Reaction-Field Calculations.

The adiabatic energy gap between the lowest singlet and triplet excited states ΔEST is a central property of thermally activated delayed fluorescence (TADF) emitters. Since these states are dominated by a charge-transfer character, causing strong orbital-relaxation and environmental effects, an accurate prediction of ΔEST is very challenging, even with modern quantum-chemical excited-state methods. Addressing this major challenge, we present an approach that combines spin-unrestricted (UKS) and restricted open-shell Kohn-Sham (ROKS) self-consistent field calculations with a polarizable-continuum model and range-separated hybrid functionals. Tests on a new representative benchmark set of 27 TADF emitters with accurately known ΔEST values termed STGABS27 reveal a robust and unprecedented performance with a mean absolute deviation of only 0.025 eV (∼0.5 kcal/mol) and few deviations greater than 0.05 eV (∼1 kcal/mol), even in electronically challenging cases. Requiring only two geometry optimizations per molecule at the ROKS/UKS level in a compact double-ζ basis, the approach is computationally efficient and can routinely be applied to molecules with more than 100 atoms.

Analytic energy gradients of spin-adapted open-shell time-dependent density functional theory.

It is now well established that the spin-adapted time-dependent density functional theory [X-TD-DFT; Li and Liu, J. Chem. Phys. 135, 194106 (2011)] for low-lying excited states of open-shell systems has very much the same accuracy as the conventional TD-DFT for low-lying excited states of closed-shell systems. In particular, this has been achieved without computational overhead over the unrestricted TD-DFT (U-TD-DFT) that usually produces heavily spin-contaminated excited states. It is shown here that the analytic energy gradients of X-TD-DFT can be obtained by just slight modifications of those of U-TD-DFT running with restricted open-shell Kohn-Sham orbitals. As such, X-TD-DFT also has no overhead over U-TD-DFT in the calculation of energy gradients of excited states of open-shell systems. Although only a few prototypical open-shell molecules are considered as showcases, it can definitely be said that X-TD-DFT can replace U-TD-DFT for geometry optimization and dynamics simulation of excited states of open-shell systems.

Excited State Orbital Optimization via Minimizing the Square of the Gradient: General Approach and Application to Singly and Doubly Excited States via Density Functional Theory.

We present a general approach to converge excited state solutions to any quantum chemistry orbital optimization process, without the risk of variational collapse. The resulting square gradient minimization (SGM) approach only requires analytic energy/Lagrangian orbital gradients and merely costs 3 times as much as ground state orbital optimization (per iteration), when implemented via a finite difference approach. SGM is applied to both single determinant ΔSCF and spin-purified restricted open-shell Kohn-Sham (ROKS) approaches to study the accuracy of orbital optimized DFT excited states. It is found that SGM can converge challenging states where the maximum overlap method (MOM) or analogues either collapse to the ground state or fail to converge. We also report that ΔSCF/ROKS predict highly accurate excitation energies for doubly excited states (which are inaccessible via TDDFT). Singly excited states obtained via ROKS are also found to be quite accurate, especially for Rydberg states that frustrate (semi)local TDDFT. Our results suggest that orbital optimized excited state DFT methods can be used to push past the limitations of TDDFT to doubly excited, charge-transfer, or Rydberg states, making them a useful tool for the practical quantum chemist's toolbox for studying excited states in large systems.

Highly Accurate Prediction of Core Spectra of Molecules at Density Functional Theory Cost: Attaining Sub-electronvolt Error from a Restricted Open-Shell Kohn-Sham Approach.

We present the use of the recently developed square gradient minimization (SGM) algorithm for excited-state orbital optimization to obtain spin-pure restricted open-shell Kohn-Sham (ROKS) energies for core excited states of molecules. The SGM algorithm is robust against variational collapse and offers a reliable route to converging orbitals for target excited states at only 2-3 times the cost of ground-state orbital optimization (per iteration). ROKS/SGM with the modern SCAN/ωB97X-V functionals is found to predict the K-edge of C, N, O, and F to a root mean squared error of ∼0.3 eV. ROKS/SGM is equally effective at predicting L-edge spectra of third period elements, provided a perturbative spin-orbit correction is employed. This high accuracy can be contrasted with traditional time-dependent density functional theory (TDDFT), which typically has greater than 10 eV error and requires translation of computed spectra to align with experiment. ROKS is computationally affordable (having the same scaling as ground-state DFT and a slightly larger prefactor) and can be applied to geometry optimizations/ab initio molecular dynamics of core excited states, as well as condensed phase simulations. ROKS can also model doubly excited/ionized states with one broken electron pair, which are beyond the ability of linear response based methods.

How Parallel Are Excited State Potential Energy Surfaces from Time-Independent and Time-Dependent DFT? A BODIPY Dye Case Study

To support the development and characterization of chromophores with targeted photophysical properties, excited-state electronic structure calculations should rapidly and accurately predict how derivatization of a chromophore will affect its excitation and emission energies. This paper examines whether a time-independent excited-state density functional theory (DFT) approach meets this need through a case study of BODIPY chromophore photophysics. A restricted open-shell Kohn-Sham (ROKS) treatment of the S1 excited state of BODIPY dyes is contrasted with linear-response time-dependent density functional theory (TDDFT). Vertical excitation energies predicted by the two approaches are remarkably different due to overestimation by TDDFT and underestimation by ROKS relative to experiment. Overall, ROKS with a standard hybrid functional provides the more accurate description of the S1 excited state of BODIPY dyes, but excitation energies computed by the two methods are strongly correlated. The two approaches also make similar predictions of shifts in the excitation energy upon functionalization of the chromophore. TDDFT and ROKS models of the S1 potential energy surface are then examined in detail for a representative BODIPY dye through molecular dynamics sampling on both model surfaces. We identify the most significant differences in the sampled surfaces and analyze these differences along selected normal modes. Differences between ROKS and TDDFT descriptions of the S1 potential energy surface for this BODIPY derivative highlight the continuing need for validation of widely used approximations in excited state DFT through experimental benchmarking and comparison to ab initio reference data.

Prediction of Excited-State Energies and Singlet-Triplet Gaps of Charge-Transfer States Using a Restricted Open-Shell Kohn-Sham Approach.

Organic molecules with charge-transfer (CT) excited states are widely used in industry and are especially attractive as candidates for fabrication of energy efficient OLEDs, as they can harvest energy from nonradiative triplets by means of thermally activated delayed fluorescence (TADF). It is therefore useful to have computational protocols for accurate estimation of their electronic spectra in order to screen candidate molecules for OLED applications. However, it is difficult to predict the photophysical properties of TADF molecules with LR-TDDFT, as semilocal LR-TDDFT is incapable of accurately modeling CT states. Herein, we study absorption energies, emission energies, zero-zero transition energies, and singlet-triplet gaps of TADF molecules using a restricted open-shell Kohn-Sham (ROKS) approach instead and discover that ROKS calculations with semilocal hybrid functionals are in good agreement with experiments-unlike TDDFT, which significantly underestimates energy gaps. We also propose a cheap computational protocol for studying excited states with large CT character that is found to give good agreement with experimental results without having to perform any excited-state geometry optimizations.

Excitation energies and Stokes shifts from a restricted open-shell Kohn-Sham approach

Restricted open-shell Kohn-Sham (ROKS) theory provides a powerful computational tool for calculating singlet excited state energies and dynamics. However, the possibility of multiple solutions to the ROKS equations - with the associated difficulty of automatically selecting the physically meaningful solution - limits its usefulness for intensive applications such as long-time Born-Oppenheimer molecular dynamics. We present an implementation of ROKS for excited states which prescribes the physically correct solution from an overlap criterion and guarantees that this solution is stationary, allowing for straightforward evaluation of nuclear gradients. The method is used to benchmark ROKS for vertical excitation energies of small and large organic dyes and for the calculation of Stokes shifts. With common density functional approximations, ROKS vertical excitation energies, and Stokes shifts show similar accuracy to those from time-dependent density functional theory and Δ-self-consistent-field approaches. Advantages of the ROKS approach for excited state structure and molecular dynamics are discussed.

Systematic Theoretical Study of the Zero-Field Splitting in Coordination Complexes of Mn(III). Density Functional Theory versus Multireference Wave Function Approaches

This paper presents a detailed evaluation of the performance of density functional theory (DFT) as well as complete active space self-consistent field (CASSCF)-based methods (CASSCF and second-order N-electron valence state perturbation theory, NEVPT2) to predict the zero-field splitting (zfs) parameters for a series of coordination complexes containing the Mn(III) ion. The physical origin of the experimentally determined zfs's was investigated by studying the different contributions to these parameters. To this end, a series of mononuclear Mn(III) complexes was chosen for which the structures have been resolved by X-ray diffraction and the zfs parameters have been accurately determined by high-field EPR spectroscopy. In a second step, small models have been constructed to allow for a systematic assessment of the factors that dominate the variations in the observed zfs parameters and to establish magnetostructural correlations. Among the tested functionals, the best predictions have been obtained with B3LYP, followed by the nonhybrid BP86 functional, which in turn is more successful than the meta-hybrid GGA functional TPSSh. For the estimation of the spin-orbit coupling (SOC) part of the zfs, it was found that the coupled perturbed SOC approach CP is more successful than the Pederson-Khanna method. Concerning the spin-spin interaction (SS), the restricted open-shell Kohn-Sham (ROKS) approach led to a slightly better agreement with the experiment than the unrestricted KS (UKS) approach. The ab initio state-averaged CASSCF (SA-CASSCF) method with a minimal active space and the most recent implementation that treats the SOC and SS contributions on an equal footing provides the best predictions for the zfs. The analysis demonstrates that the major contribution to the axial zfs parameter (D) originates from the SOC interaction but that the SS part is far from being negligible (between 10 and 20% of D). Importantly, the various excited triplet ligand field states account for roughly half of the value of D, contrary to popular ligand field models. Despite covering dynamic correlation contributions to the transition energies, NEVPT2 does not lead to large improvements in the results as the excitation energies of the Mn(III) d-d transitions are already fairly accurate at the SA-CASSCF level. For a given type of coordination sphere (e.g., elongated or compressed octahedron), the magnetic anisotropy of the Mn(III) ion, D, does not appear to be highly sensitive to the nature of the ligands, while the E/D ratio is notably affected by all octahedral distortions. Furthermore, the introduction of different halides into the coordination sphere of Mn(III) only leads to small effects on D. Nevertheless, it appears that oxygen-based ligands afford larger D values than nitrogen-based ligands.

Spin-adapted open-shell random phase approximation and time-dependent density functional theory. I. Theory

The spin-adaptation of single-reference quantum chemical methods for excited states of open-shell systems has been nontrivial. The primary reason is that the configuration space, generated by a truncated rank of excitations from only one component of a reference multiplet, is spin-incomplete. Those "missing" configurations are of higher ranks and can, in principle, be recaptured by a particular class of excitation operators. However, the resulting formalisms are then quite involved and there are situations [e.g., time-dependent density functional theory (TD-DFT) under the adiabatic approximation] that prevent one from doing so. To solve this issue, we propose here a tensor-coupling scheme that invokes all the components of a reference multiplet (i.e., a tensor reference) rather than increases the excitation ranks. A minimal spin-adapted n-tuply excited configuration space can readily be constructed by tensor products between the n-tuple tensor excitation operators and the chosen tensor reference. Further combined with the tensor equation-of-motion formalism, very compact expressions for excitation energies can be obtained. As a first application of this general idea, a spin-adapted open-shell random phase approximation is first developed. The so-called "translation rule" is then adopted to formulate a spin-adapted, restricted open-shell Kohn-Sham (ROKS)-based TD-DFT (ROKS-TD-DFT). Here, a particular symmetry structure has to be imposed on the exchange-correlation kernel. While the standard ROKS-TD-DFT can access only excited states due to singlet-coupled single excitations, i.e., only some of the singly excited states of the same spin (S(i)) as the reference, the new scheme can capture all the excited states of spin S(i)-1, S(i), or S(i)+1 due to both singlet- and triplet-coupled single excitations. The actual implementation and computation are very much like the (spin-contaminated) unrestricted Kohn-Sham-based TD-DFT. It is also shown that spin-contaminated spin-flip configuration interaction approaches can easily be spin-adapted via the tensor-coupling scheme.

Restricted open-shell Kohn-Sham theory: Simulation of the pyrrole photodissociation

The authors study the photodissociation reactions of pyrrole and N-methylpyrrole using first-principles molecular dynamics. The first excited state is described with restricted open-shell Kohn-Sham theory. They find a small barrier in the excited state potential energy surface. The possibility of energy redistribution near the Franck-Condon region leads to two different reaction channels in on-the-fly simulations on a single diabatic potential energy surface. The results are discussed in comparison with previous ab initio calculations and with experiments.

Calculation of nonadiabatic couplings with restricted open-shell Kohn-Sham density-functional theory

This paper generalizes the recently proposed approaches for calculating the derivative couplings between adiabatic states in density-functional theory (DFT) based on a Slater transition-state density to transitions such as singlet-singlet excitations, where a single-determinant ansatz is insufficient. The proposed approach is based on restricted open-shell Frank et al. [J. Chem. Phys. 108, 4060 (1998)] theory used to describe a spin-adapted Slater transition state. To treat the dependence of electron-electron interactions on the nuclear positions, variational linear-response density-functional perturbation theory is generalized to reference states with an orbital-dependent Kohn-Sham Hamiltonian and nontrivial occupation patterns. The methods proposed in this paper are not limited to the calculation of derivative coupling vectors, but can also be used for the calculation of other transition matrix elements. Moreover, they can be used to calculate the linear response of open-shell systems to arbitrary external perturbations in DFT.

Spin−Spin Contributions to the Zero-Field Splitting Tensor in Organic Triplets, Carbenes and BiradicalsA Density Functional and Ab Initio Study

An evaluation study for the direct dipolar electron spin-spin (SS) contribution to the zero-field splitting (ZFS) tensor in electron paramagnetic resonance (EPR) spectroscopy is presented. Calculations were performed on a wide variety of organic systems where the SS contribution to the ZFS dominates over the second-order spin-orbit coupling (SOC) contribution. Calculations were performed using (hybrid) density functional theory (DFT), as well as complete active space self-consistent field (CASSCF) wave functions. In the former case, our implementation is an approximation, because we use the two-particle reduced spin-density matrix of the noninteracting reference system. In the latter case, the SS contribution is approximated by a mean-field method which, nevertheless, gives accurate results, compared to the approximation free computation of the SS part in a CASSCF framework. For the case of the triplet dioxygen molecule, it was shown that restricted open-shell density functional theory (RODFT), as well as CASSCF, can provide accurate spin-spin couplings while spin-unrestricted DFT leads to much larger errors. Furthermore, 15 organic radicals, including several 1,3 and 1,5 diradicals, dinitroxide biradicals, and even a chlorophyll a model system, were examined as test cases to demonstrate the accuracy and efficiency of our approach within a DFT framework. Accurate D values with root-mean-square deviations of 0.0035 cm(-1) were obtained. Furthermore, all trends, including those due to substituent effects, were correctly reproduced. In a different set of calculations, the polyacenes benzene, naphthalene, anthracene, and tetracene were studied. Applying DFT, the absolute D values were noticeably underestimated, but it was possible to correctly reproduce the trend to smaller D values with larger size of the systems. Finally, it was demonstrated that our approach is also well-suited for the study of carbenes. The smaller organic radicals of this work were also studied, through the use of CASSCF wave functions. This was a special advantage in the case of the triplet polyacenes, where the CASSCF approach gave better results than the DFT method. In comparing spin-restricted and spin-unrestricted results, it was shown through a natural orbital analysis and comparison to high-level ab initio calculations that even small amounts of spin polarization introduced by the unrestricted calculations lead to large deviations between the unrestricted Kohn-Sham (UKS) and restricted open-shell Kohn-Sham (ROKS) approaches. It is challenging to understand why the ROKS results show much better correlation with the experimental data.

First‐Principles Simulation of the Photoreaction of a Capped Azobenzene: The Rotational Pathway is Feasible

We present first-principles molecular dynamics simulations of azobenzene and a sterically hindered derivative in the first excited state. The restricted open-shell Kohn-Sham (ROKS) approach is employed to describe the motion in the lowest excited state. The rotational pathway is observed in the molecular dynamics simulations for both azobenzene and its azacrown ether capped derivative.

Molecular dynamics in electronically excited states using time-dependent density functional theory

We describe two different implementations of time-dependent density functional theory (TDDFT) for use in excited state molecular dynamics simulations. One is based on the linear response formulation (LR-TDDFT), whereas the other uses a time propagation scheme for the electronic wave functions (P-TDDFT). Photo-induced cis–trans isomerization of C=C, C=N and N=N double bonds is investigated in three model compounds, namely the 2,4-pentadiene-1-iminium cation (PSB), formaldimine and diimide. For formaldimine and diimide, the results obtained with both schemes are in agreement with experimental data and previously reported theoretical results. Molecular dynamics simulations yield new insights into the relaxation pathways in the excited state. For PSB, which is a model system for the retinal protonated Schiff base involved in the visual process, the forces computed from the LR-TDDFT S1 surface lead to an increased bond length alternation and, consequently, to single bond rotation. On the contrary, P-TDDFT dynamics lead to a decreased bond length alternation, in agreement with CASPT2 and restricted open-shell Kohn–Sham (ROKS) calculations.

Comment on “Excitations in photoactive molecules from quantum Monte Carlo” [J. Chem. Phys. 121, 5836 (2004)

It is shown that the qualitative differences between high-level ab initio calculations and restricted open-shell Kohn-Sham (ROKS) results for the lowest singlet excited electronic state of formaldimine along a particular isomerization path found by Schautz, Buda, and Filippi [J. Chem. Phys.121, 5836 (2004)] play a minor role in molecular dynamics simulations of photoisomerization at room temperature. In fact, ROKS yields, within its well-known limitations, a good representation of the physically relevant isomerization pathway.