Global Minimum Pt13M20 (M = Ag, Au, Cu, Pd) Dodecahedral Core

Article pubs.acs.org/JPCA

Global Minimum Pt13M20 (M = Ag, Au, Cu, Pd) Dodecahedral Core− Shell Clusters Dora J. Borbón-González,† Alessandro Fortunelli,*,‡ Giovanni Barcaro,‡ Luca Sementa,‡ Roy L. Johnston,*,§ and Alvaro Posada-Amarillas*,∥ †

Departamento de Matemáticas, Universidad de Sonora, Blvd. Luis Encinas & Rosales, 83000 Hermosillo, Sonora, México CNR, Consiglio Nazionale delle Ricerche, v. G. Moruzzi 1, Pisa 56124, Italy § School of Chemistry, University of Birmingham, Edgbaston, B15 2TT Birmingham, United Kingdom ∥ Departamento de Investigación en Física, Universidad de Sonora, Blvd. Luis Encinas & Rosales, 83000 Hermosillo, Sonora, México J. Phys. Chem. A 2013.117:14261-14266. Downloaded from pubs.acs.org by UNIV OF SUNDERLAND on 10/17/18. For personal use only.

‡

S Supporting Information *

ABSTRACT: In this work, we report finding dodecahedral core−shell structures as the putative global minima of Pt13M20 (M = Ag, Au, Cu, Pd) clusters by using the basin hopping method and the many-body Gupta model potential to model interatomic interactions. These nanoparticles consist of an icosahedral 13-atom platinum core encapsulated by a 20 metal-atom shell exhibiting a dodecahedral geometry (and Ih symmetry). The interaction between the icosahedral platinum core and the dodecahedral shell is analyzed in terms of the increase in volume of the icosahedral core, and the strength and stickiness of M−Pt and M−M interactions. Low-lying metastable isomers are also obtained. Local relaxations at the DFT level are performed to verify the energetic ordering and stability of the structures predicted by the Gupta potential finding that dodecahedral core−shell structures are indeed the putative global minima for Pt13Ag20 and Pt13Pd20, whereas decahedral structures are obtained as the minimum energy configurations for Pt13Au20 and Pt13Cu20 clusters. and multiply twinned faceted nanoparticles.10−12 Thus, the production of nanoparticles exhibiting different shapes has been stimulated, and nowadays there are several different preparation methods, such that multiply faceted polyhedra, resembling Platonic and Archimedean solids in the nanometer size range, have been repeatedly reported (see e.g. ref 13). These high symmetry polyhedra are constructed from basic 3-D building blocks and, according to their symmetry elements, are characterized by their point groups. Nanoalloys are currently the subject of very active research worldwide. They are interesting because of the synergetic effect produced by combining different chemical elements,14 which allows tuning of their physical and chemical properties to develop new tailored nanomaterials. There has been significant effort to comprehend the detailed structures of bimetallic nanoparticles, and experimental and computational studies have contributed equally to the unveiling of unusual, novel morphologies, adopted as a consequence of the synthesis conditions and interatomic interactions.4 Different geometrical shapes have been obtained by varying experimental conditions, and diverse computational models have provided several different morphologies which resemble those obtained by experimentalists. Thus, it is currently accepted that small sized

1. INTRODUCTION Chemical composition, size, and atomic arrangement determine the structures and properties of metal nanoparticles.1 At the nanoscale level, a variety of structural motifs (morphologies) can be found, such as crystalline or amorphous morphologies, or morphologies with symmetries forbidden in a crystal but allowed in isolated systems, such as pentagonal axes.2 Indeed, complexity is the signature of nanoscale materials, and multiple morphologies can be obtained for a given atomic composition and nanoparticle size.3,4 In particular, bimetallic nanoparticles (nanoalloys) exhibit a vast number of structural families which, from an energy landscape perspective, represent local minima on a potential energy hypersurface,5 and the determination of their spatial atomic configuration is of crucial importance for practical applications. Such an increase in structural complexity in nanoalloys is connected with the existence of homotops,6 i.e., isomers with identical overall structures but having different chemical ordering. There is a growing interest in catalytic properties of nanoparticles, especially those composed of transition and noble metal atoms, aiming to provide practical solutions to intricate issues such as pollution control and renewable energy sources.7,8 The potential utility of nanoparticles and the availability of experimental fabrication approaches have placed morphological peculiarities at the forefront of current structural property investigations.9 Controlled synthesis methods have produced fascinating structures including cubes, rods, spirals, © 2013 American Chemical Society

