Debye sheath

Last updated

The Debye sheath (also electrostatic sheath) is a layer in a plasma which has a greater density of positive ions, and hence an overall excess positive charge, that balances an opposite negative charge on the surface of a material with which it is in contact. The thickness of such a layer is several Debye lengths thick, a value whose size depends on various characteristics of plasma (e.g. temperature, density, etc.).

Contents

A Debye sheath arises in a plasma because the electrons usually have a temperature on the order of magnitude or greater than that of the ions and are much lighter. Consequently, they are faster than the ions by at least a factor of . At the interface to a material surface, therefore, the electrons will fly out of the plasma, charging the surface negative relative to the bulk plasma. Due to Debye shielding, the scale length of the transition region will be the Debye length . As the potential increases, more and more electrons are reflected by the sheath potential. An equilibrium is finally reached when the potential difference is a few times the electron temperature.

The Debye sheath is the transition from a plasma to a solid surface. Similar physics is involved between two plasma regions that have different characteristics; the transition between these regions is known as a double layer, and features one positive, and one negative layer.

Description

Positive ion sheaths around grid wires in a thermionic gas tube, where [?] represents a positive charge (not to scale) (After Langmuir, 1929) Plasma-sheath.svg
Positive ion sheaths around grid wires in a thermionic gas tube, where represents a positive charge (not to scale) (After Langmuir, 1929)

Sheaths were first described by American physicist Irving Langmuir. In 1923 he wrote:

"Electrons are repelled from the negative electrode while positive ions are drawn towards it. Around each negative electrode there is thus a sheath of definite thickness containing only positive ions and neutral atoms. [..] Electrons are reflected from the outside surface of the sheath while all positive ions which reach the sheath are attracted to the electrode. [..] it follows directly that no change occurs in the positive ion current reaching the electrode. The electrode is in fact perfectly screened from the discharge by the positive ion sheath, and its potential cannot influence the phenomena occurring in the arc, nor the current flowing to the electrode." [1]

Langmuir and co-author Albert W. Hull further described a sheath formed in a thermionic valve:

"Figure 1 shows graphically the condition that exists in such a tube containing mercury vapor. The space between filament and plate is filled with a mixture of electrons and positive ions, in nearly equal numbers, to which has been given the name "plasma". A wire immersed in the plasma, at zero potential with respect to it, will absorb every ion and electron that strikes it. Since the electrons move about 600 times as fast as the ions, 600 times as many electrons will strike the wire as ions. If the wire is insulated it must assume such a negative potential that it receives equal numbers of electrons and ions, that is, such a potential that it repels all but 1 in 600 of the electrons headed for it."
"Suppose that this wire, which we may take to be part of a grid, is made still more negative with a view to controlling the current through the tube. It will now repel all the electrons headed for it, but will receive all the positive ions that fly toward it. There will thus be a region around the wire which contains positive ions and no electrons, as shown diagrammatically in Fig. 1. The ions are accelerated as they approach the negative wire, and there will exist a potential gradient in this sheath, as we may call it, of positive ions, such that the potential is less and less negative as we recede from the wire, and at a certain distance is equal to the potential of the plasma. This distance we define as the boundary of the sheath. Beyond this distance there is no effect due to the potential of the wire." [2]

Mathematical treatment

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 entering the sheath with a velocity , having charge opposite to the electron, conservation of energy in the sheath potential requires

,

where is the charge of the electron taken positively, i.e. x .

Ion continuity: In the steady state, the ions do not build up anywhere, so the flux is everywhere the same:

.

Boltzmann relation for the electrons: Since most of the electrons are reflected, their density is given by

.

Poisson's equation: The curvature of the electrostatic potential is related to the net charge density as follows:

.

Combining these equations and writing them in terms of the dimensionless potential, position, and ion speed,

we arrive at the sheath equation:

.

The Bohm sheath criterion

The sheath equation can be integrated once by multiplying by :

At the sheath edge (), we can define the potential to be zero () and assume that the electric field is also zero (). With these boundary conditions, the integrations yield

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 , in particular for small values. Looking at the Taylor expansion around , we see that the first term that does not vanish is the quadratic one, so that we can require

,

or

,

or

.

This inequality is known as the Bohm sheath criterion after its discoverer, David Bohm. 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 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.

The Child–Langmuir law

Although the sheath equation must generally be integrated numerically, we can find an approximate solution analytically by neglecting the 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 is much larger than unity. Then the sheath equation takes on the simple form

.

As before, we multiply by and integrate to obtain

,

or

.

This is easily integrated over ξ to yield

,

where 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 and and noting that the ion current into the wall is , we have

.

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 :

.

In recent years, the Child-Langmuir (CL) law have been revised as reported in two review papers. [3] [4]

