Modeling of photo-conversion efficiency for hydrogenated amorphous Si p-i-n structures

An analytical formalism to optimize the photoconversion efficiency η of hydrogenated amorphous silicon-based (a-Si:H) solar cells has been developed. This model allows firstly the optimization of a p⁺ -i-n sandwich in terms of carrier mobilities, thickness of the layers, doping levels, and other...

Full description

Saved in:
Bibliographic Details
Date:2007
Main Authors: Sachenko, A.V., Sokolovskyi, I.O., Kazakevitch, A., Shkrebtii, A.I., Gaspari, F.
Format: Article
Language:English
Published: Інститут фізики напівпровідників імені В.Є. Лашкарьова НАН України 2007
Series:Semiconductor Physics Quantum Electronics & Optoelectronics
Online Access:https://nasplib.isofts.kiev.ua/handle/123456789/118334
Tags: Add Tag
No Tags, Be the first to tag this record!
Journal Title:Digital Library of Periodicals of National Academy of Sciences of Ukraine
Cite this:Modeling of photo-conversion efficiency for hydrogenated amorphous Si p-i-n structures / A.V. Sachenko, I. O Sokolovskyi, A. Kazakevitch, A.I. Shkrebtii, F. Gaspari // Semiconductor Physics Quantum Electronics & Optoelectronics. — 2007. — Т. 10, № 4. — С. 60-66. — Бібліогр.: 13 назв. — англ.

Institution

