International Journal of PhotoenergyVolume 2012 (2012), Article ID 369101, 8 pagesdoi:10.1155/2012/369101

Research Article

Modelling of Light Trapping in Acidic-Textured Multicrystalline Silicon Wafers

Yang Li, Zhongtian Li, Yuebin Zhao, and Alison Lennon

School of Photovoltaic and Renewable Energy Engineering, The University of New South Wales, Sydney, NSW 2052, Australia

Received 31 October 2011; Accepted 6 November 2011

Academic Editor: Donghwan Kim

Copyright © 2012 Yang Li 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.

Abstract.

Acidic texturing has been widely used to reduce the reflection losses for silicon solar cells fabricated on multicrystalline wafers, however, there are few available models which attempt to predict the reduced reflection after texturing based on the morphology of the textured surfaces. An optical model which simulates the light trapping and scattering effects of acidic-textured surfaces based on the surface morphology is presented. The developed model was experimentally verified by reflection measurements from multicrystalline silicon wafers textured using different etching conditions. The relationship between weighted average reflection and surface morphology is demonstrated with some of the trends being explained by simulating reflection in different wavelength regions. The developed model could be embedded into solar cell simulation tools or adapted to predict optical properties of diverse surface morphologies.

1. Introduction

Surface texturing is a key process in the fabrication of silicon solar cells both in the laboratory and in commercial production. It provides enhanced light absorption with less front surface reflection, more internal reflection, and increased light path lengths inside the cell [1, 2]. A wide range of physical and chemical methods have been employed for texturing of silicon wafer-based and thin-film solar cells, including lithography [3, 4], mechanical grooving [5, 6], reactive ion texturing [7–9], laser texturing [10–12], alkaline [13–15], and acidic texturing [16–20].

Of these texturing methods, alkaline and acidic texturing have been widely adopted in the commercial production of silicon wafer-based solar cells due to their lower cost, amenability to high-throughput processing, and superior antireflection properties. Alkaline etchants etch the planes of silicon much faster than the planes. This anisotropic etching, when performed at low sodium or potassium hydroxide concentrations and temperatures of 80–90°C, results in the formation of pyramidal structures on orientated monocrystalline silicon. These pyramidal structures trap the incoming light and reduce reflection effectively [2]. However, alkaline texturing performed on multicrystalline silicon results in largely untextured regions (where crystal planes are exposed at the surface) and steps between grains due to the different grain orientations of the multicrystalline silicon [15, 21, 22]. Although acidic texturing of multicrystalline wafers results in higher average reflectance values than typically observed for alkaline texturing of monocrystalline wafers, the acidic texturing process is relatively low-cost and can be effectively performed as an in-line process in a manufacturing line. In fact, because of the ease of implementing acidic texturing compared to alkaline texturing, some silicon solar cell manufacturers select to also texture monocrystalline wafers using acidic texturing, because the difference in surface reflection between alkaline- and acidic-textured surfaces is reduced significantly after application of a silicon nitride antireflection coating and encapsulation.

Acidic texturing has also been employed to texture ZnO : Al transparent conducting oxide (TCO) layers, which serve as window layers in silicon-based thin-film solar cells [23]. It may also find application in the texturing of wafers produced using seeded casting processes [24]. The wafers that result from these new casting processes still contain grains but have large regions of monocrystalline material. A number of different methods are being developed for large-scale manufacture of this new silicon material [24–26], and because it is likely that grain orientation cannot always be carefully controlled, acidic texturing may also find a role in texturing these new wafer surfaces.

Simulation tools, such as PC1D [27], contain models that simulate the reflection and light trapping for alkaline-textured solar cells and predict the electrical performance of textured solar cell devices. These tools enable higher level solar cell production line simulators, such as the Virtual Production Line (VPL), which has been developed at the University of New South Wales (UNSW) [28, 29], to relate the effects of different alkaline texturing conditions (e.g., chemical composition, time, and temperature) on final device electrical performance. The VPL simulation software is extensively used in the teaching of photovoltaic engineering at UNSW. However, the current version of this software can only simulate the reflectance of alkaline-textured surfaces, which limits its application to current manufacturing processes which involve the use of multicrystalline silicon wafers.