Received: October 10, 2013 Revised: November 22, 2013 Published: December 4, 2013 14261

dx.doi.org/10.1021/jp410079t | J. Phys. Chem. A 2013, 117, 14261−14266

The Journal of Physical Chemistry A

Article

Table 1. Gupta Parameters Used in the Present Work param

Pt−Pt

Ag−Ag

Ag−Pt

Au−Au

Au−Pt

Cu−Cu

Cu−Pt

Pd−Pd

Pd−Pt

A (eV) ξ (eV) p q r0 (Å)

0.2975 2.695 10.612 4.004 2.7747

0.1031 1.1895 10.85 3.18 2.8921

0.175 1.79 10.73 3.59 2.833

0.2061 1.790 10.229 4.036 2.884

0.25 2.20 10.42 4.02 2.83

0.0855 1.224 10.960 2.278 2.556

0.16 1.82 10.786 3.141 2.666

0.1746 1.718 10.867 3.742 2.7485

0.23 2.2 10.74 3.87 2.76

basin hopping algorithm. Quite unexpectedly, we find that in all cases the predicted global minimum (GM) structures have a Pt13 icosahedral core (Ih13) surrounded by an M20 dodecahedral shell. To our knowledge, only once in the previous literature52 has a dodecahedral structure been identified as the global minimum using a combined EP/DFT strategy similar to that followed in the present work. The reasons for the stability of dodecahedral configurations at this size and composition are discussed below. Low-lying metastable isomers are also obtained to provide information on competing structural families at higher energy values and are used as a basis for an exploration of these nanoalloys using a higher-level theoretical approach, i.e., density functional theory (DFT). This combined approach has proven to be very effective in the elucidation of the energetics of different structural families.27 Indeed in the present case we find that only in the Ag−Pt and Pd−Pt case the dodecahedral structure is confirmed to be the putative global minimum of the system.

clusters typically adopt Platonic solid shapes while intermediate and large sized nanoparticles prefer more complex, polyhedral morphologies, such as Archimedean solids, before eventually becoming bulk-like.15 Polyhedral nanoparticles have been experimentally obtained for elemental clusters, including Ag cubes and octahedra,16 Au rhombic dodecahedra,17 and Pt cubic nanoparticles.18 Polyhedral nanoalloys have also been fabricated, such as Pd−Rh concave nanocubes,19 Au−Pd cubes, and other faceted structures.13,15,20 Theoretically, polyhedral nanoparticles have been studied by first-principles methods when the number of atoms is of the order of tens, and by semiempirical potentials (EP) for larger atomic clusters. The optical properties of polyhedral Al and Al−Pb clusters were analyzed by timedependent density functional theory (TD-DFT), showing that these properties can be tuned by controlling the size, shape, and composition, making these clusters excellent candidates for tunable optical applications.21 Montejano-Carrizales et al.15 identified small dodecahedral bimetallic Au−Pd nanoparticles (∼1−2 nm) in high-resolution electron microscopy images. They also reported DFT calculations on dodecahedral monometallic (Au) nanoparticles, generating model structures from geometrical considerations. Marville and Andreoni22 performed semiempirical studies on polyhedral clusters by minimizing the cohesive energy of transition-metal clusters for different geometry structures. Besley et al.23 predicted the stability order of polyhedral Fe clusters of selected geometries for diameters ∼3 nm, using the Murrell-Mottram EP. They predict the existence of stable rhombic dodecahedral (bcc-like) structures for medium-sized iron clusters. Yang et al. performed molecular dynamics simulations to study the thermal stability of icosahedral Pt− Au nanoparticles.24 Monte Carlo simulations have also been reported for Pt-based alloy nanoparticles25 in the size range of 3−4 nm (approximately 1654 atoms). In this work, various truncated-octahedra (TO) configurations were chosen by varying the chemical ordering and using a hybrid optimization scheme on a hypersurface defined by the embedded atom method (EAM) model potential. Core−shell TO structures were predicted for Pt−Ag and Pt−Au clusters, while an onionlike TO structure was found for Pt−Cu clusters and a mixed TO nanoalloy was established as the minimum energy configuration of Pt−Pd clusters. Recently, X-ray experiments have provided evidence for the existence of a Ag32 dodecahedral core consisting of an inner 12-atom hollow icosahedron, encapsulated by a 20-atom dodecahedron, in silver thiolated clusters.26 In contrast to the above-mentioned previous theoretical studies, in the work reported here we have performed an extensive global search of the minimum energy structures of Ptbased clusters from initial randomly generated configurations. Small 33-atom cluster structures Pt13M20, where M is Ag, Au, Cu, and Pd, are analyzed. The potential energy surface (PES) defined by the Gupta EP was explored using the Monte Carlo

