Abstract: Atmospheric pressure gliding arc (GA) discharges are gaining increasing interest for CO2 conversion and other gas conversion applications, due to their simplicity and high energy efficiency. However, they are characterized by some drawbacks, such as non-uniform gas treatment, limiting the conversion, as well as the development of a hot cathode spot, resulting in severe electrode degradation. In this work, we built a dual-vortex plasmatron, which is a GA plasma reactor with innovative electrode configuration, to solve the above problems. The design aims to improve the CO2 conversion capability of the GA reactor by elongating the arc in two directions, to increase the residence time of the gas inside the arc, and to actively cool the cathode spot by rotation of the arc and gas convection. The measured CO2 conversion and corresponding energy efficiency indeed look very promising. In addition, we developed a fluid dynamics non-thermal plasma model with argon chemistry, to study the arc behavior in the reactor and to explain the experimental results.

Abstract: Plasma catalysis, a promising technology for conversion of CO2 into value-added chemicals near room temperature, is gaining increasing interest. A dielectric barrier discharge (DBD) plasma has attracted attention due to its simple design and operation at near ambient conditions, ease to implement catalysts in the plasma zone and upscaling ability to industrial applications. To improve its main drawbacks, being relatively low conversion and energy efficiency, a packing material is used in the plasma discharge zone of the reactor, sometimes decorated by a catalytic material. Nevertheless, the extent to which different properties of the packing material influence plasma performance is still largely unexplored and unknown. In this study, the particular effect of synthesis induced differences in the morphology of a TiO2 shell covering a SiO2 core packing material on the plasma conversion of CO2 is studied. TiO2 has been successfully deposited around 1.6–1.8 mm sized SiO2 spheres by means of spray coating, starting from aqueous citratoperoxotitanate(IV) precursors. Parameters such as concentration of the Ti(IV) precursor solutions and addition of a binder were found to affect the shells’ properties and surface morphology and to have a major impact on the CO2 conversion in a packed bed DBD plasma reactor. Core-shell SiO2@TiO2 obtained from 0.25 M citratoperoxotitante(IV) precursors with the addition of a LUDOX binder showed the highest CO2 conversion 37.7% (at a space time of 70 s corresponding to an energy efficiency of 2%) and the highest energy efficiency of 4.8% (at a space time of 2.5 s corresponding to a conversion of 3%).

Abstract: We present a confined atmospheric pressure glow discharge plasma reactor, with very good performance towards dry reforming of methane, i.e., CO2 and CH4 conversion of 64 % and 94 %, respectively, at an energy cost of 3.5–4 eV/molecule (or 14–16 kJ/L). This excellent performance is among the best reported up to now for all types of plasma reactors in literature, and is due to the confinement of the plasma, which maximizes the fraction of gas passing through the active plasma region. The main product formed is syngas, with H2O and C2H2 as byproducts. We developed a quasi-1D chemical kinetics model, showing good agreement with the experimental results, which provides a thorough insight in the reaction pathways underlying the conversion of CO2 and CH4 and the formation of the different products.

Abstract: In this work, we take a crucial step towards the industrial readiness of plasma-based CO2 conversion. We present a stepwise method to study plasma reactors in series as a first approach to a recycle flow. By means of this procedure, the CO2 conversion is enhanced by a factor of 3, demonstrating that a single-pass plasma treatment performs far below the optimal capacity of the reactor. Furthermore, we explore the effect of O2 in the mixture with our flexible procedure. Addition of O2 in the mixture has a clear detrimental effect on the conversion, in agreement with other experiments in atmospheric pressure plasmas. O2 removal is however highly beneficial, demonstrating a conversion per pass that is 1.6 times higher than the standard procedure. Indeed, extracting one of the products prevents recombination reactions. Based on these insights, we discuss opportunities for further improvements, especially in the field of specialised separation techniques.

Abstract: A falling film dielectric barrier discharge (DBD) plasma reactor combined with adsorption on activated carbon textile material was optimized to minimize the formation of hazardous oxidation by-products from the treatment of persistent pesticides (alachlor, diuron and isoproturon) in water. The formation of by-products and the reaction mechanism was investigated by HPLC-TOF-MS. The maximum concentration of each by-product was at least two orders of magnitude below the initial pesticide concentration, during the first 10 min of treatment. After 30 min of treatment, the individual by-product concentrations had decreased to values of at least three orders of magnitude below the initial pesticide concentration. The proposed oxidation pathways revealed five main oxidation steps: dechlorination, dealkylation, hydroxylation, addition of a double-bonded oxygen and nitrification. The latter is one of the main oxidation mechanisms of diuron and isoproturon for air plasma treatment. To our knowledge, this is the first time that the formation of nitrificated intermediates is reported for the plasma treatment of non-phenolic compounds.

