**American Journal of Mechanical Engineering**

## Optimization of the Geometric Design of Silicon Solar Cells under Concentrated Sunlight

**Hoe-Gil Lee**^{1,}, **Singiresu S. Rao**^{2}

^{1}Department of Mechanical Engineering Shiv Nadar University NH-91, Teshsil Dadri, District Gautam Buddha Nagar U.P. 2101 314, India

^{2}Department of Mechanical and Aerospace Engineering University of Miami, Coral Gables, FL 33146, USA

Abstract | |

1. | Introduction |

2. | Theoretical Model |

3. | Deterministic Optimization of Solar Cell Performance |

4. | Conclusion |

Nomenclature | |

References | |

Appendix |

### Abstract

The optimization of silicon solar cells is presented. The desired performance characteristics of a solar cell include a high conversion efficiency and a large power output through optimal grid contact design with minimal power losses under concentrated sunlight. Both square and rectangular shapes of linear grid contact patterns are considered for silicon solar cells with constraints on design parameters. The optical and ohmic losses are noted in computing the conversion efficiency and power output, which influence the structure and size of the solar cell and the grid contact design. In the optimization of solar cell performance, the various power losses caused by grid contact that influence the solar cell size, geometry of the fingers, busbars, and spacing between the metal grid lines are also considered. The proposed methodology permits the design of the solar cell for optimal cell performance even when some of the parameters, such as the size of the cell or the geometry of the fingers, are prespecified.

**Keywords:** solar cell, optimization, geometric design parameters, genetic algorithm

**Copyright**© 2016 Science and Education Publishing. All Rights Reserved.

### Cite this article:

- Hoe-Gil Lee, Singiresu S. Rao. Optimization of the Geometric Design of Silicon Solar Cells under Concentrated Sunlight.
*American Journal of Mechanical Engineering*. Vol. 4, No. 2, 2016, pp 50-59. https://pubs.sciepub.com/ajme/4/2/2

- Lee, Hoe-Gil, and Singiresu S. Rao. "Optimization of the Geometric Design of Silicon Solar Cells under Concentrated Sunlight."
*American Journal of Mechanical Engineering*4.2 (2016): 50-59.

- Lee, H. , & Rao, S. S. (2016). Optimization of the Geometric Design of Silicon Solar Cells under Concentrated Sunlight.
*American Journal of Mechanical Engineering*,*4*(2), 50-59.

- Lee, Hoe-Gil, and Singiresu S. Rao. "Optimization of the Geometric Design of Silicon Solar Cells under Concentrated Sunlight."
*American Journal of Mechanical Engineering*4, no. 2 (2016): 50-59.

Import into BibTeX | Import into EndNote | Import into RefMan | Import into RefWorks |

### At a glance: Figures

### 1. Introduction

The purpose of a photovoltaic concentrator solar cell is to collect as much solar energy as possible from the sun’s rays. Two steps are involved in the design of a solar cell. The first one is to optimize the performance of the solar cell by adjusting its thickness using the given data, especially the doping level. The second step is to minimize the power losses for maximizing the conversion efficiency of the solar cell under concentrated sunlight. Therefore, in order to achieve the maximum conversion efficiency, it is necessary to optimize the structure of the solar cell as well as the collecting grid contact design under concentrated sunlight. Arturo (1985) ^{[1]} described a method for the optimization of the concentration factor in terms of the size and nominal efficiency (at intensity of 1sun), by assuming practical values of the specific resistance between the grid contact patterns and the semiconductor. However, factors such as the properties of materials, the metallic value of geometric grid contact factors, and the interactions among these factors were not considered; only the relationship between the conversion efficiency and concentrated sunlight (*C*) to the length of cell is taken into account. Arturo (1985) ^{[1]} did not indicate a procedure to optimize solar cells. Gessert and Coutts (1992) ^{[2]} reviewed the models and techniques utilized to design and optimize metal contacts and antireflective coatings and indicated the differences between grid metallization of cells used under electrical resistivity by using a computer program. A limitation of Gessert’s study is that it did not examine the design in terms of constraints on design variables. Liu, Li, Chen, Wang, and Yang (2010) ^{[3]} showed the influence of metal grid lines and power losses under concentrated sunlight (*C*) in the optimization of grid contact design of a solar cell by using computer simulations. Unfortunately, the study failed to establish how to obtain individual optimal design values when geometric grid contact values are fixed as constrained parameters based on the variation of metal grid properties. In addition, the thickness of the solar cell was not factored into power output, even though these thicknesses are related to sheet resistivity due to the doping level. Kulushich, Zapf-Gottwick, Baxer-Bahci, and Werner (2013) ^{[4]} presented a method to optimize the front geometric parameters with considerations of power losses. In this work, the solar cell structure and concentration ratio were not considered. In addition, although the optimization was performed by adjusting the values of the geometric grid contact parameters through a trial and error process, the resulting value of maximum conversion efficiency is expected to be lower compared to the value obtainable through the application of mathematical programming techniques.

