Lowest Singlet Research Articles

B-embedded disulfide-bridged π-conjugated compounds: structures and optical tuning.

Two novel B-embedded disulfide-bridged π-conjugated compounds (BS-CZ and BS-N) bearing different electron donor groups (phenyl carbazole and triphenylamine) have been prepared and show different optical mechanisms. The compound BS-CZ exhibits significant multiple resonance thermal activation delayed fluorescence (MR-TADF) properties with a small singlet-triplet energy gap (ΔEST = 0.16 eV) and a narrow half-peak full width (FWHM = 33 nm), while the compound BS-N shows traditional fluorescence luminescence (FL) characteristics with a larger ΔEST (0.28 eV) and FWHM (57 nm). Time-dependent density functional theory (TD-DFT) calculations show that the lowest excited singlet state (S1) of the compound BS-CZ exhibits local excited (LE) state characteristics, while the charge transfer (CT) state characteristics can be found in S1 of the compound BS-N. Considering good optical performance, the compound BS-CZ is used as an emitting layer of the organic light-emitting diode device and achieved saturated blue emission (473 nm) with a narrow FWHM (39 nm), and CIE color coordinates of (0.12, 0.21). This work provides an important strategy for the optical mechanism regulation and photoelectric applications of B-embedded disulfide-bridged π-conjugated molecules.

UV-spectrum and photodecomposition of peroxynitrous acid in the troposphere.

The UV spectrum of peroxynitrous acid, HOONO, was computed at the B3LYP/AVTZ and MCSCF/AVTZ levels using the fewest switches surface hopping algorithm. Due to large-amplitude vibrational motions of this molecule, the maxima in the simulated spectra are displaced from the positions of vertical excitations. The three lowest excited electronic singlet states, which are all repulsive, can be reached by UV absorption. The photolysis products are determined, and the photolysis rate constant is provided for the first time. We found that near the tropopause the photolysis rate constant J ≈ 6 × 10-4 s-1, exceeds that for thermal decomposition by two orders of magnitude. The photolysis lifetime is about 30 minutes. Thus, photolysis is an important process and should be included in atmospheric models.

The effect of heavy atoms on the deactivation of electronic excited states of dye molecules near the surface of metal nanoparticles.

The influence of a heavy atom and the plasmon field on the efficiency of populating the lowest triplet state (T1) and on the phosphorescence intensity has been studied for fluorescein, 2Br-fluorescein, eosin and erythrosine, which have an increasing number of substituted heavy atoms. We show that the heavy atoms affect not only the rate constant of intersystem crossing (kISC) but also the rate constant of internal conversion (kIC). The calculations show that the C-H bonds in the meso position are the primary acceptors of the excitation energy of the lowest excited electronic singlet state (S1). Substitution of the meso hydrogen atoms with I or Br leads to a smaller kIC rate constant of 1 × 108 s-1 for fluorescein to 8 × 106 s-1 for eosin. Substitution with heavy atoms also leads to a larger ISC rate constant (kISC) between the T2 and S1 states because the spin-orbit coupling matrix element 〈S1|HSO|T2〉 increases by two orders of magnitude from 0.36 cm-1 for fluorescein to 35.0 cm-1 for erythrosine. The phosphorescence rate constant increases by three orders of magnitude from 4.8 × 10 s-1 for fluorescein to 3.3 × 104 s-1 for erythrosine, which is supported by experimental data. The plasmon effect increases the intensity of the xanthene dye emissions. The intensity and the quantum yield of fluorescence increase in the series fluorescein < 2Br-fluorescein < eosin < erythrosine. The intensity of the delayed fluorescence and phosphorescence grows in the same way. The enhancement factor of the phosphorescence intensity increases from 1.8 to 5.6 in the series from fluorescein to erythrosine. The differences in the plasmon effect originate from intensity borrowing to the radiative triplet-singlet transition (T1 → S0) from the singlet-singlet transitions (Sn → S0), which is more efficient when molecules have heavy atoms in the meso position.

Fully quantal description of combined internal conversion and intersystem crossing processes in the smallest Criegee intermediate CH2OO.

