On the Potential of Silicon Intermediate Band Solar Cells

: Intermediate band solar cells (IBSCs) have an e ﬃ ciency limit of 63.2%, which is signiﬁcantly higher than the 40.7% limit for conventional single gap solar cells. In order to achieve the maximum e ﬃ ciency, the total bandgap of the cell should be in the range of ~2 eV. However, that fact does not prevent other cells based on di ﬀ erent semiconductor bandgaps from beneﬁting from the presence of an intermediate band (IB) within their bandgap. Since silicon (1.12 eV bandgap) is the dominant material in solar cell technology, it is of interest to determine the limit e ﬃ ciency of a silicon IBSC, because even a modest gain in e ﬃ ciency could trigger a large commercial interest if the IB is implemented at low cost. In this work we study the limit e ﬃ ciency of silicon-based IBSCs considering operating conditions that include the use of non-ideal photon casting between the optical transitions, di ﬀ erent light intensities and Auger recombination. The results lead to the conclusion that a silicon IBSC, operating under the conventional model in which the sub-bandgaps add to the total silicon gap, provides an e ﬃ ciency gain if operated in the medium-high concentration range. The performance of these devices is a ﬀ ected by Auger recombination only under extremely high concentrations.


Introduction
The idea behind the intermediate band solar cell (IBSC) [1,2] concept is the absorption of sub-bandgap energy photons to produce electric work. These photons are absorbed by engineering a semiconductor-like material that, in addition to the conduction and valence bands, exhibits an intermediate band (IB) within the otherwise conventional semiconductor gap.
The basic structure of an IBSC is illustrated in Figure 1, where the (IB) is shown splitting the total bandgap of the semiconductor, E G , into two sub-bandgaps, E L (the low energy sub-bandgap) and E H (the high energy sub-bandgap). The absorption of sub-bandgap photons is illustrated in Figure 1 by the arrows labeled a1 and a2. This absorption produces an increase in cell photocurrent, providing a gain in efficiency if the output voltage is not limited by any of the sub-bandgaps E L nor E H . The IB is isolated from the metal contacts by means of conventional semiconductors (called emitters) in order to assure that the output voltage is limited by E G .
(which defines the optical bandgap ) and transitions a1 take place between lower states in the IB to the conduction band (CB) (which defines the optical bandgap ) so that + . An approach to implement IB materials consists of introducing a high concentration of impurities inside the lattice of a host semiconductor. The impurity-host interaction can result in the emergence of deep-levels (DLs) inside the semiconductor bandgap, which usually act as non-radiative recombination centers. However, if the impurity concentration is high enough to produce Mott transition [4], the wave function of the DLs can be delocalized, forming an IB [5].
There is some controversy about whether sufficiently increasing the concentration of impurities can reduce the non-radiative recombination related to DLs or not, with theoretical works that support this possibility [6,7] and theoretical works that do not [8]. To test experimentally the limitations of the DLs approach well-known semiconductors such as silicon have been studied, despite its low bandgap (~1.12 eV) when compared to the optimum one (1.95 eV) [1] for IBSC operation. Several works suggest that an IBSC based on silicon could be implemented in practice by inserting impurities of transition metals such as titanium (Ti) [9][10][11], cobalt (Co) [12], selenium (Se) [13,14] and sulfur (S) [14][15][16][17]. This is supported by experimental results that prove some of the features related to the IBSC operation, such as sub-bandgap absorption [14], lifetime recovery [16][17][18] and splitting of quasi Fermi levels between the IB and the CB [19].
In spite of this controversy, since silicon is the dominant material in solar cell technology, it is important to clarify whether silicon IBSC can provide an efficiency gain, since even a modest increase in efficiency could trigger a large commercial interest. For this reason, the objective of this work is to determine the efficiency of silicon-based non-ratchet IBSCs. To this end, we will consider not only the ideal radiative limit case, but also a more realistic scenario for the case of silicon that includes Auger non-radiative recombination, unavoidable even in defect-free semiconductors. Non-ideal photon casting between the optical transitions will also be considered. The case in which the optical gaps are not constrained to + = , operating in the radiative limit, can be found in the work by Krishna and J. J. Krich [20].