Digital Library of Periodicals of National Academy of Sciences of Ukraine
id nasplib_isofts_kiev_ua-123456789-118334
record_format dspace
spelling nasplib_isofts_kiev_ua-123456789-1183342025-06-03T16:25:24Z Modeling of photo-conversion efficiency for hydrogenated amorphous Si p-i-n structures Sachenko, A.V. Sokolovskyi, I.O. Kazakevitch, A. Shkrebtii, A.I. Gaspari, F. An analytical formalism to optimize the photoconversion efficiency η of hydrogenated amorphous silicon-based (a-Si:H) solar cells has been developed. This model allows firstly the optimization of a p⁺ -i-n sandwich in terms of carrier mobilities, thickness of the layers, doping levels, and others. Second, the geometry of grid fingers that conduct the photocurrent to the bus bars and ITO/SiO₂ layers has been optimized, and the effect of non-zero incidence angles of Sun’s light has been included as well. The optimization method has been applied to typical a-Si:H solar cells. The codes allow the optimization of amorphous Si based solar cells in a wide range of parameters and are available on the e-mail request. 2007 Article Modeling of photo-conversion efficiency for hydrogenated amorphous Si p-i-n structures / A.V. Sachenko, I. O Sokolovskyi, A. Kazakevitch, A.I. Shkrebtii, F. Gaspari // Semiconductor Physics Quantum Electronics & Optoelectronics. — 2007. — Т. 10, № 4. — С. 60-66. — Бібліогр.: 13 назв. — англ. 1560-8034 PACS 72.20.Jv, 73.40.Cg, 84.60.Jt, 85.30.De https://nasplib.isofts.kiev.ua/handle/123456789/118334 en Semiconductor Physics Quantum Electronics & Optoelectronics application/pdf Інститут фізики напівпровідників імені В.Є. Лашкарьова НАН України
institution Digital Library of Periodicals of National Academy of Sciences of Ukraine
collection DSpace DC
language English
description An analytical formalism to optimize the photoconversion efficiency η of hydrogenated amorphous silicon-based (a-Si:H) solar cells has been developed. This model allows firstly the optimization of a p⁺ -i-n sandwich in terms of carrier mobilities, thickness of the layers, doping levels, and others. Second, the geometry of grid fingers that conduct the photocurrent to the bus bars and ITO/SiO₂ layers has been optimized, and the effect of non-zero incidence angles of Sun’s light has been included as well. The optimization method has been applied to typical a-Si:H solar cells. The codes allow the optimization of amorphous Si based solar cells in a wide range of parameters and are available on the e-mail request.
format Article
author Sachenko, A.V.
Sokolovskyi, I.O.
Kazakevitch, A.
Shkrebtii, A.I.
Gaspari, F.
spellingShingle Sachenko, A.V.
Sokolovskyi, I.O.
Kazakevitch, A.
Shkrebtii, A.I.
Gaspari, F.
Modeling of photo-conversion efficiency for hydrogenated amorphous Si p-i-n structures
Semiconductor Physics Quantum Electronics & Optoelectronics
author_facet Sachenko, A.V.
Sokolovskyi, I.O.
Kazakevitch, A.
Shkrebtii, A.I.
Gaspari, F.
author_sort Sachenko, A.V.
title Modeling of photo-conversion efficiency for hydrogenated amorphous Si p-i-n structures
title_short Modeling of photo-conversion efficiency for hydrogenated amorphous Si p-i-n structures
title_full Modeling of photo-conversion efficiency for hydrogenated amorphous Si p-i-n structures
title_fullStr Modeling of photo-conversion efficiency for hydrogenated amorphous Si p-i-n structures
title_full_unstemmed Modeling of photo-conversion efficiency for hydrogenated amorphous Si p-i-n structures
title_sort modeling of photo-conversion efficiency for hydrogenated amorphous si p-i-n structures
publisher Інститут фізики напівпровідників імені В.Є. Лашкарьова НАН України
publishDate 2007
url https://nasplib.isofts.kiev.ua/handle/123456789/118334
citation_txt Modeling of photo-conversion efficiency for hydrogenated amorphous Si p-i-n structures / A.V. Sachenko, I. O Sokolovskyi, A. Kazakevitch, A.I. Shkrebtii, F. Gaspari // Semiconductor Physics Quantum Electronics & Optoelectronics. — 2007. — Т. 10, № 4. — С. 60-66. — Бібліогр.: 13 назв. — англ.
series Semiconductor Physics Quantum Electronics & Optoelectronics
work_keys_str_mv AT sachenkoav modelingofphotoconversionefficiencyforhydrogenatedamorphoussipinstructures
AT sokolovskyiio modelingofphotoconversionefficiencyforhydrogenatedamorphoussipinstructures
AT kazakevitcha modelingofphotoconversionefficiencyforhydrogenatedamorphoussipinstructures
AT shkrebtiiai modelingofphotoconversionefficiencyforhydrogenatedamorphoussipinstructures
AT gasparif modelingofphotoconversionefficiencyforhydrogenatedamorphoussipinstructures
first_indexed 2025-11-24T14:44:09Z
last_indexed 2025-11-24T14:44:09Z
_version_ 1849683292289236992
fulltext Semiconductor Physics, Quantum Electronics & Optoelectronics, 2007. V. 10, N 4. P. 60-66. © 2007, V. Lashkaryov Institute of Semiconductor Physics, National Academy of Sciences of Ukraine 60 PACS 72.20.Jv, 73.40.Cg, 84.60.Jt, 85.30.De Modeling of photoconversion efficiency for hydrogenated amorphous Si p-i-n structures A.V. Sachenko1∗, I.O. Sokolovskyi1, A. Kazakevitch2, A.I. Shkrebtii2, F. Gaspari2 1V. Lashkaryov Institute of Semiconductor Physics, NAS of Ukraine, 03028 Kyiv, Ukraine 2University of Ontario Institute of Technology, Oshawa, ON, L1H 7L7, Canada ∗Corresponding author, e-mail: sachenko@inbox.ru Abstract. An analytical formalism to optimize the photoconversion efficiency η of hydrogenated amorphous silicon-based (a-Si:H) solar cells has been developed. This model allows firstly the optimization of a p+-i-n sandwich in terms of carrier mobilities, thickness of the layers, doping levels, and others. Second, the geometry of grid fingers that conduct the photocurrent to the bus bars and ITO/SiO2 layers has been optimized, and the effect of non-zero incidence angles of Sun’s light has been included as well. The optimization method has been applied to typical a-Si:H solar cells. The codes allow the optimization of amorphous Si based solar cells in a wide range of parameters and are available on the e-mail request. Keywords: photoconversion efficiency, hydrogenated amorphous silicon, a-Si:H solar cell. Manuscript received 26.11.07; accepted for publication 19.12.07; published online 31.01.08. 1. Introduction Thin film hydrogenated amorphous silicon (a-Si:H) is widely used for photovoltaic applications [1]. Amor- phous silicon-based solar cells (SC) are very promising because of low production cost, possibility of covering large uneven areas, and sufficiently high efficiency. In order to get the best possible performance of the a-Si:H solar cells, it is important to (i) produce high-quality amorphous films with p-i-n junction and (ii) optimize the films and solar cells by their parameters such as, for instance, p-, i- and n-layer thicknesses, their doping le- vels, electron and hole mobilities µn and µp, lifetimes of electrons and holes, resistance of p-, i- and n-layers, contact grid geometry, and parameters of transparent conducting and antireflecting layers. Many experimental techniques to grow a-Si:H films have been proposed recently. Many of them require a reasonably high substrate temperature to grow high-quality samples. One of the authors (F.G.) has developed the so-called dc saddle-field glow-discharge technique (see [2] and references therein). This efficient technique was used to deposit a hydrogenated thin Si film on the amorphous solar cells under consideration. It was previously demonstrated that such films can be successfully used for photovoltaic applications [3]. Regarding the theoretical models, a few approaches have been proposed to optimize the performance of solar cells of various types [1]. They range from the analytical models [4] and combined analytical and numerical approaches [4, 5] to purely numerical models using a specially developed software such as AMPS. Analytical models have the advantage of being physically intuitive and predictive on a wide range of the cell parameters and offer a possibility of the quick and accurate estimation of photoconversion efficiency including the optimization of the solar cell geometry. However, the common ap- proaches to characterize analytically hydrogenated amorphous silicon solar cells were mainly based on the models developed for a crystalline semiconductor. To extend these models to amorphous Si, the crystalline Si bandgap Eg was simply substituted by the mobility gap Eµ for an amorphous material [6]. In such approaches, the main distinction between amorphous and crystalline Si-based solar cells was due to different resistances of the bases (i-layers). However, to properly characterize a-Si:H solar cells, this is not sufficient. In addition, the Stark effect must be taken into account. For instance, the series resistance of SC built on p-i-n structures of hydrogenated amorphous silicon is primarily determined by the photoresistance of the i-layer (in contrast to crystalline Si). It is commonly accepted that this important factor determines the fill factor γ of the I-V characteristic of p-i-n structures [1]. An extra parameter responsible for γ is γm. This is the ratio of the maximum power PMPP at the so-called maximum power point Semiconductor Physics, Quantum Electronics & Optoelectronics, 2007. V. 10, N 4. P. 60-66. © 2007, V. Lashkaryov Institute of Semiconductor Physics, National Academy of Sciences of Ukraine 61 (MPP) to the product of the short-circuit current ISC and the open-circuit voltage VOC [2]. Since PMPP = ImVm (Im and Vm are the current and the voltage at MPP), γm = ImVm / ISCVOC [2]. The above theoretical models produce γm for typical hydrogenated amorphous silicon-based SC in a range of 0.81 – 0.83, which is always higher than experimental values [1]. This is due to the extra photo- voltage drops on the photoresistance of the i-layer and the layer resistance of a front contact and a contact net. This has to be included in the analytical model as well. Another problem with the above standard approaches is that they are based on the so-called diode theory of rectification which works well for the crystalline Si- based SC. For a a-Si:H based solar cell, however, the diffusion theory of rectification is more suitable because of the high defect density and the absence of a long- range order, which essentially reduces the mobility of carriers and their lifetime. Two of the authors (A.V.S. and A.I.S.) have developed analytical two-dimensional models for the optimization of crystalline solar cells [7]. The one- dimensional analytical model of photoconversion efficiency for amorphous silicon-based solar cells has been proposed by one of the author (A.V.S.) [8]. Although analytical one-dimensional photoconversion models have been successfully developed in the past [8, 9], numerical methods are being mainly used recently to optimize solar cells, especially amorphous Si-based ones [10]. In this paper, we develop an analytical model of photoconversion η for hydrogenated amorphous silicon a-Si:H based p-i-n SC under the АМ0 condition. The model allows the straightforward extension to other illumination conditions, such as, for instance, AM1.5. The model considers an interplay of various factors such as p-, i-, and n-layer doping levels and conductivity, carrier mobilities µh and µp, lifetimes of carriers, diffu- sion length, contact grid geometry, etc. Finally, the solar cells might not be optimally oriented with respect to the Sun, and the additional optimization with respect to the angle of light incidence is included in the model as well. We will demonstrate how the above factors influence the fill factor γ of the I-V characteristic. We consider the mobility gap (or electric bandgap) Eµ of amorphous Si to be 1.71 eV, although the theory is suitable for various values of this gap. A spectral peculiarity of the photo- current collection that originates from the diffusion theory of rectification [7] is taken into account as well. Theoretical results were compared with available experimental data and a good agreement has been found. 2. Theory We consider a conventional SC of the sandwich type consisting of the p-i-n structure between the frontal grid and rear electrodes (see Fig. 1). The frontal collecting grid electrode consists of parallel metal fingers connected to one another through two conductive bus bars. The grid is placed on the top of a conductive p-type ITO film of thickness d2. The area between fingers is filled with an antireflecting SiO2 film with thickness of d1. The highly doped p-layer of thickness dp is in contact with the frontal electrode on the top and the i-layer of thickness d on the bottom. The highly doped n-layer of thickness dn is located between the i-layer and the uniform Al rear electrode with reflectivity Rd (Fig. 1). The frontal grid and rear electrodes collect photocurrent. ITO film (i) serves as a wideband window; (ii) provides the transparency for a solar illumination; and (iii) collects the photocurrent between metal fingers of the grid electrode. The highly doped p and n layers create a rectifying barrier that separate electrons and holes, while the electron-hole pair generation takes place mainly in the i-layer. In the case of the single reflection from the metal rear electrode, the electron-hole pair generation function g(α,x) at a given wavelength λ is described by ,)])(2exp( )([exp),( xdddR xIxg npd α+++α−⋅+ +α−⋅⋅α=α (1) where α = α(λ) is the wavelength-dependent absorption coefficient, and I is the intensity of monochromatic light at the wavelength λ. The α(λ) dependence for hydrogenated amorphous silicon is taken from our experimental measurements (previously, the earlier results [8] were used). Since the absorption coefficient essentially depends on the photon energy, we distinguish (i) the high-energy part of the spectrum (above Eµ), (ii) intermediate energy (near Eµ), and (iii) the energy below the mobility gap Eµ, in the so- called exponential bandtail. In the present paper, we take Eµ = 1.71 eV. The analytical procedure of calculation of the hole diffusion flux at the n+-layer boundary and that of electrons in the i-layer in the plane x = xn is outlined in [5]. This gives the electron collection coefficient in the p-layer fn(α)as follows: Fig. 1. Schematic view of a solar cell under consideration. It consists of the p-i-n structure between frontal grids and the rear contact. The frontal collecting grid electrode contains parallel metal fingers connected to one another through two conductive bus bars. The frontal grid is placed on the top of the transparent conducting ITO layer. Semiconductor Physics, Quantum Electronics & Optoelectronics, 2007. V. 10, N 4. P. 60-66. © 2007, V. Lashkaryov Institute of Semiconductor Physics, National Academy of Sciences of Ukraine 62 × ⎥ ⎥ ⎦ ⎤ ⎢ ⎢ ⎣ ⎡ ⎟⎟ ⎠ ⎞ ⎜⎜ ⎝ ⎛ +⎟⎟ ⎠ ⎞ ⎜⎜ ⎝ ⎛ × −α α =α −1 02 coshsinh 1)( )( n p p n pn n L d V L d S Ln Lf ( ) ( ) ( )×α−α−⎟⎟ ⎠ ⎞ ⎜⎜ ⎝ ⎛ α−− ⎪⎩ ⎪ ⎨ ⎧ ×⎟ ⎟ ⎠ ⎞ ⎜ ⎜ ⎝ ⎛ ⎟⎟ ⎠ ⎞ ⎜⎜ ⎝ ⎛ α−−+α× pn p pn p pn dL Ln d dV Ln d dSD expsinhexp coshexp10 (2) ⎪⎭ ⎪ ⎬ ⎫ ⎥ ⎥ ⎦ ⎤ ⎢ ⎢ ⎣ ⎡ ⎟⎟ ⎠ ⎞ ⎜⎜ ⎝ ⎛ +⎟⎟ ⎠ ⎞ ⎜⎜ ⎝ ⎛ n p n n p L d V L d S coshsinh0 . The hole collection coefficient in the i-layer is ( ) ( ) ( ) ( ) ( )[ ] ( ) ( )[ ] ( )} ( ) ( ) ( ) ( ) ( )[ ] ( ) . 1 2expexp expexp )(exp112 2expexp exp exp 1 )( 22 2 1 22 −α α−α−−α−α + +⎥ ⎦ ⎤ ⎟ ⎠ ⎞ ⎜ ⎝ ⎛ −⋅−−⎢ ⎣ ⎡ ⎟ ⎠ ⎞ ⎜ ⎝ ⎛⋅+× ×+α−+−−α+ +α−α−+α−× ×⎥ ⎦ ⎤ ⎟ ⎠ ⎞ ⎜ ⎝ ⎛ −⋅−+ ⎪⎩ ⎪ ⎨ ⎧ +⎢ ⎣ ⎡ ⎟ ⎠ ⎞ ⎜ ⎝ ⎛⋅+⋅ −α α −≅α − L ddRdL L dVS L dVS ddRSRD ddRd L dVS L dVS L Lf pdp dd pddd pdp d dp (3) Here, nnn DL τ= , Vn = Dn /Ln, Dn is the electron diffusion coefficient in the p+-layer, pDL τ= , V = D /L, D is the diffusion coefficient of holes in the i- layer, S0 and Sd are surface recombination velocities at the illuminated (top) and rear surfaces, respectively. As was mentioned above, the diffusion theory of rectification has to be used for hydrogenated amorphous silicon, and the total electron-hole pair collection coefficient f (α) is [7, 9] )exp( 1 )()()( p s s pn d L L fff α− α+ α −α+α≅α . (4) Here, LS = kT / qEm, and Em is the maximum electrical field in the p+-i junction. The last term in (4) is responsible for a decrease of the collection coefficient at short wavelengths. This happens even if the surface recombination velocity S0 and the bulk recombination in the p-layer can be neglected. The calculated spectral dependence of the collection coefficients for this case is shown in Fig. 2. It is clear from the plot that the stronger the maximum value of the electric field Em, the higher is the maximum of the collection coefficient f. The increase of Em shifts the region where the collection coefficient decreases toward shorter wavelengths. 0.3 0.4 0.5 0.6 0.7 0.8 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 3 2 1 C ol le ct io n co ef fic ie nt k Wavelength λ (µm) Fig. 2. Spectral dependences of the collection coefficient in a-Si:H-based SC for different values of the maximum electric field Em equal to 105, 3·104, and 104 V/m: curves 1, 2, and 3, respectively. Equation (5) gives the short-circuit current JSC for the AM0 conditions. Approximating AM0 by the black body radiation at Тс = 5800 K, taking the mobility gap to be 1.71 eV, measuring the JSC current density in А/сm2, and combining the fundamental constants, we obtain a dimensionless convenient expression for the short-circuit current JSC in the form . 1exp cos)),(1())(( )1(2)( 1 0 4 2 3 ∫ ⎥ ⎥ ⎦ ⎤ ⎢ ⎢ ⎣ ⎡ −⎟⎟ ⎠ ⎞ ⎜⎜ ⎝ ⎛ ⋅λ ⋅ ϕ⋅ϕ−⋅α × ×−⋅π⎟⎟ ⎠ ⎞ ⎜⎜ ⎝ ⎛ λ =ϕ dz zkT hcz zRzf mq D rcJ x S S S x SC (5) Here, ϕ is the sunlight incidence angle, m is the frontal grid shadowing coefficient (m ≤ 1, the bigger m is, the smaller part of the solar cell is illuminated), RS (z, ϕ) is the reflection coefficient of the SC frontal surface, z = λ /λx, λ is the wavelength of the incident light, λx is the threshold wavelength of the photoelectric effect for hydrogenated amorphous silicon, rS is the equatorial radius of the Sun, and DS is the mean Sun- Earth distance. We used λx = 0.8 µm, which corresponds to the absorption edge of an a-Si:H film for Eµ = 1.71 eV. The spectral dependence of the reflection coefficient RS (z, ϕ) for a multilayer structure SiO2 – ITO – α-Si:H was calculated according to Berning’s approach [6] which considers the oblique incidence of the sunlight at an angle φ. To maximize the short-circuit current density JSC, the reflections from the SC should be as low as possible. For the typical refractive index of SiO2 (n = 1.45) and ITO (n = 1.96), the minimal ref- lection of the SiO2-ITO structure occurs when the SiO2 film is 0.09 µm thick, and the ITO layer is 1.5 µm thick. Semiconductor Physics, Quantum Electronics & Optoelectronics, 2007. V. 10, N 4. P. 60-66. © 2007, V. Lashkaryov Institute of Semiconductor Physics, National Academy of Sciences of Ukraine 63 The angular dependence of the short circuit current density JSC normalized to its maximum value is shown in Fig. 3. The angle ϕ corresponds to Earth’s axis in- clination with respect to the normal to the ecliptic plane. The short-circuit current density at ϕ ≈ 23° is 91 % of that at ϕ = 0, while JSC is 93 % of its maximum value for the angle-independent reflection coefficient RS. It is instructive to consider the upper limit of the short circuit current density Jlim, when all photo- generated carriers in the i-layer are collected: [ ] ∫ ⎥ ⎥ ⎦ ⎤ ⎢ ⎢ ⎣ ⎡ −⎟⎟ ⎠ ⎞ ⎜⎜ ⎝ ⎛ ⋅λ ⋅ −⋅−αα−+α−− × ×−⋅π⎟ ⎠ ⎞ ⎜ ⎝ ⎛ λ ≅ 1 0 4 2 3lim . 1exp ))(1()1))(exp(2exp(exp(1 )1(2 dz zkT hcz zRddRd qmq D rcJ x Sd S x (6) Figure 4 demonstrates Jlim and JSC dependence when the i-layer thickness d changes from 0.1 to 10 µm. As seen, Jlim monotonically increases, as d grows (curve 1). When the rear contact is not reflective (Rd = 0), the density of the short-circuit current JSC increases with d and saturates after 1 µm (curve 6). Curves from 2 to 5 located between Jlim(d) at Rd = 1 and JSC(d) at Rd = 0 correspond to Rd values of 0.8, 0.6, 0.4, and 0.2. The current density reaches its maximum when dm ≈ L and increases, as the bulk recombination velocity Sd decreases. Finally, JSC(d) approaches Jlim(d) when the diffusion length L is much larger than d, and the surface recombination S0 as well as Sd can be neglected. The standard expression for the photocurrent Jph, (or solar cell I-V characteristic) in terms of the diode ideality factor A and the saturated current density JS is as follows: ⎟⎟ ⎠ ⎞ ⎜⎜ ⎝ ⎛ −= kTA qV JJVJ SSC exp)(ph . (6) 6 11 17 23 0.90 0.92 0.94 0.96 0.98 1.00 N or m al iz ed c ur re nt d en si ty J S C Angle ϕ (degrees) Fig. 3. Normalized angular dependence of the short-current density JSC for a-Si:H-based SC. Curve 1 includes the re- flection coefficient dependence on the light incidence angle φ, while curve 2 does not consider it. 0,1 1 10 8 12 16 20 24 28 C ur re nt d en si ty J S C (m A /c m 2 ) Film thickness d (µm) 1 2 3 4 5 6 Fig. 4. Photocurrent density versus the i-layer thickness d. Curve 1 corresponds to the maximally possible photocurrent density Jlim(d), when all the light is adsorbed in the film. Curves 2-6 correspond to the photocurrent JSC(d). For curve 2, the recombination is minimal, and the field is maximal. For curve 6, Rd = 0. For curve 3, t = 10−5 s, for curves 4, 5, and 6, t = 10−6 s. For curves 5 and 6, S0 = 102 cm/s, Sd = 104 cm/s. For curves 3, 4, S0 = 102 cm/s, Sd = 10 cm/s. The open-circuit voltage VOC of the SC obtained from (6) by putting Jph = 0 is ⎟⎟ ⎠ ⎞ ⎜⎜ ⎝ ⎛ = S SC OC J J q kTAV ln . (7) The expression for the photoconversion efficiency under the AM0 condition ηAM0 can be written [2] as , )( 1 / )/ln(1 / 11 ph 0 ⎟⎟ ⎠ ⎞ ⎜⎜ ⎝ ⎛ + −×⎟⎟ ⎠ ⎞ ⎜⎜ ⎝ ⎛ −× ⎟⎟ ⎠ ⎞ ⎜⎜ ⎝ ⎛ −≅η OC SCSC OC OC OC OCSC AM V RRSJ AkTqV AkTqV AkTqVP VSJ (8) where S is the SC area, Р = 0.135 W/cm2 is the power of incident solar radiation under the AM0 condition, RSC is the component of the series layer resistance of the ITO film and the grid electrode, Rph is the component of the series resistance of photogenerated carriers in the i-layer. If SJSC RSC / VOC << 1 and SJSC Rph / VOC << 1, the last term of Eq. (8) can be rewritten as γR γσ, where γR and γσ are multipliers of the fill factor responsible for the resistance of the grid electrode and the i-layer, respectively. This product demonstrates the influence of the series resistance of SC on the fill factor of the I-V characteristic. The value of γR determines the influence of the SiO2 layer resistance and the grid electrode resistance; and γσ is related to the i-layer photore- sistance. While determining the series resistance Rph of the i- region, we consider the electron mobility much over the hole mobility. Then ∫ ∆+µ = d n xnn dx qS R 0 0 ph )( 1 , (9) Semiconductor Physics, Quantum Electronics & Optoelectronics, 2007. V. 10, N 4. P. 60-66. © 2007, V. Lashkaryov Institute of Semiconductor Physics, National Academy of Sciences of Ukraine 64 where S is the sample surface area, µn the electron mobility, n0 the equilibrium electron concentration in the i-region, and ∆n(x) the increase of the electron concentration under illumination. The ∆n(x) value is determined from the relation , )1( ))2exp()(exp( expexp)()( 22 2* 21 LD xdRxLI L xC L xCxpxn d α− α+α−+α−α + +⎟ ⎠ ⎞ ⎜ ⎝ ⎛+⎟ ⎠ ⎞ ⎜ ⎝ ⎛−=∆=∆ (10) where ∆p(x) is the solution of a differential equation for holes in the i-region for an open circuit. For AM0 conditions, the I∗ value is determined from the following formula: ∫ ⎥ ⎦ ⎤ ⎢ ⎣ ⎡ −⎟ ⎠ ⎞ ⎜ ⎝ ⎛⋅ ϕ⋅ϕ− ⋅−⋅=ϕ 1 0 4 * 1 12.3 exp cos)),(1( )1( 28.1 )( dz z z zR m q I S . (11) After the determination of the integration constants C1 and C2 from the boundary conditions for the hole flows at x = dp and x = dp + d, )()( 0 pp dxpSdxj =∆−== , (12) )()( ddxpSddxj pdp +=∆−=+= , (13) and the substitution of the values obtained in Eqs. (10) and (9), one can calculate the component of the I-V curve filling factor (1 − SJSC Rph / VOC) that is related to the series resistance of the i-region. As the next step, we consider the frontal grid consisting of parallel metal fingers connected to bus bars (see Fig. 1). To determine the grid parameters, we can neglect the resistance of the metal grid electrodes compared to the essentially higher layer resistance of the ITO film. If the inequality SJSC RSC / VOC << 1 is satis- fied, γR can be written as [11-13] ⎟⎟ ⎠ ⎞ ⎜⎜ ⎝ ⎛ =γ c c R L l l L 2 tanh 2 . (11) Here, Lc = [qµ p NpVOC / JOC]1/2 is the effective collection length, µp is the hole mobility, Np = npd, and np is the two-dimensional hole concentration (Gibbs excess) in the SiO2 film. Equation (11) is accurate when Lc >> l, where l is the distance between metal fingers. If a metal finger has width ln, and m is the relative front-side metallization, that is the shadowing of the photoactive area by the metal grid (0 < m < 1), the distance between fingers l ≅ ln(1−m) / m (see [9]). Since the open-circuit voltage VOC for AM0 is very close to that at AM1.5, the efficiency ηAM1.5 of SC at AM1.5 can be expressed in terms of the efficiency ηAM0 at AM0 as 0 5.1 05.1 1.0 135.0 AM SC AM SC AMAM J J ⋅⋅η≈η . (12) To find the peak values of 5.1AM SCJ , we carried out the numerical integration of Eq. (12) at АМ1.5, considering for simplicity that the collection coefficient equals 1. For Eµ = 1.71 eV, the value of ηAM1.5 is about 10 % higher than that of ηAM0. 3. Optimization of parameters of hydrogenated amorphous silicon solar cells Using the expressions derived above, we can finally maximize the SC performance. The parameters of SC that have to be optimized include the thicknesses and doping levels of the base (i-layer) and the emitter (p- layer), mobility gap Eµ , front-side metallization, emitter layer resistance, loss of power due to the reflection, etc. Let us carry out the optimization of the i-layer thickness taking into account that the efficiency is limited by light absorption for a thin film and by its resistance in the case of thick films. In Fig. 5, the η values versus the i-layer thickness are displayed. These curves were evaluated for different lifetimes in the i- layer, by using (8) and taking into account (9) and (10). The formation of a maximum on curves 1-6 is caused by the existing maxima of JSC(d) dependences at high values of the coefficient of reflection from the rear electrode. Another mechanism deals with the fact that the i-layer thickness expansion provides the escalation of the amount of generated pairs due to an increase in the absorption. But, on the other hand, extending the base layer enlarges the series resistance related to the i-layer photoresistance that produces the photovoltage drop. At low values of either the surface recombination velocity at the rear surface Sd or the photoelectron lifetime in the base, the maximum of η shifts to smaller values of emitter thickness, d; while larger Sd values cause the maximum to move to larger d. The calculated values of dm vary from 3·10−5 to 10−4 cm, which correlates with experimental data. As is well known, the Staebler-Wronski effect leads to a sufficient escalation of the density of states that causes a comparative reduction of the electron and hole lifetimes in the i-layer. Therefore, curves 1 and 2 can be interpreted as the efficiency of photoconversion in SC based on hydrogenated amorphous silicon according to the Staebler-Wronski effect. In this case, the efficiency maximum value reaches ~8 % that is com- parable with experimental data for high-efficient SC. Thus, from the data obtained, we can conclude that the thickness of the base (i-layer) is required to be optimized first as a key parameter. The following parameters were used in plotting the curves in Fig. 5: D = 2.5·10−3 cm2/s, Rd = 0.9; Е = 105 V/cm; RS = 0.082; JS = 10−12 A/cm2; А = 1.5; Т = 300 К; µn = 0.1 cm2/(V·s); µp = 102 cm2/(V·s); N = 1015 cm−2; m = 0.05. The values JS = 10−12 А/cm2 and А = 1.5 are taken from [1] and correspond to the parameters of high-efficient SC. Semiconductor Physics, Quantum Electronics & Optoelectronics, 2007. V. 10, N 4. P. 60-66. © 2007, V. Lashkaryov Institute of Semiconductor Physics, National Academy of Sciences of Ukraine 65 0.5 1.0 1.5 2.0 2.5 3.0 6 7 8 9 10 11 12 Ef fic ie nc y η (% ) Film thickness d (µm) 1 2 4 3 5 6 Fig. 5. Efficiency as a function of the i-layer thickness. For curves 1, 3, and 5, n0 = 1015 cm−3; for curves 2, 4, and 6, n0 = 1014 cm−3. For curves 3 and 4, S0 = 105 cm/s, Sd = 104 cm/s. For curves 1, 2, 5, and 6, S0 = Sd = 102 cm/s. For curves 5 and 6, τ = 10−7 s, and, for curves 1, 2, 3, and 4, τ = 10−6 s. The used values of µn and τn lie within the corresponding ranges realizable in amorphous hydrogenated silicon (though they should be specified), while the values µp = 102 cm2/(V⋅s) and N = 1015 cm−2 (belonging to the ITO parameters) were taken arbitrarily. Thus, the typical ITO р-layer parameters (hole mobility and concentration) are to be found in the literature. It was stated in [1, p. 387] that the doped contact р- and n-layers in SC made of hydrogenated silicon are 5– 10 nm thick. According to the more accurate data given in the report by B.A. Korevaar (Eindhoven Technical University), the thickness of a р-(n-)layer is 10 (20) nm. In our calculations, we took dp = 10 nm. Strictly speaking, one should consider (taking into account the single reflection of light from the metallized back surface) that light passes a distance 2(d + dp + dn), where dn is the heavily doped n-layer thickness. In the above equations, the quantity dn was not allowed for. This led to a change of the reflection coefficients in the р- and n- layers, since the mobility gap values in them differ from that in the i-region. However, if the conditions d >> dp, dn and αdp ≤ 1 are satisfied, the corrections due to the allowance for the distinction in α values are sufficiently small. Here, we assumed, when taking the absorption in the p-layer into account, that the coefficient of light absorption in it is the same as that in the i-region. Now let us consider the optimization of the contact grid parameters. In this case, the optimization is made with allowance for the fact that, at a small shadowing (i.e., at a big distance between the fingers), the series resistance increases, while the great shadowing (i.e., at a small distance between the fingers) leads to a decrease in the photocurrent. In Fig. 6, we show the photoconversion efficiency vs the relative degree of metallization, the finger width ln serving as a parameter. When plotting the curves in Fig. 6, we considered that the product of the hole concentration in the р-layer and the layer thickness is 1015 (1014) cm−2. An important parameter in this case is also the hole concentration (while plotting the curves in Figs. 5 and 6, we took it to be 102 cm2/(V⋅s)). One can see from the above figures that the bigger the finger width, the higher is the metallization degree, at which the maximal photoconversion efficiency is realized and the lower is the peak photoconversion efficiency. As the values of µp and N decrease, one should choose smaller values of the finger width. 4. Conclusions The primary need to optimize the base region (i-layer) thickness can be concluded from the obtained data. Efficiency as a function of metallization has a maximum. The efficiency maximum value reaches ~8 % that is comparable with experimental data for high-efficient SC. References 1. R.A. Street, Hydrogenated Amorphous Silicon. Cambridge University Press, Cambridge, 1991. 2. A.L. Fahrenbruch, R.H. Bube, Fundamentals of Solar Cells. Photovoltaic Solar Energy Conversion. Academic Press, New York, 1983. 3. Proc. 19th European Photovoltaic Solar Energy Conference, 7-11 June 2004, Paris, France, p. 1374-1377. 4. E.A. Schiff // Solar Energy Materials & Solar Cells 78, p. 567-595 (2003). 5. E.A. Schiff // Journal of Non-Crystalline Solids 352(9), p. 1087-1092 (2006). 6. T. Searle (Ed.), Properties of amorphous silicon and its alloys, INSPEC. The Institute of Electrical Engineers, London, 1998. 1 10 0 2 4 6 8 10 E ffi ci en cy η (% ) Metallization m (%) 5 20 12 3 4 Fig. 6. Efficiency as a function of metallization. L = 3·10−3, 10−2, 10−1, and 3·10−1 cm for curves 1, 2, 3, and 4, respectively. Semiconductor Physics, Quantum Electronics & Optoelectronics, 2007. V. 10, N 4. P. 60-66. © 2007, V. Lashkaryov Institute of Semiconductor Physics, National Academy of Sciences of Ukraine 66 7. A.V. Sachenko // Soviet Phys. Semicond. 19, p. 903 (1985). 8. L. Jiao, I. Chen, R.W. Collins, C.R. Wronski, and N. Hata // Appl. Phys. Lett. 72(9), p. 1057-1059 (1998). 9. A.P. Gorban, V. P. Kostylyov, A.V. Sachenko // Semiconductor Physics, Quantum Electronics & Optoelectronics 2(3), p. 26-31 (1999). 10. P.H. Berning, in: G. Hass (Ed.), Physics of Thin Films, Vol. 1. Academic Press, New York, 1963, p. 69-121. 11. A.V. Sachenko, A. I. Chkrebtii // Ukr. Fizicheskiy Zhurnal 29(12), p. 1855-1859 (1984) (in Russian). 12. A.V. Sachenko, A.I. Chkrebtii, T.V. Panichevskaya // Optoelektronika i poluprovodnikovaya tekhnika 21, p. 63-70 (1991) (in Ukrainian). 13. A.V. Sachenko, A.P. Gorban // Semiconductor Physics, Quantum Electronics & Optoelectronics 2(2), p. 42-44 (1999) (http://www.journal- spqeo.org.ua/n2_99/42_299.htm).