Abstract: Amorphous platinum clusters supported on porous carbon have been envisaged for high-performance fuel cell electrodes. For this application, it is crucial to control the morphology of the Pt layer and the Ptsubstrate interaction to maximize activity and stability. We thus investigate the morphology evolution during Pt cluster growth on a porous carbon substrate employing atomic scale molecular dynamics simulations. The simulations are based on the Pt-C interaction potential using parameters derived from density functional theory and are found to yield a Pt cluster morphology similar to that observed in low loaded fuel cell electrodes prepared by plasma sputtering. Moreover, the simulations show amorphous Pt cluster growth in agreement with X-ray diffraction and transmission electron microscopy experiments on high performance low Pt content (10 μgPt cm−2) loaded fuel cell electrodes and provide a fundamental insight in the cluster growth mechanism.

Abstract: The resistance of a TiO2 anatase (001) surface to coke formation was studied in the context of dry reforming of methane using density functional theory (DFT) calculations. As carbon atoms act as precursors for coke formation, the resistance to coke formation can be measured by the carbon coverage of the surface. This is related to the stability of different CHx (x = 0−3) species and their rate of hydrogenation and dehydrogenation on the TiO2 surface. Therefore, we studied the reaction mechanisms and their corresponding rates as a function of the temperature for the dehydrogenation of the species on the surface. We found that the stabilities of C and CH are significantly lower than those of CH3 and CH2. The hydrogenation rates of the different species are significantly higher than the dehydrogenation rates in a temperature range of 300−1000 K. Furthermore, we found that dehydrogenation of CH3, CH2, and CH will only occur at appreciable rates starting from 600, 900, and 900 K, respectively. On the basis of these results, it is clear that the anatase (001) surface has a high coke resistance, and it is thus not likely that the surface will become poisoned by coke during dry reforming of methane. As the rate limiting step in dry reforming is the dissociative adsorption of CH4, we studied an alternative approach to thermal catalysis. We found that the temperature threshold for dry reforming is at least 700 K. This threshold temperature may be lowered by the use of plasma-catalysis, where the appreciable rates of adsorption of plasma-generated CHx radicals result in bypassing the rate limiting step of the reaction.

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: Recent experiments have revealed that the macroscopic Kelvin equation remains surprisingly accurate even for nanoscale capillaries. This phenomenon was so far explained by the oscillatory behavior of the solid−liquid interfacial free energy. We here demonstrate thermodynamic and capillarity inconsistencies with this explanation. After revising the Kelvin equation, we ascribe its validity at nanoscale confinement to the effect of disjoining pressure.
To substantiate our hypothesis, we employed molecular dynamics simulations to evaluate interfacial heat transfer and wetting properties. Our assessments unveil a breakdown in a previously established proportionality between the work of adhesion and the Kapitza conductance at capillary heights below 1.3 nm, where the dominance of the work of adhesion shifts primarily from energy to entropy. Alternatively, the peak density of the initial water layer can effectively probe the work of adhesion. Unlike under bulk conditions, high confinement renders the work of adhesion entropically unfavorable.

Abstract: Background: Strong electric fields are knownto affect cell membrane permeability,which can be applied for therapeutic purposes, e.g., in cancer therapy. A synergistic enhancement of this effect may be accomplished by the presence of reactive oxygen species (ROS), as generated in cold atmospheric plasmas. Little is known about the synergy between lipid oxidation by ROS and the electric field, nor on howthis affects the cell membrane permeability.
Method: We here conduct molecular dynamics simulations to elucidate the dynamics of the permeation process under the influence of combined lipid oxidation and electroporation. A phospholipid bilayer (PLB), consisting of di-oleoyl-phosphatidylcholine molecules covered with water layers, is used as a model system for the plasma membrane.
Results and conclusions:Weshow howoxidation of the lipids in the PLB leads to an increase of the permeability of the bilayer to ROS, although the permeation free energy barriers still remain relatively high. More importantly, oxidation of the lipids results in a drop of the electric field threshold needed for pore formation (i.e., electroporation) in the PLB. The created pores in the membrane facilitate the penetration of reactive plasma species deep into the cell interior, eventually causing oxidative damage.
General significance: This study is of particular interest for plasma medicine, as plasma generates both ROS and electric fields, but it is also of more general interest for applications where strong electric fields and ROS both come into play.

