The planar sheath equation The quantitative physics of the Debye sheath is determined by four phenomena:
Energy conservation of the ions: If we assume for simplicity cold ions of mass m_\mathrm{i} entering the sheath with a velocity u_0, having charge opposite to the electron,
conservation of energy in the sheath potential requires :\frac{1}{2}m_\mathrm{i}\,u(x)^2 = \frac{1}{2}m_\mathrm{i}\,u_0^2 - e\,\varphi(x), where e is the charge of the electron taken positively, i.e. e=1.602 x 10^{-19} \mathrm{C}.
Ion continuity: In the steady state, the ions do not build up anywhere, so the flux is everywhere the same: :n_0\,u_0 = n_\mathrm{i}(x)\,u(x).
Boltzmann relation for the electrons: Since most of the electrons are reflected, their density is given by :n_\mathrm{e}(x) = n_0 \exp\Big(\frac{e\,\varphi(x)}{k_\mathrm{B}T_\mathrm{e}}\Big). '''
Poisson's equation:''' The curvature of the electrostatic potential is related to the net
charge density as follows: :\frac{d^2\varphi(x)}{dx^2} = \frac{e (n_\mathrm{e}(x)-n_\mathrm{i}(x))}{\epsilon_0} . Combining these equations and writing them in terms of the dimensionless potential, position, and ion speed, :\chi(\xi) = -\frac{e\varphi(\xi)}{k_\mathrm{B}T_\mathrm{e}} :\xi = \frac{x}{\lambda_\mathrm{D}} :\mathfrak{M} = \frac{u_\mathrm{o}}{(k_\mathrm{B}T_\mathrm{e}/m_\mathrm{i})^{1/2}} we arrive at the sheath equation: :\chi'' = \left( 1 + \frac{2\chi}{\mathfrak{M}^2} \right)^{-1/2} - e^{-\chi}.
Bohm sheath criterion The sheath equation can be integrated once by multiplying by \chi': :\int_0^\xi \chi' \chi''\,d\xi_1 = \int_0^\xi \left( 1 + \frac{2\chi}{\mathfrak{M}^2} \right)^{-1/2} \chi' \,d\xi_1 - \int_0^\xi e^{-\chi} \chi'\,d\xi_1 At the sheath edge (\xi = 0), we can define the potential to be zero (\chi = 0) and assume that the
electric field is also zero (\chi'=0). With these boundary conditions, the integrations yield :\frac{1}{2}\chi'^2 = \mathfrak{M}^2 \left[ \left( 1 + \frac{2\chi}{\mathfrak{M}^2} \right)^{1/2} - 1 \right] + e^{-\chi} - 1 This is easily rewritten as an integral in closed form, although one that can only be solved numerically. Nevertheless, an important piece of information can be derived analytically. Since the left-hand-side is a square, the right-hand-side must also be non-negative for every value of \chi, in particular for small values. Looking at the Taylor expansion around \chi = 0, we see that the first term that does not vanish is the quadratic one, so that we can require :\frac{1}{2}\chi^2\left( -\frac{1}{\mathfrak{M}^2} + 1 \right) \ge 0, or :\mathfrak{M}^2 \ge 1, or :u_0 \ge (k_\mathrm{B}T_\mathrm{e}/m_\mathrm{i})^{1/2}. This inequality is known as the
Bohm sheath criterion after its discoverer,
David Bohm who discussed it in 1949. If the ions are entering the sheath too slowly, the sheath potential will "eat" its way into the plasma to accelerate them. Ultimately a so-called
pre-sheath will develop with a potential drop on the order of (k_\mathrm{B}T_\mathrm{e}/2e) and a scale determined by the physics of the
ion source (often the same as the dimensions of the plasma). Normally the Bohm criterion will hold with equality, but there are some situations where the ions enter the sheath with
supersonic speed.
Child–Langmuir law Although the sheath equation must generally be integrated numerically, we can find an approximate solution analytically by neglecting the e^{-\chi} term. This amounts to neglecting the
electron density in the sheath, or only analyzing that part of the sheath where there are no electrons. For a "floating" surface, i.e. one that draws no net current from the plasma, this is a useful if rough approximation. For a surface biased strongly negative so that it draws the
ion saturation current, the approximation is very good. It is customary, although not strictly necessary, to further simplify the equation by assuming that 2\chi/\mathfrak{M}^2 is much larger than unity. Then the sheath equation takes on the simple form :\chi'' = \frac{\mathfrak{M}}{(2\chi)^{1/2}}. As before, we multiply by \chi' and integrate to obtain :\frac{1}{2}\chi'^2 = \mathfrak{M} (2\chi)^{1/2}, or :\chi^{-1/4}\chi' = 2^{3/4} \mathfrak{M}^{1/2}. This is easily integrated over ξ to yield :\frac{4}{3}\chi_\mathrm{w}^{3/4} = 2^{3/4} \mathfrak{M}^{1/2} d, where \chi_\mathrm{w} is the (normalized) potential at the wall (relative to the sheath edge), and
d is the thickness of the sheath. Changing back to the variables u_0 and \varphi and noting that the ion current into the wall is J=e\,n_0\,u_0, we have :J = \frac{4}{9} \left(\frac{2e}{m_i}\right)^{1/2} \frac{|\varphi_w|^{3/2}}{4\pi d^2}. This equation is known as '''Child's law
, after Clement D. Child (1868–1933), who first published it in 1911, or as the Child–Langmuir law'
, honoring as well Irving Langmuir, who discovered it independently and published in 1913. It was first used to give the space-charge-limited current in a vacuum diode with electrode spacing d''. It can also be inverted to give the thickness of the Debye sheath as a function of the
voltage drop by setting J=j_\mathrm{ion}^\mathrm{sat}: : d = \frac{2}{3} \left(\frac{2e}{m_\mathrm{i}}\right)^{1/4} \frac{|\varphi_\mathrm{w}|^{3/4}}{2\sqrt{\pi j_\mathrm{ion}^\mathrm{sat}}}. In recent years, the Child-Langmuir (CL) law have been revised as reported in two review papers. == See also ==