Abstract: The free energy surface (FES) for carbon segregation from nickel nanoparticles is obtained from advanced molecular dynamics simulations. A suitable reaction coordinate is developed that can distinguish dissolved carbon atoms from segregated dimers, chains and junctions on the nanoparticle surface. Because of the typically long segregation time scale (up to ms), metadynamics simulations along the developed reaction coordinate are used to construct FES over a wide range of temperatures and carbon concentrations. The FES revealed the relative stability of different stages in the segregation process, and free energy barriers and rates of the individual steps could then be calculated and decomposed into enthalpic and entropic contributions. As the carbon concentration in the nickel nanoparticle increases, segregated carbon becomes more stable in terms of both enthalpy and entropy. The activation free energy of the reaction also decreases with the increase of carbon concentration, which can be mainly attributed to entropic effects. These insights and the methodology developed to obtain them improve our understanding of carbon segregation process across materials science in general, and the nucleation and growth of carbon nanotube in particular.

Abstract: Plasma catalysis continues to gain increasing scientific interest, both in established fields like toxic waste abatement and emerging fields like greenhouse gas conversion into value-added chemicals. Attention is typically focused on the obtained conversion process selectivity, rates and energy efficiency. Much less attention is usually paid to the underlying mechanistic aspects of the processes that occur. In this contribution, we critically examine a number of fundamentally important nanoscale thermodynamic aspects of plasma catalysis, which are very relevant to these processes but so far have been overlooked or insufficiently covered in the plasma catalysis literature.

Abstract: In recent years, the potential use of hydrogen as a clean energy source has gained considerable attention. Especially H2 formation by Ni-catalyzed reforming of methane at elevated temperatures is an attractive process. However, a more fundamental knowledge at the atomic level is needed for a full comprehension of the reactions at the catalyst surface. In this contribution, we therefore investigate the H2 formation after CHx impacts on a Ni(1 1 1) surface in the temperature range 4001600 K, by means of reactive molecular dynamics (MD) simulations using the ReaxFF potential. While some H2 formation is already observed at the lower temperatures, substantial H2 formation is only obtained at elevated temperatures of 1400 K and above. At 1600 K, the H2 molecules are even the most frequently formed species. In direct correlation with the increasing dehydrogenation at elevated temperatures, an increased surface-to-subsurface C-diffusivity is observed as well. This study highlights the major importance of the temperature on the H2 formation.

Abstract: We performed molecular dynamics simulations to investigate the effect of lipid peroxidation products on the structural and dynamic properties of the cell membrane. Our simulations predict that the lipid order in a phospholipid bilayer, as a model system for the cell membrane, decreases upon addition of lipid peroxidation products. Eventually, when all phospholipids are oxidized, pore formation can occur. This will allow reactive species, such as reactive oxygen and nitrogen species (RONS), to enter the cell and cause oxidative damage to intracellular macromolecules, such as DNA or proteins. On the other hand, upon increasing the cholesterol fraction of lipid bilayers, the cell membrane order increases, eventually reaching a certain threshold, from which cholesterol is able to protect the membrane against pore formation. This finding is crucial for cancer treatment by plasma technology, producing a large number of RONS, as well as for other cancer treatment methods that cause an increase in the concentration of extracellular RONS. Indeed, cancer cells contain less cholesterol than their healthy counterparts. Thus, they will be more vulnerable to the consequences of lipid peroxidation, eventually enabling the penetration of RONS into the interior of the cell, giving rise to oxidative stress, inducing pro-apoptotic factors. This provides, for the first time, molecular level insight why plasma can selectively treat cancer cells, while leaving their healthy counterparts undamaged, as is indeed experimentally demonstrated.

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.

