Abstract: Plasma catalysis is an emerging new technology for the electrification and downscaling of NH3 synthesis. Increasing attention is being paid to the optimization of plasma catalysis with respect to the plasma conditions, the catalyst material, and their mutual interaction. In this work we use microkinetic models to study how the total conversion process is impacted by the combination of different plasma conditions and transition metal catalysts. We study how plasma-generated radicals and vibrationally excited N2 (present in a dielectric barrier discharge plasma) interact with the catalyst and impact the NH3 turnover frequencies (TOFs). Both filamentary and uniform plasmas are studied, based on plasma chemistry models that provided plasma phase speciation and vibrational distribution functions. The Langmuir−Hinshelwood reaction rate coefficients (i.e., adsorption reactions and subsequent reactions among adsorbates) are determined using conventional scaling relations. An additional set of Eley−Rideal reactions (i.e., direct reactions of plasma radicals with adsorbates) was added and a sensitivity analysis on the assumed reaction rate coefficients was performed. We first show the impact of different vibrational distribution functions on the catalytic dissociation of N2 and subsequent production of NH3, and we gradually include more radical reactions, to illustrate the contribution of these species and their corresponding reaction pathways. Analysis over a large range of catalysts indicates that different transition metals (metals such as Rh, Ni, Pt, and Pd) optimize the NH3TOFs depending on the population of the vibrational levels of N2. At higher concentrations of plasma-generated radicals, the NH3 TOFs become less dependent on the catalyst material, due to radical adsorptions on the more noble catalysts and Eley−Rideal reactions on the less noble catalysts.

Abstract: The 2022 Roadmap is the next update in the series of Plasma Roadmaps published by<italic>Journal of Physics</italic>D with the intent to identify important outstanding challenges in the field of low-temperature plasma (LTP) physics and technology. The format of the Roadmap is the same as the previous Roadmaps representing the visions of 41 leading experts representing 21 countries and five continents in the various sub-fields of LTP science and technology. In recognition of the evolution in the field, several new topics have been introduced or given more prominence. These new topics and emphasis highlight increased interests in plasma-enabled additive manufacturing, soft materials, electrification of chemical conversions, plasma propulsion, extreme plasma regimes, plasmas in hypersonics, data-driven plasma science and technology and the contribution of LTP to combat COVID-19. In the last few decades, LTP science and technology has made a tremendously positive impact on our society. It is our hope that this roadmap will help continue this excellent track record over the next 5–10 years.

Abstract: We present a combined study of experimental and computational work for a dielectric barrier discharge (DBD) used for CH4 conversion into H2. More specifically, we investigated the influence of N2 as an impurity (150,000 ppm) and as additive gas (199%) on the CH4 conversion and H2 yield. For this purpose, a zero-dimensional chemical kinetics model is applied to study the plasma chemistry. The calculated conversions and yields for various gas mixing ratios are compared to the obtained experimental values, and good agreement is achieved. The study reveals the significance of the View the MathML source and View the MathML source metastable states for the CH4 conversion into H2, based on a kinetic analysis of the reaction chemistry.

Abstract: Cold atmospheric plasmas (CAPs) have attracted significant interest for their potential benefits in medical applications, including cancer therapy. The therapeutic effects of CAPs are related to reactive oxygen and nitrogen species (ROS and RNS) present in the plasma. The impact of ROS has been extensively studied, but the role of RNS in CAP-treatment remains poorly understood at the molecular level. Here, we investigate the permeation of RNS and ROS across native and oxidized phospholipid bilayers (PLBs) by means of computer simulations. The results reveal significantly lower free energy barriers for RNS (i.e. NO, NO2, N2O4) and O3 compared to hydrophilic ROS, such as OH, HO2 and H2O2. This suggests that the investigated RNS and O3 can permeate more easily through both native and oxidized PLBs in comparison to hydrophilic ROS, indicating their potentially important role in plasma medicine.

