Abstract: A computer model is developed for describing argon/nitrogen glow discharges. The species taken into account in the model include electrons, Ar atoms in the ground state and in the 4s metastable levels, N2 molecules in the ground state and in six different electronically excited levels, N atoms, Ar+ ions, N+, N2+, N3+ and N4+ ions. The fast electrons are simulated with a Monte Carlo model, whereas all other species are treated in a fluid model. 74 different chemical reactions are considered in the model. The calculation results include the densities of all the different plasma species, as well as information on their production and loss processes. The effect of different N2 additions, in the range between 0.1 and 10%, is investigated.

Abstract: A Monte Carlo model is utilized for studying the behavior of electrons in the afterglow of an analytical microsecond dc pulsed glow discharge. This model uses several quantities as input data, such as electric field and potential, ion flux at the cathode, the fast argon ion and atom impact ionization rates, slow electron density, the electrical characterization of the pulse (voltage and current profiles) and temperature profile. These quantities were obtained by earlier Monte Carlo fluid calculations for a pulsed discharge. Our goal is to study the behavior of the so-called Monte Carlo electrons (i.e., those electrons created at the cathode or by ionization collisions in the plasma which are followed by using the Monte Carlo model) from their origin to the moment when they are absorbed at the cell walls or when they have lost their energy by collisions (being transferred to the group of slow electrons) in the afterglow of the pulsed discharge. The thermalization of the electrons is a phenomenon where the electron-electron Coulomb collisions acquire a special importance. Indeed, in the afterglow the cross sections of the other electron reactions taken into account in the model are very low, because of the very low electron energy. We study the electron energy distributions at several times during and after the pulse and at several positions in the plasma cell, focusing on the thermalization and on the behavior of the electrons in the afterglow. Also, the time evolution of the rates of the various collision processes, the average electron energy, the densities of Monte Carlo and slow electrons and the ionization degree are investigated.

Abstract: A model for an atmospheric pressure inductively coupled plasma (ICP) is developed which allows rather easy extension to a variable number of species and ionisation degrees. This encompasses an easy calculation of transport parameters for mixtures, ionisation and heat capacity. The ICP is modeled in an axisymmetric geometry, taking into account the gas streaming into a flowing ambient gas. A mixture of argon and helium is applied in the injector gas stream as it is often done in laser ablation ICP spectrometry. The results show a strong influence of the added helium on the center of the ICP, which is important for chemical analysis. The length of the central channel is significantly increased and the temperature inside is significantly higher than in the case of pure argon. This means that higher gas volume flow rates can be applied by addition of helium compared to the use of pure argon. This has the advantage that the gas velocity in the transport system towards the ICP can be increased, which allows shorter washout-times. Consequently, shorter measurement times can be achieved, e.g. for spatial mapping analyses in laser ablation ICP spectrometry. Furthermore, the higher temperature and the longer effective plasma length will increase the maximum size of droplets or particles injected into the ICP that are completely evaporated at the detection site. Thus, we expect an increase of the analytical performance of the ICP by helium addition to the injector gas.

Abstract: An inductively coupled plasma, connected to a mass spectrometer interface, is computationally investigated. The effect of pressure behind the sampler, injector gas flow rate, auxiliary gas flow rate, and applied power is studied. There seems to be an optimum range of injector gas flow rate for each setup which guaranties the presence and also a proper length of the central channel in the torch. Moreover, our modeling results show that for any specific purpose, it is possible to control that either only the central gas flow passes through the sampler orifice or that it is accompanied by the auxiliary gas flow. It was also found that depending on geometry, the variation of outgoing gas flow rate is much less than the variation of the injector gas flow rate and this causes a slightly higher pressure inside the torch. The general effect of increasing the applied power is a rise in the plasma temperature, which results in a higher ionization in the coil region. However, the negative effect is reducing the length of the cool central channel which is important to transfer the sample substances to the sampler. Using a proper applied power can enhance the efficiency of the system. Indeed, by changing the gas path lines, the power can control which flow (i.e., only from injector gas or also from the auxiliary gas) goes to the sampler orifice. Finally, as also reported from experiments in literature, the pressure behind the sampler has no dramatic effect on the plasma characteristics.

Abstract: In traditional capacitively coupled plasmas, the discharge can be described by an electrostatic model, in which the Poisson equation is employed to determine the electrostatic electric field. However, current plasma reactors are much larger and driven at a much higher frequency. If the excitation wavelength k in the plasma becomes comparable to the electrode radius, and the plasma skin depth d becomes comparable to the electrode spacing, the electromagnetic (EM) effects will become significant and compromise the plasma uniformity. In this regime, capacitive discharges have to be described by an EM model, i.e., the full set of Maxwells equations should be solved to address the EM effects. This paper gives an overview of the theory, simulation and experiments that have recently been carried out to understand these effects, which cause major uniformity problems in plasma processing for microelectronics and flat panel display industries. Furthermore, some methods for improving the plasma uniformity are also described and compared.

