<figure><img src="image_1.png" /></figure>

HAL
open science

# Intermetallics with sp-d orbital hybridisation: morphologies, stabilities and work-functions of In-Pd particles at the nanoscale

Alexis Front, Clovis Lapointe, Émilie Gaudry

► To cite this version:

Alexis Front, Clovis Lapointe, Émilie Gaudry. Intermetallics with sp-d orbital hybridisation: morphologies, stabilities and work-functions of In-Pd particles at the nanoscale. Nanoscale Horizons, 2024, 10.1039/D3NH00594A. cea-04622757

HAL Id: cea-04622757

https://cea.hal.science/cea-04622757

Submitted on 24 Jun 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Cite this: DOI: 00.0000/xxxxxxxxxx

Received Date

Accepted Date

DOI: 00.0000/xxxxxxxxxx

# Intermetallics with sp-d orbital hybridisation: morphologies, stabilities and work-functions of In-Pd particles at the nanoscale

Alexis Front*<sup>a,b</sup>, Clovis Lapointe<sup>a,c</sup> and Émilie Gaudry*<sup>a</sup>

The field of intermetallic catalysts, alloying a p-block and a transition metal to form a pM-TM bimetallic, is experiencing robust growth, emerging as a vibrant frontier in catalysis research. Although such materials are increasingly used in the form of nanoparticles, a precise description of their atomic arrangements at the nanoscale remains scarce. Based on the In-Pd binary as a typical pM-TM system, we performed Density Functional Theory calculations to investigate the morphologies, relative stabilities and electronic properties of 24 Å and 36 Å nanoparticles built from the $In_3Pd_2$, InPd and $InPd_3$ compounds. Wulff equilibrium structures are compared to other ordered and disordered structures. Surface energies are computed to discuss their thermodynamic stability, while work-functions are calculated to examine their electronic structures. For any compound, increasing the size leads to the stabilisation of Wulff polyhedra, which are found to offer smaller surface energies than non-crystalline and chemically disordered structures. Disordered $In_3Pd_2$ and InPd, nanoparticles show a tendency towards amorphisation, owing to repulsive short In-In bonds. Tuning nanoparticles work-functions can be achieved through the control of the surface structure and composition, by virtue of the roughly linear correlation found between the surface composition and the work-function but which nevertheless includes a certain number of outliers. This work paves the way to a rationalisation of both structural and electronic properties of pM-TM nanoparticles.

## Conceptual Insights

Intermetallic nanoparticles, combining p-block and transition metals to form pM-TM binary alloys, are of paramount importance across a broad spectrum of applications, underscoring the crucial significance of studying their nanoscale structure. Through the In-Pd system, novel insights have been uncovered, shedding new light on the structure-property relationships within this type of pM-TM intermetallic nanoparticles. Based on density functional theory calculations, regarding structures of $In_3Pd_2$, InPd and $InPd_3$ nanoparticles, within sizes ranging from 25 Å to 34 Å, we show that non crystalline morphologies are not likely, Wulff-type morphologies being thermodynamically more favourable for ordered systems, the larger the size, the better. Moreover, our computations reveal that disordered In-rich nanoparticles tend to undergo amorphisation, driven by destabilising In-In interactions, while the crystalline character of nanoparticles with a large Pd content is kept. Concerning nanoparticle properties, using an original computational method, we bring to light a roughly linear correlation between the surface composition and the work function, although it nonetheless incorporates a certain number of outliers. Amorphous particles have a lower work function compared to crystalline particles, attributable to their surface composition enriched in indium. This study establishes a foundation for systematically understanding the structural and electronic characteristics of pM-TM nanoparticles.

<footnote><sup>a</sup> Université de Lorraine, CNRS, Institut Jean Lamour, UMR 7198, Campus Artem, 2 allée André Guinier, F-54011, Nancy, France; <sup>b</sup> Department of Chemistry and Materials Science, Aalto University, 02150 Espoo, Finland; Department of Applied</footnote>
Cite this: DOI: 00.0000/xxxxxxxxxx

Received Date

Accepted Date

DOI: 00.0000/xxxxxxxxxx

# Intermetallics with sp-d orbital hybridisation: morphologies, stabilities and work-functions of In-Pd particles at the nanoscale

Alexis Front*<sup>a,b</sup>, Clovis Lapointe<sup>a,c</sup> and Émilie Gaudry*<sup>a</sup>

The field of intermetallic catalysts, alloying a p-block and a transition metal to form a pM-TM bimetallic, is experiencing robust growth, emerging as a vibrant frontier in catalysis research. Although such materials are increasingly used in the form of nanoparticles, a precise description of their atomic arrangements at the nanoscale remains scarce. Based on the In-Pd binary as a typical pM-TM system, we performed Density Functional Theory calculations to investigate the morphologies, relative stabilities and electronic properties of 24 Å and 36 Å nanoparticles built from the $In_3Pd_2$, InPd and $InPd_3$ compounds. Wulff equilibrium structures are compared to other ordered and disordered structures. Surface energies are computed to discuss their thermodynamic stability, while work-functions are calculated to examine their electronic structures. For any compound, increasing the size leads to the stabilisation of Wulff polyhedra, which are found to offer smaller surface energies than non-crystalline and chemically disordered structures. Disordered $In_3Pd_2$ and InPd, nanoparticles show a tendency towards amorphisation, owing to repulsive short In-In bonds. Tuning nanoparticles work-functions can be achieved through the control of the surface structure and composition, by virtue of the roughly linear correlation found between the surface composition and the work-function but which nevertheless includes a certain number of outliers. This work paves the way to a rationalisation of both structural and electronic properties of pM-TM nanoparticles.

## Conceptual Insights

Intermetallic nanoparticles, combining p-block and transition metals to form pM-TM binary alloys, are of paramount importance across a broad spectrum of applications, underscoring the crucial significance of studying their nanoscale structure. Through the In-Pd system, novel insights have been uncovered, shedding new light on the structure-property relationships within this type of pM-TM intermetallic nanoparticles. Based on density functional theory calculations, regarding structures of $In_3Pd_2$, InPd and $InPd_3$ nanoparticles, within sizes ranging from 25 Å to 34 Å, we show that non crystalline morphologies are not likely, Wulff-type morphologies being thermodynamically more favourable for ordered systems, the larger the size, the better. Moreover, our computations reveal that disordered In-rich nanoparticles tend to undergo amorphisation, driven by destabilising In-In interactions, while the crystalline character of nanoparticles with a large Pd content is kept. Concerning nanoparticle properties, using an original computational method, we bring to light a roughly linear correlation between the surface composition and the work function, although it nonetheless incorporates a certain number of outliers. Amorphous particles have a lower work function compared to crystalline particles, attributable to their surface composition enriched in indium. This study establishes a foundation for systematically understanding the structural and electronic characteristics of pM-TM nanoparticles.

<footnote><sup>a</sup> Université de Lorraine, CNRS, Institut Jean Lamour, UMR 7198, Campus Artem, 2 allée André Guinier, F-54011, Nancy, France; <sup>b</sup> Department of Chemistry and Materials Science, Aalto University, 02150 Espoo, Finland; Department of Applied</footnote>
# 1 Introduction

The development of technologically relevant nanoparticles (NPs) has seen unprecedented growth in recent years. Initially used by artisans, mainly for their optical properties, albeit without knowledge of their nature,¹⁻³ they are nowadays used to develop many different applications, in biomedicine,⁴⁻⁶ energy,⁷⁻⁹ and information technologies,¹⁰,¹¹ with the most widespread applications likely being in catalysis. In the latter case, mono-metallic precious and transition metal NPs, have played a major role,¹²⁻¹⁵ demonstrated, among others, by the relevance of gold NPs in CO oxidation,¹⁶ of copper-based NPs in methanol production from syngas (CO + H₂),¹⁷ of silver NPs in solar-driven photocatalytic water splitting, CO₂ reduction and degradation of organic pollutants,¹⁸ and also of platinum, iridium and ruthenium based NPs in hydrogen evolution reaction (HER), oxygen reduction reaction (ORR) and oxygen evolution reaction (OER).¹⁹ This is attributed to the ability of metals to activate chemical bonds through d-electrons or d-empty orbitals,²⁰ and also to the quantitative understanding that has been established in transition-metal catalysis over recent years.²¹,²²

Elemental transition metals however struggle to satisfy the numerous and sometimes contradictory needs of high-performance catalysts, i.e. high activity, tractable selectivity, low cost and good stability. Tailoring atomic arrangements and chemical compositions are key approaches to circumvent the previous limitations, generally achieved by alloying, faceting, nano-structuring, as well as selecting crystal structures and surface orientations.²³⁻²⁵ Thus, many bimetallic NPs have been developed over the last few years, generally by associating two precious metals or a noble metal and a first-row transition metal. The so-called “synergistic” effects between the two metals is known to greatly enhance the catalyst efficiency in many cases. Hence, Pd-Au alloys are likely among the most widely used solid-solution alloys in all fields of catalysis (Ref.²⁵ and references therein), while Pt-based nanoalloys, especially in the Pt-Ni system, are considered to be one of the most promising catalysts for the oxygen reduction reaction (ORR) – bottleneck restricting the large-scale commercialisation of fuel cell vehicles²⁶ – able to improve by a factor 4 to 10 the catalytic performance of ORR in comparison to commercial Pt/C.²⁷,²⁸

The specific performances of bimetallic catalysts is typically attributed to electronic and geometric effects, including orbital hybridisation and active site isolation. Up to now, most theoretical investigations have relied on infinite slabs, while experimental catalysts are generally synthesised in the form of nanoparticles. At the nanoscale, mixing two late transition metals leads a large variety of configurations, including core-shell, segregated, multi-shell, ordered or solid-solution alloys.²⁹ The atomic arrangements can greatly differ from the ones resulting from bulk truncation, especially when the composition range of a given phase is

large and its formation enthalpy is low. It is for instance the case of Pt-Ni nanoalloys ($|\Delta H_f| < 0.1$ eV/at.),³⁰ for which icosahedral core-shell and multishell morphologies have been predicted to be likely for 147-923 atom NPs.³¹

Compared to conventional d-d bimetallic NPs, the knowledge in systems combining transition and p-block metals – e.g. Al, Zn, Ga, Cd, In and Sn within the first periods of the periodic table – is little and is currently quickly developing. Among them, In-Pd phases have been identified as relevant catalysts for several reactions. This is especially the case for (i) methanol steam reforming, for which InₓPdᵧ/In₂O₃ aerogels exhibit excellent CO₂ selectivities (99%),³² (ii) for methanol synthesis, for which In-Pd NPs show higher CH₃OH rates (70%) and selectivities (80% at 270°C) compared to conventional Cu/ZnO/Al₂O₃,³³,³⁴ and (iii) for hydrogenation / dehydrogenation reactions, for which the addition of In to Pd improves the catalytic performances, likely attributed to site isolation.³⁵⁻³⁹ The formation of In-Pd intermetallics can also reduce the catalyst coking³⁷,³⁸ and poisoning.⁴⁰ Such performances strongly depend on the catalyst’s chemical composition and surface structure. Indeed, recent studies highlight the better catalytic performances in acetylene hydrogenation of bimetallic In-Pd at the In:Pd ratio of 0.8,³⁷ and the higher selectivity (92% vs 21% at 90 °C)⁴¹ of 1:1 InPd NPs than of InPd₃ NPs, attributed to the different type of mainly exposed facets – (110) and (111) for 1:1 InPd and InPd₃, respectively. Thus, the controlled synthesis of pM-TM (p-block metal - transition metal) NPs is crucial to design catalysts with optimal performances.

The rich structural chemistry of pM-TM systems, typically attributed to the large formation enthalpy of 1:1 bulk compounds ($|\Delta H_f| > 0.1$ eV/at. in many cases)³⁰,⁴² yields to many intermetallics in the bulk phase diagram,⁴³ and thus is expected to lead to a substantial variety of atomic arrangements at the nanoscale. Within the In-Pd system, the large compositional range of the 1:1 InPd phase, extending from 45 at.%Pd to 61.5 at.%Pd, i.e. much larger than other In-Pd binaries,⁴⁴ likely ease the experimental synthesis of InPd NPs, with sizes down to 10 nm. Other single-phase NPs have however recently been synthesised, like InPd₃ (size 5 nm)³⁹, In₃Pd₂ (size 90 nm)⁴⁵, In₇Pd₃ (size in the range 60-90 nm)³⁹,⁴⁵. Moreover, for particles sizes below 3 nm, core-shell structures have been identified, with a Pd core and a 1:1 InPd shell, when the In:Pd atomic ratio is 0.8.³⁵ Ordered 1:1 InPd particles, with an unalloyed In shell, have also been formed, with a rather high In:Pd ratio (value = 2).³⁵ However, no disordered crystalline structure has ever been observed, nor any non-crystalline structures, in contrast to the literature on d-d nanocrystals.⁴⁶⁻⁴⁸

Overall, the development of a synthetic approach for ordered intermetallic nanocrystals is challenging.⁴⁹ Thus, numerical methods are crucial to understand, explain and predict nanoalloys morphologies, structures and properties.⁵⁰,⁵¹ Current computational resources combined with Density Functional Theory (DFT) are now able to calculate key quantities, like the nanoparticles surface energy, one of the most important quantity to discuss their thermodynamic stability,³¹,⁵²,⁵³ which can be compared to experimental measurements.⁵⁴ Because chemisorption typically involves the Fermi level of the catalysts and the

physics, Aalto university, Espoo, Finland; c Université Paris-Saclay, CEA, Service de Recherches de Métallurgie Physique, 91191, Gif-sur-Yvette, France
E-mails: alexis.front@aalto.fi; emilie.gaudry@univ-lorraine.fr

† Electronic Supplementary Information (ESI) available: See DOI: 00.0000/0000000.
# 1 Introduction

The development of technologically relevant nanoparticles (NPs) has seen unprecedented growth in recent years. Initially used by artisans, mainly for their optical properties, albeit without knowledge of their nature,¹⁻³ they are nowadays used to develop many different applications, in biomedicine,⁴⁻⁶ energy,⁷⁻⁹ and information technologies,¹⁰,¹¹ with the most widespread applications likely being in catalysis. In the latter case, mono-metallic precious and transition metal NPs, have played a major role,¹²⁻¹⁵ demonstrated, among others, by the relevance of gold NPs in CO oxidation,¹⁶ of copper-based NPs in methanol production from syngas (CO + H₂),¹⁷ of silver NPs in solar-driven photocatalytic water splitting, CO₂ reduction and degradation of organic pollutants,¹⁸ and also of platinum, iridium and ruthenium based NPs in hydrogen evolution reaction (HER), oxygen reduction reaction (ORR) and oxygen evolution reaction (OER).¹⁹ This is attributed to the ability of metals to activate chemical bonds through d-electrons or d-empty orbitals,²⁰ and also to the quantitative understanding that has been established in transition-metal catalysis over recent years.²¹,²²

Elemental transition metals however struggle to satisfy the numerous and sometimes contradictory needs of high-performance catalysts, i.e. high activity, tractable selectivity, low cost and good stability. Tailoring atomic arrangements and chemical compositions are key approaches to circumvent the previous limitations, generally achieved by alloying, faceting, nano-structuring, as well as selecting crystal structures and surface orientations.²³⁻²⁵ Thus, many bimetallic NPs have been developed over the last few years, generally by associating two precious metals or a noble metal and a first-row transition metal. The so-called “synergistic” effects between the two metals is known to greatly enhance the catalyst efficiency in many cases. Hence, Pd-Au alloys are likely among the most widely used solid-solution alloys in all fields of catalysis (Ref.²⁵ and references therein), while Pt-based nanoalloys, especially in the Pt-Ni system, are considered to be one of the most promising catalysts for the oxygen reduction reaction (ORR) – bottleneck restricting the large-scale commercialisation of fuel cell vehicles²⁶ – able to improve by a factor 4 to 10 the catalytic performance of ORR in comparison to commercial Pt/C.²⁷,²⁸

The specific performances of bimetallic catalysts is typically attributed to electronic and geometric effects, including orbital hybridisation and active site isolation. Up to now, most theoretical investigations have relied on infinite slabs, while experimental catalysts are generally synthesised in the form of nanoparticles. At the nanoscale, mixing two late transition metals leads a large variety of configurations, including core-shell, segregated, multi-shell, ordered or solid-solution alloys.²⁹ The atomic arrangements can greatly differ from the ones resulting from bulk truncation, especially when the composition range of a given phase is

large and its formation enthalpy is low. It is for instance the case of Pt-Ni nanoalloys ($|\Delta H_f| < 0.1$ eV/at.),³⁰ for which icosahedral core-shell and multishell morphologies have been predicted to be likely for 147-923 atom NPs.³¹

