Abstract: Selective etching allows for obtaining carbon nanotubes with a specific chirality. While plasma-assisted etching has already been used to separate metallic tubes from their semiconducting counterparts, little is known about the nanoscale mechanisms of the etching process. We combine (reactive) molecular dynamics (MD) and force-bias Monte Carlo (tfMC) simulations to study H-etching of CNTs. In particular, during the hydrogenation and subsequent etching of both the carbon cap and the tube, they sequentially transform to different carbon nanostructures, including carbon nanosheet, nanowall, and polyyne chains, before they are completely removed from the surface of a substrate-bound Ni-nanocluster.We also found that onset of the etching process is different in the cases of the cap and the tube, although the overall etching scenario is similar in both cases. The entire hydrogenation/etching process for both cases is analysed in detail, comparing with available theoretical and experimental evidences.

Abstract: In this study, we present a computational model of a cylindrical electric probe in atmospheric pressure argon plasma. The plasma properties are varied in terms of density and electron temperature. Furthermore, results for plasmas with Maxwellian and non-Maxwellian electron energy distribution functions are also obtained and compared. The model is based on the fluid description of plasma within the COMSOL software package. The results for the ion saturation current are compared and show good agreement with existing analytical Langmuir probe theories. A strong dependence between the ion saturation current and electron transport properties was observed, and attributed to the effects of ambipolar diffusion.

Abstract: In this study we quantitatively investigate for the first time the plasma characteristics of an argon gliding arc with a 3D model. The model is validated by comparison with available experimental data from literature and a reasonable agreement is obtained for the calculated gas temperature and electron density. A complete arc cycle is modeled from initial ignition to arc decay. We investigate how the plasma characteristics, i.e., the electron temperature, gas temperature,
reduced electric field, and the densities of electrons, Ar+ and Ar2+ ions and Ar(4s) excited states, vary over one complete arc cycle, including their behavior in the discharge and post-discharge. These plasma characteristics exhibit a different evolution over one arc cycle, indicating that either the active discharge stage or the post-discharge stage can be beneficial for certain applications.

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: The formation process of a microdischarge (MD) in both μm- and nm-sized catalyst pores is simulated by a two-dimensional particle-in-cell/Monte Carlo collision model. A parallel-plate dielectric barrier discharge configuration in filamentary mode is considered in ambient air. The discharge is powered by a high voltage pulse. Our calculations reveal that a streamer can penetrate into the surface features of a porous catalyst and MDs can be formed inside both μm- and nm-sized pores, yielding ionization inside the pore. For the μm-sized pores, the ionization mainly occurs inside the pore, while for the nm-sized pores the ionization is strongest near and inside the pore. Thus, enhanced discharges near and inside the mesoporous catalyst are observed. Indeed, the maximum values of the electric field, ionization rate and electron density occur near and inside the pore. The maximum electric field and electron density inside the pore first increase when the pore size rises from 4 nm to 10 nm, and then they decrease for the 100 nm pore, due to
a more pronounced surface discharge for the smaller pores. However, the ionization rate is highest for the 100 nm pore due to the largest effective ionization region.

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: 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: A gliding arc plasma is a potential way to convert CO2 into CO and O2, due to its non-equilibrium character, but little is known about the underlying mechanisms. In this paper, a self-consistent two-dimensional (2D) gliding arc model is developed, with a detailed non-equilibrium CO2 plasma chemistry, and validated with experiments. Our calculated values of the electron number density in the plasma, the CO2 conversion and energy efficiency show reasonable agreement with the experiments, indicating that the model can provide a realistic picture of the plasma chemistry. Comparison of the results with classical thermal conversion, as well as other plasma-based technologies for CO2 conversion reported in literature, demonstrates the non-equilibrium character of the gliding arc, and indicates that the gliding arc is a promising plasma reactor for CO2 conversion. However, some process modifications should be exploited to further improve its performance. As the model provides a realistic picture of the plasma behaviour, we use it first to investigate the plasma characteristics in a whole gliding arc cycle, which is necessary to understand the underlying mechanisms. Subsequently, we perform a chemical kinetics analysis, to investigate the different pathways for CO2 loss and formation. Based on the revealed discharge properties and the underlying CO2 plasma chemistry, the model allows us to propose solutions on how to further improve the
CO2 conversion and energy efficiency by a gliding arc plasma.

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: 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: Packed bed plasma reactors (PBPRs) are gaining increasing interest for use in environmental applications, such as greenhouse gas conversion into value-added chemicals or renewable fuels and volatile pollutant removal (e.g. NOx, VOC, K), as they enhance the conversion and energy efficiency of the process compared to a non-packed reactor. However, the plasma behaviour in a PBPR is not well understood. In this paper we demonstrate, by means of a fluid model, that the discharge behaviour changes considerably when changing the size of the packing beads and their dielectric constant, while keeping the interelectrode spacing constant. At low dielectric constant, the plasma is spread out over the full discharge gap, showing significant density in the voids as well as in the connecting void channels. The electric current profile shows a strong peak during each half cycle. When the dielectric constant increases, the plasma becomes localised in the voids, with a current profile consisting of many smaller peaks during each half cycle. For large bead sizes, the shift from full gap discharge to localised discharges takes place at a higher dielectric constant than for smaller beads. Furthermore, smaller beads or beads with a lower dielectric constant require a higher breakdown voltage to cause plasma formation.

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: Microdischarges are an important type of plasma discharge that possess several unique characteristics, such as the presence of a stable glow discharge, high plasma density and intense excimer radiation, leading to several potential applications. The intense and controllable gas heating within the extremely small dimensions of microdischarges has been exploited in microthruster technologies by incorporating a micro-nozzle to generate the thrust. This kind of microthruster has a significantly improved specific impulse performance compared to conventional cold gas thrusters, and can meet the requirements arising from the emerging development and application of micro-spacecraft. In this paper, we performed a self-consistent 2D particle-in-cell simulation, with a Monte Carlo collision model, of a microdischarge operating in a prototype micro-plasma thruster with a hollow cylinder geometry and a divergent micro-nozzle. The model takes into account the thermionic electron emission including the Schottky effect, the secondary electron emission due to cathode bombardment by the plasma ions, several different collision processes, and a non-uniform argon background gas density in the cathode-anode gap. Results in the high-pressure (several hundreds of Torr), high-current (mA) operating regime showing the behavior of the plasma density, potential distribution, and energy flux towards the hollow cathode and anode are presented and discussed. In addition, the results of simulations showing the effect of different argon gas pressures, cathode material work function and discharge voltage on the operation of the microdischarge thruster are presented. Our calculated properties are compared with experimental data under similar conditions and qualitative and quantitative agreements are reached.

