Abstract: Nonthermal atmospheric pressure plasmas are gaining increasing attention for biomedical applications. However, very little fundamental information on the interaction mechanisms between the plasma species and biological cells is currently available. We investigate the interaction of important plasma species, such as OH, H2O2, O, O3, as well as O2 and H2O, with bacterial peptidoglycan by means of reactive molecular dynamics simulations, aiming for a better understanding of plasma disinfection. Our results show that OH, O, O3, and H2O2 can break structurally important bonds of peptidoglycan (i.e., CO, CN, or CC bonds), which consequently leads to the destruction of the bacterial cell wall. The mechanisms behind these breakups are, however, dependent on the impinging plasma species, and this also determines the effectiveness of the cell wall destruction.

Abstract: Device performance of Si|SiOx core-shell based nanowires critically depends on the exact control over the oxide thickness. Low-temperature plasma oxidation is a highly promising alternative to thermal oxidation allowing for improved control over the oxidation process, in particular for ultrasmall Si nanowires. We here elucidate the room temperature plasma oxidation mechanisms of ultrasmall Si nanowires using hybrid molecular dynamics / force-bias Monte Carlo simulations. We demonstrate how the oxidation and concurrent water formation mechanisms are a function of the oxidizing plasma species and we demonstrate how the resulting core-shell oxide thickness can be controlled through these species. A new mechanism of water formation is discussed in detail. The results provide a detailed atomic level explanation of the oxidation process of highly curved Si surfaces. These results point out a route toward plasma-based formation of ultrathin core-shell Si|SiOx nanowires at room temperature.

Abstract: Plasma catalysis is gaining increasing interest for various environmental applications, but the crucial question is whether plasma can be created inside catalyst pores and under which conditions. In practice, various catalytic support materials are used, with various dielectric constants. We investigate here the influence of the dielectric constant on the plasma properties inside catalyst pores and in the sheath in front of the pores, for various pore sizes. The calculations are performed by a two-dimensional fluid model for an atmospheric pressure dielectric barrier discharge in helium. The electron impact ionization rate, electron temperature, electron and ion density, as well as the potential distribution and surface charge density, are analyzed for a better understanding of the discharge behavior inside catalyst pores. The results indicate that, in a 100 μm pore, the electron impact ionization in the pore, which is characteristic for the plasma generation inside the pore, is greatly enhanced for dielectric constants below 300. Smaller pore sizes only yield enhanced ionization for smaller dielectric constants, i.e., up to εr = 200, 150, and 50 for pore sizes of 50, 30, and 10 μm. Thus, the most common catalyst supports, i.e., Al2O3 and SiO2, which have dielectric constants around εr = 8−11 and 4.2, respectively, should allow more easily that microdischarges can be formed inside catalyst pores, even for smaller pore sizes. On the other hand, ferroelectric materials with dielectric constants above 300 never seem to yield plasma enhancement inside catalyst pores, not even for 100 μm pore sizes. Furthermore, it is clear that the dielectric constant of the material has a large effect on the extent of plasma enhancement inside the catalyst pores, especially in the range between εr = 4 and εr = 200. The obtained results are explained in detail based on the surface charge density at the pore walls,
and the potential distribution and electron temperature inside and above the pores. The results obtained with this model are
important for plasma catalysis, as the production plasma species in catalyst pores might affect the catalyst properties, and thus
improve the applications of plasma catalysis.