Materials and Methods
The model to determine the radiative limiting efficiency of an IBSC was first formulated in reference [1]. This model assumed "non-overlapping" of the absorption coefficients, which means that the absorption of photons with a specific energy can only take place in one electronic transition. In particular, to achieve the limit efficiency of an IBSC the absorption coefficients have to be tailored in a way that a photon with energy E ph , high enough to be absorbed by means of several inter-band transitions, is absorbed in the electronic transition involving the highest bandgap. By imposing this condition, the thermalization losses (E loss ) resulting from absorbing photons with excessive energy (with an energy higher than the involved inter-band transition) are minimized. Therefore, when Concerning the properties of the IB, we will keep our work within the conventional framework that considers a non-ratchet IB. In the non-ratchet framework, the optical gaps are constrained to E L + E H = E G . The "photon-ratchet" intermediate band solar cell proposed by Yoshida et al. in 2012 [3] improves the conventional model by distinguishing between different energy states within the IB so that, for example, transitions a2 take place from the valence band (VB) to upper states in the IB (which defines the optical bandgap E H ) and transitions a1 take place between lower states in the IB to the conduction band (CB) (which defines the optical bandgap E L ) so that E L + E H E G .
An approach to implement IB materials consists of introducing a high concentration of impurities inside the lattice of a host semiconductor. The impurity-host interaction can result in the emergence of deep-levels (DLs) inside the semiconductor bandgap, which usually act as non-radiative recombination centers. However, if the impurity concentration is high enough to produce Mott transition [4], the wave function of the DLs can be delocalized, forming an IB [5]. There is some controversy about whether sufficiently increasing the concentration of impurities can reduce the non-radiative recombination related to DLs or not, with theoretical works that support this possibility [6,7] and theoretical works that do not [8]. To test experimentally the limitations of the DLs approach well-known semiconductors such as silicon have been studied, despite its low bandgap (~1.12 eV) when compared to the optimum one (1.95 eV) [1] for IBSC operation. Several works suggest that an IBSC based on silicon could be implemented in practice by inserting impurities of transition metals such as titanium (Ti) [9][10][11], cobalt (Co) [12], selenium (Se) [13,14] and sulfur (S) [14][15][16][17]. This is supported by experimental results that prove some of the features related to the IBSC operation, such as sub-bandgap absorption [14], lifetime recovery [16][17][18] and splitting of quasi Fermi levels between the IB and the CB [19].
In spite of this controversy, since silicon is the dominant material in solar cell technology, it is important to clarify whether silicon IBSC can provide an efficiency gain, since even a modest increase in efficiency could trigger a large commercial interest. For this reason, the objective of this work is to determine the efficiency of silicon-based non-ratchet IBSCs. To this end, we will consider not only the ideal radiative limit case, but also a more realistic scenario for the case of silicon that includes Auger non-radiative recombination, unavoidable even in defect-free semiconductors. Non-ideal photon casting between the optical transitions will also be considered. The case in which the optical gaps are not constrained to E L + E H = E G , operating in the radiative limit, can be found in the work by Krishna and J. J. Krich [20].

