Abstract: Since the discovery of graphene, a broad range of two-dimensional (2D) materials has captured the attention of the scientific communities. Materials, such as hexagonal boron nitride (hBN) and the transition metal dichalcogenides (TMDs) family, have shown promising semiconducting and insulating properties that are very appealing for the semiconductor industry. Recently, the possibility of taking advantage of the properties of 2D-based heterostructures has been investigated for low-power nanoelectronic applications. In this work, we aim at evaluating the relation between the nature of the materials used in such heterostructures and the amplitude of the layer-to-layer charge transfer induced by an external electric field, as is typically present in nanoelectronic gated devices. A broad range of combinations of TMDs, graphene, and hBN has been investigated using density functional theory. Our results show that the electric field induced charge transfer strongly depends on the nature of the 2D materials used in the van der Waals heterostructures and to a lesser extent on the relative orientation of the materials in the structure. Our findings contribute to the building of the fundamental understanding required to engineer electrostatically the doping of 2D materials and to establish the factors that drive the charge transfer mechanisms in electron tunneling-based devices. These are key ingredients for the development of 2D -based nanoelectronic devices.

Abstract: The continuous miniaturization of nanodevices, such as transistors, solar cells, and optical fibers, requires the controlled synthesis of (ultra)thin gate oxides (<10 nm), including Si gate-oxide (SiO2) with high quality at the atomic scale. Traditional thermal growth of SiO2 on planar Si surfaces, however, does not allow one to obtain such ultrathin oxide due to either the high oxygen diffusivity at high temperature or the very low sticking ability of incident oxygen at low temperature. Two recent techniques, both operative at low (room) temperature, have been put forward to overcome these obstacles: (i) hyperthermal oxidation of planar Si surfaces and (ii) thermal or plasma-assisted oxidation of nonplanar Si surfaces, including Si nanowires (SiNWs). These nanooxidation processes are, however, often difficult to study experimentally, due to the key intermediate processes taking place on the nanosecond time scale.
In this Account, these Si nano-oxidation techniques are discussed from a computational point of view and compared to both hyperthermal and thermal oxidation experiments, as well as to well-known models of thermal oxidation, including the Deal−Grove, Cabrera−Mott, and Kao models and several alternative mechanisms. In our studies, we use reactive molecular dynamics (MD) and hybrid MD/Monte Carlo simulation techniques, applying the Reax force field. The incident energy of oxygen species is chosen in the range of 1−5 eV in hyperthermal oxidation of planar Si surfaces in order to prevent energy-induced damage. It turns out that hyperthermal growth allows for two growth modes, where the ultrathin oxide thickness depends on either (1) only the kinetic energy of the incident oxygen species at a growth temperature below Ttrans = 600 K, or (2) both the incident energy and the growth temperature at a growth temperature above Ttrans. These modes are specific to such ultrathin oxides, and are not observed in traditional thermal oxidation, nor theoretically considered by already existing models. In the case of thermal or plasma-assisted oxidation of small Si nanowires, on the other hand, the thickness of the ultrathin oxide is a function of the growth temperature and the nanowire diameter. Below Ttrans, which varies with the nanowire diameter, partially oxidized SiNW are formed, whereas complete oxidation to a SiO2 nanowire occurs only above Ttrans. In both nano-oxidation processes at lower temperature (T < Ttrans), final sandwich c-Si|SiOx|a-SiO2 structures are obtained due to a competition between overcoming the energy barrier to penetrate into Si subsurface layers and the compressive stress (∼2−3 GPa) at the Si crystal/oxide interface. The overall atomic-simulation results strongly indicate that the thickness of the intermediate SiOx (x < 2) region is very limited (∼0.5 nm) and constant irrespective of oxidation parameters. Thus, control over the ultrathin SiO2 thickness with good quality is indeed possible by accurately tuning the oxidant energy, oxidation temperature and surface curvature.
In general, we discuss and put in perspective these two oxidation mechanisms for obtaining controllable ultrathin gate-oxide films, offering a new route toward the fabrication of nanodevices via selective nano-oxidation.

