A Comparative Study of Oxygen and Hydrogen Adsorption on Strained and Alloy-Supported Pt(111) Monolayers

: A comparative study of the unreacted and reacted uniaxially strained Pt(111) and the layered (111)-Pt/Ni/Pt 3 Ni and (111)-Pt/Ni/PtNi 3 surfaces has been performed using density functional theory (DFT). An in-depth study of the unreacted surfaces has been performed to evaluate the importance of geometric, magnetic and ligand effects in determining the reactivity of these different Pt surfaces. An analysis of the binding energies of oxygen and hydrogen over the high-symmetry binding positions of all surfaces has been performed. The study has shown that O and H tend to bind more strongly to the (111)-Pt/Ni/Pt 3 Ni surface and less strongly to the (111)-Pt/Ni/PtNi 3 surface compared to binding on the equivalently strained Pt(111) surfaces. Changes in the surface magnetisation of the surfaces overlaying the ferromagnetic alloys during adsorption are discussed, as well as the behaviour of the d-band centre across all surfaces, to evaluate the potential mechanisms for these differences in binding. An accompanying comparison of the accessible density functionals has been included to estimate the error in the computational binding energies.


Introduction
Understanding the mechanisms of the use of the oxygen reduction reaction (ORR) and the hydrogen evolution reaction (HER) is key to the development of modern energy storage and transport materials. The technology inherent in these materials contributes significantly to sustainable chemical production, whose energy footprint is an imperative concern when designing catalysts and networks of catalysts.
There are several different ways that catalytic surfaces can be trained towards a particular reaction sequence or energy budget. Lattice strain control [1] has been conceptually available for several decades [2] and generally involves applying biaxial or, less commonly, uniaxial strain within the surface plane of the catalyst. There are a number of possible consequences of this action, and their relative importance depends sensitively on a number of factors, including the stoichiometry and structure of the catalyst under investigation, the working conditions of the catalysts, the way that the strain is applied and the type of investigation being performed. In the latter point, the main distinction made in the current work is between investigations which are performed experimentally and those that are performed theoretically, or computationally. In these latter investigations, systems are reduced in their complexity, such as the trial systems commonly presented in density functional theory (DFT) investigations. Other computational methods can model larger systems and investigate larger scale phenomena, e.g., reaction processes around extended surface clusters, but the computational tools used, such as empirical force fields, require considerable refinement and training to produce simulations which are physically meaningful.
The current manuscript is a purely computational work. However, it is worthwhile to review some of the more important experimental phenomena observed to reinforce the importance of some of the investigations that will be performed in the current work and how applying strain can enable some of these phenomena. Fundamentally, there are three classes of phenomena: ensemble, ligand and geometric. Ensemble effects occur when ferromagnetism and its relative importance to the externally applied strain have not been fully investigated in the original work [10] or in subsequent publications but may prove significant.
Based on the current published database of catalyst performance, some systematic design protocols have been proposed. Examples of the design protocols have been based on the chemical composition, atomic coordination, and strains of the top surface layers of nanocatalysts [11]. Other examples use strain engineering and localized lattice strain [12], which focuses development on the d-band model of surface reactivity and the known high reactivity of defects compared to defect-free surfaces and facets. The importance of defects and subsurface vacancies has been shown in fundamental studies of Pt-based catalysts [13]. Further protocols have focused on ensemble effects in Pt overlayer [14] using a kinetic-mapping technique.
Whilst these broader characterizations have demonstrated some success, there is still the need to develop a fundamental understanding of Pt-based catalysts-for example, in the empirical modelling of the pure metal under strain [15]-as a characterization and potential systematic design tool. This approach has been applied, in part, in recent studies which evaluated both the strain and the ligand effects [16] by systematically geometrically constraining Pt-bearing nanoparticles but allowing adsorbates to fully relax on them. These studies identified Pt-Au-M nanoparticles (M = Cr, Mn, Co, Cu, Zn) as promising electrocatalysts for the ORR, due to their inherently weak Pt-OH binding with respect to a pure Pt nanoparticle. This fundamental approach highlights the importance of single crystal methods, e.g., [17], even amongst the literature which is increasingly dominated by nanoparticle and multiple-facet studies, and as an essential part of the important bridging studies [18,19] between the more idealized single crystal methodology and the less computationally accessible nanoparticle simulations. One recent example of single crystal studies has been a broad survey [20] of low index transition metal surfaces and their potential reactivity towards the HER. This survey concluded that under applied strain, the hydrogen atoms across a range of surfaces could migrate between preferred binding sites with a concurrent change in local valence, binding energy, and potentially significant changes to their reactivity.
The current work will focus on density functional theory investigations of single crystal systems. This approach has very recently proven successful in determining the stability of Pt skins overlying a range of intermetallic layers [21] and in the activity of Pt/PtV(111) towards the ORR [22]. The current investigations will focus on the multilayered (111)-Pt/Ni/Pt 3 Ni and (111)-Pt/Ni/PtN 3 i using a full spin-polarized approach. Recent studies [23] of the magnetic character of the ordered bulk L 10 and L 12 phases of Pt 3 Ni, PtNi and PtNi 3 using both scalar and vector relativistic density functional theory (DFT) have shown that the effect of applying compressive or tensile strain to each alloy is to decrease or increase the magnetism of the alloy, respectively. The Pt x Ni 1-x bulk alloys considered in that work [23] are part of a series of experimentally accessible alloys, which include Pt x M 1-x (M = Fe,Co), and therefore, each contain a permanently magnetised (Ni, Fe or Co) component. Comparative studies of these materials are therefore beneficial as they can highlight trends within this class of bimetallics. The response of the magnetic moment of the L 10 and L 12 ordered phases of Pt 3 M, PtM and PtM 3 (M = Fe,Co) to compressive and tensile strain [24] has shown that the magnetic moment of the Pt 3 M and PtM phases varies linearly with applied strain. The response of the PtM 3 alloys shows significant transition in the rate of change of magnetic moment at approximately zero strain. The magnetic moment changes of Pt x Fe 1-x and Pt x Co 1-x during strain arise from intraorbital charge transfer between d orbitals. These background studies demonstrate some of the characteristics of the ordering of the substrate layers considered in the current work.
The computational investigations in the current work have been prompted by the experimental observation [25] that the Pt 3 Ni(111) activity towards the ORR is 10-fold the activity seen on Pt(111), and similar observations have been made on Pt 3 M surfaces [26]. Due to the technological importance of both the HER and ORR systematic investigations, a systematic series of DFT investigations have been completed for the strained H/Pt(111) [27], O/Pt(111) and OH/Pt(111) [28] systems. These investigations have shown that for the H patterned surfaces, and for smaller unit cells of the oxygenated surfaces, applying biaxial strain in the surface (111) plane changes the lowest energy binding position of the adsorbate. This result is of fundamental importance, particularly with adsorbates that have directional bonds with the substrate, such as oxygen, because changing the valence of these adsorbates with the surface will significantly change the available valence charge of that adsorbate to other reacting molecules.
Single crystal experimental studies of the Pt 3 Ni(100), (110) and (111) [29] have demonstrated that some refinement of the single overlayer model routinely used in computational investigations in this field is necessary. After thermal processing of nanoparticle Pt 3 Ni, a conventional core-shell structure develops composed of pure metal overlayers and an alloy. However, in the case of the Pt x Ni 1-x , the Pt surface layer has been shown to lie on a purely Ni second (selvedge) layer. The current work will therefore investigate the effect of this magnetic selvedge layer. The current work will also investigate the effects of unit cell size during the oxygenation of Pt(111) under strain and will further provide a benchmark equation of state and surface energy data for the clean Pt(111) under strain. These investigations will quantitatively evaluate the density functional used and seek to establish a broad computational error bar in the binding energies and bond lengths, determined using both local density approximation (LDA) and the generalised gradient approximation (GGA).

