Abstract: Atomistic simulation methods for the quantification of free energies are in wide use. These methods operate by sampling the probability density of a system along a small set of suitable collective variables (CVs), which is, in turn, expressed in the form of a free energy surface (FES). This definition of the FES can capture the relative stability of metastable states but not that of the transition state because the barrier height is not invariant to the choice of CVs. Free energy barriers therefore cannot be consistently computed from the FES. Here, we present a simple approach to calculate the gauge correction necessary to eliminate this inconsistency. Using our procedure, the standard FES as well as its gauge-corrected counterpart can be obtained by reweighing the same simulated trajectory at little additional cost. We apply the method to a number of systems—a particle solvated in a Lennard-Jones fluid, a Diels–Alder reaction, and crystallization of liquid sodium—to demonstrate its ability to produce consistent free energy barriers that correctly capture the kinetics of chemical or physical transformations, and discuss the additional demands it puts on the chosen CVs. Because the FES can be converged at relatively short (sub-ns) time scales, a free energy-based description of reaction kinetics is a particularly attractive option to study chemical processes at more expensive quantum mechanical levels of theory.

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 work, we use density functional theory calculations to study the combined effect of external electric fields, surface morphology, and surface charge on CO2 activation over Cu(111), Cu(211), Cu(110), and Cu(001) surfaces. We observe that the binding energy of the CO2 molecule on Cu surfaces increases significantly upon increasing the applied electric field strength. In addition, rougher surfaces respond more effectively to the presence of the external electric field toward facilitating the formation of a carbonate-like CO2 structure and the transformation of the most stable adsorption mode from physisorption to chemisorption. The presence of surface charges further strengthens the electric field effect and consequently causes an improved bending of the CO2 molecule and C−O bond length elongation. On the other hand, a net charge in the absence of an externally applied electric field shows only a marginal effect on CO2 binding. The chemisorbed CO2 is more stable and further activated when the effects of an external electric field, rough surface, and surface charge are combined. These results can help to elucidate the underlying factors that control CO2 activation in heterogeneous and plasma catalysis, as well as in electrochemical processes.

Abstract: We present an approach to incorporate the effect of vibrational nonequilibrium in molecular dynamics (MD) simulations. A perturbed canonical ensemble, in which selected modes are excited to higher temperature while all others remain equilibrated at low temperature, is simulated by applying a specifically tailored bias potential. Our method can be readily applied to any (classical or quantum mechanical) MD setup at virtually no additional computational cost and allows the study of reactions of vibrationally excited molecules in nonequilibrium environments such as plasmas. In combination with enhanced sampling methods, the vibrational efficacy and mode selectivity of vibrationally stimulated reactions can then be quantified in terms of chemically relevant observables, such as reaction rates and apparent free energy barriers. We first validate our method for the prototypical hydrogen exchange reaction and then show how it can capture the effect of vibrational excitation on a symmetric SN2 reaction and radical addition on CO2.

Abstract: Carbyne is a novel material of current interest in nanotechnology. As is typically the case for nanomaterials, the growth process determines the resulting properties. While endohedral carbyne has been successfully synthesized, its catalyst and feedstock-dependent growth mechanism is still elusive. We here study the nucleation and growth mechanism of different carbon chains in a Ni-containing double walled carbon nanotube using classical molecular dynamics simulations and first-principles calculations. We find that the understanding the competitive role of the metal catalyst and the hydrocarbon is important to control the growth of 1-dimensional carbon chains, including Ni or H-terminated carbyne. Also, we find that the electronic property of the Ni-terminated carbyne can be tuned by steering the H concentration along the chain. These results suggest catalyst-containing carbon nanotubes as a possible synthesis route for carbyne formation.

