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: 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: Carbyne is a novel material of current interest in nanotechnology. As is typically the case for nanomaterials, the growth process determines the resulting properties. While endohedral carbyne has been successfully synthesized, its catalyst and feedstock-dependent growth mechanism is still elusive. We here study the nucleation and growth mechanism of different carbon chains in a Ni-containing double walled carbon nanotube using classical molecular dynamics simulations and first-principles calculations. We find that the understanding the competitive role of the metal catalyst and the hydrocarbon is important to control the growth of 1-dimensional carbon chains, including Ni or H-terminated carbyne. Also, we find that the electronic property of the Ni-terminated carbyne can be tuned by steering the H concentration along the chain. These results suggest catalyst-containing carbon nanotubes as a possible synthesis route for carbyne formation.

Abstract: The free energy surface (FES) for carbon segregation from nickel nanoparticles is obtained from advanced molecular dynamics simulations. A suitable reaction coordinate is developed that can distinguish dissolved carbon atoms from segregated dimers, chains and junctions on the nanoparticle surface. Because of the typically long segregation time scale (up to ms), metadynamics simulations along the developed reaction coordinate are used to construct FES over a wide range of temperatures and carbon concentrations. The FES revealed the relative stability of different stages in the segregation process, and free energy barriers and rates of the individual steps could then be calculated and decomposed into enthalpic and entropic contributions. As the carbon concentration in the nickel nanoparticle increases, segregated carbon becomes more stable in terms of both enthalpy and entropy. The activation free energy of the reaction also decreases with the increase of carbon concentration, which can be mainly attributed to entropic effects. These insights and the methodology developed to obtain them improve our understanding of carbon segregation process across materials science in general, and the nucleation and growth of carbon nanotube in particular.

Abstract: Using self-consistent KohnSham density-functional theory molecular dynamics simulations, we demonstrate the theoretical possibility to synthesize NiC60, the incarfullerene Ni@C60 and the heterofullerene C59Ni in an ion implantation setup. The corresponding formation mechanisms of all three complexes are elucidated as a function of the ion implantation energy and impact location, suggesting possible routes for selectively synthesizing these complexes.

Abstract: Ion irradiation of carbon nanotubes very often leads to defect formation. However, we have recently shown that Ar ion irradiation in a limited energy window of 1025 eV may enhance the initial cap nucleation process, when the carbon network is in contact with the metal nanocatalyst. Here, we employ reactive molecular dynamics simulations to demonstrate that ion irradiation in a higher energy window of 1035 eV may also heal network defects after the nucleation stage through a non-metal-mediated mechanism, when the carbon network is no longer in contact with the metal nanocatalyst. The results demonstrate the possibility of beneficially utilizing ions in e.g. plasma-enhanced chemical vapour deposition of carbon nanotubes.

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 safe and efficient use of hydrogen energy, which is in high demand worldwide today, requires efficient hydrogen storage. Despite significant advances in hydrogen storage using carbon-based nanomaterials, including carbon nanotubes (CNTs), efforts to substantially increase the storage capacity remain less effective. In this work, we demonstrate the effect of endohedral transition metal atoms on the hydrogen storage capacity of CNTs using reactive molecular dynamics simulations. We find that an increase in the volume fraction of endohedral nickel atoms leads to an increase in the concentration of physisorbed hydrogen molecules around single-walled CNTs (SWNTs) by approximately 1.6 times compared to pure SWNTs. The obtained results provide insight into the underlying mechanisms of how endohedral transition metal atoms enhance the hydrogen storage ability of SWNTs under nearly ambient conditions.

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: The application of coreshell SiSiO2 nanowires as nanoelectronic devices strongly depends on their structure, which is difficult to tune precisely. In this work, we investigate the formation of the coreshell nanowires at the atomic scale, by reactive molecular dynamics simulations. The occurrence of two temperature-dependent oxidation mechanisms of ultra-small diameter Si-NWs is demonstrated. We found that control over the Si-core radius and the SiOx (x ≤ 2) oxide shell is possible by tuning the growth temperature and the initial Si-NW diameter. Two different structures were obtained, i.e., ultrathin SiO2 silica nanowires at high temperature and Si core|ultrathin SiO2 silica nanowires at low temperature. The transition temperature is found to linearly decrease with the nanowire curvature. Finally, the interfacial stress is found to be responsible for self-limiting oxidation, depending on both the initial Si-NW radius and the oxide growth temperature. These novel insights allow us to gain control over the exact morphology and structure of the wires, as is needed for their application in nanoelectronics.