Compared to conventional d-d bimetallic NPs, the knowledge in systems combining transition and p-block metals – e.g. Al, Zn, Ga, Cd, In and Sn within the first periods of the periodic table – is little and is currently quickly developing. Among them, In-Pd phases have been identified as relevant catalysts for several reactions. This is especially the case for (i) methanol steam reforming, for which InₓPdᵧ/In₂O₃ aerogels exhibit excellent CO₂ selectivities (99%),³² (ii) for methanol synthesis, for which In-Pd NPs show higher CH₃OH rates (70%) and selectivities (80% at 270°C) compared to conventional Cu/ZnO/Al₂O₃,³³,³⁴ and (iii) for hydrogenation / dehydrogenation reactions, for which the addition of In to Pd improves the catalytic performances, likely attributed to site isolation.³⁵⁻³⁹ The formation of In-Pd intermetallics can also reduce the catalyst coking³⁷,³⁸ and poisoning.⁴⁰ Such performances strongly depend on the catalyst’s chemical composition and surface structure. Indeed, recent studies highlight the better catalytic performances in acetylene hydrogenation of bimetallic In-Pd at the In:Pd ratio of 0.8,³⁷ and the higher selectivity (92% vs 21% at 90 °C)⁴¹ of 1:1 InPd NPs than of InPd₃ NPs, attributed to the different type of mainly exposed facets – (110) and (111) for 1:1 InPd and InPd₃, respectively. Thus, the controlled synthesis of pM-TM (p-block metal - transition metal) NPs is crucial to design catalysts with optimal performances.

The rich structural chemistry of pM-TM systems, typically attributed to the large formation enthalpy of 1:1 bulk compounds ($|\Delta H_f| > 0.1$ eV/at. in many cases)³⁰,⁴² yields to many intermetallics in the bulk phase diagram,⁴³ and thus is expected to lead to a substantial variety of atomic arrangements at the nanoscale. Within the In-Pd system, the large compositional range of the 1:1 InPd phase, extending from 45 at.%Pd to 61.5 at.%Pd, i.e. much larger than other In-Pd binaries,⁴⁴ likely ease the experimental synthesis of InPd NPs, with sizes down to 10 nm. Other single-phase NPs have however recently been synthesised, like InPd₃ (size 5 nm)³⁹, In₃Pd₂ (size 90 nm)⁴⁵, In₇Pd₃ (size in the range 60-90 nm)³⁹,⁴⁵. Moreover, for particles sizes below 3 nm, core-shell structures have been identified, with a Pd core and a 1:1 InPd shell, when the In:Pd atomic ratio is 0.8.³⁵ Ordered 1:1 InPd particles, with an unalloyed In shell, have also been formed, with a rather high In:Pd ratio (value = 2).³⁵ However, no disordered crystalline structure has ever been observed, nor any non-crystalline structures, in contrast to the literature on d-d nanocrystals.⁴⁶⁻⁴⁸

Overall, the development of a synthetic approach for ordered intermetallic nanocrystals is challenging.⁴⁹ Thus, numerical methods are crucial to understand, explain and predict nanoalloys morphologies, structures and properties.⁵⁰,⁵¹ Current computational resources combined with Density Functional Theory (DFT) are now able to calculate key quantities, like the nanoparticles surface energy, one of the most important quantity to discuss their thermodynamic stability,³¹,⁵²,⁵³ which can be compared to experimental measurements.⁵⁴ Because chemisorption typically involves the Fermi level of the catalysts and the

<footnote>physics, Aalto university, Espoo, Finland; c Université Paris-Saclay, CEA, Service de Recherches de Métallurgie Physique, 91191, Gif-sur-Yvette, France
E-mails: alexis.front@aalto.fi; emilie.gaudry@univ-lorraine.fr</footnote>

<footnote>† Electronic Supplementary Information (ESI) available: See DOI: 00.0000/0000000.</footnote>
HOMOs / LUMOs (highest occupied / lowest unoccupied molecular orbitals) of the reactants, the work-function is one important feature for catalysis.⁵⁵⁻⁵⁷ In this work, by focusing on NPs with the 1:1 InPd composition, as well as on one composition in the Pd-rich side (InPd₃) and on one composition in the In-rich side (In₃Pd₂), we show that ordered structures are stable, at least when compared to non-crystalline structures. We also highlight the amorphisation of disordered structures, driven by destabilising In-In interactions. Finally, trends between the NPs work-functions and their morphologies, sizes and chemical compositions are examined, thus opening a way to design NPs with specific properties.

# 2 Results and discussion

## 2.1 The In-Pd system

The bulk In-Pd system presents a rich phase diagram, with six intermetallic compounds at room temperature:⁴⁴ In₇Pd₃ (Im₃m space group, cI40), In₃Pd₂ (P₃m1 space group, hP5), α-InPd₃ (I4/mmm space group, tI8), α-InPd₂ (Pnma space group, oP12), In₃Pd₅ (Pbam space group, oP16), InPd (Pm₃m space group, cP2). Beside the 1:1 InPd compound, we have considered in this work one In-rich compound (In₃Pd₂) and one Pd-rich compound (InPd₃). Depending on the experimental conditions, the InPd₃ compound can be synthesised with several crystal structures, including AuCu₃-type (cubic, Pm₃m space group, cP4),⁵⁸ ZrAl₃-type (tetragonal, I4/mmm space group, tl16),⁵⁸ CuInPt₂-type (tetragonal, I4/mmm space group, tP4).⁵⁹ In this work, we focus on three In-Pd compounds (Tab. 1): cubic InPd of CsCl-type and cubic InPd₃ of AuCu₃-type, in connection with the extensive investigations of nanoalloys with cubic structures.⁶⁰⁻⁶² We also consider the In₃Pd₂ phase, that crystallises in the trigonal system, for which investigations at the nanoscale are nonexistent, as far as we know.

<table><caption>Table 1 Space groups, cell parameters and formation enthalpies of the In-Pd compounds considered in this study.</caption><thead><tr><th></th><th>Space group</th><th>a / c (Å)</th><th>ΔH<sub>f</sub> (eV/at.)</th><th></th></tr></thead><tbody><tr><td rowspan="3">In<sub>3</sub>Pd<sub>2</sub></td><td rowspan="3">P3m1</td><td>4.642 / 5.586</td><td>-0.476</td><td>PBE (our work)</td></tr><tr><td>4.642 / 5.586</td><td>-0.464</td><td>PBE<sup>30</sup></td></tr><tr><td>4.535 / 5.512</td><td>-0.59</td><td>exp. <sup>59,63</sup></td></tr><tr><td rowspan="3">InPd</td><td rowspan="3">Pm&#773;3m</td><td>3.31 / -</td><td>-0.489</td><td>PBE (our work)</td></tr><tr><td>3.31 / -</td><td>-0.497</td><td>PBE<sup>30</sup></td></tr><tr><td>3.250 / -</td><td>-0.72</td><td>exp. <sup>59,63</sup></td></tr><tr><td rowspan="3">InPd<sub>3</sub></td><td rowspan="3">Pm&#773;3m</td><td>4.045 / -</td><td>-0.4396</td><td>PBE (our work)</td></tr><tr><td>4.045 / -</td><td>-0.438</td><td>PBE<sup>30</sup></td></tr><tr><td>4.027 / -</td><td>-0.554</td><td>exp. <sup>58,64</sup></td></tr></tbody></table>

The calculated electronic density of states of the three compounds is represented in Fig. 1. The low energy range (-10 eV to -5 eV) is mainly formed by In-states, with a contribution of Pd-states. As previously mentioned by Refs.⁶⁵ and⁶⁶ or GaPd₂ and In₃Pd₅, this region reflects the chemical binding in the compounds. The remaining part shows one strong contribution from Pd d-states, in a form of a band, located above 4 eV, whose bandwidth increases with the Pd content. Indeed, the isolation of Pd sites in In₃Pd₂ and InPd, more pronounced as the In concentration increases, results in a narrow d-band, the hybridisation be-

tween Pd-d states being weak. In contrast, the d-d hybridisation in InPd₃ leads to a large d-band (roughly two times as large as in In₃Pd₂ and InPd). Such differences in the electronic structure may be at the origin of the distinct behaviour of In₃Pd₂ and InPd NPs, in comparison to InPd₃ NPs.

<figure><img src="image_1.png" /><figcaption>Fig. 1 Density of states (DOS) of In₃Pd₂ (left), InPd (middle) and InPd₃ (right) bulk compounds. The total DOS of each bulk compound is given (in grey), along with the In-sp and Pd-d contributions (in brown and blue, respectively).</figcaption></figure>

## 2.2 Surface energies and morphologies of NPs based on the Wulff construction

Overall, the equilibrium shape of free nanoparticles are known to be given by the Wulff construction, assuming the bulk crystal structure, and based on the ratios between the compound surface energies.⁶⁷ Energies of several low-index surfaces (Tab. S1) have been computed as a function of the Pd chemical potential (Tab. S2) and are shown in Fig. 2. The three sub-figures correspond to energies of In₃Pd₂ (left), InPd (middle) and InPd₃ (right) low-index surfaces. A few additional orientations have been considered, but are not presented in Fig. 2 (In₃Pd₂(2\bar{1}0), InPd(102) and InPd₃(112)). In the following, we define the In- and the Pd-rich limit by the situations with Δμ<sub>Pd</sub> = 0 eV and Δμ<sub>In</sub> = 0 eV, respectively. This is valid for the three compositions. Nevertheless, it should be mentioned that the value of Δμ<sub>Pd</sub> at the In-rich limit (Δμ<sub>In</sub> = 0 eV) depends on the considered compound. It is Δμ<sub>Pd</sub> = -0.978 eV, Δμ<sub>Pd</sub> = -0.586 eV and Δμ<sub>Pd</sub> = -1.190 eV in the case of InPd, InPd₃ and In₃Pd₂, respectively.

On the whole, pure Pd-terminated surfaces are calculated to present a rather high surface energy, (label 12 for InPd in Fig. 2, not represented for InPd₃), while stable surfaces are found among mixed terminations. It is the case for InPd₃(100) (label 15 in Fig. 2), with a particularly low surface energy at the In-rich limit (γ = 0.491 J/m²) that increases at the Pd-rich limit (γ = 1.072 J/m²). The latter value is similar to the one found for InPd₃(111) (1.046 J/m²) for Δμ<sub>Pd</sub> = 0. Indeed, both InPd₃(111) and InPd₃(100) terminations contain In and Pd atoms. The InPd₃(111) topmost plane is stoichiometric and shows a higher
HOMOs / LUMOs (highest occupied / lowest unoccupied molecular orbitals) of the reactants, the work-function is one important feature for catalysis.⁵⁵⁻⁵⁷ In this work, by focusing on NPs with the 1:1 InPd composition, as well as on one composition in the Pd-rich side (InPd₃) and on one composition in the In-rich side (In₃Pd₂), we show that ordered structures are stable, at least when compared to non-crystalline structures. We also highlight the amorphisation of disordered structures, driven by destabilising In-In interactions. Finally, trends between the NPs work-functions and their morphologies, sizes and chemical compositions are examined, thus opening a way to design NPs with specific properties.

## 2 Results and discussion

### 2.1 The In-Pd system

The bulk In-Pd system presents a rich phase diagram, with six intermetallic compounds at room temperature:⁴⁴ In₇Pd₃ (Im₃m space group, cI40), In₃Pd₂ (P₃m1 space group, hP5), α-InPd₃ (I4/mmm space group, tI8), α-InPd₂ (Pnma space group, oP12), In₃Pd₅ (Pbam space group, oP16), InPd (Pm₃m space group, cP2). Beside the 1:1 InPd compound, we have considered in this work one In-rich compound (In₃Pd₂) and one Pd-rich compound (InPd₃). Depending on the experimental conditions, the InPd₃ compound can be synthesised with several crystal structures, including AuCu₃-type (cubic, Pm₃m space group, cP4),⁵⁸ ZrAl₃-type (tetragonal, I4/mmm space group, tl16),⁵⁸ CuInPt₂-type (tetragonal, I4/mmm space group, tP4).⁵⁹ In this work, we focus on three In-Pd compounds (Tab. 1): cubic InPd of CsCl-type and cubic InPd₃ of AuCu₃-type, in connection with the extensive investigations of nanoalloys with cubic structures.⁶⁰⁻⁶² We also consider the In₃Pd₂ phase, that crystallises in the trigonal system, for which investigations at the nanoscale are nonexistent, as far as we know.

<table><caption>Table 1 Space groups, cell parameters and formation enthalpies of the In-Pd compounds considered in this study.</caption><thead><tr><th></th><th>Space group</th><th>a / c (Å)</th><th>ΔH<sub>f</sub> (eV/at.)</th><th></th></tr></thead><tbody><tr><td rowspan="3">In<sub>3</sub>Pd<sub>2</sub></td><td rowspan="3">P3m1</td><td>4.642 / 5.586</td><td>-0.476</td><td>PBE (our work)</td></tr><tr><td>4.642 / 5.586</td><td>-0.464</td><td>PBE<sup>30</sup></td></tr><tr><td>4.535 / 5.512</td><td>-0.59</td><td>exp. <sup>59,63</sup></td></tr><tr><td rowspan="3">InPd</td><td rowspan="3">Pm&#773;3m</td><td>3.31 / -</td><td>-0.489</td><td>PBE (our work)</td></tr><tr><td>3.31 / -</td><td>-0.497</td><td>PBE<sup>30</sup></td></tr><tr><td>3.250 / -</td><td>-0.72</td><td>exp. <sup>59,63</sup></td></tr><tr><td rowspan="3">InPd<sub>3</sub></td><td rowspan="3">Pm&#773;3m</td><td>4.045 / -</td><td>-0.4396</td><td>PBE (our work)</td></tr><tr><td>4.045 / -</td><td>-0.438</td><td>PBE<sup>30</sup></td></tr><tr><td>4.027 / -</td><td>-0.554</td><td>exp. <sup>58,64</sup></td></tr></tbody></table>

The calculated electronic density of states of the three compounds is represented in Fig. 1. The low energy range (-10 eV to -5 eV) is mainly formed by In-states, with a contribution of Pd-states. As previously mentioned by Refs.⁶⁵ and⁶⁶ or GaPd₂ and In₃Pd₅, this region reflects the chemical binding in the compounds. The remaining part shows one strong contribution from Pd d-states, in a form of a band, located above 4 eV, whose bandwidth increases with the Pd content. Indeed, the isolation of Pd sites in In₃Pd₂ and InPd, more pronounced as the In concentration increases, results in a narrow d-band, the hybridisation be-

tween Pd-d states being weak. In contrast, the d-d hybridisation in InPd₃ leads to a large d-band (roughly two times as large as in In₃Pd₂ and InPd). Such differences in the electronic structure may be at the origin of the distinct behaviour of In₃Pd₂ and InPd NPs, in comparison to InPd₃ NPs.

<figure><img src="image_3.png" /><figcaption>Fig. 1 Density of states (DOS) of In₃Pd₂ (left), InPd (middle) and InPd₃ (right) bulk compounds. The total DOS of each bulk compound is given (in grey), along with the In-sp and Pd-d contributions (in brown and blue, respectively).</figcaption></figure>

### 2.2 Surface energies and morphologies of NPs based on the Wulff construction

Overall, the equilibrium shape of free nanoparticles are known to be given by the Wulff construction, assuming the bulk crystal structure, and based on the ratios between the compound surface energies.⁶⁷ Energies of several low-index surfaces (Tab. S1) have been computed as a function of the Pd chemical potential (Tab. S2) and are shown in Fig. 2. The three sub-figures correspond to energies of In₃Pd₂ (left), InPd (middle) and InPd₃ (right) low-index surfaces. A few additional orientations have been considered, but are not presented in Fig. 2 (In₃Pd₂(2\bar{1}0), InPd(102) and InPd₃(112)). In the following, we define the In- and the Pd-rich limit by the situations with Δμ<sub>Pd</sub> = 0 eV and Δμ<sub>In</sub> = 0 eV, respectively. This is valid for the three compositions. Nevertheless, it should be mentioned that the value of Δμ<sub>Pd</sub> at the In-rich limit (Δμ<sub>In</sub> = 0 eV) depends on the considered compound. It is Δμ<sub>Pd</sub> = -0.978 eV, Δμ<sub>Pd</sub> = -0.586 eV and Δμ<sub>Pd</sub> = -1.190 eV in the case of InPd, InPd₃ and In₃Pd₂, respectively.

On the whole, pure Pd-terminated surfaces are calculated to present a rather high surface energy, (label 12 for InPd in Fig. 2, not represented for InPd₃), while stable surfaces are found among mixed terminations. It is the case for InPd₃(100) (label 15 in Fig. 2), with a particularly low surface energy at the In-rich limit (γ = 0.491 J/m²) that increases at the Pd-rich limit (γ = 1.072 J/m²). The latter value is similar to the one found for InPd₃(111) (1.046 J/m²) for Δμ<sub>Pd</sub> = 0. Indeed, both InPd₃(111) and InPd₃(100) terminations contain In and Pd atoms. The InPd₃(111) topmost plane is stoichiometric and shows a higher
<figure><img src="image_2.png" /><figcaption>Fig. 2 Surface energies of In₃Pd₂(hkℓ) (left), InPd(hkℓ) (middle) and InPd₃(hkℓ) (right) as a function of Δμₚd = μₚd - μᵇₚd where μₚd is the chemical potential of Pd in the compound and μᵇₚd is the cohesive energy of elemental Pd bulk. Pure Pd, pure In and mixed terminations are shown with solid, dashed and dash-dotted lines, respectively. Only the most stable surfaces are represented. Colour code for surface terminations: Pd = blue and In = brown.</figcaption></figure>

atomic density than the mixed termination of InPd₃(100) – 0.15 at/Å² vs. 0.12 at/Å². But the In content of the mixed InPd₃(100) termination is larger (0.06 In at./Å², vs. 0.04 In at./Å² for the stoichiometric surface), thus leading to a stabilising effect. In- deed, the presence of In atoms on the surface is less unfavourable than that of Pd. Similar surface energies are also calculated for InPd(100), InPd(110) and InPd(111) at the In-rich limit, and for InPd(111) and InPd(110) at the Pd-rich limit. The situa- tion is rather complex for In₃Pd₂, with relative stabilities calcu- lated to rank as γ(101) < γ(100) < γ(110) < γ(001) in the In-rich and γ(001) < γ(100) < γ(110) < γ(101) in the Pd-rich side.