The aim of this work is to demonstrate how to determine grid contact design parameters and the thickness of a solar cell for any specified set of cell structure parameters. In addition, consideration of the interactions among the grid contact design parameters, cell size and structural parameters under concentrated sunlight, and constrained design parameters is emphasized. This work seeks to optimize the design of solar cells under concentrated sunlight to find the best solar cell structure and grid contact design parameters.

### 2. Theoretical Model

The performance of a solar cell can be described using the fundamental equations of semiconductor devices. There are three parts in a solar cell – the structure, cell size, and grid contact design for the optimization of solar cell performance.

The total current density in the three regions can be computed using Eq. (1) with details indicated in Appendix A.

(1) |

The maximum operating power density from a solar cell can be determined with the short-circuit current density and open-circuit voltage. The details of the computational procedure for finding the short-circuit current density and open-circuit voltage are presented in the Appendix A.

The maximum operating power density at one sun intensity can be found as

(2) |

For a sunlight concentration with intensity *C*, the equations for , and can be obtained as

(3) |

where and can be expressed by

(4) |

(5) |

where the intensity of sunlight* C* can vary in the range of 1 to 100 suns and and the ideality factor (*n*), which is chosen to lie* *between 1 and 2 for simplicity, is a measure of how closely the diode follows the ideal diode equation. When a load is connected to the diode, a current will flow in the circuit as shown in Figure 1.

Djeffal, Bendib, Arar, and Dibi (2012) ^{[5]} described an optimum grid metal design with the combined effect of the four loss mechanisms associated with the front metal grid contact design.

**Fig**

**ure**

**1**. Simple solar cell structure with grid lines and top view of contact grid structure

After calculating the theoretical performance of a solar cell, it is essential to find the solar power without considering any of power losses from the blockage of concentrated sunlight caused by the influence of geometric parameters crucial to the extraction of voltage and current characteristics. In this work, shadowing and grid contact losses are considered in the optimization of solar cell design. Thus, the total power loss can be expressed in terms of the individual fractional power losses (Appendix A):

(6) |

Arturo (1985) calculated the conversion efficiency of a solar cell under concentrated sunlight (*C*) as

(7) |

where is the incident power density at 1 sun intensity and is equal to 1 Kw/m^{2}. In this work, the solar cell performance is optimized with consideration of fractional power losses and concentrated sunlight.

### 3. Deterministic Optimization of Solar Cell Performance

The performance of a solar cell can be expressed in terms of two characteristics, namely, the conversion efficiency *,* which does not depend on the cell area, and the power output , which depends on the cell area. The goal of optimization is to obtain the best performance of the system under given restrictions. The MATLAB program, *ga*, is used for the optimization of solar cell performance. The program, *ga, *is based on genetic algorithm (GA) and, most likely, can find the global minimum of a function involving mixed-integer variables starting from a set of initial design vectors.

**3.1. Formulation of Optimization Problem**

An optimization problem can be stated in mathematical form as

(8) |

to minimize or maximize the objective function subject to the constraints

(9) |

