Theoretical Approach for the Luminescent Properties of Ir(III) Complexes to Produce Red–Green–Blue LEC Devices

With an appropriate mixture of cyclometalating and ancillary ligands, based on simple structures (commercial or easily synthesized), it has been possible to design a family of eight new Ir(III) complexes (1A, 1B, 2B, 2C, 3B, 3C, 3D and 3E) useful as luminescent materials in LEC devices. These complexes involved the use of phenylpyridines or fluorophenylpyridines as cyclometalating ligands and bipyridine or phenanthroline-type structures as ancillary ligands. The emitting properties have been evaluated from a theoretical approach through Density Functional Theory and Time-Dependent Density Functional Theory calculations, determining geometric parameters, frontier orbital energies, absorption and emission energies, injection and transport parameters of holes and electrons, and parameters associated with the radiative and non-radiative decays. With these complexes it was possible to obtain a wide range of emission colours, from deep red to blue (701–440 nm). Considering all the calculated parameters between all the complexes, it was identified that 1B was the best red, 2B was the best green, and 3D was the best blue emitter. Thus, with the mixture of these complexes, a dual host–guest system with 3D-1B and an RGB (red–green–blue) system with 3D-2B-1B are proposed, to produce white LECs.


Introduction
The concept of Solid-State Lighting (SSL) promotes energy savings and greenhouse gas reduction compared to conventional lighting (incandescent bulbs or halogen lamps) [1,2]. The most common SSL devices are LED (Light Emitting Diode) and OLED (Organic Light Emitting Diode) [3,4]. LED technology is based on high purity inorganic semiconductors, for example: AlGaAs, InGaN, GaN, and ZnSe, among others. These materials provide highly efficient and convenient point sources of light of different colors, depending on the semiconductor used [3]. Instead, OLEDs are processed in a multilayer system, using neutral organic or organometallic luminescent compounds, sandwiched between two electrodes [5][6][7]. To improve the manufacturing, costs, and performance of these SSLs, LECs (Light Emitting Electrochemical Cells) emerge as a promising alternative [8].
All SSL systems work through electron-hole recombination processes, where electrons are injected from the cathode into the luminescent material and holes are injected from the anode into the luminescent material. The recombination zone must take place in the middle of the luminescent material, producing an exciton, which if radiatively deactivated could produce light [9][10][11][12]. LEC devices contain a thin film of an ionic luminescent material between a cathode and an anode. This material is usually an ionic transition metal complex (iTMC) [13,14] and their films can be easily processed from an organic solution, using the spin-coating technique, which provides thinner and homogeneous layers [15]. Due to the ionic nature of the iTMC (unlike neutral compounds used in OLEDs), multilayer arrangements are not required in LEC because the electron-hole recombination is determined by the ionic mobility of the active molecules; therefore, LECs are a low cost alternative compare to OLEDs [16].
The cyclometalated Ir(III) complexes ([Ir(CˆN) 2 (NˆN)] + , where CˆN is a cyclometalating ligand and NˆN is an ancillary ligand) show the best performance as luminescent material in LECs [17][18][19]. The Ir atom is characterized by a high ligand-field splitting energy, thus in their complexes the metal centered (MC) excited states are less thermally accessible, avoiding the non-radiative decays [20,21]. Besides, the spin-orbit coupling (SOC) is increased in Ir-iTMC compared to complexes with metals of the 1st or 2nd transition rows, thus, a very effective singlet-to-triplet intersystem crossing occurs in Ir-iTMC, promoting radiative decay with a high quantum yield [22].
Another outstanding characteristic is that the Ir(III) complexes show emissions in a wide range of the visible spectrum, which can be modulated by the incorporation of an electron donor and/or acceptor substituents in both CˆN and NˆN ligands [22][23][24]. This can be understood since the electron density distributions of the frontier molecular orbitals (HOMO: Highest Occupied Molecular Orbital, and LUMO: Lowest Unoccupied Molecular Orbital) are similar in most of the Ir-iTMCs [23,25]. The HOMO has contributions from Ir dπ orbitals and π orbitals of the CˆN ligand, and the LUMO is composed mainly of π* orbitals centered in the NˆN ligand [22,23,26]. Consequently, the emitting triplet state (T 1 ) commonly has a 3 MLCT/ 3 LLCT (MLCT: metal-to-ligand charge transfer, LLCT: ligandto-ligand charge transfer, from CˆN to NˆN) mixed character [22,23,25,26]. The modulation of the emission energy involves, for example, the use of electron-withdrawing substituents in the CˆN ligand, decreasing the electron density on the metal, leading to the stabilization of the HOMO level. Also, electron donor substituents can be incorporated into the NˆN ligand, mainly destabilizing the LUMO. The HOMO-LUMO gap is increased with both strategies, obtaining high emission energies. Conversely, with the HOMO destabilization and LUMO stabilization, emissions at lower energies can be obtained [20,25,27].
The design of efficient Ir(III) complexes for LECs involve, in many cases, complicated synthetic procedures of the ligands, obtained at very low yields, affecting the final cost of the devices [19,28,29]. In this sense, the challenge to produce new Ir-iTMCs should involve the use of ancillary and cyclometalating ligands with favorable synthetic routes, providing specific emission colors and, of course, stable compounds.
Considering the wide range of emitting colors obtained from the Ir-iTMCs, the contributions to white light can be feasibly increased, exploring strategic blends of CˆN and NˆN ligands to produce new stable RGB Ir(III) complexes. According to the antecedents commented on, in this work are evaluated by theoretical calculations, the photophysical properties of an extended family of Ir(III) complexes (see Figure 1) with a mix of CˆN and NˆN ligands, with different electron-donating and/or -withdrawing natures and different structures to enhance the LEC performance. Therefore, the design of three series was carried out with the purpose of tuning the emitting color of the complexes, from red to blue, mainly determined by the cyclometalating ligand. In this sense, series 1 was expected to show a red emission; series 2, a green emission; and finally, series 3, a blue emission. However, some discrepancies were observed, namely, the 3E complex exhibits strongly green-shifted emissions, which will be discussed in more depth in this study. The calculations were carried out using the Density Functional Theory (DFT) level, providing a deep structural characterization of the ground and excited states, as well as a study of the absorption and emission properties, and the associated photophysical efficiency parameters, to identify the complexes with the best emitting properties. Finally, the best combinations of blue, green, and red emitters are proposed with the aim to contribute to the strategic design of white LECs.