Based on the surface energies at the In- and Pd-rich limits, and using the Wulff approach, two types of NPs are built (Fig. 3). Their morphologies strongly depend on the chemical potentials. The InPd₃ NP exhibits a cubic shape in the In-rich side and a trun- cated octahedral (TOh) morphology – typically the Wulff polyhe- dra of TM alloys – in the Pd-rich side. Because the stoichiomet- ric InPd(110) surface presents a low energy, independent of the chemical potentials, the InPd NP morphology is similar at both In- and Pd-rich limits, characterised by large (110) facets. The main difference between the two cases lies in small (100) facets ob- served at the In-rich side, due to the low InPd(100) surface energy at Δμ<sub>In</sub> = 0. The In₃Pd₂ NP displays a cylindrical shape at the Pd- rich limit, mainly due to the stability of the In₃Pd₂(001) surface for Δμ<sub>Pd</sub> = 0. In the In-rich side, the stability of the In₃Pd₂(110) surface leads to a shape with steeper facets. In the following,

Wulff Polyhedra, built using most stable facets at the In-rich and Pd-rich limits (Fig. 2, Tab. S1, section S1), are labelled WP<sub>In</sub> and WP<sub>Pd</sub>, respectively.

## 2.3 Stability of bimetallic In-Pd nanoparticles

The calculation of surface energies requires an infinite reservoir of particles, which becomes questionable when the size of the NPs is not large enough. Therefore, we build two types of NPs: on the one hand, NPs with a diameter between 22 Å and 28 Å (300 to 500 atoms, called “small” NPs in the following), and on the other hand, NPs with a diameter larger than 36 Å (made of around 1000 to 1200 atoms, called “large” NPs in the following).

Moreover, we also consider in this work other geometries, that can offer dense fcc{111} and fcc{100} facets when dealing with fcc metals, but whose relevance to the case of non fcc In-Pd NPs mostly relies on their spherical shape. We focused on icosahedra- dron (Ih), decahedron (Dh), octahedron (Oh) and cuboctahedron (Cubo) shapes (Fig. 4). These morphologies have been chosen because, for metal NPs, (i) Non-crystalline Ih and Dh geometries are typically the ones of small to intermediate sized NPs,<sup>68,69</sup> while crystalline Oh and Cubo geometries are commonly in com- petition with the Wulff polyhedron for large size NPs. Below, not to be confused with the conventional Cubo and Oh shapes of L1<sub>2</sub> bimetallic NPs (as is InPd<sub>3</sub>), we rename Oh by Oh* and Cubo by Cubo*, to indicate that * correspond to In<sub>3</sub>Pd<sub>2</sub> or InPd. Both or- dered and chemically disordered NPs are generated. ll structures
<figure><img src="image_4.png" /><figcaption>Fig. 2 Surface energies of In₃Pd₂(hkℓ) (left), InPd(hkℓ) (middle) and InPd₃(hkℓ) (right) as a function of Δμₚd = μₚd - μᵇₚd where μₚd is the chemical potential of Pd in the compound and μᵇₚd is the cohesive energy of elemental Pd bulk. Pure Pd, pure In and mixed terminations are shown with solid, dashed and dash-dotted lines, respectively. Only the most stable surfaces are represented. Colour code for surface terminations: Pd = blue and In = brown.</figcaption></figure>

atomic density than the mixed termination of InPd₃(100) – 0.15 at/Å² vs. 0.12 at/Å². But the In content of the mixed InPd₃(100) termination is larger (0.06 In at./Å², vs. 0.04 In at./Å² for the stoichiometric surface), thus leading to a stabilising effect. In- deed, the presence of In atoms on the surface is less unfavourable than that of Pd. Similar surface energies are also calculated for InPd(100), InPd(110) and InPd(111) at the In-rich limit, and for InPd(111) and InPd(110) at the Pd-rich limit. The situa- tion is rather complex for In₃Pd₂, with relative stabilities calcu- lated to rank as γ<sub>(101)</sub> < γ<sub>(100)</sub> < γ<sub>(110)</sub> < γ<sub>(001)</sub> in the In-rich and γ<sub>(001)</sub> < γ<sub>(100)</sub> < γ<sub>(110)</sub> < γ<sub>(101)</sub> in the Pd-rich side.

Based on the surface energies at the In- and Pd-rich limits, and using the Wulff approach, two types of NPs are built (Fig. 3). Their morphologies strongly depend on the chemical potentials. The InPd₃ NP exhibits a cubic shape in the In-rich side and a trun- cated octahedral (TOh) morphology – typically the Wulff polyhe- dra of TM alloys – in the Pd-rich side. Because the stoichiomet- ric InPd(110) surface presents a low energy, independent of the chemical potentials, the InPd NP morphology is similar at both In- and Pd-rich limits, characterised by large (110) facets. The main difference between the two cases lies in small (100) facets ob- served at the In-rich side, due to the low InPd(100) surface energy at Δμ<sub>In</sub> = 0. The In₃Pd₂ NP displays a cylindrical shape at the Pd- rich limit, mainly due to the stability of the In₃Pd₂(001) surface for Δμ<sub>Pd</sub> = 0. In the In-rich side, the stability of the In₃Pd₂(110) surface leads to a shape with steeper facets. In the following,

Wulff Polyhedra, built using most stable facets at the In-rich and Pd-rich limits (Fig. 2, Tab. S1, section S1), are labelled WP<sub>In</sub> and WP<sub>Pd</sub>, respectively.

## 2.3 Stability of bimetallic In-Pd nanoparticles

The calculation of surface energies requires an infinite reservoir of particles, which becomes questionable when the size of the NPs is not large enough. Therefore, we build two types of NPs: on the one hand, NPs with a diameter between 22 Å and 28 Å (300 to 500 atoms, called “small” NPs in the following), and on the other hand, NPs with a diameter larger than 36 Å (made of around 1000 to 1200 atoms, called “large” NPs in the following).

Moreover, we also consider in this work other geometries, that can offer dense fcc{111} and fcc{100} facets when dealing with fcc metals, but whose relevance to the case of non fcc In-Pd NPs mostly relies on their spherical shape. We focused on icosahed- ron (Ih), decahedron (Dh), octahedron (Oh) and cuboctahedron (Cubo) shapes (Fig. 4). These morphologies have been chosen because, for metal NPs, (i) Non-crystalline Ih and Dh geometries are typically the ones of small to intermediate sized NPs,<sup>68,69</sup> while crystalline Oh and Cubo geometries are commonly in com- petition with the Wulff polyhedron for large size NPs. Below, not to be confused with the conventional Cubo and Oh shapes of L1<sub>2</sub> bimetallic NPs (as is InPd<sub>3</sub>), we rename Oh by Oh* and Cubo by Cubo*, to indicate that * correspond to In<sub>3</sub>Pd<sub>2</sub> or InPd. Both or- dered and chemically disordered NPs are generated. ll structures
<figure><img src="image_3.png" /><figcaption>Fig. 3 Morphologies of In₃Pd₂ (left), InPd (middle) and InPd₃ (right) NPs based on the Wulff construction. The diameter of each NP is given, along with its number of atoms. The top row corresponds to the NP morphology at the In-rich side whereas the bottom row corresponds to the Pd-rich side. Colour code: Pd = blue and In = brown.</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>Fig. 4 Non-crystalline (Ih and Dh) and crystalline (Oh and Cubo) morphologies of In₃Pd₂, InPd and InPd₃ NPs. The diameter of each NP is given, along with its number of atoms. Colour code: Pd = blue and In = brown.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>Fig. 5 Surface energies of In₃Pd₂ (left), InPd (middle) and InPd₃ (right) NPs as a function of Δμ<sub>Pd</sub> = μ<sub>Pd</sub> - μ<sub>Pd</sub><sup>b</sup>. Chemically ordered and disordered NPs are shown with full and dashed lines, respectively. Colour code: Pd = blue and In = brown.</figcaption></figure>

have been relaxed. Total energies of chemically disordered WPs result from an average over nine or ten configurations, generated by random shuffling atoms in the structure. While energy differences between the different configurations can be rather high – it reaches more than 60 meV/at for Wulff shaped In₃Pd₂ NPs, the dispersion of energies after structural relaxation is less than a few meV/at. (20 meV/at for Wulff shaped In₃Pd₂ NPs). The NP stabilities have been evaluated through their surface energies, calculated as a function of the Pd chemical potential. They are shown in Fig. 5. The three sub-figures correspond to NPs built from In₃Pd₂ (left), InPd (middle) and InPd₃ (right).

In the following, we first discuss the case of small NPs (Fig. 5). At the Pd-rich limit (Δμ<sub>Pd</sub> = 0), the most stable configuration

of InPd₃ NP is WP<sub>In</sub> which exhibits only low energy mixed (100) facets, i.e. faced composed of both In and Pd. At lower Δμ<sub>Pd</sub>, this structure is progressively destabilised and the WP<sub>Pd</sub> configuration, characterised by (100) and (111) facets, becomes the most stable one in the region -0.59 eV < Δμ<sub>Pd</sub> < -0.21 eV. Interestingly, at the crossover between WP<sub>In</sub> and WP<sub>Pd</sub> (i.e. -0.21 eV < Δμ<sub>Pd</sub> < -0.19 eV), the most stable configuration is the Oh structure with only mixed (111) facets. Increasing the indium concentration to InPd leads to a situation, where WP<sub>Pd</sub>, characterised by only mixed (110) facets, thus forming a quasi-spherical NP, is the most stable configuration on almost the whole range of Pd chemical potentials. At the Pd-rich limit, there is a competition between WP<sub>Pd</sub> and WP<sub>In</sub>, which mainly exhibits (100) facets.
<figure><img src="image_5.png" /><figcaption>Fig. 3 Morphologies of In₃Pd₂ (left), InPd (middle) and InPd₃ (right) NPs based on the Wulff construction. The diameter of each NP is given, along with its number of atoms. The top row corresponds to the NP morphology at the In-rich side whereas the bottom row corresponds to the Pd-rich side. Colour code: Pd = blue and In = brown.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>Fig. 4 Non-crystalline (Ih and Dh) and crystalline (Oh and Cubo) morphologies of In₃Pd₂, InPd and InPd₃ NPs. The diameter of each NP is given, along with its number of atoms. Colour code: Pd = blue and In = brown.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>Fig. 5 Surface energies of In₃Pd₂ (left), InPd (middle) and InPd₃ (right) NPs as a function of Δμ<sub>Pd</sub> = μ<sub>Pd</sub> - μ<sub>Pd</sub><sup>b</sup>. Chemically ordered and disordered NPs are shown with full and dashed lines, respectively. Colour code: Pd = blue and In = brown.</figcaption></figure>

have been relaxed. Total energies of chemically disordered WPs result from an average over nine or ten configurations, generated by random shuffling atoms in the structure. While energy differences between the different configurations can be rather high – it reaches more than 60 meV/at for Wulff shaped In₃Pd₂ NPs, the dispersion of energies after structural relaxation is less than a few meV/at. (20 meV/at for Wulff shaped In₃Pd₂ NPs). The NP stabilities have been evaluated through their surface energies, calculated as a function of the Pd chemical potential. They are shown in Fig. 5. The three sub-figures correspond to NPs built from In₃Pd₂ (left), InPd (middle) and InPd₃ (right).

In the following, we first discuss the case of small NPs (Fig. 5). At the Pd-rich limit (Δμ<sub>Pd</sub> = 0), the most stable configuration

of InPd₃ NP is WP<sub>In</sub> which exhibits only low energy mixed (100) facets, i.e. faced composed of both In and Pd. At lower Δμ<sub>Pd</sub>, this structure is progressively destabilised and the WP<sub>Pd</sub> configuration, characterised by (100) and (111) facets, becomes the most stable one in the region -0.59 eV < Δμ<sub>Pd</sub> < -0.21 eV. Interestingly, at the crossover between WP<sub>In</sub> and WP<sub>Pd</sub> (i.e. -0.21 eV < Δμ<sub>Pd</sub> < -0.19 eV), the most stable configuration is the Oh structure with only mixed (111) facets. Increasing the indium concentration to InPd leads to a situation, where WP<sub>Pd</sub>, characterised by only mixed (110) facets, thus forming a quasi-spherical NP, is the most stable configuration on almost the whole range of Pd chemical potentials. At the Pd-rich limit, there is a competition between WP<sub>Pd</sub> and WP<sub>In</sub>, which mainly exhibits (100) facets.
At the In-rich limit, Oh* is the most stable morphology. Focusing on In₃Pd₂ NPs, the most stable configurations are WPIn at the Pd-rich side, and Cubo* and Oh* at the In-rich side. The surface energies of Cubo* and Oh* are close to each other, Cubo* becoming more favourable than Oh* at the In-rich side. Overall non-crystalline and disordered structures are not favoured on the whole range of Δμ<sub>Pd</sub>, a consequence of the large formation enthalpy of the In-Pd compounds (|ΔH<sub>f</sub>|) which act as a driving force to stabilise chemical order and crystalline structures (see section 2.4). This is rather different from the observations on transition metal nanoalloys (|ΔH<sub>f</sub>| < 0.1eV/at. in many cases), in which many morphologies distinct from polyhedral shapes might occur.<sup>70–73</sup>

<figure><img src="image_6.png" /><figcaption>Fig. 6 Radial distribution function of In-In (brown) and Pd-Pd (blue) bonds as a function of the distance from the NP centre (r). The considered NPs are small WP<sub>Pd</sub> nanoparticles (diameter 24 Å) of In<sub>3</sub>Pd<sub>2</sub> (left), InPd (middle), InPd<sub>3</sub> (right) compounds. The upper and bottom rows correspond to ordered and disordered WP<sub>Pd</sub> NPs, respectively. Colour code: Pd = blue and In = brown.</figcaption></figure>

To deepen our analysis, we computed the Radial Distribution Function (RDF) along the radius of the NP as displayed in Fig. 6. In the ordered state, NP RDFs are characterised by alternating sharp peaks, for any considered In-Pd compound, i.e. In₃Pd₂, InPd and InPd₃. The computed first neighbour distances for In-Pd (2.78 Å) in InPd NPs are found in excellent agreement with experimental data: according to X-ray diffraction and X-ray absorption spectroscopy, values are 2.81 Å and 2.79 Å, respectively, for small InPd NPs (27 Å, label Pd-In 2.0 in Ref. ³⁵). In addition, our computed value for first neighbour Pd-Pd distances (3.22 Å) is also in good agreement with the same distance in bulk InPd (3.25 Å).

In the case of disordered InPd<sub>3</sub>, the RDF exhibits wider peaks than in the ordered state. For disordered compounds with larger In concentrations, i.e. for disordered InPd and In<sub>3</sub>Pd<sub>2</sub>, only one sharp peak is still present, that of the closest neighbours. It suggests an amorphisation of the nanoparticles. This is in agreement with recent observations of amorphous disordered NP structures in the pM-TM system.<sup>74</sup> This is however quite different from the observations on transition metal nanoalloys with crystalline dis-

ordered structures, which do not present any amorphisation.<sup>75</sup>

The size of NPs is known to have a strong impact on their stability. Indeed, increasing the size might change the NP terminations and thus their Pd contents, leading to a drastic modification of their surface energies. Size effects are investigated here by comparing small (300-500 atoms) and large NPs (≈ 1000 to 1200 atoms, Fig. 7). In the case of InPd<sub>3</sub>, in the Pd-rich side, the most stable configuration is WPIn (thick yellow line, Fig. 7), i.e. the same configuration as that identified for small sizes. As μ<sub>Pd</sub> decreases, this structure is progressively destabilised in favour of the Cubo morphology (thick green line, Fig. 7). The latter morphology is then the most stable one, in the range -0.59 eV < Δμ<sub>Pd</sub> < -0.21 eV. This is consistent with the result found for smaller NPs, for which the WP<sub>Pd</sub> structure has been found to be the most stable one (thin yellow file, Fig. 5). Indeed, both Cubo and WP<sub>Pd</sub> exhibit (100) and (111) facets, the core being the same in both cases. The only difference between Cubo and WP<sub>Pd</sub> NP is their facet sizes. In the case of InPd, the most stable morphology is WP<sub>Pd</sub> for both small and large NPs (thick and thin blue lines, Fig. 7), the surface energy decreasing by around 0.05 J/m² when the number of atoms in the NP is multiplied by three, approximately. In the case of In<sub>3</sub>Pd<sub>2</sub>, the Wulff Polyhedra is the most stable structure for large NPs, i.e. WP<sub>Pd</sub> at Δμ<sub>Pd</sub> = 0 (thick blue line, Fig. 7) and WPIn at Δμ<sub>Pd</sub> = -1.19 eV with a crossover around Δμ<sub>Pd</sub> = -0.6 eV. This is a direct signature of the size effect. The Oh* and Cubo* morphologies, which are typical for fcc metals, are not favoured in InPd<sub>3</sub> NPs, at least when Δμ<sub>Pd</sub> = 0, even if bulk InPd<sub>3</sub> presents some similarities with fcc metals.

## 2.4 In-In short range interactions drive possible amorphisations of NPs

We now quantitatively investigate the amorphisation induced by chemical disorder in small In-Pd NPs. Our discussion is based on the pseudo-pairwise electrostatic interaction ($\tilde{V}_{ij}(r)$) defined as:

$$ \tilde{V}_{ij}(r) = \frac{\Delta q_i \Delta q_j}{4\pi\epsilon_0} \langle \chi_i | \hat{\mathbf{r}}_{ij}^{-1} | \chi_j \rangle, \quad (1) $$

where Δqᵢ is the effective charge of atom i, ε₀ is the vacuum permittivity, $\hat{\mathbf{r}}_{ij}$ is the distance operator between atom i and atom j and $|\chi_i\rangle$ are ad-hoc localised electronic wavefunctions basis vectors for electrons of atom i. The effective electrostatic potential ($\tilde{V}_{ij}(r)$) for a given pair of atom (ij) separated by a distance r is positive and negative when the ij interaction is repulsive and attractive, respectively. To estimate $\langle \chi_i | \hat{\mathbf{r}}_{ij}^{-1} | \chi_j \rangle$, we have performed a Crystal Orbital Overlap Populations (COOP) analysis<sup>76</sup> using the Lobster package<sup>77</sup> which is coupled to the VASP software.<sup>78–80</sup> All bonds with lengths below 5.5 Å have been considered. Then, we have computed the following quantity:

$$ \langle \tilde{V}_{ij}(r) \rangle_{\text{b,} \delta r} = \frac{1}{\int_{0}^{\infty} \mathbf{1}_{|r_b - r| \le \delta r} dr} \int_{0}^{\infty} \tilde{V}_{ij}(r) \mathbf{1}_{|r_b - r| \le \delta r} dr \quad (2) $$

where $\mathbf{1}_{|r_b - r| \le \delta r}$ is the indicator function which is equal to 1 if $|r_b - r| \le \delta r$ and 0 otherwise ($\delta r$ is set to 0.55 Å). The radial $\langle \tilde{V}_{ij}(r) \rangle$ quantity defined in Eq. 2 is an average of the $\tilde{V}_{ij}(r)$ electrostatic bound energies for ij atom pairs in the bound range dis-
At the In-rich limit, Oh* is the most stable morphology. Focusing on In₃Pd₂ NPs, the most stable configurations are WPIn at the Pd-rich side, and Cubo* and Oh* at the In-rich side. The surface energies of Cubo* and Oh* are close to each other, Cubo* becoming more favourable than Oh* at the In-rich side. Overall non-crystalline and disordered structures are not favoured on the whole range of Δμ<sub>Pd</sub>, a consequence of the large formation enthalpy of the In-Pd compounds (|ΔH<sub>f</sub>|) which act as a driving force to stabilise chemical order and crystalline structures (see section 2.4). This is rather different from the observations on transition metal nanoalloys (|ΔH<sub>f</sub>| < 0.1eV/at. in many cases), in which many morphologies distinct from polyhedral shapes might occur.<sup>70–73</sup>

<figure><img src="image_8.png" /><figcaption>Fig. 6 Radial distribution function of In-In (brown) and Pd-Pd (blue) bonds as a function of the distance from the NP centre (r). The considered NPs are small WP<sub>Pd</sub> nanoparticles (diameter 24 Å) of In<sub>3</sub>Pd<sub>2</sub> (left), InPd (middle), InPd<sub>3</sub> (right) compounds. The upper and bottom rows correspond to ordered and disordered WP<sub>Pd</sub> NPs, respectively. Colour code: Pd = blue and In = brown.</figcaption></figure>

To deepen our analysis, we computed the Radial Distribution Function (RDF) along the radius of the NP as displayed in Fig. 6. In the ordered state, NP RDFs are characterised by alternating sharp peaks, for any considered In-Pd compound, i.e. In₃Pd₂, InPd and InPd₃. The computed first neighbour distances for In-Pd (2.78 Å) in InPd NPs are found in excellent agreement with experimental data: according to X-ray diffraction and X-ray absorption spectroscopy, values are 2.81 Å and 2.79 Å, respectively, for small InPd NPs (27 Å, label Pd-In 2.0 in Ref. ³⁵). In addition, our computed value for first neighbour Pd-Pd distances (3.22 Å) is also in good agreement with the same distance in bulk InPd (3.25 Å).

In the case of disordered InPd₃, the RDF exhibits wider peaks than in the ordered state. For disordered compounds with larger In concentrations, i.e. for disordered InPd and In₃Pd₂, only one sharp peak is still present, that of the closest neighbours. It suggests an amorphisation of the nanoparticles. This is in agreement with recent observations of amorphous disordered NP structures in the pM-TM system.<sup>74</sup> This is however quite different from the observations on transition metal nanoalloys with crystalline dis-

ordered structures, which do not present any amorphisation.<sup>75</sup>

The size of NPs is known to have a strong impact on their stability. Indeed, increasing the size might change the NP terminations and thus their Pd contents, leading to a drastic modification of their surface energies. Size effects are investigated here by comparing small (300-500 atoms) and large NPs (≈ 1000 to 1200 atoms, Fig. 7). In the case of InPd₃, in the Pd-rich side, the most stable configuration is WPIn (thick yellow line, Fig. 7), i.e. the same configuration as that identified for small sizes. As μPd decreases, this structure is progressively destabilised in favour of the Cubo morphology (thick green line, Fig. 7). The latter morphology is then the most stable one, in the range -0.59 eV < ΔμPd < -0.21 eV. This is consistent with the result found for smaller NPs, for which the WPd structure has been found to be the most stable one (thin yellow file, Fig. 5). Indeed, both Cubo and WPd exhibit (100) and (111) facets, the core being the same in both cases. The only difference between Cubo and WPd NP is their facet sizes. In the case of InPd, the most stable morphology is WPd for both small and large NPs (thick and thin blue lines, Fig. 7), the surface energy decreasing by around 0.05 J/m² when the number of atoms in the NP is multiplied by three, approximately. In the case of In₃Pd₂, the Wulff Polyhedra is the most stable structure for large NPs, i.e. WPd at ΔμPd = 0 (thick blue line, Fig. 7) and WPIn at ΔμPd = -1.19 eV with a crossover around ΔμPd = -0.6 eV. This is a direct signature of the size effect. The Oh* and Cubo* morphologies, which are typical for fcc metals, are not favoured in InPd₃ NPs, at least when ΔμPd = 0, even if bulk InPd₃ presents some similarities with fcc metals.

## 2.4 In-In short range interactions drive possible amorphisations of NPs

We now quantitatively investigate the amorphisation induced by chemical disorder in small In-Pd NPs. Our discussion is based on the pseudo-pairwise electrostatic interaction ($\tilde{V}_{ij}(r)$) defined as:

$$ \tilde{V}_{ij}(r) = \frac{\Delta q_i \Delta q_j}{4\pi\epsilon_0} \langle \chi_i | \hat{\mathbf{r}}_{ij}^{-1} | \chi_j \rangle, \quad (1) $$

where Δqᵢ is the effective charge of atom i, ε₀ is the vacuum permittivity, $\hat{\mathbf{r}}_{ij}$ is the distance operator between atom i and atom j and $|\chi_i\rangle$ are ad-hoc localised electronic wavefunctions basis vectors for electrons of atom i. The effective electrostatic potential ($\tilde{V}_{ij}(r)$) for a given pair of atom (ij) separated by a distance r is positive and negative when the ij interaction is repulsive and attractive, respectively. To estimate $\langle \chi_i | \hat{\mathbf{r}}_{ij}^{-1} | \chi_j \rangle$, we have performed a Crystal Orbital Overlap Populations (COOP) analysis<sup>76</sup> using the Lobster package<sup>77</sup> which is coupled to the VASP software.<sup>78–80</sup> All bonds with lengths below 5.5 Å have been considered. Then, we have computed the following quantity:

$$ \langle \tilde{V}_{ij} \rangle (r_b, \delta r) = \frac{1}{\int_0^\infty \mathbf{1}_{|r_b - r| \le \delta r} dr} \int_0^\infty \tilde{V}_{ij}(r) \mathbf{1}_{|r_b - r| \le \delta r} dr \quad (2) $$

where $\mathbf{1}_{|r_b - r| \le \delta r}$ is the indicator function which is equal to 1 if $|r_b - r| \le \delta r$ and 0 otherwise ($\delta r$ is set to 0.55 Å). The radial $\langle \tilde{V}_{ij} \rangle$ quantity defined in Eq. 2 is an average of the $\tilde{V}_{ij}(r)$ electrostatic bound energies for ij atom pairs in the bound range dis-
<figure><img src="image_7.png" /><figcaption>Fig. 7 Surface energies of In₃Pd₂ (left), InPd (middle) and InPd₃ (right) NPs as a function of Δμₚd Thin and thick lines represent small and large NPs, respectively. Colour code: Pd = blue and In = brown.</figcaption></figure>

$$
\text{tance } [r_b - \delta r, r_b + \delta r].
$$

<figure><img src="image_8.png" /><figcaption>Fig. 8 Radial average $\tilde{V}_{ij}$ energy ($(\tilde{V}_{ij})(r_b)$) of In-In bonds (brown), Pd-Pd bonds (blue) and In-Pd bonds (grey) for small NPs (diameter 24 Å) built by the Wulff construction using In₃Pd₂ (left), InPd (middle), InPd₃ (right) bulk compounds. The upper and bottom rows display ordered and disordered WPₚd, respectively. Dots are brute data and dashed lines result from interpolation. Dot sizes are proportional to bond numbers.</figcaption></figure>

In Fig. 8 and in Fig. S3, we systematically compare, for a given
composition, averaged $\tilde{V}_{ij}$ energies for chemically ordered and
chemically disordered nanoparticles. The upper row displays or-
dered WP<sub>Pd</sub> whereas the bottom row represents disordered WP<sub>Pd</sub>.
Dots are brute data, their sizes are proportional to the number of
bonds at a bond length equal to r<sub>b</sub> and dashed lines are fitted
Pauling functions.<sup>81</sup> For all compositions (In<sub>3</sub>Pd<sub>2</sub>, InPd, InPd<sub>3</sub>),
the magnitudes of electrostatic pair energies for In-In bonds are
much larger than the ones of In-Pd and Pd-Pd. Therefore, we fo-

cus in the following on In-In bonds to investigate the stability. In
InPd₃ WPₚd NPs, most In-In bonds are attractive, suggesting that
In-In interactions tend to stabilise NPs. In InPd and In₃Pd₂ WPₚd
NPs, the In-In electrostatic pair energy is a monotonically decreas-
ing function, with a high positive maximum close to r_c = 3.0 Å
(Fig.8). Since positive values indicate repulsive interactions, the
short In-In bonds in InPd and InₚPd₂ WPₚd NPs are destabilising.
To avoid such non favourable interactions, the structure is mod-
ified and as a consequence it induces amorphisation in the cor-
responding NP structures. It is also evidenced by Fig.S3, where
differences appear between panel A (non relaxed structures) and
panel B (relaxed structures for InPd and In₃Pd₂ WPₚd NPs (bot-
tom line, left and middle columns), with a minimum appearing
for r close to 3 Å.

Looking more precisely at the region, called $\mathcal{S}$, localised be-
tween 3 Å and 4 Å, relevant for shortest In-In bonds and then
meaningful for the stability, we found that $\mathcal{S}$ is an attractive
 resp. repulsive) domain for In-In bonds in InPd$_3$ NPs (resp. in
InPd and In$_3$Pd$_2$ NPs). Based on the analysis performed in section
2.3 (Fig. 6), we have recalculated the radial distribution func-
tions, focusing on the $\mathcal{S}$ domain (Fig. 9). We have plotted RDFs
before and after the structural optimisation, i.e. we considered
both non relaxed (Fig. 9A) and relaxed (Fig. 9B) NPs. In each
panel (A & B), the upper and bottom rows displays ordered and
disordered WPd NPs, respectively. In each sub-figure, the $\mathcal{S}$ do-
main is highlighted in colour. The blue filled area represents the
attractive In-In interactions in InPd$_3$ NPs, while the red filled ar-
eas represent the repulsive In-In interactions in In$_3$Pd$_2$ and InPd
NPs.

InPd₃ NPs exhibit sharp In-In peaks within the stability domain
<figure><img src="image_9.png" /><figcaption>Fig. 7 Surface energies of In₃Pd₂ (left), InPd (middle) and InPd₃ (right) NPs as a function of Δμₚd Thin and thick lines represent small and large NPs, respectively. Colour code: Pd = blue and In = brown.</figcaption></figure>

$$
\text{tance } [r_b - \delta r, r_b + \delta r].
$$

<figure><img src="image_10.png" /><figcaption>Fig. 8 Radial average $\tilde{V}_{ij}$ energy ($(\tilde{V}_{ij})(r_b)$) of In-In bonds (brown), Pd-Pd bonds (blue) and In-Pd bonds (grey) for small NPs (diameter 24 Å) built by the Wulff construction using In₃Pd₂ (left), InPd (middle), InPd₃ (right) bulk compounds. The upper and bottom rows display ordered and disordered WPₚd, respectively. Dots are brute data and dashed lines result from interpolation. Dot sizes are proportional to bond numbers.</figcaption></figure>

In Fig. 8 and in Fig. S3, we systematically compare, for a given
composition, averaged $\tilde{V}_{ij}$ energies for chemically ordered and
chemically disordered nanoparticles. The upper row displays or-
dered WP<sub>Pd</sub> whereas the bottom row represents disordered WP<sub>Pd</sub>.
Dots are brute data, their sizes are proportional to the number of
bonds at a bond length equal to r<sub>b</sub> and dashed lines are fitted
Pauling functions.<sup>81</sup> For all compositions (In<sub>3</sub>Pd<sub>2</sub>, InPd, InPd<sub>3</sub>),
the magnitudes of electrostatic pair energies for In-In bonds are
much larger than the ones of In-Pd and Pd-Pd. Therefore, we fo-

cus in the following on In-In bonds to investigate the stability. In
InPd₃ WPₚd NPs, most In-In bonds are attractive, suggesting that
In-In interactions tend to stabilise NPs. In InPd and In₃Pd₂ WPₚd
NPs, the In-In electrostatic pair energy is a monotonically decreas-
ing function, with a high positive maximum close to r_c = 3.0 Å
(Fig.8). Since positive values indicate repulsive interactions, the
short In-In bonds in InPd and InₚPd₂ WPₚd NPs are destabilising.
To avoid such non favourable interactions, the structure is mod-
ified and as a consequence it induces amorphisation in the cor-
responding NP structures. It is also evidenced by Fig.S3, where
differences appear between panel A (non relaxed structures) and
panel B (relaxed structures for InPd and In₃Pd₂ WPₚd NPs (bot-
tom line, left and middle columns), with a minimum appearing
for r close to 3 Å.

Looking more precisely at the region, called $\mathcal{S}$, localised be-
tween 3 Å and 4 Å, relevant for shortest In-In bonds and then
meaningful for the stability, we found that $\mathcal{S}$ is an attractive
 resp. repulsive) domain for In-In bonds in InPd₃ NPs (resp. in
InPd and In₃Pd₂ NPs). Based on the analysis performed in section
2.3 (Fig. 6), we have recalculated the radial distribution func-
tions, focusing on the $\mathcal{S}$ domain (Fig. 9). We have plotted RDFs
before and after the structural optimisation, i.e. we considered
both non relaxed (Fig. 9A) and relaxed (Fig. 9B) NPs. In each
panel (A & B), the upper and bottom rows displays ordered and
disordered WP<sub>Pd</sub> NPs, respectively. In each sub-figure, the $\mathcal{S}$ do-
main is highlighted in colour. The blue filled area represents the
attractive In-In interactions in InPd₃ NPs, while the red filled ar-
eas represent the repulsive In-In interactions in In₃Pd₂ and InPd
NPs.

InPd₃ NPs exhibit sharp In-In peaks within the stability domain
<figure><img src="image_9.png" /><figcaption>Fig. 9 Radial distribution function of In-In (brown), Pd-Pd (blue) and In-Pd (grey) bounds for small Pd-rich Wulff Polyhedra nanoparticles (diameter 24 Å) of In₃Pd₂ (left), InPd (middle), InPd₃ (right) compounds. The upper row displays ordered WPₚd whereas the bottom row represents disordered WPₚd. Radial function distributions for non relaxed and relaxed nanoparticles are drawn in figure A (top) and B (bottom), respectively. Blue filled area represents the attractive domain of In-In interactions in InPd₃ NPs. Red filled areas represent the repulsive domain of In-In interactions in In₃Pd₂ and InPd NPs. (S = [3.0 Å, 4.0 Å])</figcaption></figure>

S, at r ≃ 3.0 Å and r ≃ 4.0 Å for both chemically ordered and disordered structures. After relaxation, peaks remain sharp notably for chemically disordered NPs. Such NPs are crystalline (Fig. 6, right), the crystalline character resulting from the attractive In-In interactions Fig. 6, right, plotted in brown). Both InPd and In₃Pd₂ NPs show similar behaviours, distinct from the one observed for InPd₃, with sharp and large In-In RDF peaks in the S region for chemically ordered and disordered NPs, respectively. Upon relaxation, the RDF plots do not change much for ordered NPs, at most there is a slight broadening of peaks for In-Pd NPs, while the RDF features are strongly modified in the case of disordered NPs. Based on the previous pair interactions analysis (Fig. 8), short In-In bonds (< 3 Å) are highly repulsive, i.e. much more repulsive than in the S region. Thus, bond lengths should in-