where and are the inequality and equality constraints, respectively, is the design variable, and and are the lower and upper bounds on the design variable, respectively.

**3.2. Maximization of Solar Cell Performance**

The objective is to find the optimal vector for maximization of conversion efficiency and power output through minimization of the power losses under a solar intensity factor of *C* sums. The conversion efficiency can be expressed as

(10) |

To remove the dependence of power output of the solar cell on its area, it is more common to express the short-circuit current density as in mA/cm^{2}. Thus, the conversion efficiency of the solar cell will be related to short-circuit current , open-circuit voltage , incident power density at 1 sun, concentrated sunlight of intensity of *C* suns, and the total fractional power loss . The design vector of the problem, for a rectangular solar cell, is:

(11) |

The optimization problem is solved by placing lower and upper bounds on the design variables as , with the bounds indicated in Table 1:

The constraints of the optimization problem include relationships between the heights of finger and busbar by considering the delivery to the busbars and the shading from the busbars, the ratio of width to height of finger, and the spacing (D) between the fingers and the busbars:

(12) |

(13) |

(14) |

(15) |

(16) |

(17) |

The power output .. can be calculated from the maximum operating power densit , which corresponds to the maximum operating voltage and current density including the total fractional power loss , and the size of the solar cell.

(18) |

The problem of maximization of conversion efficiency is investigated for two cases – one by maximizing the theoretical conversion efficiency of the solar cell (with no front contact material using only 8 design variables) and the other by maximizing the practical conversion efficiency of the solar cell (with front contact material which causes ohmic and optical losses). In the case of maximization of power output, the problem is addressed by including a constraint that the conversion efficiency be at least a specified percentage of the maximum conversion efficiency and is investigated for different percentages of the maximum conversion efficiencies.

**3.3. Solution of Optimization Problems: Genetic Algorithm (GA)**

Deb (2001) ^{[6]} and Rao (2009) ^{[7]} showed that the optimization problems stated above are mixed integer programming problems and are solved using Genetic Algorithm (GA). For the GA, the population size is chosen as 20 and the initial population is generated randomly within the specified lower and upper bounds of the design variables using uniform distribution for the random numbers. The fitness of the individuals (design vectors) is based on their ranks rather than their scores. The rank of an individual is its position in the scores. In other words, the rank of the fittest individual is one, that of the next fittest is two, and so on. This type of fitness scaling removes the effect of the spread of the raw scores. The selection of parents is based on uniform random number generation with each parent corresponding to the selection of the line of length proportional to its expected fitness value. The algorithm moves along the line in steps of equal size, one step for each parent. At each step, the algorithm allocates a parent from the section it lands on. The reproduction process specifies the number of individuals that are guaranteed to survive to the next generation as two. The crossover process is implemented by specifying the fraction of the next generation that crossover produces as 0.8. Mutation produces the remaining individuals in the next generation. The mutation function makes small random changes in the individuals in the population, which provide genetic diversity and enables the genetic algorithm to search in a broader space. The crossover process selects the genes where the vector is a 1 from the first parent, and the genes where the vector is a 0 from the second parent, and combines the genes to form a child. The penalty parameter for violating a constraint is taken as 10 initially and gradually increased to 100 depending on the observed rate of convergence. For the convergence of the algorithm, a function tolerance of 1e-20 and a constraint tolerance of 1e-6 (for nonlinear constraints) are used.

**3.4. Maximization of Conversion Efficiency of the Solar Cell**

The maximization of the conversion efficiency of a solar cell can be considered in two stages - theoretical and practical conversion efficiency. Figure 2 shows the variations of short-circuit current density in the emitter, base, space-charge regions (SCR) and the total current density with the thickness of the emitter over the range 0-8μm (by solving a number of optimization problems by fixing the thickness of the emitter at one specific value at a time). The thickness of the base is fixed as 250μm. The current density in the emitter is found to increase steeply from 0 to 2 μm with a slow variation beyond a value of 2μm while the current density in the base is decreased within the same range. This means that the thickness of the emitter between 2 μm and 8 μm corresponds to a small variation of the total current density between 39.88 mA and 40.01 mA. As a result, the variations in the thicknesses of the emitter and base will have a large influence on the total current density in specific ranges between 0.1 μm and 2 μm for the emitter of a solar cell. On the other hand, the range between 2 μm and 8 μm for the emitter a solar cell will have much less impact on the total current density.

