Abstract: We use microkinetic modeling to examine the potential of plasma-catalytic partial oxidation (POX) of CH4 as a promising new approach to produce oxygenates. We study how different plasma species affect POX of CH4 on the Pt(111) surface, and we discuss the associated kinetic and mechanistic changes. We discuss the effect of vibrationally excited CH4 and O2, as well as plasma-generated radicals and stable intermediates. Our results show that vibrational excitation enhances the turnover frequency (TOF) of catalytic CH4 dissociation and has good potential for improving the selectivities toward CH3OH, HCOOH, and C2 hydrocarbons. Nevertheless, when also considering plasma-generated radicals, we find that these species mainly govern the surface chemistry. Additionally, we find that plasma-generated radicals and stable intermediates enhance the TOFs of COx and oxygenates, increase the selectivity toward oxygenates, and make the formation of HCOOH more significant on Pt(111). We also briefly illustrate the potential impact of Eley−Rideal reactions that involve plasma-generated radicals. Finally, we reveal how various radicals affect the catalyst surface chemistry and we link this to the formation of different products. This allows us to make suggestions on how the plasma composition should be altered to improve the formation of desired products.

Abstract: To gain full control over the growth of carbon nanotubes (CNTs) using plasma-enhanced chemical vapor deposition (PECVD), a thorough understanding of the underlying plasma-catalyst mechanisms is required. Oxygen-containing species are often used as or added to the growth precursor gas, but these species also yield various radicals and ions, which may simultaneously etch the CNT during the growth. At present, the effect of these reactive species on the growth onset has not yet been thoroughly investigated. We here report on the etching mechanism of incipient CNT structures from OH and O radicals as derived from combined (reactive) molecular dynamics (MD) and force-bias Monte Carlo (tfMC) simulations. Our results indicate that the oxygen-containing radicals initiate a dissociation process. In particular, we show how the oxygen species weaken the interaction between the CNT and the nanocluster. As a result of this weakened interaction, the CNT closes off and dissociates from the cluster in the form of a fullerene. Beyond the specific systems studied in this work, these results are generically important in the context of PECVD-based growth of CNTs using oxygen-containing precursors.

Abstract: We present an approach to incorporate the effect of vibrational nonequilibrium in molecular dynamics (MD) simulations. A perturbed canonical ensemble, in which selected modes are excited to higher temperature while all others remain equilibrated at low temperature, is simulated by applying a specifically tailored bias potential. Our method can be readily applied to any (classical or quantum mechanical) MD setup at virtually no additional computational cost and allows the study of reactions of vibrationally excited molecules in nonequilibrium environments such as plasmas. In combination with enhanced sampling methods, the vibrational efficacy and mode selectivity of vibrationally stimulated reactions can then be quantified in terms of chemically relevant observables, such as reaction rates and apparent free energy barriers. We first validate our method for the prototypical hydrogen exchange reaction and then show how it can capture the effect of vibrational excitation on a symmetric SN2 reaction and radical addition on CO2.

Abstract: Theoretical studies on molecule-metal surface reactions have so far been limited to small surface unit cells due to computational costs. Here, for the first time molecular dynamics simulations on very large surface unit cells at the level of density functional theory are performed, allowing a direct comparison to experiments performed on a curved crystal. Specifically, the reaction of D-2 on a curved Pt crystal is investigated with a neural network potential (NNP). The developed NNP is also accurate for surface unit cells considerably larger than those that have been included in the training data, allowing dynamical simulations on very large surface unit cells that otherwise would have been intractable. Important and complex aspects of the reaction mechanism are discovered such as diffusion and a shadow effect of the step. Furthermore, conclusions from simulations on smaller surface unit cells cannot always be transfered to larger surface unit cells, limiting the applicability of theoretical studies of smaller surface unit cells to heterogeneous catalysts with small defect densities.

Abstract: We present an ab initio study of the electronic structure and of the formation energies of various point defects in BaSnO3 and SrGeO3. We show that La and Y impurities substituting Ba or Sr are shallow donors with a preferred 1 + charge state. These defects have a low formation energy within all the suitable equilibrium growth conditions considered. Oxygen vacancies behave as shallow donors as well, preferring the 2 + charge state. Their formation energies, however, are higher in most growth conditions, indicating a limited contribution to conductivity. The calculated electron effective mass in BaSnO3, with a value of 0.21 me, and the very high mobility reported recently in La-doped BaSnO3 single-crystals, suggest that remarkably low scattering rates can be achieved in the latter. In the case of SrGeO3, our results point to carrier density and mobility values in the low range for typical polycrystalline TCOs, in line with experiment.

Abstract: Our aim is to develop and optimize a model for a dielectric barrier discharge used as a chemical reactor for gas treatment. In order to determine the optimum operating conditions, we have studied the influence of the gas flow rate, reactor geometry and applied voltage parameters on the discharge characteristics. For this purpose, a two-dimensional time-dependent fluid model has been applied to an atmospheric pressure DBD in helium with nitrogen impurities, in a cylindrical geometry. The numerical model is based on the continuity and flux equations for each type of particles treated, the electron energy equation and the Poisson equation. The gas flow is incorporated in the flux equations as a source term. The set of coupled partial differential equations is solved by the so-called modified strongly implicit method. The background gas flow is numerically treated separately, assuming in the model that there is no influence of the plasma on the flow. Indeed, the gas convection velocity is calculated using the commercial code Fluent and it is used as input into the 2D fluid model. The plasma characteristics have been studied in terms of gas flow rate, applied voltage amplitude and frequency, and geometrical effects. The electric currents as a function of time for a given applied potential have been obtained, as well as the number densities and fluxes of plasma species.

