Research Article  Open Access
Enhanced Efficiencies for HighConcentration, Multijunction PV Systems by Optimizing Grid Spacing under Nonuniform Illumination
Abstract
The design of a triple junction solar cell’s front contact grid can significantly affect cell conversion efficiency under high concentration. We consider one aspect of grid design, choosing a linear grid within a distributed resistance cell model to optimize finger spacings at concentrations between 500 and 2500 suns under uniform and nonuniform illumination. Optimization for maximum efficiency under Gaussian irradiance profiles is done by SPICE analysis. Relative to the optimized uniform illumination designs, we find enhancements of 0.5% to 2% in absolute efficiencies for uniform spacing. Efficiency enhancement with nonuniform spacing under nonuniform illumination is also evaluated. Our model suggests that, at lower concentrations (<1000 suns), the penalty for using uniformly spaced fingers instead of nonuniformly spaced fingers is <0.1%. However, at a concentration of 2500 suns the penalty increases to 0.3%. Thus, relative to a uniform irradiance optimization, an absolute efficiency increase of 2.3% can be attained for an optimized nonuniform spacing given the Gaussian irradiance profile under consideration.
1. Introduction
Multijunction solar cells based on IIIV semiconductors are the most efficient of all photovoltaic technologies, with a record efficiency of 44.7% under concentration recently reported [1]. Not only does concentration permit the balance of concentrator photovoltaic system costs to offset the expensive material costs of multijunction solar cells, but it also enhances cell efficiency through the logarithmic increase in open circuit voltage [2]. The effects of concentration lead to a nonuniform profile that is not effectively flattened by secondary optic solutions [3]. Concentrator PV system optics often generate profiles that peak sharply at the center of the cell (Figure 1(a)) and follow a Gaussian distribution [4]. Peaktoaverage irradiance ratios (PAR) typically range between 2 and 10 [2], but greater PARs are not unknown [5]. The nonuniformity in irradiance leads to a nonuniformity in current production that can decrease cell efficiency [6, 7].
(a)
(b)
The decrease in cell performance becomes increasingly evident at higher concentrations, where increasing current densities drive higher series resistance () losses. The impact of these losses is a reduction in fill factor and therefore a reduction in conversion efficiency [8]. For uniform irradiance, a linear grid with even finger spacing, as shown in Figure 1(b), achieves maximum extraction efficiency regardless of such losses. This is not the case for nonuniform irradiance, where the consequent nonuniform current production within each subcell also leads to lateral currents that exacerbate carrier recombination. Continued reliance on uniform irradiance grid designs is therefore a cause for concern [6], motivating analyses of more realistic scenarios.
Simulation tools such as TCAD Sentaurus allow for valuable understanding of how solar cell performance depends on material properties and device structure [9]. However, as their complexity drives long computational times, usage of such tools is limited to simulating a small symmetrical element of the device. Therefore, electrical models designed in SPICE (simulation program with integrated circuit emphasis) have been used to characterize the spatial dependence of solar cell performance. Due to their inability to simulate the effects of distributed resistance with varying irradiance profiles, classical lumped SPICE models do not permit realistic simulation. Hence, distributed resistance models, capable of simulating actual cell areas, are required for evaluating the effects of nonuniform illumination.
In this paper, we use a twodiode equivalent circuit for each subcell in a twodimensional (2D) distributed resistance model to address nonuniform current production under various Gaussian irradiance profiles. We report on the efficiency enhancements obtained by grid spacing optimization using realistic irradiance profiles at concentrations of 500, 1000, 2000, and 2500 suns with both uniform and nonuniform grid spacing. The results are compared with those obtained from the corresponding uniform irradiance analyses.
2. Numerical Model
2.1. Device Structure and Parameter Extraction
In order to generate values for equivalent circuit parameters that act as an input to the SPICE model,  curves for each subcell are required. These can be obtained either by using experimental data from isotypes or by simulating solar cell structures using TCAD tools [9, 10]. In this work, we have designed a triplejunction solar cell (TJSC) structure using TCAD Sentaurus by Synopsys. The simulation has been performed on a symmetrical element of a latticematched TJSC with 5% shading due to an ohmic contact on a highly doped GaAs cap layer. The structure consists of three subcells composed of GaInP, InGaAs, and Ge layers on a Ge substrate [9] connected in series using lowresistance tunnel junctions [11]. Previous variants of the model have been validated with experimental results. Since a twodiode model yields more accurate results than a onediode model [12], solar cell parameters have been determined for each subcell by curvefitting the subcell  curves to where the individual diode equations are is the lightgenerated current (which is directly proportional to the solar irradiation), and are reverse saturation currents, is the electronic charge, is the Boltzmann constant, is the series resistance, is the shunt resistance, is the output voltage, and is the output current of the solar cell.
Since the values of , , and are known from the subcell  curves and from the sheet resistances of the TJSC structure, the remaining unknowns, , , , and , can be determined by curvefitting using the leastsquares method [13]. The parameters thus obtained are shown in Table 1 and have errors of less than 2%. Figure 2 shows the fitted subcell currentdensity versus voltage () curves along with the simulated () curves obtained for this structure using TCAD Sentaurus.

