Abstract: In this study, we investigated the diffusion of H-atoms to the subsurface and their further diffusion into the bulk of a Ni(111) crystal by means of density functional theory calculations in the context of thermal and plasma-assisted catalysis. The H-atoms at the surface can originate from the dissociative adsorption of H2 or CH4 molecules, determining the surface H-coverage. When a threshold H-coverage is passed, corresponding to 1.00 ML for the crystalline Ni(111) surface, the surface-bound H-atoms start to diffuse to the subsurface. A similar threshold coverage is observed for the interstitial H-coverage. Once the interstitial sites are filled up with a coverage above 1.00 ML of H, dissolution of interstitial H-atoms to the layer below the interstitial sites will be initiated. Hence, by applying a high pressure or inducing a reactive plasma and high temperature, increasing the H-flux to the surface, a large amount of hydrogen can diffuse in a crystalline metal like Ni and can be absorbed. The formation of metal hydride may modify the entire reaction kinetics of the system. Equivalently, the H-atoms in the bulk can easily go back to the surface and release a large amount of heat. In a plasma process, H-atoms are formed in the plasma, and therefore the energy barrier for dissociative adsorption is dismissed, thus allowing achievement of the threshold coverage without applying a high pressure as in a thermal process. As a result, depending on the crystal plane and type of metal, a large number of H-atoms can be dissolved (absorbed) in the metal catalyst, explaining the high efficiency of plasma-assisted catalytic reactions. Here, the mechanism of H-dissolution is established as a chemical identifier for the investigation of the reaction kinetics of a chemical process.

Abstract: The catalytic oxidative dehydrogenation of hydrocarbons is an industrially important process, in which selectivity is a key issue. We here investigate the conversion of methanol to formaldehyde on a vanadia surface employing long timescale simulations, reaching a time scale of seconds. In particular, we compare the thermal process to the case where an additional external electric field is applied, as would be the case in a direct plasma-catalysis setup. We find that the electric field influences the retention time of the molecules at the catalyst surface. These simulations provide an atomic scale insight in the thermal catalytic oxidative dehydrogenation process, and in how an external electric field may affect this process.

Abstract: In this paper, we show that the apparent delocalization of the conduction band reported from first-principles simulations for the high-mobility amorphous oxide semiconductor InGaZnO4 (a-IGZO) is an artifact induced by the periodic conditions imposed to the model. Given a sufficiently large unit-cell dimension (over 40 angstrom), the conduction band becomes localized. Such a model size is up to four times the size of commonly used models for the study of a-IGZO. This finding challenges the analyses done so far on the nature of the defects and on the interpretation of numerous electrical measurements. In particular, we re-interpret the meaning of the computed effective mass reported so far in literature. Our finding also applies to materials such as SiZnSnO, ZnSnO, InZnSnO, In2O3 or InAlZnO4 whose models have been reported to display a fully delocalized conduction band in the amorphous phase.

Abstract: The amorphous oxide semiconductor Indium-Gallium-Zinc-Oxide (a-IGZO) has gained a large technological relevance as a semiconductor for thin-film transistors in active-matrix displays. Yet, major questions remain unanswered regarding the atomic origin of threshold voltage control, doping level, hysteresis, negative bias stress (NBS), and negative bias illumination stress (NBIS). We undertake a systematic study of the effects of oxygen vacancies on the properties of a-IGZO by relating experimental observations to microscopic insights gained from first-principle simulations. It is found that the amorphous nature of the semiconductor allows unusually large atomic relaxations. In some cases, oxygen vacancies are found to behave as perfect shallow donors without the formation of structural defects. Once structural defects are formed, their transition states can vary upon charge and discharge cycles. We associate this phenomenon to a possible presence of hysteresis in the transfer curve of the devices. Under NBS, the creation of oxygen vacancies becomes energetically very stable, hence thermodynamically very likely. This generation process is correlated with the occurrence of the negative bias stress instabilities observed in a-IGZO transistors. While oxygen vacancies can therefore be related to NBS and hysteresis, it appears unlikely from our results that they are direct causes of NBIS, contrary to common belief.

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.'));