Carbon-Supported Trimetallic Catalysts (PdAuNi/C) for Borohydride Oxidation Reaction

The synthesis of palladium-based trimetallic catalysts via a facile and scalable synthesis procedure was shown to yield highly promising materials for borohydride-based fuel cells, which are attractive for use in compact environments. This, thereby, provides a route to more environmentally friendly energy storage and generation systems. Carbon-supported trimetallic catalysts were herein prepared by three different routes: using a NaBH4-ethylene glycol complex (PdAuNi/CSBEG), a NaBH4-2-propanol complex (PdAuNi/CSBIPA), and a three-step route (PdAuNi/C3-step). Notably, PdAuNi/CSBIPA yielded highly dispersed trimetallic alloy particles, as determined by XRD, EDX, ICP-OES, XPS, and TEM. The activity of the catalysts for borohydride oxidation reaction was assessed by cyclic voltammetry and RDE-based procedures, with results referenced to a Pd/C catalyst. A number of exchanged electrons close to eight was obtained for PdAuNi/C3-step and PdAuNi/CSBIPA (7.4 and 7.1, respectively), while the others, PdAuNi/CSBEG and Pd/CSBIPA, presented lower values, 2.8 and 1.2, respectively. A direct borohydride-peroxide fuel cell employing PdAuNi/CSBIPA catalyst in the anode attained a power density of 47.5 mW cm−2 at room temperature, while the elevation of temperature to 75 °C led to an approximately four-fold increase in power density to 175 mW cm−2. Trimetallic catalysts prepared via this synthesis route have significant potential for future development.


