Designing Electron-Deficient Diketone Unit Based Non-Fused Ring Acceptors with Amplified Optoelectronic Features for Highly Efficient Organic Solar Cells: A DFT Study

Organic solar cells (OSCs) made of electron-acceptor and electron-donor materials have significantly developed in the last decade, demonstrating their enormous potential in cutting-edge optoelectronic applications. Consequently, we designed seven novel non-fused ring electron acceptors (NFREAs) (BTIC-U1 to BTIC-U7) using synthesized electron-deficient diketone units and reported end-capped acceptors, a viable route for augmented optoelectronic properties. The DFT and TDDFT approaches were used to measure the power conversion efficiency (PCE), open circuit voltage (Voc), reorganization energies (λh, λe), fill factor (FF), light harvesting efficiency (LHE) and to evaluate the potential usage of proposed compounds in solar cell applications. The findings confirmed that the photovoltaic, photophysical, and electronic properties of the designed molecules BTIC-U1 to BTIC-U7 are superior to those of reference BTIC-R. The TDM analysis demonstrates a smooth flow of charge from the core to the acceptor groups. Charge transfer analysis of the BTIC-U1:PTB7-Th blend revealed orbital superposition and successful charge transfer from HOMO (PTB7-Th) to LUMO (BTIC-U1). The BTIC-U5 and BTIC-U7 outperformed the reference BTIC-R and other developed molecules in terms of PCE (23.29% and 21.18%), FF (0.901 and 0.894), normalized Voc (48.674 and 44.597), and Voc (1.261 eV and 1.155 eV). The proposed compounds enclose high electron and hole transfer mobilities, making them the ideal candidate for use with PTB7-Th film. As a result, future SM-OSC design should prioritize using these constructed molecules, which exhibit excellent optoelectronic properties, as superior scaffolds.

These findings suggest that NFREAs created by incorporating electronegativity atoms with intramolecular noncovalent interactions (NCI) can significantly contribute to improved photovoltaic performance. Because of their elevated HOMO energy levels, NFREAs containing non-fused electron-donating rings or alkyloxy-substituted benzene have typically been restricted to pairing with polymer donors such as PM6, PTQ10 having wide bandgap and lower-lying HOMO level. As a result, the energy profile of hole transfer from NFREAs to donors is unfavorable [31].
Thus, electron-deficient D moiety units are needed to lock the molecular geometry through intramolecular NCIs and achieve low-lying HOMO energy levels. Taking this into account, the need for electron-withdrawing properties and the desire for a noncovalent intramolecular interaction appear to be satisfied by electron-deficient diketone units. The O/S intramolecular interactions in thieno [3,4-c]pyrrole-4,6-dione (TPD) NFREAs were reported by Chen's group [32]. TPD core exhibited additional quinoidal character than that of commonly encountered alkyloxy-substituted benzene analogues, resulting in improved coplanarity of NFREAs skeleton and larger p-electrons delocalization. When combined with wide bandgap polymer donors, electron affinity is increased by the diketone groups' which prevents the HOMO levels from being efficiently boosted. Additionally, the molecular planarity is enhanced by the O/S intramolecular NCI, which also boosts the ICT effect. Thus, when using PBDB-T, the TPD-based device attained a PCE of 10.12%.
The 2,20-bithiophene-3,30-dicarboxyimide (BTI) unit, one of many electron-deficient diketone units, was used as the basis of many superior polymer semiconductors [33][34][35][36]. The π-electron delocalization in BTI units gives them potent electron-drawing properties. For an all-polymer solar cell, for instance, Guo's team reported a PCE of 14.3% using an acceptor polymer holding BTI entity [37]. The presence of these characteristics suggested that the BTI unit could be used to build NFREAs for efficient OSCs. In light of the preceding, Luo et al. reported a pair of straightforward electron acceptors (BTIC-4F and TPDC-4F NFREAs) comprised of BTI and TPD electron-deficient diketone cores, correspondingly, linked to the DFIC terminals via a CPDT bridge [38].
End-capped modification is a good way to design novel OSC materials with outstanding photovoltaic properties [39][40][41][42][43]. Our work suggests that building NFREAs with end-capped modifications and electron-deficient diketone units is viable for performant and economical OSCs. We also shed light on the sophisticated molecular design behind efficient yet simple electron acceptors. Using end-capped engineering, which replaces the terminal units with new electron-withdrawing groups, we created seven new simple electron acceptors (BTIC-U1 to BTIC-U7). The photovoltaic properties of proposed compounds are predicted using DFT and TDDFT computations. This theoretical perspective will allow computational scientists to utilize this suitable approach for designing novel OSC materials and experimental scientists to synthesize these extremely effective predicted molecules for solar cell applications.

Method Selection and Optimized Geometries
After the geometrical optimization of BTIC-R, its time-dependent-DFT calculations were made in the excited state, specifically the cited chloroform solvent. Chloroform was used as the solvent of choice due to its implementation on the reference molecule in the cited literature. Utilizing five distinct functionals (PBEPBE, ωB97XD, MPW1PW91, CAM-B3LYP, B3LYP) and a basis set of 6-31G*, the BTIC-R molecule's absorption measurements were analyzed. The reason for using these functional is to investigate the functional that could lead to desirable results in terms of λ max . Figure 2 demonstrates that the λ max level of the BTIC-R molecule for these functionals is 769 nm, followed by 558 nm, 713 nm, 537 nm, and 1034 nm in that order. The λ max of the BTIC-R molecule was measured experimentally to be 696 nm [38] (see Figure 2), and DFT calculations based on λ max values suggest that the MPW1PW91/6-31G* has the highest concordance with the experimentally obtained λ max . Therefore, this functional was selected to carry out the research's auxiliary theoretical calculations. As a result, the functional and basis set were applied to optimize BTIC-U1 to BTIC-U7 molecules. It's common knowledge that a molecule's structure significantly affects its optoelectronic properties. Delocalization of π-electrons allows for a more rapid and efficient charge transfer between BTIC-R and BTIC-U1 to BTIC-U7 molecules, resulting in a longer conjugation.