2. METHODOLOGY 2.1. Semiempirical Approach. We modeled the interatomic interactions by the Gupta many-body potential,28 which describes the intermetallic bonding of transition and noble metal clusters reasonably well. This is a many-body potential based on Friedel’s tight-binding model, and it consists of a twobody repulsive component and a many-body attractive term. In this interaction model, the configurational energy of a cluster is written as the sum over all the atoms of the attractive and repulsive energy components.29 Their analytic form contains a set of parameters, A, ξ, p, q, and r0, fitted to experimental properties of bulk metals and alloys, such as the cohesive energy, lattice parameters, and independent elastic constants for the reference crystal structure at 0 K; r0 denotes the nearest neighbor distance of the pure bulk elements, often taken as the average of the pure distances but it can also be taken as the experimental nearest-neighbor distance in some specific ordered bulk alloy. Values of the Gupta potential parameters describing Pt−Pt and M−M interactions are taken from the work of Cleri and Rosato.29 The heteronuclear (Pt−M) parameters utilized in this work are the average parameters which have been successfully used previously in theoretical studies of bimetallic clusters.30−34 The potential parameters are reported in Table 1. The basin hopping algorithm (BH) to explore hypersurfaces was systematically investigated by Doye and Wales as a stochastic search strategy, and applied to atomic clusters with simple interaction models.35 Though simple, the BH method has proven to be an effective strategy to study a number of systems described by empirical potentials or force fields.36−38 This method is based on a mathematical transformation of the potential energy function, which is transformed into a collection of interpenetrating staircases.35 For each Pt13M20 cluster our optimizations started from 300 randomly generated configurations. Each search consists of 5000 variable length 14262



Article

Monte Carlo steps, with an initial step size of 1.2 Å and a thermal energy of 0.086 eV (kBT). It is worth noting that different simulation conditions were used, obtaining in all cases identical final results. 2.2. DFT Calculations. DFT calculations were performed using the Quantum ESPRESSO package,39 employing a basis set of plane waves, ultrasoft pseudopotentials,40 and the Perdew−Burke−Ernzerhof (PBE) exchange-correlation (xc-) functional.41 The chosen pseudopotentials include relativistic effects for the Pt atom at the scalar relativistic level. Values of 40 and 240 Ry (1 Ry = 13.606 eV) were used as the energy cutoff for the selection of the plane-wave basis set for the description of the wave function and the electron density, respectively. Eigenvalues and eigenstates of the Kohn−Sham Hamiltonian have been calculated at the Γ point only of a cubic cell of side of approximately 20 Å, applying a Gaussian smearing technique with a broadening of the one-particle levels of 0.03 eV. The DFT local relaxations were performed on the lowestenergy structures obtained from the empirical potential search, and carrying out a DFT reoptimization according to a combined DFT-EP procedure as discussed in ref 42. Metal atom coordinates were fully relaxed until the forces were smaller than 0.05 eV Å−1, and have been checked with relaxations using a threshold of 0.01 eV Å−1 finding total energy values differing by few hundredths of an electronvolt.

3. RESULTS AND DISCUSSION Figure 1 shows the optimized Pt13M20 (M = Ag, Au, Cu, Pd) cluster structures predicted as putative global minima together with the lowest-energy isomers and the corresponding relative energy differences. The core−shell morphology is evident in these nanoalloys, as well as the Ih13 platinum core. Energetic factors favor surface segregation of the M atoms. The higher surface energy and higher cohesive energy of Pt compared to M (Ag, Au, Cu, and Pd) leads to segregation of the majority element to the outer dodecahedral shell. The superior stability of dodecahedra with respect to other possible isomers is due to the high symmetry exhibited by these structures and their chemical ordering, connected with the fact that 13−20 is a magic cluster and a magic composition for dodecahedra, achieving geometric shell closure for both elemental components. Bond counting is suggestive: one finds 132 nearest-neighbor bonds in Pt13M20 as compared to only 128 in the (larger) 34-atom truncated octahedron (TO) obtained from the perfect TO38 by removing one of the (100) facets. Even though some of the surface bonds are strained (see below), the perfect core−shell arrangement, combined with a high number of metal−metal contacts (especially the stronger Pt−Pt and M−Pt interactions), make dodecahedra prevail over other non-closed-shell competitors, such as fcc-like structures, icosahedra or decahedra. It should be noted that the 20 surface atoms (the dodecahedral shell) adopt anti-Mackay coverage (i.e., sitting over triangular faces) with respect to the icosahedral core.43 Finding a dodecahedral configuration as the GM is particularly surprising in the Pt13Cu20 case. In fact, whereas in the case of Ag−Pt44 and Au−Pt45 all energetic factors agree in favoring the Ptcore−Mshell arrangement, and this tendency is (albeit more weakly) also present in the Pd−Pt case,34 the smaller lattice parameter of Cu should favor an inverse Cucore− Ptshell segregation. Indeed, preferential occupation of core sites by Cu is what is usually found in both theory and experiment

