Boosting the oxygen evolution activity in non-stoichiometric praseodymium ferrite-based perovskites by A site substitution for alkaline electrolyser anodes SustainableEnergy&Fuels

Sustainable fossil fuel free systems are crucial for tackling climate change in the global energy market, and the identi ﬁ cation and understanding of catalysts needed to build these systems plays a vital role in their development. ABO 3 (cid:1) d perovskite oxides have been observed to be potential replacement materials for the high-performing, but low ionic conducting and economically unfavourable Pt and IrO 2 water splitting catalysts. In this work increased addition of Sr 2+ aliovalent dopant ions into the crystal lattice of Pr 1 (cid:1) x Sr x FeO 3 (cid:1) d perovskites via A site substitution was seen to drastically improve the electrocatalytic activity of the oxygen evolution reaction (OER) in alkaline environments. The undoped PrFeO 3 (cid:1) d catalyst was not catalytically active up to 1.70 V against the reversible hydrogen electrode (RHE), whilst an onset potential of 1.62 V was observed for x ¼ 0.5. Increased strontium content in Pr 1 (cid:1) x Sr x FeO 3 (cid:1) d was found to cause a reduction in the lattice parameters and crystal volume whilst retaining the orthorhombic Pbnm space group throughout all dopant levels, analysed using the Rietveld method. However, it was noted that the orthorhombic distortion was reduced as more Sr 2+ replaced Pr 3+ . The mechanism for the increased electrocatalytic activity with increased strontium is due to the increasing concentration of oxygen vacancy ( d ), leading to increased catalyst site availability, and the increased average oxidation state of Fe cations, consistent with the iodometric titration results. This results in shifting the average d shell e g electron ﬁ lling further towards unity. X-ray photoelectron spectrum of the O 1s core level also shows the presence of lattice oxide and surface hydroxide/carbonate. This work shows promise in that using the more abundant and more economically friendly material of strontium allows for improved OER catalytic activity in otherwise inactive perovskite catalyst oxides.


