# Theory of solar cells

The theory of solar cells explains the process by which light energy in photons is converted into electric current when the photons strike a suitable semiconductor device. The theoretical studies are of practical use because they predict the fundamental limits of a solar cell, and give guidance on the phenomena that contribute to losses and solar cell efficiency.

## Working explanation

* Photons in sunlight hit the solar panel and are absorbed by semi-conducting materials.
* Electrons (negatively charged) are knocked loose from their atoms as they are excited. Due to their special structure and the materials in solar cells, the electrons are only allowed to move in a single direction. The electronic structure of the materials is very important for the process to work, and often silicon incorporating small amounts of boron or phosphorus is used in different layers.
* An array of solar cells converts solar energy into a usable amount of direct current (DC) electricity.

## Photogeneration of charge carriers

When a photon hits a piece of semiconductor, one of three things can happen:
* The photon can pass straight through the semiconductor — this (generally) happens for lower energy photons.
* The photon can reflect off the surface.
* The photon can be absorbed by the semiconductor if the photon energy is higher than the band gap value. This generates an electron-hole pair and sometimes heat depending on the band structure.
When a photon is absorbed, its energy is given to an electron in the crystal lattice. Usually this electron is in the valence band. The energy given to the electron by the photon "excites" it into the conduction band where it is free to move around within the semiconductor. The network of covalent bonds that the electron was previously a part of now has one fewer electron. This is known as a hole, and it has positive charge. The presence of a missing covalent bond allows the bonded electrons of neighboring atoms to move into the "hole", leaving another hole behind, thus propagating holes throughout the lattice in the opposite direction to the movement of the negatively electrons. It can be said that photons absorbed in the semiconductor create electron-hole pairs.
A photon only needs to have energy greater than that of the band gap in order to excite an electron from the valence band into the conduction band. However, the solar frequency spectrum approximates a black body spectrum at about 5,800 K, and as such, much of the solar radiation reaching the Earth is composed of photons with energies greater than the band gap of silicon (1.12eV), which is near to the ideal value for a terrestrial solar cell (1.4eV). These higher energy photons will be absorbed by a silicon solar cell, but the difference in energy between these photons and the silicon band gap is converted into heat (via lattice vibrations — called phonons) rather than into usable electrical energy.

## The p–n junction

The most commonly known solar cell is configured as a large-area p–n junction made from silicon. As a simplification, one can imagine bringing a layer of n-type silicon into direct contact with a layer of p-type silicon. n-type doping produces mobile electrons (leaving behind positively charged donors) while p-type doping produces mobile holes (and negatively charged acceptors). In practice, p–n junctions of silicon solar cells are not made in this way, but rather by diffusing an n-type dopant into one side of a p-type wafer (or vice versa).
If a piece of p-type silicon is placed in close contact with a piece of n-type silicon, then a diffusion of electrons occurs from the region of high electron concentration (the n-type side of the junction) into the region of low electron concentration (p-type side of the junction). When the electrons diffuse into the p-type side, each one annihilates a hole, making that side net negatively charged (because now the number of mobile positive holes is now less than the number of negative acceptors). Similarly, holes diffusing to the n-type side make it more positively charged. However (in the absence of an external circuit) this diffusion current of carriers does not go on indefinitely because the charge build up on either side of the junction produces an electric field that opposes further diffusion of more charges. Eventually, an equilibrium is reached where the net current is zero, leaving a region either side of the junction where electrons and holes have diffused across the junction and annihilated each other called the depletion region because it contains practically no mobile charge carriers. It is also known as the space charge region, although space charge extends a bit further in both directions than the depletion region.
Once equilibrium is established, electron-hole pairs generated in the depletion region are separated by the electric field, with the electron attracted to the positive n-type side and holes to the negative p-type side, reducing the charge (and the electric field) built up by the diffusion just described. If the device is unconnected (or the external load is very high) then diffusion current would eventually restore the equilibrium charge by bringing the electron and hole back across the junction, but if the load connected is small enough, the electrons prefer to go around the external circuit in their attempt to restore equilibrium, doing useful work on the way.

## Charge carrier separation