Method Selection and Optimized Geometries
After the geometrical optimization of BTIC-R, its time-dependent-DFT calculations were made in the excited state, specifically the cited chloroform solvent. Chloroform was used as the solvent of choice due to its implementation on the reference molecule in the cited literature. Utilizing five distinct functionals (PBEPBE, ωB97XD, MPW1PW91, CAM-B3LYP, B3LYP) and a basis set of 6-31G*, the BTIC-R molecule's absorption measurements were analyzed. The reason for using these functional is to investigate the functional that could lead to desirable results in terms of λmax. Figure 2 demonstrates that the λmax level of the BTIC-R molecule for these functionals is 769 nm, followed by 558 nm, 713 nm, 537 nm, and 1034 nm in that order. The λmax of the BTIC-R molecule was measured experimentally to be 696 nm [38] (see Figure 2), and DFT calculations based on λmax values suggest that the MPW1PW91/6-31G* has the highest concordance with the experimentally obtained λmax. Therefore, this functional was selected to carry out the research's auxiliary theoretical calculations. As a result, the functional and basis set were applied to optimize BTIC-U1 to BTIC-U7 molecules. It's common knowledge that a molecule's structure significantly affects its optoelectronic properties. Delocalization of π-electrons allows for a more rapid and efficient charge transfer between BTIC-R and BTIC-U1 to BTIC-U7 molecules, resulting in a longer conjugation.

Quantum Mechanical Descriptors
Electronic properties, absorption and charge transmission of investigated molecules are profoundly affected by the energies of frontier molecular orbitals (LUMO and HOMO) [44][45][46][47]. The valence band, or HOMO, is where electrons are given away, while the conduction band, or LUMO, is where they are accepted [48,49]. The energy gap (Egap) is a standard identifier of photovoltaic (PV) and solar cell (SC) devices; it is a measure of how much energy is needed to split an atom's electrons apart [50]. Band gaps in organic solar cells can be narrower to increase efficiency [51]. Efficient and reliable PV devices require molecules with the smallest possible band gap value. Band gaps (Egap), HOMO, and LUMO energies for BTIC-R and BTIC-U1-U7 are shown in Table 1 below.

Quantum Mechanical Descriptors
Electronic properties, absorption and charge transmission of investigated molecules are profoundly affected by the energies of frontier molecular orbitals (LUMO and HOMO) [44][45][46][47]. The valence band, or HOMO, is where electrons are given away, while the conduction band, or LUMO, is where they are accepted [48,49]. The energy gap (E gap ) is a standard identifier of photovoltaic (PV) and solar cell (SC) devices; it is a measure of how much energy is needed to split an atom's electrons apart [50]. Band gaps in organic solar cells can be narrower to increase efficiency [51]. Efficient and reliable PV devices require molecules with the smallest possible band gap value. Band gaps (E gap ), HOMO, and LUMO energies for BTIC-R and BTIC-U1-U7 are shown in Table 1 below. In Figure 3, HOMO-LUMO plots for every molecule (built from their optimized geometries) and their corresponding E gap values are portrayed. The ground state's (HOMO) density of charges is centralized on the molecule's D and A1 parts, demonstrating that charge can be efficiently transferred from the ground to excited states. In the excited state, electron density also moves towards end-group acceptors (LUMO). In addition, in both FMOs under consideration, the molecule's predominantly planar topology is indicated by the fact that the charge density is spread out over nearly the entire compound. The HOMO energy of the BTIC-R molecule was measured to be −5.81 eV, and the LUMO energy of the same molecule was estimated to be −3.57 eV. This molecule also displayed an Egap of 2.24 eV, which is higher than all proposed compounds. The HOMO of the BTIC-U1 molecule was shown to have the lowest energy of all the molecules stud- The HOMO energy of the BTIC-R molecule was measured to be −5.81 eV, and the LUMO energy of the same molecule was estimated to be −3.57 eV. This molecule also displayed an E gap of 2.24 eV, which is higher than all proposed compounds. The HOMO of the BTIC-U1 molecule was shown to have the lowest energy of all the molecules studied, which indicates that the HOMO of the BTIC-U1 molecule is the most stable of all the molecules. Nitro groups are present in the end-capped acceptor groups of BTIC-U1 molecules; these groups have the potential to contribute to the molecules' stability.
In addition, The BTIC-U6 molecule possesses the HOMO with the second highest degree of stability, right behind the BTIC-U1 molecule. This is because the HOMO of BTIC-U6 molecule is lesser than any other molecule, coming in at −6.05 eV, and it also has a significantly higher amount of cyano units than any other molecule. In general, the order of decreasing HOMO energy levels across all investigated molecules is BTIC-U5 > BTIC-U7 > BTIC-R > BTIC-U4 > BTIC-U3 > BTIC-U2 > BTIC-U6 > BTIC-U1, and the declining array of the LUMOs for BTIC-R and BTIC-U1 to BTIC-U7 molecules is as follows: Therefore, in both of the evaluated FMOs, the energy levels present in BTIC-U1 and BTIC-U6, respectively, are incredibly low and rank second to last. Moving on, the E gap in the BTIC-U6 molecule is 2.12 eV, which is less than what is found in the BTIC-R molecule and every other designed molecule. Again, the BTIC-U1 molecule had the second smallest E gap of all the molecules tested (2.14 eV). In addition, all newly created molecules have an E gap lower than the E gap of the BTIC-R molecules. This demonstrates that the modified molecules are pretty capable of transferring charge effectively. According to the results obtained from Table 1, the decreasing order of all molecules that were studied concerning

