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: While it is well-known that the plasma-enhanced catalytic chemical vapor deposition (PECVD) of carbon nanotubes (CNTs) offers a number of advantages over thermal CVD, the influence of the various individual contributing factors is not well understood. Especially the role of ions is unclear, since ions in plasmas are generally associated with sputtering rather than with growing a material. Even so, various studies have demonstrated the beneficial effects of ion bombardment during the growth of CNTs. This review looks at the role of the ions in plasma-enhanced CNT growth as deduced from both experimental and simulation studies. Specific attention is paid to the beneficial effects of ion bombardment. Based on the available literature, it can be concluded that ions can be either beneficial or detrimental for carbon nanotube growth, depending on the exact conditions and the control over the growth process.

Abstract: We report on multi-level atomistic simulations for the interaction of reactive oxygen species (ROS) with the head groups of the phospholipid bilayer, and the subsequent effect of head group and lipid tail oxidation on the structural and dynamic properties of the cell membrane. Our simulations are validated by experiments using a cold atmospheric plasma as external ROS source. We found that plasma treatment leads to a slight initial rise in membrane rigidity, followed by a strong and persistent increase in fluidity, indicating a drop in lipid order. The latter is also revealed by our simulations. This study is important for cancer treatment by therapies producing (extracellular) ROS, such as plasma treatment. These ROS will interact with the cell membrane, first oxidizing the head groups, followed by the lipid tails. A drop in lipid order might allow them to penetrate into the cell interior (e.g., through pores created due to oxidation of the lipid tails) and cause intracellular oxidative damage, eventually leading to cell death. This work in general elucidates the underlying mechanisms of ROS interaction with the cell membrane at the atomic level.

Abstract: We present an approach to incorporate the effect of vibrational nonequilibrium in molecular dynamics (MD) simulations. A perturbed canonical ensemble, in which selected modes are excited to higher temperature while all others remain equilibrated at low temperature, is simulated by applying a specifically tailored bias potential. Our method can be readily applied to any (classical or quantum mechanical) MD setup at virtually no additional computational cost and allows the study of reactions of vibrationally excited molecules in nonequilibrium environments such as plasmas. In combination with enhanced sampling methods, the vibrational efficacy and mode selectivity of vibrationally stimulated reactions can then be quantified in terms of chemically relevant observables, such as reaction rates and apparent free energy barriers. We first validate our method for the prototypical hydrogen exchange reaction and then show how it can capture the effect of vibrational excitation on a symmetric SN2 reaction and radical addition on CO2.

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: We investigated the deposition of MoS2 multilayers on large area substrates. The pre-deposition of metal or metal oxide with subsequent sulfurization is a promising technique to achieve layered films. We distinguish a different reaction behavior in metal oxide and metallic films and investigate the effect of the temperature, the H2S/H-2 gas mixture composition, and the role of the underlying substrate on the material quality. The results of the experiments suggest a MoS2 growth mechanism consisting of two subsequent process steps. At first, the reaction of the sulfur precursor with the metal or metal oxide occurs, requiring higher temperatures in the case of metallic film compared to metal oxide. At this stage, the basal planes assemble towards the diffusion direction of the reaction educts and products. After the sulfurization reaction, the material recrystallizes and the basal planes rearrange parallel to the substrate to minimize the surface energy. Therefore, substrates with low roughness show basal plane assembly parallel to the substrate. These results indicate that the substrate character has a significant impact on the assembly of low dimensional MoS2 films.

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: van der Waals (vdW) forces play an important role in the adsorption of molecules on the surface of solids. However, the choice of the most suitable vdW functional for different systems is an essential problem which must be addressed for different systems. The lack of a systematic study on the performance of the vdW functionals in the bulk and adsorption properties of metal-oxides motivated us to examine different vdW approaches and compute the bulk and molecular adsorption properties of alpha-Cr2O3, alpha-Fe2O3, and alpha-Al2O3. For the bulk properties, we compared our results for the heat of formation, cohesive energy, lattice parameters and bond distances between the different vdW functionals and available experimental data. Next we studied the adsorption of benzene and CH3 molecules on top of different oxide surfaces. We employed different approximations to exchange and correlation within DFT, namely, the Perdew-Burke-Ernzerhof (PBE) GGA, (PBE)+U, and vdW density functionals [ DFT(vdW-DF/DF2/optPBE/optB86b/optB88)+U] as well as DFT-D2/D3(+U) methods of Grimme for the bulk calculations and optB86b-vdW(+U) and DFT-D2(+U) for the adsorption energy calculations. Our results highlight the importance of vdW interactions not only in the adsorption of molecules, but importantly also for the bulk properties. Although the vdW contribution in the adsorption of CH3 (as a chemisorption interaction) is less important compared to the adsorption of benzene (as a physisorption interaction), this contribution is not negligible. Also adsorption of benzene on ferryl/chromyl terminated surfaces shows an important chemisorption contribution in which the vdW interactions become less significant.

Abstract: Using reactive molecular dynamics simulations, we have investigated the effect of single-impact, low-energy (thermal-100 eV) bombardment of a Si(100){2 × 1} surface by atomic and molecular oxygen. Penetration probability distributions, as well as defect formation distributions, are presented as a function of the impact energy for both species. It is found that at low impact energy, defects are created chemically due to the chemisorption process in the top layers of the surface, while at high impact energy, additional defects are created by a knock-on displacement of Si. These results are of particular importance for understanding device performances of silica-based microelectronic and photovoltaic devices.

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: Using reactive molecular dynamics simulations based on the ReaxFF potential, we studied the growth mechanism of ultrathin silica (SiO2) layers during hyperthermal oxidation as a function of temperature in the range 100-1300 K. Oxidation of Si(100){2 x 1} surfaces by both atomic and molecular oxygen was investigated for hyperthermal impact energies in the range of 1 to 5 eV. Two different growth mechanisms are found, corresponding to a low temperature oxidation and a high temperature one. The transition temperature between these mechanisms is estimated to be about 700 K. Also, the initial step of the Si oxidation process is analyzed in detail. Where possible, we validated our results with experimental and ab initio data, and good agreement was obtained. 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.

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: Uniform acceptance force biased Monte Carlo (UFMC) simulations have previously been shown to be a powerful tool to simulate atomic scale processes, enabling one to follow the dynamical path during the simulation. In this contribution, we present a simple proof to demonstrate that this uniform acceptance still complies with the condition of detailed balance, on the condition that the characteristic parameter lambda = 1/2 and that the maximum allowed step size is chosen to be sufficiently small. Furthermore, the relation to Metropolis Monte Carlo (MMC) is also established, and it is shown that UFMC reduces to MMC by choosing the characteristic parameter lambda = 0 [Rao, M. et al. Mol. Phys. 1979, 37, 1773]. Finally, a simple example compares the UFMC and MMC methods.

Abstract: Understanding the fundamental dynamics in carbon nanotube (CNT) catalysts is of primary importance to understand CNT nucleation. This Letter reports on calculated bond switching (BS) rates in pure and carbon containing nickel nanoclusters. The rates are analyzed in terms of their temperature dependent spatial distribution and the mobility of the cluster atoms. The BS mechanism is found to change from vibrational to diffusional at around 900 K, with a corresponding strong increase in activation energy. Furthermore, the BS activation energy is observed to decrease as the carbon content in the cluster increases, resulting in an effective liquification of the cluster.

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.