Computational Methods
The investigations performed in this work used the Quantum Espresso density functional theory (DFT) package [30]. Spin-polarized simulations were performed using both the local density approximation (LDA) and the generalized gradient approximation (GGA). Hybrid functionals, which would nominally have required an order of magnitude more computational time, were not considered in the current work.
The reason that the two sets of calculations were performed was to estimate the computational error in the current set of investigations, which will assist experimental teams who may wish to compare their experiment results with the results presented in the current work. The reason that these two approximations were used was because the current work focusses on the bond lengths and the cohesive, surface and adsorption energies. The LDA is conventionally found to overestimate cohesive energies and underestimate bond lengths. The GGA tends to overcompensate for the LDA overbinding and consequently overestimates the bond lengths. Comparative studies of the two functionals can then at least approximately provide an upper and lower estimate of the expected values of these quantities. Furthermore, the two functionals are computationally accessible for the size of systems that are being investigated in the current work.
The Perdew-Zunger (PZ) LDA and the Perdew-Burke-Ernzerhof (PBE) GGA exchangecorrelation functionals were used to generate the norm-conserving pseudo-potentials [31] used in the current simulations. A wave-function kinetic energy cut-off of 75 Ry and a charge density/potential cut-off of 300 Ry were used throughout these investigations, and a Brillouin zone sampling of (10 × 10 × 10) and (6 × 6 × 1) was used for the bulk and surface investigations, respectively. A first-order Methfessel-Paxton smearing of 0.02 Ry is used throughout this work [32].
The cohesive energy E coh of the bulk Pt crystal is defined in Equation (1): where E bulk is the energy of the primitive bulk crystal, which contains 1 Pt atom/unit cell. E iso is the energy of a single Pt atom in the gas phase, which was calculated by setting the atom in an orthogonal box with dimensions 10 × 10 × 10 Å 3 . (1 × 1 × 1) unit cells were used to perform the bulk Pt and bulk alloy investigations. The surface energy E sur f of the clean Pt(111) slab was calculated using Equation (2): where A is the total surface area of a slab containing n atoms with a total energy E slab .
The energies E of any density of states which have been considered in this work are normalized to the Fermi level E F , i.e., E F = 0 eV. The centre E c and width w of each of the density of states P(E) are defined in terms of their first and second moment, respectively [33].
Strain σ was defined in Equation (5) in terms of general lattice constants L and L 0 , where L 0 is the lattice constant L at equilibrium.
We applied the strains biaxially throughout the current work, so the surface lattice constant along the 112 and 110 directions was changed by the same fractional amount when strain was applied to surface FCC (111) simulations. All atoms in each strained slab were allowed to fully relax, except for the central atomic layers. Each Pt(111) simulation had 6 Pt layers, whilst each alloy simulation had 7 layers of Pt 3 Ni or PtNi 3 alloy with further Ni and Pt layers-detailed later in the text-on each side of the slab. Subsequent slabs were separated by a vacuum with a width of approximately ten interplanar distances. For bulk simulations, the strains were applied identically and simultaneously along the [100], [010] and [001] directions. In the bulk case, however, no subsequent relaxation was considered.
To simulate the adsorption of atoms, the slabs were patterned symmetrically by attaching adsorbate atoms-either oxygen or hydrogen-in the same binding position on either side of a symmetric slab. The adsorption energy of E ads was then defined as where E A/S is the Kohn-Sham energy of the adsorbate patterned slab, E S is the energy of the unpatterned slab and E A is the energy of the adsorbate. The factors of 2 and 1 /2 account for the patterning by adsorbates on both faces of the slab. Symmetric slabs were used throughout the current-with one minor exception, which is discussed in the text. This means that all quantitative outcomes were calculated using slabs which were either entirely unreacted or which were reacted identically adsorbate on either face. (3 × 3) and (2 × 2) surface unit cells were used to simulate the reacted Pt(111), and the Pt/Ni/Pt 3 Ni-(111) and Pt/Ni/PtNi 3 -(111) surfaces, respectively.