There are two causes of charge carrier motion and separation in a solar cell:
* drift of carriers, driven by the electric field, with electrons being pushed one way and holes the other way
* diffusion of carriers from zones of higher carrier concentration to zones of lower carrier concentration (following a gradient of chemical potential).
These two "forces" may work one against the other at any given point in the cell. For instance, an electron moving through the junction from the p region to the n region (as in the diagram at the beginning of this article) is being pushed by the electric field against the concentration gradient. The same goes for a hole moving in the opposite direction.
It is easiest to understand how a current is generated when considering electron-hole pairs that are created in the depletion zone, which is where there is a strong electric field. The electron is pushed by this field toward the n side and the hole toward the p side. (This is opposite to the direction of current in a forward-biased diode, such as a light-emitting diode in operation.) When the pair is created outside the space charge zone, where the electric field is smaller, diffusion also acts to move the carriers, but the junction still plays a role by sweeping any electrons that reach it from the p side to the n side, and by sweeping any holes that reach it from the n side to the p side, thereby creating a concentration gradient outside the space charge zone.
In thick solar cells there is very little electric field in the active region outside the space charge zone, so the dominant mode of charge carrier separation is diffusion. In these cells the diffusion length of minority carriers (the length that photo-generated carriers can travel before they recombine) must be large compared to the cell thickness. In thin film cells (such as amorphous silicon), the diffusion length of minority carriers is usually very short due to the existence of defects, and the dominant charge separation is therefore drift, driven by the electrostatic field of the junction, which extends to the whole thickness of the cell.
Once the minority carrier enters the drift region, it is 'swept' across the junction and, at the other side of the junction, becomes a majority carrier. This reverse current is a generation current, fed both thermally and (if present) by the absorption of light. On the other hand, majority carriers are driven into the drift region by diffusion (resulting from the concentration gradient), which leads to the forward current; only the majority carriers with the highest energies (in the so-called Boltzmann tail; cf. Maxwell–Boltzmann statistics) can fully cross the drift region. Therefore, the carrier distribution in the whole device is governed by a dynamic equilibrium between reverse current and forward current.

## Connection to an external load

Ohmic metal-semiconductor contacts are made to both the n-type and p-type sides of the solar cell, and the electrodes connected to an external load. Electrons that are created on the n-type side, or created on the p-type side, "collected" by the junction and swept onto the n-type side, may travel through the wire, power the load, and continue through the wire until they reach the p-type semiconductor-metal contact. Here, they recombine with a hole that was either created as an electron-hole pair on the p-type side of the solar cell, or a hole that was swept across the junction from the n-type side after being created there.
The voltage measured is equal to the difference in the quasi Fermi levels of the majority carriers (electrons in the n-type portion and holes in the p-type portion) at the two terminals.

## Equivalent circuit of a solar cell

An equivalent circuit model of an ideal solar cell's p–n junction uses an ideal current source (whose photogenerated current $I_{\text{L}}$ increases with light intensity) in parallel with a diode (whose current $I_{\text{D}}$ represents recombination losses). To account for resistive losses, a shunt resistance $R_{\text{SH}}$ and a series resistance $R_{\text{S}}$ are added as lumped elements. The resulting output current $I_{\text{out}}$ equals the photogenerated current minus the currents through the diode and shunt resistor:
    $I_{\text{out}}=I_{\text{L}}-I_{\text{D}}-I_{\text{SH}}$
The junction voltage (across both the diode and shunt resistance) is:
    $V_{\text{j}}=V_{\text{out}}+I_{\text{out}}\,R_{\text{S}}$
where $V_{\text{out}}$ is the voltage across the output terminals. The leakage current $I_{\text{SH}}$ through the shunt resistor is proportional to the junction's voltage $V_{\text{j}}$, according to Ohm's law:
    $I_{\text{SH}}={\frac {V_{\text{j}}}{R_{\text{SH}}}}$
By the Shockley diode equation, the current diverted through the diode is:
    $I_{\text{D}}=I_{0}\left\{\exp \left[{\frac {V_{\text{j}}}{nV_{\text{T}}}}\right]-1\right\}$
where
* I$_0$, reverse saturation current
* n, diode ideality factor (1 for an ideal diode)
* q, elementary charge
* k, Boltzmann constant
* T, absolute temperature
* $V_{\text{T}}=kT/q,$ the thermal voltage. At 25 °C, $V_{\text{T}}\approx 0.0259$ volt.
Substituting these into the first equation produces the characteristic equation of a solar cell, which relates solar cell parameters to the output current and voltage:
    $I_{\text{out}}=I_{\text{L}}-I_{0}\left\{\exp \left[{\frac {V_{\text{out}}+I_{\text{out}}R_{\text{S}}}{nV_{\text{T}}}}\right]-1\right\}-{\frac {V_{\text{out}}+I_{\text{out}}R_{\text{S}}}{R_{\text{SH}}}}.$
