Open Access Article
This Open Access Article is licensed under a
Creative Commons Attribution 3.0 Unported Licence

Identifying efficiency-loss pathways in triplet–triplet annihilation upconversion systems

Abhishek Kalpattu a, Daniel E. Falvey b and John T. Fourkas *bcd
aDepartment of Materials Science and Engineering, University of Maryland, College Park, MD 20742, USA
bDepartment of Chemistry and Biochemistry, University of Maryland, College Park, MD 20742, USA. E-mail: fourkas@umd.edu
cInstitute for Physical Science and Technology, University of Maryland, College Park, MD 20742, USA
dMaryland Quantum Materials Center, University of Maryland, College Park, MD 20817, USA

Received 13th February 2025 , Accepted 25th April 2025

First published on 2nd May 2025


Abstract

Triplet–triplet annihilation upconversion (TTA-UC) systems have been studied extensively recently, and have been proposed for use in a wide range of applications. Identification of the dominant mechanisms of upconversion-efficiency loss (UEL) will assist in the development of efficient TTA-UC systems. In this work, we combine experiments and kinetic analysis to study UEL. We identify exciplex formation and reverse triplet energy transfer (TET) as the two most important UEL mechanisms in the model TTA-UC system of platinum octaethylporphyrin (PtOEP) and 9,10-diphenylanthracene (DPA). Based on spectral analysis and time-resolved photoluminescence experiments, we show that exciplex formation is a potent UEL pathway in the PtOEP–DPA system. We demonstrate that prolonged sensitizer phosphorescence arises from reverse TET from annihilator triplet states, and that the reverse TET is likely facilitated by thermal population of low-frequency vibrational states in the sensitizer and the annihilator. Additionally, we demonstrate how the rate constants for reverse TET and exciplex formation can be estimated based on knowledge of a few key parameters and the experimental value of the optimum sensitizer concentration.


Introduction

Triplet–triplet annihilation upconversion (TTA-UC) is a process in which two annihilator triplet states that were created via triplet energy transfer (TET) from photogenerated sensitizer triplet states undergo disproportionation to produce an annihilator in a fluorescent singlet state and an annihilator in the ground electronic state. TTA-UC does not require high excitation powers to achieve frequency upconversion,1–3 and can even be driven by an incoherent light source, such as a light-emitting diode (LED).4–7 TTA-UC is being pursued for a variety of applications, such as bioimaging,3,8–10 solar-cell performance enhancement,11–15 and organic LEDs.16–19 A kinetic analysis of the TTA-UC process indicates that the intensity of the upconverted light depends quadratically on irradiance under conditions in which unimolecular annihilator triplet-decay processes are favored over the bimolecular TTA process.20–22 In the converse situation, the intensity of the upconverted light depends linearly on the irradiance. The irradiance at which the quadratic and linear regimes meet is commonly referred to as the threshold irradiance (Ith).23 A higher sensitizer concentration increases the concentration of photogenerated sensitizer triplet states, thereby decreasing Ith. A smaller Ith is generally desirable, and is indicative of a TTA-UC system that performs efficiently at low irradiances. However, increasing the concentration of sensitizers has been observed to be deleterious to both the upconverted fluorescence intensity and the upconversion quantum yield.24,25 Sensitizer TTA, which results in the net loss of one sensitizer triplet for every two that are photogenerated, is one suspected culprit in the decrease of the upconversion quantum yield with increasing sensitizer concentration.26,27 However, for ideal TTA-UC systems, sensitizer TTA is often outcompeted by TET to the annihilator, and so sensitizer TTA may not be the sole driver behind the upconversion efficiency loss (UEL) that arises at high sensitizer concentrations. Back energy transfer from the emissive annihilator singlet states to the sensitizer, via long-range energy-transfer mechanisms such as Förster resonance energy transfer (FRET), has also been identified as a UEL mechanism that has the potential to become more important at high sensitizer concentrations.28,29

Isokuortii et al. studied the effects of reverse TET, i.e., from the annihilator triplet state to the sensitizer, in cases in which the annihilator triplet energy was lower than, or slightly higher than, the sensitizer triplet energy.30 These authors found that for an endothermic system in which the triplet energy of the annihilator was ∼0.08 eV lower than the triplet energy of the sensitizer, the rate constant for reverse TET from the annihilator to the sensitizer was twenty times larger than that for forward TET to the annihilator. A key finding of Isokuortii et al. is that even for a system in which the difference between the sensitizer triplet energy and the annihilator triplet energy is several times larger than the thermal energy, reverse TET can be the dominant annihilator triplet decay pathway at low irradiances. As discussed below, we believe that low-frequency vibrational modes are a key consideration in this phenomenon.

Efforts to reduce efficiency losses in the TTA-UC process, for example by minimizing the energy difference between the singlet and triplet states of the sensitizer,31–33 or by minimizing the energy difference between the annihilator and sensitizer triplet states,34,35 may inadvertently lead to UEL due to increased rate of reverse intersystem crossing (RISC) in the sensitizer in the first case, and of reverse TET in the second case. There have been numerous studies on photochemical upconversion systems in which the effect of reverse TET is manifested in the form of long-lived sensitizer emission.36–38 Meroni et al. note, in a study of a TTA-UC system comprising Pd(II) meso-tetraphenyl tetrabenzoporphine and perylene,38 that the lifetimes of the annihilator triplets and of the delayed phosphorescence that arises due to reverse TET must be equal. These authors concluded that an increase in the concentration of annihilators suppresses the effects of reverse TET, and increases the lifetime of upconverted emission. Although Meroni et al. did not study the effects of an increase in the sensitizer concentration, presumably such an increase would induce the opposite effect, i.e., enhancing the effects of reverse TET and reducing the lifetime of the upconverted emission.

Gholizadeh et al. reported the sensitizer-concentration-dependent quenching of annihilator triplets, a process they termed “dynamic quenching.”24 These authors discovered that for a TTA-UC system consisting of platinum octaethylporphyrin (PtOEP) and diphenylanthracene (DPA), the dynamic quenching rate constant was 4.6 × 106 M−1 s−1, roughly four orders of magnitude larger than the rate of intrinsic annihilator decay. Spin–orbit coupling (SOC) induced in the annihilator by the heavy central atoms common in photosensitizers, which is known as the heavy-atom effect,39,40 is one possible origin of dynamic quenching. However, Gholizadeh et al. showed that the addition of a bromine-containing compound to a solution containing PtOEP and DPA did not affect the upconverted fluorescence intensity significantly. These authors also reported an increase in the empirically determined rate of dynamic quenching upon the replacement of the platinum in PtOEP with zinc, a lighter element. An important distinction between ZnOEP and PtOEP is that the triplet energy of the former (1.78 eV) is closer to that of the annihilator, DPA (1.77 eV), than is the triplet energy of the former.41 Therefore, reverse TET may also affect the apparent rate of dynamic quenching when ZnOEP is the sensitizer. The work of Gholizadeh et al. demonstrates that the heavy-atom effect cannot be the sole origin of dynamic quenching.

Although a large body of work has been dedicated to exploring UEL mechanisms on an individual basis, we are not aware of any concerted effort to analyze the potentially important UEL mechanisms with a mass-conserving kinetic model. A mass-conserving TTA-UC model accounts accurately for the steady-state concentration of ground-state molecules that can participate in various UEL processes.42 Therefore, it would be beneficial to use such a model to analyze UEL mechanisms that are driven by interactions between ground-state and excited molecules at high irradiances. Here we perform such an analysis of mechanisms whereby an increase in the concentration of sensitizers hinders the upconversion ability of a TTA-UC system. Using PtOEP and DPA as a benchmark TTA-UC system, we also performed steady-state and time-resolved TTA-UC measurements. Based on the results of these experiments, we explore a range of UEL pathways, and use our kinetic model to determine which of these UEL mechanisms are the most influential. We employ the ratio of the sensitizer phosphorescence intensity to the upconverted fluorescence intensity, a quantity that we term the emission ratio, to identify UEL mechanisms that reduce the upconverted fluorescence intensity at steady state, while either not affecting, or even enhancing, the amount of sensitizer phosphorescence. We also use our kinetic model to analyze the effects of reverse TET on the upconversion quantum yield, the sensitizer emission decay, and the annihilator triplet lifetime. We demonstrate how reverse TET can affect the rate of dynamic quenching, as determined from the dependence of the rate of annihilator triplet decay on the sensitizer concentration. We show how the balance between enhanced upconverted emission due to an increase in sensitizer triplet concentration on the one hand, and enhanced rates of reverse TET and dynamic quenching due to a high concentration of sensitizers on the other hand, is affected by the rate constants for the latter two processes. Finally, we demonstrate how the rates of reverse TET and dynamic quenching can be ascertained independently based on determination of the sensitizer concentration that leads to the maximum upconversion quantum yield and of the experimentally-determined overall rate of sensitizer-dependent annihilator triplet decay.

Methods

Experimental

PtOEP, DPA, 1,3-diphenylisobenzofuran (DPBF), and dimethyl sulfide (DMS) were purchased from Sigma Aldrich and used as received. TTA-UC samples were prepared by mixing stock solutions of PtOEP and DPA in toluene. DMS (0.5 M) was added to each sample to scavenge the singlet oxygen that forms during photoexcitation. For optical measurements, TTA-UC samples were prepared in a 1 cm × 1 cm quartz cuvette (FireflySci). To ensure effective removal of oxygen, the TTA-UC samples were sealed with rubber septa and then purged with N2 gas for 10 min. To ensure that the DMS had removed the majority of any remaining oxygen, we only recorded data after the steady-state upconverted fluorescence intensity had stabilized. The phosphorescence lifetime of DPA in our samples (vide infra) is well in line with literature values for freeze/pump/thawed samples.43,44 Taken together, these observations suggest that the oxygen removal in our samples was sufficient.