crease to reduce energetically destabilising In-In interactions. In-deed, after relaxation, In-In RDF peaks disappeared at r < 3 Å, and shortest In-In peak here lies at r = 3.3 Å. Hence, the combination of chemical disorder and In-In bond length increase are demonstrated to drive the tendency towards amorphisation of chemically disordered InPd and In₃Pd₂ NPs. To summarise, In-In interactions destabilise the chemically disordered In₃Pd₂ and InPd NPs, leading to amorphisation. In contrast, the favourable In-In interactions observed for InPd₃ nanoparticles contribute to their stabilisation, and support their crystalline character.

## 2.5 NP work-functions, impact of the surface composition and structure

The work-function (WF), i.e. the energy barrier of an electron escaping from the surface, is a fundamental materials property, with many applications, especially those involving electron exchanges at the material surface, such as adsorption. Thus, it is directly linked to catalytic properties, as illustrated in the literature.⁸²,⁸³ Although it lacks direct structural and chemical information, it can be a powerful tool to design optimised material catalysts.⁵⁷,⁸⁴⁻⁹⁰

The first successful calculations of WFs for simple metal surfaces have been achieved on the basis of the jellium model and the slab method by Lang and Kohn:⁹¹

$$W_{\text{slab}} = \lim_{||\mathbf{r}|| \to +\infty} \langle V_{\text{LP}} \rangle(\mathbf{r}, \mathbf{n}) - \varepsilon_{\text{F}}, \quad (3)$$

where $\mathbf{r} \in \mathbb{R}^{3N}$ is the system coordinates, $\varepsilon_F$ is the Fermi level, $V_{LP}$ is the potential energy and $\langle V_{LP} \rangle(\mathbf{r}, \mathbf{n})$ is its averaged value in a plane defined by $\mathbf{n}$ (normal surface vector). The computation of $\langle V_{LP} \rangle(\mathbf{r}, \mathbf{n})$ is rather easy using slabs, but can be influenced by quantum-size effects, as discussed in the literature.⁹²⁻⁹⁴ Nevertheless, over the years, WFs have been computed for a wide range of materials, using this approach.⁹⁵⁻⁹⁹

In the case of NPs, the previous pseudo one-dimension scheme cannot be used. Thus, we propose to introduce a new methodology to compute NP WFs, based on the Lang and Kohn work:⁹¹

$$W_{\text{NP}}(\delta r) = \lim_{||\mathbf{r}|| \to +\infty} \langle V_{\text{LP}} \rangle(r, \delta r) - \varepsilon_{\text{F}}. \quad (4)$$

where $\langle V_{LP} \rangle(r, \delta r)$ is the radial averaged potential energy function (see equation S9, supporting information). We successfully validated our method by comparing computational results with experimental data obtained from gold nanoparticles (Fig. 10, Fig.S4 and supplementary S5 for more details about the procedure).¹⁰⁰ Sizes of calculated nanoparticles have been evaluated by examining the variations in $\langle V_{LP} \rangle$ as a function of the distance from their centre. Relative work functions have been considered, i.e. $\Delta W = W_{\text{NP}} - W_{\text{NP}}^{\infty}$, where $W_{\text{NP}}^{\infty}$ is the limit of the NP work function when its size tends towards infinity. This approach is adopted because the experimental measurements are conducted on supported nanoparticles, while the calculations are performed on isolated nanoparticles. The comparison of relative values serves to correct for substrate influence and potential measurement-related shifts, providing a means to align the experimental and calculated
<figure><img src="image_11.png" /><figcaption>Fig. 9 Radial distribution function of In-In (brown), Pd-Pd (blue) and In-Pd (grey) bounds for small Pd-rich Wulff Polyhedra nanoparticles (diameter 24 Å) of In₃Pd₂ (left), InPd (middle), InPd₃ (right) compounds. The upper row displays ordered WPₚd whereas the bottom row represents disordered WPₚd. Radial function distributions for non relaxed and relaxed nanoparticles are drawn in figure A (top) and B (bottom), respectively. Blue filled area represents the attractive domain of In-In interactions in InPd₃ NPs. Red filled areas represent the repulsive domain of In-In interactions in In₃Pd₂ and InPd NPs. (S = [3.0 Å, 4.0 Å])</figcaption></figure>

S, at r ≃ 3.0 Å and r ≃ 4.0 Å for both chemically ordered and disordered structures. After relaxation, peaks remain sharp notably for chemically disordered NPs. Such NPs are crystalline (Fig. 6, right), the crystalline character resulting from the attractive In-In interactions Fig. 6, right, plotted in brown). Both InPd and In₃Pd₂ NPs show similar behaviours, distinct from the one observed for InPd₃, with sharp and large In-In RDF peaks in the S region for chemically ordered and disordered NPs, respectively. Upon relaxation, the RDF plots do not change much for ordered NPs, at most there is a slight broadening of peaks for In-Pd NPs, while the RDF features are strongly modified in the case of disordered NPs. Based on the previous pair interactions analysis (Fig. 8), short In-In bonds (< 3 Å) are highly repulsive, i.e. much more repulsive than in the S region. Thus, bond lengths should in-

crease to reduce energetically destabilising In-In interactions. In-deed, after relaxation, In-In RDF peaks disappeared at r < 3 Å, and shortest In-In peak here lies at r = 3.3 Å. Hence, the combination of chemical disorder and In-In bond length increase are demonstrated to drive the tendency towards amorphisation of chemically disordered InPd and In₃Pd₂ NPs. To summarise, In-In interactions destabilise the chemically disordered In₃Pd₂ and InPd NPs, leading to amorphisation. In contrast, the favourable In-In interactions observed for InPd₃ nanoparticles contribute to their stabilisation, and support their crystalline character.

## 2.5 NP work-functions, impact of the surface composition and structure

The work-function (WF), i.e. the energy barrier of an electron escaping from the surface, is a fundamental materials property, with many applications, especially those involving electron exchanges at the material surface, such as adsorption. Thus, it is directly linked to catalytic properties, as illustrated in the literature.⁸²,⁸³ Although it lacks direct structural and chemical information, it can be a powerful tool to design optimised material catalysts.⁵⁷,⁸⁴⁻⁹⁰

The first successful calculations of WFs for simple metal surfaces have been achieved on the basis of the jellium model and the slab method by Lang and Kohn:⁹¹

$$W_{\text{slab}} = \lim_{||\mathbf{r}|| \to +\infty} \langle V_{\text{LP}} \rangle (\mathbf{r}, \mathbf{n}) - \varepsilon_{\text{F}}, \quad (3)$$

where $\mathbf{r} \in \mathbb{R}^{3N}$ is the system coordinates, $\varepsilon_F$ is the Fermi level, $V_{LP}$ is the potential energy and $\langle V_{LP} \rangle (\mathbf{r}, \mathbf{n})$ is its averaged value in a plane defined by $\mathbf{n}$ (normal surface vector). The computation of $\langle V_{LP} \rangle (\mathbf{r}, \mathbf{n})$ is rather easy using slabs, but can be influenced by quantum-size effects, as discussed in the literature.⁹²⁻⁹⁴ Nevertheless, over the years, WFs have been computed for a wide range of materials, using this approach.⁹⁵⁻⁹⁹

In the case of NPs, the previous pseudo one-dimension scheme cannot be used. Thus, we propose to introduce a new methodology to compute NP WFs, based on the Lang and Kohn work:⁹¹

$$W_{\text{NP}}(\delta r) = \lim_{||\mathbf{r}|| \to +\infty} \langle V_{\text{LP}} \rangle (\mathbf{r}, \delta r) - \varepsilon_{\text{F}}. \quad (4)$$

where $\langle V_{LP} \rangle (\mathbf{r}, \delta r)$ is the radial averaged potential energy function (see equation S9, supporting information). We successfully validated our method by comparing computational results with experimental data obtained from gold nanoparticles (Fig. 10, Fig.S4 and supplementary S5 for more details about the procedure).¹⁰⁰ Sizes of calculated nanoparticles have been evaluated by examining the variations in $\langle V_{LP} \rangle$ as a function of the distance from their centre. Relative work functions have been considered, i.e. $\Delta W = W_{\text{NP}} - W_{\text{NP}}^{\infty}$, where $W_{\text{NP}}^{\infty}$ is the limit of the NP work function when its size tends towards infinity. This approach is adopted because the experimental measurements are conducted on supported nanoparticles, while the calculations are performed on isolated nanoparticles. The comparison of relative values serves to correct for substrate influence and potential measurement-related shifts, providing a means to align the experimental and calculated
results. Although the size domains explored in the calculations (2 nm - 4 nm) and in the experimental approach (larger than 4 nm) slightly differ, the variations of the work functions with the NP sizes exhibit strong consistencies. Moreover, the trends nicely align with a Coulombic decay with the nanoparticle size, as expected by a simple electrostatic model.¹⁰¹ This strengthens the overall confidence in the reliability and applicability of our computational methodology.

<figure><img src="image_10.png" /><figcaption>Fig. 10 Direct comparison between computed and measured work-function for gold NPs. Experimental data (in blue) are extracted from Ref. ¹⁰⁰ and ab initio work-functions (in red) are calculated based on our methodology.</figcaption></figure>

The WFs computed for crystalline and amorphous In-Pd NPs are presented in Fig. 11 (values in Tab. S4). They are found of the same order of magnitude as infinite slabs, and rather different from the ones of pure In and pure Pd surfaces, thus highlighting an alloy effect, as frequently reported in the literature (for instance, see Refs ¹⁰²,¹⁰³). Our results suggest tiny size effects, since the work-functions of larger NPs are larger than those of smaller particles with similar surface compositions. According to the jellium model, the difference in work-functions is expected to be of the order of 0.07 eV.¹⁰¹,¹⁰⁴ But we observed a smaller difference (0.03 eV) between the WFs of the 369 and 1105 atom InPd WPd NPs with similar surface Pd contents (0.48 and 0.49 at.%, respectively). The chemical potential drastically impacts NP WFs: for instance, the WFs of the small In₃Pd₂ WPd and WPIn differ by roughly 0.15 eV, and reaches approximately 0.5 eV for InPd NPs. Moreover, the correlation of the WFs with the Pd surface content is clearly observed (Fig. 11), the correlation being roughly linear.

Finally, we quantitatively compare the WF difference between small ordered NPs described in Fig. 6 (upper row) and small chemically disordered NPs, with same size and composition, generated randomly. Work-functions of disordered NPs are gathered in Tab. 2, along with $\langle \Delta W_{NP} \rangle$, a parameter introduced to quantify the changes in WFs induced by disorder/amorphisation and defined by:

$$ \langle \Delta W_{NP} \rangle (N^{\text{dis}}) = \frac{1}{N^{\text{dis}}} \sum_{i=1}^{N^{\text{dis}}} (W_{\text{NP}}^{\text{ord}} - W_{\text{NP}}^{\text{dis}}(i)), \quad (5) $$

<figure><img src="image_11.png" /><figcaption>Fig. 11 Work-function (in eV) as a function of the surface Pd concentration (in at. %). Small and large symbols correspond to small and large nanoparticles, respectively. Colours give the NP morphologies while symbols give the compound ($\bigtriangledown$ In₃Pd₂, ○ InPd, △ InPd₃) from which it has been built. Data are gathered in Tab. S4.</figcaption></figure>

where $N^{\text{dis}}$ is the number of configurations considered to describe one disordered nanoparticle, while $W_{\text{NP}}^{\text{ord}}$ and $W_{\text{NP}}^{\text{dis}}(i)$ are the work-functions of ordered and disordered NPs (labelled with i), respectively. The standard error $\sigma(W_{\text{NP}}^{\text{dis}})/\sqrt{N^{\text{dis}}}$ of $\langle \Delta W_{NP} \rangle$ is calculated based on its standard deviation $\sigma(W_{\text{NP}}^{\text{dis}})$. In all cases, the computed standard error is small in comparison to $\langle \Delta W_{NP} \rangle$, which means that values of $\langle \Delta W_{NP} \rangle$ are statistically significant.

On the one hand, the change in WFs ($\langle \Delta W_{NP} \rangle$) for InPd₃ WPd NPs is tiny and slightly negative with a rather large – larger than in the other cases – standard error value. This is consistent with the previous structural analysis. Only small structural differences have been observed between ordered and disordered NPs in this case, the crystalline character of the InPd₃ WPd NP being kept when disorder is introduced. Thus, the WF modification ($\langle \Delta W_{NP} \rangle$) only comes from a surface Pd content which may slightly differ from the one in the ordered NP.

On the other hand, the change in WFs ($\langle \Delta W_{NP} \rangle$) for InPd and In₃Pd₂ WPd NPs is positive and rather large, of the order of 0.2 eV. This result is quite different from the one for InPd₃ WPd NPs, and may be induced by the structural changes investigated above, especially the tendency towards amorphisation observed for disordered WPd InPd and In₃Pd₂ NPs. Here, the change in WFs mainly comes from the modification of the NP surface structure.
results. Although the size domains explored in the calculations (2 nm - 4 nm) and in the experimental approach (larger than 4 nm) slightly differ, the variations of the work functions with the NP sizes exhibit strong consistencies. Moreover, the trends nicely align with a Coulombic decay with the nanoparticle size, as expected by a simple electrostatic model.¹⁰¹ This strengthens the overall confidence in the reliability and applicability of our computational methodology.

<figure><img src="image_12.png" /><figcaption>Fig. 10 Direct comparison between computed and measured work-function for gold NPs. Experimental data (in blue) are extracted from Ref. ¹⁰⁰ and ab initio work-functions (in red) are calculated based on our methodology.</figcaption></figure>

The WFs computed for crystalline and amorphous In-Pd NPs are presented in Fig. 11 (values in Tab. S4). They are found of the same order of magnitude as infinite slabs, and rather different from the ones of pure In and pure Pd surfaces, thus highlighting an alloy effect, as frequently reported in the literature (for instance, see Refs ¹⁰²,¹⁰³). Our results suggest tiny size effects, since the work-functions of larger NPs are larger than those of smaller particles with similar surface compositions. According to the jellium model, the difference in work-functions is expected to be of the order of 0.07 eV.¹⁰¹,¹⁰⁴ But we observed a smaller difference (0.03 eV) between the WFs of the 369 and 1105 atom InPd WPd NPs with similar surface Pd contents (0.48 and 0.49 at.%, respectively). The chemical potential drastically impacts NP WFs: for instance, the WFs of the small In₃Pd₂ WPd and WPIn differ by roughly 0.15 eV, and reaches approximately 0.5 eV for InPd NPs. Moreover, the correlation of the WFs with the Pd surface content is clearly observed (Fig. 11), the correlation being roughly linear.

Finally, we quantitatively compare the WF difference between small ordered NPs described in Fig. 6 (upper row) and small chemically disordered NPs, with same size and composition, generated randomly. Work-functions of disordered NPs are gathered in Tab. 2, along with $\langle \Delta W_{NP} \rangle$, a parameter introduced to quantify the changes in WFs induced by disorder/amorphisation and defined by:

$$ \langle \Delta W_{NP} \rangle (N^{\text{dis}}) = \frac{1}{N^{\text{dis}}} \sum_{i=1}^{N^{\text{dis}}} (W_{\text{NP}}^{\text{ord}} - W_{\text{NP}}^{\text{dis}}(i)), \quad (5) $$

<figure><img src="image_13.png" /><figcaption>Fig. 11 Work-function (in eV) as a function of the surface Pd concentration (in at. %). Small and large symbols correspond to small and large nanoparticles, respectively. Colours give the NP morphologies while symbols give the compound ($\bigtriangledown \text{In}_3\text{Pd}_2$, $\bigcirc \text{InPd}$, $\bigtriangleup \text{InPd}_3$) from which it has been built. Data are gathered in Tab. S4.</figcaption></figure>

where $N^{\text{dis}}$ is the number of configurations considered to describe one disordered nanoparticle, while $W_{\text{NP}}^{\text{ord}}$ and $W_{\text{NP}}^{\text{dis}}(i)$ are the work-functions of ordered and disordered NPs (labelled with $i$), respectively. The standard error $\sigma(W_{\text{NP}}^{\text{dis}})/\sqrt{N^{\text{dis}}}$ of $\langle \Delta W_{NP} \rangle$ is calculated based on its standard deviation $\sigma(W_{\text{NP}}^{\text{dis}})$. In all cases, the computed standard error is small in comparison to $\langle \Delta W_{NP} \rangle$, which means that values of $\langle \Delta W_{NP} \rangle$ are statistically significant.

On the one hand, the change in WFs ($\langle \Delta W_{NP} \rangle$) for InPd₃ WPd NPs is tiny and slightly negative with a rather large – larger than in the other cases – standard error value. This is consistent with the previous structural analysis. Only small structural differences have been observed between ordered and disordered NPs in this case, the crystalline character of the InPd₃ WPd NP being kept when disorder is introduced. Thus, the WF modification ($\langle \Delta W_{NP} \rangle$) only comes from a surface Pd content which may slightly differ from the one in the ordered NP.

