SPOCK.CI: A multireference spin-orbit configuration interaction method for large molecules

Unified Implementation of Relativistic Wave Function Methods: 4C-iCIPT2 as a Showcase

In parallel to the unified construction of relativistic Hamiltonians based solely on physical arguments (J. Chem. Phys. 2024, 160, 084111), a unified implementation of relativistic wave function methods is achieved here via programming techniques (e.g., template metaprogramming and polymorphism in C++). That is, once the code for constructing the Hamiltonian matrix is made ready, all the rest can be generated automatically from existing templates used for the nonrelativistic counterparts. This is facilitated by decomposing a second-quantized relativistic Hamiltonian into diagrams that are topologically the same as those required for computing the basic coupling coefficients between spin-free configuration state functions (CSF). Moreover, both time reversal and binary double point group symmetries can readily be incorporated into molecular integrals and Hamiltonian matrix elements. The latter can first be evaluated in the space of (randomly selected) spin-dependent determinants and then transformed to that of spin-dependent CSFs, thanks to simple relations in between. As a showcase, we consider here the no-pair four-component relativistic iterative configuration interaction with selection and perturbation correction (4C-iCIPT2), which is a natural extension of the spin-free iCIPT2 (J. Chem. Theory Comput. 2021, 17, 949), and can provide near-exact numerical results within the manifold of positive energy states (PES), as demonstrated by numerical examples.

Toward the improvement of vibronic spectra and non-radiative rate constants using the vertical Hessian method

For the computation of vibrationally resolved electronic spectra, various approaches can be employed. Adiabatic approaches simulate vibronic transitions using harmonic potentials of the initial and final states, while vertical approaches extrapolate the final state potential from the gradients and Hessian at the Franck-Condon point, avoiding a full exploration of the potential energy surface of the final state. Our implementation of the vertical Hessian (VH) method has been validated with a benchmark set of four small molecules, each presenting unique challenges, such as complex topologies, problematic low-frequency vibrations, or significant geometrical changes upon electronic excitation. We assess the quality of both adiabatic and vertical approaches for simulating vibronic transitions. For two types of donor-acceptor compounds with promising thermally activated delayed fluorescence properties, our computations confirm that the vertical approaches outperform the adiabatic ones. The VH method significantly reduces computational costs and yields meaningful emission spectra, where adiabatic models fail. More importantly, we pioneer the use of the VH method for the computation of rate constants for non-radiative processes, such as intersystem crossing and reverse intersystem crossing along a relaxed interpolated pathway of a donor-acceptor compound. This study highlights the potential of the VH method to advance computational vibronic spectroscopy by providing meaningful simulations of intricate decay pathway mechanisms in complex molecular systems.

Consistent Second-Order Treatment of Spin-Orbit Coupling and Dynamic Correlation in Quasidegenerate N-Electron Valence Perturbation Theory

We present a formulation and implementation of second-order quasidegenerate N-electron valence perturbation theory (QDNEVPT2) that provides a balanced and accurate description of spin-orbit coupling and dynamic correlation effects in multiconfigurational electronic states. In our approach, the energies and wave functions of electronic states are computed by treating electron repulsion and spin-orbit coupling operators as equal perturbations to the nonrelativistic complete active-space wave functions, and their contributions are incorporated fully up to the second order. The spin-orbit effects are described using the Breit-Pauli (BP) or exact two-component Douglas-Kroll-Hess (DKH) Hamiltonians within spin-orbit mean-field approximation. The resulting second-order methods (BP2- and DKH2-QDNEVPT2) are capable of treating spin-orbit coupling effects in nearly degenerate electronic states by diagonalizing an effective Hamiltonian expanded in a compact non-relativistic basis. For a variety of atoms and small molecules across the entire periodic table, we demonstrate that DKH2-QDNEVPT2 is competitive in accuracy with variational two-component relativistic theories. BP2-QDNEVPT2 shows high accuracy for the second- and third-period elements, but its performance deteriorates for heavier atoms and molecules. We also consider the first-order spin-orbit QDNEVPT2 approximations (BP1- and DKH1-QDNEVPT2), among which DKH1-QDNEVPT2 is reliable but less accurate than DKH2-QDNEVPT2. Both DKH1- and DKH2-QDNEVPT2 hold promise as efficient and accurate electronic structure methods for treating electron correlation and spin-orbit coupling in a variety of applications.

Multiexcitonic and optically bright states in subunits of pentacene crystals: A hybrid DFT/MRCI and molecular mechanics study

A hybrid quantum mechanics/molecular mechanics setup was used to model electronically excited pentacene in the crystal phase. Particularly interesting in the context of singlet fission (SF) is the energetic location of the antiferromagnetically coupled multiexcitonic singlet state, 1(TT), and the ferromagnetically coupled analog in relation to the optically bright singlet state. To provide photophysical properties of the accessible spin manifold, combined density functional theory and multi-reference configuration interaction calculations were performed on pentacene dimers and a trimer, electrostatically embedded in the crystal. The likelihood of a quintet intermediate in the SF process was estimated by computing singlet-quintet electron spin-spin couplings employing the Breit-Pauli Hamiltonian. The performance of the applied methods was assessed on the pentacene monomer. The character of the optically bright state and the energetic location of the 1(TT) state depend strongly on the relative orientation of the pentacene units. In the V-shaped dimers and in the trimer, the optically bright state is dominated by local and charge transfer (CT) excitations, with admixtures of doubly excited configurations. The CT excitations gain weight upon geometry relaxation, thus supporting a CT-mediated SF mechanism as the primary step of the SF process. For the slip-stacked dimer, the energetic order of the bright and the 1(TT) states swaps upon geometry relaxation, indicating strong nonadiabatic coupling close to the Franck-Condon region-a prerequisite for a coherent SF process. The multiexcitonic singlet, triplet, and quintet states are energetically too far apart and their spin-spin couplings are too small to bring about a noteworthy multiplicity mixing.

