Langmuir probe
From Academic Kids

A Langmuir probe is a device named after Nobel Prize winning physicist Irving Langmuir which is used to determine the electron temperature, electron density, and plasma potential. It works by inserting one or more electrodes into a plasma, with a constant or timevarying 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.
Contents 
IV characteristic of the Debye sheath
The beginning of Langmuir probe theory is the IV characteristic of the Debye sheath, that is, the current density flowing to a surface in a plasma as a function of the voltage drop across the sheath. The analysis presented here indicates how the electron temperature, electron density, and plasma potential can be derived from the IV characterisitic. In some situations a more detailed analysis can yield information on the ion density (<math>n_i<math>), the ion temperature <math>T_i<math>, or the electron energy distribution function (EEDF) or <math>f_e(v)<math>.
Ion saturation current density
Consider first an surface biased to a large negative voltage. If the voltage is large enough, essentially all electrons (and any negative ions) will be repelled. The ion velocity will satisfy the Bohm sheath criterion, which is, strictly speaking, an inequality, but which is usually marginally fulfilled. The Bohm criterion in its marginal form says that the ion velocity at the sheath edge is simply the sound speed given by
<math> c_s = \sqrt{k_B(T_e+Z\gamma_iT_i)/m_i}<math>.
The ion temperature term is often neglected, which is justified if the ions are cold. Even if the ions are known to be warm, the ion temperature is usually not known, so it is usually assumed to be simply equal to the electron temperature. In that case, consideration of finite ion temperature only results in a small numerical factor. Z is the (average) charge state of the ions, and <math>\gamma_i<math> is the adiabatic coefficient for the ions. The proper choice of <math>\gamma_i<math> is a matter of some contention. Most analyses use <math>\gamma_i=1<math>, corresponding to isothermal ions, but some kinetic theory suggests that <math>\gamma_i=3<math>, corresponding to one degree of freedom, is more appropriate. For <math>Z=1<math> and <math>T_i=T_e<math>, using the larger value results in the conclusion that the density is <math>\sqrt{2}<math> times smaller. Uncertainties of this magnitude arise several places in the analysis of Langmuir probe data and are very difficult to resolve.
The charge density of the ions depends on the charge state Z, but quasineutrality allows one to write it simply in terms of the electron density as <math>en_e<math>.
Using these results we have the current density to the surface due to the ions. The current density at large negative voltages is due solely to the ions and, except for possible sheath expansion effects, does not depend on the bias voltage, so it is referred to as the ion saturation current density and is given by
<math>j_{sat}^{ion} = en_ec_s<math>
where the plasma parameters, in particular the density, are those at the sheath edge.
Exponential electron current
As the potential drop in the Debye sheath is reduced, the more energetic electrons are able to overcome the potential barrier of the electrostatic sheath. We can model the electrons at the sheath edge with a Boltzmann distribution, i.e.,
<math>f(v_x)\,dv_x \propto e^{\frac{1}{2}m_ev_x^2/k_BT_e}<math>,
except that the high energy tail moving away from the surface is missing, because only the lower energy electrons moving toward the surface are reflected. The higher energy electrons overcome the sheath potential and are absorbed. The mean velocity of the electrons which are able to overcome the potential drop of the sheath is
<math> \langle v_e \rangle = \frac {\int_{v_{e0}}^\infty f(v_x)\,v_x\,dv_x} {\int_{v_{e0}}^\infty f(v_x)\,dv_x} <math>,
where the cutoff velocity for the upper integral is
<math>v_{e0} = \sqrt{2e\Phi_{sh}/m_e}<math>.
<math>\Phi_{sh}<math> is the potential drop across the Debye sheath, that is, the potential at the sheath edge minus the potential of the surface. For a potential drop large compared to the electron temperature, the result is
<math> \langle v_e \rangle = \sqrt{\frac{k_BT_e}{2\pi m_e}}\, e^{e\Phi_{sh}/k_BT_e} <math>.
With this expression we can write the electron contribution to the current to the probe in terms of the ion saturation current as
<math> j_{elec} = j_{ion}^{sat}\sqrt{m_i/2\pi m_e}\, e^{e\Phi_{sh}/k_BT_e} <math>,
valid as long as the electron current is not more than two or three times the ions current.
Floating potential
The total current, of course, is the sum of the ion and electron currents:
<math> j = j_{ion}^{sat} \left( 1 + \sqrt{m_i/2\pi m_e}\,e^{e\Phi_{sh}/k_BT_e} \right) <math>.
We are using the convention that current from the surface into the plasma is positive. An interesting and practical question is the potential of a surface to which no net current flows. It is easily seen from the above equation that
<math>\Phi_{fl} = (k_BT_e/e)\,(1/2)\ln(m_i/2\pi m_e)<math>.
If we introduce the reduced ion mass <math>\mu_i=m_i/m_p<math>, we can write
<math> \Phi_{fl} = (k_BT_e/e)\, ( 2.8 + 0.5\ln \mu_i ) <math>
Since the floating potential is the experimentally accessible quantity, the current (below electron saturation) is usually written as
<math> j = j_{ion}^{sat} \left( 1 + \,e^{e(V_{pr}V_{fl})/k_BT_e} \right) <math>.
Electron saturation current
When the electrode potential is equal to or greater than the plasma potential, then there is no longer a sheath to reflect electrons, so we may expect the electron current to saturate. Using the Boltzmann expression for the mean electron velocity given above with <math>v_{e0} = 0<math> and setting the ion current to zero, the electron saturation current density would be
<math> j_{elec}^{sat} = j_{ion}^{sat}\sqrt{m_i/\pi m_e} = j_{ion}^{sat} \left( 24.2 * \sqrt{\mu_i} \right) <math>
This derivation is not rigorous and may be missing important physical effects. In practice, it is often difficult and usually considered uninformative to measure the electron saturation current experimentally. When it is measured, it is found to be highly variable and generally much (a factor of three or more) lower than the value given above. Often a clear saturation is not seen at all. Understanding electron saturation is one of the most important outstanding problem of Langmuir probe theory.
Effects of the bulk plasma
Debye sheath theory explains the basic behavior of Langmuir probes, but is not complete. Merely inserting an object like a probe into a plasma changes the density, temperature, and potential at the sheath edge and perhaps everywhere. Changing the voltage on the probe will also, in general, change various palsma parameters. Such effects are less well understood than sheath physics, but they can at least in some cases be roughly accounted.
Presheath
We have noted that the Bohm criterion requires the ions to enter the Debye sheath at the sound speed. The potential drop that accelerates them to this speed is called the presheath. It has a spatial scale that depnds on the physics of the ion source but which is large compared to the Debye length and often of the order of the plasma dimensions. The magnitude of the potential drop is equal to (at least)
<math> \Phi_{pre} = \frac{\frac{1}{2}m_ic_s^2}{Ze} = k_B(T_e+Z\gamma_iT_i)/(2Ze) <math>
The acceleration of the ions also entails a decrease in the density, usually by a factor of about 2 depending on the details.
Resistivity
Collisions between ions and electrons will also affect the IV characteristic of a Langmuir probe. When an electrode is biased to any voltage other than the floating potential, the current it draws must pass through the plasma, which has a finite resistivity. The resistivity and current path can be calculated with relative ease in an unmagnetized plasma. In a magnetized plasma, the problem is much more difficult. In either case, the effect is to add a voltage drop proportional to the current drawn, which shears the characteristic. The deviation from an exponential function is usually not possible to observe directly, so that the flattening of the characteristic is usually misinterpreted as a larger plasma temperature. Looking at it from the other side, any measured IV characteristic can be interpreted as a hot plasma, where most of the voltage is dropped in the Debye sheath, or as a cold plasma, where most of the voltage is dropped in the bulk plasma. Without quantitative modeling of the bulk resistivity, Langmuir probes can only give an upper limit on the electron temperature.
Sheath expansion
It is not enough to know the current density as a function of bias voltage, since it is the absolute current which is measured. In an unmagnetized plasma, the currentcollecting area is usually taken to be the exposed surface area of the electrode. In a magnetized plasma, the projected area is taken, that is, the area of the electrode as viewed along the magnetic field. If the electrode is not shadowed by a wall or other nearby object, then the area must be doubled to account for current coming along the field from both sides. If the electrode dimensions are not small in comparison to the Debye length, then the size of the electrode is effectively increased in all directions by the sheath thickness. In a magnetized plasma, the electrode is sometimes assumed to be increased in a similar way by the ion Larmor radius. It is certainly true that the finite Larmor radius allows some ions to reach the electrode that would have otherwise gone past it, but the details of the effect have never been calculated in a fully selfconsistent way.
If we refer to the probe area including these effects as <math>A_{eff}<math> (which may be a function of the bias voltage) and make the assumptions
 <math>T_i=T_e<math>,
 <math>Z=1<math>
 <math>\gamma_i=3<math>, and
 <math>n_{e,sh}=0.5\,n_e<math>,
