Multi-mode ionization instability induced striations in RF driven He/H 2 O atmospheric pressure plasma (APP) discharges
Previously [E. Kawamura et al., Plasma Sources Sci. Technol. 25, 054009 (2016) and E. Kawamura et al., J. Phys. D: Appl. Phys. 50, 145204 (2017)], one dimensional (1D) particle-in-cell (PIC) simulations of 1 to 4mm gap, He/2%H 2 O atmospheric pressure plasma discharges showed an ionization instability resulting in bulk striations. Assuming that the ionization rate coefficient K iz is related to the root mean square electric field E by K iz ∝ E q , a striation theory showed that q<0 is a necessary condition for the instability. A local calculation yielded q>0, implying that nonlocal electron kinetics are required for the instability. Wider gaps can fit in a wider range of wavelengths, resulting in multi-mode striations. Previously, we assumed a single mean q value for each discharge, and did not calculate q for each mode separately. Here, we develop a wavelength resolved global striation model and apply it to 1D PIC simulations of 4mm gap discharges with input currents J rf = 0.04 to 0.30A/cm 2 . We first examine a base case at 0.23A/cm 2 and observe a mixture of unstable modes within a window of wavelengths λ. At shorter λ, the modes are suppressed by diffusion. At longer λ, we observe a transition to locality in which q becomes less negative with increasing λ, approaching its local positive value and stabilizing the modes. The window of unstable modes shifts to shorter λ with increasing J rf , causing the modes to be suppressed by diffusion at higher J rf . At lower J rf , the decrease in bulk plasma density with decreasing J rf suppresses the striations.