A theoretical investigation of 38-atom CuPd clusters: the effect of potential parameterisation on structure and segregation

Understanding the structure of bimetallic clusters is increasingly important due to their emerging practical applications. Herein we investigate the structure of 38-atom CuPd clusters using a genetic algorithm with cluster energies described by the semi-empirical Gupta potential. Selected clusters are then refined with density functional theory. Three different parameterisations of the Gupta potential are used and their performance assessed to understand what features of bulk and surfaces are necessary to capture for accurate description of small clusters. Three general regions of motif stability exist; for the Pd majority clusters (Pd38 to Cu4Pd34) the truncated octahedron is most stable, while for clusters of intermediate compositions (Cu5Pd33 to Cu25Pd13) a "pancake" icosahedron is most stable, and for the Cu majority clusters (Cu26Pd12 to Cu38) again the truncated octahedron is most stable. CuPd clusters tend to segregate to a Cu-core, Pd-shell structure if possible, and at higher Cu compositions, the Pd segregates to the faces of the cluster. Using multiple parameterisations of the Gupta potential ensures the full variety of possible structures is found, and improves the search for the most stable CuPd clusters.

Density functional theory study of thiophene desulfurization and conversion of desulfurization products on the Ni(111) surface and Ni55 cluster: implication for the mechanism of reactive adsorption desulfurization over Ni/ZnO catalysts

Density functional theory calculations were performed to study thiophene desulfurization and conversion of desulfurization products on the Ni(111) surface and Ni₅₅ cluster.

Adsorption and dehydrogenation of C2-C6n-alkanes over a Pt catalyst: a theoretical study on the size effects of alkane molecules and Pt substrates

Adsorption and dehydrogenation of C2-C6n-alkanes are investigated on a Pt substrate using density functional theory (DFT) calculations, and the size effects of alkane molecules and Pt substrates are discussed in detail. The Pt(111) surface and Pt55 cluster are chosen to represent large and small Pt nanoparticles, respectively. The C2-C6 straight-chain alkanes show no site preference on Pt(111) drifting over the surface, but prefer to locate along the edge sites of Pt55. Our results suggest that a linear relationship holds for the adsorption energies of n-alkanes against the chain length on Pt(111), in accordance with the experimental observations. Pt55 also exhibits a similar linear relationship for n-alkanes but with larger adsorption energies due to the low-coordinated Pt atoms at the edge site. For the two-step dehydrogenation from alkanes to alkenes, the first dehydrogenation reaction is the rate-determining step (RDS) on Pt(111), and a larger size of alkane molecule will lead to a lower dehydrogenation activity. While on Pt55, no RDS is present and the dehydrogenation activity oscillates slightly as the chain length of n-alkane increases. Generally, Pt55 involves lower energy barriers for most dehydrogenation steps compared to Pt(111), indicating that small Pt particles with more low-coordinated Pt atoms are more active towards alkane dehydrogenation. In addition, a clear BEP relationship is identified for all the dehydrogenation reactions of C2-C6n-alkanes on Pt substrates, and this linear relationship is independent of the particle size of the Pt substrate and the chain length of alkanes.

Bimetallic Pt-Pd nano-catalyst: size, shape and composition matter

The goal of this work is to understand how the phase diagram (PHAD) of the platinum-palladium (Pt-Pd) alloy changes with size and shape and how it correlates with catalytic properties. By using nano-thermodynamics, the size and shape effects on the PHAD of Pt-Pd nanoparticles were determined theoretically. The PHAD of some nanoparticles (sphere, tetrahedron, octahedron, decahedron, cube, cuboctahedron and rhombic dodecahedron) exhibits a congruent melting point that becomes more and more pronounced when the size decreases. At the right of the congruent melting point i.e. close to the Pt-rich side, the coexistence region exhibits a contraction while an expansion is noticed at larger palladium concentrations. From the Gibbs free energy analysis, the stability of all the considered shapes has been determined versus temperature and composition. Furthermore, the surface segregation was also calculated and it is shown that the surface segregation is reversed at very small sizes. Indeed, below a critical size, Pd does not segregate anymore at the surface like it normally does for larger nanoparticles; but Pt does. The critical size range has been determined for each considered shape; and within this range Pt and Pd co-exist at the surface. Finally, the most catalytically active shapes are predicted to be the tetrahedron and the cube in agreement with the available experimental data and other theoretical results.