Materials and Methods
The model to determine the radiative limiting efficiency of an IBSC was first formulated in reference [1]. This model assumed "non-overlapping" of the absorption coefficients, which means Energies 2020, 13, 3044 3 of 10 that the absorption of photons with a specific energy can only take place in one electronic transition. In particular, to achieve the limit efficiency of an IBSC the absorption coefficients have to be tailored in a way that a photon with energy E ph , high enough to be absorbed by means of several inter-band transitions, is absorbed in the electronic transition involving the highest bandgap. By imposing this condition, the thermalization losses (E loss ) resulting from absorbing photons with excessive energy (with an energy higher than the involved inter-band transition) are minimized. Therefore, when investigating the maximum efficiency of IBSCs, the condition of non-overlapping absorption coefficients is imposed and the values of all the bandgaps are optimized in order to seek this maximum. However, when the value of all the bandgaps is fixed and determined by the fabrication process, as it would be in the case of IBSCs based on silicon, an increase in the thermalization losses resulting from overlapping absorption coefficients (illustrated in Figure 2d) can be compensated by a better current matching in the IBSC [20]. investigating the maximum efficiency of IBSCs, the condition of non-overlapping absorption coefficients is imposed and the values of all the bandgaps are optimized in order to seek this maximum. However, when the value of all the bandgaps is fixed and determined by the fabrication process, as it would be in the case of IBSCs based on silicon, an increase in the thermalization losses resulting from overlapping absorption coefficients (illustrated in Figure 2d) can be compensated by a better current matching in the IBSC [20]. The mathematical model for IBSCs with overlapping absorption coefficients was described in reference [21]. We use this model, based on the detailed balance, to study the limiting efficiency of silicon IBSCs, assuming, for simplification purposes, the sun is a black body at 6000 K (TS) and the cell and its surroundings are black bodies at 300 K (TC). In this context, the general equation describing the current-voltage ( ) characteristic is given by: where is the photon energy, and factors and , which can be considered equivalent to electron and hole absorptivity, are expressed as The mathematical model for IBSCs with overlapping absorption coefficients was described in reference [21]. We use this model, based on the detailed balance, to study the limiting efficiency of silicon IBSCs, assuming, for simplification purposes, the sun is a black body at 6000 K (T S ) and the cell and its surroundings are black bodies at 300 K (T C ). In this context, the general equation describing the current-voltage (J − V) characteristic is given by: where E is the photon energy, and factors M e and M h , which can be considered equivalent to electron and hole absorptivity, are expressed as α XY is the photon absorption coefficient related to the transition from the X band to the Y band; e is the electron charge; F abs is the photon flux absorbed by the cell (which equals the incoming photon flux); F em is the photon flux emitted by the cell; and J B,e and J B,h are loss factors. Those loss factors are associated with the irreversibility due to photons emitted as a consequence of a radiative transition from band Y to band X being reabsorbed, not by means of an absorption from band X to band Y but as a transition in another band.
As explained in [21], when the absorption coefficients overlap, the loss factors J B,e and J B,h depend also on the cell thickness and can be considered arbitrarily small when the cell thickness tends to zero. In the case of silicon, the use of light trapping techniques [22,23] makes it plausible to consider J B,e = J B,h = 0 as a limiting case to calculate the cell limiting efficiency. The absorption enhancement obtained with these techniques allows assuming that every photon with energies higher than E L is absorbed, so that an absorptivity equal to 1 can be supposed. To simplify the study, we shall also assume that the absorption coefficients α IC and α VI can overlap but they do not overlap with α VC . This hypothesis can be considered as an approximation describing the working case in which photons capable of causing transitions from the valence band to the conduction band are absorbed in the front emitter of the cell (as illustrated in Figure 1) leaving for the IB material the task of dealing with the absorption of below bandgap photons. Under this hypothesis, M e and M h are given by: and, Equation (1) is expanded as where X is the solar concentration, sin 2 θ S = 1/46050 corresponds to the semi-angle of the sun disc seen from Earth and describes the flux of photons per unit of etendùe at temperature T and chemical potential µ for photons with energy between E 1 and E 2 . In Equation (7) c is the speed of light in vacuum, h is the Planck constant and k B is the Boltzmann constant. For a given voltage (V) the following condition must be fulfilled: where µ CV , µ CI and µ IV are the chemical potentials of photons involved in transitions between the VB and the CB; the IB and the CB; and the VB and the IB, respectively. The presence of the factor "−1" in Equation (7), which ultimately derives from the inclusion of stimulated emission in the model [24], imposes that the value of those chemical potentials, which equals the splitting of the quasi-Fermi levels, cannot be larger than the bandgap involved in the electronic transitions (that is, µ CV < E G , µ IV < E H , and µ CI < E L ). Assuming that the IB has a high density of states, the quasi-Fermi level of the IB can be fixed within this energy band (illustrated in Figure 1), enabling both empty states to receive electrons from the VB as filled states to supply electrons to the CB. Equation (1) assumes that only radiative recombination takes place between bands. However, it is well known that silicon cells are affected by non-radiative recombination like Auger recombination, which is the predominant recombination mechanism when silicon cells are illuminated with concentrated light. For this reason, it is worth including this non-radiative recombination mechanism in order to know whether the efficiency potential of an IBSC based on silicon still remains under more realistic circumstances.
In order to include Auger recombination, Equation (1) must be modified as follows: where J Au is the Auger recombination current. This current is calculated assuming flat quasi-Fermi levels, which leads to the following equations: In this expression, W is the thickness of the cell; n, p and n i are electron, holes and intrinsic concentration, respectively; and C e and C h are Auger coefficients [25].
With the conditions mentioned above, Equation (6) is numerically solved in order to find µ CI and µ IV and, therefore, the current voltage-characteristic J(V). Once this current-voltage characteristic is solved, the maximum power point and the efficiency of IBSCs can be calculated.

