Open Access Article
Federica
Fina†
ab,
Caterina
Bellatreccia†
ab,
Xia
Wu
c,
Pier Giorgio
Cozzi
ab,
Alessandro
Troisi
c,
Sergei
Vinogradov
de and
Paola
Ceroni
*ab
aAlma Mater Studiorum-Dipartimento di Chimica “G. Ciamician”, Università di Bologna, Via Gobetti 83, Bologna 40129, Italy. E-mail: paola.ceroni@unibo.it
bCenter for Chemical Catalysis - C3, Alma Mater Studiorum, Università di Bologna, Via Gobetti 83, Bologna 40129, Itlay
cDepartment of Chemistry, University of Liverpool, Liverpool L69 7ZD, UK
dDepartment of Biochemistry and Biophysics, Perelman School of Medicine, University of Pennsylvania, Philadelphia, USA
eDepartment of Chemistry, School of Arts and Sciences University of Pennsylvania, Philadelphia, USA
First published on 7th October 2025
Cyanoarene chromophores exhibiting thermally activated delayed fluorescence (TADF) are increasingly used in photoredox catalysis. At high concentrations of organic substrates, which are typically employed in preparative synthesis, the primary photoinduced electron transfer (PeT) steps in the photocatalytic processes can involve both singlet (S1) and triplet (T1) excited states of TADF chromophores, despite very short lifetimes (nanoseconds) of the former. However, the difference between the reactivities of these states is not well understood, while being critically important for the photocatalytic process. In this work, three representative TADF chromophores were examined in reductive and oxidative PeT quenching reactions. First, using kinetic simulations, we assert that Stern–Volmer quenching plots based on the experimentally measured prompt and delayed fluorescence lifetimes, but not integrated intensities, yield accurate bimolecular rate constants for the PeT quenching reactions involving S1 and T1 excited states. Secondly, experimental measurements of prompt and delayed fluorescence reveal significantly higher quenching constants for reductive quenching of S1 compared to T1 states, while for oxidative quenching the rate constants are nearly equal. Electronic structure calculations provide insight into the difference between the PeT rates for reductive quenching, suggesting that it might stem from the different spatial hole–electron distributions in S1 and T1 states. Taken together, our findings bring crucial information about the photocatalytic process involving TADF chromophores that should aid the design of the next-generation of TADF photocatalysts.
The widespread use of TADF chromophores in photocatalysis and, more specifically, in photoredox catalysis, stems from: (i) availability via easy synthetic procedure; (ii) tunable redox properties by modification of either the core or the acceptor units;15 (iii) close energetic proximity of S1 and T1 states, so that only a small fraction of energy is lost when S1 transitions to T1via isc, while, at high substrate concentrations, both states can participate in the photoinduced electron transfer (PeT) reactions.16–19
Nevertheless, despite the interest in using TADF chromophores as photocatalysts, the quenching of their S1 and T1 states by PeT is not well understood, while being critically important for the photocatalytic process. For example, recently we observed that reductive quenching of the prompt fluorescence of 4DPAIPN (Fig. 1, kqp = 2.7 × 109 M−1 s−1) by PeT from 1,3-dimethyl-2-phenyl-2,3-dihydro-1H-benzimidazole (BIH) is characterized by a more than 10-fold higher quenching constant than quenching of the delayed fluorescence (kqd = 1.5 × 108 M−1 s−1).20 Similarly, we observed faster quenching of S1 compared to T1 in the case of the quenching of another TADF chromophore, namely 9,9’-(sulfonylbis(pyrimidine-5,2-diyl))bis(3,6-di-tert-butyl-9Hcarbazole) pDTCz-DPmS by four different electron donors.21 Since the energies of the S1 and T1 states of 4DPAIPN are close (vide infra), the driving forces (ΔG) of the PeT processes involving these states should be similar, and, therefore, one would expect only a small difference between the corresponding quenching rate constants. The objective of the present study was to understand the origin of this anomalous behaviour by analysing other isophthalonitrile-based TADF chromophores (Fig. 1) and other quenchers to determine whether the difference persists across a broader set of compounds and to rationalize the experimental findings with the help of kinetic simulations and computational analysis.
![]() | (1) |
![]() | (2) |
![]() | (3) |
![]() | (4) |
We performed kinetic simulations by numerical integration of the rate equations of a TADF system. In these simulations, a specific combination of kinetic constants of Fig. 2 (input parameters) yields simulated traces of S1 and T1, which can be fitted with biexponential functions to obtain τ1S, τ2S, τ1T and τ2T (output values, see SI for more details). The time constants obtained by the fitting of the simulated S1 trace (τ1S, τ2S) correspond to the lifetimes obtained by fitting of the fluorescence decay in a real experiment. As expected, the simulated decay of S1 in absence of the quencher is always in agreement with the exact solution, meaning that τ1S and τ2S match with τp, τd calculated viaeqn (3) and (4) based on the input parameters (Table S1). In addition to that, our simulated data show that, under any combination of kinetic constants, the longer lifetime obtained by the fitting of S1 trace (τ2S) matches with the decay lifetime of the T1 trace (τ2T), demonstrating that the delayed fluorescence lifetime mirrors T1 decay. The importance of this finding stems from the fact that for most common TADF organic photocatalysts T1– >S0 phosphorescence, which would be considered a direct signature of the T1 state, is not observable at ambient temperatures and measuring evolution of the triplet state by transient absorption experiments23–25 is not straightforward to set-up and/or interpret. Instead, the analysis of delayed fluorescence is simpler and more accessible from the experimental point of view.
In the presence of the quencher, it is necessary to assess how the quenching constants kqp and kqd, experimentally determined by prompt and delayed fluorescence quenching, can be correlated to the quenching constants of S1 and T1, namely kqS and kqT. To address this issue, we run kinetic simulations introducing the quenching terms in the differential equations of S1 and T1. Each rate constant shown in Fig. 2 was varied over three orders of magnitude in the simulations (see Table S1). For each combination of kinetic constants, the concentration of the quencher [Q] was varied to generate the respective [S1] decays, which were then fitted (using the non-linear least squares method) with a bi-exponential equation:
| [S1] = A1 × e(−t/τ1S)+A2 × e(−t/τ2S) | (5) |
The extracted values of τ1S (shorter lifetime) and τ2S (longer lifetime) were used to generate simulated Stern–Volmer plots, from which the respective output quenching constants kqp and kqd were obtained by linear fitting (Fig. 3a, b and Table S1). The performed simulations unambiguously showed that the quenching constants kqp and kqd, obtained by the numerical treatment described above, match the constants kqS and kqT that were used generate the decays (Fig. 3a, b and SI).
The key findings of our modelling, which can be translated to real experiments, are the following: (i) for any combination of rate constants and quencher concentration, the delayed lifetime of S1, easily measurable by time-resolved fluorescence, matches with the decay lifetime of T1; (ii) the rate constants obtained from Stern–Volmer experiments for prompt and delayed fluorescence lifetimes (kqp and kqd) are the quenching constants of the singlet (kqS) and the triplet states (kqT), respectively.
The physical meaning of these results is that when the rates of prompt and delayed fluorescence differ by three or more orders of magnitude, S1 state does not get populated to any significant extent during the emission of delayed fluorescence, but is visited transiently, merely providing a radiative deactivation channel for the T1 state. Thus, all chemistry involving S1 state occurs during the initial phase right after the irradiation pulse (prompt fluorescence regime). At later times, chemical reactions are carried out almost exclusively by T1 state, whose evolution is mirrored by delayed fluorescence. Thus, the difference between kqp and kqd for isophthalonitrile-based TADF molecules, which was reported earlier,20,21 reflects the difference between the PeT quenching rate constants for the respective S1 and T1 states (kqS and kqT). Considering the importance of TADF fluorophores for the field of photocatalysis, understanding of the physical origin of this difference is crucial for the development of new optimized catalytic systems.
Notably, Stern–Volmer plots constructed using integrated intensities, as opposed to lifetimes, are highly non-linear (Fig. 3c), which is expected for biexponential decays, whose integrals do not correlate with the underlying individual rate constants. Limburgh has proposed a formula to obtain kqS and kqT from the fitting of the non-linear trend of intensities.26 Nevertheless, this approach requires the evaluation of 3 parameters (kisckrisc, kqS and kqT) from a single set of experimental data, which leads to high uncertainty on the determined parameters. On the other hand, the Stern–Volmer analysis on emission lifetimes allows to obtain 2 parameters (kqS, kqT) from 2 distinct sets of data (τp and τd), ensuring higher precision in the evaluation of the quenching constants. Therefore, using emission lifetimes for determination of kqS and kqT should be preferred over using integrated intensities. Unfortunately, in the majority of the literature reports on TADF chromophores in photoredox catalysis, Stern–Volmer analyses have been carried out using linear fitting of luminescence intensities, and, therefore, the extracted quenching constants are inaccurate.
| *TADF + Q → TADF˙− + Q˙+ | (6) |
In the case of 4DPAIPN and 3DPAFIPN, the quenching constant for S1 state, kqS, extracted from the prompt fluorescence measurements, was found to be ca. 3 orders of magnitude higher than that for T1 state, kqT, obtained using delayed fluorescence. For 4CzIPN, kqS was 3 times higher than kqT (Table 1). To rationalise these results, we calculated the free energies (ΔG) for the PeT processes (6) involving S1 and T1 states using the ground-state redox potentials, measured in dichloromethane (DCM), and the spectroscopic energy E00 estimated from the fluorescence and phosphorescence spectra registered at 77 K (see Fig. S2–S5):
| ΔG = F[E(TADF/TADF˙−) − E(Q˙+/Q)] + E00(TADF/*TADF) | (7) |
It is worth noting that the analysed TADF chromophores in dichloromethane solution display a chemically and electrochemically reversible electron transfer process in the cathodic region (Fig. S6), while DIPEA is reported to undergo chemically irreversible electron transfer process in the anodic region and significantly different redox potentials are reported in the literature also in the same solvent.28–30 We assumed a value of E(DIPEA˙+/DIPEA) = 0.94 V vs. SCE, based on one of the literature reports,28 but uncertainty of this value is reflected in the calculated ΔG.
As expected, when the PeT driving force was positive, like in the case of case of T1 state of 4DPAIPN, the quenching constant (kqT) was small. However, the driving force is clearly not the only factor that defines the difference between the quenching rate constants for the studied PeT reactions. For example, the ΔG values for the processes involving S1 of 4DPAIPN and T1 of 3DPAFIPN, are very similar, while the respective quenching constants are different by 3 orders of magnitude.
To eliminate the uncertainty in the ΔG values stemming from the irreversibility of DIPEA oxidation, the measurements were carried out using ferrocene derivatives as electron donors (Fig. 4) that, in contrast to DIPEA, display a chemically reversible anodic electron transfer process (Fig. S7) with redox potentials such that the PeT processes involving either S1 or T1 states of all of the studied TADF sensitizers are exergonic. The absorption spectra of the photosensitizers and ferrocene-based quenchers correspond to the sum of the spectra of the isolated compounds (see Fig. S14), and both prompt and delayed fluorescence decays can be individually fitted with a monoexponential function (Fig. S15–S24). Therefore, there is no evidence of static quenching and the kinetic constants determined are solely attributed to dynamic quenching. The measured quenching constants were found to be in the range of 108–1010 M−1 s−1, and the quenching constants for S1 were consistently higher than for T1 (Table 2).
![]() | ||
| Fig. 4 Structures of (a) ferrocene derivatives (electron donors) and (b) naphthoquinone derivatives (electron acceptors) used as quenchers in PET reactions. | ||
| Photocatalyst | Quencher | E/V vs. SCE | State | ΔG/eV | k q/109 M−1 s−1 |
|---|---|---|---|---|---|
| 4DPAIPN | Fc | 0.39 | S1 | −0.48 | 9.4 |
| T1 | −0.32 | 2.0 | |||
| FcCHO | −0.16 | S1 | −0.18 | 8.5 | |
| T1 | −0.02 | 1.5 | |||
| DmFc | 0.69 | S1 | −1.03 | 12 | |
| T1 | −0.87 | 9.0 | |||
| 4CzIPN | Fc | 0.39 | S1 | −1.18 | 14 |
| T1 | −1.10 | 11 | |||
| FcCHO | −0.16 | S1 | −0.88 | 15 | |
| T1 | −0.79 | 15 | |||
| DmFc | 0.69 | S1 | −1.73 | 7.9 | |
| T1 | −1.65 | 2.2 | |||
| 3DPAFIPN | Fc | 0.39 | S1 | −0.77 | 3.5 |
| T1 | −0.51 | 0.81 | |||
| FcCHO | −0.16 | S1 | −0.46 | 10 | |
| T1 | −0.20 | 3.8 | |||
| DmFc | 0.69 | S1 | −1.31 | 4.5 | |
| T1 | −1.05 | 1.0 |
It should be noted that the same trend is also observed in DCM (Fig. S15 and S18), suggesting that the polarity of the solvent does not significantly affect the quenching dynamics.
As expected, the differences became smaller and eventually vanished (for example, see the data for 4CzIPN and FcCHO quencher) as the constants approached the diffusion limit, taken as an approximation as the self-diffusion of toluene (kd = 1.2 × 1010 M−1 s−1).31 The Rehm–Weller plot (Fig. S28), constructed assuming the same electronic and nuclear factors (νN, κel, λ, see SI for more details) for the PC/Q couples in Table 2 did not reveal any correlation between the quenching rates and the PeT driving forces (ΔG) in the series. Therefore, the difference between the quenching constants for S1 and T1 states of the same photocatalyst cannot be explained solely on the basis of ΔG of the respective reactions, i.e. S1 state being a stronger oxidant than T1 state.
It is important to mention that while the cage escape of the geminate radical pair TADF˙−…Q˙+ and radical-pair recombination are important for the outcome of the overall catalytic cycle,32,33 spin selectivity on these processes do not influence the rates of the PeT reactions, since they occur after the quenching process.
Next, we examined oxidative PeT reactions involving one of the TADF chromophores, 4CzIPN. Such reactions follow the equation:
| *TADF + Q → TADF˙+ + Q˙− | (8) |
Three naphthoquinone derivatives (Fig. 4) were selected as quenchers based on the established chemical reversibility of their cathodic electron transfer process and considering the exergonicity of processes (8) involving either S1 or T1 states of 4CzIPN.
Remarkably, unlike in the case of reductive quenching, the values of kqS and kqT (Table 3) for reactions with naphthoquinones were found to be very similar (<10% difference, Fig. 5 and Table 3).
| PC | Q | Ea/V vs. SCE | ΔGb/eV | k q/109 M−1 s−1 | |
|---|---|---|---|---|---|
| a Redox potentials of the quinone derivatives in dichloromethane.34 b ΔG is calculated by considering E(4CzIPN˙+/4CzIPN) = 1.52 V vs. SCE.27 | |||||
| 4CzIPN | 2,3-CINQ | −0.38 | S1 | −1.31 | 4.7 |
| T1 | −1.05 | 4.5 | |||
| NQ | −0.69 | S1 | −0.77 | 5.6 | |
| T1 | −0.51 | 5.0 | |||
| 2-MetNQ | −0.73 | S1 | −0.46 | 3.7 | |
| T1 | −0.20 | 3.5 | |||
![]() | ||
| Fig. 5 Rate constants of reductive and oxidative quenching processes involving S1 and T1 states of selected TADF chromophores measured in toluene solution (values reported in Tables 2 and 3). | ||
The key new analysis to explain the observations in this paper is the hole and electron distributions in the adiabatic S1 and T1 states presented in Fig. 6, with a visual decomposition of the molecule into “core” (orange) and “shell” (green) fragments. For all three TADF chromophores there is a highly relevant change between the location of the hole density in S1 and T1. The hole density is much larger in the core fragment for T1 (32–48%) with respect to S1 (4–16%) and, as a consequence, we expect a reduced rate to accept an electron in T1via reductive quenching because the core fragment is not in contact with the quencher. Conversely, the electron density is fairly similar between T1 and S1 across the three molecules consistent with a similar rate of oxidative quenching for singlet and triplet. It should be noted that this situation is relatively uncommon as simpler chromophores often have the lowest singlet and triplet excited states sharing a similar HOMO–LUMO configuration, with almost identical hole and electron density. On the other hand, in this case, different combinations of intramolecular charge transfer transitions are involved in S1 and T1. This observation is robust not only across the three molecules considered but it is also not sensitive to changes in the solvent model and is seen also in vertical rather than adiabatic transitions (see SI for further detail).
![]() | ||
| Fig. 6 Computed electron and hole density distributions in adiabatic S1 and T1 states for 4CzIPN, 4DPAIPN, and 3DPAFIPN. SUM refers to total density of hole or electron in the core fragment. | ||
The key findings of the present study are as follows: (i) the decay of the delayed fluorescence in TADF fluorophores mirrors the dynamics of T1 state, so that the kinetics of T1 state can be measured using fluorescence spectroscopy rather than transient absorption methods; (ii) Stern Volmer plots based on the prompt and delayed fluorescence lifetimes yield the quenching rate constants for the reactions of S1 and T1 states, respectively; (iii) for the TADF chromophores studied, reductive quenching involving S1 state is considerably faster compared to that involving T1 state, while the energy difference between S1 and T1 do not fully account for the observed difference; (iv) electronic structure calculations suggest that the difference between the quenching rates for S1 and T1 states is related to the difference between the hole–electron density distributions in these states. In S1 state the hole is more peripherally located, facilitating reactions in which the photocatalyst serves as an electron acceptor (reductive quenching).
Future studies will aim to determine whether the observed effects are present in other classes of TADF chromophores, such as Cu(I) complexes or Zn(II) porphyrins, to define the scope and limitations of the present findings. Furthermore, cage escape yields for the radical pair (e.g., TADF˙−…Q˙+) generated by quenching S1 or T1 excited states will be investigated. S1 is generally expected to have lower cage-escape yield based on spin considerations. The initial spin state of the radical pair matches that of the original excited state (S1 or T1). Fast spin-allowed singlet charge recombination is expected to effectively compete with cage escape,37,38 whereas the latter should outcompete the much slower spin-forbidden triplet charge recombination. Indeed, it has been recently demonstrated that the overall efficiency of the reaction depends critically on which excited state undergoes reductive quenching and, as expected, the triplet was found to be more productive.38 This observation highlights the importance of promoting efficient quenching from the triplet state while minimizing quenching from the singlet. From this perspective, the insights gained from our study are especially significant, as they clarify the underlying factors responsible for the relative difference in quenching rates for the singlet and triplet excited state. Overall, these findings integrate and complement existing strategies for the rational engineering of efficient TADF photocatalysts15,18 by opening up new perspectives, including the design of sensitizers with hole densities in the T1 state located at the periphery to facilitate reductive electron transfer processes, ultimately enhancing the efficiency and selectivity of the catalytic process.
Supplementary information is available. See DOI: https://doi.org/10.1039/d5sc04948b.
Footnote |
| † Authors contributed equally. |
| This journal is © The Royal Society of Chemistry 2025 |