Although the chemical process for acidic texturing can be simulated [29], a model is required to predict the reflectance as a function of wavelength and hence the electrical performance of the textured solar cell device. The light-trapping properties of the acidic-textured surface morphology have not been extensively studied. Earlier studies have considered only the front-surface reflection properties of these surfaces [21, 22]. However, this paper describes a more comprehensive model which considers the reflection from the entire solar power spectrum, bulk absorption, and the contribution of the rear surface to light trapping. These enhancements to the model become critical as thinner wafers are increasingly used in the manufacturing of cells. Acidic texturing experiments were performed and the measured reflectance results were compared to those predicted by the model. The results suggest that the model can predict the reflection for a wide range of wavelengths. Further work is required to integrate this reflection model in solar cell simulators such as PC1D or higher level simulators such as VPL.

2. Model of Acidic-Textured Surfaces

Multicrystalline silicon surfaces textured using the UKN recipe for acidic texturing [18–20] are typified by an irregular array of concavities, or etching pits, as shown in Figure 1. These pits in the textured surface are generated by the etching of surface defects in the crystal structure caused by sawing damage. Unlike in alkaline texturing where an initial saw-damage etching process must be performed before the texturing, in acidic texturing the saw damage is removed during the texturing thus simplifying the manufacturing process. However, the weighted average reflectance (WAR) values achieved by alkaline texturing are significantly lower than those achieved with acidic texturing. Furthermore, if acidic texturing times are used that are longer than optimal, then the textured concavities become larger and flatter and the resulting surface is significantly more reflective [30].

Figure 1: A scanning electron microscope (SEM) image of an acidic-textured multicrystalline silicon wafer surface showing the formation of randomly arranged etching pits on the surface. The surface was textured for 2 mins at 10°C in a solution comprising an HF : HNO₃ : H₂O ratio of 15 : 50 : 35.

To simplify the simulation of the reflection based on surface morphology described in this paper, a spherical cap surface was used to model the etching pit structure of the acidic-textured surface. The front-surface reflection, bulk absorption, rear surface reflection, and further internal reflections of light reaching each part of the etching pit were estimated, with the reflection of the wafer being estimated as the weighted average of all components.

2.1. Front Surface Scattering

The spherical model is shown schematically in Figure 2. The variables and denote the centre and radius of the sphere cap, and and donate the diameter and height of the sphere cap, respectively. The sphere cap is divided into small annular surface elements, , each element having a subtended angle to . The kth-order reflection at the mth annular surface element is denoted using .

Figure 2: Schematic depicting a cross-section of a spherical cap, and the associated variables, used to model the etch pits of an acidic-textured surface.

It is assumed that incident light is directed perpendicularly to the silicon surface, and hence also perpendicular to the base of the sphere cap. Considering the incident light reaching annular surface element with an incident angle of and a corresponding refracted angle of , the primary reflection at the mth annular surface element is given by where and represent the reflection of the transverse-electric and transverse-magnetic waves of the incident light, respectively. Light rays having a large incident angle can be reflected more than once as shown in Figure 2. The reflected ray intersects the sphere again at annular surface element with the incident angle, and refracted angle , which are equal to and , respectively, before being reflected away from the surface. Thus, the front surface reflection for the mth annular surface element is given by where donates the highest reflection order and is given by The expression in (3) is derived from internal angle sum of the polygon, which consists of the intersection points and as its vertices.

2.2. Bulk Absorption

Considering the same incident light ray that reached discussed above, the attenuation factor after being absorbed inside the wafer is given by where denotes the actual light path length, is the thickness of the wafer, is the absorption coefficient, and is the refracted angle of light reaching annular surface element referring to the normal direction to the wafer, as shown in Figure 3. This expression shows that the light path is increased due to the textured surface and hence the probability of absorption by the silicon wafer is increased for the light’s first pass through the wafer.

Figure 3: Schematic showing the path of a representative transmitted light ray, scattered within the wafer.

2.3. Rear Surface Reflection