Computational Details
All calculations were performed using the DFT approach with the B3LYP functional [41,42], which proved to be the most suitable to reproduce the absorption and emission properties of the [Ir(ppy) 2 (bpy)] + complex, which was used as the reference complex for the DFT benchmark (see Table S1 of Supporting Information (SI)). The LANL2DZ [43] basis set and quasi-relativistic pseudopotential were adopted for the Iridium atom, while the 6-31G(d,p) [44,45] basis set was used for the other atoms. The optimized structures for all complexes correspond to the energy minima, according to the vibrational frequencies (real values). The Time-Dependent DFT (TD-DFT) methodology was used to calculate the first 50 and 6 triplet excited states, respectively. The triplet (T 1 ) excited states were optimized using TD-DFT gradients (TD-DFT optimization). Emission energies were estimated as the vertical energy difference between the total energies of the relaxed triplet state and ground state at the optimized triplet geometry. The continuous solvent effect was incorporated by the IEF-PCM model (polarized continuum model), [46] using dichloromethane as the solvent. The optimized structures of 3 MC states were also determined; for this purpose unrestricted triplet optimization calculations were performed on the bases of the geometry of the distorted excited triplet state, where the six coordinate bond lengths around the metal atom (Ir-C CˆN , Ir-N CˆN and Ir-N NˆN ) are all gradually elongated (up to about 0.8 Å) and spin density distribution calculations confirm that the metal-centered excited states are obtained, as described in the literature [47,48]. The ionization potentials (IP), the electron affinities (EA), and the hole/electron reorganization energies were obtained by the differences between the total energies of the molecular system in its fundamental state and with one more electron or one less electron [49][50][51]. All the calculations were performed in the Gaussian16 program package [52] and the wavefunction analyses were performed in the Multiwfn 3.4 code [53].