A quantal description of nuclear motion using coupled fifteen-state potential energy and spin-orbit coupling surfaces for studying the photodissociation of CH2OO to H2CO(X1A1) + O1D and H2CO(X1A1) + O3P channels is presented. For the evaluation of surfaces, multireference electronic wave functions are employed. For the fully quantal description of the nuclear motion, we diabatize the PESs of the two and four lowest excited singlet and triplet states, respectively, within the three sets of vibronically coupled states, i.e. (B1A', C1A'), (a3A', b3A') and (a3A'', b3A''), employing the diabatization by ansatz method. This yields three different adiabatic-to-diabatic mixing angles, which are used for the diabatization of the spin-orbit coupling surfaces and allow to investigate simultaneously the internal conversion and intersystem crossing processes in CH2OO using a nuclear quantum dynamical approach for the first time. Our calculation predicts the presence of a weak spectral band with irregular and discrete structures, which originates from the role of spin-orbit couplings. This band of the spectrum is mainly located between the minimum energy of the a3A' state and the onset of the B1A' ← X1A' spectral band. Furthermore, sizable SOC between the B1A' and a3A'' states mixes them via intersystem crossing. Due to the vibronic interaction between the a3A'' and a3A' states, the molecule finally relaxes to the a3A' state and is dissociated to H2CO(X1A1) and O3P products.

Theoretical elucidation of the effect of the linkage and orientation of carbazole and naphthalenediimide on the TADF and RTP propensities.

Thermally activated delayed fluorescence (TADF) and room temperature phosphorescence (RTP) are most promising processes for harvesting triplet excitons in organic light-emitting diodes. In this work, the effect of the linkage between the carbazole (Cz) donor (D) and the naphthalenediimide (NDI) acceptor (A) on the TADF and RTP propensities is elucidated using density functional theory computations employing D-A, D-A-D, D-π-A, and D-π-A-π-D structural designs. The effects of the dihedral angle between the donor and acceptor units on the energy difference between the singlet and triplet excited states (ΔEST), the spin-orbit coupling (SOC) constants, and radiative (kr), intersystem crossing (kISC) and reverse intersystem crossing (kRISC) rates are unravelled. The molecules possessing a direct linkage between Cz and NDI exhibit large ΔEST values due to substantial overlap between the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO). However, the insertion of a phenyl spacer between the Cz and NDI units led to disjoint HOMO and LUMO and consequently resulted in a small ΔEST. Furthermore, the presence of two donors with and without a phenyl spacer on NDI resulted in a high-lying triplet state (T2) that is energetically lower than the lowest singlet excited state (S1), hence providing additional channels to the TADF and RTP processes. Also, the orientation of Cz and NDI in the ortho-positions of the phenyl unit resulted in a T1 state with dominant LE character which led to moderate spin-orbit coupling constants and highest kr rates compared to the analogous meta- and para-linked derivatives. Thus, the ortho-derivatives possessed small ΔEST, charge transfer dominated S1, joint holes and electrons for the T1 state, characteristic of local excitation, high SOC, and promising rISC and kr rates. Overall, the phenyl linked derivatives possess TADF characteristics, while the directly linked analogues show RTP propensity.

Sulphur-atom positional engineering in perylenimide: structure-property relationships and H-aggregation directed type-I photodynamic therapy.