The rear surface is modelled using an inverted sphere cap as shown in Figure 4. The incident light reaching on the front surface is scattered with an angle of and has the possibility of reaching the convex sphere surface at the rear of the textured wafer. To evaluate the angle distribution function (ADF) of the rear incident angles, the convex surface of the rear etch pit is divided into annular surface elements perpendicular to the incident light so, for example, light rays reaching different points on pth annular surface element would have the same incident angle of . The ADF is determined by the projection area of the annular surface elements to the incident direction as represented by where is the projected area of unshaded part of pth rear annular surface element, is the shading percentage of it induced by the adjacent etching pit derived from Figure 4, is the radius of the sphere cap on the rear surface, and is a geometric factor determining the fraction of annular surface element of the sphere cap exposed in the incident direction. The rear surface reflection of the light, which entered via the mth annular front surface element, is the average of the light reflected from the convex surface of the weighted by the ADF as represented by Based on the primary reflection and absorption calculated above and the assumption that these values remain the same during internal reflection, the overall reflection of incident light reaching at the front surface of a free standing acidic-textured bare (i.e., uncoated) silicon wafer can be expressed as Figure 5 depicts the dependence of on . The value of for wavelengths shorter than 1000 nm does not change with increasing values of due to high absorption coefficient of this light and hence the high probability that this light is absorbed before it reaches the rear surface. However, for longer wavelengths of incident light, first increases and then decreases with increasing , with a very abrupt decrease being predicted for wavelengths of ≥1200 nm. This increase in as increases for longer wavelength light is due to the light, which is reflected from the rear surface, ultimately escaping from the wafer before absorption. This reflected loss of initially captured light is more likely for longer wavelength light, which has lower absorption coefficient and is more likely to encounter internal surfaces during its light path through the wafer. It is also enhanced for surfaces which are characterised by lower values where light is coupled less obliquely into the wafer.

Figure 4: Schematic showing parallel light rays incident on the rear acidic-textured surface.

Figure 5: Reflection, , as a function of rear surface reflection, , for different wavelengths of incident light. The values of and wafer thickness used in simulation were 0.25 and 200 μm respectively.

The overall reflection of the etching pit, , is the average reflection of all annular surface elements weighted by their projection area, as given by where and are the projected areas of the small annular surface element and sphere cap, respectively, and is the radius of the sphere cap on the front surface, which is not necessarily the same as due to processes like rear etching. Figure 6 shows the reflection calculated from (8) as a function of wavelength and ratio. The calculations use light absorption coefficients and refractive indices at 300 K provided in [31].

Figure 6: Reflection predicted by the model as a function of wavelength and surface morphology indictor . The wafer thickness used in the simulations was 200 μm.

3. Experimental

Texturing experiments were performed to examine the accuracy of the developed model. In order to show the relationship between reflection and surface morphology, multicrystalline silicon wafers of thickness 200 μm were textured in solutions comprising HF, HNO₃, and deionised water to achieve a range of etching pit scales by changing the etching conditions (e.g., chemical composition of the solution and reaction temperature).

The dimensions of the spherical cap (i.e., and ) were estimated from SEM cross-sectional images taken from different locations through the wafer. For each wafer surface, values of and were estimated as the average of ~20 measurement points taken from five different SEM cross-sectional images. The reflection from the front surface of wafers, which had been textured on both sides, was measured for a wavelength range of 250 nm to 1300 nm using a Perkin Elmer spectrophotometer. It was estimated that <1% of the light that was transmitted through the wafer was reflected back into the wafer. Weighted average reflection values were calculated by weighting the measured reflection against the intensity of direct radiation in the AM1.5G spectrum.

4. Results and Discussion

Figure 7 shows SEM images of some of the different wafer surface morphologies that were observed after acidic texturing. It demonstrates how different etch pit geometries can be achieved using different etching conditions, with higher HF : HNO₃ ratios and lower texturing temperatures resulting in deeper pits which are characterised by larger values. Figure 8 shows a comparison between simulated and experimental reflection of two textured wafers which had different values over the wavelength range of 250 nm to 1300 nm. The simulated reflection data predicts the main trends in the experimental reflection curves, however, some inaccuracy in the simulated data is to be expected due to the simplified assumptions of the model (e.g., assumptions that the surface is characterised by a single value). Furthermore, mismatches between regions selected for reflection measurement and SEM imaging may contribute to discrepancies between simulated and experimental reflection values. Values of greater than 0.5 were not considered in this analysis because these surfaces involve more coupling of surrounding pits and possible porous silicon formation than surface morphologies that are typically observed with UKN texturing.

(a)

(b)

(c)

(d)

Figure 7: Scanning electron microscope images of the surface and cross sections of multicrystalline silicon wafer surfaces textured under two different conditions: (a,b) HF : HNO₃ : H₂O = 15 : 50 : 35, 25°C, 2 min; (c,d) HF : HNO₃ : H₂O = 24 : 40 : 36, 10°C, 2 min.