Abstract: A two-dimensional hybrid Monte Carlofluid model, incorporating a full-wave solution of Maxwell's equations, is employed to describe the behavior of high frequency (HF) and very high frequency capacitively coupled plasmas (CCPs), operating both at single frequency (SF) and dual frequency (DF) in a CF4/O2 gas mixture. First, the authors investigate the plasma composition, and the simulations reveal that besides CF4 and O2, also COF2, CF3, and CO2 are important neutral species, and CF+3 and F− are the most important positive and negative ions. Second, by comparing the results of the model with and without taking into account the electromagnetic effects for a SF CCP, it is clear that the electromagnetic effects are important, both at 27 and 60 MHz, because they affect the absolute values of the calculation results and also (to some extent) the spatial profiles, which accordingly affects the uniformity in plasma processing. In order to improve the plasma radial uniformity, which is important for the etch process, a low frequency (LF) source is added to the discharge. Therefore, in the major part of the paper, the plasma uniformity is investigated for both SF and DF CCPs, operating at a HF of 27 and 60 MHz and a LF of 2 MHz. For this purpose, the authors measure the etch rates as a function of position on the wafer in a wide range of LF powers, and the authors compare them with the calculated fluxes toward the wafer of the plasma species playing a role in the etch process, to explain the trends in the measured etch rate profiles. It is found that at a HF of 60 MHz, the uniformity of the etch rate is effectively improved by adding a LF power of 2 MHz and 300 W, while its absolute value increases by about 50%, thus a high etch rate with a uniform distribution is observed under this condition.

Abstract: To get better performance from inductively coupled plasma (ICP)-based methods, it is informative to study the properties of the ICP under different conditions. Annemie Bogaerts and Maryam Aghaei at the University of Antwerp, Belgium, are using computational modeling to examine how various properties of the ICP, such as gas flow path lines and velocity, temperature changes, and ionization effects, are affected by numerous factors, such as the gas flow rates of injector and auxiliary gas, applied power, and even the very presence of a mass spectrometry (MS) sampler. They have also applied their models to study particle transport through the ICP. Using their developed model, it is now possible to predict optimum conditions for specific analyses. Bogaerts and Aghaei spoke to us about this work.

Abstract: In this paper, the effects of prolonged bias-enhanced nucleation (prolonged BEN) on the growth mechanisms of diamond are investigated by molecular dynamics (MD) and combined MD-Metropolis Monte Carlo (MD-MMC) simulations. First, cumulative impacts of CxHy+ and Hx+ on an a-C:H/nanodiamond composite were simulated; second, nonconsecutive impacts of the dominant ions were simulated in order to understand the observed phenomena in more detail. As stated in the existing literature, the growth of diamond structures during prolonged BEN is a process that takes place below the surface of the growing film. The investigation of the penetration behavior of CxHy+ and Hx+ species shows that the carbon-containing ions remain trapped within this amorphous phase where they dominate mechanisms like precipitation of sp3 carbon clusters. The H+ ions, however, penetrate into the crystalline phase at high bias voltages (>100 V), destroying the perfect diamond structure. The experimentally measured reduction of grain sizes at high bias voltage, reported in the literature, might thus be related to penetrating H+ ions. Furthermore, the CxHy+ ions are found to be the most efficient sputtering agents, preventing the build up of defective material.

Abstract: Catalysis science and technology increased efforts recently to progress beyond conventional “thermal” catalysis and face the challenges of net-zero emissions and electrification of production. Nevertheless, a better gaps and opportunities analysis is necessary. This review analyses four emerging areas of unconventional or less- conventional catalysis which share the common aspect of using directly renewable energy sources: (i) plasma catalysis, (ii) catalysis for flow chemistry and process intensification, (iii) application of electromagnetic (EM) fields to modulate catalytic activity and (iv) nanoscale generation at the catalyst interface of a strong local EM by plasmonic effect. Plasma catalysis has demonstrated synergistic effects, where the outcome is higher than the sum of both processes alone. Still, the underlying mechanisms are complex, and synergy is not always obtained. There is a crucial need for a better understanding to (i) design catalysts tailored to the plasma environment, (ii) design plasma reactors with optimal transport of plasma species to the catalyst surface, and (iii) tune the plasma conditions so they work in optimal synergy with the catalyst. Microfluidic reactors (flow chemistry) is another emerging sector leading to the intensification of catalytic syntheses, particularly in organic chemistry. New unconventional catalysts must be designed to exploit in full the novel possibilities. With a focus on (a) continuous-flow photocatalysis, (b) electrochemical flow catalysis, (c) microwave flow catalysis and (d) ultra­ sound flow activation, a series of examples are discussed, with also indications on scale-up and process indus­ trialisation. The third area discussed regards the effect on catalytic performances of applying oriented EM fields spanning several orders of magnitude. Under well-defined conditions, gas breakdown and, in some cases, plasma formation generates activated gas phase species. The EM field-driven chemical conversion processes depend further on structured electric/magnetic catalysts, which shape the EM field in strength and direction. Different effects influencing chemical conversion have been reported, including reduced activation energy, surface charging, hot spot generation, and selective local heating. The last topic discussed is complementary to the third, focusing on the possibility of tuning the photo- and electro-catalytic properties by creating a strong localised electrical field with a plasmonic effect. The novel possibilities of hot carriers generated by the plasmonic effect are also discussed. This review thus aims to stimulate the reader to make new, creative catalysis to address the challenges of reaching a carbon-neutral world.

Abstract: In recent years, the potential use of hydrogen as a clean energy source has gained considerable attention. Especially H2 formation by Ni-catalyzed reforming of methane at elevated temperatures is an attractive process. However, a more fundamental knowledge at the atomic level is needed for a full comprehension of the reactions at the catalyst surface. In this contribution, we therefore investigate the H2 formation after CHx impacts on a Ni(1 1 1) surface in the temperature range 4001600 K, by means of reactive molecular dynamics (MD) simulations using the ReaxFF potential. While some H2 formation is already observed at the lower temperatures, substantial H2 formation is only obtained at elevated temperatures of 1400 K and above. At 1600 K, the H2 molecules are even the most frequently formed species. In direct correlation with the increasing dehydrogenation at elevated temperatures, an increased surface-to-subsurface C-diffusivity is observed as well. This study highlights the major importance of the temperature on the H2 formation.