To maximize the practical conversion efficiency of a solar cell, the grid variations of the contact design variables, namely, the length of the cell, the width, height and number of the fingers, and the width, height and number of the busbars are to be considered. These variables lead to optical and ohmic losses that cause a reduction in the theoretical conversion efficiency. Two shapes- square and rectangular- are considered for the maximization of the practical efficiency of the solar cell. In the case of the square solar cell configuration, the width and height of the cell will be same; as such only ten design variables are considered by eliminating in the design vector of Eq. (11). A parametric study is conducted to find the influence of the concentrated sunlight on the open-circuit voltage of the solar cell. For this, the optimization problem stated in section 3.2 is solved several times by using only the first nine variables in the design vector, Eq. (11), with the value of the concentrated sunlight fixed at a different value each time.

**Fig**

**ure**

**2**

**.**Relationships between short-circuit current density and thicknesses of emitter

**3.4.1. Theoretical Conversion Efficiency of the Solar Cell Under a Concentration Ratio of One**

The efficiency of a solar cell depends on a range of factors, such as materials, dopant concentrations ( and ), thickness of a solar cell ( and ), recombination velocity, and the variations associated with the fabrication process. An efficient solar cell consists of a thin emitter formed by low energy, low recombination velocity at the front surface, and enough base thickness for proper absorption from a wide range of wavelengths. However, if the thickness is thin enough, there is a limited recombination velocity of minority carriers toward the p-n junction. The emitter allows a solar cell to generate power. Heavily doped materials are the most important factors to determine the cell performance, because those minority-carrier diffusion coefficients ( and ), minority-carrier lifetime, and minority diffusion lengths ( and ) influence the open-circuit voltage and short-circuit current . Once the design factors of a solar cell’s performance are determined, there are a number of potential methods that could be utilized based on a study of the performance characteristics at various concentrations of sunlight. The optimal design factor of cell thickness reduces material costs because it results in a thin solar cell, which is competitively priced to have significant impact on the size of large-scale power systems. If a solar cell has a smaller thickness, material is saved. At the same time, its performance is maximized. The objective function for the maximization of the theoretical conversion efficiency of a solar cell can be expressed as

(19) |

All the constraints stated in section 3.2, except those associated with the grid contact materials, are considered in the solution process. Note that the first three design variables in Eq. (11), which relate to the grid contact materials, are excluded from the design vector. The theoretical conversion efficiency is determined to find the characteristics of the solar cell using the values of structural parameters and pre-specified values of the cell thickness and concentrated sunlight. The thickness of base is assumed as 250 μm, same as the value recommended in the literature under a concentrated sunlight of 1 sun intensity. The input design factors used in this study are given by : dopant concentration of emitter =2×10^{17}cm^{-3}, dopant concentration of base = 5×10^{1}^{6}cm^{-3}, recombination velocity of front surface = 1 ×10^{2}cm/s , and recombination velocity of back surface = 1 ×10^{4}cm/s. The charateristics of the solar cell at optimal design, including the optimal values of the thicknesses of the emitter and base and the optimal conversion efficiency, are shown in Table 2. Next, a parametric study is conducted to find the influence of emitter thickness on the short-circuit current density.

**3.4.2. Maximization of Practical Conversion Efficiency Considering All Ten Design Variables in Eq. (11)**

