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 growth mechanism and chirality formation of a single-walled carbon nanotube (SWNT) on a surface-bound nickel nanocluster are investigated by hybrid reactive molecular dynamics/force-biased Monte Carlo simulations. The validity of the interatomic potential used, the so-called ReaxFF potential, for simulating catalytic SWNT growth is demonstrated. The SWNT growth process was found to be in agreement with previous studies and observed to proceed through a number of distinct steps, viz., the dissolution of carbon in the metallic particle, the surface segregation of carbon with the formation of aggregated carbon clusters on the surface, the formation of graphitic islands that grow into SWNT caps, and finally continued growth of the SWNT. Moreover, it is clearly illustrated in the present study that during the growth process, the carbon network is continuously restructured by a metal-mediated process, thereby healing many topological defects. It is also found that a cap can nucleate and disappear again, which was not observed in previous simulations. Encapsulation of the nanoparticle is observed to be prevented by the carbon network migrating as a whole over the cluster surface. Finally, for the first time, the chirality of the growing SWNT cap is observed to change from (11,0) over (9,3) to (7,7). It is demonstrated that this change in chirality is due to the metal-mediated restructuring process.

Abstract: In this paper, we quantitatively investigate with atom probe tomography, the effect of temperature on the interfacial transition layer suboxide species due to the thermal oxidation of silicon. The chemistry at the interface was measured with atomic scale resolution, and the changes in chemistry and intermixing at the interface were identified on a nanometer scale. We find an increase of suboxide (SiOx) concentration relative to SiO2 and increased oxygen ingress with elevated temperatures. Our experimental findings are in agreement with reactive force field molecular dynamics simulations. This work demonstrates the direct comparison between atom probe derived chemical profiles and atomistic-scale simulations for transitional interfacial layer of suboxides as a function of temperature.

Abstract: In order to optimize the plasma-synthesis and modification process of carbon nanomaterials for applications such as nanoelectronics and energy storage, a deeper understanding of fundamental hydrogengraphite/graphene interactions is required. Atomistic simulations by Molecular Dynamics have proven to be indispensable to illuminate these phenomena. However, severe time-scale limitations restrict them to very fast processes such as reflection, while slow thermal processes such as surface diffusion and molecular desorption are commonly inaccessible. In this work, we could however reach these thermal processes for the first time at time-scales and surface temperatures (1000 K) similar to high-flux plasma exposure experiments during the simulation of multilayer graphene etching by 5 eV H ions. This was achieved by applying the Collective Variable-Driven Hyperdynamics biasing technique, which extended the inter-impact time over a range of six orders of magnitude, down to a more realistic ion-flux of 1023m
2s
1. The results show that this not only causes a strong shift from predominant ion-to thermally induced interactions, but also significantly affects the hydrogen uptake and surface evolution. This study thus elucidates H ion-graphite/graphene interaction mechanisms and stresses the importance of including long time-scales in atomistic simulations at high surface temperatures to understand the dynamics of the ion-surface system.

Abstract: The possibility of confinement effects inside a carbon nanotube provides new application opportunities, e.g., growth of novel carbon nanostructures. However, the understanding the precise role of catalystfeedstock in the nanostructure synthesis is still elusive. In our simulation-based study, we investigate the Ni-catalyzed growth mechanism of encapsulated carbon nanostructures, viz. double-wall carbon nanotube and graphene nanoribbon, from carbon and hydrocarbon growth precursors, respectively. Specifically, we find that the tube and ribbon growth is determined by a catalyst-vs-feedstock competition effect. We compare our results, i.e., growth mechanism and structure morphology with all available theoretical and experimental data. Our calculations show that all encapsulated nanostructures contain metal (catalyst) atoms and such structures are less stable than their pure counterparts. Therefore, we study the purification mechanism of these structures. In general, this study opens a possible route to the controllable synthesis of tubular and planar carbon nanostructures for today’s nanotechnology.

Abstract: The interaction of thermal and hyperthermal Ni ions with gas-phase C60 fullerene was investigated at two temperatures with classical molecular dynamics simulations using a recently developed interatomic many-body potential. The interaction between Ni and C60 is characterized in terms of the NiC60 binding sites, complex formation, and the collision and temperature induced deformation of the C60 cage structure. The simulations show how ion implantation theoretically allows the synthesis of both endohedral Ni@C60 and exohedral NiC60 metallofullerene complexes.

Abstract: Low-temperature growth of carbon nanotubes (CNTs) has been claimed to provide a route towards chiral-selective growth, enabling a host of applications. In this contribution, we employ reactive molecular dynamics simulations to demonstrate how plasma-based deposition allows such low-temperature growth. We first show how ion bombardment during the growth affects the carbon dissolution and precipitation process. We then continue to demonstrate how a narrow ion energy window allows CNT growth at 500 K. Finally, we also show how CNTs in contrast cannot be grown in thermal CVD at this low temperature, but only at high temperature, in agreement with experimental data. (C) 2013 Elsevier Ltd. All rights reserved.

