Abstract: The electron behaviour in an Ar/CF4 inductively coupled plasma is investigated by a Langmuir probe and a hybrid model. The simulated and measured results include electron density, temperature and electron energy distribution function for different values of Ar/CF4 ratio, coil power and gas pressure. The hybrid plasma equipment model simulations show qualitative agreement with experiment. The effect of F2 electron attachment on the electron behaviour is explored by comparing two sets of data based on different F atom boundary conditions. It is demonstrated that electron attachment at F2 molecules is responsible for the depletion of low-energy electrons, causing a density decrease as well as a temperature increase when CF4 is added to an Ar plasma.

Abstract: A combined Monte Carlo (MC)/analytical surface model is developed to study the plasma processes occurring during the reactive sputter deposition of TiOx thin films. This model describes the important plasma species with a MC approach (i.e. electrons, Ar+ ions, {\rm O}_2

Abstract: In this paper, an O2 inductively coupled plasma used for plasma enhanced atomic layer deposition of Al2O3 thin films is investigated by means of modeling. This work intends to provide more information about basic plasma properties such as species densities and species fluxes to the substrate as a function of power and pressure, which might be hard to measure experimentally. For this purpose, a hybrid model developed by Kushner et al is applied to calculate the plasma characteristics in the reactor volume for different chamber pressures ranging from 1 to 10 mTorr and different coil powers ranging from 50 to 500 W. Density profiles of the various oxygen containing plasma species are reported as well as fluxes to the substrate under various operating conditions. Furthermore, different orientations of the substrate, which can be placed vertically or horizontally in the reactor, are taken into account. In addition, special attention is paid to the recombination process of atomic oxygen on the different reactor walls under the stated operating conditions. From this work it can be concluded that the plasma properties change significantly in different locations of the reactor. The plasma density near the cylindrical coil is high, while it is almost negligible in the neighborhood of the substrate. Ion and excited species fluxes to the substrate are found to be very low and negligible. Finally, the orientation of the substrate has a minor effect on the flux of O2, while it has a significant effect on the flux of O. In the horizontal configuration, the flux of atomic oxygen can be up to one order of magnitude lower than in the vertical configuration.

Abstract: Non-equilibrium inductively coupled plasmas (ICPs) operating in hydrogen are of significant interest for applications including large-area materials processing. Increasing control of spatial gas heating, which drives the formation of neutral species density gradients and the rate of gas-temperature-dependent reactions, is critical. In this study, we use 2D fluid-kinetic simulations with the Hybrid Plasma Equipment Model to investigate the spatially resolved production of atomic hydrogen in a low-pressure planar ICP operating in pure hydrogen (10-20 Pa or 0.075-0.15 Torr, 300 W). The reaction set incorporates self-consistent calculation of the spatially resolved gas temperature and 14 vibrationally excited states. We find that the formation of neutral-gas density gradients, which result from spatially non-uniform electrical power deposition at constant pressure, can drive significant variations in the vibrational distribution function and density of atomic hydrogen when gas heating is spatially resolved. This highlights the significance of spatial gas heating on the production of reactive species in relatively high-power-density plasma processing sources.

Abstract: Since the onset of pattern transfer technologies for chip manufacturing, various strategies have been developed to circumvent or overcome aspect ratio dependent etching (ARDE). These methods have, however, their own limitations in terms of etch non-idealities, throughput or costs. Moreover, they have mainly been optimized for individual in-device features and die-scale patterns, while occasionally ending up with poor patterning of metrology marks, affecting the alignment and overlay in lithography. Obtaining a better understanding of the underlying mechanisms of ARDE and how to mitigate them therefore remains a relevant challenge to date, for both marks and advanced nodes. In this work, we accordingly assessed the neutral transport mechanisms in ARDE by means of experiments and multiscale modeling for SiO₂etching with CHF₃/Ar and CF₄/Ar plasmas. The experiments revealed a local maximum in the etch rate for an aspect ratio around unity, i.e. the simultaneous occurrence of regular and inverse reactive ion etching lag for a given etch condition. We were able to reproduce this ARDE trend in the simulations without taking into account charging effects and the polymer layer thickness, suggesting shadowing and diffuse reflection of neutrals as the primary underlying mechanisms. Subsequently, we explored four methods with the simulations to regulate ARDE, by varying the incident plasma species fluxes, the amount of polymer deposition, the ion energy and angular distribution and the initial hardmask sidewall angle, for which the latter was found to be promising in particular. Although our study focusses on feature dimensions characteristic to metrology marks and back-end-of-the-line integration, the obtained insights have a broader relevance, e.g. to the patterning of advanced nodes. Additionally, this work supports the insight that physisorption may be more important in plasma etching at room temperature than originally thought, in line with other recent studies, a topic on which we recommend further research.

Abstract: Self-organized striated structure has been observed experimentally and numerically in CF4 plasmas in radio-frequency capacitively coupled plasmas recently (Liu et al 2016 Phys. Rev. Lett. 116 255002). In this work, the striated structure is investigated in a capacitively coupled oxygen discharge with the introduction of the effect from the secondary electron emission, based on a particle-in-cell/Monte Carlo collision model. As we know, the transport of positive and negative ions plays a key role in the formation of striations in electronegative gases, for which, the electronegativity needs to be large enough. As the secondary electron emission increases, electrons in the sheaths gradually contribute more ionization to the discharge. Meanwhile, the increase of the electron density, especially in the plasma bulk, leads to an increased electrical conductivity and a reduced bulk electric field, which would shield the ions' mobility. These changes result in enlarged striation gaps. And then, with more emitted electrons, obvious disruption of the striations is observed accompanied with a transition of electron heating mode. Due to the weakened field, the impact ionization in the plasma bulk is attenuated, compared with the enhanced ionization caused by secondary electrons. This would lead to the electron heating mode transition from striated (STR) mode to gamma-mode. Besides, our investigation further reveals that gamma-mode is more likely to dominate the discharge under high gas pressures or driving voltages.

Abstract: Plasma generation inside catalyst pores is of utmost importance for plasma catalysis, as the existence of plasma species inside the pores affects the active surface area of the catalyst available to the plasma species for catalytic reactions. In this paper, the electric field enhancement, and thus the plasma production inside catalyst pores with different pore shapes is studied with a two-dimensional fluid model. The results indicate that the electric field will be significantly enhanced near tip-like structures. In a conical pore with small opening, the strongest electric field appears at the opening and bottom corners of the pore, giving rise to a prominent ionization rate throughout the pore. For a cylindrical pore, the electric field is only enhanced at the bottom corners of the pore, with lower absolute value, and thus the ionization rate inside the pore is only slightly enhanced. Finally, in a conical pore with large opening, the electric field is characterized by a maximum at the bottom of the pore, yielding a similar behavior for the ionization rate. These results demonstrate that the shape of the pore has a significantly influence on the electric field enhancement, and thus modifies the plasma properties.

Abstract: This work explains the need for plasma models, introduces arguments for choosing the type of model that better fits the purpose of each study, and presents the basics of the most common nonequilibrium low-temperature plasma models and the information available from each one, along with an extensive list of references for complementary in-depth reading. The paper presents the following models, organised according to the level of multi-dimensional description of the plasma: kinetic models, based on either a statistical particle-in-cell/Monte-Carlo approach or the solution to the Boltzmann equation (in the latter case, special focus is given to the description of the electron kinetics); multi-fluid models, based on the solution to the hydrodynamic equations; global (spatially-average) models, based on the solution to the particle and energy rate-balance equations for the main plasma species, usually including a very complete reaction chemistry; mesoscopic models for plasma–surface interaction, adopting either a deterministic approach or a stochastic dynamical Monte-Carlo approach. For each plasma model, the paper puts forward the physics context, introduces the fundamental equations, presents advantages and limitations, also from a numerical perspective, and illustrates its application with some examples. Whenever pertinent, the interconnection between models is also discussed, in view of multi-scale hybrid approaches.

Abstract: A gliding arc plasmatron (GAP) is very promising for CO2 conversion into value-added chemicals, but to further improve this important application, a better understanding of the arc behavior is indispensable. Therefore, we study here for the first time the dynamic arc behavior of the GAP by means of a high-speed camera, for different reactor configurations and in a wide range of operating conditions. This allows us to provide a complete image of the behavior of the gliding arc. More specifically, the arc body shape, diameter, movement and rotation speed are analyzed and discussed. Clearly, the arc movement and shape relies on a number of factors, such as gas turbulence, outlet diameter, electrode surface, gas contraction and buoyance force. Furthermore, we also compare the experimentally measured arc movement to a state-of-the-art 3D-plasma model, which predicts the plasma movement and rotation speed with very good accuracy, to gain further insight in the underlying mechanisms. Finally, we correlate the arc dynamics with the CO2 conversion and energy efficiency, at exactly the same conditions, to explain the effect of these parameters on the CO2 conversion process. This work is important for understanding and optimizing the GAP for CO2 conversion.

Abstract: Packed bed plasma reactors (PBPRs) are gaining increasing interest for use in environmental applications, such as greenhouse gas conversion into value-added chemicals or renewable fuels and volatile pollutant removal (e.g. NOx, VOC, K), as they enhance the conversion and energy efficiency of the process compared to a non-packed reactor. However, the plasma behaviour in a PBPR is not well understood. In this paper we demonstrate, by means of a fluid model, that the discharge behaviour changes considerably when changing the size of the packing beads and their dielectric constant, while keeping the interelectrode spacing constant. At low dielectric constant, the plasma is spread out over the full discharge gap, showing significant density in the voids as well as in the connecting void channels. The electric current profile shows a strong peak during each half cycle. When the dielectric constant increases, the plasma becomes localised in the voids, with a current profile consisting of many smaller peaks during each half cycle. For large bead sizes, the shift from full gap discharge to localised discharges takes place at a higher dielectric constant than for smaller beads. Furthermore, smaller beads or beads with a lower dielectric constant require a higher breakdown voltage to cause plasma formation.

Abstract: In recent years there has been growing interest in the use of plasma technology for CO2 conversion. To improve this application, a good insight into the underlying mechanisms is of great importance. This can be obtained from modeling the detailed plasma chemistry in order to understand the chemical reaction pathways leading to CO2 conversion (either in pure form or mixed with another gas). Moreover, in practice, several plasma reactor types are being investigated for CO2 conversion, so in addition it is essential to be able to model these reactor geometries so that their design can be improved, and the most energy efficient CO2 conversion can be achieved. Modeling the detailed plasma chemistry of CO2 conversion in complex reactors is, however, very time-consuming. This problem can be overcome by using a combination of two different types of model: 0D chemical reaction kinetics models are very suitable for describing the detailed plasma chemistry, while the characteristic features of different reactor geometries can be studied by 2D or 3D fluid models. In the first instance the latter can be developed in argon or helium with a simple chemistry to limit the calculation time; however, the ultimate aim is to implement the more complex CO2 chemistry in these models. In the present paper, examples will be given of both the 0D plasma chemistry models and the 2D and 3D fluid models for the most common plasma reactors used for CO2 conversion in order to emphasize the complementarity of both approaches. Furthermore, based on the modeling insights, the paper discusses the possibilities and limitations of plasma-based CO2 conversion in different types of plasma reactors, as well as what is needed to make further progress in this field.

Abstract: One of the most challenging and recurring problems when modeling plasmas is the lack of data on the key atomic and molecular reactions that drive plasma processes. Even when there are data for some reactions, complete and validated datasets of chemistries are rarely available. This hinders research on plasma processes and curbs development of industrial applications. The QDB project aims to address this problem by providing a platform for provision, exchange, and validation of chemistry datasets. A new data model developed for QDB is presented. QDB collates published data on both electron scattering and heavy-particle reactions. These data are formed into reaction sets, which are then validated against experimental data where possible. This process produces both complete chemistry sets and identifies key reactions that are currently unreported in the literature. Gaps in the datasets can be filled using established theoretical methods. Initial validated chemistry sets for SF6/CF4/O2 and SF6/CF4/N2/H2 are presented as examples.