Abstract: Using reactive molecular dynamics simulations, the melting behavior of nickelcarbon 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 nickelcarbon 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: Structural defects inevitably appear during the nucleation event that determines the structure and properties of single-walled carbon nanotubes. By combining ion bombardment experiments with atomistic simulations we reveal that ion bombardment in a suitable energy range allows these defects to be healed resulting in an enhanced nucleation of the carbon nanotube cap. The enhanced growth of the nanotube cap is explained by a nonthermal ion-induced graphene network restructuring mechanism.

Abstract: Natural deep eutectic solvents (NADES) are mixtures of naturally derived compounds with a significantly decreased melting point due to the specific interactions among the constituents. NADES have benign properties (low volatility, flammability, toxicity, cost) and tailorable physicochemical properties (by altering the type and molar ratio of constituents), hence they are often considered as a green alternative to common organic solvents. Modeling the relation between their composition and properties is crucial though, both for understanding and predicting their behavior. Several efforts were done to this end, yet this review aims at structuring the present knowledge as an outline for future research. First, we reviewed the key properties of NADES and relate them to their structure based on the available experimental data. Second, we reviewed available modeling methods applicable to NADES. At the molecular level, density functional theory and molecular dynamics allow interpreting density differences and vibrational spectra, and computation of interaction energies. Additionally, properties at the level of the bulk media can be explained and predicted by semi-empirical methods based on ab initio methods (COSMO-RS) and equation of state models (PC-SAFT). Finally, methods based on large datasets are discussed; models based on group contribution methods and machine learning. A combination of bulk media and dataset modeling allows qualitative prediction and interpretation of phase equilibria properties on the one hand, and quantitative prediction of melting point, density, viscosity, surface tension and refractive indices on the other hand. In our view, multiscale modeling, combining the molecular and macroscale methods, will strongly enhance the predictability of NADES properties and their interaction with solutes, yielding truly tailorable solvents to accommodate (bio)chemical reactions.

Abstract: Atomistic simulations can in principle provide an unbiased description of all mechanisms, intermediates, and products of complex chemical processes. However, due to the severe time scale limitation of conventional simulation techniques, unrealistically high simulation temperatures are usually applied, which are a poor approximation of most practically relevant low-temperature applications. In this work, we demonstrate the direct observation at the atomic scale of the pyrolysis and oxidation of n-dodecane at temperatures as low as 700 K through the use of a novel simulation technique, collective variable-driven hyperdynamics (CVHD). A simulated timescale of up to 39 seconds is reached. Product compositions and dominant mechanisms are found to be strongly temperature-dependent, and are consistent with experiments and kinetic models. These simulations provide a first atomic-level look at the full dynamics of the complicated fuel combustion process at industrially relevant temperatures and time scales, unattainable by conventional molecular dynamics simulations.

Abstract: tWe investigate microdischarge formation inside catalyst pores by a two-dimensional fluid model forvarious pore sizes in the m-range and for various applied voltages. Indeed, this is a poorly understoodphenomenon in plasma catalysis. The calculations are performed for a dielectric barrier discharge inhelium, at atmospheric pressure. The electron and ion densities, electron temperature, electric field andpotential, as well as the electron impact ionization and excitation rate and the densities of excited plasmaspecies, are examined for a better understanding of the characteristics of the plasma inside a pore. Theresults indicate that the pore size and the applied voltage are critical parameters for the formation of amicrodischarge inside a pore. At an applied voltage of 20 kV, our calculations reveal that the ionizationmainly takes place inside the pore, and the electron density shows a significant increase near and inthe pore for pore sizes larger than 200m, whereas the effect of the pore on the total ion density isevident even for 10m pores. When the pore size is fixed at 30m, the presence of the pore has nosignificant influence on the plasma properties at an applied voltage of 2 kV. Upon increasing the voltage,the ionization process is enhanced due to the strong electric field and high electron temperature, andthe ion density shows a remarkable increase near and in the pore for voltages above 10 kV. These resultsindicate that the plasma species can be formed inside pores of structured catalysts (in the m range),and they may interact with the catalyst surface, and affect the plasma catalytic process.