Study of the Contributions of Donor and Acceptor Photoexcitations to Open Circuit Voltage in Bulk Heterojunction Organic Solar Cells

One of the key parameters in determining the power conversion efficiency (PCE) of bulk heterojunction (BHJ) organic solar cells (OSCs) is the open circuit voltage (VOC). The processes of exciting the donor and acceptor materials individually in a BHJ OSC are investigated and are found to produce two different expressions for VOCs. Using the contributions of electron and hole quasi-Fermi levels and charge carrier concentrations, the two different VOC expressions are derived as functions of the energetics of the donor and acceptor materials and the photo-generated charge carrier concentrations, and calculated for a set of donor-acceptor blends. The simultaneous excitation of both the donor and acceptor materials is also considered and the corresponding VOC, which is different from the above two, is derived. The VOC calculated from the photoexcitation of the donor is found to be somewhat comparable with that obtained from the photoexcitation of the acceptor in most combinations of the donor and acceptor materials considered here. It is also found that the VOC calculated from the simultaneous excitations of donor and acceptor in BHJ OSCs is also comparable with the other two VOCs. All three VOCs thus derived produce similar results and agree reasonably well with the measured values. All three VOCs depend linearly on the concentration of the photoexcited charge carriers and hence incident light intensity, which agrees with experimental results. The outcomes of this study are expected to help in finding materials that may produce higher VOC and hence enhanced PCE in BHJ OSCs.