An innovative design strategy of placing sulfur (S)-atoms within the pendant functional groups and at carbonyl positions in conventional perylenimide (PNI-O) has been demonstrated to investigate the condensed state structure-property relationship and potential photodynamic therapy (PDT) application. Incorporation of simply S-atoms at the peri-functionalized perylenimide (RPNI-O) leads to an aggregation-induced enhanced emission luminogen (AIEEgen), 2-hexyl-8-(thianthren-1-yl)-1H-benzo[5,10]anthra[2,1,9-def]isoquinoline-1,3(2H)-dione (API), which achieves a remarkable photoluminescence quantum yield (Φ PL) of 0.85 in aqueous environments and established novel AIE mechanisms. Additionally, substitution of the S-atom at the carbonyl position in RPNI-O leads to thioperylenimides (RPNI-S): 2-hexyl-8-phenyl-1H-benzo[5,10]anthra[2,1,9-def]isoquinoline-1,3(2H)-dithione (PPIS), 8-([2,2'-bithiophen]-5-yl)-2-hexyl-1H-benzo[5,10]anthra[2,1,9-def]isoquinoline-1,3(2H)-dithione (THPIS), and 2-hexyl-8-(thianthren-1-yl)-1H-benzo[5,10]anthra[2,1,9-def]isoquinoline-1,3(2H)-dithion (APIS), with distinct photophysical properties (enlarged spin-orbit coupling (SOC) and Φ PL ≈ 0.00), and developed diverse potent photosensitizers (PSs). The present work provides a novel SOC enhancement mechanism via pronounced H-aggregation. Surprisingly, the lowest singlet oxygen quantum yield (Φ Δ) and theoretical calculation suggest the specific type-I PDT for RPNI-S. Interestingly, RPNI-S efficiently produces superoxide (O2˙-) due to its remarkably lower Gibbs free energy (ΔG) values (THPIS: -40.83 kcal mol-1). The non-toxic and heavy-atom free very specific thio-based PPIS and THPIS PSs showed selective and efficient PDT under normoxia, as a rare example.

Thermally activated delayed fluorescence tetradentate ligand-containing gold(III) complexes with preferential molecular orientation and their application in organic light-emitting devices.

A new class of thermally activated delayed fluorescence (TADF) pyridine-/pyrazine-containing tetradentate C^C^N^N gold(III) complexes have been designed and synthesized. Displaying photoluminescence quantum yields (PLQYs) of up to 0.77 in solid-state thin films, these complexes showed at-least a six-fold increase in the radiative decay rate constant (kr) in toluene upon increasing temperature from 210 to 360 K. Using variable-temperature (VT) ultrafast transient absorption (TA) spectroscopy, the reverse intersystem crossing (RISC) processes were directly observed and the activation parameters were determined, in line with the results of the Boltzmann two-level model fittings, in which the energy separation values between the lowest-lying singlet excited state (S1) and the lowest-lying triplet excited state (T1), ΔE(S1-T1), of these complexes were estimated to be in the range of 0.16-0.18 eV. Through strategic modification of the position of the electron-donating -tBu substituent in the cyclometalating ligand, the permanent dipole moments (PDMs) of these tetradentate gold(III) emitters could be manipulated to enhance their horizontal alignment in the emitting layer of organic light-emitting devices (OLEDs). Consequently, the resulting vacuum-deposited OLEDs demonstrated a 30% increase in the theoretical out-coupling efficiency (ηout), as well as promising electroluminescence (EL) performance with maximum external quantum efficiencies (EQEs) of up to 15.7%.

Optimizing the energy level alignment for achieving record-breaking efficiency in hot exciton deep red OLEDs.

To effectively compete with the quenching process in long-wavelength regions like deep red (DR) and near-infrared (NIR), rapid radiative decay is urgently needed to address the challenges posed by the "energy gap law". Herein, we confirmed that it is crucial for hot exciton emitters to attain a narrow energy gap (ΔES1-T2) between the lowest singlet excited (S1) state and second triplet excited (T2) state, while ensuring that T2 slightly exceeds S1 in the energy level. Two proofs-of-concept of hot exciton DR emitters, namely αT-IPD and βT-IPD, were successfully designed and synthesized by coupling electron-acceptors N,N-diphenylnaphthalen-2-amine (αTPA) and N,N-diphenylnaphthalen-1-amine (βTPA) with an electron-withdrawing unit 5-(4-(tert-butyl) phenyl)-5H-pyrazino[2,3-b]indole-2,3-dicarbonitrile (IPD). Both emitters exhibited a narrow ΔES1-T2, with T2 being slightly higher than S1. Additionally, both emitters showed significantly large ΔET2-T1. Moreover, due to their aggregation-induced emission characteristics, J-aggregated packing modes, moderate strength intermolecular CN⋯H-C and C-H⋯π interactions, and unique, comparatively large center-to-center distances among trimers in the crystalline state, both αT-IPD and βT-IPD emitters exhibited remarkable photoluminescence quantum yields of 68.5% and 73.5%, respectively, in non-doped films. Remarkably, the corresponding non-doped DR-OLED based on βT-IPD achieved a maximum external quantum efficiency of 15.5% at an emission peak wavelength of 667 nm, representing the highest reported value for hot exciton DR-OLEDs.

Nature of ultrafast dynamics in the lowest-lying singlet excited state of [Ru(bpy)3]2.

This work presents mechanisms to rationalize the nature of ultrafast photochemical and photophysical processes on the first singlet metal-ligand charge transfer state (1MLCT1) of the [Ru(bpy)3]2+ complex. The 1MLCT1 state is the lowest-lying singlet excited state and the most important intermediate in the early evolution of photoexcited [Ru(bpy)3]2+*. The results obtained from simple but interpretable theoretical models show that the 1MLCT1 state can be very quickly formed via both direct photo-excitation and internal conversions and then can efficiently relax to its equilibrium geometry in ca. 5 fs. The interligand electron transfer (ILET) on the potential energy surface of the 1MLCT1 state is also extremely fast, with a rate constant of ca. 1.38 × 1013 s-1. The ultrafast ILET implies that the excited electron can dynamically delocalize over the three bpy ligands, despite the fact that the excited electron may be localized on either one of the three ligands at the equilibrium geometries of the three symmetric equivalent minima. Since rapid ILET essentially suggests delocalization, the long-standing controversy in inorganic photophysics-whether the excited electron is localized or delocalized-may therefore be calmed down to some extent.

Anthracene derivatives with strong spin-orbit coupling and efficient high-lying reverse intersystem crossing beyond the El-Sayed rule.

The attention to materials with hot exciton channel and triplet-triplet fusion (TTF) mediated high-lying reverse intersystem crossing (hRISC) has been raised for their ability to convert non-emissive 'dark' triplets into radiative singlet excitons. This spin conversion process results in high exciton utilization efficiency (EUE) that exceeds the theoretical limits. Notably, it is known that such spin conversion processes from the high-lying excited triplet to the singlet state are facilitated by the orthogonal orbital transition effect governed by the El-Sayed's rule. In this study, an anthracene derivative with indenoquinoline substituent 7,7-dimethyl-9-(10-(4-(naphthalen-1-yl)phenyl)anthracen-9-yl)-7H-indeno[1,2-f]quinoline (2MIQ-NPA) was synthesized and analyzed to investigate whether the hRISC process occurs in these molecules, even when the El-Sayed's rule is not followed. The hRISC channels of the emitter were fully unraveled through DFT calculations and experiments, which were quantitatively subdivided using transient electroluminescence measurements. The results showed that 2MIQ-NPA, which does not follow the El-Sayed's rule and has a relatively strong spin-orbit coupling matrix element of 0.116 cm-1 between the high-lying triplet state of T4 and the lowest singlet state of S1, effectively converted triplet excitons into singlet excitons with an EUE of 64.3%, contributed by a direct hot exciton channel of 19.2% and a TTF-mediated hot exciton channel of 15.1%. Despite the low outcoupling efficiency, the non-doped device with 2MIQ-NPA achieved an excellent device performance with an external quantum efficiency of 7.0%.

Computational design of boron-free triangular molecules with inverted singlet-triplet energy gap.

A novel, computationally designed, class of triangular-shape organic molecules with an inverted singlet-triplet (IST) energy gap is investigated with ab initio electronic structure methods. The considered molecular systems are cyclic oligomers and their common feature is electronic conjugation along the molecular rim. Vertical excitation energies from the electronic ground state to the lowest singlet and triplet excited states were computed, as well as vertical emission energies from these states to the ground state. The results underscore the significance of optimizing excited-state geometries to accurately describe the optoelectronic properties of IST molecules, in particular with respect to their application in OLEDs.

Axially-chiral boramidine for detailed (chir)optical studies.

The inclusion of boron atoms into chiral π-conjugated systems is an effective strategy to unlock unique chiroptical properties. Herein, the preparation and characterization of a configurationally stable axially-chiral boramidine are reported, showcasing absorption in the UV domain, deep-blue fluorescence (Φ up to 94%), and ca. |10-3| gabs and glum values. Detailed photophysical studies and quantum-chemical calculations clearly elucidate the deactivation pathways of the emissive state to triplet excited states, involving increased spin-orbit coupling between the lowest singlet excited state and an upper triplet state.

Accidental triplet harvesting in donor-acceptor dyads with low spin-orbit coupling.

We present an accidental mechanism for efficient intersystem crossing (ISC) between singlet and triplet states with low spin-orbit coupling (SOC) in molecules having donor-acceptor (D-A) moieties separated by a Sigma bond. Our study shows that SOC between the lowest singlet excited state and the higher-lying triplet states, together with nuclear motion-driven coupling of this triplet state with lower-lying triplet states during the free rotation about a Sigma bond, is one of the possible ways to achieve the experimentally observed ISC rate for a class of D-A type photoredox catalysts. This mechanism is found to be the dominant contributor to the ISC process with the corresponding rate reaching a maximum at a dihedral angle in the range of 72°-78° between the D-A moieties of 10-(naphthalen-1-yl)-3,7-diphenyl-10H-phenoxazine and other molecules included in the study. We have further demonstrated that the same mechanism is operative in a specific spirobis[anthracene]dione molecule, where the D and A moieties are interlocked near to the optimal dihedral angle, indicating the plausible effectiveness of the proposed mechanism. The present finding is expected to have implications in strategies for the synthesis of new generations of triplet-harvesting organic molecules.

Stimulating Phonon Bottleneck Effect in Organic Semiconductors by Charge-Transfer-Mediated J-Aggregation.

Hot carriers rapidly lose kinetic energies on a subpicosecond time scale, posing significant limitations on semiconductors' photon-conversion efficiencies. To slow the hot carrier cooling, the phonon bottleneck effect is constructed prevalently in quantum-confined structures with discrete energy levels. However, the maximum energy separation (ΔEES) between the energy levels is in a range of several hundred meV, leading to unsatisfactory cooling time. To address this, we design a novel organic semiconductor capable of forming intermolecular charge transfer (CT) in J-aggregates, where the lowest singlet excited state (S1) splits into two states due to the significant interplay between the Coulomb interaction and intermolecular CT coupling. The ΔEES between the two states can be adjusted up to 1.02 eV, and an extremely slow carrier cooling process of ∼72.3 ps was observed by femtosecond transient absorption spectroscopy. Moreover, the phonon bottleneck effect was identified in organic materials for the first time, and CT-mediated J-aggregation with short-range interactions was found to be the key to achieving large ΔEES. The significantly prolonged carrier cooling time, compared to <100 fs in the isolated molecule (10-6 M), highlights the potential of organic molecules with diversified aggregation structures in achieving long-lived hot carriers. These findings provide valuable insights into the intrinsic photophysics of electron-phonon scattering in organic semiconductors.

Bimolecular photoinduced symmetry-breaking charge separation of perylene in solution

Photoinduced symmetry-breaking charge separation (SB-CS) results in the generation of charge carriers through electron transfer between two identical molecules, after photoexcitation of one of them. It is usually studied in systems where the two reacting moieties are covalently linked. Examples of photoinduced bimolecular SB-CS with organic molecules yielding free ions remain scarce due to solubility or aggregation issues at the high concentrations needed to study this diffusion-assisted process. Here we investigate the excited-state dynamics of perylene (Pe) at high concentrations in solvents of varying polarity. Transient absorption spectroscopy on the subnanosecond to microsecond timescales reveal that self-quenching of Pe in the lowest singlet excited state leads to excimer formation in all solvents used. Additionally, bimolecular SB-CS, resulting in the generation of free ions, occurs concurrently to excimer formation in polar media, with a relative efficiency that increases with the polarity of the solvent. Moreover, we show that SB-CS is most efficient in room-temperature ionic liquids due to a charge-shielding effect leading to a larger escape of ions and due to the high viscosity that disfavours excimer formation.Graphical abstract

Roles of Singlet Fission in the Photosensitization of Silicon Phthalocyanine.

The roles of singlet fission in the triplet generation of silicon phthalocyanine (SiPc), a compound analogous to the IRDye700DX photosensitizer used in near-infrared photoimmunotherapy, are investigated by considering the energetical relation between the excitations of this compound. These excitations are obtained through spin-flip long-range corrected time-dependent density functional theory calculations. To initiate singlet fission, chromophores must meet two conditions: (1) near-degenerate low-lying singlet and quintet (triplet-triplet) excitations with a considerable energy gap of the lowest singlet and triplet excited states and (2) moderate π-stacking energy of chromophores, which is higher than but not far from the solvation energy, to facilitate the dissociation and generation of triplet-state chromophores. The present calculations demonstrate that SiPc satisfies both of these conditions after the formation of π-stacking irrespective of the presence of an axial ligand(s), suggesting that singlet fission plays a crucial role in the triplet generation process, although intersystem crossing occurs simultaneously at a very slow rate.

Anti-Kasha triplet energy transfer and excitation wavelength dependent persistent luminescence from host-guest doping systems

Anti-Kasha’s process in organic luminogens has attracted many attentions since its discovery. However, only limited examples of anti-Kasha’s rule have been reported and anti-Kasha triplet energy transfer (ET) is even less-touched. Benefiting from anti-Kasha’s rule, this work provided an efficient strategy to realize excitation wavelength dependent (Ex-De) afterglow in a host-guest system. The host has almost imperceptible RTP upon 365 nm excitation and guest is totally RTP inactive, while the doping host-guest system exhibits Ex-De afterglow with improved quantum yields. Anti-Kasha triplet ET process is realized from the higher excited triplet state T2 of host to the lowest excited singlet state S1 of the aggregated/unimolecular guest. ET efficiency in the doping system could be tuned by simply changing its processing methods to guide host and guest to adopt denser or looser intermolecular packing. The strategy of anti-Kasha triplet ET endows the host-guest doping system with multiple stimuli-responsive properties, including Ex-De afterglow, mechano-, and thermal-triggered afterglow behaviors. The corresponding applications of these properties are also realized in multiple information anti-counterfeiting and display.

Heptazine, Cyclazine, and Related Compounds: Chemically-Accurate Estimates of the Inverted Singlet-Triplet Gap.

Molecules that violate Hund's rule and exhibit an inverted gap between the lowest singlet S1 and triplet T1 excited states have attracted considerable attention due to their potential applications in optoelectronics. Among these molecules, the triangular-shaped heptazine, and its derivatives, have been in the limelight. However, conflicting reports have arisen regarding the relative energies of S1 and T1. Here, we employ highly accurate levels of theory, such as CC3, to not only resolve the debate concerning the sign but also quantify the magnitude of the S1-T1 gap. We also determined the 0-0 energies to evaluate the significance of the vertical approximation. In addition, we compute reference S1-T1 gaps for a series of 10 related molecules. This enables us to benchmark lower-order methods for future applications in larger systems within the same family of compounds. This contribution can serve as a foundation for the design of triangular-shaped molecules with enhanced photophysical properties.

Inverted Lowest Singlet and Triplet Excitation Energy Ordering of Graphitic Carbon Nitride Flakes.

In organic light-emitting diodes (OLEDs), only 25% of electrically generated excitons are in a singlet state, S1, and the remaining 75% are in a triplet state, T1. In thermally activated delayed fluorescence (TADF) chromophores the transition from the nonradiative T1 state to the radiative S1 state can be thermally activated, which improves the efficiency of OLEDs. Chromophores with inverted energy ordering of S1 and T1 states, S1 < T1, are superior to TADF chromophores, thanks to the absence of an energy barrier for the transition from T1 to S1. We benchmark the performance of time-dependent density functional theory using different exchange-correlation functionals and find that scaled long-range corrected double-hybrid functionals correctly predict the inverted singlet-triplet gaps of N-substituted phenalene derivatives. We then show that the inverted energy ordering of S1 and T1 is an intrinsic property of graphitic carbon nitride flakes. A design strategy of new chromophores with inverted singlet-triplet gaps is proposed. The color of emitted light can be fine-tuned through flake size and amine substitution on flake vertices.

Computational Design of Multiple Resonance-Type BN Molecules for Inverted Singlet and Triplet Excited States.

A computational design of linearly extended multiple resonance (MR)-type BN molecules based on DABNA-1 is proposed herein in the quest to find potential candidates that exhibit a negative singlet-triplet gap (ΔEST) and a large oscillator strength value. The impact of a proper account of the electron correlation in the lowest singlet and triplet excited states is systematically investigated by using double-hybrid functionals within the TD-DFT framework, as well as wavefunction-based methods (EOM-CCSD and SCS-CC2), since this contribution plays an essential role in driving the magnitude of the ΔEST in MR-TADF and inverted singlet-triplet gap compounds. Our results point out a gradual reduction of the ΔEST gap with respect to the increasing sum of the number of B and N atoms, reaching negative ΔEST values for some molecules as a function of their size. The double-hybrid functionals reproduce the gap with only slight deviation compared to available experimental data for DABNA-1, ν-DABNA, and mDBCz and nicely agree with high-level quantum mechanical methods (e.g., EOM-CCSD and SCS-CC2). Larger oscillator strengths are found compared to the azaphenalene-type molecules, also exhibiting the inversion of their singlet and triplet excited states. We hope this study can serve as a motivation for further design of the molecules showing negative ΔEST based on boron- and nitrogen-doped polyaromatic hydrocarbons.