Abstract: Hydrogen–graphite interactions are relevant to a wide variety of applications, ranging from astrophysics to fusion devices and nano-electronics. In order to shed light on these interactions, atomistic simulation using Molecular Dynamics (MD) has been shown to be an invaluable tool. It suffers, however, from severe timescale
limitations. In this work we apply the recently developed Collective Variable-Driven Hyperdynamics (CVHD) method to hydrogen etching of graphite for varying inter-impact times up to a realistic value of 1 ms, which corresponds to a flux of
1020 m2 s1. The results show that the erosion yield, hydrogen surface coverage and species distribution are significantly affected by the time between impacts. This can be explained by the higher probability of C–C bond breaking due to the prolonged exposure to thermal stress and the subsequent transition from ion- to thermal-induced etching. This latter regime of thermal-induced etching – chemical erosion – is here accessed for the first time using atomistic simulations. In conclusion, this study demonstrates that accounting for long time-scales significantly affects ion bombardment simulations and should not be neglected in a wide range of conditions, in contrast to what is typically assumed.

Abstract: CO2 conversion into value-added chemicals and fuels is considered as one of the great challenges of the 21st century. Due to the limitations of the traditional thermal approaches, several novel technologies are being developed. One promising approach in this field, which has received little attention to date, is plasma
technology. Its advantages include mild operating conditions, easy upscaling, and gas activation by energetic electrons instead of heat. This allows thermodynamically difficult reactions, such as CO2 splitting and the dry reformation of methane, to occur with reasonable energy cost. In this review, after exploring the traditional thermal approaches, we have provided a brief overview of the fierce competition between various novel approaches in a quest to find the most effective and efficient CO2 conversion technology. This is needed to critically assess whether plasma technology can be successful in an already crowded arena. The following questions need to be answered in this regard: are there key advantages to using plasma technology over other novel approaches, and if so, what is the flip side to the use of this technology? Can plasma technology be successful on its own, or can synergies be achieved by combining it with other technologies? To answer
these specific questions and to evaluate the potentials and limitations of plasma technology in general, this review presents the current state-of-the-art and a critical assessment of plasma-based CO2 conversion, as well as the future challenges for its practical implementation.

Abstract: We report on the thin film resistivity of several platinum-group metals (Ru, Pd, Ir, and Pt). Platinum-group thin films show comparable or lower resistivities than Cu for film thicknesses below about 5 nm due to a weaker thickness dependence of the resistivity. Based on experimentally determined mean linear distances between grain boundaries as well as ab initio calculations of the electron mean free path, the data for Ru, Ir, and Cu were modeled within the semiclassical Mayadas-Shatzkes model [Phys. Rev. B 1, 1382 (1970)] to assess the combined contributions of surface and grain boundary scattering to the resistivity. For Ru, the modeling results indicated that surface scattering was strongly dependent on the surrounding material with nearly specular scattering at interfaces with SiO2 or air but with diffuse scattering at interfaces with TaN. The dependence of the thin film resistivity on the mean free path is also discussed within the Mayadas-Shatzkes model in consideration of the experimental findings. Published by AIP Publishing.

Abstract: Ablation-controlled plasmas have been used in a range of technical applications where local thermodynamic equilibrium (LTE) is often violated near the wall due to the strong cooling effect caused by the ablation of wall materials. The thermodynamic and transport properties of ablated polytetrafluoroethylene (PTFE) vapor, which determine the flowing plasma behavior in such applications, are calculated based on a two-temperature model at atmospheric pressure. To our knowledge, no data for PTFE have been reported in the literature. The species composition and thermodynamic properties are numerically determined using the two-temperature Saha equation and the Guldberg-Waage equation according to van de Sanden et al's derivation. The transport coefficients, including viscosity, thermal conductivity and electrical conductivity, are calculated with the most recent collision interaction potentials using Devoto's electron and heavy-particle decoupling approach but expanded to the third-order approximation (second-order for viscosity) in the frame of the Chapman-Enskog method. Results are computed for different degrees of thermal non-equilibrium, i.e. the ratio of electron to heavy-particle temperatures, from 1 to 10, with electron temperature ranging from 300 to 40 000 K. Plasma transport properties in the LTE state obtained from the present work are compared with existing published results and the causes for the discrepancy analyzed. The two-temperature plasma properties calculated in the present work enable the modeling of wall ablation-controlled plasma processes.