Introduction
The increasing effects of climate change are fuelling a drive to develop environmentally friendly renewable energy technologies in order to facilitate a move away from fossil fuel dependency. [1][2][3][4][5] Hydrogen, apart from being the most abundant element known in the universe, is an ideal fuel source for renewable energy applications. Hydrogen is known to have one of the highest relative specic energy constants (120 MJ kg À1 ) when compared to fuels such as liqueed natural gas (54.4 MJ kg À1 ), propane (49.6 MJ kg À1 ), automotive gasoline (46.4 MJ kg À1 ), and ethanol (29.6 MJ kg À1 ). 6 Therefore, the development of technologies that facilitate clean hydrogen generation is an important step in clean energy development.
Water is one of the most abundant resources on the planet, is environmentally neutral, and is a good source of hydrogen. The electrochemical splitting of water into hydrogen and oxygen using electrolysers is a key technology in the development of clean energy, with the oxygen evolution reaction (OER) and the hydrogen evolution reaction (HER) taking place at the anode and cathode respectively. However, although water may be an abundant source of hydrogen, the splitting reaction is rate limited by the kinetically sluggish nature of the oxygen evolution reaction. [7][8][9] This sluggishness has led to energy efficiency rates of as low as 53%. [10][11][12][13] Hence, one of the keys to practical hydrogen production lies in overcoming the sluggishness of the OER.
Current state-of-the-art catalysts employed in the electrochemical splitting of water are the noble metal oxides IrO 2 and RuO 2 for the anodic half-cell reaction (eqn (1)), for the OER in alkaline solutions. 7,13 These materials are some of the highest performing anodes, requiring overpotentials of around 470 mV to achieve a current density of 10 mA cm À2 in 0.1 M KOH electrolyte. [13][14][15] These materials however, are in low abundance, demonstrate poor ionic conductivity of oxygen species leading to high internal resistance manifesting as a limitation to the reaction rate, and utilise high cost materials. It is for this reason that there is a drive to replace these materials with high performing and more economically sustainable materials.
Perovskite oxides with the formula ABO 3Àd (A ¼ rare earth metal and B ¼ transition metal), synthesised using various routes, [16][17][18][19][20] have been well researched as potential alkaline electrolyser anodes due to their relatively low cost, high mixed ion-electron conductivity and high surface oxygen exchange coefficient. 21,22 Perovskite oxides are suitable OER catalysts as the materials can overcome certain barriers when facilitating the complex oxidation reaction. The slow nature of the 4-electron oxidation process per mole of oxygen at the catalyst surface in alkaline media (eqn (1) 7,23,24 ) has been suggested by Rossmeisl et al. to be directly related to the binding strength of oxygen intermediates with the transition metal catalyst surface. More specically, catalysts that form strong bonds with oxygen are limited by HOO* formation, whereas catalysts which form weak bonds with oxygen are limited by the O* formation, where * indicates an active catalyst surface site. 25 This is represented as a volcano type relationship between oxygen-bonding strength and catalytic activity. Suntivich et al. expanded on this theory, discovering an activity descriptor for catalytic activity relating to the lling of the transition metal higher energy d shell e g s* antibonding orbital, with an e g lling of single unit occupancy giving the highest oxygen species to transition metal catalyst interaction in octahedral complexes. 7,26 The understanding of the lling of the transition metal d shell e g s* antibonding orbital, and the interaction with the 2p p bonding orbital of the oxygen species has led to developments of rst row transition metal based perovskite oxide catalysts such as Ba 0.5 Sr 0.5 Co 0.8 Fe 0.2 O 3Àd , [27][28][29] LaNiO 3Àd , 30-32 and LaFeO 3Àd 33 for use as electrolyser anodes. Perovskite oxide materials have even been shown to have a comparable OER performance with IrO 2 catalysts. [34][35][36] Perovskite oxides also offer considerable control over their crystal structure with relatively small variations to the synthesis process. The exact science behind how tuning the crystal structure of a perovskite to orchestrate higher catalytic activity is oen overlooked in the literature. However, a full understanding of how tuning simple perovskite oxides is critical in the design of more complex catalytic systems and could hold the key to unlocking feasible noble metal free OER catalysts. One such structural tuning is the introduction of oxygen vacancies within the lattice, thought to increase the OER catalytic activity. [37][38][39][40][41][42][43] Doping perovskite oxide A sites through the substitution of a trivalent rare earth (RE) metal by a divalent alkaline earth (AE) metal will cause a charge imbalance throughout the material. This imbalance is rectied by the introduction of oxygen vacancies and/or increasing the oxidation state of the transition metals in order to retain overall charge neutrality. In the case of strontium doped praseodymium ferrites (PSFO, Pr 1Àx Sr x FeO 3Àd ), increased strontium concentration (x) leads to increased oxygen vacancy concentration (d). 44 The proposed mechanism for oxygen vacancy (V) manufacture arising from substituting RE 3+ with AE 2+ denoted by the Kröger-Vink notation is shown in eqn (2).
where Â, 0 and c represent neutral, single negative and single positive charges respectively. The presence of a vacancy in the perovskite lattice is thought to increase the availability of the transition metal catalytic sites at the triple phase boundary, leading to increased OER performance. 45,46 Electrolyte selection is important in electrolyser systems. In order to lessen high overpotentials due to the electrolyte, the water splitting reaction is best carried out under strongly acidic or alkaline conditions due to the lower solution resistance provided via the higher availability of ions contained in these solutions. 47 In order to avoid pH gradients from forming due to the high pH, the electrodes are separated by an ion conducting membrane in electrolyser setups. 48 In acidic media, transition metal oxides are associated with energy loss due to corrosion and high overpotentials at the anode. 25,[49][50][51][52] Therefore, in order to extend the lifetime of the catalyst and to reduce overpotentials, alkaline environments are preferential when utilising transition metal based catalysts, which will allow for a replacement of noble metal based systems, increasing economic and environmental viability.
In this paper we investigate the effects of different doping levels of strontium into praseodymium ferrite oxide perovskite structures on the crystal structure, surface electronic states and the electrochemical performance, specically of the oxygen evolution reaction.

