Synergistic end-capped engineering on non-fused thiophene ring-based acceptors to enhance the photovoltaic properties of organic solar cells

In this study, a series of non-fused thiophene ring-based small molecular acceptors (4T1–4T7) of A-D-A type are developed by the replacement of the end-groups of the 4TR molecule. The optoelectronic characteristics of the 4TR and 4T1–4T7 molecules are investigated employing the MPW1PW91 functional with the 6-31G (d,p) basis set, and solvent-state computations are studied using the TD-SCF. All the parameters estimated in this research are improved to a substantial level for the developed molecules as compared to the 4TR molecule, e.g. all the newly developed molecules have shown a red shift in their maximum absorption (λmax) and a reduced bandgap compared to the 4TR molecule, with ranges of 646 nm to 692 nm (in chlorobenzene solvent) and 2.34 eV to 2.47 eV, respectively. The reorganization energies of electron and hole mobility for almost all developed molecules are smaller than those for the 4TR molecule, with ranges of 0.00766–0.01034 eV and 0.01324–0.01447 eV, respectively. Hence, all the modified chromophores exhibit better charge capabilities than the 4TR molecule. The charge mobility of almost all the developed molecules is improved because of their reduced reorganization energies. The 4T2 molecule has minimum RE values for both electrons (0.00766) and holes (0.01324). The VOC values of all acceptor molecules are calculated with respect to the PTB7-Th donor. An elevation in VOC and FF values is exhibited by the 4T5 and 4T7 molecules. As a result, these end-capped engineered molecules should be proposed for the future manufacturing of highly efficient organic solar cells.