Molecular Geometries in the Ground States and Lowest Excited Triplet States
The different Ir(III) complexes studied are shown in Figure 1, with the numbering of some key atoms. According to this nomenclature, the selected geometric parameters of the ground states (S 0 ) and the excited triplet states (T 1 , T 2 or T 3 ) are summarized in Table 1. All the data of the determined geometric parameters are shown in Table S2 of the Supplementary Materials. The calculated data are in agreement with the bond lengths and angles reported in literature for similar cyclometalated Ir(III) complexes [15,19,54,55]. In the ground state it is possible observe that the all complexes have a distorted octahedral geometry around the Ir metallic center, and exhibit similar geometric parameters between them, namely, for the Ir-C and Ir-N bonds of the CˆN ligand, an average length of 2.02 and 2.08 Å is found, respectively; while the Ir-N NˆN bond lengths are in the range of 2.19 and 2.32 Å. This is significantly longer than those for CˆN, which is attributed to the strong trans-effect of the C donor in the CˆN ligands, in agreement with the reported literature for the Ir(III) complexes with similar characteristics [56][57][58][59]. Furthermore, the bond angles involving the metallic center are also very similar between all complexes; the C 1 -Ir-N 4 angle is~96 • and the C 1 -Ir-N 3 angle is around 172-179 • for all the complexes. Finally, the C 1 -Ir-C 2 angle is range between 82-89 • .
Comparing the geometric parameters of the triplet excited states with respect to the parameters of the S 0 structures, marginal variations are found in the distances between the metal center and the ligands; some significant variations were found for Ir-N NˆN (Ir-N 2 and Ir-N 3 ) in almost all complexes, where a shortening is identified in the excited states.

Frontier Molecular Orbitals Analysis
To understand the absorption and emission properties of the studied complexes, it is necessary to analyze their electronic structures in the ground state, namely: the electron distributions of HOMO and LUMO, the energies of the frontier orbitals, and HOMO-LUMO gaps (∆HL). The contributions of molecular fragments to each molecular orbital (including energies) are summarized in Tables S3 and S4    The HOMO orbital is distributed between the d-orbital of the metal center (33-41%) and the phenyl ring orbitals of the CˆN ligand (57-62%), and the LUMO orbital being exclusively on the NˆN ligand (~95%) in all complexes. This electron density distribution agrees with the literature data for similar Ir(III) complexes [15,25,55].
Series 1 displays small values of ∆HL; 1A = 3.00 eV and 1B = 3.39 eV. With the increase of an aromatic ring fused in the framework of the A ancillary ligand, it is incrementing the π-accepting character of this ligand. The increment of conjugation in the ancillary ligands causes an enhancement of the electron withdrawing effect, therefore, the LUMO orbitals are stabilized, reducing the ∆HL of 1A compared to 1B [60,61]. Thus, since 1A shows the lowest ∆HL of the complexes studied, the absorption and emission energy of this complex is expected to be the most red-shifted. Series 2 has intermediate ∆HL values; 2B = 3.65 eV and 2C = 3.61 eV. This variation, with respect to series 1, is mainly attributed to the energetic stabilization that the HOMO orbital undergoes due to the electron-withdrawing character of the fluorine atoms present in the CˆN ligand [62][63][64][65]; therefore, series 2 is expected to present emissions in the green region. This trend has been observed in similar Ir(III) complexes, where the presence of the electron-withdrawing substituents present in the aromatic rings of the CˆN ligands promote greater stabilization of the HOMO orbital, demonstrated, for example, by cyclic voltammetry studies, where the oxidation peak associated with the CˆN ligand and the d orbitals of the metal is shifted to higher potentials, with respect to the oxidation of the analogous complexes with the CˆN ligands without the electron withdrawing substituents [66].
Finally, the largest ∆HL values are for the series 3; 3B = 4.02 eV and 3C = 3D = 4.01 eV, except for 3E = 3.51 eV. The introduction of the N onto the aromatic ring of the CˆN ligand, plus the presence of the fluorine atoms, has a significant effect on the HOMO energy stabilization, [23] increased in the HOMO-LUMO energy gap compared to the complexes of the series 1 and 2, which probably will result in blue-shifted emissions. In the case of the 3E complex, its ∆HL is considerably reduced due to the stabilization of the LUMO by the presence of the E ligand, with high electron delocalization that increases their electron acceptor character, as has been described in analogous Ir(III) complexes with a 2,2 -biquinoline as NˆN ancillary ligand [38,67].
In summary, we observe that the presence of the fluorine atoms in the CˆN ligand (in series 2), as well as the additional incorporation of the nitrogen atom in the CˆN ligand skeleton (series 3), has a significant influence in the blue-shift of the lowest absorption bands compared to series 1, which allows the prediction that such complexes present the expected emission at a higher energy. On the other hand, the remarkable participation of MLCT in the region of interest could be indicative of an efficient intersystem crossing throughout the three series, which anticipates a good emission performance.

