The unexpected effect of vacancies and wrinkling on the electronic properties of MoS2 layers.

We report a combined experimental/theoretical approach to study the connection of S-vacancies and wrinkling on MoS2 layers, and how this feature produces significant changes in the electronic structure and reactivity of this 2D material. The MoS2 material, when used as a catalyst in operative conditions, was found to be mainly composed of thin and short 1-5 layer sheets instead of a poorly crystalline structure, as it was previously assumed. Notably wrinkled structures with S-vacancies were also found through transmission electron microscopy. Atomistic simulations revealed a natural connection between sulfur-vacancies, wrinkling and folding. Density functional calculations further revealed that such curved structures present a lower electronic band-gap and a higher reactivity towards thiophene compared to the planar MoS2 counterpart.

Interaction of Water with the Gypsum (010) Surface: Structure and Dynamics from Nonlinear Vibrational Spectroscopy and Ab Initio Molecular Dynamics.

Water-mineral interfaces are important for several environmental, industrial, biological, and geological processes. Gypsum, CaSO4·2H2O, is a widespread mineral of high technological, medical, and environmental relevance, but little is known about its surface structure and its interaction with water. A molecular-level understanding of gypsum/water interface is given here by a combined experimental/theoretical study. We investigate the structure and dynamics of water adsorbed from vapor on the gypsum (010) single-crystal surface at room temperature, combining sum-frequency generation (SFG) vibrational spectroscopy experiments and ab initio molecular dynamics (AIMD) simulations. The SFG spectra of gypsum at low relative humidity (RH) show an anisotropic arrangement of structural water molecules and the presence of dangling OH groups. The AIMD simulations allow a detailed assignment of the SFG spectra and show that the cleaved (010) surface rearranges to have only 25% of the OH groups pointing away from the surface. At higher RHs, the first adsorbed water layer binds to these OH groups and forms an anisotropic arrangement, but with the amount of free OH groups significantly suppressed and without any significant diffusion. Upon adsorption of a second water layer, although the topmost layer of molecules is more disordered and dynamic than the previous one, its structure is still influenced by the gypsum surface underneath because it has a much reduced amount of free OH groups with respect to the free surface of water, and a slower surface diffusion with respect to bulk water. The theoretical results corroborate the experimental ones and provide an accurate atomic characterization of the surface structure.

Bimetallic Ag-Pt Sub-nanometer Supported Clusters as Highly Efficient and Robust Oxidation Catalysts.

A combined experimental and theoretical investigation of Ag-Pt sub-nanometer clusters as heterogeneous catalysts in the CO→CO2 reaction (COox) is presented. Ag9 Pt2 and Ag9 Pt3 clusters are size-selected in the gas phase, deposited on an ultrathin amorphous alumina support, and tested as catalysts experimentally under realistic conditions and by first-principles simulations at realistic coverage. In situ GISAXS/TPRx demonstrates that the clusters do not sinter or deactivate even after prolonged exposure to reactants at high temperature, and present comparable, extremely high COox catalytic efficiency. Such high activity and stability are ascribed to a synergic role of Ag and Pt in ultranano-aggregates, in which Pt anchors the clusters to the support and binds and activates two CO molecules, while Ag binds and activates O2 , and Ag/Pt surface proximity disfavors poisoning by CO or oxidized species.

Reactivity of atomically dispersed Pt(2+) species towards H2: model Pt-CeO2 fuel cell catalyst.

The reactivity of atomically dispersed Pt(2+) species on the surface of nanostructured CeO2 films and the mechanism of H2 activation on these sites have been investigated by means of synchrotron radiation photoelectron spectroscopy and resonant photoemission spectroscopy in combination with density functional calculations. Isolated Pt(2+) sites are found to be inactive towards H2 dissociation due to high activation energy required for H-H bond scission. Trace amounts of metallic Pt are necessary to initiate H2 dissociation on Pt-CeO2 films. H2 dissociation triggers the reduction of Ce(4+) cations which, in turn, is coupled with the reduction of Pt(2+) species. The mechanism of Pt(2+) reduction involves reverse oxygen spillover and formation of oxygen vacancies on Pt-CeO2 films. Our calculations suggest the existence of a threshold concentration of oxygen vacancies associated with the onset of Pt(2+) reduction.

Ligand/cluster/support catalytic complexes in heterogeneous ultrananocatalysis: NO oxidation on Ag3/MgO(100).