Figure 8: Simulated and measured reflection values graphed as a function of wavelength for two samples with different values. Reflection values for are shown in red with the solid line and (red triangle) symbols denoting the simulated and measured values, respectively. Reflection values for are shown in black with the solid line and () symbols denoting the simulated and measured values, respectively.

The relationship between the WAR and was also investigated, due to the usefulness of WAR values in comparing different texturing methods. Figure 9 shows how the WAR varies with values. The WAR does not vary significantly with increasing until an value of ~0.2, after which it begins to decrease and subsequently plateaus at a value of 0.24. The decrease in WAR, which occurs at ~0.2 in Figure 9, is due to the onset of higher-order reflection (i.e., ). The critical value of can be calculated from (3) to be 0.2071, after which the reflection is reduced due to the reduced front surface reflection.

Figure 9: Simulated (solid line) and measured (◆) WAR values, computed using the wavelength range of 300 nm to 1200 nm, as a function of .

The contribution of the rear surface reflection to the increased reflection predicted by the model and depicted in Figure 5 suggests that, at least for longer wavelengths when is small and light has a greater probability of reaching the rear surface because of the shorter light path in the wafer, some increase in reflection is expected. By graphing the simulated WAR of the short (550–600 nm) and long wavelengths (1150–1200 nm) of light separately, as shown in Figure 10, this predicted increase in reflection is demonstrated for only the longer wavelength light.

Figure 10: Simulated and measured WAR values for short-wavelength region (550 nm to 600 nm) and long-wavelength region (1150 nm to1200 nm) graphed as a function of . The black solid line and (◆) values represent the simulated and measured values for the short wavelength region and the red solid line and (red triangle) symbols represent the simulated and measured values for the long-wavelength region.

However, although the model predicts the WAR reasonably well for the short wavelength light where the reflection is largely determined by the front surface, it appears to only successfully predict the reflection of long wavelength light for values greater than ~0.25. For the cases where is higher than the critical value of high-order reflection, the WAR of both short- and long-wavelength light decreases as increases because the front surface reflection is effectively reduced by higher-order reflections within etching pits. However, as becomes even larger, higher-order reflection only occurs on a smaller surface of the etching pit, thus the value of the WAR plateaus as the value of approaches 0.5. For the region where is lower than the critical value, the WAR of short wavelength light is substantially unchanged because there is no significant higher-order reflection arising from the front surface. The short wavelength light has a high probability of absorption before reaching the rear surface without requiring a lengthened light path due to front surface scattering.

Increased reflection for long wavelength light at low values is not evident in the experimental data. This is most likely due to higher-order reflection not being considered in the rear surface reflection model described in Section 2.3. This can result in underestimation of the light path length in the wafer and hence higher predicted WAR values for low values. However, the model appears to accurately predict the WAR for larger values, because the longer light paths caused by the increased scattering increases the probability that the light is absorbed in the bulk and therefore reduces the sensitivity of the model to the rear surface reflection as shown in Figure 5. For small values, light is coupled into the wafer at less oblique angles and hence experiences shorter light paths in the wafer and is less likely to be internally reflected.

5. Conclusions

An optical model was presented for the modelling of light trapping in acidic-textured silicon wafers. Sphere caps were used to model the etching pits that form during the texturing when the commonly-used UKN texturing process is used. The developed model was used to predict the reflection from bare acidic-textured multicrystalline wafers on the basis of the surface morphology which was varied by altering the texturing conditions. The model was shown to predict the reflection for a wide range of wavelengths. The model predicted that, for shorter wavelength light, the WAR was constant for small values and then began to decrease at where higher-order reflection effects became significant resulting in increased scattering of light in the wafer. For longer wavelength light, where the absorption coefficient is significantly reduced, the model predicted a small increase in WAR for small values due to the shorter paths experienced by light that is coupled less obliquely into the wafer. This results in a greater chance of initially captured light escaping from the wafer on subsequent encounters with wafer surfaces. However, this increase in WAR for low values was not observed experimentally. It is proposed that this was due to the fact that the rear surface reflection model used did not take into account the effects of higher-order reflection at the convex surface of the etching pits and thus, most likely, underestimated the length of the internally reflected light path.

