When electrodes are biased above the plasma potential, electrons accelerated through the associated electron sheath can dramatically increase the ionization rate of neutrals near the electrode surface. It has previously been observed that if the ionization rate is great enough, a double layer separates a luminous high-potential plasma attached to the electrode surface (called an anode spot or fireball) from the bulk plasma. Here, results of the first 2D particle-in-cell simulations of anode spot formation are presented along with a theoretical model describing the formation process. It is found that ionization leads to the build-up of an ion-rich layer adjacent to the electrode, forming a narrow potential well near the electrode surface that traps electrons born from ionization. It is shown that anode spot onset occurs when a quasineutral region is established in the potential well and the density in this region becomes large enough to violate the steady-state Langmuir condition, which is a balance between electron and ion fluxes across the double layer. A model for steady-state properties of the anode spot is also presented, which predicts values for the anode spot size, double layer potential drop, and form of the sheath at the electrode by considering particle, power, and current balance. These predictions are found to be consistent with the presented simulation and previous experiments.
Low pressure anode spots are a discharge phenomenon that occurs near positively biased electrodes at pressures less than a few hundred millitorr.1 When the potential across an electron sheath is large, electrons can gain enough energy to ionize neutral species by electron impact ionization. If the ionization rate is sufficient, a positive space charge layer will develop at the electrode next to a negative space charge layer adjacent to the plasma. These adjacent positive and negative space charge layers are a type of anode double layer2 called an anode glow.3 The anode glow potential profile is sketched in Fig. 1 and is shown in the photograph in Fig. 2(a). In experiments with sufficient ionization, it has been observed that an increase in bias leads to the rapid formation of an anode spot.3–6 The main characteristic of an anode spot is a double layer detached from the electrode separating a luminous high-potential plasma from the bulk plasma. An anode spot is shown in the photograph in Fig. 2(b). Potential profiles typical of these situations are shown in Fig. 1. Anode spots are just one of the several forms of sheath structure that are possible near positively biased electrodes4 and are particularly of interest due to their application in electron temperature control,7 dust confinement,8 plasma contactors,9 as an ion source,10 and as a platform to study plasma self organization.11,12 In this paper, the anode spot onset is studied using particle-in-cell (PIC) simulations for the first time and a theory for their onset and equilibrium properties is presented.
Stationary double layers are expected to satisfy a flux density balance criterion called the Langmuir condition13
where Γe = neVe and Γi = niVi are the electron and ion flux densities and Ve and Vi are the electron and ion flow velocities at the double layer sheath edge. In 1972, Block14 noted that the Langmuir condition is satisfied in the reference frame of the double layer. Two decades later, Song, Merlino, and D'Angelo15 derived a modified Langmuir condition by considering the momentum and continuity equation in the frame of a double layer moving with velocity UDL with respect to the lab frame. The condition in the lab frame is
where nHi is the density of the plasma at the double layer's high potential side. Here, the double layer motion results from an imbalance in the electron and ion fluxes crossing the double layer. Using this newly derived condition, the stability of the anode spot steady-state was studied.16
In this paper, the modified Langmuir condition is used to derive a density criterion between the high and low potential sides of the double layer required for spot onset. This density criterion is then used to form a model for the critical electrode bias required for spot onset. When a flux imbalance in Eq. (1) is present, the double layer is predicted to move away from the electrode forming an anode spot. The model, informed by particle-in-cell (PIC) simulations, also describes a mechanism for establishing a quasineutral region with ni ≈ ne and E ≈ 0 at the high potential side. Once sufficient ion density has built up next to the electrode, a potential well for electrons forms. This allows trapping of electrons near the electrode and the formation of a quasineutral spot plasma. Similar potential wells have been predicted for the high pressure anode spot phenomenon.17 Previous work had suggested that an instability was needed for spot onset.18 The model presented here suggests a significantly different picture that does not include instability as a feature of spot onset.
This paper also presents a model for steady-state properties of the anode spot. In the model, the balance of particle creation and loss, power, and current loss within the spot are used to determine the anode spot diameter, double layer potential drop, and form of the sheath between the electrode and spot plasma as a function of electron temperature. This model couples these different properties through the energy dependence of the electron impact ionization cross section of neutrals. Prior descriptions of the anode spot length scale L have used equations of the form
where σI is the electron impact ionization cross section and nn is the neutral density.3,19 In these descriptions, a constant estimate of the cross section was assumed instead of considering its energy dependence.
This paper presents the first study of anode spots using particle-in-cell simulations. The simulations presented in Sec. II are used to formulate a theory for the spot onset in Sec. III A. Following this, the theory of Secs. III B and III C predicts steady-state properties of the anode spot. Finally, properties of the ion presheath leading up to the high potential entrance of the double layer are calculated in Sec. III D, showing that, under typical low temperature plasma conditions, the presheath length is determined by the electron impact ionization rate. Concluding remarks are given in Sec. IV.
A. Simulation setup
Simulations were performed using the electrostatic PIC code Aleph.20 Aleph solves for the electric field on an unstructured mesh, and includes three velocity components for particles. The simulations incorporated collisions between charged and neutral particles using the Direct Simulation Monte-Carlo (DSMC) method.21
The 2D simulation domain shown in Fig. 3 was chosen to resemble experiments,22 though at reduced spatial dimensions of 7.5 cm × 7.5 cm. The boundary conditions for particles were one reflecting boundary (left wall) and three absorbing boundaries (right, top, and bottom walls). The absorbing walls had a ϕ = 0 V Dirichlet boundary condition for the electric field except for a small 0.25 cm electrode embedded in the lower boundary. The electrode bias was increased linearly from 0 V at t = 0 to 40 V at t = 9 μs. At t = 9 μs, the electrode bias was stepped to a fixed value and held constant for the remainder of the simulation (see Fig. 3). The final electrode bias was varied over several simulations to determine that required for anode spot formation. The reflecting boundary had a Neumann boundary condition, which, in combination with the reflecting boundary condition for particles, allows the simulation to represent half of a symmetric 15 cm × 7.5 cm domain. The domain was discretized using an unstructured triangular mesh with an element size of approximately 0.0116 cm that resolved the Debye length across the domain throughout the duration of the simulation. A time step of 10 ps was chosen to resolve a Courant-Friedrichs-Lewy-like (CFL) condition for all particles, preventing particles from crossing more than one cell length before a velocity update and ensuring that particle trajectories were resolved.
A helium plasma was continuously generated at a rate of 4.7 × 109 cm−3 μs−1 in a 0.25 cm by 6.5 cm source region 6.25 cm above the electrode. This was designed to model a source chamber in the experiments described in Ref. 22. Electrons and ions were sourced with temperatures of 4 eV and 86 meV, respectively. Particles left the source region by ambipolar diffusion and filled the domain. After expansion of the plasma from the source region, the electron temperature in the bulk plasma was approximately 2.4 eV. A neutral helium background with density 7 × 1015 cm−3 at 24 meV (200 mTorr at 273 K) was present in the simulation. Electrons, ions, and neutral particles had weights of 2000, 8000, and 2 × 1010, respectively. Two populations of electrons were tracked to better understand the effects of ionization within the anode spot. Electrons which have participated in an ionization event are denoted and those which have not are denoted (bulk electrons). The included ionization interactions for these species are
Elastic electron neutral collisions with the helium background were also included. Ion-neutral elastic and charge exchange collisions are not expected to play a significant role in the qualitative behavior of the anode spot since this region is dominated by ionization (see Sec. III D for justification of this assumption). The cross sections used for collisions in the simulations and for calculations in Secs. III A, III C, and III D were obtained from the Phelps database on LXcat.23
Differences between simulated and actual anode spots can be expected since a 2D simulation is used to represent a 3D anode spot. However, the qualitative behavior of the anode spot is expected to be the same since the dimensionality does not change the ionization or particle and current loss mechanisms. For the simulation, a minimum of two dimensions is needed to represent the problem since electrodes cannot be biased above the plasma potential in 1D without current emission from the boundary. Three dimensional simulations were not attempted since the 2D case was already computationally expensive. Nevertheless, changes in the rates of ionization and particle or current loss may be expected between 2D and 3D representations. In particular, this may lead to differences in the plasma density and average electron temperature within the anode spot.
B. Anode spot simulation
Simulations were carried out for final electrode biases of 40 V, 45 V, 46 V, 47 V, 48 V, and 50 V to determine the critical electrode bias for spot onset. For each simulation, the plasma potential at t = 9 μs, just before the electrode bias was increased, was 6.4 V. Anode spots were observed for the simulation with a 48 V electrode bias, but not with a 47 V bias, indicating that the critical bias with respect to the plasma potential is in the range 40.6–41.6 V. In this section, the simulation with a 50 V bias will be used to study the spot onset.
Figure 4 shows profiles of the ion, bulk electron, and ionization electron densities along with the electrostatic potential and electric field along the symmetry axis (reflecting boundary) in time throughout the 50 V simulation. Early in the simulation, particles were observed to fill the domain by ambipolar diffusion from the source region. Initially, a large portion of the lower half of the domain was electron rich. After approximately 4 μs, the ions had enough time to traverse the domain and encounter the electric field at the electrode. Once these ions were reflected by the field, an electron sheath formed. This can be seen in the ion density and electric field in Fig. 4. Between 5 μs and 9 μs, waves propagating toward the electron sheath can be seen in the ion density and electric field. These are ion acoustic waves excited by the differential flow between ions in the plasma and electrons accelerating to their thermal speed in the electron presheath.24
At 9 μs, the electrode bias was set to 50 V. Figure 4 shows a dramatic increase in ion density within the sheath immediately following this jump in applied bias. Although the mean free path for electron impact ionization is much longer than the sheath scale (∼50 cm for 26 eV electrons), significant build-up of positive space charge results from the large incoming flux of energetic electrons and the relative immobility of ions compared with electrons. The effect of this increase can be seen by comparing the first two panels of Fig. 5. The increase in ion density is not immediately accompanied by an increase in neI since the sheath field quickly accelerates ionization-born electrons to the electrode. At t = 9.5 μs, the build-up of ions within the sheath has resulted in a change in concavity of the potential profile. The panel of Fig. 5 marked t = 10.2 μs shows that a further increase in ion density results in a change in the sign of the electric field, which is indicated by the maximum in the potential just off of the electrode surface accompanied by a small increase in neI. Low energy electrons resulting from electron impact ionization interactions have temperature of approximately 1 eV and are trapped at the maximum by the ∼1 V drop to the electrode. Once this electron trapping occurs, the electron density increases rapidly until quasineutrality is established. This is followed by expansion leading to the formation of the anode spot plasma shown in the densities and electric field plotted in Fig. 4 and the panel marked t = 11 μs in Fig. 5.
By t = 13 μs, expansion of the anode spot slows and the double layer potential settles at approximately 27.2 V, which is 2.6 V above the potential needed for electron impact ionization of helium. This is similar to experiments which have measured the double layer potential to be a few volts above the ionization energy of the neutral gas.19,25,26 The plasma potential within the spot settles at around 51 V, which is 1 V above the electrode potential. The scale length of the spot at this time (13 μs) is 0.5 cm, indicated by the 2D geometry of the zero electric field area shown in Fig. 9.
The duration of the simulation was limited to ∼13 μs due to its computational expense (100 000 CPU hours). For this reason, longer simulations which could assess the occasionally observed ∼100 μs period anode spot oscillations6,27,28 were not attempted.
C. Velocity distribution functions
The velocity distribution functions (VDFs) were examined at t = 10.8 μs by taking histograms of particle velocities in four 2D regions inside the anode spot. Ion data were obtained in a region 0.05 cm × 0.05 cm, and electron data in a 0.1 cm × 0.1 cm region (see Fig. 6). Different sized regions were used due to the different macroparticle weight of electrons and ions. The corresponding VDFs are shown in Fig. 7. Section II of Ref. 29 provides further details regarding the calculation of VDFs along the reflecting boundary.
Figure 6 shows that at y = 0.7 cm, a location within the double layer, electrons have a peak density near vy= –180 cm/μs corresponding to the ≈10 eV energy gained from the double layer electric field between 1 cm and 0.7 cm. These electrons increase in kinetic energy to nearly 26 eV at the maximum potential near y = 0.1–0.2 cm. This corresponds to the half-ring-like distribution with a peak density near v ≈ 300 cm/μs. The geometry of the VDF within the spot is due to the curvature of the surface of the double layer and the multitude of directions from which electrons enter the spot. Neutral scattering of double layer accelerated electrons is not expected to contribute to the isotropization of the VDF since the mean free path for elastic scattering is approximately 5 cm, an order of magnitude larger than the mean free path for ∼eV electrons and significantly larger than the anode spot size. In addition, electron thermalization due to streaming instabilities caused by the accelerated electrons is expected to be negligible. As Fig. 6 shows, the beam population density is very small compared with the trapped electron density. As a result, we do not expect significant electron-electron beam instabilities to form, and as expected, high frequency oscillations of the electric filed consistent with instability are not observed in the simulation.
Electrons born from ionization () within the spot have a nearly Maxwellian VDF. Few of these electrons are present between y = 0.5 cm and y = 0.7 cm since they are trapped by the strong double layer electric field. The VDFs also show that this population is much more dense than that of electrons from the bulk. This can also be seen in the density profiles of Fig. 6.
Ions born within the spot are well described by a Maxwellian distribution with no flow shift at y = 0.2 cm. However, very close to the electrode (y = 0.1 cm), ions are flow shifted due to the ion presheath leading up to the electrode surface. Within the double layer, ions are accelerated toward the bulk plasma as indicated by their velocity and decreased density.
D. Current collection
One common feature observed in experiments is a jump in the current collected by the electrode when a spot forms.3,19,25 It has previously been suggested that the cause of this increase is an increase of electron collecting area resulting from the relatively large surface area of the double layer compared with the electrode.6 Figure 8 shows that a current increase was also observed in the simulations. Most of the current collected at the electrode is due to the bulk plasma electrons () accelerated to the electrode by the double layer. There is also a ≈ 20% contribution from electrons formed by ionization after spot onset. This confirms that the majority of the increased current collection after onset is due to the increased collection area for electrons from the bulk plasma. In Fig. 8, the current component due to ions born from ionization collected at the walls of the domain is also shown. This indicates that approximately 5%–10% of the ion current collected at the walls after spot onset is due to ions born in the spot and accelerated to the bulk plasma by the double layer. Once the anode spot has formed, all of the bulk ion current is collected by the walls and all the bulk electron current is lost to the electrode. This situation, known as global non-ambipolar flow,26 occurs near electrodes that are small enough to be biased above the plasma potential and large enough to collect the entirety of the electron current. In this case, it is the effective collection area of the anode spot that establishes this flow scenario. After onset, the current shows oscillations on a microsecond timescale, possibly related to the increased fluctuations in the ion density and electric field shown in Fig. 4.
Additional evidence for the increase in electron current collection area can be found by integrating test electron trajectories starting at locations outside the double layer. The integration of the trajectories used the velocity Verlet method along with the electric field from a single time step of the simulation at t = 10.8 μs. The electric field within the anode spot had small fluctuations on the timescale of an electron crossing the spot, but these contributed to a negligible change in the trajectories, so the time dependence of the electric field was neglected. Selected particle trajectories are plotted over the electric field in Fig. 9. The initial velocities for the particles are typical of those found in a flow-shifted distribution at the sheath edge. Trajectories starting between x = 0 cm and x ≈ 0.4 cm impact the electrode, which extends from 0 to 0.25 cm on the lower boundary. Particles starting further to the right typically impact other parts of the lower boundary. These trajectories show that the anode spot increases the current collection area.
A. Anode spot onset
Section II B showed simulation results demonstrating that anode spot onset is a consequence of the build-up of positive space-charge within the sheath due to electron impact ionization of neutral atoms. Initially, electrons born from ionization were quickly accelerated to the electrode by the sheath electric field, but when the positive space-charge grew large enough, the potential profile developed a local maximum off the electrode surface. This local maximum is a potential well for electrons. Low energy electrons born from ionization within this well are trapped, allowing the electron density to increase. It was observed that this trapping led to the formation of a quasineutral region where E ≈ 0 and ne ≈ ni.
Once quasineutrality is established, the double layer is predicted to move if there is an imbalance in the Langmuir condition, resulting in the double layer velocity given in Eq. (2). To determine when the double layer will move outward, resulting in the expansion of the anode spot plasma, the flux densities entering the double layer are considered. Bulk electrons are accelerated to their thermal speed () in an electron presheath,24,30 leading up to the double layer at the low potential side. This flux density is
where the subscript B indicates values in the bulk plasma. Ions enter the high potential side of the double layer near the electrode at their sound speed due to the Bohm criterion, resulting in the flux density
where the subscript Hi indicates values at the high potential side of the double layer. Inserting Eqs. (6) and (7) into the modified Langmuir condition Eq. (2) and requiring UDL > 0 as a condition for anode spot onset31 leads to
Prior to the formation of the potential well, it may be thought that the modified Langmuir condition, Eq. (2), suggests UDL < 0 so the double layer will move toward the electrode. However, Eq. (2) is only valid when there is a region with E ≈ 0 at both the high and low potential sides of the double layer. Thus, quasineutrality must first be established on the high potential side, creating a region with E ≈ 0, before an imbalance in the static Langmuir condition implies expansion of the double layer.
Previous attempts to model spot onset solved an integral form of Poisson's equation.18,32,33 In these models, the integral of the charge density within the sheath, including contributions from electron impact ionization, was formulated as a function of the sheath electric field
Models based on Eq. (9) do not provide solutions which include a potential well for electrons because it requires that the electric field becomes multivalued as a function of the potential.
For typical experimental conditions where Te,B/Te,Hi ≃ 1, Eq. (8) predicts that the double layer will move when the high potential density exceeds the bulk density. Figure 4 shows that when this occurs, the double layer begins to move outward, which is in agreement with the condition in Eq. (8). After 10.7 μs, the double layer position indicated by the electric field in Fig. 4 shows that its speed has decreased. At first glance, this may seem to contradict Eq. (2) since the peak density within the center of spot has increased significantly even though the double layer velocity has not. However, a careful inspection of the density and potential profiles at t = 10.8 μs in Fig. 6 reveals that the plasma density at the high potential sheath edge is slightly less than half the peak density within the spot. The high potential density was still large enough for the double layer to continue to slowly move outward. The theory predicts that when the densities at the entrance to the double layer balance the Langmuir condition, its motion stops. The simulation time was not long enough for a steady-state configuration to be observed.
To predict when an anode spot will form, it is desirable to connect the density criterion of Eq. (8) with the electrode bias. In experiments, anode spot onset is observed when the electrode exceeds some critical bias. The critical bias can be predicted by considering the conditions required for spot onset. Two conditions are required: (1) a quasineutral region is established at the high potential side of the double layer so that the assumptions of the modified Langmuir condition are satisfied and the double layer is free to move, and (2) the ion flux leaving the high potential side leads to an imbalance in the Langmuir condition resulting in motion of the double layer and expansion of the anode spot plasma. Electrons born from ionization near the maximum of the potential profile are trapped if their energy is less than the depth of the well. A potential well for electrons is a potential hill for ions, therefore, for a deep enough well, the electron density increases faster than the ion density. This occurs if the depth of the well is greater than the magnitude of the floating potential. Once the electron density increases to match the ion density at the bottom of the well, condition (1) is satisfied. If ions are born faster than they are lost, the plasma density within the electron potential well will increase until condition (2) is satisfied. The statement that the ion birth rate is greater than the ion loss rate can be expressed as
where Asheath is the cross-sectional area of the sheath prior to the double layer expansion, cs,Hi is the ion sound speed at the high potential side, and zsheath is the sheath thickness. The numerical factor is due to half of the ions born near the electron well being lost in the electrode direction and the other half being lost in the plasma direction.
To determine the critical bias, the condition of exact balance in Eq. (10) is considered. Combining this with the steady-state Langmuir condition of Eq. (1) and the high potential ion flux in Eq. (7), the threshold condition becomes
This is similar to the equation used by Song et al.19 to estimate the sheath thickness of the anode glow prior to spot onset. While the use of the Langmuir condition for planar double layers is not strictly valid for the analysis of spherical double layers such as anode spots, the initial form of the double layer is approximately planar when attached to the electrode prior to spot formation. Using this condition along with an independent estimate of the electron sheath thickness24
Eq. (11) can be solved numerically for the critical bias Δϕc given the energy-dependent electron impact ionization cross section σI(eΔϕc). The critical bias for a helium plasma with an electron temperature of 2.4 eV is shown in Fig. 10 for several values of the bulk plasma density.
An interesting observation from experiments is that increasing the size of the electrode decreases the critical bias.3,7 In deriving Eq. (11), it was assumed that the sheath and electrode have the same surface area. This is appropriate for the simulations of Sec. II where the initial approximately planar electron sheath is embedded within the ion sheath of the surrounding walls. However, for an isolated electrode, the thickness of the sheath increases its surface area. If the sheath surface area is modeled as a disk with diameter D and finite thickness zsheath, the ratio of sheath area to electrode area is
Using this, Eq. (11) is modified to include the sheath area that was unaccounted for
Figure 11 shows a solution of Eq. (14) for an argon plasma along with data from the experiments of Ref. 3 conducted with 1 cm and 5.5 cm diameter electrodes. The plasma potential in the experiment was approximately 5 V, which was subtracted from the data points in Fig. 11, although an independent measurement was not available at each pressure. The predicted values of the critical bias are found to be in good agreement with the experimental data.
Section II B presents the results of a series of simulations used to determine the critical bias for a plasma with a 200 mTorr helium background, Te = 2.4 eV, and a bulk density of 3 × 108 cm−3. The critical bias for these conditions was determined to be 41.6 V above the plasma potential. This is approximately 10 V more than the critical bias estimated from Fig. 10. One possibility for the agreement with ∼1 mTorr experiments in Fig. 11 and discrepancy for ∼200 mTorr simulations is that elastic collisions of electrons result in a thicker sheath than that described by Eq. (12).34 A thicker sheath provides a greater loss area for ions formed by ionization; hence, a greater ionization rate may be needed to satisfy Eq. (14).
B. Anode spot current balance
As all steady-state plasmas, the anode spot maintains quasineutrality by the equal loss of electron and ion currents. Baalrud, Herskowitz, and Longmier developed a theory for the form of the sheath at a positive electrode using global current balance arguments.26 In their work, the current loss to a positively biased electrode of area AE in a plasma chamber with wall area AW was considered. The form of the sheath was predicted to be either an electron sheath, ion sheath, or a sheath with a non-monotonic potential, depending on the ratio of surface area of the positively biased electrode to the wall area AE/AW. When the electrode area was small, an electron sheath was present since it would not significantly modify the balance of global current loss. However, when the electrode area was large, the plasma potential locked to a few volts above the electrode potential, preventing electrons from being lost at a faster rate than ions. This behavior has been observed in both experiment and simulation.22,35
Similar arguments can be applied to the anode spot by considering the current lost through a double layer of surface area AS, which is the analog of AW for the quasineutral spot plasma, and the current lost to the electrode of area AE. For this situation, the form of the sheath between the electrode and spot plasma is expected to be either an ion sheath, electron sheath, or a non-monotonic electron sheath with virtual cathode depending on the area ratio AS/AE. These possibilities are shown in Fig. 12. Figure 4 shows that nearly all of the electrons in the spot are those formed by ionization of neutrals.36 Due to this observation, the analysis below will assume that the quasineutrality condition within the anode spot is between the ions and electrons produced by ionization interactions in Eqs. (4) and (5), i.e., neI ≈ ni.
1. Ion sheath
If the sheath between the spot plasma and electrode is an ion sheath, the electron current collected is , where ΔϕI = ϕS – ϕE, ϕS is the anode spot plasma potential, and is the average speed of electrons born from ionization assuming a Maxwellian VDF with temperature TeI. The ion current lost to the electrode and double layer is Ii = 0.6ecsni(AE + AS). Balancing the two currents gives the potential within the spot
Here, TeS = (TeIneI + TeBneB)/(neI + neB) is the total electron temperature within the spot. Both TeI and TeS appear in this equation since the total electron temperature TeS determines the Bohm speed for ions and TeI determines the loss rate for ionization electrons. Although neI ≫ neB, TeI ≪ TeB, so both components must be considered in the total electron temperature TeS. For Eq. (15) to be consistent with the assumption that an ion sheath is present (ΔϕI > 0), the area ratio must satisfy
2. Electron sheath
If the sheath between the spot plasma and electrode is an electron sheath, the electron current lost to the electrode is IeI = AEeneIvTeS, where . The ion current lost to the electrode and double layer is , where ΔϕE = ϕE – ϕS, and is the average speed for Maxwellian distributed ions with temperature Ti. The value of ΔϕE can be determined by balancing the two currents and using the quasineutrality condition resulting in
For this to be consistent with the assumptions of an electron sheath at the electrode (ΔϕE > 0)
3. Electron sheath with virtual cathode
For area ratios between the electron and ion sheath cases, the potential may have a non-monotonic profile such as a virtual cathode which limits the electron flux to the electrode.37
In experiments, larger anode spots at low pressure were observed to have electron sheath potential profiles both with and without virtual cathodes.3,19 An exact comparison between these experiments and the theory of this section is not possible since the electron temperature was not reported. For smaller anode spots at higher pressures, the potential profile is expected to include an ion sheath at the electrode. In the 2D simulation, the anode spot area is approximately half the circumference of a circle resulting in AS/AE ≈ 3. For this area ratio, the sheath at the electrode is expected to be an ion sheath. This is consistent with the observed potential profile shown in Fig. 5.
C. Particle and power balance and anode spot size
In this subsection, the steady-state size of the anode spot is predicted by considering conservation of particle number and power in a global model of the spot plasma. This global model ties the potential across the double layer to the ionization rate within the spot, a feature that was not present in previous estimations of the anode spot size. Once the size and double layer potential are known, the form of the sheath at the electrode can be predicted using the theory of Sec. III B.
Balancing the volume ionization rate with the ion loss rate, and neglecting the loss of particles to the electrode38
Here, nn is the neutral gas density, is the rate constant for impact ionization by electrons accelerated from the bulk, nB(ΔϕDL) is the density of the bulk electrons accelerated into the spot by ΔϕDL, and VS is the anode spot volume. A length scale can be defined by relating the spot surface area and volume as VS = ASL. Equation (19) can be written in terms of L
Using the approximation , and assuming that the cross section does not vary over the thermal width of the electron VDF
where the second equality results from the use of the Langmuir condition, returning the result of Eq. (3). The use of the Langmuir condition for planar double layers is applied to this situation assuming that the double layer thickness is small compared with the radius of curvature (or radius for a spherical anode spot) of the anode spot surface.16 The cross section in Eq. (21) is strongly dependent on the energy of electrons gained while passing the double layer. The effect of this energy dependence was not previously considered. Instead, estimates based on a constant cross section were used to determine the length scale.3,19 The length scale as a function of double layer energy is plotted from Eq. (21) in Fig. 13 for helium neutral background pressures of 50, 100, and 200 mTorr. To self consistently predict the size of the anode spot, a constraint on the double layer potential is needed.
In steady-state, the power entering and leaving the spot plasma will also balance. Most of the bulk electrons entering the spot pass through without undergoing an ionizing collision and exit at the electrode. However, a small fraction ionize neutrals, creating the spot plasma. For the non-ionizing electron component, the power entering and leaving the anode spot is in balance. The power input associated with the ionizing population is the number of ionization events per unit time multiplied by the residual energy left over after ionization of the neutral atom that goes into the kinetic energy of resulting particles. The power lost is the number of particles leaving the spot per unit time multiplied by the amount of energy that they carry with them. The resulting balance equation is
where ΓeI and Γi are the fluxes of electrons and ions leaving through their respective loss areas AeI and Ai at energy and , and is the threshold energy for electron impact ionization of neutrals. Making the approximation by assuming that the thermal width of electrons does not appreciably contribute to the ionization rate constant, Eq. (22) is
This can be evaluated for each of the three different anode spot potential profile configurations discussed in Sec. III B.
When the sheath at the anode is an electron sheath, the loss area for ions is the double layer surface area AS and the loss area for electrons is AE. Particle balance requires that the loss of electrons to the anode balance the loss of ions to the double layer resulting in
When the sheath at the anode is an ion sheath, ions are lost to the double layer and the electrode resulting in Ai = AS + AE. Once again, the ion flux is Γi = nScs and the energy is . The electron loss only occurs to the anode since ΔϕDL is much greater than the difference in the potential between the spot plasma and the anode. Particle balance results in . In this case, an electron lost to the anode is lost at the energy determined by the average speed of electrons, . Using these considerations, the double layer potential is
Note that the use of different temperatures TeI and TeS is due to the total electron temperature TeS controlling the Bohm speed, while TeI controls the thermal flux of the particles of the dominant electron population in the spot ( electrons).
Finally, consider the case where the sheath at the anode is an electron sheath with a virtual cathode; here, the collection areas are the same as in the previous case. The flux of ions to the double layer is given by the same relation used in the electron sheath case, while the flux of electrons is reduced by the dip potential of the virtual cathode and is . Using particle balance along with for ions and , the relation in Eq. (23) reduces to
The predictions of this section can be compared with the anode spot from the simulation in Sec. II B. Although these simulations never reached steady-state, the anode spot size varied slowly beginning around 10.7 μs after the initial double layer motion around 10.2 μs (see the electric field magnitude in Fig. 4). Using the available particle velocity data from Fig. 7, the temperatures can be calculated allowing a comparison. See Sec. II of Ref. 29 for details of the temperature calculation. From these data, the total electron temperature within the anode spot due to both and components is TeS = 3.8 eV, while the temperature for the component is TeI = 1.15 eV. At this time, the 2D electric field shows that the area ratio of the double layer surface to electrode surface is roughly AS/AE ≈ 3. Inserting these quantities into Eq. (26) results in V, which is close to the value V from the double layer potential in Fig. 6.
D. Ion presheath in the anode spot
At the high potential side of the double layer, ions leave the anode spot plasma entering into a positive space charge region. From the point of view of an ion within the spot plasma, this is similar to an ion entering an ion sheath. The Bohm criterion applies. An ion must enter the non-neutral region at a velocity exceeding its sound speed. This velocity is attained in a presheath region leading up to the sheath edge. It was previously suggested that the anode spot presheath length was determined by the ion-neutral collision mean free path,3 although in some cases this would result in a presheath significantly longer than the observed size of the anode spot.
The ion presheath length can be predicted by considering the rate constants for processes in the presheath. If the rate constant for ion-neutral collisions is greater than that for electron impact ionization, the presheath length would be determined by the ion-neutral mean free path. However, if the rate constant for electron impact ionization is larger, the presheath is dominated by ionization. In the latter case, the presheath length is predicted to be half of the plasma length, in this case, half of the anode spot size L. This is one of the main results described by the various formulations of the Tonks-Langmuir presheath models.39,40 In this section, the rate constants for these two processes are estimated, showing that the presheath is ionization dominated for typical plasma parameters.
The rate constant is defined as averaged over the distribution function of the incident and background particles denoted by subscripts 1 and 2, respectively. This is
where σ is the interaction cross section and the velocity distribution functions are normalized to unity. If the characteristic velocities of the background particles are much less than those of the incident particles, then and the v2 integral can be evaluated.
Equation (28) is evaluated assuming that the particles accelerated by the presheath have a flow shifted Maxwellian distribution, , where U is the flow shift and is the thermal speed. Writing the integration variable in terms of energy by using , and evaluating the integral in cylindrical coordinates aligned with the flow direction
The ratio of the rate constants for electron impact ionization and the sum of ion-neutral elastic and charge exchange collisions is considered, calculated using cross sections obtained from the Phelps database on LXcat,23 to determine which processes is dominant. Because of the strong dependence on energy, the rate constant for ion neutral collisions is calculated using Eq. (29) with a flow of cs/2 which is typical of a presheath ion. The ionization rate constant is estimated as
The ratio of these rate constants is shown in Fig. 14. The rate constant for electron impact ionization is several times that due to ion-neutral scattering for room temperature ions and values of eΔϕ ∼ 25 eV. Equations (25)–(27) predict that the double layer potential is approximately an electron temperature above the ionization energy of 24.6 eV. Since the temperatures encountered are typically , this predicts that the presheath length will be determined by ionization processes for a helium neutral gas background. Similar results are obtained for argon plasmas.
In this paper, the anode spot was studied for the first time using PIC simulations. These simulations demonstrated that electron impact ionization within the sheath results in a positive space charge layer adjacent to the electrode. With sufficient ionization, the positive space charge forms a potential well which traps low energy electrons formed by ionization in front of the electrode. This electron trapping results in an increase in electron density and formation of a quasineutral plasma.
A model for the spot onset was formulated based on observations of the simulated anode spot. The main feature of this model is that an imbalance in flux densities crossing the double layer leads to its motion and the expansion of the anode spot plasma. Using estimates of the sheath ionization rate, the value of the electrode bias relative to the plasma potential was tied to the sheath ionization. Predictions of the critical bias for spot onset were found to be in agreement with past experiments for different electrode sizes and plasma conditions, allowing an experimental test of the spot onset model.
Steady-state properties were predicted based on an analysis of current, power, and particle balance of the spot plasma. Maintenance of quasineutrality within the spot dictates the form of the sheath between the anode spot plasma and the electrode, determining how particles are lost from the anode spot. In the model, balance of the total ionization rate and particle loss rate determines the anode spot size as a function of the energy of electrons entering the spot from the bulk. The size is determined once the double layer potential is known. Balance of power lost from and deposited into the spot plasma sets this potential in the model. The predicted energy gain of an electron crossing the double layer potential is predicted to be a few electron volts above the ionization energy for typical experimental conditions. This result is consistent with several anode spot experiments.
The authors thank Nathaniel Shaffer for his comments on the manuscript and Robert Merlino and Sina Jahanbakhsh for useful discussion. This research was supported by the Office of Fusion Energy Science at the U.S. Department of Energy under Contract No. DE-AC04-94SL85000. The first author was also supported by the U.S. Department of Energy, Office of Science, Office of Workforce Development for Teachers and Scientists, Office of Science Graduate Student Research (SCGSR) program. The SCGSR program is administered by the Oak Ridge Institute for Science and Education for the DOE under Contract No. DE-AC05-06OR23100.