Introduction
Experts are attempting to make the transfer from fossil fuels to sustainable and innite energy sources in order to combat global warming and pollution. 1,2 Geothermal, solar, wind, tidal, and biofuels are some examples of renewable resources. Solar energy is the only abundant source, since it delivers 3.8 million EJ (EJ ¼ 1018 J) of energy each year which is entirely unrestricted and pollutant-free. 3,4 Solar cells utilize solar energy for the generation of electricity based on the photovoltaic effect; this is illustrated by the fact that when a photon of a particular energy is absorbed, free holes as well as electrons are generated in inorganic solar cells (SCs) or in the joint form in organic SCs as excitons, which create electricity as they go toward the opposite electrodes. [5][6][7] The rst generation of inorganic SCs were presented in 1954, but because of several drawbacks, such as short lifetimes, expensive synthetic approaches, and brittleness, secondgeneration solar cells have replaced them. 8,9 Cadmium telluride (CdTe) and amorphous silicon lms as well as copper indium gallium selenide (CIG), etc., 9 are thin-lm photovoltaic cells that have inexpensive production approaches, are less hazardous, and have higher visible spectrum absorption; however, they have low power conversion efficiency (PCE), i.e., 7-10%. 10 Meanwhile, organic solar cells (OSCs), which are thirdgeneration solar devices, are under the consideration of researchers because of their efficient properties, such as their transparent morphology, light weight, maximum absorption, and low cost. 11,12 Bulk heterojunctions are another modern approach toward photovoltaic cells in which a blended form of donor and acceptor material is present in the active layer, 13 such as the bulk heterojunction of a P 3 HT donor and PC 61 BM and PC 71 BM acceptors in fullerene-based OSCs. 14,15 In order to attain effective dissociation of excitons as well as transfer of charge, the donor and acceptor areas should possess their own consistent phase structures, correct phase sizes, maximum phase purities, and organized molecular packing and orientation. The optimization of morphology has long been an important topic in OSC research. Clearly, the active layer's shape is not only inuenced by the molecular aggregation qualities of discrete donors and acceptors but also by their intramolecular interactions. Today, the majority of highly efficient OSCs depend on techniques of device optimization such as additives, solvent annealing, thermal annealing and thermal spin-coating procedures. 16 Fullerene derivatives are considered good acceptor choices for organic photovoltaic (OPV) cells due to their electron-withdrawing efficiency, effective charge mobility and 3-D structures. 17 Fullerene-based acceptors have affected the OPV industry though some of their imperfections, such as nontunable energy levels, structural instability, expensive production methods, and absorption of shorter wavelengths. 18 Nonfullerene acceptors (NFAs) exclude the downsides of fullerenebased OSCs by their excellent properties, such as exible structure, tunable energy levels, cost-efficient production approaches and transparent morphology. 17,19,20 A variety of modications have been performed, such as alkyl chain engineering, interfacial alteration of donor-acceptor moieties as well as substitution of the terminal ends, to tune the optoelectronic properties and increase the efficiency of nonfullerene-based OSCs. [21][22][23] Fused ring electron-acceptors (FREAs) with an A-D-A molecular structure were described with PCEs equivalent or superior to those of the corresponding OSCs using fullerene acceptors since 2015. 24 OSCs have progressed to a new level in their modication, and a number of highly efficient FREAs, such as IDIC, 25 EICO-4F, 26 IHIC, 27 ITIC, 24 BT-IC, 28 and DTDP-IC, 29 have emerged. In spite of their remarkable efficiency, FREAs have complicated chemical structures that need lengthy syntheses and many purication procedures, resulting in poor yields and high production costs. 30 Low-cost acceptor materials should be made with a minimum of synthetic stages and simple purifying methods using readily accessible starting materials. 31,32 Using widely accessible derivatives of thiophene as precursors, non-fused electron acceptors (FNEAs) were produced based on tetrathiophene groups with very high yields in two stages. 30 In the present research, a non-fused tetrathiophene-based acceptor molecule with A-D-A symmetry is taken as a reference molecule (4TR). The reference molecule basically consists of four non-fused thiophene rings that are further linked to 2-(5,6-diuoro-2-methylene-3-oxo-indan-1-ylidene)-malononitrile accepting moieties at the end sites. 30 We developed seven new molecules (4T1-4T7) by the modication of the terminal acceptor moieties and replaced them with new acceptor moieties, i.e. 1-dicyanomethylene-2-methylene-3-oxo-indan-5,6dicarboxylic acid dimethyl ester (4T1), 1-dicyanomethylene-2methylene-3-oxo-indan-5,6-dicarbonitrile (4T2), 6-cyano-1dicyanomethylene-2-methylene-3-oxo-indan-5-carboxylic acid methyl ester (4T3), 2-(1-chloro-5-methylene-6-oxo-5,6-dihydrocyclopenta[c]thiophen-4-ylidene)-malononitrile (4T4), 2-(2methylene-3-oxo-2,3-dihydro-cyclopenta[b]naphthalen-1ylidene)-malononitrile (4T5), 2-(6,7-diuoro-2-methylene-3-oxo-2,3-dihydro-cyclopenta[b]naphthalen-1-ylidene)-malononitrile (4T6) and 2-(2-bromo-6-methylene-7-oxo-6,7-dihydro-1-thia-sindacen-5-ylidene)-malononitrile (4T7). ChemDraw geometries

