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: 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: 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: The interaction between cold atmospheric pressure plasma and liquids is receiving increasing attention for various applications. In particular, the use of plasma-treated liquids (PTL) for biomedical applications is of growing importance, in particular for sterilization and cancer treatment. However, insight into the
underlying mechanisms of plasma–liquid interactions is still scarce. Here, we present a 2D fluid dynamics model for the interaction between a plasma jet and liquid water. Our results indicate that the formed reactive species originate from either the gas phase (with further solvation) or are formed at the liquid interface. A clear increase in the aqueous density of H2O2, HNO2/NO2- and NO3-
is observed as a function of time, while the densities of O3, HO2/O2- and ONOOH/ONOO- are found to quickly reach a maximum due to chemical reactions in solution. The trends observed in our model correlate well with experimental observations from the literature.

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: 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: 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: CO2 conversion into value-added chemicals and fuels is considered as one of the great challenges of the 21st century. Due to the limitations of the traditional thermal approaches, several novel technologies are being developed. One promising approach in this field, which has received little attention to date, is plasma
technology. Its advantages include mild operating conditions, easy upscaling, and gas activation by energetic electrons instead of heat. This allows thermodynamically difficult reactions, such as CO2 splitting and the dry reformation of methane, to occur with reasonable energy cost. In this review, after exploring the traditional thermal approaches, we have provided a brief overview of the fierce competition between various novel approaches in a quest to find the most effective and efficient CO2 conversion technology. This is needed to critically assess whether plasma technology can be successful in an already crowded arena. The following questions need to be answered in this regard: are there key advantages to using plasma technology over other novel approaches, and if so, what is the flip side to the use of this technology? Can plasma technology be successful on its own, or can synergies be achieved by combining it with other technologies? To answer
these specific questions and to evaluate the potentials and limitations of plasma technology in general, this review presents the current state-of-the-art and a critical assessment of plasma-based CO2 conversion, as well as the future challenges for its practical implementation.

Abstract: Dry reforming of methane has gained significant interest over the years. A novel reforming technique with great potential is plasma technology. One of its drawbacks, however, is energy consumption. Therefore, we performed an extensive computational study, supported by experiments, aiming to identify the influence of the operating parameters (gas mixture, power, residence time and frequency) of a dielectric barrier discharge plasma on the conversion and energy efficiency, and to investigate which of these parameters lead to the most promising results and whether these are eventually sufficient for industrial implementation. The best results, in terms of both energy efficiency and conversion, are obtained at a specific energy input (SEI) of 100 J cm−3, a 1090 CH4CO2 ratio, 10 Hz, a residence time of 1 ms, resulting in a total conversion of 84% and an energy efficiency of 8.5%. In general, increasing the CO2 content in the gas mixture leads to a higher conversion and energy efficiency. The SEI couples the effect of the power and residence time, and increasing the SEI always results in a higher conversion, but somewhat lower energy efficiencies. The effect of the frequency is more complicated: we observed that the product of frequency (f) and residence time (τ), being a measure for the total number of micro-discharge filaments which the gas molecules experience when passing through the reactor, was critical. For most cases, a higher number of filaments yields higher values for conversion and energy efficiency. To benchmark our model predictions, we also give an overview of measured conversions and energy efficiencies reported in the literature, to indicate the potential for improvement compared to the state-of-the art. Finally, we identify the limitations as well as the benefits and future possibilities of plasma technology.

Abstract: This paper discusses our recent results on plasma-based CO2 conversion, obtained by a combination of experiments and modeling, for a dielectric barrier discharge (DBD), a microwave plasma and a packed bed DBD reactor. The results illustrate that plasma technology is quite promising for CO2 conversion, but more research is needed to better understand the underlying mechanisms and to further improve the capabilities.

Abstract: Gold nanocluster properties exhibit unique size-dependence. In this contribution, we employ reactive molecular dynamics simulations to calculate the size- and temperature-dependent surface energies, strain energies and atomic displacements for icosahedral, cuboctahedral, truncated octahedral and decahedral Au-nanoclusters. The calculations demonstrate that the surface energy decreases with increasing cluster size at 0 K but increases with size at higher temperatures. The calculated melting curves as a function of cluster size demonstrate the Gibbs-Thomson effect. Atomic displacements and strain are found to strongly depend on the cluster size and both are found to increase with increasing cluster size. These results are of importance for understanding the size-and temperature-dependent surface processes on gold nanoclusters.

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: In recent years, the ability of cold atmospheric pressure plasmas (CAPS) to selectively induce cell death in cancer cells has been widely established. This selectivity has been assigned to the reactive oxygen and nitrogen species (RONS) created in CAPs. To provide new insights in the search for an explanation
for the observed selectivity, we calculate the transfer free energy of multiple ROS across membranes containing a varying amount of cholesterol. The cholesterol fraction is investigated as a selectivity parameter because membranes of cancer cells are known to contain lower fractions of cholesterol compared to healthy cells. We find that cholesterol has a significant effect on the permeation of
reactive species across a membrane. Indeed, depending on the specific reactive species, an increasing cholesterol fraction can lead to (i) an increase of the transfer free energy barrier height and width, (ii) the formation of a local free energy minimum in the center of the membrane and (iii) the creation of extra free energy barriers due to the bulky sterol rings. In the context of plasma oncology, these observations suggest that the increased ingress of RONS in cancer cells can be explained by the decreased cholesterol fraction of their cell membrane.

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.