An alternative derivation produces an equation similar in appearance, but with $V_{\text{out}}$ on the left-hand side. The two alternatives are identities; that is, they yield precisely the same results.
Since the parameters I$_0$, n, R$_S$, and R$_{SH}$ cannot be measured directly, the most common application of the characteristic equation is nonlinear regression to extract the values of these parameters on the basis of their combined effect on solar cell behavior.
When R$_S$ is not zero, the above equation does not give $I_{\text{out}}$ directly, but it can then be solved using the Lambert W function:
    $I_{\text{out}}={\frac {(I_{\text{L}}+I_{0})-V_{\text{out}}/R_{\text{SH}}}{1+R_{\text{S}}/R_{\text{SH}}}}-{\frac {nV_{\text{T}}}{R_{\text{S}}}}W\left({\frac {I_{0}R_{\text{S}}}{nV_{\text{T}}(1+R_{\text{S}}/R_{\text{SH}})}}\exp \left({\frac {V_{\text{out}}}{nV_{\text{T}}}}\left(1-{\frac {R_{\text{S}}}{R_{\text{S}}+R_{\text{SH}}}}\right)+{\frac {(I_{\text{L}}+I_{0})R_{\text{S}}}{nV_{\text{T}}(1+R_{\text{S}}/R_{\text{SH}})}}\right)\right)$
When an external load is used with the cell, its resistance can simply be added to R$_S$ and $V_{\text{out}}$ set to zero in order to find the current.
When $I_{0}R_{\text{S}}/nV_{\text{T}}$ is small, we can use the approximation $x^{-1}W\left(xy\right)\to y$ as $x\to 0$ to produce something much easier to work with
    $I_{\text{out}}={\frac {(I_{\text{L}}+I_{0})-V_{\text{out}}/R_{\text{SH}}}{1+R_{\text{S}}/R_{\text{SH}}}}-{\frac {I_{0}}{(1+R_{\text{S}}/R_{\text{SH}})}}\exp \left({\frac {V_{\text{out}}}{nV_{\text{T}}}}\left(1-{\frac {R_{\text{S}}}{R_{\text{S}}+R_{\text{SH}}}}\right)+{\frac {(I_{\text{L}}+I_{0})R_{\text{S}}}{nV_{\text{T}}(1+R_{\text{S}}/R_{\text{SH}})}}\right).$
Several further simplifications are now possible, such as when $R_{\text{S}}\ll R_{\text{SH}}$ which leads to
    $I_{\text{out}}=I_{\text{L}}+I_{0}-{\frac {V_{\text{out}}}{R_{\text{SH}}}}-I_{0}\exp \left({\frac {V_{\text{out}}+(I_{\text{L}}+I_{0})R_{\text{S}}}{nV_{\text{T}}}}\right).$
When the current generated by the PV is large compared with the current in the shunt, i.e. $I_{\text{L}}\gg V_{\text{out}}/R_{\text{SH}}$ (because the shunt resistance is large) there is an analytical solution for $V_{\text{out}}$ for any $I_{\text{out}}$ less than $I_{\text{L}}+I_{0}$:
    $V_{\text{out}}=nV_{\text{T}}\ln \left({\frac {I_{\text{L}}-I_{\text{out}}}{I_{0}}}+1\right)-(I_{\text{L}}+I_{0})R_{\text{S}}.$
Otherwise one can solve for $V_{\text{out}}$ using the Lambert W function:
    $V=(I_{\text{L}}+I_{0})R_{\text{SH}}-I(R_{\text{S}}+R_{\text{SH}})-nV_{\text{T}}W\left({\frac {I_{0}R_{\text{SH}}}{nV_{\text{T}}}}\exp \left({\frac {(I_{\text{L}}+I_{0}-I)R_{\text{SH}}}{nV_{\text{T}}}}\right)\right)$
However, when R$_{SH}$ is large it's better to solve the original equation numerically.
The general form of the solution is a curve with $I_{\text{out}}$ decreasing as $V_{\text{out}}$ increases (see graphs lower down). The slope at small or negative $V_{\text{out}}$ (where the W function is near zero) approaches $-1/(R_{\text{S}}+R_{\text{SH}})$, whereas the slope at high $V_{\text{out}}$ approaches $-1/R_{\text{S}}$. Therefore for high optimum output power $P_{\text{out}}=I_{\text{out}}V_{\text{out}}$, it is desirable to have $R_{\text{SH}}$ large and $R_{\text{S}}$ should be small.

### Open-circuit voltage and short-circuit current

When the cell is operated at open circuit, $I_{\text{out}}$= 0 and the voltage across the output terminals is defined as the open-circuit voltage. Assuming the shunt resistance is high enough to neglect the final term of the characteristic equation, the open-circuit voltage V$_{OC}$ is:
    $V_{\text{OC}}\approx {\frac {nkT}{q}}\ln \left({\frac {I_{\text{L}}}{I_{0}}}+1\right).$
