Abstract: Unlocking the enormous technological potential of carbon nanotubes strongly depends on our ability to specifically produce metallic or semiconducting tubes. While selective etching of both has already been demonstrated, the underlying reasons, however, remain elusive as yet. We here present computational and experimental evidence on the operative mechanisms at the atomic scale. We demonstrate that during the adsorption of H atoms and their coalescence, the adsorbed ortho hydrogen pairs on single-walled carbon nanotubes induce higher shear stresses than axial stresses, leading to the elongation of HC–CH bonds as a function of their alignment with the tube chirality vector, which we denote as the γ-angle. As a result, the C–C cleavage occurs more rapidly in nanotubes containing ortho H-pairs with a small γ-angle. This phenomenon can explain the selective etching of small-diameter semiconductor nanotubes with a similar curvature. Both theoretical and experimental results strongly indicate the important role of the γ-angle in the selective etching mechanisms of carbon nanotubes, in addition to the nanotube curvature and metallicity effects and lead us to clearly understand the onset of selective synthesis/removal of CNT-based materials.

Abstract: Heterostructures of low-dimensional semiconducting materials, such as transition metal dichalcogenides (MX2), are promising building blocks for future electronic and optoelectronic devices. The patterning of one MX2 material on top of another one is challenging due to their structural similarity. This prevents an intrinsic etch stop when conventional anisotropic dry etching processes are used. An alternative approach consist in a two-step process, where a sacrificial silicon layer is pre-patterned with a low damage plasma process, stopping on the underlying MoS2 film. The pre-patterned layer is used as sacrificial template for the formation of the top WS2 film. This study describes the optimization of a cyclic Ar/Cl-2 atomic layer etch process applied to etch silicon on top of MoS2, with minimal damage, followed by a selective conversion of the patterned Si into WS2. The impact of the Si atomic layer etch towards the MoS2 is evaluated: in the ion energy range used for this study, MoS2 removal occurs in the over-etch step over 1-2 layers, leading to the appearance of MoOx but without significant lattice distortions to the remaining layers. The combination of Si atomic layer etch, on top of MoS2, and subsequent Si-to-WS2 selective conversion, allows to create a WS2/MoS2 heterostructure, with clear Raman signals and horizontal lattice alignment. These results demonstrate a scalable, transfer free method to achieve horizontally individually patterned heterostacks and open the route towards wafer-level processing of 2D materials.

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: 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: Natural deep eutectic solvents (NADES) represent a green alternative to common organic solvents in the biochemical industry due to their benign behavior and tailorable properties, in particular as media for enzymatic reactions. However, to fully exploit their potential in enzymatic reactions, there is a need for a more fundamental understanding of how these neoteric solvents influence the course of these reac-tions. Thus, the aim of this study is to investigate the influence of NADES with various molar composi-tions on the stability and structure of enzymes, applying molecular dynamics simulations. This can help to better understand the effect of individual compounds of NADES, in addition to eutectic mixtures. More specifically, we simulate the behavior of Candida antarctica lipase B (CALB) enzyme in NADES com-posed of choline chloride with either urea, ethylene glycol or glycerol. Hereto, we monitor the NADES microstructure, the general stability of the enzyme and changes in the structure of its active sites and sur-face residues. Our simulations show that none of the studied NADES systems significantly disrupt the microstructure of the solvent or the stability of the CALB enzyme within the time scales of the simula-tions. The enzyme preserves its initial structure, size and intra-chain hydrogen bonds in all investigated compositions and, for the first time reported, also in NADES with increased hydrogen bond donating com-pound ratios. As the main novelty, our results indicate that, in addition to the composition, the molar ratio can be an additional variable to fine-tune the physicochemical properties of NADES without altering the enzyme characteristics. These findings could facilitate the development and application of task -tailored NADES media for biocatalytic processes. (c) 2022 Elsevier B.V. All rights reserved.

Abstract: Plasma catalysis is an emerging new technology for the electrification and downscaling of NH3 synthesis. Increasing attention is being paid to the optimization of plasma catalysis with respect to the plasma conditions, the catalyst material, and their mutual interaction. In this work we use microkinetic models to study how the total conversion process is impacted by the combination of different plasma conditions and transition metal catalysts. We study how plasma-generated radicals and vibrationally excited N2 (present in a dielectric barrier discharge plasma) interact with the catalyst and impact the NH3 turnover frequencies (TOFs). Both filamentary and uniform plasmas are studied, based on plasma chemistry models that provided plasma phase speciation and vibrational distribution functions. The Langmuir−Hinshelwood reaction rate coefficients (i.e., adsorption reactions and subsequent reactions among adsorbates) are determined using conventional scaling relations. An additional set of Eley−Rideal reactions (i.e., direct reactions of plasma radicals with adsorbates) was added and a sensitivity analysis on the assumed reaction rate coefficients was performed. We first show the impact of different vibrational distribution functions on the catalytic dissociation of N2 and subsequent production of NH3, and we gradually include more radical reactions, to illustrate the contribution of these species and their corresponding reaction pathways. Analysis over a large range of catalysts indicates that different transition metals (metals such as Rh, Ni, Pt, and Pd) optimize the NH3TOFs depending on the population of the vibrational levels of N2. At higher concentrations of plasma-generated radicals, the NH3 TOFs become less dependent on the catalyst material, due to radical adsorptions on the more noble catalysts and Eley−Rideal reactions on the less noble catalysts.