Abstract: Heterostructures of low-dimensional semiconducting materials, such as transition metal dichalcogenides (MX2), are promising building blocks for future electronic and optoelectronic devices. The patterning of one MX2 material on top of another one is challenging due to their structural similarity. This prevents an intrinsic etch stop when conventional anisotropic dry etching processes are used. An alternative approach consist in a two-step process, where a sacrificial silicon layer is pre-patterned with a low damage plasma process, stopping on the underlying MoS2 film. The pre-patterned layer is used as sacrificial template for the formation of the top WS2 film. This study describes the optimization of a cyclic Ar/Cl-2 atomic layer etch process applied to etch silicon on top of MoS2, with minimal damage, followed by a selective conversion of the patterned Si into WS2. The impact of the Si atomic layer etch towards the MoS2 is evaluated: in the ion energy range used for this study, MoS2 removal occurs in the over-etch step over 1-2 layers, leading to the appearance of MoOx but without significant lattice distortions to the remaining layers. The combination of Si atomic layer etch, on top of MoS2, and subsequent Si-to-WS2 selective conversion, allows to create a WS2/MoS2 heterostructure, with clear Raman signals and horizontal lattice alignment. These results demonstrate a scalable, transfer free method to achieve horizontally individually patterned heterostacks and open the route towards wafer-level processing of 2D materials.

Abstract: Optimization of catalytic materials for a given application is greatly constrained by linear scaling relations. Recently, however, it has been demonstrated that it is possible to reversibly modulate the chemisorption of molecules on nanomaterials by charging (i.e., injection or removal of electrons) and hence reversibly and selectively modify catalytic activity beyond structure−activity correlations. The fundamental physical relation between the properties of the material, the charging process, and the chemisorption energy, however, remains unclear, and a systematic exploration and optimization of charge-switchable sorbent materials is not yet possible. Using hybrid DFT calculations of CO2 chemisorption on hexagonal boron nitride nanosheets with several types of defects and dopants, we here reveal the existence of fundamental correlations between the electron affinity of a material and charge-induced chemisorption, show how defect engineering can be used to modulate the strength and efficiency of the adsorption process, and demonstrate that excess electrons stabilize many topological defects. We then show how these insights could be exploited in the development of new electrocatalytic materials and the synthesis of doped nanomaterials. Moreover, we demonstrate that calculated chemical properties of charged materials are highly sensitive to the employed computational methodology because of the self-interaction error, which underlines the theoretical challenge posed by such systems.

Abstract: Atomic scale simulations have been shown to be a powerful tool for elucidating the growth mechanisms of carbon nanotubes. The growth picture is however not entirely clear yet due to the gap between current simulations and real experiments. We here simulate for the first time the nucleation and subsequent growth of single-wall carbon nanotubes (SWNTs) from oxygen-containing hydrocarbon feedstocks using the hybrid Molecular Dynamics/Monte Carlo technique. The underlying nucleation mechanisms of Ni-catalysed SWNT growth are discussed in detail. Specifically, we find that as a function of the feedstock, different carbon fractions may emerge as the main growth species, due to a competition between the feedstock decomposition, its rehydroxylation and its contribution to etching of the growing SWNT. This study provides a further understanding of the feedstock effects in SWNT growth in comparison with available experimental evidence as well as with<italic>ab initio</italic>and other simulation data, thereby reducing the simulation–experiment gap.

Abstract: The integration of van der Waals materials in nanoelectronic devices requires the deposition of few-layered MX2 films with excellent quality crystals covering a large area. In recent years, astonishing progress in the monolayer growth of WS2 and MoS2 was demonstrated, but multilayer growth resulted often in separated triangular or hexagonal islands. These polycrystalline films cannot fully employ the specific MX2 properties since they are not connected in-plane to the other domains. To coalesce separated islands, ultrahigh-temperature postdeposition anneals in H2S are applied, which are not compatible with bare silicon substrates. Starting from the deposition of stoichiometric short-ordered films, the present work studies different options for subsequent high-temperature annealing in an inert atmosphere to form crystalline films with large grains from stoichiometric films with small grains. The rapid thermal annealing, performed over a few seconds, is compared to excimer laser annealing in the nanosecond range, which are both able to crystallize the thin WS2. The WS2 recrystallization temperature can be lowered using metallic crystallization promoters (Co and Ni). The best result is obtained using a Co cap, due to the circumvention of Co and S binary phase formation below the eutectic temperature. The recrystallization above a critical temperature is accompanied by sulfur loss and 3D regrowth. These undesired effects can be suppressed by the application of a dielectric capping layer prior to annealing. A SiO2 cap can suppress the sulfur loss successfully during annealing and reveals improved material quality in comparison to noncapped films Published by the AVS.

