Abstract: Atomic scale simulations of the nucleation and growth of carbon nanotubes is essential for understanding their growth mechanism. In spite of over twenty years of simulation efforts in this area, limited progress has so far been made on addressing the role of the hydrocarbon growth precursor. Here we report on atomic scale simulations of cap nucleation of single-walled carbon nanotubes from hydrocarbon precursors. The presented mechanism emphasizes the important role of hydrogen in the nucleation process, and is discussed in relation to previously presented mechanisms. In particular, the role of hydrogen in the appearance of unstable carbon structures during in situ experimental observations as well as the initial stage of multi-walled carbon nanotube growth is discussed. The results are in good agreement with available experimental and quantum-mechanical results, and provide a basic understanding of the incubation and nucleation stages of hydrocarbon-based CNT growth at the atomic level.

Abstract: The cleaning of two-dimensional (2D) materials is an essential step in the fabrication of future devices, leveraging their unique physical, optical, and chemical properties. Part of these emerging 2D materials are transition metal dichalcogenides (TMDs). So far there is limited understanding of the cleaning of “monolayer” TMD materials. In this study, we report on the use of downstream H₂plasma to clean the surface of monolayer WS₂grown by MOCVD. We demonstrate that high-temperature processing is essential, allowing to maximize the removal rate of polymers and to mitigate damage caused to the WS₂in the form of sulfur vacancies. We show that low temperature in situ carbonyl sulfide (OCS) soak is an efficient way to resulfurize the material, besides high-temperature H₂S annealing. The cleaning processes and mechanisms elucidated in this work are tested on back-gated field-effect transistors, confirming that transport properties of WS₂devices can be maintained by the combination of H₂plasma cleaning and OCS restoration. The low-damage plasma cleaning based on H₂and OCS is very reproducible, fast (completed in a few minutes) and uses a 300 mm industrial plasma etch system qualified for standard semiconductor pilot production. This process is, therefore, expected to enable the industrial scale-up of 2D-based devices, co-integrated with silicon technology.

Abstract: This study reports on the possible effects of OH radical impact on the transmembrane domain 6 of P-glycoprotein, TM6, which plays a crucial role in drug binding in human cells. For the first time, we employ molecular dynamics (MD) simulations based on the self-consistent charge density functional tight binding (SCC-DFTB) method to elucidate the potential sites of fragmentation and mutation in this domain upon impact of OH radicals, and to obtain fundamental information about the underlying reaction mechanisms. Furthermore, we apply non-reactive MD simulations to investigate the long-term effect of this mutation, with possible implications for drug binding. Our simulations indicate that the interaction of OH radicals with TM6 might lead to the breaking of C-C and C-N peptide bonds, which eventually cause fragmentation of TM6. Moreover, according to our simulations, the OH radicals can yield mutation in the aromatic ring of phenylalanine in TM6, which in turn affects its structure. As TM6 plays an important role in the binding of a range of cytotoxic drugs with P-glycoprotein, any changes in its structure are likely to affect the response of the tumor cell in chemotherapy. This is crucial for cancer therapies based on reactive oxygen species, such as plasma treatment.

Abstract: Modeling the oxidation process of silicon nanowires through reactive force field based molecular dynamics simulations suggests that the formation of Si epoxide defects occurs both at the Si/SiOx interface and at the nanowire surface, whereas for flat surfaces, this defect is experimentally observed to occur only at the interface as a result of stress. In this paper, we argue that the increasing curvature stabilizes the defect at the nanowire surface, as suggested by our density functional theory calculations. The latter can have important consequences for the opto-electronic properties of thin silicon nanowires, since the epoxide induces an electronic state within the band gap. Removing the epoxide defect by hydrogenation is expected to be possible but becomes increasingly difficult with a reduction of the diameter of the nanowires.

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: We performed molecular dynamics simulations to investigate the effect of lipid peroxidation products on the structural and dynamic properties of the cell membrane. Our simulations predict that the lipid order in a phospholipid bilayer, as a model system for the cell membrane, decreases upon addition of lipid peroxidation products. Eventually, when all phospholipids are oxidized, pore formation can occur. This will allow reactive species, such as reactive oxygen and nitrogen species (RONS), to enter the cell and cause oxidative damage to intracellular macromolecules, such as DNA or proteins. On the other hand, upon increasing the cholesterol fraction of lipid bilayers, the cell membrane order increases, eventually reaching a certain threshold, from which cholesterol is able to protect the membrane against pore formation. This finding is crucial for cancer treatment by plasma technology, producing a large number of RONS, as well as for other cancer treatment methods that cause an increase in the concentration of extracellular RONS. Indeed, cancer cells contain less cholesterol than their healthy counterparts. Thus, they will be more vulnerable to the consequences of lipid peroxidation, eventually enabling the penetration of RONS into the interior of the cell, giving rise to oxidative stress, inducing pro-apoptotic factors. This provides, for the first time, molecular level insight why plasma can selectively treat cancer cells, while leaving their healthy counterparts undamaged, as is indeed experimentally demonstrated.