See also

Footnotes

  1. Langmuir, Irving, "Positive Ion Currents from the Positive Column of Mercury Arcs" (1923) Science, Volume 58, Issue 1502, pp. 290-291
  2. Albert W. Hull and Irving Langmuir, "Control of an Arc Discharge by Means of a Grid", Proc Natl Acad Sci USA. 1929 March 15; 15(3): 218–225
  3. P. Zhang, A. Valfells, L. K. Ang, J. W. Luginsland and Y. Y. Lau (2017). "100 years of the physics of diodes". Applied Physics Reviews . 4 (1): 011304. Bibcode:2017ApPRv...4a1304Z. doi: 10.1063/1.4978231 .{{cite journal}}: CS1 maint: multiple names: authors list (link)
  4. P Zhang, Y. S. Ang, A. L. Garner, A. Valfells, J. L. Luginsland, and L. K. Ang (2021). "Space–charge limited current in nanodiodes: Ballistic, collisional, and dynamical effects". Journal of Applied Physics . 129 (10): 100902. Bibcode:2021JAP...129j0902Z. doi:10.1063/5.0042355. hdl: 20.500.11815/2643 . S2CID   233643434.{{cite journal}}: CS1 maint: multiple names: authors list (link)

Related Research Articles

<span class="mw-page-title-main">Permittivity</span> Measure of the electric polarizability of a dielectric material

In electromagnetism, the absolute permittivity, often simply called permittivity and denoted by the Greek letter ε (epsilon), is a measure of the electric polarizability of a dielectric material. A material with high permittivity polarizes more in response to an applied electric field than a material with low permittivity, thereby storing more energy in the material. In electrostatics, the permittivity plays an important role in determining the capacitance of a capacitor.

A Langmuir probe is a device used to determine the electron temperature, electron density, and electric potential of a plasma. It works by inserting one or more electrodes into a plasma, with a constant or time-varying electric potential between the various electrodes or between them and the surrounding vessel. The measured currents and potentials in this system allow the determination of the physical properties of the plasma.

In quantum chemistry and molecular physics, the Born–Oppenheimer (BO) approximation is the best-known mathematical approximation in molecular dynamics. Specifically, it is the assumption that the wave functions of atomic nuclei and electrons in a molecule can be treated separately, based on the fact that the nuclei are much heavier than the electrons. Due to the larger relative mass of a nucleus compared to an electron, the coordinates of the nuclei in a system are approximated as fixed, while the coordinates of the electrons are dynamic. The approach is named after Max Born and his 23-year-old graduate student J. Robert Oppenheimer, the latter of whom proposed it in 1927 during a period of intense ferment in the development of quantum mechanics.

In physics, screening is the damping of electric fields caused by the presence of mobile charge carriers. It is an important part of the behavior of charge-carrying fluids, such as ionized gases, electrolytes, and charge carriers in electronic conductors . In a fluid, with a given permittivity ε, composed of electrically charged constituent particles, each pair of particles interact through the Coulomb force as

Plasma oscillations, also known as Langmuir waves, are rapid oscillations of the electron density in conducting media such as plasmas or metals in the ultraviolet region. The oscillations can be described as an instability in the dielectric function of a free electron gas. The frequency depends only weakly on the wavelength of the oscillation. The quasiparticle resulting from the quantization of these oscillations is the plasmon.

In plasmas and electrolytes, the Debye length, is a measure of a charge carrier's net electrostatic effect in a solution and how far its electrostatic effect persists. With each Debye length the charges are increasingly electrically screened and the electric potential decreases in magnitude by 1/e. A Debye sphere is a volume whose radius is the Debye length. Debye length is an important parameter in plasma physics, electrolytes, and colloids. The corresponding Debye screening wave vector for particles of density , charge at a temperature is given by in Gaussian units. Expressions in MKS units will be given below. The analogous quantities at very low temperatures are known as the Thomas–Fermi length and the Thomas–Fermi wave vector. They are of interest in describing the behaviour of electrons in metals at room temperature.

The Vlasov equation is a differential equation describing time evolution of the distribution function of plasma consisting of charged particles with long-range interaction, such as the Coulomb interaction. The equation was first suggested for the description of plasma by Anatoly Vlasov in 1938 and later discussed by him in detail in a monograph.

<span class="mw-page-title-main">Debye–Hückel equation</span> Electrochemical equation

The chemists Peter Debye and Erich Hückel noticed that solutions that contain ionic solutes do not behave ideally even at very low concentrations. So, while the concentration of the solutes is fundamental to the calculation of the dynamics of a solution, they theorized that an extra factor that they termed gamma is necessary to the calculation of the activities of the solution. Hence they developed the Debye–Hückel equation and Debye–Hückel limiting law. The activity is only proportional to the concentration and is altered by a factor known as the activity coefficient . This factor takes into account the interaction energy of ions in solution.