Abstract: The mechanisms of plasma in medicine are broadly attributed to plasma-derived reactive oxygen and nitrogen species (RONS). In order to exert any intracellular effects, these plasma-derived RONS must first traverse a major barrier in the cell membrane. The cell membrane lipid composition, and thereby the magnitude of this barrier, is highly variable between cells depending on type and state (e.g. it is widely accepted that healthy and cancerous cells have different membrane lipid compositions). In this study, we investigate how plasma-derived RONS interactions with lipid membrane components can potentially be exploited in the future for treatment of diseases. We couple phospholipid vesicle experiments, used as simple cell models, with molecular dynamics (MD) simulations of the lipid membrane to provide new insights into how the interplay between phospholipids and cholesterol may influence the response of healthy and diseased cell membranes to plasma-derived RONS. We focus on the (i) lipid tail saturation degree, (ii) lipid head group type, and (iii) membrane cholesterol fraction. Using encapsulated molecular probes, we study the influence of the above membrane components on the ingress of RONS into the vesicles, and subsequent DNA damage. Our results indicate that all of the above membrane components can enhance or suppress RONS uptake, depending on their relative concentration within the membrane. Further, we show that higher RONS uptake into the vesicles does not always correlate with increased DNA damage, which is attributed to ROS reactivity and lifetime. The MD simulations indicate the multifactorial chemical and physical processes at play, including (i) lipid oxidation, (ii) lipid packing, and (iii) lipid rafts formation. The methods and findings presented here provide a platform of knowledge that could be leveraged in the development of therapies relying on the action of plasma, in which the cell membrane and oxidative stress response in cells is targeted.

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: In this tutorial review paper, we illustrate how computer modeling can contribute to a better insight in inductively coupled plasma-mass spectrometry (ICP-MS). We start with a brief overview on previous efforts, studying the fundamentals of the ICP and ICP-MS, with main focus on previous modeling activities. Subsequently, we explain in detail the model that we developed in previous years, and we show typical calculation results, illustrating the plasma characteristics, gas flow patterns and the sample transport, evaporation and ionization. We also present the effect of various experimental parameters, such as operating conditions, geometrical aspects and sample characteristics, to illustrate how modeling can help to elucidate the optimal conditions for improved analytical performance.

Abstract: We studied the transport of copper droplets through an inductively coupled plasma, connected to the sampling cone of a mass spectrometer, by means of a computational model. The sample droplets are followed until they become evaporated. They are inserted as liquid particles from the central inlet and the effects of injection position (i.e. “on” and “off” axis), droplet diameter, as well as mass loading flow rate are investigated. It is shown that more “on-axis” injection of the droplets leads to a more straight path line, so that the droplets move less in the radial direction and are evaporated more on the central axis, enabling a better sample transfer efficiency to the sampler cone. Furthermore, there are optimum ranges of diameters and flow rates, which guarantee the proper position of evaporation along the torch, i.e. not too early, so that the sample can get lost in the torch, and not too late, which reduces the chance of becoming ionized before reaching the sampler.

Abstract: We present here some of our modeling efforts for PECVD growth of nanostructured carbon materials with focus on amorphous hydrogenated carbon. Experimental data from an expanding thermal plasma setup were used as input for the simulations. Attention was focused both on the film growth mechanism, as well as on the hydrocarbon reaction mechanisms during growth of the films. It is found that the reaction mechanisms and sticking coefficients are dependent on the specific surface sites, and the structural properties of the growth radicals. The film growth results are in correspondence with the experiment. Furthermore, it is found that thin a-C:H films can be densified using an additional H-flux towards the substrate.

Abstract: The aim of this work consists of the evaluation of atmospheric pressure dielectric barrier discharges for the conversion of greenhouse gases into useful compounds. Therefore, pure CO2 feed flows are administered to the discharge zone at varying discharge frequency, power input, gas temperature and feed flow rates, aiming at the formation of CO and O2. The discharge obtained in CO2 is characterized as a filamentary mode with a microdischarge zone in each half cycle of the applied voltage. It is shown that the most important parameter affecting the CO2-conversion levels is the gas flow rate. At low flow rates, both the conversion and the CO-yield are significantly higher. In addition, also an increase in the gas temperature and the power input give rise to higher conversion levels, although the effect on the CO-yield is limited. The optimum discharge frequency depends on the power input level and it cannot be unambiguously stated that higher frequencies give rise to increased conversion levels. A maximum CO2 conversion of 30% is achieved at a flow rate of 0.05 L min−1, a power density of 14.75 W cm−3 and a frequency of 60 kHz. The most energy efficient conversions are achieved at a flow rate of 0.2 L min−1, a power density of 11 W cm−3 and a discharge frequency of 30 kHz.

Abstract: A fluid model is self-consistently established to investigate the harmonic effects in an inductively coupled plasma, where the electromagnetic field is solved by the finite difference time domain technique. The spatiotemporal distribution of harmonic current density, harmonic potential, and other plasma quantities, such as radio frequency power deposition, plasma density, and electron temperature, have been investigated. Distinct differences in current density have been observed when calculated with and without Lorentz force, which indicates that the nonlinear Lorentz force plays an important role in the harmonic effects, especially at low frequencies. Moreover, the even harmonics are larger than the odd harmonics both in the current density and the potential. Finally, the dependence of various plasma quantities with and without the Lorentz force on various driving frequencies is also examined. It is shown that the deposited power density decreases and the depth of penetration increases slightly because of the Lorentz force. The electron density increases distinctly while the electron temperature remains almost the same when the Lorentz force is taken into account.