Abstract: tWe investigated the plasma-assisted catalytic reactions for the production of value-added chemicalsfrom Ni-catalyzed plasma dry reforming of methane by means of density functional theory (DFT). Weinspected many activation barriers, from the early stage of adsorption of the major chemical fragmentsderived fromCH4andCO2molecules up to the formation of value-added chemicals at the surface, focusingon the formation of methanol, as well as the hydrogenation of C1and C2hydrocarbon fragments. Theactivation barrier calculations show that the presence of surface-bound H atoms and in some cases alsoremaining chemical fragments at the surface facilitates the formation of products. This implies that thehydrogenation of a chemical fragment on the hydrogenated crystalline surface is energetically favouredcompared to the simple hydrogenation of the chemical fragment at the bare Ni(111) surface. Indeed, thepresence of hydrogen modifies the electronic structure of the surface and the course of the reactions.We therefore conclude that surface-bound H atoms, and to some extent also the remaining chemicalfragments at the crystalline surface, induce the following effects: they facilitate associative desorption ofmethanol and ethane by increasing the rate of H-transfer to the adsorbed fragments while they impedehydrogenation of ethylene to ethane, thus promoting again the desorption of ethylene. Overall, they thusfacilitate the catalytic conversion of the formed fragments from CH4and CO2, into value-added chemicals.Finally, we believe that the retention of methane fragments, especially CH3, in the presence of surface-boundHatoms (as observed here for Ni) can be regarded as an identifier for the proper choice of a catalystfor the production of value-added chemicals.

Abstract: A hybrid model, called the hybrid plasma equipment model (HPEM), was used to study an inductively coupled plasma in gas mixtures of H2 or NH3 with CH4 or C2H2 used for the synthesis of carbon nanotubes or carbon nanofibres (CNTs/CNFs). The plasma properties are discussed for different gas mixtures at low and moderate pressures, and the growth precursors for CNTs/CNFs are analysed. It is found that C2H2, C2H4 and C2H6 are the predominant molecules in CH4 containing plasmas besides the feedstock gas, and serve as carbon sources for CNT/CNF formation. On the other hand, long-chain hydrocarbons are observed in C2H2-containing plasmas. Furthermore, the background gases CH4 and C2H2 show a different decomposition rate with H2 or NH3 addition at moderate pressures.

Abstract: Using reactive molecular dynamics simulations, the melting behavior of nickel-carbon nanoclusters is examined. The phase diagrams of icosahedral and Wulff polyhedron clusters are determined using both the Lindemann index and the potential energy. Formulae are derived for calculating the equilibrium constants and the solid and liquid fractions during a phase transition, allowing more rational determination of the melting temperature with respect to the arbitrary Lindemann value. These results give more insight into the properties of nickel-carbon nanoclusters in general and can specifically be very useful for a better understanding of the synthesis of carbon nanotubes using the catalytic chemical vapor deposition method.

Abstract: Dry reforming of methane (DRM) in a gliding arc plasmatron is studied for different CH4 fractions in the mixture. The CO2 and CH4 conversions reach their highest values of approximately 18 and 10%, respectively, at 25% CH4 in the gas mixture, corresponding to an overall energy cost of 10 kJ L@1 (or 2.5 eV per molecule) and an energy efficiency of 66%. CO and H2 are the major products, with the formation of smaller fractions of C2Hx (x=2, 4, or 6) compounds and H2O. A chemical kinetics model is used to investigate the underlying chemical processes. The calculated CO2 and CH4 conversion and the energy efficiency are in good agreement with the experimental data. The model calculations reveal that the reaction of CO2 (mainly at vibrationally excited levels) with H radicals is mainly responsible for
the CO2 conversion, especially at higher CH4 fractions in the mixture, which explains why the CO2 conversion increases with increasing CH4 fraction. The main process responsible for CH4 conversion is the reaction with OH radicals. The excellent energy efficiency can be explained by the non-equilibrium character of the plasma, in which the electrons mainly activate the gas molecules, and by the important role of the vibrational kinetics of CO2. The results demonstrate that a gliding arc plasmatron is very promising for DRM.

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: The initial stage of nanoparticle formation and growth in radiofrequency acetylene (C2H2) plasmas is investigated by means of a self-consistent one-dimensional fluid model. A detailed chemical kinetic scheme, containing electron impact, ion-neutral, and neutral-neutral reactions, has been developed in order to predict the underlying dust growth mechanisms and the most important dust precursors. The model considers 41 different species (neutrals, radicals, ions, and electrons) describing hydrocarbons (CnHm) containing up to 12 carbon atoms. Possible routes for particle growth are discussed. Both positive and negative ion reaction pathways are considered, as consecutive anion- and cation-molecule reactions seem to lead to a fast build up of the carbon skeleton.

