Solar cells rely on photogeneration of charge carriers in p-n junctions and their transport and subsequent recombination in the quasineutral regions. A number of basic issues concerning the physics of the operation of solar cells still remain obscure. This paper discusses some of those unsolved basic problems. In conventional solar cells, recombination of photogenerated charge carriers plays a major limiting role in the cell efficiency. High quality thin-film solar cells may overcome this limit if the minority diffusion lengths become large as compared to the cell dimensions, but, strikingly, the conventional model fails to describe the cell electric behavior under these conditions. A new formulation of the basic equations describing charge carrier transport in the cell along with a set of boundary conditions is presented. An analytical closed-form solution is obtained under a linear approximation. In the new framework given, the calculation of the open-circuit voltage of the solar cell diode does not lead to unphysical results.
1. Introduction
Charge carrier transport underlies the electrical behavior of any semiconductor device and, in particular, of solar cells. Despite efforts made over the years to correctly model charge-carrier transport in semiconductors, some important questions still remain unanswered. These questions need to be addressed in order to correctly model present and future devices.
One of these open questions is how to model carrier recombination. The mathematical expression routinely used to model the recombination rate [1] is basically incorrect. It has recently been demonstrated that it contradicts Maxwell’s equations, and hence a new corrected model must be developed [2, 3]. This problem was partly addressed in previous works (see, e.g., [4]). Recombination is a key feature when describing carrier transport in solar cells because it strongly affects the electrical response of the semiconductor at all levels of external excitation. A remark that must be made about the need for a correct modeling of recombination is that in devices operating under a strong excitation regime (the operation of solar cells lies in this regime), the importance of a correct formulation of the recombination terms is even more important.
Since in general the set of Poisson and transport equations cannot be solved analytically, some simplifications must be introduced to obtain a closed-form solution. One approximation commonly used to solve this system of equations is the assumption of quasineutrality (QN) [5]. The use of the QN approximation is acceptable if the sample’s and the diffusion lengths are both larger than the Debye length. Although QN has been routinely used in semiconductor device modeling for many years, the role of space charge in the formation of the current-voltage (I-V) characteristic in a semiconductor is still controversial.
A very important question is the choice of the boundary conditions used when solving the carrier-transport equations. It should be noted that the expressions commonly used are valid only for semiconductor devices operating in open-circuit conditions (see, e.g., [6]). Since in normal operation a current flows at the terminals, the widespread use of boundary conditions for open-circuit conditions is incorrect. For closed-circuit conditions, a different set of boundary conditions needs to be derived. This problem has only been addressed in the last few years [7–9].
All the abovementioned issues need to be addressed when modeling any semiconductor devices. Thin-film solar cells are one of the strongest technologies in the steadily growing photovoltaic market [10]. In this paper we will address a very basic problem encountered when modeling cells with large values of the minority diffusion lengths.
2. Problem Statement
The commonly accepted electrical model for the current density-voltage characteristic (J-V) of an ideal p-n junction is [11]J=J0(eV/VT-1),
where J is the diode current density, V is the bias voltage, VT is the thermal voltage, and J0 the reverse-saturation current density that for a long abrupt diode isJ0=VT(σpnLDn+σnpLDp),σnp (σpn) is the electron (hole) conductivity of minority carriers in the p-side (n-side) of a junction diode; LDn and LDp are, respectively, the hole and electron minority carrier diffusion lengths of minority carriers in n- and p- sides defined as
LDn=Dpτp,LDp=Dnτn,Dn (Dp) is the electron (hole) diffusivity and τn (τp) is the electron (hole) minority carrier lifetime.
Equation (1) may be rewritten for a solar cell under illumination as [12]-J=J0(eV/VT-1)+JI,
where JI is the photocurrent and J is additionally related to the generated voltage V through the external resistive load (R) in the circuit and the solar cell area (S) as shown in Figure 1:J=VRS.
Solar cell diode with an external resistor load.
For the sake of simplicity, let us examine the photocurrent generated in a slab of semiconductor of thickness L with incident light at x=0. JI may be obtained as the integral of the light intensity across the device I(x) as|JI|=qαβ∫0LI(x)dx,
where q is the absolute value of the electron charge, α is the absorption coefficient of light, and β is the quantum efficiency,β=α-α0α.
The light intensity inside the semiconductor slab will decrease following the law:I=I0e-αx,
where I0 is the light intensity (number of incident by photons per unit time and unit area) at the plane of the incident light (x=0) and the carrier generation rate due to photon-absorption can be written as g=(α-α0)I0e-αx.
The set of (4)–(9) may be extended or augmented with other equations for a finer modeling of a real technology [13], but this set constitutes the fundamental framework to describe the electrical behaviour of any diode-based solar cell. Let us assume a one-dimensional problem for a single-diode solar cell such as the one given in Figure 1. According to the values of diffusion lengths in Si, in a c-Si thin-film solar cell, or in a-Si ones of nanometer dimensions [14], the following inequality will hold:ln,p≪LDn,p.
From (10) it follows that LDn,p→∞ and τn,p→∞; therefore, recombination is negligible and according to (2) J0→0 and, moreover, from (4), the basic model leads to the unphysical result (see Figure 2(b)): J≈-JI, i.e., the p-n junction behaves as a current-source that translates photo-excitation into current at any value of the diode-voltage independently of the p-n junction properties.
J-V characteristic of an ideal solar cell with the short- and open-circuit expressions for the axis-intercept points Voc and JSH, respectively, (a). J-V characteristic in the limiting case of a thin-film cell ln,p≪LDn,p (b).
From the above it follows that the framework of basic equations conventionally used to model solar cells fails in high-quality thin-film solar cells [15]. Regardless of its apparent simplicity, since this framework constitutes the foundations of any electrical modeling of a solar cell diode, the discussed problem will persist in any other model that one can build; therefore, it is crucial to identify the origin of the problem and modify the basic model to build a new framework free from unphysical errors. In this paper we will present a new model and we will analytically show in a simple case that this new framework does not fail to correctly describe the limiting case of a thin-film solar cell.
3. New Model for Thin-Film Solar Cells
Let us assume once again that inequality of (10) holds and, accordingly, the solar cell recombination is negligible. The macroscopic description of the transport of nonequilibrium charge is done with the continuity equations for the electron (jn) and hole (jp) current densities [3]:∂n∂t=gn+1edivjn-Rn,∂p∂t=gp-1edivjp-Rp.
In a stationary case and one-dimensional problem along the x-direction, normal to the p-n junction plane as in Figure 1, (11) becomesdjndx=q(Rn-gn),djpdx=-q(Rp-gp).
Since volume recombination is negligible Rn=Rp=0 and, accordingly, gn=gp=g(x), that is, the recombination is band-to-band and the continuity equations reduce todjndx=-qg(x),djpdx=qg(x).
Additionally, it can be shown that in absence of recombination and under band-to-band generation, the density of charge in any impurity level remains constant [15], under all the above conditions the QN approximation reduces to δn=δp. In a general case, for a nondegenerate semiconductor, in a linear approximation δn and δp can be written in terms of the variations in their respective chemical potentials (μn and μp) as [16]δn=n0T0δμn,δp=p0T0δμp,
where the subscript 0 stand for equilibrium value of the magnitude, and using the QN approximation the following relationship between the deviations from equilibrium of the chemical potentials holdsδμp=n0p0δμn.
In order, to further simplify the mathematical problem, while ensuring no loss of generality of the physical one, we will assume a step-like profile for the light absorption instead of the model in (8):I={I0-lp<x<0,00<x<ln.
Using (9) and (16), we may write the continuity equations, (13), asdjnpdx=-qg0,djppdx=qg0,djnndx=0,djpndx=0.
Superscripts in current densities in (17) indicate the doping type of the region. The current densities may be written as [3]jnn=-σnn(dδφndx-1qdδμnndx),jpn=-σpn(dδφndx+1qdδμpndx),jnp=-σnp(dδφpdx-1qdδμnpdx),jpp=-σpp(dδφpdx+1qdδμppdx),
where we introduced the electric potential (φ). To solve (17) with (14), (15), and (18), we need to determine and impose enough boundary conditions at the semiconductor interfaces [8]. We will assume ideal metallic (ohmic) contacts placed at x=-lp and at x=ln; therefore, the excess of carries are null:δnp(-lp)=δpp(-lp)=0,δnn(ln)=δpn(ln)=0.
And, according to (19), we obtain the following four boundary conditions for the chemical potentials at the metal-semiconductor contacts:δμnp(-lp)=δμpp(-lp)=0,δμnn(ln)=δμpn(ln)=0.
Two additional boundary conditions may be imposed on the electric potential at the same metal-semiconductor contacts:δφp(-lp)=0,δφn(ln)=V.
At the interface of the two semiconductor regions, since there is no recombination, we may write two additional boundary conditions at x=0 [8]:jnn(0)=jnp(0),jpn(0)=jpp(0).
Two more boundary conditions areδφp(0)-δμnp(0)q=δφn(0)-δμnn(0)q,δφp(0)+δμpp(0)q=δφn(0)+δμpn(0)q.
Using the obtained boundary conditions from the continuity of the Fermi quasilevels it is straightforward to solve the linear system of equations for the current densities (18) and obtain the solutionJ(lnσpnlpσip)+V(lpσnp+lnσpn)+qg02(lp2lnσnpσpn)=0.
In the third term of (24), we may identify the photocurrent density created under illuminationJI=qg02lp.
Since we are under a small voltage-/current-linear approximation with the aim of obtaining an analytical closed form for the J-V (24), we must compare the obtained equation with a linearization of (4):JI+V(J0VT)=-J.
By rewriting (24) in a similar form we obtainJI+V(lpσnp+lnσpn)(σnpσpnlpln)=-J.
This is a main result of this work, identifying the second terms in the left hand side of (26) and (27), we obtain a new expression for the current density J0:J0=(lpσnp+lnσpn)(σnpσpnlpln)VT.
This new expression does not exhibit an unphysical behavior (J0→∞) in thin-film solar cells (Figure 2(b)). Figure 3 gives the two axis-intercept points (Voc and JSH) of the J-V characteristic, it can be trivially verified that the x-intercept point (Voc) will remain at a finite value independently of the minority carriers diffusion lengths at each side of the junction and the diode dimensions. Since the analytical model presented in this paper was obtained by a linearization of the full model and some other approximations have been used to simplify the mathematical problem, the large-signal J-V is not properly obtained (Figure 3), but the framework of transport equations and boundary conditions can be solved numerically to study any real solar cell diode.
J-V characteristic of an ideal thin-film solar cell with the new expression for the x axis-intercept point Voc.
4. Conclusions
Solar cells rely on photogeneration of charge carriers in p-n junctions and their subsequent recombination in the quasi-neutral regions. A number of basic issues concerning the physics of the operation of solar cells still remain obscure. In this paper we briefly discussed some of those unsolved basic problems: need for not unphysical models of the carrier recombination and revision of the QN concept. In conventional solar cells recombination of photogenerated charge carriers plays a major limiting role in the cell efficiency. High quality thin-film solar cells may overcome this limit if the minority diffusion lengths become large as compared to the cell dimensions, but, strikingly, the conventional model fails to describe the cell electric behavior under these conditions. A new formulation of the basic equations describing charge carrier transport in the cell along with a set of boundary conditions is presented. An analytical closed-form solution is obtained under a linear approximation. In the new framework given, the calculation of the open-circuit voltage of the solar cell diode does not lead to unphysical results.
Acknowledgments
J. E. Velazquez-Perez wants to thank Ministerio de Ciencia e Innovación (MICINN) and Junta de Castilla y León for partial financial support under Grants TEC2008-02281 and SA061A09. Yu. G. Gurevich wants to thank CONACYT-México for partial financial support.
LandsbergP. T.1991Cambridge, UKCambridge UniversityVolovichevI. N.LogvinovG. N.TitovO. Y.GurevichY. G.Recombination and lifetimes of charge carriers in semiconductors2004958449444962-s2.0-234256258110.1063/1.1669074GurevichYu. G.Velazquez-PerezJ. E.Espejo-LopezG.VolovichevI. N.TitovO. Y.Transport of nonequilibrium carriers in bipolar semiconductors2007101282-s2.0-3384775864510.1063/1.2424502023705GurevichY. G.VolovichevI. N.Forgotten mechanism of nonlinearity in the theory of hot electrons19996011771577172-s2.0-0000174488VolovichevI. N.Velazquez-PerezJ. E.GurevichYu. G.Transport boundary condition for semiconductor structures20085211170317092-s2.0-5504914111010.1016/j.sse.2008.04.037SeegerK.1985Berlin, GermanySpringerTitovO. Yu.GiraldoJ.GurevichYu. G.Boundary conditions in an electric current contact2002801731082-s2.0-003719333110.1063/1.1473875VolovichevI. N.Velazquez-PerezJ. E.GurevichY. G.Transport boundary conditions for solar cells20099316102-s2.0-5594911610710.1016/j.solmat.2008.01.001GurevichY. G.LogvinovG. N.VelazquezJ. E.TitovO. Y.Transport and recombination in solar cells: new perspectives20079115-16140814112-s2.0-3444728084210.1016/j.solmat.2007.03.006GreenM. A.Thin-film solar cells: review of materials, technologies and commercial status200718supplement 1S15S192-s2.0-3454760597210.1007/s10854-007-9177-9SinghJ.2000New York, NY, USAJohn Wiley & SonsLuqueA.HegedusS.2003Chichester, UKJohn Wiley & Sons2009Mountain View, Calif, USASynopsysBergmannR. B.Crystalline Si thin-film solar cells: a review19996921871942-s2.0-003259326710.1007/s003390050989Velázquez-PérezJ. E.GurevichY. G.Transport and recombination in solar cells, new boundary conditions and perspectivesProceedings of the 25th European Photovoltaic Solar Energy Conference and Exhibition/5th World Conference on Photovoltaic Energy ConversionSeptember 2010Valencia, SpainGurevichYu. G.LogvinovG. N.VolovichevI. N.EspejoG.TitovO. Yu.MeriutsA.The role of non-equilibrium carriers in the formation of thermo-E.M.F. in bipolar semiconductors200223112782932-s2.0-0035994516