Computational procedure
The Gaussian 09 (ref. 33) and Gaussview 6.0 (ref. 34) soware programs were used to perform quantum chemistry structure simulations based on density functional theory (DFT) 35 as well as geometric visualization of the reference 4TR and the designed molecules (4T1-4T7), respectively. To compute the molar absorptivity of reference 4TR, at restricted spin, 6-3IG (d, p) basis set, four hybrid functionals: B3LYP, 36 CAM-B3LYP, 37 MPW1PW91, 38 and wB97XD were used to optimize the ground state. 39 The molar absorptivity of the 4TR molecule was also estimated by employing the ZINDO/S method. 40 To model the absorption spectrum of 4TR, time-dependent self-consistent eld (TD-SCF) computations in the gas phase were performed aer geometry optimization. 41 The integral equation formalism polarizable continuum model (IEFPCM) was employed to study the inuence of a solvent (chlorobenzene). 42 Swizard soware 43 was used to calculate the maximum absorbance value (l max ) for 4TR and the developed molecules (4T1-4T7). The plots of the molar absorption coefficient values (l max ) were drawn using Origin 6.0 soware. 44 The l max achieved using the modied Perdew-Wang 1-parameter (MPW1PW91) functional was consistent with the results provided in the literature for the reference molecule 4TR, which clearly offers us a strong indication that it is suitable to calculate the modied molecules employing the nominated combination of functional and basis set.
The densities of states (DOS) of 4TR and the developed molecules (4T1-4T7) were explored using PyMOlyze-1.1 soware to determine the roles of different portions of the molecules, such as the donor and acceptor portions, in the absorption. 45 To show transition density matrix (TDM) results at the level of theory applied in electronic excited state computations, the Multiwfn-Multifunctional wavefunction analyzer was utilized. 46 The reorganization energies (RE) in semiconducting devices are closely related to the charge (electron and hole) mobilities. The entire RE is distributed into two classes: the outer area RE (l ext ) and the inner region RE (l int ). l int deals only with rapid changes in the internal geometry, 47 and l ext is only related to the external environment; hence, we can ignore it because of its minimal impact on the present research. The reorganization energies for the electrons (l e ) and holes (l h ) were estimated by eqn (1) and (2). 48 E 0 À and E 0 + represent the energies of anions and cations' neutral molecules. E + 0 and E À 0 are the ground-state energies of the anionic and cationic species, respectively. E À and E + are the anion and cation energies obtained by neutral molecule anion and cation geometrical optimization, while the neutral molecule with single point energy calculated in its ground energy state is E 0 .

Method selection and structural optimization
The geometry of the reference molecule 4TR was optimized by employing for four disparate functionals, namely, B3LYP, CAM-B3LYP, MPW1PW91 and wB97XD, with the basis set 6-31G (d,p) for the measurement of maximum absorption. The l max values of the 4TR molecule for these four functionals are 695, 497, 642, and 478 nm, respectively. The absorbance of the 4TR molecule was also estimated using the ZINDO/S method (one of the best methods for transition energy states, which is also a very effective and low-cost method 49 ), and the observed l max of this method for the reference molecule is 688 nm. The UV-visible spectra of the 4TR molecule with the abovementioned ve methods are shown in Fig. 2, and the estimated l max for all these methods are given in Table 1. According to the literature, the 4TR molecule exhibits signicant absorbance in the 500-700 nm wavelength spectrum, with l max around 600 nm; 30 this demonstrated the closest association with the MPW1PW91 functional and 6-31G (d,p) basis set, indicating a superior choice of DFT technique aer thorough consideration. As a result, this functional was selected to carry out more theoretical computations. Therefore, all the developed molecules (4T1-4T7) were optimized with the selected pair of functional and basis set, as shown in Fig. 3. Extended conjugation can be seen from the core to the terminal part of the molecules, indicating efficient charge transfer possibilities.
The bond lengths (L c-c ) as well as the dihedral angles (q ) of 4TR and 4T1-4T7 are listed in Table 2. It was found that the optimized structures of all the examined molecules (4TR, 4T1-4T7) possess comparable geometries, as they all have comparable bond lengths and dihedral angles. The C-C single bond length is 1.54Å and the C]C bond length is 1.33Å; the bond lengths of all the examined molecules ranged from 1.40 to 1.41 A, which clearly indicates enhanced charge transport characteristics and enables the delocalization of p-electrons via conjugation. The dihedral angles of all the examined molecules have a range of 0.254-0.594 , which means that all the molecules possess a planar geometry without any twist between the non-fused thiophene groups and terminal acceptors and also without any groups causing hindrance.

