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.

Abstract: Because of the unique properties of plasma technology, its use in gas conversion applications is gaining significant interest around the globe. Plasma-based CO2 and CH4 conversion has become a major research area. Many investigations have already been performed regarding the single-component gases, that is, CO2 splitting and CH4 reforming, as well as for two-component mixtures, that is, dry reforming of methane
(CO2/CH4), partial oxidation of methane (CH4/O2), artificial photosynthesis (CO2/H2O), CO2 hydrogenation (CO2/H2), and even first steps toward the influence of N2 impurities have been taken, that is, CO2/N2 and CH4/N2. In this Feature Article we briefly discuss the advances made in literature for these different steps from a plasma chemistry modeling point of view. Subsequently, we present a comprehensive plasma chemistry set, combining the knowledge gathered in this field so far and supported with extensive experimental data. This set can be used for chemical kinetics plasma modeling for all possible combinations of CO2, CH4, N2, O2, and H2O to investigate the bigger picture of the underlying plasmachemical pathways for these mixtures in a dielectric barrier discharge plasma. This is extremely valuable
for the optimization of existing plasma-based CO2 conversion and CH4 reforming processes as well as for investigating the influence of N2, O2, and H2O on these processes and even to support plasma-based multireforming processes.

Abstract: Surface charging is an often overlooked factor in many plasma-surface interactions and in particular in plasma catalysis. In this study, we investigate the effect of excess electrons induced by a plasma on the adsorption properties of CO2 on titania-supported Cu-5 and Ni-5 clusters using spin-polarized and dispersion-corrected density functional theory calculations. The effect of excess electrons on the adsorption of Ni and Cu pentamers as well as on CO2 adsorption on a pristine anatase TiO2(101) slab is studied. Our results indicate that adding plasma-induced excess electrons to the system leads to further stabilization of the bent CO2 structure. Also, dissociation of CO2 on charged clusters is energetically more favorable than on neutral clusters. We hypothesize that surface charge is a plausible cause for the synergistic effects sometimes observed in plasma catalysis.

Abstract: In this contribution, we present an overview of the various techniques for combining atomistic molecular dynamics with Monte Carlo simulations, mainly in the context of condensed matter systems, as well as a brief summary of the main accelerated dynamics techniques. Special attention is given to the force bias Monte Carlo technique and its combination with molecular dynamics, in view of promising recent developments, including a definable timescale. Various examples of the application of combined molecular dynamics / Monte Carlo simulations are given, in order to demonstrate the enhanced simulation efficiency with respect to either pure molecular dynamics or Monte Carlo.

Abstract: A seriesof [TiN/TiAlN]nmultilayer coatingswith different bilayer numbers n=5, 10, 25, 50, and 100 were deposited on stainless steel substrate AISI 304 by a lateral rotating cathode arc technique in a flowing nitrogen atmosphere. The composition and microstructure of the coatings have been analyzed by using energy dispersive X-ray spectroscopy, X-ray diffraction (XRD), and conventional and high-resolution transmission electron microscopy (HRTEM). XRD analysis shows that the preferential orientation growth along the (111) direction is reduced in the multilayer coatings. TEM analysis reveals that the grain size of the coatings decreases with increasing bilayer number. HRTEMimaging of the multilayer coatings shows a high density misfit dislocation between the TiN and TiAlN layers. The cross-plane thermal conductivity of the coatings was measured by a pulsed photothermal reflectance technique. With increasing bilayer number, the multilayer coatings' thermal conductivity decreases gradually. This reduction of thermal conductivity can be ascribed to increased phonon scattering due to the disruption of columnar structure, reduced preferential orientation, decreased grain size of the coatings and present misfit dislocations at the interfaces.

Abstract: Pesticides are emerging contaminants frequently detected in the aquatic environment. In this work, a novel approach combining activated carbon adsorption, oxygen plasma treatment and ozonation was studied for the removal of the persistent chlorinated pesticide alachlor. A comparison was made between the removal efficiency and energy consumption for two different reactor operation modes: batchrecirculation and single-pass mode. The kinetics study revealed that the insufficient removal of alachlor by adsorption was significantly improved in terms of degradation efficiency and energy consumption when combined with the plasma treatment. The best efficiency (ca. 80% removal with an energy cost of 19.4 kWh mÀ3) was found for the single-pass operational mode of the reactor. In the batch-recirculating process, a complete elimination of alachlor by plasma treatment was observed after 30 min of treatment. Analysis of the reactive species induced by plasma in aqueous solutions showed that the decomposition of alachlor mainly occurred through a radical oxidation mechanism, with a minor contribution of long-living oxidants (O3, H2O2). Investigation of the alachlor oxidation pathways revealed six different oxidation mechanisms, including the loss of aromaticity which was never before reported for plasma-assisted degradation of aromatic pesticides. It was revealed that the removal rate and energy cost could be further improved with more than 50% by additional O3 gas bubbling in the solution reservoir.