Determination of the energy landscape of Pd12Pt1 using a combined genetic algorithm and threshold energy method3

In this work we present a thorough exploration of the potential energy surface (PES) of Pd–Pt bimetallic nanoparticles at the specific composition Pd12Pt1, using the combination of a genetic algorithm and the threshold method for global optimization and exploration of the barrier structure, employing the semiempirical Gupta many-body potential for modeling the interatomic interactions. The structural and energetic analysis of Pd12Pt1 nanoparticles, including binding energies (Eb), symmetries and commonneighbor analysis (CNA) allowed us to identify a large set of representative structures of local minima, with an icosahedral motif found to be the putative global minimum for Pd12Pt1. A detailed study of the icosahedral motif was carried out by an exhaustive exploration of low energy isomers, in order to understand qualitatively structural interconversion. 2-D tree (disconnectivity) graphs are plotted to map the structures of minima on the PES of Pd12Pt1. DFT calculations were performed on representative structures to establish the energetic hierarchy and structural stability.


Introduction
3][4][5][6][7] At the nanoscale, metallic nanoalloys have a high surface/volume ratio implying a large number of surface atoms, a vital property which makes them ideal as catalysts, as many reactions take place at nanoparticle surface sites.6][7][8][9][10] Compared to monometallic clusters, bimetallic nanoalloys such as Pd-Pt show a wider range of structural possibilities.It is known that the preparation conditions of the bimetallic nanoparticles can directly influence their size, shape and segregation properties. 5,11Using a combination of experimental techniques (transmission electron microscopy (TEM), x-ray absorption spectroscopy (XAS) and low energy ion scattering (LEIS)), Renouprez, Rousset et al. demonstrated that their Pd-Pt nanoparticles (with diameters of 1-5 nm) were truncated octahedra with a Pt-rich core surrounded by a Pd-rich shell. 8,9azin et al. have studied supported Pd-Pt nanoalloys deposited on c-alumina, using a combination of XAS, TEM and volumetric H 2 O 2 titration.They showed that their Pd-Pt particles, with sizes ,1 nm, had ''cherry-like'' (spherical) structures with Pt atoms occupying core sites and with Pd atoms decorating the surface. 10Habas et al. have also reported the characterization of Pt core Pd shell nanoparticles, obtained by epitaxial seeded growth. 12On the other hand, D'Souza and Sampath have reported the preparation of Pd-Pt nanoparticles by the sol-gel process.These bimetallic nanoparticles (with sizes 1-3 nm) were characterized by X-ray photoelectron spectroscopy (XPS) and CO adsorption measurements, the results of which indicated that the particles have a Pd core and a Pt shell. 13Sanchez et al. have shown, using high angle annular dark field scanning transmission electron microscopy (HAADF-STEM), that Pt core Pd shell , Pd core Pt shell or alloyed Pd-Pt configurations can be generated by sequential reduction or coreduction of Pd and Pt salts. 14rom a theoretical point of view, a challenging problem in cluster science is to determine the most stable isomer of a particular cluster for a given size (number of atoms) and composition as this involves optimizing both the geometry and chemical ordering (degree of segregation or mixing).6][17] For the case of bimetallic (binary) nanoparticles, the problem increases dramatically, as for each geometrical isomer there is a combinatorial dependence of the number of homotops (inequivalent permutational isomers) on the composition. 18Global optimization techniques, [19][20][21][22][23] such as genetic algorithms, [24][25][26][27] basin hopping Monte Carlo [28][29][30] and the threshold method 31,32 have been developed to overcome this problem.This has allowed the exploration of large regions of the PES of monometallic and bimetallic nanoparticles (both in the gas-phase or supported on an oxide substrate) and the successful identification of their corresponding GM structures within a vast number of geometrical isomers and homotops.
Since the pioneering work of Goldstein, 33 it has been clear that many of the intrinsic physical and chemical properties of complex systems (e.g.4][35] An accurate PES description is thus needed in order to locate the putative GM configuration as well as the pathways connecting them via the corresponding transition states (TS) to other local-minimum configurations.
In this work we present a thorough exploration of the PES for the bimetallic cluster Pd 12 Pt 1 , using a combination of a genetic algorithm and the threshold method for global optimization and exploration of the barrier structure, employing the semi-empirical Gupta many-body potential for modeling the interatomic interactions.Structural and energetic analysis of Pd 12 Pt 1 nanoparticles include the calculation of binding energies (E b ), common-neighbor analysis (CNA) indices and their normalized relative abundance (RA) 36 in order to understand plausible structural interconversion.It is worth mentioning that the CNA technique gives a threedimensional representation of the atomic configuration around a selected bonded pair of atoms (root pair) which facilitate the structural classification in terms of a set of three numbers (indices) that provide the order type and corresponding abundance. 37Tree graphs 32,[38][39][40][41] are plotted to map both the representative structures of the minima on the PES of Pd 12 Pt 1 and a variety of archetypical multiply twinned (MT) structures found in our searches.Prominent structural isomers were found in our searches, ranging for example, from helical (right and left-handed) configurations and slightly distorted biplanar (BP) structures to icosahedral motifs (Ih).The minimum energy configuration in the case of the BP structure has the Pt atom located at the center of the hexagonal plane, while in the helical structures the minimum energy configuration has the Pt atom in the second position of the 5-atom strand.Disconnectivity graphs show the homotops energetical ordering of selected MT structures, rationalizing the different tetrahedral growth types in terms of some basic information provided by the energy landscape analysis.Energy landscape analysis including the determination of the barrier structure is a novel approach for bimetallic clusters, and has been utilized to generate a representation of a partial PES region, 61 thus yielding a deeper understanding of the entire potential energy surface.
Finally, the relative stability of putative minima is tested against more accurate methods, namely density functional theory (DFT). 42,43Full relaxations are carried out on selected configurations and the results are discussed in terms of the energetic ordering and structural properties of these clusters, including the novel helical and biplanar Pd 12 Pt 1 nanoparticles.

