Abstract: We studied the evolution of an isolated pulsed plasma in a vortex flow stabilised microwave (MW) discharge in N2 at 25 mbar via the combination of 0D kinetics modelling, iCCD imaging and laser scattering diagnostics. Quenching of electronically excited N2 results in fast gas heating and the onset of a thermal-ionisation instability, contracting the discharge volume. The onset of a thermal-ionisation instability driven by vibrational excitation pathways is found to facilitate significantly higher N2 conversion (i.e. dissociation to atomic N2 ) compared to pre-instability conditions, emphasizing the potential utility of this dynamic in future fixation applications. The instability onset is found to be instigated by super-elastic heating of the electron energy distribution tail via vibrationally excited N2 . Radial contraction of the discharge to the skin depth is found to occur post instability, while the axial elongation is found to be temporarily contracted during the thermal instability onset. An increase in power reflection during the thermal instability onset eventually limits the destabilising effects of exothermic electronically excited N2 quenching. Translational and vibrational temperature reach a quasi-non-equilibrium after the discharge contraction, with translational temperatures reaching ∼1200 K at the pulse end, while vibrational temperatures are found in near equilibrium with the electron energy (1 eV, or ∼11 600 K). This first description of the importance of electronically excited N2 quenching in thermal instabilities gives an additional fundamental understanding of N2 plasma behaviour in pulsed MW context, and thereby brings the eventual implementation of this novel N2 fixation method one step closer.

Abstract: The limitation in time and length scale is a major issue of molecular dynamics (MD) simulation. Although several methods have been developed to extend the MD time scale, their performance usually deteriorates with increasing system size. Therefore, an acceleration method which is applicable to large systems is required to bridge the gap between the MD simulations and target phenomena. In this study, an accelerated MD method for large system is developed based on the collective variable-driven hyperdynamics (CVHD) method [K.M. Bal and E.C. Neyts, 2015]. The key idea is to run CVHD in parallel with rate control and accelerate multiple possible events simultaneously. Using this novel method, carbon diffusion in bcc-iron bicrystal with grain boundary is examined as an application for practical materials. Carbon atoms reaching at the grain boundary are trapped whereas carbon atoms in the bulk region diffuse randomly, and both dynamic regimes can be simultaneously accelerated with the parallel CVHD technique.

Abstract: A hybrid model is developed for describing the effects of oxygen addition to argon glow discharges. The species taken into account in the model include Ar atoms in the ground state and the metastable level, O2 gas molecules in the ground state and two metastable levels, O atoms in the ground state and one metastable level, O3 molecules, Ar+, O+, O2+ and O− ions, as well as the electrons. The hybrid model consists of a Monte Carlo model for electrons and fluid models for the other plasma species. In total, 87 different reactions between the various plasma species are taken into account. Calculation results include the species densities and the importance of their production and loss processes, as well as the dissociation degree of oxygen. The effect of different O2 additions on these calculation results, as well as on the sputtering rates, is discussed.

Abstract: A computer model is developed for describing argon/nitrogen glow discharges. The species taken into account in the model include electrons, Ar atoms in the ground state and in the 4s metastable levels, N2 molecules in the ground state and in six different electronically excited levels, N atoms, Ar+ ions, N+, N2+, N3+ and N4+ ions. The fast electrons are simulated with a Monte Carlo model, whereas all other species are treated in a fluid model. 74 different chemical reactions are considered in the model. The calculation results include the densities of all the different plasma species, as well as information on their production and loss processes. The effect of different N2 additions, in the range between 0.1 and 10%, is investigated.