Abstract: We explored the potential of plasma-based In-Situ Resource Utilization (ISRU) for Mars through the conversion of Martian atmosphere (~96% CO2, 2% N2, and 2% Ar) into life-sustaining chemicals. As the Martian surface pressure is about 1% of the Earth’s surface pressure, it is an ideal environment for plasma-based gas conversion using microwave reactors. At 1000 W and 10 Ln/min (normal liters per minute), we produced ~76 g/h of O2 and ~3 g/h of NOx using a 2.45 GHz waveguided reactor at 25 mbar, which is ~3.5 times Mars ambient pressure. The energy cost required to produce O2 was ~0.013 kWh/g, which is very promising compared to recently concluded MOXIE experiments on the Mars surface. This marks a crucial step towards realizing the extension of human exploration.

Abstract: We explain the basic operation of a nanowire pinch-off FET and graphene nanoribbon tunnelFET. For the nanowire pinch-off FET we construct an analytical model to obtain the threshold voltage as a function of radius and doping density. We use the gradual channel approximation to calculate the current-voltage characteristics of this device and we show that the nanowire pinch-off FET has a subthreshold slope of 60 mV/dec and good ION and ION/IOFF ratios. For the graphene nanoribbon tunnelFET we show that an improved analytical model yields more realistic results for the transmission probability and hence the tunneling current. The first simulation results for the graphene nanoribbon tunnelFET show promising subthreshold slopes.

Abstract: We evaluate the importance of quantum effects on the sticking of H2 on Al(110) for conditions that are close to those of molecular beam experiments that have been done on this system. Calculations with the quasi-classical trajectory (QCT) method and with quantum dynamics (QD) are performed using a model in which only motion in the six molecular degrees of freedom is allowed. The potential energy surface used has a minimum barrier height close to the value recently obtained with the quantum Monte Carlo method. Monte Carlo averaging over the initial rovibrational states allowed the QD calculations to be done with an order of magnitude smaller computational expense. The sticking probability curve computed with QD is shifted to lower energies relative to the QCT curve by 0.21 to 0.05 kcal/mol, with the highest shift obtained for the lowest incidence energy. Quantum effects are therefore expected to play a small role in calculations that would evaluate the accuracy of electronic structure methods for determining the minimum barrier height to dissociative chemisorption for H2 + Al(110) on the basis of the standard procedure for comparing results of theory with molecular beam experiments.

Abstract: We discuss the synthesis of carbon nanotubes (CNTs) and graphene by catalytic chemical vapour deposition (CCVD) and plasma-enhanced CVD (PECVD), summarising the state-of-the-art understanding of mechanisms controlling their growth rate, chiral angle, number of layers (walls), diameter, length and quality (defects), before presenting a new model for 2D nucleation of a graphene sheet from amorphous carbon on a nickel surface. Although many groups have modelled this process using a variety of techniques, we ask whether there are any complementary ideas emerging from the different proposed growth mechanisms, and whether different modelling techniques can give the same answers for a given mechanism. Subsequently, by comparing the results of tight-binding, semi-empirical molecular orbital theory and reactive bond order force field calculations, we demonstrate that graphene on crystalline Ni(111) is thermodynamically stable with respect to the corresponding amorphous metal and carbon structures. Finally, we show in principle how a complementary heterogeneous nucleation step may play a key role in the transformation from amorphous carbon to graphene on the metal surface. We conclude that achieving the conditions under which this complementary crystallisation process can occur may be a promising method to gain better control over the growth processes of both graphene from flat metal surfaces and CNTs from catalyst nanoparticles.

Abstract: We discuss the selection and assessment of alternative metals by a combination of ab initio computation of electronic properties, experimental resistivity assessments, and calibrated line resistance models. Pt-group metals as well as Nb are identified as the most promising elements, with Ru showing the best combination of material properties and process maturity. An experimental assessment of the resistivity of Ru, Ir, and Co lines down to similar to 30 nm(2) is then used to devise compact models for line and via resistance that can be compared to Cu predictions. The main advantage of alternative metals originates from the possibility for barrierless metallization.

