Theoretical design of D-π-A system new dyes candidate for DSSC application.

Currently, dye-sensitized solar cells (DSSCs) are one of the energy technologies that has piqued the interest of researchers, due to their distinct characteristics such as excellent air stability, ease of synthesis and photovoltaic properties interesting. This work aims to study the optoelectronic properties and photovoltaic of six organic dyes based on phenothiazine (PTZ). The effects of bridging core modifications of recently synthesized PSB-4(R) molecule on structural, photovoltaic, electronic, and optical properties of D1-D6 are studied. Using the method Density Functional Theory (DFT) level of the B3LYP (Becke three-parameter Lee-Yang-Parr) exchange correlation functional with 6-31G (d, p) and time-dependent DFT (TD-DFT). According to the obtained results, optoelectronic properties and photovoltaic of the dyes, we can suggest that these designed molecules are better sensitizers as a candidate for the production of dye solar cells (DSSCs). This theoretical study paves the way for chemists to synthesize more efficient sensitizers for applications in dye solar cells.

Introduction

The strong gn class="Chemical">rowth of developing countries has reflected an increase in global energy consumption estimated at double current demand if policies, in terms of both energy saving and efficiency, are not more proactive. One of the most promising solutions for the future energy of humanity is the development of other forms (of energy) known as "renewable energies". This type of energy has immense advantage of being the natural origin, inexhaustible and non-polluting, since they do not emit greenhouse gases nor radioactive waste. Although solar cells made of silicon achieve yields that vary between 15% and 25.4% [1]. Their manufacturing costs and weight are high. This is an obstacle to their massive use by individuals and manufacturers. New alternatives to reduce the cost of these cells are (DSSCs) or “Grätzel cell”, which is a photoelectrochemical device that converts sunlight into electricity [2, 3]. With a relatively low cost [4, 5, 6, 7] and design that results in a wide range of molecular structures that can be accessed and tunable color [8, 9], their current yield is over 13 % [10]. Iodide/triiodide, electrolytes, and then a platinum counter electrode [11, 12, 13] compared to traditional P–N junction solar cells, they have considerably lower production costs. The most effective organic dyes (DSSCs) have the molecular architecture donor-bridge-acceptor (D–π–A); donors (D) are triarylamine or carbazole. In the literature, we can find many up-to-date examples of various D-π-A compounds based on PTZ framework with various linkers and significant photovoltaic performance [14, 15, 16, 17, 18]. In our study, the phenothiazine group (PTZ) is chosen as a donor while the groups derived from cyanoacrylic acid are considered an acceptor. This push-pull system facilitates the transfer of intramolecular load from the donor group to the acceptor through the π-conjugated units. In the past few years, several research groups provide an effort to improve the properties of D-π-A dyes [19, 20]. The principle of a dye solar cell operation is the absorption of the photosensitizer to the surface of the semiconductor. As light (in the visible spectrum) is absorbed, the dye is excited in its first excited state (S∗) and then oxidized as electrons are transferred into the conduction band (CB) of TiO2, the carboxylic group (COOH) at the end of the dye rings: the COOH forms a bond with the surface of the TiO2, network, giving it a proton. The oxidized dye is then diminished by a reducing agent and regenerated against the counter-electrode [21]. The separation of charges in (DSSC) is done by a process of electron transfer from the dye to TiO2, as well as the process of transporting holes from the oxidized dye to the electrolyte. The mechanism of electron transfer is strongly dependent on the electronic structure of the absorbent dye molecule and the adequacy of energy levels between the excited state of the dye and the CB of the TiO2. However, the separation of charges is directly related to the positioning of the energy levels between the molecule of the dye and other nanoparticles. The excited state S∗ of the dye is greater than the limit of the TiO2 conduction band; and (level 1) of this dye is lower than the chemical potential of the redox iodide/tri-iodide pair (I−/I3-)in the electrolyte [22]. These two are now an energetic driving force for the separation of electrons and holes. The semiconductor nanoparticle array functions not only as a large substrate surface for the dye molecules, but also as a transport medium for the electrons injected from them.