Abstract: A hybrid model, called the hybrid plasma equipment model (HPEM), is used to describe the plasma chemistry in an inductively coupled plasma, operating in a gas mixture of C2H2 with either H2 or NH3, as typically used for carbon nanotube (CNT) growth. Two-dimensional profiles of power density, electron temperature and density, gas temperature, and densities of some plasma species are plotted and analyzed. Besides, the fluxes of the various plasma species towards the substrate (where the CNTs can be grown), as well as the decomposition rates of the feedstock gases (C2H2, NH3 and H2), are calculated as a function of the C2H2 fraction in both gas mixtures.

Abstract: Dielectric barrier discharges (DBDs) typically operate in the filamentary regime and thus exhibit great spatial and temporal non-uniformity. In order to optimize DBDs for various applications, such as in plasma catalysis, more fundamental insight is needed. Here, we consider how the millions of microdischarges, characteristic for a DBD, influence individual gas molecules. We use a Monte Carlo approach to determine the number of microdischarges to which a single molecule would be exposed, by means of particle tracing simulations through a full-scale packed bed DBD reactor, as well as an empty DBD reactor. We find that the fraction of microdischarges to which the molecules are exposed can be approximated as the microdischarge volume over the entire reactor gas volume. The use of this concept provides good agreement between a plasma-catalytic kinetics model and experiments for plasma-catalytic NH3 synthesis. We also show that the concept of the fraction of microdischarges indicates the efficiency by which the plasma power is transferred to the gas molecules. This generalised concept is also applicable for other spatially and temporally non-uniform plasmas.

Abstract: Packed bed dielectric barrier discharge (PB-DBD) plasma reactors are very promising for various plasma catalysis applications, but the exact mechanisms of plasma-catalyst interaction are far from understood, because the plasma discharge and catalyst/packing properties are mutually dependent. To better understand the effect of packing dielectric material on the electrical plasma properties, we study here a single bead DBD plasma reactor operating in dry air, with beads of different dielectric constant and for different applied voltages, by means of fluid modelling validated by optical imaging experiments. Our study reveals that the plasma in the single bead DBD reactor can manifest itself in two different modalities, i.e. (a) polar discharges at the bead poles in contact with the electrodes, and (b) a streamer discharge caused by surface ionization waves, which bridges the gas gap. Beads with high dielectric constant result in localised electric field enhancement and hence yield a reduction of the applied voltage required for plasma production. At low applied voltage, the discharge appears as polar discharges between the bead and the electrodes, and upon higher voltage it undergoes a transition into a bridging streamer discharge. The transition voltage to the streamer mode rises for beads with higher dielectric constant. These observations are important for plasma catalysis applications. A higher dielectric constant yields a higher electric field and thus higher average electron energy and density, giving rise to more reactive species, but it also yields a confined discharge near the contact points of packing beads, limiting the interaction area between the catalyst and the active plasma species. In addition, our model reveals that the dielectric bead behaves as a capacitor and traps charges, which can explain the significant occurrence of partial discharging in PB-DBDs and non-parallelogram shaped Lissajous plots. Hence, equivalent circuit modelling of PB-DBDs should take into account the role of packing beads in charge trapping as a capacitor.

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: In this paper, some modelling approaches to describe direct current (dc) magnetron discharges developed in our research groups will be presented, including an analytical model, Monte Carlo simulations for the electrons and for the sputtered atoms, a hybrid Monte Carlo-fluid model and particle-in-cell-Monte Carlo collision simulations. The strengths and limitations of the various modelling approaches will be explained, and some characteristic simulation results will be illustrated. Furthermore, some other simulation methods related to the magnetron device will be briefly explained, more specifically for calculating the magnetic field distribution inside the discharge, and for describing the (reactive) sputtering.

Abstract: Plasma medicine is a rapidly evolving multidisciplinary field at the intersection of chemistry, biochemistry, physics, biology, medicine and bioengineering. It holds great potential in medical, health care, dentistry, surgical, food treatment and other applications. This multidisciplinary nature and variety of possible applications come along with an inherent and intrinsic complexity. Advancing plasma medicine to the stage that it becomes an everyday tool in its respective fields requires a fundamental understanding of the basic processes, which is lacking so far. However, some major advances have already been made through detailed experiments over the last 15 years. Complementary, computer simulations may provide insight that is difficultif not impossibleto obtain through experiments. In this review, we aim to provide an overview of the various simulations that have been carried out in the context of plasma medicine so far, or that are relevant for plasma medicine. We focus our attention mostly on atomistic simulations dealing with plasmabiomolecule interactions. We also provide a perspective and tentative list of opportunities for future modelling studies that are likely to further advance the field.

Abstract: A numerical parameter study has been performed for a cylindrical atmospheric pressure dielectric barrier discharge (DBD) in helium with nitrogen impurities using a two-dimensional time-dependent fluid model. The calculated electric currents and gap voltages as a function of time for a given applied potential are presented, as well as the number densities of the various plasma species. This study shows that for the geometry under consideration the applied voltage parameters have a large impact on the electric current profiles and that the discharge current is always determined by the electron and ion conduction currents while the displacement current is nearly negligible. A relative broadening of the current profiles (compared with the duration of the half cycle of the applied voltage) with an increase in the applied frequency is obtained. Nearly sinusoidal current wave forms, usually typical for radio frequency DBDs, are observed while still operating at the frequencies of tens of kilohertz. For the setup under investigation, the Townsend mode of the DBD is observed in the entire range of applied voltage amplitudes and frequencies. It is shown that the average power density dissipated in the discharge increases with rising applied voltage and frequency. An increase in applied voltage frequency leads to an increase in the electron density and a decrease in electron energy, while increasing the voltage amplitude has the opposite effect.