Introduction
Fuel cells (FCs) are electrochemical energy conversion devices with superior performance compared to conventional internal combustion engines, as they are not limited by the Carnot cycle [1]. Additionally, FCs are more environmentally friendly than conventional energy generation processes since they generally do not generate carbon dioxide (with few exceptions, such as direct alcohol/formic acid FCs) and do not rely on the consumption of fossil resources. Consequently, FC research and development has grown linearly with the worldwide energy demand and is projected to increase as FCs replace conventional heat engines in various applications [2]. However, many challenges need to be overcome if mass uptake of fuel cells is to be realized, including issues of safety, fuel source, and ensuring high efficiencies at low operating temperatures.
Low-temperature FCs targeted for use in portable applications and, particularly, for use in the transportation area have recently undergone rapid development [3]. They are BOR NaBH 4 Many authors have discussed the advantages, disadvantages, and state-of-the-art in DBFCs, including postulating the replacement of O 2 by hydrogen peroxide (H 2 O 2 ) [7][8][9][10]. The use of H 2 O 2 as an oxidant has become more popular in recent years, in particular as this results in an entirely liquid FC, thereby facilitating its use in applications where a compact design is required. When H 2 O 2 is used as the oxidant, the catholyte consists of an acidic solution of H 2 O 2 in order to avoid its decomposition into H 2 O and O 2 . Consequently, the DBFC in which the hydrogen peroxide reduction reaction (HPRR) takes place at the cathode and the BOR at the anode is known as the direct borohydride-peroxide fuel cell (DBPFC). It can generate electricity through Equations (1) and (4), with Equation (5) The majority of published research on DBFCs discusses the need for effective and low-cost anode and cathode materials in order to facilitate the commercialization of DBFCs. This is a result of the high operating costs of the DBFC, which is correlated with the performance, efficiency, and cost of the electrocatalysts [5,[8][9][10][11]. In fact, anodic catalysts for DBFCs must have high catalytic activity for BOR and be inert towards BH 4 − hydrolysis, a chemical side reaction where BH 4 − is consumed, generating H 2 and BO 2 − . In addition, DBFC anodes need to exhibit chemical stability in a basic solution. Common anodic catalysts for DBFCs include noble metal catalysts such as Au [12][13][14][15][16], Pt [17][18][19][20], and their corresponding alloy materials [21][22][23][24], as well as some nonprecious metal catalysts, e.g., Ni [25,26], Zn [27], and Co [28][29][30]. In particular, Au has been extensively studied, as it was originally thought to be capable of inhibiting borohydride hydrolysis; however, it exhibits rather slow kinetics for borohydride oxidation [31,32]. Studies involving Au, Pd, and Pt bulk electrodes for BOR revealed that surface poisoning by reaction intermediates, including BH ads or BH 3 , was significantly different for the three noble metals. Results have shown that this effect was more severe for Pt and also that the Au surface was the Nanomaterials 2021, 11, 1441 3 of 24 least affected by poisoning [33]. Despite this observation, Pt-based catalysts typically demonstrate faster electrode kinetics in BOR [24,34,35]; however, limited resources hinder their application as DBFC anodic catalysts. Pd is a promising alternative to Pt, presenting similar chemical properties and reactivity. Pd is considerably more abundant in the Earth's crust than Pt (ca. fifty times) and hence presents fewer challenges regarding resource scarcity [36][37][38]. Therefore, Pd has been tested in FCs both as a cocatalyst with Pt and as a Pt-free catalyst [37,[39][40][41]. Decreasing the required quantity of Pd by mixing it with non-noble metals in a well-designed catalytic structure capable of efficient borohydride oxidation can decrease the cost of membrane electrode assemblies (MEAs) and boost the commercialization of DBFCs. It was previously found that co-doping can improve the catalytic performance of noble metal catalysts [42][43][44]. For example, a Pd-Co/C nanocatalyst has been demonstrated as an electrocatalyst for formic acid oxidation and methanol oxidation [45,46].
Previous studies reported Pd nanocatalysts and Pd nanoalloys as anodic catalysts for BOR in DBFCs [47,48]. Furthermore, Simões et al. [49] analyzed how Pd, Au, and Pd-Au carbon-supported nanoparticles with different compositions performed during BOR in alkaline media. Pd x Au 1-x /C electrocatalyst exhibited an open circuit potential (OCP) lower than Pd and Au did. Merino-Jimenez et al. [50] evaluated a titanium-supported and carbon nanofiber-covered Pd-Ir alloy for BOR. The authors concluded that alloying Pd with Ir improves Pd electroactivity towards BOR. This enhancement was associated with the favoring of BOR (direct oxidation) and suppression of hydrogen evolution reaction (HER). Furthermore, Grimmer et al. [51] studied BOR at Pd/C by cyclic voltammetry, rotating disc electrode (RDE) measurements, chronoamperometry coupled to NMR spectrometry, and in situ Fourier transform infrared spectroscopy (FTIR). RDE results yielded n = 4 (at 0.4 V vs. reversible hydrogen electrode (RHE)) and n = 8 (at 0.8 V vs. RHE). Thus, at 0.8 V vs. RHE, BOR is 100% direct, while at 0.4 V, two different pathways are proposed to operate: (i) the total oxidation of BH 4 − to BO 2 − and two H 2 molecules and (ii) a four-electron electrooxidation with the production of BH 2 − species. The authors also suggested that at low potentials, Pd-BH 2 OH species block surface sites, while at higher potentials, this intermediate is fully oxidized (n = 8). Additionally, Pd has been shown to be relatively unharmed by surface poisoning, preserving good BOR kinetics even in highly concentrated BH 4 − solutions relevant for DBPFC operation [33]. Recently, by screening selected metal combinations (e.g., Pd with Co, Au, Ti) over a wide range of compositions, bimetallic and trimetallic Pd-based materials with optimum compositions that exhibited good activity for BOR were identified [51][52][53][54][55][56]. According to multiple published reports, Ni presents several benefits when added as a cocatalyst to C-supported Pd. First, it is an Earth-abundant metal, and therefore its use tackles issues of resource scarcity and reduces cost. Secondly, it is oxyphilic and capable of generating the oxygen species necessary for the redox reaction. Finally, the presence of Ni can exert an electronic effect over Pd, altering and potentially enhancing the catalytic properties of the latter [57][58][59].
In this work, investigations on the oxidation of BH 4 − in alkaline media at different working conditions are presented for a selection of trimetallic carbon-supported catalysts based on Pd (PdAuNi/C), prepared by a facile and scalable synthetic route. There is a consensus that adding a second metal as a cocatalyst to Pd is beneficial for its electrocatalytic performance. This is due to the electronic, geometric, and bifunctional effects that the second metal can provide. Previous works have considered multiplying the catalytic benefits by adding two metals to Pd instead of one [56,[60][61][62][63][64][65][66]. It is hypothesized that adding two metals, such as Au and Ni, each of which is known to be a beneficial addition to Pd individually [67][68][69], would significantly enhance Pd electrocatalytic performance. A small number of previous studies have previously investigated PdAuNi electrocatalysts [67][68][69][70][71][72]. Au [12,[73][74][75][76][77] and Ni [25,[78][79][80], like Pd, have been proven to possess good catalytic activity towards BOR, which makes their incorporation in a catalyst additionally attractive. Specifically, in the case of Ni electrodes, the formation of β-NiOOH has been pointed out as essential for the oxidation of BH 4 − [26]. Furthermore, to maximize the potential for BOR, a high-surface inert, mesoporous, and conductive support material is necessary to anchor the trimetallic nanoparticles and prevent their coalescence. Vulcan carbon (XC-72R) is commonly applied for this application due to its high surface area, mesoporous structure, and affordable cost. It was therefore chosen to support the PdAuNi and Pd nanoparticles. There are various chemical and physical synthetic routes in the literature to prepare dispersed metal nanoparticles [81]. The synthesis protocol is thought to play a crucial role and contribute to the catalyst physicochemical and catalytic characteristics. Any chosen method should balance between fine-tuned catalytic properties and ease of preparation, leading potentially to facile mass production [82]. NaBH 4 metal reduction is a very common means to prepare C-supported catalysts due to its easy procedure and quality-produced catalysts. Moreover, there is more than one NaBH 4 protocol to prepare Pd-supported nanoparticles. This work aims to synthesize three carbon-based trimetallic catalysts based on different borohydride reduction protocols. The physicochemical and morphological properties of those materials were determined and compared. Their performance as electrocatalysts for BOR was analyzed and discussed considering the different synthesis routes employed and the final properties obtained for each material.

Catalyst Preparation
Three borohydride reduction protocols using (i) a NaBH 4 -ethylene glycol complex, (ii) a NaBH 4 -2-propanol complex, and (iii) a three-step route were deployed to synthesize carbon-supported PdAuNi catalysts denoted herein as PdAuNi/C SBEG , PdAuNi/C SBIPA , and PdAuNi/C 3-step , respectively. The theoretical metal loading based on the quantity of metal salts used in the synthesis was 12 wt.% for PdAuNi/C SBIPA and 20 wt.% for both PdAuNi/C SBEG and PdAuNi/C 3-step . Although electrocatalysts for laboratory studies usually load the carbon support with 20 to 40 wt.% of metal, the use of lower loadings is more representative of industrial catalyst loadings which are typically~5 wt.% [83].
The SBEG synthesis protocol uses a NaBH 4 -ethylene glycol complex, following previous work [36,84,85]. Vulcan carbon (202 mg) was first dispersed in ethylene glycol (10 mL). The metallic precursor solution containing PdCl 2 (28 mg mL −1 ), AuCl 3 (55 mg mL −1 ), and NiCl 2 (36.7 mg mL −1 ) was then suspended in ethylene glycol (10 mL). The two suspensions were subsequently mixed. Next, a solution of NaBH 4 (20 mL) was slowly added to the metal and carbon mixture under magnetic stirring. The mixture was heated to 35 • C and kept at that temperature for 3.5 h. The mixture was then left to cool and washed copiously with deionized water until neutralization. Following this, it was dried in a vacuum oven at 80 • C for 2 h.
The second reduction protocol (SBIPA) applied NaBH 4 following methods described in previous reports [86][87][88][89]. The Vulcan carbon and metal precursors were sonicated in a mixture of 2-propanol/water (50 vol.%). KBr was added to stabilize the metal nanoparticles. The mechanism of this is that the larger Br − ions substitute for Cl − ions in solution, stabilizing Pd particles and preventing further particle agglomeration [90]. Unlike an organic surfactant such as PVP, KBr can be washed off easily after synthesis from the particle surface. The KBr/metal atomic ratio was 1.5. The mixture was stirred for 10 min, followed by the addition of NaBH 4 solution (0.5 M, 15 mL). Stirring then continued for another 30 min. Following this, the sample was washed under vacuum filtration. Finally, the wet powder was dried at 80 • C in a vacuum oven overnight.
The third protocol to prepare PdAuNi/C (three-step) aimed to compare the stepwise reduction of metals on the carbon surface to the co-reduction procedure practiced in the first two protocols. Consecutive reduction of different metals can lead to the formation of core@shell structures, such as Ru@Pd/C [91], Au@Pd [92,93], and Cu@Pd [94]. The NaBH 4 -2-propanol reduction mixture was once more used to prepare this catalyst. The procedure started with sonicating Vulcan carbon (202.4 mg) in a mixture of 2-propanol/water (50 vol.%). Then, NiCl 2 was added and kept in the solution under stirring for 1 h. This was longer than the reduction time required for Au and Pd reduction time, as Ni is more oxyphilic and less reducible [95,96]. Subsequently, AuCl 3 was added with stirring for 30 min. Finally, PdCl 2 was added with stirring for 30 min. To illustrate the impact of adding both Au and Ni onto Pd, a monometallic C-supported Pd catalyst (Pd/C SBIPA ) was synthesized using the SBIPA procedure.

Physical Characterization
The crystal structure of the prepared catalysts was investigated using X-ray diffraction (XRD) with a Bruker D2 Phaser operating using Cu K α radiation at 30 kV, 10 mA, and 12 • min −1 scan rate. To evaluate the crystallite size, the Scherrer equation was applied using the Pd (111) peak details. To quantify the metallic species in each catalyst powder, inductively coupled plasma optical emission spectrometry (ICP-OES) was applied using a Spectro Ciros Vision spectrometer by Spectro Analytical Instruments Inc. (Kleve, Germany). The samples were first digested in HNO 3 (2%), followed by 10-fold dilution in aquaregia (8%). The surface-weighted chemical composition of the catalysts was examined through energy-dispersive X-ray spectroscopy (EDX) attached to a JEOL JSM 6010LA scanning electron microscope (SEM) by JEOL Ltd. (Tokyo, Japan). Each catalyst surface was examined twice applying two different accelerating voltages: 10 kV and 20 kV [97] enabling studying the composition at two different depths. The catalyst morphology was examined by transmission electron microscopy (TEM) employing a Phillips C100 microscope (Hillsboro, OR, USA) operating at 100 kV applying LaB 6 filament. The TEM samples were prepared by suspending 5 mg of each sample powder in 2 mL of ethanol. Then, the mixture was sonicated for 1 h. A spray atomizer was then used to take 10 µL of the suspension to paint the holey carbon-coated Cu grid. The grid was left to dry for 24 h. The surface chemistry of the prepared electrocatalysts was studied by X-ray photoelectron spectroscopy (XPS) employing a Thermo Scientific K-Alpha+ spectrometer (Waltham, MA, USA). An Al-X ray source (72 W) and 400-µm 2 were applied. The pass energy to record the data was 150 eV for survey scans and 40 eV for high-resolution ones. The survey scan step size was 1 eV, but that of the high-resolution was 0.1 eV. Low-energy electrons and argon ions were used to neutralize the charge. CasaXPS from Casa Software Ltd. (Teignmouth, UK) was used to analyze the data and employed a Shirley-type background and Scofield cross-sections with an energy dependence of −0.6.

Electrochemical Evaluation
Electrochemical characterization was carried out in a conventional three-electrode setup of 125 mL. An ALS 2325 bipotentiostat (ALS Co., Ltd, Japan) combined with an RRDE-3A apparatus was used for fundamental experiments. Pt served as a counter electrode, and a saturated calomel electrode (SCE) as a reference electrode. All potential values within this paper were converted and are presented relative to RHE.
BOR measurements were carried out using cyclic and linear scan voltammetry. A 2 M NaOH solution was employed as the supporting electrolyte for all fundamental experiments. Cyclic voltammetry (CV) measurements were carried out in the potential range from OCP to 1.2 V applying scan rates from 5 to 1000 mV s −1 . Furthermore, the temperature effect was analyzed in the 25-65 • C temperature range. The BH 4 − concentration influence was studied in the 0.01-0.12 M range. Studies on the effect of temperature and concentration were conducted at a scan rate of 50 mV s −1 . RDE studies were carried out by linear scan voltammetry (LSV) at 10 mV s −1 by applying different rotation rates (0-2400 rpm). Apart from the concentration study, all experiments were conducted in a 0.03 M NaBH 4 in 2 M NaOH solution. Experiments were run at 25 • C and 0 rpm unless otherwise stated.

Fuel Cell Testing
A laboratory-scale direct borohydride peroxide fuel cell (DBPFC, V = 100 mL for both anodic and cathodic compartments) was assembled and tested. The two compartments were separated by a Nafion ® 117 (DuPont, Wilmington, DE, USA) cation-exchange membrane. A PAR 273A (Princeton Applied Research, Inc., Oak Ridge, TN, USA) potentiostat using PowerSuite software (PowerCV+ PowerSTEP) was employed for control of the experiments. The anolyte and catholyte were 1 M NaBH 4 in 4 M NaOH and 5 M H 2 O 2 in 1.5 M HCl solutions, respectively. The chosen compositions are proven to be optimal in previous studies [99,100]. PdAuNi/C SBIPA was employed as the anodic catalyst (1 cm 2 , 0.100 mg cm −2 catalyst loading, 0.012 mg cm −2 metal loading). The catalytic ink used for the anode was prepared as described in Section 2.4, with a GC tip used as substrate electrode. The cathode was a Pt mesh with a high surface area (Johnson Matthey, A = 50 cm 2 ) to ensure that fuel cell performance was not limited by the cathodic process. Interelectrode distance was the lowest possible to minimize ohmic resistance. DBPFC performance was analyzed in the 25-65 • C temperature range. Cell polarization curves were recorded, and the corresponding power density curves were plotted. Figure 1 shows the XRD patterns of Pd/C and the three trimetallic catalysts. All catalysts showed one broad peak at 25 • , typical of graphitic Vulcan carbon. For the monometallic Pd/C, the pattern showed four peaks that are characteristic of Pd (111), (200), (220), and (311) faces at 39.8 • , 46.5 • , 68.3 • , and 82.1 • , respectively (PDF#46-1043). However, for both trimetallic catalysts prepared by co-reduction (SBEG and SBIPA), there was a small shift in the peak positions to lower diffraction angles. This shift is ascribed to alloying among the three metals. Pd and Au, in particular, are well-known to exhibit remarkable alloying behavior [95,[101][102][103][104][105][106][107][108]. Although Pd and Ni are known to form a solid solution for their bulk alloys, their potential for formulating a nanoalloy is lower [96,[107][108][109]. It is not uncommon to see distinctive Ni(OH) 2 peaks at 60.1 • , which is indicative of bulk Ni species segregation. In the present study, the shift in the two co-reduced trimetallic samples was suggestive of a nanoalloy formation among the metals. However, the appearance of Ni(OH) 2 peaks at 33.4 • and 60.1 • in the case of PdAuNi/C SBEG implied there was some Ni segregation in the bulk of that catalyst. Additionally, these Ni(OH) 2 peaks were also present in the three-step catalyst, although smaller than that of PdAuNi/C SBEG . not uncommon to see distinctive Ni(OH)2 peaks at 60.1°, which is indicative of bulk Ni species segregation. In the present study, the shift in the two co-reduced trimetallic samples was suggestive of a nanoalloy formation among the metals. However, the appearance of Ni(OH)2 peaks at 33.4° and 60.1° in the case of PdAuNi/CSBEG implied there was some Ni segregation in the bulk of that catalyst. Additionally, these Ni(OH)2 peaks were also present in the three-step catalyst, although smaller than that of PdAuNi/CSBEG. In contrast to PdAuN/C3-step and PdAuNi/CSBEG, the catalyst synthesized via the SBIPA protocol did not show any individual Ni peaks. This indicated a higher potential of mixing of the three metals when the SBIPA protocol was used. KBr addition may also have contributed due to the replacement of the smaller Cl − anion with Br − during synthesis [110].

Physical Characterization Results
Unlike the SBEG and SBIPA catalysts, the XRD peaks of PdAuNi/C3-step were closer to pure Au (PDF#04-0784) than pure Pd (PDF#46-1043). This suggested the formation of Au@Pd (core@shell) structure as a consequence of the stepwise reduction process as previously reported [92,93]. Another important observation is that the XRD peak shift to lower diffraction angles was suggestive of the expansion of the Pd crystal lattice by the incorporation of the larger Au atoms into the Pd lattice. Ni is, theoretically, expected to contract the lattice, not to expand it. However, previous studies have found that adding Ni to Pd does not seem to displace the Pd diffraction peaks and therefore Ni does not necessarily contract the Pd lattice [59,111,112]. The peak shift in the case of PdAuNi/CSBEG and PdAuNi/CSBIPA is therefore consistent with lattice expansion and tensile strain that could upshift the d-band center and decrease the adsorbate-Pd bond strength [105].
The XRD Pd(111) peak of the four catalysts was analyzed to extract the peak broadening and diffraction angle values of each sample, and subsequently to determine their structural parameters: particle size, interplanar distance, and lattice constants ( Table 1). The detailed analysis is explained in section "1. XRD calculation" of supplementary information (SI). While the interplanar distance and lattice constant were the same for Pd/CSBIPA, PdAuNi/CSBIPA, PdAuNi/CSBEG (2.26 and 3.92 Å, respectively), it has increased in the case of PdAuNi/C3-step to 2.33 and 4.03 Å, respectively. It is noteworthy that the smallest particle size and highest broadening (FWHM) were those of PdAuNi/CSBIPA, followed by the monometallic Pd/CSBIPA, PdAuNi/CSBEG, and, finally, PdAuNi/C3-step. Factors other than crystallite size can also contribute to peak broadening, e.g., the instrument employed and the lattice strain. However, the Scherrer equation (Equation (S3)) is known to provide a reliable approximation of the crystallite size. The large particle size of PdAuNi/C3-step could be attributed to its Au-like structure, which presents larger lattice and interplanar distances. In contrast to PdAuN/C 3-step and PdAuNi/C SBEG , the catalyst synthesized via the SBIPA protocol did not show any individual Ni peaks. This indicated a higher potential of mixing of the three metals when the SBIPA protocol was used. KBr addition may also have contributed due to the replacement of the smaller Cl − anion with Br − during synthesis [110].
Unlike the SBEG and SBIPA catalysts, the XRD peaks of PdAuNi/C 3-step were closer to pure Au (PDF#04-0784) than pure Pd (PDF#46-1043). This suggested the formation of Au@Pd (core@shell) structure as a consequence of the stepwise reduction process as previously reported [92,93]. Another important observation is that the XRD peak shift to lower diffraction angles was suggestive of the expansion of the Pd crystal lattice by the incorporation of the larger Au atoms into the Pd lattice. Ni is, theoretically, expected to contract the lattice, not to expand it. However, previous studies have found that adding Ni to Pd does not seem to displace the Pd diffraction peaks and therefore Ni does not necessarily contract the Pd lattice [59,111,112]. The peak shift in the case of PdAuNi/C SBEG and PdAuNi/C SBIPA is therefore consistent with lattice expansion and tensile strain that could upshift the d-band center and decrease the adsorbate-Pd bond strength [105].
The XRD Pd(111) peak of the four catalysts was analyzed to extract the peak broadening and diffraction angle values of each sample, and subsequently to determine their structural parameters: particle size, interplanar distance, and lattice constants ( Table 1). The detailed analysis is explained in section "1. XRD calculation" of supplementary information (SI). While the interplanar distance and lattice constant were the same for Pd/C SBIPA , PdAuNi/C SBIPA , PdAuNi/C SBEG (2.26 and 3.92 Å, respectively), it has increased in the case of PdAuNi/C 3-step to 2.33 and 4.03 Å, respectively. It is noteworthy that the smallest particle size and highest broadening (FWHM) were those of PdAuNi/C SBIPA , followed by the monometallic Pd/C SBIPA , PdAuNi/C SBEG , and, finally, PdAuNi/C 3-step . Factors other than crystallite size can also contribute to peak broadening, e.g., the instrument employed and the lattice strain. However, the Scherrer equation (Equation (S3)) is known to provide a reliable approximation of the crystallite size. The large particle size of PdAuNi/C 3-step could be attributed to its Au-like structure, which presents larger lattice and interplanar distances. The bulk elemental composition of the synthesized catalysts was investigated by ICP-OES, as illustrated in Table 2. For Pd/C SBIPA , the identified Pd loading was 12 wt.%, equivalent to the theoretical value. For PdAuNi/C SBIPA , PdAuNi/C SBEG , and PdAuNi/C 3-step , Pd loadings of 9 wt.%, 16 wt.%, and 22 wt.% were identified, respectively. While ICP-OES determines the whole composition of the catalysts, the catalytically active sites are not dispersed throughout the bulk of the material. Instead, interaction with reactants takes place at surface sites. Therefore, the surface-weighted technique of EDX was used to gain insights into the composition of the catalysts in the catalytically relevant surface region. Figure 2 shows EDX metal elemental maps of PdAuNi/C SBEG . While the Pd (A) and Au (B) maps show similar and relatively even distribution of both metals, the Ni map (C) shows three areas that were more densely populated with Ni species, indicating that Ni was not homogeneously distributed at the surface. The regions of high Ni density correspond to regions of low density in the C map (D) and in the O map (E). This may suggest that some bulk segregation of Ni(OH)2 had occurred, as supported by XRD data where a Ni(OH)2 peak was observed, Figure 1. Table 3 shows the quantitative composition data extracted from EDX in both weight and atomic percent for the three trimetallic catalysts investigated. This aimed to highlight (E). This may suggest that some bulk segregation of Ni(OH) 2 had occurred, as supported by XRD data where a Ni(OH) 2 peak was observed, Figure 1. Table 3 shows the quantitative composition data extracted from EDX in both weight and atomic percent for the three trimetallic catalysts investigated. This aimed to highlight the different surface metal segregation by applying two accelerating voltages: 10 kV and 20 kV, the latter providing greater surface penetration. Unlike Pd and Au, the atomic concentration of Ni at 10 kV was larger than that at 20 kV. This further confirms that Ni segregation to the catalyst surface occurred. EDX analyses revealed higher metal loading than either ICP-OES or the theoretical compositions, as the synthesis resulted in a heterogeneous composition throughout the carbon support, with metal preferentially located towards the surface, beneficial for catalysis.  Figure 3 shows the EDX maps of PdAuNi/C SBIPA . Unlike PdAuNi/C SBEG , the Au (A), Pd (B), and Ni (C) maps show that the three metals were co-located with high-density regions of one metal corresponding to high-density regions of the others. Areas of low metal density corresponded to density regions in the C map (D). This is anticipated, as at regions of low metal density, the bare C support was exposed at the surface. Although to a lesser extent than PdAuNi/C SBEG , the quantitative data for PdAuNi/C SBIPA (Table 3) indicated a slightly higher concentration of Ni at 10 kV (2.90 at.%) as compared to that at 20 kV (2.52 at.%), perhaps indicating some surface segregation of Ni for this material as well. Overall, PdAuNi/C SBIPA showed relatively homogeneous metal distribution throughout the surface region, indicative of good mixing characteristics of the three metals in this case. This was further supported by the absence of any distinctive Ni peaks in this catalyst XRD pattern. Similar to PdAuNi/C SBIPA , the EDX analysis of PdAuNi/C 3-step revealed that Pd, Au, Ni, and C were co-located with each other ( Figure S1), suggesting good mixing and little Ni segregation in this material. Figure 4 shows representative TEM micrographs obtained for the Pd/C SBIPA , PdAuNi/C SBIPA , PdAuNi/C SBEG , and PdAuNi/C 3-step electrocatalysts, with additional TEM images being included in SI (Figures S2-S5). Examining the micrographs, the presence of highly dispersed metal nanoparticles on the carbon surface could be inferred. The very small particles indicated metal nanoparticles, while the large grey particles (40-60 nm) indicated the carbon aggregates. However, some particle aggregation and agglomeration could also be seen. This could be due to the synthesis method and/or being a monometallic system, which can promote particle growth and agglomeration more than a bimetallic or trimetallic one [56]. Comparing the trimetallic PdAuNi/C SBIPA catalyst with the monometallic one prepared by applying the same protocol suggests that the latter explanation may apply in this case. kV (2.52 at.%), perhaps indicating some surface segregation of Ni for this material as well. Overall, PdAuNi/CSBIPA showed relatively homogeneous metal distribution throughout the surface region, indicative of good mixing characteristics of the three metals in this case. This was further supported by the absence of any distinctive Ni peaks in this catalyst XRD pattern. Similar to PdAuNi/CSBIPA, the EDX analysis of PdAuNi/C3-step revealed that Pd, Au, Ni, and C were co-located with each other (Figure S1), suggesting good mixing and little Ni segregation in this material.   highly dispersed metal nanoparticles on the carbon surface could be inferred. The very small particles indicated metal nanoparticles, while the large grey particles (40-60 nm) indicated the carbon aggregates. However, some particle aggregation and agglomeration could also be seen. This could be due to the synthesis method and/or being a monometallic system, which can promote particle growth and agglomeration more than a bimetallic or trimetallic one [56]. Comparing the trimetallic PdAuNi/CSBIPA catalyst with the monometallic one prepared by applying the same protocol suggests that the latter explanation may apply in this case. PdAuNi/CSBIPA particles (1-6 nm, consistent with XRD analysis, Table 1) were very well dispersed, and less particle aggregation could be seen as compared to Pd/C. The reverse situation was noted for the PdAuNi/CSBEG, which showed significant particle aggregation. Above those few carbon particles, high metal particle aggregation could be seen. This could be promoted by the NaBH4-ethylene glycol reduction complex at 40 °C. Additionally, as discussed above with reference to the EDX map ( Figure 3C), the high Ni content in this catalyst might contribute to this agglomeration, noting that Ni exists only in the form of Ni(OH)2 ( Figure 5). The preparation method may also play a role in the promotion of agglomeration. For PdAuNi/CSBEG, the metal precursors were initially mixed in ethylene glycol, followed by the slow addition of NaBH4. Once the NaBH4 addition commenced, the orange/brown solution turned black instantly, alongside the formation of bubbles. PdAuNi particles had PdAuNi/C SBIPA particles (1-6 nm, consistent with XRD analysis, Table 1) were very well dispersed, and less particle aggregation could be seen as compared to Pd/C. The reverse situation was noted for the PdAuNi/C SBEG , which showed significant particle aggregation. Above those few carbon particles, high metal particle aggregation could be seen. This could be promoted by the NaBH 4 -ethylene glycol reduction complex at 40 • C. Additionally, as discussed above with reference to the EDX map ( Figure 3C), the high Ni content in this catalyst might contribute to this agglomeration, noting that Ni exists only in the form of Ni(OH) 2 ( Figure 5). PdAuNi/CSBEG has furthermore achieved the highest shift of Pd 3d to the left, which was also the broadest from the four-sample Pd 3d peaks as reported for the addition of Ni to Pd [108]. Examining the Ni 2p peaks of the three trimetallic samples, two high-bindingenergy satellite peaks were adjoined to the main peaks of Ni, which indicates Ni multielectron excitations [59,96,111,112]. Additionally, in the case of PdAuNi/CSBEG and PdAuNi/C3-step samples, the single pattern of two Ni 2p peaks suggests the presence of only one oxide species, Ni(OH)2, in contrast to the findings of Ramulifho et al. [113]. In terms of co-reduction, Au is the easiest of the three metals to reduce, followed by Pd. On the contrary, Ni is the slowest to reduce [95,114]. Therefore, it may not be incorrect to assume that the surface of the trimetallic catalysts (especially PdAuNi/CSBEG) was richest in Ni and poor in Au, which was evident from the XPS and EDX measurements. On the other hand, the core would be highly populated with Au and lacking a substantial Ni proportion.

Borohydride Oxidation Reaction Studies
The catalytic performance of the three PdAuNi/C electrocatalysts towards the BOR was initially investigated by scanning comparative CVs in 2 M NaOH with 0.03 M NaBH4 and in pure 2 M NaOH. Along with evaluating the activity of the synthesized trimetallic catalysts, the response of monometallic Pd/C was also evaluated in order to better understand the role played by Au and Ni. The CVs of three catalysts as well as of monometallic Pd/C depicted in Figure 6A-C show the absence of any relevant anodic peaks in the BOR potential region in NaOH solution, implying that the current generated during CVs recorded in the presence of BH4 − originated from its oxidation. Based on the shape of the CVs and similar OCP values of ca. −0.12 V, it could be assumed that the studied catalysts promoted similar reaction mechanisms. Similar OCP values were reported in other recent studies on BOR at Pd-containing materials [53,54]. Within the scanned potential range, CVs of three different catalysts showed two oxidation peaks in the positive scan, while an additional sharp oxidation peak was evident in the reverse scan. The first oxidation peak appeared as a wide oxidation hump at ca. 0.29 V for the PdAuNi/CSBEG and PdAuNi/CSBIPA catalysts, while in the case of PdAuNi/C3-step, this peak was more pronounced and positioned at more positive potential (0.38 V). It is believed that this peak can be attributed to H2 oxidation, which originates from the borohydride hydrolysis reaction, and to the initial The preparation method may also play a role in the promotion of agglomeration. For PdAuNi/C SBEG , the metal precursors were initially mixed in ethylene glycol, followed by the slow addition of NaBH 4 . Once the NaBH 4 addition commenced, the orange/brown solution turned black instantly, alongside the formation of bubbles. PdAuNi particles had most likely been reduced and segregated before the carbon support was added. For PdAuNi/C 3-step , on the other hand, Ni was initially reduced on the carbon surface (for 30 min), followed by Au (for 30 min), and finally Pd (for 30 min). As is clear from the XRD and EDX measurements, Ni segregation occurred once again, further verified by the observation of Ni(OH) 2 peaks, albeit in lower intensity than for PdAuNi/C SBEG . Additionally, Au had most likely formed a core@shell structure with Pd, since the XRD pattern was more representative of Au than of Pd. Although to a lesser degree than PdAuNi/C SBEG , significant particle aggregation was observed, and particle size was the largest in comparison to the other catalysts. Following PdAuNi/C 3-step , PdAuNi/C SBEG particle size was larger than that of both Pd/C SBIPA and PdAuNi/C SBIPA . The PdAuNi particle sizes obtained were slightly smaller than previously reported for trimetallic systems [68,69]. Figure 5 shows the XPS measured data for the Pd 3d peak positions of the four catalysts. The molar Pd:Au:Ni ratio of PdAuNi/C SBIPA equaled 64:17:19, while it was 5:1:94 and 58:6:36 for PdAuNi/C SBEG and PdAuNi/C 3-step , respectively. This showed the abundant Pd presence on the surface of PdAuNi/C 3-step and PdAuNi/C SBIPA samples, while PdAuNi/C SBEG contained mainly Ni (94 at.%) and only 6 at.% Pd and Au combined. Furthermore, it is noteworthy that the highest oxygen content (5.19 at.%) existed in PdAuNi/C SBEG , likely a consequence of the presence of significant Ni(OH) 2 (2.19 at.%). The Pd 3d peak of Pd/C SBIPA (Table S1 in SI) could be deconvoluted into a high-energy band corresponding to Pd 3d 3/2 (340.73 eV) and a low-energy band corresponding to Pd 3d 5/2 (335.43 eV). Those bands were shifted approximately 0.40 eV to higher binding energy compared to pure Pd, due to the interaction with the C support. Moreover, the Pd 3d binding energies of PdAuNi/C SBEG and PdAuNi/C SBIPA samples were shifted~0.5 eV and 0.15 eV, respectively, to lower binding energies. On the contrary, Pd 3d of PdAuNi/C 3-step was shifted 0.05 eV to binding energy higher than that of the Pd/C SBIPA sample. Pd 3d peak positions are shown in Figure 5. Ni 2p and Au 4f were also shifted in the three trimetallic samples compared to pure Au and Ni, the binding energies of which were 84 eV and 855.6 eV (Ni(OH) 2 ), respectively. While the Ni peaks shifted towards higher binding energy compared to pure Ni, the reverse trend was observed for Au. The Ni 2p shift was estimated to be +0.3 eV, +0.7 eV, and +0.8 eV for PdAuNi/C SBIPA , PdAuNi/C SBEG , and PdAuNi/C 3-step , respectively. The Au 4f shift was −0.3 eV, −0.4 eV, and −0.2 eV for PdAuNi/C SBIPA , PdAuNi/C SBEG , and PdAuNi/C 3-step , respectively. The shift in the binding energy was indicative of a change in the adsorption behavior of the respective catalysts. The presence of Au and Ni potentially explains why PdO was not detected in the trimetallic sample, unlike in the monometallic one, which contained 20 wt.% of Pd as PdO. Adding both Au and Ni enhanced Pd air stability. The full XPS surveys are illustrated in Figure S6. The individual elemental (Pd 3d, Au 4f, and Ni 2p) XPS peaks are available in Table S1, which confirms that Pd and Au existed in the metallic state in the trimetallic samples, Pd existed in metallic and oxide form in the monometallic sample, and Ni predominantly existed as Ni(OH) 2 . The Ni 2p peak in the PdAuNi/C SBIPA sample contained some Ni metal alongside Ni(OH) 2 , unlike PdAuNi/C SBEG and PdAuNi/C 3-step . This further exemplifies the good mixing of the three metals obtained by the SBIPA protocol, which was corroborated by the XRD and EDX data.
PdAuNi/C SBEG has furthermore achieved the highest shift of Pd 3d to the left, which was also the broadest from the four-sample Pd 3d peaks as reported for the addition of Ni to Pd [108]. Examining the Ni 2p peaks of the three trimetallic samples, two highbinding-energy satellite peaks were adjoined to the main peaks of Ni, which indicates Ni multi-electron excitations [59,96,111,112]. Additionally, in the case of PdAuNi/C SBEG and PdAuNi/C 3-step samples, the single pattern of two Ni 2p peaks suggests the presence of only one oxide species, Ni(OH) 2 , in contrast to the findings of Ramulifho et al. [113]. In terms of co-reduction, Au is the easiest of the three metals to reduce, followed by Pd. On the contrary, Ni is the slowest to reduce [95,114]. Therefore, it may not be incorrect to assume that the surface of the trimetallic catalysts (especially PdAuNi/C SBEG ) was richest in Ni and poor in Au, which was evident from the XPS and EDX measurements. On the other hand, the core would be highly populated with Au and lacking a substantial Ni proportion.

Borohydride Oxidation Reaction Studies
The catalytic performance of the three PdAuNi/C electrocatalysts towards the BOR was initially investigated by scanning comparative CVs in 2 M NaOH with 0.03 M NaBH 4 and in pure 2 M NaOH. Along with evaluating the activity of the synthesized trimetallic catalysts, the response of monometallic Pd/C was also evaluated in order to better understand the role played by Au and Ni. The CVs of three catalysts as well as of monometallic Pd/C depicted in Figure 6A-C show the absence of any relevant anodic peaks in the BOR potential region in NaOH solution, implying that the current generated during CVs recorded in the presence of BH 4 − originated from its oxidation. Based on the shape of the CVs and similar OCP values of ca. −0.12 V, it could be assumed that the studied catalysts promoted similar reaction mechanisms. Similar OCP values were reported in other recent studies on BOR at Pd-containing materials [53,54]. Within the scanned potential range, CVs of three different catalysts showed two oxidation peaks in the positive scan, while an additional sharp oxidation peak was evident in the reverse scan. The first oxidation peak appeared as a wide oxidation hump at ca. 0.29 V for the PdAuNi/C SBEG and PdAuNi/C SBIPA catalysts, while in the case of PdAuNi/C 3-step , this peak was more pronounced and positioned at more positive potential (0.38 V). It is believed that this peak can be attributed to H 2 oxidation, which originates from the borohydride hydrolysis reaction, and to the initial oxidation of adsorbed BH 3 OH − [54]. The second anodic peak at~0.8 V for all examined samples can be considered to originate from electrooxidation of BH 4 − as well as BH 3 OH − [48,54]. PdAuNi/C SBIPA gave the best catalytic response, reaching the peak current density of 46 mA cm −2 at 50 mV s −1 , a value considerably higher than those reached employing carbon-supported Pd-Cu [115], Pd-Au [116], and Au-Ni [117] alloy nanoparticles. Additionally, the oxidation peak on the cathodic scan in the potential window from 0.5 to 0.8 V for all three studied samples was most likely related to the oxidation of BH 3 OH − [118], which was formed during the anodic scan and stayed adsorbed on the oxidized electrode surface until it was reactivated by reduction of the surface oxides. BOR kinetic parameters are summarized in Table 4 for each catalyst. For all samples, an upward trend of current density with increasing polarization rate, as well as a small shift in the position of the characteristic BH4 − oxidation peak to more positive potentials were observed. Figure 7A illustrates this behavior for PdAuNi/CSBIPA as the best material. This allowed for the determination of the charge transfer coefficient (α) values using the expression describing the peak potential dependence on the logarithm of polarization rate (Ep vs. ln ν) for irreversible processes [12,121]. The values were found to be relatively similar for all studied catalysts (Table 4), with PdAuNi/CSBIPA showing the highest α value. High values indicate that BOR at examined electrocatalysts was an irreversible and diffusion-controlled process. The obtained values are comparable to the ones reported for catalysts containing Pd, Au, and Ni [12,25,47,79,100,121].
In the case of PdAuNi/CSBIPA ( Figure 7B), a significant increase in anodic current den- The catalytic performance of the synthesized materials was further compared with the response obtained for monometallic Pd/C SBIPA . Figure 6D reveals the inferior catalytic activity of Pd/C SBIPA . Introducing Au and Ni into Pd provided more active sites where BH 4 − adsorption could take place and as such, trimetallic catalysts were expected to exhibit higher activity towards BOR. Improved catalytic properties can be explained by the produced synergistic effect of combining Pd with Au, a metal that is relatively inert towards hydrolysis, and Ni, a 3d transition metal, and the resultant electronic effect. In terms of electronic interaction, as XPS results suggest, the addition of Ni may have influenced the atomic band structure of the other two metals. This can further lead to the reduced adsorption energy of BOR intermediates, thus preserving those active sites for BH 4 − adsorption and further encouraging BOR activity [119].
The PdAuNi/C SBIPA highest catalytic activity most likely resulted from the absence of segregation of Ni. This enabled an improved degree of alloying in comparison with the other two studied catalysts, as confirmed by XRD results. Additionally, the enhanced catalytic activity of PdAuNi/C SBIPA could also be associated with the better dispersed multi-metallic particles of smaller average size, confirmed by TEM analysis and the re-spective EDX elemental maps. Smaller particle size and interparticle distance produced larger surface area and, therefore, were deemed beneficial for catalytic reactions. However, caution must be practiced as Antolini et al. [120] reviewed that decreasing the particle size below 5 nm and interparticle distance below 20 nm adversely impacts the reduction of oxygen, as well as reactions of ethanol and methanol oxidation on Pt surfaces. The situation would be worse for a highly aggregated particle system due to exceedingly smaller interparticle distance. That is a likely reason for the lower PdAuNi/C SBEG performance in the present study. The CVs shown in Figure 6 are also available in specific current density (mA mg metal −1 ) in SI ( Figure S7). BOR kinetic parameters are summarized in Table 4 for each catalyst. For all samples, an upward trend of current density with increasing polarization rate, as well as a small shift in the position of the characteristic BH 4 − oxidation peak to more positive potentials were observed. Figure 7A illustrates this behavior for PdAuNi/C SBIPA as the best material. This allowed for the determination of the charge transfer coefficient (α) values using the expression describing the peak potential dependence on the logarithm of polarization rate (E p vs. ln ν) for irreversible processes [12,121]. The values were found to be relatively similar for all studied catalysts (Table 4), with PdAuNi/C SBIPA showing the highest α value. High values indicate that BOR at examined electrocatalysts was an irreversible and diffusion-controlled process. The obtained values are comparable to the ones reported for catalysts containing Pd, Au, and Ni [12,25,47,79,100,121].
In the case of PdAuNi/C SBIPA ( Figure 7B), a significant increase in anodic current densities when recording the CVs starting from 0.01 M to 0.12 M NaBH 4 concentration was evident. The CVs revealed that the anodic current density value amounted to 8.5 mA cm −2 for the lowest concentration (0.01 M), while a value nine times higher was achieved for the highest NaBH 4 concentration. The observed trend corresponded to enhanced mass transfer, and this behavior agrees with previous reports employing Au-, Pd-, and/or Ni-based electrocatalysts [48,78,79,122,123]. In addition, the oxidation process started at lower potentials but reached its maximum at more positive potential values. The order of reaction (β) values was obtained using ln j vs. ln c plot slope and ranged from 0.91 for PdAuNi/C SBEG to 1.12 for PdAuNi/C 3-step .
CVs obtained at elevated temperatures also exhibited higher current densities for all studied electrocatalysts. This behavior was expected, since higher temperatures lead simultaneously to a decrease in solution viscosity, resulting in a higher BH 4 − diffusion coefficient, and faster electron transfer kinetics. By constructing the Arrhenius plots [121], depicted in the inset of Figure 7C, and applying the Arrhenius equation, the values of the apparent activation energy, E a app , were estimated. PdAuNi/C SBIPA showed the lowest value of 16.7 kJ mol −1 , while the highest value of this parameter was observed for PdAuNi/C SBEG (34 kJ mol −1 ). Similar values were reported for mono-and bimetallic noble metal or/and transition metal-containing materials. For instance, BOR at carbon-supported Pd catalyst (prepared using different synthesis methods and types of carbon support) [53,124] has been shown to proceed with activation energy ranging from as low as 10 kJ mol −1 to 26 kJ mol −1 . Moreover, E a app of BOR at a commercially available Pt/C [22] catalyst was reported to have a value of 34 kJ mol −1 , while alloying with transition metals (Co and Ni) resulted in lower values (25 and 20 kJ mol −1 , respectively).
to have a value of 34 kJ mol −1 , while alloying with transition metals (Co and Ni) resulted in lower values (25 and 20 kJ mol −1 , respectively).
The reaction kinetics were additionally examined using RDE LSVs, and polarization curves were obtained; the corresponding j −1 vs. ω −1/2 plot for the currents taken at 0.9 V (inset) is shown in Figure 7D. From the slope of j −1 vs. ω −1/2 plots and using the Koutecky-Levich equation, n values were obtained [125]. The shift of half-wave potential with an increase of the electrode rotation speed confirmed the irreversibility of the reaction.  The reaction kinetics were additionally examined using RDE LSVs, and polarization curves were obtained; the corresponding j −1 vs. ω −1/2 plot for the currents taken at 0.9 V (inset) is shown in Figure 7D. From the slope of j −1 vs. ω −1/2 plots and using the Koutecky-Levich equation, n values were obtained [125]. The shift of half-wave potential with an increase of the electrode rotation speed confirmed the irreversibility of the reaction.
The n values for PdAuNi/C 3-step and PdAuNi/C SBIPA electrocatalysts (7.4 and 7.1, respectively) were determined to nearly match the theoretical value of 8, while significantly lower values were obtained for the other two materials (2.8 and 1.2 for PdAuNi/C SBEG and Pd/C SBIPA , respectively). Although in terms of current densities, PdAuNi/C 3-step showed significantly lower values than PdAuNi/C SBIPA , the highest n value could be explained by Au@Pd (core@shell) structure formation previously identified in the XRD analysis. Such a structure makes the behavior of the catalyst more Au-like due to the core position of those nanoparticles and as such, less active toward the BH 4 − hydrolysis reaction [12,100,126]. The lower faradaic efficiency of PdAuNi/C SBEG and Pd/C SBIPA was most likely a consequence of the higher hydrolysis rates at those two electrocatalysts. Obtained n values for PdAuNi/C SBEG and Pd/C SBIPA are similar to the ones obtained for other Pd-containing materials reported in previous studies. For instance, Martins et al. [54] reported n values of 2 and 5.6 during BOR at Pd nanoparticles with two different types of biobased carbon supports. For other Pd-based catalysts, depending on NaBH 4 concentration, n values between 2 and 6 [53,118] were obtained, indicating hydrolysis occurrence in parallel with the oxidation of BH 4 − . Song et al. [127] recently reported high BOR current densities for a ternary CuPdNi electrode, but at high overpotentials and an n value of 4.9, far from those obtained for the two best catalysts in this work. Alloying Pd with Au resulted in n~6, as reported by Simões et al. [57], also lower than values than obtained herein for PdAuNi/C SBEG and Pd/C SBIPA . Figure 7 is also available in specific current density ( Figure S8 of SI). Figure 8 (available in SI in specific current density, Figure S9) depicts the polarization and power density curves of a single cell with PdAuNi/C SBIPA anode operating at two different temperatures, 25 and 75 • C, along with anode and cathode potentials measured during the cell operation and experimentally measured cell voltage. The cathode potential measured in zero current regime (open-circuit conditions) at 25 • C was almost identical to the equilibrium potential of H 2 O 2 reduction given in Equation (4) (E 0 = 1.78 V), suggesting direct H 2 O 2 reduction. As for the cell operating at elevated temperature (75 • C), this value was slightly less positive, as expected, implying H 2 O 2 disproportionation was followed by O 2 reduction to a certain extent [18,21]. On the other hand, the anode potential measured during the experiment was significantly higher than the equilibrium potential for BOR described by Equation (1), implying that the anodic process also involved side reactions, such as oxidation of H and BH 3 OH − generated during spontaneous hydrolysis of BH 4 − [128]. Moreover, from Figure 8B,C, a small difference between cell voltage, E cell , and the calculated potential difference of the half-reactions, E c -E a , could be observed, which revealed minor resistance of system components.

Fuel Cell Testing
ing materials reported in previous studies. For instance, Martins et al. [54] reported n values of 2 and 5.6 during BOR at Pd nanoparticles with two different types of biobased carbon supports. For other Pd-based catalysts, depending on NaBH4 concentration, n values between 2 and 6 [53,118] were obtained, indicating hydrolysis occurrence in parallel with the oxidation of BH4 − . Song et al. [127] recently reported high BOR current densities for a ternary CuPdNi electrode, but at high overpotentials and an n value of 4.9, far from those obtained for the two best catalysts in this work. Alloying Pd with Au resulted in n ~6, as reported by Simões et al. [57], also lower than values than obtained herein for PdAuNi/CSBEG and Pd/CSBIPA. Figure 7 is also available in specific current density ( Figure  S8 of SI). Figure 8 (available in SI in specific current density, Figure S9) depicts the polarization and power density curves of a single cell with PdAuNi/CSBIPA anode operating at two different temperatures, 25 and 75 °C, along with anode and cathode potentials measured during the cell operation and experimentally measured cell voltage. The cathode potential measured in zero current regime (open-circuit conditions) at 25 °C was almost identical to the equilibrium potential of H2O2 reduction given in Equation (4) (E 0 = 1.78 V), suggesting direct H2O2 reduction. As for the cell operating at elevated temperature (75 °C), this value was slightly less positive, as expected, implying H2O2 disproportionation was followed by O2 reduction to a certain extent [18,21]. On the other hand, the anode potential measured during the experiment was significantly higher than the equilibrium potential for BOR described by Equation (1), implying that the anodic process also involved side reactions, such as oxidation of H and BH3OH − generated during spontaneous hydrolysis of BH4 − [128]. Moreover, from Figure 8B,C, a small difference between cell voltage, Ecell, and the calculated potential difference of the half-reactions, Ec-Ea, could be observed, which revealed minor resistance of system components.  Considering the pH values of electrolytes, the theoretical voltage value of such fuel cells was expected to be close to 3 V. However, the cell polarization curves show an OCV value of~1.9 V at both temperatures. This fluctuation in OCV value was mainly triggered by the mixed potential occurrence on the anodic side and it agrees with the reports from the literature for DBPFCs studies performed under conditions similar to the ones reported herein [18,125,128]. After the initial voltage decay due to activation losses, cell voltage continued to decrease linearly with current density, while mass transport effects were not evident. The power density value reached its maximum of 47.5 mW cm −2 at 0.6 V and 78 mA cm −2 . As expected, at an elevated temperature of 75 • C, the kinetics of both fuel cell reactions (BOR and HPRR) was faster, and a significant improvement in cell performance was displayed. The maximum value of power density that was obtained at this temperature was almost four times higher than that at 25 • C, amounting to 175 mW cm −2 . This value was exhibited at a cell voltage value of 0.7 V, while current density reached 247 mA cm −2 . The comparison of the DBPFC results obtained herein with the ones from the literature is shown in Table 5.  [128] PdAuNi/C SBIPA exhibited the best performance among the presented electrocatalysts, even when compared with the typical (and expensive) electrocatalysts exclusively composed of noble metals (e.g., Pd [128], Au/C [129]). It shows that the alloying effect, especially between noble and low-cost metals, is one of the best ways to decrease the price of electrocatalysts while increasing their catalytic activity. In fact, the highest activity of PdAuNi/C SBIPA towards BOR, i.e., the most pronounced synergistic effect, might be due to the complete mixing of three metals by the SBIPA procedure, as indicated by XRD analysis. Furthermore, XRD and XPS analysis indicated the decrease of adsorbate-Pd bond strength that could lead to faster release and subsequent oxidation of the reaction intermediates. Thus, PdAuNi/C SBIPA was demonstrated to be 1.5 and 1.6 times more active (in terms of current density) for BOR than Ni@Pd/PANI [130] and PtNi/C [131] electrocatalysts, which reflects the advantage of using a trimetallic alloy, probably due to the mentioned synergic effect. Carbon-supported Au x Ni y Cu z trimetallic electrocatalysts [119] exhibited better results than PdAuNi/C SBIPA at high Au amounts (Au 2 NiCu). For the other compositions, PdAuNi/C SBIPA was demonstrated to be a more effective anode, showing that the affinity between the metals is crucial to promote a notable alloying effect.

Conclusions
PdAuNi carbon-supported trimetallic catalysts were successfully synthesized by three different routes. A PdAuNi/C SBIPA catalyst, synthesized via reduction using a NaBH 4 -2-propanol complex, presented the best alloying effect, while the other two materials showed the presence of segregated Ni species, as confirmed by XRD, TEM, and XPS analyses. The results obtained from CV and LSV measurements revealed that the catalysts performed well towards BOR. The PdAuNi/C SBIPA catalyst exhibited the highest activity, evidenced by the high number of exchanged electrons amounting to 7.1, and the lowest apparent activation energy (ca. 16.7 kJ mol −1 ). Furthermore, DBPFC operating at two different temperatures (25 and 75 • C) employing a PdAuNi/C SBIPA anode showed a maximum power density of 47.5 mW cm −2 at 25 • C, with a significant improvement with a temperature increase and peak power density of 175 mW cm −2 at 75 • C. Considering the economic benefits of alloying and using support materials, and the results obtained herein, it can be concluded that PdAuNi/C SBIPA is an interesting candidate for the application in DBPFC as anode material.