Abstract: Mannosylerythritol lipids (MELs) are a promising group of biosurfactants due to their high fermentation yield, selfassembly and biological activity. During fermentation by Pseudozyma aphidis, a mixture of MELs with different levels of acylation is formed, of which the fully deacetylated form is the most valuable. In order to reduce the environmental impact of deacetylation, an enzymatic process using natural deep eutectic solvents (NADES) has been developed. We tested the deacetylation of a purified MELs mixture with immobilized Candida antarctica lipase B enzyme and 2-ethylhexanol as co-substrate in 140 h reactions with different NADES. We identified hydrophobic NADES systems with similar yields and kinetics as in pure 2-ethylhexanol solvent. Our results indicate that deacetylation of MELs mixtures in NADES as a solvent is possible with yields comparable to pure co-substrate and that hydrophobic NADES without carboxylic acid compounds facilitate the reaction to the greatest extent.

Abstract: There is currently a growing interest in the realisation and optimization of hybrid plasma/catalyst systems for a multitude of applications, ranging from nanotechnology to environmental chemistry. In spite of this interest, there is, however, a lack in fundamental understanding of the underlying processes in such systems. While a lot of experimental research is already being carried out to gain this understanding, only recently the first simulations have appeared in the literature. In this contribution, an overview is presented on atomic scale simulations of plasma catalytic processes as carried out in our group. In particular, this contribution focusses on plasma-assisted catalyzed carbon nanostructure growth, and plasma catalysis for greenhouse gas conversion. Attention is paid to what can routinely be done, and where challenges persist.

Abstract: Molecular dynamics simulations using improved TersoffBrenner potential parameters were performed to investigate Cl+ etching of a {2×1} reconstructed Si(100) surface. Steady-state Si etching accompanying the Cl coverage of the surface is observed. Furthermore, a steady-state chlorinated reaction layer is formed. The thickness of this reaction layer is found to increase with increasing energy. The stoichiometry of SiClx species in the reaction layer is found to be SiCl:SiCl2:SiCl3 = 1.0:0.14:0.008 at 50 eV. These results are in excellent agreement with available experimental data. While elemental Si products are created by physical sputtering, most SiClx (0<x<4) etch products are produced by chemical-enhanced physical sputtering.

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 hyperdynamics method is a powerful tool to simulate slow processes at the atomic level. However, the construction of an optimal hyperdynamics potential is a task that is far from trivial. Here, we propose a generally applicable implementation of the hyperdynamics algorithm, borrowing two concepts from metadynamics. First, the use of a collective variable (CV) to represent the accelerated dynamics gives the method a very large flexibility and simplicity. Second, a metadynamics procedure can be used to construct a suitable history-dependent bias potential on-the-fly, effectively turning the algorithm into a self-learning accelerated molecular dynamics method. This collective variable-driven hyperdynamics (CVHD) method has a modular design: both the local system properties on which the bias is based, as well as the characteristics of the biasing method itself, can be chosen to match the needs of the considered system. As a result, system-specific details are abstracted from the biasing algorithm itself, making it extremely versatile and transparent. The method is tested on three model systems: diffusion on the Cu(001) surface and nickel-catalyzed methane decomposition, as examples of reactive processes with a bond-length-based CV, and the folding of a long polymer-like chain, using a set of dihedral angles as a CV. Boost factors up to 109, corresponding to a time scale of seconds, could be obtained while still accurately reproducing correct dynamics.

Abstract: We investigated the deposition of MoS2 multilayers on large area substrates. The pre-deposition of metal or metal oxide with subsequent sulfurization is a promising technique to achieve layered films. We distinguish a different reaction behavior in metal oxide and metallic films and investigate the effect of the temperature, the H2S/H-2 gas mixture composition, and the role of the underlying substrate on the material quality. The results of the experiments suggest a MoS2 growth mechanism consisting of two subsequent process steps. At first, the reaction of the sulfur precursor with the metal or metal oxide occurs, requiring higher temperatures in the case of metallic film compared to metal oxide. At this stage, the basal planes assemble towards the diffusion direction of the reaction educts and products. After the sulfurization reaction, the material recrystallizes and the basal planes rearrange parallel to the substrate to minimize the surface energy. Therefore, substrates with low roughness show basal plane assembly parallel to the substrate. These results indicate that the substrate character has a significant impact on the assembly of low dimensional MoS2 films.

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