Introduction
Organic solar cells (OSCs) based on bulk heterojunction (BHJ) structure are promising candidates for the generation of clean, affordable electricity owing to their low fabrication cost compared to inorganic solar cells [1,2]. Recent research efforts into the development of new donor-acceptor blends have led to a significant improvement in the power conversion efficiency (PCE) of OSCs, above 11% [3][4][5][6], which is gradually moving OSCs closer to commercialization. However, poor material properties such as low charge carrier mobilities, low dielectric constants, and poor band offsets have led to the low PCE of OSCs relative to their silicon counterparts [1]. Hence, optimization of these parameters is essential for improving the performance of OSCs. Considerable attention has recently been focused on exploring strategies for the further optimization of the PCE, guided by a thorough understanding of the fundamental mechanisms that govern the photovoltaic performance of OSCs [7]. The following four processes distinguish BHJ OSCs from their inorganic counterparts: (i) photon absorption and exciton formation; (ii) exciton diffusion to the donor-acceptor (D-A) interface; (iii) exciton dissociation into free charge carriers at the D-A interface; and (iv) charge transport and Electronics 2017, 6, 75 2 of 14 collection at their respective electrodes [8,9]. These processes should be efficient enough to achieve an enhanced PCE of BHJ OSCs.
One of the important parameters of BHJ OSCs is the open-circuit voltage V OC , which is the voltage at which the net current flowing through the device is zero. It is established that in BHJ OSCs, V OC depends on the energetic difference between the lowest unoccupied molecular orbital (LUMO) of the acceptor material (E A LUMO ) and the highest occupied molecular orbital (HOMO) of the donor material (E D HOMO ), called the donor-acceptor effective energy gap (E DA ) [10][11][12][13][14]. Until now, the most explored approach to increasing V OC of BHJ OSCs has been to increase the energy of the LUMO of the fullerene acceptor material or decrease the energy of the HOMO of the donor material, or both [10][11][12][13][14][15]. Scharber and co-workers [10] empirically defined V OC as: where q is the elementary charge and ∆E loss is an empirical value denoting the energy losses occurring in transporting the charge carriers to the electrodes [16][17][18]. The validity of Equation (1) has been established for a number of donor-acceptor (fullerene) blends with a voltage loss ∆E loss ranging from 0.3−0.6 eV [19,20]. The question that may be asked is: what are the possible causes of the voltage loss ∆E loss and how can it be reduced? The origin of ∆E loss is still being debated, with several possible mechanisms such as bimolecular recombination, coulombic interactions, energetic disorder, etc., [12,[16][17][18]21,22]. An important loss mechanism in BHJ OSCs is the bimolecular recombination of free charge carriers with its rate R given by [23,24]: where γ = γ pre q ε o ε r (µ e + µ h ) is the bimolecular recombination coefficient (m 3 s −1 ): γ pre is the dimensionless reduction prefactor, ε 0 is the permittivity of free space, ε r is the dielectric constant, and µ e (µ h ) is the electron (hole) mobility (m 2 /Vs) in the active organic layer.
By minimizing the nonradiative bimolecular recombination in BHJ OSCs [25], one can minimize the voltage loss and eventually increase V OC . Recent studies [26,27] have shown that V OC can be determined more accurately by measuring the bimolecular recombination rate. These studies have also considered the influence of light intensity on V OC by varying the carrier concentration due to photoexcitation. Lange et al. [28] studied the dependence of V OC on the charge carrier concentration in polymer:fullerene blends using bias amplified charge extraction (BACE) measurements. In addition, as stated above, the energy loss also occurs due to the energetic disorder, which reduces the V OC . A few recent studies have shown the effect of energetic disorder on V OC using a Gaussian density of states (DOS) [29][30][31]. Garcia et al. [29] developed a model to predict the V OC values by introducing DOS distributions of the HOMO in donor and the LUMO in acceptor materials. In their study, the energetic disorder is modeled by the Gaussian DOS with the mean energy equal to E A LUMO and standard deviation σ e for the acceptor material and corresponding E D HOMO and σ h for the donor material. In other studies [16,25,[32][33][34][35][36], the role of energetics with respect to charge carrier dynamics and coulombic interactions have been analyzed. Credgington and Durrant [32] studied the dependence of V OC on the generation and recombination dynamics; they studied the dependence of V OC on the charge carrier concentration and carrier recombination lifetime. Vandewal et al. [25,[33][34][35] derived a relationship between V OC and the energy of the charge transfer (CT) exciton state E CT by incorporating the energetic disorders. However, other works [16,36] have derived different expressions relating V OC and E CT in OSCs which exclude the influence of energetic disorder.
Despite numerous attempts to study V OC and the associated energetic loss mechanisms in BHJ OSCs, there appears to be few attempts being made to understand whether it depends on the concentration of photoexcited charge carriers in the donor and acceptor materials. This may be due to the fact that as the two materials, donor and acceptor, are blended together in a BHJ structure, they cannot be separately excited experimentally. However, in order to understand the contributions of photoexcitation of donor and acceptor materials to the performance of a BHJ OSC, one may seek to assess the performance considering the two cases separately: (1) when absorption and photoexcitation occur only in the donor and (2) when these occur only in the acceptor. Such a study can only be carried out theoretically, as presented here. When an exciton is excited in the donor in a BHJ OSC, its electron can be transferred from the donor LUMO to the acceptor LUMO at the D-A interface because the acceptor LUMO lies at a lower energy (see Figure 1a). In this process, an exciton becomes a CT exciton as its electron is transferred to the acceptor and its hole remains in the donor, as shown in Figure 1c. Likewise, when an exciton is excited in the acceptor, its hole can be transferred from the acceptor HOMO to the donor HOMO, being at a lower hole energy (see Figure 1b), which also creates a CT exciton with an electron in the acceptor LUMO and a hole in the donor HOMO, as shown in Figure 1c.
OC CT by incorporating the energetic disorders. However, other works [16,36] have derived different expressions relating OC V and CT E in OSCs which exclude the influence of energetic disorder.
Despite numerous attempts to study OC V and the associated energetic loss mechanisms in BHJ OSCs, there appears to be few attempts being made to understand whether it depends on the concentration of photoexcited charge carriers in the donor and acceptor materials. This may be due to the fact that as the two materials, donor and acceptor, are blended together in a BHJ structure, they cannot be separately excited experimentally. However, in order to understand the contributions of photoexcitation of donor and acceptor materials to the performance of a BHJ OSC, one may seek to assess the performance considering the two cases separately: (1) when absorption and photoexcitation occur only in the donor and (2) when these occur only in the acceptor. Such a study can only be carried out theoretically, as presented here. When an exciton is excited in the donor in a BHJ OSC, its electron can be transferred from the donor LUMO to the acceptor LUMO at the D-A interface because the acceptor LUMO lies at a lower energy (see Figure 1a). In this process, an exciton becomes a CT exciton as its electron is transferred to the acceptor and its hole remains in the donor, as shown in Figure 1c. Likewise, when an exciton is excited in the acceptor, its hole can be transferred from the acceptor HOMO to the donor HOMO, being at a lower hole energy (see Figure 1b), which also creates a CT exciton with an electron in the acceptor LUMO and a hole in the donor HOMO, as shown in Figure 1c. This has prompted several researchers to investigate the above two processes in more detail [37][38][39][40][41][42][43][44][45][46]. Recently, two different external quantum efficiencies (EQE) due to different donor and acceptor absorption bands have been reported [37]. The difference in the absorption bands can be attributed to the difference in the energy levels involved. The times of transfer of electrons from the donor LUMO to the acceptor LUMO and that of holes from the HOMO of the acceptor to the HOMO of the donor have been measured [38,39] and have been found to be different; the hole transfer was found to be faster than the electron transfer. This may be expected to lead to generating different concentrations This has prompted several researchers to investigate the above two processes in more detail [37][38][39][40][41][42][43][44][45][46]. Recently, two different external quantum efficiencies (EQE) due to different donor and acceptor absorption bands have been reported [37]. The difference in the absorption bands can be attributed to the difference in the energy levels involved. The times of transfer of electrons from the donor LUMO to the acceptor LUMO and that of holes from the HOMO of the acceptor to the HOMO of the donor have been measured [38,39] and have been found to be different; the hole transfer was found to be faster than the electron transfer. This may be expected to lead to generating different concentrations n and p of electrons and holes. Different charge generations by exciting the donor and acceptor individually have also been reported in some donor-acceptor combinations [40], and different internal quantum efficiencies (IQE) have also been observed [41]. Some studies have demonstrated that BHJ OSCs based on thienothiophene-substituted diketopyrrolopyrrole (DPP) polymers (PDPP2TT-T) possess an EQE of 0.63-0.78 within the wavelength range of acceptor-fullerene absorption (400-650 nm), whereas a lower EQE of 0.35-0.50 is achieved within the range of donor-polymer absorption (650-900 nm) [41,42]. Such a difference may not be attributed to a difference in absorption of photons in the different materials, but may be due to the difference in charge generation. In another study, Baulin et al. [43] compared the dynamics of CT states generated through electron-transfer and hole-transfer exciton-dissociation pathways in three BHJ OSCs, and they found that the dynamics of CT state recombination are very similar for the two charge-generation processes, implying that the nature of the generated CT states is independent of the pathways of CT formation. They did, however, observe that the generation of CT states and free charge carriers greatly depends on the actual LUMO energy offset ∆E LUMO = E D LUMO − E A LUMO (electron transfer) and the HOMO energy offset ∆E HOMO = E D HOMO − E A HOMO (hole transfer) of the blended system [44]. These energy offsets are independent of the pathway of the formation of CT excitons. Singh et al. [45,46] recently studied the contribution of exciton generation in the donor and acceptor materials and its dissociation to the performance of BHJ OSCs. They found that the absorption and dissociation rates as well as exciton diffusion lengths are comparable whether excitons are generated in the donor or acceptor.
In this paper, we investigate in detail the difference in the open circuit voltage from the photoexcited donor and photoexcited acceptor separately in the following seven different polymer:fullerene blends: [3, 4-b]thiophenediyl] (PTB7), mono-DPP, and bis-DPP blended with different derivatives of methano-fullerene[6,6]-phenyl C 61 ,-butyric acid methyl ester (PCBM). We derived two different expressions for V OC in BHJ OSCs. The first expression, which is a result of photoexcitation of the donor, depends primarily on the energetic difference between the hole quasi-Fermi level (E h F ) and the HOMO of the donor material (E h F − E D HOMO ) and the electron concentration n. The second expression, which is a result of photoexcitation of the acceptor, depends primarily on the energetic difference between the LUMO and the electron quasi-Fermi level (E e F ) of the acceptor material (E A LUMO − E e F ) and the hole concentration p. Thus, the two open circuit voltages calculated in each selected D-A blend are presented here and their influences on the performance of BHJ OSCs are discussed. Following this, the simultaneous photoexcitation of both the donor and acceptor materials is also considered and the corresponding V OC , which is different from the above two, is derived, and the calculated V OC s are compared with their corresponding measured values and a reasonably good agreement is found between them.

Theoretical Formalism
Under no illumination condition of equilibrium, assuming the condition of the open circuit voltage in a BHJ OSC, the electron and hole current densities are assumed to be independently zero, and under this condition the quasi-Fermi levels in the donor and acceptor materials are iso-energetic at energy E F throughout the device. Once such a BHJ OSC is illuminated, charge carriers are continuously generated in the active layer. As the concentration of charge carriers increases, the quasi-Fermi level E F splits into the hole and electron quasi-Fermi levels located, respectively, at energies E h F in the donor material and at E e F in the acceptor, as shown in Figure 2. Under the condition of illumination, V OC in any BHJ OSC is given by [47]: Under the condition of illumination, OC V in any BHJ OSC is given by [47]:

Photoexcitation of the Donor
We consider the photoexcitation of only the donor material here, where an excited exciton at the D-A interface forms a CT exciton by transferring its electron to the acceptor's LUMO, being at a lower energy [48]. The electron concentration ) (n in the LUMO of the fullerene-acceptor material and the hole concentration ) ( p in the HOMO of the polymer-donor material are, respectively, given by the Boltzmann distribution as [49]: assume that p n ≠ [39,40], and then divide Equation (6) by 2 n to obtain: The Fermi energies E e F and E h F may be expected to be different for excitations of the donor and acceptor individually. Although the difference may not be expected to be significant, this may lead to two different values of V OC according to Equation (3).

Photoexcitation of the Donor
We consider the photoexcitation of only the donor material here, where an excited exciton at the D-A interface forms a CT exciton by transferring its electron to the acceptor's LUMO, being at a lower energy [48]. The electron concentration (n) in the LUMO of the fullerene-acceptor material and the hole concentration (p) in the HOMO of the polymer-donor material are, respectively, given by the Boltzmann distribution as [49]: where is the effective density of states for the LUMO (HOMO) of the acceptor (donor) material, T is the temperature, and k B is the Boltzmann constant. Under the conditions of V OC and bimolecular recombination, the generation rate of photoexcited electron and hole pairs G is equal to the bimolecular recombination rate given in Equation (2), which gives: Assuming a non-optimized BHJ OSC with unbalanced charge carrier mobilities (µ e = µ h ), we assume that n = p [39,40], and then divide Equation (6) by n 2 to obtain: Electronics 2017, 6, 75 6 of 14 Dividing Equation (5) by Equation (4) and using it in Equation (7) gives: Substituting E A LUMO = E D HOMO + E DA into Equation (8) and using Equation (3), we get: Assuming N C = N V [38,41], and rearranging Equation (9), we obtain V OC related to the excitation of donor only as: where:

Photoexcitation of the Acceptor
Next, we consider the photoexcitation of only the acceptor material, where the excited exciton at the D-A interface forms a CT exciton by transferring its hole to the donor's HOMO, being at a lower energy [48]. This gives a different open circuit voltage, which is derived as follows.
Here, we divide both sides of Equation (6) by p 2 to obtain: Dividing Equation (4) by Equation (5) and using it in Equation (12) gives: Substituting E D HOMO = E A LUMO − E DA into Equation (13) and using Equation (3), we get: Assuming N C = N V [38,41] and rearranging Equation (14), we obtain V OC for this case as: where: We have thus derived two different expressions for the V OC ; the first one in Equation (10) is obtained when the donor material is excited and the second in Equation (15) is obtained when the acceptor is excited.

Simultaneous Photoexcitation of the Donor and Acceptor
The theoretical developments presented in Sections 2.1 and 2.2, where only the donor or acceptor is excited, are mostly applicable in layered structures of donor and acceptor organic solar cells. However, in a BHJ OSC where the donor and acceptor are blended together in the active layer, both the donor and acceptor may be excited simultaneously. In this case, we obtain a different V OC as follows. Multiplying Equations (4) and (5), we obtain: Substituting −E DA = E D HOMO − E A LUMO and using Equation (3) in Equation (17) we thus get: where: It may be noted that, unlike V OC s obtained in Equations (10) and (15), which depend on the electron (n ) and hole (p) concentrations, respectively, Equation (18) depends on both n and p. The dependence of V OC in Equation (18) on √ np is presented mainly to compare it with the results derived in Equations (10) and (15).

Results
For calculating the V OC using Equations (10), (15), and (18), we require some material-dependent input parameters. For this, as stated above, we chose seven donor-acceptor blends with known input parameters available in the literature, as listed in Table 1. Table 1. List of input parameters required for calculating V OC from Equations (10), (15), and (18): LUMO − E D HOMO )/q, generation rate of bound polaron pairs (G) and bimolecular recombination coefficient (γ), electron (hole) mobility µ e (µ h ), and dielectric constant ε r of each donor-acceptor (D-A) blend.