Frontier molecular orbitals (FMOs) analysis
Frontier molecular orbitals (FMOs, also known as the HOMO and LUMO) are important parameters of quantum chemistry that are studied to quantify transfer of charge, charge mobility, and molecular correlation. 50,51 The valence band HOMO provides electrons, whereas the conduction band LUMO receives them. The energy gap (E g ) shows the need for dissociation of excitons and is thus a distinguishing feature of SCs and other PV devices. 50,52 The electrical energy levels of an organic solar device are inversely related to the PCE, i.e., a smaller band gap results in greater efficiency. By the modication of the parent molecule 4TR, such as replacement of the terminal acceptor groups in this case, the HOMO-LUMO gap (E g ) can be improved. The charge distribution patterns in the FMOs of the reference molecule with all altered molecules at the MPW1PW91/6-31G (d,p) level are presented in Fig. 4.
In reference 4TR, the bonding orbitals have a calculated energy of À5.97 eV, while the antibonding orbitals have a calculated value of À3.49 eV with a bandgap of 2.47 eV; this is greater than the E g values of all the developed molecules (4T1-4T7), because the 4T1-4T7 molecules have efficient electronwithdrawing acceptor parts at the end sites of the molecules, and therefore these modied molecules possess shorter bandgaps than the 4TR molecule. In all the investigated molecules (4TR, 4T1-4T7), the charge density is mainly dominated by the non-fused thiophene core part of the molecules at the ground state, and the charge shis to the terminal acceptor groups of the molecules in the excited state. The band gaps of the reference and created compounds, as well as their HOMO and LUMO energies, are listed in Table 3. The HOMO and LUMO energy levels for all the examined molecules have , respectively. All the compounds developed exhibit a narrow bandgaps than the 4TR; only the 4T7 molecule has a comparable bandgap with the 4TR molecule. The declining order of E g for all the developed molecules is In general, molecules with robust electron-withdrawing moieties possess a smaller bandgap as well as greater light absorption, encouraging charge mobility. 53 The results illustrate that 4T2 possesses a smaller bandgap of 2.34 eV than all other examined molecules; the reason is the presence of a robust electronwithdrawing acceptor group at the terminal position, i.e., 1dicyanomethylene-2-methylene-3-oxo-indan-5,6-dicarbonitrile, which indicates the greatest mobility of charge from the HOMO to the LUMO.

Ionization potential and electron affinity
Ionization potential (IP) and electron affinity (EA) are suitable parameters to investigate charge transfer capabilities. 54 Charge transfer is improved when chromophores have a low IP and a high EA, because electron-donating groups destabilize the   56 The IP and EA values of all the considered molecules were calculated by using eqn (3) and (4), as suggested by Koopman's theorem, 57,58 and they are given in Table 3.

Absorption spectra
The electrical specications of chromophores can be estimated by employing quantum absorbance spectra. Chromophores are stimulated when they absorb photons of a specic energy relating to their respective band gaps. The absorbance spectrum (l max ), excitation energy (E x ), oscillator strength (f), assignment and dipole moment (D) of all the considered molecules were investigated by employing MPW1PW91/6-31G (d,p) in both the gas and solvent (chlorobenzene) forms and are given in Tables 4  and 5, respectively. The UV-visible spectra of molecules 4TR and 4T1-4T7 are represented in Fig. 5. Oscillator strength (f) is dimensionless, and it is critical in identifying the optical characteristics of photovoltaic cells and in calculating the radiation intensity emitted as a consequence of electrical excitation between two levels of energy. E x is the amount of energy that is mandatory for the possible transition; therefore, the increment of the oscillator strength (f), lower excitation energy, and wide absorption spectrum at the red shi are all predicted to result in efficient ICT. A considerable red shi was observed in the absorption of all the molecules compared to that of the 4TR molecule. The increasing absorption order of all the explored molecules in the gaseous form is 4TR < 4T7 < 4T1 ¼ 4T5 < 4T4 < 4T6 < 4T3 < 4T2, and the increasing sequence in the solvent (chlorobenzene) phase is 4TR < 4T7 < 4T5 < 4T4 < 4T6 ¼ 4T1 < 4T3 < 4T2. It is clear from the results that 4T2 exhibits the l max of 634 nm in the gas phase and 692 nm in chlorobenzene. In chlorobenzene solvent, all of the molecules displayed a bathochromic shi because the polar excited state was stabilized using a polar solvent. The LUMO produced a contraction in the E g and a rise in the wavelength because of the reverse connection between energy and wavelength. 4

