Structural, Electronic, and Charge Transport Properties of New Materials based on 2-(5-Mercapto-1,3,4-Oxadiazol-2-yl) Phenol for Organic Solar Cells and Light Emitting Diodes by DFT and TD-DFT

)is work reports on the density functional theory (DFT) and its time-dependent extension (TD-DFT) study of the structural, electronic, and charge transport properties of 2-(5-mercapto-1,3,4-oxadiazol-2-yl) phenol (MODP) and some of its transitionM complexes (M� Fe, Co, Cu, Ni, Zn, Pd, Pt). Reorganization energy, integral charge transfer, mobility, open circuit voltage, and electronic properties of these compounds have been calculated by employing the global hybrid functional PBE0 in conjunction with the Karlsruhe basis set def2-TZVP. Results show that MODP and its transition metal complexes are good electron donors for organic solar cells (OSC) owing to their relatively higher HOMO and LUMO energies compared to the prototypical (6, 6)-phenylC61-butyric acid methyl ester (PCBM). Energy gaps ranging between 2.502 and 4.455 eV, energy driving forces (∆EL-L) ranging between 2.08 and 2.44 eV, and large open circuit voltages (VOC) ranging from 1.12 to 2.05 eV were obtained. )e results also revealed that MODP and its Pd(II) and Pt(II) complexes could serve as ambipolar charge transport materials owing to their very small reorganization energies, integral charge transfers, high rate charge transfers, and mobilities. All studied molecules showed OSC donor and hole/electron transport characteristics required by organic light-emitting diodes (OLEDs). Based on these results, new ways for designing charge transport materials for OLEDs as well as donor materials in OSCs are proposed.