and ignore the effects of
 bulk resistivity, and
 electron saturation,
then the IV characteristic becomes
<math> I = I_{ion}^{sat}e^{(V_{pr}V_{fl})/(k_BT_e/e)} <math>,
where
<math> I_{ion}^{sat} = en_e\sqrt{k_BT_e/m_i}\,A_{eff} <math>.
Magnetized plasmas
The theory of Langmuir probes is much more complex when the plasma is magnetized. The simplest extension of the unmagnetized case is simply to use the projected area rather than the surface area of the electrode. For a long cylinder far from other surfaces, this reduces the effective area by a factor of π/2 = 1.57. As mentioned before, it might be necessary to increase the radius by about the thermal ion Larmor radius, but not above the effective area for the unmagnetized case.
The use of the projected area seems to be closely tied with the existence of a magnetic sheath. Its scale is the ion Larmor radius at the sound speed, which is normally between the scales of the Debye sheath and the presheath. The Bohm criterion for ions entering the magnetic sheath applies to the motion along the field, while at the entrance to the Debye sheath it applies to the motion normal to the surface. This results in a reduction of the density by the sine of the angle between the field and the surface. The associated increase in the Debye length must be taken into account when considering ion nonsaturation due to sheath effects.
Especially interesting and difficult to understand is the role of crossfield currents. Naively, one would expect the current to flow parallel to the magnetic field along a flux tube. In many geometries, this flux tube will end at a surface in a distant part of the device, and this spot should itself exhibit an IV characteristic. The net result would be the measurement of a doubleprobe characteristic, in other words, electron saturation current equal to the ion saturation current.
When this picture is considered in detail, it is seen that the flux tube must charge up and the surrounding plasma must spin around it. The flow of current into or out of the flux tube must be associated with a force that slow down this spinning. Candidate forces are viscosity, friction with neutrals, and inertial forces associated with plasma flows, either steady or fluctuating. It is not known which force is strongest in practice, and in fact it is generally difficult to find any force that is powerful enough to explain the characteristics actually measured.
It is also likely that the magnetic field plays a decisive role in determining the level of electron saturation, but no quantitative theory is as yet available.
Electrode configurations
Once one has a theory of the IV characteristic of an electrode, one can proceed to measure it and then fit the data with the theoretical curve to extract the plasma parameters. The straightforward was to do this is to sweep the voltage on a single electrode, but, for a number of reasons, configurations using multiple electrodes and/or exploring only a part of the characteristic are used in practice.
Single probe
The most straightforward way to measure the IV characteristic of a plasma is with a single probe, consisting of one electrode biased with a voltage ramp relative to the vessel. The advantages are simplicity of the electrode and redundancy of information, i.e. one can check whether the IV characteristic has the expected form. Potentially additional information can be extracted from details of the characteristic. The disadvantages are more complex biasing and measurement electronics and a poor time resolution. If fluctuations are present (as they always are) and the sweep is slower than the fluctuation frequency (as it usually is), then the IV is is the average current as a function of voltage, which may result in systematic errors if it is analyzed as though it were an instantaneous IV. The ideal situation is to sweep the voltage at a frequenc above the fluctuation frequency but still below the ion cyclotron frequency. This, however, requires sophisticated electronics and a great deal of care.
Double probe
An electrode can be biased relative to a second electrode, rather than to the vessel. The theory is similar to that of a single probe, except that the current is limited to the ion saturation current for both positive and negative voltages. In particular, if <math>V_{bias}<math> is the voltage applied between two identical electrodes, we can write
<math> I = I_{ion}^{sat} \left( 1 + \,e^{e(V_2V_{fl})/k_BT_e} \right) = I_{ion}^{sat} \left( 1 + \,e^{e(V_1V_{fl})/k_BT_e} \right) <math>,
which can be rewritten using <math>V_{bias}=V_2V_1<math> as
<math> I = I_{ion}^{sat} \tanh\left( \frac{1}{2}\,\frac{eV_{bias}}{k_BT_e} \right) <math>.
One advantage of the double probe is that that neither electrode is ever very far above floating, so the theoretical uncertainties at large electron currents are avoided. If it is desired to sample more of the exponential electron portion if the characteristic, and asymmetric double probe may be used, with one electrode larger than the other. The characteristic in this case is still a hyperbolic tangent, but shifted vertically. Another advantage is that there is no reference to the vessel, so it is too some extent immune to the disturbances in an rf plasma. On the other hand, it shares the limitations of a single probe concerning complicated electronics and poor time resolution. In addition, the second electrode not only complicates the system, but it makes it susceptible to disturbance by gradients in the plasma.
Triple probe
An elegant electrode configuration is the triple probe, consisting of two electrodes biased with a fixed voltage and a third which is floating. The bias voltage is chosen to be a few time the electron temperature so that the negative electrode draws the ion saturation current, which, like the floating potential, is directly measured. The positive tip draws a positive current equal in magnitude to the (normally negative) ion saturation current, so we can write
<math> I_{ion}^{sat} = I_{ion}^{sat} \left( 1 + \,e^{e(V_{pos}V_{fl})/k_BT_e} \right) <math>,
or
<math> (V_{pos}V_{fl}) = \ln 2 \, (k_BT_e/e) <math>.
That is, the voltage difference between the positive and floating electrodes is proportional to the electron temperature. (This was especially important in the sixties and seventies before sophisticated data processing became widely available.)
More sophisticated analysis of triple probe data can take into account such factors as incomplete saturation, nonsaturation, unequal areas.
Triple probes have the advantage of simple biasing electronics (no sweeping required), simple data analysis, excellent time resolution, and insensitivity to potential fluctuations (whether imposed by an rf source or inherent fluctuations). Like double probes, they are sensitive to gradients in plasma parameters.
Special arrangements
Arrangements with four (tetra probe) or five (penta probe) have sometimes been used, but the advantage over triple probes has never been entirely convincing.
A pinplate probe consists of a small electrode directly in front of a large electrode, the idea being that the voltage sweep of the large probe can perturb the plasma potential at the sheath edge and thereby aggravate the difficulty of interpreting the IV characteristic. The floating potential of the small electrode can be used to correct for changes in potential at the sheath edge of the large probe. Experimental results from this arrangement look promising, but experimental complexity and residual difficulties in the interpretation have prevented this configuration from becoming standard.
Various geometries have been proposed for use as ion temperature probes, for example, two cylindrical tips that rotate past each other im a magnetized plasma. Since shadowing effects depend on the ion Larmor radius, the results can be interpreted in terms of ion temperature. The ion temperature is an important quantity that is very difficult to measure. Unfortunately, it is also very difficult to analyze such probes in a fully selfconsistent way.
To measure fluctuations in plasma parameters, arrays of electrodes are use, usually one but occasionally twodimensional. A typical array has a spacing of 1 mm and a total of 16 or 32 electrodes. A simpler arrangement to measure fluctuations is a negativelybiased electrode flanked by two floating electrodes. The ionsaturation current is taken as a surrogate for the density and the floating potential as a surrogate for the plasma potential. This allows a rough measurement of the turbulent particle flux
<math> \Phi_{turb} = \langle \tilde{n}_e \tilde{v}_{E\times B} \rangle \propto \langle \tilde{I}_{ion}^{sat} ( \tilde{V}_{fl,2}  \tilde{V}_{fl,1} ) \rangle <math>
Practical considerations
For laboratory and technical plasmas, the electrodes are most commonly tungsten wires several thousandths of an inch thick, because they have a high melting point but can be made small enough not to perturb the plasma. For fusion plasmas, graphite electrodes with dimensions from 1 to 10 mm are usually used because they can withstand the highest power loads (They not only withstand high temperature, when they do overheat, they sublimate rather than melt.), and the atoms that do get into the plasma result in less radiation due to the low atomic number of carbon. The electrode surface exposed to the plasma must be defined, e.g. by insulating all but the tip of a wire electrode. If there can be significant deposition of conducting materials (metals or graphite), then the insulator should be separated from the electrode by a meander to prevent shortcircuiting.
In a magnetized plasma, it appears to be best to choose a probe size a few times larger than the ion Larmor radius. A point of contention is whether it is better to use proud probes, where the angle between the magnetic field and the surface at least 15°, or flushmounted probes, which are imbedded in the plasmafacing components and generally have an angle of 1 to 5 °. Many plasma physicists feel more comfortable with proud probes, which have a longer tradition and possibly are less perturbed by electron saturation effects, although this is disputed. Fluchmounted probes, on the other hand, being part of the wall, are less pertabative. Knowledge of the field angle is necessary with proud probes to determine the fluxes to the wall, whereas it is necessary with flushmounted probes to determine the density.
In very hot and dense plasmas, as found in fusion research, it is often necessary to limit the thermal load to the probe by limiting the exposure time. A reciprocating probe is mounted on an arm that is moved into and back out of the plasma, usually in about one second by means of either a pneumatic drive or an electromagnetic drive using the ambient magnetic field. Popup probes are similar, but the electodes rest behind a shield and are only moved the few millimeters necessary to bring them into the plasma near the wall.
A Langmuir probe can be purchased off the shelf for on the order of 15,000 U.S. Dollars, or they can be built by an experienced researcher and/or technician. When working at frequencies under 100 MHz, it is advisable to use blocking filters, and take necessary grounding precautions.