Understanding the role of physical processes contributing to breakdown is critical for many applications in which breakdown is undesirable, such as capacitors, and applications in which controlled breakdown is intended, such as plasma medicine, lightning protection, and materials processing. The electron emission from the cathode is a critical source of electrons which then undergo impact ionization to produce electrical breakdown. In this study, the role of secondary electron yields due to photons (γ ph) and ions (γ i) in direct current breakdown is investigated using a particle-in-cell direct simulation Monte Carlo model. The plasma studied is a one-dimensional discharge in 50 Torr of pure helium with a platinum cathode, gap size of 1.15 cm, and voltages of 1.2-1.8 kV. The current traces are compared with experimental measurements. Larger values of γ ph generally result in a faster breakdown, while larger values of γ i result in a larger maximum current. The 58.4 nm photons emitted from He(21P) are the primary source of electrons at the cathode before the cathode fall is developed. Of the values of γ ph and γ i investigated, those which provide the best agreement with the experimental current measurements are γ ph = 0.005 and γ i = 0.01. These values are significantly lower than those in the literature for pristine platinum or for a graphitic carbon film which we speculate may cover the platinum. This difference is in part due to the limitations of a one-dimensional model but may also indicate surface conditions and exposure to a plasma can have a significant effect on the secondary electron yields. The effects of applied voltage and the current produced by a UV diode which was used to initiate the discharge, are also discussed.
This paper describes the verification and validation (V&V) framework developed for the stochastic Particle-in-Cell, Direct Simulation Monte Carlo code Aleph. An ideal framework for V&V from the viewpoint of the authors is described where a physics problem is defined, and relevant physics models and parameters to the defined problem are assessed and captured in a Phenomena Identification and Ranking Table (PIRT). Furthermore, numerous V&V examples guided by the PIRT for a simple gas discharge are shown to demonstrate the V&V process applied to a real-world simulation tool with the overall goal to demonstrably increase the confidence in the results for the simulation tool and its predictive capability. Although many examples are provided here to demonstrate elements of the framework, the primary goal of this work is to introduce this framework and not to provide a fully complete implementation, which would be a much longer document. Comparisons and contrasts are made to more usual approaches to V&V, and techniques new to the low-temperature plasma community are introduced. Specific challenges relating to the sufficiency of available data (e.g., cross sections), the limits of ad hoc validation approaches, the additional difficulty of utilizing a stochastic simulation tool, and the extreme cost of formal validation are discussed.
Carbone, Emile; Graef, Wouter; Hagelaar, Gerjan; Boer, Daan; Hopkins, Matthew M.; Stephens, Jacob C.; Yee, Benjamin T.; Pancheshnyi, Sergey; Van Dijk, Jan; Pitchford, Leanne
Technologies based on non-equilibrium, low-temperature plasmas are ubiquitous in today’s society. Plasma modeling plays an essential role in their understanding, development and optimization. An accurate description of electron and ion collisions with neutrals and their transport is required to correctly describe plasma properties as a function of external parameters. LXCat is an open-access, web-based platform for storing, exchangig and manipulating data needed for modeling the electron and ion components of non-equilibrium, low-temperature plasmas. The data types supported by LXCat are electron- and ion-scattering cross-sections with neutrals (total and differential), interaction potentials, oscillator strengths, and electron- and ion-swarm/transport parameters. Online tools allow users to identify and compare the data through plotting routines, and use the data to generate swarm parameters and reaction rates with the integrated electron Boltzmann solver. In this review, the historical evolution of the project and some perspectives on its future are discussed together with a tutorial review for using data from LXCat.
Helium is frequently used as a working medium for the generation of plasmas and is capable of energetic photon emissions. These energetic photon emissions are often attributed to the formation of helium excimer and subsequent photon emission. When the plasma device is exposed to another gas, such as nitrogen, this energetic photon emission can cause photoionization and further ionization wave penetration into the additional gas. Often ignored are the helium resonance emissions that are assumed to be radiation trapped and therefore not pertinent to photoionization. Here, experimental evidence for the presence of helium atomic emission in a pulsed discharge at ten's of Torr is shown. Simulations of a discharge in similar conditions agree with the experimental measurements. In this context, the role of atomic and molecular helium light emission on photoionization of molecular nitrogen in an ionization wave is studied using a kinetic modeling approach that accounts for radiation dynamics in a developing low-temperature plasma. Three different mixtures of helium at a total pressure of 250 Torr are studied in simulation. Photoionization of the nitrogen molecule by vacuum ultraviolet helium emission is used as the only seed source ahead of the ionization front. It is found that even though radiation trapped, the atomic helium emission lines are the significant source of photoionization of nitrogen. The significant effect of radiation trapped photon emission on ionization wave dynamics demonstrates the need to consider these radiation dynamics in plasma reactors where self-absorbed radiation is ignored.
The boundary regions of low-temperature plasmas are known to be susceptible to kinetic instabilities, which can affect the energies and fluxes of particles directed at the material boundary. For example, both the ion acoustic instability as well as an instability near the electron plasma frequency have been observed. Particle-in-cell (PIC) simulation is a tool that, alongside experiments, can capture the effects these instabilities have on the particle distribution functions. Ultimately, simulations can determine under what conditions these effects are significant by comparing to theoretical predictions and explore conditions unamenable to experiments.
Clark, Raimi; Brooks, William; Hopkins, Matthew M.; Mankowski, John; Stephens, Jacob; Neuber, Andreas
Early light emission provides information about the dominant mechanisms culminating in vacuum surface flashover (anode-initiated vs. cathode-initiated) for particular geometries. From experimental evidence gathered elsewhere, for the case of an insulator oriented at 45° with respect to the anode, anode-initiated flashover is believed to dominate since the field at the anode triple point is roughly three times that of the cathode. Similar to previous work performed on cathode-initiated flashover, light emission from the voltage rise through the impedance collapse is collected into two optical fibers focused on light emanating from the insulator in regions near the anode and cathode. The optical fibers are either connected to PMTs for spectrally integrated localized light intensity information or to a spectrograph used in conjunction with an ICCD camera. Challenges associated with localizing the flashover for optical diagnostics and incorporating the optical diagnostics into the high-field environment are discussed. Initial results for cross-linked polystyrene (Rexolite 1422) support the premise that flashover is initiated from the anode for these geometries, as early light from the anode leads cathode light up to photocathode saturation. Early spectroscopy results show promise for future characterization of the spatio-temporal development of emission from desorbed gas species across the insulator surface and identification of bulk insulator involvement if it occurs.
3D Particle-In-Cell Direct Simulation Monte Carlo (PIC-DSMC) simulations of cm-sized devices cannot resolve atomic-scale (nm) surface features and thus one must generate micron-scale models for an effective “local” work function, field enhancement factor, and emission area. Here we report on development of a stochastic effective model based on atomic-scale characterization of as-built electrode surfaces. Representative probability density distributions of the work function and geometric field enhancement factor (beta) for a sputter-deposited Pt surface are generated from atomic-scale surface characterization using Scanning Tunneling Microscopy (STM), Atomic Force Microscopy (AFM), and Photoemission Electron Microscopy (PEEM). In the micron-scale model every simulated PIC-DSMC surface element draws work functions and betas for many independent “atomic emitters”. During the simulation the field emitted current from an element is computed by summing each “atomic emitter's” current. This model has reasonable agreement with measured micron-scale emitted currents across a range of electric field values.
PIC MCC simulation results on the breakdown in the pulse discharge in helium at pressure of 100 Torr and voltage of U=3.25 kV are presented. The delay of the breakdown development is studied with different initial densities of plasma and excited helium atoms, which corresponds to various discharge operation frequencies. It is shown that for high concentration of excited atoms the photoemission determines the breakdown delay time. In opposite case of low excited atoms density, the ion-electron emission plays a key role in the breakdown development. The photoemission from the cathode is set with a flux of the photons with Doppler shift over the frequency. These photons are generated in reactions between exited atoms and fast atoms. A wide distribution of breakdown delay time was observed in different runs and analyzed.
The purpose of this paper is to characterize the need for improved predictive capabilities in low-temperature plasma (LTP) science, and to identify possible means of accomplishing this. While these means may constitute an initiative of their own, we consider these ideas to have widespread importance to discovery plasma science. Therefore, it is our hope that these ideas are more generally incorporated in future work.
Modern computational validation efforts rely on comparison of known experimental quantities such as current, voltage, particle densities, and other plasma properties with the same values determined through simulation. A discrete photon approach for radiation transport was recently incorporated into a particle-in-cell/direct simulation Monte Carlo code. As a result, spatially and temporally resolved synthetic spectra may be generated even for non-equilibrium plasmas. The generation of this synthetic spectra lends itself to potentially new validation opportunities. In this work, initial comparisons of synthetic spectra are made with experimentally gathered optical emission spectroscopy. A custom test apparatus was constructed that contains a 0.5 cm gap distance parallel plane discharge in ultra high purity helium gas (99.9999%) at a pressure of 75 Torr. Plasma generation is initiated with the application of a fast rise-time, 100 ns full-width half maximum, 2.0 kV voltage pulse. Transient electrical diagnostics are captured along with time-resolved emission spectra. A one-dimensional simulation is run under the same conditions and compared against the experiment to determine if sufficient physics are included to model the discharge. To sync the current measurements from experiment and simulation, significant effort was undertaken to understand the kinetic scheme required to reproduce the observed features. Additionally, the role of the helium molecule excimer emission and atomic helium resonance emission on photocurrent from the cathode are studied to understand which effect dominates photo-feedback processes. Results indicate that during discharge development, atomic helium resonance emission dominates the photo-flux at the cathode even though it is strongly self-absorbed. A comparison between the experiment and simulation demonstrates that the simulation reproduces observed features in the experimental discharge current waveform. Furthermore, the synthesized spectra from the kinetic method produces more favorable agreement with the experimental data than a simple local thermodynamic equilibrium calculation and is a first step towards using spectra generated from a kinetic method in validation procedures. The results of this study produced a detailed compilation of important helium plasma chemistry reactions for simulating transient helium plasma discharges.
A fully resolved kinetic model (particle-in-cell and direct simulation Monte Carlo for particle/photon collisions) of a near atmospheric pressure ionization wave is presented here. Fully resolving the required numerical spatial (sub-μm) and temporal scales (tens of fs) for atmospheric pressure discharges in three-dimensions is still a challenging task on modern super computers. To keep the overall problem tractable, the total number of elements are reduced by only simulating a 10° wedge rather than a full 360° geometry. The ionization wave is generated in a needle-plane configuration with a gap size of 250 μm and a background of nitrogen and helium gas. A voltage of 1500 V is applied to the anode and an initial electron and ion density of 109 cm-3 is seeded in a region near the anode electrode tip and extending towards the cathode. As these initial electrons are swept away, photoionization and photoemission create new electrons and allow the ionization front to propagate towards the cathode. Results from the 90% N2, 10% He discharge indicate that photoionization has minimal impact on plasma formation processes and cathode photoemission is the dominant mechanism for new electrons. In the 90% He, 10% N2 discharge case, however, photoionization likely has an impact as the observed locations of photoionization occur far enough away from the ionization front to allow for sufficient avalanche processes that contribute to the propagation of the ionization wave. Additionally, the electron energy distribution functions in the 90% He, 10% N2 case indicate that there is less energy loss to the low lying molecular N2 electronic states as well as the vibrational and rotational modes. This leads to higher electron energies and faster plasma development times of ∼0.4 ns for the 90% He, 10% N2 case, and ∼1.5 ns for the 90% N2, 10% He case. In addition to analysis of the ionization wave results, the overall challenges associated with simulations near atmospheric pressure discharges in three-dimensions are discussed, including the limitations of the 10° wedge that produces, at least qualitatively, minimal 3D effects.