Emission Properties
The luminescent properties of the complexes in this study were determined based on the minimum energy structures of the lowest triplet excited states, obtained using the TD-DFT formalism. The main results are listed in Table 3.  In series 2, it is observed that the fluorine atoms in the cyclometalating ligand cause a displaced emission at high energy in complexes 2B and 2C (~2.26 eV), compared to series 1, as we expected, since this strategy has been widely used to favor blue-shifted emissions [64,65]. For both 2B and 2C the emitter state found corresponds to T 2 . The emission of 2B originates from the LUMO→HOMO-1 transition with a 3 MLCT and 3 ILCT character, while, for 2C, it is predominantly contributed to by the LUMO+1(LUMO+2)→HOMO transitions and its emissive state shows a mixed 3 MLCT/ 3 LLCT/ 3 LC character.
In series 3, the addition of fluorine and nitrogen atoms to the CˆN ligand blue-shifts the emission further with respect to series 2, i.e., the emission energies obtained were 2.76, 2.73 and 2.82 eV, for 3B, 3C and 3D, respectively. The exception is 3E, which displays a significantly lower emission energy compared to the rest of the complexes, shifting its emission to the yellow region of the visible spectrum; this behavior is attributed to the increased conjugation of the E ancillary ligand, which stabilizes the energy of LUMO, as commented previously. The emitting state for 3B, 3D and 3E corresponds to T 2 , while for 3C, it is T 3 . Two trends are observed in the nature of the emitting excited state: the complexes 3B and 3D can be described with 3 MLCT/ 3 LLCT mixed characters and for 3C and 3E, a mixed 3 MLCT/ 3 LLCT/ 3 LC character. In all these complexes, the deactivation pathway originates from LUMO towards HOMO (HOMO-1, HOMO-2).
Regarding the emission energies found, the designed complexes showed a significant emission color tuning from deep red to blue, i.e., the complexes 1A and 1B could be exhibiting red emissions; for 2B, 2C and 3E the emissions ranged at green, and finally, emissions at blue should be displayed for 3B, 3C and 3D.