2.2. 2Dimensional Distributed Resistance Model
Equivalent circuit representation of solar cells is useful for CPV system analysis [14–17]. In this work, a 2D distributed resistance model for a TJSC is used [13]. This model, designed using LTSPICE, is an improvement over the 2D model described by Nishioka et al. [14] as it includes equivalent circuits for the directly illuminated regions and for the shaded areas underneath the metal contact. The inclusion of equivalent circuits for the shaded areas is essential to analyze the effects of nonuniform illumination on cell efficiency [16].
Figure 3 presents a detailed circuit diagram for the SPICE model depicting illuminated and shaded equivalent circuits that comprise the repeated functional block. Since finger spacing is variable (given constant finger width), the number of functional blocks varies with spacing for each simulation and is given by where is the total cell area, is the finger spacing, and is the finger width. Narrower spacing requires a larger number of functional blocks to represent the physical extent of the TJSC being modeled.
The linear grid means a 2D model will suffice to describe the TJSC. The model incorporates lateral resistance components connecting the functional blocks across the rows (), tunnel junction resistances () connecting the subcells in series, and contact resistances () for the electrodes. The number of function blocks for a fixed cell dimension determines the lateral resistance values [13, 18]. Electrode resistances are assumed to be negligible, finger widths are fixed at 5 μm, and the active cell area is equal to 1 cm^{2}. The model therefore accounts for shading losses due to the fingers and ohmic losses due to series resistances. Bus bars are assumed to be outside the active cell area and are not included in the model. Shading losses and series resistances due to bus bars are neglected for our simulations.
3. Optimization
In order to optimize finger spacing, the solar cell conversion efficiency is assumed to be a unimodal objective function with one maximum. A derivativefree optimization approach is adopted to maximize cell efficiency. In order to reduce the number of iterations and for a faster convergence, a Goldensection search algorithm is used. Finger spacings range from 50 μm to 300 μm, for a shading range of ~2% to 10%, which is typical for solar cell designs. Convergence of the optimization process is ensured by requiring the absolute change in efficiency to be less than 0.01%; spacing increments are 3 μm. Eight to ten iterations are sufficient to achieve convergence. Section 3.1 below describes the optimization process for an ideal uniform irradiance distribution; Section 3.2 describes the process for a nonuniform irradiance distribution.
3.1. Under Uniform Illumination
Figure 4 shows efficiency as a function of finger spacing at concentrations of 500, 1000, 2000, and 2500 suns under a uniform illumination profile (1 sun = 900 W/m^{2}). The relationship between the shortcircuit current and the concentration is assumed to be linear [2]. This implies that the shortcircuit current for this cell at 1000 suns is 1000 times that of the shortcircuit current under one sun operation. The figure reveals that the efficiency of the cell reaches a maximum value at each concentration for a specific finger spacing.
The optimal finger spacing at 500 suns is found to be 120 μm. However, at 1000 suns the optimal finger spacing is reduced to 95 μm. It is evident that this cell can perform almost optimally, without much loss in efficiency (<0.2%), at both 500 and 1000 suns if the fingers are spaced from 95 to 120 μm. At 2000 and 2500 suns, however, the cell needs to be designed with more closely spaced fingers (optimal spacing is 75 μm and 72 μm resp.) in order to get the best possible efficiency. A similar trend is reported in an NREL technical report [19]. The highest efficiencies obtained with this cell under uniform illumination are 37.9% at 500 suns, 38.1% at 1000 suns, 38.0% at 2000 suns, and 37.9% at 2500 suns.
An increase in the slope of the efficiency versus spacing curve with concentration indicates that the penalty for using a nonoptimal spacing is greater at higher concentrations (2000 and 2500 suns) than at lower concentrations (500 and 1000 suns). If the cell is designed with a 120 μm spacing and operated under uniform illumination at a concentration of 2500 suns, a 1.3% absolute loss in efficiency is observed as compared to that obtained with the optimal spacing.
3.2. Uniform Spacing under Nonuniform Illumination
To model nonuniform illumination, three Gaussian profiles with different peak to average ratios (PARs) are used to represent various types of concentrating optics. Efficiency optimization is performed for PARs of 1.6, 2.6, and 5.3 at concentrations of 500, 1000, 2000, and 2500 suns across the range of finger spacings. A PAR of 1.6 approximates an ideal concentrating system, while a PAR of 5.3 may represent a system with only a primary Fresnel lens element and no secondary homogenizer. Note that, while the total intensity of each irradiance profile is kept constant, chromatic aberration effects are ignored in the optimization. This implies that current mismatch between top, middle, and bottom subcells is neglected. Figure 5 shows how the efficiency variation with finger spacing depends on the illumination profile; uniform illumination has also been included, for comparison. The three Gaussian irradiance distributions are shown in the inset. The optimal finger spacing is seen to decrease with both increasing nonuniformity and increasing concentration, and the efficiency is likewise reduced. Figure 5(a), at 500 suns, shows that, for the highest PAR, the efficiency at the optimal uniform spacing of 120 μm is about 0.5% absolute lower than it would be at its optimal 90 μm spacing; for Figure 5(b), at 1000 suns, the efficiency at the optimal uniform spacing of 98 μm is about 0.9% absolute lower than it would be at its optimal 69 μm spacing; and for Figure 5(c), at 2000 suns, the efficiency at the optimal uniform spacing of 75 μm is about 1.5% absolute lower than it would be at its optimal 50 μm. The results clearly illustrate that an increase in concentration and in nonuniformity has an increasingly significant effect on the conversion efficiency of the solar cell.
(a)
(b)
(c)
3.3. Nonuniform Finger Spacing under Nonuniform Illumination
In order to further improve the efficiency of the cell exposed to a nonuniform irradiance profile, a nonuniform grid spacing is investigated to ascertain if unevenly spaced fingers enhance cell efficiency. To reduce computational complexity, the cell is virtually divided into a number of segments of equal area. Spacing is uniform within a segment but differs between segments. For example, consider a solar cell virtually divided into ten segments (say to ). Since the segments closest to the center of the cell will receive higher illumination, the fingers should be more narrowly spaced at the center of the cell than at the edges. Also, since Gaussian irradiance profiles have been used to represent nonuniform illumination, a symmetric spacing distribution can be used. For example, the spacing for the two outermost segments, and , is the same.
Simulations show a negligible change in efficiency (<0.01%) when the number of segments is increased beyond ten. Therefore, ten segments are used for our model. Each segment consists of a number of functional blocks given by where is the segment number for the symmetric left half of the cell, is the slope of the linear profile, and is a constant.
By varying the values of and , the number of functional blocks, and hence the finger spacing, is varied. Changing varies the spacing between segments more drastically than changing . In order to implement this in SPICE, recalculation of resistance values for each segment, and for intermediate areas between segment transitions, is required. Maximum slope () and constant () values are set to 10 and maximum efficiencies are obtained as the values of and are varied between 1 and 10, yielding optimal configurations for all irradiance profiles and concentrations. Figure 6 shows results for a Gaussian irradiance profile of PAR = 5.3 at 1000 suns and 2500 suns. At 1000 suns, a maximum efficiency of 36.9% is obtained at and . In this case the spacing varies from 195 μm at the cell edges to 72 μm at the cell center. At 2500 suns, = 36.4% is obtained at and . In this case the spacing varies from 137 μm at the cell edges to 48 μm at the cell center.
(a)
(b)
Figure 7 shows the optimal spacing profile as a function of segment position under nonuniform irradiance profile with a PAR = 5.3 at concentrations of 500 and 2000 suns. It is evident that the spacings become narrower toward the center of the cell. For 500 suns, the spacing varies from 328 μm at the cell edges to 86 μm at the cell center; and for 2000 suns, it varies from 162 μm at the cell edges to 51 μm at the cell center.
Similar simulations have been carried out for illumination profiles with PAR = 1.6 and PAR = 2.6. Negligible change in maximum efficiency was seen for the nearoptimal profile (PAR = 1.6). Table 2 shows the maximum efficiency for Gaussian profiles having a PAR of 2.6 and 5.3 at various concentrations.