Upconversion data were collected at room temperature using a custom-designed setup that is shown schematically in Fig. S1 (ESI). A 532-nm, continuous-wave (CW) laser (Coherent, Verdi V12) served as the excitation source. For time-resolved emission experiments, TTA-UC samples were excited by a train of square pulses that was generated by an acousto-optic modulator. A function generator (WaveTek, Model 29) and a radio-frequency generator (IntraAction, Model ME) were used to control the repetition rate, amplitude, and duty cycle of the pulse train. Emission from the TTA-UC samples was projected onto a forward-biased Si photodiode (NewFocus, 2001-FS) after passing through a 532-nm notch filter and a bandpass filter designed to isolate either the phosphorescence from PtOEP or the upconverted fluorescence from DPA. The time-resolved emission decay was detected by the Si photodiode and collected by a digital oscilloscope (LeCroy 9400A).

For steady-state TTA-UC measurements, samples were excited by a spatially-filtered, 532-nm CW beam. The excitation power was controlled with a half-wave plate and a polarizing beam cube. The half-wave plate was mounted in a computer-controlled rotation stage (Aerotech, Ensemble Epaq). The excitation power was measured by directing a portion of the spatially-filtered excitation beam through a mechanical chopper (Maxxon DC Motors) and onto a calibrated Si photodiode (Thorlabs, DET36A2). Steady-state emission from the TTA-UC sample was collected and projected through a mechanical chopper (Maxxon DC Motors) onto a matched Si photodiode. Upconverted light was isolated with a ±20 nm bandpass filter with its passband centered at 435 nm. Two digital lock-in amplifiers (Stanford Research Systems, DSP-SR810) were used to amplify the photodiode signals and to reject noise. Steady-state fluorescence measurements at various excitation powers were acquired in an automated fashion using a LabView program that rotated the half-wave plate by a 0.5° increment following each 2-s fluorescence-measurement interval. An entire set of upconverted fluorescence data typically consisted of data points collected at 32 different irradiances over a period of less than 70 s. Where needed, additional measurement intervals were added and a longer signal averaging time was used, resulting in an extended data-acquisition period of 135 s.

For spectral measurements, emission from the TTA-UC sample was coupled into a fiber-optic cable that was connected to a photospectrometer (Ocean Optics USB2000). A 532-nm notch filter was used to reject scattered excitation light. Emission ratios were calculated using the ratio of the “area under the curve” (AUC) of the emission spectra between 600 nm and 700 nm and the AUC of the emission spectra between 380 nm and 530 nm. These calculations were performed with MATLAB. To obtain emission ratio measurements at different excitation powers, a LabView program was used to rotate the half-wave plate by a 0.5° increment following each 4-s spectral measurement interval.

Kinetic modeling

The time-dependent concentrations of sensitizers in the lowest triplet state ([3S*]), annihilators in the emissive excited singlet state ([1A*]), annihilators in the lowest triplet state ([3A*]), and annihilators in a higher triplet state ([3A**]) were computed numerically using MATLAB. A simulation covering dynamics over a period of 400 μs was performed with a fixed irradiation profile that consisted of a 200 μs, constant-amplitude, 2 W cm−2 pulse, followed by a 200 μs “dark” period. [3S*]t, [1A*]t, [3A*]t, and [3A**]t were computed with a time step of 1 ns. We note that the sensitizer singlet state was not modelled explicitly. Because ISC in most TTA-UC sensitizers, including PtOEP, takes place on a time scale that is far shorter than the time step of our simulations, the sensitizer triplet state was treated as being populated directly and losslessly via absorption of light from the ground state.

We now develop equations for [3S*]t, [1A*]t, [3A*]t, and [3A**]t. A detailed derivation of these equations is given in the ESI. For [3S*]t, we find

 
image file: d5cp00578g-t1.tif(1)
where dt is the simulation time step. [S]0 and [A]0 represent the initial concentrations of sensitizers and annihilators in the TTA-UC system, respectively, and kex, kFRET, kST, kph, and kSq are the rate constants for excitation of ground-state sensitizers, FRET, intrinsic sensitizer triplet decay, sensitizer phosphorescence decay, and sensitizer triplet self-quenching, respectively. Additionally, kSTTA is the rate constant for TTA between two sensitizers, and kFsens and kRsens are the rate constants for triplet energy transfer from the sensitizer to the annihilator and from the annihilator to the sensitizer, respectively. For [3A*]t we find
 
image file: d5cp00578g-t2.tif(2)
where kATTA is the rate constant for TTA between two annihilators, kAT is the rate constant for intrinsic annihilator triplet quenching, kIC is the rate constant for internal conversion from the Tn state to the T1 state of the annihilator, and kAq is the rate constant for dynamic quenching. For [1A*]t we find
 
image file: d5cp00578g-t3.tif(3)
where kfl and kANR are the rate constants for the radiative and non-radiative decay of annihilator singlets, respectively, and kRISC is the rate constant for reverse ISC from the Tn state to the excited singlet state of the annihilator. For [3A**]t we find
 
[3A**]t ≈ [3A**]t−dt + dt(0.75kATTA[3A*]t−dt2 − (kIC + kRISC)[3A**]t−dt).(4)
Finally, the upconverted fluorescence rate, Fl(t), and the sensitizer phosphorescence rate, Ph(t), were calculated using
 
Fl(t) = kfl[1A*]t(5)
and
 
Ph(t) = kph[3S*]t,(6)
respectively. The initial conditions [3S*]0, [3A*]0, [1A*]0, [3A**]0 = 0 were used when solving these equations. The smaller the time step dt, the more accurate the results for [3S*]t, [3A*]t, [1A*]t, [3A**]t. A time step of 1 ns was chosen as a compromise between accuracy and efficiency.

The effective sensitizer and annihilator triplet lifetimes were calculated from the slopes of logarithms of the phosphorescence and fluorescence decays immediately after the excitation ceased,

 
image file: d5cp00578g-t4.tif(7)
and
 
image file: d5cp00578g-t5.tif(8)
respectively. The factor of ½ in eqn (8) arises from the fact that upconverted fluorescence depends quadratically on the population of annihilator triplets.

To model the formation and decay of radiative excimers and exciplexes explicitly, the concentrations of two additional species were introduced. The time evolution of the concentration of sensitizer excimers ([S*S]t) is given as

 
[S*S]t ≈ [S*S]t−dt + dt(kSq([S]0 − [3S*]t−dt)[3S*]t−dtkexcimerph·[S*S]t−dt),(9)
where kexcimerph is the rate constant for the radiative decay of sensitizer excimers. No non-radiative pathway was considered for the decay of sensitizer excimers. The rate of sensitizer excimer phosphorescence is then given by
 
Phexcimer(t) = kexcimerph[S*S]t.(10)

The time evolution of the concentration of exciplexes ([A*S]t) is given as

 
[A*S]t ≈ [A*S]t−dt + dt(δ·kAq([S]0 − [3S*]t−dt)[3A*]t−dtkexph[A*S]t−dt),(11)
where δ is the fraction of dynamic quenching events that result in the formation of sensitizer–annihilator exciplexes and kexph is the rate constant for the radiative decay of sensitizer–annihilator exciplexes. No non-radiative pathway was considered for the decay of sensitizer–annihilator exciplexes.

The steady-state phosphorescence and upconverted fluorescence rates were calculated by starting with the kinetic equations

 
image file: d5cp00578g-t6.tif(12)
 
image file: d5cp00578g-t7.tif(13)
 
image file: d5cp00578g-t8.tif(14)
and
 
image file: d5cp00578g-t9.tif(15)
eqn (12)–(15) are solved under steady-state conditions to obtain the steady-state concentration of sensitizer triplets ([3S*]SS), annihilator triplets ([3A*]SS), annihilator singlets ([1A*]SS), and higher annihilator triplet states ([3A**]SS). The rates of phosphorescence and upconverted fluorescence at steady state are given by
 
FlSS = kfl[1A*]SS(16)
and
 
PhSS = kph[3S*]SS,(17)
respectively. The emission ratio at steady state is given by
 
image file: d5cp00578g-t10.tif(18)
The quantum yield of delayed phosphorescence from reverse TET (ΦRph) is given by
 
image file: d5cp00578g-t11.tif(19)
[3A*]t, [3S*]t, [S]t, and [A]t were calculated numerically. Simulated irradiation of the TTA-UC sample occurred during the first 0.2 ms. Therefore, the total number of photons consumed as a result of photoexcitation is given by
 
image file: d5cp00578g-t12.tif(20)

Density–functional–theory calculations

Density–functional–theory (DFT) calculations on DPA and DPBF in their lowest triplet states were performed using (u)B3LYP-d3/def2SVP, i.e., the unrestricted B3LYP functional using Grimme's empirical dispersion correction45 (d3) and the basis set of Weigend and Ahlrichs46 (def2SVP). The vibrational frequencies obtained from the calculations on DPA were used to calculate the average energy in each mode via the equipartition theorem, which dictates that, for each mode,
 
Evib〉 = ħω(1/(exp(ħω/kBT) − 1)).(21)

Because the goal of these calculations was to determine the order of magnitude of the thermal vibrational energy, we did not include an anharmonic correction to the mode frequencies (which would increase the thermal vibrational energy) or attempt to use a more accurate method for treating the triplet state.

Results and discussion

