Theoretical study of structure sensitivity on ceria-supported single platinum atoms and its influence on carbon monoxide adsorption.

Density functional theory (DFT) calculations explore the stability of a single platinum atom on various flat, stepped, and defective ceria surfaces, in the context of single-atom catalysts (SACs) for the water-gas shift (WGS) reaction. The adsorption properties and diffusion kinetics of the metal strongly depend on the support termination with large stability on metastable and stepped CeO2(100) and (210) surfaces where the diffusion of the platinum atom is hindered. At the opposite, the more stable CeO2(111) and (110) terminations weakly bind the platinum atom and can promote the growth of metallic clusters thanks to fast diffusion kinetics. The adsorption of carbon monoxide on the single platinum atom supported on the various ceria terminations is also sensitive to the surface structure. Carbon monoxide weakly binds to the single platinum atom supported on reduced CeO2(111) and (211) terminations. The desorption of the CO2 formed during the WGS reaction is thus facilitated on the latter terminations. A vibrational analysis underlines the significant changes in the calculated scaled anharmonic CO stretching frequency on these catalysts.

Identifying the Structure of Supported Metal Catalysts Using Vibrational Fingerprints from Ab Initio Nanoscale Models.

Identifying active sites of supported noble metal nanocatalysts remains challenging, since their size and shape undergo changes depending on the support, temperature, and gas mixture composition. Herein, the anharmonic infrared spectrum of adsorbed CO is simulated using density functional theory (DFT) to gain insight into the nature of Pd nanoparticles (NPs) supported on ceria. The authors systematically determine how the simulated infrared spectra are affected by CO coverage, NP size (0.5-1.5 nm), NP morphology (octahedral, icosahedral), and metal-support contact angle, by exploring a diversity of realistic models inspired by ab initio molecular dynamics. The simulated spectra are then used as a spectroscopic fingerprint to characterize nanoparticles in a real catalyst, by comparison with in-situ diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS) experiments. Truncated octahedral NPs with an acute Pd-ceria angle reproduce most of the measurements. In particular, the authors isolate features characteristic of CO adsorbed at the metal-support interface appearing at low frequencies, both seen in simulation and experiment. This work illustrates the strong need for realistic models to provide a robust description of the active sites, especially at the interface of supported metal nanocatalysts, which can be highly dynamic and evolve considerably during reaction.

Coordination of Ethylamine on Small Silver Clusters: Structural and Topological (ELF, QTAIM) Analyses.

Amine ligands are expected to drive the organization of metallic centers as well as the chemical reactivity of silver clusters early growing during the very first steps of the synthesis of silver nanoparticles via an organometallic route. Density functional theory (DFT) computational studies have been performed to characterize the structure, the atomic charge distribution, and the planar two-dimensional (2D)/three-dimensional (3D) relative stability of small-size silver clusters (Agn, 2 ≤ n ≤ 7), with or without an ethylamine (EA) ligand coordinated to the Ag clusters. The transition from 2D to 3D structures is shifted from n = 7 to 6 in the presence of one EA coordinating ligand, and it is explained from the analysis of the Ag-N and Ag-Ag bond energies. For fully EA saturated silver clusters (Agn-EAn), the effect on the 2D/3D transition is even more pronounced with a shift between n = 4 and 5. Subsequent electron localization function (ELF) and quantum theory of atoms in molecules (QTAIM) topological analyses allow for the fine characterization of the dative Ag-N and metallic Ag-Ag bonds, both in nature and in strength. Electron transfer from ethylamine to the coordinated silver atoms induces an increase of the polarization of the metallic core.

Revealing Size Dependent Structural Transitions in Supported Gold Nanoparticles in Hydrogen at Atmospheric Pressure.