Dipole moment
The dipole moment (D), which concerns the electron density and polarization of organic molecules, is a crucial aspect that adequately validates the photovoltaic activity of OSCs. The dipole moment determines the solubility of engineered molecules in polar organic solvents. 59 The dipole moment describes the efficient ow of electrons and holes between the donor and acceptor components. A large value of D indicates enhanced crystallinity and opens the possibility of electron-hole transit by increasing atom self-assembly. 2 The dipole moment values of molecules 4TR and 4T1-4T7 in both the gaseous and solvent phases are given in Tables 4 and 5, respectively. All the examined molecules possess the same increasing sequence of D in both the gaseous as well as the solvent phase, i.e., 4T2 < 4T4 < 4T3 < 4TR < 4T7 < 4T6 < 4T5 < 4T1. In the present case, the 4T2 molecule has the smallest dipole moment both in the gaseous and solvent forms, clearly indicating its lessened charge transfer and reduced solubility in solvent, while the 4T1 molecule possesses the greatest dipole moment in both phases because it has two different and robust electron-withdrawing groups (COOCH 3 and CN) at the edges of the acceptors, indicating efficient solubility in solvent and procient charge transfer.

Density of states (DOS)
DOS investigations are performed to further validate the roles of each fragment of a molecule, i.e., donor and acceptor, in the charge mobility of the molecule, as it is a helpful aspect for determining the arrangement of the electronic density distribution on the FMOs based on the Mulliken charge distribution theory. 60,61 DOS computations of all the examined molecules were executed by employing the MPW1PW91/6-31G (d, p) method, and the graphs were drawn using PyMOlyze 1.1 soware. These graphs have the energy (eV) on the x-axis and the relative intensity on the y-axis. For the analysis of the contributions of specic fragments of the molecules to the FMOs, each molecule was segmented into two parts, i.e., donor and acceptor. In the graphs of 4TR and 4T1-4T7, red, black and green lines represent the involvements of the donor and acceptor and the total contribution of the moieties to increasing the FMOs, respectively, as represented in Fig. 6, and the contribution data is listed in Table 6.
The results from the DOS analysis of all the molecules (4TR and 4T1-4T7) illustrated that all the molecules have comparable contributions of donor and acceptor, i.e., the donor fragment makes the major contribution in the ground state and the acceptor moieties are dominant in the excited state; this evidently signies the conjugation and facilitating the charge transfer from the donor to the acceptor moieties of the molecules, which will possibly improve the overall efficacy of organic SCs. These results also support the FMO analysis for reference and the designed molecules as stated above, as shown by Fig. 4. The results illustrate that the electron-rich donor core majorly  participates in the HOMO, and the charge shis from the core donor to the terminal acceptor part, which is the LUMO; moreover, the participation of the acceptor moieties is greater in almost all the newly developed molecules compared to the 4TR molecule, proving the efficient charge transfer and the greater efficiency of the molecules.

Reorganization energy
The reorganization energy (RE), which quanties the movement of charge from the donor to the acceptor fragments and is coupled to the hole and electron migration charge, is the driving factor in developing competent materials for OSCs. 62 The charge mobility is inversely related to the RE. The smaller the RE value, the more efficient the charge transfer that will take place. 63 The RE is affected by many factors, including the geometry of cations and anions. The electron mobility of acceptor and donor groups depends on their anionic and cationic geometries, respectively. The RE values of electrons and holes for all the examined molecules were calculated using eqn (1) and (2), and the values are listed in Table 7.
All the modied molecules (4T1-4T7) possess smaller electron mobility energies compared with 4TR, demonstrating that they have good electron mobilities and are hence excellent electron transporters. The decreasing sequence for l e is 4TR > 4T7 > 4T1 > 4T4 > 4T5 > 4T6 > 4T3 > 4T2. The hole (l h ) mobility value for the 4TR molecule is 0.01398 eV, and it is estimated that all the modied molecules except 4T1 have smaller values. The computationally calculated hole (l h ) mobility values have It was determined from the computationally analysed results that 4T2 has the lowest RE values for both electrons and holes; hence, it possesses the maximum electron and hole mobility. Overall, the results demonstrate that the modication of the 4TR molecule by new acceptor groups has a positive inuence on the charge transfer capabilities of the molecules.