Abstract: Porous materials are commonly used in microelectronics, as they can meet the demand for continuously shrinking electronic feature dimensions. However, they are facing severe challenges in plasma etching, due to plasma induced damage. In this paper, we present both the plasma characteristics and surface processing during the etching of porous materials. We explain how the damage occurs in the porous material during plasma etching for a wide range of chuck temperatures and the responsible mechanism for plasma damage-free etching at cryogenic temperature, by a combination of experiments and numerical modeling.

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: A packed bed dielectric barrier discharge (DBD) is widely used for plasma catalysis applications, but the exact plasma characteristics in between the packing beads are far from understood. Therefore, we study here these plasma characteristics by means of fluid modelling and experimental observations using ICCD imaging, for packing materials with different dielectric constants. Our study reveals that a packed bed DBD reactor in dry air at atmospheric pressure may show three types of discharges, i.e. positive restrikes, filamentary microdischarges, which can also be localized between two packing beads, and surface discharges (so-called surface ionization
waves). Restrikes between the dielectric surfaces result in the formation of filamentary microdischarges, while surface charging creates electric field components parallel to the dielectric surfaces, leading to the formation of surface discharges. A transition in discharge mode occurs from surface discharges to local filamentary discharges between the packing beads when the dielectric constant of the packing rises from 5 to 1000. This may have implications for the efficiency of plasma catalytic gas treatment, because the catalyst activation may be limited by constraining the discharge to the contact points of the beads. The production of reactive species occurs most in the positive restrikes, the surface discharges and the local microdischarges in between the beads, and is less significant in the longer filamentary microdischarges. The faster streamer propagation and discharge development with higher dielectric constant of the packing beads leads to a faster production of reactive species. This study is of great interest for plasma catalysis, where packing beads with different dielectric constants are often used as supports for the catalytic materials. It allows us to better understand how different packing materials can influence the performance of packed bed plasma reactors for environmental applications.

Abstract: This work explains the need for plasma models, introduces arguments for choosing the type of model that better fits the purpose of each study, and presents the basics of the most common nonequilibrium low-temperature plasma models and the information available from each one, along with an extensive list of references for complementary in-depth reading. The paper presents the following models, organised according to the level of multi-dimensional description of the plasma: kinetic models, based on either a statistical particle-in-cell/Monte-Carlo approach or the solution to the Boltzmann equation (in the latter case, special focus is given to the description of the electron kinetics); multi-fluid models, based on the solution to the hydrodynamic equations; global (spatially-average) models, based on the solution to the particle and energy rate-balance equations for the main plasma species, usually including a very complete reaction chemistry; mesoscopic models for plasma–surface interaction, adopting either a deterministic approach or a stochastic dynamical Monte-Carlo approach. For each plasma model, the paper puts forward the physics context, introduces the fundamental equations, presents advantages and limitations, also from a numerical perspective, and illustrates its application with some examples. Whenever pertinent, the interconnection between models is also discussed, in view of multi-scale hybrid approaches.

Abstract: Understanding the nature and effect of the multitude of plasma–surface interactions in plasma catalysis is a crucial requirement for further process development and improvement. A particularly intriguing and rather unique property of a plasma-catalytic setup is the ability of the plasma to modify the electronic structure, and hence chemical properties, of the catalyst through charging, i.e. the absorption of excess electrons. In this work, we develop a quantum chemical model based on density functional theory to study excess negative surface charges in a heterogeneous catalyst exposed to a plasma. This method is specifically applied to investigate plasma-catalytic CO2 activation on supported M/Al2O3 (M=Ti, Ni, Cu) single atom catalysts. We find that (1) the presence of a negative surface charge dramatically improves the reductive power of the catalyst, strongly promoting the splitting of CO2 to CO and oxygen, and (2) the relative activity of the investigated transition metals is also changed upon charging, suggesting that controlled surface charging is a powerful additional parameter to tune catalyst activity and selectivity. These results strongly point to plasma-induced surface charging of the catalyst as an important factor contributing to the plasma-catalyst synergistic effects frequently reported for plasma catalysis.