Abstract: In order to optimize the plasma-synthesis and modification process of carbon nanomaterials for applications such as nanoelectronics and energy storage, a deeper understanding of fundamental hydrogengraphite/graphene interactions is required. Atomistic simulations by Molecular Dynamics have proven to be indispensable to illuminate these phenomena. However, severe time-scale limitations restrict them to very fast processes such as reflection, while slow thermal processes such as surface diffusion and molecular desorption are commonly inaccessible. In this work, we could however reach these thermal processes for the first time at time-scales and surface temperatures (1000 K) similar to high-flux plasma exposure experiments during the simulation of multilayer graphene etching by 5 eV H ions. This was achieved by applying the Collective Variable-Driven Hyperdynamics biasing technique, which extended the inter-impact time over a range of six orders of magnitude, down to a more realistic ion-flux of 1023m
2s
1. The results show that this not only causes a strong shift from predominant ion-to thermally induced interactions, but also significantly affects the hydrogen uptake and surface evolution. This study thus elucidates H ion-graphite/graphene interaction mechanisms and stresses the importance of including long time-scales in atomistic simulations at high surface temperatures to understand the dynamics of the ion-surface system.

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: 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: In this study, the healing of N-vacancy boron carbonitride nanosheet (NV-BC2NNS) and nanotube (NV-BC2NNT) by NO molecule is studied by means of density functional theory calculations. Two different N-vacancies are considered in each of these structures in which the vacancy site is surrounded by either three B-atoms (NB) or by two B- and one C-atom (NBC). By means of the healed BC2NNS and BC2NNT as a support, the removal of two toxic gas molecules (NO and CO) are applicable. It should be noted that the obtained energy barriers of both healing and oxidizing processes are significantly lower than those of graphene, carbon nanotubes or boron nitride nanostructures. Also, at the end of the oxidation process, the pure BC2NNS or BC2NNT is obtained without any additional defects. Therefore, by using this method, we can considerably purify the defective BC2NNS/BC2NNT. Moreover, according to the thermochemistry calculations we can further confirm that the healing process of the NV-BC2NNS and NV-BC2NNT by NO are feasible at room temperature. So, we can claim that this study could be very helpful in both purifying the defective BC2NNS/BC2NNT while in the same effort removing toxic NO and CO gases.

Abstract: Catalyst preparation with plasmas is increasingly attracting interest. A plasma is a partially ionized gas, consisting of electrons, ions, molecules, radicals, photons, and excited species, which are all active species for catalyst preparation and treatment. Under the influence of plasma, nucleation and crystal growth in catalyst preparation can be very different from those in the conventional thermal approach. Some thermodynamically unfavorable reactions can easily take place with plasmas. Compounds such as sulfides, nitrides, and phosphides that are produced under harsh conditions can be synthesized by plasma under mild conditions. Plasmas can produce catalysts with smaller particle sizes and controllable structure. Plasma is also a facile tool for reduction, oxidation, doping, etching, coating, alloy formation, surface treatment, and surface cleaning in a simple and direct way. A rapid and convenient plasma template removal has thus been established for zeolite synthesis. It can operate at room temperature and allows the catalyst preparation on temperature-sensitive supporting materials. Plasma is typically effective for the production of various catalysts on metallic substrates. In addition, plasma-prepared transition-metal catalysts show enhanced low-temperature activity with improved stability. This provides a useful model catalyst for further improvement of industrial catalysts. In this review, we aim to summarize the recent advances in catalyst preparation with plasmas. The present understanding of plasma-based catalyst preparation is discussed. The challenges and future development are addressed.

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.