Abstract: Plasma-catalytic CO, hydrogenation for methanol production is gaining increasing interest, but our understanding of its reaction mechanism remains primitive. We present a combined experimental/computational study on plasma-catalytic CO, hydrogenation to CH3OH over a size-selected Cu/gamma-Al2O3 catalyst. Our experiments demonstrate a synergistic effect between the Cu/gamma-Al2O3 catalyst and the CO2/H-2 plasma, achieving a CO2 conversion of 10% at 4 wt % Cu loading and a CH3OH selectivity near 50% further rising to 65% with H2O addition (for a H2O/CO2 ratio of 1). Furthermore, the energy consumption for CH3OH production was more than 20 times lower than with plasma only. We carried out density functional theory calculations over a Cu-13/gamma-Al2O3 model, which reveal that the interfacial sites of the Cu-13 cluster and gamma-Al2O3 support show a bifunctional effect: they not only activate the CO2 molecules but also strongly adsorb key intermediates to promote their hydrogenation further. Reactive plasma species can regulate the catalyst surface reactions via the Eley-Rideal (E-R) mechanism, which accelerates the hydrogenation process and promotes the generation of the key intermediates. H2O can promote the CH3OH desorption by competitive adsorption over the Cu-13/gamma-Al2O3 surface. This study provides new insights into CO2 hydrogenation through plasma catalysis, and it provides inspiration for the conversion of some other small molecules (CH4, N-2, CO, etc.) by plasma catalysis using supported-metal clusters.

Abstract: Plasma-surface interactions are in general highly complex due to the interplay of many concurrent processes. Molecular dynamics simulations provide insight in some of these processes, subject to the accessible time and length scales, and the availability of suitable force fields. In this introductory tutorial-style review, we aim to describe the current capabilities and limitations of molecular dynamics simulations in this field, restricting ourselves to low-temperature nonthermal plasmas. Attention is paid to the simulation of the various fundamental processes occurring, including sputtering, etching, implantation, and deposition, as well as to what extent the basic plasma components can be accounted for, including ground state and excited species, electric fields, ions, photons, and electrons. A number of examples is provided, giving an bird’s eye overview of the current state of the field.

Abstract: Purpose Natural deep eutectic solvents (NADES) represent a green alternative to conventional organic solvents as reaction medium, offering more benign properties. To efficiently design NADES for biocatalysis, a better understanding of their effect on these reactions is needed. We hypothesize that this effect can be described by separately considering (1) the solvent interactions with the substrates, (2) the solvent viscosities and (3) the enzyme stability in NADES. Methods We investigated the effect of substrate solvation and viscosity on the reaction rate; and the stability of the enzyme in NADES. To this end, we monitored the conversion over time of the transesterification of vinyl laurate with 1- butanol by the lipase enzyme Candida antarctica B in NADES of different compounds and molar ratios. Results The initial reaction rate is higher in most NADES ( varying between 1.14 and 15.07 mu mol min(-1) mg(-1)) than in the reference n-hexane (4.0 mu mol min(-1) mg(-1))), but no clear relationship between viscosity and initial reaction rate was found. The increased reaction rate is most likely related to the solvation of the substrate due to a change in the activation energy of the reaction or a change in the conformation of the substrate. The enzyme retained part of its activity after the first 2 h of reaction (on average 20 % of the substrate reacted in the 2-24 h period). Enzyme incubation in ethylene glycol-based NADES resulted in a reduced reaction rate ( 15.07 vs. 3.34 mu mol min(-1) mg(-1)), but this may also be due to slow dissolution of the substrate. Conclusions The effect of viscosity seems to be marginal next to the effect of solvation and possible enzyme-NADES interaction. The enzyme retains some of its activity during the 24-hour measurements, but the enzyme incubation experiments did not yield accurate, comparable values. [GRAPHICS] .

Abstract: Recent experiments have revealed that the macroscopic Kelvin equation remains surprisingly accurate even for nanoscale capillaries. This phenomenon was so far explained by the oscillatory behavior of the solid−liquid interfacial free energy. We here demonstrate thermodynamic and capillarity inconsistencies with this explanation. After revising the Kelvin equation, we ascribe its validity at nanoscale confinement to the effect of disjoining pressure.
To substantiate our hypothesis, we employed molecular dynamics simulations to evaluate interfacial heat transfer and wetting properties. Our assessments unveil a breakdown in a previously established proportionality between the work of adhesion and the Kapitza conductance at capillary heights below 1.3 nm, where the dominance of the work of adhesion shifts primarily from energy to entropy. Alternatively, the peak density of the initial water layer can effectively probe the work of adhesion. Unlike under bulk conditions, high confinement renders the work of adhesion entropically unfavorable.