Similarly, when the cell is operated at short circuit, $V_{\text{out}}$ = 0 and the current $I_{\text{SC}}$ through the terminals is defined as the short-circuit current. It can be shown that for a high-quality solar cell (low R$_S$ and I$_0$, and high R$_{SH}$) the short-circuit current is:
    $I_{\text{SC}}\approx I_{\text{L}}.$
It is not possible to extract any power from the device when operating at either open circuit or short circuit conditions.

### Effect of physical size

The values of I$_L$, I$_0$, R$_S$, and R$_{SH}$ are dependent upon the physical size of the solar cell. In comparing otherwise identical cells, a cell with twice the junction area of another will, in principle, have double the I$_L$ and I$_0$ because it has twice the area where photocurrent is generated and across which diode current can flow. By the same argument, it will also have half the R$_S$ of the series resistance related to vertical current flow; however, for large-area silicon solar cells, the scaling of the series resistance encountered by lateral current flow is not easily predictable since it will depend crucially on the grid design (it is not clear what "otherwise identical" means in this respect). Depending on the shunt type, the larger cell may also have half the R$_{SH}$ because it has twice the area where shunts may occur; on the other hand, if shunts occur mainly at the perimeter, then R$_{SH}$ will decrease according to the change in circumference, not area.
Since the changes in the currents are the dominating ones and are balancing each other, the open-circuit voltage is practically the same; V$_{OC}$ starts to depend on the cell size only if R$_{SH}$ becomes too low. To account for the dominance of the currents, the characteristic equation is frequently written in terms of current density, or current produced per unit cell area:
    $J=J_{\text{L}}-J_{0}\left\{\exp \left[{\frac {q(V_{\text{out}}+Jr_{\text{S}})}{nkT}}\right]-1\right\}-{\frac {V_{\text{out}}+Jr_{\text{S}}}{r_{\text{SH}}}}$
where
* J, current density (ampere/cm$^2$)
* J$_L$, photogenerated current density (ampere/cm$^2$)
* J$_0$, reverse saturation current density (ampere/cm$^2$)
* r$_S$, specific series resistance (Ω·cm$^2$)
* r$_{SH}$, specific shunt resistance (Ω·cm$^2$).
This formulation has several advantages. One is that since cell characteristics are referenced to a common cross-sectional area they may be compared for cells of different physical dimensions. While this is of limited benefit in a manufacturing setting, where all cells tend to be the same size, it is useful in research and in comparing cells between manufacturers. Another advantage is that the density equation naturally scales the parameter values to similar orders of magnitude, which can make numerical extraction of them simpler and more accurate even with naive solution methods.
There are practical limitations of this formulation. For instance, certain parasitic effects grow in importance as cell sizes shrink and can affect the extracted parameter values. Recombination and contamination of the junction tend to be greatest at the perimeter of the cell, so very small cells may exhibit higher values of J$_0$ or lower values of R$_{SH}$ than larger cells that are otherwise identical. In such cases, comparisons between cells must be made cautiously and with these effects in mind.
This approach should only be used for comparing solar cells with comparable layout. For instance, a comparison between primarily quadratical solar cells like typical crystalline silicon solar cells and narrow but long solar cells like typical thin film solar cells can lead to wrong assumptions caused by the different kinds of current paths and therefore the influence of, for instance, a distributed series resistance contribution to r$_S$. Macro-architecture of the solar cells could result in different surface areas being placed in any fixed volume - particularly for thin film solar cells and flexible solar cells which may allow for highly convoluted folded structures. If volume is the binding constraint, then efficiency density based on surface area may be of less relevance.

### Transparent conducting electrodes

Transparent conducting electrodes are essential components of solar cells. It is either a continuous film of indium tin oxide or a conducting wire network, in which wires are charge collectors while voids between wires are transparent for light. An optimum density of wire network is essential for the maximum solar cell performance as higher wire density blocks the light transmittance while lower wire density leads to high recombination losses due to more distance traveled by the charge carriers.

### Cell temperature