Material synthesis
Pr 1Àx Sr x FeO 3Àd (x ¼ 0, 0.1, 0.2, 0.3, 0.4, 0.5) powders were prepared by mixing stoichiometric amounts of Pr(NO 3 ) 3 $6H 2 O (99.9% Alfa Aesar), Sr(NO 3 ) 2 (99%, Alfa Aesar), and Fe(NO 3 ) 3 -$9H 2 O (min. 98%, Alfa Aesar) in deionised water to give an overall concentration of 0.02 M Pr 1Àx Sr x FeO 3Àd . Citric acid monohydrate (min. 99.5%, Alfa Aesar) was added as a chelating agent with a molar ratio of 1 : 10 metal to citrate complex. The solution was stirred overnight and then slowly heated until a gel-like substance had formed. This gel was then transferred to an alumina crucible before sintering in a muffle furnace at 800 C for 6 hours under static air, with a heating rate of 2 C min À1 .

Characterisation methods
The crystallinity of the powders was analysed using laboratory powder X-ray diffraction, using a PANalytical X'Pert Pro MPD, tted with an X'Celerator detector. The diffraction patterns were acquired by irradiating powder samples with Cu Ka X-ray radiation (l average ¼ 1.5418Å). The samples were scanned over a range of 5-120 in 2q with a step size of 0.0334 and nominal time-per-step of around 6 seconds. Fixed divergence and antiscatter slits of 1 2 and 1 respectively were used together with a beam mask of 10 mm. The resulting diffraction patterns were rened using the Rietveld method of least squares 53 using TOPAS-Academic V6. 54,55 The starting model used the Pbnm space group with the atomic positions of undoped PrFeO 3Àd (ICSD code 27274), 56 with the Sr dopant sharing the Pr crystallographic site. Fundamental parameter peak shape prole 57 and X-ray form factors of Pr 3+ , Sr 2+ , Fe 3+ and O 2À were used. For each data set, the Chebyshev polynomial and 1/x background, a zero error, unit cell parameters, scale factors, and crystal size/strain were rened sequentially. The sum of occupancies of Pr 3+ and Sr 2+ was constrained to unity. Isotropic temperature displacement factors for all atoms were rened independently. The cation positions were rened, followed by the oxygen positions.
The oxygen content of the perovskite oxides was calculated from the oxidation state of iron, which was determined by iodometric titration. $20 mg of perovskite powder was dissolved in 20 ml of hydrochloric acid (12 M) by continuous stirring in a nitrogen gas ushed ask. An excess of 0.1 M KI was added once the powder had fully dissolved. 0.005 M Na 2 S 2 O 3 solution was used to titrate against the iodine generated. A 1% starch solution, as an indicator, was added towards the end of titration.
The electrocatalytic activity of Pr 1Àx Sr x FeO 3Àd perovskite oxides for the oxygen evolution reaction was measured in a three-electrode cell setup. Perovskite inks were prepared by mixing 2 mg of perovskite catalyst, 0.4 ml deionised H 2 O, 0.1 ml propan-2-ol, and 60 ml of Naon 117 (5% w/v). Homogeneous mixing was assured by ultrasonicating the ink for at least 20 minutes. 3 ml of the ink was drop cast onto a 3 mm diameter gold electrode (99.95% purity, Alvatek). 10.7 mg of catalyst was loaded onto the electrode. Cyclic voltammetry (CV), electrical impedance spectroscopy (EIS), and chronoamperometry measurements were carried out in a glass cell using an Ag/AgCl leakless electrode and platinum mesh as reference and counter electrodes respectively. A solution of nitrogen saturated 0.1 M KOH (Alfa Aesar, 99.98%) in deionised water was used as the electrolyte. [58][59][60] The Ag/AgCl reference electrode was calibrated against the reversible hydrogen electrode (RHE) in the same electrolyte solution giving a conversion value of E RHE ¼ E Ag/AgCl + 0.98 V. Potential and currents were applied and measured respectively using a Gamry Interface 5000E potentiostat. Cyclic voltammetry experiments were carried out in a range of 1.30-1.80 V RHE , using a scan rate of 5 mV s À1 . 3 cycles were measured aer activating the material by running chronoamperometry at 1.70 V RHE for 30 minutes. To account for resistances arising from the electrolyte solution, and internal connections, ohmic drop (iR) compensation was applied to the CV data post testing using resistance data gathered from EIS measurements. All current densities were normalised through the area of the electrode (0.07 cm 2 ). EIS scans were recorded from 10 5 to 10 À1 Hz at 1.70 V RHE . The performance of the material as a catalyst at electrolyser anodes was tested using a 300 mm thick titanium bre felt gas diffusion layer (GDL, Bekaert S.A. Fibre Technologies) with a 78% porosity using 20 mm thick titanium bres used as a support for the PSFO material. The mesh was cut into an area of 2 cm 2 . The anode was prepared by spraying the catalyst material directly onto the titanium mesh. The spray consisted of 10 wt% polystyrene-b-poly (ethylene/butylene)-bpolystyrene base ionomer, made with polymer synthesis as outlined by Gupta et al., 61 with 2 mg cm À2 loading of the PSFO perovskite samples, with a small amount of 5 wt% polytetra-uoroethylene (PTFE) binder.
The cathode for the hydrogen evolution reaction (HER) was constructed using a non-wet proofed carbon GDL with a microporous layer (MPL, Freudenberg), cut into a 13 mm diameter circle. For the HER catalyst, 0.4 mg cm À2 20 wt% Pt on carbon was used. The catalyst ink contained 20 wt% Pt on carbon, 28 wt% ionomer, 5 wt% PTFE binder, and propan-2-ol. This was sprayed directly onto the carbon GDL/MPL paper.
The radiation graed anion exchange membrane (AEM) was synthesised also according to Gupta et al., 61 using a low-density polyethylene (LDPE) with vinylbenzyl chloride (VBC) as the gra monomer. 0.1 M NaOH was used as the electrolyte. EIS and CV scans were measured at 20 C, 40 C, 60 C, and 80 C.
X-ray photoelectron spectroscopy (XPS) analysis was performed using a Kratos Axis SUPRA XPS tted with a monochromated Al Ka X-ray source (1486.7 eV), a spherical sector analyser, 3 multichannel resistive plates, and 128 channel delay line detectors. Ultraviolet photoelectron spectroscopy (UPS) measurements were performed using a He(I) lamp (21.22 eV) and an emission current of 25 mA. All spectra were recorded at 150 W and a spot size of 700 mm Â 300 mm. Survey scans were recorded at a pass energy of 160 eV, and high-resolution scans recorded at a pass energy of 20 eV. UPS scans were recorded with a pass energy of 5 eV. Electronic charge neutralization was achieved using an electron ood gun. The parameters used for the lament current, charge balance, and lament bias were 0.27 A, 3.3 V, and 3.8 V respectively. All sample data were recorded at a pressure below 10 À9 Torr and room temperature. Data was analysed using CasaXPS v2.3.19PR1.0. Peaks were t with a Shirley background prior to component analysis. Lineshapes of LA (1.53243) were used to t components. Charge calibration was achieved using adventitious carbon at 284.8 eV.