Figure 1. Gupta potential global minima structures and low-lying isomers of Pt13M20 (M = Cu, Ag, Au, Pd) clusters. The binding energy difference (in eV) with respect to the most stable structure is given below each structural motif. Dark (blue) spheres represent Pt atoms, while gray (orange, light blue, yellow, and cherry) spheres represent Cu, Ag, Au, and Pd atoms.

for larger clusters and extended surfaces; however, the presence of ligands such as CO can drive Cu to the surface.46,47 The predicted Cu surface segregation is thus a peculiarity of the nanoscale character of these Pt13Cu20 clusters. These remarkable dodecahedral structures, previously reported on the basis of experiment, but so far never found as global minima in the theoretical literature, exhibit some peculiarities according to the shell atom types; i.e., the chemical nature and the symmetry of the dodecahedral cage modify some of the structural characteristics of the icosahedral Pt core, through a core−shell interaction that symmetrically changes the core’s volume. Consider, for example, the central Pt atom distance to the neighboring Pt atoms (r0): Table 2 lists the nearest-neighbor distances for the clusters considered in this study (the bare Pt13 cluster distance values are presented here as a reference).48 As a result of the core−shell interaction, the Pt13 cluster in M20 slightly expands, depending on the M−Pt bond strength, which is larger for Pd−Pt and Au−Pt than for Ag−Pt and Cu−Pt. The other two geometrical parameters reported in Table 2 are the M−Pt and M−M distances which are clearly in competition: for a given radial distance of the M atoms, the shorter D(M−Pt) the larger D(M−M). This competition is realized in different ways in the systems here considered: in the Cu−Pt case, e.g., one finds very short 14263



Article

Table 2. Selected Nearest-Neighbor Distances (D) for the Global Minima Dodecahedral Structures of Pt13M20 (M = Cu, Ag, Au, Pd=) Clusters Considered in This Worka D

Pt13

Pt−Ag

Pt−Cu

Pt−Pd

Pt−Au

m

r0 M−Pt M−M

2.566

2.626 2.711 3.054

2.627 2.588 2.944

2.668 2.599 2.964

2.670 2.669 3.028

12 60 30

a

r0 is the radial distance of the 12 Ih Pt atoms, M−Pt and M−M are the closest M−Pt and M−M distances, respectively. Multiplicity (m) indicates the distance’s frequency of occurrence. All distances are in Å.48.

D(Cu−Pt), but quite elongated D(Cu−Cu). The basic quantities which determine this behavior are the strength and the ‘stickiness’ of the various metal−metal interactions, i.e., the relative easiness with which a metal−metal bond is strained.49 The greater the strength and stickiness of M−Pt with respect to M−M, the shorter D(M−Pt) and the longer D(M−M). Stickiness in general increases going down a column of the periodic table, and it is large for third-row transition and noble metals and in particular for Au and Pt, but it is small for Cu. Ag stickiness is intermediate between Cu and Pt whereas Pd stickiness is more or less similar to that of Pt depending on the empirical potential parametrization. Indeed we find strongly elongated Cu−Cu distances: D(Cu−Cu) = 2.944 Å compared to 2.544 Å in the bulk. In contrast, D(Ag−Ag) = 3.054 Å is much less elongated with respect to the bulk value of 2.889 Å. Similar behavior is observed for distances D(Pd−Pd) and D(Au−Au), whose relationship to bulk values are 2.964 Å with respect to 2.75 Å for D(Pd−Pd), and 3.028 Å with respect to 2.88 Å in the case of D(Au−Au). In addition to global minima, our PES exploration also provides isomeric structures. This is relevant because two types of effects are found when the level of theory includes electronic interactions: energetic reordering or structural distortions. The lower symmetry of the chemical ordering of the isomeric structures obtained in this work suggests structurally dissimilar configurations as possible competing motifs to be considered in, for example, DFT calculations. Indeed, the binding energy/ atom (−Eclust/N) difference between the lowest energy and first isomer structures of, e.g., the Pt13Cu20 cluster is rather small at the empirical potential theory level, which emphasizes the importance of carrying out first principles calculations. A detailed study of Pt13M20 (M = Cu, Ag, Au, Pd) clusters by the DFT method was then conducted with the aim of clarifying the energetics of these nanoparticles at a higher theoretical level. As mentioned in section 2.2, several structural motifs were singled out in the empirical potential global optimization search and were reoptimized at the DFT level. Results are shown in Figure 2, where we present only the lowest energy and first isomeric structures obtained after DFT optimization. Coordinates and the corresponding geometric characterization for these and other highest energy isomers are available in Supporting Information. By inspecting this figure, we observe that the dodecahedral structure is confirmed as the lowestenergy structure at the DFT level for Pt13Ag20 cluster, exhibiting a large energy difference of ∼1 eV between this structure and the first higher-energy structure (first isomer), which makes the dodecahedral structure a very stable one. Dodecahedral geometry was also obtained as the putative global minimum for the Pt13Pd20 cluster, for which the energy difference between the dodecahedral structure and the first