Abstract: We developed a microkinetic model to reveal the effects of plasma-generated radicals, intermediates, and vibrationally excited species on the catalytic hydrogenation of CO2 to CH3OH on a Cu(111) surface. As a benchmark, we first present the mechanisms of thermal catalytic CH3OH formation. Our model predicts that the reverse water-gas shift reaction followed by CO hydrogenation, together with the formate path, mainly contribute to CH3OH formation in thermal catalysis. Adding plasma-generated radicals and intermediates results in a higher CH3OH turnover frequency (TOF) by six to seven orders of magnitude, showing the potential of plasma-catalytic CO2 hydrogenation into CH3OH, in accordance with the literature. In addition, CO2 vibrational excitation further increases the CH3OH TOF, but the effect is limited due to relatively low vibrational temperatures under typical plasma catalysis conditions. The predicted increase in CH3OH formation by plasma catalysis is mainly attributed to the increased importance of the formate path. In addition, the conversion of plasma-generated CO to HCO* and subsequent HCOO* or H2CO* formation contribute to CH3OH formation. Both pathways bypass the HCOO* formation from CO2, which is the main bottleneck in the process. Hence, our model points toward the important role of CO, but also O, OH, and H radicals, as they influence the reactions that consume CO2 and CO. In addition, our model reveals that the H pressure should not be smaller than ca. half of the O pressure in the plasma as this would cause O* poisoning, which would result in very small product TOFs. Thus, plasma conditions should be targeted with a high CO and H content as this is favorable for CH3OH formation, while the O content should be minimized.

Abstract: We derive an analytical model for the electrostatics and the drive current in a silicon nanowire operating in JFET mode. We show that there exists a range of nanowire radii and doping densities for which the nanowire JFET satisfies reasonable device characteristics. For thin nanowires we have developed a self-consistent quantum mechanical model to obtain the electronic structure.

Abstract: We apply an implicit particle-in-cell Monte-Carlo (PIC-MC) method to study a radio-frequency argon microdischarge at steady state in the glow discharge limit, in which the microdischarge is sustained by secondary electron emission from the electrodes. The plasma density, electron energy distribution function (EEDF), and electron temperature are calculated in a wide range of operating conditions, including driving voltage, microdischarge gap, and pressure. Also, the effect of gap size scaling (in the range of 50-1000 μm) on the plasma sustaining voltage and peak electron density at atmospheric pressure is examined, which has not been explored before. In our simulations, three different EEDFs, i.e., a so-called three temperature hybrid mode, a two temperature α mode, and a two temperature γ mode distribution, are identified at different gaps and voltages. The maximum sustaining voltage to avoid a transition from the glow mode to an arc is predicted, as well as the minimum sustaining voltage for a steady glow discharge. Our calculations elucidate that secondary electrons play an essential role in sustaining the discharge, and as a result the relationship between breakdown voltage and gap spacing is far away from the Paschen law at atmospheric pressure.

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: Vibrational excitation represents an efficient channel to drive the dissociation of CO₂in a non-thermal plasma. Its viability is investigated in low-pressure pulsed discharges, with the intention of selectively exciting the asymmetric stretching mode, leading to stepwise excitation up to the dissociation limit of the molecule. Gas heating is crucial for the attainability of this process, since the efficiency of vibration–translation (V–T) relaxation strongly depends on temperature, creating a feedback mechanism that can ultimately thermalize the discharge. Indeed, recent experiments demonstrated that the timeframe of V–T non-equilibrium is limited to a few milliseconds at ca. 6 mbar, and shrinks to the<italic>μ</italic>s-scale at 100 mbar. With the aim of backtracking the origin of gas heating in pure CO₂plasma, we perform a kinetic study to describe the energy transfers under typical non-thermal plasma conditions. The validation of our kinetic scheme with pulsed glow discharge experiments enables to depict the gas heating dynamics. In particular, we pinpoint the role of vibration–vibration–translation relaxation in redistributing the energy from asymmetric to symmetric levels of CO₂, and the importance of collisional quenching of CO₂electronic states in triggering the heating feedback mechanism in the sub-millisecond scale. This latter finding represents a novelty for the modelling of low-pressure pulsed discharges and we suggest that more attention should be paid to it in future studies. Additionally, O atoms convert vibrational energy into heat, speeding up the feedback loop. The efficiency of these heating pathways, even at relatively low gas temperature and pressure, underpins the lifetime of V–T non-equilibrium and suggests a redefinition of the optimal conditions to exploit the ‘ladder-climbing’ mechanism in CO₂discharges.

