Spectral response measurements of multijunction solar cells with low shunt resistance and breakdown voltages

Spectral response measurements of germanium‐based triple‐junction solar cells were performed under a variety of light and voltage bias conditions. Two of the three junctions exhibited voltage and light bias‐dependent artifacts in their measured responses, complicating the true spectral response of these junctions. To obtain more insight into the observed phenomena, a set of current‐voltage measurement combinations were also performed on the solar cells under identical illumination conditions, and the data were used in the context of a diode‐based analytical model to calculate and predict the spectral response behavior of each junction as a function of voltage. The analysis revealed that both low shunt resistance and low breakdown voltages in two of the three junctions influenced the measured quantum efficiency of all three junctions. The data and the modeling suggest that combination of current‐voltage measurements under various light bias sources can reveal important information about the spectral response behavior in multijunction solar cells.


Introduction
With significant advances in high-efficiency multijunction solar cell (MJSC) technologies, the number of subcells within each structure has steadily increased over the years with the goal of maximizing the operating voltage while preserving the other performance parameters such as the short circuit current (I sc ) and the fill factor [1][2][3][4][5][6]. Recently, five-junction solar cells based on semiconductor bonded technology have demonstrated confirmed non-concentrator terrestrial air mass 1.5 global spectrum power conversion efficiency of ≈ 39% [7] and concentrator 4-junction cell efficiencies of 46% [8]. With increased complexities in design and architecture of these types of cells, the task of performing electrical characterization and spectral response (SR) measurements of the cells will also become more involved. Although there are well-established procedures for performing current-voltage (I-V) or spectral response measurements in MJSCs [9,10], it has been shown previously that certain factors such as low shunt resistance [11][12][13][14][15][16][17], low reverse breakdown voltage [10,18] or luminescence coupling [12,[19][20][21][22][23][24][25] can cause artifacts in the spectral response or the quantum efficiency (QE) of the device.
Luminescence coupling is usually a factor between two adjacent junctions and its effects generally appear as an increased QE signal in the wavelength range corresponding to energies above the bandgap of the current-limited junction [24][25][26]. It has also been reported that luminescence coupling has a light bias intensity dependence, resulting in a decrease in QE of the junction of interest with an increase in the light bias intensity of the junction above it [12]. This reduction factor is substantially dependent on the LC coupling coefficients between the junctions. When considering artifacts originating mostly from shunt and breakdown voltage issues, most of the work reported to-date in QE measurements of MJSCs, particularly 3-junction cells based on a germanium (Ge) bottom cell (BC), has been focused on cases where only a single junction has been found to behave nonideally. If the nonideality is due to a low breakdown voltage only, then an applied voltage bias of enough magnitude is able to fully recover the QE of the junction of interest without much dependence on the light bias intensity [18].

Spectral response measurements of multijunction solar cells with low shunt resistance and breakdown voltages
If the junction of interest suffers from a low shunt resistance, increasing the light bias intensity of the nonlimiting junctions and application of a voltage bias is generally needed to achieve a good measurement of the QE signal [14]. There are even reported cases where a combination of luminescence coupling and low shunt resistance affect the QE response of the Ge junction [12].
Although generalized opto-electronic models have been proposed that combine and explain a variety of simultaneous effects in MJSCs measurements [21], it is still very illuminating to consider simpler analytical models to explain the influence of the most severe artifacts in certain devices. Our aim has particularly been focused on cells where two junctions revealed both low shunt resistance and low breakdown voltage effects and when LC has not been a major factor in causing artifacts in the QE response of the junctions. Even though several previous works have been dedicated to this subject and some variations of it [10,12,[14][15][16][17][18], we found that significant insight can be obtained into the behavior of each junction when we first investigated the I-V response of one of these types of cells under several combinations of light sources where each light source is a high-power light emitting diode (LED) with a unique wavelength emission peak suitable to light bias a single junction. For a triple junction cell, this approach involves eight combinations of measurements (3 single-LED illuminations, 3 double LED combinations, one triple LED illumination, and one with no illumination). These data can then be used within the context of a straight-forward two-diode equivalent circuit (EC) model to determine a set of parameters that will fit all eight data combinations at once. This technique eliminates the need for separate iso-cell measurements for parameter estimates because there are sufficient data to allow for obtaining good fits to the data. After the EC model parameters are determined, the results are used to predict SR and QE of each junction with similar levels of light bias as a function of voltage and compare the findings to the experimental results. We generally find very good agreement between the predicted model results and the experimental data. When more than one junction suffers from nonideal shunt and breakdown voltage effects generally caused during manufacturing or nonoptimal material processing, the ideal diode-like behavior expected from a wellfabricated p-n junction is nonexistent. In these cases, the QE measurement of the intended junction will reveal artifacts originating from all of the nonideal junctions, in a combined way. Furthermore, in almost all cases studied, the magnitude of the measured QE itself, as measured across the two terminals of the entire stack, will be less than the individual junctions' actual QE response, even under the best voltage bias conditions.

Experimental
For this study, the samples were monolithic, two-terminal triple junction GaInP/GaAs/Ge solar cells of commercial grade with a nominal area of 2.3 cm 2 and air mass (AM) 1.5 G 1-sun efficiency of 27% at 25°C [4]. The currentvoltage and the spectral response data reported here are representative of multiple solar cells tested within the batch received. Combinations of current-voltage measurements were performed with a source-meter electrometer and three high-power LEDs with projector optics illuminating the entire area of the cell with bright illumination beams of good uniformity. The spectral irradiance of all three LEDs at the incident plane of illumination, was measured by a calibrated spectroradiometer. The irradiance has a Gaussian form, with the three-center wavelengths at 520 nm, 741 nm, and 1077 nm. All three LEDs were operated with battery sources close to their rated source current values for better control over signal-to-noise ratios when performing differential QE measurements. The reported irradiance values for each LED were obtained by integrating the spectral irradiance curves over wavelength. For QE, a dual-light source monochromator set up was used, utilizing the principle of differential spectral response (DSR) measurement as described previously [27] with calibrated monitor detectors from 300 nm to 1900 nm, allowing for simultaneous measurements of the cell photocurrent and the incident monochromatic beam power. The monochromator beam was chopped by a mechanical chopper, hence generating a small AC photocurrent that is collected via a transimpedance preamplifier with the ability to separate AC and DC components of the signal and built-in voltage bias capability. The AC signal of the cell was measured by a lock-in amplifier. The voltage bias was applied in the forward bias direction of the cell operation from 0 V to values slightly higher than the open circuit voltage, V oc , of the cell under the illumination conditions to verify modeling results. Although application of higher voltage amounts is possible, it can lead to significant instability in the recorded AC lock-in signal. The monochromator measurement uncertainty for spectral response measurements under low voltage bias is generally better than 1%. This includes the uncertainty of the reference photodiode and the transfer calibration onto the monitor diode, the repeatability and reproducibility measurements and the standard deviation of the collected data. However, at bias voltages approaching the V oc , it can reach 5% or more, mostly due to signal instability.

The Modeling Framework
In this section, we first describe the mathematical model used to fit the composite I-V curve data under each unique lighting condition. Then, by use of a simplified analysis, we demonstrate how the true value of the spectral response of a given junction is affected by various artifacts.

Mathematical model for each junction
We consider a typical two-diode model for each junction [28], with its corresponding shunt and series resistances. The implicit current density versus voltage J-V curves for each junction i = 1,2,3 can be written in the forward bias direction: and in reverse bias direction: where J i ph is the photocurrent density of the ith junction, J i 01 is the saturation current density, J i 02 is the saturation current density for the depletion region, A is the device area, n i is the ideality factor, R i sh is the shunt resistance, R i s is the series resistance, n i m is the Miller index and V i bd is the breakdown voltage. In equation (2), a multiplicative factor is added to account for the cell's avalanche breakdown in the reverse bias [28]. We define the differential spectral responsivity of each junction as the difference in its short circuit current between the two chopper states of closed and open, r, divided by the power of the incident monochromatic beam. The difference between these two current densities is simply: Note that in reality, in a MJSC, only the voltage across the entire cell stack can be fixed, say at 0 V, leaving each junction under a different operating voltage. However, this definition is consistent with single-cell measurements where the short circuit current is the quantity of interest for spectral response measurements. For a very small value of r, which is the case for differential monochromator-based measurements, this change is equivalent to a small displacement of the entire curve, J i (V) in an amount equal to r.
To calculate the composite J(V) curve for the entire cell, which constitutes the series sum of the individual junction J i (V), first the best fit parameters for each implicit equation corresponding to each junction are selected and the J i versus V functional form is calculated. Inverting these functions and obtaining V(J i ) curves then allow for summation of the voltages at fixed current values. Finally, this composite curve can be inverted back to the J(V) form. As it will be described in Discussion and Results, the objective of this modeling was to obtain the best fit parameters from as many combination I-V curves as possible so that a robust estimate of the voltage bias dependence of the spectral responsivity of each junction can be made. Here, by best fit, we mean the best visual fit of the model to the experimental data, that is, no least square calculations were employed.

Factors limiting the spectral response measurements
The mathematical model given by equations (1) and (2) allows for determination of the relationship between the measured spectral response quantity, which is influenced by the whole composite curve and measured at the terminals, and the real increase r of the zero-voltage current for the junction of interest. We have observed that in cases where one or more junctions have nonideal diode properties, application of a voltage bias or light bias will not fully recover the entire value r, as desired. Although this concept can be mathematically explained by use of equations (1) and (2), here we present a simplified case that demonstrates this point well. Figure 1 shows an example of an MJSC case with a well-behaved GaAs junction but with faulty Ge and GaInP junctions with low shunt resistances and low breakdown voltages. The plotted J-V curves are representative of the sort of behavior we have seen in our measurement and modeling work. This plot represents a case where the intention is to extract the SR of the GaInP junction by light biasing the Ge junction (with the 1077 nm LED) and the GaAs junction (with the 741 nm LED), while light modulating the GaInP junction with the much weaker 520 nm monochromatic light. For the sake of clarity, the J-V curve corresponding to the open chopper state is drawn in an expanded view with an exaggerated r value, but in reality, r is many orders of magnitude smaller than the other two junctions' light bias-induced short circuit current values. Given this, the composite J-V curve can be constructed for both chopper states, showing that points The difference ΔJ between these two current densities is the amplitude of the AC current density that is actually being measured: J chopper open   ΔJ is always lower than the change r in the short circuit current density of the junction being measured, but gets close to it only as we tend to certain values for the slopes m of the J versus V curves. This behavior can be easily understood by observing the following limit: Equation (12) indicates that the measurement of the spectral response of one junction, such as the GaInP junction, can be affected by the J-V characteristics of the other junctions as well. Furthermore, it can be seen from equation (13) that having a higher intensity light bias on the Ge junction would make m Ge larger, therefore increasing the ΔJ closer to the r limit, in essence increasing the measured spectral response of the GaInP junction. The section on voltage bias dependence of spectral response discusses this issue in detail.
The responsivity of a junction is defined as: where P is the monochromatic beam power. However, the measurement gives the following value: Therefore, the ratio of these two quantities, R m ∕R describes a factor that effectively reduces the spectral response from its maximum expected value: In the section on the voltage bias dependence of spectral response, it will be shown that this ratio is close to 1 for well-behaved junctions and lower than 1 for bad cases of low shunt resistance or breakdown voltages. This ratio can be extracted from combinations of I-V measurements at a representative wavelength for each junction, prior to performing any spectral response measurements.

Discussion and Results
Performing a set of I-V measurements of the cell under all possible combinations of the three LED light sources prior to performing QE measurements can significantly help in determining: (1) which junctions are likely to be affected by artifacts, (2) what level of light and voltage bias should be used to obtain the best measurement results, and finally (3) how the best QE measurement results compare with model predictions for the maximum QE that each junction is able to produce internally.

Combinations of I-V measurements
In order to use the mathematical model described in Mathematical model for each junction to predict the SR behavior, the values of all the parameters in equations (1) and (2) need to be determined as accurately as possible. This was accomplished by performing eight combinations of I-V measurements on the solar cell, with three unique LEDs appropriate for generating a photocurrent signal in each of the junctions. These were the 520 nm LED for excitation of the GaInP junction, the 741 nm LED for the GaAs junction, and the 1077 nm LED for the Ge junction. The best fit parameters for each junction are chosen so that the calculated composite I-V curve gives the best fit to the measured cell's I-V curve for each of the eight unique measurements. Figure 2 Note that if the cell has either one or no nonideal junctions, it will suffice to perform just one I-V curve measurement with all three LEDs on to obtain a good estimate of the parameters. However, in a case such as this where two junctions have poor diode qualities, data from additional combinations of I-V measurements significantly improve the estimate of the parameters. The modeling results here indicate that both the GaInP junction and the Ge junction suffer from low shunt resistance and low breakdown voltages [14,18].
The two-diode model seems to provide a reasonably good fit to the various composite I-V curves. However, even a single-diode model may be sufficient for obtaining reasonable fits to the data. The data fitting, especially in the forward bias direction, could potentially be improved further with additional circuit elements and fit parameters, but those additions will add complexities to the simple model presented here. Table 1 reports the values of the fit parameters for each junction. As expected, the shunt and breakdown voltage values for the GaInP and the Ge junctions are lower than what has typically been reported for solar cells [14] The photocurrent density values J ph correspond to cases where the junctions are illuminated by their respective LEDs, that is, J ph would be zero for an LED off-state for a particular junction. The Miller indices are set to 3, as suggested elsewhere [28].

Voltage bias dependence of spectral response
Once the I-V curve parameters are determined, the measured spectral response behavior as a function of voltage can be calculated and compared to the measured data. This calculation is based on the modeled junction I-Vs and not based on equation (12) which was derived with a simpler explanation of the observed effects. Furthermore, the important ratio R m ∕R can also be plotted as a function of V b to see the discrepancy between the measured SR of a junction versus the predicted response from the modeled individual I-V curve. Since the SR is a differential measurement, the mechanical motion of the chopper results in a pair of I-V composite curves, one for when the chopper blocks the monochromatic light and one for when it unblocks it. The SR calculations are straightforward. Taking the GaInP junction as an example, the individual I-V curves for GaAs and Ge junctions are simulated under the 741 nm and the 1077 nm LEDs 'on', and the I-V curve for GaInP junction is simulated under dark (520 nm LED 'off'). A second simulation of this case is run again, with the singular difference that a small photocurrent value is picked for the GaInP junction to simulate its response under the open chopper scenario. This value is many orders of magnitude lower than the other two photocurrent values, just as is the case in actual measurements. From this result, the change due to photogenerated current in the GaInP junction at 0 V, r, is calculated from the simulated subcell J-V curve. This of course is proportional to the real responsivity of the junction. However, the actual quantity that is measured, R m , is proportional to ΔJ, the difference in the composite current between the two chopper states of on and off, as previously shown in Figure 1.
ΔJ is also obtained from the simulation as a function of bias voltage. Plotting the calculated ratio ΔJ∕r as a function of voltage is very revealing; for example, it can show whether a true, artifact-free spectral response for a given junction can be obtained from a measurement. LED. The left Y-axis represents the measured and calculated responsivities and the right Y-axis simultaneously shows the ratio R m ∕R, obtained purely from the calculations. Remembering that a R m ∕R = 1 ratio means that the measured responsivity is equal to the actual junction responsivity, that is, photocurrent generated in the junction is the photocurrent measured at the terminals, two important observations can be made here: First, the diode-based model predicts the V b dependence of the responsivity relatively well, predicting a peak response that can be achieved by an application of V b ≈ 1.8 V in all three cases with the exception that both combinations of the 1077 and 741 nm LEDs are needed to obtain the maximum response. Note Table 1. The model fit parameters obtained from fitting the measured composite I-V curves to the model based on equations (1) and (2). Parameter n is higher than the typical value of 2, most likely due to Ohmic or recombination current distributions [28].  that beyond the peak value, the responsivity should drop to 0, as the data also suggest, though the DSR measurement becomes very unstable at higher V b . The second observation is that in Figure 3A, R m ∕R is reaching a maximum value at 0.92 for the GaInP junction, meaning that no matter what voltage bias is used, 8% of the GaInPgenerated photocurrent can never be recovered at the external terminals. Figure 3D shows the V b dependence of the responsivities for the GaAs junction, with the monochromator wavelength at 741 nm, and both the light bias 520 and 1077 nm LEDs on. The data and the model are in good agreement here, showing no significant V b dependence until the responsivity drops at higher voltages. This behavior is consistent with the high breakdown voltage and high shunt resistance of the GaAs junction. Yet, the R m ∕R ratio for this junction is at 0.96, indicating that the poor behavior of the other two junctions slightly affects the maximum measured responsivity of an otherwise ideal junction. The other two light bias cases for this junction were similar and are therefore not shown.
Finally, Figure 3E-G show the V b dependence of the responsivities for the Ge junction, with the monochromator wavelength at 1077 nm and the light bias scenarios: (E) both 520 and 741 nm LEDs on, (F) with only 741 nm LED on, and (G) with only 520 nm LED on. Although there are small discrepancies between the modeled and measured results, the findings are similar to the GaInP case in that a peak response is achieved at a voltage bias higher than 0 V, typically around 2 V. One notable exception is that unlike the GaInP junction measurements, the maximum R m ∕R for the Ge junction is 0.55, indicating that a much smaller photogenerated current value is extracted at the terminals for this junction. A portion of this 45% photocurrent loss can be recovered for the Ge junction if the 520 nm LED intensity were increased.
We believe that the discrepancies observed between the model and the data, particularly at the low voltage values for Figure 3A and B and the shifts in the peak positions for Figure 3E and F, originate from imperfections of the double-diode I-V model in capturing all the relevant physics inside these imperfect devices. Since we are using the parameters strictly obtained from I-V curve fits, we do not make any alterations or adjustments to the individual R m ∕R ratio plots afterwards. Actual spectral response measurements are differential lock-in based measurements and furthermore, the monochromatic light intensity is very low and hence it is possible that there are nonlinearity or second-order effects present between the monochromator "on" and "off" states that we have not captured in our simulations; that is, it is possible the shape of the I-V curve for the chopper-on and chopper-off states change slightly. There is no way for us to know this, because our instruments are not sensitive to this while there is a large light bias signal present on the other junctions. The purpose of this section is to show that it is possible to understand the voltage bias behavior but the exact source of discrepancy between the model and the data is beyond the scope of this work.

Complete wavelength measurements
Now that the voltage dependence of the SR has been described at fixed wavelengths within the regions of responsivity of each junction, we now turn our attention to the full wavelength measurements for all three junctions. Figure 4A shows the measured spectral responses of all three junctions under the best voltage bias conditions and light bias conditions corresponding to LED light intensities of 116 W/m 2 , 228 W/m 2 and 31 W/m 2 for the 520 nm, 741 nm, and the 1077 nm LEDs, respectively. Notice that the nonzero photocurrent response for the GaInP top cell in the region > 700 nm is likely due to some leaked photocurrent from the GaAs middle cell due to the low shunt resistance values of the Ge and the GaInP subcells. However, the nonzero response (0.03 A/W) for the Ge junction for < 860 nm is related to LC from the middle junction. The calculated R m ∕R ratios for each junction from the previous section are used as correction factors for each corresponding curve. The corrected spectral response results are plotted as in Figure 4B. The largest correction factor of 1/0.55 for the Ge junction substantially increases the responsivity of that junction, to values expected from a well-behaved Ge junction. In Figure 4C, the corrected QE curves are also shown for reference.
Furthermore, the analytical model allows one to predict the magnitude of the light bias needed to obtain the best QE measurement for the affected junctions. Focusing on the Ge junction's response, we know that the poor diode behavior of the GaInP junction is partially to blame for Ge's low SR values. Applying a stronger light bias to the GaInP junction should shift the operating voltage point for that junction farther to the right towards a higher shunt resistance region, that is, m GaInP → ∞, causing ΔJ → r Ge , hence approaching the corrected SR curve. To show this trend, the 520 nm LED intensity was progressively increased from the initial value of 116 W/m 2 to 147 W/m 2 and 243 W/m 2 , while the 741 nm LED was kept fixed at 228 W/m 2 , partially to maintain the effect of LC between the GaAs and the Ge junctions at a fixed level. To calculate the model predicted R m ∕R ratios, a linear relationship was assumed between the 520 nm LED light intensity and the photogenerated current in the GaInP junction as a result. The predicted R m ∕R ratios for the Ge junction should rise from 0.55 to 0.63 and finally to 0.81 for the three 520 nm LED light intensity settings. Figure 5A and B show the results of these three measurements for the Ge junction, and the corrected Ge responsivities based on these calculated ratios. The plots indicate good agreement between the model and the measurements and show that all three distinct curves can be collapsed onto one curve, within the context of the simple analytical model. Since the 741 nm bias light was kept constant for all three measurements and only the 520 nm light was changed, the luminescence coupling between the middle and the bottom junctions were mostly unchanged as evidenced by the plateau below 870 nm. So the shunt resistance effect is the dominant mechanism affecting the QE response of these devices.

Conclusions
Combinations of I-V measurements under various types of light bias sources were used within the context of an analytical two-diode I-V curve model to predict the voltage bias and light bias dependence of the differential spectral response measurements in a triple junction solar cell. The solar cells studied for this work were Ge based, with both the GaInP and the Ge junctions presenting severe cases of low shunt resistance and low reverse bias breakdown voltages. We have shown that these conditions can significantly affect the spectral response measurements of all three junctions, with the Ge junction affected the most. Our work suggests that it is valuable to perform these sets of current voltage measurements and subsequent analysis prior to any spectral response measurements so that one can choose a priori how to set up, or what to expect from, the DSR measurements.