Results and Discussion
This section is structured in the following way: first, the bulk Pt and unreacted Pt(111) surface are investigated comparatively using the LDA and GGA exchange correlation functionals before investigating the reactivity of the Pt(111) surface towards oxygen and hydrogen. A similar methodology is then performed to investigate the more complex Pt/Ni/Pt 3 Ni-(111) and Pt/Ni/PtNi 3 -(111) surfaces before investigating their reactivity towards oxygen and hydrogen. The effects of adsorption on these surfaces are discussed with some particular focus on the magnetisation of the various surface, the origins of this magnetisation and the effects of adsorption on the magnetisation.

Unreacted Bulk Pt and Pt(111)
The orthogonal unit cell of the bulk FCC Pt crystal is shown in Figure 1.
lation functionals before investigating the reactivity of the Pt(111) surface towards oxygen and hydrogen. A similar methodology is then performed to investigate the more complex Pt/Ni/Pt3Ni-(111) and Pt/Ni/PtNi3-(111) surfaces before investigating their reactivity towards oxygen and hydrogen. The effects of adsorption on these surfaces are discussed with some particular focus on the magnetisation of the various surface, the origins of this magnetisation and the effects of adsorption on the magnetisation.

Unreacted Bulk Pt and Pt(111)
The orthogonal unit cell of the bulk FCC Pt crystal is shown in Figure 1. To determine the equilibrium lattice parameter a, the cohesive energy Ecoh of the unit cell was calculated at different strains σ. The resulting equations of state are shown in Figure 2 and were calculated using the LDA and GGA approximations. The equilibrium lattice spacing occurs at the minimum of these curves. To determine the equilibrium lattice parameter a, the cohesive energy E coh of the unit cell was calculated at different strains σ. The resulting equations of state are shown in Figure 2 and were calculated using the LDA and GGA approximations. The equilibrium lattice spacing occurs at the minimum of these curves. lation functionals before investigating the reactivity of the Pt(111) surface towards oxygen and hydrogen. A similar methodology is then performed to investigate the more complex Pt/Ni/Pt3Ni-(111) and Pt/Ni/PtNi3-(111) surfaces before investigating their reactivity towards oxygen and hydrogen. The effects of adsorption on these surfaces are discussed with some particular focus on the magnetisation of the various surface, the origins of this magnetisation and the effects of adsorption on the magnetisation.

