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: Carbon nanotubes (CNTs) are nowadays routinely grown in a thermal CVD setup. State-of-the-art plasma-enhanced CVD (PECVD) growth, however, offers advantages over thermal CVD. A lower growth temperature and the growth of aligned freestanding single-walled CNTs (SWNTs) makes the technique very attractive. The atomic scale growth mechanisms of PECVD CNT growth, however, remain currently entirely unexplored. In this contribution, we employed molecular dynamics simulations to focus on the effect of applying an electric field on the SWNT growth process, as one of the effects coming into play in PECVD. Using sufficiently strong fields results in (a) alignment of the growing SWNTs, (b) a better ordering of the carbon network, and (c) a higher growth rate relative to thermal growth rate. We suggest that these effects are due to the small charge transfer occurring in the Ni/C system. These simulations constitute the first study of PECVD growth of SWNTs on the atomic level.

Abstract: In this contribution, we present an overview of the various techniques for combining atomistic molecular dynamics with Monte Carlo simulations, mainly in the context of condensed matter systems, as well as a brief summary of the main accelerated dynamics techniques. Special attention is given to the force bias Monte Carlo technique and its combination with molecular dynamics, in view of promising recent developments, including a definable timescale. Various examples of the application of combined molecular dynamics / Monte Carlo simulations are given, in order to demonstrate the enhanced simulation efficiency with respect to either pure molecular dynamics or Monte Carlo.

Abstract: Although the heterogeneity of pyrolyzed Fe???N???C materials is known and has been reported previously, the atomic structure of the active sites and their detailed reaction mechanisms are still unknown. Here, we identified two pyridinic Fe???N4-like centers with different local C coordinates, i.e., FeN4C8 and FeN4C10, and studied their electrocatalytic activity for the nitrogen reduction reaction (NRR) based on density functional theory (DFT) calculations. We also discovered the influence of the adsorption of NH2 as a functional ligand on catalyst performance on the NRR. We confirmed that the NRR selectivity of the studied catalysts is essentially governed either by the local C coordination or by the dynamic structure associated with the FeII/FeIII. Our investigations indicate that the proposed traditional pyridinic FeN4C10 has higher catalytic activity and selectivity for the NRR than the robust FeN4C8 catalyst, while it may have outstanding activity for promoting other (electro)catalytic reactions. <comment>Superscript/Subscript Available</comment

Abstract: Immobilization of single metal atoms on a solid host opens numerous possibilities for catalyst designs. If that host is a two-dimensional sheet, sheet curvature becomes a design parameter potentially complementary to host and metal composition. Here, we use a combination of density functional theory calculations and microkinetic modeling to compare the mechanisms and kinetics of formic acid decomposition and formation, chosen for their relevance as a potential hydrogen storage medium, over single Ru atoms anchored to pyridinic nitrogen in a planar graphene flake (RuN4-G) and curved carbon nanotube (RuN4-CNT). Activation barriers are lowered and the predicted turnover frequencies are increased over RuN4-CNT relative to RuN4-CNT. The results highlight the potential of curvature control as a means to achieve high performance and robust catalysts.

Abstract: Optimization of catalytic materials for a given application is greatly constrained by linear scaling relations. Recently, however, it has been demonstrated that it is possible to reversibly modulate the chemisorption of molecules on nanomaterials by charging (i.e., injection or removal of electrons) and hence reversibly and selectively modify catalytic activity beyond structure−activity correlations. The fundamental physical relation between the properties of the material, the charging process, and the chemisorption energy, however, remains unclear, and a systematic exploration and optimization of charge-switchable sorbent materials is not yet possible. Using hybrid DFT calculations of CO2 chemisorption on hexagonal boron nitride nanosheets with several types of defects and dopants, we here reveal the existence of fundamental correlations between the electron affinity of a material and charge-induced chemisorption, show how defect engineering can be used to modulate the strength and efficiency of the adsorption process, and demonstrate that excess electrons stabilize many topological defects. We then show how these insights could be exploited in the development of new electrocatalytic materials and the synthesis of doped nanomaterials. Moreover, we demonstrate that calculated chemical properties of charged materials are highly sensitive to the employed computational methodology because of the self-interaction error, which underlines the theoretical challenge posed by such systems.

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: To gain insight into the nature of the adhesion mechanism between hydroxyapatite (HA) and rutile (rTiO(2)), the mutual affinity between their surfaces was systematically studied using density functional theory (DFT). We calculated both bulk and surface properties of HA and rTiO(2), and explored the interfacial bonding mechanism of amorphous HA (aHA) surface onto amorphous as well as stoichiometric and nonstoichiometric crystalline rTiO(2). Formation energies of bridging and subbridging oxygen vacancies considered in the rTiO(2)(110) surface were evaluated and compared with other theoretical and experimental results. The interfacial interaction was evaluated through the work of adhesion. For the aHA/rTiO(2)(110) interfaces, the work of adhesion is found to depend strongly on the chemical environment of the rTiO(2)(110) surface. Electronic analysis indicates that the charge transfer is very small in the case of interface formation between aHA and crystalline rTiO(2)(110). In contrast, significant charge transfer occurs between aHA and amorphous rTiO(2) (aTiO(2)) slabs during the formation of the interface. Charge density difference (CDD) analysis indicates that the dominant interactions in the interface have significant covalent character, and in particular the Ti-O and Ca-O bonds. Thus, the obtained results reveal that the aHA/aTiO(2) interface shows a more preferable interaction and is thermodynamically more stable than other interfaces. These results are particularly important for improving the long-term stability of HA-based implants.

