Abstract: The resistance of a TiO2 anatase (001) surface to coke formation was studied in the context of dry reforming of methane using density functional theory (DFT) calculations. As carbon atoms act as precursors for coke formation, the resistance to coke formation can be measured by the carbon coverage of the surface. This is related to the stability of different CHx (x = 0−3) species and their rate of hydrogenation and dehydrogenation on the TiO2 surface. Therefore, we studied the reaction mechanisms and their corresponding rates as a function of the temperature for the dehydrogenation of the species on the surface. We found that the stabilities of C and CH are significantly lower than those of CH3 and CH2. The hydrogenation rates of the different species are significantly higher than the dehydrogenation rates in a temperature range of 300−1000 K. Furthermore, we found that dehydrogenation of CH3, CH2, and CH will only occur at appreciable rates starting from 600, 900, and 900 K, respectively. On the basis of these results, it is clear that the anatase (001) surface has a high coke resistance, and it is thus not likely that the surface will become poisoned by coke during dry reforming of methane. As the rate limiting step in dry reforming is the dissociative adsorption of CH4, we studied an alternative approach to thermal catalysis. We found that the temperature threshold for dry reforming is at least 700 K. This threshold temperature may be lowered by the use of plasma-catalysis, where the appreciable rates of adsorption of plasma-generated CHx radicals result in bypassing the rate limiting step of the reaction.

Abstract: An inductively coupled plasma (ICP) is analyzed by means of experiments and numerical simulation. Important plasma properties are analyzed, namely, the effective temperature inside the central channel and the mean flow velocity inside the plasma. Furthermore, the effect of torches with different injector diameters is studied by the model. The temperature inside the central channel is determined from the end-on collected line-to-background ratio in dependence of the injector gas flow rates. Within the limits of 3% deviation, the results of the simulation and the experiments are in good agreement in the range of flow rates relevant for the analysis of relatively large droplets, i.e., 50 μm. The deviation increases for higher gas flow rates but stays below 6% for all flow rates studied. The velocity of the gas inside the coil region was determined by side-on analyte emission measurements with single monodisperse droplet introduction and by the analysis of the injector gas path lines in the simulation. In the downstream region significantly higher velocities were found than in the upstream region in both the simulation and the experiment. The quantitative values show good agreement in the downstream region. In the upstream region, deviations were found in the absolute values which can be attributed to the flow conditions in that region and because the methods used for velocity determination are not fully consistent. Eddy structures are found in the simulated flow lines. These affect strongly the way taken by the path lines of the injector gas and they can explain the very long analytical signals found in the experiments at low flow rates. Simulations were performed for different injector diameters in order to find conditions where good analyte transport and optimum signals can be expected. The results clearly show the existence of a transition flow rate which marks the lower limit for effective analyte transport conditions through the plasma. A rule-of-thumb equation was extracted from the results from which the transition flow rate can be estimated for different injector diameters and different injector gas compositions.

Abstract: We have investigated the plasma processes in an atmospheric pressure glow discharge (APGD) in He used for analytical spectrometry by means of fluid and Monte Carlo (MC) simulations. Typical results include the potential and electric field distributions in the plasma, the density profiles of the various plasma species throughout the discharge, the mean electron energy, as well as the rates of the various collision processes in the plasma, and the relative importance of the different production and loss rates for the various species. The similarities and differences with low-pressure glow discharges are discussed. The main differences are a very small cathode dark space region and a large positive column as well as the dominant role of molecular ions. Some characteristic features of the APGD, such as the occurrence of the different spatial zones in the discharge, are illustrated, with links to experimental observations.

Abstract: Using reactive molecular dynamics simulations based on the ReaxFF potential, we studied the growth mechanism of ultrathin silica (SiO2) layers during hyperthermal oxidation as a function of temperature in the range 100-1300 K. Oxidation of Si(100){2 x 1} surfaces by both atomic and molecular oxygen was investigated for hyperthermal impact energies in the range of 1 to 5 eV. Two different growth mechanisms are found, corresponding to a low temperature oxidation and a high temperature one. The transition temperature between these mechanisms is estimated to be about 700 K. Also, the initial step of the Si oxidation process is analyzed in detail. Where possible, we validated our results with experimental and ab initio data, and good agreement was obtained. This study is important for the fabrication of silica-based devices in the micro- and nanoelectronics industry and, more specifically, for the fabrication of metal oxide semiconductor devices.

Abstract: A number of recent computational material design studies based on density functional theory (DFT) calculations have put forward a new class of materials with electrically switchable chemical characteristics that can be exploited in the development of tunable gas storage and electrocatalytic applications. We find systematic flaws in almost every computational study of gas adsorption on polarized or charged surfaces, stemming from an improper and unreproducible treatment of periodicity, leading to very large errors of up to 3 eV in some cases. Two simple corrective procedures that lead to consistent results are proposed, constituting a crucial course correction to the research in the field.

Abstract: We report the formation of nanobubbles on graphene with a radius of the order of 1 nm, using ultralow energy implantation of noble gas ions (He, Ne, Ar) into graphene grown on a Pt(111) surface. We show that the universal scaling of the aspect ratio, which has previously been established for larger bubbles, breaks down when the bubble radius approaches 1 nm, resulting in much larger aspect ratios. Moreover, we observe that the bubble stability and aspect ratio depend on the substrate onto which the graphene is grown (bubbles are stable for Pt but not for Cu) and trapped element. We interpret these dependencies in terms of the atomic compressibility of the noble gas as well as of the adhesion energies between graphene, the substrate, and trapped atoms.