Unreacted Bulk Pt and Pt(111)
The orthogonal unit cell of the bulk FCC Pt crystal is shown in Figure 1. To determine the equilibrium lattice parameter a, the cohesive energy Ecoh of the unit cell was calculated at different strains σ. The resulting equations of state are shown in Figure 2 and were calculated using the LDA and GGA approximations. The equilibrium lattice spacing occurs at the minimum of these curves.  Table 1 summarizes the bulk lattice parameters, equilibrium cohesive energies and magnetic moments for both the bulk Pt system (shown in Figure 1) and for the bulk Ni system, which has the same FCC structure as Pt. It can be seen from Table 1 that the LDA tends to underestimate the bond lengths compared to the experimental values. To compensate for this, the GGA estimates a lower cohesive energy compared to the LDA cohesive energies, and this can be seen in Figure 1. This lower estimation of the energy predicts longer equilibrium bond lengths for the GGA in Table 1. Figure 3 shows the (a) Pt(111) surface and (b) the same slab in cross-sectional view. When strain is applied to the surface, the lattice constant along both the orthogonal 7 of 23 112 and 110 directions was changed by the same fractional amount. Once this strain was applied, the atoms within the slab could relax freely, except for the central two layers which were constrained. Table 1. The bulk lattice parameters a, cohesive energies E coh and magnetic moments µ Pt and µ Ni for the bulk Ni and Pt metals. Cohesive energies are at the minima in the equations of state shown in Figure 2. The 'Theory/Experiment' column shows the exchange-correlation functional used (LDA/GGA) in the case of theoretical estimates of each of the physical parameter or shows the experimental value ('Exptl') where available.

Metal
Theory/Experiment  The resulting surface energies and layer spacings are shown in Figures 4 and 5, respectively. Esurf is overestimated by the LDA in the same way that the LDA cohesive energies were overestimated in Figure 2, with the GGA reducing the overestimation in both cases.   The results and analysis in this section have so far focused on the energetic and th structural character of the bulk and unreacted pure metal surfaces. In order to develop th description of these systems, the projected density of states is now considered. Followin Equations (3) and (4)     The results and analysis in this section have so far focused on the energetic and t structural character of the bulk and unreacted pure metal surfaces. In order to develop t description of these systems, the projected density of states is now considered. Followin Equations (3) and (4) in the 'Computational Methods' section, the key metrics which w be used in the current work are the centre Ec and width w of the density of states. Figu 6 shows explicitly how these quantities develop for bulk Pt. The centre of both the d an  The results and analysis in this section have so far focused on the energetic and the structural character of the bulk and unreacted pure metal surfaces. In order to develop the description of these systems, the projected density of states is now considered. Following Equations (3) and (4) in the 'Computational Methods' section, the key metrics which will be used in the current work are the centre E c and width w of the density of states. Figure 6 shows explicitly how these quantities develop for bulk Pt. The centre of both the d and the sp bands approaches the Fermi level as the strain σ increases, i.e., as the atoms in the bulk become increasingly separated. Accompanying this shift in E c , the width w narrows. These two phenomena can be interpreted as the atoms approaching a more gas-like state as σ increases and are indicative that over the range of strains considered for the bulk Pt phase, a structural or electronic phase transition might not ordinarily be expected.  Figure 7 shows the surface projected density of states (PDOS) in a strained P slab at strains which correspond to the surface lattice constants of (111)-Pt3Ni and and shows little evidence that the surface Pt atoms become magnetised simply d strain alone.  Figure 7 shows the surface projected density of states (PDOS) in a strained Pt(111) slab at strains which correspond to the surface lattice constants of (111)-Pt 3 Ni and PtNi 3 and shows little evidence that the surface Pt atoms become magnetised simply due to strain alone.
The layer-wise variation in the centre E c and width w of the density of states under strain is summarized in Figure 8 which shows that the surface Pt E c approach the Fermi level as strain becomes increasingly tensile and the distance between the atoms becomes greater. This effect is largely independent of the functional. For the deeper level atoms, the effect is only apparent for the more directional d bonds. The sp band is less directional, containing the diffuse s electrons which are less sensitive to the changes in the lattice constant. This is evidenced in the small gradient of the second and deeper layer atoms in Figure 8b. The layer-wise variation in the centre Ec and width w of the density of states und strain is summarized in Figure 8 which shows that the surface Pt Ec approach the Ferm level as strain becomes increasingly tensile and the distance between the atoms becom greater. This effect is largely independent of the functional. For the deeper level atom the effect is only apparent for the more directional d bonds. The sp band is less directiona