Abstract: To gain insight into the nature of the adhesion mechanism between hydroxyapatite (HA) and rutile (rTiO(2)), the mutual affinity between their surfaces was systematically studied using density functional theory (DFT). We calculated both bulk and surface properties of HA and rTiO(2), and explored the interfacial bonding mechanism of amorphous HA (aHA) surface onto amorphous as well as stoichiometric and nonstoichiometric crystalline rTiO(2). Formation energies of bridging and subbridging oxygen vacancies considered in the rTiO(2)(110) surface were evaluated and compared with other theoretical and experimental results. The interfacial interaction was evaluated through the work of adhesion. For the aHA/rTiO(2)(110) interfaces, the work of adhesion is found to depend strongly on the chemical environment of the rTiO(2)(110) surface. Electronic analysis indicates that the charge transfer is very small in the case of interface formation between aHA and crystalline rTiO(2)(110). In contrast, significant charge transfer occurs between aHA and amorphous rTiO(2) (aTiO(2)) slabs during the formation of the interface. Charge density difference (CDD) analysis indicates that the dominant interactions in the interface have significant covalent character, and in particular the Ti-O and Ca-O bonds. Thus, the obtained results reveal that the aHA/aTiO(2) interface shows a more preferable interaction and is thermodynamically more stable than other interfaces. These results are particularly important for improving the long-term stability of HA-based implants.

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: Surface charging is an often overlooked factor in many plasma-surface interactions and in particular in plasma catalysis. In this study, we investigate the effect of excess electrons induced by a plasma on the adsorption properties of CO2 on titania-supported Cu-5 and Ni-5 clusters using spin-polarized and dispersion-corrected density functional theory calculations. The effect of excess electrons on the adsorption of Ni and Cu pentamers as well as on CO2 adsorption on a pristine anatase TiO2(101) slab is studied. Our results indicate that adding plasma-induced excess electrons to the system leads to further stabilization of the bent CO2 structure. Also, dissociation of CO2 on charged clusters is energetically more favorable than on neutral clusters. We hypothesize that surface charge is a plausible cause for the synergistic effects sometimes observed in plasma catalysis.

Abstract: Although the heterogeneity of pyrolyzed Fe???N???C materials is known and has been reported previously, the atomic structure of the active sites and their detailed reaction mechanisms are still unknown. Here, we identified two pyridinic Fe???N4-like centers with different local C coordinates, i.e., FeN4C8 and FeN4C10, and studied their electrocatalytic activity for the nitrogen reduction reaction (NRR) based on density functional theory (DFT) calculations. We also discovered the influence of the adsorption of NH2 as a functional ligand on catalyst performance on the NRR. We confirmed that the NRR selectivity of the studied catalysts is essentially governed either by the local C coordination or by the dynamic structure associated with the FeII/FeIII. Our investigations indicate that the proposed traditional pyridinic FeN4C10 has higher catalytic activity and selectivity for the NRR than the robust FeN4C8 catalyst, while it may have outstanding activity for promoting other (electro)catalytic reactions. <comment>Superscript/Subscript Available</comment

Abstract: Immobilization of single metal atoms on a solid host opens numerous possibilities for catalyst designs. If that host is a two-dimensional sheet, sheet curvature becomes a design parameter potentially complementary to host and metal composition. Here, we use a combination of density functional theory calculations and microkinetic modeling to compare the mechanisms and kinetics of formic acid decomposition and formation, chosen for their relevance as a potential hydrogen storage medium, over single Ru atoms anchored to pyridinic nitrogen in a planar graphene flake (RuN4-G) and curved carbon nanotube (RuN4-CNT). Activation barriers are lowered and the predicted turnover frequencies are increased over RuN4-CNT relative to RuN4-CNT. The results highlight the potential of curvature control as a means to achieve high performance and robust catalysts.

Abstract: Natural deep eutectic solvents (NADES) represent a green alternative to common organic solvents in the biochemical industry due to their benign behavior and tailorable properties, in particular as media for enzymatic reactions. However, to fully exploit their potential in enzymatic reactions, there is a need for a more fundamental understanding of how these neoteric solvents influence the course of these reac-tions. Thus, the aim of this study is to investigate the influence of NADES with various molar composi-tions on the stability and structure of enzymes, applying molecular dynamics simulations. This can help to better understand the effect of individual compounds of NADES, in addition to eutectic mixtures. More specifically, we simulate the behavior of Candida antarctica lipase B (CALB) enzyme in NADES com-posed of choline chloride with either urea, ethylene glycol or glycerol. Hereto, we monitor the NADES microstructure, the general stability of the enzyme and changes in the structure of its active sites and sur-face residues. Our simulations show that none of the studied NADES systems significantly disrupt the microstructure of the solvent or the stability of the CALB enzyme within the time scales of the simula-tions. The enzyme preserves its initial structure, size and intra-chain hydrogen bonds in all investigated compositions and, for the first time reported, also in NADES with increased hydrogen bond donating com-pound ratios. As the main novelty, our results indicate that, in addition to the composition, the molar ratio can be an additional variable to fine-tune the physicochemical properties of NADES without altering the enzyme characteristics. These findings could facilitate the development and application of task -tailored NADES media for biocatalytic processes. (c) 2022 Elsevier B.V. All rights reserved.