In another won class="Chemical">rk, we investigate the design of six light-absorbing dyes named Di (i = 1–6) with a donor–π–bridge–acceptor (D-π-A) structure. The optoelectronic properties are calculated using DFT approaches. Then, photovoltaic applications of these designed molecules are related to molecule PSB-4(R) through computational studies.

In this computational study, we have designed six new molecules (D1-D6) after modification in reference molecule synthesized. From the study, the molecules designed exhibit broad absorption bands in the visible region. The calculated parameters include the energy gap, the density of states (DOS), reorganization energies, electron injection driving force (ΔG inject), light-harvesting efficiency (LHE) curve, open-circuit voltage (Voc), natural bond orbital (NBO) analyses, (TDM) analysis, exciting state time, and frontier molecular orbital FMOs. According to the obtained results from these parameters, our work is useful to effectively guide efforts to synthesize these designed molecules (D1-D6) in the discovery of highly efficient (DSSCs).

Theoretical methodology

In the current study, we present a theoretical study of optical and electronic properties of the six molecules and the PSB-4(R) (Figure 1). This study has been performed without symmetry restriction on the ground state using DFT at the B3LYP level with the standard 6-31G (d, p) basis in the gas phase (Figure 2). To evaluate the properties, we made use of the Gaussian 09 package [23], with Gauss view 06 for visualization [24]. The spectrum of absorption of designed molecules has been computed through the TD-B3LYP/6-31G (d, p) basis set from the optimization structures. DFT with the Becke three-parameter exchange functional and Lee-Yang-Parr functional (B3LYP) [25, 26, 27, 28] and 6-31G (d, p) [29] have been employed. Electronic transitions (vertical excitation spectra including wavelengths and oscillator strength f), LHE, and electron injection driving force (ΔGinject) of designed molecules are calculated using TD-DFT [30] with method-B3LYP (CAM-B3LYP) [31]. The charge transfer upon electronic transitions has been examined by charge density difference (CDD) and means of the (TDM) implanted in Multiwfn 3.1 [32].

Figure 1

Chemical structure for the PSB-4(R) and designed molecules Di (i = 1–6).

Figure 2

Ground-state optimized geometries of the compounds (D1-D6) and R obtained by B3LYP/6-31G (d, p).

Chemical structun class="Chemical">re for the PSB-4(R) and designed molecules Di (i = 1–6).

Ground-state optimized geometn class="Chemical">ries of the compounds (D1-D6) and R obtained by B3LYP/6-31G (d, p).

Frontier molecular analysis and electronic properties

Analysis of electron density distn class="Chemical">ribution over the frontier molecular orbitals across different parts of the studied dyes indicates intramolecular charge transfer (ICT) of the ground state dyes have been performed (Figure 3). The analysis shows that the HOMO's of molecules D1, D2, D3, and D4 are localized in the donor groups phenothiazine (PTZ), while LUMO's are spread on the π-spacers and the central acceptor unit. For the designed molecules D5, D6 and reference (R), the HOMO electrons are largely spread on these molecules with maximum density on donor groups, and with less localization on the acceptor group. The LUMOs are primarily localized at the level of the cyanoacrylic acid units, with maximum density on π-spacers groups for D5 and D6. In contrast, the D1-D6 dyes exhibit a high LUMO charge density. This indicates the formation of a semiconductor electronic dye coupling resulted from high electron injection. These outcomes expose that the π-bridge groups substituted have clear effects on the distribution of (FMOs).

Figure 3

(FMOs) diagram of the molecular orbitals HOMO and LUMO of studied dyes D1-D6 and (R) at B3LYP/6-31G (d, p) level of theory.

(FMOs) diagram of the moleculan class="Chemical">r orbitals HOMO and LUMO of studied dyes D1-D6 and (R) at B3LYP/6-31G (d, p) level of theory.