Abstract: In this study, we investigated the diffusion of H-atoms to the subsurface and their further diffusion into the bulk of a Ni(111) crystal by means of density functional theory calculations in the context of thermal and plasma-assisted catalysis. The H-atoms at the surface can originate from the dissociative adsorption of H2 or CH4 molecules, determining the surface H-coverage. When a threshold H-coverage is passed, corresponding to 1.00 ML for the crystalline Ni(111) surface, the surface-bound H-atoms start to diffuse to the subsurface. A similar threshold coverage is observed for the interstitial H-coverage. Once the interstitial sites are filled up with a coverage above 1.00 ML of H, dissolution of interstitial H-atoms to the layer below the interstitial sites will be initiated. Hence, by applying a high pressure or inducing a reactive plasma and high temperature, increasing the H-flux to the surface, a large amount of hydrogen can diffuse in a crystalline metal like Ni and can be absorbed. The formation of metal hydride may modify the entire reaction kinetics of the system. Equivalently, the H-atoms in the bulk can easily go back to the surface and release a large amount of heat. In a plasma process, H-atoms are formed in the plasma, and therefore the energy barrier for dissociative adsorption is dismissed, thus allowing achievement of the threshold coverage without applying a high pressure as in a thermal process. As a result, depending on the crystal plane and type of metal, a large number of H-atoms can be dissolved (absorbed) in the metal catalyst, explaining the high efficiency of plasma-assisted catalytic reactions. Here, the mechanism of H-dissolution is established as a chemical identifier for the investigation of the reaction kinetics of a chemical process.

Abstract: In this contribution, we present an overview of the various techniques for combining atomistic molecular dynamics with Monte Carlo simulations, mainly in the context of condensed matter systems, as well as a brief summary of the main accelerated dynamics techniques. Special attention is given to the force bias Monte Carlo technique and its combination with molecular dynamics, in view of promising recent developments, including a definable timescale. Various examples of the application of combined molecular dynamics / Monte Carlo simulations are given, in order to demonstrate the enhanced simulation efficiency with respect to either pure molecular dynamics or Monte Carlo.

Abstract: The unique plasma-specific features and physical phenomena in the organization of nanoscale soild-state systems in a broad range of elemental composition, structure, and dimensionality are critically reviewed. These effects lead to the possibility to localize and control energy and matter at nanoscales and to produce self-organized nano-solids with highly unusual and superior properties. A unifying conceptual framework based on the control of production, transport, and self-organization of precursor species is introduced and a variety of plasma-specific non-equilibrium and kinetics-driven phenomena across the many temporal and spatial scales is explained. When the plasma is localized to micrometer and nanometer dimensions, new emergent phenomena arise. The examples range from semiconducting quantum dots and nanowires, chirality control of single-walled carbon nanotubes, ultra-fine manipulation of graphenes, nano-diamond, and organic matter to nano-plasma effects and nano-plasmas of different states of matter.

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: Plasma-catalytic CO, hydrogenation for methanol production is gaining increasing interest, but our understanding of its reaction mechanism remains primitive. We present a combined experimental/computational study on plasma-catalytic CO, hydrogenation to CH3OH over a size-selected Cu/gamma-Al2O3 catalyst. Our experiments demonstrate a synergistic effect between the Cu/gamma-Al2O3 catalyst and the CO2/H-2 plasma, achieving a CO2 conversion of 10% at 4 wt % Cu loading and a CH3OH selectivity near 50% further rising to 65% with H2O addition (for a H2O/CO2 ratio of 1). Furthermore, the energy consumption for CH3OH production was more than 20 times lower than with plasma only. We carried out density functional theory calculations over a Cu-13/gamma-Al2O3 model, which reveal that the interfacial sites of the Cu-13 cluster and gamma-Al2O3 support show a bifunctional effect: they not only activate the CO2 molecules but also strongly adsorb key intermediates to promote their hydrogenation further. Reactive plasma species can regulate the catalyst surface reactions via the Eley-Rideal (E-R) mechanism, which accelerates the hydrogenation process and promotes the generation of the key intermediates. H2O can promote the CH3OH desorption by competitive adsorption over the Cu-13/gamma-Al2O3 surface. This study provides new insights into CO2 hydrogenation through plasma catalysis, and it provides inspiration for the conversion of some other small molecules (CH4, N-2, CO, etc.) by plasma catalysis using supported-metal clusters.

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: 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: Plasma catalysis holds great promise for environmental applications, provided that the process viability can be maximized in terms of energy efficiency and product selectivity. This requires a fundamental understanding of the various processes taking place and especially the mutual interactions between plasma and catalyst. In this review, we therefore first examine the various effects of the plasma on the catalyst and of the catalyst on the plasma that have been described in the literature. Most of these studies are purely experimental. The urgently needed fundamental understanding of the mechanisms underpinning plasma catalysis, however, may also be obtained through modelling and simulation. Therefore, we also provide here an overview of the modelling efforts that have been developed already, on both the atomistic and the macroscale, and we identify the data that can be obtained with these models to illustrate how modelling and simulation may contribute to this field. Last but not least, we also identify future modelling opportunities to obtain a more complete understanding of the various underlying plasma catalytic effects, which is needed to provide a comprehensive picture of plasma catalysis.

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.