Abstract: In this paper, a simulation method is described to predict the etching behaviour of Cl2/O2/Ar inductively coupled plasmas on a Si substrate, as used in shallow trench isolation for the production of electronic devices. The hybrid plasma equipment model (HPEM) developed by Kushner et al is applied to calculate the plasma characteristics in the reactor chamber and two additional Monte Carlo simulations are performed to predict the fluxes, angles and energy of the plasma species bombarding the Si substrate, as well as the resulting surface processes such as etching and deposition. The simulations are performed for a wide variety of operating conditions such as gas composition, chamber pressure, power deposition and substrate bias. It is predicted by the simulations that when the fraction of oxygen in the gas mixture is too high, the oxidation of the Si substrate is superior to the etching of Si by chlorine species, resulting in an etch rate close to zero as is also observed in the experiments.

Abstract: In this work, we extensively investigate the structural and electronic properties of van der Waals heterostructures (HTs) constructed by MoS${2}$/$BC₃$, MoS${2}$/$C₃N$, MoS${2}$/$C₃N₄$, MoS${2}$/$C₄N₃$ and those using Janus MoSSe instead of MoS$2$ by performing density functional theory calculations. The electronic band structure calculations and the corresponding partial density of states reveal that the significant changes are driven by quite strong layer-layer interaction between the constitutive layers. Our results show that although all monolayers are semiconductors as free-standing layers, the MoS${2}$/$C₃N$ and MoS${2}$/$C₄N₃$ bilayer HTs display metallic behavior as a consequence of transfer of charge carriers between two constituent layers. In addition, it is found that in MoSSe/$C₃N$ bilayer HT, the degree of metallicity is affected by the interface chalcogen atom type when Se atoms are facing to $C₃N$ layer, the overlap of the bands around the Fermi level is smaller. Moreover, the half-metallic magnetic $C₄N₃$ is shown to form magnetic half-metallic trilayer HT with MoS$2$ independent of the stacking sequence, i.e. whether it is sandwiched or two $C₄N₃$ layer encapsulate MoS$2$ layer. We further analyze the trilayer HTs in which MoS$2$ is encapsulated by two different monolayers and it is revealed that at least with one magnetic monolayer, it is possible to construct a magnetic trilayer. While the trilayer of $C₄N₃$/MoS${2}$/$BC₃$ and $C₄N₃$/MoS${2}$/$C₃N₄$ exhibit half-metallic characteristics, $C₄N₃$/MoS${_2}$/$C₃$N possesses a magnetic metallic ground state. Overall, our results reveal that holly structures of BCN crystals are suitable for heterostructure formation even over van der Waals type interaction which significantly changes electronic nature of the constituent layers.

Abstract: A two-dimensional fluid model, including the full set of Maxwell equations, has been developed and applied to investigate the effect of a phase shift between two power sources on the radial uniformity of several plasma characteristics in a hydrogen capacitively coupled plasma. This study was carried out at various frequencies in the range 13.56200 MHz. When the frequency is low, at 13.56 MHz, the plasma density is characterized by an off-axis peak when both power sources are in-phase (phgr = 0), and the best radial uniformity is obtained at phgr = π. This trend can be explained because the radial nonuniformity caused by the electrostatic edge effect can be effectively suppressed by the phase-shift effect at a phase difference equal to π. When the frequency rises to 60 MHz, the plasma density profiles shift smoothly from edge-peaked over uniform to centre-peaked as the phase difference increases, due to the pronounced standing-wave effect, and the best radial uniformity is reached at phgr = 0.3π. At a frequency of 100 MHz, a similar behaviour is observed, except that the maximum of the plasma density moves again towards the radial edge at the reverse-phase case (phgr = π), because of the dominant skin effect. When the frequency is 200 MHz, the bulk plasma density increases significantly with increasing phase-shift values, and a better uniformity is obtained at phgr = 0.4π. This is because the density in the centre increases faster than at the radial edge as the phase difference rises, due to the increasing power deposition Pz in the centre and the decreasing power density Pr at the radial edge. As the phase difference increases to π, the maximum near the radial edge becomes obvious again. This is because the skin effect has a predominant influence on the plasma density under this condition, resulting in a high density at the radial edge. Moreover, the axial ion flux increases monotonically with phase difference, and exhibits similar profiles to the plasma density. The calculation results illustrate that the radial uniformity of the various plasma characteristics is strongly dependent on the applied frequency and the phase shift between both power sources, which is important to realize, for controlling the uniformity of the plasma etch and deposition processes.

Abstract: The atmospheric pressure dielectric barrier discharge in helium is a pulsed discharge in nature. If during the electrical current pulse a glow discharge is reached, then this pulse will last only a few microseconds in operating periods of sinusoidal voltage with lengths of about 10 to 100 µs. In this paper we demonstrate that right before a glow discharge is reached, the discharge very closely resembles the commonly assumed Townsend discharge structure, but actually contains some significant differing features and hence should not be considered as a Townsend discharge. In order to clarify this, we present calculation results of high time and space resolution of the pulse formation. The results indicate that indeed a maximum of ionization is formed at the anode, but that the level of ionization remains high and that the electric field at that time is significantly disturbed. Our results also show where this intermediate structure comes from.