A summary of the kinetic processes that constitute our basic mass-conserving TTA-UC model is shown in Fig. 1. The basic model does not consider any UEL mechanisms other than the intrinsic decays of sensitizer and annihilator triplets and the non-radiative decay of annihilator singlets. The rate equations for [3S*]t, [3A*]t, [1A*]t, and [3A**]t were presented previously.42 Although our mass-conserving TTA-UC model provides valuable mechanistic insights into the TTA-UC process, the model does not predict that UEL is correlated with the sensitizer concentration.
image file: d5cp00578g-f1.tif
Fig. 1 A schematic illustrating the processes that are considered in the basic TTA-UC model.

Fig. 2a and b show that the experimental upconverted fluorescence intensity at any given value of the irradiance increases when the sensitizer concentration goes from 2.5 μM to 10 μM. However, at higher sensitizer concentrations, the upconverted fluorescence intensity at any given value of the irradiance decreases with increasing [S]0. A high concentration of sensitizers therefore appears to limit the upconversion yield, presumably through mechanisms such as quenching of annihilator triplets, expenditure of sensitizer triplets through sensitizer TTA, parasitic absorption of upconverted emission, and FRET.


image file: d5cp00578g-f2.tif
Fig. 2 Experimental TTA-UC data for PtOEP–DPA samples. (a) The dependence of the upconverted fluorescence intensity on irradiance. The data (squares) were fit with the basic mass-conserving TTA-UC model. (b) The dependence of the upconverted fluorescence intensity at specific excitation powers on the sensitizer concentration. (c) The dependence of the emission ratio on the irradiance. (d) Emission spectra of PtOEP-DPA samples normalized to the emission intensity at 445 nm. The peaks that originate from upconverted fluorescence are indicated by *, and the peaks that originate from sensitizer phosphorescence are indicated by #. The arrow indicates scattered excitation light and sample emission around 532 nm, which were blocked using a notch filter. The concentration of the annihilator was 1 mM.

The emission ratio, eqn (18), allows us to differentiate between loss mechanisms that affect the yields of both sensitizer triplets and annihilator singlets upon photoexcitation on the one hand, and loss mechanisms that only impact the yield of annihilator singlets on the other. In the basic TTA-UC model, the expression for the emission ratio can be expanded as

 
image file: d5cp00578g-t13.tif(22)
In the low-irradiance regime, [3S*]SS is small, so the dominant term within the radicand in eqn (20) is kAT. It is also the case when [3S*]SS is small that
 
image file: d5cp00578g-t14.tif(23)
Under these conditions eqn (22) can be simplified to
 
image file: d5cp00578g-t15.tif(24)

Therefore, at low irradiances, the emission ratio decreases as [3S*]SS increases. Because [3S*]SS is proportional to I, the emission ratio is inversely proportional to I at low irradiances. As I increases, [1A*]SS becomes proportional to I until, at high enough irradiance, [1A*]SS ceases to increase. When [1A*]SS and [3S*]SS are both proportional to I, the emission ratio is independent of I. This trend in decreasing emission ratios with increasing I is consistent with the experimental TTA-UC data in Fig. 2c, which show that emission ratios decrease monotonically with I for all samples. Moreover, the rate of this decrease, i.e., the slope in Fig. 2c, becomes smaller with increasing I until the emission ratios are almost independent of I. An increase in [S]0 will result in an increase in [3S*]SS, all other things being equal. Therefore, according to eqn (24), we can expect that samples with higher sensitizer concentrations should exhibit lower emission ratios. This expectation is not consistent with the data in Fig. 2c. Our experiments show that when [S]0 is within two orders of magnitude of [A]0, the emission ratio increases with increasing [S]0. We posit that this behavior is indicative of the existence of UEL mechanisms that quench excited annihilator molecules selectively while leaving the excited sensitizer population unaffected. Another possibility is that these UEL mechanisms enhance the yield of sensitizer triplets at the expense of excited annihilator molecules.

Fig. 2d shows emission spectra of TTA-UC samples with varying sensitizer concentrations. The spectra are normalized to the emission intensity at 445 nm. With increasing sensitizer concentration, the phosphorescence intensity increases relative to the intensity of the upconverted emission at 445 nm. Furthermore, an additional emission feature appears in the near-infrared (NIR) between 750 nm and 800 nm at higher sensitizer concentrations. This feature is not observed in the absence of the annihilator (Fig. S2a, ESI). Therefore, we can rule out sensitizer-excimer emission and aggregation-induced emission as the cause of this NIR emission band. Radiative decay of exciplexes is a more likely explanation for the appearance of this feature. Indeed, we find that the integrated intensity of the NIR emission band has a dynamic dependence on excitation power, which is akin to the relationship between the upconverted fluorescence intensity and the excitation power (Fig. S2b, ESI). Because the relationship between the integrated intensity of the NIR emission band and the excitation power is linear, rather than quadratic, at low excitation powers, we conclude that the concentration of the NIR-emissive species is directly proportional to the concentration of annihilator triplets. The data in Fig. S2a (ESI) also reveal the presence of a secondary sensitizer emission peak at 680 nm that is suppressed upon the addition of the annihilator. Upon close inspection (Fig. S2c, ESI), it is apparent that in the absence of the annihilator, the intensity of this secondary sensitizer-emission peak increases with increasing sensitizer concentration. We therefore believe that the sensitizer-concentration dependence of the emission peak at 680 nm is likely indicative of the formation of PtOEP triplet excimers, which have a broader, and slightly red-shifted, emission profile as compared to that of the monomer (at ∼650 nm).

Based on the above, we identify sensitizer-triplet excimer formation and quenching of annihilator triplets by ground-state sensitizer molecules to form exciplexes as two key UEL mechanisms in the PtOEP–DPA system. Additionally, because sensitizer phosphorescence is enhanced at the expense of upconverted fluorescence, we also consider the possibilities that FRET between the annihilator singlet state and the sensitizer singlet state and/or reverse TET are important UEL mechanisms. We additionally assess the effects of sensitizer TTA, a UEL mechanism that could be influential under circumstances in which the ratio between the concentrations of sensitizers and annihilators is high.

Fig. 3 shows a version of the Jablonski diagram from Fig. 1 in which the above-mentioned additional UEL mechanisms are included. Because the emissive annihilator singlet excited state is typically substantially higher in energy than the lowest sensitizer singlet excited state, FRET is assumed to be a unidirectional process from the annihilator to the sensitizer. Sensitizer TTA refers to the disproportionation of two sensitizer triplets, a process that ultimately results in the loss of one triplet, assuming a lossless ISC process should the annihilation generate an excited singlet state. Sensitizer self-quenching refers to the sensitizer-triplet excimer-formation process. The rate and quantum yield of the radiative decay of sensitizer triplet excimers may differ from those of an isolated sensitizer in the triplet state. Reverse TET refers to the Dexter energy-transfer process from the annihilator triplet to the sensitizer triplet. Finally, dynamic quenching is a catch-all term for processes in which annihilator triplets are quenched by sensitizer molecules, including the formation of an exciplex between an annihilator triplet and a sensitizer molecule and SOC induced by heavy-metal-based sensitizer compounds. The list of UEL mechanisms that we consider here is not exhaustive. For example, the reabsorption of upconverted light by the sensitizer may lead to outcoupling losses in TTA-UC systems,47 and a concomitant reduction in overall upconversion efficiency.7 Functionally, this UEL pathway is similar to FRET, in that a sensitizer singlet state is generated in conjunction with the loss of an annihilator singlet state. The total rate at which upconverted light is extracted from the TTA-UC system can be expressed as βextkfl[1A*]SS. βext is a branching ratio that determines the proportion of upconverted light that can be extracted from a TTA-UC system, and is given by image file: d5cp00578g-t16.tif where kext is the rate constant for the external extraction of upconverted photons from the TTA-UC system and image file: d5cp00578g-t17.tif is the rate constant for the reabsorption of upconverted light by the sensitizer. When image file: d5cp00578g-t18.tif is large and kext is small, βext becomes small and the reabsorption of upconverted light can emerge as a highly influential UEL mechanism. In reality, kext depends on a variety of factors, including the geometry of system, and therefore it is non-trivial to obtain a definitive estimate for βext. For those reasons, here we have opted to not study the impact of the reabsorption of upconverted light by the sensitizer on emission ratios or upconversion quantum yields. However, we do offer further insights into this parasitic reabsorption process in the ESI. Other potential UEL mechanisms, such as absorption of the upconverted light or of the excitation light by annihilators in the T1 state and annihilator triplet-excimer formation, are not considered here, because the magnitudes of these mechanisms do not depend on the concentration of sensitizers.


image file: d5cp00578g-f3.tif
Fig. 3 A schematic illustrating the processes that are considered in the enhanced TTA-UC model, with the UEL pathways depicted as red arrows.

Analytical, steady-state modelling was conducted to determine which of the aforementioned UEL mechanisms could replicate an upward trend in the emission ratio with increasing sensitizer concentration. The emission ratio at steady state was calculated from eqn (18). Fig. 4a shows that when kFRET, kSTTA, kSq, kRsens, and kAq are all set to 0, the emission ratio at irradiances below 100 mW cm−2 is inversely proportional to [S]0. A similar trend in the emission ratio is observed when kFRET = 2 × 109 M−1 s−1 and kSTTA, kSq, kRsens, and kAq are zero (Fig. 4b). It is apparent that if FRET is to be competitive with the radiative annihilator singlet-decay pathway, it must be the case that kFRET[S] ≈ kFRET[S]0 > kfl. Typically, [S]0 is orders of magnitude smaller than 1 M. Therefore, kFRET must be several orders of magnitude larger than kfl, which is usually on the order of 108 s−1, for FRET to become an influential loss pathway. An inverse relationship between the emission ratio and [S]0 persists at moderate-to-low irradiances in instances in which sensitizer TTA or sensitizer self-quenching are the only UEL loss pathways, barring intrinsic quenching effects (Fig. 4c and d, respectively). This result is not surprising, because both sensitizer TTA and sensitizer self-quenching decrease the upconversion efficiency by reducing the availability of sensitizer triplets.