The Gupta potential
In this study, the intermetallic interaction is modeled by the Gupta many-body potential, 44 which is a phenomenological potential energy function that has been shown to describe transition and noble metal cluster structures reasonably accurately. 45The Gupta potential is an example of a manybody potential which is based on the second-moment approximation of the tight binding density of states.The potential consists of a two-body repulsive component and a many-body attractive component.The homonuclear parameters are fitted to the cohesive energy, lattice parameter and independent elastic constants for the bulk metals in the corresponding crystal structure at 0 K.The heteronuclear (Pd-Pt) parameters used here are the average parameters introduced by Massen et al. 46 which have been used previously. 26At the Gupta potential level of theory, the binding energy (E b ) per atom of a cluster in our work is calculated as the total energy of the cluster (E Total ), divided by the total number of atoms (N).Thus, negative values imply improved structural stability.

Parameters of the ab initio DFT calculations
8][49] A gradient corrected exchangecorrelation functional is employed as proposed by Perdew, Burke and Ernzerhof (PBE). 50The one-electron Kohn-Sham orbitals are expanded in a localized numerical basis set.A double numerical basis set is used together with polarization functions (dnp).A real space cutoff of 4.9 Å is used for the basis functions.A pseudo potential 51 which also accounts for the relativistic effects is employed for Pd and Pt to describe the interaction between the valence electrons and the nuclei together with the inner shell electrons.The states treated in the valence shell are 4s 2 4p 6 4d 10 and 5s 2 5p 6 5d 9 6s 1 for Pd and Pt, respectively.The Kohn-Sham equations are solved self-consistently using an integration technique based on weighted overlapping spheres centered at each atom.The direct Coulomb potential is obtained by projection of the charge density onto angular dependent weighting functions also centered at each atom.The Poisson equation is, thereafter, solved by one-dimensional integration.3][54][55] Structures are considered to be relaxed when changes in the energy, the largest element of the gradient and the atomic displacements are smaller than 10 25 Ha (i.e.2.72 6 10 24 eV), 0.05 eV Å 21 and 0.005 Å, respectively.All calculations are performed spin unrestricted.The relative stability of the different Pd 12 Pt 1 cluster geometries is analyzed at the DFT level by calculating their corresponding average binding energy per atom: where N denotes the total number of atoms in the bimetallic cluster.