The inaccuracies of the rear surface reflection model may not be critical for practical purposes because the rear surfaces of many acidic-textured silicon solar cells are partially planarised by a rear surface etching process used to isolate the junction to the front surface of the cell and then alloyed with aluminium to form a rear-surface field and rear-cell electrode. Consequently, the rear surface reflection model, developed here for bare wafers, will need to be further developed to more accurately represent the light trapping properties of final devices.

Finally, the model has the potential to be adapted for use in solar cell simulators, such as PC1D, or higher-level simulators such as the VPL where it can be used to relate acidic texturing processing conditions to final device efficiency and hence enhance the educational and training experiences for future PV engineers. The accuracy of the model could be enhanced by an improved characterisation of the surface morphology of the textured surfaces which are formed using different processing parameters, and further characterisation of the reflectance from final devices.

Acknowledgment

The authors would like to thank Dr. Deming Zhu of UNSW Electron Microscopy Unit for assistance with scanning electron microscopy.

References

E. Yablonovitch and G. D. Cody, “Intensity enhancement in textured optical sheets for solar cells,” IEEE Transactions on Electron Devices, vol. 29, no. 2, pp. 300–305, 1982.

P. Campbell and M. A. Green, “Light trapping properties of pyramidally textured surfaces,” Journal of Applied Physics, vol. 62, no. 1, pp. 243–249, 1987.

A. Wang, J. Zhao, and M. Green, “24-percent efficient silicon solar cells,” Applied Physics Letters, vol. 57, pp. 602–604, 1990.

J. Zhao, A. Wang, and P. Campbell, “A 19.8% efficient honeycomb multicrystalline silicon solar cell with improved light trapping,” IEEE Transactions on Electron Devices, vol. 46, no. 10, pp. 1978–1983, 1999.

G. Willeke, H. Nussbaumer, H. Bender, and E. Bucher, “A simple and effective light trapping technique for polycrystalline silicon solar cells,” Solar Energy Materials and Solar Cells, vol. 26, no. 4, pp. 345–356, 1992.

C. Gerhards, C. Marckmann, R. Tolle, et al., “Mechanically V-textured low cost multicrystalline silicon solar cells with a novel printing metallization,” in Proceedings of the of the 26th IEEE Photovoltaic Specialists Conference, pp. 43–46, 1997.

Y. Inomata, K. Fukui, and K. Shirasawa, “Surface texturing of large area multicrystalline silicon solar cells using reactive ion etching method,” Solar Energy Materials and Solar Cells, vol. 48, no. 1–4, pp. 237–242, 1997.

W. A. Nositschka, C. Beneking, O. Voigt, and H. Kurz, “Texturisation of multicrystalline silicon wafers for solar cells by reactive ion etching through colloidal masks,” Solar Energy Materials and Solar Cells, vol. 76, no. 2, pp. 155–166, 2003.

D. S. Ruby, S. H. Zaidi, S. Narayanan, B. M. Damiani, and A. Rohatgi, “Rie-texturing of multicrystalline silicon solar cells,” Solar Energy Materials and Solar Cells, vol. 74, no. 1–4, pp. 133–137, 2002.

J. C. Zolper, S. Narayanan, S. R. Wenham, and M. A. Green, “16.7% efficient, laser textured, buried contact polycrystalline silicon solar cell,” Applied Physics Letters, vol. 55, no. 22, pp. 2363–2365, 1989.

M. Abbott and J. Cotter, “Optical and electrical properties of laser texturing for high-efficiency solar cells,” Progress in Photovoltaics, vol. 14, no. 3, pp. 225–235, 2006.

L. A. Dobrzański, A. Drygała, K. Gołombek, P. Panek, E. Bielańska, and P. Zieba, “Laser surface treatment of multicrystalline silicon for enhancing optical properties,” Journal of Materials Processing Technology, vol. 201, no. 1–3, pp. 291–296, 2008.

H. Seidel, L. Csepregi, A. Heuberger, and H. Baumgaertel, “Anisotropic etching of crystalline silicon in alkaline solutions. I. Orientation dependence and behavior of passivation layers,” Journal of the Electrochemical Society, vol. 137, no. 11, pp. 3612–3626, 1990.

Z. Xi, D. Yang, W. Dan, C. Jun, X. Li, and D. Que, “Investigation of texturization for monocrystalline silicon solar cells with different kinds of alkaline,” Renewable Energy, vol. 29, no. 13, pp. 2101–2107, 2004.