image file: d5cp00578g-f4.tif
Fig. 4 Plots, based on kinetic modeling, of the dependence of the emission ratio on the irradiance for TTA-UC systems with UEL mechanisms including: (a) intrinsic quenching effects only; (b) FRET and intrinsic quenching effects only; (c) sensitizer TTA and intrinsic quenching effects only; (d) sensitizer self-quenching and intrinsic quenching effects only; (e) reverse TET and intrinsic quenching effects only; and (f) dynamic and intrinsic quenching effects only.

Reverse TET is an interesting case. Fig. 4e demonstrates that when kRsens = 2 × 109 M−1 s−1 and kSTTA, kSq, kFRET and kAq are 0, the emission ratio decreases with an increase in [S]0 from 10 μM to 100 μM. However, a further increase in [S]0 from 100 μM to 1 mM results in an increase in the emission ratio. There are two reasons why reverse TET is more likely to enhance the rate of sensitizer phosphorescence at the expense of upconverted fluorescence in TTA-UC systems than is FRET. First, reverse TET is a competitive pathway for the depletion of annihilator triplets, as long as kRsens[S] > kAT and kRsens[S] > kATTA[3A*]. The rate constant kAT for annihilators such as DPA has been shown to be as low as 200 s−1. Therefore, at low irradiances reverse TET can easily outcompete intrinsic annihilator triplet-quenching effects and TTA. At high irradiances, reverse TET can continue to outpace TTA, as long as the value of kRsens is on par with that of kATTA and [3A*] does not exceed [S]. Second, the upconverted fluorescence rate depends quadratically on [3A*]. Therefore, mechanisms that reduce the population of annihilator triplet states will have a stronger impact on emission ratios than do mechanisms that reduce the population of annihilator singlet states, such as FRET.

Finally, as seen in Fig. 4f, dynamic quenching on its own, when kAq = 2 × 109 M−1 s−1 and kSTTA, kSq, kFRET and kRsens are 0, can reverse the decreasing trend in emission ratios with increasing [S]0. An increase in [S]0 from 100 μM to 1 mM elicits a larger increase in the emission ratio than does an increase from 10 μM to 100 μM. In a similar vein to reverse TET, dynamic quenching becomes a dominant depletion mechanism for annihilator triplets when kAq[S] > kAT and kAq[S] > kATTA[3A*]. Unlike reverse TET, however, dynamic quenching does not generate an excited molecule that can participate in the upconversion process or can phosphoresce. Therefore, when compared to reverse TET, dynamic quenching has a stronger effect on suppressing the rate of upconverted fluorescence, but does not enhance sensitizer phosphorescence. Of the five UEL mechanisms discussed here, only reverse TET and dynamic quenching are consistent with the experimentally observed trends in emission ratios with increasing [S]0.

We first consider the behavior when only reverse TET is added to the basic model. In this situation, the effective annihilator-triplet lifetime, which is half the value of the slope of the log–log plots of normalized upconverted emission at time scales longer than ∼50 μs, decreases with increasing sensitizer concentration (Fig. 5a). On the other hand, under the influence of reverse TET alone, the phosphorescence lifetime is initially dominated by the rapid forward TET process, and only on time scales longer than ∼50 μs does the radiative decay of sensitizer triplets that have been repopulated via reverse TET appear as a second exponential component (Fig. 5b). The form of the slow decay component in Fig. 5b mimics the form of the upconverted-emission decays in Fig. 5a.


image file: d5cp00578g-f5.tif
Fig. 5 Transient decays derived from kinetic models including: (a) upconverted fluorescence in the presence of reverse TET only; (b) phosphorescence in the presence of reverse TET only; (c) upconverted fluorescence in the presence of dynamic quenching only; (d) phosphorescence in the presence of dynamic quenching only; (e) upconverted fluorescence in the presence of both reverse TET and dynamic quenching; and (f) phosphorescence in the presence of both reverse TET and dynamic quenching.

Fig. 5c and d show the results of kinetic modeling of the sensitizer phosphorescence and upconverted emission for systems in which dynamic quenching is the only UEL mechanism. Due to the dominance of dynamic quenching over TTA, no bimolecular decay component is visible in the upconverted-emission decays. However, the effective lifetime of annihilator triplets decreases rapidly with increasing sensitizer concentration. On the other hand, there is no change in the sensitizer-phosphorescence decays with increasing sensitizer concentration.

Fig. 5e demonstrates that when the effects of reverse TET and dynamic quenching are combined, the upconverted-emission decay at long times is dominated by the effects of dynamic quenching. The sensitizer-phosphorescence decays under the influence of both reverse TET and dynamic quenching resemble biexponential decays at all three sensitizer concentrations studied here (Fig. 5f). The form of the slow decay component in this case mimics that of the upconverted-emission decay curves in Fig. 5e.

The effects of reverse TET and dynamic quenching on the effective lifetime of the annihilator are most prominent when kAT < kRsens[S] < kFsens[A] or kAT < kAq[S] < kFsens[A]. Fig. S3a (ESI) shows the dependence of the effective annihilator triplet lifetime on log(kRsens). This function has an inverse sigmoid shape. With increasing kRsens, the repopulation of annihilator triplets via forward TET becomes the rate-limiting step in the quenching of annihilator triplets via reverse TET, as shown in Fig. S3b (ESI). Therefore, the upper asymptote in Fig. S3a and b (ESI) is limited by the rate of intrinsic annihilator-triplet decay, and the lower asymptote is limited by the rate of forward TET. A similar relationship exists between both the effective sensitizer and annihilator triplet lifetimes and log([S]0) (Fig. S3c and d, ESI). Note that because the slow decay of the phosphorescence is a direct result of reverse TET, the effective lifetime of the sensitizer triplets is identical to that of annihilator triplets (i.e., τA = τS). Experimentally, a straightforward method of establishing the influence of reverse TET and dynamic quenching is to determine the dependence of τA−1 on [S]0. Fig. S4a (ESI) illustrates the dependence of τA−1 on [S]0 under different UEL mechanisms. When dynamic quenching is the sole UEL mechanism, τA−1 is linearly dependent on [S]0, and the rate constant for dynamic quenching can be determined from the slope of the curve. However, with the inclusion of reverse TET, τA−1 tapers off with increasing [S]0. Fig. S4b (ESI) shows that when reverse TET is the only UEL mechanism, this tapering-off effect of τA−1 is more prominent for smaller values of [A]0. With increasing [A]0/[S]0, the rate of forward TET is enhanced, thus decreasing τA−1 and diminishing the overall impact of reverse TET, in agreement with the experimental results of Meroni et al.38 Meroni et al. also found that the improvement in the overall efficiency of triplet sensitization with an increased annihilator concentration reaches a plateau. However, our results indicate that the critical parameter is not [S]0, but rather [A]0/[S]0. Indeed, for a given set of rate parameters, we find that a near universal relationship exists between τA−1 and log([A]0/[S]0), indicating that the effects of reverse TET on τA−1 are negated once [A]0/[S]0 exceeds a certain value. Fig. S5a and b (ESI) show that in a TTA-UC system that possesses identical kRsens and kFsens of 2 × 109 M−1 s−1, the effects of reverse TET on τS−1 and τA−1 are largely negated when the annihilator concentration is roughly two orders of magnitude larger than the sensitizer concentration.

We next compare the dependence of τA on log([S]0) when reverse TET is the only UEL mechanism and when both reverse TET and dynamic quenching are present (Fig. S6a and b, respectively, ESI). These figures also show the dependence of ΦRph, the quantum yield for delayed phosphorescence from reverse TET measured over a 400-μs period, on log([S]0). ΦRph is the difference in the quantum yield for phosphorescence (i.e., the number of emissive sensitizer triplets divided by the number of absorbed photons) with and without reverse TET. Fig. S6a (ESI) shows that the growth in ΦRph with increasing [S]0 begins to accelerate when [S]0 > 10−6 M. The growth in ΦRph becomes linear with respect to log([S]0) when log([S]0) is ∼10−5 M. With a further increase in [S]0, ΦRph reaches an asymptote. On the other hand, the decrease in τA with increasing [S]0 accelerates when kRsens[S]0kAT, a region in which the growth in ΦRph with log([S]0) has already begun to saturate. Fig. S6b (ESI) demonstrates that dynamic quenching causes τA to begin to decrease at approximately the value of [S]0 at which ΦRph begins to increase, and that the behaviors of these two quantities roughly mirror one another. Once again, we find that dynamic quenching has a stronger effect on τA than does reverse TET. Moreover, dynamic quenching noticeably reduces ΦRph. When kRsens = 2 × 109 M−1 s−1, kAq = 0 M−1 s−1, and [S]0 = 1 mM, ΦRph is almost 90%, indicating that the vast majority of sensitizer phosphorescence originates from resensitized triplet states. The total phosphorescence quantum yield in this case is ∼19.0%. On the other hand, when kRsens and [S]0 are maintained at their original values but kAq is increased to 2 × 109 M−1 s−1, ΦRph drops to just below 50%, indicating that a smaller proportion of the sensitizer phosphorescence can be attributed to reverse TET in this case.