Structural Control of Highly Efficient Thermally Activated Delayed Fluorescence in Carbene Zinc(II) Dithiolates

Luminescent metal complexes based on earth abundant elements are a valuable target to substitute 4d/5d transition metal complexes as triplet emitters in advanced photonic applications. Whereas CuI complexes have been thoroughly investigated in the last two decades for this purpose, no structure-property-relationships for efficient luminescence involving triplet excited states from ZnII complexes are established. Herein, we report on the design of monomeric carbene zinc(II) dithiolates (CZT) featuring a donor-acceptor-motif that leads to highly efficient thermally activated delayed fluorescence (TADF) with for ZnII compounds unprecedented radiative rate constants kTADF =1.2×106  s-1 at 297 K. Our high-level DFT/MRCI calculations revealed that the relative orientation of the ligands involved in the ligand-to-ligand charge transfer (1/3 LLCT) states is paramount to control the TADF process. Specifically, a dihedral angle of 36-40° leads to very efficient reverse intersystem-crossing (rISC) on the order of 109  s-1 due to spin-orbit coupling (SOC) mediated by the sulfur atoms in combination with a small ΔES1-T1 of ca. 56 meV.

Synthesis, Structural Characterization, and Phosphorescence Properties of Trigonal Zn(II) Carbene Complexes

