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: 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: 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: 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: 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: The adsorption of C and CHx radicals on anatase (001) was studied using DFT within the generalized gradient approximation using the Perde-Burke-Ernzerhof (PBE) functional. We have studied the influence of oxygen vacancies in and at the surface on the adsorption properties of the radicals. For the oxygen vacancies in anatase (001), the most stable vacancy is located at the surface. For this vacancy, the maximal adsorption strength of C and CH decreases compared to the adsorption on the stoichiometric surface, but it increases for CH2 and CH3. If an oxygen vacancy is present in the first subsurface layer, the maximal adsorption strength increases for C, CH, CH2, and CH3. When the vacancy is present in the next subsurface layer, we find that only the CH3 adsorption is enhanced, while the maximal adsorption energies for the other radical species decrease. Not only does the precise location of the oxygen vacancy determine the maximal adsorption interaction, it also influences the adsorption strengths of the radicals at different surface configurations. This determines the probability of finding a certain adsorption configuration at the surface, which in turn influences the possible surface reactions. We find that C preferentially adsorbs far away from the oxygen vacancy, while CH2 and CH3 adsorb preferentially at the oxygen vacancy site. A fraction of CH partially adsorbs at the oxygen vacancy, and another fraction adsorbs further away from the vacancy.

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: New & beta;-stabilizedTi-based alloys are highly promising forbone implants, thanks in part to their low elasticity. The natureof this elasticity, however, is as yet unknown. We here present combinedfirst-principles DFT calculations and experiments on the microstructure,structural stability, mechanical characteristics, and electronic structureto elucidate this origin. Our results suggest that the studied & beta;Ti-35Nb-7Zr-5Ta wt % (TNZT) alloy manufacturedby the electron-beam powder bed fusion (E-PBF) method has homogeneousmechanical properties (H = 2.01 & PLUSMN; 0.22 GPa and E = 69.48 & PLUSMN; 0.03 GPa) along the building direction,which is dictated by the crystallographic texture and microstructuremorphologies. The analysis of the structural and electronic properties,as the main factors dominating the chemical bonding mechanism, indicatesthat TNZT has a mixture of strong metallic and weak covalent bonding.Our calculations demonstrate that the softening in the Cauchy pressure(C & PRIME; = 98.00 GPa) and elastic constant C ̅ ( 44 ) = 23.84 GPa is the originof the low elasticity of TNZT. Moreover, the nature of this softeningphenomenon can be related to the weakness of the second and thirdneighbor bonds in comparison with the first neighbor bonds in theTNZT. Thus, the obtained results indicate that a carefully designedTNZT alloy can be an excellent candidate for the manufacturing oforthopedic internal fixation devices. In addition, the current findingscan be used as guidance not only for predicting the mechanical propertiesbut also the nature of elastic characteristics of the newly developedalloys with yet unknown properties.

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: 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: 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: The recent advantages of the fabrication of artificial nanochannels enabled new research on the molecular transport, permeance, and selectivity of various gases and molecules. However, the physisorption/chemisorption of the unwanted molecules (usually hydrocarbons) inside nanochannels results in the alteration of the functionality of the nanochannels. We investigated contamination due to hydrocarbon molecules, nanochannels made of graphene, hexagonal boron nitride, BC2N, and molybdenum disulfide using molecular dynamics simulations. We found that for a certain size of nanochannel (i.e., h = 0.7 nm), as a result of the anomalous hydrophilic nature of nanochannels made of graphene, the hydrocarbons are fully adsorbed in the nanochannel, giving rise to full uptake. An increasing temperature plays an important role in unclogging, while pressure does not have a significant role. The results of our pioneering work contribute to a better understanding and highlight the important factors in alleviating the contamination and unclogging of nanochannels, which are in good agreement with the results of recent experiments.