Now that we have established that an increase in the sensitizer concentration can lead to an increase in τA, and thereby limit the upconverted fluorescence yield and the upconversion quantum efficiency, we turn to the dependence of the upconverted quantum yield, ΦUC, on the sensitizer concentration. ΦUC is calculated under the steady-state approximation for a specific value of I, and is defined as

 
image file: d5cp00578g-t19.tif(25)
In Fig. 6a–c, we demonstrate that ΦUC possesses a Gaussian-like distribution as a function of log([S]0). When [S]0 is small, ΦUC increases with increasing [S]0, because a larger value of [S]0 leads to larger values of [3S*]SS and [3A*]SS. A higher value of [3A*]SS implies that the bimolecular TTA process is favored over all other annihilator-triplet-depletion mechanisms, thereby increasing ΦUC. However, with further increases in [S]0, the effects of an increase in [3A*]SS are outweighed by an increase in the rates at which annihilator triplets undergo dynamic quenching and reverse TET. For pure dynamic quenching, an increase in [S]0 is detrimental to ΦUC when kAq[S]0 > kAT.


image file: d5cp00578g-f6.tif
Fig. 6 Dependence of ΦUC on log([S]0) for TTA-UC systems with UEL mechanisms including: (a) dynamic quenching effects only, varying kAq; (b) reverse TET only, varying kRsens; and (c) both dynamic quenching and reverse TET, varying kAq. (d) [S]0,ideal for TTA-UC systems with both dynamic quenching and reverse TET, varying kAq and kRsens. Ascending numeric labels for kAq from 1 through 7 indicate increasing values, by one order of magnitude each, from 2 × 103 M−1 s−1 to 2 × 109 M−1 s−1, and the ascending numeric labels for kRsens from 1 through 4 indicate increasing values, by one order of magnitude each, from 2 × 104 M−1 s−1 to 2 × 107 M−1 s−1. The dashed lines in (a)–(c) represent the concentrations at which the maximum value for ΦUC is found.

When a TTA-UC system possesses a large enough pool of ground-state annihilator molecules that annihilator mass conservation can be ignored, d[3A*]/dt can be expressed as

 
image file: d5cp00578g-t20.tif(26)
When the unimolecular decay terms are dominant, we can approximate that
 
image file: d5cp00578g-t21.tif(27)
At steady state, [3A*]SS can then be expressed as
 
image file: d5cp00578g-t22.tif(28)
and the rate of fluorescence becomes
 
image file: d5cp00578g-t23.tif(29)
When the bimolecular decay term is dominant, we can approximate that
 
image file: d5cp00578g-t24.tif(30)
[3A*]SS can then be expressed as
 
image file: d5cp00578g-t25.tif(31)
and the rate of fluorescence is
 
image file: d5cp00578g-t26.tif(32)
Setting eqn (29) equal to eqn (32), we find that
 
image file: d5cp00578g-t27.tif(33)
ΦUC for any TTA-UC system is higher when the bimolecular TTA process dominates over unimolecular annihilator triplet-decay processes. Therefore, a system with a lower Ith exhibits a higher ΦUC. According to eqn (33), the critical condition for an increase in [S]0 to be detrimental to ΦUC is that image file: d5cp00578g-t28.tif. The value of [S]0 that leads to the largest ΦUC, for a TTA-UC system that is affected by dynamic quenching, which we denote [S]0,ideal, is
 
image file: d5cp00578g-t29.tif(34)
The ideal concentration of ground-state sensitizers can therefore be estimated from the value of kAT/kAq. Fig. 6c shows that with increasing kAq, the peak of a plot of ΦUC as a function of [S]0 (dashed lines) shifts towards smaller values of [S]0. The magnitude of the peak value of ΦUC also decreases with increasing kAq.

In the case of reverse TET, determining the value of [S]0,ideal becomes more complicated. Once again ignoring the effects of annihilator mass conservation, d[3A*]/dt can be expressed as

 
image file: d5cp00578g-t30.tif(35)
Under the steady-state approximation, eqn (35) yields a quadratic expression for [3A*] in the low-irradiance regime (in which the unimolecular decay term dominates), and a cubic expression for [3A*] in the high-irradiance regime (in which the bimolecular decay term dominates). Although Ith can be obtained analytically from eqn (35), no analytical solution exists for [S]0,ideal. Fig. 6b demonstrates that when kRsens is below 2 × 109 M−1 s−1 and [A]0 = 1 mM, ΦUC continues to increase with increasing [S]0, even as [A]0/[S]0 reaches unity. When kRsens = 2 × 109 M−1 s−1 and [A]0 remains unchanged, a peak value of ΦUC is reached when [S]0 ∼ 3.1 × 103 M. When kRsens is increased further to 2 × 109 M−1 s−1 and [A]0 remains unchanged, [S]0,ideal is near 3.2 × 10−4 M. Fig. 6c and d illustrate the combination of the effects of reverse TET and dynamic quenching on [S]0,ideal. In Fig. 6c, kRsens was fixed at 2 × 109 M−1 s−1 and kAq was varied from 2 × 107 M−1 s−1 to 2 × 1010 M−1 s−1. Dynamic quenching has its strongest effect on [S]0,ideal when kRq is on the same order of magnitude as kAq.

The intrinsic annihilator triplet-quenching rate has a prominent effect on ΦUC and [S]0,ideal. In the case of dynamic quenching, [S]0,ideal is given roughly by kAT/kAq. Fig. S7a (ESI) shows that when kAT is set at 2 × 102 s−1 and kAq is fixed at 2 × 109 M−1 s−1, [S]0,ideal is precisely 10−7 M. When kAT is increased by an order of magnitude, [S]0,ideal becomes 10−6 M. In contrast, with increasing kAT, ΦUC decreases even as [S]0,ideal increases. In the case of reverse TET, [S]0,ideal increases proportionally with kAT, ΦUC increases with increasing kAT, and there is a change in the shape of the distribution of ΦUC. When kAT is small, the distribution in ΦUC is broader than when kAT is large. Fig. S7b (ESI) shows that when kAT is set at 2 × 102 s−1 and kRsens is fixed at 2 × 109 M−1 s−1, there is a less than a 1% drop in ΦUC as [S]0 increases by one order of magnitude from [S]0,ideal. On the other hand, when kAT is set at 2 × 103 s−1, there is a nearly 2% drop in ΦUC as [S]0 increases by an order of magnitude from [S]0,ideal. The implication of this finding is that for systems that possess smaller values of kAT, the impact on ΦUC due to deviations from [S]0,ideal is greatly reduced compared to what happens in systems with larger values of kAT. ΦUC increases with increasing irradiance. However, [S]0,ideal is largely independent of I. Fig. S7c (ESI) demonstrates that even under varying conditions, [S]0,ideal is independent of I.

Spin statistics dictate that for every excited singlet state that is populated via TTA, a higher triplet state and a quintuplet state should be populated as well. Typically, the quintuplet states are energetically inaccessible to the TTA pair. However, reverse ISC from a higher triplet state to the excited singlet state is possible. The effects of RISC can be captured by a branching ratio, βRISC, which is given by

 
image file: d5cp00578g-t31.tif(36)

Here, the rate constant kRISC governs the rate at which RISC occurs between the higher triplet state and the emissive excited singlet state.

The maximum theoretical value for ΦUC, ΦUC,max, is dependent on βRISCvia

 
image file: d5cp00578g-t32.tif(37)
We investigated whether delayed upconverted emission via reverse ISC from higher triplet states could alter our analysis of [S]0,ideal. Because RISC between higher triplet states and excited singlet states in the annihilator affects only the total yield of singlets following a TTA event, βRISC would not affect the impact of reverse TET or dynamic quenching on τA, FSS, or ΦUC. Indeed, in our analysis of the effects of reverse TET and dynamic quenching, we find that βRISC affects the magnitude of ΦUC, but does not affect [S]0,ideal (Fig. S7d, ESI).

To assess the impacts of dynamic quenching and reverse TET on the model PtOEP–DPA upconversion system, we conducted transient photoluminescence experiments. The upconverted fluorescence decay of the PtOEP/DPA system in deaerated toluene is shown in Fig. 7a. Annihilator-singlet decay occurs on a time scale of several ns, and does not limit the rate of the fluorescence decay. Therefore, the data in Fig. 7a must reflect the normalized logarithmic decay of [3A*]2. A function of the form

 
image file: d5cp00578g-t33.tif(38)
was used to fit the fluorescence decays. In eqn (38), α is a proportionality constant and β is a constant that describes the fraction of annihilator triplets that decay through bimolecular channels. Eqn (38) was adapted from Edhborg et al., whose β′ is identical to β.48 However, to obtain the best fit to our data, β′ must be slightly larger than β. We believe that the need for using a different value in the denominator arises from the fact that in the approximation used by Edhborg et al., no new annihilator triplet states are generated after irradiation ceases, i.e. their kinetic equation contains only loss terms for annihilator triplet states. However, there remain excited sensitizers when irradiation ceases that can generate additional annihilator triplet states via TET. The full kinetic expression does not have an analytical solution, but we find that the partially ad hoceqn (38) provides better fits to our data than does the expression of Edhborg et al.


image file: d5cp00578g-f7.tif
Fig. 7 Transient photoluminescence data. (a) Normalized upconverted fluorescence decay for TTA-UC samples with varying sensitizer concentrations. The data were fit to eqn (38). (b) Values of τA−1 extracted from the data in (a). The solid line is the best linear fit to the data. (c) Normalized phosphorescence decays fit to eqn (39). (d) Values of τS−1 extracted from the data in (c). The solid line is the best linear fit to the data. (e) Normalized phosphorescence decays fit to eqn (40). (f) Values of τA−1 extracted from (a) and (e). The concentration of the annihilator was 1 mM.