The enhancement of the catalytic activity of gold nanoparticles with their decreasing size is often attributed to the increasing proportion of low-coordinated surface sites. This correlation is based on the paradigmatic picture of working gold nanoparticles as perfect crystal forms having complete and static outer surface layers whatever their size. This picture is incomplete as catalysts can dynamically change their structure according to the reaction conditions and as such changes can be eventually size-dependent. In this work, using aberration-corrected environmental electron microscopy, size-dependent crystal structure and morphological evolution in gold nanoparticles exposed to hydrogen at atmospheric pressure, with loss of the face-centered cubic crystal structure of gold for particle size below 4 nm, are revealed for the first time. Theoretical calculations highlight the role of mobile gold atoms in the observed symmetry changes and particle reshaping in the critical size regime. An unprecedented stable surface molecular structure of hydrogenated gold decorating a highly distorted core is identified. By combining atomic scale in situ observations and modeling of nanoparticle structure under relevant reaction conditions, this work provides a fundamental understanding of the size-dependent reactivity of gold nanoparticles with a precise picture of their surface at working conditions.

Capturing Solvation Effects at a Liquid/Nanoparticle Interface by Ab Initio Molecular Dynamics: Pt201 Immersed in Water.

Solvation can substantially modify the adsorption properties of heterogeneous catalysts. Although essential for achieving realistic theoretical models, assessing such solvent effects over nanoparticles is challenging from a computational standpoint due to the complexity of those liquid/metal interfaces. This effect is investigated by ab initio molecular dynamics simulations at 350 K of a large platinum nanoparticle immersed in liquid water. The first solvation layer contains twice as much physisorbed water molecules above the terraces, than chemisorbed ones located only at edges and corners. The solvent stabilizes the binding energy of chemisorbates: 66% of the total gain comes from interactions with physisorbed molecules and 34% from the influence of bulk liquid.

Coverage-dependent thermodynamic analysis of the formation of water and hydrogen peroxide on a platinum model catalyst.

Understanding the selectivity of the oxygen reduction reaction, especially the formation of water versus hydrogen peroxide in fuel cells, is an ongoing challenge in electrochemistry, surface science and catalysis. In this study, we propose a comprehensive thermodynamic analysis of the reaction intermediates for the formation of water on Pt(111). Density functional theory calculations of all the elementary steps linking hydroxyl and hydroperoxyl surface species with water and hydrogen peroxide have been performed at low (1/12 ML, ML = monolayer) and high (1/4 ML) coverages. The reaction energy variation for the two competing elementary events (molecular oxygen dissociation and hydroperoxyl formation) is strongly coverage-dependent. For the direct dissociation, an increase is observed at low coverage with respect to the usual high coverage picture. The stability of the reaction intermediates is investigated from thermodynamic diagrams. At 353 K and a total pressure of 1 atm, water and hydroxyl surface species are expected to compete for adsorption on Pt(111).

Fast prediction of adsorption properties for platinum nanocatalysts with generalized coordination numbers.

Platinum is a prominent catalyst for a multiplicity of reactions because of its high activity and stability. As Pt nanoparticles are normally used to maximize catalyst utilization and to minimize catalyst loading, it is important to rationalize and predict catalytic activity trends in nanoparticles in simple terms, while being able to compare these trends with those of extended surfaces. The trends in the adsorption energies of small oxygen- and hydrogen-containing adsorbates on Pt nanoparticles of various sizes and on extended surfaces were analyzed through DFT calculations by making use of the generalized coordination numbers of the surface sites. This simple and predictive descriptor links the geometric arrangement of a surface to its adsorption properties. It generates linear adsorption-energy trends, captures finite-size effects, and provides more accurate descriptions than d-band centers and usual coordination numbers. Unlike electronic-structure descriptors, which require knowledge of the densities of states, it is calculated manually. Finally, it was shown that an approximate equivalence exists between generalized coordination numbers and d-band centers.

Influence of air exposure on structural isomers of silver nanoparticles.

Up to date, the influence of ambient air exposure on the energetics and stability of silver clusters has rarely been investigated and compared to clusters in vacuum. Silver clusters up to 3000 atoms in size, on an amorphous carbon film, have been exposed to ambient air and investigated by atomic-resolution imaging in the aberration-corrected Scanning Transmission Electron Microscope. Ordered structures comprise more than half the population, the rest are amorphous. Here, we show that the most common ordered isomer structures is the icosahedron. These results contrast with the published behaviour of silver clusters protected from atmospheric exposure, where the predominant ordered isomer is face-centred cubic. We propose that the formation of surface oxide or sulphide species resulting from air exposure can account for this deviation in stable isomer. This interpretation is consistent with density functional theory calculations based on silver nanoclusters, in the size range 147-201 atoms, on which methanethiol molecules are adsorbed. An understanding of the effects of ambient exposure on the atomic structure and therefore functional properties of nanoparticles is highly relevant to their real-world performance and applications.