In the present work we explore via first-principles simulations whether the ligand/cluster/support catalytic complex generated by CO oxidation over silver trimers deposited on the regular MgO(100) surface - i.e. a Ag3/carbonate or Ag3(CO3)/MgO(100) species - can be used as a catalyst in a different reaction: the selective oxidation of NO to NO2 (or NOox). The Ag3(CO3)/MgO(100) complex is first shown to be reasonably stable at room temperature in terms of both disaggregation and sintering, and that it can be generated from Ag3 adsorbed onto an oxygen vacancy defect of the regular MgO(100) surface under oxidation conditions. It is then found that the Ag3(CO3)/MgO(100) species transforms under NOox conditions into an even more complex aggregate, a mixed carbonate/double-nitrite Ag3(CO3)(NO2)2/MgO(100) species, which can then act as an efficient catalyst of NOox. It is noteworthy that under NOox reaction conditions a different ligand/cluster/support catalytic complex is formed with respect to the original COox one. These findings prove the diversity of the catalytic chemistry of subnanometer (or ultranano) metal clusters deposited on oxide substrates, associated with the formation of many different ligand/cluster/support aggregates, the vast amount of combinatorial possibilities thus opening, and the need for computational approaches to perform systematic structural and stoichiometric searches in order to cope with such a multiform diversity.

Communication: Striking dependence of diffusion kinetics in Ag-Cu nanoalloys upon composition and quantum effects.

The kinetics of elemental inter-diffusion in Ag-Cu nanoalloys of 32-34 atoms around 80:20 composition is theoretically investigated by combining analytic-potential and first-principles calculations. An extremely varied behavior is found, with transformation times ranging from tens of nanoseconds to weeks at room temperature in a narrow interval of size and composition, also depending on quantum effects in magic clusters. Predictions are consistent with time-of-flight experiments and suggest their interpretation in a new light.

Adsorption-induced restructuring and early stages of carbon-nanotube growth on Ni nanoparticles.

Carbon adsorption on various Ni surfaces is investigated as a function of coverage via a combination of first-principles simulations and field emission microscope experiments. It is found that carbon can be efficiently stored as subsurface carbides, but with different energetics on differently oriented surfaces depending on their compactness and density of adsorption sites. In the resulting morphological reshaping, {113} facets are predicted to grow at the expense of {111} and {100} facets, in excellent agreement with experimental observations. Moreover, at high coverage on the {113} surface the carbon adsorption energy passes through a maximum after which a structural crossover is realized such that carbon atoms tend to ascend to the surface to form one-dimensional chains (which are the precursors of graphitic nanostructures). This rationalizes the experimental observation of an incubation time between carbon storage and the beginning of catalytic growth, and provides insight into the early stages (nucleation mechanism) of carbon nanotubes on Ni nanoparticles.

Modelling the metal-on-top effect for Pd clusters on the MgO{100} substrate.

We introduce a novel empirical model for the adhesion of Pd clusters on the MgO{100} substrate. The new model corrects the known bias of previous models toward structures with large interfaces with the substrate due to the failure to account for the so-called "metal-on-top" effect, i.e., the enhancement of the adhesion due to the presence of other metal atoms on top of those which are directly in contact with the substrate. The new model is parametrised using density-functional theory calculations on MgO-supported Pd clusters with sizes up to 80 atoms. The proposed potential is continuous with respect to spatial coordinates and can therefore be used directly in molecular dynamics simulations.

Ordered arrays of size-selected oxide nanoparticles.

A bottom-up approach to produce a long-range ordered superlattice of monodisperse and isomorphic metal-oxide nanoparticles (NP) supported onto an oxide substrate is demonstrated. The synthetic strategy consists of self-assembling metallic NP on an ultrathin nanopatterned aluminum oxide template followed by a morphology-conserving oxidation process, and is exemplified in the case of Ni, but is generally applicable to a wide range of metallic systems. Both fully oxidized and core-shell metal-metal-oxide particles are synthesized, up to 3-4 nm in diameter, and characterized via spectroscopic and theoretical tools. This opens up a new avenue for probing unit and ensemble effects on the properties of oxide materials in the nanoscale regime.

Interface effects on the magnetism of CoPt-supported nanostructures.

The magnetism of CoPt nanostructures supported on the MgO(100) surface is investigated via first-principles simulations using 1D models. Nanostructures with L1(0) chemical ordering and cube-on-cube epitaxy are predicted to possess large magnetic moments and easy magnetization axis perpendicular to the surface. However, their magnetic anisotropy energy is roughly halved with respect to the bulk alloy due to a peculiar mixing of particle and support electronic states. The general factors at play in determining this behavior and the implications of these findings are discussed in view of designing room-temperature magnetic bits.