The most important fitting parameter in eqn (38) is τA. Unsurprisingly, we find the bimolecular component of the upconverted-fluorescence decay to be dominant when the initial concentration of the sensitizer is high. More interestingly, we find that an increase in the initial concentration of PtOEP from 2.5 μM to 273 μM is accompanied by a steady decrease in τA−1. When the experimentally-extracted values of τA−1 are plotted as a function of [S]0 (Fig. 7b), a linear relationship between τA−1 and [S]0 is apparent. The slope of the best linear fit in Fig. 7b is 3.6 × 107 M−1 s−1, and the y-intercept is 312 s−1. From the slope of the linear fit in Fig. 7b, we can estimate that kRsens + kAq = 3.6 × 107 M−1 s−1, and from the y-intercept of the best linear fit in Fig. 7b, we estimate that kAT = 312 s−1. The experimentally determined value of kAT is likely an underestimate, given that τA−1 in the limit [S]0 → 0 is smaller than kAT according to our simulations. Our result for kRsens + kAq is roughly an order of magnitude larger than the rate constant for dynamic quenching that was obtained by the Schmidt group.24

To determine whether reverse TET takes place in the PtOEP–DPA system, we measured the sensitizer-phosphorescence decay of the samples for which the upconversion data are shown in Fig. 7a. It is evident that the phosphorescence decays in Fig. 7c are biphasic. There is a rapid decay of sensitizer phosphorescence that we ascribe to forward TET, as well as a slower, non-exponential decay component. The decay curves in Fig. 7a and c bear a strong degree of similarity. We fit the phosphorescence decays to

 
image file: d5cp00578g-t34.tif(39)
where k1 is a rate constant corresponding to sensitizer triplet decay, which arises primarily from forward TET. The inverse lifetime of the long-lived component of phosphorescence decay, which we denote τP−1, depends linearly on [S]0, in analogy with the relationship between τA−1 and [S]0. However, the values of τA−1 extracted from the experimental data are nearly twice as large as the corresponding values of τP−1. When the phosphorescence-decay curves are fit instead with an equation of the form
 
image file: d5cp00578g-t35.tif(40)
the experimentally extracted values of τP−1 are substantially closer to those of τA−1 (Fig. 7e and f). Therefore, the slow component in the phosphorescence decays must represent a static portion of the population of the annihilator triplet states, rather than the population of the annihilator singlet states. Accordingly, FRET can be ruled out as the origin of the slow decay component in the sensitizer-phosphorescence decays. Hence, reverse TET is the most probable source of this slow decay, because the rate at which reverse TET occurs must depend linearly on the population of annihilator triplets. Furthermore, the rate of forward TET is typically expected to be higher than the rate at which the sensitizer triplet manifold is repopulated via reverse TET. Therefore, the rate-limiting step in prolonged sensitizer emission should be the lifetime of the annihilator triplets, in agreement with our experimental observations.

It is natural to question why reverse TET from DPA to PtOEP might occur in the first place, given that it is generally accepted that an energy difference of ∼0.14 eV exists between the energies of the triplet states of the two molecules (∼1.91 eV and ∼1.77 eV for PtOEP and DPA, respectively).41,49–51 The triplet energy of PtOEP is commonly determined from a phosphorescence peak at ∼650 nm.18,51–53 The triplet energy of DPA that is cited most often was determined from direct singlet–triplet absorption experiments in chloroform.49 Given the uncertainty inherent in these measurements, it is certainly possible that the energy difference between the triplet states of PtOEP and DPA is smaller than the commonly accepted value.

The use of pure electronic energies, as opposed to vibronic energies, to determine the likelihood of triplet-energy transfer can be misleading. For instance, it has been shown that triplet sensitization with a large energy gap can lead to vibrational excitation in the triplet state of the acceptor molecule.54 Detailed balance implies that the opposite process can also occur, i.e., a vibrationally excited “acceptor” molecule may be more likely to undergo reverse TET to a “donor” molecule. To assess this possibility, we performed DFT calculations on the lowest triplet state of DPA to determine the frequencies of its vibrational modes. Accordingly, we can estimate the average amount of vibrational energy at room temperature. Given the relatively large size of DPA, it is not surprising that this average energy comes out to be 0.46 eV, which is more than three times the gap between the triplet energies of PtOEP and DPA. PtOEP is also a large molecule that has floppy pendant groups, and so should have an even larger average thermal vibrational energy than does DPA. This factor also works in favor of reverse TET, in that vibrations that are populated in the electronic ground state of PtOEP need not be populated following reverse TET, effectively lowering the energy gap between the singlet and triplet states in this molecule. The rate of Dexter transfer between two species is related to the integral of the product of the normalized singlet–triplet absorption spectrum of the acceptor (PtOEP for reverse TET in our case) and the phosphorescence spectrum of the donor (DPA for reverse TET in our case),55 so we can think of the thermally-populated vibrations as creating hot bands in these spectra. Therefore, only thermally-populated, Franck–Condon-active modes of the donor and acceptor will help to promote reverse TET. Although it is not possible to make a quantitative estimate of the rate of reverse TET based on these arguments, it is clear that even if a small fraction of the thermally-populated vibrations in PtOEP and DPA are involved in a reverse TET transition, the 0.14 eV “barrier” can easily be overcome.

Changes in the conformation of the sensitizer or the annihilator that result from thermal fluctuations can also alter the driving force for TET substantially.56,57 For TTA-UC systems in which the triplet energy gap between the sensitizer and annihilator is small, endothermic reverse TET can be facilitated by small conformational changes that decrease the S0–T1 energy gap. Computational work by Zapata et al. showed that in-plane and out-of-plane distortions in the geometry of porphyrin molecules can alter the S0–T1 gap by as much as 0.15 eV at room temperature.58 Moreover, Gray et al. showed that at room temperature, the S0–T1 gap of DPA can fluctuate by ±1kBT as the orientations of the molecules’ two phenyl rings change.41 In our case, such conformational changes may induce small changes in the S0–T1 gaps of DPA and PtOEP molecules, thereby reducing the energy barrier to reverse TET by several kBT and subsequently increasing the rate constant for reverse TET by an order of magnitude or more.

Finally, it is important to understand that even when the rate constant for reverse TET is several orders of magnitude smaller than the diffusion-limited value of ∼109 M−1 s−1, ΦRph may be large enough for the effects of reverse TET to be observed in the form of a slow decay component in sensitizer phosphorescence-decay curves. Critically, ΦRph values increase exponentially as the magnitude of the product kRsens[S] begins to exceed the rate constant for intrinsic annihilator triplet decay. We highlight this behavior in Fig. S6a (ESI), in which ΦRph remains close to 0 when kRsens[S] ≪ kTA, but rises exponentially with increases in [S]. As [S] rises to 10−6 M, ΦRph grows to a value of 0.1, indicating that more than 10% of total sensitizer phosphorescence now arises from the decay of triplet states that are populated via reverse TET. In Fig. S6 (ESI), kRsens is ∼109 M−1 s−1, and kTA is ∼104 s−1. In a more experimentally relevant scenario, kTA values may be as small as ∼102 s−1. Therefore, ΦRph values well in excess of 10% can easily be achieved even when reverse TET is energetically unfavorable and kRsens is ∼107 M−1 s−1. When the effects of dynamic quenching are factored in, ΦRph is smaller, yet remains above 10% for sensitizer concentrations in excess of 10−6 M (Fig. S6b, ESI). In our case, we find that even if reverse TET only accounts for a fraction of the value of kRsens + kAq (3.6 × 107 M−1 s−1) that we obtained from transient photoluminescence studies, ΦRph should be large enough to lead to the observation of delayed sensitizer phosphorescence.

We also consider the possibility that the slow component in the sensitizer phosphorescence-decay curves arises from exciplex emission. The rate of formation of sensitizer–annihilator exciplexes, a process that falls under the umbrella of dynamic quenching, is linearly dependent on the concentration of annihilator triplets. In a scenario in which the overall rate constant for the radiative and nonradiative decay of the sensitizer–annihilator exciplexes exceeds the rate at which the exciplexes are formed, the combined decay of sensitizer phosphorescence and exciplex emission would be similar to the decay we observe experimentally. In Fig. S8a (ESI), decays for the sum of sensitizer phosphorescence and exciplex emission were calculated for a variety of rate constants for the radiative decay of sensitizer–annihilator exciplexes (kexph). Reverse TET was not included in the model, and the emission quantum yield of sensitizer–annihilator exciplexes (Φex) was assumed to be unity. When kexph is 200 s−1, a strong accumulation in the population of exciplexes is observed at long times, as evidenced by the growth in total emission even 0.2 ms after the excitation ceases. Only when the value of kexph is sufficiently large do the total emission decays resemble the phosphorescence decays observed experimentally.

