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: It remains a real challenge to control the selectivity of the electrocatalytic CO2 reduction (eCO(2)R) reaction to valuable chemicals and fuels. Most of the electrocatalysts are made of non-renewable metal resources, which hampers their large-scale implementation. Here, we report the preparation of bimetallic copper-lead (CuPb) electrocatalysts from industrial metallurgical waste. The metal ions were extracted from the metallurgical waste through simple chemical treatment with ammonium chloride, and CuxPby electrocatalysts with tunable compositions were fabricated through electrodeposition at varying cathodic potentials. X-ray spectroscopy techniques showed that the pristine electrocatalysts consist of Cu-0, Cu1+ and Pb2+ domains, and no evidence for alloy formation was found. We found a volcano-shape relationship between eCO(2)R selectivity toward two electron products, such as CO, and the elemental ratio of Cu and Pb. A maximum Faradaic efficiency towards CO was found for Cu9.00Pb1.00, which was four times higher than that of pure Cu, under the same electrocatalytic conditions. In situ Raman spectroscopy revealed that the optimal amount of Pb effectively improved the reducibility of the pristine Cu1+ and Pb2+ domains to metallic Cu and Pb, which boosted the selectivity towards CO by synergistic effects. This work provides a framework of thinking to design and tune the selectivity of bimetallic electrocatalysts for CO2 reduction through valorization of metallurgical waste.

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: Group II-VI semiconductor honeycomb monolayers have a noncentrosymmetric crystal structure and therefore are expected to be important for nano piezoelectric device applications. This motivated us to perform first principles calculations based on density functional theory to unveil the piezoelectric properties (i.e., piezoelectric stress (e(11)) and piezoelectric strain (d(11)) coefficients) of these monolayer materials with chemical formula MX (where M = Be, Mg, Ca, Sr, Ba, Zr, Cd and X = S, Se, Te). We found that these two-dimensional materials have peculiar piezoelectric properties with d(11) coefficients 1 order of magnitude larger than those of commercially utilized bulk materials. A clear trend in their piezoelectric properties emerges, which

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: 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 study the oxygen thermal stability of two binary
systems. The larger particles are magnetic amorphous Co (7.2 nm) or
Fe3O4 (7.5 nm) nanocrystals, whereas the smaller ones (3.7 nm) are
Au nanocrystals. The nanocrystal ordering as well as the choice of the
magnetic nanoparticles very much influence the stability of the binary
system. A perfect crystalline structure is obtained with the Fe3O4/Au
binary supracrystals. For the Co/Au binary system, oxidation of Co
results in the chemical transformation from Co to CoO, where the size
of the amorphous Co nanoparticles increases from 7.2 to 9.8 nm in
diameter. During the volume expansion of the Co nanoparticles, Au
nanoparticles within the binary assemblies coalesce and are at the
origin of the instability of the binary nanoparticle supracrystals. On the
other hand, for the Fe3O4/Au binary system, the oxidation of Fe3O4 to
γ-Fe2O3 does not lead to a size change of the nanoparticles, which
maintains the stability of the binary nanoparticle supracrystals. A similar behavior is observed for an AlB2-type Co−Ag binary
system: The crystalline structure is maintained, whereas in disordered assemblies, coalescence of Ag nanocrystals is observed.

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: 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: Two-dimensional polyaniline with structural unit C3N is an indirect semiconductor with 0.4 eV band gap, which has attracted a lot of interest because of its unusual electronic, optoelectronic, thermal, and mechanical properties useful for various applications. Adsorption of adatoms is an effective method to improve and tune the properties of C3N. Using first-principles calculations, we investigated the adsorption of adatoms, including H, O, S, F, Cl, B, C, Si, N, P, Al, Li, Na, K, Be, Mg, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, and Zn, on C3N. Depending on the adatom size and the number of valence electrons, they may induce metallic, half-metallic, semiconducting, and ferromagnetic-metallic behavior. In addition, we investigate the effects of an electrical field, charging, and strain on C3N and found how the electronic and magnetic properties are modified. Semi- and full hydrogenation are studied. From the mechanical and thermal stability of C3N monolayer, we found it to be a hard material that can withstand large strain. From our calculations, we gained novel insights into the properties of C3N demonstrating its unique electronic and magnetic properties that can be useful for semiconducting, nanosensor, and catalytic applications.