Absorption Spectrum
The spectral analysis of BTIC-R and BTIC-U1 to BTIC-U7 molecules was performed using the previously determined MPW1PW91/6-31G*. Tables 2 and 3 present the assessed statistics for the solvent (chloroform) and gaseous phases, respectively.  As seen in Figure 4, the absorbance spectrum of each of the compounds that are the subject of this investigation spans the 350 nm to 1100 nm interval in both tested phases. The maximum wavelength (λ max ) of the BTIC-R molecule is 670 nm when it is in the gaseous phase, but it is 713 nm in chloroform (CHCl 3 ). Molecules in the BTIC-U1 to BTIC-U7 series have λ max values between 680 and 707 nm in the gaseous phase and between 723 and 761 nm in the CHCl 3 solvent. Compared to the BTIC-R molecule, all modified molecules exhibited a bathochromic shift in their λ max in solvent and gaseous phases, as shown by this UV-visible spectral analysis. And this bathochromic shift displayed by BTIC-U1 to BTIC-U7 molecules is between 18-37 nm in the gas and 10-48 nm in the solvent phase concerning the BTIC-R molecule.
Molecules 2023, 28, x FOR PEER REVIEW 8 of 26 As seen in Figure 4, the absorbance spectrum of each of the compounds that are the subject of this investigation spans the 350 nm to 1100 nm interval in both tested phases. The maximum wavelength (λmax) of the BTIC-R molecule is 670 nm when it is in the gaseous phase, but it is 713 nm in chloroform (CHCl3). Molecules in the BTIC-U1 to BTIC-U7 series have λmax values between 680 and 707 nm in the gaseous phase and between 723 and 761 nm in the CHCl3 solvent. Compared to the BTIC-R molecule, all modified molecules exhibited a bathochromic shift in their λmax in solvent and gaseous phases, as shown by this UV-visible spectral analysis. And this bathochromic shift displayed by BTIC-U1 to BTIC-U7 molecules is between 18-37 nm in the gas and 10-48 nm in the solvent phase concerning the BTIC-R molecule. In both the gas phase and the CHCl3, the increasing λmax sequence of BTIC-R and BTIC-U1 to BTIC-U7 molecules is as follows: BTIC-R < BTIC-U5 < BTIC-U4 < BTIC-U7 < BTIC-U3 < BTIC-U2 < BTIC-U1 < BTIC-U6 and BTIC-R < BTIC-U5 < BTIC-U7 < BTIC-U4 < BTIC-U3 < BTIC-U2 < BTIC-U6 < BTIC-U1, sequentially. All studied compounds in chlo-  In both the gas phase and the CHCl 3 , the increasing λ max sequence of BTIC-R and BTIC-U1 to BTIC-U7 molecules is as follows: BTIC-R < BTIC-U5 < BTIC-U4 < BTIC-U7 < BTIC-U3 < BTIC-U2 < BTIC-U1 < BTIC-U6 and BTIC-R < BTIC-U5 < BTIC-U7 < BTIC-U4 < BTIC-U3 < BTIC-U2 < BTIC-U6 < BTIC-U1, sequentially. All studied compounds in chloroform solvent exhibited a bathochromic shift relative to the gas phase due to the stabilisation of the polar excited state by a polar solvent. The shorter E gap of all the molecules resulted in a longer wavelength as a result of the inverse correlation between wavelength and energy.
The findings indicate that BTIC-U1 has a λ max of 761 nm in a solvent medium, pointing to a notable end-capped acceptor that may have affected its λ max . The oscillator strength (f) is a dimensionless factor that determines the photovoltaic cell's optical properties and the calculation of the strength of the radiation generated upon the electrical excitation of two energy levels. The excitation energy of OSC is also crucial to enhancing their photovoltaic or optoelectronic properties. All the molecules we've designed have low excitation energies relative to the standard BTIC-R molecule in the solvent phase, making them ideal for photovoltaic or optoelectronic applications. Tables 2 and 3 show data from the gas and solvent phases, highlighting significant HOMO→LUMO (≤95%) transitions or assignments for various molecules. These results apply to all BTIC-U1 to BTIC-U7 designed molecules and the BTIC-R reference molecules. This quality is also present in molecules engineered to have superior optoelectronic properties. Figure 4 depicts the simulated absorption spectra of molecules BTIC-U1 through BTIC-U7, and the reference molecule BTIC-R attained in chloroform and gas phases, illustrating each molecule's broader and more intense peak. The tabulated results in Tables 2 and 3 show that the λ max values for BTIC-U1 to BTIC-U7 are significantly impacted by the end-capped unit modifications, resulting in a 92% HOMO-LUMO transition, low excitation energies and a red shift in the UV/Visible spectrum. Thus, all of these molecules (BTIC-U1 to BTIC-U7) are superior to BTIC-R, the standard molecule, in their optoelectronic properties. Improving the photovoltaic capabilities and adjusting the photophysical characteristics of proposed compounds is easily accomplished through end-capped unit modification. That's why the newly designed molecules BTIC-U1 through BTIC-U7 have better optoelectronic properties than BTIC-R.

Light Harvesting Efficiency (LHE)
Every component that goes into the SCs needs to generate a charge following the light collection process, also referred to as the LHE. This LHE affects the device's solar efficiency, which in turn is affected by the oscillator strength due to its effect on short-circuit current generation. The LHE was determined for each of the compounds that were investigated by applying the Equation (1).
where the LHE is denoted by η λ and f denotes oscillator strength. In the following Table 4, the LHE values that were computed for molecules of BTIC-R and BTIC-U1 to BTIC-U7 in both the gas phase and the chloroform solvent are listed. The BTIC-U5 to BTIC-U7 molecules in the solvent phase displayed higher LHE values than BTIC-R due to the more significant oscillator strengths possessed by these molecules ( Figure 5). Because of the greater oscillation strength of BTIC-U6, LHE is at its highest level there. The findings indicate that the LHE of these small molecules is significantly affected by the terminal acceptor moieties.