Abstract: Despite significant progress in single walled carbon nanotube (SWCNT) production by plasma enhanced chemical vapor deposition (PECVD), the growth mechanism in this method is not clearly understood. We employ reactive molecular dynamics simulations to investigate how plasma-based deposition allows growth at low temperature. We first investigate the SWCNT growth mechanism at low and high temperatures under conditions similar to thermal CVD and PECVD. We then show how ion bombardment during the nucleation stage increases the carbon solubility in the catalyst at low temperature. Finally, we demonstrate how moderate energy ions sputter amorphous carbon allowing for SWCNT growth at 500 K. (C) 2013 Elsevier B. V. All rights reserved.

Abstract: Atomistic simulation methods for the quantification of free energies are in wide use. These methods operate by sampling the probability density of a system along a small set of suitable collective variables (CVs), which is, in turn, expressed in the form of a free energy surface (FES). This definition of the FES can capture the relative stability of metastable states but not that of the transition state because the barrier height is not invariant to the choice of CVs. Free energy barriers therefore cannot be consistently computed from the FES. Here, we present a simple approach to calculate the gauge correction necessary to eliminate this inconsistency. Using our procedure, the standard FES as well as its gauge-corrected counterpart can be obtained by reweighing the same simulated trajectory at little additional cost. We apply the method to a number of systems—a particle solvated in a Lennard-Jones fluid, a Diels–Alder reaction, and crystallization of liquid sodium—to demonstrate its ability to produce consistent free energy barriers that correctly capture the kinetics of chemical or physical transformations, and discuss the additional demands it puts on the chosen CVs. Because the FES can be converged at relatively short (sub-ns) time scales, a free energy-based description of reaction kinetics is a particularly attractive option to study chemical processes at more expensive quantum mechanical levels of theory.

Abstract: Uniform-acceptance force-bias Monte Carlo (fbMC) methods have been shown to be a powerful technique to access longer timescales in atomistic simulations allowing, for example, phase transitions and growth. Recently, a new fbMC method, the time-stamped force-bias Monte Carlo (tfMC) method, was derived with inclusion of an estimated effective timescale; this timescale, however, does not seem able to explain some of the successes the method. In this contribution, we therefore explicitly quantify the effective timescale tfMC is able to access for a variety of systems, namely a simple single-particle, one-dimensional model system, the Lennard-Jones liquid, an adatom on the Cu(100) surface, a silicon crystal with point defects and a highly defected graphene sheet, in order to gain new insights into the mechanisms by which tfMC operates. It is found that considerable boosts, up to three orders of magnitude compared to molecular dynamics, can be achieved for solid state systems by lowering of the apparent activation barrier of occurring processes, while not requiring any system-specific input or modifications of the method. We furthermore address the pitfalls of using the method as a replacement or complement of molecular dynamics simulations, its ability to explicitly describe correct dynamics and reaction mechanisms, and the association of timescales to MC simulations in general.

Abstract: Experiments show that the etch rate of Si with SF6 inductively coupled plasma (ICP) is significantly influenced by the absolute gas flow rate in the range of 50–600 sccm, with a maximum at around 200 sccm. Therefore, we numerically investigate the effects of the gas flow rate on the bulk plasma properties and on the etch rate, to obtain more insight in the underlying reasons of this effect. A hybrid Monte Carlo—fluid model is applied to simulate an SF6 ICP. It is found that the etch rate is influenced by two simultaneous effects: (i) the residence time of the gas and (ii) the temperature profile of the plasma in the ICP volume, resulting indeed in a maximum etch rate at 200 sccm.

Abstract: The rate, selectivity and efficiency of plasma-based conversion processes is strongly affected by nonequilibrium phenomena. High concentrations of vibrationally excited molecules are such a plasma-induced effect. It is frequently assumed that vibrationally excited molecules are important in plasma catalysis because their presence lowers the apparent activation energy of dissociative chemisorption reactions and thus increases the conversion rate. A detailed atomic-level understanding of vibrationally stimulated catalytic reactions in the context of plasma catalysis is however lacking. Here, we couple a recently developed statistical model of a plasma-induced vibrational nonequilibrium to molecular dynamics simulations, enhanced sampling methods, and machine learning techniques. We quantify the impact of a vibrational nonequilibrium on the dissociative chemisorption barrier of H2 and CH4 on nickel catalysts over a wide range of vibrational temperatures. We investigate the effect of surface structure and compare the role of different vibrational modes of methane in the dissociation process. For low vibrational temperatures, very high vibrational efficacies are found, and energy in bend vibrations appears to dominate the dissociation of methane. The relative impact of vibrational nonequilibrium is much higher on terrace sites than on surface steps. We then show how our simulations can help to interpret recent experimental results, and suggest new paths to a better understanding of plasma catalysis.

Abstract: Plasma medicine has been claimed to provide a novel route to heal wounds and regenerate skin, although very little is currently known about the elementary processes taking place. We carried out a series of ReaxFF-based reactive molecular dynamics simulations to investigate the interaction of O and OH radicals with lipids, more specifically with α-linolenic acid as a model for the free fatty acids present in the upper skin layer. Our calculations predict that the O and OH radicals most typically abstract a H atom from the fatty acids, which can lead to the formation of a conjugated double bond, but also to the incorporation of alcohol or aldehyde groups, thereby increasing the hydrophilic character of the fatty acids and changing the general lipid composition of the skin. Within the limitations of the investigated model, no formation of possibly toxic products was observed.