Abstract: Recently, Sr-doped LaCrO3 has been experimentally introduced as a new p-type transparent conducting oxide. It is demonstrated that substituting Sr for La results in inducing p-type conductivity in LaCrO3. Performing first principles calculations we study the electronic structure and formation energy of various point defects in LaCrO3. Our results for the formation energies show that in addition to Sr, two more divalent defects, Ca and Ba, substituting for La in LaCrO3, behave as shallow acceptors in line with previous experimental reports. We further demonstrate that under oxygen-poor growth conditions, these shallow acceptors will be compensated by intrinsic donor-like defects (an oxygen vacancy and Cr on an oxygen site), but in the oxygen-rich growth regime the shallow acceptors have the lowest formation energies between all considered defects and will lead to p-type conductivity.

Abstract: Structural defects inevitably appear during the nucleation event that determines the structure and properties of single-walled carbon nanotubes. By combining ion bombardment experiments with atomistic simulations we reveal that ion bombardment in a suitable energy range allows these defects to be healed resulting in an enhanced nucleation of the carbon nanotube cap. The enhanced growth of the nanotube cap is explained by a nonthermal ion-induced graphene network restructuring mechanism.

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: The adsorption, dissociation, and diffusion of CO2 on the anatase (001) surface was studied using DFT by means of the generalized gradient approximation using the Perdew−Burcke−Ernzerhof (PBE)-functional and applying corrections for long-range dispersion interactions. Different stable adsorption configurations were identified for the fully oxidized surface. The most stable adsorption configuration is the monodentated carbonate-like structure. Small energy barriers were identified for the conversion of a physisorbed to a chemisorbed configuration.
CO2 dissociation is found to be unfeasible on the stoichiometric surface. The introduction of oxygen vacancy defects gives rise to new highly stable adsorption configurations with a stronger activation of the C−O bonds. This leads to the possibility of exothermic dissociation of CO2 with barriers up to 22.2 kcal/mol,
corresponding to chemical lifetimes of less than 4 s at 300 K. These reactions cause a CO molecule to be formed, which will easily desorb, and the reduced surface to become oxidized. It is clear that oxygen vacancy defects play a key role in the catalytic activity of an anatase (001) surface. Oxygen vacancies play an important role in the dissociation of CO2 on the anatase (001) surface, and will play a significant role in complex problems, such as the catalytic conversion of CO2 to value-added chemicals.

Abstract: The application of coreshell SiSiO2 nanowires as nanoelectronic devices strongly depends on their structure, which is difficult to tune precisely. In this work, we investigate the formation of the coreshell nanowires at the atomic scale, by reactive molecular dynamics simulations. The occurrence of two temperature-dependent oxidation mechanisms of ultra-small diameter Si-NWs is demonstrated. We found that control over the Si-core radius and the SiOx (x ≤ 2) oxide shell is possible by tuning the growth temperature and the initial Si-NW diameter. Two different structures were obtained, i.e., ultrathin SiO2 silica nanowires at high temperature and Si core|ultrathin SiO2 silica nanowires at low temperature. The transition temperature is found to linearly decrease with the nanowire curvature. Finally, the interfacial stress is found to be responsible for self-limiting oxidation, depending on both the initial Si-NW radius and the oxide growth temperature. These novel insights allow us to gain control over the exact morphology and structure of the wires, as is needed for their application in nanoelectronics.