Abstract: <script type='text/javascript'>document.write(unpmarked('The crystal and defect structure of SnS crystals grown using chemical vapour deposition for application in electronic devices are investigated. The structural analysis shows the presence of two distinct crystal morphologies, that is thin flakes with lateral sizes up to 50 m and nanometer scale thickness, and much thicker but smaller crystallites. Both show similar Raman response associated with SnS. The structural analysis with transmission electron microscopy shows that the flakes are single crystals of -SnS with [010] normal to the substrate. Parallel with the surface of the flakes, lamellae with varying thickness of a new SnS phase are observed. High-resolution transmission electron microscopy (TEM), scanning transmission electron microscopy (STEM), first-principles simulations (DFT) and nanobeam diffraction (NBD) techniques are employed to characterise this phase in detail. DFT results suggest that the phase is a strain stabilised \u0027 one grown epitaxially on the -SnS crystals. TEM analysis shows that the crystallites are also -SnS with generally the [010] direction orthogonal to the substrate. Contrary to the flakes the crystallites consist of two to four grains which are tilted up to 15 degrees relative to the substrate. The various grain boundary structures and twin relations are discussed. Under high-dose electron irradiation, the SnS structure is reduced and -Sn formed. It is shown that this damage only occurs for SnS in direct contact with SiO2. Lay description SnS is a p-type semiconductor, which has attracted significant interest for electronic devices due to its unique properties, low-toxicity and abundance of Sn in nature. Although in the past it has been most extensively studied as the absorber material in solar cells, it has recently garnered interest for application as a p-type two-dimensional semiconductor in nanoelectronic devices due to its anisotropic layered structure similar to the better known phosphorene. Tin sulphide can take the form of several phases and the electronic properties of the material depend strongly on its crystal structure. It is therefore crucial to study the crystal structure of the material in order to predict the electronic properties and gain insight into the growth mechanism. In this work, SnS crystals deposited using a chemical vapour deposition technique are investigated extensively for their crystal and defect structure using transmission electron microscopy (TEM) and related techniques. We find the presence of two distinct crystal morphologies, that is thin flakes with lateral sizes up to 50 m and nm scale thickness, and much thicker but smaller crystallites. The flakes are single crystals of -SnS and contain lamellae with varying thickness of a different phase which appear to be -SnS at first glance. High-resolution scanning transmission electron microscopy is used to characterise these lamellae where the annular bright field (ABF) mode better reveals the position of the sulphur columns. The sulphur columns in the lamellae are found to be shifted relative to the -SnS structure which indicates the formation of a new phase which is a distorted version of the phase which we tentatively refer to as \u0027-SnS. Simulations based on density functional theory (DFT) are used to model the interface and a similar shift of sulphur columns in the -SnS layer is observed which takes place as a result of strong interaction at the interface between the two phases resulting in strain transfer. Nanobeam electron diffraction (NBD) is used to map the lattice mismatch in the thickness of the flakes which reveals good in-plane matching and some expansion out-of-plane in the lamellae. Contrary to the flakes the crystallites are made solely of -SnS and consist of two to four grains which are tilted up to 15 degrees relative to the substrate. The various grain boundary structures and twin relations are discussed. At high electron doses, SnS is reduced to -Sn, however the damage occurs only for SnS in direct contact with SiO2.'));

Abstract: It is generally known that antidiabetic activity is associated with an increased level of glucose uptake in adipocytes and skeletal muscle cells. However, the role of exogenous reactive oxygen and nitrogen species (RONS) in muscle development and more importantly in glucose uptake is largely unknown. We investigate the effect of RONS generated by cold atmospheric plasma (CAP) in glucose uptake. We show that the glucose uptake is significantly enhanced in differentiated L6 skeletal muscle cells after CAP treatment. We also observe a significant increase of the intracellular Ca++ and ROS level, without causing toxicity. One of the possible reasons for an elevated level of glucose uptake as well as intracellular ROS and Ca++ ions is probably the increased oxidative stress leading to glucose transport.