Abstract: Nitrogen oxyanions and oxyacids are important agents in atmospheric chemistry and medical biology. Although their chemical behavior in solution is relatively well understood, they may behave very differently at the water/air interface of atmospheric aerosols or at the membrane/water interface of cells. Here, we developed a fully classical model for molecular dynamics simulations of NO3−, NO2−, HNO3, and HNO2 in the framework of the GROMOS 53A6 and 54A7 force field versions. The model successfully accounted for the poorly structured solvation shell and ion pairing tendency of NO3−. Accurate pure-liquid properties and hydration free energies were obtained for the oxyacids. Simulations at the water/air interface showed a local enrichment of HNO3 and depletion of NO3−. The effect was discussed in light of earlier spectroscopic data and ab initio calculations, suggesting that HNO3 behaves as a weaker acid at the surface of water. Our model will hopefully allow for efficient and accurate simulations of nitrogen oxyanions and oxyacids in solution and at microheterogeneous interface environments.

Abstract: Cyrogenic etching of silicon is envisaged to enable better control over plasma processing in the microelectronics industry, albeit little is known about the fundamental differences compared to the room temperature process. We here present molecular dynamics simulations carried out to obtain sticking probabilities, thermal desorption rates, surface diffusion speeds, and sputter yields of F, F2, Si, SiF, SiF2, SiF3, SiF4, and the corresponding ions on Si(100) and on SiF13 surfaces, both at cryogenic and near room temperature. The different surface behavior during conventional etching and cryoetching is discussed. F2 is found to be relatively reactive compared to other species like SiF03. Thermal desorption occurs at a significantly lower rate under cryogenic conditions, which results in an accumulation of physisorbed species. Moreover, ion incorporation is often observed for ions with energies of 30400 eV, which results in a relatively low net sputter yield. The obtained results suggest that the actual etching of Si, under both cryogenic and near room temperature conditions, is based on the complete conversion of the Si surface to physisorbed SiF4, followed by subsequent sputtering of these molecules, instead of direct sputtering of the SiF03 surface.

Abstract: Continued miniaturization of metal-oxide-semiconductor field-effect transistors (MOSFETs) requires an ever-decreasing thickness of the gate oxide. The structure of ultrathin silicon oxide films, however, critically depends on the oxidation mechanism. Using reactive atomistic simulations, we here demonstrate how the oxidation mechanism in hyperthermal oxidation of such structures may be controlled by the oxidation temperature and the oxidant energy. Specifically, we study the interaction of hyperthermal oxygen with energies of 15 eV with thin SiOx (x ≤ 2) films with a native oxide thickness of about 10 Å. We analyze the oxygen penetration depth probability and compare with results of the hyperthermal oxidation of a bare Si(100){2 × 1} (c-Si) surface. The temperature-dependent oxidation mechanisms are discussed in detail. Our results demonstrate that, at low (i.e., room) temperature, the penetrated oxygen mostly resides in the oxide region rather than at the SiOx|c-Si interface. However, at higher temperatures, starting at around 700 K, oxygen atoms are found to penetrate and to diffuse through the oxide layer followed by reaction at the c-Si boundary. We demonstrate that hyperthermal oxidation resembles thermal oxidation, which can be described by the DealGrove model at high temperatures. Furthermore, defect creation mechanisms that occur during the oxidation process are also analyzed. This study is useful for the fabrication of ultrathin silicon oxide gate oxides for metal-oxide-semiconductor devices as it links parameters that can be straightforwardly controlled in experiment (oxygen temperature, velocity) with the silicon oxide structure.

Abstract: Molecular dynamics simulations were used to investigate the size-dependent melting mechanism of nickel nanoclusters of various sizes. The melting process was monitored by the caloric curve, the overall cluster Lindemann index, and the atomic Lindemann index. Size-dependent melting temperatures were determined, and the correct linear dependence on inverse diameter was recovered. We found that the melting mechanism gradually changes from dynamic coexistence melting to surface melting with increasing cluster size. These findings are of importance in better understanding carbon nanotube growth by catalytic chemical vapor deposition as the phase state of the catalyst nanoparticle codetermines the growth mechanism.

Abstract: We present a computational study for the conversion of CH4 and CO2 into value-added chemicals, i.e., the so-called dry reforming of methane, in a dielectric barrier discharge reactor. A zero-dimensional chemical kinetics model is applied to study the plasma chemistry in a 1:1 CH4/CO2 mixture. The calculations are first performed for one microdischarge pulse and its afterglow, to study in detail the chemical pathways of the conversion. Subsequently, long time-scale simulations are carried out, corresponding to real residence times in the plasma, assuming a large number of consecutive microdischarge pulses, to mimic the conditions of the filamentary discharge regime in a dielectric barrier discharge (DBD) reactor. The conversion of CH4 and CO2 as well as the selectivity of the formed products and the energy cost and energy efficiency of the process are calculated and compared to experiments for a range of different powers and gas flows, and reasonable agreement is reached.

Abstract: The adsorption probability and reaction behavior of CHx plasma species on various nickel catalyst surfaces is investigated by means of reactive molecular dynamics (MD) simulations using the ReaxFF potential. Such catalysts are used in the reforming of hydrocarbons and in the growth of carbon nanotubes, and further insight in the underlying mechanisms of these processes is needed to increase their applicability. Single and consecutive impacts of CHx radicals (x={1,2,3}) were performed on four different Ni surfaces, at a temperature of 400 K. The adsorption probability is shown to be related to the number of free electrons, i.e. a higher number leads to more adsorptions, and the steric hindrance caused by the hydrogen atoms bonded to the impacting CHx species. Furthermore, some of the CH bonds break after adsorption, which generally leads to diffusion of the hydrogen atom over the surface. Additionally, these adsorbed H-atoms can be used in reactions to form new molecules, such as CH4 and C2Hx, although this is dependent on the precise morphology of the surface. New molecules are also formed by subtraction of H-atoms from adsorbed radicals, leading to occasional formation of H2 and C2Hx molecules.