Results and Discussion
The efficiency limit of silicon-based IBSCs is calculated with the model and considerations described above. To this end the J(V) characteristic is numerically solved, with Equation (6) for the radiative limit and Equation (9) when including Auger recombination. Notice that, in Equation (6), the J(V) characteristic depends on µ CI but not on µ IV . This is because the total current density (J) is defined in terms of the electron current density. J can also be defined as a function of the hole current density with an equation similar to Equation (6), but replacing E L with E H , γ with 1−γ and µ CI with µ IV . By imposing that electron and hole current density must be equal, the values for µ CI and µ IV can be obtained, and consequently the J(V) characteristic, the maximum power point and the efficiency of IBSCs can be calculated. The results are obtained as a function of the E L , swiping the value of this parameter from 0 eV to half of the silicon bandgap, since symmetry with respect to the mid-gap is expected.
Results are plotted in Figures 3 and 4. In those figures the curves labeled "Non-overlapping" correspond to absorption coefficients that do not overlap (γ = 0). Curves labeled "Overlapping" correspond to absorption coefficients α IC and α VI that do overlap (0 < γ ≤ 1), so for each E L , γ is optimized for maximum efficiency. Curves labeled as "Auger" plot the results obtained using Equation (9). For the calculus, cell thicknesses of W = 200 µm and W = 50 µm are assumed and we still consider J B,e = 0 and J B,h = 0 (neglecting losses related to inefficient photon recycling) in order to provide an upper bound of the limit efficiency. The study includes several solar concentrations: One sun (Figure 3), 100 suns and the maximum concentration of 46,050 suns ( Figure 4).

One-Sun Case
Silicon solar cells are mostly used in flat modules operating under irradiances in the one-sun range. Therefore, studying the limiting efficiency of silicon IBSCs under these circumstances is our first goal. From the results plotted in Figure 3, it can be noticed that at one sun there is no efficiency gain for a non-ratched IBSC over the single gap solar cell. The limit efficiency of a silicon non-ratched IBSC is not higher than the limit efficiency of a single gap silicon solar cell, and only when the IB is located at the middle of the gap and the absorption coefficients do not overlap does it equals the limit efficiency of a single gap silicon cell. This is because under those circumstances the IB does not contribute at all, neither to the enhancement of the photocurrent nor to the radiative recombination, so the operation of an IBSC is equivalent to the one of a single gap cell. In this work the efficiency of the silicon cell is not improved by the inclusion of an IB due to the restrictions imposed on the slitting of the quasi-Fermi level. However, it must be pointed out that imposing different conditions like absorption edges whose value adds up to more than the bandgap involved [20] or different energy band diagrams can still enable an efficiency gain [26].
Energies 2020, 13, x FOR PEER REVIEW 6 of 11 imposed on the slitting of the quasi-Fermi level. However, it must be pointed out that imposing different conditions like absorption edges whose value adds up to more than the bandgap involved [20] or different energy band diagrams can still enable an efficiency gain [26]. Results from Figure 3 corroborate what we advanced previously regarding that the efficiency limit of IBSCs with non-optimum bandgaps might be higher with overlapped absorption coefficients. In particular, this is true when 0 < EL/EG ≤ 0.16 (for γ = 0.99) and 0.33 ≤ EL/EG < 0.5 (for 0.11 ≤ γ ≤ 0.5). The reason is an enhancement in the current generation due to a better match between the carriers generated in inter-band transitions from the IB to the CB and from the VB to the IB, as a consequence of the overlapping of the corresponding absorption coefficients (αIC and αVI).

Figure 3.
Limit efficiency of silicon IBSC at one sun for different IB positions with and without overlapping and with and without Auger recombination (discussed in Section 3.3). The value of γ is optimized to maximize the efficiency at every EL and the limit efficiency of a conventional silicon solar cell is shown for comparison.

