Reverse Intersystem Crossing: Mechanisms & Kinetics

• Reverse intersystem crossing (RISC) is the non-radiative transition from triplet to singlet states, critical for enabling TADF in organic emitters.
• It leverages spin–orbit coupling and vibronic interactions to overcome the singlet–triplet gap, with efficiency determined by molecular symmetry and environmental effects.
• Design strategies focus on minimizing ΔEST and optimizing donor–acceptor architectures to boost rISC rates up to the order of 10⁷ s⁻¹ for advanced OLED applications.

Reverse intersystem crossing (RISC) denotes the non-radiative up-conversion of molecular excited-state population from a triplet ($T$) to a singlet ($S$) manifold, enabling repopulation of bright singlet states from energetically proximate triplets. This mechanism is foundational to thermally activated delayed fluorescence (TADF), fluorescence via higher triplets (FvHT), and related up-conversion phenomena utilized in organic light-emitting diodes (OLEDs) and photofunctional organic materials. RISC is driven by spin-orbit coupling (SOC), vibronic effects, and can be modulated by environmental dielectric response, nuclear motion, and molecular symmetry. Efficiency of RISC is determined by the relative singlet-triplet gap ($\Delta E_\mathrm{ST}$), SOC magnitude, reorganization energy, and the detailed kinetics of competing non-radiative processes.

1. Fundamental Mechanisms, Definitions, and Energy-Level Structure

Ordinary intersystem crossing (ISC) represents the spin-forbidden relaxation from an excited singlet (typically $S_1$) to a triplet ($T_1$). RISC is the reverse process, thermally or vibronically assisted, resulting in $T_1 \rightarrow S_1$ transfer. In photophysical and optoelectronic contexts, these transitions are governed by the spin–orbit coupling operator, and are typically slower than allowed radiative transitions due to their spin-forbidden nature.

Key energy-level parameters are as follows:

• Singlet and triplet states with strong charge-transfer (CT) character: $^1CT$ or $S_1$, $^3CT$ or $T_1$.
• The singlet–triplet gap: $S$0, which needs to be $S$1 eV for efficient RISC at room temperature (Gillett et al., 2021).
• In higher triplet mechanisms such as FvHT, RISC can occur from $S$2 (energy gap $S$3 meV), followed by internal conversion (IC) to $S$4 (Sato et al., 2016).

A schematic pathway for a four-state system:

• $S$5

For symmetric systems (SC-TADF, iST), specialized selection rules and orbital structures allow direct or barrierless RISC or even reversal of Hund’s rule ordering (INVEST systems where $S$6) (Karak et al., 2024).

2. Kinetic Rate Expressions and Quantum Theories

The rate of RISC is described under various frameworks, depending on the level of electronic-vibrational interaction and thermal fluctuation considered.

• Semiclassical Marcus Theory (applicable to charge-transfer systems) (Gillett et al., 2021):

$S$7

where $S$8 is the SOC matrix element, $S$9 is total reorganization energy, $\Delta E_\mathrm{ST}$0.

In the high-temperature (Arrhenius) limit:

$\Delta E_\mathrm{ST}$1

• Fermi’s Golden Rule (Vibronic Coupling Picture) (Sato et al., 2016, Karak et al., 2024):

$\Delta E_\mathrm{ST}$2

where $\Delta E_\mathrm{ST}$3 is the Franck–Condon–weighted density of states.

• Generating-Function and Wigner Averaging (phase-space formalism) (Karak et al., 2024):

$\Delta E_\mathrm{ST}$4

incorporating nuclear coordinate sampling and dynamical modulations.

Numerical rates highlight the effect of the medium: in a weakly polar solvent (toluene, $\Delta E_\mathrm{ST}$5), $\Delta E_\mathrm{ST}$6 s$\Delta E_\mathrm{ST}$7, versus $\Delta E_\mathrm{ST}$8 s$\Delta E_\mathrm{ST}$9 in vacuum for TXO-TPA (Gillett et al., 2021). In INVEST emitters, $S_1$0 s$S_1$1 (Karak et al., 2024).

3. Role of Molecular and Environmental Structure

Donor–Acceptor Architectures and Dipole Moments

Efficient TADF RISC is realized in molecules exhibiting strong D–A separation, minimizing orbital overlap to keep $S_1$2 small, with partial local-exciton (LE) character preserved to ensure sufficient SOC (Gillett et al., 2021).

Large changes in dipole moment $S_1$3 upon excitation favor strong stabilization by polar media. Conjugated molecules like TXO-TPA show an environment-induced reduction of $S_1$4 by $S_1$5 eV, facilitating RISC (Gillett et al., 2021). Table 1 summarizes key variables involved in environmental tuning:

Variable	Effect on RISC	Characteristic Value/Scale
$S_1$6	Activation barrier	$S_1$7 eV (optimal)
$S_1$8	Dielectric stabilization	$S_1$9 D (effective tuning)
$T_1$0	Outer sphere (solvent)	$T_1$1 meV (toluene/TXO-TPA)
$T_1$2	Host dielectric const.	$T_1$3

Vibrational Modes and Vibronic Coupling