Abstract: Utilizing ammonia as a storage medium for hydrogen is currently receiving increased attention. A possible method to retrieve the hydrogen is by plasma-catalytic decomposition. In this work, we combined an experimental study, using a dielectric barrier discharge plasma reactor, with a plasma kinetic model, to get insights into the decomposition mechanism. The experimental results revealed a similar effect on the ammonia conversion when changing the flow rate and power, where increasing the specific energy input (higher power or lower flow rate) gave an increased conversion. A conversion as high as 82 % was achieved at a specific energy input of 18 kJ/Nl. Furthermore, when changing the discharge volume from 31 to 10 cm3, a change in the plasma distribution factor from 0.2 to 0.1 was needed in the model to best describe the conversions of the experiments. This means that a smaller plasma volume caused a higher transfer of energy through micro-discharges (non-uniform plasma), which was found to promote the decomposition of ammonia. These results indicate that it is the collisions between NH3 and the high-energy electrons that initiate the decomposition. Moreover, the rate of ammonia destruction was found by the model to be in the order of 1022 molecules/(cm3 s) during the micro-discharges, which is 5 to 6 orders of magnitude higher than in the afterglows. A considerable re-formation of ammonia was found to take place in the afterglows, limiting the overall conversion. In addition, the model revealed that implementation of packing material in the plasma introduced high concentrations of surface-bound hydrogen atoms, which introduced an additional ammonia re-formation pathway through an Eley-Rideal reaction with gas phase NH2. Furthermore, a more uniform plasma is predicted in the presence of MgAl2O4, which leads to a lower average electron energy during micro-discharges and a lower conversion (37 %) at a comparable residence time for the plasma alone (51 %).

Abstract: Using reactive molecular dynamics simulations, we have investigated the effect of single-impact, low-energy (thermal-100 eV) bombardment of a Si(100){2 × 1} surface by atomic and molecular oxygen. Penetration probability distributions, as well as defect formation distributions, are presented as a function of the impact energy for both species. It is found that at low impact energy, defects are created chemically due to the chemisorption process in the top layers of the surface, while at high impact energy, additional defects are created by a knock-on displacement of Si. These results are of particular importance for understanding device performances of silica-based microelectronic and photovoltaic devices.

Abstract: Using reactive molecular dynamics simulations, the melting behavior of nickel-carbon nanoclusters is examined. The phase diagrams of icosahedral and Wulff polyhedron clusters are determined using both the Lindemann index and the potential energy. Formulae are derived for calculating the equilibrium constants and the solid and liquid fractions during a phase transition, allowing more rational determination of the melting temperature with respect to the arbitrary Lindemann value. These results give more insight into the properties of nickel-carbon nanoclusters in general and can specifically be very useful for a better understanding of the synthesis of carbon nanotubes using the catalytic chemical vapor deposition method.

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: Using a molecular dynamics model the crystallinity of MgxAlyOz thin films with a variation in the stoichiometry of the thin film is studied at operating conditions similar to the experimental operating conditions of a dual magnetron sputter deposition system. The films are deposited on a crystalline or amorphous substrate. The Mg metal content in the film ranged from 100% (i.e. MgO film) to 0% (i.e. Al2O3 film). The radial distribution function and density of the films are calculated. The results are compared with x-ray diffraction and transmission electron microscopy analyses of experimentally deposited thin films by the dual magnetron reactive sputtering process. Both simulation and experimental results show that the structure of the MgAlO film varies from crystalline to amorphous when the Mg concentration decreases. It seems that the crystalline MgAlO films have a MgO structure with Al atoms in between.