The genetic algorithm
Genetic algorithms (GA) are global search optimization techniques, based on the concepts of natural selection, which have been widely used to obtain minimum energy structures of atomic and molecular clusters. 56Initially, a collection of individuals (clusters) is generated at random and then energyminimized to produce an initial population, where each member of the population corresponds to a local minimum on the PES.Secondly, a fitness function is used by the GA in order to determine which individuals will survive from one generation to the next.The fitness value depends on the cluster's total energy such that low-energy clusters will have high fitness and vice versa.Subsequently, individuals are selected to take part in crossover (''mating'') based on their fitness.A cluster is accepted for mating if its fitness value is greater than a randomly generated number between 0 and 1.In this way, parent structures are combined in order to generate offspring.Crossover moves are performed using the "cut-and-splice" method introduced by Deaven and Ho, 57 which involves randomly rotating two parent clusters, cutting both horizontally in one or two positions and then combining complementary fragments. 58As in nature, mutation is introduced in order to maintain population diversity.In our optimizations, each cluster has the same probability of being mutated. 58The specific details of the GA we have adopted 56 and the descriptions of previous applications can be found elsewhere. 26,46,593 The threshold energy method 2.3.1 Description of method.In order to survey the PES of Pd 12 Pt 1 , we have used the threshold method.31,32 This is a Monte Carlo based approach which makes use of random walks below a sequence of prescribed energy lids combined with stochastic quenches.This method requires a low-energy configuration usually obtained from other optimization techniques (in this case the GA) initial input.Starting from the GM or a local minimum, random walks are carried out through the configurational space, under the constraint that a given energy E lid threshold is not exceeded all moves that fulfill this criterion are accepted.After a certain number of steps the system is relaxed into one neighboring minimum.One usually performs several stochastic quenches from the same stopping point in order to determine whether the walker is inside a basin or in a transition region.By carrying out a number of Monte Carlo runs at different E lid values, energy barriers separating the putative GM from other (higher-energy) minimum structures can be estimated.Such a set of threshold runs is repeated starting from many local minima on the landscape.In this way, the iterative procedure allows us to carefully explore the PES landscape.
A set of simulation parameters selected to facilitate an extensive exploration of the PES determine the success of finding minimum energy structures.Their values depend on the system's size, the number of chemical elements involved, the PES function complexity, etc.In this study, the parameters needed in the G42{ code were obtained as a result of running several tests before finding the simulation conditions which allowed us to have the most accurate view of the PES belonging to the Pd 12 Pt 1 nanoalloy.
2.3.2Application to the specific system under investigation.Particular details of the PES exploration are as follows: a ladder of 30 energy lids (E lid ) was constructed in order to carry out each of the runs; the energy gap between successive lids was 0.1 eV/atom.The first lid was placed at 23.5 eV/atom, an energy value slightly above that corresponding to the icosahedral structure, whereas the final lid was placed very close to the cluster's dissociation energy, at 20.6 eV/atom.Within this energy window, the system was allowed to evolve in order to access several local minima.For each run below every energy lid, a different seed was used, releasing a set of random walkers for each seed.The random walks in each run accept movements limited by the value of the energy lid, a procedure which will eventually give access to other local minima by crossing energy barriers.Three different seeds were used for each energy lid, 10 random walkers for each seed giving each walker 2.5 6 10 5 Monte Carlo steps and, at the stopping point, 20 stochastic quenches of 2.5 6 10 5 steps each, to find out if the random walk had led to another local minimum.Local relaxations were carried out on each new structure in order to reach the bottom of the local basin.The L-BFGS-B algorithm of Zhu et al. 60 was used at the end of each quench to efficiently reach the bottom of the basin.
The threshold runs provided 18 000 configuration files which were analyzed taking into account structural and energetic criteria (CNA indices and total energy).We found only 5016 different configurations (minima) which permit us, at least in principle, to gain insight on the Pd 12 Pt 1 cluster PES, i.e., on the physical and chemical properties of this nanoalloy.It was necessary to analyze the morphology of each of the different configurations to elucidate the existence of structural or compositional variations.This procedure yielded 79 structures and made it possible to build the disconnectivity graph for the energy region where these structures were found computationally.For this, 25 structures out of the 79 were { G42 is a modular code consisting of a large package of routines for the exploration of energy landscapes.chosen as representatives of the cluster's energy landscape, based on geometrical and compositional considerations.
After identifying the different structures with the lowest local energy, we constructed their disconnectivity graphs starting in the lowest energy minimum selected for each isomer and employing a set of 10 energy thresholds with a spacing between energy lids of 0.01 eV/atom.Reaching the top energy lid allows us to determine the corresponding energy barrier for each of the isomers found in the global exploration.Considering the stochastic nature of the TM and the fact that the number of steps performed by the walkers is finite, there is always a possibility that some minima cannot be reached even if there exists a trajectory under some energy lid.However, the occurrence of a structural transition from a given structure to another can be explored systematically by means of an appropriate procedure.
In order to extract further information about the equilibrium properties of the isomers with the lowest energy values, additional threshold runs were carried out by defining a particular energy exploration region for selected structures.These target structures were chosen mainly because we expected to find sets of homotops (e.g. the tetrahelices), or due to their apparent polytetrahedral (isotropic or anisotropic) growth behavior.Ten energy lids were generated, starting from the corresponding local minimum, giving a spacing of 0.01 eV/ atom between energy lids.This procedure allowed us to distinguish the corresponding energy barriers between different local minima distributed throughout the hypersurface, obtained from the G42 runs.In this way, several new local minimum configurations of Pd 12 Pt 1 clusters were located, and used as starting configurations in new explorative PES searches (with identical computer simulation conditions), performed to obtain a deeper understanding around the selected local minima structures.This methodology permitted us previously to identify a set of enantiomeric helical structures, 61 detected in an energy region where axial growth of fused tetrahedra (sharing faces or edges) are local minima.In the present work, we extend the analysis choosing a set of computational parameters appropriate for a global PES exploration.

Results and discussion
Table 1 (in the ESI3) lists 25 representatives of the structurally diverse minima found after performing an exhaustive exploration of the Pd 12 Pt 1 cluster PES.The putative GM corresponds to an icosahedral structure, in which the Pt atom occupies the core position (Ih).The CNA indices and their normalized RA provide an indication of the type and degree of order as the energy increases.Several types of pairs characterize different types of order, 62 for instance, an icosahedral packing arrangement is defined by indices 555, whereas 421 and 422 are typical of fcc and hcp order.In our results, the lowest energy icosahedral structure exhibits only 555 (y29%) and 322 (y71%) indices.This structure is an example of isotropic tetrahedral growth, obtained by joining tetrahedra with no preferred stacking.The existence of non-equilateral triangular faces, defined by 322 indices, may be associated with the tetrahedral coordination of this structure while 555 indices define the presence of pentagonal planes as characteristic geometrical elements of this structure.
Relevant higher energy structures are: a fragment of a waistcapped Ino truncated decahedron (Ih-Ino); a fcc BP structure with a small distortion due to the presence of the Pt atom, which was confirmed through local reoptimizations of all the homotops after changing the Pt atom to Pd; linear stacking of tetrahedra forming chiral pairs of right-handed and lefthanded helices (HR and HL); fragments of a double icosahedron (MT-Ih2 and MT-Ih3); a structure composed of two edge-connected pentagonal bipyramids, sharing a capped atom over this edge (S6); a polytetrahedral elongated cluster exhibiting two coupled square bipyramids, resembling a monoclinic structure (S13) and some other structures showing both isotropic and anisotropic polytetrahedral growth.The latter are elongated structures appearing in high energy regions, usually for energy values above 23.365eV/atom.These elongated structures are formed by multiple, twisted or straight strands of atoms.Amorphous structures were also found in our search, as well as a large number of other distorted icosahedral, decahedral and biplanar structures.These have been omitted as they only constituted y1% of the minima found during our PES search.Fig. 1 displays the multiple funnel palm-tree type 16 energy landscape obtained from exploring the semi-empirical potential energy hypersurface.
Once an initial exhaustive global exploration of the Pd 12 Pt 1 cluster PES had been carried out, we then focused on the lowest energy structure, i.e. the putative icosahedral GM structure (see Table 1, 3 Ih), and used this configuration as the initial structure for additional threshold method runs.Fig. 2 shows the different icosahedron-related local-minima found after carrying out the new PES exploration.Four different icosahedron-based configurations were located, and  Our results show that the energy of low-lying structures varies with the number of Pt-Pd bonds, energetically favoring those structures that have a higher number of Pt-Pd bonds (see Table 23).The effect of the varying number of Pt-Pd bonds can be seen by comparing the structures of the homotops Ih 1 and Ih 12 .In Ih 1 we found the maximum number of Pt-Pd bonds (12).On the other hand, Ih A tree or disconnectivity graph 32,[38][39][40][41] was constructed to gain a better picture of the Pd 12 Pt 1 complex hypersurface.In order to construct these diagrams, we grouped all those configurations having identical total energy and CNA indices, i.e. those with essentially the same structure up to permutational symmetry operations, into one basin (see Fig. 1).Fig. 2 shows the tree graph arising from the refined exploration centered on the GM icosahedral structure (Ih 1 ), thus giving a more precise representation of the unique energy landscape of the Pd 12 Pt 1 cluster in terms of structures Ih 1 to Ih 12 within the energy range 23.5 eV/atom to 23.3 eV/atom, which correspond to the lower region of the PES.In this figure, a single transition state connects local minima, but we note that potential structural interconversions, even in this simplistic representation, are affected by the physical nature of the barrier, energetic or entropic.
When the random walker performs a finite number of steps and moves a finite distance, there is a certain probability that the system will not reach a particular target minimum from a given starting minimum, even if a (short enough) pathway below the energy lid between the two neighboring minima exists.The reason for this is that in high-dimensional spaces the likelihood to find such a fast low-energy route between two minima is not only controlled by energetic barriers but also by entropic and/or kinetic barriers taking the width of the transition regions and the internal structure of the basins into account.All these aspects are subsumed under the concept of a generalized barrier discussed elsewhere 63 that can be quantified via the logarithm of the inverse transition probability obtained by counting the number of successful transitions between the minima under consideration.In Fig. 2 we show the 12 different Pd 12 Pt 1 cluster structures (Ih 1 -Ih 12 ), all of which correspond to local minima on the PES.Structural interconversion may be rationalized as a process that transforms one structure into another by crossing the energy barrier connecting the two minima, promoted by physical or chemical processes.
One particularly interesting feature is the grouping of minima on the landscape according to their underlying structure (i.e., modulo the homotopy §).Both the barriers and the energy differences among the helical structures and also among the biplanar structures are quite low compared with the typical differences among structurally distinct isomers.While this is not too surprising as far as the energies are concerned, the rather small barriers, of about 0.05 eV/ atom, are somewhat surprising, since the transformations involve complex atom exchange paths.However, the fact that most atoms can be considered surface atoms is a likely cause for such low-energy transition routes between members of a structural group.As a consequence, we would expect that homotopically related isomers can relatively easily transform into each other, at least in the case of small intermetallic clusters.
Our computational surface energy analysis allowed us to find other relevant structures discussed recently in the scientific literature, 61,64-73 some of which are important for current or potential practical technological applications.
Besides biological examples, in diverse synthesis procedures based on chemical or physical methods, structures with helical symmetry appear as a result of experimental conditions and forces governing interactions between atoms.This motivated an investigation of the geometric structures of some of the elongated Pd 12 Pt 1 clusters existing on the PES, which we published very recently, 61 describing a subset of enantiomeric helical structures (''Bernal spirals'' [64][65][66][67] ) resembling those found experimentally for Au-Ag nanowires. 68Computational restriction means that we have been unable to carry out an exhaustive PES exploration, so several structures were not found during the search but have been subsequently generated (based on symmetry) and then locally minimized.Current results, shown in Fig. 3 as a tree graph, present a set of 14 multitwinned tetrahelices, 7 right-handed and 7 left-handed.The atomic distribution on strands and the helix symmetry for the Pd 12 Pt 1 nanoalloy reduce the number of inequivalent homotops to 7; therefore our current searches provided all the expected right-handed structures and out of 7 of the lefthanded ones.The left-handed structure HL 5 was not reached even though the energy threshold between the lids used in this work had a value of 0.01 eV/atom, i.e., half the spacing used previously. 61Additional structural and energetical details are given in Table 3 (see the ESI 3 ), where the missing left-handed structure was built by hand and locally optimized to provide complete information.Notice that the difference in the barrier heights of the left-and right-handed modifications is due to sampling effects, and the best estimate of the energy barrier is the lower of the two values.71]73 Most of these studies agree with the idea that 13-atom Pd and Pt clusters adopt the BP structure as the ground state configuration, nevertheless, to date there are only a few very recent studies reporting BP structures for bimetallic systems 72 by any of the current theoretical methods to search for the lowest energy structures, i.e., using DFT or semi-empirical potentials.Even though the Gupta potential tends to give 3-D ground state structures, our global exploration allows us to have a better idea of the structural families belonging to a particular PES.
The results reported here also include those for BP Pd 12 Pt 1 structures, which are not the lowest energy structures (see Fig. 1 and Table 4 in the ESI3) but correspond to local minima on the Gupta potential energy hypersurface.The parent Pd 13 BP structure has C 3v symmetry and consists of a plane of 6 atoms forming a triangle and a plane of 7 atoms forming a centered hexagon.Upon substituting one Pd atom by Pt, 4 homotops can be generated for BP Pd 12 Pt 1 (which we have named BP 1 -BP 4 , in order of increasing energy): the disconnectivity diagram for the BP isomers Pd 12 Pt 1 is shown in Fig. 4. The homotops can be described as follows: BP 1 -the Pt lies in the center of the hexagonal plane; BP 2 -the Pt lies in the inner triangle of the triangular plane; BP 3 -the Pt lies on the periphery of the hexagonal plane; BP 4 -the Pt lies in the outer triangle of the triangular plane.Homotops BP 1 , BP 2 and BP 4 are achiral (having C s symmetry), as the Pt atom lies on a mirror plane of symmetry.However, as the Pt atom in homotop BP 3 does not lie on a mirror plane of the parent BP structure, this homotop is chiral, existing as the enantiomeric pair BP 3 9 and BP 3 99 which confirms that chirality is recurrent in bimetallic clusters. 61,74inally, we verified the structural stability and established the ranking of the Pd 12 Pt 1 clusters by energy using more accurate first principles methods (see Section 2.2).For this, a  12) compared to the BP1 structure (9).In terms of the binding energy, the Ih 1 and BP 1 structures have the same value to two decimal places: 2.24 eV.The helical HR 1 and HL 1 structures are found to be rather high in total energy (1.13 eV) compared to the Ih 1 cluster, with a calculated E b value of 2.16 eV for both cases.The relaxed structures and energies are shown in Fig. 5.
Our results show that, even at the higher level of theory, the energetic ordering predicted by the semi-empirical potential and global optimization methods is preserved, and that the conclusions arising from our global searches remain valid.Similar energetic trends have recently been reported by Ferrando et al. for larger Pd-Pt clusters, comparing semiempirical Gupta potential minima and DFT optimized configurations. 75In their work, the putative Gupta minimum became a high-energy isomer at the DFT level; thus highlighting the need to test the reliability of the potential against electronic structure methods even though the computational cost can be rather high.It is important to emphasize that the results reported in ref. 75 correspond to searches for low-lying minima, rather than mapping the energy landscape or calculating the energy barriers to isomerization.Here, we describe a different, insightful process which provides a good description of the full potential energy surface and which allows us to realize e.g. that homotopically related structures are separated by comparatively low barriers.In addition, based on our DFT calculations, we predict that the lowest energy structures reported in this work for these bimetallic clusters can co-exist under appropriate experimental conditions (e.g.temperature, pressure and support type). 68