A second experimentally relevant scenario is that Φex < 1, but the sum of the radiative and nonradiative sensitizer–annihilator exciplex decay processes is sufficiently large to avoid an accumulation of exciplexes on the sub-ms time scale. However, a small Φex would decrease the contribution of exciplex emission. Exciplex emission from DPA–PtOEP samples, as characterized by emission peaks in the 750 nm to 800 nm range, was suppressed completely by the addition of DPBF (Fig. S8b, ESI). DFT calculations revealed that the energy of the T1 state of DPBF lies in the vicinity of 1.43 eV, which is well below the energy of DPA's T1 state. Therefore, it is not surprising that the addition of DPBF to a sample containing DPA and PtOEP rapidly quenches the triplet states of DPA, thereby suppressing exciplex formation and subsequent emission. However, a slow bimolecular component in the sensitizer–phosphorescence decay is still observed for a TTA-UC sample containing DPBF, DPA, and PtOEP, despite the complete absence of exciplex emission from this sample (Fig. S8c, ESI). For reference, no slow component could be identified in the sensitizer-phosphorescence decay for a TTA-UC sample containing only DPBF and PtOEP. Therefore, from Fig. S8b and c (ESI) it is evident that the slow component that we observe in the phosphorescence decays in Fig. 7 cannot be attributed solely to exciplex emission, and must instead arise in large part from reverse TET.

We can conceive of only one scenario involving excimers in which the radiative decay of sensitizer triplet excimers results in the nonexponential decay of sensitizer emission. In the presence of the annihilator, the rate of forward TET would need to exceed the rate of sensitizer-excimer formation, and of all other sensitizer triplet-depletion mechanisms. Therefore, a biexponential decay of sensitizer emission, without any bimolecular decay features, would result when the rate of excimer decay, including all radiative and nonradiative pathways, is less than the rate of forward TET (Fig. S8d, ESI). However, in this scenario, the lifetime of the delayed excimer emission would be independent of the concentration of sensitizer molecules. Therefore, we conclude that sensitizer excimer emission cannot be the origin of the slow component observed in the phosphorescence decays.

Having determined the sum of the rate constants for the sensitizer-dependent quenching processes from the slope of the linear fit to Fig. 7b, we now turn to determining the rate constants for reverse TET and dynamic quenching. We demonstrate that with knowledge of [S]0,ideal and kRsens + kAq, we can estimate kRsens and kAq. Fig. 8a shows the expected value of ΦUC as a function of [S]0 for two cases, one in which kAq is equal to the slope of the linear fit to Fig. 7b (3.6 × 107 M−1 s−1) and kRsens is zero, and another in which kAq is zero and kRsens is equal to the slope of the linear fit to Fig. 7b (3.6 × 107 M−1 s−1). In the first case, we find that [S]0,ideal ∼ 10−5 M, which is close to the lower limit of [S]0,ideal that is given by kAT/kAq = 8.64 × 10−6 M. In the second case, [S]0,ideal is roughly 3.1 × 10−4 M. Therefore, for the given set of system parameters, including the experimentally obtained values of kAT and kAq + kRsens, we find that 8.64 × 10−6 M < [S]0,ideal < 3.1 × 10−4 M. Fig. 8b shows [S]0,ideal for different combinations of kRsens and kAq, with the constraint that kAq + kRsens = 3.6 × 107 M−1 s−1. Experimentally, [S]0,ideal can be estimated from the dependence of ΦUC on [S]0. Here, ΦUC can be obtained experimentally at any value of the irradiance, as long as the irradiance at which ΦUC is measured is unchanged for all samples with different [S]0. As an example, if [S]0,ideal were to be estimated as 1 × 10−5 M, kAq ∼ 2.5 × 107 M−1 s−1, which would mean that kRsens ∼ 1.1 × 107 M−1 s−1. With kRsens = 1.1 × 107 M−1 s−1, kAq = 2.5 × 107 M−1 s−1 and kAT = 312 s−1, we calculate ΦRph to be 0.67% when [S]0 = 10−5 M. When [S]0 = 10−4 M, ΦRph increases to 1.92%. Finally, when [S]0 is increased to be on par with [A]0, 1 mM, ΦRph increases to 5.15%.


image file: d5cp00578g-f8.tif
Fig. 8 Determining the rate constants for reverse TET and dynamic quenching. (a) Dependence of ΦUC on log([S]0) for TTA-UC systems with different values of kRsens and kAq. (b) Dependence of [S]0,ideal on kAq; kAq + kRsens was constrained to be 3.6 × 107 M−1 s−1 and kAT was set to the experimentally-determined value of 312 s−1.

Conclusions

We have used new experimental approaches for the rapid collection of steady-state and time-resolved TTA-UC data to explore several possible mechanisms whereby an increased concentration of sensitizers might lead to UEL. The observation of an enhancement in sensitizer phosphorescence with increasing concentration of sensitizers, at the expense of upconverted fluorescence, led to the identification of reverse TET and dynamic quenching as the two most influential UEL mechanisms.

Spectral measurements on PtOEP–DPA samples revealed the emergence of a new NIR emission peak with increasing concentration of sensitizers; this peak was absent when no annihilator was present. Unlike the intensity of the principal phosphorescence peak centered at ∼650 nm, the intensity of the NIR emission peak is not proportional to the excitation power. Therefore, we believe that the NIR emission results from exciplex formation between annihilator triplets and sensitizers. These exciplexes constitute a potential source of dynamic quenching.

Transient photoluminescence measurements of PtOEP–DPA samples in toluene indicated that the effective lifetime of annihilator triplets decreases with increasing concentration of sensitizers, with an overall sensitizer-dependent-quenching rate constant of 3.6 × 107 M−1 s−1. Reverse TET was identified as the probable cause of the majority of the prolonged sensitizer emission, which had an effective lifetime on par with that of the annihilator triplets. Although at first glance reverse TET might appear unlikely to be an efficient UEL mechanism in the PtOEP–DPA system, we believe that thermal vibrational excitation in DPA and PtOEP, and conformation-dependent S0–T1 energy gaps make the reverse TET process considerably more likely than is implied by the energies of the triplet-state origins in the equilibrium geometries of the molecules. The same phenomenon can occur in any TTA-UC system in which one or more of the components has a substantial number of low-frequency vibrational modes, and so we believe that reverse TET is likely to be an important UEL mechanism in many TTA-UC systems.

We also demonstrated that the rate constants for reverse TET and dynamic quenching can be determined given knowledge of [S]0,ideal, kRsens + kAq, and kAT. We showed that [S]0,ideal can be determined from plots of the upconversion quantum yield as a function of the sensitizer concentration at constant irradiance, and that kRsens + kAq and kAT can be determined from transient TTA-UC measurements performed at different sensitizer concentrations. This methodology should be applicable to a wide range of TTA-UC systems.

It will be of great interest to perform further experiments that can assess the role that thermal vibrational excitation plays in reverse TET in TTA-UC systems. One obvious approach is to perform temperature-dependent experiments in which the role of reverse TET can be assessed. However, a potentially confounding factor is such experiments is that other experimental parameters that are relevant to TTA-UC, such as the diffusion rate, are also sensitive to temperature.

Data availability

The data sets generated and/or analyzed in this study are available from the first/corresponding author on reasonable request.

Conflicts of interest

There are no conficts of interest to declare.

Acknowledgements

This work was supported by the National Science Foundation, grant CHE-1800491 (J. T. F. and A. K.). A. K. would also like to thank the University of Maryland's Graduate School and the Department of Materials Science and Engineering for supporting this research through a Summer Graduate Research Fellowship.