Abstract: Ultrathin two-dimensional Janus-type platinum dichalcogenide crystals formed by two different atoms at opposite surfaces are investigated by performing state-of-the-art density functional theory calculations. First, it is shown that single-layer PtX2 structures (where X = S, Se, or Te) crystallize into the dynamically stable IT phase and are indirect band gap semiconductors. It is also found that the substitutional chalcogen doping in all PtX2 structures is favorable via replacement of surface atoms with a smaller chalcogen atom, and such a process leads to the formation of Janus-type platinum dichalcogenides (XPtY, where X and Y stand for S, Se, or Te) which are novel single-layer crystals. While all Janus structures are indirect band gap semiconductors as their binary analogues, their Raman spectra show distinctive features that stem from the broken out-of-plane symmetry. In addition, it is revealed that the construction of Janus crystals enhances the piezoelectric constants of PtX2 crystals significantly both in the in plane and in the out-of-plane directions. Moreover, it is shown that vertically stacked van der Waals heterostructures of binary and ternary (Janus) platinum dichalcogenides offer a wide range of electronic features by forming bilayer heterojunctions of type-I, type-II, and type-III, respectively. Our findings reveal that Janus-type ultrathin platinum dichalcogenide crystals are quite promising materials for optoelectronic device applications.

Abstract: We present a first-principles investigation on the stability, electronic structure, and mechanical response of ultrathin heterostructures composed of single layers of InSe and SiGe. First, by performing total energy optimization and phonon calculations, we show that single layers of InSe and SiGe can form dynamically stable heterostructures in 12 different stacking types. Valence and conduction band edges of the heterobilayers form a type-I heterojunction having a tiny band gap ranging between 0.09 and 0.48 eV. Calculations on elastic-stiffness tensor reveal that two mechanically soft single layers form a heterostructure which is stiffer than the constituent layers because of relatively strong interlayer interaction. Moreover, phonon analysis shows that the bilayer heterostructure has highly Raman active modes at 205.3 and 43.7 cm(-1), stemming from the out-of-plane interlayer mode and layer breathing mode, respectively. Our results show that, as a stable type-I heterojunction, ultrathin heterobilayer of InSe/SiGe holds promise for nanoscale device applications.

Abstract: The synthesis of gold nanorods requires the presence of symmetry-breaking and shape-directing additives, among which bromide ions and quaternary ammonium surfactants have been reported as essential. As a result, hexadecyltrimethylammonium bromide (CTAB) has been selected as the most efficient surfactant to direct anisotropic growth. One of the difficulties arising from this selection is the low solubility of CTAB in water at room temperature, and therefore the seeded growth of gold nanorods is usually performed at 25 degrees C or above, which has restricted so far the analysis of kinetic effects derived from lower temperatures. We report a systematic study of the synthesis of gold nanorods from pentatwinned seeds using hexadecyltrimethylammonium chloride (CTAC) as the principal surfactant and a low concentration of bromide as shape-directing agent. Under these conditions, the synthesis can be performed at temperatures as low as 8 degrees C, and the corresponding kinetic effects can be studied, resulting in temperature-controlled aspect ratio tunability.

Abstract: The reduction of dimensionality is a very effective way to achieve appealing properties in two-dimensional materials (2DMs). First-principles calculations can greatly facilitate the prediction of 2DM properties and find possible approaches to enhance their performance. We employed first-principles calculations to gain insight into the impact of different types of point defects (vacancies and substitutional dopants) on the electronic and magnetic properties of a ZnSnN2 (ZSN) monolayer. We show that Zn, Sn, and N + Zn vacancy-defected structures are p-type conducting, while the defected ZSN with a N vacancy is n-type conducting. For substitutional dopants, we found that all doped structures are thermally and energetically stable. The most stable structure is found to be B-doping at the Zn site. The highest work function value (5.0 eV) has been obtained for Be substitution at the Sn site. Li-doping (at the Zn site) and Be-doping (at the Sn site) are p-type conducting, while B-doping (at the Zn site) is n-type conducting. We found that the considered ZSN monolayer-based structures with point defects are magnetic, except those with the N vacancy defects and Be-doped structures. The ab initio molecular dynamics simulations confirm that all substitutionally doped and defected structures are thermally stable. Thus, our results highlight the possibility of tuning the magnetism in ZnSnN2 monolayers through defect engineering.

Abstract: Direct experimental measurement of thermal expansion coefficient without substrate effects is a challenging task for two-dimensional (2D) materials, and its accurate estimation with large-scale ab initio molecular dynamics is computationally very expensive. Machine learning-based interatomic potentials trained with ab initio data have been successfully used in molecular dynamics simulations to decrease the computational cost without compromising the accuracy. In this study, we investigated using Gaussian approximation potentials to reproduce the density functional theory-level accuracy for graphene within both lattice dynamical and molecular dynamical methods, and to extend their applicability to larger length and time scales. Two such potentials are considered, GAP17 and GAP20. GAP17, which was trained with pristine graphene structures, is found to give closer results to density functional theory calculations at different scales. Further vibrational and structural analyses verify that the same conclusions can be deduced with density functional theory level in terms of the reasoning of the thermal expansion behavior, and the negative thermal expansion behavior is associated with long-range out-of-plane phonon vibrations. Thus, it is argued that the enabled larger system sizes by machine learning potentials may even enhance the accuracy compared to small-size-limited ab initio molecular dynamics.