How far the chemistry of self-assembled monolayers on gold surfaces affects their work function?

Self-assembled monolayers composed of various long-chain aliphatic molecules and different tail functional groups have been synthesized on the Au(111) surface and characterized by Kelvin probe force microscopy and ultraviolet photoelectron spectroscopy. Carboxy, amino, thio and methyl terminal groups have been considered in the design of self-assembled monolayers with different aliphatic chain lengths (from C6 to C16). Work function measurements by Kelvin probe force microscopy have been carried out under a controlled and room atmosphere. Remarkably, a reduction of the relative humidity from 40% to 3% has induced a work function shift of up to 0.3 eV. As expected, the changes of the chain length of the aliphatic moiety and of the tail group have a significant impact on the tuning of the measured work function (3.90 eV for dodecanethiol versus 4.57 eV for mercaptohexadecylamine). Surprisingly, the change of the net dipole moment of the tail group (sign and amplitude) does not dominate the work function variations. In contrast, the change of the chain length and the possibility of the tail group to form a complex hydrogen bond network between molecules lead to significant modulations of the work function. In order to interpret these original findings, density functional theory models of equivalent self-assembled monolayers adsorbed on the Au(111) surface have been developed at an unprecedented level of description with large supercells including simultaneously 27 co-adsorbed molecules and weak van der Waals interactions between them. Such large systems have allowed the theoretical modeling of complex hydrogen bond networks between molecules when possible (carboxy tail group). The comparison between computed and measured work functions shows a striking agreement, thus allowing the disentanglement of the previously mentioned competing effects. This consistency between experiment and theory will help in designing the electronic properties of self-assembled monolayers in the context of molecular electronics and organic transistors.

Growth of a Pt film on non-reduced ceria: a density functional theory study.

The growth of platinum on non-reduced CeO(2) (111) surface is studied by means of calculations based on the density functional theory. Particles of increasing size are formed on the oxide surface by incorporating the platinum atoms one by one until multilayer films are obtained. The main conclusion is that platinum atoms tend to maximize the number of metallic bonds and to approach the situation of the bulk, hence preferring films to particles, particles to isolated atoms, and a three-dimensional growth to a two-dimensional one. The supported particles and the films exhibit a contraction of the Pt-Pt distances, with respect to those of the Pt bulk, in order to match the ceria lattice. The density of states projected on the film surface platinum atoms shows important differences in shape and energy (lower d-band center) compared to the Pt(111) reference surface, which could be the major reason for the observed changes in catalytic reactivity when deposited particles are compared with single crystal surfaces.

Nature of adhesion of condensed organic films on platinum by first-principles simulations.

Understanding the nature of the adhesion of an organic liquid on a metal surface is of paramount importance for elucidating the stability and chemical reactivity at these complex interfaces. However, to date, the morphology, layering and chemical properties at organic liquid metal interfaces have been rarely known. Using semi-empirical dispersion corrected density functional theory calculations and ab initio molecular dynamics simulations, we show that carbon tetrachloride and ethanol films confined to a platinum surface alter their intrinsic properties and exhibit interfacial reactivity. A few interface carbon tetrachloride (ethanol) molecules adsorb dissociatively (molecularly) on platinum thanks to the surrounding medium. The adsorption strength of the interfacial molecules is consequently increased in the condensed phase as compared to the gas phase. This remarkable effect is rationalized by an interaction energy decomposition model and an electrostatic potential analysis.

Importance of Defective and Nonsymmetric Structures in Silver Nanoparticles.

