Nguyen Le Phuoca,
Arun Kumarb,
Aliaksandr Baidakbc,
Alexander S. Romanov*b and
Mikko Linnolahti
*a
aDepartment of Chemistry and Sustainable Technology, University of Eastern Finland, Joensuu Campus, FI-80101 Joensuu, Finland. E-mail: mikko.linnolahti@uef.fi
bDepartment of Chemistry, The University of Manchester, Oxford Rd, Manchester M13 9PL, UK. E-mail: alexander.romanov@manchester.ac.uk
cDalton Cumbrian Facility, The University of Manchester, West Lakes Science Park, Moor Row, CA24 3HA, UK
First published on 30th July 2025
Carbene–metal–amides (CMAs) have emerged as promising thermally activated delayed fluorescence (TADF) emitters for organic light-emitting diodes (OLEDs). Here, we present a comprehensive computational investigation of amide ligand effects on CMA photophysics, examining ca. 70 complexes through density functional theory and its time-dependent extension. Our systematic analysis reveals how structural modifications influence key parameters, including HOMO–LUMO overlap, singlet–triplet energy gaps, oscillator strengths, and metal–ligand bond energies. We demonstrate emission tunability across the visible spectrum through strategic modification of carbazole, indole, carboline, and guanidine-based amides. The computational screening identified promising candidates balancing TADF efficiency with molecular stability. Guided by these predictions, we synthesized two rationally designed complexes with contrasting excited state alignments, and tested their photo- and radioluminescence performance. This combined theoretical–experiment approach establishes clear structure–property relationships for CMA design and demonstrates the effectiveness of computational screening in accelerating OLED material development.
Since the introduction of the archetypical CMA1 in 2017,3 extensive studies have been reported, exploring structural modifications to enhance photophysical properties and OLED performance. On the carbene side, variations in cyclic (alkyl)(amino)carbene (CAAC) structure, including ring size alterations and substituent modifications, have demonstrated significant effects on emission wavelengths and quantum yields. Changes to the carbene π-system through benzannulation have enabled access to deep red emission, while modified ring sizes have influenced emission lifetimes and gaps between the highest occupied molecular (HOMO) and lowest unoccupied molecular (LUMO) orbitals.4–6 On the amide size, systematic modification of the carbazole core through electron-donating and -withdrawing substituents have achieved emission tuning from yellow to deep blue, with some derivatives showing near-unity quantum yields and sub-microsecond delayed fluorescence lifetimes.4,7
Although these experimental studies have demonstrated the versatility of CMA complexes as OLED emitters, a comprehensive understanding of the structure–property relationships governing their photophysical behavior remains elusive. Building on our previous investigation of CAAC effects in CMA emitters,24 we now present a systematic computational study examining how the amide structural modifications influence key photophysical parameters.
Using density functional theory (DFT)25,26 and time-dependent DFT (TD-DFT),27 we analyze diverse amide-structural motifs including modified carbazoles, expanded ring systems, and heterocyclic derivatives. Our investigation examines key parameters such as frontier orbital energies, singlet–triplet gaps, oscillator strengths, and bond dissociation energies. The computational analysis guides the synthesis of two rationally designed CMA complexes, demonstrating bright emission in agreement with theoretical predictions. By establishing clear structure–property relationships validated through targeted synthesis, this work provides design principles for high-performance CMA-based OLED materials and demonstrates the effectiveness of computational screening in accelerating emitter development.
![]() | ||
Fig. 1 Computed structures and properties of e-2-1 complex: (a) optimized ground state S0 geometry; (b) dipole moments in S0 and (c) S1@S0 state; (d) the HOMO and (e) LUMO orbitals; (f) optimized S1; (g) optimized T1; (h) vertical excitation energies indicated by the upward arrows; (i) fluorescence (left) and phosphorescence (right) energies indicated by the downward arrows at optimized S1 and T1 geometries.24 |
The ground state electronic structure reveals characteristic donor–acceptor behavior. The HOMO is predominantly localized on the carbazole amide with minimal (2.8%) gold contribution, while the LUMO distributes across the Au-CAAC moiety with 10.8% metal character (Fig. 1d and e). This orbital arrangement yields a HOMO–LUMO overlap integral of 0.35, favorable for TADF.39 The complex exhibits a substantial ground state dipole moment (10.7 D) oriented along the C–Au–N axis toward the carbene ligand. Notably, upon vertical excitation to S1@S0, the dipole moment not only decreases in magnitude (to 3.0 D) but also reverses direction (Fig. 1b and c). This dramatic reversal change reflects the charge-transfer (CT) nature of the excitation, as electron density shifts from the amide to the carbene moiety and may significantly influence the interaction between the complex and its environment in OLED devices.40 Vertical excitations analysis shows that both S0–S1 and S0–T1 transitions are dominated by HOMO-to-LUMO character (>94%), confirming their CT nature. The calculated energies (S1: 3.17 eV, T1: 2.88 eV) yield a ΔEST of 0.29 eV. Notably, the lowest triplet excitation (3LE) localized to the amide ligand (3LE(A)) corresponds to the second lowest excited triplet (T2) state at 3.35 eV, positioned above the CT states as required for efficient TADF.41
Geometry optimization of the excited states reveals significant structural and energetic changes from the ground state. While S0 and T1 geometries maintain a near-coplanarity between amide and carbene ligands, the S1 state prefers a perpendicular arrangement that lies 0.17 eV lower than the coplanar configuration. This structural reorganization has profound implications for the photophysical properties. The ΔEST gap decreases from 0.29 eV in the vertical excitation to 0.11 eV when determined from optimized geometries, potentially enhancing RISC efficiency.42 Simultaneously, the S1–S0 oscillator strength undergoes a dramatic reduction, falling from the ground state S0–S1 oscillator strength of 0.1952 to zero in the fully relaxed S1 geometry. The resulting emission energies are 2.25 eV for fluorescence (S1–S0@S1) and 2.34 eV for phosphorescence (T1–S0@T1), reflecting significant Stokes shifts43 (Fig. 1h and i). The vanishing oscillator strength in the relaxed S1 geometry, when considered alongside the experimentally observed high PLQY, suggests that emission likely occurs from a higher-energy near-coplanar configuration. The impressive TADF performance of the complex can thus be attributed to its ability to dynamically access geometries that strike a crucial balance between efficient RISC, enabled by the reduced ΔEST in the perpendicular configuration, and strong radiative transitions from near-coplanar geometries.
Analysis of the structural stability of the complex through bond dissociation energy calculations reveals robust metal–ligand interactions. The Au–C bond exhibits higher strength (403.2 kJ mol−1) compared to the Au–N bond (375.5 kJ mol−1), with the difference (ΔE = 27.7 kJ mol−1) identifying the metal–amide interface as the more vulnerable point for potential degradation. These values indicate overall strong metal–ligand coordination, crucial for maintaining structural integrity under typical OLED operating conditions.44
![]() | ||
Fig. 2 Structural variations of EtCAAC–Au–amide complexes investigated in this study, categorized by amide modification type. |
1. Carbazole modifications: we systematically functionalized the carbazole's six-membered rings at α, β, γ, and δ positions with diverse substituents including: electron-donating groups (OH, OCH3), electron-withdrawing groups (CF3, CN), alkyl groups (CH3, tBu), and halogens (Cl, Br). This comprehensive substitution strategy was motivated by experimental findings showing that tBu- and CF3-substituted carbazoles yield CMA emitters with tunable emission from yellow to deep blue.7
2. Ring expansions: we investigated the effect of expanding the carbazole's five-membered ring using heteroatoms (N, O, S) and alkyl groups. Previous work has demonstrated that such modifications enable broad spectral tuning, with PLQYs reaching 89% and successful implementation in yellow, sky-blue, and warm-white OLEDs.8
3. Indole and carboline derivatives: following reports of efficient blue emission from azacarbazole-based CMAs,9 we explored structurally related heterocycles. This series includes indole and carboline cores with established synthetic accessibility (e.g. skatole,45 indoline,46 harmane,47 pinoline48), allowing for potential experimental validation.
4. Guanidine amides: recent demonstrations of exceptional performance of rigid benzoguanidine-based CMAs10 prompted investigation of related structures. We examined both simple guanidine cores (G1–G6) and extended π-systems (G7–G13), targeting the promising deep-blue emission (424–466 nm), unity PLQY, and sub-microsecond delayed fluorescence demonstrated experimentally.10
Our computational analysis encompassed 70 complexes, providing systematic insights into how amide structural variations influence key PL parameters including excited state energies, HOMO–LUMO characteristics, and charge-transfer properties. This comprehensive investigation establishes design principles for optimizing CMA emitters, as detailed in subsequent sections.
The gold atom's participation in frontier molecular orbitals varies significantly across the investigated complexes. For the HOMO, primarily localized on the amide ligand, shows Au contributions ranging from 0.0% to 5.4%. The LUMO, distributed over the Au–EtCAAC moiety, typically exhibits larger Au character (10.0%–15.6%). The reference complex e-2-1 shows moderate Au participation of 2.8% and 10.8% in HOMO to LUMO, respectively.
Systematic modification of the carbazole structure reveals distinct trends in metal orbital contributions. Substitution at the α position generally leads to higher Au participation in frontier orbitals, with complex α_tBu showing the highest LUMO contribution (15.6%) across all studied compounds. The effect of β-substitution depends strongly on electronic character: electron-donating groups at this position (as in β_OH and β_OCH3) completely eliminate Au contribution to HOMO (0.0%), while electron-withdrawing substituents lead to moderate Au participation (2.8–3.3%). Substitution at γ- and δ-positions typically results in lower Au contribution to HOMO (γ: 2.4–3.2%, δ: 2.0–2.9%). Ring expansion of the carbazole core maintains moderate HOMO contributions (2.4–4.5%) while showing relatively consistent LUMO participation (11.1%–13.6%). Both indole and carboline derivatives show similar patterns of metal participation, with LUMO contributions (10.2–11.6%) comparable to the reference complex e-2-1. Guanidine-based amides exhibit the widest range of HOMO Au contributions (1.1–5.0%) while maintaining typical LUMO participation (10.3–12.2%). Notably, across all modifications, the variation in metal contribution to HOMO versus the consistency in LUMO reflects the fundamental electronic structure of CMAs: while the HOMO is predominantly localized on the modified amide ligand, the LUMO distributes across the unchanged Au–EtCAAC moiety.
The HOMO–LUMO overlap integral, crucial for TADF efficiency,49 varies from 0.15 to 0.45 across our dataset. The reference complex e-2-1 exhibits a moderate overlap integral of 0.35. Substitution effects are most pronounced at the β position, where electron-donating groups significantly reduce orbital overlap: complexes β_OH and β_OCH3 show the lowest values in our dataset (0.15 and 0.17, respectively). Ring-expanded complexes maintain overlap integrals (0.34–0.39) similar to the reference compound. Among heterocyclic variants, complex I4 exhibits the highest overlap integral (0.45) in our dataset, while most carboline derivatives show values comparable to e-2-1. Guanidine-based amides (G1–G13) show moderate to low overlap integrals (0.28–0.37). As amide modifications primarily affect the HOMO while leaving the Au-EtCAAC-centered LUMO largely unchanged, the observed variations in HOMO–LUMO overlap reflect the extent to which different amide structures alter the spatial and energetic characteristics of the HOMO.
To assess potential device stability, we calculated Au–N(amide) and Au–C(carbene) bond dissociation energies. The reference complex e-2-1 exhibits Au–N and Au–C bond energies of 375.5 and 403.2 kJ mol−1, respectively. Across the series, the Au–N bond shows greater variability (274.7–411.7 kJ mol−1) compared to the Au–C bond (373.7–424.6 kJ mol−1), reflecting its greater sensitivity to amide modifications. Electron-withdrawing substituents generally strengthen both metal–ligand bonds, exemplified by complex β_CF3 with the highest Au–N bond energy of (411.7 kJ mol−1) among carbazole derivatives. This strengthening can be attributed to the reduction of electron density on the amide nitrogen, enhancing the electrostatic interaction with the gold center. Conversely, electron-donating groups tend to weaken this interaction by increasing electron density on the nitrogen. Ring expansions significantly weaken the Au–N bond (274.7–365.2 kJ mol−1) while maintaining relatively strong Au–C bonds (377.0–416.2 kJ mol−1). This weakening is likely due to electronic effects rather than steric factors, with different heteroatoms or groups in the expanded ring altering the electronic structure of the entire π-system. Indole derivatives (I1–I4) show moderate Au–N bond strengths (318.4–378.0 kJ mol−1), while carboline-based complexes (C1–C10) maintain stronger metal–ligand bonds (Au–N: 348.0–403.9 kJ mol−1, Au–C: 395.8–412.2 kJ mol−1). Guanidine-based amides (G1–G13) exhibit lower Au–N bond energies (299.2–375.1 kJ mol−1) compared to the reference complex, which could be explained by the varying degrees of electron delocalization within the guanidine moiety. The markedly different responses of Au–N versus Au–C bonds to amide modifications reflect the localized nature of structural changes: while Au–N interaction is directly influenced by modifications to the amide, the Au–C bond to the unchanged carbene remains relatively stable.
Our computational results reveal systematic tuning of excited state energies through structural modification (Fig. 3). S1 energies span from 2.38 eV in the nitrogen-expanded carbazole complex (R_N) to 3.91 eV in the α-substituted complex (α_CF3), corresponding to emission wavelengths from 521 nm to 317 nm before accounting for Stokes shift effects. T1 states follow a similar pattern, ranging from 2.04 eV (R_N) to 3.32 eV (α_CF3), with S1/T1 energy ratios remaining notably consistent between 1.06–1.22 across the series.
Another excited state that warrants careful consideration relative to S1 and T1 states is T2. For efficient TADF, the optimal energetic hierarchy is T2 > S1 > T1. This ordering ensures that thermal equilibration occurs primarily between CT states S1 and T1, without interference from higher triplet manifolds. If T2 is positioned energetically below S1, it may serve as a competing decay channel, thereby disrupting the desired S1 ↔ T1 equilibrium essential for TADF. In such cases, excitons in the S1 state can undergo competing ISC to the lower-lying T2 state, removing population from the TADF-active S1 state. Additionally, thermal activation from T1 to S1 becomes less favorable when T2 provides an alternative, lower-energy triplet state. Furthermore, T2 states often exhibits distinct orbital character, potentially introducing additional non-radiative decay pathway. In the systems under investigation, T2 state energies (2.79–3.94 eV) typically lie above S1, thereby supporting favorable conditions for TADF. However, several compounds show T2 states below their respective S1 energies. For instance, in γ_OH, T2 (2.89 eV) lies below S1 (2.95 eV), and similar inversion occurs with many other substitution patterns. This unfavorable energy ordering could compromise TADF efficiency. The occurrence of such cases, particularly with electron-withdrawing substituents, highlights a key molecular design challenge: modifications targeting blue emission through raised S1 energies must simultaneously maintain appropriate excited state ordering.
Analysis of other structural modifications reveals distinct trends in excited state ordering. Ring-expanded complexes maintain favorable T2 > S1 ordering, with R_N and R_N(CH3) exhibiting the lowest S1 energies in our dataset (2.38–2.39 eV), though the series extends to higher energies up to 3.31 eV in R_SO2. Indole derivatives (I1–I4) similarly preserve T2 above S1 with substantial gaps (T2–S1 ≥ 0.19 eV), though at typically higher S1 energies (2.87–3.25 eV). Carboline derivatives, despite their promise as blue emitters with S1 energies up to 3.48 eV, show several cases of unfavorable energy ordering, particularly in C7 and C8, where T2 lies below S1 by 0.10–0.12 eV. In contrast, guanidine-based complexes (G1–G13) consistently maintain T2 above S1 across a wide range of S1 energies (2.64–3.56 eV), suggesting this modification strategy may be particularly robust for maintaining appropriate state ordering while tuning emission wavelength.
The TADF-critical ΔEST gap and S0–S1 oscillator strength show systematic variation with structural modification (Fig. 4). Among carbazole derivatives, α-substitution produces the widest range of ΔEST values, from α_tBu (0.17 eV) to α_CF3 (0.59 eV). Ring expansions generally maintain moderate ΔEST gaps (0.26–0.35 eV), while guanidine derivatives, particularly G9–11 and G13, achieve consistently small gaps (0.20–0.21 eV). These variations correlate with HOMO–LUMO overlap integrals: complexes with smaller overlaps typically show reduced ΔEST gaps but often at the cost of diminished oscillator strength. Notably, β_OH and β_OCH3 show minimal overlaps (0.15–0.17) and moderate ΔEST (0.29–0.30 eV), with reasonably high oscillator strengths (≈0.19), but their nearly degenerate S1 and T2 states (ΔE ≈ 0.03 eV) may complicate TADF performance. Complex γ_tBu maintains moderate overlap integral (0.35) and achieves high oscillator strength (≈0.22) while keeping ΔEST manageable (0.29 eV), with favorable T2–S1 gap of 0.13 eV. However, despite promising ΔEST (0.31 eV) and high oscillator strength (≈0.23), carboline derivative C5 shows unfavorable state ordering with T2 (3.22 eV) below S1 (3.29 eV). These results highlight the multiple requirements for efficient TADF emitter design: optimal complexes must balance moderate ΔEST gaps for efficient RISC, maintain sufficient oscillator strength for strong radiative decay, and crucially, preserve favorable excite state ordering with T2 above S1.
Electrostatic dipole moments provide insight into potential reorganization energies and TADF efficiency.40 Ground state dipole moments in all complexes align along the C(carbene)–Au–N(amide) axis towards the carbene, mirroring e-2-1 (10.7 D), but with magnitudes varying widely (4.1–17.0 D). Notably, guanidine-based complexes (G1–G6) show the smallest S0 dipole moments (4.1–5.4 D), while some γ- and δ-substituted carbazoles exhibit the largest (e.g. δ_CF3 with 17.0 D; γ_Cl, γ_Br, and γ_CF3 with 14.2–16.0 D). Upon S0–S1 excitation, most complexes exhibit a reversal in dipole moment direction, consistent with CT character, and often a significant change in magnitude. The S1@S0 dipole moments range from −8.0 D to 4.9 D. Large changes are observed in complexes like (α_CF3) (β_OCH3) (from 8.0 D to −7.2 D) and γ_Br (from 14.5 D to −0.8 D). These substantial reorganizations of charge distribution may impact TADF performance through enhanced non-radiative decay or altered solid-state packing interactions.40
Carbazole α-substitution offers the widest tunability but with important trade-offs for TADF design. Electron-donating groups (OH, OCH3, tBu) maintain favorable HOMO–LUMO overlap and reduce ΔEST gaps, supporting efficient RISC, while enhancing oscillator strength for bright emission. However, electron-withdrawing groups (CF3, Cl, Br) achieve substantial blue-shifts (up to 3.91 eV) but at the cost of increased ΔEST gaps (up to 0.59 eV) and reduced oscillator strength. The α-position shows the strongest electronic coupling to the metal center, with α_tBu exhibiting the highest Au contribution to frontier orbitals (15.6% LUMO). For blue emission applications, α-substitution requires careful balance between achieving target energies and maintaining TADF efficiency.
β-substitution provides unique electronic control through dramatic HOMO–LUMO decoupling. Electron-donating groups (β_OH, β_OCH3) create the most spatially separated donor–acceptor systems in our dataset (overlap integrals 0.15–0.17), potentially enabling very small ΔEST gaps while maintaining reasonable oscillator strengths (ca. 0.19). However, their nearly degenerate S1 and T2 states may complicate TADF dynamics. Electron-withdrawing groups offer substantial blue-shifts (up to 3.59 eV with β_CF3) while preserving electronic coupling, making this position ideal for high-energy blue emitters. Uniquely among all positions, β-substitution universally strengthens metal–ligand bonds, suggesting enhanced operational stability for device applications.
γ-substitution emerges as the “sweet spot” for TADF optimization, offering emission tuning without compromising key electronic parameters. Unlike other positions, γ-substitution maintains consistent HOMO–LUMO overlap (0.34–0.36) and favorable ΔEST gaps (0.27–0.31 eV) across all substituents, while uniquely preserving or enhancing oscillator strength. Complex γ_tBu exemplifies this balance with high oscillator strength (0.22), moderate ΔEST (0.29 eV), and favorable T2–S1 ordering. The position enables systematic emission tuning (2.95–3.55 eV) while maintaining TADF-critical parameters, making it ideal for fine-tuning emission color without sacrificing efficiency.
δ-substitution offers similar benefits to γ-substitution but with enhanced ΔEST control. Electron-donating groups achieve the smallest ΔEST gaps in the carbazole series (δ_OH: 0.24 eV), making this position attractive for maximizing RISC efficiency. However, electron-withdrawing groups produce detrimental large gaps (δ_CF3: 0.42 eV) with reduced oscillator strength, limiting blue emission applications. The position provides a viable alternative to γ-substitution when minimal ΔEST is prioritized over oscillator strength, though the latter trade-off must be carefully considered for overall TADF performance.
Ring expansion of the five-membered carbazole ring provides the most dramatic red-shifts in our dataset but with significant stability trade-offs. Nitrogen incorporation (R_N, R_N(CH3)) achieves the lowest S1 energies (2.38–2.39 eV) for accessing orange-red emission, while maintaining favorable ΔEST gaps for TADF. However, these systems suffer from substantially weakened Au–N bonds (274.7–277.1 kJ mol−1 vs. 375.5 kJ mol−1 reference) and reduced oscillator strength, potentially compromising device stability and brightness. The R_SO2 system offers a compromise, approaching reference S1 energy while providing the smallest ΔEST gap (0.26 eV) in the series. Ring expansion is most suitable for red-shifted emission applications where the stability trade-offs are acceptable.
Indole derivatives split into two distinct categories with contrasting design implications. Complexes I1–I3 maintain TADF-favorable properties similar to the reference (ΔEST: 0.27–0.29 eV, S1: 3.15–3.25 eV) but with reduced oscillator strength, making them suitable for blue emission with moderate brightness. In stark contrast, I4 represents an outlier with the highest HOMO–LUMO overlap in our dataset (0.45) and enhanced oscillator strength (0.25), but suffers from a prohibitively large ΔEST gap (0.52 eV) that likely prevents efficient TADF. The indole scaffold offers a reliable platform for blue emitters (I1–I3), while I4 serves as a cautionary example of how structural modifications can disrupt TADF balance.
Carboline derivatives demonstrate how nitrogen position critically affects TADF performance, creating two distinct design pathways. The C5–C6 series achieves excellent oscillator strength (0.20–0.23) with moderate ΔEST gaps, making them promising blue emitters despite slightly unfavorable T2 ordering. Conversely, C7–C8 suffer from poor oscillator strength (0.08–0.09) and large ΔEST gaps (0.44–0.45 eV), effectively eliminating TADF capability. The carboline scaffold offers high-performance blue emission potential (C5, C6) but requires careful nitrogen positioning to avoid TADF-detrimental configurations. Enhanced metal–ligand bonding across the series suggests improved stability compared to other heterocyclic modifications.
Guanidine amides reveal how π-conjugation extent fundamentally determines emitter performance, creating two complementary design strategies. Single π-bond systems (G1–G6) target deep-blue emission (3.33–3.56 eV) but suffer from large ΔEST gaps (up to 0.41 eV) and poor oscillator strength, limiting TADF efficiency. Extended conjugation systems (G7–G13) unlock exceptional TADF performance with small ΔEST gaps (0.20–0.21 eV) and tunable emission (2.64–3.27 eV), with G13 demonstrating optimal balance (ΔEST: 0.21 eV, oscillator strength: 0.18). The guanidine platform offers a unique design lever – π-extension – to transition from stable blue phosphors (G1–G6) to high-performance TADF emitters (G7–G13).
Based on this framework, computational screening identifies five leading candidates that exemplify optimal design principles. Complex γ_tBu demonstrates the γ-position advantage with balanced parameters across all metrics: moderate ΔEST (0.29 eV), high oscillator strength (0.22), and favorable T2–S1 ordering. α_OH achieves the highest oscillator strength in our dataset (0.25) while maintaining blue emission capability (3.21 eV). Carboline C6 balances strong emission properties (oscillator strength: 0.20) with enhanced stability, though T2 positioning requires consideration. Guanidine complexes G12 and G13 represent the success of the π-extension strategy, combining small ΔEST gaps (0.21–0.27 eV) with robust T2 > S1 ordering essential for uncompromised TADF.
To validate these computational predictions and understand the emission mechanisms, excited state geometry optimizations were performed for all five candidates. The optimized structures (Fig. 5) reveal distinct photophysical behaviors that explain their relative performance. Complexes γ_tBu, C6, G12, and G13 exhibit the expected TADF-favorable rotation from coplanar S0/T1 geometries to perpendicular S1 configuration, reducing ΔEST to 0.11–0.20 eV (Table 2). This rotation mechanism enables thermal S1 ↔ T1 equilibration essential for TADF, though requiring geometric reorganization (0.06–0.17 eV) to achieve maximum oscillator strength for emission. In contrast, α_OH demonstrates why α-substitution requires careful design consideration: steric constraints from the proximal hydroxyl group prevent S1 rotation, maintaining coplanar geometry and resulting in a large ΔEST (0.26 eV) that impedes RISC. All complexes maintain beneficial 3LE positioning (3.13–3.37 eV) above CT states, confirming appropriate excited state ordering for TADF applications.
α_OH | γ_tBu | C6 | G12 | G13 | ||
---|---|---|---|---|---|---|
a Vertical S1 and T1 excitations.b Optimized S1 and T1 geometries.c Optimized T1 geometry.d Carbene and amide fixed coplanar. | ||||||
Vertical S1 excitation energy (eV) | 3.21 | 3.09 | 3.31 | 3.16 | 3.17 | |
Vertical T1 excitation energy (eV) | 2.91 | 2.80 | 3.02 | 2.89 | 2.96 | |
ΔEST![]() |
0.30 | 0.29 | 0.29 | 0.27 | 0.21 | |
Energy relative to optimized S0 (eV) | Optimized coplanar S1 | 2.84 | 2.81 | 3.02 | 2.79 | 2.84 |
Optimized rotated S1 | 2.87 | 2.64 | 2.95 | 2.66 | 2.78 | |
Optimized coplanar T1 | 2.58 | 2.53 | 2.75 | 2.52 | 2.66 | |
ΔEST![]() |
0.26 | 0.11 | 0.20 | 0.14 | 0.12 | |
3LE energyc (eV) | 3.15 | 3.13 | 3.14 | 3.28 | 3.37 | |
Maximum S1–S0 oscillator strengthd | 0.0820 | 0.1670 | 0.1518 | 0.1136 | 0.1322 | |
Fluorescence (S1–S0@S1) (eV) | 2.13 | 2.16 | 2.32 | 2.05 | 2.31 | |
Phosphorescence (T1–S0@T1) (eV) | 2.18 | 2.24 | 2.45 | 2.08 | 2.31 | |
1CT and 3CT SOCME (cm−1) | 100.79 | 3.77 | 4.03 | 6.20 | 3.65 |
Spin–orbit coupling analysis provides additional mechanistic insight into these performance differences. The calculated SOCMEs between 1CT and 3CT states reveal significant variations across the series (3.65–100.79 cm−1, Tables 2 and S3). Complex α_OH exhibits the highest SOCME (100.79 cm−1), indicating rapid ISC from singlet to triplet state. However, its large ΔEST (0.26 eV) and inability to adopt the perpendicular S1 geometry hinder efficient RISC, negatively impacting TADF performance. In contrast, γ_tBu, C6, G12, and G13 show significantly smaller SOCMEs (3.65–6.20 cm−1) combined with favorable small ΔEST gaps (0.11–0.20 eV) that enable efficient RISC. Notably, G13 achieves optimal balance with the smallest SOCME (3.65 cm−1) and ΔEST gaps (0.12 eV), consistent with its superior experimental TADF performance (vide infra). These results suggest that moderate SOC strength is preferable for CMA emitters, providing sufficient ISC/RISC rates while minimizing competing non-radiative decay pathways.
The molecular structure for CMA complexes C5 and G13 were confirmed by single crystal X-ray diffraction (Fig. 6). Both C5 and G13 have a similar linear geometry for the gold atom with negligible deviation up to 2° from the ideal 180° angle for linear geometry. Unlike complex e-2-124 having a C2v symmetry for carbazolide amide donor moiety, complexes C5 and G13 possess unsymmetrically substituted amide donor. This results in crystallization of the complex C5 with two independent molecules A and B in the unit cell which are different by conformation of the ethyl-group for the EtCAAC–carbene π-acceptor ligand and relative orientation with the amide π-donor ligand. For instance, molecule B possesses methyl C-atom C25 pointed towards the center of the 2,6-diisopropylaniline moiety of carbene whereas C25 positioned in the space between of two ethyl groups of the carbene for molecule A. Complex G13 has only one independent molecules in the unit cell. Both Au–C and Au–N bond lengths for C5 and G13 are up to 0.02 Å longer compared to benchmark complex e-2-1, indicating greater separation between donor and acceptor ligands for the new CMA materials and influencing their photophysical behaviour, vide infra. The torsion angle (C1–N1–N2–C23, Fig. 6) between carbene and carbazole ligands varies in the range from 3.0(1) to 10.7(1)° which is close for complex e-2-1 (ca. 15.1(1)°), indicating near co-planar orientation of the ligands. Analysis of the intermolecular interactions between neighbouring molecules of complexes C5 and G13 reveals only weak C–H⋯π and C–H⋯N contacts resulting in formation of the 3D-network.
The redox behaviour of complexes C5 and G13 was analysed in THF solution with profiles shown in Fig. S6 and data collected in Table S5. Both complexes exhibit a quasi-reversible, carbene ligand-centred reduction process at E1/2 value at −2.80 V. All oxidation processes are centered on the amide and irreversible with Ep values at ca. +0.46 V (Table S5). This results in close similar values for LUMO at ca. 2.70 eV for both complexes but slightly more stabilised HOMO energy level at −5.66 eV for C5 compared to −5.56 eV for G13. The UV-vis absorption spectra in solvents with different polarity (Fig. S5, Table S4) show two major absorption bands: high-energy (ca. 300 nm) and low-energy absorption (ca. 360 nm) band. Complex C5 exhibits only well-resolved and narrow absorption bands for both high- and low-energy absorption bands with no solvatochromism effect. We ascribed both bands to π–π* transition localised on the amide-donor ligand, indicating absence of the low energy CT band. Unlike C5, complex G13 possesses a broad L(M)L CT low-energy absorption band, which experience 25 nm blue-shift from methylcyclohexane (MCH) to dichloromethane (DCM, Fig. S5, Table S4). Such a negative solvatochromism for complex G13 is characteristic for the CMA materials having bright TADF emission mechanism, vide infra.
The PL spectra were measured in neat crystals and polystyrene (PS) matrix at 295 K and 77 K (Table 3 and Fig. 7). Both complexes C5 and G13 emit blue light at 451 and 460 nm, respectively, which is up to 31 nm blue-shift compared with the e-2-1 in PS matrix24 (484 nm, Table 3). The PL profile for harmane complex C5 possesses poorly resolved vibronic structure (Fig. 7a top) which is becoming apparent upon cooling to 77 K (Fig. 7a bottom). Unlike C5, PS films of benzoguanidinato complex G13 possess a broad CT profile at all temperatures from 295 K to 16 K (Fig. 7b) thus demonstrating the PL behaviour similar to complex e-2-1. The CT-state energy for complexes C5 and G13 was estimated at ca. 3.06 eV from the blue onset of the PL profile in PS matrix at 295 K (Table 3). The excited state lifetime of benzoguanidinato CMA complex G13 is 0.38 μs which is three orders of magnitude shorter compared to harmane complex C5 (up to 298 μs) and nearly two-time shorter compared to carbazolide CMA complex e-2-1 (ca. 1.0 μs). Complex G13 possesses a radiative rate of 2.6 × 106 s−1 which is two-orders of magnitude faster compared to C5 (2 × 104 s−1, Table 3) thanks to a short excited state lifetime and near unity PLQY for G13 and only 40% PLQY for C5. Such drastic differences in PL performance clearly indicate a different luminescence mechanism governing the PL behaviour of the complexes C5 and G13.
λem (nm) | τ (μs) | Φa (%) | kr![]() |
knr![]() |
CT/3LE(A)/ΔEST![]() |
λem (nm, 77 K) | τ (μs, 77 K) | |
---|---|---|---|---|---|---|---|---|
a Quantum yields determined using an integrating sphere.b Radiative rate constant kr = Φ/τ.c Nonradiative constant knr = (1 − Φ)/τ.d CT/3LE(A) energies based on the onset values of the emission spectra blue edge in MeTHF glasses at 77 K and 295 K. | ||||||||
MeTHF frozen glass | ||||||||
C5 | — | — | — | — | — | —/2.87 | 450 | 2522 |
G13 | — | — | — | — | — | —/3.33 | 402 | 67 (90%), 1100 (10%) |
Polystyrene Matrix (1 wt%) | ||||||||
e-2-1 | 484 | 1.07 | 82 | 7.66 | 1.68 | 2.92/2.95/−0.03 | 458 | 48.0 (63%), 188.3 (32%) |
C5 | 451 | 40 (30%), 298 (70%) | 40 | 0.02 | 0.03 | 3.06/2.87/+0.19 | 449 | 401 (44%), 1162 (32%) |
G13 | 460 | 0.38 | 99 | 26 | 0.26 | 3.07/3.33/-0.26 | 448 | 55 |
To reveal the energy of the 3LE state, we performed the PL spectra for both complexes in 105M MeTHF frozen glass at 77 K (Fig. 7c). Both complexes possess a vibronically resolved PL profile indicative of a dominant phosphorescence from the 3LE(A) at 2.87 eV for C5 and 3.33 eV for G13. Both complexes possess the excited state lifetime in millisecond range supporting assignment of the phosphorescence at 77 K and in line with our previous observations for CMA complexes emitting phosphorescence from the locally excited triplet 3LE(A).10 The energy difference between CT and 3LE(A) states, ΔEST, is larger for compound G13 (−0.26 eV) compared with C5 (+0.19 eV), while 3LE(A) state is higher in energy compared with CT state for G13 (Table 3). Such experimental observation unequivocally corroborates our theoretical calculations above. For instance, complex G13 has the same excited CT and LE-states orientation with the reference compound e-2-1 (see Fig. 1), where triplet 3LE state is a T2 state laying above the CT singlet S1 and triplet T1. Theoretical calculations predict the energy of the triplet 3LE(A) to be 3.37 eV and energy gap of −0.21 eV (Table 2) which is in excellent agreement with the experimental value of 3.33 eV for 3LE(A) and −0.26 eV for the energy gap in complex G13. Therefore, theory and experiment enabled us to assign TADF emission mechanism for the complex G13. We also performed a variable-temperature lifetime measurement to estimate the activation energy barrier, ΔEa, for G13 in PS film. On cooling to 16 K, the PL profile remains broad and unstructured (Fig. 7b and d) while the emissive lifetime increases from 0.38 μs at 295 K to 55 μs at 16 K (Table 3, Fig. 7d). The ΔEa for G13 was estimated as 49.1 meV by an Arrhenius plot analysis of the varied temperature data with lnkTADF vs. 1/T, supporting high radiative rate of TADF.
Theoretical calculations for complex C5 predicted that both CT and 3LE(A) triplets are below in energy compared to CT singlet S1 state. This may result in a mixed character emission from CT and 3LE states simultaneously which is in a good agreement with the experimental observation for PS film of C5 – a rather broad CT emission profile with features of vibronic structure even at room temperature (Fig. 7a, Table 3).
Spectacular photophysical characteristics for complex G13 (unity PLQY and short excited state lifetime of 380 ns) spurred our interest to test the applied potential of G13 in radioluminescence experiment. This is important for development of the advanced scintillators for detecting high-energy X-ray radiation (radioluminescence, RL). Crystals of G13 were irradiated in an X-ray cabinet at 350 kV (max) and 11.4 mA current, with emission spectra measured using a fibre-optic cable connected to a fluorimeter. Complex G13 exhibits a RL emission spectra similar to its PL profile (Fig. 8). Radiostability was assessed by monitoring emission intensity at the compound's emission maxima every 10 seconds during constant irradiation at 16.357 Gy min−1 (total dose ca. 980 Gy). Complex G13 shows excellent stability with LT75 = 1 h (where LT75 represents the time required for RL intensity to decrease by 25% under constant 350 kV X-ray radiation). The radioluminescence stability for complex G13 sample is comparable with the industrial standard BGO (bismuth germanium oxide, powder sample) as shown in Fig. 8.
Our computational screening reveals remarkable emission tunability across the visible spectrum (2.38–3.91 eV), with systematic trends in how structural modifications affect emission energy. Electron-withdrawing substituents on carbazole consistently blue-shift emission, while electron-donating groups show position-dependent effects: ring expansions and extended π-conjugation provide routes to red-shifted emission, while minimal π-systems enable access to blue regions.
The predictive power of our computational approach is validated through the successful synthesis of two rationally designed complexes, C5 and G13, which exhibit bright phosphorescence and TADF emission, respectively, as predicted by their calculated electronic structures. These complexes, strategically selected for their contrasting excited-state alignments, demonstrate how computational insight can guide the targeted development of high-performance emitters. Complex G13 shows superior deep-blue TADF luminescence with radiative rates up to 2.6 × 106 s−1 enabling excellent radioluminescence stability under constant 350 kV X-ray radiation.
This combined theoretical–experimental study establishes a robust framework for CMA emitter design, providing clear guidelines for optimizing TADF performance through structural modification. The validated computational methodology offers an efficient pathway for screening candidate structures before synthesis, accelerating the development of next-generation OLED materials. These insights advance the broader goal of achieving highly efficient, stable, and color-tunable displays and lighting devices based on TADF technology.
Calculated parameters and experimental details, coordinates of optimized structures, cif-files and check CIF for complexes C5 and G13. See DOI: https://doi.org/10.1039/d5qi01245g.
CCDC 2454948 and 2454949 contain the supplementary crystallographic data for this paper.50a,b
This journal is © the Partner Organisations 2025 |