For lower concentrations and lower PAR values, the penalty for using uniform spacing is minimal. A similar result has been reported in [20] at 1000 suns. However, at higher concentrations, a 0.3% absolute efficiency gain is observed. The effect of using a nonuniform spacing is more pronounced as the irradiance PAR increases.
4. Conclusion
Our results indicate that, under uniform illumination, the finger spacing needs to be optimized for the system concentration. Under nonuniform illumination, the effect of finger spacing on solar cell efficiency can be pronounced. For uniform spacing at 2500 suns, an efficiency enhancement of 2% absolute is obtained when optimizing using a Gaussian profile of PAR = 5.3 rather than a flat profile. For nonuniform spacing an additional enhancement in absolute efficiency of 0.3% was found. At lower concentrations, the penalty for using evenly spaced fingers becomes negligible. The need to optimize grid designs based on nonuniform illumination is thus clearly demonstrated by our optimization of front contact finger spacing on multijunction solar cells using more realistic irradiance profiles across a range of concentrations.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The authors would like to thank the National Science and Engineering Research Council of Canada, the Canadian Foundation for Innovation, Canada Research Chair Program, Ontario Research Fund, and Cyrium Technologies Inc. for supporting this work and CMC Microsystems for providing access to Sentaurus software.
References
 F. Dimroth, M. Grave, P. Beutel et al., “Wafer bonded fourjunction GaInP/GaAs//GaInAsP/GaInAs concentrator solar cells with 44.7% efficiency,” Progress in Photovoltaics: Research and Applications, vol. 22, pp. 277–282, 2014. View at: Publisher Site  Google Scholar
 C. Algora, “The importance of the very high concentration in thirdgeneration solar cells,” in Next Generation Photovoltaics, pp. 108–136, Institute of Physics Publishing, London, UK, 2003. View at: Google Scholar
 L. Fu, R. Leutz, and H. P. Annen, “Secondary optics for Fresnel lens solar concentrators,” in 7th Nonimaging Optics: Efficient Design for Illumination and Solar Concentration, vol. 7785 of Proceedings of SPIE, San Diego, Calif, USA, August 2010. View at: Publisher Site  Google Scholar
 R. Herrero, M. Victoria, C. Domínguez, S. Askins, I. Antón, and G. Sala, “Concentration photovoltaic optical system irradiance distribution measurements and its effect on multijunction solar cells,” Progress in Photovoltaics: Research and Applications, vol. 20, no. 4, pp. 423–430, 2012. View at: Publisher Site  Google Scholar
 C. Algora and V. Diaz, “Influence of series resistance on guidelines for manufacture of concentrator ponn GaAs solar cells,” Progress in Photovoltaics: Research and Applications, vol. 8, pp. 211–225, 2000. View at: Google Scholar
 H. Baig, K. C. Heasman, and T. K. Mallick, “Nonuniform illumination in concentrating solar cells,” Renewable and Sustainable Energy Reviews, vol. 16, no. 8, pp. 5890–5909, 2012. View at: Publisher Site  Google Scholar
 H. Baig, N. Sarmah, K. C. Heasman, and T. K. Mallick, “Numerical modelling and experimental validation of a low concentrating photovoltaic system,” Solar Energy Materials and Solar Cells, vol. 113, pp. 201–219, 2013. View at: Publisher Site  Google Scholar
 J. Nelson, The Physics of Solar Cells, Imperial College Press, London, UK, 2003.
 A. W. Walker, Bandgap engineering of multijunction solar cells using nanostructures for enhanced performance under concentration [Ph.D. dissertation], Dept. of Physics, University of Ottawa, Ottawa, Canada, 2013.
 A. W. Walker, O. Thériault, M. M. Wilkins, J. F. Wheeldon, and K. Hinzer, “Tunneljunctionlimited multijunction solar cell performance over concentration,” IEEE Journal of Selected Topics in Quantum Electronics, vol. 19, no. 5, pp. 1–8, 2013. View at: Publisher Site  Google Scholar
 J. F. Wheeldon, C. E. Valdivia, A. W. Walker et al., “Performance comparison of AlGaAs, GaAs and InGaP tunnel junctions for concentrated multijunction solar cells,” Progress in Photovoltaics: Research and Applications, vol. 19, no. 4, pp. 442–452, 2011. View at: Publisher Site  Google Scholar
 G. Segev, G. Mittelman, and A. Kribus, “Equivalent circuit models for triplejunction concentrator solar cells,” Solar Energy Materials and Solar Cells, vol. 98, pp. 57–65, 2012. View at: Publisher Site  Google Scholar
 P. Sharma, A. Walker, J. Wheeldon, H. Schriemer, and K. Hinzer, “Optimization of finger spacing for concentrator photovoltaic cells under nonuniform illumination using SPICE,” in Photonics North Conference, vol. 8915 of Proceedings of SPIE, Ottawa, Canada, June 2013. View at: Publisher Site  Google Scholar
 K. Nishioka, T. Takamoto, T. Agui, M. Kaneiwa, Y. Uraoka, and T. Fuyuki, “Evaluation of InGaP/InGaAs/Ge triplejunction solar cell under concentrated light by simulation program with integrated circuit emphasis,” Japanese Journal of Applied Physics, vol. 43, no. 3, part 1, pp. 882–889, 2004, 10.1143/JJAP.43.882. View at: Google Scholar
 A. Zekry and A. Y. AlMazroo, “A distributed SPICEmodel of a solar cell,” IEEE Transactions on Electron Devices, vol. 43, no. 5, pp. 691–700, 1996. View at: Publisher Site  Google Scholar
 B. Galiana, C. Algora, and I. ReyStolle, “Comparison of 1D and 3D analysis of the front contact influence on GaAs concentrator solar cell performance,” Solar Energy Materials and Solar Cells, vol. 90, no. 16, pp. 2589–2604, 2006. View at: Publisher Site  Google Scholar
 P. Rodrigo, E. F. Fernández, F. Almonacid, and P. J. PérezHigueras, “Models for the electrical characterization of high concentration photovoltaic cells and modules: a review,” Renewable and Sustainable Energy Reviews, vol. 26, pp. 752–760, 2013. View at: Publisher Site  Google Scholar
 K. Nishioka, T. Takamoto, T. Agui, M. Kaneiwa, Y. Uraoka, and T. Fuyuki, “Evaluation of InGaP/InGaAs/Ge triplejunction solar cell and optimization of solar cell's structure focusing on series resistance for highefficiency concentrator photovoltaic systems,” Solar Energy Materials and Solar Cells, vol. 90, no. 9, pp. 1308–1321, 2006. View at: Publisher Site  Google Scholar
 R. King, “Ultrahigh efficiency multijunction cell and receiver module, Phase 1B: high performance PV exploring and accelerating ultimate pathways,” Tech. Rep. NREL/SR52047602, National Renewable Energy Laboratory, Boulder, Colo, USA, 2010. View at: Google Scholar
 I. Garcia, C. Algora, I. ReyStolle, and B. Galiana, “Study of nonuniform light profiles on high concentration IIIV solar cells using quasi3D distributed models,” in Proceedings of the 33rd IEEE Photovoltaic Specialists Conference (PVSC '08), pp. 1–6, San Diego, Calif, USA, May 2008. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2014 Pratibha Sharma et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.