To illustrate the split of the excited state exciton and chan class="Chemical">rge separation, it is useful to investigate the (EDDM) between the excited states and the ground ones. The electron density difference maps (Figure 4) reveal that a load transfer has been performed between the excited states and the ground ones for all the studied dyes. This is because of the electronic charge movement from the phenothiazine (PTZ) to the cyanoacrylic acid. The (EDDM) plots show the decrease of electron density localized on the main electron phenothiazine (PTZ) donor and the π-spacer segments, Additionally, the electron density on the π-linker and the cyanoacrylic acid units increased. For all of the dyes tested. The electron density increment zone (purple) is primarily focused on the acceptor region, with a minor contribution from the π-spacer; however, the electron density (blue) is primarily focused on the donor group, with a minor contribution from the π-spacer for dyes D1-D6. These results indicate that the (ICT) between the excited states and the ground-state have been performed for all the designed dyes and an absorption shift towards the red and a reduction band gap upon light absorption process have been observed, which is essential for an efficient dye-sensitized solar cells.

Figure 4

Electron density difference maps between the excited states and the ground ones of dyes Di (i = 1–6) and the reference R.

The EHOMO, ELUMO, and HOMO–n class="Chemical">LUMO energy gap of all designed molecules and the molecule (R) have been computed in the gas phase from the ground state optimized geometry. Their values are presented in Table 1. The theoretically calculated values of the HOMO and LUMO energies of the reference compound are -4.67 eV and -2.68 eV, while the associated energy gap is 1.99 eV. The experimental HOMO and LUMO energy levels of the reference molecule are -5.07 eV and -2.88 eV, along with the energy band gap (2.19 eV) [33]. It is noted that the theoretically obtained results are in good agreement with those obtained experimentally. These results reflect that the choice of the methods of density theory (DFT) and functional theory of density as a function of time (TD-DFT), cited above, is a good option for predicting the photovoltaic, electronic, and optical properties of the newly designed molecules. The calculated values of the associated energy Egap of the reference molecule show a slight deviation in comparison to those of the experimental values, since our calculation model does not consider the solid-state packing effect of the compounds.

Table 1

The energy of HOMOs, LUMOs, and Egap at B3LYP/6–31 G (d, p) level of theory for the compounds Di (i = 1–6) and R. In the parenthesis, the experimental values are given.

Dye	EHOMO	ELUMO	Egap
R	-5.01 (-5.07)	-2.88 (-2.88)	2.13 (2.19)
D1	-5.20	-3.33	1.87
D2	-5.15	-3.40	1.76
D3	-5.20	-3.25	1.96
D4	-5.16	-3.41	1.75
D5	-5.16	-3.73	1.43
D6	-5.03	-3.50	1.53

The energy of HOMOs, n class="Chemical">LUMOs, and Egap at B3LYP/6–31 G (d, p) level of theory for the compounds Di (i = 1–6) and R. In the parenthesis, the experimental values are given.

The NBO analysis of all the designed compounds (D1-D6) and refen class="Chemical">rence (R) at B3LYP/6-31G (d,p) level of theory.

Results and discussion

Structure of studied dyes

In our won class="Chemical">rk, the experimentally synthesized PSB-4(R) is taken as a reference. The structure of the reference molecule consists of phenothiazine (PTZ) as donor and cyanoacrylic acid as electron acceptor [25]. We have replaced the terminal 4,8-bis(5-isopropylthiophen-2-yl)-2,6-dimethylthieno [2,3-f]benzofuran group of the reference molecule with different reported π-spacers (D1-D6) moieties. The chemical structures of the six dyes (Figure 1). By carrying out the above-mentioned π-spacers (D1-D6) modifications in the reference molecule, six distinct dyes are obtained.

It is noticed that all our novel dyes show a nan class="Chemical">rrow energy gap compared to the reference molecule. This means that all the designed molecules exhibit easy charge transfer phenomenon from the ground state to the excited one. In comparison to the reference, this transfer favors efficient electron injection at the conduction band TiO2. D5 and D6 show the lowest energy gap value among all other designed molecules. This is due to higher extended conjugation in the π-linker segments inserted and the presence of different chromophore units (–C=O, –C=N, and –C–S–C–), in order to make higher aromatic stability through resonance so large charge transfer phenomenon which can be done from the ground state to the excited one. That's why D5 and D6 are recommended as the most effective candidates of the series for applications in organic dye solar cells.