RISC efficiency is further enhanced by tuning vibrational modes that mediate interstate coupling. For TXO-TPA in toluene, impulsive Raman measurements reveal vibrational modes at 412 and 813 cm$T_1$4 as fingerprints of the fully relaxed CT product state, supporting fast reorganization (Gillett et al., 2021).

Off-diagonal vibronic coupling constants (VCCs) control the relative rates of IC and RISC in higher triplet channels. In BD1, VCCs support ultrafast IC from $T_1$5 to $T_1$6 and suppress IC from $T_1$7 to lower triplets, channeling population into the RISC-allowed $T_1$8 transition (Sato et al., 2016).

4. Extended RISC Pathways and Generalized Frameworks

Beyond $T_1$9 RISC characteristic of standard TADF, higher-energy mechanisms such as FvHT involve conversion from $T_1 \rightarrow S_1$0, followed by S$T_1 \rightarrow S_1$1S$T_1 \rightarrow S_1$2 internal conversion. This process is facilitated by:

• Small energy gap $T_1 \rightarrow S_1$3 meV (sufficient for thermal activation at room temperature).
• Symmetry-allowed SOC for $T_1 \rightarrow S_1$4, but symmetry- or overlap-suppressed transitions for $T_1 \rightarrow S_1$5 and $T_1 \rightarrow S_1$6.
• Pseudo-degenerate frontier orbitals enabling specific construction/cancellation of transition densities, optimizing up-conversion (Sato et al., 2016).

SC-TADF and iST systems (where $T_1 \rightarrow S_1$7 at equilibrium) exemplify alternative frameworks that also leverage RISC, with the unique case of INVEST systems allowing barrierless or “downhill” RISC due to negative $T_1 \rightarrow S_1$8 (Karak et al., 2024). All cases are subsumed by the “fluorescence via RISC (FvRISC)” superordinate classification.

5. Dielectric, Dynamic, and Nuclear Effects

The environment impacts RISC both by shifting $T_1 \rightarrow S_1$9 and by modulating the reorganization energy. For polar solvents, the CT state stabilization ($^1CT$00.3 eV) narrows $^1CT$1 substantially (Gillett et al., 2021). Marcus-type outer-sphere calculations yield:

$^1CT$2

where $^1CT$3 is an effective radius.

Explicit QM/MM molecular dynamics show large thermal fluctuations in the S$^1CT$4–T$^1CT$5 energy gap (standard deviation $^1CT$60.29 eV), requiring ensemble or phase-space models for accurate rISC rate predictions (Gillett et al., 2021, Karak et al., 2024).

Wigner phase-space sampling reveals that, even in INVEST systems, the fraction of configurations supporting $^1CT$7 varies with nuclear geometry, and out-of-plane dihedral (puckering) motion directly modulates the local $^1CT$8 (Karak et al., 2024). This underlies the observed weak temperature dependence and robustness of barrierless RISC in such systems.

6. Quantitative Comparisons and Design Guidelines

Computed and experimentally measured rates for rISC and ISC for various systems are summarized in Table 2:

Material/Env.	$^1CT$9 (s$S_1$0)	Activation Energy $S_1$1 (eV)
TXO-TPA (vac)	$S_1$2	$S_1$3
TXO-TPA (tol)	$S_1$4	$S_1$5
INVEST 1 (calc)	$S_1$6	$S_1$7
INVEST 2 (calc)	$S_1$8	$S_1$9

RISC rates can thus be boosted by 2–3 orders of magnitude through environment and molecular engineering.

Design recommendations include (Gillett et al., 2021, Karak et al., 2024, Sato et al., 2016):

• Maximize $^3CT$0 and D–A separation for strong CT character and dielectric tunability.
• Minimize $^3CT$1 over the thermal distribution (not just at equilibrium).
• Retain partial LE character for SOC compatibility.
• Employ host matrices of moderate polarity ($^3CT$2) with dynamic reorganization capacity.
• Target vibronic modes in the 400–800 cm$^3CT$3 range for effective spin–vibronic coupling.
• Optimize double-excitation and multiresonance character to lower $^3CT$4 in INVEST-type emitters.
• Favor rigid architectures that accommodate the key normal modes influencing the singlet–triplet gap.

7. Limitations, Open Problems, and Synthesis Implications

Limitations in current theoretical and experimental approaches include the commonly neglected fast (optical) component of solvent response, incomplete treatment of higher-lying triplet and singlet states except via indirect mechanisms, and classical descriptions for most vibrational contributions except for key modes.

Theoretical models such as DFT (PBEh-3c) with QM/MM forces are benchmarked by higher-level methods (e.g., CC2), which confirm the accuracy of trends but introduce systematic absolute energy offsets.

A central theme affirmed by phase-space studies (Karak et al., 2024) is that RISC efficiency is not dictated solely by static ground-state properties, but by a thermodynamically weighted ensemble of molecular configurations in which the singlet–triplet energy gap can be periodically inverted or minimized by nuclear motion. Design of next-generation rISC emitters thus necessitates a dynamic, multidisciplinary approach integrating molecular electronic structure, environment, vibrational dynamics, and device-level considerations.