Concentration Case
The situation changes when concentration is used. At 100 suns the limiting efficiency of silicon IBSC surpasses the limiting efficiency of a single gap silicon solar cell (Figure 4 top). At this concentration, the maximum limiting efficiency of the IBSC reaches 41.8% when the IB is located at an optimum position (EL/EG = 0.31), being significantly higher than the 34.82% of a silicon cell without IB. Overlapping absorption coefficients provide a higher limiting efficiency that the non-overlapping ones when EL/EG > 0.31 and 0 ≤ γ ≤ 0.5. However, those efficiencies do not surpass the maximum efficiency obtained with non-overlapping absorption coefficients and an optimum IB. This is different to what is reported in reference [20], where the highest efficiency is obtained with overlapping absorption coefficients (also at one sun) but for an IB with optical gaps not constrained to + = , as anticipated in the introduction.
At maximum concentration (46050 suns), the obtained results show a similar behavior as in the case of 100 suns. The silicon IBSC increases its limiting efficiency with respect to the cell without IB, but in this case for every value of EL/EG (Figure 4 bottom). Silicon IBSC obtains a maximum efficiency of 54.6% compared with 40.7% for bare silicon. The highest efficiency is obtained when EL/EG = 0.32. As in the previous case, for EL/EG above 0.32, overlapping absorption coefficients provide higher efficiency (for 0 ≤ γ ≤ 0.5) than the non-overlapping ones.  Figure 3 corroborate what we advanced previously regarding that the efficiency limit of IBSCs with non-optimum bandgaps might be higher with overlapped absorption coefficients. In particular, this is true when 0 < E L /E G ≤ 0.16 (for γ = 0.99) and 0.33 ≤ E L /E G < 0.5 (for 0.11 ≤ γ ≤ 0.5). The reason is an enhancement in the current generation due to a better match between the carriers generated in inter-band transitions from the IB to the CB and from the VB to the IB, as a consequence of the overlapping of the corresponding absorption coefficients (α IC and α VI ).

Concentration Case
The situation changes when concentration is used. At 100 suns the limiting efficiency of silicon IBSC surpasses the limiting efficiency of a single gap silicon solar cell (Figure 4 top). At this concentration, the maximum limiting efficiency of the IBSC reaches 41.8% when the IB is located at an optimum position (E L /E G = 0.31), being significantly higher than the 34.82% of a silicon cell without IB. Overlapping absorption coefficients provide a higher limiting efficiency that the non-overlapping ones when E L /E G > 0.31 and 0 ≤ γ ≤ 0.5. However, those efficiencies do not surpass the maximum efficiency obtained with non-overlapping absorption coefficients and an optimum IB. This is different to what is reported in reference [20], where the highest efficiency is obtained with overlapping absorption coefficients Energies 2020, 13, 3044 7 of 10 (also at one sun) but for an IB with optical gaps not constrained to E L + E H = E G , as anticipated in the introduction.
At maximum concentration (46050 suns), the obtained results show a similar behavior as in the case of 100 suns. The silicon IBSC increases its limiting efficiency with respect to the cell without IB, but in this case for every value of E L /E G (Figure 4 bottom). Silicon IBSC obtains a maximum efficiency of 54.6% compared with 40.7% for bare silicon. The highest efficiency is obtained when E L /E G = 0.32. As in the previous case, for E L /E G above 0.32, overlapping absorption coefficients provide higher efficiency (for 0 ≤ γ ≤ 0.5) than the non-overlapping ones.
An important conclusion worth emphasizing from these results is that IBSCs with overlapping absorption coefficients exhibit higher limiting efficiency when E L is higher than the optimum value. This result draws attention to the fact that, when dealing with non-optimum bandgaps, an IBSC in which overlap exists can give more efficiency than the case in which overlap does not exist. Therefore, the fact that the condition of non-overlapping is required for attaining the IBSC with the absolute maximum efficiency must not mislead researchers to pursue the same condition when non-ideal bandgaps are used. In addition, these results imply that pursuing silicon IBSCs with optical gaps constrained to E L + E H = E G only makes sense if used under concentrated light.
Energies 2020, 13, x FOR PEER REVIEW 7 of 11 An important conclusion worth emphasizing from these results is that IBSCs with overlapping absorption coefficients exhibit higher limiting efficiency when EL is higher than the optimum value. This result draws attention to the fact that, when dealing with non-optimum bandgaps, an IBSC in which overlap exists can give more efficiency than the case in which overlap does not exist. Therefore, the fact that the condition of non-overlapping is required for attaining the IBSC with the absolute maximum efficiency must not mislead researchers to pursue the same condition when non-ideal bandgaps are used. In addition, these results imply that pursuing silicon IBSCs with optical gaps constrained to + = only makes sense if used under concentrated light. An important conclusion worth emphasizing from these results is that IBSCs with overlapping absorption coefficients exhibit higher limiting efficiency when EL is higher than the optimum value. This result draws attention to the fact that, when dealing with non-optimum bandgaps, an IBSC in which overlap exists can give more efficiency than the case in which overlap does not exist. Therefore, the fact that the condition of non-overlapping is required for attaining the IBSC with the absolute maximum efficiency must not mislead researchers to pursue the same condition when non-ideal bandgaps are used. In addition, these results imply that pursuing silicon IBSCs with optical gaps constrained to + = only makes sense if used under concentrated light.