Crystal structure
Powder X-ray diffraction patterns for all PSFO samples indicate pure perovskite phase with minimal secondary impurity phases within the detection limit of the XRD instrument (Fig. 1). The x ¼ 0.4 and x ¼ 0.5 material have a small amount (<3 wt%) of the Pr 2 FeO 4+d Ruddlesden-Popper phase present. The Rietveld ts of the XRD patterns together with the Bragg markers for both perovskite and Ruddlesden-Popper phases are collated in Fig. S1-S6. † As the concentration of the divalent strontium species (Sr 2+ ) increases, the peaks shi towards a higher 2q value which results in the contraction of the perovskite lattice parameters ( Fig. 2a and b) but the overall crystal structure  where the peak splits are less apparent. This indicates the reduction of the orthorhombic distortion, leading to a higher degree of crystal symmetry indicated by the ratio of b/a lattice parameters approaching unity with increased Sr 2+ (Fig. 2c). 44,63,64 This reduction of orthorhombic distortion leading closer to an ideal cubic perovskite increases the electrical conductivity of the material. 65 This is explained by a wider charge carrier bandwidth and overlap between the Fe 3d and O 2p orbitals, which is due to the changes in the buckling angle of the Fe-O-Fe bonds, and the bond lengths of the Fe-O environment. 66,67 The contraction of lattice parameters and crystal volume as strontium doping increases is shown in Fig. 2. Overall, the a lattice parameter stayed at 5.48 to 5.49Å, whereas the b and c lattice parameters decreased in line with increasing Sr 2+ content. The overall crystal volume of the perovskite decreased with increased addition of Sr 2+ . Despite the larger ionic size of Sr 2+ (1.31Å, in 9-coordination) 68 occupying the Pr 3+ site (1.179Å ), the decreasing trend of the unit cell volume can be explained by an increasing average oxidation state of Fe, obtained from the iodometric titration, and as a result, a decrease in the coordination number of Fe-O, from 6 to 5, due to a higher oxygen vacancy content. This is consistent with a slight shi of the X-ray photoelectron spectrum towards a lower binding energy of Fe 2p as x increases, which will be explained in the next section. Lattice parameters and crystal volume are shown in Table 1. The reduction of the crystal lattice volume as Sr 2+ doping concentration increases is consistent with what has been observed in similar perovskite oxide materials. 69,70 Bond valence sum (BVS) analysis 71 which gives an estimation of the bond strength in ionic compounds, was calculated taking into consideration the oxygen vacancy concentration determined by iodometric titration using eqn (3).
where s i is the bond valence sum, R 0 and B are previously determined parameters, 71 and R is the bond length. 72 For determining the BVS of the Fe cation, 1.76Å and 0.35 were used  (5) 1.966 (7) 1.958 (4) 1.970 (10) 3.19 (7) 0.515 (7) 1.962 (7) 1.976 (4) 1.955 (10) 3.19 (7) for R 0 and B respectively. The oxygen vacancy weighted bond valence sum of Fe cations and Fe-O bond lengths for each material are listed in Table 2. Interestingly, Fig. 3a and b indicate that for lower Sr 2+ concentration (x ¼ 0-0.2) it is the axial Fe-O bonds that appear to have a higher contribution to the reduction in the average overall bond length, whereas for Sr 2+ concentrations of 0.3-0.5 the equatorial bond lengths have a greater contribution. This observation is worthy of note, but the reasoning behind this change is beyond the scope of this investigation. In line with what is hypothesised in eqn (2), Fig. 3d shows that as Sr 2+ concentration increases, so does the concentration of oxygen vacancies. The reduction of the oxygen stoichiometry (3 À d) whilst retaining the Pbnm space group was thought to lead to higher active sites of the B site Fe catalyst surface throughout the crystal lattice to the oxygen species to be oxidised, leading to an increase in OER catalytic activity.
The presence of Fe 4+ was ruled out in this system as this was not observed in either the iodometric titration or the XPS measurements. The higher overall oxidation state will lead to a d shell degenerate e g orbital lling closer to unity, which as previously discussed 7,26 will lead to higher OER catalyst activity.
This combined with the increased availability of the transition metal species through vacancies in the oxygen site were hypothesised to be the reason behind the increase in the OER electrocatalytic performance of the higher Sr 2+ doped PSFO material. This increase in the oxidation state is in good agreement with the recent literature regarding studies of the perovskite oxide La 1Àx Sr x FeO 3Àd for OER catalysis, where Fe oxidation states as high as +4 have been observed in high Sr 2+ doping concentrations through X-ray absorption spectroscopy of the Fe L 2,3 edge. 73,74 This increase of the overall oxidation state will induce changes in the electronic properties of the material, with increasing oxidation states affecting the density of states (DOS) near the Fermi level, thus altering the binding strength between the Fe 3d and oxygen species 2p orbitals which as previously discussed is in line with the activity descriptor set out by Suntivich et al. 58,75