Introduction
Organic optoelectronic devices like organic light-emitting diodes (OLEDs) and organic solar cells (OSCs) have attracted tremendous research attention in recent years [1][2][3][4]. OLED devices find applications in display and lighting [5,6], while OSC devices are most promising for electrical energy applications [7]. Organic electronics based on small π-conjugated molecules and polymers are ideal for optoelectronic applications due to their absorption, light emission, and charge transport properties [4][5][6][7][8]. Although enormous advances have been made in the search and development of new materials, the low efficiency and thermal stability of small molecule-based organic electronics materials still remain a problem for their application and has limited their marketing. It is therefore necessary to develop new multifunctional organic materials with high-efficiency for OLEDs and OSCs, which are more stable and capable of functioning as better light emitters, absorbers, and charge transporters. Such materials serve as efficient light emitters in OLEDs or donor materials capable of converting abundant sunlight energy directly into electricity for OSCs and as charge transport materials, or both. OLEDs and OSCs are devices whose components comprise an emitting layer (EL) made out of a film of organic compound, a hole transport layer (HTL), and an electron transport layer (ETL), all located between two electrodes [2,9]. In OLED devices, N,Ndiphenyl-N,N-bis (3-methylphenyl) (1,1-biphenyl)-4,4-diamine (TPD) is often used as a good material for hole transport because of its efficient hole injection and high hole mobility properties [10], whereas tris (8-hydroxyquinolinato) aluminum (Alq3) is commonly used as a material for electron transport thanks to its Highest Occupied Molecular Orbital (HOMO) and Lowest Unoccupied Molecular Orbital (LUMO) energy levels, as well as its effective hole blocking ability [11].
Recently, several studies on compounds with a D-π-A framework have been carried out [1,8,12]. e presence of an extended delocalization of π-electrons between a donor (D) and an acceptor (A) is a key element in the design of OLEDs and OSCs [1,13]. is has stimulated great interest in the synthesis of π-conjugated systems comprising fullerene derivatives [14]. For instance, (6, 6)-phenyl-C 61 -butyric acid methyl ester (PCBM), a fullerene derivative, plays a crucial role in the synthesis of OSCs because of its stability, capacity to accept electrons, and its proper frontier molecular orbital (FMO) energy levels [8,14]. In OSC materials, the active layer generally consists of two components (i.e., a donor and an acceptor material), which are mounted in a two-layered structure or in the form of a mixture. It is worth mentioning that such an arrangement is observed in PCBM and its derivatives, which are mostly used as electron acceptor materials [8,14].
is arrangement is also seen in poly(3-hexylthiophene) (P3HT) and its derivatives, which are generally used as electron donors in OSCs [15,16]. For the design of OSCs and OLEDs materials, knowledge of the HOMO and LUMO energy levels of the compounds under investigation is of prime importance [1,17]. ese energy levels are also necessary to elucidate the propensity of charge transfer [17,18].
Within the framework of the design of OLEDs and OSCs, 1,3,4-oxadiazole and their derivatives serve as important materials thanks to their π-conjugated nature, high quantum efficiency, high photoluminescence, and their chemical and thermal stability [9,19,20]. ese compounds have proven to be good electron transporters (ET) and important hole blockers (HB) in organic optoelectronic devices [20,21]. Bhava and co-workers [22] synthesized and characterized Cu(II), Co(II), Ni(II), and Zn(II) complexes based on 2-(5-Mercapto-1,3,4-Oxadiazol-2-yl) Phenol (MODP) (see Figure 1) with the aim of searching for new antimicrobial agents based on the 1,3,4-oxadiazoles skeleton. In spite of the fact that these oxadiazole-based molecules and their metal complexes have been studied both experimentally and theoretically, very little attention has been focused on the structural, electronic, and charge transport properties of MODP and its transition metal complexes. It is clear from Figure 1 that MODP is a heterocyclic compound with promising push-pull potentials necessary for charge transport and photovoltaic properties of interest in OLEDs and OSCs. In addition to the metal ions that have been investigated experimentally in conjunction with MODP, Pd(II) and Pt(II) metal ions are also studied herein, since several works have shown that complexes bearing these ions are good materials for OLEDs [23,24].
In order to shed more light on the applicability of the studied molecules as OLED and OSC materials, this work sets out to study the geometric, electronic, and charge transport properties of MODP, as a bidentate ligand, along with its Cu(II), Ni(II), Co(II), Zn(II), Fe(II), Pd(II), and Pt(II) complexes. In order to achieve this goal, the HOMO and LUMO energies, (E HOMO and E LUMO ), the E HOMO − E HOMO gaps(Egap), the open circuit voltage (V oc ), the energy driving force (∆EL-L), the absorption spectra, reorganization energies (λ), integral charge transfer (V), and the λ max of electronic absorption have been computed and discussed. e density functional theory (DFT) method (and its time-dependent extension (TD-DFT)) has been used in this work since it is good compromise between computational time and accuracy and also takes into account the electron correlation [25,26].