Abstract: The impact of the scaling of the channel length on the performances of metal-oxide-semiconductor field effect transistors, based on two-dimensional (2D) channel materials, is theoretically investigated, using density functional theory combined with the non-equilibrium Green's function method. It is found that the scaling of the channel length below 10nm leads to strong device performance degradations. Our simulations reveal that this degradation is essentially due to the tunneling current flowing between the source and the drain in these aggressively scaled devices. It is shown that this electron tunneling process is modulated by the effective mass of the 2D channel material, and sets the limit of the scaling in future transistor designs. (C) 2016 AIP Publishing LLC.

Abstract: The discharge properties of the plasma bulk flow near the surface of heated compound-materials strongly affects the kinetic layer parameters modeled and manifested in the Knudsen layer. This paper extends the widely used two-layer kinetic ablation model to the ablation controlled non-equilibrium discharge due to the fact that the local thermodynamic equilibrium (LTE) approximation is often violated as a result of the interaction between the plasma and solid walls. Modifications to the governing set of equations, to account for this effect, are derived and presented by assuming that the temperature of the electrons deviates from that of the heavy particles. The ablation characteristics of one typical material, polytetrafluoroethylene (PTFE) are calculated with this improved model. The internal degrees of freedom as well as the average particle mass and specific heat ratio of the polyatomic vapor, which strongly depends on the temperature, pressure and plasma non-equilibrium degree and plays a crucial role in the accurate determination of the ablation behavior by this model, are also taken into account. Our assessment showed the significance of including such modifications related to the non-equilibrium effect in the study of vaporization of heated compound materials in ablation controlled arcs. Additionally, a two-temperature magneto-hydrodynamic (MHD) model accounting for the thermal non-equilibrium occurring near the wall surface is developed and applied into an ablation-dominated discharge for an electro-thermal chemical launch device. Special attention is paid to the interaction between the non-equilibrium plasma and the solid propellant surface. Both the mass exchange process caused by the wall ablation and plasma species deposition as well as the associated momentum and energy exchange processes are taken into account. A detailed comparison of the results of the non-equilibrium model with those of an equilibrium model is presented. The non-equilibrium results show a non-equilibrium region near the plasma-wall interaction region and this indicates the need for the consideration of the influence of the possible departure from LTE in the plasma bulk on the determination of ablation rate.