The optimization problem stated in section 3.2 is solved by placing lower and upper bounds on all the design variables. In reference to the work of Sharan (1987) ^{[8]}, Gessert (1992) ^{[2]} and Djeffal (2012) ^{[5]}, the design data are assumed as = 2 ×10^{17}cm^{-3}, = 5 ×10^{1}^{6}cm^{-3} (dopant concentrations of emitter and base), = 1 *cm/s** *, = 1 ×10^{2}*cm/s* (recombination velocities) to determine the thickness of the solar cell by specifying grid contact resistance, sheet resistance, and metal resistivity for silicon as = 3 ×10^{-3}Ω cm^{2}, = 100 Ω/cm^{2}, and = 1.6 ×10^{-6}Ω·*cm*, respectively. The solution of the optimization problem yields the optimal values of the solar cell thickness (i.e., emitter and base thicknesses), grid contact design parameters, and the intensity of concentrated sunlight.

If the cell size is permitted to vary along with other design variables during optimization, the optimal values of the remaining design factors are expected to be different compared to the values found in the previous cases. This study is conducted in this section. The optimal values of the design variables, the value of the geometric grid lines, short-circuit current, open-circuit voltage, and the optimal theoretical and practical efficiencies of the solar cell can be obtained as shown in Table 3.

As can be seen from Table 3, the length of the cell and width of the finger approached their respective lower bound values. The maximum practical efficiencies of the square and rectangular cells with minimal power losses are found to be 20.28 % and 20.54 %, respectively. The corresponding theoretical efficiencies are found to be 22.80 % in both types of cells.

**3.4.3. Optimization with Pre-specified Cell Size and Width of Fingers**

The practical conversion efficiency of the solar cell is maximized by specifying the area of the cell as 11cm^{2} for a square cell and 32cm^{2} for a rectangular cell, and the width of fingers as 80 μm. Thus there will be only eight design variables left in Eq. (11) in the optimization problem. The other data are same as the ones given in Table 2. The results obtained from the solution of the optimization problem (using the program, *ga*) are shown in Table 4.

It is noticed that the specification of the geometric parameters of the cell area and the width of a finger did not affect the thickness of the solar cell, short-circuit current and open-circuit voltage. The practical conversion efficiency is decreased from 20.28 % to 18.17 % and 20.54 % to 19.17 %, while the fractional power losses are increased from 11.09 % to 18.88 % and 9.93 % to 13.32 % for the square and rectangular cells, respectively, due to the optical and resistive losses. This shows that the change in the conversion efficiency is not proportional to the change in fractional power losses due to concentrated sunlight and interactions among the various geometric parameters.

**3.5. Maximization of Power Output**

The conversion efficiency of the solar cell is independent of the power output. It is only a factor when the solar cell is designed for maximum power output within a specified area. Thus, in order to find the maximum power output of the solar cell, a minimum constraint of 80 % of the maximum conversion efficiencies (in both square and rectangular cells) becomes necessary to prevent unrealistic conversion efficiencies resulting from the intensity of sunlight. Table 5 shows the values of design variables and other outputs corresponding to maximum power output with a constraint of at least 80 % of the maximum conversion efficiency in square and rectangular cells.

In the case of a square cell, the maximum power output is 0.4770 W and the conversion efficiency is 19.08 %, which is 94 % of the maximum conversion efficiency. Since the power output is related the size of the solar cell and conversion efficiency, which is associated with an increase in individual power losses, the total power loss in a square cell for maximum power output is 9.27 %. The dominant power loss of 5.39 % is the fractional power loss due to shadowing from busbars and fingers (sunlight blockage). Thus, the maximum power output is given by the product of the power density, 0.0191 and the area of the solar cell, 25cm^{2} as 0.4770 W, for a square cell.

In the case of a rectangular cell, the maximum power output is found to be 0.4770 W with the conversion efficiency of 19.08 %, nearing 93 % of the maximum conversion efficiency (20.54 %). The difference in total power loss between maximum conversion efficiency and maximum power output is 1.93 %, an indication that the total power loss in maximum power output is higher as a result of differing sizes of the solar cell. The main fractional power loss (5.56%) is caused by shadowing from fingers and busbars because the number of fingers and busbars has increased from 12 to 41 and from 3 to 7, respectively. In addition, the widths of the finger and busbar associated with the power losses have increased from 20.03 μm to 37.64 μm and from 100.21 μm to 188.21 μm, respectively. With consideration of cell structure, cell size, grid contact design, and conversion efficiency, the power output of 0.4770 W can be found as the product of density, 0.0191 , and the area of the solar cell, 25 cm^{2} for a rectangular cell. This indicates that the relationship between the conversion efficiency and power output plays a role in the optimization of solar cell performance.