Figure 2. DFT relaxed structures and low-lying isomers of Pt13M20 (M = Cu, Ag, Au, Pd) clusters. The binding energy difference (in eV) with respect to the most stable structure is given below each structural motif. Color coding is similar to that in Figure 1.

isomer is nevertheless smaller: 0.24 eV. In contrast, in the case of the Pt13Cu20 and Pt13Au20 systems, the putative global minimum structures are Marks decahedra (Dh-Mk) with an energy difference of 0.18 eV (Cu−Pt) and 0.10 eV (Au−Pt) with respect to the first isomer. In both cases, the dodecahedral structure is higher in energy. With reference to the ground state structure of Pt13Cu20 and Pt13Au20 clusters, the corresponding dodecahedral structures have an energy difference of 0.30 and 1.57 eV, respectively. The reasons for the disagreement between empirical potential and DFT predictions are different in the Au−Pt and Cu−Pt cases. In the former case, the strong directionality effects which play a decisive role in Au−Pt interactions are not well accounted for by the isotropic Gupta potential.28 In the latter case, instead, in addition to slightly stronger directionality effects DFT predicts a stickier Cu−Cu bonding which does not survive the large elongation imposed by the anti-Mackay arrangement of the surface shell (see Table 2), explaining the energy destabilization of 0.30 eV thus found. In any case, the present results suggest that in appropriate conditions of size and composition dodecahedral arrangements may be the thermodynamically favored global minima of the system (such as in the present Ag−Pt and Pd−Pt examples) or low-energy isomers (such as in the present Cu−Pt example), thus providing a possible rationalization of the direct identification of the bimetallic dodecahedral nanoparticles in, for example, HAADF-STEM experiments.15,50 An additional 14264



Article