Dipole Moment
There are few properties more crucial than crystallinity and solubility, and both are profoundly affected by dipole moment (µ). These characteristics play an essential role in establishing polarization situation in the required solvent for efficient OSCs compounds. Continuous charge transfer is enabled by the excellent polar solubility, superior crystallinity and tight molecular assembling of molecules with planar and ordered geometries containing significant dipole moments. A molecule's solubility and crystallinity in polar substances increase with increasing dipole moment, which improves charge transmission. In general, molecules with zero dipole moments cannot be dissolved by organic polar solvents like chloroform.

Dipole Moment
There are few properties more crucial than crystallinity and solubility, and both are profoundly affected by dipole moment (μ). These characteristics play an essential role in establishing polarization situation in the required solvent for efficient OSCs compounds. Continuous charge transfer is enabled by the excellent polar solubility, superior crystallinity and tight molecular assembling of molecules with planar and ordered geometries containing significant dipole moments. A molecule's solubility and crystallinity in polar substances increase with increasing dipole moment, which improves charge transmission. In general, molecules with zero dipole moments cannot be dissolved by organic polar solvents like chloroform.
However, this relationship is not constant because the charge transfers properties and solubility of each molecule depend on its specific molecular arrangement. Both the solvent and gaseous states are represented in Table 5 and Figure 6 by the μ values of the molecules BTIC-R and BTIC-U1 to BTIC-U7. In the solvent state, the order is, i.e., BTIC-U3 < BTIC-R < BTIC-U2 < BTIC-U7 < BTIC-U4 < BTIC-U5 < BTIC-U1 < BTIC-U6. The μ of the BTIC-U3 molecule was lower than the BTIC-R molecule in the solvent phase. However, it was greater than BTIC-R for the other molecules. Regarding solubility and polarity, BTIC-U6 had a huge dipole moment, possibly due to the cyano groups at its ends.

Molecules
μ (Gaseous Phase) μ (Solvent Phase) However, this relationship is not constant because the charge transfers properties and solubility of each molecule depend on its specific molecular arrangement. Both the solvent and gaseous states are represented in Table 5 and Figure 6 by the µ values of the molecules BTIC-R and BTIC-U1 to BTIC-U7. In the solvent state, the order is, i.e., BTIC-U3 < BTIC-R < BTIC-U2 < BTIC-U7 < BTIC-U4 < BTIC-U5 < BTIC-U1 < BTIC-U6. The µ of the BTIC-U3 molecule was lower than the BTIC-R molecule in the solvent phase. However, it was greater than BTIC-R for the other molecules. Regarding solubility and polarity, BTIC-U6 had a huge dipole moment, possibly due to the cyano groups at its ends.
Compared to the reference molecule, the incremental exchange of the µ was found in the range of 1.532-5.582 Debye (when measured in gas) and 0.866-7.307 Debye (when measured in solvent), as shown in Table 5. These findings validate the simulated materials' applicability for organic solar cells by proving their claim that crystalline behavior becomes more stable in the solvent. In-depth comparisons suggest that BTIC-U6, with its maximum µ value of 7.307 Debye, is the more crystalline material. An in-depth view of comparative analysis prompts BTIC-U6 as a more crystalline material because of its highest dipole moment values of 7.307 Debye. Additionally, all modeled entities have changed to have more dipole moments, guaranteeing accurate modeling for highly efficient polymer solar cells.  Compared to the reference molecule, the incremental exchange of the μ was found in the range of 1.532-5.582 Debye (when measured in gas) and 0.866-7.307 Debye (when measured in solvent), as shown in Table 5. These findings validate the simulated materials' applicability for organic solar cells by proving their claim that crystalline behavior becomes more stable in the solvent. In-depth comparisons suggest that BTIC-U6, with its maximum μ value of 7.307 Debye, is the more crystalline material. An in-depth view of comparative analysis prompts BTIC-U6 as a more crystalline material because of its highest dipole moment values of 7.307 Debye. Additionally, all modeled entities have changed to have more dipole moments, guaranteeing accurate modeling for highly efficient polymer solar cells.

Density of State
Studies of DOS are necessary for understanding how the total and partial DOS in the charge mobility of molecules validate the utilities of individual units (acceptor1, donor, and acceptor2) and the inclusive performance of investigated molecules. When analyzing the FMOs regarding the Mulliken charge density, it is essential to determine the order in which the FMOs are arranged. Calculations of the DOS for each of the molecules under investigation were carried out by employing the MPW1PW91/6-31G* methodology. The DOS plots have an x-axis that shows the total amount of energy and a y-axis that shows the relative intensity of that value. For this investigation, studied compounds are cut up into three different units-acceptor1, donor, and acceptor2-so that the researchers could examine the contributions made by each fragment to the FMOs. In Figure 7, black, blue, pink, and red lines represent the total contribution of the molecules (BTIC-R and BTIC-U1 to BTIC-U7) as a whole, contributions of the donor moieties, acceptor2 and acceptor1, re-