Computational Details
All quantum chemical calculations herein were carried out using the ORCA 4.0.1 program [27]. In all calculations, the medium-sized numerical quadrature grid5 and the tight selfconsistent field (SCF) convergence criterion were utilized. e initial geometries of the molecules studied were prepared using the Avogadro visualization program [28]. Geometry optimization and frequency calculations were done with the aid of the low-cost composite electronic structure B97-3c method, which is very effective especially for geometry optimization of metal containing systems [29,30]. e B97-3c approach is used for molecular geometries and to exhibit excellent performance for non-covalent interaction energies of small and large molecules [29,31]. To confirm all optimized geometries as stable structures (i.e., minima on the potential surface energy), vibrational frequencies were computed at the same level of theory as that used for geometry optimization, and no imaginary frequencies were obtained.
Based on the B97-3c optimized geometries, single-point (SP) calculations were carried out using the hybrid-GGA functional PBE0 [30] along with the Karlsruhe basis set def2-TZVP [27,32]. e PBE0 functional [30] has been used in this study as it provides good energies for transition metals [33,34]. PBE0 has been combined to def2-TZVP because it is suitable for energy calculations [35]. e Grimme's atompairwise dispersion correction with Becke-Johnson damping dubbed (D3BJ) [36] was used in all calculations to account for long-range dispersion interactions. To speed up the calculations though with a minimal loss in accuracy, the Resolution-of-Identity (RI-J) approximation [37] was employed together with the Chain-of-Spheres (COSX) approximation [38] resulting in the RIJCOSX approximation. With the aid of the TD-DFT method [39], some electronic absorption properties (excited-state energies, maximum absorption wavelengths, λ max , and oscillation strengths (f)) of the investigated compounds were calculated based on the B97-3c optimized ground state geometries in the gaseous phase and in Dimethylsulfoxide (DMSO) at the CAM-B3LYP/def2-TZVP level of theory. e range-separated functional Coulomb-attenuating method applied to Becke's 3 Lee-Yang-Parameter (CAM-B3LYP) [39] was used because it is suitable for electronic spectra studies [39]. Implicit solvent effects on the absorption spectra of the compounds studied were included using the conductor-like polarizable model (CPCM) as implemented in the Orca 4.0.1 program [27,34].
Given that the performance of OSCs is influenced by the FMO energy level [40], two major FMO-based performance parameters were computed in this work. ese include the energy driving forces (∆E L − L ), which correspond to the energy difference between the LUMO energy levels of the donor (the oxadiazoles studied herein) and the acceptor (PCBM) and the open circuit voltage (V OC ). V OC is a significant parameter for studying organic solar cells [3] (see Figure 2). eoretically, V OC has been calculated as shown in the following equation [41,42]: where V OC is the open circuit voltage, E Donor HOMO is the HOMO energy of the donor, E Acceptor LUMO is the LUMO energy of the acceptor, and 0.3 is an empirical factor.
HOMO energy determines the ability of a molecular species to give out an electron, while the LUMO energy indicates its ability to accept electrons [43]. Studies have shown that the difference in energy between the HOMO of the donor and the LUMO of the acceptor significantly affects the V OC value [40,41].
∆E L − L is the determining factor of an exciton's dissociation efficiency at the donor/acceptor interface in OSCs materials as shown in Figure 2. ∆E L − L is calculated according to equation (2). It is worth noting that the value of ∆E L − L must be at least 0.3 eV; otherwise, dissociation will be unlikely [15].
V OC and ∆E L − L of the investigated molecules were obtained by carrying out Single Point calculations on the B97-3coptimized geometries at the RIJCOSX-PBE0-D3(BJ)/def2-TZVP level of theory, using equations (1) and (2). LUMOD and LUMOA are the LUMOs of the donor and acceptor, respectively. In a similar manner, HOMOD and HOMOA are the HOMOs of the donor and acceptor, and Egap D and Egap A are the energy gaps of the donor and acceptor, respectively. Marcus' theory, as summarized in equation (3), has been widely used to describe charge transfer and rate charge transfer speed [44,45].