Relative Stability of Small Silver, Platinum, and Palladium Doped Gold Cluster Cations

The stability patterns of single silver, platinum, and palladium atom doped gold cluster cations, MAuN−1+ (M = Ag, Pt, Pd; N = 3–6), are investigated by a combination of photofragmentation experiments and density functional theory calculations. The mass spectra of the photofragmented clusters reveal an odd-even pattern in the abundances of AgAuN−1+, with local maxima for clusters containing an even number of valence electrons, similarly to pure AuN+. The odd-even pattern, however, disappears upon Pt and Pd doping. Computed dissociation energies agree well with the experimental findings for the different doped clusters. The effect of Ag, Pt, and Pd doping is discussed on the basis of an analysis of the density of states of the N = 3–5 clusters. Whereas Ag delocalizes its 5s valence electron in all sizes, this process is size-specific for Pt and Pd.

Towards operando computational modeling in heterogeneous catalysis

An increased synergy between experimental and theoretical investigations in heterogeneous catalysis has become apparent during the last decade. Experimental work has extended from ultra-high vacuum and low temperature towards operando conditions. These developments have motivated the computational community to move from standard descriptive computational models, based on inspection of the potential energy surface at 0 K and low reactant concentrations (0 K/UHV model), to more realistic conditions. The transition from 0 K/UHV to operando models has been backed by significant developments in computer hardware and software over the past few decades. New methodological developments, designed to overcome part of the gap between 0 K/UHV and operando conditions, include (i) global optimization techniques, (ii) ab initio constrained thermodynamics, (iii) biased molecular dynamics, (iv) microkinetic models of reaction networks and (v) machine learning approaches. The importance of the transition is highlighted by discussing how the molecular level picture of catalytic sites and the associated reaction mechanisms changes when the chemical environment, pressure and temperature effects are correctly accounted for in molecular simulations. It is the purpose of this review to discuss each method on an equal footing, and to draw connections between methods, particularly where they may be applied in combination.

New molecular insights into the stability of Ni–Pd hollow nanoparticles

Molecular dynamics simulations were used to investigate the structural evolution and thermal behavior of Ni–Pd hollow nanoparticles.

Using structural diversity to tune the catalytic performance of Pt nanoparticle ensembles

While reducing the size, and restricting shape of nanocatalysts can improve performance, monodispersed samples are not necessarily ideal.

Size-dependent subnanometer Pd cluster (Pd4, Pd6, and Pd17) water oxidation electrocatalysis

Water oxidation is a key catalytic step for electrical fuel generation. Recently, significant progress has been made in synthesizing electrocatalytic materials with reduced overpotentials and increased turnover rates, both key parameters enabling commercial use in electrolysis or solar to fuels applications. The complexity of both the catalytic materials and the water oxidation reaction makes understanding the catalytic site critical to improving the process. Here we study water oxidation in alkaline conditions using size-selected clusters of Pd to probe the relationship between cluster size and the water oxidation reaction. We find that Pd4 shows no reaction, while Pd6 and Pd17 deposited clusters are among the most active (in terms of turnover rate per Pd atom) catalysts known. Theoretical calculations suggest that this striking difference may be a demonstration that bridging Pd-Pd sites (which are only present in three-dimensional clusters) are active for the oxygen evolution reaction in Pd6O6. The ability to experimentally synthesize size-specific clusters allows direct comparison to this theory. The support electrode for these investigations is ultrananocrystalline diamond (UNCD). This material is thin enough to be electrically conducting and is chemically/electrochemically very stable. Even under the harsh experimental conditions (basic, high potential) typically employed for water oxidation catalysts, UNCD demonstrates a very wide potential electrochemical working window and shows only minor evidence of reaction. The system (soft-landed Pd4, Pd6, or Pd17 clusters on a UNCD Si-coated electrode) shows stable electrochemical potentials over several cycles, and synchrotron studies of the electrodes show no evidence for evolution or dissolution of either the electrode material or the clusters.

