Abstract: Bimetallic nanoparticles play a crucial role in various applications. A better understanding of their properties would facilitate these applications and possibly even enable chirality-specific growth of carbon nanotubes (CNTs). We here examine the stabilities of NiFe, NiGa, and FeGa nanoparticles and the effect of carbon dissolved in NiFe nanoparticles through density functional theory (DFT) calculations and Born Oppenheimer molecular dynamics (BOMD) simulations. We establish that nanoparticles with more Fe in the core and more Ga on the surface are more stable and compare these results with well-known properties such as surface energy and atom size. Furthermore, we find that the nanoparticles become more stable with increasing carbon content, both at 0 K and at 700 K. These results provide a basis for further research into the chirality-specific growth of CNT's.

Abstract: Bacteria can be inactivated through various physical and chemical means, and these have always been the focus of extensive research. To further improve the methodology for these ends, two types of plasma systems were investigated: nano-second pulsed plasma (NPP) as liquid discharge plasma and an Argon gas-feeding dielectric barrier discharge (Ar-DBD) as a form of surface plasma. To understand the sterilizing action of these two different plasma sources, we performed experiments with Staphylococcus aureus (S. aureus) bacteria (wild type) and multidrug resistant bacteria (Penicillum-resistant, Methicillin-resistant and Gentamicin-resistant). We observed that both plasma sources can inactivate both the wild type and multidrug-resistant bacteria to a good extent. Moreover, we observed a change in the surface morphology, gene expression and β-lactamase activity. Furthermore, we used X-ray photoelectron spectroscopy to investigate the variation in functional groups (C-H/C-C, C-OH and C=O) of the peptidoglycan (PG) resulting from exposure to plasma species. To obtain atomic scale insight in the plasma-cell interactions and support our experimental observations, we have performed molecular dynamics simulations to study the effects of plasma species, such as OH, H2O2, O, O3, as well as O2 and H2O, on the dissociation/formation of above mentioned functional groups in PG.

Abstract: Plasma catalysis is an emerging new technology for the electrification and downscaling of NH3 synthesis. Increasing attention is being paid to the optimization of plasma catalysis with respect to the plasma conditions, the catalyst material, and their mutual interaction. In this work we use microkinetic models to study how the total conversion process is impacted by the combination of different plasma conditions and transition metal catalysts. We study how plasma-generated radicals and vibrationally excited N2 (present in a dielectric barrier discharge plasma) interact with the catalyst and impact the NH3 turnover frequencies (TOFs). Both filamentary and uniform plasmas are studied, based on plasma chemistry models that provided plasma phase speciation and vibrational distribution functions. The Langmuir−Hinshelwood reaction rate coefficients (i.e., adsorption reactions and subsequent reactions among adsorbates) are determined using conventional scaling relations. An additional set of Eley−Rideal reactions (i.e., direct reactions of plasma radicals with adsorbates) was added and a sensitivity analysis on the assumed reaction rate coefficients was performed. We first show the impact of different vibrational distribution functions on the catalytic dissociation of N2 and subsequent production of NH3, and we gradually include more radical reactions, to illustrate the contribution of these species and their corresponding reaction pathways. Analysis over a large range of catalysts indicates that different transition metals (metals such as Rh, Ni, Pt, and Pd) optimize the NH3TOFs depending on the population of the vibrational levels of N2. At higher concentrations of plasma-generated radicals, the NH3 TOFs become less dependent on the catalyst material, due to radical adsorptions on the more noble catalysts and Eley−Rideal reactions on the less noble catalysts.

Abstract: In this paper, an overview is given of modeling investigations carried out in our research group for a better understanding of plasmas used for medical, environmental and nano applications. The focus is both on modeling the plasma chemistry and the plasma-surface interactions. The plasma chemistry provides the densities and fluxes of the important plasma species. This information can be used as input when modeling the plasma-surface interactions. The combination of plasma simulations and plasma – surface interaction simulations provides a more comprehensive understanding of the underlying processes for these applications.

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 this paper, we quantitatively investigate with atom probe tomography, the effect of temperature on the interfacial transition layer suboxide species due to the thermal oxidation of silicon. The chemistry at the interface was measured with atomic scale resolution, and the changes in chemistry and intermixing at the interface were identified on a nanometer scale. We find an increase of suboxide (SiOx) concentration relative to SiO2 and increased oxygen ingress with elevated temperatures. Our experimental findings are in agreement with reactive force field molecular dynamics simulations. This work demonstrates the direct comparison between atom probe derived chemical profiles and atomistic-scale simulations for transitional interfacial layer of suboxides as a function of temperature.

Abstract: Graphene oxide (GO) is one of the most promising candidates for next generation of atomically thin membranes. Nevertheless, one of the major issues for real world application of GO membranes is their undesirable swelling in an aqueous environment. Recently, we demonstrated that generation of H3O+ and OH- ions (e.g., with an external electric field) in the interlayer gallery could impart aqueous stability to the layered GO membranes (A. Gogoi, ACS Appl. Mater. Interfaces, 2022, 14, 34946). This, however, compromises the water flux through the membrane. In this study, we report on reducing the GO nanosheets as a solution to this issue. With the reduction of the GO nanosheets, the water flux through the layered GO membrane initially increases and then decreases again beyond a certain degree of reduction. Here, two key factors are at play. Firstly, the instability of the H-bond network between water molecules and the GO nanosheets, which increases the water flux. Secondly, the pore size reduction in the interlayer gallery of the membranes, which decreases the water flux. We also observe a significant improvement in the salt rejection of the membranes, due to the dissociation of water molecules in the interlayer gallery. In particular, for the case of 10% water dissociation, the water flux through the membranes can be enhanced without altering its selectivity. This is an encouraging observation as it breaks the traditional tradeoff between water flux and salt rejection of a membrane.

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: Density functional theory is widely used to gain insights into molecule−metal surface reaction networks, which is important for a better understanding of catalysis. However, it is well-known that generalized gradient approximation (GGA)
density functionals (DFs), most often used for the study of reaction networks, struggle to correctly describe both gas-phase molecules and metal surfaces. Also, GGA DFs typically underestimate reaction barriers due to an underestimation of the selfinteraction energy. Screened hybrid GGA DFs have been shown to reduce this problem but are currently intractable for wide usage. In this work, we use a more affordable meta-GGA (mGGA) DF in combination with a nonlocal correlation DF for the first time to study and gain new insights into a catalytically important surface
reaction network, namely, CO2 hydrogenation on Cu. We show that the mGGA DF used, namely, rMS-RPBEl-rVV10, outperforms typical GGA DFs by providing similar or better predictions for metals and molecules, as well as molecule−metal surface adsorption
and activation energies. Hence, it is a better choice for constructing molecule−metal surface reaction networks.

Abstract: We report the formation of nanobubbles on graphene with a radius of the order of 1 nm, using ultralow energy implantation of noble gas ions (He, Ne, Ar) into graphene grown on a Pt(111) surface. We show that the universal scaling of the aspect ratio, which has previously been established for larger bubbles, breaks down when the bubble radius approaches 1 nm, resulting in much larger aspect ratios. Moreover, we observe that the bubble stability and aspect ratio depend on the substrate onto which the graphene is grown (bubbles are stable for Pt but not for Cu) and trapped element. We interpret these dependencies in terms of the atomic compressibility of the noble gas as well as of the adhesion energies between graphene, the substrate, and trapped atoms.

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: 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.