On the other hand, the change in WFs ($\langle \Delta W_{NP} \rangle$) for InPd and In₃Pd₂ WPd NPs is positive and rather large, of the order of 0.2 eV. This result is quite different from the one for InPd₃ WPd NPs, and may be induced by the structural changes investigated above, especially the tendency towards amorphisation observed for disordered WPd InPd and In₃Pd₂ NPs. Here, the change in WFs mainly comes from the modification of the NP surface structure.
<table><caption>Table 2 WF differences between (chemically) ordered and disordered NPs, along with the standard error deviation.</caption><thead><tr><th>NP</th><th>N<sup>dis</sup></th><th>⟨ΔW<sub>NP</sub>⟩(N<sup>dis</sup>) (eV)</th><th>σ(W<sup>dis</sup><sub>NP</sub>)/√N<sup>dis</sup> (eV)</th></tr></thead><tbody><tr><td>In<sub>3</sub>Pd<sub>2</sub> WP<sub>Pd</sub></td><td>9</td><td>0.1842</td><td>0.005</td></tr><tr><td>InPd WP<sub>Pd</sub></td><td>9</td><td>0.2468</td><td>0.007</td></tr><tr><td>InPd<sub>3</sub> WP<sub>Pd</sub></td><td>10</td><td>-0.0431</td><td>0.010</td></tr></tbody></table>

### 3 Conclusions

This study contributes to a better understanding of structure- properties in NPs, through the investigation of In-Pd NPs built from bulk In<sub>3</sub>Pd<sub>2</sub>, InPd and InPd<sub>3</sub> intermetallics. According to surface energies, ordered NPs likely present crystalline forms, the Wulff Polyhedron being the most stable NP shape in the range 300-1200 atoms. A noticeable tendency towards amorphisation is observed when disorder is introduced in NPs built from In<sub>3</sub>Pd<sub>2</sub> and InPd, driven by non favourable In-In interactions. Tuning the NP surface structure and composition modifies the NP work- function, with variations of the order of 0.5 eV. Non crystalline morphologies, such as Ih or Dh, are not likely. This original be- haviour compared to that observed in metals and TM-TM alloys, likely stems from the rather elevated formation enthalpy of pM- TM intermetallics.

The current study is limited to three In-Pd phases with rather small crystal cells, and already gives important insights. More complex phases do exist in the phase diagram of In-Pd, thus chal- lenging the interplay between bulk structures and 3-dimensional NP morphologies built from them. Indeed, in addition to the selection of the termination plane discussed in the present pa- per, several questions arise when considering more complex com- pounds. Actually, the approach based on the Wulff construction should be adapted in the case of facets whose structure is only a fraction of the periodic surface. More generally, when the crystal cell is much larger than the NP size, one may expect that the dif- ferent part of the crystal are not equivalent towards nanoparticle formation, as, for instance, the clusters generally used to rep- resent the structure of such complex intermetallics.<sup>105</sup> Finally, because most heterogeneous reactions catalysed by nanoparticles occur at relatively high temperatures, molecular dynamic simu- lations could be useful to identify dynamical modulations of NP morphologies during reactions.<sup>106</sup>

### 4 Computational details

Electronic structure spin polarised calculations have been per- formed for static NPs with the plane wave Vienna ab initio simu- lation package (VASP),<sup>78–80</sup> using the projector augmented wave (PAW) method,<sup>107,108</sup> the generalised gradient approximation (GGA-PBE),<sup>109,110</sup> and an energy cutoff set to 450 eV. Ten va- lence electrons have been explicitly treated for Pd (5s<sup>1</sup>4d<sup>9</sup>) and three for In (5s<sup>2</sup>5p<sup>1</sup>).

Total energies have been minimised until the energy differences became less than 10<sup>-5</sup> eV between two electronic cycles. Atomic structures have been relaxed till the Hellmann-Feynman forces were as low as 0.01 eV/Å. Structures were plotted using the Ovito software<sup>111</sup>.

For each compositions, bulk-, slab- and nanoparticle-type cal-

culations have been performed. The different systems have been built with ASE<sup>112</sup> combined to MPInterfaces<sup>113</sup> and pymatgen<sup>114</sup>. The Brillouin zones of the different systems have been sampled with Γ-centred Monkhorst-Pack k-point meshes. In- finite surfaces have been modelled with symmetric slabs (Tab. S1), separated by a void thickness larger than 20 Å. Nanoparticles have been computed in large cubic simulation boxes (parameter equal to 50 Å). Only one k-point has been used for the calcula- tions with nanoparticles. The stabilities of the considered systems have been computed using formation enthalpies for bulk and NP ysstems (Tabs. 1 and S6) and surface energies as a function of the Pd chemical potential for slabs (see supporting information, sec- tion S1).<sup>115</sup> For NPs, historically, stabilities are evaluated through the Δ descriptor:

$$ \Delta = \frac{E_{\text{tot}}^{\text{NP}} - \sum_{i=1}^{M} N_i E_i^{\text{caw}}}{(\sum_{i=1}^{M} N_i)^{2/3}} \quad (6) $$

where $E_{\text{tot}}^{\text{NP}}$ is the total energy of the NP composed of M elements, $E_i^{\text{car}}$ is the bulk cohesive energy of the i species and $N_i$ is the num- ber of atom i in the nanoparticle. The Δ descriptor approximates the number of surface atoms by $(\sum_{i=1}^{M} N_i)^{2/3}$, which is not accu- rate for small NPs. In this work, stabilities of NPs are evaluated through surface energies ($\gamma_{NP}$):

$$ \gamma_{\text{NP}} = \frac{1}{A_{\text{NP}}} (E_{\text{tot}}^{\text{NP}} - n_{\text{Pd}}\mu_{\text{Pd}} - n_{\text{In}}\mu_{\text{In}}) \quad (7) $$

where $A_{NP}$ is the NP surface area, computed with Ovito<sup>111</sup> based on the Gaussian density method.<sup>116</sup> The quantities $\mu_{Pd}$ and $\mu_{In}$ (resp. $n_{Pd}$ and $n_{In}$) are the chemical potentials (resp. number of atoms) of Pd and In in the nanoparticle. The chemical potentials of Pd and In are not independent, since:

$$ (x + y)\Delta H_f(\text{In}_x\text{Pd}_y) = x(\mu_{\text{Pd}} - \mu_{\text{Pd}}^b) + y(\mu_{\text{In}} - \mu_{\text{In}}^b) \quad (8) $$

Assuming an equilibrium between the surface and the underlying NP core (assimilated to bulk), the chemical potentials are con- fined in a given range, i.e.

$$ \frac{(x+y)\Delta H_f}{y} \le \Delta\mu_{\text{Pd}} \le 0 \quad (9) $$

where $\Delta\mu_{Pd} = \mu_{Pd} - \mu_{Pd}^b$. In the following, we will discuss NP stabilities at the In-rich side, corresponding to $\Delta\mu_{Pd} = \frac{(x+y)\Delta H_f}{y}$ ($\Delta\mu_{In} = 0$), and at the Pd-rich side, corresponding to $\Delta\mu_{Pd} = 0$.

### Acknowledgements

This work was granted access to the High Performance Comput- ing (HPC) resources of IDRIS and TGCC under the allocation 99642 attributed by GENCI. HPC resources were also partially provided by the EXPLOR centre hosted by the Université de Lor- raine (project 2017M4XXX0108). One of us thanks Dr. David Loffreda for fruitful discussions. This work is supported by the European Integrated Centre for the Development of New Metal- lic Alloys and Compounds. Some of us acknowledge financial support through the COMETE project (CONception in silico de Matériaux pour l'Environnement et l'énergie) cofunded by the
<table><caption>Table 2 WF differences between (chemically) ordered and disordered NPs, along with the standard error deviation.</caption><thead><tr><th>NP</th><th>N<sup>dis</sup></th><th>⟨ΔW<sub>NP</sub>⟩(N<sup>dis</sup>) (eV)</th><th>σ(W<sup>dis</sup><sub>NP</sub>)/√N<sup>dis</sup> (eV)</th></tr></thead><tbody><tr><td>In<sub>3</sub>Pd<sub>2</sub> WP<sub>Pd</sub></td><td>9</td><td>0.1842</td><td>0.005</td></tr><tr><td>InPd WP<sub>Pd</sub></td><td>9</td><td>0.2468</td><td>0.007</td></tr><tr><td>InPd<sub>3</sub> WP<sub>Pd</sub></td><td>10</td><td>-0.0431</td><td>0.010</td></tr></tbody></table>

### 3 Conclusions

This study contributes to a better understanding of structure- properties in NPs, through the investigation of In-Pd NPs built from bulk In<sub>3</sub>Pd<sub>2</sub>, InPd and InPd<sub>3</sub> intermetallics. According to surface energies, ordered NPs likely present crystalline forms, the Wulff Polyhedron being the most stable NP shape in the range 300-1200 atoms. A noticeable tendency towards amorphisation is observed when disorder is introduced in NPs built from In<sub>3</sub>Pd<sub>2</sub> and InPd, driven by non favourable In-In interactions. Tuning the NP surface structure and composition modifies the NP work- function, with variations of the order of 0.5 eV. Non crystalline morphologies, such as Ih or Dh, are not likely. This original be- haviour compared to that observed in metals and TM-TM alloys, likely stems from the rather elevated formation enthalpy of pM- TM intermetallics.

The current study is limited to three In-Pd phases with rather small crystal cells, and already gives important insights. More complex phases do exist in the phase diagram of In-Pd, thus chal- lenging the interplay between bulk structures and 3-dimensional NP morphologies built from them. Indeed, in addition to the selection of the termination plane discussed in the present pa- per, several questions arise when considering more complex com- pounds. Actually, the approach based on the Wulff construction should be adapted in the case of facets whose structure is only a fraction of the periodic surface. More generally, when the crystal cell is much larger than the NP size, one may expect that the dif- ferent part of the crystal are not equivalent towards nanoparticle formation, as, for instance, the clusters generally used to rep- resent the structure of such complex intermetallics.<sup>105</sup> Finally, because most heterogeneous reactions catalysed by nanoparticles occur at relatively high temperatures, molecular dynamic simu- lations could be useful to identify dynamical modulations of NP morphologies during reactions.<sup>106</sup>

### 4 Computational details

Electronic structure spin polarised calculations have been per- formed for static NPs with the plane wave Vienna ab initio simu- lation package (VASP),<sup>78–80</sup> using the projector augmented wave (PAW) method,<sup>107,108</sup> the generalised gradient approximation (GGA-PBE),<sup>109,110</sup> and an energy cutoff set to 450 eV. Ten va- lence electrons have been explicitly treated for Pd (5s<sup>1</sup>4d<sup>9</sup>) and three for In (5s<sup>2</sup>5p<sup>1</sup>).

Total energies have been minimised until the energy differences became less than 10<sup>-5</sup> eV between two electronic cycles. Atomic structures have been relaxed till the Hellmann-Feynman forces were as low as 0.01 eV/Å. Structures were plotted using the Ovito software<sup>111</sup>.

For each compositions, bulk-, slab- and nanoparticle-type cal-

culations have been performed. The different systems have been built with ASE<sup>112</sup> combined to MPInterfaces<sup>113</sup> and pymatgen<sup>114</sup>. The Brillouin zones of the different systems have been sampled with Γ-centred Monkhorst-Pack k-point meshes. In- finite surfaces have been modelled with symmetric slabs (Tab. S1), separated by a void thickness larger than 20 Å. Nanoparticles have been computed in large cubic simulation boxes (parameter equal to 50 Å). Only one k-point has been used for the calcula- tions with nanoparticles. The stabilities of the considered systems have been computed using formation enthalpies for bulk and NP ysstems (Tabs. 1 and S6) and surface energies as a function of the Pd chemical potential for slabs (see supporting information, sec- tion S1).<sup>115</sup> For NPs, historically, stabilities are evaluated through the Δ descriptor:

$$ \Delta = \frac{E_{\text{tot}}^{\text{NP}} - \sum_{i=1}^{M} N_i E_i^{\text{caw}}}{(\sum_{i=1}^{M} N_i)^{2/3}} \quad (6) $$

where $E_{\text{tot}}^{\text{NP}}$ is the total energy of the NP composed of M elements, $E_i^{\text{car}}$ is the bulk cohesive energy of the i species and $N_i$ is the num- ber of atom i in the nanoparticle. The Δ descriptor approximates the number of surface atoms by $(\sum_{i=1}^{M} N_i)^{2/3}$, which is not accu- rate for small NPs. In this work, stabilities of NPs are evaluated through surface energies ($\gamma_{\text{NP}}$):

$$ \gamma_{\text{NP}} = \frac{1}{A_{\text{NP}}} (E_{\text{tot}}^{\text{NP}} - n_{\text{Pd}}\mu_{\text{Pd}} - n_{\text{In}}\mu_{\text{In}}) \quad (7) $$

where $A_{\text{NP}}$ is the NP surface area, computed with Ovito<sup>111</sup> based on the Gaussian density method.<sup>116</sup> The quantities $\mu_{\text{Pd}}$ and $\mu_{\text{In}}$ (resp. $n_{\text{Pd}}$ and $n_{\text{In}}$) are the chemical potentials (resp. number of atoms) of Pd and In in the nanoparticle. The chemical potentials of Pd and In are not independent, since:

$$ (x + y)\Delta H_f(\text{In}_x\text{Pd}_y) = x(\mu_{\text{Pd}} - \mu_{\text{Pd}}^b) + y(\mu_{\text{In}} - \mu_{\text{In}}^b) \quad (8) $$

Assuming an equilibrium between the surface and the underlying NP core (assimilated to bulk), the chemical potentials are con- fined in a given range, i.e.

$$ \frac{(x+y)\Delta H_f}{y} \le \Delta\mu_{\text{Pd}} \le 0 \quad (9) $$

where $\Delta\mu_{\text{Pd}} = \mu_{\text{Pd}} - \mu_{\text{Pd}}^b$. In the following, we will discuss NP stabilities at the In-rich side, corresponding to $\Delta\mu_{\text{Pd}} = \frac{(x+y)\Delta H_f}{y}$ ($\Delta\mu_{\text{In}} = 0$), and at the Pd-rich side, corresponding to $\Delta\mu_{\text{Pd}} = 0$.

### Acknowledgements

This work was granted access to the High Performance Comput- ing (HPC) resources of IDRIS and TGCC under the allocation 99642 attributed by GENCI. HPC resources were also partially provided by the EXPLOR centre hosted by the Université de Lor- raine (project 2017M4XXX0108). One of us thanks Dr. David Loffreda for fruitful discussions. This work is supported by the European Integrated Centre for the Development of New Metal- lic Alloys and Compounds. Some of us acknowledge financial support through the COMETE project (CONception in silico de Matériaux pour l'Environnement et l'énergie) cofunded by the
European Union under the program FEDER-FSE Lorraine et Mas-
sif des Vosges 2014-2020. One of us thanks financial support from
"Lorraine Université d'Excellence".

Notes and references

1. I. Freestone, N. Meeks, M. Sax and C. Higgitt, *Gold Bulletin*, **2007**, *40*, 270.

2. J. Jeevanandam, A. Barhoum, Y. Chan, A. Dufresne and M. Danquah, *Beilstein J Nanotechnol.*, 2018, **9**, 1050.

3 *Fundamentals and Applications of Nano Silicon in Plasmonics and Fullerines*, ed. M. Nayfeh, Elsevier, 2018, pp. 497–518.

4 E. C. Dreaden, A. M. Alkilany, X. Huang, C. J. Murphy and M. A. El-Sayed, *Chem. Soc. Rev.*, 2012, **41**, 2740.

5 C. Shasha and K. Krishnan, *Adv. Mater.*, 2021, **33**, 1904131.

6 S. Campora and G. Ghersi, *Nanotechnology Reviews*, 2022, **11**, 2595.

7 H. You, S. Yang, B. Ding and H. Yang, *Chem. Soc. Rev.*, 2013, **42**, 2880.

8 D. J. de Aberasturi, A. B. Serrano-Montes and L. M. Liz-Marzan, *Advanced Optical Materials*, 2015, **3**, 602.

9 A. Kumar, P. Choudhary, A. Kumar, P. H. C. Camargo and V. Krishn, *small*, 2022, **18**, 2101638.

10 J. A. Darr, J. Zhang, N. M. Makwana and X. Weng, *Chemical Reviews*, 2017, **117**, 11125.

11 L. Wang, M. H. Kafshgari and M. Meunier, *Advanced Functional Materials*, 2020, **30**, 2005400.

12 J. Grunes, J. Zhu and G. Somorjai, *Chem. Commun.*, 2003, **18**, 2257.

13 T. Bligaard and J. Norskov, in *Chemical Bonding at Surfaces and Interfaces*, ed. A. Nilsson, L. G. Pettersson and J. K. Norskov, Elsevier, 2008, ch. Chapter 4 - Heterogeneous Catalysis, p. 255.

14 J. R. Ludwig and C. S. Schindler, *Chem*, 2017, **2**, 313.

15 D. Wang and D. Astruc, *Chem. Soc. Rev.*, 2017, **46**, 816.

16 M. H. ns T. Kobayashi, H. Sano and N. Yamada, *Chem. Lett.*, 1987, **16**, 405.

17 S. Vukojevic, O. Trapp, J.-D. Grunwaldt, C. Kiener and F. Schüth, *Angew. Chem.*, 2005, **117**, 8192.

18 X. Liu, J. Iocozzia, Y. Wang, X. Cui, Y. Chen, S. Zhao and Z. Li, *Energy Environ. Sci.*, 2017, **10**, 402.

19 Z. W. Seh, J. Kibsgaard, C. F. Dickens, I. Chorkendorff, J. K. Norskov and T. F. Jaramillo, *Science*, 2017, **355**, eaad4998.

20 L. Liu and A. Corma, *Chemical Reviews*, 2018, **118**, 4981.

21 A. Medford, A. Vojvodic, J. S. Hummelshoj, J. Voss, F. Abild-Pedersen, F. Studt, T. Bligaard, A. Nilsson and J. K. Norskov, *J. Catal.*, 2015, **328**, 36.