Phosphorescence Quantum Efficiency
From the photoluminescence quantum yield (Φ) we can quantify the efficiency of the emission process. The Φ is determined by the competition between the radiative rate constant (k r ) and the non-radiative rate constant (k nr ), according to the following relation: [68,69] Φ = k r /(k r + k nr ) The k r is usually expressed as: [63] k r = γ( Ψ S1 |H S0 |Ψ Tn 2 µ S1 2 )/(∆E(S 1 − T n ) 2 ); with γ = (16π 3 10 6 n 3 E emi 3 )/(3hε o ) (2) where Ψ S1 |H S0 |Ψ Tn is the spin-orbit coupling (SOC) matrix element between the S 1 state and the emitting state (T 1 , T 2 or T 3 ), µ S1 is the transition electronic dipole moment in S 0 →S 1 transition, ∆E(S 1 − T n ) is the energy gap between the S 1 state and the emitting triplet state, and n, E emi , h and ε 0 are the refractive index of the medium, emission energy, Planck's constant, and the permittivity in vacuum, respectively. This expression is applicable to coordination compounds with a heavy metal center since the radiative rate is directly proportional to the SOC matrix element related to the emitting triplet and singlet states, and inversely proportional to the degree of mix between them (∆E(S 1 − T n )). In these type of complexes, large intersystem crossing rates (ISC) predominate, so the fluorescence rate can be considered null, as well as some non-radiative pathways (inverse ISC, internal conversion, conversion from MLCT to MC) [70]. The SOC effects can be elucidated from the metal contribution in the emitting state (% 3 MLCT), which were calculated from the sum of the 3 MLCT contributions from each monoexcitation (according to their respective orbitals involved), considering the corresponding configuration coefficient, as has been reported in the literature [71,72]. The results are summarized in Table 4. Table 4. Metal-ligand charge transfer character ( 3 MLCT, %), transition electric dipole moment (µ S1 , D) and energy gaps between the S 1 and T n states (∆E(S 1 − T n ), eV) of studied complexes. According to Equation (2), higher values of % 3 MLCT and µ S1 increase the k r value and, conversely, higher values of ∆E(S 1 − T n ) decrease the k r . We observe that, in all complexes, the % 3 MLCT is high enough, therefore, k r is favored in all cases. However, the higher values are obtained for the complexes of series 1 and series 3 (except for 3E), which suggests that the % 3 MLCT is affected by the different CˆN ligands used. Interestingly, these complexes showed structural changes in the excited state (T 1 , T 2 or T 3 ) that would favor the metal-ligand interaction, namely, a shortening of the length of Ir-N NˆN ligand bonds (Ir-N 2 and Ir-N 3 ).
For µ S1 , the variation of the CˆN and NˆN ligands do not show a clear trend. The largest values are shown by 1B, 2B, 2C, and 3C (between 1.26 to 1.00 D), so higher k r values are expected in these complexes.
On the other hand, a minimal difference between the S 1 and the emitting state, ∆E(S 1 -T n ), is favorable for enhancing the intersystem crossing (ISC) rate, implying an increase in k r . The results show significantly lower values for complexes 3D, 3B, 1A, and 1B, and precisely, in these complexes, the highest values of % 3 MLCT (~27% to 33%) are observed, as expected since a high SOC (% 3 MLTC) promotes an effective ISC.
In summary, considering all the parameters that determine the k r value, it is observed that the 1A, 1B, 2B, 2C, 3B, and 3D complexes would present characteristics that favor the radiative deactivation processes.
Equation (1) shows that another determining factor in the efficiency of the emission process corresponds to the k nr . The population of metal-centered ( 3 MC, d-d*) triplet excited states is one of the most important deactivation pathways of the phosphorescence, namely, the 3 MLCT (π-π*) excited state can be rapidly converted to a short-lived 3 MC (d-d*) state, from which no emission occurs; to avoid this situation, a splitting in energy is required between the 3 MC and 3 MLCT states [25,73,74].
The metal-ligand bond lengths and bond angles in 3 MC states are listed in Table S5 (see Supplementary Materials), along with the spin density distributions determined for the 3 MC state (see Figure S7) where the metal-centered character is observed. The results show that, in almost all complexes, the main structural change in the 3 MC states correspond to the distance between the metal and the N atoms located in the axial position (Ir-N 1 /Ir-N 4 ), which increase up to 0.6 Å, with respect to ground state. However, complexes 1A and 3C showed a relatively small structural distortion between the 3 MC and S 0 states, which would explain the lower energies obtained for these 3 MC states.
In general, it is observed that, in almost all complexes, the relative energy position of the 3 MC state lies above the emitting triplet state 3 MLCT, except for complexes 3B and 3C, as shown in Figure 4. This is explained due to the electron withdrawing character of the cyclometalating ligand (3), therefore, it strongly influences the increase in the energy of the 3 MLCT excited state, which, in some cases, causes it to exceed the energy of the 3 MC state. High energies of the 3 MLCT excited states have been identified for analogous Ir(III) complexes with this cyclometalating ligand [75][76][77][78]. This effect is almost experienced by 3D ( 3 MLCT energy close to 3 MC energy) and, in the case of 3E, the effect is not observed since the E ligand (electron withdrawing by resonance) compensates for the effect of the cyclometalating ligand. With regard to the behavior of the 3B and 3D complexes, a similar tendency should be expected, since the aliphatic bonds and substituents in the 4 and 4 positions should not influence the electronic properties. However, to elucidate why, in one case, the 3 MC (3B) is stabilized and in the other, the 3 MLCT (3D) is destabilized, we calculated the root-meansquare deviation (RMSD), determining the difference between two sets of the coordinate values, in this case between the 3 MLCT and 3 MC excited states (see Supplementary Materials Figure S8). For the 3B complex, a higher value (0.668 Å) was found compared to 3D (0.494 Å), which would explain the higher energy gap in 3B, probably due to the high mobility of the substituents in the 4,4 positions, yielding a 3 MC state that is more stabilized. Whereas in 3D, the lower value of RMSD indicates less distortion (in fact, both states have very similar energies), where the 3 MLCT state is most stable.
The computed adiabatic energy differences between the 3 1A, 1B, 2B, 2C, 3B, 3C, 3D, and 3E, respectively. High ∆E( 3 MC-3 MLCT) values of iTMCs, in the range of 0.26-0.60 eV, have shown high electroluminescent performance in LEC devices [79]. Consequently, the complexes 1B, 2B, and 3E would lead to a lower k nr and would present a high performance in their application in LEC. While the probability of populating the 3 MC states increases in complexes 1A, 2C, and 3D, which would result in a more favored k nr , as the emission of 3B and 3C comes from very high 3 MLCT states (blue-shifted emission), [21] the conversion from the emitting state to the 3 MC state is easily completed, thereby increasing the probability of non-radiative decay to the ground state and consequently, a significantly larger k nr is obtained. Note that the 3D blue emitting state also displays a high 3 MLCT state.

