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: 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: Plasma–liquid interactions represent a growing interdisciplinary area of research involving plasma science, fluid dynamics, heat and mass transfer, photolysis, multiphase chemistry and aerosol science. This review provides an assessment of the state-of-the-art of this multidisciplinary area and identifies the key research challenges. The developments in diagnostics, modeling and further extensions of cross section and reaction rate databases that are necessary to address these challenges are discussed. The review focusses on nonequilibrium plasmas.

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: The chemical complexity of non-equilibrium plasmas poses a challenge for plasma modeling because of the computational load. This paper presents a dimension reduction method for such chemically complex plasmas based on principal component analysis (PCA). PCA is used to identify a low-dimensional manifold in chemical state space that is described by a small number of parameters: the principal components. Reduction is obtained since continuity equations only need to be solved for these principal components and not for all the species. Application of the presented method to a CO2 plasma model including state-to-state vibrational kinetics of CO2 and CO demonstrates the potential of the PCA method for dimension reduction. A manifold described by only two principal components is able to predict the CO2 to CO conversion at varying ionization degrees very accurately.

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: The formation process of a microdischarge (MD) in both μm- and nm-sized catalyst pores is simulated by a two-dimensional particle-in-cell/Monte Carlo collision model. A parallel-plate dielectric barrier discharge configuration in filamentary mode is considered in ambient air. The discharge is powered by a high voltage pulse. Our calculations reveal that a streamer can penetrate into the surface features of a porous catalyst and MDs can be formed inside both μm- and nm-sized pores, yielding ionization inside the pore. For the μm-sized pores, the ionization mainly occurs inside the pore, while for the nm-sized pores the ionization is strongest near and inside the pore. Thus, enhanced discharges near and inside the mesoporous catalyst are observed. Indeed, the maximum values of the electric field, ionization rate and electron density occur near and inside the pore. The maximum electric field and electron density inside the pore first increase when the pore size rises from 4 nm to 10 nm, and then they decrease for the 100 nm pore, due to
a more pronounced surface discharge for the smaller pores. However, the ionization rate is highest for the 100 nm pore due to the largest effective ionization region.

Abstract: This work examines the properties of a dielectric barrier discharge (DBD) reactor, built for CO2 decomposition, by means of electrical characterization, optical emission spectroscopy and gas chromatography. The discharge, formed in an electronegative gas (such as CO2, but also O2), exhibits clearly different electrical characteristics, depending on the surface conductivity of the reactor walls. An asymmetric current waveform is observed in the metaldielectric (MD) configuration, with sparse high-current pulses in the positive half-cycle (HC) and a more uniform regime in the negative HC. This indicates that the discharge is operating in two alternating regimes with rather different properties. At high CO2 conversion regimes, a conductive coating is deposited on the dielectric. This so-called coated MD configuration yields a symmetric current waveform, with current peaks in both the positive and negative HCs. In a double-dielectric (DD) configuration, the current waveform is also symmetric, but without current peaks in both the positive and negative HC. Finally, the DD configuration with conductive coating on the inner surface of the outer dielectric, i.e. so-called coated DD, yields again an asymmetric current waveform, with current peaks in the negative HC. These different electrical characteristics are related to the presence of the conductive coating on the dielectric wall of the reactor and can be explained by an increase of the local barrier capacitance available for charge transfer. The different discharge regimes affect the CO2 conversion, more specifically, the CO2 conversion is lowest in the clean DD configuration. It is somewhat higher in the coated DD configuration, and still higher in the MD configuration. The clean and coated MD configuration, however, gave similar CO2 conversion. These results indicate that the conductivity of the dielectric reactor walls can highly promote the development of the high-amplitude discharge current pulses and subsequently the CO2 conversion.

Abstract: Understanding the nature and effect of the multitude of plasma–surface interactions in plasma catalysis is a crucial requirement for further process development and improvement. A particularly intriguing and rather unique property of a plasma-catalytic setup is the ability of the plasma to modify the electronic structure, and hence chemical properties, of the catalyst through charging, i.e. the absorption of excess electrons. In this work, we develop a quantum chemical model based on density functional theory to study excess negative surface charges in a heterogeneous catalyst exposed to a plasma. This method is specifically applied to investigate plasma-catalytic CO2 activation on supported M/Al2O3 (M=Ti, Ni, Cu) single atom catalysts. We find that (1) the presence of a negative surface charge dramatically improves the reductive power of the catalyst, strongly promoting the splitting of CO2 to CO and oxygen, and (2) the relative activity of the investigated transition metals is also changed upon charging, suggesting that controlled surface charging is a powerful additional parameter to tune catalyst activity and selectivity. These results strongly point to plasma-induced surface charging of the catalyst as an important factor contributing to the plasma-catalyst synergistic effects frequently reported for plasma catalysis.

Abstract: We investigate discharge contraction in a microwave plasma at sub-atmospheric pressure, operating in CO2 and CO2/CH4 mixtures. The rise of the electron number density with plasma contraction intensifies the gas heating in the core of the plasma. This, in turn, initiates fast core-periphery transport and defines the rate of thermal chemistry over plasma chemistry. In this context, power concentration describes the overall mechanism including plasma contraction and chemical kinetics. In a complex chemistry such as dry reforming of methane, transport of reactive species is essential to define the performance of the reactor and achieve the desired outputs. Thus, we couple experimental observations and thermodynamic calculations for model validation and understanding of reactor performance. Adding CH4 alters the thermodynamic properties of the mixture, especially the reactive component of the heat conductivity. The increase in reactive heat conductivity increases the pressure at which plasma contraction occurs, because higher rates of gas heating are required to reach the same temperature. In addition, we suggest that the predominance of heat conduction over convection is a key condition to observe the effect of heat conductivity on gas temperature.

Abstract: We present a zero-dimensional kinetic model of CO2 splitting in non-equilibrium plasmas. The model includes a description of the CO2 vibrational kinetics (25 vibrational levels up to the dissociation limit of the molecule), taking into account state-specific VT and VV relaxation reactions and the effect of vibrational excitation on other chemical reactions. The model is applied to study the reaction kinetics of CO2 splitting in an atmospheric-pressure dielectric barrier discharge (DBD) and in a moderate-pressure microwave discharge. The model results are in qualitative agreement with published experimental works. We show that the CO2 conversion and its energy efficiency are very different in these two types of discharges, which reflects the important dissociation mechanisms involved. In the microwave discharge, excitation of the vibrational levels promotes efficient dissociation when the specific energy input is higher than a critical value (2.0 eV/molecule under the conditions examined). The calculated energy efficiency of the process has a maximum of 23%. In the DBD, vibrationally excited levels do not contribute significantly to the dissociation of CO2 and the calculated energy efficiency of the process is much lower (5%).