Active Layer Blend
, as well as the carrier concentrations n and p in each blend to calculate V OC from Equations (10), (15), and (18), which are also listed in Tables 2 and 3. However, it may be noted that both n and p are only available for PTB7:PC 71 BM [62] as given in Tables 2 and 3. In mono-DPP:PC 71 BM and bis-DPP:PC 71 BM, only the electron concentration n is known in the literature [63]. Consequently, to find the hole concentration p in these two blends, we used Equation (7) and calculated p; these n and p thus obtained are listed in Tables 2 and 3. We first used the n and p values determined above for PTB7:PC 71   OC .  OC .   Tables 2 and 3). That is to say, using Equations (4) and (5) with N C = N V = 2.5 × 10 25 m −3 known for MDMO-PPV:PC 61 BM [17,57], N C = N V = 2 × 10 26 m −3 for P3HT:PC 60 BM [1], N C = N V = 5 × 10 25 m −3 for PCDTBT:PC 71 BM [52], and N C = N V = 2.5 × 10 25 m −3 for OC 1 C I0 -PPV:PC 61 BM [1,59], in addition to the corresponding average charge carrier concentrations n = 1.8 × 10 22 m −3 and p = 0.95 × 10 22 m −3 , we calculated (E A LUMO − E e F ) and (E h F − E D HOMO ), respectively, in the other four blends. Using these parameters in Equations (10) and (15), the V OC is calculated in all seven blends for the photoexcitation of the donor and acceptor as given in Tables 2 and 3, respectively. The calculated values of V OC from Equation (10) when the donor is excited in most selected materials are somewhat comparable with those calculated from Equation (15) when the acceptor is excited, as also shown in Figure 3. In addition, considering the case of simultaneous photoexcitation of donor and acceptor materials, the open circuit voltage in each blend is calculated from Equation (18) as given in Table 4. All the three V OC s obtained from the asynchronous and simultaneous photoexcitations of donor and acceptor are found to be comparable and agree well with the corresponding measured values [51,54,59,61,64,65].