DFT study of the structures and energetics of 98-atom AuPd clusters

The energetics, structures and segregation of 98-atom AuPd nanoclusters are investigated using a genetic algorithm global optimization technique with the Gupta empirical potential (comparing three different potential parameterisations) followed by local minimizations using Density Functional Theory (DFT) calculations. A shell optimization program algorithm is employed in order to study the energetics of the highly symmetric Leary Tetrahedron (LT) structure and optimization of the chemical ordering of a number of structural motifs is carried out using the Basin Hopping Monte Carlo approach. Although one of the empirical potentials is found to favour the LT structure, it is shown that Marks Decahedral and mixed FCC-HCP motifs are lowest in energy at the DFT level.

Magic pairs and structural transitions in binary metallic clusters

Structures and binding energies for bimetallic clusters consisting of a large variety of atomic species are obtained for all atomic sizes N≤40 and all concentrations, using an interatomic potential derived within a quasi-classical description. It is found that increasing the difference between the two types of atoms leads to a gradual disappearance of the well-known homo-atomic geometric magic numbers and the appearance of magic pairs corresponding to the number of atoms of each atomic species in binary nanostructures with higher stability. This change is accompanied by structural transitions and ground-state↔isomer inversions, induced by changes in composition or concentration. There is a clear tendency towards phase separation, the core-shell radial segregation being predominant (energetically favored) in this model.

Tailoring of Pd-Pt bimetallic clusters with high stability for oxygen reduction reaction

The composition-dependent equilibrium structure and thermal stability of Pd-Pt clusters with the size of 55 atoms, and CO, O, OH, and O(2) adsorption on these clusters have been studied using molecular simulation based on the Gupta empirical potential and density functional theory (DFT) calculations. It is found that Pd(43)Pt(12) with a three-shell onionlike structure (TS-cluster) exhibits the highest relative stability in both DFT and Gupta levels and also the highest melting point at the Gupta level among these Pd-Pt clusters. In addition, the Pd(43)Pt(12) TS-cluster possesses the weakest CO, O, OH, and O(2) adsorption strength, compared to the Pt(55), Pd(55), and Pd(13)Pt(42) clusters, indicating good catalytic activities toward the oxygen reduction reaction (ORR) among these Pd-Pt clusters considered. We expect that this kind of DFT-guided strategy by controlling the composition could provide a simple way for possibly searching new electrocatalysts.

Structures and energetics of 98 atom Pd-Pt nanoalloys: potential stability of the Leary tetrahedron for bimetallic nanoparticles

The energetics of 98 atom bimetallic Pd-Pt clusters are studied using a combination of: a genetic algorithm technique (to explore vast areas of the configurational space); a basin-hopping atom-exchange routine (to search for lowest-energy homotops at fixed composition); and a shell optimisation approach (to search for high symmetry isomers). The interatomic interactions between Pd and Pt are modelled by the Gupta many-body empirical potential. For most compositions, the putative global minima are found to have structures based on defective Marks decahedra, but in the composition range from Pd46Pt52 to Pd63Pt35, the Leary tetrahedron (LT)--a structure previously identified for 98 atom Lennard-Jones clusters--is consistently found as the most stable structure. Based on the excess energy stability criterion, Pd56Pt42 represents the most stable cluster across the entire composition range. This structure, a Td-symmetry LT, exhibits multi-layer segregation with an innermost core of Pd atoms, an intermediate layer of Pt atoms and an outermost Pd surface shell (Pd-Pt-Pd). The stability of the Leary tetrahedron is compared against other low-energy competing structural motifs: the Marks decahedron (Dh-M), a "quasi" tetrahedron (a closed-packed structure) and two other closed-packed structures. The stability of LT structures is rationalized in terms of their spherical shape and the large number of nearest neighbours.