Journal of Chemistry
where k B is the Boltzmann constant, T is the temperature, λ is the reorganization energy, h is the Planck constant, ΔG o is the standard Gibbs free energy for the change in the process, and V is the element of electron coupling matrices between two adjacent species. e reorganization energy (λ) and the transfer integral (V) are two key parameters that govern charge transport properties [46]. λ is the parameter, which determines the charge transfer properties of materials and is comprised of the intermolecular reorganization energy λ ext and the intramolecular reorganization energy λ int [47]. Attention is generally focused on λ int because environmental relaxations are ignored [48]. In general, the lower the λ values, the higher the charge transfer rate. In this paper, we focus on the internal reorganization energy exclusively because of the low dielectric constant of medium in OSCs materials. λ int is the sum of the electron (λ e ) and the hole (λ h ) reorganization energies, respectively, and evaluated according to equations (4) and (5) [46]: where E + 0 and E − 0 are, respectively, the cationic and anionic energies of a given molecule, which are computed using the optimized geometries of an electrically neutral structure; E + + and E − − are the energies of the cationic and anionic optimized geometries; and E 0 + , E 0 − and E 0 0 are energies of the neutral geometries computed using the cationic, anionic, and neutral geometries, respectively. In this study, the reorganization energies were calculated at the RIJCOSX-PBE0-D3(BJ)/def2-TZVP level of theory, based on the Single Point energies of the neutral and ionic structures of the studied molecules.
V represents the orbital coupling of neighboring molecules [24]. Dimers comprising monomers were constructed from the optimized geometries of the individual monomers using the B97-3c method. A 4.0Å distance between the monomer units of each dimer was used here, which corresponds to the minimal distance where the overlap between the molecular wavefunctions of the donor and acceptor systems could be neglected [49,50]. e simplest way to evaluate the charge transfer integral in a molecular dimer is by using the method of dividing the boundary orbital energy levels [50]. is method was used herein based on Koopmans' theorem, whereby the HOMO and HOMO +1 , LUMO and LUMO +1 energy levels were calculated. Using Marcus's theory, the charge transfer integrals of the electron (V e ) and the hole (V h ) were obtained according to equations (6) and (7) [49].
where V e and V h are the electron transfer and hole integrals, respectively, E H is the HOMO energy, E H−1 is the HOMO-1 energy, E L is the LUMO energy, and E L+1 is the LUMO+1 energy. In order to compute full charge transfer of the electron and hole, the PBE0-D3(BJ)/def2-TZVP method was used. e manufacture of high-performant OLED devices requires materials with good charge transport properties and high chemical stability, as well as a considerable shelf life [51]. Stability is an important factor that often limits the functionality of charge transport and luminescent materials [2]. From molecular orbitals perspectives, the prediction of the stability of the molecules studied was via chemical hardness (η), which was calculated according to equation (8) [51,52].
where µ is the chemical potential, N is the number of electrons. IP and EA represent the ionization potential and electronic affinity, respectively. e adiabatic IP of these molecules is defined as the energy difference between the cationic and neutral species, while the EA is the energy difference between the neutral and anionic species, respectively.
ese values were calculated according to equations (9) and (10): where E(M 0 ), E(M + ) and E(M − ) are the energies of the neutral, cationic, and anionic forms of the investigated molecules, respectively. In this work, these energies were calculated at the PBE0-D3(BJ)/def2-TZVP level of theory. e stability of transition metal complexes has also been studied based on the interaction energy (E int ) as described by equation (11). where is the energy of the complex, E MODP is the energy of the ligand, E M 2+ is the energy of metal ion, and E H 2 O is the energy of the water molecule.

Stability Analysis of Spin States.
Transition metal ions with incompletely filled d-orbitals can exhibit high and low  Table 1. e results show that apart from the Co(II) complex, the energies of the high spin complexes are lower than those of the low spin counterparts. As a result, only the triplet state for the Ni(II), Pd(II), and Pt(II) complexes, the doublet state for the Co(II) complex, and the quintet state for the Fe(II) complex have further been considered throughout this work.

Infrared (IR) Vibrational Analysis.
In order to confirm the suitability of the DFT/B97-3c level of theory employed in the computation of IR frequencies, computed frequencies have been compared with their experimentally obtained counterparts from [22]. e spectra of MODP and its complexes are presented in Figure S1 of the accompanying electronic supplementary material (ESM). For better comparison of the theoretical and experimental IR frequencies, correlation equations (12) and (13) have been established based on the vibrational modes of the MODP and the complexes, respectively.
where ] cal and ] exp represent the calculated and experimental wavenumbers, respectively. e high correlation coefficients (R 2 � 0.998) for the MODP and (R 2 � 0.949) depict a good linear correlation between the calculated (unscaled) and experimental IR frequencies.
is confirms the suitability of the level of theory employed in these calculations.
Presented in Table 2 are some experimental and calculated gas phase IR frequencies of MODP and its M(II) complexes along with their probable assignments.
e results show that the IR spectra calculated for MODP and its transition metal complexes showed bands e metal-ligand vibration frequencies of these complexes are in agreement with the experimental values in [22]. e IR spectra of the complexes indicate that MODP behaves as a bidentate ligand and is bound to the metal ion via the N-atom of the azomethine group (C � N) and the O-atom of the phenolic group (OH) (after deprotonation).

