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: A two-dimensional hybrid Monte Carlofluid model is applied to study HBr/He inductively coupled plasmas used for etching of Si. Complete sets of gas-phase and surface reactions are presented and the effects of the gas mixing ratio on the plasma characteristics and on the etch rates are discussed. A comparison with experimentally measured etch rates is made to validate the modelling results. The etch rate in the HBr plasma is found to be quite low under the investigated conditions compared to typical etch rates of Si with F- or Cl-containing gases. This allows for a higher control and fine-tuning of the etch rate when creating ultra-small features. Our calculations predict a higher electron temperature at higher He fraction, because the electrons do not lose their energy so efficiently in vibrational and rotational excitations. As a consequence, electron impact ionization and dissociation become more important, yielding higher densities of ions, electrons and H atoms. This results in more pronounced sputtering of the surface. Nevertheless, the overall etch rate decreases upon increasing He fraction, suggesting that chemical etching is still the determining factor for the overall etch rate.

Abstract: Plasma modelling is an exciting subject in which virtually all physical disciplines are represented. Plasma models combine the electromagnetic, statistical and fluid dynamical theories that have their roots in the 19th century with the modern insights concerning the structure of matter that were developed throughout the 20th century. The present cluster issue consists of 20 invited contributions, which are representative of the state of the art in plasma modelling and numerical simulation. These contributions provide an in-depth discussion of the major theories and modelling and simulation strategies, and their applications to contemporary plasma-based technologies. In this editorial review, we introduce and complement those papers by providing a bird's eye perspective on plasma modelling and discussing the historical context in which it has surfaced.

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: 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: Using a molecular dynamics model the crystallinity of MgxAlyOz thin films with a variation in the stoichiometry of the thin film is studied at operating conditions similar to the experimental operating conditions of a dual magnetron sputter deposition system. The films are deposited on a crystalline or amorphous substrate. The Mg metal content in the film ranged from 100% (i.e. MgO film) to 0% (i.e. Al2O3 film). The radial distribution function and density of the films are calculated. The results are compared with x-ray diffraction and transmission electron microscopy analyses of experimentally deposited thin films by the dual magnetron reactive sputtering process. Both simulation and experimental results show that the structure of the MgAlO film varies from crystalline to amorphous when the Mg concentration decreases. It seems that the crystalline MgAlO films have a MgO structure with Al atoms in between.

Abstract: This paper focuses on the investigation of CF4 decomposition in a low-pressure inductively coupled plasma by means of a global model. The influence of O2 on the CF4 decomposition process is studied for conditions used in semiconductor manufacturing processes. The model is applied for different powers and O2 contents ranging between 2% and 98% in the CF4/O2 gas mixture. The model includes the reaction mechanisms in the gas phase coupled with the surface reactions and sticking probabilities of the species at the walls. The calculation results are first compared with experimental results from the literature (for the electron density, temperature and F atom density) at a specific power, in the entire range of CF4/O2 gas mixture ratios, and the obtained agreements indicate the validity of the model. The main products of the gas mixture, obtained from this model, include CO, CO2 and COF2 together with a low fraction of F2. The most effective reactions for the formation and loss of the various species in this process are also determined in detail. Decomposition of CF4 produces mostly CF3 and F radicals. These radicals also contribute to the backward reactions, forming again CF4. This study reveals that the maximum decomposition efficiency of CF4 is achieved at a CF4/O2 ratio equal to 1, at the applied power of 300 W.

Abstract: High density polyethylene surfaces were exposed to the atmospheric post-discharge of a radiofrequency plasma torch supplied in helium and oxygen. Dynamic water contact angle measurements were performed to evaluate changes in surface hydrophilicity and angle resolved x-ray photoelectron spectroscopy was carried out to identify the functional groups responsible for wettability changes and to study their subsurface depth profiles, up to 9 nm in depth. The reactions leading to the formation of CO, C = O and OC = O groups were simulated by molecular dynamics. These simulations demonstrate that impinging oxygen atoms do not react immediately upon impact but rather remain at or close to the surface before eventually reacting. The simulations also explain the release of gaseous species in the ambient environment as well as the ejection of low molecular weight oxidized materials from the surface.