Abstract: Plasma technology is gaining increasing interest for CO2 conversion, but to maximize the energy efficiency, it is important to track the different energy transfers taking place in the plasma. In this paper, we study these mechanisms by a 0D chemical kinetics model, including the vibrational kinetics, for different conditions of reduced electric field, gas temperature and ionization degree, at a pressure of 100 mbar. Our model predicts a maximum conversion and energy efficiency of 32% and 47%, respectively, at conditions that are particularly beneficial for energy efficient CO2 conversion, i.e. a low reduced electric field (10 Td) and a low gas temperature (300 K). We study the effect of the efficiency by which the vibrational energy is used to dissociate CO2, as well as of the activation energy of the reaction CO2+O→CO+O2, to elucidate the theoretical limitations to the energy
efficiency. Our model reveals that these parameters are mainly responsible for the limitations in the energy efficiency. By varying these parameters, we can reach a maximum conversion and energy efficiency of 86%. Finally, we derive an empirical formula to estimate the maximum possible energy efficiency that can be reached under the assumptions of the model.

Abstract: The deposition of ultra-thin tungsten films and their related 2D chalcogen compounds on large area dielectric substrates by gas phase reactions is challenging. The lack of nucleation sites complicates the adsorption of W-related precursors and subsequent sulfurization usually requires high temperatures. We propose here a technique in which a thin solid amorphous silicon film is used as reductant for the gas phase precursor WF6 leading to the conversion to metallic W. The selectivity of the W conversion towards the underlying dielectric surfaces is demonstrated. The role of the Si surface preparation, the conversion temperature, and Si thickness on the formation process is investigated. Further, the in situ conversion of the metallic tungsten into thin stoichiometric WS2 is achieved by a cyclic approach based on WF6 and H2S pulses at the moderate temperature of 450 1C, which is much lower than usual oxide sulfurization processes.

Abstract: Epitaxially grown ultra-thin Si layers are often used to passivate Ge surfaces in the high-k gate module of (strained) Ge FinFET and Gate All Around devices. We use Si4H10 as Si precursor as it enables epitaxial Si growth at temperatures down to 330 degrees. C-V characteristics of blanket capacitors made on Ge virtual substrates point to the presence of an optimal Si thickness. In case of compressively strained Ge fin structures, the Si growth results in non-uniform and high strain levels in the strained Ge fin. These strain levels have been calculated for different shapes of the Ge fin and in function of the grown Si thickness. The high strain is the driving force for potential (unwanted) Ge surface reflow during Si deposition. The Ge surface reflow is strongly affected by the strength of the H-passivation during Si-capping and can be avoided by carefully selected process conditions. (C) The Author(s) 2018. Published by ECS.

Abstract: Although plasma catalysis is gaining increasing interest for various environmental applications, the underlying mechanisms are still far from understood. For instance, it is not yet clear whether and how plasma streamers can propagate in catalyst pores, and what is the minimum pore size to make this happen. As this is crucial information to ensure good plasma-catalyst interaction, we study here the mechanism of plasma streamer propagation in a catalyst pore, by means of a twodimensional particle-in-cell/Monte Carlo collision model, for various pore diameters in the nm range to μm-range. The so-called Debye length is an important criterion for plasma penetration into catalyst pores, i.e. a plasma streamer can penetrate into pores when their diameter is larger than the Debye length. The Debye length is typically in the order of a few 100 nm up to 1 μm at the conditions under study, depending on electron density and temperature in the plasma streamer. For pores in the range of ∼50 nm, plasma can thus only penetrate to some extent and at
very short times, i.e. at the beginning of a micro-discharge, before the actual plasma streamer reaches the catalyst surface and a sheath is formed in front of the surface. We can make plasma streamers penetrate into smaller pores (down to ca. 500 nm at the conditions under study) by increasing the applied voltage, which yields a higher plasma density, and thus reduces the Debye length. Our simulations also reveal that the plasma streamers induce surface charging of the catalyst pore sidewalls, causing discharge enhancement inside the pore, depending on pore diameter and depth.