The single-objective optimization problems for maximum conversion efficiency provide a solar cell design having a maximum conversion efficiency of 20.28 % and 20.54 % for square and rectangular cells, respectively. In this section, the maximum power output of a solar cell is investigated through sensitivity analysis with respect to applying different constraint values on the maximum conversion efficiency (formulated as the ratio) in the range of 70 % to 100 %.

The difference between the lowest and highest power output of a solar cell is mainly caused by changes in power density, collected amount of sunlight, and total fractional power loss. An increase in the thicknesses of a solar cell emitter and base has the influence of improving the conversion efficiency . The widths of fingers and busbars contribute to conversion efficiency and power density since they prevent collection of the proper amount of sunlight. Beyond a value of 90 % for minimum permissible value of the maximum conversion efficiency, the cell areas have steeply decreased from 25.00 cm^{2} to 0.66 cm^{2} for a square cell and from 25.00 cm^{2} to 1.23 cm^{2} for a rectangular cell. This indicates that an increase in the conversion efficiency combined with an increase in the power output occur due to changes in the geometric design and the number of fingers and busbars; these are also associated with total power loss and power density.

Table 6 shows the variations of conversion efficiency, total power loss, power density, cell area, and power output with respect to different minimum permissible values of constraint on conversion efficiency in finding the maximum power output. The relationship between the conversion efficiency and power output has been observed through variations in power density and cell area associated with geometric design parameters, which indicated that the amount of power output has increased even though the power density has decreased by the reduction of conversion efficiency. Thus, the area of a solar cell contributes to an increase in the total power output in both types of cells. This indicates that the behaviors of conversion efficiency and power output are approximately opposite; hence, a compromise solution is to be found in a practical solar cell design.

#### Table 6. Variations of conversion efficiency, power density, cell area, and power output with respect to maximum conversion efficiency ratio

### 4. Conclusion

The optical and ohmic losses are the major causes that reduce the efficiency of a solar cell in the practical conversion of solar energy to electricity. To maximize the conversion efficiency and power output from sunlight to electricity, the parameters of the solar cell structure and the mechanism of collecting solar energy are determined using mathematical programming approaches. Solar cell modeling has revealed that, in optimizing solar cell performance, a correlation exists between the cell structure and geometric design parameters for maximizing conversion efficiency of the solar cell. The solar cell size, power density, and total power loss influence the maximizing power output. Heavy doping levels reduce the open-circuit voltage and short-circuit current of a solar cell due to the emitter layer with resistive properties, while close spacing between fingers results in high power losses. This is the most important reason why optimal geometric design factors should be found for a solar cell for maximum conversion efficiency. The design of solar cells under specified constraints on design parameters, finds the optimal intensity values of *C* suns, the values of the solar cell structure, such as the thicknesses of emitter and base and grid contact design parameters when the conversion efficiency is maximized. Even when the areas of square and rectangular solar cells are restricted to lie in any specific ranges, such as, 0.5 × 0.5 cm^{2} to 5 × 5 cm^{2} and 0.5 to 5 × 0.5 to 5 cm^{2}, respectively, the values of the optimal geometric design parameter can be determined through optimization. It is found that, in some cases, higher total power losses do not always correspond to maximum conversion efficiency on account of concentrated sunlight. Therefore, the correlation of cell structure parameters, grid contact design, cell size, conversion efficiency, and power output under intensity of sunlight should be deliberated by using optimum design procedure for the solar cell.

### Nomenclature

= Intensity of sunlight (an integer)

= Spacing between the busbars (cm)

= Spacing between the fingers (cm)

= Minority electron diffusion coefficient (cm^{2}/s)

= Minority hole diffusion coefficient (cm^{2}/s)

=Fractional power loss of the resistivity of the busbars (%)