Active Layer Blend
The open circuit voltages V OC s calculated from Equations (10) and (15) are plotted as a function of the effective band gap E DA /q as shown in Figure 3, which illustrates that V OC calculated from Equation (10) is slightly higher than that from Equation (15). [51,54,59,61,64,65].
The open circuit voltages s OC V calculated from Equations (10) and (15) are plotted as a function of the effective band gap q E DA / as shown in Figure 3, which illustrates that OC V calculated from Equation (10) is slightly higher than that from Equation (15).  Table 4 for comparison.    Tables 2 and 3, respectively, voltage offset ( E DA/q − V DA OC ) upon simultaneous excitation of both the donor and acceptor materials, and the corresponding measured open circuit voltage V meas.
OC . Finally, as stated previously, in a BHJ OSC the donor and acceptor materials are blended together. Consequently, both the donor and acceptor are excited simultaneously, resulting in one open circuit voltage as derived in Equation (18), which is used to calculate the V DA OC for simultaneous photoexcitations of both the donor and acceptor in each of the selected blends, as listed in Table 4.

Active Layer Blend
The measured values of open circuit voltage V meas.
OC are also given in Table 4 for comparison.

Discussions
A comprehensive study of the open circuit voltage in BHJ OSCs has been carried out. Three different expressions for V OC are derived; the first in Equation (10) is obtained when the donor material in a BHJ OSC is excited, the second in Equation (15) is obtained when the acceptor material is excited, and the third in Equation (18) is obtained when both the donor and acceptor materials are excited simultaneously. As stated previously, when the photoexcitation occurs in the donor, an electron is transferred from the LUMO of the donor to the LUMO of the acceptor, and conversely, when photoexcitation occurs in the acceptor, a hole is transferred from the HOMO of the acceptor to the HOMO of the donor [48]. The photoexcitation of the donor and acceptor individually produces comparable V OC values from Equations (10) and (15), as shown in Figure 3. Moreover, according to Figure 3, the linear correlation between V OC and E DA /q agrees well with previous studies [20,30,32,33,64]. The third open circuit voltage V DA OC calculated in Equation (18) from the simultaneous excitations of both the donor and acceptor, as given in Table 4, is also comparable with the other two V OC s. All three V OC s thus derived produce comparable results and agree reasonably well with the experimental values (see Tables 2-4). However, the calculated V OC values in Tables 2-4 are obtained from slightly different energetic parameters; as a result, we obtained three different sets of V OC s for the photoexcitation of the donor, photoexcitation of the acceptor, and simultaneous photoexcitation of both. However, experimentally, there is only one set of measured values of V OC to compare the theoretical results with. As a result, the calculated range of V OC for the materials considered here is 0.50-0.96 V, which is slightly different from the range of experimental values (0.52-0.85 V). Thus, as stated above, this discrepancy may be attributed to the different energies (E h F − E D HOMO )/q used in the calculation of V OC from Equations (10), (15), and (18). In addition, the two V OC s in Equations (10) and (15) depend linearly on ln n and ln p, respectively, and this agrees well with previously measured results [28,63]. This also holds for V OC derived from Equation (18), which depends on ln( np). By setting n = p in Equations (10), (15), and (18), we found that all three V OC s increase linearly with ln n and have the same slope equal to k B T q / log(exp(1)) ≈ 2.3 × 2k B T/q ≈ 119 mV. This implies that the materials with higher photo-generated charge carrier concentration will have higher V OC in comparison with materials with lower carrier concentration.
It is also worth noting that the dependence of V OC on the incident light intensity has been studied [17,24]. Substituting np = G γ ∝ I abs = I o [1 − exp(−αt)] into Equation (18), where I abs is the absorbed light intensity, I 0 is the incident light intensity, α is the absorption coefficient, and t is the thickness of the active layer, we get: The V OC in Equation (20) depends linearly on ln I 0 with a slope of k B T q , which agrees perfectly well with the plotted V OC as a function of ln I 0 in [17,24].
Orlowski et al. [66] used the energetics of quasi-Fermi levels to study V OC in a heterojunction solar cell. They found that for a particular donor (p-type)-acceptor (n-type) heterojunction, the electron concentration (n) in the conduction band (LUMO) of the acceptor and the hole concentration (p) in the valence band (HOMO) of the donor, coupled with their respective quasi-Fermi levels, can create open circuit voltages like two independent cells. This agrees very well with our concept of deriving two different V OC s in Equations (10) and (15), although both the expressions produce similar results.
Furthermore, we find that the V D OC lies 0.34-0.49 V below the E DA /q in column 5 of Table 2 when the donor is excited, comparable to 0.34-0.48 V below the E DA /q in column 5 of Table 3 when the acceptor is excited. This implies that both the donor and acceptor photoexcitations can contribute to the achievable open circuit voltage as well as the photovoltaic performance of BHJ OSCs; this agrees with our earlier results [45]. The voltage offsets ( E DA/q − V j OC ) (where j = D, A or DA) of the various blends listed in Tables 2-4, ranging from 0.34-0.50 V, agree well with the experimental values measured at room temperature in the range of 0.25-0.48 V for small molecule OSCs [67], 0.32 V for polymer:fullerene solar cells of indacenoedithiophene (IDT) polymer [68], 0.34-0.44 V for a range of donor-acceptor blends [69], 0.38 V for OC 1 C 10 -PPV:PCBM solar cells [12], and 0.30-0.60 V for BHJ OSCs [20]. Other low voltage offsets ( E DA/q − V j OC ) that have been reported recently are 0.25 V for an evaporated bilayer OSC [70], 0.23 V and 0.26 V for diketopyrrolopyrrole-thieno[2,3-f]benzofuran (DTD):PC 60 BM and DTD:naphthalene diimide acceptor-polymer (N2200) systems, respectively [71], and 0.34-0.40 V for BHJ OSCs [72], with which our calculated values also somewhat agree in the range.
Overall, it may be deduced from our results that separate and simultaneous photoexcitations of the donor and acceptor produce comparable V OC s, and the highest loss in V OC is found to be contributed by the energetics of the donor and acceptor materials. Therefore, to achieve a high V OC , it is necessary to choose higher efficiency polymer-donor and acceptor materials, as well as control the donor and acceptor material energetics, which may lead to improvements in the PCEs of BHJ OSCs [32].

Conclusions
In summary, we have derived two different expressions for the calculation of the V OC in BHJ OSCs; one is used when the exciton is excited in the donor material and the other is employed when the exciton is excited in the acceptor material. We have also derived an expression for V OC when both the donor and acceptor materials are excited simultaneously, and its calculated values agree better with their measured values. All three V OC s obtained from separate and simultaneous excitations are found to be comparable, implying that the contributions of both donor and acceptor photoexcitations are comparable in BHJ OSCs. Also, it is found that materials that generate higher photo-generated charge carrier concentrations and have lower energetic distances of (E h F − E D HOMO ) or (E A LUMO − E e F ) may lead to higher V OC . Invention and application of such materials in the fabrication of BHJ OSCs may lead to higher V OC s and hence enhanced PCEs.
Author Contributions: Both the authors have contributed to the publication of this paper. Douglas Yeboah carried out all the calculations. With the guidance from Jai Singh, the paper was initially drafted by Yeboah and then was finalized by them together.

Conflicts of Interest:
The authors declare no conflict of interest.