Charge Transfer Properties: IP and EA
The good performance of LEC devices is determined by the appropriate injection and transfer of holes and electrons, which can be evaluated by the ionization potential (IP), electron affinity (EA), and reorganization energy (λ). The IP and EA are obtained as: [49][50][51] where E N , E N+1 , E N-1 are the total energies of the molecular system in its ground state, with one electron more and one electron less, respectively. Small IP values indicate easy hole injection from the anode to the HOMO of the iTMC, and large EA values can be related to a favored electron injection from the cathode to the LUMO of the iTMC [25,80,81]. As shown in Table 5, the IP values gradually increase in the following order: Series 1 < Series 2 < Series 3, which is consistent with the HOMO energy levels (see Section 3.2). Therefore, the CˆN ligand largely determines the HOMO level, establishing that ligand 1 causes a lower hole injection energy barrier compared to ligands 2 and 3. In this sense, series 1 could display the best hole injection performance. Related to the EA, the highest values are found in 1A and 3E, showing agreement with the lowest energies of the LUMO levels, according to the electron acceptor character of the NˆN ligand. Consequently, these three complexes will have an enhanced electron injection ability compared to other complexes.
On the other hand, the reorganization energy can be used to estimate the charge transport rate and balance between holes (λ h ) and electrons (λ e ) according to the following expression: [23,82,83] λ h = IP − HEP (5) HEP (EEP) is the hole (electron) extraction potential and is determined as the vertical energy difference between the ground state and the relaxed state with one electron less (more), using the geometry with one electron less (more), respectively. Generally, a low reorganization energy (λ h , λ e ) is necessary for an efficient charge transport process and in this respect, the full series 3 appears to be more efficient in the hole transport (λ h = 0.12 to 0.13 eV) while complexes 1A, 2C, 3C, and 3E showed greater efficiency in electron transport. However, in all complexes, the hole transport performance is favored over the electron transport ability, due to the higher values obtained of λ e (0.32 to 0.46 eV), with respect to λ h (0.12 to 0.19 eV).
In addition, we determined the difference between λ h and λ e , and it was observed that the complexes 1A, 2C, 3C, and 3E showed less discrepancy (∆λ < 0.21), indicating that the balance of electron and hole transfer could be easily achieved in the emitting layer of the LEC devices.