Abstract: Doped Cr2O3 has been shown to be a p-type transparent conducting oxide (TCO). Its conductivity, however, is low. As for most p-type TCOs, the main problem is the high effective hole mass due to flat valence bands. We use first-principles methods to investigate whether one can increase the valence band dispersion (i.e. reduce the hole mass) by anion alloying with sulfur, while keeping the band gap large enough for transparency. The alloying concentrations considered are given by Cr(4)SxO(6-x), with x = 1-5. To be able to describe the electronic properties of these materials accurately, we first study Cr2O3, examining critically the accuracy of different density functionals and methods, including PBE, PBE+U, HSE06, as well as perturbative approaches within the GW approximation. Our results demonstrate that Cr4S2O4 has an optical band gap of 3.08 eV and an effective hole mass of 1.8 m(e). This suggests Cr4S2O4 as a new p-type TCO host candidate.

Abstract: Recently, Sr-doped LaCrO3 has been experimentally introduced as a new p-type transparent conducting oxide. It is demonstrated that substituting Sr for La results in inducing p-type conductivity in LaCrO3. Performing first principles calculations we study the electronic structure and formation energy of various point defects in LaCrO3. Our results for the formation energies show that in addition to Sr, two more divalent defects, Ca and Ba, substituting for La in LaCrO3, behave as shallow acceptors in line with previous experimental reports. We further demonstrate that under oxygen-poor growth conditions, these shallow acceptors will be compensated by intrinsic donor-like defects (an oxygen vacancy and Cr on an oxygen site), but in the oxygen-rich growth regime the shallow acceptors have the lowest formation energies between all considered defects and will lead to p-type conductivity.

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: Atomic scale simulations have been shown to be a powerful tool for elucidating the growth mechanisms of carbon nanotubes. The growth picture is however not entirely clear yet due to the gap between current simulations and real experiments. We here simulate for the first time the nucleation and subsequent growth of single-wall carbon nanotubes (SWNTs) from oxygen-containing hydrocarbon feedstocks using the hybrid Molecular Dynamics/Monte Carlo technique. The underlying nucleation mechanisms of Ni-catalysed SWNT growth are discussed in detail. Specifically, we find that as a function of the feedstock, different carbon fractions may emerge as the main growth species, due to a competition between the feedstock decomposition, its rehydroxylation and its contribution to etching of the growing SWNT. This study provides a further understanding of the feedstock effects in SWNT growth in comparison with available experimental evidence as well as with<italic>ab initio</italic>and other simulation data, thereby reducing the simulation–experiment gap.

Abstract: While spectroscopic data on small hydrocarbons in interstellar media in combination with crossed molecular beam (CMB) experiments have provided a wealth of information on astrochemically relevant species, much of the underlying mechanistic pathways of their formation remain elusive. Therefore, in this work, the chemical reaction mechanisms of C(³P_J) + C₆H₆and C⁺(²P) + C₆H₆systems using the quantum mechanical molecular dynamics (QMMD) technique at the PBE0-D3(BJ) level of theory is investigated, mimicking a CMB experiment. Both the dynamics of the reactions as well as the electronic structure for the purpose of the reaction network are evaluated. The method is validated for the first reaction by comparison to the available experimental data. The reaction scheme for the C(³P_J) + C₆H₆system covers the literature data,<italic>e.g.</italic>the major products are the 1,2-didehydrocycloheptatrienyl radical (C₇H₅) and benzocyclopropenyl radical (C₆H₅–CH), and it reveals the existence of less common pathways for the first time. The chemistry of the C⁺(²P_J) + C₆H₆system is found to be much richer, and we have found that this is because of more exothermic reactions in this system in comparison to those in the C(³P_J) + C₆H₆system. Moreover, using the QMMD simulation, a number of reaction paths have been revealed that produce three distinct classes of reaction products with different ring sizes. All in all, at all the collision energies and orientations, the major product is the heptagon molecular ion for the ionic system. It is also revealed that the collision orientation has a dominant effect on the reaction products in both systems, while the collision energy mostly affects the charged system. These simulations both prove the applicability of this approach to simulate crossed molecular beams, and provide fundamental information on reactions relevant for the interstellar medium.

Abstract: Molecular dynamics simulations are used to investigate the effect of an AFM tip when indenting graphene nano bubbles filled by a noble gas (i.e. He, Ne and Ar) up to the breaking point. The failure points resemble those of viral shells as described by the Foppl-von Karman (FvK) dimensionless number defined in the context of elasticity theory of thin shells. At room temperature, He gas inside the bubbles is found to be in the liquid state while Ne and Ar atoms are in the solid state although the pressure inside the nano bubble is below the melting pressure of the bulk. The trapped gases are under higher hydrostatic pressure at low temperatures than at room temperature.

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: We calculate bubble nucleation rates in a Lennard-Jones fluid through explicit molecular dynamics simulations. Our approach-based on a recent free energy method (dubbed reweighted Jarzynski sampling), transition state theory, and a simple recrossing correction-allows us to probe a fairly wide range of rates in several superheated and cavitation regimes in a consistent manner. Rate predictions from this approach bridge disparate independent literature studies on the same model system. As such, we find that rate predictions based on classical nucleation theory, direct brute force molecular dynamics simulations, and seeding are consistent with our approach and one another. Published rates derived from forward flux sampling simulations are, however, found to be outliers. This study serves two purposes: First, we validate the reliability of common modeling techniques and extrapolation approaches on a paradigmatic problem in materials science and chemical physics. Second, we further test our highly generic recipe for rate calculations, and establish its applicability to nucleation processes.

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.