Abstract: Using density functional theory in combination with the Green’s functional formalism, we study the effect of surface functionalization on the electronic transport properties of 1D carbon allotrope—carbyne. We found that both hydrogenation and fluorination result in structural changes and semiconducting to metallic transition. Consequently, the current in the functionalization systems increases significantly due to strong delocalization of electronic states along the carbon chain. We also study the electronic transport in partially hydrogenated carbyne and interface structures consisting of pristine and functionalized carbyne. In the latter case, current rectification is obtained in the system with rectification ratio up to 50%. These findings can be useful for developing carbyne-based structures with tunable electronic transport 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: The cleaning of two-dimensional (2D) materials is an essential step in the fabrication of future devices, leveraging their unique physical, optical, and chemical properties. Part of these emerging 2D materials are transition metal dichalcogenides (TMDs). So far there is limited understanding of the cleaning of “monolayer” TMD materials. In this study, we report on the use of downstream H₂plasma to clean the surface of monolayer WS₂grown by MOCVD. We demonstrate that high-temperature processing is essential, allowing to maximize the removal rate of polymers and to mitigate damage caused to the WS₂in the form of sulfur vacancies. We show that low temperature in situ carbonyl sulfide (OCS) soak is an efficient way to resulfurize the material, besides high-temperature H₂S annealing. The cleaning processes and mechanisms elucidated in this work are tested on back-gated field-effect transistors, confirming that transport properties of WS₂devices can be maintained by the combination of H₂plasma cleaning and OCS restoration. The low-damage plasma cleaning based on H₂and OCS is very reproducible, fast (completed in a few minutes) and uses a 300 mm industrial plasma etch system qualified for standard semiconductor pilot production. This process is, therefore, expected to enable the industrial scale-up of 2D-based devices, co-integrated with silicon technology.

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: 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: Metal-catalyzed growth mechanisms of carbon nanotubes (CNTs) were studied by hybrid molecular dynamics−Monte Carlo simulations using a recently developed ReaxFF reactive force field. Using this novel approach, including relaxation effects, a CNT with definable chirality is obtained, and a step-by-step atomistic description of the nucleation process is presented. Both root and tip growth mechanisms are observed. The importance of the relaxation of the network is highlighted by the observed healing of defects.

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: 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: 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: Thermal-catalytic gas processing is integral to many current industrial processes. Ever-increasing demands on conversion and energy efficiencies are a strong driving force for the development of alternative approaches. Similarly, synthesis of several functional materials (such as nanowires and nanotubes) demands special processing conditions. Plasma catalysis provides such an alternative, where the catalytic process is complemented by the use of plasmas that activate the source gas. This combination is often observed to result in a synergy between plasma and catalyst. This Review introduces the current state-of-the-art in plasma catalysis, including numerous examples where plasma catalysis has demonstrated its benefits or shows future potential, including CO2 conversion, hydrocarbon reforming, synthesis of nanomaterials, ammonia production, and abatement of toxic waste gases. The underlying mechanisms governing these applications, as resulting from the interaction between the plasma and the catalyst, render the process highly complex, and little is known about the factors leading to the often-observed synergy. This Review critically examines the catalytic mechanisms relevant to each specific application.

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: Reactive molecular dynamics simulations are performed to study the interaction of reactive oxygen species, such as OH, HO2 and H2O2, with the endotoxic biomolecule lipid A of the gram-negative bacterium Escherichia coli. It is found that the aforementioned plasma species can destroy the lipid A, which consequently results in reducing its toxic activity. All bond dissociation events are initiated by hydrogen-abstraction reactions. However, the mechanisms behind these dissociations are dependent on the impinging plasma species, i.e. a clear difference is observed in the mechanisms upon impact of HO2 radicals and H2O2 molecules on one hand and OH radicals on the other hand. Our simulation results are in good agreement with experimental observations.

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.