Abstract: Based on first-principles calculations, we study the structural stability of the H and T phases of monolayer MoS2 upon Li doping. Our calculations demonstrate that it is possible to stabilize a distorted T phase of MoS2 over the H phase through adsorption of Li atoms on the MoS2 surface. Through molecular dynamics and phonon calculations, we show that the T phase of MoS2 is dynamically unstable and undergoes considerable distortions. The type of distortion depends on the concentration of adsorbed Li atoms and changes from zigzag-like to diamond-like when increasing the Li doping. There exists a substantial energy barrier to transform the stable H phase to the distorted T phases, which is considerably reduced by increasing the concentration of Li atoms. We show that it is necessary that the Li atoms adsorb on both sides of the MoS2 monolayer to reduce the barrier sufficiently. Two processes are examined that allow for such two-sided adsorption, namely, penetration through the MoS2 layer and diffusion over the MoS2 surface. We show that while there is only a small barrier of 0.24 eV for surface diffusion, the amount of energy needed to pass through a pure MoS2 layer is of the order of similar or equal to 2 eV. However, when the MoS2 layer is covered with Li atoms the amount of energy that Li atoms should gain to penetrate the layer is drastically reduced and penetration becomes feasible.

Abstract: The quantum electrostatic heterostructure (QEH) model allows for efficient computation of the dielectric screening properties of layered van der Waals (vdW)-bonded heterostructures in terms of the dielectric functions of the individual two-dimensional (2D) layers. Here, we extend the QEH model by including (1) contributions to the dielectric function from infrared active phonons in the 2D layers, (2) screening from homogeneous bulk substrates, and (3) intraband screening from free carriers in doped 2D semiconductor layers. We demonstrate the potential of the extended QEH model by calculating the dispersion of coupled phonons in multilayer stacks of hexagonal boron-nitride (hBN), the strong hybridization of plasmons and optical phonons in graphene/hBN heterostructures, the effect of substrate screening on the exciton series of monolayer MoS2, and the properties of hyperbolic plasmons in a doped phosphorene sheet. The new QEH code is distributed as a Python package with a simple command line interface and a comprehensive library of dielectric building blocks for the most common 2D materials, providing an efficient open platform for dielectric modeling of realistic vdW heterostructures.

Abstract: To gain insight into the nature of the adhesion mechanism between hydroxyapatite (HA) and rutile (rTiO(2)), the mutual affinity between their surfaces was systematically studied using density functional theory (DFT). We calculated both bulk and surface properties of HA and rTiO(2), and explored the interfacial bonding mechanism of amorphous HA (aHA) surface onto amorphous as well as stoichiometric and nonstoichiometric crystalline rTiO(2). Formation energies of bridging and subbridging oxygen vacancies considered in the rTiO(2)(110) surface were evaluated and compared with other theoretical and experimental results. The interfacial interaction was evaluated through the work of adhesion. For the aHA/rTiO(2)(110) interfaces, the work of adhesion is found to depend strongly on the chemical environment of the rTiO(2)(110) surface. Electronic analysis indicates that the charge transfer is very small in the case of interface formation between aHA and crystalline rTiO(2)(110). In contrast, significant charge transfer occurs between aHA and amorphous rTiO(2) (aTiO(2)) slabs during the formation of the interface. Charge density difference (CDD) analysis indicates that the dominant interactions in the interface have significant covalent character, and in particular the Ti-O and Ca-O bonds. Thus, the obtained results reveal that the aHA/aTiO(2) interface shows a more preferable interaction and is thermodynamically more stable than other interfaces. These results are particularly important for improving the long-term stability of HA-based implants.

Abstract: By means of chemical kinetics modeling, it is possible to elucidate the main dissociation mechanisms of CO2 in a gliding arc plasmatron (GAP). We obtain good agreement between the calculated and experimental conversions and energy efficiencies, indicating that the model can indeed be used to study the underlying mechanisms. The calculations predict that vibration-induced dissociation is the main dissociation mechanism of CO2, but it occurs mainly from the lowest vibrational levels because of fast thermalization of the vibrational distribution. Based on these findings, we propose ideas for improving the performance of the GAP, but testing of these ideas in the simulations reveals that they do not always lead to significant enhancement, because of other side effects, thus illustrating the complexity of the process. Nevertheless, the model allows more insight into the underlying mechanisms to be obtained and limitations to be identified.