The energy level diagram of the designed compounds D1, D2, D3, D4, D5, D6 and the reference compound R is displayed in Figure 5, along with CB and VB of TiO2 and the regeneration of the oxidized dyes by I3-/I - redox couple. It is seen that the LUMO values of the compounds D1 to D6 lie sufficiently above the CB edge of TiO2, providing a required driving force for successful injection of the electrons from the LUMO (excited-state) of the dye to the CB TiO2. It is also worth noting that the HOMO energy level of all engineered molecules is lower than the redox energy level of the I3-/I- redox shuttle. So, the impede back-electron transfer from the conduction band of TiO2 to the redox energy level of the redox shuttle (the recombination process) cannot be performed.

Figure 5

Energy level plot of the designed compounds D1 to D6 and R.

Energy level plot of the designed compounds D1 to D6 and n class="Chemical">R.

Milliken's population analysis

In orden class="Chemical">r to learn more, density distribution over the (FMOs) [34], a density of states (DOS) investigation has also been accomplished by DFT/B3LYP/6-31G (d, p) method. The outcomes from the DOS involvement (Figure 6) support the evidence revealed by the FMOs in (Figure 3). The density of states of the spectrum represents that HOMO density significantly spreads on the donor moiety and slightly on the π-linker group, while LUMO is contributed by the acceptor and π-linker moieties. The density of the state's spectrum indicates that the effect of the bridge group on distribution pattern is seen in DOS spectra. The major contribution of bridge groups in D5 and D6 towards FMOs is seen compared to other bridge groups. This analysis also predicts that there is a significant delocalization of electron density. Also, a significant amount of charge transfer occurs from the donor group to acceptor regions.

Figure 6

Density of states for the studied dyes D1-D6 and (R) at B3LYP/6-31G (d,p) level of theory.

Density of states for the studied dyes D1-D6 and (n class="Chemical">R) at B3LYP/6-31G (d,p) level of theory.

NBO analysis

In molecular systems, the (NBO) analysis is a cn class="Chemical">ritical parameter for examining charge transfer or conjugative interaction. It also provides a convenient platform for researching intramolecular and intermolecular bonding, as well as interactions between bonds [35].

To compare the pan class="Chemical">rtial charges on the donor moiety (phenothiazine), group π-spacers (bridge), and acceptor moiety (cyanoacrylic acid); the calculation of NBO has been carried out using the B3LYP/6-31G (d, p) method in the gas phase (Table 2). It is noted that all the molecules have positive charges that are situated in the donor part and the π-spacers part, while the negative charges are located in the acceptor part. Contrarily, the positive charges in the bridge inserted demonstrate that this unit behaves as a donor. The most positive charge on the donor of the compound D5 and the π-spacer for D6 is the indication of being the most effective electron donor, while the highest negative (NBO) charge for all the studied dyes is located on the acceptor moiety which acts as the best acceptor moiety. The Natural bond orbital analysis demonstrates the charge-separated state generated in the dye molecule and proves that the charge transfer can be produced between donor and acceptor moieties.

Table 2

The NBO analysis of all the designed compounds (D1-D6) and reference (R) at B3LYP/6-31G (d,p) level of theory.

dyes	DonatingGroup	π-spacer	Acceptor Group
R	0.019	0.108	-0.127
D1	0.045	0.058	-0.109
D2	0.045	0.083	-0.123
D3	0.037	0.083	-0.120
D4	0.035	0.076	-0.114
D5	0.064	0.067	-0.131
D6	0.045	0.110	-0.155

Simulated Uv-visible absorption spectra

The UV-Visible absorption spectra of the newly designed photovoltaic materials (D1-D6) and the molecule R are computed at the TD-CAM-B3LYP functional with the 6-31G (d,p) basis set as stated in the theoretical methodology. The spectrum of absorption for the molecules (D1-D6) set has been computed in the gas phase (vacuum) from the fully optimized 3-D structure. The corresponding results are tabulated in Table 3 and illustrated in Figure 7. The results obtained with TD-CAM-B3LYP level of theory are very similar to the reference molecules experimental ultraviolet spectrum [25]. This slight difference between the experimental values and theoretical ones may be attributed to the fact that the theoretical calculation has been carried out in the gas phase, while these experimental values are obtained from the solid-state molecules of dyes. According to Figure 7, the intense absorption peaks are obtained in the visible region. These results are a good argument of the use of organic materials in solar cell devices.