Abstract: An extensive CO2 plasmamodel is presented that is relevant for the production of ‘‘solar fuels.’’ It is based on reaction rate coefficients fromrigorously reviewed literature, and is augmented with reactionrate coefficients that are obtained fromscaling laws.The input data set,which is suitable for usage with the plasma simulation software Plasimo (https://plasimo.phys.tue.nl/), is available via the Plasimo and publisher’s websites.1 The correctness of this model implementation has been established by independent ZDPlasKin implementation (http://www.zdplaskin.
laplace.univ-tlse.fr/), to verify that the results agree. Results of these ‘‘global models’’ are presented for a DBD plasma reactor.

Abstract: An Ar plasma sustained by a surfaguide wave launcher is investigated at intermediate pressure (200–2667 Pa). Two 2D self-consistent models (quasi-neutral and plasma bulk-sheath) are developed and benchmarked. The complete set of electromagnetic and fluid equations and the boundary conditions are presented. The transformation of fluid equations from a local reference frame, that is, moving with plasma or when the gas flow is zero, to a laboratory reference frame, that is,
accounting for the gas flow, is discussed. The pressure range is extended down to 80 Pa by experimental measurements. The electron temperature decreases with pressure. The electron density depends linearly on power, and changes its behavior with pressure depending on the product of pressure and radial plasma size.

Abstract: The modelling of a gliding arc discharge (GAD) is studied by means of the quasineutral (QN) plasma modelling approach. The model is first evaluated for reliability and proper description of a gliding arc discharge at atmospheric pressure, by comparing with a more elaborate non-quasineutral (NQN) plasma model in two different geometries – a 2D axisymmetric and a Cartesian geometry. The NQN model is considered as a reference, since it provides a continuous self-consistent plasma description, including the near electrode regions. In general, the results of the QN model agree very well with those obtained from the NQN model. The small differences between both models are attributed to the approximations in the derivation of the QN model. The use of the QN model provides a substantial reduction of the computation time compared to the NQN model, which is crucial for the development of more complex models in three dimensions or with complicated chemistries. The latter is illustrated for (i) a reverse vortex flow(RVF) GAD in argon, and (ii) a GAD in CO2. The RVF discharge is modelled in three dimensions and the effect of the turbulent heat transport on the plasma and gas characteristics is
discussed. The GAD model in CO2 is in a 1D geometry with axial symmetry and provides results for the time evolution of the electron, gas and vibrational temperature of CO2, as well as for the molar fractions of the different species.

Abstract: The conversion of atmospheric nitrogen into valuable compounds, that is, so-called nitrogen fixation, is gaining increased interest, owing to the essential role in the nitrogen cycle of the biosphere. Plasma technology, and more specifically gliding arc plasma, has great potential in this area, but little is known about the underlying mechanisms. Therefore, we developed a detailed chemical kinetics model for a pulsed-power gliding-arc reactor operating at atmospheric pressure for nitrogen oxide synthesis. Experiments are performed to validate the model and reasonable agreement is reached between the calculated and measured NO and NO2 yields and the corresponding energy efficiency for NOx formation for different N2/O2 ratios, indicating that the model can provide a realistic picture of the plasma chemistry. Therefore, we can use the model to investigate the reaction pathways for the formation and loss of NOx. The results indicate that vibrational excitation of N2 in the gliding arc contributes significantly to activating the N2 molecules, and leads to an energy efficient way of NOx production, compared to the thermal process. Based on the underlying chemistry, the model allows us to propose solutions on how to further improve the NOx formation by gliding arc technology. Although the energy efficiency of the gliding-arc-based nitrogen fixation process at the present stage is not comparable to the world-scale HaberBosch process, we believe our study helps us to come up with more realistic scenarios of entering a cutting-edge innovation in new business cases for the decentralised production of fertilisers for agriculture, in which low-temperature plasma technology might play an important role.