Abstract: Device performance of Si|SiOx core-shell based nanowires critically depends on the exact control over the oxide thickness. Low-temperature plasma oxidation is a highly promising alternative to thermal oxidation allowing for improved control over the oxidation process, in particular for ultrasmall Si nanowires. We here elucidate the room temperature plasma oxidation mechanisms of ultrasmall Si nanowires using hybrid molecular dynamics / force-bias Monte Carlo simulations. We demonstrate how the oxidation and concurrent water formation mechanisms are a function of the oxidizing plasma species and we demonstrate how the resulting core-shell oxide thickness can be controlled through these species. A new mechanism of water formation is discussed in detail. The results provide a detailed atomic level explanation of the oxidation process of highly curved Si surfaces. These results point out a route toward plasma-based formation of ultrathin core-shell Si|SiOx nanowires at room temperature.

Abstract: Plasma catalysis is gaining increasing interest for various environmental applications, but the crucial question is whether plasma can be created inside catalyst pores and under which conditions. In practice, various catalytic support materials are used, with various dielectric constants. We investigate here the influence of the dielectric constant on the plasma properties inside catalyst pores and in the sheath in front of the pores, for various pore sizes. The calculations are performed by a two-dimensional fluid model for an atmospheric pressure dielectric barrier discharge in helium. The electron impact ionization rate, electron temperature, electron and ion density, as well as the potential distribution and surface charge density, are analyzed for a better understanding of the discharge behavior inside catalyst pores. The results indicate that, in a 100 μm pore, the electron impact ionization in the pore, which is characteristic for the plasma generation inside the pore, is greatly enhanced for dielectric constants below 300. Smaller pore sizes only yield enhanced ionization for smaller dielectric constants, i.e., up to εr = 200, 150, and 50 for pore sizes of 50, 30, and 10 μm. Thus, the most common catalyst supports, i.e., Al2O3 and SiO2, which have dielectric constants around εr = 8−11 and 4.2, respectively, should allow more easily that microdischarges can be formed inside catalyst pores, even for smaller pore sizes. On the other hand, ferroelectric materials with dielectric constants above 300 never seem to yield plasma enhancement inside catalyst pores, not even for 100 μm pore sizes. Furthermore, it is clear that the dielectric constant of the material has a large effect on the extent of plasma enhancement inside the catalyst pores, especially in the range between εr = 4 and εr = 200. The obtained results are explained in detail based on the surface charge density at the pore walls,
and the potential distribution and electron temperature inside and above the pores. The results obtained with this model are
important for plasma catalysis, as the production plasma species in catalyst pores might affect the catalyst properties, and thus
improve the applications of plasma catalysis.

Abstract: Molecular dynamics (MD) simulations were performed to provide atomic scale insight in the initial interaction between hydroxyl radicals (OH) and peptide systems in solution. These OH radicals are representative reactive oxygen species produced by cold atmospheric plasmas. The use of plasma for biomedical applications is gaining increasing interest, but the fundamental mechanisms behind the plasma modifications still remain largely elusive. This study helps to gain more insight in the underlying mechanisms of plasma medicine but is also more generally applicable to peptide oxidation, of interest for other applications. Combining both reactive and nonreactive MD simulations, we are able to elucidate the reactivity of the amino acids inside the peptide systems and their effect on their structure up to 1 μs. Additionally, experiments were performed, treating the simulated peptides with a plasma jet. The computational results presented here correlate well with the obtained experimental data and highlight the importance of the chemical environment for the reactivity of the individual amino acids, so that specific amino acids are attacked in higher numbers than expected. Furthermore, the long time scale simulations suggest that a single oxidation has an effect on the 3D conformation due to an increase in hydrophilicity and intra- and intermolecular interactions.

Abstract: The gliding arc plasmatron (GAP) is a highly efficient atmospheric plasma source, which is very promising for CO2 conversion applications. To understand its operation principles and to improve its application, we present here comprehensive modeling results, obtained by means of computational fluid dynamics simulations and plasma modeling. Because of the complexity of the CO2 plasma, a full 3D plasma model would be computationally impractical. Therefore, we combine a 3D turbulent gas flow model with a 2D plasma and gas heating model in order to calculate the plasma parameters and CO2 conversion characteristics. In addition, a complete 3D gas flow and plasma model with simplified argon chemistry is used to evaluate the gliding arc evolution in space and time. The calculated values are compared with experimental data from literature as much as possible in order to validate the model. The insights obtained in this study are very helpful for improving the application of CO2 conversion, as they allow us to identify the limiting factors in the performance, based on which solutions can be provided on how to further improve the capabilities of CO2 conversion in the GAP.