Abstract: In the last twenty years, numerical modeling has become an indispensable part of magnetism research. It has become a standard tool for both the exploration of new systems and for the interpretation of experimental data. In the last five years, the capabilities of micromagnetic modeling have dramatically increased due to the deployment of graphical processing units (GPU), which have sped up calculations to a factor of 200. This has enabled many studies which were previously unfeasible. In this topical review, we give an overview of this modeling approach and show how it has contributed to the forefront of current magnetism research.

Abstract: Microdischarges are an important type of plasma discharge that possess several unique characteristics, such as the presence of a stable glow discharge, high plasma density and intense excimer radiation, leading to several potential applications. The intense and controllable gas heating within the extremely small dimensions of microdischarges has been exploited in microthruster technologies by incorporating a micro-nozzle to generate the thrust. This kind of microthruster has a significantly improved specific impulse performance compared to conventional cold gas thrusters, and can meet the requirements arising from the emerging development and application of micro-spacecraft. In this paper, we performed a self-consistent 2D particle-in-cell simulation, with a Monte Carlo collision model, of a microdischarge operating in a prototype micro-plasma thruster with a hollow cylinder geometry and a divergent micro-nozzle. The model takes into account the thermionic electron emission including the Schottky effect, the secondary electron emission due to cathode bombardment by the plasma ions, several different collision processes, and a non-uniform argon background gas density in the cathode-anode gap. Results in the high-pressure (several hundreds of Torr), high-current (mA) operating regime showing the behavior of the plasma density, potential distribution, and energy flux towards the hollow cathode and anode are presented and discussed. In addition, the results of simulations showing the effect of different argon gas pressures, cathode material work function and discharge voltage on the operation of the microdischarge thruster are presented. Our calculated properties are compared with experimental data under similar conditions and qualitative and quantitative agreements are reached.

Abstract: In this review paper, an overview is given of different modelling efforts for plasmas used for the formation and growth of nanostructured materials. This includes both the plasma chemistry, providing information on the precursors for nanostructure formation, as well as the growth processes itself. We limit ourselves to carbon (and silicon) nanostructures. Examples of the plasma modelling comprise nanoparticle formation in silane and hydrocarbon plasmas, as well as the plasma chemistry giving rise to carbon nanostructure formation, such as (ultra)nanocrystalline diamond ((U)NCD) and carbon nanotubes (CNTs). The second part of the paper deals with the simulation of the (plasma-based) growth mechanisms of the same carbon nanostructures, i.e. (U)NCD and CNTs, both by mechanistic modelling and detailed atomistic simulations.

Abstract: Dielectric barrier discharges (DBDs) typically operate in the filamentary regime and thus exhibit great spatial and temporal non-uniformity. In order to optimize DBDs for various applications, such as in plasma catalysis, more fundamental insight is needed. Here, we consider how the millions of microdischarges, characteristic for a DBD, influence individual gas molecules. We use a Monte Carlo approach to determine the number of microdischarges to which a single molecule would be exposed, by means of particle tracing simulations through a full-scale packed bed DBD reactor, as well as an empty DBD reactor. We find that the fraction of microdischarges to which the molecules are exposed can be approximated as the microdischarge volume over the entire reactor gas volume. The use of this concept provides good agreement between a plasma-catalytic kinetics model and experiments for plasma-catalytic NH3 synthesis. We also show that the concept of the fraction of microdischarges indicates the efficiency by which the plasma power is transferred to the gas molecules. This generalised concept is also applicable for other spatially and temporally non-uniform plasmas.

Abstract: In this paper, some modelling approaches to describe direct current (dc) magnetron discharges developed in our research groups will be presented, including an analytical model, Monte Carlo simulations for the electrons and for the sputtered atoms, a hybrid Monte Carlo-fluid model and particle-in-cell-Monte Carlo collision simulations. The strengths and limitations of the various modelling approaches will be explained, and some characteristic simulation results will be illustrated. Furthermore, some other simulation methods related to the magnetron device will be briefly explained, more specifically for calculating the magnetic field distribution inside the discharge, and for describing the (reactive) sputtering.

Abstract: The application of atmospheric pressure plasmas in medicine is increasingly gaining attention in recent years, although very little is currently known about the plasma-induced processes occurring on the surface of living organisms. It is known that most bio-organisms, including bacteria, are coated by a liquid film surrounding them, and there might be many interactions between plasma species and the liquid layer before the plasma species reach the surface of the bio-organisms. Therefore, it is essential to study the behaviour of the reactive species in a liquid film, in order to determine whether these species can travel through this layer and reach the biomolecules, or whether new species are formed along the way. In this work, we investigate the interaction of reactive oxygen species (i.e. O, OH, HO2 and H2O2) with water, which is assumed as a simple model system for the liquid layer surrounding biomolecules. Our computational investigations show that OH, HO2 and H2O2 can travel deep into the liquid layer and are hence in principle able to reach the bio-organism. Furthermore, O, OH and HO2 radicals react with water molecules through hydrogen-abstraction reactions, whereas no H-abstraction reaction takes place in the case of H2O2. This study is important to gain insight into the fundamental operating mechanisms in plasma medicine, in general, and the interaction mechanisms of plasma species with a liquid film, in particular.