Abstract: To gain full control over the growth of carbon nanotubes (CNTs) using plasma-enhanced chemical vapor deposition (PECVD), a thorough understanding of the underlying plasma-catalyst mechanisms is required. Oxygen-containing species are often used as or added to the growth precursor gas, but these species also yield various radicals and ions, which may simultaneously etch the CNT during the growth. At present, the effect of these reactive species on the growth onset has not yet been thoroughly investigated. We here report on the etching mechanism of incipient CNT structures from OH and O radicals as derived from combined (reactive) molecular dynamics (MD) and force-bias Monte Carlo (tfMC) simulations. Our results indicate that the oxygen-containing radicals initiate a dissociation process. In particular, we show how the oxygen species weaken the interaction between the CNT and the nanocluster. As a result of this weakened interaction, the CNT closes off and dissociates from the cluster in the form of a fullerene. Beyond the specific systems studied in this work, these results are generically important in the context of PECVD-based growth of CNTs using oxygen-containing precursors.

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 recent energy crisis revealed that there is a strong need to replace hydrocarbon-fueled industrial nitrogen fixation processes by alternative, more sustainable methods. In light of this, plasma-based nitrogen fixation remains one of the most promising options, considering both theoretical and experimental aspects. Lately, plasma interacting with water has received considerable attention in nitrogen fixation applications as it can trigger a unique gas- and liquid-phase chemistry. Within this context, a critical exploration of plasma-assisted nitrogen fixation with or without water presence is of great interest with an emphasis on energy costs, particularly in plasma reactors which have potential for large-scale industrial application. In this work, the presence of water in a multi-pin plasma system on nitrogen oxidation is experimentally investigated by comparing two pulsed negative DC voltage plasmas in metal–metal and metal–liquid electrode configurations. The plasma setups are designed to create similar plasma properties, including plasma power and discharge regime in both configurations. The system energy cost is calculated, considering nitrogen-containing species generated in gas and liquid phases as measured by a gas analyzer, nitrate sensor, and a colorimetry method. The energy cost profile as a function of specific energy input showed a strong dependency on the plasma operational frequency and the gas flow rate, as a result of different plasma operation regimes and initiated reverse processes. More importantly, the presence of the plasma/liquid interface increased the energy cost up to 14 ± 8%. Overall, the results showed that the presence of water in the reaction zone has a negative impact on the nitrogen fixation process.

Abstract: Water permeation between stacked layers of hBN sheets forming 2D nanochannels is investigated using large-scale ab initio-quality molecular dynamics simulations. A high-dimensional neural network potential trained on density-functional theory calculations is employed. We simulate water in van der Waals nanocapillaries and study the impact of nanometric confinement on the structure and dynamics of water using both equilibrium and nonequilibrium methods. At an interlayer distance of 10.2 A confinement induces a first-order phase transition resulting in a well-defined AA-stacked bilayer of hexagonal ice. In contrast, for h < 9 A, the 2D water monolayer consists of a mixture of different locally ordered patterns of squares, pentagons, and hexagons. We found a significant change in the transport properties of confined water, particularly for monolayer water where the water-solid friction coefficient decreases to half and the diffusion coefficient increases by a factor of 4 as compared to bulk water. Accordingly, the slip-velocity is found to increase under confinement and we found that the overall permeation is dominated by monolayer water adjacent to the hBN membranes at extreme confinements. We conclude that monolayer water in addition to bilayer ice has a major contribution to water transport through 2D nanochannels.

Abstract: We investigate the theoretical possibility of single layer graphene formation on a nickel surface at different substrate temperatures under far-from-equilibrium high precursor flux conditions, employing state-of-the-art hybrid reactive molecular dynamics/uniform acceptance force bias Monte Carlo simulations. It is predicted that under these conditions, the formation of a single layer graphene-like film may proceed through a combined depositionsegregation mechanism on a nickel substrate, rather than by pure surface segregation as is typically observed for metals with high carbon solubility. At 900 K and above, nearly continuous graphene layers are obtained. These simulations suggest that single layer graphene deposition is theoretically possible on Ni under high flux conditions.

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: Nitrogen disproportionation i.e. its simultaneous conversion to compounds of higher (NOx) and lower (NH3) oxidation states in a N-2 DC plasma-driven electrolysis process with a plasma cathode is investigated. This type of plasma-liquid interaction exhibits a growing interest for many applications, in particular nitrogen fixation where it represents a green alternative to the Haber-Bosch process. Optical emission spectroscopy, FTIR and electrochemical sensing systems are used to characterize the gas phase physico-chemistry while the liquid phase is analyzed via ionic chromatography and colorimetric assays. Experiments suggest that lowering the discharge current enhances nitrogen reduction and facilitates the transfer of nitrogen compounds to the liquid phase. Large amounts of water vapor appear to impact the gas discharge physico-chemistry and to favor the vibrational excitation of N-2, a key parameter for an energy-efficient nitrogen fixation.