Abstract: In this review, we discuss the progress of emerging dry processes for nanoscale fabrication. Experts in the fields of plasma processing have contributed to addressing the increasingly challenging demands in achieving atomic-level control of material selectivity and physicochemical reactions involving ion bombardment. The discussion encompasses major challenges shared across the plasma science and technology community. Focus is placed on advances in the development of fabrication technologies for emerging materials, especially metallic and intermetallic compounds and multiferroic, and two-dimensional (2D) materials, as well as state-of-the-art techniques used in nanoscale semiconductor manufacturing with a brief summary of future challenges.

Abstract: Low-temperature plasmas are gaining a lot of interest for environmental and energy applications. A large research field in these applications is the conversion of CO2 into chemicals and fuels. Since CO2 is a very stable molecule, a key performance indicator for the research on plasma-based CO2 conversion is the energy efficiency. Until now, the energy efficiency in atmospheric plasma reactors is quite low, and therefore we employ here a novel type of plasma reactor, the gliding arc plasmatron (GAP). This paper provides a detailed experimental and computational study of the CO2 conversion, as well as the energy cost and efficiency in a GAP. A comparison with thermal conversion, other plasma types and other novel CO2 conversion technologies is made to find out whether this novel plasma reactor can provide a significant contribution to the much-needed efficient conversion of CO2. From these comparisons it becomes evident that our results are less than a factor of two away from being cost competitive and already outperform several other new technologies. Furthermore, we indicate how the performance of the GAP can still be improved by further exploiting its non-equilibrium character. Hence, it is clear that the GAP is very promising for CO2 conversion.

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: For the energy industry, a process that is able to transform methane—being the prime component of natural gas—efficiently into a liquid product would be equivalent to a goose with golden eggs. As such it is no surprise that research efforts in this field already date back to the nineteen hundreds. Plasma technology can be considered to be a novel player in this field, but nevertheless one with great potential. Over the past decades this technology has evolved from sole hydrogen production, over indirect methane liquefaction to eventually direct plasma-assisted methane liquefaction processes. An overview of this evolution and these processes is presented, from which it becomes clear that the near future probably lies with the direct two phase plasma-assisted methane liquefaction and the far future with the direct oxidative methane liquefaction.

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: Journal of Physics D: Applied Physics published the first Plasma Roadmap in 2012
consisting of the individual perspectives of 16 leading experts in the various sub-fields of low temperature plasma science and technology. The 2017 Plasma Roadmap is the first update of a planned series of periodic updates of the Plasma Roadmap. The continuously growing interdisciplinary nature of the low temperature plasma field and its equally broad range of applications are making it increasingly difficult to identify major challenges that encompass all of the many sub-fields and applications. This intellectual diversity is ultimately a strength of the field. The current state of the art for the 19 sub-fields addressed in this roadmap demonstrates the enviable track record of the low temperature plasma field in the development of plasmas as an enabling technology for a vast range of technologies that underpin our modern society. At the same time, the many important scientific and technological challenges shared in this roadmap show that the path forward is not only scientifically rich but has the potential to make wide and far reaching contributions to many societal challenges.

Abstract: Since the advent of graphene, other 2D materials have garnered interest; notably the single element materials silicene, germanene, and stanene. Weinvestigate the ballistic current-voltage (I-V) characteristics of armchair silicene and stanene armchair nanoribbons (AXNRs with X = Si, Sn) using a combination of density functional theory and non-equilibrium Green's functions. The impact of out-of-plane electric field and in-plane uniaxial strain on the ribbon geometries, electronic structure, and (I-V)s are considered and contrasted with graphene. Since silicene and stanene are sp(2)/sp(3) buckled layers, the electronic structure can be tuned by an electric field that breaks the sublattice symmetry, an effect absent in graphene. This decreases the current by similar to 50% for Sn, since it has the largest buckling. Uniaxial straining of the ballistic channel affects the AXNR electronic structure in multiple ways: it changes the bandgap and associated effective carrier mass, and creates a local buckling distortion at the lead-channel interface which induces a interface dipole. Due to the increasing sp(3) hybridization character with increasing element mass, large reconstructions rectify the strained systems, an effect absent in sp(2) bonded graphene. This results in a smaller strain effect on the current: a decrease of 20% for Sn at 15% tensile strain compared to a similar to 75% decrease for C.

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.