Light harvesting efficiency (LHE)
The capability of a compound to generate charge aer light collection or to induce the conducting band is also a vital factor to consider for any component used in a solar cell. Eqn (5) was used to determine the LHE. 64 In the equation above, LHE is symbolised by h l and f is the oscillator strength. The oscillator strength of the solvent phase is used for the attribution of LHE. The computationally calculated values of LHE for 4TR and 4T1-4T7 are listed in Table 8. Among all the molecules, 4T7 has the maximum LHE value (0.9987). According to the following eqn (6), 65 LHE has a direct relationship with the short circuit current (J SC ): where LHE stands for the light harvesting efficiency, J SC for the short circuit current, h collect for the determined charge assembly, and f injected for the electron injection in the above equation. Because the above equation links LHE with J SC , the changes in LHE were shown to have a direct effect on J SC .

Molecular electrostatic potential (MEP)
MEP illustrates the 3-D charge dispersal taking place in a molecule and indicates the electron-rich and electron-poor sites of the molecule. 66 MEPs are colorful maps that indicate the existence of electron-withdrawing moieties, lone pairs and electrons on the molecular surface. The colored MEP maps of 4TR and 4T1-4T7 are shown in Fig. 7. Red color on the maps species the negative region where excessive electron density is present, green color shows the neutral sites, and blue color indicates the positive region, i.e., electron-decient sites. The maps of the investigated molecules show red color on the oxygen and nitrogen atoms and a small amount of red color on the uorine atoms present on the terminal acceptor groups, indicating strong electron density on these sites. Meanwhile, the core of non-fused thiophene rings with carbon atoms in the vicinity presents blue color on the MEP maps, representing an extreme deciency of electrons on those particular positions of the molecules.

Transition density matrix (TDM) and exciton binding energy
TDM gives qualitative information regarding the twodimensional movement of electrons, the dynamics of the transition as well as the electronic excitation of related hole-   67 A DFT model was used to compute the TDM energy with the MPW1PW91/6-31G (d,p) method, and then the schemes were plotted using g Multiwfn_3.7 (Multiwavefunction); the x-axis and the le y-axis represent the total numbers of atoms present in that particular molecule, but these numbers do not include hydrogen atoms because the system eliminates them by default because of their minute involvement in excitation. The y-axis (right side) shows the values of the electron density coefficient. In all the examined molecules (4TR and 4T1-4T7), the electron density is dispersed over all the molecules, mainly on their donor region, as shown in Fig. 8. Electronic consistency spreads diagonally and off-diagonally, but more dominantly diagonally, from donor to acceptor. The scheme illustrated effectual charge dispersal from the donor to the acceptor moiety of the molecule, indicating the efficient transfer of electrons from the donor to the acceptor by successive conjugation. The declining sequence of the interaction coefficients for 4TR and 4T1-4T7 is 4T4 > 4TR , and the values are given in Table 9. 4T2 has the lowest value of the interaction coefficient, demonstrating efficient charge transfer from the donor to the acceptor part of the molecule. Exciton binding energy (E b ) is a useful method for generating coulombic forces among charge carriers in conjugated molecules as a consequence of photoexcitation. There is less coulombic interaction between the electron-hole pairs of molecules with lower binding energies, and they move more quickly toward their corresponding electrodes. E b was calculated by using eqn (7). 68 Here, E g is the bandgap and E x is the rst excitation energy. The E b values of all the examined molecules were estimated in both gas and solvent phases and are listed in Table 9. Among all the investigated molecules, 4T4 has the smallest E b value in the gaseous and solvent phases; therefore, its exciton easily diffuses into independent charge carriers, making it a unique option for higher current charge density. Polar solvents interact and strongly bond with excitons; therefore, E b has a higher value in chlorobenzene than in the gaseous form.