Conclusions
In this work we have established the appropriateness of using combined optimization (genetic algorithm and the threshold method) schemes in order to explore the complex PES of a bimetallic Pd 12 Pt 1 nanoalloy, as modeled by the Gupta manybody potential.A thorough analysis of the cluster PES allowed us to locate and identify a set of representative structures, where an icosahedral structure was unambiguously found as the corresponding global minimum configuration, for this specific cluster size and composition.The knowledge of global and local minima is a preliminary step in the study of more complex theoretical and experimental analysis of sub-nanoparticles, such as those identified in this research.The construction of tree-type disconnectivity graphs facilitates an understanding of the possible structural interconversion between icosahedron-based geometrical isomers and other multitwinned conformations of the Pd 12 Pt 1 cluster.Our results provide a full energetic description of these local minima structures, their place on the energy landscapes, as well as those configurations closely connected to the icosahedral global minimum structure.The energetics of the low-lying structures predicted by the Gupta potential are essentially confirmed through DFT computations, additionally finding that icosahedral and biplanar structures are also likely to coexist under the same physical conditions.Considering that the modeling of the interatomic interactions between the constituent atoms is a complex task, having realistic model potentials coupled with high-level theoretical calculations will facilitate the PES exploration of larger Pd-Pt bimetallic nanoparticles (also displaying different concentrations of Pd and Pt atoms); thus allowing the prediction of their corresponding ground-state configurations with a good degree of certainty.In this context, the combination of theoretical approaches and state-of-the-art computational tools become crucial in the development of novel tailor-made materials at the nanoscale, underpinning experimental studies.In view of the increasing relevance of sub-nanoparticles in technologically important processes such as heterogeneous catalysis, we expect that our theoretical results will stimulate the realization of experiments on systems such as the one studied here, with new techniques enabling finite temperature measurements to be carried out on real systems, as well as other theoretical studies modeling free and supported binary nanoparticles.