Density of State
Studies of DOS are necessary for understanding how the total and partial DOS in the charge mobility of molecules validate the utilities of individual units (acceptor1, donor, and acceptor2) and the inclusive performance of investigated molecules. When analyzing the FMOs regarding the Mulliken charge density, it is essential to determine the order in which the FMOs are arranged. Calculations of the DOS for each of the molecules under investigation were carried out by employing the MPW1PW91/6-31G* methodology. The DOS plots have an x-axis that shows the total amount of energy and a y-axis that shows the relative intensity of that value. For this investigation, studied compounds are cut up into three different units-acceptor1, donor, and acceptor2-so that the researchers could examine the contributions made by each fragment to the FMOs. In Figure 7, black, blue, pink, and red lines represent the total contribution of the molecules (BTIC-R and BTIC-U1 to BTIC-U7) as a whole, contributions of the donor moieties, acceptor2 and acceptor1, respectively.
In addition, the numerical contribution that each distinct sector of BTIC-R and BTIC-U1 makes to BTIC-U7 molecules is mentioned in Table 6. The donor group is an essential factor in the overall elevation of the ground state for all of the recently developed molecules and the BTIC-R molecule (HOMO). The contributions of both acceptor2 and acceptor1 to the HOMO energy level are similar but modest. However, acceptor1 contributes slightly more to the HOMO energy level than acceptor2.
Acceptor2 is responsible for most of the excited state's (LUMO) contribution. Based on these findings, it can be deduced that charges will be transferred via subsequent conjugation from the ground state of the donor fragment (electron-rich) part of molecules to the A2 unit in their excited state. Moreover, these results are consistent with the FMO analysis of BTIC-R and BTIC-U1 to BTIC-7, as shown in Figure 7. In addition, the numerical contribution that each distinct sector of BTIC-R and BTIC-U1 makes to BTIC-U7 molecules is mentioned in Table 6. The donor group is an essential factor in the overall elevation of the ground state for all of the recently developed molecules and the BTIC-R molecule (HOMO). The contributions of both acceptor2 and accep-tor1 to the HOMO energy level are similar but modest. However, acceptor1 contributes slightly more to the HOMO energy level than acceptor2.

Electrostatic Potential (ESP)
ESP shows the distribution of charges throughout a molecule in three dimensions and draws attention to the various locations within the molecule where electrons can be found. An ESP analysis was carried out on the molecules under consideration, in order to make a prediction regarding the reactive potential of molecular frameworks.
ESP maps depict the three-dimensional distribution of electronegative substances, lone pairs and electrons readily available for nucleophilic action. Red indicates an area with a high concentration of electrons on the ESP maps. Green represents a neutral region, and blue represents an area with a low concentration of electrons. Figure 8 displays the BTIC-R and BTIC-U1-BTIC-U7 ESP color maps. Because of the high electron density above these positions, ESP maps show the N and O atoms in the outer acceptor molecular regions as deep red. Likewise, the oxygen atoms in the core's acceptor region show up in red, making the existence of electron density in those regions obvious. Donor parts with methyl groups and thiophene rings are shown to have severe electron depletion in these ESP maps, as indicated by their bluish coloration. Because of the high electron density above these positions, ESP maps show the N and O atoms in the outer acceptor molecular regions as deep red. Likewise, the oxygen atoms in the core's acceptor region show up in red, making the existence of electron density in those regions obvious. Donor parts with methyl groups and thiophene rings are shown to have severe electron depletion in these ESP maps, as indicated by their bluish coloration.

Analysis of Charge Mobility
By computing the RE of BTIC-U1 to BTIC-U7 and BTIC-R molecules at the MPW1PW91 functional, we could examine the charge transfer efficiency between fragments of an electron D and A. The RE, which quantifies the potential for the transfer of charge from D to A parts, and is related to the electron and hole charge mobility, is a critical factor in the research and development of adequate materials for OSCs. According to Marcus' theory, there is an inverse connection between RE and the movement of these charges (electrons and holes). As a result, the reduction in the value of the RE will cause the charge transfer to be completed with greater effectiveness.
The internal reorganization energy of organic solar cell molecules has a substantial impact on the device's charge mobility, CT, and efficiency. Higher reorganization energy can lead to slower charge transfer rates, lower charge mobility, and decreased device efficiency. On the other hand, lower reorganization energy can result in faster charge transfer rates, higher charge mobility, and increased device efficiency. The internal reorganization energy of the organic solar cell molecule can also affect the stability of the device. If the reorganization energy is too high, the device may be unstable and prone to degradation. Therefore, optimizing the internal RE of the organic solar cell molecule is necessary to optimize the device's performance.
The RE is affected by various factors, including the geometry of cation and anion arrangements. The hole and electron exchange between acceptor and donor units is indicated by the cationic and anionic geometries, respectively. By solving Equations (6) and (7), we calculated the electron and hole RE values for each of the molecules under investigation; these values are summarized in Table 7 and Figure 9. For BTIC-R molecule, the RE of λ e and λ h are 0.1914 eV and 0.2108 eV, respectively. Except for BTIC-U4, all of the modified molecules improve upon the mobility of electrons compared to the reference BTIC-R by having lower RE energies for electrons. Based on their λ e values, the molecules can be arranged as follows: BTIC-U6 < BTIC-U1 < BTIC-U2 < BTIC-U3 < BTIC-U7 < BTIC-U5 < BTIC-R < BTIC-U4.
Based on these calculations, all other hypothetical molecules have a smaller value of λ e than BTIC-U4. In addition, BTIC-U6 conjugation framework and strong elctron withdrawing group in acceptor regions give it a minimum value of λ e (0.1217 eV), meaning it has a better rate of electron flow and can perform extraordinarily well in the process of charge transfer. Regarding acceptor molecules, the BTIC-U6 molecule has been shown to be the most effective Molecules with lower λ h such asBTIC-U2, BTIC-U3, BTIC-U5, BTIC-U6, and BTIC-U7, compared to BTIC-R, have a greater hole mobility than BTIC-R. In order of increasing hole moblility, molecules can be arrangedas follows: BTIC-U5 < BTIC-U6 < BTIC-U7 < BTIC-U2 < BTIC-U3 < BTIC-R < BTIC-U1 < BTIC-U4. To maximize hole mobility, the end-group acceptor effectively contributes to reducing the hole RE of the BTIC-U5 and BTIC-U6 molecules. In light of these findings, it is clear that the mobility of holes formed by exciton dissociations is maximized in BTIC-U5 and BTIC-U6. But it's easy to see that every molecule possesses λ e < λ h . It demonstrates that electrons in these molecules are more likely to be mobile and generate current than holes are. All compounds' lower λ e values relative to their λ h values suggest that they have improved acceptor capacities in modern organic photovoltaic cells. In addition, the low values of λ e and λ h for BRIC-U6 suggest that it may be the most stable molecule. In a nutshell, the improved performance of the newly proposed molecules in terms of internal reorganization energy can be interpreted to apply to the vast majority of them. As a result, they hold greater promise for enhancing charge mobility and enhancing performance. the reorganization energy is too high, the device may be unstable and prone to degradation. Therefore, optimizing the internal RE of the organic solar cell molecule is necessary to optimize the device's performance.
The RE is affected by various factors, including the geometry of cation and anion arrangements. The hole and electron exchange between acceptor and donor units is indicated by the cationic and anionic geometries, respectively. By solving Equations (6) and (7), we calculated the electron and hole RE values for each of the molecules under investigation; these values are summarized in Table 7 and Figure 9.