Abstract: A chemical kinetics model is developed for a CO2/N2 microwave plasma, focusing especially on the vibrational levels of both CO2 and N2. The model is used to calculate the CO2 and N2 conversion as well as the energy efficiency of CO2 conversion for different power densities and for N2 fractions in the CO2/N2 gas mixture ranging from 0 to 90%. The calculation results are compared with measurements, and agreements within 23% and 33% are generally found for the CO2 conversion and N2 conversion, respectively. To explain the observed trends, the destruction and formation processes of both CO2 and N2 are analyzed, as well as the vibrational distribution functions of both CO2 and N2. The results indicate that N2 contributes in populating the lower asymmetric levels of CO2, leading to a higher absolute CO2 conversion upon increasing N2 fraction. However, the effective CO2 conversion drops because there is less CO2 initially present in the gas mixture; thus, the energy efficiency also drops with rising N2 fraction.

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 structural, electronic, and magnetic properties of pristine, defective, and oxidized monolayer TiS3 are investigated using first-principles calculations in the framework of density functional theory. We found that a single layer of TiS3 is a direct band gap semiconductor, and the bonding nature of the crystal is fundamentally different from other transition metal chalcogenides. The negatively charged surfaces of single layer TiS3 makes this crystal a promising material for lubrication applications. The formation energies of possible vacancies, i.e. S, Ti, TiS, and double S, are investigated via total energy optimization calculations. We found that the formation of a single S vacancy was the most likely one among the considered vacancy types. While a single S vacancy results in a nonmagnetic, semiconducting character with an enhanced band gap, other vacancy types induce metallic behavior with spin polarization of 0.3-0.8 mu(B). The reactivity of pristine and defective TiS3 crystals against oxidation was investigated using conjugate gradient calculations where we considered the interaction with atomic O, O-2, and O-3. While O-2 has the lowest binding energy with 0.05-0.07 eV, O-3 forms strong bonds stable even at moderate temperatures. The strong interaction (3.9-4.0 eV) between atomic O and TiS3 results in dissociative adsorption of some O-containing molecules. In addition, the presence of S-vacancies enhances the reactivity of the surface with atomic O, whereas it had a negative effect on the reactivity with O-2 and O-3 molecules.

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: To determine and control the spatial confinement of charge carriers is of importance for nanoscale optoelectronic device applications. Using first-principles calculations, we investigate the tunability of band alignment and Charge localization in lateral and combined lateral vertical heterostructures of MoS2 and WS2. First, we Show that a type-II to type-I band alignment transition takes place when tensile strain is applied on the WS2 region. This band alignment transition is a result of the different response of the band edge states with strain and is caused by their different wave function characters. Then we show that the presence of the grain boundary introduces localized in-gap states. The boundary at the armchair interface significantly modifies the charge distribution of the valence band maximum (VBM) state, whereas in a heterostructure with tilt grain domains both conducation band maximum (CBM) and VBM are found to be localized around the grain boundary. We also found that the thickness of the constituents in a lateral heterostructure also determines how the electrons and holes are confined. Creating combined lateral vertical heterostructures of MOS2/WS2 provides another way cif tuning the charge confinement. These results provide possible ways to tune the carrier confinement in MoS2/WS2 heterostructures, which are interesting for its practical: applications in the future.

Abstract: Using reactive molecular dynamics simulations by means of the ReaxFF potential, we studied the growth mechanism of ultrathin silica (SiO2) layers during hyperthermal oxidation at room temperature. Oxidation of Si(100){2 × 1} surfaces by both atomic and molecular oxygen was investigated in the energy range 15 eV. The oxidation mechanism, which differs from thermal oxidation, is discussed. In the case of oxidation by molecular O2, silica is quickly formed and the thickness of the formed layers remains limited compared to oxidation by atomic oxygen. The Si/SiO2 interfaces are analyzed in terms of partial charges and angle distributions. The obtained structures of the ultrathin SiO2 films are amorphous, including some intrinsic defects. This study is important for the fabrication of silica-based devices in the micro- and nanoelectronics industry, and more specifically for the fabrication of metal oxide semiconductor devices.

Abstract: Using reactive molecular dynamics simulations based on the ReaxFF potential, we studied the growth mechanism of ultrathin silica (SiO2) layers during hyperthermal oxidation as a function of temperature in the range 100-1300 K. Oxidation of Si(100){2 x 1} surfaces by both atomic and molecular oxygen was investigated for hyperthermal impact energies in the range of 1 to 5 eV. Two different growth mechanisms are found, corresponding to a low temperature oxidation and a high temperature one. The transition temperature between these mechanisms is estimated to be about 700 K. Also, the initial step of the Si oxidation process is analyzed in detail. Where possible, we validated our results with experimental and ab initio data, and good agreement was obtained. This study is important for the fabrication of silica-based devices in the micro- and nanoelectronics industry and, more specifically, for the fabrication of metal oxide semiconductor devices.