Fig. 1
Fig. 1 Tree (disconnectivity) graph showing the most abundant structures on the Pd 12 Pt 1 potential energy hypersurface.The black horizontal lines indicate energy lid positions; the vertical red lines give the magnitudes of the energy barriers, while the blue vertical lines indicate the last energy lid where this structure was found.
several homotops were also identified, giving a total of 12 structures.These structures can be described starting with the lowest energy structure as follows: icosahedron (Ih 1 ); truncated icosahedron with one waist-down capping Pd atom (Ih 2 ); truncated icosahedron with one waist-up capping Pd atom (Ih 3 ); truncated icosahedron with one top face-capping Pd atom (Ih 4 ) (it should be noted that structures Ih 1 to Ih 4 all have the Pt atom in the central site of the parent icosahedron); homotop of Ih 1 with the Pt atom in the shell; homotop of Ih 2 with Pt at the apex (Ih 6 ); homotops of Ih 3 with Pt in a pentagonal plane (Ih 8 ), at the top (Ih 10 ) and in the waist-up capping position (Ih 11 ); 3 homotops of Ih 4 with Pt at the top (Ih 7 ), in the upper pentagonal plane (Ih 9 ), and as a capping atom (I 12 ).Despite the vast PES analysis carried out by our global optimization strategies, we are aware that many other possible homotops were not reached by our search methods.Consider, for example, the right-handed helical structures where all the 7 homotops (considering symmetry) are found just by changing the Pt atom position.The obtained configurations are grouped together in Table 2 (ESI 3 ) for clarity, where we give additional structural and energetic details.

Fig. 2
Fig. 2 Tree (disconnectivity) graph indicating the possible interconversion pathways between icosahedral-like structures.Structural interconversion between homotops is in the energy range explored by the random walk.The color code is the same as in Fig. 1.

Fig. 3
Fig. 3 Tree (disconnectivity) graph exhibiting the energy landscape region where helical enantiomeric structures are abundant.The color code is the same as in Fig. 1.Fig. 4 Possible homotops BP 1 -BP 4 (and enantiomers BP 3 9 and BP 3 99) of the biplanar structure for the Pd 12 Pt 1 cluster.The tree (disconnectivity) graph shows the structural distribution according to the energy of each configuration.The color code is the same as in Fig. 1.

Fig. 5
Fig. 5 Energetical ordering of selected lowest energy configurations, according to semiempirical (Gupta potential) and DFT calculations.The binding energy (E b ) is in eV/atom.