Structures of gas-phase Ag-Pd nanoclusters: a computational study.

Gas-phase Ag-Pd clusters in the size range of 38-100 atoms are studied via a combined density-functional/empirical-potential (DF-EP) approach. Many-body EPs describing Pd-Pd, Ag-Ag, and Ag-Pd interactions are reparametrized and used in thorough global optimization searches at sizes N=38, 60, and 100 and compositions 25%, 50%, and 75%. The results are analyzed in terms of structural families, whose lowest-energy isomers are reoptimized at the DF level to investigate the crossover among structural motifs. It is found that the reparametrized EPs show a better qualitative and quantitative agreement with DF results when compared to the original potentials taken from literature: Both methods agree on which is the lowest-energy isomer at each size and composition, and the energy differences in the various isomers are in good qualitative agreement, especially for 60- and 100-atom clusters. The reparametrized potentials should thus be applicable to large clusters, where DF calculations are not feasible any more.

Creating single-atom Pt-ceria catalysts by surface step decoration.

Single-atom catalysts maximize the utilization of supported precious metals by exposing every single metal atom to reactants. To avoid sintering and deactivation at realistic reaction conditions, single metal atoms are stabilized by specific adsorption sites on catalyst substrates. Here we show by combining photoelectron spectroscopy, scanning tunnelling microscopy and density functional theory calculations that Pt single atoms on ceria are stabilized by the most ubiquitous defects on solid surfaces--monoatomic step edges. Pt segregation at steps leads to stable dispersions of single Pt(2+) ions in planar PtO4 moieties incorporating excess O atoms and contributing to oxygen storage capacity of ceria. We experimentally control the step density on our samples, to maximize the coverage of monodispersed Pt(2+) and demonstrate that step engineering and step decoration represent effective strategies for understanding and design of new single-atom catalysts.

The atomistic origin of the extraordinary oxygen reduction activity of Pt3Ni7 fuel cell catalysts.

Recently Debe et al. reported that Pt3Ni7 leads to extraordinary Oxygen Reduction Reaction (ORR) activity. However, several reports show that hardly any Ni remains in the layers of the catalysts close to the surface ("Pt-skin effect"). This paradox that Ni is essential to the high catalytic activity with the peak ORR activity at Pt3Ni7 while little or no Ni remains close to the surface is explained here using large-scale first-principles-based simulations. We make the radical assumption that processing Pt-Ni catalysts under ORR conditions would leach out all Ni accessible to the solvent. To simulate this process we use the ReaxFF reactive force field, starting with random alloy particles ranging from 50% Ni to 90% Ni and containing up to ∼300 000 atoms, deleting the Ni atoms, and equilibrating the resulting structures. We find that the Pt3Ni7 case and a final particle radius around 7.5 nm lead to internal voids in communication with the exterior, doubling the external surface footprint, in fair agreement with experiment. Then we examine the surface character of these nanoporous systems and find that a prominent feature in the surface of the de-alloyed particles is a rhombic structure involving 4 surface atoms which is crystalline-like but under-coordinated. Using density-functional theory, we calculate the energy barriers of ORR steps on Pt nanoporous catalysts, focusing on the Oad-hydration reaction (Oad + H2Oad → OHad + OHad) but including the barriers of O2 dissociation (O2ad → Oad + Oad) and water formation (OHad + Had → H2Oad). We find that the reaction barrier for the Oad-hydration rate-determining-step is reduced significantly on the de-alloyed surface sites compared to Pt(111). Moreover we find that these active sites are prevalent on the surface of particles de-alloyed from a Pt-Ni 30 : 70 initial composition. These simulations explain the peak in surface reactivity at Pt3Ni7, and provide a rational guide to use for further optimization of improved catalytic and nanoporous materials.

Metal tungstates at the ultimate two-dimensional limit: fabrication of a CuWO4 nanophase.

Metal tungstates (with general formula MWO4) are functional materials with a high potential for a diverse set of applications ranging from low-dimensional magnetism to chemical sensing and photoelectrocatalytic water oxidation. For high level applications, nanoscale control of film growth is necessary, as well as a deeper understanding and characterization of materials properties at reduced dimensionality. We succeeded in fabricating and characterizing a two-dimensional (2-D) copper tungstate (CuWO4). For the first time, the atomic structure of an ultrathin ternary oxide is fully unveiled. It corresponds to a CuWO4 monolayer arranged in three sublayers with stacking O-W-O/Cu from the interface. The resulting bidimensional structure forms a robust framework with localized regions of anisotropic flexibility. Electronically it displays a reduced band gap and increased density of states close to the Fermi level with respect to the bulk compound. These unique features open a way for new applications in the field of photo- and electrocatalysis, while the proposed synthesis method represents a radically new and general approach toward the fabrication of 2-D ternary oxides.