Abstract: In this work, the effects of N-doping into the Co-doped single vacancy (Co-SV-G) and di-vacancy graphene flake (Co-dV-G) are investigated and compared toward direct oxidation of methane to methanol (DOMM) employing two different oxidants (N2O and O-2) using density functional theory (DFT) calculation. We found that DOMM on CoN3-G utilizing the N2O molecule as oxygen-donor proceeds via a two-step reaction with low activation energies. In addition, we found that although CoN3-G might be a good catalyst for methane conversion, it can also catalyze the oxidation of methanol to CO2 and H2O due to the required low activation barriers. Moreover, the adsorption behaviors of CHx (x = 0-4) species and dehydrogenation of CHx (x = 1-4) species on CoN3-G are investigated. We concluded that CoN3-G can be used as an efficient catalyst for DOMM and N-2O reduction at ambient conditions which may serve as a guide for fabricating effective C/N catalysts in energy-related devices.

Abstract: The deposition of ultra-thin tungsten films and their related 2D chalcogen compounds on large area dielectric substrates by gas phase reactions is challenging. The lack of nucleation sites complicates the adsorption of W-related precursors and subsequent sulfurization usually requires high temperatures. We propose here a technique in which a thin solid amorphous silicon film is used as reductant for the gas phase precursor WF6 leading to the conversion to metallic W. The selectivity of the W conversion towards the underlying dielectric surfaces is demonstrated. The role of the Si surface preparation, the conversion temperature, and Si thickness on the formation process is investigated. Further, the in situ conversion of the metallic tungsten into thin stoichiometric WS2 is achieved by a cyclic approach based on WF6 and H2S pulses at the moderate temperature of 450 1C, which is much lower than usual oxide sulfurization processes.

Abstract: The hyperdynamics method is a powerful tool to simulate slow processes at the atomic level. However, the construction of an optimal hyperdynamics potential is a task that is far from trivial. Here, we propose a generally applicable implementation of the hyperdynamics algorithm, borrowing two concepts from metadynamics. First, the use of a collective variable (CV) to represent the accelerated dynamics gives the method a very large flexibility and simplicity. Second, a metadynamics procedure can be used to construct a suitable history-dependent bias potential on-the-fly, effectively turning the algorithm into a self-learning accelerated molecular dynamics method. This collective variable-driven hyperdynamics (CVHD) method has a modular design: both the local system properties on which the bias is based, as well as the characteristics of the biasing method itself, can be chosen to match the needs of the considered system. As a result, system-specific details are abstracted from the biasing algorithm itself, making it extremely versatile and transparent. The method is tested on three model systems: diffusion on the Cu(001) surface and nickel-catalyzed methane decomposition, as examples of reactive processes with a bond-length-based CV, and the folding of a long polymer-like chain, using a set of dihedral angles as a CV. Boost factors up to 109, corresponding to a time scale of seconds, could be obtained while still accurately reproducing correct dynamics.

Abstract: Magnetron sputtering is a widely used physical vapor deposition technique for deposition and formation of nanocatalyst thin films and clusters. Nevertheless, so far only few studies investigated this formation process at the fundamental level. We here review atomic scale molecular dynamics simulations aimed at elucidating the nanocatalyst growth process through magnetron sputtering. We first introduce the basic magnetron sputtering background and machinery of molecular dynamics simulations, and then describe the studies conducted in this field so far. We also present a perspective view on how the field may be developed further.

Abstract: The adsorption of C and CHx radicals on anatase (001) was studied using DFT within the generalized gradient approximation using the Perde-Burke-Ernzerhof (PBE) functional. We have studied the influence of oxygen vacancies in and at the surface on the adsorption properties of the radicals. For the oxygen vacancies in anatase (001), the most stable vacancy is located at the surface. For this vacancy, the maximal adsorption strength of C and CH decreases compared to the adsorption on the stoichiometric surface, but it increases for CH2 and CH3. If an oxygen vacancy is present in the first subsurface layer, the maximal adsorption strength increases for C, CH, CH2, and CH3. When the vacancy is present in the next subsurface layer, we find that only the CH3 adsorption is enhanced, while the maximal adsorption energies for the other radical species decrease. Not only does the precise location of the oxygen vacancy determine the maximal adsorption interaction, it also influences the adsorption strengths of the radicals at different surface configurations. This determines the probability of finding a certain adsorption configuration at the surface, which in turn influences the possible surface reactions. We find that C preferentially adsorbs far away from the oxygen vacancy, while CH2 and CH3 adsorb preferentially at the oxygen vacancy site. A fraction of CH partially adsorbs at the oxygen vacancy, and another fraction adsorbs further away from the vacancy.

Abstract: Cyrogenic etching of silicon is envisaged to enable better control over plasma processing in the microelectronics industry, albeit little is known about the fundamental differences compared to the room temperature process. We here present molecular dynamics simulations carried out to obtain sticking probabilities, thermal desorption rates, surface diffusion speeds, and sputter yields of F, F2, Si, SiF, SiF2, SiF3, SiF4, and the corresponding ions on Si(100) and on SiF13 surfaces, both at cryogenic and near room temperature. The different surface behavior during conventional etching and cryoetching is discussed. F2 is found to be relatively reactive compared to other species like SiF03. Thermal desorption occurs at a significantly lower rate under cryogenic conditions, which results in an accumulation of physisorbed species. Moreover, ion incorporation is often observed for ions with energies of 30400 eV, which results in a relatively low net sputter yield. The obtained results suggest that the actual etching of Si, under both cryogenic and near room temperature conditions, is based on the complete conversion of the Si surface to physisorbed SiF4, followed by subsequent sputtering of these molecules, instead of direct sputtering of the SiF03 surface.