One of the guiding principles in modern chemical dynamics and spectroscopy is that the motion of the nuclei in a molecule is slow compared to that of its electrons. This is justified by the large disparity between the mass of an electron, and the typical mass of a nucleus and leads to the Born–Oppenheimer approximation and the idea that the structure and dynamics of a chemical species are largely determined by nuclear motion on potential energy surfaces.

<span class="mw-page-title-main">Electrodynamic tether</span> Long conducting wires which can act as electrical motors or generators

Electrodynamic tethers (EDTs) are long conducting wires, such as one deployed from a tether satellite, which can operate on electromagnetic principles as generators, by converting their kinetic energy to electrical energy, or as motors, converting electrical energy to kinetic energy. Electric potential is generated across a conductive tether by its motion through a planet's magnetic field.

In physics, the Einstein relation is a previously unexpected connection revealed independently by William Sutherland in 1904, Albert Einstein in 1905, and by Marian Smoluchowski in 1906 in their works on Brownian motion. The more general form of the equation in the classical case is

A dusty plasma is a plasma containing micrometer (10−6) to nanometer (10−9) sized particles suspended in it. Dust particles are charged and the plasma and particles behave as a plasma. Dust particles may form larger particles resulting in "grain plasmas". Due to the additional complexity of studying plasmas with charged dust particles, dusty plasmas are also known as complex plasmas.

Plasma parameters define various characteristics of a plasma, an electrically conductive collection of charged and neutral particles of various species that responds collectively to electromagnetic forces. Such particle systems can be studied statistically, i.e., their behaviour can be described based on a limited number of global parameters instead of tracking each particle separately.

The Poisson–Boltzmann equation describes the distribution of the electric potential in solution in the direction normal to a charged surface. This distribution is important to determine how the electrostatic interactions will affect the molecules in solution. The Poisson–Boltzmann equation is derived via mean-field assumptions. From the Poisson–Boltzmann equation many other equations have been derived with a number of different assumptions.

The Orr–Sommerfeld equation, in fluid dynamics, is an eigenvalue equation describing the linear two-dimensional modes of disturbance to a viscous parallel flow. The solution to the Navier–Stokes equations for a parallel, laminar flow can become unstable if certain conditions on the flow are satisfied, and the Orr–Sommerfeld equation determines precisely what the conditions for hydrodynamic stability are.

In classical mechanics, a Liouville dynamical system is an exactly solvable dynamical system in which the kinetic energy T and potential energy V can be expressed in terms of the s generalized coordinates q as follows:

In mathematics, the spectral theory of ordinary differential equations is the part of spectral theory concerned with the determination of the spectrum and eigenfunction expansion associated with a linear ordinary differential equation. In his dissertation, Hermann Weyl generalized the classical Sturm–Liouville theory on a finite closed interval to second order differential operators with singularities at the endpoints of the interval, possibly semi-infinite or infinite. Unlike the classical case, the spectrum may no longer consist of just a countable set of eigenvalues, but may also contain a continuous part. In this case the eigenfunction expansion involves an integral over the continuous part with respect to a spectral measure, given by the Titchmarsh–Kodaira formula. The theory was put in its final simplified form for singular differential equations of even degree by Kodaira and others, using von Neumann's spectral theorem. It has had important applications in quantum mechanics, operator theory and harmonic analysis on semisimple Lie groups.

In mathematics, the Plancherel theorem for spherical functions is an important result in the representation theory of semisimple Lie groups, due in its final form to Harish-Chandra. It is a natural generalisation in non-commutative harmonic analysis of the Plancherel formula and Fourier inversion formula in the representation theory of the group of real numbers in classical harmonic analysis and has a similarly close interconnection with the theory of differential equations. It is the special case for zonal spherical functions of the general Plancherel theorem for semisimple Lie groups, also proved by Harish-Chandra. The Plancherel theorem gives the eigenfunction expansion of radial functions for the Laplacian operator on the associated symmetric space X; it also gives the direct integral decomposition into irreducible representations of the regular representation on L2(X). In the case of hyperbolic space, these expansions were known from prior results of Mehler, Weyl and Fock.

The optical metric was defined by German theoretical physicist Walter Gordon in 1923 to study the geometrical optics in curved space-time filled with moving dielectric materials.

Taylor–Maccoll flow refers to the steady flow behind a conical shock wave that is attached to a solid cone. The flow is named after G. I. Taylor and J. W. Maccoll, whom described the flow in 1933, guided by an earlier work of Theodore von Kármán.