Geometries and FMO analysis in S0 ground state
The optimized ground state geometries of the eleven molecules 1a-d, 2a-d and 3a-c indicate large torsion angles between the acceptor and the donor fragments. The dihedral angles j between the DMA unit and the connected phenyl are evaluated to be approximately 82° for 1a-1d, 80° for 2a-2d and 87° for 3a-3c. As depicted in Figure S1, the electron density distribution of the frontier molecular orbitals (FMOs) for all derivatives 1a-d, 2a-d and 3a-c in the S0 state shows a p-character. It is noted that the lowest unoccupied molecular orbital (LUMO) is located mainly on the diphenyl-diazine core, whereas the highest occupied molecular orbital (HOMO) concentrated on the dimethyl-acridine donor fragment. As a result, with just a very small spatial orbital overlap, the HOMO and LUMO are almost completely separated which is demonstrated by the compositions (%) of the HOMOs and LUMOs (Table S1). In Table 1, the low calculated values of donor/acceptor compositions (D/A) ratio of LUMO in the S0 state are also provided to indicate the effective HOMO-LUMO separation in the all designed molecules. This great spatial separation made possible by the orthogonal geometry between the dimethylacridine donor units and the biphenyl-diazine accepting core, enabling thermally activated delayed fluorescence.
As depicted in Figure 1, the calculated electrostatic potential maps indicate clearly that the negative region (red color) is mainly localized on the two nitrogen atoms of the diazine core. From Table 1, the LUMO energy levels of 1a-d, 2a-d and 3a-c are varied from -1.69 to -2.10 eV, while the HOMO energies vary between -1.69 and -2.10 eV. Furthermore, the studied molecules display a large energy gap DEH-L ranging from 2.83 to 3.35 eV, which allows to expect blue light emission [25,26].
<Insert Table 1 here>
<Insert Figure 1 here>
Charge injection and transport properties
The reorganization energies for electron (λe) and hole (λh) of the studied molecules were predicted from the single point energy calculations at the B3LYP/6-31G(d,p) level based on the B3LYP/6-31G(d,p) optimized neutral, cationic, and anionic geometries. The computed values of reorganization energies are summarized in Table 1. The nature of diazine core and the donor linking position seem to have significant effect on the reorganization energy for both electrons (0.29-0.61eV) and holes (0.12-0.21eV). The obtained results given in Table 1 show that the λh value is the lowest when the two phenyl bridges are linked at the 3,5-positions or 3,6-positions of pyridazine (2a and 3a) and 2,5-positions of pyrazine (3c), whereas the λe value is the lowest when the two phenyl units are linked at the 4,6-positions of pyrimidine (2b). Furthermore, the reorganization energies for hole transport (λh) of all the studied compounds are comparatively smaller than the reorganization energies for electron transport (λe), which reveals that the hole transport performance of these compounds is better than the electron transport ability. Moreover, the calculated λh values of 1a-3c are smaller than that of N,N’-diphenyl-N,N’-bis(3-methlphenyl)-(1,1’-biphenyl)-4,4’-diamine (TPD), which is a typical hole transport material (λh = 0.292 eV) [55], showing relevance of the studied derivatives as p-type material for organic light emitting devices. According to Table 1, the IP value is the lowest and the EA value is the highest when the two phenyl units are linked at 2,5-positions of pyrazine (3c). Thus, the molecule 3c would have a greater capacity for both electron or hole injection.
Exciton binding energies
The computed values of exciton binding energies of the investigated molecules are given in Table 2. The all investigated molecules 1a-d, 2a-d, 3a-c have appreciable exciton binding energies ranging from 0.39 to 0.54 eV, due to the effective attraction in the electron-hole pair, suggesting that when used in OLEDs, the excitons are formed regardless of the diazine kind and regardless of the linking positions of the DMA units.
<Insert Table 2 here>
Singlet-triplet energy gap
The S1 and T1 energy levels of the studied materials were calculated using TD-DFT at B3LYP/6-31G(d,p) level of theory (See Table 2). The eleven studied molecules exhibit high S1 energy levels above 2.47 eV, indicating that they are blue materials and their T1 energy levels are calculated to be in the range 2.38-2.82 eV. The calculated ΔEST gaps for the investigated molecules range from 4.2 to 135.1 eV indicating that they exhibit TADF properties and demonstrating that the singlet-triplet energy gap is significantly controlled by adjusting the linking position of the two dimethyl-acridine electron-donor units with the diphenyldiazine acceptor core. From Table 2, it turns out that, among the designed molecules, 1c (where the phenyl bridges are linked at the 4,5-positions of pyrimidine), 2b and 2c (where the phenyl bridges are linked at the 4,6- and 2,3-positions of pyrimidine) and 3b (where the phenyl bridges are linked at the 2,5-positions of pyrimidine) are characterized by the smaller DEST values which suggests more efficient reverse intersystem crossing and consequently results in thermally activated delayed fluorescence.
Emission spectra
The values of the variation of dipole moment with respect to ground state for the first singlet excited S1 state (DmS1-S0) of the studied molecules are listed in Table 2. It can be seen that all molecules exhibit significant change in dipole moment between fundamental and first excited states which is equivalent to important geometric distortion in the S1 state [56]. The large DmS1-S0 values suggest that the emission of these compounds originates from a state which is much polar than the ground-state [57]. Further, the large values of relaxation energies (λS1−S0~2.0 eV) in the singlet excited state obtained for the all studied compounds and listed in Table 2, suggest significant geometrical deformations in the S1 state and indicate that compounds 2c and 3c with lowest λS1−S0 values possess slower radiative rates and thus extended excited state lifetimes. Indeed, the rate of the non-radiative internal conversion from the S1 to S0 is directly correlated with the relaxation energy [58], which quantifies the strength of the electron-vibration interactions in the S1 state. According to the Einstein spontaneous emission relationship: [59], typically large λS1-S0 values correspond to large rates for the internal conversion from the S1 to S0 state. On the other hand, from the values of the emission wavelength from the singlet state listed in Table 2, it can be seen that the all molecules show light emission in the region from 480 nm to 540 nm in toluene solvent and can be considered as potential TADF blue emitters.
Aromaticity of the diazine core
In order to examine the effect of aromatic diazine core on the optoelectronic properties of the designed D-p-diazine-p-D framework molecules, we have calculated the topological HOMA (harmonic oscillator model for aromaticity) aromaticity indices [60] of the diazine ring. The HOMA indices are computed using Multiwfn program [44] and the obtained results are presented in Table 3. The obtained HOMA values indicate that among the modeled molecules, 1c, 2b, 2c and 3b are characterized by the more aromatic diazine core which is consistent with their lowest values of the singlet-triplet ΔEST energy gap. We conclude that the TADF properties of the studied molecules are as much affected by the aromaticity of the diazine core whether it is pyridazine, pyrimidine or pyrazine as by the linking positions of the donor groups.
<Insert Table 3 here>
Natural Population Analysis (NPA)
As part of the natural bonding orbital (NBO) method [61], the net natural charges of the two nitrogen atoms of the diazine core calculated at B3LYP/6-31G(d,p) and HF/6-31G(d,p) [62] levels of theory are given in Table 3. The negative charge values obtained clearly validate the electron-withdrawing capacity of the diazine ring. It is clear that the nitrogen charge is affected both by the linking positions of the donor groups and by the diazine kind, it varies in the following order : pyridazine < pyrazine < pyrimidine. The Table 3 shows that compounds 1c, 2b, 2c and 3b are characterized by the more electron-deficient diazine core which is consistent with their lowest values of the singlet-triplet energy gap.
NTO analysis in S1 and T1 excited states
The natural transition orbitals (NTOs) describing the S1 and T1 excited states of the designed molecules are examined in order to better characterize the nature of the electronic transitions and to learn more about the charge transfer properties within these molecules. Figure 2 displays the hole-particle pairs of NTOs for the first singlet and triplet excited states. The highest occupied NTO (HONTO) is represented by the hole, whereas the lowest unoccupied NTO (LUNTO) is represented by the particle. The weight (ν) of the corresponding NTOs pair is also given. According to Figure 2, the NTOs for the lowest singlet excitation of the all investigated molecules showed good agreement with their HOMO and LUMO distributions (Figure S1). For the S1 and T1 excited states in the eleven emitters, the hole is mainly localized on the dimethyl-acridine donor moiety, whereas the particle is distributed on the diphenyl-diazine fragment. There is practically no overlap between holes and particles, which denotes that the S1 and T1 states have charge transfer (CT) type excitations. On the other hand, according to the obtained results, by adjusting the linking position of the two donor moieties for the same diazine kind (1a,1b, 2a and 3a for pyridazine, 1c, 2b, 2c and 3b for pyrimidine, 1d, 2d and 3c for pyrazine) lead to the decreasing of S1 and T1 energy levels (see Table 2). Furthermore, the hole and particle overlap integral (Sh/e) [41], the hole-electron centroid distance (D) [42] of the singlet and triplet states are calculated for a quantitative comparison and listed in Table 4 along with the (Dr) [43] index values. The integral of overlap of hole-electron distribution is a measure of spatial separation of hole and electron and it indicates how much the hole and the particle overlap. A smaller value of Sh/e index denotes a greater hole-electron separation, which increases the possibility of the rISC process [63]. The distance D between centroid of hole and electron is a measure of CT length. In general, longer centroid distance of the orbitals indicates more CT transition component. We have also calculated the charge transfer index (Dr) [43] to assess the degree of charge transfer. The degree of charge transfer is greater the higher the value of Dr. From the obtained results (Table 4), the Sh/e values are closer to zero implying that the CT character is dominant in the studied molecules. Additionally, all of the molecules have large Dr indices (larger than the threshold value 2.0 Å) [43] leading to strong CT character type of the electronic excitation. It is clear from Table 4 that in general the decrease in ΔEST values (from 135 to 4.2 eV) is consistent with the decrease in Sh/e values (from 0.114 in S1 and 0.326 in T1 for 1a to 0.027 in S1 and 0.029 in T1 for 3c) which increases the possibility of the triplet's up-conversion in these compounds.
<Insert Figure 2 here>
<Insert Table 4 here>