The sterically demanding N-heterocyclic carbene ITr (N,N'-bis(triphenylmethyl)imidazolylidene) was employed for the preparation of novel trigonal zinc(II) complexes of the type [ZnX2(ITr)] [X = Cl (1), Br (2), and I (3)], for which the low coordination mode was confirmed in both solution and solid state. Because of the atypical coordination geometry, the reactivity of 1-3 was studied in detail using partial or exhaustive halide exchange and halide abstraction reactions to access [ZnLCl(ITr)] [L = carbazolate (4), 3,6-di-tert-butyl-carbazolate (5), phenoxazine (6), and phenothiazine (7)], [Zn(bdt)(ITr)] (bdt = benzene-1,2-dithiolate) (8), and cationic [Zn(μ2-X)(ITr)]2[B(C6F5)4]2 [X = Cl (9), Br (10), and I (11)], all of which were isolated and structurally characterized. Importantly, for all complexes 4-11, the trigonal coordination environment of the ZnII ion is maintained, demonstrating a highly stabilizing effect due to the steric demand of the ITr ligand, which protects the metal center from further ligand association. In addition, complexes 1-3 and 8-11 show long-lived luminescence from triplet excited states in the solid state at room temperature, according to our photophysical studies. Our quantum chemical density functional theory/multireference configuration interaction (DFT/MRCI) calculations reveal that the phosphorescence of 8 originates from a locally excited triplet state on the bdt ligand. They further suggest that the phenyl substituents of ITr are photochemically not innocent but can coordinate to the electron-deficient metal center of this trigonal complex in the excited state.

SOiCISCF: Combining SOiCI and iCISCF for Variational Treatment of Spin-Orbit Coupling

It has recently been shown that the SOiCI approach [Zhang, N.; J. Phys.: Condens. Matter 2022, 34, 224007], in conjunction with the spin-separated exact two-component relativistic Hamiltonian, can provide very accurate fine structures of systems containing heavy elements by treating electron correlation and spin-orbit coupling (SOC) on an equal footing. Nonetheless, orbital relaxations/polarizations induced by SOC are not yet fully accounted for due to the use of scalar relativistic orbitals. This issue can be resolved by further optimizing the still real-valued orbitals self-consistently in the presence of SOC, as done in the spin-orbit coupled CASSCF approach [Ganyushin, D.; et al. J. Chem. Phys. 2013, 138, 104113] but with the iCISCF algorithm [Guo, Y.; J. Chem. Theory Comput. 2021, 17, 7545-7561] for large active spaces. The resulting SOiCISCF employs both double group and time reversal symmetries for computational efficiency and the assignment of target states. The fine structures of p-block elements are taken as showcases to reveal the efficacy of SOiCISCF.

Polarity-Tunable Luminescence and Intersystem Crossing of a Zinc(II) Diimine Dithiolate Complex

Combined density functional theory and multireference configuration interaction methods including spin-orbit interactions have been employed to investigate the photophysical properties and deactivation pathways of a zinc diimine dithiolate complex involving the phenanthroline derivative bathocuproine and the dianionic dithiosquarate as chelating ligands. Zn(batho)(dtsq) is one of the few luminescent zinc complexes for which triplet emission had been reported in the solid state [Gronlund, P. Inorg. Chim. Acta 1995, 234, 13-18]. Because of the high dipole moment of the complex in the electronic ground state, ligand-to-ligand charge-transfer (LLCT) states experience strong hypsochromic shifts in polar media, while ligand-centered (LC) states are nearly unaffected. Rate constants for the thermally activated upconversion of the TLLCT population to the SLLCT state are promising due to a small singlet-triplet energy gap and the participation of the sulfur in the electronic excitation, but the TLLCT state is not the lowest-lying excited triplet state in ethanol solution. In addition to the TLLCT electronic structure, TLC(batho)' and TLC(dtsq) ππ* excitations form minima on the T1 potential energy surface. The SLLCT luminescence is expected to be quenched at the nanosecond time scale by the dark TLC(dtsq)ππ* state. Moreover, a TLC(dtsq)σπ* state has been identified, which leads to degradation of the compound. In mildly polar media, the dark triplet LC states are energetically inaccessible and the lowest excited singlet and triplet states clearly exhibit an LLCT character. However, their mutual spin-orbit coupling is reduced to the extent that reverse intersystem crossing is not very likely at room temperature. While Zn(diimine)(dithiolate) complexes continue to be perceived as an interesting substance class with potential application as emitters in electroluminescent devices, the particular Zn(batho)(dtsq) complex is not considered suitable for that purpose.

A relativistic configuration interaction method with general expansions and initial applications to electronic g-factors

A new implementation of the orbital-based two-component relativistic configuration interaction approach is reported and applied to calculations of the electronic g-shifts of three diatomic radicals: AlO, HgF, and PdH. The new implementation augments efficient routines for the calculation of nonrelativistic Hamiltonians with new vectorized routines for the calculation of the action of the one-electron spin-orbit operator and allows efficient calculations for the expansion of generalized active space type. The program makes full use of double group as well as time-reversal symmetry. Particle-hole reorganization of the operators is used to improve the efficiency for expansions with nearly fully occupied orbital spaces. The flexibility of the algorithm and program is used to investigate the convergence of electronic g-shifts for the three diatomic radicals as functions of the active space, states included in the orbital optimization, and excitation levels. It was possible to converge to the valence limits within a few percent using expansions containing up to quadruple excitations. However, when excitations from the core orbitals were added, it was not possible to demonstrate convergence to within a few percent with expansions containing at most 10 × 109 determinants.

Finding Design Principles of OLED Emitters through Theoretical Investigations of Zn(II) Carbene Complexes

In this work, Zn(II) carbene complexes carrying a dianionic 1,2-dithiolbenzene (dtb) or 1,2-diolbenzene (dob) ligand were investigated regarding their suitability as organic light-emitting diode (OLED) emitter. For the optimization of the complexes, density functional-based methods were used and frequency analyses verified the obtained structures as minima. All calculations were carried out including a polarizable continuum model to mimic solvent-solute interactions. Multireference configuration interaction methods were used to determine excitation energies, spin-orbit couplings, and luminescence properties. Rate constants of spin-allowed and spin-forbidden transitions were calculated according to a Fermi golden rule expression. Using carbene ligands with varying σ-donor and π-acceptor strengths, the luminescence is found to be tunable from yellow to orange/red to deep red/near-infrared. The calculated intersystem crossing (ISC) time constants indicate thermally activated delayed fluorescence (TADF) to be the main decay channel. In contrast to many d¹⁰ coinage metal complexes, a parallel orientation of dtb or dob and the carbene ligand is found to be highly favorable. For the complexes with a cyclic (alkyl)(amino) carbene (CAAC) or cyclic (amino)(aryl) carbene (CAArC) ligand, the S₁ and T₁ states have ligand-to-ligand charge-transfer (LLCT) character and are energetically close. The complex with a classical N-heterocyclic carbene (NHC) ligand has S₁ and T₁ states with mixed ligand-to-metal charge-transfer (LMCT)/LLCT character and is a very rare example in which the zinc ion contributes to the excitation.

The dominant nature of Herzberg-Teller terms in the photophysical description of naphthalene compared to anthracene and tetracene

The first order and second order corrected photoluminescence quantum yields are computed and compared to experiment for naphthalene in this manuscript discussing negative results. Results for anthracene and tetracene are recalled from previous work (Manian et al. in J Chem Phys 155:054108, 2021), and the results for all three polyacenes are juxtaposed to each other. While at the Franck-Condon point, each of the three noted polyacenes were found to possess a quantum yield near unity. Following the consideration of Herzberg-Teller effects, quantum yields stabilised for anthracene and tetracene to 0.19 and 0.08, respectively. Conversely, the second order corrected quantum yield for naphthalene was found to be 0.91. Analysis of this result showed that while the predicted non-radiative pathways correlate well with what should be expected, the approximation used to calculate second order corrected fluorescence, which yielded very positive results for many other molecular systems, here is unable to account for strong second order contributions, resulting in a grossly overestimated rate of fluorescence. However, substitution of an experimental radiative rate results in a quantum yield of 0.33. This work extols the importance of Herzberg-Teller terms in photophysical descriptions of chromophores, and highlights those cases in which a treatment beyond the above approximation is required.

A first principles examination of phosphorescence

This paper explores phosphorescence from a first principles standpoint, and examines the intricacies involved in calculating the spin-forbidden T 1 → S 0 transition dipole moment, to highlight that the mechanism is not as complicated to compute as it seems. Using gas phase acridine as a case study, we break down the formalism required to compute the phosphorescent spectra within both the Franck-Condon and Herzberg-Teller regimes by coupling the first triplet excited state up to the S 4 and T 4 states. Despite the first singlet excited state appearing as an L b state and not of nπ* character, the second order corrected rate constant was found to be 0.402 s-1, comparing well with experimental phosphorescent lifetimes of acridine derivatives. In showing only certain states are required to accurately describe the matrix elements as well as how to find these states, our calculations suggest that the nπ* state only weakly couples to the T 1 state. This suggest its importance hinges on its ability to quench fluorescence and exalt non-radiative mechanisms rather than its contribution to the transition dipole moment. A followup investigation into the T 1 → S 0 transition dipole moment's growth as a function of its coupling to other electronic states highlights that terms dominating the matrix element arise entirely from the inclusion of states with strong spin-orbit coupling terms. This means that while the expansion of the transition dipole moment can extend to include an infinite number of electronic states, only certain states need to be included.

Ultra-Long Lived Luminescent Triplet Excited States in Cyclic (Alkyl)(amino)carbene Complexes of Zn(II) Halides

The high element abundance and d10 electron configuration make ZnII -based compounds attractive candidates for the development of novel photoactive molecules. Although a large library of purely fluorescent compounds exists, emission involving triplet excited states is a rare phenomenon for zinc complexes. We have investigated the photophysical and -chemical properties of a series of dimeric and monomeric ZnII halide complexes bearing a cyclic (alkyl)(amino)carbene (cAAC) as chromophore unit. Specifically, [(cAAC)XZn(μ-X)2 ZnX(cAAC)] (X=Cl (1), Br (2), I (3)) and [ZnX2 (cAAC)(NCMe)] (X=Br (4), I (5)) were isolated and fully characterized, showing intense visible light photoluminescence under UV irradiation at 297 K and fast photo-induced transformation. At 77 K, the compounds exhibit improved stability allowing to record ultra-long lifetimes in the millisecond regime. DFT/MRCI calculations confirm that the emission stems from 3 XCT/LEcAAC states and indicate the phototransformation to be related to asymmetric distortion of the complexes by cAAC ligand rotation. This study enhances our understanding of the excited state properties for future development and application of new classes of ZnII phosphorescent complexes.

Vibronic effects accelerate the intersystem crossing processes of the through-space charge transfer states in the triptycene bridged acridine-triazine donor-acceptor molecule TpAT-tFFO

Quantum chemical studies employing combined density functional and multireference configuration interaction methods suggest five excited electronic states to be involved in the prompt and delayed fluorescence emission of TpAT-tFFO. Three of them, a pair of singlet and triplet charge transfer (CT) states (S1 and T1) and a locally excited (LE) triplet state (T3), can be associated with the (Me → N) conformer, the other two CT-type states (S2 and T2) form the lowest excited singlet and triplet states of the (Me → Ph) conformer. The two conformers, which differ in essence by the shearing angle of the face-to-face aligned donor and acceptor moieties, are easily interconverted in the electronic ground state whereas the reorganization energy is substantial in the excited singlet state, thus explaining the two experimentally observed time constants of prompt fluorescence emission. Forward and reverse intersystem crossing between the singlet and triplet CT states is mediated by vibronic spin-orbit interactions involving the LE T3 state. Low-frequency vibrational modes altering the distance and alignment of the donor and acceptor π-systems tune the S1 and T3 states (likewise S2 and T3) into and out of resonance. The enhancement of intersystem crossing due to the interplay of vibronic and spin-orbit coupling is considered a general feature of organic through-space charge-transfer thermally activated delayed fluorescence emitters.

SOiCI and iCISO: combining iterative configuration interaction with spin-orbit coupling in two ways

The near-exact iCIPT2 approach for strongly correlated systems of electrons, which stems from the combination of iterative configuration interaction (iCI, an exact solver of full CI) with configuration selection for static correlation and second-order perturbation theory (PT2) for dynamic correlation, is extended to the relativistic domain. In the spirit of spin separation, relativistic effects are treated in two steps: scalar relativity is treated by the infinite-order, spin-free part of the exact two-component (X2C) relativistic Hamiltonian, whereas spin-orbit coupling (SOC) is treated by the first-order, Douglas-Kroll-Hess-like SOC operator derived from the same X2C Hamiltonian. Two possible combinations of iCIPT2 with SOC are considered, i.e., SOiCI and iCISO. The former treats SOC and electron correlation on an equal footing, whereas the latter treats SOC in the spirit of state interaction, by constructing and diagonalizing an effective spin-orbit Hamiltonian matrix in a small number of correlated scalar states. Both double group and time reversal symmetries are incorporated to simplify the computation. Pilot applications reveal that SOiCI is very accurate for the spin-orbit splitting (SOS) of heavy atoms, whereas the computationally very cheap iCISO can safely be applied to the SOS of light atoms and even of systems containing heavy atoms when SOC is largely quenched by ligand fields.

Exciton Dynamics of a Diketo-Pyrrolopyrrole Core for All Low-Lying Electronic Excited States Using Density Functional Theory-Based Methods

Ab initio treatments of interexcited state internal conversion (IC) are more often than not missing from exciton dynamic descriptions, because of their inherent complexity. Here, we define "interexcited state IC" as a same-spin nonradiative transition between states i and j, where i ≠ j ≠ 0. Competing directly with multiexciton processes such as singlet fission or triplet photoupconversion, inclusion of this mechanism in the narrative of molecular photophysics would allow for strategic synthesis of chromophores for more efficient photon-harvesting applications. Herein, we present a robust formalism which can model these rates using density functional theory (DFT)-based methods within the Franck-Condon and Herzberg-Teller regime. Using an unsubstituted diketo-pyrrolopyrrole (DPP) core as a case study, we illustrate the exciton dynamics along the first four excited states for both singlet and triplet manifolds, showing ultrafast same-spin transfer mechanisms due to all excited states, excluding the first triplet level, being in close energetic proximity (within 0.8 eV of each other). The resulting electron same-spin rates outcompete the electron spin-flipping intersystem crossing (ISC) rates, with excitons firmly obeying Kasha's rule as they cascade down from the high-lying excited states toward the lower states. Furthermore, we calculated that only the first singlet excited state displayed a reasonable probability of triplet exciton generation, of ∼40%, with a near-zero chance of the exciton reverting to the singlet manifold once the electron-hole pair are of parallel spin.

Removing the Deadwood from DFT/MRCI Wave Functions: The p-DFT/MRCI Method

The combined density functional theory and multireference configuration interaction (DFT/MRCI) method is a powerful tool for the calculation of excited electronic states of large molecules. There exists, however, a large amount of superfluous configurations in a typical DFT/MRCI wave function. We show that this deadwood may be effectively removed using a simple configuration pruning algorithm based on second-order Epstein-Nesbet perturbation theory. The resulting method, which we denote p-DFT/MRCI, is shown to result in orders of magnitude saving in computational timings, while retaining the accuracy of the original DFT/MRCI method.

Modeling radiative and non-radiative pathways at both the Franck-Condon and Herzberg-Teller approximation level

Here, we present a concise model that can predict the photoluminescent properties of a given compound from first principles, both within and beyond the Franck-Condon approximation. The formalism required to compute fluorescence, Internal Conversion (IC), and Inter-System Crossing (ISC) is discussed. The IC mechanism, in particular, is a difficult pathway to compute due to difficulties associated with the computation of required bosonic configurations and non-adiabatic coupling elements. Here, we offer a discussion and breakdown on how to model these pathways at the Density Functional Theory (DFT) level with respect to its computational implementation, strengths, and current limitations. The model is then used to compute the photoluminescent quantum yield (PLQY) of a number of small but important compounds: anthracene, tetracene, pentacene, diketo-pyrrolo-pyrrole (DPP), and Perylene Diimide (PDI) within a polarizable continuum model. Rate constants for fluorescence, IC, and ISC compare well for the most part with respect to experiment, despite triplet energies being overestimated to a degree. The resulting PLQYs are promising with respect to the level of theory being DFT. While we obtained a positive result for PDI within the Franck-Condon limit, the other systems require a second order correction. Recomputing quantum yields with Herzberg-Teller terms yields PLQYs of 0.19, 0.08, 0.04, 0.70, and 0.99 for anthracene, tetracene, pentacene, DPP, and PDI, respectively. Based on these results, we are confident that the presented methodology is sound with respect to the level of quantum chemistry and presents an important stepping stone in the search for a tool to predict the properties of larger coupled systems.

Internal conversion of singlet and triplet states employing numerical DFT/MRCI derivative couplings: Implementation, tests, and application to xanthone

We present an efficient implementation of nonadiabatic coupling matrix elements (NACMEs) for density functional theory/multireference configuration interaction (DFT/MRCI) wave functions of singlet and triplet multiplicity and an extension of the Vibes program that allows us to determine rate constants for internal conversion (IC) in addition to intersystem crossing (ISC) nonradiative transitions. Following the suggestion of Plasser et al. [J. Chem. Theory Comput. 12, 1207 (2016)], the derivative couplings are computed as finite differences of wave function overlaps. Several measures have been taken to speed up the calculation of the NACMEs. Schur's determinant complement is employed to build up the determinant of the full matrix of spin-blocked orbital overlaps from precomputed spin factors with fixed orbital occupation. Test calculations on formaldehyde, pyrazine, and xanthone show that the mutual excitation level of the configurations at the reference and displaced geometries can be restricted to 1. In combination with a cutoff parameter of t_norm = 10^-8 for the DFT/MRCI wave function expansion, this approximation leads to substantial savings of cpu time without essential loss of precision. With regard to applications, the photoexcitation decay kinetics of xanthone in apolar media and in aqueous solution is in the focus of the present work. The results of our computational study substantiate the conjecture that S₁ T₂ reverse ISC outcompetes the T₂ ↝ T₁ IC in aqueous solution, thus explaining the occurrence of delayed fluorescence in addition to prompt fluorescence.

Intersystem crossing processes in the 2CzPN emitter: a DFT/MRCI study including vibrational spin-orbit interactions

Multireference quantum chemical calculations were performed in order to investigate the (reverse) intersystem crossing ((R)ISC) mechanisms of 4,5-di(9H-carbazol-9-yl)-phthalonitrile (2CzPN). A combination of density funcional theory (DFT) and multireference configuration interaction methods (MRCI) was used. The excellent agreement of the computed absorption spectrum with available experimental absorption spectra lends confidence to the chosen computational protocol. Vertically, two triplet excited states (T₁ and T₂) are found below the S₁ state. At the excited state minima, the calculated adiabatic energies locate only the T₁ state below the S₁ state. The enhanced charge transfer (CT) character of the geometrically relaxed excited states causes their mutual (direct) spin-orbit coupling (SOC) interaction to be low. Contributions of vibronic SOC to the (R)ISC probability, evaluated by a Herzberg-Teller-like procedure for a temperature of 300 K, are small but not negligible. For ISC, the S₁→ T₁ channel is the fastest (8 × 10⁶ s^-1), while the S₁→ T₂ channel is found to be thermally activated (9 × 10⁴ s^-1) and less efficient when proceeding from the adiabatic S₁ state. Our calculations also reveal, however, a barrierless S₁→ T₂ ISC pathway near the Franck-Condon region. RISC is found to essentially proceed via the T₁→ S₁ channel, with a rate constant of (3 × 10⁴ s^-1) if our adiabatic singlet-triplet energy gap in vacuum (ΔE_ST = 0.12 eV) is employed. Shifting the potentials to match two experimentally reported singlet-triplet energy gaps in toluene (ΔE_ST = 0.21 and 0.31 eV, respectively) leads to a drastic reduction of the computed rate constant by up to 4 orders of magnitude. The T₂ state is not expected to play a major role in mediating triplet-singlet transitions in 2CzPN unless it is directly populated by hot excitons. No indication for a strong vibronic coupling of the T₂ and T₁ potentials is found, which could help overcome the negative exponential dependence of the RISC rate constant on the magnitude of the energy gap.

Behind the scenes of spin-forbidden decay pathways in transition metal complexes

The interpretation of the ultrafast photophysics of transition metal complexes following photo-absorption is quite involved as the heavy metal center leads to a complicated and entangled singlet-triplet manifold. This opens up multiple pathways for deactivation, often with competitive rates. As a result, intersystem crossing (ISC) and phosphorescence are commonly observed in transition metal complexes. A detailed understanding of such an excited-state structure and dynamics calls for state-of-the-art experimental and theoretical methodologies. In this review, we delve into the inability of non-relativistic quantum theory to describe spin-forbidden transitions, which can be overcome by taking into account spin-orbit coupling, whose importance grows with increasing atomic number. We present the quantum chemical theory of phosphorescence and ISC together with illustrative examples. Finally, a few applications are highlighted, bridging the gap between theoretical studies and experimental applications, such as photofunctional materials.

Effect of spin-orbit coupling on strong field ionization simulated with time-dependent configuration interaction

Time-dependent configuration interaction with a complex absorbing potential has been used to simulate strong field ionization by intense laser fields. Because spin-orbit coupling changes the energies of the ground and excited states, it can affect the strong field ionization rate for molecules containing heavy atoms. Configuration interaction with single excitations (CIS) has been employed for strong field ionization of closed shell systems. Single and double excitation configuration interaction with ionization (CISD-IP) has been used to treat ionization of degenerate states of cations on an equal footing. The CISD-IP wavefunction consists of ionizing single (one hole) and double (two hole/one particle) excitations from the neutral atom. Spin-orbit coupling has been implemented using an effective one electron spin-orbit coupling operator. The effective nuclear charge in the spin-orbit coupling operator has been optimized for Ar⁺, Kr⁺, Xe⁺, HX⁺ (X = Cl, Br, and I). Spin-orbit effects on angular dependence of the strong field ionization have been studied for HX and HX⁺. The effects of spin-orbit coupling are largest for ionization from the π orbitals of HX⁺. In a static field, oscillations are seen between the ²Π_3/2 and ²Π_1/2 states of HX⁺. For ionization of HX⁺ by a two cycle circularly polarized pulse, a single peak is seen when the maximum in the carrier envelope is perpendicular to the molecular axis and two peaks are seen when it is parallel to the axis. This is the result of the greater ionization rate for the π orbitals than for the σ orbitals.

Understanding the luminescence properties of Cu(i) complexes: a quantum chemical perusal

Electronic structures and excited-state properties of Cu(i) complexes with varying coordination numbers have been investigated by means of advanced quantum chemical methods. The computational protocol employs density functional-based methods for geometry optimizations and vibrational analyses including solvent effects through continuum models. Excitation energies, spin-orbit couplings and luminescence properties are evaluated using multireference configuration interaction methods. Rate constants of spin-allowed and spin-forbidden transitions have been determined according to the Fermi golden rule. The computational results for the 4-coordinate (DPEPhos)Cu(PyrTet), the 3-coordinate [IPr-Cu-Py₂]⁺, and the linear CAAC^Me₂-Cu-Cl complexes agree well with experimental absorption and emission wavelengths, intersystem crossing (ISC) time constants, and radiative lifetimes in liquid solution. Spectral shifts on the ligand-to-ligand charge transfer (LLCT) and metal-to-ligand charge transfer (MLCT) transitions caused by the polarity of the environment are well represented by the continuum models whereas the shifts caused by pseudo-Jahn-Teller distortions in the MLCT states are too pronounced in comparison to solid-state data. Systematic variation of the ligands in linear Cu(i) carbene complexes shows that only those complexes with S₁ and T₁ states of LLCT character possess sufficiently small singlet-triplet energy gaps ΔE_ST to enable thermally activated delayed fluorescence (TADF). Complexes whose S₁ and T₁ wavefunctions are dominated by MLCT excitations tend to emit phosphorescence instead. Unlike the situation in metal-free TADF emitters, the presence of low-lying locally excited triplet states does not promote ISC. These states rather hold the danger of trapping the excitation with nonradiative deactivation being the major deactivation channel.

Visible Light-Induced Homolytic Cleavage of Perfluoroalkyl Iodides Mediated by Phosphines

In an effort to explain the experimentally observed variation of the photocatalytic activity of t Bu 3 P, n Bu 3 P and (MeO) 3 P in the blue-light regime [Helmecke et al., Org. Lett. 21 (2019) 7823], we have explored the absorption characteristics of several phosphine- and phosphite-IC 4 F 9 adducts by means of relativistic density functional theory and multireference configuration interaction methods. Based on the results of these computational and complementary experimental studies, we offer an explanation for the broad tailing of the absorption of t Bu 3 P-IC 4 F 9 and (MeO) 3 P-IC 4 F 9 into the visible-light region. Larger coordinate displacements of the ground and excited singlet potential energy wells in n Bu 3 P-IC 4 F 9 , in particular with regard to the P-I-C bending angle, reduce the Franck-Condon factors and thus the absorption probability compared to t Bu 3 P-IC 4 F 9 . Spectroscopic and computational evaluation of conformationally flexible and locked phosphites suggests that the reactivity of (MeO) 3 P may be the result of oxygen lone-pair participation and concomitant broadening of absorption. The proposed mechanism for the phosphine-catalyzed homolytic C-I cleavage of perfluorobutane iodide involves S₁ ← S₀ absorption of the adduct followed by intersystem crossing to the photochemically active T 1 state.

On the photophysical properties of IrIII, PtII, and PdII (phenylpyrazole) (phenyldipyrrin) complexes

The absorption and emission characteristics of (ppz)2(dipy)IrIII, (ppz)(dipy)PtII and (ppz)(dipy)PdII, where ppz stands for phenylpyrazole and dipy for a phenyl meso-substituted dipyrrin ligand, have been investigated by means of combined density functional theory and multireference configuration interaction including scalar relativistic and spin-orbit coupling effects. These results were compared with experimental spectra. The complexes exhibit a high density of low-lying electronically excited states originating from ligand-centered (LC) and metal-to-ligand charge transfer (MLCT) states involving the dipyrrin ligand. In addition, metal-centered (MC) states are found to be low-lying in the Pd complex. In all three cases, the first strong absorption band and the phosphorescence emission band stem from LC excitations on the dipyrrin ligand with small MLCT contributions. The MLCT states show more pronounced relaxation effects than the LC states, with the consequence that the first excited state with predominant singlet multiplicity is of SMLCT/LC type in the heavier Ir and Pt complexes. Substantial spin-orbit coupling between SMLCT/LC and TLC enables fast and efficient intersystem crossing (ISC) and a high triplet quantum yield. Phosphorescence rate constants are rather small in accord with the dominant LC character of the transitions. Out-of-plane distortion promotes nonradiative decay of the excited state population via the MC states thus explaining the lower phosphorescence quantum yield of the Pt complex. The spectral properties of the Pd complex are different in many aspects. Optimization of the S1 state yields a dipyrrin intraligand charge transfer (ILCT) state with highly distorted nuclear arrangement in the butterfly conformers leading to nonradiative deactivation. In contrast, the primarily excited SLC state and the SMLCT/LC state of the twist conformer have nearly equal adiabatic excitation energies. The lack of a driving force toward the SMLCT/LC minimum, the high fluorescence rate constant of the bright SLC state and its moderately efficient ISC to the triplet manifold explain the experimentally observed dual emission of the Pd complex at room temperature.

DFT/MRCI assessment of the excited-state interplay in a coumarin-schiff Mg2+ fluorescent sensor

Fluorescent sensors with selectivity and sensitivity to metal ions are an active field in supramolecular chemistry for biochemical, analytical, and environmental problems. Mg²⁺ is one of the most abundant divalent ions in the cell, and it plays a critical role in many biological processes. Coumarin-based sensors are widely used as desirable fluorophore and binding moieties showing a remarkable sensitivity and fluorometric enhancement for Mg²⁺ . In this work, density functional theory/multireference configuration interaction (DFT/MRCI) calculations were performed in order to understand the sensing behavior of the organic fluorescent sensor 7-hydroxy-4-methyl-8-((2-(pyridin-2-yl)hydrazono)methyl)-2H-chromen-2-one (PyHC) in ethanol to solvated Mg²⁺ ions. The computed optical properties reproduce well-reported experimental data. Our results suggest that after photoexcitation of the free PyHC, a photo-induced electron transfer (PET) mechanism may compete with the fluorescence decay to the ground state. In contrast, this PET channel is no longer available in the complex with Mg²⁺ making the emissive decay more efficient. © 2019 Wiley Periodicals, Inc.

Structure-Emission Property Relationships in Cyclometalated Pt(II) β-Diketonate Complexes

Extending the ligand π-system of phosphorescent (C^∧C*) or (C^∧N) cyclometalated platinum(II) β-diketonate complexes can lead to large and seemingly abrupt variations of the photophysical properties such as triplet quantum yields and phosphorescence lifetimes. Quantum chemical studies using methods including elements from density functional theory (DFT) and multireference configuration interaction (MRCI) as well as spin-orbit coupling (SOC) provide a rationale for these observations. In the Franck-Condon region, the first excited singlet states (S₁) of these complexes are characterized by mixed metal-to-ligand charge-transfer (MLCT) and ligand-centered (LC) excitations. With increasing extension of the effective π-system, the lowest-lying triplet state yields more and more LC character, thus leading to a decrease of the phosphorescence rate constant. The ability to undergo efficient intersystem crossing from S₁ to T₁ is not diminished as the S₁ state largely retains its character. In the N-heterocyclic carbene (NHC) complexes investigated here, at least two triplet states are found energetically below the S₁ state. Out-of-plane distortion enhances the probability for nonradiative decay of the triplet population. In the smaller compounds emitting in the violet or blue spectral region, the phosphorescent state is separated from the lowest-lying dark metal-centered (MC) triplet state by a small barrier only, explaining their experimentally observed low photoluminescence quantum yields in liquid solution. The semiempirical DFT/MRCI-R2018 Hamiltonian employed in our studies proves well-suited for investigating the absorption and emission properties of these platinum(II) complexes. Generally, good agreement is observed between our calculated data and the experimental findings.

Computer-Aided Design of Luminescent Linear N-Heterocyclic Carbene Copper(I) Pyridine Complexes

Multireference configuration interaction methods including spin-orbit interactions have been employed to investigate the photophysical properties of various linear N-heterocyclic carbene (NHC) copper(I) pyridine complexes with the aim of designing performant thermally activated delayed fluorescence (TADF) emitters for use in organic-light-emitting diodes. Our theoretical results indicate that this structural motif is very favorable for generating excited triplet states with high quantum yield. The first excited singlet (S_MLCT) and corresponding triplet state (T_MLCT) are characterized by d_σ → π_Py metal-to-ligand charge-transfer (MLCT) excitations. Efficient intersystem crossing (ISC) and reverse ISC (rISC) between these states is mediated by a near-degenerate second triplet state (T_MLCT/LC) with large d_π → π_Py contributions. Spin-vibronic coupling is strong and is expected to play a major role in the (r)ISC processes. The calculations reveal, however, that the luminescence is effectively quenched by locally excited triplet states if the NHC ligand carries two diisopropylphenyl (DIPP) substituents. When DIPP is replaced with 1-adamantyl residues, this quenching process is suppressed and TADF in the UV spectral regime is predicted to proceed at a rate of about 1/μs. The introduction of +I substituents on the carbene and -M substituents on the pyridine allows tuning of the emission wavelength from the UV to the blue-green or green spectral region.

Rotationally Assisted Spin-State Inversion in Carbene-Metal-Amides Is an Artifact

High-level quantum chemical calculations, presented in this Letter, show that the unusual luminescence properties of the high-performance thermally activated delayed fluorescence emitter CMA1, observed recently by Di et al. [ Science 2017, 356, 159-163 ], can be explained without resorting to the rotationally assisted spin-state inversion mechanism proposed by these authors. Multiconfiguration and relativistic effects lead to fast and efficient thermal equilibration of the excited singlet and triplet populations of this linear gold complex even for coplanar orientations of the ligands. The calculations predict S₁ ⇝ T₁ intersystem crossing outcompetes the submicrosecond fluorescence by 2 orders of magnitude, thus quenching prompt fluorescence. The significant time- and environment-dependent shifts of the CMA1 luminescence, observed in experiment, are attributed to effects of (hindered) solvent reorganization.

Theory and Calculation of the Phosphorescence Phenomenon

Phosphorescence is a phenomenon of delayed luminescence that corresponds to the radiative decay of the molecular triplet state. As a general property of molecules, phosphorescence represents a cornerstone problem of chemical physics due to the spin prohibition of the underlying triplet-singlet emission and because its analysis embraces a deep knowledge of electronic molecular structure. Phosphorescence is the simplest physical process which provides an example of spin-forbidden transformation with a characteristic spin selectivity and magnetic field dependence, being the model also for more complicated chemical reactions and for spin catalysis applications. The bridging of the spin prohibition in phosphorescence is commonly analyzed by perturbation theory, which considers the intensity borrowing from spin-allowed electronic transitions. In this review, we highlight the basic theoretical principles and computational aspects for the estimation of various phosphorescence parameters, like intensity, radiative rate constant, lifetime, polarization, zero-field splitting, and spin sublevel population. Qualitative aspects of the phosphorescence phenomenon are discussed in terms of concepts like structure-activity relationships, donor-acceptor interactions, vibronic activity, and the role of spin-orbit coupling under charge-transfer perturbations. We illustrate the theory and principles of computational phosphorescence by highlighting studies of classical examples like molecular nitrogen and oxygen, benzene, naphthalene and their azaderivatives, porphyrins, as well as by reviewing current research on systems like electrophosphorescent transition metal complexes, nucleobases, and amino acids. We furthermore discuss modern studies of phosphorescence that cover topics of applied relevance, like the design of novel photofunctional materials for organic light-emitting diodes (OLEDs), photovoltaic cells, chemical sensors, and bioimaging.

Combining the spin-separated exact two-component relativistic Hamiltonian with the equation-of-motion coupled-cluster method for the treatment of spin–orbit splittings of light and heavy elements

An accurate and efficient treatment of spin–orbit splittings has been achieved by combining the sf-X2C+soc-DKH1 Hamiltonian with the equation-of-motion coupled-cluster method.

Effect of Conjugation Pathway in Metal-Free Room-Temperature Dual Singlet-Triplet Emitters for Organic Light-Emitting Diodes

Metal-free dual singlet-triplet organic light-emitting diode (OLED) emitters can provide direct insight into spin statistics, spin correlations and spin relaxation phenomena, through a comparison of fluorescence to phosphorescence intensity. Remarkably, such materials can also function at room temperature, exhibiting phosphorescence lifetimes of several milliseconds. Using electroluminescence, quantum chemistry, and electron paramagnetic resonance spectroscopy, we investigate the effect of the conjugation pathway on radiative and nonradiative relaxation of the triplet state in phenazine-based compounds and demonstrate that the contribution of the phenazine nπ* excited state is crucial to enabling phosphorescence.

Spin-Free CC2 Implementation of Induced Transitions between Singlet Ground and Triplet Excited States

In most organic molecules, phosphorescence has its origin in transitions from triplet exited states to the singlet ground state, which are spin-forbidden in nonrelativistic quantum mechanics. A sufficiently accurate description of phosphorescence lifetimes for molecules that contain only light elements can be achieved by treating the spin-orbit coupling (SOC) with perturbation theory (PT). We present an efficient implementation of this approach for the approximate coupled cluster singles and doubles model CC2 in combination with the resolution-of-the-identity approximation for the electron repulsion integrals. The induced oscillator strengths and phosphorescence lifetimes from SOC-PT are computed within the response theory framework. In contrast to previous work, we employ an explicitly spin-coupled basis for singlet and triplet operators. Thereby, a spin-orbital treatment can be entirely avoided for closed-shell molecules. For compounds containing only light elements, the phosphorescence lifetimes obtained with SOC-PT-CC2 are in good agreement with those of exact two-component (X2C) CC2, whereas the calculations are roughly 12 times faster than with X2C. Phosphorescence lifetimes computed for two thioketones with the SOC-PT-CC2 approach agree very well with reference results from experiment and are similar to those obtained with multireference spin-orbit configuration interaction and with X2C-CC2. An application to phosphorescent emitters for metal-free organic light-emitting diodes (OLEDs) with almost 60 atoms and more than 1800 basis functions demonstrates how the approach extends the applicability of coupled cluster methods for studying phosphorescence. The results indicate that other decay channels like vibrational relaxation may become important in such systems if lifetimes are large.