Abstract: The deposition of ultra-thin tungsten films and their related 2D chalcogen compounds on large area dielectric substrates by gas phase reactions is challenging. The lack of nucleation sites complicates the adsorption of W-related precursors and subsequent sulfurization usually requires high temperatures. We propose here a technique in which a thin solid amorphous silicon film is used as reductant for the gas phase precursor WF6 leading to the conversion to metallic W. The selectivity of the W conversion towards the underlying dielectric surfaces is demonstrated. The role of the Si surface preparation, the conversion temperature, and Si thickness on the formation process is investigated. Further, the in situ conversion of the metallic tungsten into thin stoichiometric WS2 is achieved by a cyclic approach based on WF6 and H2S pulses at the moderate temperature of 450 1C, which is much lower than usual oxide sulfurization processes.

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: Uniform-acceptance force-bias Monte Carlo (fbMC) methods have been shown to be a powerful technique to access longer timescales in atomistic simulations allowing, for example, phase transitions and growth. Recently, a new fbMC method, the time-stamped force-bias Monte Carlo (tfMC) method, was derived with inclusion of an estimated effective timescale; this timescale, however, does not seem able to explain some of the successes the method. In this contribution, we therefore explicitly quantify the effective timescale tfMC is able to access for a variety of systems, namely a simple single-particle, one-dimensional model system, the Lennard-Jones liquid, an adatom on the Cu(100) surface, a silicon crystal with point defects and a highly defected graphene sheet, in order to gain new insights into the mechanisms by which tfMC operates. It is found that considerable boosts, up to three orders of magnitude compared to molecular dynamics, can be achieved for solid state systems by lowering of the apparent activation barrier of occurring processes, while not requiring any system-specific input or modifications of the method. We furthermore address the pitfalls of using the method as a replacement or complement of molecular dynamics simulations, its ability to explicitly describe correct dynamics and reaction mechanisms, and the association of timescales to MC simulations in general.

Abstract: We calculate bubble nucleation rates in a Lennard-Jones fluid through explicit molecular dynamics simulations. Our approach-based on a recent free energy method (dubbed reweighted Jarzynski sampling), transition state theory, and a simple recrossing correction-allows us to probe a fairly wide range of rates in several superheated and cavitation regimes in a consistent manner. Rate predictions from this approach bridge disparate independent literature studies on the same model system. As such, we find that rate predictions based on classical nucleation theory, direct brute force molecular dynamics simulations, and seeding are consistent with our approach and one another. Published rates derived from forward flux sampling simulations are, however, found to be outliers. This study serves two purposes: First, we validate the reliability of common modeling techniques and extrapolation approaches on a paradigmatic problem in materials science and chemical physics. Second, we further test our highly generic recipe for rate calculations, and establish its applicability to nucleation processes.

Abstract: The growth mechanism and chirality formation of a single-walled carbon nanotube (SWNT) on a surface-bound nickel nanocluster are investigated by hybrid reactive molecular dynamics/force-biased Monte Carlo simulations. The validity of the interatomic potential used, the so-called ReaxFF potential, for simulating catalytic SWNT growth is demonstrated. The SWNT growth process was found to be in agreement with previous studies and observed to proceed through a number of distinct steps, viz., the dissolution of carbon in the metallic particle, the surface segregation of carbon with the formation of aggregated carbon clusters on the surface, the formation of graphitic islands that grow into SWNT caps, and finally continued growth of the SWNT. Moreover, it is clearly illustrated in the present study that during the growth process, the carbon network is continuously restructured by a metal-mediated process, thereby healing many topological defects. It is also found that a cap can nucleate and disappear again, which was not observed in previous simulations. Encapsulation of the nanoparticle is observed to be prevented by the carbon network migrating as a whole over the cluster surface. Finally, for the first time, the chirality of the growing SWNT cap is observed to change from (11,0) over (9,3) to (7,7). It is demonstrated that this change in chirality is due to the metal-mediated restructuring process.