Auger Recombination
The discussions in Sections 3.1 and 3.2 have been focused on the results obtained for the radiative limit. When Auger recombination is included non-radiative losses reduce the efficiency of silicon solar cells without IB (case in which E L = 0) by approximately four points for all the considered irradiances: 1100 and 46,050 suns. However, the reduction in efficiency due to Auger recombination is less severe in IBSCs. Indeed, at 1 and 100 suns, the efficiency in silicon IBSC remains intact for E L higher than the optimum value in which the maximum efficiency is reached. All these results are shown in Figures 3 and 4 and the most remarkable ones are summarized in Table 1.
The reason why Auger recombination affects IBSCs less is because the presence of the IB in these cells already introduces an additional radiative recombination that, at one sun, is higher than the recombination introduced by Auger and, at 100 suns, is comparable. This is illustrated in Figure 5, where higher recombination rates are manifested as higher reductions in the output voltage, since this is the parameter more sensitive to that (the current density is not). At maximum concentration, the radiative recombination introduced by the IB is reduced enough to obtain no significant reduction in the output voltage (black triangles in Figure 5). So, the performance of IBSCs can be affected by the additional recombination introduced by Auger at that concentration (blue symbols in Figure 5).
Energies 2020, 13, x FOR PEER REVIEW 8 of 11 Section 3.3). The value of γ is optimized to maximize the efficiency at every EL and the limit efficiency of a conventional silicon solar cell is shown for comparison.

Auger Recombination
The discussions in Sections 3.1 and 3.2 have been focused on the results obtained for the radiative limit. When Auger recombination is included non-radiative losses reduce the efficiency of silicon solar cells without IB (case in which EL = 0) by approximately four points for all the considered irradiances: 1100 and 46,050 suns. However, the reduction in efficiency due to Auger recombination is less severe in IBSCs. Indeed, at 1 and 100 suns, the efficiency in silicon IBSC remains intact for EL higher than the optimum value in which the maximum efficiency is reached. All these results are shown in Figures 3 and 4 and the most remarkable ones are summarized in Table 1. 3) The reason why Auger recombination affects IBSCs less is because the presence of the IB in these cells already introduces an additional radiative recombination that, at one sun, is higher than the recombination introduced by Auger and, at 100 suns, is comparable. This is illustrated in Figure  5, where higher recombination rates are manifested as higher reductions in the output voltage, since this is the parameter more sensitive to that (the current density is not). At maximum concentration, the radiative recombination introduced by the IB is reduced enough to obtain no significant reduction in the output voltage (black triangles in Figure 5). So, the performance of IBSCs can be affected by the additional recombination introduced by Auger at that concentration (blue symbols in Figure 5).   Figure 5 is also useful to understand why the inclusion of an IB does not always provide a gain in efficiency, despite providing an increase in current density (black stars in Figure 5). This is the case of silicon IBSC at one sun, where the reduction in power due to the relative reduction in voltage is Energies 2020, 13, 3044 9 of 10 slightly higher than the increase in power due to the relative increase in current density. When the concentration increases from one sun to 46,050 suns, the efficiency increases mainly due to the increase in voltage, which approaches the voltage of a conventional cell without IB. The current density at the maximum power point also shows an increase with concentration, which is also related to the fact that under concentration higher FF are obtained (in the ideal cases with zero series resistances).

Conclusions
In this work the efficiency of silicon-based IBSCs under different concentrations of sunlight is studied, including for some realistic parameters such as overlapping in absorption coefficients and non-radiative Auger recombination. Results reveal that silicon non-ratchet IBSCs can surpass the limit efficiency of conventional silicon solar cells only under concentrated sunlight. The maximum efficiency of silicon IBSCs is achieved for an optimum bandgap distribution (E H and E L ), assuming selective absorption coefficients. However, if non-optimum bandgaps are considered, overlapping absorption coefficients can lead to higher efficiencies than selective ones. The performance of silicon IBSC is less affected by Auger recombination than conventional silicon cells except if extremely high concentration is used to illuminate the devices.

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