Scanning transmission electron microscopy experiments indicate that face-centered cubic (FCC) is the predominant ordered structure for Ag309 ± 7 nanoclusters, synthesized in vacuum. Historically, experiments do not present a consensus on the morphology at these sizes, whereas theoretical studies find the icosahedral symmetry for Ag309 and the decahedral shape for nearby sizes. We employ density functional theory calculations to rationalize these observations, considering both regular and defective Ag nanoparticles (281-321 atoms). The change of stability induced by the presence of defects, symmetry loss, and change of number of atoms is evaluated by the nanoparticle surface energy, which was measured previously. FCC and decahedral symmetries are found to be more favorable than icosahedral, consistent with our measurements of clusters protected from extended atmospheric exposure. In addition, an energy-free descriptor, surface atomic density, is proposed and qualitatively reproduces the surface energy data. Nonsymmetric and defective structures may be preferred over perfectly regular ones within a given size range.

Effect of the Ligand Binding Strength on the Morphology of Functionalized Gold Nanoparticles.

Functionalized gold nanoparticles are investigated by density functional theory calculations in the context of cancer radiotherapy. Several typical experimental shapes, including nanostars, nanospheres, and nanorods, are modeled by optimizing Au clusters covered by organic monolayers composed of hydrated short-chain polyethylene glycol (PEG) ligands. The PEGylation stabilizes significantly the stellation of decahedral Au54 by deforming significantly its geometry at the spikes. The higher stability of the PEG molecules adsorbed on this stellated nanocluster with respect to the more spherical icosahedral Au55 and truncated octahedral Au79 leads to a larger energy cost to desorb them and thus a weaker propensity for the starred nanoparticle to exchange ligands with the cell membrane, in agreement with experiments. These results open interesting possibilities for advancing our understanding of the cellular uptake of gold nanoparticles.

Restructuring of the Pt3Sn(111) surfaces induced by atomic and molecular oxygen from first principles.

The surface restructuring of Pt(3)Sn(111) induced by oxygen chemisorption is examined by means of density-functional theory calculations. Molecular and atomic oxygen chemisorption is investigated on the two available terminations of the bulk alloy--(2 x 2) and (square root(3) x square root(3))R30 degrees--these two surfaces differing by the tin content and the nature of chemical sites. An extensive geometric, energetic, and vibrational analysis is performed including the influence of oxygen coverage in the case of atomic adsorption. For molecular adsorption, regular structures have been obtained for both surfaces with a clear effect of tin on the stability of the adsorption forms. In contrast, for atomic adsorption, two oxygen chemical states are found. In particular, a peculiar surface restructuring, involving the formation of a network of SnO(2) species, appears for large oxygen coverage. However the two terminations present discrepancies for the restructuring mechanism all along the oxygen coverage increase. All these results are supported by a systematic vibrational analysis.

A Water Solvation Shell Can Transform Gold Metastable Nanoparticles in the Fluxional Regime.

Solvated gold nanoparticles have been modeled in the fluxional regime by density functional theory including dispersion forces for an extensive set of conventional morphologies. The study of isolated adsorption of one water molecule shows that the most stable adsorption forms are similar (corners and edges) regardless of the nanoparticle shape and size, although the adsorption strength differs significantly (0.15 eV). When a complete and explicit water solvation shell interacts with gold nanoclusters, metastable in vacuum and presenting a predominance of (100) square facets (ino-decahedra Au55 and Au147), these nanoparticles are found unstable and transform into the closest morphologies exhibiting mainly (111) triangular facets and symmetries. The corresponding adsorption strength per water molecule becomes independent of shape and size and is enhanced by the formation of two hydrogen bonds on average. For applications in radiotherapy, this study suggests that the shapes of small gold nanoparticles should be homogenized by interacting with the biological environment.

Why conclusions from platinum model surfaces do not necessarily lead to enhanced nanoparticle catalysts for the oxygen reduction reaction.

