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: A one-dimensional, implicit particle-in-cell Monte Carlo collision model is used to simulate the plasma kinetic properties at a steady state in a parallel-plate direct current argon glow microdischarge under various operating conditions, such as driving voltage (301000 V) and gap size (101000 µm) at atmospheric pressure. First, a comparison between rf and dc modes is shown for the same pressure, driving voltage and gap spacing. Furthermore, the effect of gap size scaling (in the range of 101000 µm) on the breakdown voltage, peak electron density and peak electron current density at the breakdown voltage is examined. The breakdown voltage is lower than 150 V in all gaps considered. The microdischarge is found to have a neutral bulk plasma region and a cathode sheath region with size varying with the applied voltage and the discharge gap. In our calculations, the electron and ion densities are of the order of 10181023 m−3, which is in the glow discharge limit, as the ionization degree is lower than 1% . The electron energy distribution function shows a two-energy group distribution at a gap of 10 µm and a three-energy group distribution at larger gaps such as 200 µm and 1000 µm, emphasizing the importance of the gap spacing in dc microdischarges.

Abstract: Very recently, the 2D form of MoSi2N4 has been successfully fabricated (Hong et al 2020 Science 369 670). Motivated by these recent experimental results, we investigate the structural, mechanical, thermal, electronic and optical properties of the MoSi2N4 monolayer. The mechanical study confirms the stability of the MoSi2N4 monolayer. The Young's modulus decreases by similar to 30%, while the Poisson's ratio increases by similar to 30% compared to the corresponding values of graphene. In addition, the MoSi2N4 monolayer's work function is very similar to that of phosphorene and MoS2 monolayers. The electronic structure shows that the MoSi2N4 monolayer is an indirect semiconductor with a band gaps of 1.79 (2.35) eV using the GGA (HSE06) functional. The thermoelectric performance of the MoSi2N4 monolayer has been revealed and a figure of merit slightly larger than unity at high temperatures is calculated. The optical analysis shows that the first absorption peak for in-plane polarization is located in the visible range of the spectrum, therefore, the MoSi2N4 monolayer is a promising candidate for advanced optoelectronic nanodevices. In summary, the fascinating MoSi2N4 monoloayer is a promising 2D material for many applications due to its unique physical properties.

Abstract: In this paper, the initial mechanisms of nanoparticle formation and growth in radiofrequency acetylene (C2H2) plasmas are investigated by means of a comprehensive self-consistent one-dimensional (1D) fluid model. This model is an extension of the 1D fluid model, developed earlier by De Bleecker et al. Based on the comparison of our previous results with available experimental data for acetylene plasmas in the literature, some new mechanisms for negative ion formation and growth are proposed. Possible routes are considered for the formation of larger (linear and branched) hydrocarbons C2nH2 (n = 3, 4, 5), which contribute to the generation of C2nH− anions (n = 3, 4, 5) due to dissociative electron attachment. Moreover, the vinylidene anion (H2CC−) and higher anions (n = 24) are found to be important plasma species.

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: A hybrid model, called the hybrid plasma equipment model, was used to study inductively coupled SF6 plasmas used for Si etching applications. The plasma properties such as number densities of electrons, positive and negative ions, and neutrals are calculated under typical etching conditions. The electron kinetics is analysed by means of the electron energy probability function. The plasma chemistry taking place in pure SF6 and in an Ar/SF6 mixture is also discussed, and finally the effect of the argon fraction on the plasma properties is investigated.

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 this paper, an overview of the fabrication and properties of high-quality La0.67Sr0.33MnO3 (LSMO) thin films is given. A high-quality LSMO film combines a smooth surface morphology with a large magnetization and a small residual resistivity, while avoiding precipitates and surface segregation. In the literature, typically only a few of these issues are adressed. We therefore present a thorough characterization of our films, which were grown by pulsed laser deposition. The films were characterized with reflection high energy electron diffraction, atomic force microscopy, x-ray diffraction, magnetization and transport measurements, x-ray photoelectron spectroscopy and scanning transmission electron microscopy. The films have a saturation magnetization of 4.0 µB/Mn, a Curie temperature of 350 K and a residual resistivity of 60 µΩ cm. These results indicate that high-quality films, combining both large magnetization and small residual resistivity, were realized. A comparison between different samples presented in the literature shows that focussing on a single property is insufficient for the optimization of the deposition process. For high-quality films, all properties have to be adressed. For LSMO devices, the thin-film quality is crucial for the device performance. Therefore, this research is important for the application of LSMO in devices.

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.