Photovoltaic performance
Open circuit voltage (V OC ) is an inuential aspect to estimate the maximum efficacy of an OSC. The supreme current which can be achieved from an electrical equipment while the input voltage is zero is called V OC . 69 Certain factors that impact V OC include the light intensity, temperature of the solar cell, and charge mobility as well as the difference in energy between the acceptor's LUMO and the donor's HOMO. 70 In this study, the HOMO of PTB7-Th as a donor molecule was compared to the LUMOs of all the acceptor molecules (4T1-4T7), and V OC was estimated computationally. PTB7-Th is a competent donor, with À5.20 eV energy of its HOMO and À3.60 energy of its LUMO as reported in the literature. 71 The given eqn (8) was employed to calculate the V OC values of our investigated molecules by forming their complexes with the PTB7-Th donor.
In the preceding equation, e is the molecule charge, which is 1; 0.3 is the commonly used number for intersurface charge generation. Table 10 shows the computed V OC values of all the studied compounds. The HOMO and LUMO energy levels of the PTB7-Th donor and examined acceptor molecules with the calculated V OC values are shown in Fig. 9. The increasing sequence of V OC for the reference and modied molecules is 4T2 < 4T3 < 4T4 < 4TR ¼ 4T1 < 4T6 < 4T7 < 4T5. According to the given increasing order, the 4T5 molecule possesses the greatest V OC value; therefore, it can be utilized to boost the power PCEs of OSCs.

Fill factor (FF)
The ll factor (FF) has a direct link to the PCE of an SC as another critical component. It is a parameter that is mostly determined by the donor-acceptor V OC . Using eqn (9), FF was calculated: 72 in the above equation, and e is the standard charge, which is constantly 1. K B is the Boltzmann constant, and it has a value of 8.61733034 Â 10 À5 electron volts per kelvin. T is the temperature, which is constant (300 K). The computationally estimated normalized V OC and FF values for molecules 4TR and 4T1-4T7 are listed in Table 10. The range of normalized V OC is from 58.02 eV to 83.14 eV, and FF ranges from 0.9140 to 0.9354 with the same increasing sequence as V OC , i.e., The power conversion efficiency (PCE) is a method for gathering all the performance characteristics of a solar cell in a single value to verify that a photovoltaic substance is efficient enough for effective application. When analysing PCE, J SC , FF, and V OC are all critical factors to consider, as shown in the given eqn (10). 73  PCE is directly related to the aspects of FF, V OC , and J SC . P in is a constant value that represents the strength of incoming light hitting the cell interface. Theoretically, the FF and V OC values, and J SC , which is directly connected to LHE, were determined. The 4T5, 4T6 and 4T7 molecules possess better V OC , normalized V OC , FF and LHE values as compared to the reference molecule 4TR; therefore, it can be estimated that these molecules have improved PCEs as compared to the 4TR molecule.

Conclusion
In the present research, we designed a series of non-thiophene ring-based acceptor-donor-acceptor type molecules (4T1-4T7) by replacing the end-capped acceptor groups of the 4TR molecule. The optoelectronic characteristics of molecules 4TR and 4T1-4T7, such as the maximum absorption (l max ), frontier molecular orbitals (FMOs), density of states (DOS), dipole moment, transition density matrix (TDM), reorganization energy, LHE, and V OC , were investigated using the MPW1PW91 functional with the 6-31G (d,p) basis set, and the solvent state computations were studied using TD-SCF. The l max values of all the developed molecules (4T1-4T7) were shied toward higher wavelengths, ranging from 646 to 692 nm, as compared to the l max of the reference 4TR (642 nm). Except for 4T7, all the developed molecules have shorter bandgaps compared to the reference molecule, ranging from 2.47 to 2.34 eV; the 4T2 molecules possesses the lowest E g value of 2.34 eV. Almost all the developed molecules possess better charge mobility capacities because of their lower RE values as compared to the reference molecule, ranging from 0.01034 to 0.00766 eV for electron mobility and 0.01447 to 0.01324 eV for hole mobility; 4T2 has the minimum RE values for both electrons (0.00766) and holes (0.01324). Among all the examined molecules, 4T4 has the smallest E b value in the gaseous phase (0.34 eV) and the solvent phase (0.51 eV); therefore, its excitons easily diffuse into independent charge carriers. 4T2 has the lowest interaction coefficient (0.688), demonstrating efficient transfer of charge from the donor to the acceptor moieties of the molecule. The V OC values of all the studied acceptor compounds were estimated by forming their complexes with PTB7-Th as a donor; the values ranged from 1.50 to 2.12 eV. The normalized V OC and FF values have ranges of 58.02-83.14 eV and 0.9140-0.9354,