Table 3

Computed maximum absorption wavelengths (λmax, in nm), f, and major percentage contribution for all designed dyes D1 to D6 and molecule R using CAM -B3LYP/6-31G (d,p) method in the gas phase.

Dyes	Calcu.λmax (nm)	Eex (eV)	f	MO/character	(%)
R	463.35	2.6740	1.5586	HOMO→LUMO	(63%)
	383.12	2.6740	0.4411	HOMO-2 →LUMO	(35%)
	361.22	3.2340	0.7634	HOMO-3→LUMO	(35%)
D1	514.50	2.4081	1.5103	HOMO →LUMO	(63%)
	389.32	3.1824	0.0977	HOMO-1 →LUMO	(58%)
	336.58	3.6810	0.1855	HOMO-1→LUMO+1	(35%)
D2	569.26	2.1765	1.7884	HOMO→LUMO	(66%)
	403.57	3.0701	0.0138	HOMO-1 →LUMO	(59%)
	360.80	3.4340	0.0586	HOMO →LUMO+1	(46%)
D3	511.21	2.0289	1.3967	HOMO →LUMO	(63%)
	379.54	2.0724	0.0246	HOMO-1 →LUMO	(55%)
	347.73	2.8230	0.3830	HOMO-1 →LUMO+1	(39%)
D4	537.59	2.6740	1.4992	HOMO →LUMO	(61%)
	404.05	2.0151	0.0643	HOMO-1 →LUMO	(49%)
	365.79	2.9178	0.2026	HOMO-1 →LUMO+1	(33%)
D5	729.15	1.6992	0.9254	HOMO → LUMO	(85%)
	478.02	2.5919	0.0384	HOMO-1 →LUMO	(77%)
	403.36	3.0717	0.0011	HOMO-8 →LUMO	(81%)
D6	732.38	1.6918	0.8290	HOMO → LUMO	(92%)
	446.12	2.7773	0.2193	HOMO-1 →LUMO	(71%)
	420.73	2.9449	0.6509	HOMO →LUMO+1	(72%)

Figure 7

Stimulated absorption spectra of molecules D1-D6 and R in the gas phase using CAM-B3LYP/6-31G (d,p) method.

Computed maximum absorption wavelengths (λmax, in nm), f, and majon class="Chemical">r percentage contribution for all designed dyes D1 to D6 and molecule R using CAM -B3LYP/6-31G (d,p) method in the gas phase.

Stimulated absorption spectn class="Chemical">ra of molecules D1-D6 and R in the gas phase using CAM-B3LYP/6-31G (d,p) method.

