1
|
Lashkaripour A, Park W, Mazaherifar M, Choi CH. Addressing the High-Throughput Screening Challenges of Inverted Singlet-Triplet Materials by MRSF-TDDFT. J Chem Theory Comput 2025; 21:5661-5668. [PMID: 40392574 DOI: 10.1021/acs.jctc.5c00451] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/22/2025]
Abstract
A new computational protocol utilizing mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT) and the DTCAM-STG exchange-correlation functional has been developed to identify materials with inverted singlet-triplet (INVEST) energy levels. This protocol surpasses existing quantum chemical methods in both accuracy and computational efficiency for predicting ΔEST, addressing challenges in high-throughput screening for INVEST materials. Based on this approach, novel heptazine derivatives have been proposed, which are anticipated to outperform currently known systems. Additionally, dynamic spin polarization (DSP) has been identified as a critical factor influencing INVEST phenomena. This insight provides a foundation for new design principles, enabling the discovery of materials with superior performance compared to existing alternatives.
Collapse
Affiliation(s)
- Alireza Lashkaripour
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Woojin Park
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Mohsen Mazaherifar
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Cheol Ho Choi
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| |
Collapse
|
2
|
Rishi V, Abou Taka A, Hratchian HP, McCaslin LM. Quantifying Design Principles for Light-Emitting Materials with Inverted Singlet-Triplet Energy Gaps. J Phys Chem Lett 2025:5213-5220. [PMID: 40387203 DOI: 10.1021/acs.jpclett.5c00827] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/20/2025]
Abstract
Molecular engineering of organic emitter molecules with inverted singlet-triplet energy gaps (INVEST) has emerged as a powerful approach for enhancing fluorescence efficiency through triplet harvesting. In these unique materials, the first excited singlet state (S1) lies below the lowest triplet state (T1), enabling efficient reverse intersystem crossing. Previous computational studies have focused on accurately calculating the inverted energy gap and establishing qualitative structure-property relationships. Here, we present quantitative relationships that link the molecular structure to the S1-T1 energy gap, ΔEST, by introducing a benchmark set of 15 heptazine-based INVEST molecules (HEPTA-INVEST15). We identify a strong linear correlation (R2 > 0.94) between ΔEST and both the degree of intramolecular charge transfer and the deviation from a single-excitation character, as quantified by %R1 values and transition density matrix norms. These trends persist across our expanded set of 44 mono-, di-, and tri-substituted heptazines (HEPTA-INVEST44), underscoring the generality of our findings. Notably, strongly electron-donating groups, such as -NH2, minimize the magnitude of inverted gaps in mono-substituted heptazines yet produce the most negative ΔEST in certain tri-substituted derivatives, a result arising from competing resonance effects and excited-state aromaticity. Although ΔEST shows no clear correlation with Hammett parameters, our results reveal that physically meaningful, computable descriptors offer a mechanistic foundation for the future data-driven design of INVEST emitters. These findings pave the way for machine-learning approaches that connect the molecular structure to ΔEST without requiring high-level excited-state calculations.
Collapse
Affiliation(s)
- Varun Rishi
- Sandia National Laboratories, Livermore, California 94550, United States
| | - Ali Abou Taka
- Sandia National Laboratories, Livermore, California 94550, United States
| | - Hrant P Hratchian
- University of California, Merced, Merced, California 95340, United States
| | - Laura M McCaslin
- Sandia National Laboratories, Livermore, California 94550, United States
| |
Collapse
|
3
|
Kumar H, Kumar P, Kumar A, Mahadevan A, Dhingra K, Kailasam K, Venkataramani S. Tripodal Heptazine Core-Based C 3-Symmetric Multi-Azo(Hetero)Arenes: Photoswitching, Supramolecular Sol-Gel Behavior, and Base Sensing Prospects. Chemistry 2025; 31:e202404757. [PMID: 40127109 DOI: 10.1002/chem.202404757] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/27/2024] [Revised: 03/21/2025] [Accepted: 03/24/2025] [Indexed: 03/26/2025]
Abstract
Electron-deficient heptazine-based molecular and polymeric materials are a fast-growing area of research with their exciting electronic and optical properties, finding potential applications in catalysis, optoelectronics etc. However, the continuous thirst for the design of heptazine-core-based molecular systems has grown in interest in the search for a new class of materials. Herein, we have synthesized two heptazine core-based C3-symmetric tripodal molecular systems decorated with three azobenzene and phenylazo-3,5-dimethylisoxazole as photoswitchable units connected through ─NH linker. The photoswitching characteristics and thermal half-lives of photoswitched states were studied using UV-vis and nuclear magnetic resonance (NMR) spectroscopy. In addition, for azobenzene derivative 1, intriguing properties such as reversible sonication-induced gelation and base sensing ability have been demonstrated and supported by computations.
Collapse
Affiliation(s)
- Himanshu Kumar
- Department of Chemical Sciences, Indian Institute of Science Education and Research (IISER) Mohali, Sector 81, SAS Nagar, Knowledge City, Manauli, Punjab, India
| | - Pravesh Kumar
- Advanced Functional Nanomaterials, Institute of Nano Science and Technology, (INST), Knowledge city, Sector 81, SAS Nagar, Manauli, Punjab, India
| | - Amit Kumar
- Advanced Functional Nanomaterials, Institute of Nano Science and Technology, (INST), Knowledge city, Sector 81, SAS Nagar, Manauli, Punjab, India
| | - Anjali Mahadevan
- Department of Chemical Sciences, Indian Institute of Science Education and Research (IISER) Mohali, Sector 81, SAS Nagar, Knowledge City, Manauli, Punjab, India
| | - Kirti Dhingra
- Advanced Functional Nanomaterials, Institute of Nano Science and Technology, (INST), Knowledge city, Sector 81, SAS Nagar, Manauli, Punjab, India
| | - Kamalakannan Kailasam
- Advanced Functional Nanomaterials, Institute of Nano Science and Technology, (INST), Knowledge city, Sector 81, SAS Nagar, Manauli, Punjab, India
| | - Sugumar Venkataramani
- Department of Chemical Sciences, Indian Institute of Science Education and Research (IISER) Mohali, Sector 81, SAS Nagar, Knowledge City, Manauli, Punjab, India
| |
Collapse
|
4
|
Pu YJ. Computational Design of Molecules Having Less Overlapping HOMO and LUMO in the Same Plane. J Phys Chem A 2025; 129:3466-3473. [PMID: 40179312 DOI: 10.1021/acs.jpca.5c00539] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/05/2025]
Abstract
The negative energy difference between singlet and triplet excited states (ΔEST) is currently attracting significant attention; however, molecular designs remain largely confined to azaphenalene structures, as reported by Leupin and Wirz in 1980. To show negative ΔEST, a maximally separated highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) arrangement is crucial to minimizing the exchange interaction in the excited state. We revisited the electronic structure of cyclazine, consisting of cyclododecahexaene ([12]annulene) and a central nitrogen atom. The 12 π-electrons of the peripheral cyclic oligoene play an important role in achieving the less overlapping HOMO and LUMO arrangement, and the bridging by the nitrogen atom inside produces the energy difference between HOMO and LUMO while maintaining a stable planar structure. Based on these insights, we designed a set of 10 molecules in which the number of π-electrons (N) in the peripheral cyclic oligoene is 16, 20, and 24, satisfying N = 4·n (n = 4, 5, 6), and a further set of 11 molecules in which N in the peripheral cyclic oligoene is extended to 14, 18, 22, and 26, satisfying N = 4·n + 2 (n = 3, 4, 5, 6). HOMO, LUMO, exchange interaction (K), and ΔEST were calculated using configuration interaction singles, TD-DFT, and equation of motion coupled-cluster singles and doubles (EOM-CCSD), with the structure optimized without any symmetry constraint. Among the molecular structures with N = 4·n, only the molecules without bond alternation exhibit less overlapping HOMO and LUMO and a small K and ΔEST. In contrast, among the molecular structures with N = 4·n + 2, none of the molecules exhibit less overlapping HOMO and LUMO arrangement. The molecules with both N = 4·n and no bond alternation show negative ΔEST in the EOM-CCSD calculation. The findings of this study will pave the way for broader molecular designs of molecules exhibiting negative ΔEST, where a less overlap of HOMO and LUMO is essential.
Collapse
Affiliation(s)
- Yong-Jin Pu
- RIKEN Center for Emergent Matter Science (CEMS), Wako, Saitama 351-0198, Japan
| |
Collapse
|
5
|
Wrigley L, Schlenker CW. Singlet-Triplet Inversion. Annu Rev Phys Chem 2025; 76:329-355. [PMID: 40258242 DOI: 10.1146/annurev-physchem-082423-120310] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/23/2025]
Abstract
The inversion of singlet and triplet states is a rare phenomenon, where, in opposition to Hund's first rule, singlet electronic states are stabilized relative to their triplet counterparts. The recent discovery of organic molecules exhibiting this inversion presents exciting new technological opportunities, such as addressing stability issues in organic light-emitting diodes (OLEDs). In this review, we describe fundamental molecular properties that can yield singlet-triplet inversion, generally ascribed to a phenomenon known as dynamic spin polarization. We discuss the systems in which singlet-triplet inversion was theoretically proposed, experimentally verified, and first implemented in an OLED device. We highlight key insights from the extensive computational work being carried out to understand the intricacies of these systems. Finally, we consider the outlook for future inverted singlet-triplet (IST) emitters.
Collapse
Affiliation(s)
- Liam Wrigley
- Department of Chemistry, University of Washington, Seattle, Washington, USA;
| | - Cody W Schlenker
- Department of Chemistry, University of Washington, Seattle, Washington, USA;
| |
Collapse
|
6
|
Kunze L, Hansen A, Grimme S, Mewes JM. The Best of Both Worlds: ΔDFT Describes Multiresonance TADF Emitters with Wave-Function Accuracy at Density-Functional Cost. J Phys Chem Lett 2025; 16:1114-1125. [PMID: 39846377 DOI: 10.1021/acs.jpclett.4c03192] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2025]
Abstract
With their narrow-band emission, high quantum yield, and good chemical stability, multiresonance thermally activated delayed fluorescence (MR-TADF) emitters are promising materials for OLED technology. However, accurately modeling key properties, such as the singlet-triplet (ST) energy gap and fluorescence energy, remains challenging. While time-dependent density functional theory (TD-DFT), the workhorse of computational materials science, suffers from fundamental issues, wave function-based coupled-cluster (CC) approaches, like approximate CC of second-order (CC2), are accurate but suffer from high computational cost and unfavorable scaling with system size. This work demonstrates that a state-specific ΔDFT approach based on unrestricted Kohn-Sham (ΔUKS) combines the best of both worlds: on a diverse benchmark set of 35 MR-TADF emitters, ΔUKS performs as good as or better than CC2, recovering experimental ST gaps with a mean absolute deviation (MAD) of 0.03 eV at a small fraction of the computational cost of CC2. When combined with a tuned range-separated LC-ωPBE functional, the excellent performance extends to fluorescence energies and ST gaps of MR- and donor-acceptor TADF emitters and even molecules with an inverted ST gap (INVEST), rendering this approach a jack of all trades for organic electronics.
Collapse
Affiliation(s)
- Lukas Kunze
- Mulliken Center for Theoretical Chemistry, Clausius Institute for Physical and Theoretical Chemistry, Rheinische Friedrich-Wilhelms Universität Bonn, Beringstraße 4, 53115 Bonn, Germany
| | - Andreas Hansen
- Mulliken Center for Theoretical Chemistry, Clausius Institute for Physical and Theoretical Chemistry, Rheinische Friedrich-Wilhelms Universität Bonn, Beringstraße 4, 53115 Bonn, Germany
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Clausius Institute for Physical and Theoretical Chemistry, Rheinische Friedrich-Wilhelms Universität Bonn, Beringstraße 4, 53115 Bonn, Germany
| | - Jan-Michael Mewes
- Mulliken Center for Theoretical Chemistry, Clausius Institute for Physical and Theoretical Chemistry, Rheinische Friedrich-Wilhelms Universität Bonn, Beringstraße 4, 53115 Bonn, Germany
- beeOLED GmbH, Niedersedlitzer Str. 75c, 01257 Dresden, Germany
| |
Collapse
|
7
|
Chanda S, Saha S, Sen S. Benchmark computations of nearly degenerate singlet and triplet states of N-heterocyclic chromophores. II. Density-based methods. J Chem Phys 2025; 162:024111. [PMID: 39774883 DOI: 10.1063/5.0238105] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/09/2024] [Accepted: 12/23/2024] [Indexed: 01/11/2025] Open
Abstract
In this paper, we demonstrate the performance of several density-based methods in predicting the inversion of S1 and T1 states of a few N-heterocyclic triangulene based fused ring molecules (popularly known as INVEST molecules) with an eye to identify a well performing but cost-effective preliminary screening method. Both conventional linear-response time-dependent density functional theory (LR-TDDFT) and ΔSCF methods (namely maximum overlap method, square-gradient minimization method, and restricted open-shell Kohn-Sham) are considered for excited state computations using exchange-correlation (XC) functionals from different rungs of Jacob's ladder. A well-justified systematism is observed in the performance of the functionals when compared against fully internally contracted multireference configuration interaction singles and doubles and/or equation of motion coupled-cluster singles and doubles (EOM-CCSD), with the most important feature being the capture of spin-polarization in the presence of correlation. A set of functionals with the least mean absolute error is proposed for both the approaches, LR-TDDFT and ΔSCF, which can be more cost-effective alternatives for computations on synthesizable larger derivatives of the templates studied here. We have based our findings on extensive studies of three cyclazine-based molecular templates, with additional studies on a set of six related templates. Previous benchmark studies for subsets of the functionals were conducted against the domain-based local pair natural orbital-similarity transformed EOM-CCSD (STEOM-CCSD), which resulted in an inadequate evaluation due to deficiencies in the benchmark theory. The role of exact-exchange, spin-contamination, and spin-polarization in the context of DFT comes to the forefront in our studies and supports the numerical evaluation of XC functionals for these applications. Suitable connections are drawn to two and three state exciton models, which identify the minimal physics governing the interactions in these molecules.
Collapse
Affiliation(s)
- Shamik Chanda
- Department of Chemical Sciences, Indian Institute of Science Education and Research (IISER) Kolkata, Nadia, Mohanpur 741246, WB, India
| | - Subhasish Saha
- Department of Chemical Sciences, Indian Institute of Science Education and Research (IISER) Kolkata, Nadia, Mohanpur 741246, WB, India
- MLR Institute of Technology, Dundigal, Hyderabad 500043, Telangana, India
| | - Sangita Sen
- Department of Chemical Sciences, Indian Institute of Science Education and Research (IISER) Kolkata, Nadia, Mohanpur 741246, WB, India
| |
Collapse
|
8
|
Valverde D, Ser CT, Ricci G, Jorner K, Pollice R, Aspuru-Guzik A, Olivier Y. Computational Investigations of the Detailed Mechanism of Reverse Intersystem Crossing in Inverted Singlet-Triplet Gap Molecules. ACS APPLIED MATERIALS & INTERFACES 2024; 16:66991-67001. [PMID: 38728616 DOI: 10.1021/acsami.4c04347] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/12/2024]
Abstract
Inverted singlet-triplet gap (INVEST) materials have promising photophysical properties for optoelectronic applications due to an inversion of their lowest singlet (S1) and triplet (T1) excited states. This results in an exothermic reverse intersystem crossing (rISC) process that potentially enhances triplet harvesting, compared to thermally activated delayed fluorescence (TADF) emitters with endothermic rISCs. However, the processes and phenomena that facilitate conversion between excited states for INVEST materials are underexplored. We investigate the complex potential energy surfaces (PESs) of the excited states of three heavily studied azaphenalene INVEST compounds, namely, cyclazine, pentazine, and heptazine using two state-of-the-art computational methodologies, namely, RMS-CASPT2 and SCS-ADC(2) methods. Our findings suggest that ISC and rISC processes take place directly between the S1 and T1 electronic states in all three compounds through a minimum-energy crossing point (MECP) with an activation energy barrier between 0.11 to 0.58 eV above the S1 state for ISC and between 0.06 and 0.36 eV above the T1 state for rISC. We predict that higher-lying triplet states are not populated, since the crossing point structures to these states are not energetically accessible. Furthermore, the conical intersection (CI) between the ground and S1 states is high in energy for all compounds (between 0.4 to 2.0 eV) which makes nonradiative decay back to the ground state a relatively slow process. We demonstrate that the spin-orbit coupling (SOC) driving the S1-T1 conversion is enhanced by vibronic coupling with higher-lying singlet and triplet states possessing vibrational modes of proper symmetry. We also rationalize that the experimentally observed anti-Kasha emission of cyclazine is due to the energetically inaccessible CI between the bright S2 and the dark S1 states, hindering internal conversion. Finally, we show that SCS-ADC(2) is able to qualitatively reproduce excited state features, but consistently overpredict relative energies of excited state structural minima compared to RMS-CASPT2. The identification of these excited state features elaborates design rules for new INVEST emitters with improved emission quantum yields.
Collapse
Affiliation(s)
- Danillo Valverde
- Laboratory for Computational Modeling of Functional Materials, Namur Institute of Structured Matter, Université de Namur, Rue de Bruxelles, 61, 5000 Namur, Belgium
| | - Cher Tian Ser
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto, 80 St. George St, Toronto, Ontario M5S 3H4, Canada
- Department of Computer Science, University of Toronto, 40 St. George St, Toronto, Ontario M5S 2E4, Canada
| | - Gaetano Ricci
- Laboratory for Computational Modeling of Functional Materials, Namur Institute of Structured Matter, Université de Namur, Rue de Bruxelles, 61, 5000 Namur, Belgium
| | - Kjell Jorner
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto, 80 St. George St, Toronto, Ontario M5S 3H4, Canada
- Department of Computer Science, University of Toronto, 40 St. George St, Toronto, Ontario M5S 2E4, Canada
- Department of Chemistry and Chemical Engineering, Chalmers University of Technology, 412 96 Kemigård, Sweden
| | - Robert Pollice
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto, 80 St. George St, Toronto, Ontario M5S 3H4, Canada
- Department of Computer Science, University of Toronto, 40 St. George St, Toronto, Ontario M5S 2E4, Canada
| | - Alán Aspuru-Guzik
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto, 80 St. George St, Toronto, Ontario M5S 3H4, Canada
- Department of Computer Science, University of Toronto, 40 St. George St, Toronto, Ontario M5S 2E4, Canada
- Department of Chemical Engineering & Applied Chemistry, University of Toronto, 200 College St., Ontario M5S 3E5, Canada
- Department of Materials Science & Engineering, University of Toronto, 184 College St., Ontario M5S 3E4, Canada
- Vector Institute for Artificial Intelligence, 661 University Ave., Suite 710, Toronto, Ontario M5G 1M1, Canada
- Acceleration Consortium, 700 University Ave., Toronto, Ontario M5G 1X6, Canada
| | - Yoann Olivier
- Laboratory for Computational Modeling of Functional Materials, Namur Institute of Structured Matter, Université de Namur, Rue de Bruxelles, 61, 5000 Namur, Belgium
| |
Collapse
|
9
|
Landi A, Ricci G, Olivier Y, Capobianco A, Peluso A. Toward Efficient Modeling of Nonradiative Decay in Extended INVEST: Overcoming Computational Challenges in Quantum Dynamics Simulations. J Phys Chem Lett 2024; 15:11042-11050. [PMID: 39470168 DOI: 10.1021/acs.jpclett.4c02713] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/30/2024]
Abstract
In recent years, an increasing number of fully organic molecules capable of thermally activated delayed fluorescence (TADF) have been reported, often with very small or even inverted singlet-triplet (INVEST) energy gaps. These molecules typically exhibit complex photophysics due to the close energy levels of multiple singlet and triplet states, which create various transition pathways toward emission. A predictive model for the rates of these transitions is thus essential for assessing the suitability of new materials for light-emitting devices. Quantum Dynamics (QD) calculations are ideal for this purpose, as they include quantum effects, without the limitations of first-order perturbative approaches, also allowing taking into account more than two electronic states at once. However, the huge computational demands of QD methodologies, especially for large molecules, currently limit their use as a standard tool. To address this problem, we here employ a strategy that allows us to include almost the whole set of the vibrational coordinates by selecting the key elements of the Hilbert space that significantly impact dynamics, thereby hugely reducing the computational burden. Application of this protocol to two relatively large INVEST molecules reveals that internal conversion in these systems is very fast, making indirect emissive pathways a possible channel for the population of the S1 state. More importantly, this study demonstrates that the dynamics can be accurately described even with a significantly reduced vibrational space, thus allowing quantum dynamics calculations that yield accurate transition rates in a few minutes of computational time.
Collapse
Affiliation(s)
- Alessandro Landi
- Dipartimento di Chimica e Biologia "Adolfo Zambelli", Università di Salerno, Via Giovanni Paolo II, I-84084 Fisciano (SA), Italy
| | - Gaetano Ricci
- Laboratory for Computational Modeling of Functional Materials, Namur Institute of Structured Matter, Universitè de Namur, Rue de Bruxelles 61, B-5000 Namur, Belgium
| | - Yoann Olivier
- Laboratory for Computational Modeling of Functional Materials, Namur Institute of Structured Matter, Universitè de Namur, Rue de Bruxelles 61, B-5000 Namur, Belgium
| | - Amedeo Capobianco
- Dipartimento di Chimica e Biologia "Adolfo Zambelli", Università di Salerno, Via Giovanni Paolo II, I-84084 Fisciano (SA), Italy
| | - Andrea Peluso
- Dipartimento di Chimica e Biologia "Adolfo Zambelli", Università di Salerno, Via Giovanni Paolo II, I-84084 Fisciano (SA), Italy
| |
Collapse
|
10
|
Chanda S, Sen S. Benchmark computations of nearly degenerate singlet and triplet states of N-heterocyclic chromophores. I. Wavefunction-based methods. J Chem Phys 2024; 161:174117. [PMID: 39503472 DOI: 10.1063/5.0225537] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/25/2024] [Accepted: 08/07/2024] [Indexed: 11/08/2024] Open
Abstract
In this paper, we investigate the role of electron correlation in predicting the S1-S0 and T1-S0 excitation energies and, hence, the singlet-triplet gap (ΔEST) in a set of cyclazines, which act as templates for potential candidates for fifth generation organic light emitting diode materials. This issue has recently garnered much interest with the focus being on the inversion of the ΔEST, although experiments have indicated near degenerate levels with both positive and negative being within the experimental error bar [J. Am. Chem. Soc. 102, 6068 (1980), J. Am. Chem. Soc. 108, 17(1986)]. We have carried out a systematic and exhaustive study of various excited state electronic structure methodologies and identified the strengths and shortcomings of the various approaches and approximations in view of this challenging case. We have found that near degeneracy can be achieved either with a proper balance of static and dynamic correlation in multireference theories or with state-specific orbital corrections, including its coupling with correlation. The role of spin contamination is also discussed. Eventually, this paper seeks to produce benchmark numbers for establishing cost-effective theories, which can then be used for screening derivatives of these templates with desirable optical and structural properties. Additionally, we would like to point out that the use of domain-based local pair natural orbital-similarity transformed EOM-coupled cluster singles and doubles as the benchmark for ΔEST [as used in J. Phys. Chem. A 126(8), 1378 (2022), Chem. Phys. Lett. 779, 138827 (2021)] is not a suitable benchmark for these classes of molecules.
Collapse
Affiliation(s)
- Shamik Chanda
- Department of Chemical Sciences, Indian Institute of Science Education and Research (IISER) Kolkata Nadia, Mohanpur 741246, West Bengal, India
| | - Sangita Sen
- Department of Chemical Sciences, Indian Institute of Science Education and Research (IISER) Kolkata Nadia, Mohanpur 741246, West Bengal, India
| |
Collapse
|
11
|
Majumdar A, Jindal K, Das S, Ramakrishnan R. Influence of pseudo-Jahn-Teller activity on the singlet-triplet gap of azaphenalenes. Phys Chem Chem Phys 2024; 26:26723-26733. [PMID: 39403735 DOI: 10.1039/d4cp02761b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/01/2024]
Abstract
We analyze the possibility of symmetry-lowering induced by pseudo-Jahn-Teller interactions in six previously studied azaphenalenes that are known to have their first excited singlet state (S1) lower in energy than the triplet state (T1). The primary aim of this study is to explore whether Hund's rule violation is observed in these molecules when their structures are distorted from C2v or D3h point group symmetries by vibronic coupling. Along two interatomic distances connecting these point groups to their subgroups Cs or C3h, we relaxed the other internal degrees of freedom and calculated two-dimensional potential energy subsurfaces. The many-body perturbation theory (MP2) suggests that the high-symmetry structures are the energy minima for all six systems. However, single-point energy calculations using the coupled-cluster method (CCSD(T)) indicate symmetry lowering in four cases. The singlet-triplet energy gap plotted on the potential energy surface also shows variations when deviating from high-symmetry structures. A full geometry optimization at the CCSD(T) level with the cc-pVTZ basis set reveals that the D3h structure of cyclazine (1AP) is a saddle point, connecting two equivalent minima of C3h symmetry undergoing rapid automerization. The combined effects of symmetry lowering and high-level corrections result in a nearly zero singlet-triplet gap for the C3h structure of cyclazine. Azaphenalenes containing nitrogen atoms at electron-deficient sites - 2AP, 3AP, and 4AP - exhibit more pronounced in-plane structural distortion; the effect is captured by the long-range exchange-interaction corrected DFT method, ωB97XD. Excited state calculations of these systems indicate that in their low-symmetry energy minima, T1 is indeed lower in energy than S1, upholding the validity of Hund's rule. Jahn-Teller analysis predicts the symmetries of the in-plane distortion vibrational modes as or B2: C2v → Cs agreeing with the vibrational frequencies of the saddle-points.
Collapse
Affiliation(s)
- Atreyee Majumdar
- Tata Institute of Fundamental Research, Hyderabad 500046, India.
| | - Komal Jindal
- Tata Institute of Fundamental Research, Hyderabad 500046, India.
| | - Surajit Das
- Tata Institute of Fundamental Research, Hyderabad 500046, India.
| | | |
Collapse
|
12
|
Michalak P, Lesiuk M. Rank-Reduced Equation-of-Motion Coupled Cluster Triples: an Accurate and Affordable Way of Calculating Electronic Excitation Energies. J Chem Theory Comput 2024; 20:8970-8983. [PMID: 39347964 PMCID: PMC11500429 DOI: 10.1021/acs.jctc.4c00959] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2024] [Revised: 09/09/2024] [Accepted: 09/11/2024] [Indexed: 10/01/2024]
Abstract
In the present work, we report an implementation of the rank-reduced equation-of-motion coupled cluster method with approximate triple excitations (RR-EOM-CC3). The proposed variant relies on tensor decomposition techniques in order to alleviate the high cost of computing and manipulating the triply excited amplitudes. In the RR-EOM-CC3 method, both ground-state and excited-state triple-excitation amplitudes are compressed according to the Tucker-3 format. This enables factorization of the working equations such that the formal scaling of the method is reduced to N6, where N is the system size. An additional advantage of our method is the fact that the accuracy can be strictly controlled by proper choice of two parameters defining sizes of triple-excitation subspaces in the Tucker decomposition for the ground and excited states. Optimal strategies of selecting these parameters are discussed. The developed method has been tested in a series of calculations of electronic excitation energies and compared to its canonical EOM-CC3 counterpart. Errors several times smaller than the inherent error of the canonical EOM-CC3 method (in comparison to FCI) are straightforward to achieve. This conclusion holds both for valence states dominated by single excitations and for states with pronounced doubly excited character. Taking advantage of the decreased scaling, we demonstrate substantial computational costs reductions (in comparison with the canonical EOM-CC3) in the case of two large molecules - l-proline and heptazine. This illustrates the usefulness of the RR-EOM-CC3 method for accurate determination of excitation energies of large molecules.
Collapse
Affiliation(s)
- Piotr Michalak
- Faculty of Chemistry, University of Warsaw, Pasteura 1, Warsaw 02-093, Poland
| | - Michał Lesiuk
- Faculty of Chemistry, University of Warsaw, Pasteura 1, Warsaw 02-093, Poland
| |
Collapse
|
13
|
Bedogni M, Di Maiolo F. Singlet-Triplet Inversion in Triangular Boron Carbon Nitrides. J Chem Theory Comput 2024; 20:8634-8643. [PMID: 39264103 DOI: 10.1021/acs.jctc.4c00706] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/13/2024]
Abstract
The discovery of singlet-triplet (ST) inversion in some π-conjugated triangle-shaped boron carbon nitrides is a remarkable breakthrough that defies Hund's first rule. Deeply rooted in strong electron-electron interactions, ST inversion has garnered significant interest due to its potential to revolutionize triplet harvesting in organic LEDs. Using the well-established Pariser-Parr-Pople model for correlated electrons in π-conjugated systems, we employ a combination of CISDT and restricted active space configuration interaction calculations to investigate the photophysics of several triangular boron carbon nitrides. Our findings reveal that ST inversion in these systems is primarily driven by a network of alternating electron-donor and electron-acceptor groups in the molecular rim, rather than by the triangular molecular structure itself.
Collapse
Affiliation(s)
- Matteo Bedogni
- Department of Chemistry, Life Science and Environmental Sustainability, Università di Parma, 43124 Parma, Italy
| | - Francesco Di Maiolo
- Department of Chemistry, Life Science and Environmental Sustainability, Università di Parma, 43124 Parma, Italy
| |
Collapse
|
14
|
Knysh I, Lipparini F, Blondel A, Duchemin I, Blase X, Loos PF, Jacquemin D. Reference CC3 Excitation Energies for Organic Chromophores: Benchmarking TD-DFT, BSE/ GW, and Wave Function Methods. J Chem Theory Comput 2024. [PMID: 39237472 DOI: 10.1021/acs.jctc.4c00906] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/07/2024]
Abstract
To expand the QUEST database of highly accurate vertical transition energies, we consider a series of large organic chromogens ubiquitous in dye chemistry, such as anthraquinone, azobenzene, BODIPY, and naphthalimide. We compute, at the CC3 level of theory, the singlet and triplet vertical transition energies associated with the low-lying excited states. This leads to a collection of more than 120 new highly accurate excitation energies. For several singlet transitions, we have been able to determine CCSDT transition energies with a compact basis set, finding minimal deviations from the CC3 values for most states. Subsequently, we employ these reference values to benchmark a series of lower-order wave function approaches, including the popular ADC(2) and CC2 schemes, as well as time-dependent density-functional theory (TD-DFT), both with and without applying the Tamm-Dancoff approximation (TDA). At the TD-DFT level, we evaluate a large panel of global, range-separated, local, and double hybrid functionals. Additionally, we assess the performance of the Bethe-Salpeter equation (BSE) formalism relying on both G0W0 and evGW quasiparticle energies evaluated from various starting points. It turns out that CC2 and ADC(2.5) are the most accurate models among those with respective O ( N 5 ) and O ( N 6 ) scalings with system size. In contrast, CCSD does not outperform CC2. The best performing exchange-correlation functionals include BMK, M06-2X, M06-SX, CAM-B3LYP, ωB97X-D, and LH20t, with average deviations of approximately 0.20 eV or slightly below. Errors on vertical excitation energies can be further reduced by considering double hybrids. Both SOS-ωB88PP86 and SOS-ωPBEPP86 exhibit particularly attractive performances with overall quality on par with CC2, whereas PBE0-DH and PBE-QIDH are only slightly less efficient. BSE/evGW calculations based on Kohn-Sham starting points have been found to be particularly effective for singlet transitions, but much less for their triplet counterparts.
Collapse
Affiliation(s)
- Iryna Knysh
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, Via Moruzzi 3, 56124 Pisa, Italy
| | - Aymeric Blondel
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| | - Ivan Duchemin
- Université Grenoble Alpes, CEA, IRIG-MEM-L Sim, 38054 Grenoble, France
| | - Xavier Blase
- Université Grenoble Alpes, CNRS, Institut NEEL, F-38042 Grenoble, France
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques, Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Denis Jacquemin
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
- Institut Universitaire de France (IUF), F-75005 Paris, France
| |
Collapse
|
15
|
Okumura R, Tanaka H, Shizu K, Fukushima S, Yasuda Y, Kaji H. Development of an Organic Emitter Exhibiting Reverse Intersystem Crossing Faster than Intersystem Crossing. Angew Chem Int Ed Engl 2024; 63:e202409670. [PMID: 38943493 DOI: 10.1002/anie.202409670] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/22/2024] [Revised: 06/25/2024] [Accepted: 06/27/2024] [Indexed: 07/01/2024]
Abstract
In thermally activated delayed fluorescence (TADF)-based organic light-emitting diodes (OLEDs), acceleration of reverse intersystem crossing (RISC) and suppression of intersystem crossing (ISC) are demanded to shorten a lifetime of triplet excitons. As a system realizing RISC faster than ISC, inverted singlet-triplet excited states (iST) with a negative energy difference (ΔEST) between the lowest excited singlet and the lowest triplet states have been gathering much attention recently. Here, we have focused on an asymmetric hexa-azaphenalene (A6AP) core to obtain a new insight into iST. Based on A6AP, we have newly designed A6AP-Cz with the calculated ΔEST of -44 meV. The experimental studies of a synthesized A6AP-Cz revealed that the lifetime of delayed fluorescence (τDF) was only 54 ns, which was the shortest among all organic materials. The rate constant of RISC (kRISC=1.9×107 s-1) was greater than that of ISC (kISC=1.0×107 s-1). The negative ΔEST of A6AP-Cz was experimentally confirmed from 1) the kRISC and kISC (-45 meV) and 2) the temperature-dependent τDF. 3) The onsets of fluorescence and phosphorescence spectra at 77 K also supported the evidence of negative ΔEST (-73 meV). This study demonstrated the potential of A6AP as an iST core for the first time.
Collapse
Affiliation(s)
- Ryosuke Okumura
- Institute for Chemical Research, Kyoto University, Uji, Kyoto, 611-0011, Japan
| | - Hiroyuki Tanaka
- Institute for Chemical Research, Kyoto University, Uji, Kyoto, 611-0011, Japan
| | - Katsuyuki Shizu
- Institute for Chemical Research, Kyoto University, Uji, Kyoto, 611-0011, Japan
| | - Shohei Fukushima
- Institute for Chemical Research, Kyoto University, Uji, Kyoto, 611-0011, Japan
| | - Yuka Yasuda
- Institute for Chemical Research, Kyoto University, Uji, Kyoto, 611-0011, Japan
| | - Hironori Kaji
- Institute for Chemical Research, Kyoto University, Uji, Kyoto, 611-0011, Japan
| |
Collapse
|
16
|
Karak P, Manna P, Banerjee A, Ruud K, Chakrabarti S. Reverse Intersystem Crossing Dynamics in Vibronically Modulated Inverted Singlet-Triplet Gap System: A Wigner Phase Space Study. J Phys Chem Lett 2024; 15:7603-7609. [PMID: 39028946 DOI: 10.1021/acs.jpclett.4c01437] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/21/2024]
Abstract
We inspect the origin of the inverted singlet-triplet gap (INVEST) and slow change in the reverse intersystem crossing (rISC) rate with temperature, as recently observed. A Wigner phase space study reveals that, though INVEST is found at equilibrium geometry, variation in the exchange interaction and the doubles-excitation for other geometries in the harmonic region leads to non-INVEST behavior. This highlights the importance of nuclear degrees of freedom for the INVEST phenomenon, and in this case, geometric puckering of the studied molecule determines INVEST and the associated rISC dynamics.
Collapse
Affiliation(s)
- Pijush Karak
- Department of Chemistry, University of Calcutta, 92 A. P. C. Road, Kolkata 700009, West Bengal, India
| | - Pradipta Manna
- Department of Chemistry, University of Calcutta, 92 A. P. C. Road, Kolkata 700009, West Bengal, India
| | - Ambar Banerjee
- Research Institute for Sustainable Energy (RISE), TCG Centres for Research and Education in Science and Technology (TCG-CREST), Kolkata 700091, India
| | - Kenneth Ruud
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Tromsø - The Arctic University of Norway, 9037 Tromsø, Norway
- Norwegian Defence Research Establishment, P.O. Box 25, 2027 Kjeller, Norway
| | - Swapan Chakrabarti
- Department of Chemistry, University of Calcutta, 92 A. P. C. Road, Kolkata 700009, West Bengal, India
| |
Collapse
|
17
|
Kunze L, Froitzheim T, Hansen A, Grimme S, Mewes JM. ΔDFT Predicts Inverted Singlet-Triplet Gaps with Chemical Accuracy at a Fraction of the Cost of Wave Function-Based Approaches. J Phys Chem Lett 2024:8065-8077. [PMID: 39083761 DOI: 10.1021/acs.jpclett.4c01649] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/02/2024]
Abstract
Efficient OLEDs need to quickly convert singlet and triplet excitons into photons. Molecules with an inverted singlet-triplet energy gap (INVEST) are promising candidates for this task. However, typical INVEST molecules have drawbacks like too low oscillator strengths and excitation energies. High-throughput screening could identify suitable INVEST molecules, but existing methods are problematic: The workhorse method TD-DFT cannot reproduce gap inversion, while wave function-based methods are too slow. This study proposes a state-specific method based on unrestricted Kohn-Sham DFT with common hybrid functionals. Tuned on the new INVEST15 benchmark set, this method achieves an error of less than 1 kcal/mol, which is traced back to error cancellation between spin contamination and dynamic correlation. Applied to the larger and structurally diverse NAH159 set in a black-box fashion, the method maintains a small error (1.2 kcal/mol) and accurately predicts gap signs in 83% of cases, confirming its robustness and suitability for screening workflows.
Collapse
Affiliation(s)
- Lukas Kunze
- Mulliken Center for Theoretical Chemistry, Clausius Institute for Physical and Theoretical Chemistry, Rheinische Friedrich-Wilhelms Universität Bonn, Beringstraße 4, 53115 Bonn, Germany
| | - Thomas Froitzheim
- Mulliken Center for Theoretical Chemistry, Clausius Institute for Physical and Theoretical Chemistry, Rheinische Friedrich-Wilhelms Universität Bonn, Beringstraße 4, 53115 Bonn, Germany
| | - Andreas Hansen
- Mulliken Center for Theoretical Chemistry, Clausius Institute for Physical and Theoretical Chemistry, Rheinische Friedrich-Wilhelms Universität Bonn, Beringstraße 4, 53115 Bonn, Germany
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Clausius Institute for Physical and Theoretical Chemistry, Rheinische Friedrich-Wilhelms Universität Bonn, Beringstraße 4, 53115 Bonn, Germany
| | - Jan-Michael Mewes
- Mulliken Center for Theoretical Chemistry, Clausius Institute for Physical and Theoretical Chemistry, Rheinische Friedrich-Wilhelms Universität Bonn, Beringstraße 4, 53115 Bonn, Germany
- beeOLED GmbH, Niedersedlitzer Str. 75c, 01257 Dresden, Germany
| |
Collapse
|
18
|
Duszka MW, Rode MF, Sobolewski AL. Computational design of boron-free triangular molecules with inverted singlet-triplet energy gap. Phys Chem Chem Phys 2024; 26:19130-19137. [PMID: 38973243 DOI: 10.1039/d4cp01658k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/09/2024]
Abstract
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.
Collapse
Affiliation(s)
| | - Michał F Rode
- Institute of Physics, Polish Academy of Sciences, Warsaw, Poland.
| | | |
Collapse
|
19
|
Majumdar A, Ramakrishnan R. Resilience of Hund's rule in the chemical space of small organic molecules. Phys Chem Chem Phys 2024; 26:14505-14513. [PMID: 38741560 DOI: 10.1039/d4cp00886c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/16/2024]
Abstract
We embark on a quest to identify small molecules in the chemical space that can potentially violate Hund's rule. Utilizing twelve TDDFT approximations and the ADC(2) many-body method, we report the energies of S1 and T1 excited states of 12 880 closed-shell organic molecules within the bigQM7ω dataset with up to 7 CONF atoms. In this comprehensive dataset, none of the molecules, in their minimum energy geometry, exhibit a negative S1-T1 energy gap at the ADC(2) level while several molecules display values <0.1 eV. The spin-component-scaled double-hybrid method, SCS-PBE-QIDH, demonstrates the best agreement with ADC(2). Yet, at this level, a few molecules with a strained sp3-N center turn out as false-positives with the S1 state lower in energy than T1. We investigate a prototypical cage molecule with an energy gap <-0.2 eV, which a closer examination revealed as another false positive. We conclude that in the chemical space of small closed-shell organic molecules, it is possible to identify geometric and electronic structural features giving rise to S1-T1 degeneracy; still, there is no evidence of a negative gap. We share the dataset generated for this study as a module, to facilitate seamless molecular discovery through data mining.
Collapse
Affiliation(s)
- Atreyee Majumdar
- Tata Institute of Fundamental Research, Hyderabad 500046, India.
| | | |
Collapse
|
20
|
Jorner K, Pollice R, Lavigne C, Aspuru-Guzik A. Ultrafast Computational Screening of Molecules with Inverted Singlet-Triplet Energy Gaps Using the Pariser-Parr-Pople Semiempirical Quantum Chemistry Method. J Phys Chem A 2024; 128:2445-2456. [PMID: 38485448 PMCID: PMC10983003 DOI: 10.1021/acs.jpca.3c06357] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2023] [Revised: 02/23/2024] [Accepted: 02/23/2024] [Indexed: 04/04/2024]
Abstract
Molecules with an inverted energy gap between their first singlet and triplet excited states have promising applications in the next generation of organic light-emitting diode (OLED) materials. Unfortunately, such molecules are rare, and only a handful of examples are currently known. High-throughput virtual screening could assist in finding novel classes of these molecules, but current efforts are hampered by the high computational cost of the required quantum chemical methods. We present a method based on the semiempirical Pariser-Parr-Pople theory augmented by perturbation theory and show that it reproduces inverted gaps at a fraction of the cost of currently employed excited-state calculations. Our study paves the way for ultrahigh-throughput virtual screening and inverse design to accelerate the discovery and development of this new generation of OLED materials.
Collapse
Affiliation(s)
- Kjell Jorner
- Institute
of Chemical and Bioengineering, Department of Chemistry and Applied
Biosciences, ETH Zürich, Vladimir-Prelog-Weg 1, Zürich CH-8093, Switzerland
- Department
of Chemistry and Chemical Engineering, Chalmers
University of Technology, Kemigården 4, Gothenburg SE-41258, Sweden
- Chemical
Physics Theory Group, Department of Chemistry, University of Toronto, 80 St. George Street, Toronto M5S 3H6, Canada
- Department
of Computer Science, University of Toronto, 40 St. George Street, Toronto M5S 2E4, Canada
| | - Robert Pollice
- Chemical
Physics Theory Group, Department of Chemistry, University of Toronto, 80 St. George Street, Toronto M5S 3H6, Canada
- Department
of Computer Science, University of Toronto, 40 St. George Street, Toronto M5S 2E4, Canada
- Stratingh
Institute for Chemistry, University of Groningen, Nijenborgh 4, Groningen 9747, AG, The Netherlands
| | - Cyrille Lavigne
- Chemical
Physics Theory Group, Department of Chemistry, University of Toronto, 80 St. George Street, Toronto M5S 3H6, Canada
- Department
of Computer Science, University of Toronto, 40 St. George Street, Toronto M5S 2E4, Canada
| | - Alán Aspuru-Guzik
- Chemical
Physics Theory Group, Department of Chemistry, University of Toronto, 80 St. George Street, Toronto M5S 3H6, Canada
- Department
of Computer Science, University of Toronto, 40 St. George Street, Toronto M5S 2E4, Canada
- Department
of Chemical Engineering & Applied Chemistry, University of Toronto, 200 College Street, Toronto M5S 3E5, Canada
- Department
of Materials Science & Engineering, University of Toronto, 184 College Street, Toronto M5S 3E4, Canada
- Vector
Institute for Artificial Intelligence, 661 University Ave. Suite 710, Toronto M5G 1M1, Canada
- Lebovic
Fellow, Canadian Institute for Advanced
Research (CIFAR), 661
University Avenue, Toronto M5G 1M1, Canada
- Acceleration
Consortium, University of Toronto, 700 University Avenue, Toronto M5G 1Z5, Canada
| |
Collapse
|
21
|
Loos PF, Lipparini F, Jacquemin D. Heptazine, Cyclazine, and Related Compounds: Chemically-Accurate Estimates of the Inverted Singlet-Triplet Gap. J Phys Chem Lett 2023; 14:11069-11075. [PMID: 38048474 DOI: 10.1021/acs.jpclett.3c03042] [Citation(s) in RCA: 29] [Impact Index Per Article: 14.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/06/2023]
Abstract
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.
Collapse
Affiliation(s)
- Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques, Université de Toulouse, CNRS, UPS, 31400 Toulouse, France
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, Via Moruzzi 3, 56124 Pisa, Italy
| | - Denis Jacquemin
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
- Institut Universitaire de France, 75005 Paris, France
| |
Collapse
|
22
|
Wang X, Wang A, Zhao M, Marom N. Inverted Lowest Singlet and Triplet Excitation Energy Ordering of Graphitic Carbon Nitride Flakes. J Phys Chem Lett 2023:10910-10919. [PMID: 38033187 DOI: 10.1021/acs.jpclett.3c02835] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/02/2023]
Abstract
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.
Collapse
Affiliation(s)
- Xiaopeng Wang
- School of Foundational Education, University of Health and Rehabilitation Sciences, Qingdao 266114, China
- Qingdao Institute for Theoretical and Computational Sciences, Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao, Shandong 266237, P. R. China
| | - Aizhu Wang
- Collaborative Innovation Center of Technology and Equipment for Biological Diagnosis and Therapy in Universities of Shandong, Institute for Advanced Interdisciplinary Research (iAIR), University of Jinan, Jinan, 250022, P. R. China
| | - Mingwen Zhao
- School of Physics and State Key Laboratory of Crystal Materials, Shandong University, Jinan, Shandong 250100, China
| | - Noa Marom
- Department of Materials Science and Engineering, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, United States
- Department of Physics, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, United States
- Department of Chemistry, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, United States
| |
Collapse
|