Abstract: The reproducible temporal separation of ion signals generated from a single multi-element droplet, observed in previous studies, was investigated in detail in this work using an ICPTOFMS with high temporal resolution. It was shown that the signal peak intensities of individual elements temporally shift relative to each other only for droplets moving through the plasma off-axis. The magnitude of these shifts correlated with the vaporization temperatures of the analytes and depended on the radial position of the droplets as well as on the thermal properties and velocity profiles of the carrier gases of the ICP. The occurrence of the signal shifting was explained by a spatial separation of analytes already present in the vapor phase in the ICP from a yet unvaporized residue of the droplet. This separation is most likely driven by anisotropic diffusion of vaporized analytes towards the plasma axis and a radial velocity gradient. The proposed explanation is supported by modeling of the gas velocities inside the ICP and imaging of the atomic and ionic emissions produced from single droplets, whose patterns were sloping towards the center of the torch. The effects observed in these studies are important not only for the fundamental understanding of analyteplasma interactions but have also a direct impact on the signal intensities and stability.

Abstract: The electron bounce resonance heating (BRH) in dual-frequency capacitively coupled plasmas operated in oxygen is studied by different experimental methods and a particle-in-cell/Monte Carlo collision (PIC/MCC) simulation, and compared with the electropositive argon discharge. In comparison with argon, the experimental results show that in an oxygen discharge the resonance peaks in positive-ion density and light intensity tend to occur at larger electrode gaps. Moreover, at electrode gaps L > 2.5 cm, the positive-ion (and electron) density and the light emission drop monotonically in the oxygen discharge upon increasing L, whereas they rise (after an initial drop) in the argon case. At resonance gap the electronegativity reaches its maximum due to the BRH. All these experimental observations are explained by PIC/MCC simulations, which show that in the oxygen discharge the bulk electric field becomes quite strong and is out of phase with the sheath field. Therefore, it retards the resonance electrons when traversing the bulk, resulting in a suppressed BRH. Both experiment and simulation results show that this effect becomes more pronounced at lower high-frequency power, when the discharge mode changes from electropositive to electronegative. In a pure oxygen discharge, the BRH is suppressed with increasing pressure and almost diminishes at 12 Pa. Finally, the driving frequency significantly affects the BRH, because it determines the phase relation between bulk electric field and sheath electric field.

Abstract: In this work particle-in-cell/Monte Carlo collision simulations are performed to study the heating mechanism and plasma characteristics in direct current (dc) superposed radio-frequency (RF) capacitively coupled plasmas, operated both in single-frequency (SF) and dual-frequency (DF) regimes. An RF (60/2 MHz) source is applied on the bottom electrode to sustain the discharge, and a dc source is fixed on the top electrode. The heating mechanism appears to be very different in dc superposed SF and DF discharges. When only a single source of 60 MHz is applied, the plasma bulk region is reduced by the dc source, thus the ionization rate and hence the electron density decrease with rising dc voltage. However, when a DF source of 60 and 2 MHz is applied, the electron density can increase upon addition of a dc voltage, depending on the gap length and applied dc voltage. This is explained from the spatiotemporal ionization rates in the DF discharge. In fact, a completely different behavior is observed for the ionization rate in the two half-periods of the LF source. In the first LF half-period, the situation resembles the dc superposed SF discharge, and the reduced plasma bulk region due to the negative dc bias results in a very small effective discharge area and a low ionization rate. On the other hand, in the second half-period, the negative dc bias is to some extent counteracted by the LF voltage, and the sheath close to the dc electrode becomes particularly thin. Consequently, the amplitude of the high-frequency sheath oscillations at the top electrode is largely enhanced, while the LF sheath at the bottom electrode is in its expanding phase and can thus well confine the high-energy electrons. Therefore, the ionization rate increases considerably in this second LF half-period. Furthermore, in addition to the comparison between SF and DF discharges and the effect of gap length and dc voltage, the effect of secondary electrons is examined.

Abstract: Carbon nanotubes (CNTs) are nowadays routinely grown in a thermal CVD setup. State-of-the-art plasma-enhanced CVD (PECVD) growth, however, offers advantages over thermal CVD. A lower growth temperature and the growth of aligned freestanding single-walled CNTs (SWNTs) makes the technique very attractive. The atomic scale growth mechanisms of PECVD CNT growth, however, remain currently entirely unexplored. In this contribution, we employed molecular dynamics simulations to focus on the effect of applying an electric field on the SWNT growth process, as one of the effects coming into play in PECVD. Using sufficiently strong fields results in (a) alignment of the growing SWNTs, (b) a better ordering of the carbon network, and (c) a higher growth rate relative to thermal growth rate. We suggest that these effects are due to the small charge transfer occurring in the Ni/C system. These simulations constitute the first study of PECVD growth of SWNTs on the atomic level.