Surface electronic structure and chemistry
The binding energies (BE) obtained from XPS were corrected by calibrating the scans so that the dominant low binding energy C 1s peak was at 284.8 eV. High resolution scans of Pr 3d in Fig. 4a show a common lineshape characteristic of Pr 3+ . 76 Fig. 4a shows two distinct typical peaks at around 933 and 953 eV, which are attributed to Pr 3d 5/2 and Pr 3d 3/2 respectively. 77,78 The difference of $20 eV is consistent with the literature values for the spin orbital splitting of the 3d 5/2 and 3d 3/2 praseodymium orbitals. 76,79 Unfortunately, the determination between Pr +3 and +4 is challenging due to the complex nature of nal state effects in such systems. For example, it has been reported 76,79 that a satellite peak at 967 eV, arising from 3d4f 1 nal states, is exclusive of Pr 4+ from studies of PrO 2 , which is absent in the scans presented in Fig. 4a. This could indicate that the praseodymium is present in a +3 oxidation state as would be expected in ABO 3Àd perovskite oxides. However, there is overlap between this satellite peak and the oxygen O KLL Auger peak ($973 eV) which makes satellite peak determination difficult. 80 However, careful analysis of the Pr 3d 5/2 lineshape and high energy nal state effects afford an appropriate methodology for praseodymium oxidation state determination. Analysis of the 3d 5/2 photoemission reveals the presence of one photoemission peak (m) and one shake-down satellite (s) at lower binding energy arising from well screened 4f 3 nal states. The ratio between these two peaks varies between Pr 3+ and Pr 4+ , with Pr 2 O 3 reporting a higher ratio of m : s than PrO 2 . 76 Deconvolution of our experimental data reports a constant m : s ratio of 3.05 AE 0.05 a.u. which is consistent with the theoretically determined ratio for Pr 3+ . 81 As the value for x increases, praseodymium is replaced with strontium. As expected, intensities of the peaks found in photoelectron scans of the Sr 3d core orbital increase with increasing values of x (Fig. 4b). Peak models were t using an Sr 3d doublet separation of 1.76 eV and dening two species, SrO at 132 eV and SrCO 3 at 133.8 eV (Fig. 5). 82 It was determined that as the strontium content increased in the perovskite, the ratio between SrO and SrCO 3 also increased in favour of the oxide with the exception of x ¼ 0.4, which reported an anomalous quantity of carbonate in relation to the rest of the series. It is of note that no carbonate species was detected within detection limits of the XRD scans. This is indicative that the carbonate species is located mainly at the surface of the catalyst which will play a critical role in blocking the surface catalytic reaction. 83 It is reasoned that this carbonate interference is the primary factor in explaining the anomalous electrochemical result seen for the x ¼ 0.4 material, which is further elaborated in the next section. This performance was replicated in experimental repeats which suggests that the Pr 0.6 Sr 0.4 FeO 3Àd (x ¼ 0.4) is prone to carbonate formation at this specic stoichiometry. The mechanism behind this specic carbonate formation is unclear, especially as this phenomenon is not seen in the surrounding stoichiometries (x ¼ 0.3 and 0.5), which presents an opportunity for future study.
The presence of both oxide and carbonate environments is further conrmed from O 1s scans (Fig. 4c), where there are two distinct peaks. Contributions from the lattice oxide at 528 eV are accompanied by contributions at a higher binding energy of around 531 eV, which is attributed to a mixed contribution of hydroxide and/or carbonate environments, with the deconvoluted scans presented in Fig. S7. † Unfortunately, the binding energies of hydroxide and carbonate environments are very close together and therefore difficult to distinguish from each other from the O 1s spectra. 84 Fig. 4d shows Fe 2p spectra across the series. The lineshape is characteristic of Fe(III) and was t using a model developed via a combination of hematite (a-Fe 2 O 3 ) reference material and the work of Grosvenor et al. 85 From these ts we can see an increase in binding energy of the principle photoemission rising from 709.7 eV (x ¼ 0) to 710 eV (x ¼ 0.5), indicating an increased localised charge (>d + ) among the iron centres. The shi in the Fe 2p core orbital towards a higher binding energy  with increased Sr 2+ content is consistent with the iodometric titration results.
Ultraviolet photoelectron spectroscopy was used to determine the workfunction for each of the perovskite catalysts. The workfunction decreased as the concentration of Sr 2+ cations increased on the A site of the perovskite lattice ( Fig. S8a and b †). This decrease is related to the energy requirement to ionise the perovskite lattice, or to move electrons throughout the material. The increased electron movement within the crystal lattice will aid in the increased electrocatalytic performance, as discussed in the next section.