Redox processes at a nanostructured interface under strong electric fields.

Manipulation of chemistry and film growth via external electric fields is a longstanding goal in surface science. Numerous systems have been predicted to show such effects but experimental evidence is sparse. Here we demonstrate in a custom-designed UHV apparatus that the application of spatially extended, homogeneous, very high (>1 V nm(-1)) DC-fields not only changes the system energetics but triggers dynamic processes which become important much before static contributions appreciably modify the potential energy landscape. We take a well characterized ultrathin NiO film on a Ag(100) support as a proof-of-principle test case, and show how it gets reduced to supported Ni clusters under fields exceeding the threshold of +0.9 V nm(-1). Using an effective model, we trace the observed interfacial redox process down to a dissociative electron attachment resonant mechanism. The proposed approach can be easily implemented and generally applied to a wide range of interfacial systems, thus opening new opportunities for the manipulation of film growth and reaction processes at solid surfaces under strong external fields.

Computational approaches to the chemical conversion of carbon dioxide.

The conversion of CO2 into fuels and chemicals is viewed as an attractive route for controlling the atmospheric concentration and recycling of this greenhouse gas, but its industrial application is limited by the low selectivity and activity of the current catalysts. Theoretical modeling, in particular density functional theory (DFT) simulations, provides a powerful and effective tool to discover chemical reaction mechanisms and design new catalysts for the chemical conversion of CO2, overcoming the repetitious and time/labor consuming trial-and-error experimental processes. In this article we give a comprehensive survey of recent advances on mechanism determination by DFT calculations for the catalytic hydrogenation of CO2 into CO, CH4, CH3OH, and HCOOH, and CO2 methanation, as well as the photo- and electrochemical reduction of CO2. DFT-guided design procedures of new catalytic systems are also reviewed, and challenges and perspectives in this field are outlined.

A first-principles theoretical approach to heterogeneous nanocatalysis.

A theoretical approach to heterogeneous catalysis by sub-nanometre supported metal clusters and alloys is presented and discussed. Its goal is to perform a computational sampling of the reaction paths in nanocatalysis via a global search in the phase space of structures and stoichiometry combined with filtering which takes into account the given experimental conditions (catalytically relevant temperature and reactant pressure), and corresponds to an incremental exploration of the disconnectivity diagram of the system. The approach is implemented and applied to the study of propylene partial oxidation by Ag(3) supported on MgO(100). First-principles density-functional theory calculations coupled with a Reactive Global Optimization algorithm are performed, finding that: (1) the presence of an oxide support drastically changes the potential energy landscape of the system with respect to the gas phase, favoring configurations which interact positively with the electrostatic field generated by the surface; (2) the reaction energy barriers for the various mechanisms are crucial in the competition between thermodynamically and kinetically favored reaction products; (3) a topological database of structures and saddle points is produced which has general validity and can serve for future studies or for deriving general trends; (4) the MgO(100) surface captures some major features of the effect of an oxide support and appears to be a good model of a simple oxide substrate; (5) strong cooperative effects are found in the co-adsorption of O(2) and other ligands on small metal clusters. The proposed approach appears as a viable route to advance the role of predictive computational science in the field of heterogeneous nanocatalysis.

Work Function of Oxide Ultrathin Films on the Ag(100) Surface.

Theoretical calculations of the work function of monolayer (ML) and bilayer (BL) oxide films on the Ag(100) surface are reported and analyzed as a function of the nature of the oxide for first-row transition metals. The contributions due to charge compression, charge transfer and rumpling are singled out. It is found that the presence of empty d-orbitals in the oxide metal can entail a charge flow from the Ag(100) surface to the oxide film which counteracts the decrease in the work function due to charge compression. This flow can also depend on the thickness of the film and be reduced in passing from ML to BL systems. A regular trend is observed along first-row transition metals, exhibiting a maximum for CuO, in which the charge flow to the oxide is so strong as to reverse the direction of rumpling. A simple protocol to estimate separately the contribution due to charge compression is discussed, and the difference between the work function of the bare metal surface and a Pauling-like electronegativity of the free oxide slabs is used as a descriptor quantity to predict the direction of charge transfer.