= Fractional power loss of contact resistance (%)

= Fractional power loss of the resistivity of the fingers (%)

= Total fractional power losses (%)

= Fractional power loss of shadowing (%)

= Fractional power loss of sheet resistance (%)

= Height of cell (cm)

= Current density of base (mA/ cm^{2})

= Current density of emitter (mA/ cm^{2})

= Light-generated current density (mA/ cm^{2})

= Maximum operating current density (mA/ cm^{2})

= Current density of space-charge region (mA/ cm^{2})

= Saturation current density (mA/ cm^{2})

= Short-circuit density (mA/ cm^{2})

= Boltzmann constant (8.617×10^{-5}eV/K)

= Minority electron diffusion length

= Minority hole diffusion length

= Current transfer length

= Length of cell (cm)

= Acceptor concentration (cm^{3})

= Donor concentration (cm^{3})

= Intrinsic carrier concentration (cm^{3})

= Photon flux

= Maximum operating power density (W/ cm^{2})

= Input incident power density (W/ m^{2})

= Power output (W)

= Electron charge (1.602×)

= Reflection coefficient of the anti-reflective coating

=Sheet resistance (Ω/ cm^{2})

= Recombination velocity of the front surface (cm/s)

= Recombination velocity of the back surface (cm/s)

= Width of the space charge region

= Temperature (K)

= Thickness of the emitter region

= Thickness of the base region

= Open-circuit voltage (mV)

= Maximum operating voltage (mV)

= Width of cell (cm)

= Absorption coefficient (cm^{-1})

= Conversion efficiency of a solar cell (%)

= Maximum conversion efficiency of a solar cell (%)

= Minority carrier lifetime in the emitter region

= Minority carrier lifetime in the base region

= Shrinkage of the energy gap (eV)

= Contact resistance (Ω·cm^{2})

= Metal resistivity (Ω·cm)

= Wavelength

### References

[1] | Arturo, M.A, 1985, Optimum concentration factor for silicon solar cells, Solar Cells, Vol. 14, No. 1, pp. 43-49. | ||

In article | View Article | ||

[2] | Gessert, T.A and Coutts, T.J., 1992, Grid metallization and antireflection coating optimization for concentrator and one sun photovoltaic solar cells, Science & Technology of Materials, Vol. 10, No. 4, pp. 2013-2024. | ||

In article | View Article | ||

[3] | Liu, W., Li, Y., Chen, J., Chen, Y., Wang, X., and Yang, F., 2010, Optimization of grid design for solar cells, Journal of Semiconductors, Vol.31, No. 1, pp. 0140061-4. | ||

In article | |||

[4] | Kulushich, G.,Zapf-Gottwick, R., Baxer-Bahci,B, and Werner, J.H., 2013, 18.4% Efficient Grid Optimized Cells With 100 Ω/sq Emitter, IEEE, Vol. 3, No. 1, pp. 254-260. | ||

In article | |||

[5] | Djeffal, F., Bendib, T., Arar, D., Dibi, Z., 2012, An optimized metal grid design to improve the solar cell performance under solar concentration using multiobjective computation, Materials Science and Engineering, Vol. 178, No. 9, pp. 574-579. | ||

In article | View Article | ||

[6] | Deb, K., 2001, Multi-Objective Optimization Using Evolutionary Algorithms, Wiley, Chichester. | ||

In article | |||

[7] | Rao, Singiresu S., 2009, Engineering Optimization Theory and Practice, 4th Ed., Wiley, Hoboken. | ||

In article | |||

[8] | Sharan, S. N., Mathur, S. S., Kandpal, T. C., 1987, Optimum concentration ratio for an actively cooled solar concentrator photovoltaic system with a fin-type absorber, Energy Convers. Mgmt, Vol. 27 No. 4, pp 351-354. | ||

In article | View Article | ||

[9] | Jain, S.C, Heasell, E.L and Roulston, D.J., 1987, Recent advances in the physics of silicon P-N junction solar cells including their transient response, Prog. Quant. Electr., Vol.11, No. 2, pp. 105-204. | ||