Abstract: The behavior of hydrocarbon species at step edges of diamond terraces is investigated by means of combined molecular dynamics−Metropolis Monte Carlo simulations. The results show that the formation of ballas-like diamond films (like UNCD) and well-faceted diamond films (like NCD) can be related to the gas phase concentrations of CxHy in a new manner: Species that have high concentrations above the growing UNCD films suppress the extension of step edges through defect formation. The species that are present above the growing NCD film, however, enhance the extension of diamond terraces, which is believed to result in well-faceted diamond films. Furthermore, it is shown that, during UNCD growth, CxHy species with x ≥ 2 play an important role, in contrast to the currently adopted CVD diamond growth mechanism. Finally, the probabilities for the extension of the diamond (100) terrace are much higher than those for the diamond (111) terrace, which is in full agreement with the experimental observation that diamond (100) facets are more favored than diamond (111) facets during CVD diamond growth.

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: The continuous miniaturization of nanodevices, such as transistors, solar cells, and optical fibers, requires the controlled synthesis of (ultra)thin gate oxides (<10 nm), including Si gate-oxide (SiO2) with high quality at the atomic scale. Traditional thermal growth of SiO2 on planar Si surfaces, however, does not allow one to obtain such ultrathin oxide due to either the high oxygen diffusivity at high temperature or the very low sticking ability of incident oxygen at low temperature. Two recent techniques, both operative at low (room) temperature, have been put forward to overcome these obstacles: (i) hyperthermal oxidation of planar Si surfaces and (ii) thermal or plasma-assisted oxidation of nonplanar Si surfaces, including Si nanowires (SiNWs). These nanooxidation processes are, however, often difficult to study experimentally, due to the key intermediate processes taking place on the nanosecond time scale.
In this Account, these Si nano-oxidation techniques are discussed from a computational point of view and compared to both hyperthermal and thermal oxidation experiments, as well as to well-known models of thermal oxidation, including the Deal−Grove, Cabrera−Mott, and Kao models and several alternative mechanisms. In our studies, we use reactive molecular dynamics (MD) and hybrid MD/Monte Carlo simulation techniques, applying the Reax force field. The incident energy of oxygen species is chosen in the range of 1−5 eV in hyperthermal oxidation of planar Si surfaces in order to prevent energy-induced damage. It turns out that hyperthermal growth allows for two growth modes, where the ultrathin oxide thickness depends on either (1) only the kinetic energy of the incident oxygen species at a growth temperature below Ttrans = 600 K, or (2) both the incident energy and the growth temperature at a growth temperature above Ttrans. These modes are specific to such ultrathin oxides, and are not observed in traditional thermal oxidation, nor theoretically considered by already existing models. In the case of thermal or plasma-assisted oxidation of small Si nanowires, on the other hand, the thickness of the ultrathin oxide is a function of the growth temperature and the nanowire diameter. Below Ttrans, which varies with the nanowire diameter, partially oxidized SiNW are formed, whereas complete oxidation to a SiO2 nanowire occurs only above Ttrans. In both nano-oxidation processes at lower temperature (T < Ttrans), final sandwich c-Si|SiOx|a-SiO2 structures are obtained due to a competition between overcoming the energy barrier to penetrate into Si subsurface layers and the compressive stress (∼2−3 GPa) at the Si crystal/oxide interface. The overall atomic-simulation results strongly indicate that the thickness of the intermediate SiOx (x < 2) region is very limited (∼0.5 nm) and constant irrespective of oxidation parameters. Thus, control over the ultrathin SiO2 thickness with good quality is indeed possible by accurately tuning the oxidant energy, oxidation temperature and surface curvature.
In general, we discuss and put in perspective these two oxidation mechanisms for obtaining controllable ultrathin gate-oxide films, offering a new route toward the fabrication of nanodevices via selective nano-oxidation.

Abstract: The ever increasing global production and dispersion of methane requires novel chemistry to transform it into easily condensable energy carriers that can be integrated into the chemical infrastructure. In this context, single atom catalysts have attracted considerable interest due to their outstanding catalytic activity. We here use density functional theory (DFT) computations to compare the reaction and activation energies of M and MN4 embedded graphene (M = Ni and Si) on the methane-to-methanol conversion near room temperature. Thermodynamically, conversion of methane to methanol is energetically favorable at ambient conditions. Both singlet and triplet spin state of the studied systems are considered in all of the calculations. The DFT results show that the barriers are significantly lower when the complexes are in the triplet state than in the singlet state. In particular, Si-G with the preferred spin multiplicity of triplet seems to be viable catalysts for methane oxidation thanks to the corresponding lower energy barriers and higher stability of the obtained configurations. Our results provide insights into the nature of methane conversion and may serve as guidance for fabricating cost-effective graphene-based single atom catalysts.

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