References

  1. Y. C. Simon and C. Weder, J. Mater. Chem., 2012, 22, 20817–20830 RSC.
  2. Y. C. Simon, S. Bai, M. K. Sing, H. Dietsch, M. Achermann and C. Weder, Macromol. Rapid Commun., 2012, 33, 498–502 CrossRef CAS PubMed.
  3. Q. Dou, L. Jiang, D. Kai, C. Owh and X. J. Loh, Drug Discovery Today, 2017, 22, 1400–1411 CrossRef CAS PubMed.
  4. R. W. MacQueen, Y. Y. Cheng, A. N. Danos, K. Lips and T. W. Schmidt, RSC Adv., 2014, 4, 52749–52756 RSC.
  5. X. Yu, X. Cao, X. Chen, N. Ayres and P. Zhang, Chem. Commun., 2015, 51, 588–591 RSC.
  6. D. K. Limberg, J.-H. Kang and R. C. Hayward, J. Am. Chem. Soc., 2022, 144, 5226–5232 CrossRef CAS PubMed.
  7. L. Zeng, L. Huang, J. Han and G. Han, Acc. Chem. Res., 2022, 55, 2604–2615 CrossRef CAS PubMed.
  8. C. Wohnhaas, V. Mailänder, M. Dröge, M. A. Filatov, D. Busko, Y. Avlasevich, S. Baluschev, T. Miteva, K. Landfester and A. Turshatov, Macromol. Biosci., 2013, 13, 1422–1430 CrossRef CAS PubMed.
  9. Q. Liu, M. Xu, T. Yang, B. Tian, X. Zhang and F. Li, ACS Appl. Mater. Interfaces, 2018, 10, 9883–9888 CrossRef CAS PubMed.
  10. W. Lin, J. Li, H. Feng, F. Qi and L. Huang, J. Anal. Test., 2023, 7, 327–344 CrossRef.
  11. Y. Y. Cheng, A. Nattestad, T. F. Schulze, R. W. MacQueen, B. Fückel, K. Lips, G. G. Wallace, T. Khoury, M. J. Crossley and T. W. Schmidt, Chem. Sci., 2016, 7, 559–568 RSC.
  12. T. Dilbeck and K. Hanson, J. Phys. Chem. Lett., 2018, 9, 5810–5821 CrossRef CAS PubMed.
  13. M. Kinoshita, Y. Sasaki, S. Amemori, N. Harada, Z. Hu, Z. Liu, L. K. Ono, Y. Qi, N. Yanai and N. Kimizuka, ChemPhotoChem, 2020, 4, 5271–5278 CrossRef CAS.
  14. T. F. Schulze, J. Czolk, Y.-Y. Cheng, B. Fückel, R. W. MacQueen, T. Khoury, M. J. Crossley, B. Stannowski, K. Lips, U. Lemmer, A. Colsmann and T. W. Schmidt, J. Phys. Chem. C, 2012, 116, 22794–22801 CrossRef CAS.
  15. W. Sheng, J. Yang, X. Li, G. Liu, Z. Lin, J. Long, S. Xiao, L. Tan and Y. Chen, Energy Environ. Sci., 2021, 14, 3532–3541 RSC.
  16. C.-H. Chen, Y.-S. Li, S.-C. Fang, B.-Y. Lin, C.-Y. Li, Y.-C. Liao, D.-G. Chen, Y.-R. Chen, Y.-C. Kung, C.-C. Wu, Y. Lin, W.-Y. Hung, T.-L. Chiu, C.-F. Lin, E. Y. Li, T.-F. Guo, J.-H. Lee, K.-T. Wong and P.-T. Chou, Adv. Photonics Res., 2023, 4, 2200204 CrossRef CAS.
  17. D. Di, L. Yang, J. M. Richter, L. Meraldi, R. M. Altamimi, A. Y. Alyamani, D. Credgington, K. P. Musselman, J. L. MacManus-Driscoll and R. H. Friend, Adv. Mater., 2017, 29, 1605987 CrossRef PubMed.
  18. L. Xing, Z.-L. Zhu, J. He, Z. Qiu, Z. Yang, D. Lin, W.-C. Chen, Q. Yang, S. Ji, Y. Huo and C.-S. Lee, Chem. Eng. J., 2021, 421, 127748 CrossRef CAS.
  19. L. Yang, X. W. Chua, Z. Yang, X. Ding, Y. Yu, A. Suwardi, M. Zhao, K. L. Ke, B. Ehrler and D. Di, Nanoscale Adv., 2022, 4, 1318–1323 RSC.
  20. Y. Murakami and K. Kamada, Phys. Chem. Chem. Phys., 2021, 23, 18268–18282 RSC.
  21. A. Monguzzi, R. Tubino, S. Hoseinkhani, M. Campione and F. Meinardi, Phys. Chem. Chem. Phys., 2012, 14, 4322–4332 RSC.
  22. E. A. Josie, C. Yuen Yap, K. Tony, G. C. R. C. Raphaël, N. J. Ekins-Daukes, J. C. Maxwell and W. S. Timothy, J. Phys.: Conf. Ser., 2009, 185, 012002 CrossRef.
  23. A. Monguzzi, J. Mezyk, F. Scotognella, R. Tubino and F. Meinardi, Phys. Rev. B: Condens. Matter Mater. Phys., 2008, 78, 195112 CrossRef.
  24. E. M. Gholizadeh, L. Frazer, R. W. MacQueen, J. K. Gallaher and T. W. Schmidt, Phys. Chem. Chem. Phys., 2018, 20, 19500–19506 RSC.
  25. T. Kashino, R. Haruki, M. Uji, N. Harada, M. Hosoyamada, N. Yanai and N. Kimizuka, ACS Appl. Mater. Interfaces, 2022, 14, 22771–22780 CrossRef CAS PubMed.
  26. S. P. Hill, T. Dilbeck, E. Baduell and K. Hanson, ACS Energy Lett., 2016, 1, 3–8 CrossRef CAS.
  27. A. Nattestad, Y. Y. Cheng, R. W. MacQueen, G. G. Wallace and T. W. Schmidt, J. Visualized Exp., 2014, e52028 Search PubMed.
  28. D. Beery, T. W. Schmidt and K. Hanson, ACS Appl. Mater. Interfaces, 2021, 13, 32601–32605 CrossRef CAS PubMed.
  29. V. Gray, K. Moth-Poulsen, B. Albinsson and M. Abrahamsson, Coord. Chem. Rev., 2018, 362, 54–71 CrossRef CAS.
  30. J. Isokuortti, S. R. Allu, A. Efimov, E. Vuorimaa-Laukkanen, N. V. Tkachenko, S. A. Vinogradov, T. Laaksonen and N. A. Durandin, J. Phys. Chem. Lett., 2020, 11, 318–324 CrossRef CAS PubMed.
  31. J. Peng, X. Guo, X. Jiang, D. Zhao and Y. Ma, Chem. Sci., 2016, 7, 1233–1237 RSC.
  32. M. Yang, S. Sheykhi, Y. Zhang, C. Milsmann and F. N. Castellano, Chem. Sci., 2021, 12, 9069–9077 RSC.
  33. Y. Zhao, Y. Wu, W. Chen, R. Zhang, G. Hong, J. Tian, H. Wang, D. Zheng, C. Wu, X. Jiang, X. Huo, L. Sun, W. Deng, K. Han and F. Song, Adv. Opt. Mater., 2022, 10, 2102275 CrossRef CAS.
  34. N. Yanai, M. Kozue, S. Amemori, R. Kabe, C. Adachi and N. Kimizuka, J. Mater. Chem. C, 2016, 4, 6447–6451 RSC.
  35. N. Nishimura, V. Gray, J. R. Allardice, Z. Zhang, A. Pershin, D. Beljonne and A. Rao, ACS Mater. Lett., 2019, 1, 660–664 CrossRef CAS.
  36. C. Wang, F. Reichenauer, W. R. Kitzmann, C. Kerzig, K. Heinze and U. Resch-Genger, Angew. Chem., Int. Ed., 2022, 61, e202202238 CrossRef CAS PubMed.
  37. X. Qiao, H. Liu, S. Xiao, P. Yuan, J. Yao, Y. Chen, P. Lu and D. Ma, Adv. Opt. Mater., 2022, 10, 2200074 CrossRef CAS.
  38. D. Meroni, A. Monguzzi and F. Meinardi, J. Chem. Phys., 2020, 153, 114302 CrossRef CAS PubMed.
  39. S. P. McGlynn, R. Sunseri and N. Christodouleas, J. Chem. Phys., 2004, 37, 1818–1824 CrossRef.
  40. M. N. Berberan Santos, Phys. Chem. Commun., 2000, 3, 18–23 RSC.
  41. V. Gray, A. Dreos, P. Erhart, B. Albinsson, K. Moth-Poulsen and M. Abrahamsson, Phys. Chem. Chem. Phys., 2017, 19, 10931–10939 RSC.
  42. A. Kalpattu, T. Dilbeck, K. Hanson and J. T. Fourkas, Phys. Chem. Chem. Phys., 2022, 24, 28174–28190 RSC.
  43. A. Olesund, V. Gray, J. Mårtensson and B. Albinsson, J. Am. Chem. Soc., 2021, 143, 5745–5754 CrossRef CAS PubMed.
  44. P. Baronas, J. Lekavičius, M. Majdecki, J. L. Elholm, K. Kazlauskas, P. Gaweł and K. Moth-Poulsen, ACS Cent. Sci., 2025, 11, 413–421 CrossRef CAS PubMed.
  45. S. Grimme, J. Antony, S. Ehrlich and H. Krieg, J. Chem. Phys., 2010, 132, 154104 CrossRef PubMed.
  46. F. Weigend and R. Ahlrichs, Phys. Chem. Chem. Phys., 2005, 7, 3297–3305 RSC.
  47. T. F. Schulze and T. W. Schmidt, Energy Environ. Sci., 2015, 8, 103–125 RSC.
  48. F. Edhborg, A. Olesund and B. Albinsson, Photochem. Photobiol. Sci., 2022, 21, 1143–1158 CrossRef CAS PubMed.
  49. J. S. Brinen and J. G. Koren, Chem. Phys. Lett., 1968, 2, 671–672 CrossRef CAS.
  50. Y. V. Aulin, M. van Sebille, M. Moes and F. C. Grozema, RSC Adv., 2015, 5, 107896–107903 Search PubMed.
  51. X. Cao, B. Hu and P. Zhang, J. Phys. Chem. Lett., 2013, 4, 2334–2338 CrossRef CAS.
  52. F. Laquai, C. Im, A. Kadashchuk and H. Bässler, Chem. Phys. Lett., 2003, 375, 286–291 CrossRef CAS.
  53. K. Guo, C. Lin, Y. Wu, S. Xiao, X. Qiao, D. Yang, Y. Dai, Q. Sun, J. Chen, D. Hu, L. Ying, Y. Ma and D. Ma, Adv. Opt. Mater., 2023, 11, 2202988 CrossRef CAS.
  54. K.-Y. Kuan and D. A. Singleton, J. Am. Chem. Soc., 2020, 142, 19885–19888 CrossRef CAS PubMed.
  55. C. B. Murphy, Y. Zhang, T. Troxler, V. Ferry, J. J. Martin and W. E. Jones, J. Phys. Chem. B, 2004, 108, 1537–1543 CrossRef CAS.
  56. F. Zapata, M. Marazzi, O. Castaño, A. U. Acuña and L. M. Frutos, J. Chem. Phys., 2014, 140, 034102 Search PubMed.
  57. J. Saltiel, J. L. Charlton and W. B. Mueller, J. Am. Chem. Soc., 1979, 101, 1347–1348 CrossRef CAS.
  58. F. Zapata, M. Nucci, O. Castaño, M. Marazzi and L. M. Frutos, J. Chem. Theory Comput., 2021, 17, 5429–5439 CrossRef CAS PubMed.

Footnote

Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d5cp00578g

This journal is © the Owner Societies 2025
Click here to see how this site uses Cookies. View our privacy policy here.