Temperature affects the characteristic equation in two ways: directly, via T in the exponential term, and indirectly via its effect on I$_0$ (strictly speaking, temperature affects all of the terms, but these two far more significantly than the others). While increasing T reduces the magnitude of the exponent in the characteristic equation, the value of I$_0$ increases exponentially with T. The net effect is to reduce V$_{OC}$ (the open-circuit voltage) linearly with increasing temperature. The magnitude of this reduction is inversely proportional to V$_{OC}$; that is, cells with higher values of V$_{OC}$ suffer smaller reductions in voltage with increasing temperature. For most crystalline silicon solar cells the change in V$_{OC}$ with temperature is about −0.50%/°C, though the rate for the highest-efficiency crystalline silicon cells is around −0.35%/°C. By way of comparison, the rate for amorphous silicon solar cells is −0.20 to −0.30%/°C, depending on how the cell is made.
The amount of photogenerated current I$_L$ increases slightly with increasing temperature because of an increase in the number of thermally generated carriers in the cell. This effect is slight, however: about 0.065%/°C for crystalline silicon cells and 0.09% for amorphous silicon cells.
The overall effect of temperature on cell efficiency can be computed using these factors in combination with the characteristic equation. However, since the change in voltage is much stronger than the change in current, the overall effect on efficiency tends to be similar to that on voltage. Most crystalline silicon solar cells decline in efficiency by 0.50%/°C and most amorphous cells decline by 0.15−0.25%/°C. The figure above shows I-V curves that might typically be seen for a crystalline silicon solar cell at various temperatures.

### Series resistance

As series resistance increases, the voltage drop between the junction voltage and the terminal voltage becomes greater for the same current. The result is that the current-controlled portion of the I-V curve begins to sag toward the origin, producing a significant decrease in $V_{\text{out}}$ and a slight reduction in I$_{SC}$, the short-circuit current. Very high values of R$_S$ will also produce a significant reduction in I$_{SC}$; in these regimes, series resistance dominates and the behavior of the solar cell resembles that of a resistor. These effects are shown for crystalline silicon solar cells in the I-V curves displayed in the figure to the right.
Power lost through the series resistance is $I_{\text{out}}^{2}R_{\text{S}}$. During illumination when $I_{\text{D}}$ and $I_{\text{SH}}$ are small relative to photocurrent $I_{\text{L}}$, power loss also increases quadratically with $I_{\text{L}}$. Series resistance losses are therefore most important at high illumination intensities.

### Shunt resistance

As shunt resistance decreases, the current diverted through the shunt resistor increases for a given level of junction voltage. The result is that the voltage-controlled portion of the I-V curve begins to sag far from the origin, producing a significant decrease in $I_{\text{out}}$ and a slight reduction in V$_{OC}$. Very low values of R$_{SH}$ will produce a significant reduction in V$_{OC}$. Much as in the case of a high series resistance, a badly shunted solar cell will take on operating characteristics similar to those of a resistor. These effects are shown for crystalline silicon solar cells in the I-V curves displayed in the figure to the right.

### Reverse saturation current

If one assumes infinite shunt resistance, the characteristic equation can be solved for V$_{OC}$:
    $V_{\text{OC}}={\frac {kT}{q}}\ln \left({\frac {I_{\text{SC}}}{I_{0}}}+1\right).$
Thus, an increase in I$_0$ produces a reduction in V$_{OC}$ proportional to the inverse of the logarithm of the increase. This explains mathematically the reason for the reduction in V$_{OC}$ that accompanies increases in temperature described above. The effect of reverse saturation current on the I-V curve of a crystalline silicon solar cell are shown in the figure to the right. Physically, reverse saturation current is a measure of the "leakage" of carriers across the p–n junction in reverse bias. This leakage is a result of carrier recombination in the neutral regions on either side of the junction.

### Ideality factor

The ideality factor (also called the emissivity factor) is a fitting parameter that describes how closely the diode's behavior matches that predicted by theory, which assumes the p–n junction of the diode is an infinite plane and no recombination occurs within the space-charge region. A perfect match to theory is indicated when n = 1. When recombination in the space-charge region dominate other recombination, however, n = 2. The effect of changing ideality factor independently of all other parameters is shown for a crystalline silicon solar cell in the I-V curves displayed in the figure to the right.
Most solar cells, which are quite large compared to conventional diodes, well approximate an infinite plane and will usually exhibit near-ideal behavior under standard test conditions (n ≈ 1). Under certain operating conditions, however, device operation may be dominated by recombination in the space-charge region. This is characterized by a significant increase in I$_0$ as well as an increase in ideality factor to n ≈ 2. The latter tends to increase solar cell output voltage while the former acts to erode it. The net effect, therefore, is a combination of the increase in voltage shown for increasing n in the figure to the right and the decrease in voltage shown for increasing I$_0$ in the figure above. Typically, I$_0$ is the more significant factor and the result is a reduction in voltage.
Sometimes, the ideality factor is observed to be greater than 2, which is generally attributed to the presence of Schottky diode or heterojunction in the solar cell. The presence of a heterojunction offset reduces the collection efficiency of the solar cell and may contribute to low fill-factor.

### Other models

While the above model is most common, other models have been proposed, like the d1MxP discrete model.