Abstract: Natural deep eutectic solvents (NADES) are mixtures of naturally derived compounds with a significantly decreased melting point due to the specific interactions among the constituents. NADES have benign properties (low volatility, flammability, toxicity, cost) and tailorable physicochemical properties (by altering the type and molar ratio of constituents), hence they are often considered as a green alternative to common organic solvents. Modeling the relation between their composition and properties is crucial though, both for understanding and predicting their behavior. Several efforts were done to this end, yet this review aims at structuring the present knowledge as an outline for future research. First, we reviewed the key properties of NADES and relate them to their structure based on the available experimental data. Second, we reviewed available modeling methods applicable to NADES. At the molecular level, density functional theory and molecular dynamics allow interpreting density differences and vibrational spectra, and computation of interaction energies. Additionally, properties at the level of the bulk media can be explained and predicted by semi-empirical methods based on ab initio methods (COSMO-RS) and equation of state models (PC-SAFT). Finally, methods based on large datasets are discussed; models based on group contribution methods and machine learning. A combination of bulk media and dataset modeling allows qualitative prediction and interpretation of phase equilibria properties on the one hand, and quantitative prediction of melting point, density, viscosity, surface tension and refractive indices on the other hand. In our view, multiscale modeling, combining the molecular and macroscale methods, will strongly enhance the predictability of NADES properties and their interaction with solutes, yielding truly tailorable solvents to accommodate (bio)chemical reactions.

Abstract: In this paper we employ first-principles calculations to investigate the effect of substitutional Si doping in the amorphous calcium-phosphate (a-HAP) structure on the work of adhesion, integral charge transfer, charge density difference and theoretical tensile strengths between an a-HAP coating and amorphous titanium dioxide (a-TiO2) substrate systemically. Our calculations demonstrate that substitution of a P atom by a Si atom in a-HAP (a-Si-HAP) with the creation of OH-vacancies as charge compensation results in a significant increase of the bonding strength of the coating to the substrate. The work of adhesion of the optimized Si-doped interfaces reaches a value of up to -2.52 J m(-2), which is significantly higher than for the stoichiometric a-HAP/a-TiO2. Charge density difference analysis indicates that the dominant interactions at the interface have significant covalent character, and in particular two Ti-O and three Ca-O bonds are formed for a-Si-HAP/a-TiO2 and one Ti-O and three Ca-O bonds for a-HAP/a-TiO2. From the stress-strain curve, the Young's modulus of a-Si-HAP/a-TiO2 is calculated to be about 25% higher than that of the a-HAP/a-TiO2, and the yielding stress is about 2 times greater than that of the undoped model. Our calculations therefore demonstrate that the presence of Si in the a-HAP structure strongly alters not only the bioactivity and resorption rates, but also the mechanical properties of the a-HAP/a-TiO2 interface. The results presented here provide an important theoretical insight into the nature of the chemical bonding at the a-HAP/a-TiO2 interface, and are particularly significant for the practical medical applications of HAP-based biomaterials.

Abstract: The possibility of application of nanomaterials is determined by our ability to control the properties of the materials, which are ultimately determined by their structure and hence their growth processes. We employ hybrid molecular dynamics / Monte Carlo (MD/MC) simulations to explore the growth of SWNTs and graphene on nickel as a catalyst, with the specific goal of unraveling the growth mechanisms. While the general observations are in agreement with the literature, we find a number of interesting phenomena to be operative which are crucial for the growth, and which are not accessible by MD simulations alone due to the associated time scale. Specifically, we observe metal mediated healing and restructuring processes to take place, reorganizing the carbon network during the initial nucleation step. In the case of carbon nanotube growth, this leads to the growth of tubes with a determinable chirality. In the case of graphene formation, we find that graphene is only formed at temperatures above 700 K. These results are of importance for understanding the growth mechanisms of these carbon nanomaterials on the fundamental level.

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: The safe and efficient use of hydrogen energy, which is in high demand worldwide today, requires efficient hydrogen storage. Despite significant advances in hydrogen storage using carbon-based nanomaterials, including carbon nanotubes (CNTs), efforts to substantially increase the storage capacity remain less effective. In this work, we demonstrate the effect of endohedral transition metal atoms on the hydrogen storage capacity of CNTs using reactive molecular dynamics simulations. We find that an increase in the volume fraction of endohedral nickel atoms leads to an increase in the concentration of physisorbed hydrogen molecules around single-walled CNTs (SWNTs) by approximately 1.6 times compared to pure SWNTs. The obtained results provide insight into the underlying mechanisms of how endohedral transition metal atoms enhance the hydrogen storage ability of SWNTs under nearly ambient conditions.

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.