Reacted Pt(111)
In the previous section, the energetic, structural, magnetic and electronic character of the Pt(111) surface under strain was presented. In the current section, the consequences of this character will be demonstrated in the response of the surface to oxygen and hydrogen. Figure 9 shows the high-symmetry binding positions on the Pt(111) surface together with the (3 × 3) unit cell used in the current simulations. Figure 10 shows the adsorption energy under strain σ for (3 × 3)-O/Pt(111) and (3 × 3)-H/Pt(111). The oxygen atoms are seen to preferentially occupy the FCC sites for all strains. This observation is different from previous simulations [28], which showed that for compressive strains (σ < 0.97), the oxygen atoms preferentially occupy a bridge site. However, these earlier studies used a smaller (2 × 2) surface unit cell and observed a reconstruction of the surface, particularly

Reacted Pt(111)
In the previous section, the energetic, structural, magnetic and electronic character of the Pt(111) surface under strain was presented. In the current section, the consequences of this character will be demonstrated in the response of the surface to oxygen and hydrogen. Figure 9 shows the high-symmetry binding positions on the Pt(111) surface together with the (3 × 3) unit cell used in the current simulations. Figure 10 shows the adsorption energy under strain σ for (3 × 3)-O/Pt(111) and (3 × 3)-H/Pt(111). The oxygen atoms are seen to preferentially occupy the FCC sites for all strains. This observation is different from previous simulations [28], which showed that for compressive strains (σ < 0.97), the oxygen atoms preferentially occupy a bridge site. However, these earlier studies used a smaller (2 × 2) surface unit cell and observed a reconstruction of the surface, particularly for the related case of hydroxyl adsorption. This reconstruction does lead to a lower energy binding position-the values of E ads in the earlier work are noted to be of the order of 0.5 eV lower than those in the current work. Consequentially, oxygen adsorption on strained Pt(111) should be considered in low coverage (current work) and high coverage regimes. The interactions in the high coverage regime are not fully characterized, but the current evidence does suggest that lateral O-Pt interactions are significant. for the related case of hydroxyl adsorption. This reconstruction does lead to a lower energy binding position-the values of Eads in the earlier work are noted to be of the order of 0.5 eV lower than those in the current work. Consequentially, oxygen adsorption on strained Pt(111) should be considered in low coverage (current work) and high coverage regimes. The interactions in the high coverage regime are not fully characterized, but the current evidence does suggest that lateral O-Pt interactions are significant.  To provide a direct quantitative comparison between oxygen and hydrogen adsorption on the strained Pt(111) surface considered in the current section and adsorption on the alloybased surfaces considered later in this work, Table 2 shows the adsorption energy E ads for the unstrained (σ = 1) and compressively strained surfaces. The results in this table are consistent with those seen earlier in the current work-for example, in Figures 2 and 4-where the LDA is seen to overbind and consequently predict larger E ads . The results in Table 2, more importantly, can be used to discuss the error in this energy, which arises from the choice of functional. Generally, oxygen is more tightly bound to the Pt(111) surface independent of the strain state.   To provide a direct quantitative comparison between oxygen and hydrogen adsorption on the strained Pt(111) surface considered in the current section and adsorption on the alloy-based surfaces considered later in this work, Table 2 shows the adsorption energy Eads for the unstrained (σ = 1) and compressively strained surfaces. The results in this table are consistent with those seen earlier in the current work-for example, in Figures 2  and 4-where the LDA is seen to overbind and consequently predict larger Eads. The results in Table 2, more importantly, can be used to discuss the error in this energy, which arises from the choice of functional. Generally, oxygen is more tightly bound to the Pt(111) surface independent of the strain state.