Abstract: The rate, selectivity and efficiency of plasma-based conversion processes is strongly affected by nonequilibrium phenomena. High concentrations of vibrationally excited molecules are such a plasma-induced effect. It is frequently assumed that vibrationally excited molecules are important in plasma catalysis because their presence lowers the apparent activation energy of dissociative chemisorption reactions and thus increases the conversion rate. A detailed atomic-level understanding of vibrationally stimulated catalytic reactions in the context of plasma catalysis is however lacking. Here, we couple a recently developed statistical model of a plasma-induced vibrational nonequilibrium to molecular dynamics simulations, enhanced sampling methods, and machine learning techniques. We quantify the impact of a vibrational nonequilibrium on the dissociative chemisorption barrier of H2 and CH4 on nickel catalysts over a wide range of vibrational temperatures. We investigate the effect of surface structure and compare the role of different vibrational modes of methane in the dissociation process. For low vibrational temperatures, very high vibrational efficacies are found, and energy in bend vibrations appears to dominate the dissociation of methane. The relative impact of vibrational nonequilibrium is much higher on terrace sites than on surface steps. We then show how our simulations can help to interpret recent experimental results, and suggest new paths to a better understanding of plasma catalysis.

Abstract: In this work, the role of vibrationally excited HBr (HBr(vib)) is computationally investigated for a HBr/He inductively coupled plasma applied for Si etching. It is found that at least 50% of all dissociations of HBr occur through HBr(vib). This additional dissociation pathway through HBr(vib) makes the plasma significantly more atomic. It also results in a slightly higher electron temperature (i.e. about 0.2 eV higher compared to simulation results where HBr(vib) is not included), as well as a higher gas temperature (i.e. about 50 K higher than without including HBr(vib)), due to the enhanced Franck–Condon heating through HBr(vib) dissociation,
at the conditions investigated. Most importantly, the calculated etch rate with HBr(vib) included in the model is a factor 3 higher than in the case without HBr(vib), due to the higher fluxes of etching species (i.e. H and Br), while the chemical composition of the wafer surface shows no significant difference. Our calculations clearly show the importance of including HBr(vib) for accurate modeling of HBr-containing plasmas.

Abstract: Both geometrical and electrical asymmetry effects in capacitive argon discharges are investigated using a two-dimensional particle-in-cell coupled with Monte Carlo collision model. When changing the ratio of the top and bottom electrode surface areas and the phase shift between the two applied harmonics, the induced self-bias was found to develop separately. By adjusting the ratio between the high and low harmonic amplitudes, the electrical asymmetry effect at a fixed phase shift can be substantially optimized. However, the self-bias caused by the geometrical asymmetry hardly changed. Moreover, the separate control of these two asymmetry effects can also be demonstrated from their power absorption profiles. Both the axial and radial plasma density distributions can be modulated by the electrical asymmetry effect.

Abstract: As electron energy-loss spectroscopy (EELS) and x-ray absorption spectroscopy (XAS) probe the same transitions from coreshell states to unoccupied states above the Fermi energy, it should always be possible to apply the two techniques to the same physical phenomena, such as magnetic dichroism, and obtain the same information. Indeed, the similarity in the expression of the electron and x-ray cross-sections had been already exploited to prove the equivalence of x-ray magnetic linear dichroism and anisotropy in EELS, by noting that the polarization vector of a photon plays the same role as the momentum transfer in electron scattering. Recently, the same was proven true for x-ray magnetic circular dichroism (XMCD) by establishing a new TEM technique called EMCD (electron energy-loss magnetic chiral dichroism) (Schattschneider P et al 2006 Nature 441 4868), which makes use of special electron scattering conditions to force the absorption of a circularly polarized virtual photon. The intrinsic advantage of EMCD over XMCD is the high spatial resolution of electron microscopes, which are readily available. Among the particular obstacles in EMCD that do not exist for synchrotron radiation, is the notoriously low signal and the very particular scattering conditions necessary to observe a chiral dichroic signal. In spite of that, impressive progress has been made in recent years. The signal strength could be considerably increased, and some innovations such as using a convergent beam have been introduced. EMCD has evolved into several techniques, which make full use of the versatility of the TEM and energy filtering, spectroscopy or STEM conditions (Rubino S 2007 Magnetic circular dichroism in the transmission electron microscope PhD Thesis Vienna University of Technology, Vienna, Austria).

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.