(3) Borbón-González, D. J.; R. Pacheco-Contreras, R.; PosadaAmarillas, A.; Schön, J. C.; Johnston, R. L.; Montejano-Carrizales, J. M. Structural Insights into 19-Atom Pd/Pt Nanoparticles: A Computational Perspective. J. Phys. Chem. C 2009, 113, 15904−15908. (4) Pacheco-Contreras, R.; Borbón-González, D. J.; Dessens-Félix, M.; Paz-Borbón, L. O.; Johnston, R. L.; Schön, J. C.; Jansen, M.; Posada-Amarillas, A. Determination of the Energy Landscape of Pd12Pt1 Using a Combined Genetic Algorithm and Threshold Method. RSC Adv. 2013, 3, 11571−11579. (5) Wales, D. J. Energy Landscapes; Cambridge University Press: Cambridge, U.K., 2003. (6) Jellinek, J.; Krissinel, E. B. NinAlm Alloy Clusters: Analysis of Structural Forms and Their Energy Ordering. Chem. Phys. Lett. 1996, 258, 283−292. (7) Haruta, M.; Yamada, N.; Kobayashi, T.; Iijima, S. Gold Catalysts Prepared by Coprecipitation for Low-Temperature Oxidation of Hydrogen and of Carbon Monoxide. J. Catal. 1989, 115, 301−309. (8) Nørskov, J. K.; Bligaard, T.; Rossmaisl, J.; Christensen, C. H. Towards the Computational Design of Solid Catalysts. Nat. Chem. 2009, 1, 37−46. (9) Wang, Y.; Xie, S.; Liu, J.; Park, J.; Huang, C. Z.; Xia, Y. ShapeControlled Synthesis of Palladium Nanocrystals: A Mechanistic Understanding of the Evolution from Octahedrons to Tetrahedrons. Nano Lett. 2013, 13, 2276−2281. (10) Huang, X.; Li, Y.; Li, Y.; Zhou, H.; Duan, X.; Huang, Y. Synthesis of PtPd Bimetal Nanocrystals with Controllable Shape, Composition, and Their Tunable Catalytic Properties. Nano Lett. 2012, 12, 4265−4270. (11) Garcia-Gutierrez, D. I.; Gutierrez-Wing, C. E.; Giovanetti, L.; Ramallo-López, J. M.; Requejo, F. G.; Jose-Yacaman, M. Temperature Effect on the Synthesis of Au−Pt Bimetallic Nanoparticles. J. Phys. Chem. B 2005, 109, 3813−3821. (12) Velázquez-Salazar, J. J.; Esparza, R.; Mejía-Rosales, S. J.; EstradaSalas, R.; Ponce, A.; Deepak, F. L.; Castro-Guerrero, C.; José-Yacamán, M. Experimental Evidence of Icosahedral and Decahedral Packing in One-Dimensional Nanostructures. ACS Nano 2011, 5, 6272−6278. (13) Quan, Z. W.; Wang, Y. X.; Fang, J. Y. High-Index Faceted Noble Metal Nanocrystals. Acc. Chem. Res. 2013, 46, 191−202. (14) Bazin, D.; Guillaume, D.; Pichon, Ch.; Uzio, D.; Lopez, S. Structure and Size of Bimetallic Palladium-Platinum Clusters in an Hydrotreatment. Oil Gas Sci. Technol. 2005, 60, 801−813. (15) Montejano-Carrizales, J. M.; Rodríguez-López, J. L.; Pal, U.; Miki-Yoshida, M.; José-Yacamán, M. The Completion of the Platonic Atomic Polyhedra: The Dodecahedron. Small 2006, 2, 351−351. (16) Zeng, J.; Zheng, Y.; Rycenga, M.; Tao, J.; Li, Z. Y.; Zhang, Q.; Zhu, Y.; Xia, Y. Controlling the Shapes of Silver Nanocrystals with Different Capping Agents. J. Am. Chem. Soc. 2010, 132, 8552−8553. (17) Jeong, G. H.; Kim, M.; Lee, Y. W.; Choi, W.; Oh, W. T.; Park, Q.-H.; Han, S. W. Polyhedral Au Nanocrystals Exclusively Bound by {110} Facets: The Rhombic Dodecahedron. J. Am. Chem. Soc. 2009, 131, 1672−1673. (18) Shen, Z.; Matsuki, Y.; Shimoda, T. Preparation of Large Thermally Stable Platinum Nanocubes by Using Solvent-Thermal Reaction. Chem. Commun. 2010, 46, 8606−8608. (19) Xie, S.; Lu, N.; Xie, Z. X.; Wang, J.; Kim, M. J.; Xia, Y. Synthesis of Pd-Rh Core−Frame Concave Nanocubes and Their Conversion to Rh Cubic Nanoframes by Selective Etching of the Pd Cores. Angew. Chem., Int. Ed. 2012, 51, 10266−10270. (20) Kaiser, J.; et al. Catalytic Activity of Nanoalloys from Gold and Palladium. Phys. Chem. Chem. Phys. 2012, 14, 6487−6495. (21) Xie, R.-H.; Bryant, G. W.; Zhao, J.; Kar, T.; Smith, V. H., Jr. Tunable Optical Properties of Icosahedral, Dodecahedral, and Tetrahedral Clusters. Phys. Rev. B 2005, 71, 125422−1−125422−5. (22) Marvillet, L.; Andreoni, W. Size Dependence of the Structural Properties of Transition-Metal Aggregates from an Empirical Interatomic Potential Scheme. J. Phys. Chem. 1987, 91, 2645−2649. (23) Besley, N. A.; Johnston, R. L.; Stace, A. J.; Uppenbrink, J. Theoretical Study of the Structures and Stabilities of Iron Clusters. THEOCHEM 1995, 341, 75−90.

analysis of the Pt13M20 clusters’ electronic properties revealed no electronic shell closure for any of the systems studied here.