Abstract: Device performance of Si|SiOx core-shell based nanowires critically depends on the exact control over the oxide thickness. Low-temperature plasma oxidation is a highly promising alternative to thermal oxidation allowing for improved control over the oxidation process, in particular for ultrasmall Si nanowires. We here elucidate the room temperature plasma oxidation mechanisms of ultrasmall Si nanowires using hybrid molecular dynamics / force-bias Monte Carlo simulations. We demonstrate how the oxidation and concurrent water formation mechanisms are a function of the oxidizing plasma species and we demonstrate how the resulting core-shell oxide thickness can be controlled through these species. A new mechanism of water formation is discussed in detail. The results provide a detailed atomic level explanation of the oxidation process of highly curved Si surfaces. These results point out a route toward plasma-based formation of ultrathin core-shell Si|SiOx nanowires at room temperature.

Abstract: Plasma catalysis is gaining increasing interest for various environmental applications, but the crucial question is whether plasma can be created inside catalyst pores and under which conditions. In practice, various catalytic support materials are used, with various dielectric constants. We investigate here the influence of the dielectric constant on the plasma properties inside catalyst pores and in the sheath in front of the pores, for various pore sizes. The calculations are performed by a two-dimensional fluid model for an atmospheric pressure dielectric barrier discharge in helium. The electron impact ionization rate, electron temperature, electron and ion density, as well as the potential distribution and surface charge density, are analyzed for a better understanding of the discharge behavior inside catalyst pores. The results indicate that, in a 100 μm pore, the electron impact ionization in the pore, which is characteristic for the plasma generation inside the pore, is greatly enhanced for dielectric constants below 300. Smaller pore sizes only yield enhanced ionization for smaller dielectric constants, i.e., up to εr = 200, 150, and 50 for pore sizes of 50, 30, and 10 μm. Thus, the most common catalyst supports, i.e., Al2O3 and SiO2, which have dielectric constants around εr = 8−11 and 4.2, respectively, should allow more easily that microdischarges can be formed inside catalyst pores, even for smaller pore sizes. On the other hand, ferroelectric materials with dielectric constants above 300 never seem to yield plasma enhancement inside catalyst pores, not even for 100 μm pore sizes. Furthermore, it is clear that the dielectric constant of the material has a large effect on the extent of plasma enhancement inside the catalyst pores, especially in the range between εr = 4 and εr = 200. The obtained results are explained in detail based on the surface charge density at the pore walls,
and the potential distribution and electron temperature inside and above the pores. The results obtained with this model are
important for plasma catalysis, as the production plasma species in catalyst pores might affect the catalyst properties, and thus
improve the applications of plasma catalysis.

Abstract: A method has been developed for damage-free cesium (Cs) encapsulation within single-walled carbon nanotubes (SWNTs) with fine position selectivity. Precise energy tuning of Cs-ion irradiation revealed that there is a clear energy window (2060 eV) for the efficient encapsulation of Cs through the hexagonal network of SWNT sidewalls without causing significant damage. This minimum energy threshold of Cs-ion encapsulation (∼20 eV) matches well with the value obtained by ab initio simulation (∼22 eV). Furthermore, position-selective Cs encapsulation was carried out, resulting in the successful formation of pn-junction SWNT thin films with excellent environmental stability.

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: Nonthermal atmospheric pressure plasmas are gaining increasing attention for biomedical applications. However, very little fundamental information on the interaction mechanisms between the plasma species and biological cells is currently available. We investigate the interaction of important plasma species, such as OH, H2O2, O, O3, as well as O2 and H2O, with bacterial peptidoglycan by means of reactive molecular dynamics simulations, aiming for a better understanding of plasma disinfection. Our results show that OH, O, O3, and H2O2 can break structurally important bonds of peptidoglycan (i.e., CO, CN, or CC bonds), which consequently leads to the destruction of the bacterial cell wall. The mechanisms behind these breakups are, however, dependent on the impinging plasma species, and this also determines the effectiveness of the cell wall destruction.

Abstract: The exact structure and properties of the Si vertical bar SiO2 interface are very important in microelectronics and photovoltaic devices such as metal-oxide-semiconductor field-effect transistors (MOSFETs) and solar cells. Whereas Si vertical bar SiO2 structures are traditionally produced by thermal oxidation, hyperthermal oxidation shows a number of promising advantages. However, the Si vertical bar SiO2 interface induced in hyperthermal Si oxidation has not been properly investigated yet. Therefore, in this work, the interface morphology and interfacial stresses during hyperthermal oxidation at room temperature are studied using reactive molecular dynamics simulations based on the ReaxFF potential. Interface thickness and roughness, as well as the bond length and bond angle distributions in the interface are discussed and compared with other models developed for the interfaces induced by traditional thermal oxidation. The formation of a compressive stress is observed. This compressive stress, which at the interface amounts about 2 GPa, significantly slows down the inward silica growth. This value is close to the experimental value in the Si vertical bar SiO2 interface obtained in traditional thermal oxidation.

Abstract: Using reactive molecular dynamics simulations by means of the ReaxFF potential, we studied the growth mechanism of ultrathin silica (SiO2) layers during hyperthermal oxidation at room temperature. Oxidation of Si(100){2 × 1} surfaces by both atomic and molecular oxygen was investigated in the energy range 15 eV. The oxidation mechanism, which differs from thermal oxidation, is discussed. In the case of oxidation by molecular O2, silica is quickly formed and the thickness of the formed layers remains limited compared to oxidation by atomic oxygen. The Si/SiO2 interfaces are analyzed in terms of partial charges and angle distributions. The obtained structures of the ultrathin SiO2 films are amorphous, including some intrinsic defects. 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.