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: Controlling and steering the growth of single walled carbon nanotubes is often believed to require controlling of the nucleation stage. Yet, little is known about the microscopic mechanisms governing the nucleation from hydrocarbon molecules. Specifically, we address here the dehydrogenation of hydrocarbon molecules and the formation of all-carbon graphitic islands on metallic nanoclusters from hydrocarbon molecules under conditions typical for carbon nanotube growth. Employing reactive molecular dynamics simulations, we demonstrate for the first time that the formation of a graphitic network occurs through the intermediate formation of vertically oriented, not fully dehydrogenated graphitic islands. Upon dehydrogenation of these vertical graphenes, the islands curve over the surface, thereby forming a carbon network covering the nanoparticle. The results indicate that controlling the extent of dehydrogenation offers an additional parameter to control the nucleation of carbon nanotubes.

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: 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: 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: 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: Hydrogen–graphite interactions are relevant to a wide variety of applications, ranging from astrophysics to fusion devices and nano-electronics. In order to shed light on these interactions, atomistic simulation using Molecular Dynamics (MD) has been shown to be an invaluable tool. It suffers, however, from severe timescale
limitations. In this work we apply the recently developed Collective Variable-Driven Hyperdynamics (CVHD) method to hydrogen etching of graphite for varying inter-impact times up to a realistic value of 1 ms, which corresponds to a flux of
1020 m2 s1. The results show that the erosion yield, hydrogen surface coverage and species distribution are significantly affected by the time between impacts. This can be explained by the higher probability of C–C bond breaking due to the prolonged exposure to thermal stress and the subsequent transition from ion- to thermal-induced etching. This latter regime of thermal-induced etching – chemical erosion – is here accessed for the first time using atomistic simulations. In conclusion, this study demonstrates that accounting for long time-scales significantly affects ion bombardment simulations and should not be neglected in a wide range of conditions, in contrast to what is typically assumed.

Abstract: Methane reforming by plasma catalysis is a complex process that is far from understood. It requires a multidisciplinary approach which ideally takes into account all effects from the plasma on the catalyst, and vice versa. In this contribution, we focus on the interactions of CHx (x = {1,2,3}) radicals that are created in the plasma with several nickel catalyst surfaces. To this end, we perform reactive molecular dynamics simulations, based on the ReaxFF potential, in a wide temperature range of 4001600 K. First, we focus on the H2 formation as a function of temperature and surface structure. We observe that substantial H2 formation is obtained at 1400 K and above, while the role of the surface structure seems limited. Indeed, in the initial stage, the type of nickel surface influences the CH bond breaking efficiency of adsorbed radicals; however, the continuous carbon diffusion into the surface gradually diminishes the surface crystallinity and therefore reduces the effect of surface structure on the H2 formation probability. Furthermore, we have also investigated to what extent the species adsorbed on the catalyst surface can participate in surface reactions more in general, for the various surface structures and as a function of temperature. These results are part of the ongoing research on the methane reforming by plasma catalysis, a highly interesting yet complex alternative to conventional reforming processes.

Abstract: tWe investigated the plasma-assisted catalytic reactions for the production of value-added chemicalsfrom Ni-catalyzed plasma dry reforming of methane by means of density functional theory (DFT). Weinspected many activation barriers, from the early stage of adsorption of the major chemical fragmentsderived fromCH4andCO2molecules up to the formation of value-added chemicals at the surface, focusingon the formation of methanol, as well as the hydrogenation of C1and C2hydrocarbon fragments. Theactivation barrier calculations show that the presence of surface-bound H atoms and in some cases alsoremaining chemical fragments at the surface facilitates the formation of products. This implies that thehydrogenation of a chemical fragment on the hydrogenated crystalline surface is energetically favouredcompared to the simple hydrogenation of the chemical fragment at the bare Ni(111) surface. Indeed, thepresence of hydrogen modifies the electronic structure of the surface and the course of the reactions.We therefore conclude that surface-bound H atoms, and to some extent also the remaining chemicalfragments at the crystalline surface, induce the following effects: they facilitate associative desorption ofmethanol and ethane by increasing the rate of H-transfer to the adsorbed fragments while they impedehydrogenation of ethylene to ethane, thus promoting again the desorption of ethylene. Overall, they thusfacilitate the catalytic conversion of the formed fragments from CH4and CO2, into value-added chemicals.Finally, we believe that the retention of methane fragments, especially CH3, in the presence of surface-boundHatoms (as observed here for Ni) can be regarded as an identifier for the proper choice of a catalystfor the production of value-added chemicals.

Abstract: Over the past decade, graphene oxide (GO) has emerged as a promising membrane material with superior separation performance and intriguing mechanical/chemical stability. However, its practical implementation remains very challenging primarily because of its undesirable swelling in an aqueous environment. Here, we demonstrated that dissociation of water molecules into H3O+ and OH- ions inside the interlayer gallery of a layered GO membrane can strongly affect its stability and performance. We reveal that H3O+ and OH- ions form clusters inside the GO laminates that impede the permeance of water and salt ions through the membrane. Dynamics of those clusters is sensitive to an external ac electric field, which can be used to tailor the membrane performance. The presence of H3O+ and OH- ions also leads to increased stability of the hydrogen bond (H-bond) network among the water molecules and the GO layers, which further reduces water permeance through the membrane, while crucially imparting stability to the layered GO membrane against undesirable swelling. KEYWORDS: layered graphene-oxide membrane, aqueous stability, H3O+ and OH- ions, external electric field, molecular dynamics

Abstract: Industrial production of HCN from NH3 and CH4 not only uses precious Pt or Pt−Rh catalysts but also requires extremely high temperatures (∼1600 K). From an energetic, operational, and safety perspective, a drastic decrease in temperature is highly desirable. Here, we report ammonia reforming of methane for the production of HCN and H2 at 673 K by the combination of CH4/NH3 plasma and a supported Cu/silicalite-1 catalyst. 30% CH4 conversion has been achieved with 79% HCN selectivity. Catalyst characterization and plasma diagnostics reveal that the excellent reaction performance is attributed to metallic Cu active sites. In addition, we propose a possible reaction pathway, viz. E-R reactions with N, NH, NH2, and CH radicals produced in the plasma, for the production of HCN, based on density functional theory calculations. Importantly, the Cu/silicalite-1 catalyst costs less than 5% of the commercial Pt mesh catalyst.

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: Although recent studies have advanced the understanding of pyrolyzed
Fe−N−C materials as oxygen reduction reaction (ORR) catalysts, the atomic and
electronic structures of the active sites and their detailed reaction mechanisms still remain unknown. Here, based on first-principles density functional theory (DFT) computations, we discuss the electronic structures of three FeN4 catalytic centers with different local topologies of the surrounding C atoms with a focus on unraveling the mechanism of their ORR activity in acidic electrolytes. Our study brings back a forgotten, synthesized pyridinic Fe−N coordinate to the community’s attention, demonstrating that this catalyst can exhibit excellent activity for promoting direct four-electron ORR through the addition of a fifth ligand such as −NH2, −OH, and −SO4. We also identify sites with good stability properties through the combined use of our DFT calculations and Mössbauer spectroscopy data.

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: Metal-catalyzed growth mechanisms of carbon nanotubes (CNTs) were studied by hybrid molecular dynamics−Monte Carlo simulations using a recently developed ReaxFF reactive force field. Using this novel approach, including relaxation effects, a CNT with definable chirality is obtained, and a step-by-step atomistic description of the nucleation process is presented. Both root and tip growth mechanisms are observed. The importance of the relaxation of the network is highlighted by the observed healing of defects.