Chemical Reactivity Parameters
Chemical potential (µ), chemical hardness (η), chemical softness (S), and electrophilicity (ω) are all studied under the umbrella term "chemical-reactivity" to examine and compare the chemical behavior of BTIC-R and proposed molecules (BTIC-U1 to BTIC-U7). The term "electrophilicity" was coined to describe the phenomenon of an increased rate of electron flow in response to a decreased amount of energy. Table 8 shows that the values of all molecules, except for BTIC-U5, are higher than they were for BTIC-R, resulting in a good saturation of the designed structures with electrons. Compared to BTIC-R, the designed molecules have higher electrophilicity values, indicating that they function as acceptors. The chemical potential allows an acceptor to gain a net negative charge. To quantify the propensity of electrons to jump nuclei, the µ value is calculated for molecules with the BTIC-R and BTIC-U1 to BTIC-U7 hybrid electronic structure. The more significant negative potential of molecules in the range of BTIC-U1 to BTIC-U7 results in a marked increase in their reactivity. Chemical hardness, denoted by the formula η = 1/S, measures how well a material resists breakdown when subjected to chemical or physical stress. When an electron is added to the proposed structures, Table 8 below shows how much energy is lost (ensuring a negative chemical potential and keeping the hardness positive). The ability of the structures to transfer charge is explained by their pliability. The resistance to electron allocation decreases as S increases, so a material with a higher S value will have a softer surface (as shown in Figure 10).

Device Performance
The open circuit voltage (Voc) of any solar instrument can be measured to determine the instrument's photovoltaic activity, which is an essential step in deciding how the mechanism operates. When no other load is applied to an optical device, the VOC value represents the complete voltage range made available by the device. The voltage output of a photovoltaic cell is affected by different factors, including the temperature of the solar device, charge transfer, and incident light, amongst others. The valence band of the donor material is typically coupled with the conduction band of the acceptor material, which results in optimal voltage. In order to achieve the desired increased Voc values, the donor compound must have a lesser HOMO level, while the acceptor unit must have an upper LUMO level.
To achieve a high PCE, a photovoltaic system requires a high fill factor, which rises as Voc rises. Voc and intermolecular energy gaps are inextricably linked for predicting the efficiency of a solar cell. Donor polymer PTB7-Th's HOMO level is blended with BTIC-R LUMO level, and small acceptor molecules were claimed in this investigation, enhancing Voc (BTIC-U1 to BTIC-U7). PTB7-Th has a HOMO energy of −5.01 eV and a LUMO energy of −3.60 eV, suggesting it is a solid donor. In this study, Voc was determined statistically using Equation (2).

Device Performance
The open circuit voltage (Voc) of any solar instrument can be measured to determine the instrument's photovoltaic activity, which is an essential step in deciding how the mechanism operates. When no other load is applied to an optical device, the Voc value represents the complete voltage range made available by the device. The voltage output of a photovoltaic cell is affected by different factors, including the temperature of the solar device, charge transfer, and incident light, amongst others. The valence band of the donor material is typically coupled with the conduction band of the acceptor material, which results in optimal voltage. In order to achieve the desired increased Voc values, the donor compound must have a lesser HOMO level, while the acceptor unit must have an upper LUMO level.
To achieve a high PCE, a photovoltaic system requires a high fill factor, which rises as Voc rises. Voc and intermolecular energy gaps are inextricably linked for predicting the efficiency of a solar cell. Donor polymer PTB7-Th's HOMO level is blended with BTIC-R LUMO level, and small acceptor molecules were claimed in this investigation, enhancing Voc (BTIC-U1 to BTIC-U7). PTB7-Th has a HOMO energy of −5.01 eV and a LUMO energy of −3.60 eV, suggesting it is a solid donor. In this study, Voc was determined statistically using Equation (2).
In the equation, e stands for the molecular charge, which has a value of 1, whereas the inter-surface charge is a separate factor that typically has a value of 0.3. It was discovered that the Voc of BTIC-R is 1.143 eV, and the range of Voc values for BTIC-U1 to BTIC-U7 molecules was found to be 0.698-1.261 eV. When contrasted with the BTIC-R molecule, the BTIC-U5 and BTIC-U7 molecules have a Voc higher than the BTIC-R molecule. The theoretically calculated Voc levels for all of the compounds studied relative to PTB7-Th are depicted in Figure 11, whereas the statistically calculated data are provided in Table 9.
In the equation, e stands for the molecular charge, which has a value of 1, whereas the inter-surface charge is a separate factor that typically has a value of 0.3. It was discovered that the VOC of BTIC-R is 1.143 eV, and the range of VOC values for BTIC-U1 to BTIC-U7 molecules was found to be 0.698-1.261 eV. When contrasted with the BTIC-R molecule, the BTIC-U5 and BTIC-U7 molecules have a Voc higher than the BTIC-R molecule. The theoretically calculated Voc levels for all of the compounds studied relative to PTB7-Th are depicted in Figure 11, whereas the statistically calculated data are provided in Table  9.   The order of the molecules BTIC-R and BTIC-U1 to BTIC-U7 increased as follows: BTIC-U1 > BTIC-U6 > BTIC-U2 > BTIC-U3 > BTIC-U4 > BTIC-R > BTIC-U7 > BTIC-U5 According to this investigation, the BTIC-U5 molecule has the highest Voc value; as a result, it has the potential to be utilized in an effort to improve the PCE of OSCs.
To calculate the PCE of a PV system, one of the most crucial parameters is the fill factor (FF), which is directly related to the PCE. The Voc present at the interface between the acceptor molecule and the donor molecule significantly influences this aspect. To compute the FF of all of the molecules that we looked at, we used the equation shown 5 below.
The normalized Voc is eVoc K B T , where the standard charge (e) is always 1, K B stands for the Boltzmann constant, and the letter T stands for a constant temperature of 300 K. The Voc values for the compounds BTIC-R and BTIC-U1 to BTIC-U7 that were computed using computational methods are listed in Table 9 along with the normalized Voc and FF. The BTIC-R molecule has a normalized Voc of 44.146 and an FF of 0.893. The normalized Voc (48.674 and 44.597) and FF (0.901 and 0.894) of the BTIC-U5 and BTIC-U7 molecules are higher than those of the BTIC-R molecule due to their more excellent Voc. The results of the study show that these small molecules have amazing photoelectronic properties that may be useful in real-world applications.
The PCE, which sums up all of a solar cell's performance metrics, can determine whether a photovoltaic material is practical enough to use. A molecule's PCE is directly affected by the short circuit voltage (Jsc), the FF, and Voc. Whereas the strength of the radiation absorbed by the cell has the contrary effect. Equation (4) effectively describes this relationship.
By comparing each designed molecule to a standard reference molecule, we can see their Voc and FF values. It has been found that as Voc rises, so do FF and PCE. As can be seen in Table 9, BTIC-U5 and BTIC-U7 have the highest PCE and can be relied upon to improve the efficiency of solar devices.