4. CONCLUSIONS We obtain the lowest-energy structures and low-lying isomers of free Pt13M20 clusters (M = Ag, Au, Cu, Pd) through a thorough global optimization approach using the basin hopping method. We find that highly symmetrical structures are energetically favored as putative global minima at the empirical potential level for this cluster size, consisting of a M20 dodecahedral cage encapsulating an Ih13 platinum core. Because of the symmetry, the core−shell interaction is isotropic and causes an expansion of the Ih13 Pt core depending on the nature of the monolayer shell, in the following order: Pt−Ag ≈ Pt−Cu < Pt−Pd ≈ Pt−Au. The radial position of the surface shell is determined by the competition between M−Pt and M−M bonds. The magic character of these dodecahedral motifs in terms of both structural and compositional shell closure, with the more cohesive element (Pt) segregating into an Ih core, explains the superior stability of dodecahedra with respect to other competing motifs. To the best of our knowledge, this is the first time that a dodecahedral structure has been predicted as the global minimum by a systematic PES exploration rather than enforced by a biased model, thus making a first step in rationalizing the common occurrence of this and other type of structures in experiments, such as the complete 45-atom antiMackay icosahedron,51 which results from capping each of the pentagonal faces of the 33-atom dodecahedral structure. DFT calculations confirm the energetic preference for dodecahedral structures in bimetallic Pt13Ag20 and Pt13Pd20 clusters, and their character of low-energy isomers for Pt13Cu20, thus providing a possible rationalization of the experimental observation of nanoparticles with dodecahedral symmetry.15 No electronic shell closure was found for any of the Pt-based clusters studied here.

■

ASSOCIATED CONTENT

S Supporting Information *

Coordinates of the relaxed structures of the ground state configurations and higher energy isomers for the Pt-based clusters analyzed in this article. This material is available free of charge via the Internet at http://pubs.acs.org.

■

AUTHOR INFORMATION

Corresponding Authors

*E-mail: [email protected]. *E-mail: [email protected]. *E-mail: [email protected]. Notes

The authors declare no competing financial interest.

■

ACKNOWLEDGMENTS A.P.-A. acknowledges CONACYT-México for financial support through Project 180424. A.F. acknowledges financial support by the ERC-AG SEPON project.

■

REFERENCES

(1) Johnston, R. L. Atomic and Molecular Clusters; Taylor and Francis, London, 2002. (2) Ferrando, R.; Johnston, R. L.; Jellinek, J. Nanoalloys: From Theory to Applications of Alloy Clusters and Nanoparticles. Chem. Rev. 2008, 108, 845−910. 14265



Article

(48) Note: Assuming (0, 0, zPt) as the Cartesian coordinates of one of the 12 Ih Pt atoms and (0, yM, zM) as those of one of the 20 dodecahedral M atoms, and defining r0, D(M−Pt), and D(M−M) as the basic nearest-neighbor distances reported in Table 2, one has the following relationships: zPt = r0, D(M−M) = 2yM sin(36°), zM = zPt + (D(M−Pt)2 − y2M)1/2, implying 2 sin(36°)D(M−Pt) ≥ D(M−M). (49) Baletto, F.; Ferrando, R.; Fortunelli, A.; Montalenti, F.; Mottet, C. Crossover Among Structural Motifs in Transition and Noble-Metal Clusters. J. Chem. Phys. 2002, 116, 3856−3863. (50) Wang, Z. W.; Palmer, R. E. Direct Atomic Imaging and Dynamical Fluctuations of the Tetrahedral Au20 Cluster. Nanoscale 2012, 4, 4947−4949. (51) Rossi, G.; Rapallo, A.; Mottet, C.; Fortunelli, A.; Baletto, F.; Ferrando, R. Magic Polyicosahedral Core-Shell Clusters. Phys. Rev. Lett. 2004, 93, 105503−1−105503−4. (52) Damianos, K.; Solokha, P.; Ferrando, R. Core-Shell and Matryoshka Structures in MgNi Nanoalloys: A Computational Study. RSC Adv. 2013, 3, 9419−9430.