Abstract: A thorough understanding of the skyrmion motion through nanotracks is a prerequisite to realize the full potential of spintronic applications like the skyrmion racetrack memory. One of the challenges is to place the data, i.e. skyrmions, on discrete fixed positions, e.g. below a read or write head. In the domain-wall racetrack memory, one proposed solution to this problem was patterning the nanotrack with notches. Following this approach, this paper reports on the skyrmion mobility through a nanotrack with periodic notches (constrictions) made using variations in the chiral Dzyaloshinskii-Moriya interaction. We observe that such notches induce a coupling between the mobility and the skyrmion breathing mode, which manifests itself as velocity-dependent oscillations of the skyrmion diameter and plateaus in which the velocity is independent of the driving force. Despite the fact that domain walls are far more rigid objects than skyrmions, we were able to perform an analogous study and, surprisingly, found even larger plateaus of constant velocity. For both systems it is straightforward to tune the velocity at these plateaus by changing the design of the notched nanotrack geometry, e.g. by varying the distance between the notches. Therefore, the notch-induced coupling between the excited modes and the mobility could offer a strategy to stabilize the velocity against unwanted perturbations in racetrack-like applications. In the last part of the paper we focus on the low-current mobility regimes, whose very rich dynamics at nonzero temperatures are very similar to the operating principle of recently developed probabilistic logic devices. This proves that the mobility of nanomagnetic structures through a periodically modulated track is not only interesting from a fundamental point of view, but has a future in many spintronic applications.

Abstract: A hybrid model, called the hybrid plasma equipment model (HPEM), was used to study an inductively coupled plasma in gas mixtures of H2 or NH3 with CH4 or C2H2 used for the synthesis of carbon nanotubes or carbon nanofibres (CNTs/CNFs). The plasma properties are discussed for different gas mixtures at low and moderate pressures, and the growth precursors for CNTs/CNFs are analysed. It is found that C2H2, C2H4 and C2H6 are the predominant molecules in CH4 containing plasmas besides the feedstock gas, and serve as carbon sources for CNT/CNF formation. On the other hand, long-chain hydrocarbons are observed in C2H2-containing plasmas. Furthermore, the background gases CH4 and C2H2 show a different decomposition rate with H2 or NH3 addition at moderate pressures.

Abstract: In this work, a time modulated RF driven DBD-like atmospheric pressure plasma jet in Ar + 2%O2, operating at a time averaged power of 6.5 W is investigated. Spatially resolved ozone densities and gas temperatures are obtained by UV absorption and Rayleigh scattering, respectively. Significant gas heating in the core of the plasma up to 700 K is found and at the position of this increased gas temperature a depletion of the ozone density is found. The production and destruction reactions of O3 in the jet effluent as a function of the distance from the nozzle are obtained from a zero-dimensional chemical kinetics model in plug flow mode which considers relevant air chemistry due to air entrainment in the jet fluent. A comparison of the measurements and the models show that the depletion of O3 in the core of the plasma is mainly caused by an enhanced destruction of O3 due to a large atomic oxygen density.

Abstract: Packed bed dielectric barrier discharge (PB-DBD) plasma reactors are very promising for various plasma catalysis applications, but the exact mechanisms of plasma-catalyst interaction are far from understood, because the plasma discharge and catalyst/packing properties are mutually dependent. To better understand the effect of packing dielectric material on the electrical plasma properties, we study here a single bead DBD plasma reactor operating in dry air, with beads of different dielectric constant and for different applied voltages, by means of fluid modelling validated by optical imaging experiments. Our study reveals that the plasma in the single bead DBD reactor can manifest itself in two different modalities, i.e. (a) polar discharges at the bead poles in contact with the electrodes, and (b) a streamer discharge caused by surface ionization waves, which bridges the gas gap. Beads with high dielectric constant result in localised electric field enhancement and hence yield a reduction of the applied voltage required for plasma production. At low applied voltage, the discharge appears as polar discharges between the bead and the electrodes, and upon higher voltage it undergoes a transition into a bridging streamer discharge. The transition voltage to the streamer mode rises for beads with higher dielectric constant. These observations are important for plasma catalysis applications. A higher dielectric constant yields a higher electric field and thus higher average electron energy and density, giving rise to more reactive species, but it also yields a confined discharge near the contact points of packing beads, limiting the interaction area between the catalyst and the active plasma species. In addition, our model reveals that the dielectric bead behaves as a capacitor and traps charges, which can explain the significant occurrence of partial discharging in PB-DBDs and non-parallelogram shaped Lissajous plots. Hence, equivalent circuit modelling of PB-DBDs should take into account the role of packing beads in charge trapping as a capacitor.