Transition Density Matrix (TDM) and Binding Energy (Eb)
In conjugated chemical systems, TDM analysis seems to be needed to figure out where exciton (hole electron pair) transition will happen amongst the D-A regions. This analysis could look at exciton mobility during the absorption and emission process in the excited state, fundamental charge locations, and charge transition prediction. TDM plots are generally utilized to evaluate electronic properties, such as the level of delocalization, and the effects of resonance, and to gain insight into the charge mobility inside a molecule. Due to its negligible charge transfer properties, hydrogen's contribution to the transition is typically disregarded in TDM analysis. All atoms, with the exception of hydrogen, are represented by numbers on the left y-and bottom x-axes of the matrix, and the charge density coefficient is shown by the colored bars that progress from blue to red along the right y-axis. In order to determine the transition pathway, the molecule was split into three sections, which were labeled donor (D), acceptor (A), and acceptor1 (A1) in Figure 12. The charge density is found to be efficiently distributed in both the off-diagonal and diagonal directions throughout the molecules of all the compounds, with the diagonal direction responsible for the mainstream charge distribution. It was discovered that the flow of charge was effectively moving from donor to acceptor and acceptor1, which suggests that the whole molecule is enduring sequential conjugation. sections, which were labeled donor (D), acceptor (A), and acceptor1 (A1) in Figure 12. The charge density is found to be efficiently distributed in both the off-diagonal and diagonal directions throughout the molecules of all the compounds, with the diagonal direction responsible for the mainstream charge distribution. It was discovered that the flow of charge was effectively moving from donor to acceptor and acceptor1, which suggests that the whole molecule is enduring sequential conjugation.  The E b of the excitons is added aspect that should be considered in this context. By measuring the binding energy, the cell's operational efficiency, soar cell potential for holeelectron(exciton) dissociation, and its electrical properties can be estimated. The coulombic force between the hole and electron, denoted by E b , is the determining factor in their interactions. Less coulombic interaction exists between electrons and holes when E b is low and vice versa. The E b values in Table 10 were determined using the following Equation (5). E x in the above equation is the solvent phase's first excitation energy, and E gap is the HOMO-LUMO gap. In the chloroform solvent, the BTIC-U1, BTIC-U4, and BTIC-U5 contain higher E b than BTIC-R, whereas the BTIC-U2, BTIC-U3, BTIC-U6, and BTICU7 have lower E b in comparison to the BTIC-R molecules. It is important to point out that the E b of each molecule is comparable to the molecule that serves as the reference. Excitons present in molecules with a lower E b value are able to easily diffuse into free charge carriers, enabling them a good option for boosting the charge density of an electric current. The device's efficiency will increase as these (electron and hole) charge carriers move more quickly to their corresponding electrodes and produce more current. So, the BTIC-U6 molecule appears to have the lowest E b among all molecules, making it the most competent candidate.

Charge Transfer Analysis
The investigation of the interaction between BTIC-U1 and the well-known donor-type polymer PTB7-Th is being done so that we can evaluate the charge transfer (CT) capabilities of the molecules that we have designed. We have completed an in-depth analysis of the complex depicted in Figure 13.
lombic force between the hole and electron, denoted by Eb, is the determining factor in their interactions. Less coulombic interaction exists between electrons and holes when Eb is low and vice versa. The Eb values in Table 10 were determined using the following Equation (5).
Ex in the above equation is the solvent phase's first excitation energy, and Egap is the HOMO-LUMO gap. In the chloroform solvent, the BTIC-U1, BTIC-U4, and BTIC-U5 contain higher Eb than BTIC-R, whereas the BTIC-U2, BTIC-U3, BTIC-U6, and BTICU7 have lower Eb in comparison to the BTIC-R molecules. It is important to point out that the Eb of each molecule is comparable to the molecule that serves as the reference. Excitons present in molecules with a lower Eb value are able to easily diffuse into free charge carriers, enabling them a good option for boosting the charge density of an electric current. The device's efficiency will increase as these (electron and hole) charge carriers move more quickly to their corresponding electrodes and produce more current. So, the BTIC-U6 molecule appears to have the lowest Eb among all molecules, making it the most competent candidate.