Table 3 presents the recorded maxima absorption values of all the studied molecules. The value of absorption maxima of all designed molecules (D1-D6), the reference R, the oscillator strength (ƒ), excitation energies (Eex), and percentage electron transport contributions (%ETC) are theoretically calculated using the method TD-CAM-B3LYP/6-31G (d, p). The theoretical results reveal that the absorption spectra (calculated λmax) of the reference molecule is 484 nm, which is very close to its experimental value 518 nm [33], ensuring a better choice of the used theoretical methods [14, 17]. The molecules absorption spectra show the presence of a large peak between 400 and 1100 nm and narrower peaks in the 250–450 nm wavelength region characterized as a typical π→π∗ transition (HOMO to LUMO). These results suggest that all compounds D1 to D6 have only one band in the visible region (λmax > 400 nm). Among all studied dyes, the highest λmax is exhibited by D5 and D6. This may be due to the increase in the length of conjugation through the π-spacer part. So, the change of the 4,8-bis(5-isopropylthiophen-2-yl)-2,6-dimethylthieno [2,3-f] benzofuran unit by the different π-spacers particularly for molecules D5 and D6 increases the maximum absorption wavelength. Notably, compared to D5 and D6, the compound D4 has a maximum absorption shifted by about 200 nm towards shorter wavelength than D5 and D6. This may be due to a result of the more electron deficient nature of the naphtho [1,2-c:5,6-c']bis ([1,2,5]thiadiazole) acceptor, their planar and rigid backbones usually indicate a close intermolecular interaction [1,2,5]thiadiazole for high-performance polymer solar cells. J. Am. Chem. Soc.. 2011 ">36]. This shows that the designed molecules have excellent absorption properties as compared to the reference molecule.

Overall efficiencies

The perfon class="Chemical">rmance of DSSC solar cells relates to broader coverage of the sunlight into electricity. The proposed compounds (D1-D6) act as sensitizers in (DSSCs) photovoltaic parameters such as the open-circuit voltage describe these devices (Voc), which can be calculated according to Eq. (1) [37]:

For the designed molecules, the obtained results from the Voc range from 0.27 to 0.75eV (Table 3). These values are sufficient for a better efficient injection of electrons from ELUMO of the dye to the CB of TiO2. The light-harvesting efficiency (LHE (λ)) can be determined by Eq. (2) [38, 39]:

The more negative value of ΔGinject is attn class="Chemical">ributed to a higher driving force for electrons injected into the semiconductor substrate from the dyes excited states. ΔGinject can be calculated from the difference between excited state oxidation potential of the dye molecules CB of the TiO2. It is determined by Eq. (3) [40, 41]:where Edye∗ denotes the potential oxidation energy of the dye in their excited state. It can be obtained from Eq. (4) [42]:where is the ground state oxidation potential of the dye. is the lowest absorption energy associated with the photo-induced (ICT).

Another electn class="Chemical">rochemical parameter can be optimized to obtain a higher photoelectric conversion. The free energy of the regeneration of the dye (ΔGreg) can be calculated by Eq. (5) [43]:whereis defined as the ground state oxidation potential and is equal to the negative value of the HOMO energy of the isolated dyes, and I−/I3- is the iodide/tri-iodide redox potential (−4.8 eV) [44]. Table.4 shows that the findings driving force (ΔG inject) of the designed dyes are negative. ΔG inject is shown to be -1.2, -1.03, -0.83, -1.51, -0.54, and -0.67, for Di (i = 1–6), respectively. This predicts spontaneous electron injection from the excited dye molecules to the CB edge of the TiO2. It is also noted that the designed molecules demonstrate thermodynamically favorable regeneration ability as all ΔGreg values are negative. This means that the dyes excited state is lower than its potential of the I−/I−3 redox couple.

Table 4

The calculated values of (ELUMO (eV), EHOMO (eV), and VOC (eV), ΔGinject, ΔGreg, and LHE for the dyes D1 to D6 and reference R calculated at DFT/B3LYP/6-31G (d,p).

dyes	EHOMO	ELUMO	E00	Edyes	Edyes∗	ΔGinject	ΔGreg	LHE	Voc
R	-5.01	-2.88	2.67	5.01	2.34	-1.66	-0.21	0.97	1.12
D1	-5.20	-3.33	2.67	5.20	2.80	-1.2	-0.40	0.97	0.67
D2	-5.15	-3.40	2.14	5.15	2.97	-1.03	-0.35	0.98	0.60
D3	-5.20	-3.25	2.03	5.20	3.17	-0.83	-0.40	0.96	0.75
D3	-5.16	-3.41	2.02	5.16	2.49	-1.51	-0.36	0.97	0.59
D5	-5.16	-3.73	1.69	5.16	3.46	-0.54	-0.36	0.88	0.27
D6	-5.03	-3.50	2.02	5.03	3.33	-0.67	-0.23	0.85	0.50
TiO2	-7.2	-4

Table.4 shows that the oscillator strengths are 1.5103, 1.7884, 1.3967, 1.4992, 0.9254, and 0.8290 for D1, D2, D3, D4, D5, and D6, respectively. Therefore, the calculated LHE is in the order: D2>D1 = D4 = R > D3> D5> D6. This shows that our entire novel designed molecules have the highest solar light utilization efficiency compared to reference molecule to obtain the higher Jsc.

The reon class="Chemical">rganization energy λtotal has great importance to evaluate the energy penalty resulted from the molecular arrangement that occurs during photoexcitation. Therefore, the calculation of λtotal is also important to analyze the kinetics of electron injection from the LUMO (excited state) of the dye to the CB of TiO2. Marcus semiclassical theory (Eq. (6)) depicts the rate of electron transfer (KET), which is dependent on the λtotal. The relation between them is given in Eq. (6) [45]:where KB denotes Boltzmann constant; ℏ is the Planks' constant, V represents the electronic coupling constant between initial states and final ones, and T is the absolute temperature. This shows that the electron transfer rate constant KET depends only on λtotal which is a sum of the electron (λe) and hole (λh) reorganization energies whose values can be estimated from (Figure 8) [46] and calculated by Eqs. (7) and (8) respectively [47]:where e and denote the energy of the neutral molecule in the optimized anionic (cationic) geometry, represents the energy of the neutral molecule at the ground state, and and are the energy of anionic and cationic calculated from the optimized anion (cation) geometry, respectively. For low reorganization energies, charge transfer is fast and efficient.

Figure 8

Energies of the neutral compounds at the ground state and the cationic (anionic) systems.

Energies of the neutn class="Chemical">ral compounds at the ground state and the cationic (anionic) systems.

The rate of chan class="Chemical">rge transfer increases with the decreasing of the reorganization energy λtotal. According to the results obtained all the dyes Di (i = 1–6) exhibit low reorganizational energy. Therefore, the dyes (D1-D6) have greater electron transportability between the cyanoacrylic acid (acceptor) and phenothiazine (PTZ) (donor) moieties. The calculated λtotal of the compounds are closer and more increased in order: D4< R < D5 < D6 < D3< D2< D1. Therefore, there is an almost equal load transfer capacity. D4, D5, and D6 possess the lowest λtotal. The lower dye regeneration driving force and high negative NBO charge on the cyanoacrylic acid moiety affect the photovoltaic performance of D4, D5 and D6. The dye D1 has the largest λtotal. As a result, dyes D4, D5, and D6 exhibit a favorable Jsc that can be used for future efficient OSCs.

Excited-state lifetime (t) has an important influence on the chan class="Chemical">rge transfer properties of the material. In comparison to the reference molecule, the dyes are engineered to have a longer excited state lifetime. Therefore, the cationic form of the dye molecule remains the most favorable state for efficient charge transfer. The lifetime can be calculated by Eq. (9) [48]:where f denotes the oscillator strength of the excited state and Eex represents the excitation energy, Table 5 shows the lifetime of the excited states with the order of D6 > D5 >D3 > D2 > D1> R = D4, which means that the insertion of novel π-spacers group increases conjugated length. This is helpful to increase the lifetime of excited states compared to the introduction of 4,8-bis(5-isopropylthiophen-2-yl)-2,6 dimethylthieno [2,3-f] benzofuran unit.

Table 5

Reorganization energies, λtotal values (eV) and excited-state lifetime t (ns) of the designed molecules (D1-D6) and R.

dyes	λe	λh	λtotal	t (ns)
R	0.32	0.25	0.57	0.14
D1	0.33	0.31	0.64	0.17
D2	0.32	0.31	0.63	0.18
D3	0.31	0.31	0.62	0.26
D4	0.26	0.29	0.55	0.14
D5	0.25	0.33	0.59	0.56
D6	0.27	0.32	0.60	0.62

Reon class="Chemical">rganization energies, λtotal values (eV) and excited-state lifetime t (ns) of the designed molecules (D1-D6) and R.

Transition density matrix (TDM)

The (TDMs) analysis is a technique for deciphen class="Chemical">ring and understanding electronic transformation processes in dye-sensitized solar cells [49, 50, 51]. In fact, this parameter provides information on the distribution of the interactions between donor group and acceptor, and also on the linked electron-hole pairs and identify their delocalization and coherence lengths [52]. The method DFT with 6-31G (d, p) level has been used to examine the TDM diagrams of all molecules. To eliminate the complications, we divided the dyes into three parts: Donor (D), Acceptor (A), and Bridge (B), (Figure 9).

Figure 9

Transition Density Matrix (TDM) of the compounds D1-D6 and R, Donor (D), Acceptor (A) and Bridge (B) calculated at DFT/B3LYP/6-31G (d,p).

Transition Density Matn class="Chemical">rix (TDM) of the compounds D1-D6 and R, Donor (D), Acceptor (A) and Bridge (B) calculated at DFT/B3LYP/6-31G (d,p).

According to the TDM diagrams, it can be seen that the designed molecules D1-D6 have a uniform distribution of electrons throughout the molecule. Coherence successfully passes from the donor to the acceptor unit with an effective contribution to the bridge unit, which acts as a facilitator for the transfer of electrons without trapping them. The order of interaction coefficient of the studied molecule between the donor and acceptor parts is D4 < D1< D2< D3< R < D5< D6. This order of transition density matrixes analysis provides evidence that newly designed dyes D5 and D6 show better results than the other molecules.

Exciton binding energy (Eb), which can be used to judge the exciton ability to dissociate, is a cn class="Chemical">rucial parameter for estimating the optoelectronic properties of (DSSCs) [52, 53, 54, 55]. For the dyes (D1-D6), the energy Egap is defined as the difference between the energies of HOMO and LUMO orbitals, while the optical gap energy is the first excited state energy. Binding energy Eb of the electron-hole pair of designed molecules D1 to D6 can be calculated with the help of Eq. (10) [56, 57, 58].

According to the results expressed in Table 6, the designed molecule D4 has lower exciton binding energy (Eb) than the other molecules. This means that all molecules exhibit an important charge transfer rate. Consequently, a high current charge density (Jsc) is compared to the other molecules. Besides, the sequence of the Eb values for the designed dyes is D4 < D1< D2< D3< R < D5< D6. These results are in line with previous findings (TDMs) analysis.

Table 6

Calculated energy gap Egap, first singlet excitation energy (Eopt), and the exciton binding energies (Eb).

dyes	Egap (eV)	Eopt (eV)	Eb (eV)
R	2.13	1.93	0.20
D1	1.87	1.77	0.10
D2	1.76	1.62	0.14
D3	1.96	1.79	0.17
D4	1.75	1.71	0.04
D5	1.43	1.12	0.31
D6	1.53	1.15	0.38

Calculated energy gap n class="Gene">Egap, first singlet excitation energy (Eopt), and the exciton binding energies (Eb).

Conclusion

In this work, we have designed and studied new six dyes D-π-A type and studied fon class="Chemical">r photovoltaic applications by employing the DFT functional and TD-DFT with a 6–31 G (d, p) basis. The designed molecules contain phenothiazine (PTZ) and cyanoacrylic acid by substituting different π-spacer groups at the peripheral sites of PSB-4(R) to examine the bridge effect units on the photovoltaic performance of solar energy driven devices. The designed molecules show a reduction in the energy gap (1.43–1.96 eV). The maximum absorption (732 nm and 729) has been observed for the molecules D6 and D5 with the lowest Egap (1.53 eV and 1.43eV) due to the insertion effect of the different π-spacers. Besides, the molecules D5 and D6 have the highest excited-state lifetime (t) compared to the other molecules. In addition, our designed dyes are good for electron mobility, because the binding energy (Eb) of the-designed molecules (D1 to D6) is lower than the synthesized reference molecule, which is in agreement with the results of the transition density matrix (TDM). Based on this study, the modifications π-spacer groups of the molecule reference affect the performance of designed dyes. In a nutshell, the newly designed dyes D1-D6, especially D5 and D6, should be synthesised in order to create highly effective dye-sensitized solar cells (DSSCs).

Declarations

Author contribution statement

R. Kacimi: Pen class="Chemical">rformed the experiments; Analyzed and interpreted the data; Wrote the paper.

M. Raftani, T. Abn class="Chemical">ram, A. Azaid, H. Ziyat: Analyzed and interpreted the data.

M. Bouachrine: Conceived and designed the expen class="Chemical">riments; Wrote the paper.

L. Bejjit, M. N. Bennani: Contributed n class="Chemical">reagents, materials, analysis tools or data.

Funding statement

This resean class="Chemical">rch did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.

Data availability statement

Data included in article/supp. maten class="Chemical">rial/referenced in article.

Declaration of interests statement

The authors declare no conflict of interest.

Additional information

No additional information is available fon class="Chemical">r this paper.