Abstract: A modelling study is performed of high-density low-pressure inductively coupled CF4/CHF3/H2/Cl2/O2/HBr plasmas under different gas mixing ratios. A reaction set describing the complete plasma chemistry is presented and discussed. The gas fraction of each component in this mixture is varied to investigate the sensitivity of the plasma properties, like electron density, plasma potential and species densities, towards the gas mixing ratios. This research is of great interest for microelectronics applications because these gases are often combined in two (or more)-component mixtures, and mixing gases or changing the fraction of a gas can sometimes yield unwanted reaction products or unexpected changes in the overall plasma properties due to the increased chemical complexity of the system. Increasing the CF4 fraction produces more F atoms for chemical etching as expected, but also more prominently lowers the density of Cl atoms, resulting in an actual drop in the etch rate under certain conditions. Furthermore, CF4 decreases the free electron density when mixed with Cl2. However, depending on the other gas components, CF4 gas can also sometimes enhance free electron density. This is the case when HBr is added to the mixture. The addition of H2 to the gas mixture will lower the sputtering process, not only due to the lower overall positive ion density at higher H2 fractions, but also because more H+, H2 + and H3 + are present and they have very low sputter yields. In contrast, a larger Cl2 fraction results in more chemical etching but also in less physical sputtering due to a smaller abundance of positive ions. Increasing the O2 fraction in the plasma will always lower the etch rate due to more oxidation of the wafer surface and due to a lower plasma density. However, it is also observed that the density of F atoms can actually increase with rising O2 gas fraction. This is relevant to note because the exact balance between fluorination and oxidation is important for fine-tuning the overall etch rate and for control of the sidewall profile. Finally, HBr is often used as a chemical etcher, but when mixed with F- or Cl-containing gases, HBr creates the same diluting effects as Ar or He, because a
higher fraction results in less chemical etching but more (physical) sputtering.

Abstract: Experiments show that the etch rate of Si with SF6 inductively coupled plasma (ICP) is significantly influenced by the absolute gas flow rate in the range of 50–600 sccm, with a maximum at around 200 sccm. Therefore, we numerically investigate the effects of the gas flow rate on the bulk plasma properties and on the etch rate, to obtain more insight in the underlying reasons of this effect. A hybrid Monte Carlo—fluid model is applied to simulate an SF6 ICP. It is found that the etch rate is influenced by two simultaneous effects: (i) the residence time of the gas and (ii) the temperature profile of the plasma in the ICP volume, resulting indeed in a maximum etch rate at 200 sccm.

Abstract: Here, we present routes to increase CO2 conversion into CO using an atmospheric pressure dielectric-barrier discharge. The change in conversion as a function of simple plasma parameters, such as power, flow rate, but also frequency, on-and-off power pulse, thickness and the chemical nature of the dielectric, wall and gas temperature, are described. By means of an in-depth electrical characterization of the discharge (effective plasma voltage, dielectric voltage, plasma current, number and lifetime of the microdischarges), combined with infrared analysis of the walls of the reactor, optical emission spectroscopy for the gas temperature, and mass spectrometry for the CO2 conversion, we propose a global interpretation of the effect of all the experimental parameters on the conversion and efficiency of the reaction.

Abstract: Plasma-enhanced chemical deposition (PECVD) of carbon nanotubes has already been shown to allow chirality control to some extent. In PECVD, however, etching may occur simultaneously with the growth, and the occurrence of intermediate processes further significantly complicates the growth process.
We here employ a computational approach with experimental support to study the plasma-based formation of Ni nanoclusters, Ni-catalyzed CNT growth and subsequent etching processes, in order to understand the underpinning nanoscale mechanisms. We find that hydrogen is the dominant factor in both the re-structuring of a Ni film and the subsequent appearance of Ni nanoclusters, as well as in the CNT nucleation and etching processes. The obtained results are compared with available theoretical and experimental studies and provide a deeper understanding of the occurring nanoscale mechanisms in plasma-assisted CNT nucleation and growth.