(24) Yang, Z.; Yang, X.; Xu, Z. Molecular Dynamics Simulation of the Melting Behavior of Pt−Au Nanoparticles with Core−Shell Structure. J. Phys. Chem. C 2008, 112, 4937−4947. (25) Yun, K.; Cho, Y.-H.; Cha, P.-R.; Lee, J.; Nam, H.-S.; Oh, J. S.; Choi, J.-H.; Lee, S.-C. Monte Carlo Simulations of the Structure of PtBased Bimetallic Nanoparticles. Acta Mater. 2012, 60, 4908−4916. (26) Desireddy, A.; et al. Ultrastable Silver Nanoparticles. Nature 2013, 501, 399−402. (27) Dessens-Félix, M.; Pacheco-Contreras, R.; Barcaro, G.; Sementa, L.; Fortunelli, A.; Posada-Amarillas, A. Structural Motifs of Bimetallic Pt101‑xAux Nanoclusters. J. Phys. Chem. C 2013, 117, 20976−20974. (28) Gupta, R. P. Lattice Relaxation at a Metal Surface. Phys. Rev. B 1981, 23, 6265−6270. (29) Cleri, F.; Rosato, V. Tight-Binding Potentials for Transition Metals and Alloys. Phys. Rev. B 1993, 48, 22−33. (30) Guerrero-Jordan, J. Computational Modelling of Nanoalloys. M. Sc. Thesis, Universidad de Sonora, 2012. (31) Paz-Borbón, L. O.; Johnston, R. L.; Barcaro, G.; Fortunelli, A. Structural Motifs, Mixing, and Segregation Effects in 38-Atom Binary Clusters. J. Chem. Phys. 2008, 128, 134517−1−134517−12. (32) Logsdail, A.; Paz-Borbón, L. O.; Johnston, R. L. Structures and Stabilities of Platinum-Gold Nanoclusters. J. Comput. Theor. Nanosci. 2009, 6, 857−866. (33) Massen, C.; Mortimer-Jones, T. V.; Johnston, R. L. Geometries and Segregation Properties of Platinum−Palladium Nanoalloy Clusters. J. Chem. Soc., Dalton Trans. 2002, 4375−4388. (34) Barcaro, G.; Fortunelli, A.; Polak, M.; Rubinovich, L. Patchy Multishell Segregation in Pd−Pt Alloy Nanoparticles. Nano Lett. 2011, 11, 1766−1769. (35) Wales, D. J.; Doye, J. P. K. Global Optimization by BasinHopping and the Lowest Energy Structures of Lennard-Jones Clusters Containing up to 110 Atoms. J. Phys. Chem. A 1997, 101, 5111−5116. (36) Wales, D. J.; Head-Gordon, T. Evolution of the Potential Energy Landscape with Static Pulling Force for Two Model Proteins. J. Phys. Chem. B 2012, 116, 8394−8411. (37) Molayem, M.; Grigoryan, V. G.; Springborg, M. Theoretical Determination of the Most Stable Structures of NimAgnBimetallic Nanoalloys. J. Phys. Chem. C 2011, 115, 7179−7192. (38) Ferrando, R.; Rossi, G.; Nita, F.; Barcaro, G.; Fortunelli, A. Interface-Stabilized Phases of Metal-on-Oxide Nanodots. ACS Nano 2008, 2, 1849−1856. (39) Giannozzi, P.; et al. QUANTUM ESPRESSO: A Modular and Open-Source Software Project for Quantum Simulations of Materials. J. Phys.: Condens. Matter. 2009, 21, 395502-1−395502-19. (40) Vanderbilt, D. Soft Self-Consistent Pseudopotentials in a Generalized Eigenvalue Formalismo. Phys. Rev. B 1990, 41, 7892− 7895. (41) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865−3868. (42) Ferrando, R.; Fortunelli, A.; Johnston, R. L. Searching for the Optimum Structures of Alloy Nanoclusters. Phys. Chem. Chem. Phys. 2008, 10, 640−649. (43) Baletto, F.; Ferrando, R. Structural Properties of Nanoclusters: Energetic, Thermodynamic, and Kinetic Effects. Rev. Mod. Phys. 2005, 77, 371−423. (44) Negreiros, F. R.; Taherkhani, F.; Parsafar, G.; Caro, A.; Fortunelli, A. Kinetics of Chemical Ordering in a Ag-Pt Nanoalloy Particle Via First-Principles Simulations. J. Chem. Phys. 2012, 137, 194302-1−194302-7. (45) Wanjala, B. N.; Luo, J.; Fang, B.; Mott, D.; Zhong, C.-J. GoldPlatinum Nanoparticles: Alloying and Phase Segregation. J. Mater. Chem. 2011, 21, 4012−4020. (46) West, P. S.; Johnston, R. L.; Barcaro, G.; Fortunelli, A. The Effect of CO and H Chemisorption on the Chemical Ordering of Bimetallic Clusters. J. Phys. Chem. C 2010, 114, 19678−19686. (47) Andersson, K. J.; Calle-Vallejo, F.; Rossmeisl, J.; Chorkendorff, I. Adsorption-Driven Surface Segregation of the Less Reactive Alloy Component. J. Am. Chem. Soc. 2009, 131, 2404−2407. 14266


Global Minimum Pt13M20 (M = Ag, Au, Cu, Pd) Dodecahedral Core

Recommend Documents