Experiments on model surfaces commonly help in identifying the structural sensitivity of catalytic reactions. Nevertheless, their conclusions do not frequently lead to devising superior "real-world" catalysts. For instance, this is true for single-crystal platinum electrodes and the oxygen reduction reaction (ORR), an important reaction for sustainable energy conversion. Pt(111) is substantially enhanced by steps, reaching a maximum at short terrace lengths of 3-4 atoms. Conversely, regular platinum nanoparticles with similar undercoordinated defects are less active than Pt(111) and their activity increases alongside the terrace-to-defect ratio. We show here that a model to design ORR active sites on extended surfaces can also be used to solve this apparent contradiction and provide accurate design rules for nanoparticles. Essentially, only surfaces and nanostructures with concave defects can surpass the activity of Pt(111), whereas convex defects are inactive. Importantly, only the latter are present in regular nanoparticles, which is why we design various concave nanoparticles with high activities.

Promoter effect of hydration on the nucleation of nanoparticles: direct observation for gold and copper on rutile TiO2 (110).

Direct observation of the promoting effect of hydration on the nucleation of gold and copper nanoparticles supported on partially reduced rutile TiO2 (110) is achieved by combined scanning tunneling microscopy experiments and density functional theory calculations. The experiments show a clear difference between the two metals. Gold nanoparticles grow at the vicinity of the surface hydroxyl domains, whereas the nucleation of copper is not substantially affected by hydration. The nucleation of gold on surface oxygen vacancies is observed although this is not the only preferential site. Theoretical calculations of the coadsorbed phases of gold, copper and hydroxyl species on stoichiometric and reduced TiO2 (110) surfaces under relevant conditions of temperature and pressure support the experimental interpretation. Surface hydration tends to stabilize significantly gold adsorption on the stoichiometric support, while its influence on copper adsorption is not pronounced. The theoretical analysis shows that the early stages of the nucleation on hydrated stoichiometric surfaces correspond to mono-hydroxylated metallic species co-chemisorbed with hydroxyl species, whereas those on hydrated reduced surfaces are metallic atoms bound to oxygen vacancies and weakly perturbed by surface hydration.

Formation of a high coverage (3 x 3) NO phase on Pd(111) at elevated pressures: interplay between kinetic and thermodynamic accessibility.

Using in situ polarization modulation infrared reflection absorption spectroscopy and density functional theory calculations, a new high-coverage monomeric NO adsorption state on Pd(111) was observed and proposed to have a (3 x 3)-7NO structure. Formation of this high coverage NO phase was found to take place only at elevated pressure and temperature conditions showing that some of the accessible thermodynamic equilibrium states at elevated temperatures and pressures are thermodynamically unfavorable or kinetically hindered at lower temperatures and pressures. Our results emphasize the danger of extrapolating results from traditional surface science experiments performed under ultrahigh vacuum to elevated temperature and pressure conditions encountered in heterogeneous catalysis.

Introducing structural sensitivity into adsorption-energy scaling relations by means of coordination numbers.

The search for improved heterogeneous catalysts is an important but difficult task. Scaling relations between the adsorption energies of reaction intermediates greatly facilitate the computational design of catalysts. However, this methodology does not currently incorporate structure sensitivity and hence cannot describe adequately the overall activity of realistic catalyst particles and extended surfaces with several facets, edges and apices. Here, we generalize scaling relations by examining twelve different low-index, stepped and kinked surfaces of nine transition metals. This allows us to quantify the effect of the adsorption-site geometry on these relations, ensures a full prediction of their parameters, and helps in identifying intrinsic thermodynamic restrictions to the performance of catalysts. The resulting fully predictable, structure-sensitive scaling relations are a step towards the long-sought rational design of multifaceted catalytic particles. Such a design can now target not only the chemical nature of active materials but also the actual geometry of their active sites.

Finding optimal surface sites on heterogeneous catalysts by counting nearest neighbors.

A good heterogeneous catalyst for a given chemical reaction very often has only one specific type of surface site that is catalytically active. Widespread methodologies such as Sabatier-type activity plots determine optimal adsorption energies to maximize catalytic activity, but these are difficult to use as guidelines to devise new catalysts. We introduce "coordination-activity plots" that predict the geometric structure of optimal active sites. The method is illustrated on the oxygen reduction reaction catalyzed by platinum. Sites with the same number of first-nearest neighbors as (111) terraces but with an increased number of second-nearest neighbors are predicted to have superior catalytic activity. We used this rationale to create highly active sites on platinum (111), without alloying and using three different affordable experimental methods.