Unreacted Pt x Ni 1-x and Pt/Ni/Pt x Ni 1-x
Thus far, the results in the current work have focused on pure Pt in either bulk or slab form. This section now investigates the unreacted Pt 3 Ni and PtNi 3 alloys in both bulk and slab form. Figure 11 shows the bulk unit cells of these alloys.  Table 3 shows the equilibrium lattice parameters, cohesive energies and magnetic moments for both unit cells. The Pt atoms can be seen to have developed a permanent magnetic moment in the presence of Ni, which is particularly noticeable for the PtNi3 alloy.  [34] 3.645 Using the bulk lattice parameters presented in Table 3, the surfaces of unreacted (111)-Pt/Ni/Pt3Ni and (111)-Pt/Ni/PtNi3 were fully relaxed. Based on the experimental studies of these systems [29], a Ni selvedge layer is set beneath the surface Pt layer. However, because of the complexity of the bulk alloy phase, the surface and selvedge layers are comprised of non-equivalent atoms. This is show in Figure 12a,b, where a subscript notation is adopted to identify the different types of surface and second layer Pt and Ni atom.  Table 3 shows the equilibrium lattice parameters, cohesive energies and magnetic moments for both unit cells. The Pt atoms can be seen to have developed a permanent magnetic moment in the presence of Ni, which is particularly noticeable for the PtNi 3 alloy. Using the bulk lattice parameters presented in Table 3, the surfaces of unreacted (111)-Pt/Ni/Pt 3 Ni and (111)-Pt/Ni/PtNi 3 were fully relaxed. Based on the experimental studies of these systems [29], a Ni selvedge layer is set beneath the surface Pt layer. However, because of the complexity of the bulk alloy phase, the surface and selvedge layers are comprised of non-equivalent atoms. This is show in Figure 12a,b, where a subscript notation is adopted to identify the different types of surface and second layer Pt and Ni atom.  Table 4 summarizes the interlayer spacing for the alloy systems. The variation in height between the different surface Pt and second layer Ni species is averaged across each layer.  Figure 12 of the unreacted (111)-Pt/Ni/Pt 3 Ni and (111)-Pt/Ni/PtNi 3 surfaces. The 'Functional' column shows the exchangecorrelation functional used (LDA/GGA). All lengths are given in units of Å.  Figure 13 shows the planar-averaged spin density ρ ↑ − ρ ↓ for the clean (111)-Pt/Ni/Pt 3 Ni and (111)-Pt/Ni/PtNi 3 slabs. Each slab is composed of 7 layers of alloy bounded by Ni selvedge and Pt surface layers, and the averaging is taken parallel to the (111) plane.