RGB Systems to Produce White LECs
According to the emission energies determined for each complex studied, in Table 6, in the first line, the complexes were organized from the reddest emitter to the bluest emitter, from left to right, respectively. It can be observed that this order differs with respect to the original tendency designed, since it was expected that the type of cyclometalating ligand chosen was the main factor determining the emission energy. In this sense, the 3E complex is out of the original tendency since their ancillary ligands promote a strong stabilization of LUMO, which shifted their emissions at lower energies. Then, in the following lines of the table, an arbitrary qualification has been used for each studied complex with respect to the values of k r , k nr , IP, EA, and ∆λ. This qualification arises from the parameters determined in Sections 3.5 and 3.6, based on the minimum and maximum values determined, and also considering key values from the literature. Then, the +++ symbol is assigned to a very favourable value, ++ to a favourable value, and + is slightly favourable. In the case of k nr , a 0 value is added to qualify the parameter with an unfavourable value. The extension of the criteria assigned can be observed in the Supporting Information. By the analysis of k r and k nr kinetic parameters, that describe the efficiency of the intrinsic radiative process of the complex, it is observed that, for red emitter complexes, the best combination of these parameters is obtained in complex 1B, therefore, this is the best candidate to be a good red emitter in an RGB system. Then, if IP, EA, and ∆λ parameters are considered, 1A could be an appropriate red emitter. Next, in terms of the evaluation for the green emitter complexes, the best performance in the kinetic parameters of the phosphorescent processes is found for 2B, therefore being the best candidate as a green emitter. Then, if the charge transport and injection parameters are considered, plus its balance, 3E and finally 2C could be proposed.
Finally, for the emitter complexes in the blue range, only 3D could act as an appropriate emitter for an RGB system, since the other blue complexes have a very unfavourable k nr . This behaviour in 3C and 3B is due to the higher energy of the 3 MLCT states with respect to 3 MC, promoting this 3 MC state as the phosphorescence deactivator, which is known to prefer non-radiative paths.
Consequently, according to this analysis, a proposal for dual host-guest system with blue and red emitters to produce white LECs could be obtained by mixing 3D-1B, as well as 3D-1A.
In the case of an RGB system, the proposal would be to assemble the 3D-2B-1B complexes.