22 J. Norskov, F. Studt, F. Abild-Pedersen and T. Bligaard, *Fundamental Concepts in Heterogeneous Catalysis*, John Wiley & Sons Inc., 2014.

23 A. Vojvodic and J. K. Nørskov, *National Science Review*, 2015, **2**, 140.

24 M. Armbruester, *Sci Technol Adv Mater.*, 2020, **21**, 303.

25 Y. Nakaya and S. Furukawa, *Chem. Rev.*, 2023, **123**, 5859.

26 B. G. Pollet, S. S. Kocha and I. Staffell, *Current Opinion in*

*Electrochemistry*, 2019, **16**, 90.

27 V. Stamenkovic, T. Schmidt, P. Ross and N. Markovic, *Journal of Electroanalytical Chemistry*, 2003, **554**, 191.

28 V. R. Stamenkovic, B. Fowler, B. S. Mun, G. Wang, P. N. Ross, C. A. Lucas and N. M. Markovic, *Science*, 2007, **315**, 493.

29 R. Ferrando, J. Jellinek and R. L. Johnston, *Chem. Rev.*, 2008, **108**, 845.

30 A. Jain, S. P. O. G. Hautier, W. Chen, W. D. Richards, W. Dacek, S. Cholia, D. Gunter, D. Skinner, G. Ceder and K. A. Persson, *APL Mater.*, 2013, **1**, 011002.

31 J. E. Montoya-Cardona, A. Salichon, N. Tarrat, E. Gaudry and D. Loffreda, *J. Phys. Chem. C*, 2023, **127**, 18043.

32 N. Köwitsch, L. Thoni, B. Klemmed, A. Benad, P. Paciok, M. Heggen, I. Köwitsch, M. Mehring, A. Eychmüller and M. Armbrüster, *ACS Catalysis*, 2021, **11**, 304–312.

33 A. Garcia-Trenco, A. Regoutz, E. R. White, D. J. Payne, M. S. Shaffer and C. K. Williams, *Applied Catalysis B: Environmental*, 2018, **220**, 9.

34 P. Wu and B. Yang, *Catal. Sci. Technol.*, 2019, **9**, 6102.

35 Z. Wu, E. Wegener, H.-T. Tseng, J. R. Gallagher, J. W. Harris, R. E. Diaz, Y. Ren, F. H. Ribeiro and J. T. Miller, *Catal. Sci. Technol.*, 2016, **6**, 6965.

36 Y. Luo, S. Alarcón Villaseca, M. Friedrich, D. Teschner, A. Knop-Gericke and M. Armbrüster, *Journal of Catalysis*, 2016, **338**, 265–272.

37 Y. Cao, Z. Sui, Y. Zhu, X. Zhou and D. Chen, *ACS Catalysis*, 2017, **7**, 7835–7846.

38 P. Markov, A. Bukhtiyarov, I. Mashkovsky, N. S. Smirnova, I. P. Prosvirin, Z. S. Vinokurov, M. A. Panafidin, G. N. Baeva, Y. V. Zubavichus, V. I. Bukhtiyarov and A. Y. Stakheev, *Kinet Catal*, 2019, **60**, 842.

39 S. Chen, X. Huang, D. Schild, D. Wang, C. Kübel and S. Behrens, *Nanoscale*, 2022, **14**, 17661.

40 D. Pavesi, F. S. M. Ali, D. Anastasiadou, T. K. M. Figueiredo, G.-J. M. Gruter, M. T. M. Koper and K. J. P. Schouten, *Catal. Sci. Technol.*, 2020, **10**, 4264.

41 Q. Feng and S. Zhao and Y. Wang and J. Dong and W. Chen and D. He and D. Wang and J. Yang and Y. Zhu and H. Zhu and L. Gu and Z. Li and Y. Liu and R. Yu and J. Li and Y. Li, *J. Am. Chem. Soc.*, 2017, **139**, 7294–7301.

42 C. Wolverton and V. Ozolins, *Phys. Rev. B*, 2006, **73**, 144104.

43 R. Ferro and A. Saccone, *Intermetallic Chemistry*, Pergamon
Elsevier, Oxford, UK, 2008.

44 C. White and H. Okamoto, Phase Diagrams of Indium Alloys and their Engineering Applications, 1992.

45 M. Neumann, D. Teschner, A. Knop-Gericke, W. Reschetilowski and M. Armbruester, *J. Catal.*, 2016, **340**, 49.

46 J. Pirart, A. Front, D. Rapetti, C. Andreazza-Vignolle, P. Andreazza, C. Mottet and R. Ferrando, *Nature Communications*, 2019, **10**, 1982.

47 F. Tournus, K. Sato, T. Epicier, T. J. Konno and V. Dupuis, *Phys. Rev. Lett.*, 2013, **110**, 055501.

48 J. Penuelas, P. Andreazza, C. Andreazza-Vignolle, H. C. N.
European Union under the program FEDER-FSE Lorraine et Mas-
sif des Vosges 2014-2020. One of us thanks financial support from
"Lorraine Université d'Excellence".

Notes and references

1. I. Freestone, N. Meeks, M. Sax and C. Higgitt, *Gold Bulletin*, **2007**, *40*, 270.

2 J. Jeevanandam, A. Barhoum, Y. Chan, A. Dufresne and M. Danquah, *Beilstein J Nanotechnol.*, 2018, **9**, 1050.

3 *Fundamentals and Applications of Nano Silicon in Plasmonics and Fullerines*, ed. M. Nayfeh, Elsevier, 2018, pp. 497–518.

4 E. C. Dreaden, A. M. Alkilany, X. Huang, C. J. Murphy and M. A. El-Sayed, *Chem. Soc. Rev.*, 2012, **41**, 2740.

5 C. Shasha and K. Krishnan, *Adv. Mater.*, 2021, **33**, 1904131.

6 S. Campora and G. Ghersi, *Nanotechnology Reviews*, 2022, **11**, 2595.

7 H. You, S. Yang, B. Ding and H. Yang, *Chem. Soc. Rev.*, 2013, **42**, 2880.

8 D. J. de Aberasturi, A. B. Serrano-Montes and L. M. Liz-Marzan, *Advanced Optical Materials*, 2015, **3**, 602.

9 A. Kumar, P. Choudhary, A. Kumar, P. H. C. Camargo and V. Krishn, *small*, 2022, **18**, 2101638.

10 J. A. Darr, J. Zhang, N. M. Makwana and X. Weng, *Chemical Reviews*, 2017, **117**, 11125.

11 L. Wang, M. H. Kafshgari and M. Meunier, *Advanced Functional Materials*, 2020, **30**, 2005400.

12 J. Grunes, J. Zhu and G. Somorjai, *Chem. Commun.*, 2003, **18**, 2257.

13 T. Bligaard and J. Norskov, in *Chemical Bonding at Surfaces and Interfaces*, ed. A. Nilsson, L. G. Pettersson and J. K. Norskov, Elsevier, 2008, ch. Chapter 4 - Heterogeneous Catalysis, p. 255.

14 J. R. Ludwig and C. S. Schindler, *Chem*, 2017, **2**, 313.

15 D. Wang and D. Astruc, *Chem. Soc. Rev.*, 2017, **46**, 816.

16 M. H. ns T. Kobayashi, H. Sano and N. Yamada, *Chem. Lett.*, 1987, **16**, 405.

17 S. Vukojevic, O. Trapp, J.-D. Grunwaldt, C. Kiener and F. Schüth, *Angew. Chem.*, 2005, **117**, 8192.

18 X. Liu, J. Iocozzia, Y. Wang, X. Cui, Y. Chen, S. Zhao and Z. Li, *Energy Environ. Sci.*, 2017, **10**, 402.

19 Z. W. Seh, J. Kibsgaard, C. F. Dickens, I. Chorkendorff, J. K. Norskov and T. F. Jaramillo, *Science*, 2017, **355**, eaad4998.

20 L. Liu and A. Corma, *Chemical Reviews*, 2018, **118**, 4981.

21 A. Medford, A. Vojvodic, J. S. Hummelshoj, J. Voss, F. Abild-Pedersen, F. Studt, T. Bligaard, A. Nilsson and J. K. Norskov, *J. Catal.*, 2015, **328**, 36.

22 J. Norskov, F. Studt, F. Abild-Pedersen and T. Bligaard, *Fundamental Concepts in Heterogeneous Catalysis*, John Wiley & Sons Inc., 2014.

23 A. Vojvodic and J. K. Nørskov, *National Science Review*, 2015, **2**, 140.

24 M. Armbruester, *Sci Technol Adv Mater.*, 2020, **21**, 303.

25 Y. Nakaya and S. Furukawa, *Chem. Rev.*, 2023, **123**, 5859.

26 B. G. Pollet, S. S. Kocha and I. Staffell, *Current Opinion in*

*Electrochemistry*, 2019, **16**, 90.

27 V. Stamenkovic, T. Schmidt, P. Ross and N. Markovic, *Journal of Electroanalytical Chemistry*, 2003, **554**, 191.

28 V. R. Stamenkovic, B. Fowler, B. S. Mun, G. Wang, P. N. Ross, C. A. Lucas and N. M. Markovic, *Science*, 2007, **315**, 493.

29 R. Ferrando, J. Jellinek and R. L. Johnston, *Chem. Rev.*, 2008, **108**, 845.

30 A. Jain, S. P. O. G. Hautier, W. Chen, W. D. Richards, W. Dacek, S. Cholia, D. Gunter, D. Skinner, G. Ceder and K. A. Persson, *APL Mater.*, 2013, **1**, 011002.

31 J. E. Montoya-Cardona, A. Salichon, N. Tarrat, E. Gaudry and D. Loffreda, *J. Phys. Chem. C*, 2023, **127**, 18043.

32 N. Köwitsch, L. Thoni, B. Klemmed, A. Benad, P. Paciok, M. Heggen, I. Köwitsch, M. Mehring, A. Eychmüller and M. Armbrüster, *ACS Catalysis*, 2021, **11**, 304–312.

33 A. Garcia-Trenco, A. Regoutz, E. R. White, D. J. Payne, M. S. Shaffer and C. K. Williams, *Applied Catalysis B: Environmental*, 2018, **220**, 9.

34 P. Wu and B. Yang, *Catal. Sci. Technol.*, 2019, **9**, 6102.

35 Z. Wu, E. Wegener, H.-T. Tseng, J. R. Gallagher, J. W. Harris, R. E. Diaz, Y. Ren, F. H. Ribeiro and J. T. Miller, *Catal. Sci. Technol.*, 2016, **6**, 6965.

36 Y. Luo, S. Alarcón Villaseca, M. Friedrich, D. Teschner, A. Knop-Gericke and M. Armbrüster, *Journal of Catalysis*, 2016, **338**, 265–272.

37 Y. Cao, Z. Sui, Y. Zhu, X. Zhou and D. Chen, *ACS Catalysis*, 2017, **7**, 7835–7846.

38 P. Markov, A. Bukhtiyarov, I. Mashkovsky, N. S. Smirnova, I. P. Prosvirin, Z. S. Vinokurov, M. A. Panafidin, G. N. Baeva, Y. V. Zubavichus, V. I. Bukhtiyarov and A. Y. Stakheev, *Kinet Catal*, 2019, **60**, 842.

39 S. Chen, X. Huang, D. Schild, D. Wang, C. Kübel and S. Behrens, *Nanoscale*, 2022, **14**, 17661.

40 D. Pavesi, F. S. M. Ali, D. Anastasiadou, T. K. M. Figueiredo, G.-J. M. Gruter, M. T. M. Koper and K. J. P. Schouten, *Catal. Sci. Technol.*, 2020, **10**, 4264.

41 Q. Feng and S. Zhao and Y. Wang and J. Dong and W. Chen and D. He and D. Wang and J. Yang and Y. Zhu and H. Zhu and L. Gu and Z. Li and Y. Liu and R. Yu and J. Li and Y. Li, *J. Am. Chem. Soc.*, 2017, **139**, 7294–7301.

42 C. Wolverton and V. Ozolins, *Phys. Rev. B*, 2006, **73**, 144104.

43 R. Ferro and A. Saccone, *Intermetallic Chemistry*, Pergamon
Elsevier, Oxford, UK, 2008.

44 C. White and H. Okamoto, Phase Diagrams of Indium Alloys and their Engineering Applications, 1992.

45 M. Neumann, D. Teschner, A. Knop-Gericke, W. Reschetilowski and M. Armbruester, *J. Catal.*, 2016, **340**, 49.

46 J. Pirart, A. Front, D. Rapetti, C. Andreazza-Vignolle, P. Andreazza, C. Mottet and R. Ferrando, *Nature Communications*, 2019, **10**, 1982.

47 F. Tournus, K. Sato, T. Epicier, T. J. Konno and V. Dupuis, *Phys. Rev. Lett.*, 2013, **110**, 055501.

48 J. Penuelas, P. Andreazza, C. Andreazza-Vignolle, H. C. N.
Tolentino, M. De Santis and C. Mottet, *Phys. Rev. Lett.*, 2008, **100**, 115502.

49 M. Zhou, C. Li and J. Fang, *Chemical Reviews*, 2021, **121**, 736.

50 P. Makkar and N. N. Ghosh, *RSC Adv.*, 2021, **11**, 27897.

51 H. Demir and L. C. Grabow, *ACS Applied Nano Materials*, 2020, **3**, 6127.

52 D. Vollath, F. Fischer and D. Holec, *Beilstein J Nanotechnol.*, 2018, **9**, 2265.

53 H. Amara, J. Nelayah, J. Creuze, A. Chmielewski, D. Alloyeau, C. Ricolleau and B. Legrand, *Phys. Rev. B*, 2022, **105**, 165403.

54 A. Chmielewski, J. Nelayah, H. Amara, J. Creuze, D. Alloyeau, G. Wang and C. Ricolleau, *Phys. Rev. Lett.*, 2018, **120**, 025901.

55 L. Lin, R. Jacobs, T. Ma, D. Chen, J. Booske and D. Morgan, *Phys. Rev. Applied*, 2023, **19**, 037001.

56 F. Calle-Vallejo, M. T. M. Koper and A. S. Bandarenka, *Chem. Soc. Rev.*, 2013, **42**, 5210.

57 S. Trasatti, *J. Electroanal. Chem.*, 1972, **39**, 163.

58 H. Kohlmann, *Journal of Solid State Chemistry*, 2010, **183**, 367.

59 I. Harris, M. Norman and A. W. Bryant, *Journal of the Less-Common Metals*, 1968, **16**, 427.

60 C. Garnero, M. Lepesant, C. Garcia-Marcelot, Y. Shin, C. Meny, P. Farger, B. Warot-Fonrose, R. Arenal, G. Viau, K. Soulantica, P. Fau, P. Poveda, L.-M. Lacroix and B. Chaudret, *Nano Letters*, 2019, **19**, 1379–1386.

61 A. Front and C. Mottet, *J. Phys. Chem. C*, 2021, **125**, 16358–16365.

62 Y. Yang, C.-C. Chen, M. C. Scott, C. Ophus, R. Xu, A. Pryor, L. Wu, F. Sun, W. Theis, J. Zhou, M. Eisenbach, P. R. C. Kent, R. F. Sabirianov, H. Zeng, P. Ercius and J. Miao, *Nature*, 2017, **542**, 75–79.

63 S. Amore, S. Delsante, N. Parodi and G. Borzone, *Thermochimica Acta*, 2009, **481**, 1.

64 A. Ciccioli, G. Balducci, G. Gigli, L. Perring and F. Bussy, *Intermetallics*, 2000, **8**, 195.

65 A. Ormeci, E. Gaudry, M. Armbrüster and Y. Grin, *Chemistry-Open*, 2022, **11**.

66 N. Roy, A. Chakrabarty, B. Koley, T. Saha-Dasgupta and P. P. Jana, *Journal of Solid State Chemistry*, 2020, **290**, 121567.

67 G. Wulff, *Zeitschrift für Kristallographie - Crystalline Materials*, 1901, **34**, 449–530.

68 H. Wang, S. Zhou, K. D. Gilroy, Z. Cai and Y. Xia, *Nano Today*, 2017, **15**, 121–144.

69 G. Casillas, J. J. Velazquez-Salazar and M. Jose-Yacaman, *J. Phys. Chem. C*, 2012, **116**, 8844.

70 A. Front and C. Mottet, *J. Phys. Chem. C*, 2022, 2.

71 A. Front and C. Mottet, *Theoretical Chemistry Accounts*, 2022, **141**, 2.

72 K. Loza, M. Heggen and M. Epple, *Adv. Funct. Mater.*, 2020, **30**, 1909260.

73 B. Farkas and N. H. de Leeuw, *Nanotechnology*, 2020, **31**,

195711.

74 H. Yang, K. Wang, Z. Tang, Z. Liu and S. Chen, *Journal of Catalysis*, 2020, **382**, 181–191.

75 D. Alloyeau, C. Ricolleau, C. Mottet, T. Oikawa, C. Langlois, Y. Le Bouar, N. Braidy and A. Loiseau, *Nature Materials*, 2009, **8**, 940.

76 R. Dronskowski and P. E. Bloechl, *The Journal of Physical Chemistry*, 1993, **97**, 8617–8624.

77 S. Maintz, V. L. Deringer, A. L. Tchougréeff and R. Dronskowski, *Journal of Computational Chemistry*, 2016, **37**, 1030–1035.

