Yosuke
Tani
*ab,
Kiyoshi
Miyata
*c,
Erika
Ou
a,
Yuya
Oshima
a,
Mao
Komura
a,
Morihisa
Terasaki
a,
Shuji
Kimura
c,
Takumi
Ehara
c,
Koki
Kubo
c,
Ken
Onda
c and
Takuji
Ogawa
a
aDepartment of Chemistry, Graduate School of Science, Osaka University, Machikaneyama 1-1, Toyonaka, Osaka 560-0043, Japan. E-mail: y-tani@chem.sci.osaka-u.ac.jp
bInnovative Catalysis Science Division, Institute for Open and Transdisciplinary Research Initiatives (ICS-OTRI), Osaka University, Suita, Osaka 560-8531, Japan
cDepartment of Chemistry, Faculty of Science, Kyushu University, 744 Motooka, Nishi, Fukuoka 819-0395, Japan. E-mail: kmiyata@chem.kyushu-univ.jp
First published on 3rd June 2024
We report metal-free organic 1,2-diketones that exhibit fast and highly efficient room-temperature phosphorescence (RTP) with high colour purity under various conditions, including solutions. RTP quantum yields reached 38.2% in solution under Ar, 54% in a polymer matrix in air, and 50% in crystalline solids in air. Moreover, the narrowband RTP consistently dominated the steady-state emission, regardless of the molecular environment. Detailed mechanistic studies using ultrafast spectroscopy, single-crystal X-ray structure analysis, and theoretical calculations revealed picosecond intersystem crossing (ISC) followed by RTP from a planar conformation. Notably, the phosphorescence rate constant kp was unambiguously established as ∼5000 s−1, which is comparable to that of platinum porphyrins (representative heavy-metal phosphor). This inherently large kp enabled the high-efficiency RTP across diverse molecular environments, thus complementing the streamlined persistent RTP approach. The mechanism behind the photofunction has been elucidated as follows: (1) the large kp is due to efficient intensity borrowing of the T1 state from the bright S3 state, (2) the rapid ISC occurs from the S1 to the T3 state because these states are nearly isoenergetic and have a considerable spin–orbit coupling, and (3) the narrowband emission results from the minimal geometry change between the T1 and S0 states. Such mechanistic understanding based on molecular orbitals, as well as the structure-RTP property relationship study, highlighted design principles embodied by the diketone planar conformer. The fast RTP strategy enables development of organic phosphors with emissions independent of environmental conditions, thereby offering alternatives to precious-metal based phosphors.
Organic RTP must overcome several challenges stemming from the absence of heavy-metal atoms. The most severe is the inherently small phosphorescence transition probability (i.e., rate constant kp), causing poor RTP quantum yields (Φp). kp can be 104–105 s−1 for heavy-metal complexes,3 and only ∼10 s−1 or less for conventional organic compounds.1a Introducing heavy atoms, such as bromine, iodine, selenium, and tellurium, in conjunction with carbonyl functionalities, is a classical approach for enhancing kp.4 Despite its long history, however, an organic kp over ∼100 s−1 and/or a Φp over 1% in solutions are rarely observed (Fig. 1b, blue filled circles; Fig. S1 and Table S1†).5 The use of thiocarbonyls could provide a large kp in some cases (Fig. 1b, blue open circles).6 However, they are associated with high (photo)chemical reactivities and instability,7 thereby limiting their high Φp in perfluoroalkane solvents.6b The only practical way to improve Φp has been to reduce the nonradiative decay rate constant (knr). This approach has provided solid-state phosphors with RTP lasting for subseconds after excitation cessation (persistent RTP or afterglow).8 However, it inevitably restricts the organic phosphors to rigid media, such as crystals,9 or tailored host–guest systems with strong intermolecular interactions.10
A promising direction for environmentally independent organic RTP could be the significant increase in kps. Complementary to the persistent RTP, fast RTP could have significant potential use in solutions or any non-rigid molecular environment, including biological conditions.11 Moreover, because fast RTP can avoid efficiency loss caused by the accumulation of the triplet excitons, it also could be used in optoelectronics such as OLEDs.12 However, the quest for fast RTP is facing difficulty so far. In 2022, a selenium-containing molecule (phenoxaselenine) was designed as a candidate for fast RTP; however, its RTP was weak even in a polymer matrix, and only the photoluminescence spectrum and intrinsic phosphorescence rate (ca. 4000 s−1) at 78 K were reported.13 Quite recently, organic ionic crystals were reported to exhibit fast RTP with impressive kp values up to ∼105 s−1.14 The authors ascribed the large kp to the external heavy atom effect of two appropriately arranged iodide counter anions. Nonetheless, such effects relied on the crystal packing, and no RTP was observed in solutions. Thus, molecular design for a large kp that enables RTP in solution is needed.
Low colour purity is another organic RTP challenge. Considering the luminescence of the same quantum efficiency, high-colour-purity emission has a sharp and intense spectral peak that is visible, distinguishable, and vivid. These are essential features for OLED displays and bioimaging.15 However, organic RTP is usually broad with 80–120 nm full-width at half-maxima (FWHMs) (Fig. 1b, Table S1†). Moreover, because of an inefficient intersystem crossing (ISC) from singlet to triplet states, organic phosphors often exhibit both RTP and fluorescence, further impairing colour purity. As a long-known exception, 1,2-diketones, such as biacetyl and benzil, exhibit narrowband RTP that consists of one main peak with a small FWHM value accompanied by weak vibronic bands.16 However, the RTP is feeble in solution. Thus, simultaneously achieving a high Φp and colour purity in organic RTP remains a tremendous challenge.
Using heteroaromatic 1,2-diketones, we previously developed a series of solid-state mechanoresponsive RTP materials,17 solvent-free liquid RTP materials,18 and photoresponsive RTP crystals (Fig. S3†).19 They exhibited RTP in non-rigid molecular environments, such as amorphous solids or solvent-free liquids, wherein conventional metal-free organic compounds rarely show RTP. Thus, our previous results not only demonstrated unique advantages of the 1,2-diketone-based materials in applications, but also implied large kps. However, the fundamental molecular RTP properties remain elusive due to the complexity of the condensed materials. While we have briefly reported some solution-phase RTP properties of thienyl diketones 1a and 1b,17a further investigations were necessary to evaluate kp quantitatively (Fig. 1a and S2†). Typically, kp is derived from experimental Φp, ISC quantum yield ϕISC, and RTP lifetime τp, expressed as kp = Φp/(ϕISC × τp). However, ϕISC was not evaluated in our prior research as we did not focus on kp. In addition, while we determined Φp and τp of 1b, the previous protocol did not assure consistent degassing for those measurements (i.e., the extent of oxygen quenching of RTP would be different; Fig. S2†). Consequently, these values cannot be employed to establish kp. Therefore, the mechanism governing the expected fast RTP remained unexplored, not to mention the absence of the molecular design principle for substantial kp.
Here, we disclose the inherent molecular RTP properties of the thienyl diketone derivatives 1, exhibiting high-efficiency narrowband RTP based on an exceptionally large kps of ∼5000 s−1 (Fig. 1). In particular, 1a exhibited 38.2% Φp in cyclohexane. To our knowledge, this is the highest efficiency for metal-free organic molecules in common solvents. Relative to benzil, the kp was increased by a factor of ∼100 from 39 s−1 to 5300 s−1, enabling a 100-time increase in Φp from 0.31% to 38% (Fig. 1a).16c Moreover, high Φp with a narrow ∼30 nm FWHM was observed in solution (Fig. 1b). The high efficiency and colour purity were also observed in various conventional polymer matrices and crystalline solids, indicating that the RTP was an inherent molecular feature. The RTP mechanism and origin were revealed through single-crystal X-ray structure analysis, time-resolved photoluminescence (TRPL) spectroscopy, transient absorption (TA) spectroscopy, time-resolved infrared (TRIR) spectroscopy, and quantum-chemical calculations. Furthermore, structure-RTP property relationship study considering molecular- and natural transition orbitals provided design principles for fast organic RTP. Our work demonstrates the potential of metal-free organic molecular materials to exhibit fast RTP comparable to precious-metal complexes. Complementary to persistent RTP, organic fast RTP is promising not only for application in solution but also in functional condensed materials as we have witnessed in our previous studies (Fig. S3†).
To quantify the rapid ISC, we conducted femtosecond TA spectrum (fsTAS) measurement with ∼100 fs time resolution (Fig. 2c). The broad transient absorption over 450–650 nm changed to sharp spectra with a peak at 580 nm in a ∼10 ps timescale. We globally analysed the TAS by using a sequential model assuming two excited species, resulting in successful fitting with the time constant of the transition estimated to be 7.9 ps (Fig. 2d and e and S14†). We also conducted TAS measurements over nanosecond timescales (Fig. S13†), and observed that the TAS shape was unchanged up to the microsecond timescale. The long-lived spectral component was the lowest triplet excited state (T1 state), populated via ISC with a time constant of <10 ps for 1a in solution (Fig. 2f). Because the ISC outcompeted other relaxation processes, we can reasonably assume that the ISC quantum yield (ϕISC) was unity.
The total photoluminescence quantum yield can be expressed as ΦPL = Φf + Φp, where Φf is the fluorescence quantum yield. Given that ϕISC ∼ 1, phosphorescence dominates the total emission, and Φp ∼ ΦPL = 0.38. To the best of our knowledge,21–23 this is a record-breaking RTP efficiency for metal-free organic molecules in common solvents. In addition, with ϕISC ∼ 1 confirmed, Φp = ϕISCkpτp ∼ kpτp. For Φp = 0.38 and τp = 72.7 μs, kp was derived to be 5300 s−1, which greatly exceeded those of previously reported organic phosphors, except for thiocarbonyl compounds (Fig. 1a and Table S1†). Because of this exceptionally high kp, the RTP of 1a in solution was visible even in air (Φp = 2.6%, Fig. 2a and b, S10, Table S2, and ESI Movie 1†). We note that the estimation of kp in the literature often assumed ϕISC = 1 and Φp ∼ ΦPL without evaluating the time constant of ISC, and employed Φp and τp that were determined without assuring consistent degassing. This would be practical for qualitative purposes, but may cause significant errors in the kp value. Experimental evaluation of ϕISC and assuring consistent degassing are essential for the quantitative evaluation of kp.
The emission colour purity is also outstanding. Fluorescence was not discernible in the steady-state photoluminescence, and the FWHM of the RTP was only 29 nm for 1a (Fig. 2a and S10†). The RTP contained two other vibronic bands at around 615 and 680 nm; nonetheless, the first band at 560 nm covered 71% of the whole spectral area (Table S3†). This is distinct from other organic phosphors, which have broad and structure-less spectra with a 80–120 nm FWHM (Fig. 1c and Table S1†). Although there is a room for further improvement, the narrowband RTP with weak vibronic bands and the fluorescence-free nature represented a high colour purity emission; its coordinates (0.45, 0.54) are located almost at the edge of the chromaticity diagram of Commission internationale de l'éclairage (CIE) 1931 (Fig. 3b).
We introduced triisopropylsilyl (1b) and tributylsilyl groups (1c) onto the thienyl diketone core (Fig. 1a) and examined their properties in cyclohexane solution. The silylated derivatives 1b and 1c provided identical emission spectra with the maxima slightly red-shifted from 560 nm to 568 nm compared to those of 1a (Fig. 3a and S10†). The steady-state emissions were fluorescence-free and maintained a small FWHM (32 nm). This narrowband yellow emission had a CIE 1931 coordinate of (0.49, 051) (Fig. 3b).2c We also conducted fsTAS measurement for 1b (Fig. S15†) to analyse ISC; the main feature of the TAS was largely the same as that for 1a and the ISC time constants were 0.63 and 2.3 ps. The time constant is even faster than that of 1a and we can assume unity ϕISC for the silylated diketones. Φp for 1b and 1c was less than that for 1a, but still 24.6 and 17%, respectively (Fig. 3c). The decrease in Φp could be attributed to the doubled knr values, from 8500 s−1 for 1a to 15000 s−1 and 22000 s−1 for 1b and 1c, respectively. This was most likely because of the nonradiative decay accelerated by molecular motion involving the silyl groups. In contrast, kp values were only slightly changed, 5000 and 4500 s−1 for 1b and 1c. These results indicated that the large kp of 1 basically originated from the thienyl diketone core.
Moreover, phosphorescent films with Φp values of 38–18% were obtained by doping 1b into other five conventional polymers (Fig. 3f and S23, S24 and Table S4†). The emission maxima and FWHMs were almost unchanged, with colour purities independent of the polymers. Thus, RTP properties of the films were derived from inherent characteristics of 1.
We would like to emphasize that achieving efficient organic RTP in solution and in polymer matrices is challenging because nonradiative decay is not suppressed under these conditions. In such cases, the fast RTP (large kp) has a distinct advantage over persistent RTP with small kp. The same can be said for the RTP in non-crystalline aggregate states, such as the amorphous solid state and solvent-free liquid state. It becomes clear that our previous studies of mechanoresponsive RTP materials,17 solvent-free liquid RTP materials,18 and photoresponsive RTP crystals19 (Fig. S3†) demonstrated the usefulness and applications of fast RTP.
In the present work, we found that crystals of 1c exhibited yellow RTP with a 50% Φp in air (Fig. 3g and h). The RTP decayed as a single exponential with a τp of 79 μs, and kp was 6000 s−1 (Fig. S25†). Most importantly, the PL spectrum was almost identical to that of 1c in solution, indicating emission from the same conformer (Fig. 3g). Single-crystal X-ray structure analysis revealed the TP conformation of 1c, with the thienyl diketone core spatially separated from neighbouring molecules by tributylsilyl substituents (Fig. 3i and S32†; no π-stacking, and the interplanar distance was 4.45 Å). Thus, the RTP from the 1c crystal, as well as 1 in solution, was unambiguously assigned to the monomer emission from its TP conformer. Notably, 1c represents a rare example of metal-free organic molecules exhibiting remarkable Φp in solution, in a polymer matrix, and in a crystal, with almost identical PL spectra.
To experimentally elucidate conformation dynamics in the excited state, TRIR spectroscopy was performed. Although we could not perform the TRIR measurement on 1a because of poor solubility, we managed to conduct the TRIR measurement on 1b in cyclohexane using a 267 nm wavelength pump light (Fig. S16†). The spectra converged to sharp spectra over 0.74 ps and 2.67 ps timescales, consistent with the time constants extracted from fsTAS of 1b (0.63 and 2.3 ps, Fig. S15†). Note that these time constants were comparable to those observed for 1a in the fsTAS measurement. This indicated that the excited-state dynamics were primarily affected by the thienyl diketone core, while the silyl groups had minor effects. The two timescales were attributed to the conformation change and ISC. More importantly, the converged TRIR spectra were consistent with simulated spectra, assuming a TP conformation in the T1 state (Fig. S17†). Thus, both the conformation changes and the ISC were ultrafast, and the RTP is confirmed to originate from the TP conformer.
(1) |
(2) |
(3) |
TDDFT calculations for the T1-optimized TP geometry of 1a yielded kp = 5400 s−1 (for n = 1–6), which excellently reproduced the experimental value of 5300 s−1. Remarkably, the T1 state strongly coupled with the S3 state (Fig. 5a); a large SOC matrix element 〈S3|HSO|T1〉 = 167 cm−1 and a reasonable ΔES3–T1 = 1.74 eV = 14000 cm−1 provided a large mixing coefficient 〈S3|HSO|T1〉/ΔES3–T1 of 1.19 × 10−2. This coefficient indicates that the T1 state can borrow 1.19% of the transition dipole moment between the S3 and S0 states (i.e., μS3–S0), which is quite large for phosphorescence. Furthermore, the |μS3–S0| was as large as 4.68 D, thus realizing a large |mp3| = 56 × 10−3D. The contribution from mp3 corresponded to 93% of the total kp value and is obviously the source of the large kp (Fig. 6b, vide infra). Similar results were obtained for 1b, confirming that the origin of the substantial kp lies in the thienyl diketone core (Table S5†).
The natural transition orbitals (NTOs) of the S3 and T1 states provided an intuitive understanding of the large μS3–S0 and 〈S3|HSO|T1〉 (Fig. 5b). First, the S3 state (S3–S0 transition) had a pure (π,π*) configuration, where the NTOs delocalized over the entire molecule in a planar geometry, which was favourable for large μ. On the other hand, the T1 state had a pure (n,π*) configuration, i.e., it was not contaminated with (π,π*) character. The pure 3(n,π*) configuration enabled SOC with the pure 1(π,π*) state, consistent with El Sayed's rule.4a The pure (n,π*)/(π,π*) characters are noteworthy, because excited states with mixed (n,π*) and (π,π*) configurations diminish SOC.30 Most importantly, the hole NTOs of the S3 and T1 states around the heavy atoms, especially Br, were mutually perpendicular. This is ideal because the spin–orbit Hamiltonian HSO involves the orbital angular-momentum operator that rotates the orbitals by 90°.13,28a Thus, for the hole-constituting atomic orbital of Br, the in-plane σ-symmetric n orbital of T1 states (Fig. 5b bottom-right) became π-symmetric, producing a significant spatial integral with the out-of-plane p orbital of the S3 state (Fig. 5b bottom-left). The large spatial integral near the heavy nuclei boosts the heavy-atom effect, resulting in a large 〈S3|HSO|T1〉 = 167 cm−1. Thus, 1 embodies an ideal electronic structure for intensity-borrowing that enhances kp: (1) a planar π-system that produces large μ; (2) carbonyl groups in the π-plane that joined the pure 1(π,π*) and pure 3(n,π*) states; and (3) heavy atoms with electrons conjugated to both the π and σ-/n-systems to boost 〈Sn|HSO|T1〉.
The superior Φp of 1 was attributed to its large kp, which was more than one order of magnitude larger than the estimated kps of 2 and 3 (Fig. 6a and Table S2†). TDDFT calculations reasonably reproduced this trend; the calculated kp for 1a, 2a, and 3a was 5400, 9, and 71 s−1, respectively (Fig. S27, S30 and S31†). The small kp of aldehyde 2 was because of the mismatched 〈Sn|HSO|T1〉 and μSn–S0 (Fig. 6b; see Table S4† for 1b, 2b, and 3b). Thus, the T1 state of 2 had a (π,π*) configuration (Fig. 6e), which coupled strongly to (n,π*) Sn states (e.g., n = 5; 〈S5|HSO|T1〉 was as large as 247 cm−1 for 2a). However, the (n,π*) configuration diminishes μSn–S0 (e.g., |μS5–S0| was only 0.22 D for 2a). Consequently, intensity-borrowing was insufficient, resulting in a small kp. Hence, a large 〈Sn|HSO|T1〉 (or a large mixing coefficient, 〈Sn|HSO|T1〉/ΔESn–T1) itself is insufficient and the combination with a large μSn–S0 in the same n-th singlet state is essential for achieving a huge kp.31 To meet this requirement, the T1 state should have a (n,π*) configuration. Indeed, the kp of Br-free diketone 3a was larger than that of Br-containing aldehyde 2a, because the T1 state of 3a had a (n,π*) configuration. Thus, the (n,π*) T1 state prefers coupling with (π,π*) Sn states, whose μSn–S0 can be sizable (e.g., n = 3, |μS3–S0| = 5.16 D for 3a; Fig. 6b). However, SOC between these states was not as effective as that for 1a because of the absence of the heavy atom, Br (e.g., 〈S3|HSO|T1〉 = 23 cm−1 for 3a). Therefore, intensity borrowing was less effective, yielding a lower kp than the brominated diketone 1a.
The side-on spatial arrangement of C–Br and CO was crucial for the heavy-atom effect because of better mutual alignment of the in-plane Br and O p orbitals. This resulted in large hole coefficients in the T1 state (Fig. 5b).32 However, the arrangement is energetically unfavourable due to electronic repulsion, because Br and O in the TP geometry of 1 were within the sum of their van der Waals radii according to the X-ray structure (Fig. S32†).17a,33,34 Indeed, DFT calculations for the aldehyde 2a indicated that the side-on conformation was less stable (by 3 kcal mol−1) for both S0 and T1 states than the conformation with oxygen pointing away from the Br (Fig. S28†).32 Interestingly, the TP geometry of 1 involved two-fold intramolecular S⋯O noncovalent chalcogen bonding interactions,35 with interaction energies as large as 7.74 kcal mol−1, according to the natural bonding orbital analysis of 1a.17a The chalcogen bonds stabilized the TP conformer and increased kp by forcibly fixing Br in the vicinity of carbonyl oxygens.
The bandwidths could be strongly correlated with molecular geometry changes in the ground and excited states, which are quantified using the Huang–Rhys (HR) factor obtained from a Franck–Condon analysis of the spectra (Fig. 6c, S18, and Table S3†).15,36 The HR factor of diketone 1a (1b) was 0.51 (0.50), while that of aldehyde 2a (2b) was 2.22 (2.39). The smaller HR factors in 1 corresponded to smaller geometry changes in the T1-to-S0 transition, which were corroborated by DFT calculations. The results indicated that although diketone 1a and aldehyde 2a retained their planarity in the transition, bond length changes were small for 1a and large for 2a.
We visualized these bond length changes using elongation ratios, defined as (RT1 − RS1)/RS1, where RT1 and RS1 represent bond lengths in the T1-and S0-optimized geometries; a positive value indicates that the bond is longer in the T1 state (Fig. 6d). The bond lengths within the thiophene ring of diketone 1a exhibited minimal changes. Conversely, the absolute values of the elongation ratios were large in the case of aldehyde 2a (e.g., the bond length between 3C and 4C changed by 8.4%).
These considerable bond length changes in 2a were attributed to the (π,π*) character of the T1 state, as evident in the NTOs (Fig. 6e). The electronic transition from bonding to antibonding orbitals significantly affects the bond strength. In contrast, the NTOs of 1a were mainly localized on the 1,2-dicarbonyl moiety, and the electronic transition was nonbonding to antibonding, which led to minor changes in the bond lengths (Fig. 6d). Moreover, the transition in the diketones 1 had no net charge-transfer character because of the centrosymmetric geometry. This was in contrast to most aromatic (mono)carbonyl compounds, where 3(n,π*) states have a charge-transfer character from the carbonyl to the aromatic ring. These features of diketones minimized molecular geometry changes during phosphorescence emission. Therefore, the 1,2-diketone-based phosphor would be a promising platform for efficient narrowband RTP.
1. The planar geometry with a delocalized π-system
2. σ-Symmetric n-orbitals perpendicular to the π-system
3. Heavy atoms conjugated with both π- and n-electron systems
4. (n,π*) configuration of the T1 state
Fulfilling these points would be a promising design principle for molecules with large μSn–S0 and 〈Sn|HSO|T1〉 in the same Sn states, and results in a kp leap. In addition, the centrosymmetric nonbonding-to-antibonding electronic transition of the 1,2-diketone skeleton was the key to the narrowband emission. The fast RTP (large kp) is potentially advantageous for various applications, such as phosphorescence OLEDs and lasers, triplet-to-singlet conversion via Förster resonance energy transfer, bioimaging, and theranostics. However, such applications have been monopolized by precious-metal phosphors. We believe that the mechanistic elucidation paves the way for developing fast organic RTP, which could expand applications of metal-free organic materials.
Footnote |
† Electronic supplementary information (ESI) available: Experimental procedures, analytical data, physicochemical properties, computational details, and CIF files for the single-crystal X-ray structure analysis. CCDC 2269866. For ESI and crystallographic data in CIF or other electronic format see DOI: https://doi.org/10.1039/d4sc02841d |
This journal is © The Royal Society of Chemistry 2024 |