Abstract: Device performance of Si|SiOx core-shell based nanowires critically depends on the exact control over the oxide thickness. Low-temperature plasma oxidation is a highly promising alternative to thermal oxidation allowing for improved control over the oxidation process, in particular for ultrasmall Si nanowires. We here elucidate the room temperature plasma oxidation mechanisms of ultrasmall Si nanowires using hybrid molecular dynamics / force-bias Monte Carlo simulations. We demonstrate how the oxidation and concurrent water formation mechanisms are a function of the oxidizing plasma species and we demonstrate how the resulting core-shell oxide thickness can be controlled through these species. A new mechanism of water formation is discussed in detail. The results provide a detailed atomic level explanation of the oxidation process of highly curved Si surfaces. These results point out a route toward plasma-based formation of ultrathin core-shell Si|SiOx nanowires at room temperature.

Abstract: Semiconductor monolayer transition metal dichalcogenides (TMDs) have brought a new paradigm by introducing optically addressable valley degree of freedom. Concomitantly, their high flexibility constitutes a unique platform that links optics to mechanics via valleytronics. With the intention to expedite the research in this direction, we investigated ten TMDs, namely MoS2, MoSe2, MoTe2, WS2, WSe2, WTe2, MoSSe, MoSeTe, WSSe, and WSeTe, which particularly includes their so-called janus types (JTMDs). First, we obtained their electronic band structures using regular and hybrid density functional theory (DFT) calculations in the presence of the spin-orbit coupling and biaxial or uniaxial strain. Our DFT results indicated that against the expectations based on their reported piezoelectric behavior, JTMDs typically interpolated between the standard band properties of the constituent TMDs without producing a novel feature. Next, by fitting to our DFT data we generated both spinless and spinful k center dot p parameter sets which are quite accurate over the K valley where the optical activity occurs. As an important application of this parametrization, we considered the circular and linear dichroism under strain. Among the studied (J)TMDs, WTe2 stood out with its largest linear dichroism under uniaxial strain because of its narrower band gap and large K valley uniaxial deformation potential. This led us to suggest WTe2 monolayer membranes for optical polarization-based strain measurements, or conversely, as strain tunable optical polarizers.

Abstract: Anatase nanosheets with exposed {001} facets
have gained increasing interest for photocatalytic applications. To
fully understand the structure-to-activity relation, combined
experimental and computational methods have been exploited.
Anatase nanosheets were prepared under hydrothermal conditions
in the presence of fluorine ions. High resolution scanning
transmission electron microscopy was used to fully characterize
the synthesized material, confirming the TiO2 nanosheet
morphology. Moreover, the surface structure and composition
of a single nanosheet could be determined by annular bright-field
scanning transmission electron microscopy (ABF-STEM) and
STEM electron energy loss spectroscopy (STEM-EELS). The photocatalytic activity was tested for the decomposition of organic
dyes rhodamine 6G and methyl orange and compared to a reference TiO2 anatase sample. The anatase nanosheets with exposed
{001} facets revealed a significantly lower photocatalytic activity compared to the reference. In order to understand the
mechanism for the catalytic performance, and to investigate the role of the presence of F−, light-induced electron paramagnetic
resonance (EPR) experiments were performed. The EPR results are in agreement with TEM, proving the presence of Ti3+
species close to the surface of the sample and allowing the analysis of the photoinduced formation of paramagnetic species.
Further, ab initio calculations of the anisotropic effective mass of electrons and electron holes in anatase show a very high effective
mass of electrons in the [001] direction, having a negative impact on the mobility of electrons toward the {001} surface and thus
the photocatalysis. Finally, motivated by the experimental results that indicate the presence of fluorine atoms at the surface, we
performed ab initio calculations to determine the position of the band edges in anatase slabs with different terminations of the
{001} surface. The presence of fluorine atoms near the surface is shown to strongly shift down the band edges, which indicates
another reason why it can be expected that the prepared samples with a large amount of {001} surface, but with fluorine atoms
near the surface, show only a low photocatalytic activity.

Abstract: We use microkinetic modeling to examine the potential of plasma-catalytic partial oxidation (POX) of CH4 as a promising new approach to produce oxygenates. We study how different plasma species affect POX of CH4 on the Pt(111) surface, and we discuss the associated kinetic and mechanistic changes. We discuss the effect of vibrationally excited CH4 and O2, as well as plasma-generated radicals and stable intermediates. Our results show that vibrational excitation enhances the turnover frequency (TOF) of catalytic CH4 dissociation and has good potential for improving the selectivities toward CH3OH, HCOOH, and C2 hydrocarbons. Nevertheless, when also considering plasma-generated radicals, we find that these species mainly govern the surface chemistry. Additionally, we find that plasma-generated radicals and stable intermediates enhance the TOFs of COx and oxygenates, increase the selectivity toward oxygenates, and make the formation of HCOOH more significant on Pt(111). We also briefly illustrate the potential impact of Eley−Rideal reactions that involve plasma-generated radicals. Finally, we reveal how various radicals affect the catalyst surface chemistry and we link this to the formation of different products. This allows us to make suggestions on how the plasma composition should be altered to improve the formation of desired products.