### Poisson's Equation

$\frac{dE}{dx}=\frac{\rho}{\epsilon}=\frac{q}{\epsilon}(p(x)-n(x)-{N}_{A}^{-}+{N}_{D}^{+})$

where **E** is the electric field, *ρ* is the charge density and *ε* is the material permittivity. This equation gives the basic relationship between charge and electric field strength. In semiconductors we divide the charge up into four components: hole density,* p*, electron density,* n*, acceptor atom density,* N _{A}* and donor atom density,

*N*.

_{D}For the ideal diode derivation *N _{A }*is assumed constant in the p-region and zero in the n-region. Similarly,

*N*is assumed constant in the n-region and zero in the p-region.

_{D}### Transport Equations

${J}_{n}=q{\mu}_{n}n\mathit{E}+q{D}_{n}\frac{dn}{dx},{J}_{p}=q{\mu}_{p}p\mathit{E}+q{D}_{p}\frac{dp}{dx}$

where *J _{n}* is the electron current density,

*μ*is the electron mobility and

_{n}*D*is the electron diffusivity. Similarly,

_{n}*J*is the hole current density,

_{p}*μ*is the hole mobility and

_{p}*D*is the hole diffusivity.

_{p}*q*is the electronic charge and

**E**is the electric field. Note that in this section we use

**E**, $\hat{E}$ and $\overrightarrow{E}$ interchangably for the electric field. Its a mistake and we should be consistent and stick to

**E**.

The transport equations describe how carriers move, i.e. the flow of carriers or current. Its easier to use the current density,* J*, in A/cm^{2} rather than the absolute current, *I*, in A since we are not concerned with the area of the device. *I = J × Area* and for a 1 cm^{2} device *J* and *I* are equal.

The first term in each equation is for drift and the second term is for diffusion.

### Continuity Equations

The continuity equation keep track of all the carriers in terms of movement, generation and recombination. They are sometimes temed 'book keeping' equations since they make sure that every carrier is accounted for.

**General Conditions **

where *U *is the recombinaton rate and *G *is the generation rate.

Solar cells operate in steady state and we are not concerned with transients or switching times. Under thermal equilibrium and steady state conditions the carrier concetrations do not change with time so that:

$\frac{dn}{dt}=\frac{dp}{dt}=0$

Rearranging the equations above leads to:

### Summary

We now have the five basic equations to solve:

$\frac{dE}{dx}=\frac{\rho}{\epsilon}=\frac{q}{\epsilon}(p(x)-n(x)-{N}_{A}^{-}+{N}_{D}^{+})$

${J}_{n}=q{\mu}_{n}n\mathit{E}+q{D}_{n}\frac{dn}{dx},{J}_{p}=q{\mu}_{p}p\mathit{E}+q{D}_{p}\frac{dp}{dx}$

,

The equations are readily solved using numerical approaches[1], and there are many device simulators that perform this task. By making a few approximations it is also possible to solve the equations closed form as outlined on the following pages.

- 1. , “Numerical modeling of textured silicon solar cells using PC-1D”, Electron Devices, IEEE Transactions on, vol. 37, pp. 337 -343, 1990.