In situ investigation of Pt(100--x)Au(x) and Pt(100--y)Sn(y) nanoalloys

Dispersed sols of 1-10 nm sized Pt(100--x)Au(x) and Pt(100--y)Sn(y) nanoalloys have been prepared separately at various x and y above and below the miscibility limit in the bulk metals. Pt(100--x)Au(x) was derived from trisodium citrate reduction of aqueous solutions of H2PtCl6 and HAuCl4. Pt(100--y)-Sn(y) was produced by (i) complexing Sn2+ with glucose at 323 K at pH > 7, (ii) neutralising this with H2PtC16 addition and (iii) reducing the bimetallic precursor with glucose on raising the temperature to 373 K. For Pt(100--x)Au(x (where both metals were zero-valent) as x increased the average size of nanoalloy particles increased. These particles adsorbed onto graphite, where the extent of hydrogen chemisorption at 298 K decreased by 67% at 9 at % Au. Pt/SnO2 nanoparticles (<3 nm in size) were adsorbed onto alumina. The Pt interacted with and catalysed the reduction of SnO2, with some Pt(100--y)Sn(y) nanoalloy formation at about 673 K which even in the bulk occurs over a wider range of compositions than Pt-Au) and enhanced H2 chemisorption at 17-33 at % Sn. Nevertheless some Sn must remain in a positive oxidation state on the alumina surface. The ratio of rates of 2MP/ 3MP formation from MCP and n-hexane may be informative in chemically fingerprinting (and revealing fundamental differences in) these nanoalloy surfaces. The reasons for this are seen in terms of the surface structures on these two types of nanoalloy particles (i.e. the availability of contiguous asymmetric pairs of active surface atoms *, which, as expected, is found to pass through a maximum or decrease beyond specific values of x and y).

Energetic, electronic, and thermal effects on structural properties of Ag-Au nanoalloys

Using a genetic algorithm global optimization approach combined with density functional theory calculations, a search has been made for the lowest energies of (AgAu)(m) nanoalloys with 20-150 atoms (diameters of 1.0-2.0 nm). A total of 31 decahedra, 35 icosahedra, and 2 close-packed motifs are identified in two icosahedral windows and one Marks-decahedral window. These structural motifs have twinned, capped, defective, and distorted atomic packing compared to classical clusters, such as the icosahedron. The magic numbers, atomic ordering, electronic structure, and melting behavior are further studied, and a new poly-nanocrystalline decahedral motif, Ag(44)Au(44), is found to have high structural, electronic, and thermal stability. Our results show that alloying can lead to a remarkable stabilization of local order and provide a comprehensive model for the structures and properties of Ag-Au nanoalloys.

Electronic and Structural Shell Closure in AgCu and AuCu Nanoclusters

The structures of AgCu clusters containing 40 atoms are investigated. The most promising structural families (fcc clusters, capped decahedra, and two types of capped polyicosahedra) are singled out by means of global optimization techniques within an atom-atom potential model. Then, representative clusters of each family are relaxed by means of density-functional methods. It is shown that, for a large majority of compositions, a complex interplay of geometric and electronic shell-closure effects stabilizes a specific polyicosahedral family, whose clusters are much lower in energy and present large HOMO-LUMO gaps. Within this family, geometric and quantum effects concur to favor magic structures associated with core-shell chemical ordering and high symmetry, so that these clusters are very promising from the point of view of their optical properties. Our results also suggest a natural growth pathway of AgCu clusters through high-stability polyicosahedral structures. Results for AuCu clusters of the same size are reported for comparison, showing that the interplay of the different effects is highly material specific.