Geometric Parameters.
e optimized geometries of MODP and its Cu(II), Ni(II), Co(II), Zn(II), Fe(II), Pd(II), and Pt(II) complexes are displayed in Figure 3, and their optimized Cartesian coordinates are provided in Tables S1-S16 of the accompanying electronic supplementary material (ESM). Figure 3 shows that an intramolecular hydrogen bond (H-bond) is formed between the hydroxyl group and one of the oxadiazole ring nitrogen in MODP, which certainly stabilizes the molecule. is H-bond also hinders free rotation about the C 6 -C 7 bond, thus maintaining the molecule in a somewhat planar structure that promotes π-conjugation.
Selected geometric parameters (bond lengths, bond angles, and dihedral angles) are listed in Table 3.
From the dihedral angles, MODP is improved upon metal ions complexation. From the metal-ligand bond Species ] exp represents experimental wavenumbers by [22]. ] cal represents wavenumbers calculated in this work.
Journal of Chemistry lengths, we conclude that there is effective metal-ligand binding in all complexes studied. e results from this table also show that the bond angles around the central metal ion in the complexes investigated are around 90°, with exception of the Cu(II) complex where large discrepancies are observed in some cases. is indicates a high planarity of MODP in its Co(II), Ni(II), Zn(II), Fe(II), Pd(II), and Pt(II) complexes, as well as octahedral geometries around the central metal ions. e Cu(II) complex of MODP apparently adopts a trigonal bipyramidal geometry around the Cu(II) ion, since one of the ancillary H 2 O ligands is clearly out of the coordination sphere and is rather attached to one of the MODP ligands through an H-bond. is may be attributed to the fact that the Cu(II) ion is not always stable in octahedral geometries, and the complexes of this ion often have a square or tetrahedral planar structure. is arrangement comes from the stretching of the octahedron, along the z-axis (Jahn-Teller effect). e Cu(II) complex is Jahn-Teller distorted because the bonds to the central Cu(II) ion along the z-axis (M 15 -O 16   by important factors such as the HOMO and LUMO energy levels of the donor and acceptor [18]. For effective charge transfer to occur, the HOMO and LUMO energies of efficient charge transporters must be higher than that of (6, 6)-Phenyl-C 61 -butyric acid methyl ester (PCBM) (HOMO � −6.1 eV and LUMO � −3.75 eV) and lower than that of Poly(3-hexylthiophene) (P3HT) (HOMO � −4.65 eV and LUMO � −2.13 eV) [16,42]. P3HT and PCBM are reference electron donor and acceptor materials, respectively. In order to assess the electron transfer possibilities of the molecules studied, their HOMO and LUMO energy levels were compared with those of PCBM (see Table 4). PCBM was chosen in this study as electron acceptor because   of its stability, its capacity to accept electrons, and its proper energy level [53] as earlier mentioned. e results reveal that all the studied molecules have higher HOMO and LUMO energies than PCBM. is indicates the ability of molecules to inject electrons into PCBM. e results also show that complexation modifies the HOMO and LUMO energies and the energy gaps. Complexation increases the HOMO energy values by a factor of about 0.50 eV and the LUMO energy values by about a factor of 0.25 eV, except for the LUMO energy of the Pt(II) complex (which is slightly reduced by 0.03 eV). Accordingly, the investigated molecules can be used as donor materials in organic solar cells. Moreover, photo-excited electron transfer from the studied molecules to PCBM may be sufficiently effective for applications in photovoltaic devices [40,54]. Also, the Egap values of the molecules are in the range 3.50-4.46 eV and decrease in the order: Consequently, metal complexation decreases the Egap of MODP. It is also observed that the Pt(II) complex has the smallest Egap (3.50 eV), while that of the Zn(II) complex (4.41 eV) is the highest.
is implies that the Pt(II) complex is the most reactive and least stable, while the Zn(II) complex is the most stable and least reactive. Consequently, the HOMO-LUMO electron transport is easier in the Pt (II) complex than in the other complexes. e frontier molecular orbital (FMOs) energy gaps, Egap, of MODP and its complexes are displayed in Figure 4

Open Circuit Voltage (V OC ).
One of the parameters that influence the efficiency of solar cells is Voc. In organic solar cells, Voc is a function of the HOMO energy of the donor and the LUMO energy of the acceptor, which are important in determining whether an effective transfer will take place between the donor and the acceptor [40]. e Voc values of the studied molecules (donor), as well as those of the prototypical donor material P3HT and the acceptor PCBM, are presented in Table 5 We also note that complexation reduces the Voc values of the complexes, suggesting that the transfer of electrons is easier from the compounds studied to PCBM. e results in Table 5 also show that the differences in the LUMO energy levels of the studied compounds and the acceptor range from 2.08 to 2.44 eV and are all higher than 0.3 eV.
is suggests that the transfer of photo-excited electrons from the molecules to PCBM is sufficient to be used in photovoltaic devices [54]. is assures an electron transfer between the donors and the acceptor. For all the molecules studied, the E L−L values obtained with respect to the PCBM acceptors are higher than those obtained with respect to the P3HT/PCBM pair. e E L−L values of all the molecules studied are higher than those of V OC and the binding energy (0.2 to 1.0 eV). erefore, MODP and its complexes are promising donors with respect to the PCBM acceptor. As a result, these compounds can be used as solar cells.

UV-Visible Absorption Spectra.
e electronic transitions within a material can greatly influence its photovoltaic applications. It is desirable that donor materials should have intense electronic absorption compared to those of the acceptor materials. In fact, the maximum absorption wavelength (λ max ) of a suitable material must satisfy the absorption property requirement, according to which, λ max should not exceed 920 nm [55]. For each compound, 10 excited states were calculated. e calculated absorption wavelengths (λ cal ), oscillator strengths (f), electronic transitions, the molecular orbital (MO) contributions to electronic transitions, and excitation energies (E) of the investigated molecules are presented in Table 6. All the parameters in Table 6 were obtained in the excited state where the oscillation strength is the greatest and corresponds to the most intense band of the UV-Vis spectrum.
Interestingly, the λ cal of all molecules studied is below 500 nm in both DMSO and the gas phase. λ cal of the molecules in DMSO is in the order: MODP < [Zn (MODP) 2 2 ]. From this trend, we note increases in the wavelengths after complexation in all complexes in both the gas phase and in DMSO solvent. is bathochromic effect can be explained by an increase in the number of double bonds in the complexes compared to MODP, since conjugated systems are important for studying electronic and absorption properties of compounds. Generally, light absorbance of a molecule is improved by higher number of conjugated double bonds [56]. e UV-vis absorption spectra of the ligand and its transition metal complexes as studied in gas phase and in DMSO are shown in Figure 5. e absorption spectra (see   Table 5. ese values correspond to the MO contribution H-3 ⟶ L + 2 and H ⟶ L, respectively, in DMSO solvent. ese values increase in the gas phase. For the Cu(II), Zn(II), Ni(II), and Fe(II) complexes in DMSO solvent, their absorption spectra as illustrated in Figure 5 exhibit maximum wavelength bands as follows: λ � 491.6 nm, λ � 302.2 nm, λ � 339.2 nm, λ � 325.7 nm corresponding to the oscillator strengths: f � 0.0760, f � 0.6186, f � 0.0003, f � 0.0576, respectively. ese bands are due to the MO contribution Table 6: Absorption spectra data calculated by TD-DFT methods for species at CAMB3LYP/def2-TZVP level of theory in the gas phase and in DMSO.   Table 5. An increase in these values is observed in the gas phase. e absorption wavelengths obtained in the gas phase and in DMSO show that the MODP and its complexes can be used in organic solar cells.

Reorganization Energies and Charge Transport.
e charge transport properties of the investigated molecules have been predicted from computed reorganization energies of the electron (λ e ) and the hole (λ h ). In Table 7, the reorganization energies of the molecules studied are compared with those of the prototypical molecules, N-N′-diphenyl-N-N′-bis (3-methylphenyl-(1,1′-biphyl)-4,4′-diamine (TPD) molecules (λ h � 0.290 eV) and Tris(8-hydroxyquinolinato) aluminum (Alq3) (λ e � 0.276 eV) [10,11] for hole and electron transfer, respectively. e results show that the λ e values of MODP and its Zn(II), Fe(II), Pd(II), and Pt(II) complexes are smaller than those of Tris (8-hydroxyquinolato) aluminum (Alq3) (λ e � 0.276 eV) indicating that the charge transfer rate of these molecules may be higher than that of Alq3. Accordingly, MODP and its Zn(II), Fe(II), Pd(II), and Pt(II) complexes may serve as good electron transporters in OLEDs owing to their small reorganization energies. It is clear from Table 7 that the coordination of MODP with the Zn(II), Fe(II), Pd(II), and Pt(II) ions increases its λ e , whereas coordination with Cu(II), Co(II), and Ni(II) reduces its λ e . It can also be inferred from this table that the λ h values of MODP and its Cu (II), Co (II), Ni (II), Pd(II), and Pt(II) complexes are smaller than those of TPD (λ h � 0.290 eV), showing that their charge transfer rates could be higher than those of TPD. erefore, MODP and its Cu (II), Co (II), Ni (II), Pd (II), and Pt (II) complexes are promising hole transporters for use in OLEDs. Again, MODP (λ e � 0.187 eV and λ h � 0.154 eV), the Pd(II) complex (λ e � 0.142 eV and λ h � 0.161 eV), and the Pt(II) (λ e � 0,121 eV and λ h � 0,257 eV) complexes are excellent hole and electron transporters, respectively. erefore, these materials are likely to function as ambipolar charge transporters (electron and hole) under the proper operating conditions for OLEDs and OSCs based on the smaller reorganization energies.

Charge Transfer Integral and Charge Transport.
e values of charge transfer integral for electrons (V e ) and holes (V h ), as calculated at the PBE0/def2-TZVP/D3BJ level theory, are presented in Table 8.
e results show that the V e and V h values of MODP are smaller than those of its complexes (except for the V e value of the Zn(II) complex). e smaller values of the V e and V h of MODP are confirmed by the small values of the reorganization energies lower than the prototype (see Table 7). MODP can function as an ambipolar electron/hole transport material. It is clear from

Charge Transport Rate and Stability.
e calculated values of chemical hardness (η) and interaction energies, ΔE int , of molecules studied herein are listed in Table 9.
It is clear from Table 9 and Figure 6 that the coordination of MODP with the Zn(II), Fe(II), Pd(II), and Pt(II) ions   A careful inspection of same Table 9 and Figure 6 clearly shows that the η value of MODP is very close to those of the complexes. It is also clear from Table 9 that the ligand has the largest η, which indicates that MODP is the most kinetically stable and least reactive molecule studied. e complexation does not significantly affect the stability of these molecules. Among the complexes, the Co(II) complex is the most kinetically stable, and this stability is apparently confirmed by its negative E int value. e relatively high η values of the Cu(II), Zn(II), Co(II), Ni(II), Fe(II), Pd(II), and Pt(II) complexes indicate their chemical stabilities. Moreover, the stabilities of these complexes are indicated by their negative E int values, which show that their formation from the ligand and the various metal ions are energetically favorable. In conclusion, metal coordination causes a slight decrease in the chemical hardness of MODP.

e Charge Transfer Rates and Mobility. According to
Marcus' theory, the rate of charge transport is important for determining the mobility of electrons and holes in a compound. e charge mobility is an important parameter to evaluate the performance of OLED and OSC materials [57,58]. Mobility can be approximately calculated by Einstein's relation [44,49]: With e electronic charge, d is the transport distance from the molecular center to center in a dimer, K B Boltzmann constant k (e/h) charge transfer rate, and T temperature at 298 k. e µ (e/hole) values and the charge transfer rate of the electron (k e ) and hole (k h ) of the studied molecules are presented in Table 10.
It is clear from  Figure 6: IP, EA energies, and η of the molecules studied at the PBE0-D3(BJ)/def2-TZVP level of theory. respectively.
ese results show that an increase in the charge transfer rate of electron leads to an increase in the mobility of electrons.
It can also be seen from Table 10 that an increase in the value of the charge transfer rate of the hole leads to an increase in the mobility of the hole. e results also show that complexation increases the values of ku and μu for all complexes studied except for those of the Fe(II) and Cu(II) complexes, which reduces these values. e Co(II) complex has a high charge transfer rate (1.153 × 1016 s −1 ) and mobility (358.742 cm 2 V −1 s −1 ). Based on these results, the complexes of Pd(II), Pt(II), and Fe(II) have high values of electron charge transfer rate and electron mobility, and those of Co(II) and Cu(II) have high values of charge transfer rate and hole mobility. us, these molecules are good candidates for the design of OLED and OSC devices in optoelectronics.

Conclusion
In this work, a theoretical study of the geometries, electronic properties, charge transfer, and absorption of the ligand and its metal complexes was carried out using the DFT and TD-DFT methods. e possibility of their application as components of organic solar cells and organic light-emitting diodes has been discussed. All optimized complexes have been found to adopt octahedral geometries around the central metal ion, except the Cu(II) complex, which adopts a trigonal bipyramidal geometry. e calculated absorption maxima are in the range 179-492 nm in the gas phase and 238-492 nm in DMSO. e results also reveal that MODP and the studied complexes can be used as donor materials in OSCs, thanks to their higher HOMO and LUMO energy values than those of PCBM, along with their high energy gaps, high open circuit voltage (V oc ), and their intense absorption bands. For all molecules studied, the values of E L−L obtained with PCBM acceptors are higher than 0.3 eV, which suggests that the transfer of photo-excited electrons from the molecules to PCBM is sufficiently efficient to be used in photovoltaic devices. e results also showed that MODP and its Pd(II) and Pt(II) complexes are good electron (0.187, 0.142, and 0.121 eV, respectively) and hole transporters (0.154, 0.161, and 0.257 eV, respectively) due to their very small reorganization energies, integral charge transfers, high charge transfer rates, and mobilities. erefore, these materials can be good transport materials for the manufacture of OLEDs and OSCs. Furthermore, the Zn(II) and Fe(II) complexes of MODP are found to be very good electron transporters, while the Cu(II), Co(II), and Ni(II) complexes are likely to be good hole carriers. e results of chemical hardness (η) and interaction energies ΔE int have shown that MODP and its Co(II) complex are amongst the most kinetically stable molecules, while all complexes have been found to be thermodynamically stable materials, and complexation does not significantly affect the stability of these molecules. erefore, MODP and its complexes are suitable charge transport materials (electron and hole) that can be used in the manufacture of organic light-emitting diodes (OLEDs) and organic solar cells (OSCs).

Data Availability
Data is available upon request.  Figure S1: IR spectrum of MODP and its complexes calculated at the B97-3c level in gas phase. Table S1. Cartesian coordinates of the gas phase optimized geometry of MODP19. Table S2: Cartesian coordinates in DMSO optimized geometry of MODP19. Table S3: Cartesian coordinates of the gas phase optimized geometry of Cu(II) complex 43. Table S4: Cartesian coordinates in DMSO optimized geometry of Cu(II) complex 43. Table S5: Cartesian coordinates of the gas phase optimized geometry of Ni(II) complex 43. Table S6: Cartesian coordinates of the gas phase optimized geometry of Ni(II) complex 43. Table S7: Cartesian coordinates of the gas phase optimized geometry of Co(II) complex 43. Table S8. Cartesian coordinates in DMSO optimized geometry of Co(II) complex 43. Table S9: Cartesian coordinates of the gas phase optimized geometry of Zn(II) complex 43. Table S10: Cartesian coordinates in DMSO phase optimized geometry of Zn(II) complex 43.  Table S11: Cartesian coordinates of the gas phase optimized geometry of Fe(II) complex 43. Table S12: Cartesian coordinates in DMSO phase optimized geometry of Fe(II) complex 43. Table S13: Cartesian coordinates of the gas phase optimized geometry of Pd(II) complex 43. Table S14: Cartesian coordinates in DMSO phase optimized geometry of Pd(II) complex 43. Table S15: Cartesian coordinates of the gas phase optimized geometry of Pt(II) complex 43.