Charge Transfer Analysis
The investigation of the interaction between BTIC-U1 and the well-known donortype polymer PTB7-Th is being done so that we can evaluate the charge transfer (CT) capabilities of the molecules that we have designed. We have completed an in-depth analysis of the complex depicted in Figure 13.  Significant dipole moment, subtle reorganization energy, low excitation energy, and sizeable bathochromic shift led us to choose BTIC-U1 over all other fabricated molecules. As a result, the BTIC-U1 molecule is selected as the best acceptor molecule. PTB7-Th, on the other hand, is a commonly used donor polymer in CT analysis. At the MPW1PW91/6-31G* functional, the composite made up of PTB7-Th and BTIC-U1 is optimized. Figure 14 depicts the BTIC-U1: PTB7-Th complex HOMO and LUMO pictorially. The figure clearly demonstrates that the donor polymer is aligned in a precisely parallel direction to the acceptor BTIC-U1. When undergoing CT or excitation, this orientation is optimal. The findings suggest that the density of charge in the HOMO state is located on the electron-rich donor moiety. 31G* functional, the composite made up of PTB7-Th and BTIC-U1 is optimized. Figure 14 depicts the BTIC-U1: PTB7-Th complex HOMO and LUMO pictorially. The figure clearly demonstrates that the donor polymer is aligned in a precisely parallel direction to the acceptor BTIC-U1. When undergoing CT or excitation, this orientation is optimal. The findings suggest that the density of charge in the HOMO state is located on the electron-rich donor moiety. In contrast, the density of electrons in the LUMO state is distributed across the acceptor molecule of BTIC-U1. This distribution proved that electrons were transferred from the donor molecule to the acceptor atom. This change, which shows the superposition of orbitals between molecules, is a crucial confirmation of charge transfer between different moieties. Therefore, this acceptor group has been shown to be helpful for the effective molecular design of organic solar devices.

Computational Methodology
The latest method for studying quantum mechanics in computing is called density functional theory, or DFT. Due to this, all molecules' geometric calculations were carried out using Gaussian 09 [52], and all three-dimensional molecular assemblies were created and viewed using GaussView 5.0 [53]. First, the reference structure (BTIC-R) [38] was optimized using the 6-31G* basis set and four exchange-correlation functionals: MPW1PW91 [54], CAM-B3LYP [55], B3LYP [56], ωB97XD [57], and PBEPBE [58]. The most promising of these functionals was then chosen after careful consideration, and subsequent analysis was performed on that one. TD-SCF simulations in solvent (chloroform) state was used to select the functional by comparing the resulting absorption spectra to the experimental spectra of BTIC-R [38]. Based on a comparison of the computed and experimental values of λmax, we concluded that the Modified Perdew Wang 1parameter (MPW1PW91) functional [54] is most coherent with the latter, thus providing a solid indication that this functional and the basis set will calculate modified molecules' photovoltaic properties. In contrast, the density of electrons in the LUMO state is distributed across the acceptor molecule of BTIC-U1. This distribution proved that electrons were transferred from the donor molecule to the acceptor atom. This change, which shows the superposition of orbitals between molecules, is a crucial confirmation of charge transfer between different moieties. Therefore, this acceptor group has been shown to be helpful for the effective molecular design of organic solar devices.

Computational Methodology
The latest method for studying quantum mechanics in computing is called density functional theory, or DFT. Due to this, all molecules' geometric calculations were carried out using Gaussian 09 [52], and all three-dimensional molecular assemblies were created and viewed using GaussView 5.0 [53]. First, the reference structure (BTIC-R) [38] was optimized using the 6-31G* basis set and four exchange-correlation functionals: MPW1PW91 [54], CAM-B3LYP [55], B3LYP [56], ωB97XD [57], and PBEPBE [58]. The most promising of these functionals was then chosen after careful consideration, and subsequent analysis was performed on that one. TD-SCF simulations in solvent (chloroform) state was used to select the functional by comparing the resulting absorption spectra to the experimental spectra of BTIC-R [38]. Based on a comparison of the computed and experimental values of λ max , we concluded that the Modified Perdew Wang 1parameter (MPW1PW91) functional [54] is most coherent with the latter, thus providing a solid indication that this functional and the basis set will calculate modified molecules' photovoltaic properties.
At this point, the absorption maximum results were processed using GaussSum [59], and the spectral representation was displayed using the Origin 2022 program. The density of states (DOS) and transition density matrix data of the BTIC-R and BTIC-U1 to BTIC-U7 molecules were analyzed, and the roles of various fragments of the molecules were determined with the help of Multiwfn 3.7 [60], PyMOlyze 2.0 [59] softwares.
The Marcus theory can be used to analyze the reorganization energy (RE), an essential factor for determining the number of intramolecular or intermolecular charge movements. Our investigations, however, are primarily concerned with intramolecular charge transfer (ICT). The overall RE comprises both the internal and external RE. Likewise, instances of external RE include polarization changes during charge transfer and abrupt changes in the outer atmosphere. Variability in the structures of molecules is involved in the internal RE. Environmental conditions affect external RE, making it difficult to calculate [61]. In this particular study, we concentrated solely on RE performed internally. We measured the anions (λ − ), and cations (λ + ) mobility, which together makes up the internal RE.
Optimized neutral molecular geometries yield cation and anion energies E + 0 and E − 0 . E + and E − are optimized cation and anion geometries. Neutral molecule energies E 0 − and E 0 + are calculated using optimal anion and cation structures. E 0 is the optimised neutral molecule's ground state point charge [62].