Alloy
on-top, FCC and HCP surface binding positions and three non-equivalent 2-fold bridge sites. Table 4 summarizes the interlayer spacing for the alloy systems. The varia height between the different surface Pt and second layer Ni species is averaged each layer. Table 4. The interlayer spacing d12, d23, d34, d45 and d56 defined in Figure 12 of the unreacte Pt/Ni/Pt3Ni and (111)-Pt/Ni/PtNi3 surfaces. The 'Functional' column shows the exchange-cor functional used (LDA/GGA). All lengths are given in units of Å.  Figure 13 shows the planar-averaged spin density 〈 ↑ − ↓ 〉 for the clean Pt/Ni/Pt3Ni and (111)-Pt/Ni/PtNi3 slabs. Each slab is composed of 7 layers o bounded by Ni selvedge and Pt surface layers, and the averaging is taken parallel (111) plane. Quantitatively, the Pt surface layer has an average spin polarization of 0.006 both slabs shown in Figure 13. The Ni second layer has an average spin polariza 0.051 Å −1 and 0.061 Å −1 for the Pt3Ni and PtNi3 alloys, respectively. The values for atoms are consistent with the ferromagnetic moment that the Ni atoms carry. The in magnetisation in the surface Pt layers is significant, particularly when it is compa the negligible moment carried by the surface atoms in a pure Pt(111) slab, which w cussed earlier in this manuscript. The consequences of this induced moment in the s Pt layer for the adsorption of oxygen and hydrogen and the effects of adsorption moment will be discussed in the next section. Quantitatively, the Pt surface layer has an average spin polarization of 0.006 Å −1 for both slabs shown in Figure 13. The Ni second layer has an average spin polarization of 0.051 Å −1 and 0.061 Å −1 for the Pt 3 Ni and PtNi 3 alloys, respectively. The values for the Ni atoms are consistent with the ferromagnetic moment that the Ni atoms carry. The induced magnetisation in the surface Pt layers is significant, particularly when it is compared to the negligible moment carried by the surface atoms in a pure Pt(111) slab, which was discussed earlier in this manuscript. The consequences of this induced moment in the surface Pt layer for the adsorption of oxygen and hydrogen and the effects of adsorption on this moment will be discussed in the next section.   Table 5 shows the adsorption energies Eads at the high symmetry sites for oxyge hydrogen. These energies can be compared directly with the adsorption energies f strained Pt(111) surface presented in Table 2. For oxygen adsorption, the adsorpti ergies are between 0.176 and 1.081 eV greater on the (111)-Pt/Ni/Pt3Ni surface than  Table 5 shows the adsorption energies E ads at the high symmetry sites for oxygen and hydrogen. These energies can be compared directly with the adsorption energies for the strained Pt(111) surface presented in Table 2. For oxygen adsorption, the adsorption energies are between 0.176 and 1.081 eV greater on the (111)-Pt/Ni/Pt 3 Ni surface than on the strained Pt(111) surface, across both functionals. This shows that a ligand effect exists on the (111)-Pt/Ni/Pt 3 Ni surface which causes the oxygen atoms to bind more strongly to that surface than to the strained Pt(111) surface. Conversely, on the (111)-Pt/Ni/PtNi 3 surface, the adsorption energies are 0.863 to 1.874 eV lower than on the equivalently strained Pt(111) surface. The ligand effect does have some depth sensitivity. It does not depend on the nearest neighbour interactions between the adsorbate and the surface Pt atoms and does extend even beyond the second nearest-neighbour interactions involving the Ni atoms in the second (selvedge) layer. For hydrogen adsorption, a similar set of energy differences can be determined. The hydrogen atoms bind between 0.320 and 0.636 eV more strongly to the (111)-Pt/Ni/Pt 3 Ni surface than to the strained Pt(111) surface, and between 0.897 and 1.437 eV less strongly to the (111)-Pt/Ni/PtNi 3 surface than to the strained Pt(111) surface.    Table 6 compares the nearest-neighbour O-Pt and H-Pt bond lengths obtained using both the LDA and the GGA functionals. The LDA predicts shorter bond lengths than the GGA in almost all cases. This is indicative of a functional which is overestimating the bonding and therefore predicting shorter bond lengths than might be observed experimentally. Conversely the GGA tends to underestimate the binding and therefore predicts generally longer bond lengths. For these reasons, experimentalists might anticipate measuring a bond length between the LDA and GGA values presented in this table.