J. D. Hylton, A. R. Burgers, and W. C. Sinke, “Alkaline etching for reflectance reduction in multicrystalline silicon solar cells,” Journal of the Electrochemical Society, vol. 151, no. 6, pp. G408–G427, 2004.

R. Einhaus, E. Vazsonyi, J. Szlufcik, J. Nijs, and R. Mertens, “Isotropic texturing of multicrystalline silicon wafers with acidic texturing solutions,” in Proceedings of the IEEE 26th Photovoltaic Specialists Conference, pp. 167–170, October 1997.

K. Kim, S. K. Dhungel, S. Jung, D. Mangalaraj, and J. Yi, “Texturing of large area multi-crystalline silicon wafers through different chemical approaches for solar cell fabrication,” Solar Energy Materials and Solar Cells, vol. 92, no. 8, pp. 960–968, 2008.

A. Hauser, I. Melnyk, E. Wefringhaus, F. Delahaye, G. Vilsmeier, and P. Fath, “Acidic texturisation of mc-Si using a high throughput in-line prototype system with no th organic chemistry,” in Proceedings of the 19th European Photovoltaic Solar Energy Conference, pp. 1094–1097, 2004.

S. Mathijssen, S. Braun, I. Melnyk, et al., “Survey of acid texturing and new innovative acid processes for mc solar wafers,” in Proceedings of the 24th European Photovoltaic Solar Energy Conference, pp. 1964–1967, 2009.

A. Hauser, I. Melnyk, P. Fath, S. Narayanan, S. Roberts, and T. M. Bruton, “A simplified process for isotropic texturing of MC-SI,” in Proceddings of the 3rd World Conference on Photovoltaic Energy Conversion, vol. 2, pp. 1447–1450, May 2003.

Z. Xi, D. Yang, W. Dan, C. Jun, X. Li, and D. Que, “Texturization of cast multicrystalline silicon for solar cells,” Semiconductor Science and Technology, vol. 19, no. 3, pp. 485–489, 2004.

Y. Nishimoto, T. Ishihara, and K. Namba, “Investigation of acidic texturization for multicrystalline silicon solar cells,” Journal of the Electrochemical Society, vol. 146, no. 2, pp. 457–461, 1999.

J. I. Owen, J. Hüpkes, H. Zhu, E. Bunte, and S. E. Pust, “Novel etch process to tune crater size on magnetron sputtered ZnO:Al,” Physica Status Solidi A, vol. 208, no. 1, pp. 109–113, 2011.

H. Zhang, L. Zheng, X. Ma, B. Zhao, C. Wang, and F. Xu, “Nucleation and bulk growth control for high efficiency silicon ingot casting,” Journal of Crystal Growth, vol. >318, no. 1, pp. 283–287, 2011.

K. Fujiwara, W. Pan, N. Usami et al., “Growth of structure-controlled polycrystalline silicon ingots for solar cells by casting,” Acta Materialia, vol. 54, no. 12, pp. 3191–3197, 2006.

N. Stoddard, B. Wu, I. Witting et al., “Casting single crystal silicon: novel defect profiles from BP solar's mono2 ™ wafers,” Solid State Phenomena, vol. 131–133, pp. 1–8, 2008.

D. A. Clugston and P. A. Basore, “PC1D version 5: 32-bit solar cell modeling on personal computers,” in Proceedings of the IEEE 26th Photovoltaic Specialists Conference, pp. 207–210, October 1997.

S. R. Wenham and A. Bruce, “Virtual production line for the manufacturing of screen-printed solar cells,” in Proceedings of the 19th European Photovoltaic Solar Energy Conference, 2004.

Z. Li, Y. Li, M. J. Fu, and A. Lennon, “Simulation of acidic texturing for the virtual production line software,” in Proceedings of the 21st International Photovoltaic Science and Engineering Conference, 2011.

Y. T. Cheng, J. J. Ho, W. J. Lee et al., “Investigation of low-cost surface processing techniques for large-size multicrystalline silicon solar cells,” International Journal of Photoenergy, vol. 2010, Article ID 268035, 6 pages, 2010.

M. A. Green and M. J. Keevers, “Optical properties of intrinsic silicon at 300 K,” Progress in Photovoltaics, vol. 3, no. 3, pp. 189–192, 1995.