78 G. Kresse and J. Hafner, *Physical Review B*, 1993, **47**, 558–561.

79 G. Kresse and J. Furthmüller, *Computational Materials Science*, 1996, **6**, 15–50.

80 G. Kresse and J. Furthmüller, *Physical Review B*, 1996, **54**, 11169–11186.

81 L. Pauling, *The Nature of the Chemical Bond—An Introduction to Modern Structural Chemistry*, Cornell University Press, New-York, 1960.

82 T. Liu, C. Xi, C. Dong, C. Cheng, J. Qin, S. Hu, H. Liu and X.-W. Du, *The Journal of Physical Chemistry C*, 2019, **123**, 28319–28326.

83 A. M. Ismail, E. Csapó and C. Janáky, *Electrochimica Acta*, 2019, **313**, 171–178.

84 X. Luo, R. Jiang, Z. Ma, T. Yang, H. Liu, W. Wu, C. Dong and X.-W. Du, *Phys. Chem. Chem. Phys.*, 2022, **24**, 9188.

85 C. Vayenas, S. Bebelis and S. Ladas, *Nature*, 1990, **343**, 625.

86 B. Caglar, A. C. Kizilkaya, J. W. Niemantsverdriet and C. J. Weststrate, *Catalysis, Structure and Reactivity*, 2018, **4**, 1.

87 C. Beasley, M. K. Gnanamani, E. Santillan-Jimenez, M. Martinelli, W. D. Shafer, S. D. Hopps, N. Wanninayake and D.-Y. Kim, *ChemistrySelect*, 2020, **5**, 1013.

88 X. Wang, G. Zhang, L. Yang, E. Sharman and J. Jiang, *WIREs Comput Mol Sci.*, 2018, **8**, e1369.

89 A. Zeradjanin, A. Vimalanandan, G. Polymeros, A. A. Topalov, K. J. J. Mayrhofer and M. Rohwerder, *Phys. Chem. Chem. Phys.*, 2017, **19**, 17019.

90 Y. Zhou and G. Zhang and Q. Ji and W. Zhang and J. Zhang and H. Liu, and J., *Environmental Science and Technology*, 2019, **53**, 11383.

91 N. D. Lang and W. Kohn, *Phys. Rev. B*, 1971, **3**, 1215–1223.

92 N. E. Singh-Miller and N. Marzari, *Phys. Rev. B*, 2009, **80**, 235407.

93 S. De Waele, K. Lejaeghere, M. Sluydts and S. Cottenier, *Phys. Rev. B*, 2016, **94**, 235418.

94 J. C. Boettger, *Phys. Rev. B*, 1996, **53**, 13133.

95 J. Wang and S.-Q. Wang, *Surface Science*, 2014, **630**, 216–224.

96 D.-P. Ji, Q. Zhu and S.-Q. Wang, *Surface Science*, 2016, **651**, 137–146.

97 S.-H. Yoo, N. Siemer, M. Todorova, D. Marx and J. Neugebauer, *The Journal of Physical Chemistry C*, 2019, **123**, 5495–5506.
Tolentino, M. De Santis and C. Mottet, *Phys. Rev. Lett.*, 2008, **100**, 115502.

49 M. Zhou, C. Li and J. Fang, *Chemical Reviews*, 2021, **121**, 736.

50 P. Makkar and N. N. Ghosh, *RSC Adv.*, 2021, **11**, 27897.

51 H. Demir and L. C. Grabow, *ACS Applied Nano Materials*, 2020, **3**, 6127.

52 D. Vollath, F. Fischer and D. Holec, *Beilstein J Nanotechnol.*, 2018, **9**, 2265.

53 H. Amara, J. Nelayah, J. Creuze, A. Chmielewski, D. Alloyeau, C. Ricolleau and B. Legrand, *Phys. Rev. B*, 2022, **105**, 165403.

54 A. Chmielewski, J. Nelayah, H. Amara, J. Creuze, D. Alloyeau, G. Wang and C. Ricolleau, *Phys. Rev. Lett.*, 2018, **120**, 025901.

55 L. Lin, R. Jacobs, T. Ma, D. Chen, J. Booske and D. Morgan, *Phys. Rev. Applied*, 2023, **19**, 037001.

56 F. Calle-Vallejo, M. T. M. Koper and A. S. Bandarenka, *Chem. Soc. Rev.*, 2013, **42**, 5210.

57 S. Trasatti, *J. Electroanal. Chem.*, 1972, **39**, 163.

58 H. Kohlmann, *Journal of Solid State Chemistry*, 2010, **183**, 367.

59 I. Harris, M. Norman and A. W. Bryant, *Journal of the Less-Common Metals*, 1968, **16**, 427.

60 C. Garnero, M. Lepesant, C. Garcia-Marcelot, Y. Shin, C. Meny, P. Farger, B. Warot-Fonrose, R. Arenal, G. Viau, K. Soulantica, P. Fau, P. Poveda, L.-M. Lacroix and B. Chaudret, *Nano Letters*, 2019, **19**, 1379–1386.

61 A. Front and C. Mottet, *J. Phys. Chem. C*, 2021, **125**, 16358–16365.

62 Y. Yang, C.-C. Chen, M. C. Scott, C. Ophus, R. Xu, A. Pryor, L. Wu, F. Sun, W. Theis, J. Zhou, M. Eisenbach, P. R. C. Kent, R. F. Sabirianov, H. Zeng, P. Ercius and J. Miao, *Nature*, 2017, **542**, 75–79.

63 S. Amore, S. Delsante, N. Parodi and G. Borzone, *Thermochimica Acta*, 2009, **481**, 1.

64 A. Ciccioli, G. Balducci, G. Gigli, L. Perring and F. Bussy, *Intermetallics*, 2000, **8**, 195.

65 A. Ormeci, E. Gaudry, M. Armbrüster and Y. Grin, *Chemistry-Open*, 2022, **11**.

66 N. Roy, A. Chakrabarty, B. Koley, T. Saha-Dasgupta and P. P. Jana, *Journal of Solid State Chemistry*, 2020, **290**, 121567.

67 G. Wulff, *Zeitschrift für Kristallographie - Crystalline Materials*, 1901, **34**, 449–530.

68 H. Wang, S. Zhou, K. D. Gilroy, Z. Cai and Y. Xia, *Nano Today*, 2017, **15**, 121–144.

69 G. Casillas, J. J. Velazquez-Salazar and M. Jose-Yacaman, *J. Phys. Chem. C*, 2012, **116**, 8844.

70 A. Front and C. Mottet, *J. Phys. Chem. C*, 2022, 2.

71 A. Front and C. Mottet, *Theoretical Chemistry Accounts*, 2022, **141**, 2.

72 K. Loza, M. Heggen and M. Epple, *Adv. Funct. Mater.*, 2020, **30**, 1909260.

73 B. Farkas and N. H. de Leeuw, *Nanotechnology*, 2020, **31**,

195711.

74 H. Yang, K. Wang, Z. Tang, Z. Liu and S. Chen, *Journal of Catalysis*, 2020, **382**, 181–191.

75 D. Alloyeau, C. Ricolleau, C. Mottet, T. Oikawa, C. Langlois, Y. Le Bouar, N. Braidy and A. Loiseau, *Nature Materials*, 2009, **8**, 940.

76 R. Dronskowski and P. E. Bloechl, *The Journal of Physical Chemistry*, 1993, **97**, 8617–8624.

77 S. Maintz, V. L. Deringer, A. L. Tchougréeff and R. Dronskowski, *Journal of Computational Chemistry*, 2016, **37**, 1030–1035.

78 G. Kresse and J. Hafner, *Physical Review B*, 1993, **47**, 558–561.

79 G. Kresse and J. Furthmüller, *Computational Materials Science*, 1996, **6**, 15–50.

80 G. Kresse and J. Furthmüller, *Physical Review B*, 1996, **54**, 11169–11186.

81 L. Pauling, *The Nature of the Chemical Bond—An Introduction to Modern Structural Chemistry*, Cornell University Press, New-York, 1960.

82 T. Liu, C. Xi, C. Dong, C. Cheng, J. Qin, S. Hu, H. Liu and X.-W. Du, *The Journal of Physical Chemistry C*, 2019, **123**, 28319–28326.

83 A. M. Ismail, E. Csapó and C. Janáky, *Electrochimica Acta*, 2019, **313**, 171–178.

84 X. Luo, R. Jiang, Z. Ma, T. Yang, H. Liu, W. Wu, C. Dong and X.-W. Du, *Phys. Chem. Chem. Phys.*, 2022, **24**, 9188.

85 C. Vayenas, S. Bebelis and S. Ladas, *Nature*, 1990, **343**, 625.

86 B. Caglar, A. C. Kizilkaya, J. W. Niemantsverdriet and C. J. Weststrate, *Catalysis, Structure and Reactivity*, 2018, **4**, 1.

87 C. Beasley, M. K. Gnanamani, E. Santillan-Jimenez, M. Martinelli, W. D. Shafer, S. D. Hopps, N. Wanninayake and D.-Y. Kim, *ChemistrySelect*, 2020, **5**, 1013.

88 X. Wang, G. Zhang, L. Yang, E. Sharman and J. Jiang, *WIREs Comput Mol Sci.*, 2018, **8**, e1369.

89 A. Zeradjanin, A. Vimalanandan, G. Polymeros, A. A. Topalov, K. J. J. Mayrhofer and M. Rohwerder, *Phys. Chem. Chem. Phys.*, 2017, **19**, 17019.

90 Y. Zhou and G. Zhang and Q. Ji and W. Zhang and J. Zhang and H. Liu, and J., *Environmental Science and Technology*, 2019, **53**, 11383.

91 N. D. Lang and W. Kohn, *Phys. Rev. B*, 1971, **3**, 1215–1223.

92 N. E. Singh-Miller and N. Marzari, *Phys. Rev. B*, 2009, **80**, 235407.

93 S. De Waele, K. Lejaeghere, M. Sluydts and S. Cottenier, *Phys. Rev. B*, 2016, **94**, 235418.

94 J. C. Boettger, *Phys. Rev. B*, 1996, **53**, 13133.

95 J. Wang and S.-Q. Wang, *Surface Science*, 2014, **630**, 216–224.

96 D.-P. Ji, Q. Zhu and S.-Q. Wang, *Surface Science*, 2016, **651**, 137–146.

97 S.-H. Yoo, N. Siemer, M. Todorova, D. Marx and J. Neugebauer, *The Journal of Physical Chemistry C*, 2019, **123**, 5495–5506.
98 R. Tran, X.-G. Li, J. H. Montoya, D. Winston, K. A. Persson and S. P. Ong, *Surface Science*, 2019, **687**, 48–55.

99 T. T. Dorini, F. Brix, C. Chatelier, A. Kokalj and E. Gaudry, *Nanoscale*, 2021, **13**, 10771–10779.

100 Y. Zhang, O. Pluchery, L. Caillard, A.-F. Lamic-Humblot, S. Casale, Y. J. Chabal and M. Salmeron, *Nano Letters*, 2015, **15**, 51.

101 E. Kalered, N. Brenning, I. Pilch, L. Caillault, T. Minea and L. Ojamae, *Physics of Plasmas*, 2017, **24**, cddc.

102 S. C. Fain and J. M. McDavid, *Phys. Rev. B*, 1974, **9**, 5099.

103 P. D. Swartzentruber, M. J. Detisch and T. J. Balk, *J. Vac. Sci. Technol.*, 2015, **A 33**, 021405.

104 K. F. Wojciechowski, *Phys. Rev. B*, 1999, **60**, 9202.

105 T. G. Akhmetshina, V. A. Blatov, D. M. Proserpio and A. P. Shevchenko, *Accounts of Chemical Research*, 2018, **51**, 21.

106 X. Liu, X. Wen and R. Hoffmann, *ACS Catalysis*, 2018, **8**, 3365.

107 P. E. Blöchl, *Physical Review B*, 1994, **50**, 17953–17979.

108 G. Kresse and D. Joubert, *Physical Review B*, 1999, **59**, 1758–1775.

109 J. P. Perdew, K. Burke and M. Ernzerhof, *Phys. Rev. Lett.*, 1996, **77**, 3865.

110 J. P. Perdew, K. Burke and M. Ernzerhof, *Phys. Rev. Lett.*, 1997, **78**, 1396.

111 A. Stukowski, *Modelling and Simulation in Materials Science and Engineering*, 2009, **18**, 015012.

112 A. H. Larsen, J. J. Mortensen, J. Blomqvist, I. E. Castelli, R. Christensen, M. Dułak, J. Friis, M. N. Groves, B. Hammer, C. Hargus, E. D. Hermes, P. C. Jennings, P. B. Jensen, J. Kermode, J. R. Kitchin, E. L. Kolsbjerg, J. Kubal, K. Kaasbjerg, S. Lysgaard, J. B. Maronsson, T. Maxson, T. Olsen, L. Pastewka, A. Peterson, C. Rostgaard, J. Schiøtz, O. Schütt, M. Strange, K. S. Thygesen, T. Vegge, L. Vilhelmsen, M. Walter, Z. Zeng and K. W. Jacobsen, *Journal of Physics: Condensed Matter*, 2017, **29**, 273002.

113 K. Mathew, A. K. Singh, J. J. Gabriel, K. Choudhary, S. B. Sinnott, A. V. Davydov, F. Tavazza and R. G. Hennig, *Computational Materials Science*, 2016, **122**, 183–190.

114 S. P. Ong, W. D. Richards, A. Jain, G. Hautier, M. Kocher, S. Cholia, D. Gunter, V. L. Chevrier, K. A. Persson and G. Ceder, *Computational Materials Science*, 2013, **68**, 314–319.

115 Émilie Gaudry, in *Comprehensive Inorganic Chemistry III (Third Edition)*, ed. J. Reedijk and K. R. Poeppelmeier, Elsevier, Oxford, Third Edition edn., 2023, pp. 74–104.

116 M. Krone, J. Stone, T. Ertl and K. Schulten, *EuroVis - Short Papers*, 2012.
98 R. Tran, X.-G. Li, J. H. Montoya, D. Winston, K. A. Persson and S. P. Ong, *Surface Science*, 2019, **687**, 48–55.

99 T. T. Dorini, F. Brix, C. Chatelier, A. Kokalj and E. Gaudry, *Nanoscale*, 2021, **13**, 10771–10779.

100 Y. Zhang, O. Pluchery, L. Caillard, A.-F. Lamic-Humblot, S. Casale, Y. J. Chabal and M. Salmeron, *Nano Letters*, 2015, **15**, 51.

101 E. Kalered, N. Brenning, I. Pilch, L. Caillault, T. Minea and L. Ojamae, *Physics of Plasmas*, 2017, **24**, cddc.

102 S. C. Fain and J. M. McDavid, *Phys. Rev. B*, 1974, **9**, 5099.

103 P. D. Swartzentruber, M. J. Detisch and T. J. Balk, *J. Vac. Sci. Technol.*, 2015, **A 33**, 021405.

104 K. F. Wojciechowski, *Phys. Rev. B*, 1999, **60**, 9202.

105 T. G. Akhmetshina, V. A. Blatov, D. M. Proserpio and A. P. Shevchenko, *Accounts of Chemical Research*, 2018, **51**, 21.

106 X. Liu, X. Wen and R. Hoffmann, *ACS Catalysis*, 2018, **8**, 3365.

107 P. E. Blöchl, *Physical Review B*, 1994, **50**, 17953–17979.

108 G. Kresse and D. Joubert, *Physical Review B*, 1999, **59**, 1758–1775.

109 J. P. Perdew, K. Burke and M. Ernzerhof, *Phys. Rev. Lett.*, 1996, **77**, 3865.

110 J. P. Perdew, K. Burke and M. Ernzerhof, *Phys. Rev. Lett.*, 1997, **78**, 1396.

111 A. Stukowski, *Modelling and Simulation in Materials Science and Engineering*, 2009, **18**, 015012.

112 A. H. Larsen, J. J. Mortensen, J. Blomqvist, I. E. Castelli, R. Christensen, M. Dułak, J. Friis, M. N. Groves, B. Hammer, C. Hargus, E. D. Hermes, P. C. Jennings, P. B. Jensen, J. Kermode, J. R. Kitchin, E. L. Kolsbjerg, J. Kubal, K. Kaasbjerg, S. Lysgaard, J. B. Maronsson, T. Maxson, T. Olsen, L. Pastewka, A. Peterson, C. Rostgaard, J. Schiøtz, O. Schütt, M. Strange, K. S. Thygesen, T. Vegge, L. Vilhelmsen, M. Walter, Z. Zeng and K. W. Jacobsen, *Journal of Physics: Condensed Matter*, 2017, **29**, 273002.

113 K. Mathew, A. K. Singh, J. J. Gabriel, K. Choudhary, S. B. Sinnott, A. V. Davydov, F. Tavazza and R. G. Hennig, *Computational Materials Science*, 2016, **122**, 183–190.

114 S. P. Ong, W. D. Richards, A. Jain, G. Hautier, M. Kocher, S. Cholia, D. Gunter, V. L. Chevrier, K. A. Persson and G. Ceder, *Computational Materials Science*, 2013, **68**, 314–319.

115 Émilie Gaudry, in *Comprehensive Inorganic Chemistry III (Third Edition)*, ed. J. Reedijk and K. R. Poeppelmeier, Elsevier, Oxford, Third Edition edn., 2023, pp. 74–104.

116 M. Krone, J. Stone, T. Ertl and K. Schulten, *EuroVis - Short Papers*, 2012.