Electrochemistry
Cyclic voltammetry (CV) and electrochemical impedance spectroscopy (EIS) scan results are in line with what is suggested by the changes in crystal structure. Material with a higher Sr 2+ concentration demonstrated considerably higher electrochemical performance than that of unaltered PrFeO 3Àd . iR corrected CV plots for the PSFO series are presented in Fig. 6. The undoped PrFeO 3Àd catalyst is still OER inactive up to 1.70 V against the reversible hydrogen electrode (RHE).
The oxygen vacancies in the Sr 2+ doped PSFO catalysts play a crucial role in the increase in the oxygen evolution electrocatalytic activity. When a vacancy is introduced to the oxygen site of the crystal lattice, the coordination number of the Fe-O is reduced from 6 to 5. The 5 coordinated Fe-O unit will have a greater availability of the Fe 3d orbital with the 2p orbital of the adsorbing OH À species, creating an intermediate OH À -FeO 5 octahedral unit. The overlap between the Fe 3d and O 2p orbitals leads to increased catalytic activity as discussed in previous sections. 86 The concentration of lattice oxygen vacancies increases with Sr 2+ concentration (Fig. 3d), allowing for more of this OH À -FeO 5 octahedral unit to be present at one time in the catalyst, explaining the high performance seen in the x ¼ 0.5 catalyst. The improved OER performance of perovskite oxide catalysts through oxygen vacancy creation via A site doping is also observed in similar systems, which show comparable performance to the PSFO catalysts documented in this study. 35,73,74,87,88 Fig. 6 shows iR corrected CV scans where the maximum current density response comes from the x ¼ 0.5 sample. Table 3 shows the OER performance of similar rare earth single perovskite oxides in alkaline media, compared with the x ¼ 0.5 material from this study.
There is an anomalous result in the cyclic voltammogram for x ¼ 0.4 (Fig. 6). This can be explained by the high concentration of carbonate present in the sample as determined by X-ray catalyst active sites, whilst also reducing the hydrophobicity of the electrode backing layer which can lead to structural degradation and electrode ooding. 83 The onset potential for the OER, taken as the potential required to generate 1 mA cm À2 of the current response, reduced in line with Sr 2+ concentration, except for x ¼ 0.4. The onset potential for the doped samples ranges from 1.67 V RHE (overpotential, h ¼ 440 mV) for x ¼ 0.1 to 1.62 V RHE (h ¼ 390 mV) for x ¼ 0.5. The improved performance with regard to the onset potential is a great success of this material design as the undoped PrFeO 3Àd (x ¼ 0) material did not produce a current response of 1 mA cm À2 even at a potential of 1.70 V RHE in halfcell testing. Lowering onset potentials is a critical metric in designing materials for electrolyser anodes due to the cost factor and durability associated with long-term operation of electrolyser systems. Fig. S9 † illustrates chronoamperometry results for the PSFO catalyst series measured at 1.70 V RHE in 0.1 M KOH for 30 minutes. As with the other electrochemical testing, the Pr 0.5 -Sr 0.5 FeO 3Àd (x ¼ 0.5) material elicited the highest performance with regard to current density and showed good stability over the chronoamperometry testing conditions. The 'steps' seen in the chronoamperometry measurements are most likely the result of an insulating layer of oxygen gas bubbles forming on the electrode surface reducing the current density, before the bubbles are released and contact between the electrolyte and electrode is restored. From this testing it can be seen that Sr 2+ Fig. 9 Cyclic voltammograms using Pt/C as a cathode material, PSFO x ¼ 0, 0.2, and 0.5 as anode material, LDPE with VBA as a radiation grafted anion exchange membrane, and an electrolyte of 0.1 M NaOH. Cycles measured at temperatures of (a) 20 C, (b) 40 C, (c) 60 C, and (d) 80 C. incorporation into the perovskite structure does not seem to cause any negative effects relating to material stability. Fig. 7 highlights the huge difference in the electrical impedance response between the undoped sample with the doped ones, with the highest doped x ¼ 0.5 material having the lowest polarisation resistance, R p . These measurements were taken at 1.70 V RHE as this is in the active electrochemical region as shown in Fig. 6. All tested materials showed a semi-circle behaviour, with similar values for the total ohmic resistance, R s .
The undoped PrFeO 3Àd has a Tafel slope of $108 mV dec À1 , when compared to around 50 mV dec À1 for the doped catalyst perovskite oxides (Fig. 8). The 60 mV dec À1 shi is indicative of a possible change in the rate limiting reaction. The 4-electron OER reaction per mole of oxygen consists of 4 reactions (eqn (4)- (7)) on the transition metal surface (M), each of which can act as the rate determining step, Kinetic simulations undertaken by Shinagawa et al. (2015) determined that for electrocatalysis of the oxygen evolution reaction where the Tafel slope is around 120 mV dec À1 , as is the case for the x ¼ 0 material, the reaction is rate limited by the formation of the metal oxide bond (eqn (4)). The reductive tendency in the Tafel slope is seen when Sr 2+ is introduced to the unit cell tends towards the OER rate limiting step to be shied to the formation of the peroxide bond on the metal surface (eqn (6)), or the desorption of O 2 from the metal surface (eqn (7)). 93 The Tafel slope values seen in this study are consistent with those found in the literature for similar perovskite materials. 30,94-96 Fig. 9 shows forward cyclic voltammetry measurements for x ¼ 0, 0.2, and 0.5 at 20 C, 40 C, 60 C, and 80 C in the full cell conguration. As expected, the current response of the material increases with increased temperatures. Electrolyser tests conrm electrochemical results from half-cell testing (Fig. 6), in that the addition of Sr 2+ cations increases the electrocatalytic performance of the PSFO catalysts with regards to the oxygen evolution reaction. At an operating temperature of 80 C the overpotential (h) required to elicit a current response of 10 mA cm À2 is 442 mV, 353 mV, and 335 mV for x ¼ 0, 0.2, and 0.5 respectively; and for 150 mA cm À2 the overpotential required is 752 mV and 680 mV for x ¼ 0.2 and 0.5 respectively.