Conclusions
A detailed investigation is reported on the geometrical and electronic structures, emission properties, charge injection/transport abilities, and phosphorescence efficiency of eight new Ir(III) complexes (classified into 3 series) using DFT and TD-DFT methods.
The designed complexes showed a significant emission color tuning from deep red to blue, with emissions ranging from 440 to 701 nm. The results showed that the Ir-NˆN ligand bond lengths are shortened for complexes of series 1 and 3, which has a direct impact on the metal-ligand interaction, leading to more involvement of the metal in the triplet excited state (% 3 MLCT). The analysis of quantum efficiency showed that the 1A, 1B, 2B, 2C, 3B, and 3D complexes would present characteristics that favor the radiative deactivation processes, while 1B, 2B, and 3E complexes would lead to a lower k nr and would present a high performance in their application in LEC. In relation to transport and injection parameters, the complexes with the best balance between hole and electron injection are 1A, 2C, 3C, and 3E (∆λ < 0.21), however, the full Series 3 stands out for being more efficient in the hole transport, and 1A, 2C, 3C, and 3E, for presenting greater efficiency in the transport of electrons. Finally, we have proposed a host-guest dual system using a mix of 3D and 1B, and also an RGB system based on 3D-2B-1B, to produce white LECs.
With this family of complexes, corresponding to newly designed structures based on a mix of CˆN and NˆN ligands, commercially available and/or easy to synthesize, an important contribution of new Ir-iTMCs is provided from a theoretical approach, which is totally experimentally viable.
Supplementary Materials: The following supporting information can be downloaded at: https:// www.mdpi.com/article/10.3390/molecules27092623/s1. Section S1. Benchmark study with other DFT functionals; Table S1. Absorption and Emission wavelengths from ground and triplet excited states, respectively, for [Ir(ppy) 2 (bpy)] + (in dichloromethane); Section S2. Geometry parameters of the ground states and the triplet states; Table S2. Optimized geometric parameters of all complexes under study in the S 0 and T 1 states, determined at the B3LYP/6-31G(d)-LANL2DZ level of theory; Section S3. Energies of the molecular orbitals from HOMO-2 to LUMO+2, in the ground state; Table S3. Energy of frontier molecular orbitals (eV) and ∆HL of the ground state (S 0 ); Section S4. Contribution to molecular orbitals from HOMO-2 to LUMO+2, in the ground state; Table S4. Contribution to molecular orbitals (%) of all complexes calculated from HOMO-2 to LUMO+2 in the ground state (S 0 ); Section S5. Surface orbitals from HOMO and LUMO in the ground state; Figure S1. Surface of frontier molecular orbitals HOMO and LUMO (S 0 ) of series 1; Figure S2. Surface of frontier molecular orbitals HOMO and LUMO (S 0 ) of series 2; Figure S3. Surface of frontier molecular orbitals HOMO and LUMO (S 0 ) of series 3; Section S6. Molecular orbitals of the triplet excited state; Figure S4. Frontier molecular orbitals involved in the radiative deactivation of the lowest-lying triplet excited state of series 1; Figure S5. Frontier molecular orbitals involved in the radiative deactivation of the lowest-lying triplet excited state of series 2; Figure S6. Frontier molecular orbitals involved in the radiative deactivation of the lowest-lying triplet excited state of series 3; S7. Geometry parameters of 3 MC state; Table S5. Selected geometric parameters of all complexes under study in the 3 MC state, determined at the B3LYP/6-31G(d)-LANL2DZ level of theory; Section S8. Spin density of 3 MC state; Figure S7. Spin density distribution in the optimized 3 MC state for representative complexes; Section S9. Root-mean-square deviation (RMSD) between the 3 MLCT and 3 MC states; Figure S8. Superimposed structures of the 3 MLCT states (blue) and 3 MC states (orange) for 3B and 3D, and their RMSD values; Section S10. Qualification scales of kr, knr, IP, EA and ∆λ to choose RGB systems; Table S6. Photophysics and charge transport parameters to determine the best RGB systems; Table S7. Building scales to quantify the photophysics and charge transport parameters to determine the best RGB systems.