Abstract: Secondary electron emission is an important process that plays a significant role in several plasma‐related applications. As measuring the secondary electron yield experimentally is very challenging, quantitative modelling of this process to obtain reliable yield data is critical as input for higher‐scale simulations. Here, we build upon our previous work combining density functional theory calculations with a model originally developed by Hagstrum to extend its application to metallic surfaces. As plasmonic effects play a much more important role in the secondary electron emission mechanism for metals, we introduce an approach based on Poisson point processes to include both surface and bulk plasmon excitations to the process. The resulting model is able to reproduce the yield spectra of several available experimental results quite well but requires the introduction of global fitting parameters, which describe the strength of the plasmon interactions. Finally, we use an in‐house developed workflow to calculate the electron yield for a list of elemental surfaces spanning the periodic table to produce an extensive data set for the community and compare our results with more simplified approaches from the literature.

Abstract: technologies have been expanding, and one of the most exciting and rapidly growing
applications is in biology and medicine. Most biomedical studies with DBD plasma systems are performed in vitro, which include cells grown on the surface of plastic well plates, or in vivo, which include animal research models (e.g. mice, pigs). Since many DBD systems use the biological target as the secondary electrode for direct plasma generation and treatment, they are sensitive to the surface properties of the target, and thus can be altered based on the in vitro or in vivo system used. This could consequently affect biological response from plasma treatment. Therefore, in this study, we investigated the DBD plasma behavior both in vitro (i.e. 96-well flat bottom plates, 96-well U-bottom plates, and 24-well flat bottom plates), and in vivo (i.e. mouse skin). Intensified charge coupled device (ICCD) imaging was performed and the plasma discharges were visually distinguishable between the different systems. The geometry of the wells did not affect DBD plasma generation for low application distances (≤ 2 mm), but differentially affected plasma uniformity on the bottom of the well at greater distances. Since DBD plasma treatment in vitro is rarely performed in dry wells for plasma medicine experiments, the effect of well wetness was also investigated. In all in vitro cases, the uniformity of the DBD plasma was affected when comparing wet versus dry wells, with the plasma in the wide-bottom wells appearing the most similar to plasma generated on mouse skin. Interestingly, based on quantification of ICCD images, the DBD plasma intensity per surface area demonstrated an exponential one-phase decay with increasing application distance, regardless of the in vitro or in vivo system. This trend is similar to that of the energy per pulse of plasma, which is used to determine the total plasma treatment energy for biological systems. Optical emission spectroscopy performed on the plasma revealed similar trends in radical species generation between the plastic well plates and mouse skin. Therefore, taken together, DBD plasma intensity per surface area may be a valuable parameter to be used as a simple method for in situ monitoring during biological treatment and active plasma treatment control, which can be applied for in vitro and in vivo systems.

Abstract: The use of microfluidic devices in the field of plasma-liquid interaction can unlock unique possibilities to investigate the effects of plasma-generated reactive species for environmental and biomedical applications. So far, very little simulation work has been performed on microfluidic devices in contact with a plasma source. We report on the modelling and computational simulation of physical and chemical processes taking place in a novel plasma-microfluidic platform. The main production and transport pathways of reactive species both in plasma and liquid are modelled by a novel modelling approach that combines 0D chemical kinetics and 2D transport mechanisms. This combined approach, applicable to systems where the transport of chemical species occurs in unidirectional flows at high Péclet numbers, decreases calculation times considerably compared to regular 2D simulations. It takes advantage of the low computational time of the 0D reaction models while providing spatial information through multiple plug-flow simulations to yield a quasi-2D model. The gas and liquid flow profiles are simulated entirely in 2D, together with the chemical reactions and transport of key chemical species. The model correctly predicts increased transport of hydrogen peroxide into the liquid when the microfluidic opening is placed inside the plasma effluent region, as opposed to inside the plasma region itself. Furthermore, the modelled hydrogen peroxide production and transport in the microfluidic liquid differs by less than 50% compared with experimental results. To explain this discrepancy, the limits of the 0D–2D combined approach are discussed.

Abstract: Non-equilibrium inductively coupled plasmas (ICPs) operating in hydrogen are of significant interest for applications including large-area materials processing. Increasing control of spatial gas heating, which drives the formation of neutral species density gradients and the rate of gas-temperature-dependent reactions, is critical. In this study, we use 2D fluid-kinetic simulations with the Hybrid Plasma Equipment Model to investigate the spatially resolved production of atomic hydrogen in a low-pressure planar ICP operating in pure hydrogen (10-20 Pa or 0.075-0.15 Torr, 300 W). The reaction set incorporates self-consistent calculation of the spatially resolved gas temperature and 14 vibrationally excited states. We find that the formation of neutral-gas density gradients, which result from spatially non-uniform electrical power deposition at constant pressure, can drive significant variations in the vibrational distribution function and density of atomic hydrogen when gas heating is spatially resolved. This highlights the significance of spatial gas heating on the production of reactive species in relatively high-power-density plasma processing sources.

Abstract: The development of compact, low power, charge-neutral propulsion sources is of significant recent interest due to the rising application of micro-scale satellite platforms. Among such sources, radio frequency (rf) electrothermal microthrusters present an attractive option due to their scalability, reliability, and tunable control of power coupling to the propellant. For micropropulsion applications, where available power is limited, it is of particular importance to understand how electrical power can be transferred to the propellant efficiently, a process that is underpinned by the plasma sheath dynamics. In this work, two-dimensional fluid/Monte Carlo simulations are employed to investigate the effects of applied voltage frequency on the electron, ion, and neutral heating in an rf capacitively coupled plasma microthruster operating in argon. Variations in the electron and argon ion densities and power deposition, and their consequent effect on neutral-gas heating, are investigated with relation to the phase-averaged and phase-resolved sheath dynamics for rf voltage frequencies of 6-108 MHz at 450 V. Driving voltage frequencies above 40.68 MHz exhibit enhanced volumetric ionization from bulk electrons at the expense of the ion heating efficiency. Lower driving voltage frequencies below 13.56 MHz exhibit more efficient ionization due to secondary electrons and an increasing fraction of rf power deposition into ions. Thermal efficiencies are improved by a factor of 2.5 at 6 MHz as compared to the more traditional 13.56 MHz, indicating a favorable operating regime for low power satellite applications.