Conclusions
Praseodymium strontium ferrite perovskite oxides with the formula of Pr 1Àx Sr x FeO 3Àd (x ¼ 0, 0.1, 0.2, 0.3, 0.4, and 0.5) have been successfully synthesised by the wet chemistry citrate sol gel method, and have demonstrated increased OER electrocatalytic activity for higher values of x. The mechanism behind the increased activity due to the addition of Sr 2+ via A site substitution of Pr 3+ to the perovskite lattice is attributed to increased oxygen vacancy concentration allowing for a higher availability of the B site transition metal catalyst surface and a shi in the average oxidation state of the B site iron ion leading to a d shell e g electron orbital lling closer to the well-dened activity descriptor of unity, all of which allows for higher electron-ion conductivity and catalyst activity. This improvement is especially impressive as the undoped PrFeO 3Àd (x ¼ 0) material was demonstrated to be essentially inactive towards oxygen evolution at 1.70 V RHE in 0.1 M KOH electrolyte in 3 electrode setups, and requiring relatively high overpotentials and operating temperatures in order to start becoming active in electrolyser cells, compared to the performance of the A site doped material. This relatively small change to the stoichiometry of the material, whilst retaining the crystal structure of perovskites, demonstrates how the OER performance can be further enhanced with the implementation of more environmentally and economically attractive materials.
Author's contributions S. W. and S. S. P. conceived and designed this work. S. W. performed the XRD experiments and analysed the crystal structure of materials. S. W., G. G. and M. M. performed the electrical measurements and analysis of the results. M. A. I. performed the X-ray photoelectron spectroscopy and ultraviolet photoelectron spectroscopy measurements and analyses. All authors co-wrote and reviewed the manuscript.

Conflicts of interest
There are no conicts of interest to declare.