In article | View Article | ||

[10] | Singal, C.M, 1981, Optimum cell size for concentrated-sunlight silicon solar cells, Solar Cells, Vol. 3, No. 1, pp. 9-16. | ||

In article | View Article | ||

[11] | Liou, J.J. and Wong, W.W., 1992, Comparison and optimization of the performance of Si and GaAs solar cells, Solar Energy Materials and Solar Cells, Vol. 28, No. 1, pp. 9-28. | ||

In article | View Article | ||

[12] | Shabana, M. M., Saleh M.B., Soliman M.M., 1988, Optimization of grid design for solar cells at different illumination levels, Solar cells, Vol. 26, No. 3, pp.177-187. | ||

In article | View Article | ||

[13] | Geeves, G. K., Harrison, H. B., 1982, Obtaining the specific contact resistance from transmission line model measurements, IEEE, Vol. 3, No. 5, pp. 111-113. | ||

In article | |||

### Appendix

**Appendix A: Determination of the maximum operating power density (P**

_{m}) and total fraction power loss (F_{sum})Calculations of the short-circuit current density and open circuit voltage for the maximum operating power and the total fractional power losses from ohmic resistance losses are the essential factors to optimize a solar cell based on the geometric design parameters.

**• Calculation of short-circuit current density (**

*J*_{sc}) and open circuit voltage (*V*_{oc})The total current density can be expressed as

(A-1) |

Where the individual value of , , and can be refereed by the publication of Jain, Heasell, and Roulston (1986) ^{[9]} and Singal (1980) ^{[10]}.

Total current density can be calculated as

(A-2) |

(A-3) |

(A-4) |

(A-5) |

The photon flux density and Si absorption coefficient were described by Liou and Wong (1992) ^{[11]} under AM1.5 global normal sun condition can be approximated with two linear curves as

(A-6) |

(A-7) |

Also, the absorption coefficient of silicon material can be divided into 4 sections to obtain the data based on particular wavelength range, and is given by

(A-8) |

The short-circuit current density is due to the generation and collection of light-generated carriers and can be expressed as

(A-9) |

where the reverse saturation current density can be computed using the following equations, this is also known as the diode equations

(A-10) |

(A-11) |

Also, the open-circuit voltage can be found as

(A-12) |

Thus, the maximum operating power at one sun can be found as

(A-13) |

where the maximum current and voltage are given by

(A-14) |

**• Calculation of total fractional power loss (F**

_{sum})The calculation of total fractional power loss was explained by Shabana, Saleh, and Soliman (1988) ^{[12]} and can be expressed as

(A-15) |

where n is the number of individual terms of fractional power loss. While calculating the total fractional power loss , the weighted fractional losses for the individual terms are calculated, with the weights being the power contributions from the respective terms.

The total fractional power loss can also be expressed in terms of the individual fractional power losses as

(A-16) |

The resistance of the sheet can be expressed in a differential form as

(A-17) |

Thus the power loss due to sheet resistance can be calculated as

(A-18) |

The power generated can be expressed as

(A-19) |

The fractional power loss of sheet resistance is given by

(A-20) |

Under normal circumstances, the contact resistivity can be written, using the concept of transfer length , as

(A-21) |

The specific contact resistance was described by Harrison and Reeves (1982) ^{[13]} and the power loss of contact resistivity can be found as

(A-22) |

Thus, the fractional contact loss is given by

(A-23) |

The top of a solar cell has a series of arranged fingers intended to collect current. The corresponding resistive loss is given by

(A-24) |

Because of symmetry, the equation is applied precisely at the midway along the length of finger to obtain

(A-25) |

Thus, the fractional power loss of a finger can be expressed as

(A-26) |

The ratio of width to thickness of a contact should be within the limits of the recommended aspect ratio, which is 0.23~0.25. Also, the fractional power loss of busbars is given by

(A-27) |

The fractional power loss of shadowing depends on the size and number of grid lines ( and ) because it prevents light from entering a solar cell. The fractional power loss of shadowing is given by

(A-28) |