Alloy
The next step in characterizing the oxygen and hydrogen adsorption is to look in more detail at the density of states of the adsorbate and substrate for the alloy-bearing systems, in the same way that these quantities were analysed for the Pt(111) systems earlier in this work. Table 7 shows the first and second moments of the surface Pt atoms for these systems. As the surface Pt atoms in the strained Pt(111) are not magnetised, both the centre and width of the spin up and down components are identical, so only one value is quoted for each for Pt(111). For both alloys, the width of the spin polarized d-band orbital is larger than the equivalent width for the strained Pt(111). However, the behaviour of the centre of the d-band is less clearly defined. For the Pt/Ni/Pt 3 Ni, the E c;Pt are close to the E c↓ values for both the LDA and GGA simulations. This result is surprising within the conventional d-band model, as it was noted earlier in this section that comparisons of the adsorption energy between the Pt/Ni/Pt 3 Ni and strained Pt(111) varied significantly with the density functional. For Pt/Ni/PtNi 3 , the situations appear more consistent with the d-band model, as there is some variation in the E c;Pt and the E c↑ and E c↓ for both the LDA and GGA simulations, with the E c;Pt changing alignment with the functional, which is more consistent with the d-band model.    A further observation is that the magnetisation of the surface Pt atoms is generally reduced by both oxygen and hydrogen adsorption. This is shown qualitatively in Figure 15, where the surface average spin moment of the Pt atoms in the reacted surface is less than that on the clean surface. Table 7. The spin-resolved first and second moments of the d-band of the surface Pt atoms in the Pt3Ni and PtNi3 alloy-bearing systems and in the surface Pt atoms in strained Pt(111). The up and down arrows denote the spin up and down components of the surface Pt atoms in the alloy systems, and 'Pt' denotes the surface Pt atoms in strained Pt(111). The number in brackets after the Ec;Pt and is the strain σ at which the quantity was determined. The 'Functional' column shows the exchange-correlation functional used (LDA/GGA). All the band centers are relative to the Fermi energy; all the energies and widths are in eV. A further observation is that the magnetisation of the surface Pt atoms is generally reduced by both oxygen and hydrogen adsorption. This is shown qualitatively in Figure  15, where the surface average spin moment of the Pt atoms in the reacted surface is less than that on the clean surface.  Figure 14 in the upper 'reacted' surfaces, whilst the lower surface is unreacted. This asymmetry is simply intended to highlight the change in the magnetic character of the respective surface Pt layers during adsorption, and quantitative results from these simulations were not subsequently used in this work. The averaging was performed parallel to the (111) plane. The labels show the regions of bulk Pt3Ni and PtNi3, selvedge Ni and surface Pt.

Alloy Functional
The observation that the magnetisation of the surface Pt atoms is reduced by both oxygen and hydrogen adsorption is presented quantitatively in Table 8. Comparing the reacted and unreacted surfaces, the adsorption generally reduces the surface Pt atoms  Figure 14 in the upper 'reacted' surfaces, whilst the lower surface is unreacted. This asymmetry is simply intended to highlight the change in the magnetic character of the respective surface Pt layers during adsorption, and quantitative results from these simulations were not subsequently used in this work. The averaging was performed parallel to the (111) plane. The labels show the regions of bulk Pt 3 Ni and PtNi 3 , selvedge Ni and surface Pt.
The observation that the magnetisation of the surface Pt atoms is reduced by both oxygen and hydrogen adsorption is presented quantitatively in Table 8. Comparing the reacted and unreacted surfaces, the adsorption generally reduces the surface Pt atoms moment by up to 75%. This is not the case for on-top adsorption, however, where the magnetic moment is seen to increase.

Conclusions
The reacted (111)-Pt/Ni/Pt 3 Ni and (111)-Pt/Ni/PtNi 3 surfaces have been investigated using density functional theory (DFT). These investigations have focused on the binding energy, bond lengths, magnetisation and electronic density of states of adsorbates and the surface Pt atoms during the incorporation of oxygen and hydrogen into the high symmetry ad-sites.
The current work has shown that the surface Pt atoms in the unreacted (111)-Pt/Ni/Pt x Ni 1-x slabs are permanently magnetised. In comparison with appropriately strained clean Pt(111), this magnetisation has been seen to be due to a ligand effect between the surface Pt and the selvedge Ni atoms, instead of being seen as a geometric effect due to the strain.
The adsorption of oxygen and hydrogen into the high symmetry sites quenches this magnetisation. Comparing the adsorption energies of the reacted (111)-Pt/Ni/Pt x Ni 1-x slabs and the reacted, strained Pt(111) surface, it has been shown that this quenching of the magnetisation significantly increases the binding energy of the adsorbate. This effect has been shown to be predicted with both the LDA and the GGA.
A survey of the behaviour of the clean Pt(111) and bulk Pt systems has been included in the current work. This survey has quantified the shifting in the d-and sp-band centres and changes to the width of these bands under strain and has explained this behaviour using a simple bulk-to-gas phase model. The correlation between the adsorption energy of H/Pt(111) and O/Pt(111) under strain and the position of the centre of the d-band has been qualitatively discussed.
Funding: This work was supported by the author's membership to the UK's HEC Materials Chemistry Consortium, which is funded by the EPSRC (EP/L000202), and used the ARCHER UK National Supercomputing Service (http://www.archer.ac.uk).