Materials synthesis
Chemicals. All chemicals were used as received without further purification. Platinum (II) acetylacetonate (Pt(acac)2, 97%), copper (II) acetylacetonate (Cu(acac)2, 97%), nickel (II) acetylacetonate (Ni(acac)2, 95%), cobalt (II) acetylacetonate (Co(acac)2, 97%), iron (III) acetylacetonate (Fe(acac)3, 98%), manganese (II) acetylacetonate (Mn(acac)2, 97%), chromium (III) acetylacetonate (Cr(acac)3, 98%), tungsten hexacarbonyl (W(CO)6, 99%), Benzoic acid (99%), phenylboronic acid, perfluoroterephthalonitrile and pipemidic acid were purchased from Shanghai Macklin Biochemical Co., Ltd., China. Aluminum (III) acetylacetonate (Al(acac)3, 98%), molybdenum (II) acetylacetonate (Mo(acac)2, 97%), and N, N-dimethylformamide (DMF, 99.5%) were purchased from Shanghai Aladdin Biochemical Technology Co., Ltd., China. Iridium (III) acetylacetonate (Ir(acac)3, 97%) was purchased from Shanghai Sigma-Aldrich Trading Co., Ltd., China. Ketjen Black ECP-600JD was purchased from Guangdong Canrd New Energy Technology Co., Ltd., China. Ethanol absolute (C2H5OH, 99.7%), acetone (C3H6O, 99.5%), perchloric acid (HClO4, 70%), hydrochloric acid (HCl, 36%) and nitric acid (HNO3, 65%) were purchased from Sinopharm Chemical Reagent Co., Ltd. China. DI water (18.2 MΩ⋅cm) used in all experiments was prepared by passing through an ultra-pure purification system.
Synthesis of nonmetal-hybridized interstitial high-entropy Pt alloy catalysts. Relatively abundant metallic elements (Cu, Ni, Co, Mn, Fe, Cr, V, Ti, Mo, W, Al, Zn) were preferentially investigated. In a typical synthesis process of PtCuNiCoN, 10 mg of platinum (II) acetylacetonate, 6.7 mg of copper (II) acetylacetonate, 6.5 mg of nickel (II) acetylacetonate, 6.5 mg of cobalt (II) acetylacetonate, 150 mg of benzoic acid, 30 mg of pipemidic acid, and 10 mg of carbon black (Ketjen black ECP-600JD) were added in a vial containing 5 mL of N, N-dimethylformamide (DMF). Then the suspension was sonicated for 10 min. The resulting mixture in the vial was subsequently placed in an oil bath which was heated to 160oC, and maintained for 12 h. After that, the reactants in the vial were cooled to room temperature and subjected to centrifugation and being rinsed 3 times in mixed solutions of ethanol and acetone. The black powder retrieved from centrifuge tubes were dried in a vacuum oven overnight. The obtained material was then annealed under flowing H2/Ar (10% H2) gas at 400oC (5oC /min) for 2 h and cooled down to room temperature. Finally, an acid wash in 0.2 M HNO3 for 15 minutes had been applied to the annealed material and then rinsed three times with deionized water. After drying overnight, the pristine sample of N-hybridized interstitial high-entropy Pt alloy catalysts have been prepared.
A series of catalysts with different compositions have been synthesized using the same method (described above) by adjusting the quantities of corresponding metal precursors. Similarly, the catalysts with incorporated boron or fluorine were prepared in the same way as for incorporating nitrogen, where a boron (phenylboronic acid) or fluorine source (perfluoroterephthalonitrile) have been used instead. Of note, before MEA tests, the pristine catalyst samples were further washed by 0.2 M HNO3 at 85oC in an Ar atmosphere within a sealed three-neck flask for 8 h and then dried overnight at 60oC in vacuo.
Synthesis of (high-entropy) Pt alloy catalysts without interstitial nitrogen. The (high-entropy) Pt alloy catalysts without N-incorporation (such as, PtCuNiCoFe, PtCuNiFeMn, PtCuNiCo and so on) were also prepared following the same procedure as described above, but simply with no addition of nonmetal sources during the solvothermal reaction.
Transmission electron microscopy (TEM) and high-resolution transmission electron microscopy (HRTEM) images were obtained with a JEM 2100 microscope. High-angle annular dark field scanning transmission electron microscopy (HAADF-STEM) images and energy dispersive X-ray spectroscopy (EDS) elemental mapping analysis were carried out by using FEI Themis Z with 300 kV acceleration voltage. Electron energy loss spectroscopy (EELS) was taken on the same microscope equipped with a GIF Quantum 1065, both EEL spectra in Extended Data Fig. 1e being processed by superimposing the signal from multiple test areas (nanoparticle areas or dimmer single atom areas) and being slightly smoothed. The annular bright field (ABF) STEM images were acquired on an aberration-corrected JEOL JEM-ARM300F2 microscope with camera length at 10 cm (operating at 300 kV). The integrated differential phase contrast (iDPC)-STEM images were acquired in a ThermoFisher Scientific Spectra 300 aberration-corrected TEM operating at 300 kV. A convergence semi-angle of 25 mrad and current of 5 pA were used for imaging.
X-ray diffraction (XRD) patterns were collected by using a DX-2700B powder diffractometer (Dandong Haoyuan Instrument Co., Ltd.), operated at 40 kV and 30 mA, using a Cu-K radiation source (λ = 1.5405 Å). X-ray photoelectron spectroscopy (XPS) was carried out on an America Thermo Scientific K-Alpha surface analysis system employing a monochromatic Al-K X-ray source. The actual metal contents were determined by the inductively coupled plasma-optical emission spectrometer (ICP-OES, Perkin Elmer Optima 8300).
To prepare the APT samples, an ethanol suspension with the nanoparticles was dropped onto an aluminum sample holder. After drying, a 150 nm Cr layer was deposited on the surface with a Leica EM ACE600 sputter coater to embed the nanoparticles. APT specimens were prepared via a standard lift-out procedure in a FEI G4 CX focused ion beam/scanning electron microscope (FIB/SEM). The APT experiments were carried out on a Cameca LEAP 5000 XR in a laser mode at a specimen temperature of 65 K, a pulse rate of 125 kHz, a 0.3% detection rate and a laser pulse energy of 30 pJ. Reconstruction and analysis of the APT data were performed using the commercial AP Suite 6.3.0.90 software.
X-ray absorption fine structure (XAFS) spectroscopy was carried out at synchrotron-based light source (BSRF 1W1B) and using the Rapid XAFS 2M (Anhui Absorption Spectroscopy Analysis Instrument Co., Ltd.) by transmission (or fluorescence mode) at 20 kV and 30 mA, and the Si (533) spherically bent crystal analyzer with a radius of curvature of 500 mm was used for Co, the Si (771) spherically bent crystal analyzer with a radius of curvature of 500 mm was used for Pt, the Si (553) spherically bent crystal analyzer with a radius of curvature of 500 mm was used for Cu, the Si (551) spherically bent crystal analyzer with a radius of curvature of 500 mm was used for Ni.
Data reduction, data analysis, and EXAFS fitting were performed and analyzed with the Athena and Artemis programs of the Demeter data analysis packages57 that utilizes the FEFF6 program58 to fit the EXAFS data. The energy calibration of the sample was conducted through standard Pt foil, which as a reference was simultaneously measured. A linear function was subtracted from the pre-edge region, then the edge jump was normalized using Athena software. The χ(k) data were isolated by subtracting a smooth, third-order polynomial approximating the absorption background of an isolated atom. The k3-weighted χ(k) data were Fourier transformed after applying a HanFeng window function (Δk = 1.0). For EXAFS modeling, the global amplitude EXAFS (CN, R, σ2 and ΔE0) were obtained by nonlinear fitting, with least-squares refinement, of the EXAFS equation to the Fourier-transformed data in R-space using Artemis software. EXAFS of the Pt foil were fitted and the obtained amplitude reduction factor S02 value (0.772) was set fixed in the EXAFS analysis to determine the coordination numbers (CNs) in the Pt-N, Pt-Cu/Ni/Co and Pt-Pt scattering path in the samples. S02 was fixed to 0.772, according to the experimental EXAFS fit of Pt foil by fixing CN as the known crystallographic value. Fitting range: 3.0 ≤ k (/Å) ≤ 13.0 and 1.0 ≤ R (Å) ≤ 3.2 (Pt foil); 3.0 ≤ k (/Å) ≤ 12.0 and 1.2 ≤ R (Å) ≤ 2.1 (PtO2); 3.0 ≤ k (/Å) ≤ 12.0 and 1.0 ≤ R (Å) ≤ 3.3 (PtCuNiCo); 3.0 ≤ k (/Å) ≤ 12.0 and 1.0 ≤ R (Å) ≤ 3.3 (PtCuNiCoN). A reasonable range of EXAFS fitting parameters: 0.700 < Ѕ02 < 1.000; CN > 0; σ2 > 0 Å2; |ΔE0| < 10 eV; R factor < 0.02.
Electrochemical measurements
Rotating disk electrode (RDE) measurements. All electrochemical measurements were performed at an electrochemical workstation (Vertex. 100 mA. EIS Bistat) at room temperature, using a three-electrode electrochemical setup with a rotating disk electrode (RDE) system. A working electrode of glassy carbon (5 mm inner diameter, 0.196 cm2), a counter electrode of platinum wire and a reference electrode of Ag/AgCl were used for the tests. In this work, all potentials were quoted with respect to a reversible hydrogen electrode (RHE). The RHE calibration was performed in a highly H2-saturated 0.1 M HClO4 electrolyte with a platinum disk as the working electrode.
For the preparation of RDE test electrode, 1.5 mg catalyst was mixed with 990 µL ethanol and 10 µL Nafion solution (5 wt%) and the mixture was ultrasonicated for 120 minutes to form the catalyst ink. Prior to the catalyst ink deposition, the GC disk was polished with 50 nm alumina powder and rinsed generously with high purity ethanol. Then the ink was deposited on glass carbon electrode to provide a uniform catalyst layer (CL). For comparison, the RDE electrode of commercial Pt/C (TKK 20%) was prepared with the same method.
Cyclic voltammetry (CV) measurement was typically carried out in the potential range of 0.05–1.1 V (vs. RHE) at a scan rate of 50 mV s− 1 in a N2-saturated 0.1 M HClO4 solution. Linear sweep voltammetry (LSV) curves were tested in an O2-saturated 0.1 M HClO4 electrolyte at a rotation speed of 1600 rpm and a scan rate of 10 mV s− 1. The activation test was operated with 150 cycles of CV in an O2-saturated atmosphere. Accelerated durability test was performed with 30,000 cycles of CV in the voltage range of 0.6–1.1 V. After the activation or durability test, the catalysts could be collected by sonicating the GCE in ethanol for further structural and compositional analysis.
The electrochemical active surface area (ECSA) of the catalyst was derived from the average charges associated with the adsorption and desorption of hydrogen between 0.05 and 0.4 VRHE, taking a value of 210 µC cm− 2 for monolayer adsorption of hydrogen on the Pt surface. The MA was defined as the SA multiplied by the ECSA. For the calculation of the ECSA, we used the loading, as measured by ICP-OES.
Membrane electrode assembly (MEA) fabrication and single fuel cell test. The activity and durability of nitrogen-hybridized interstitial HEA catalysts as cathodes were evaluated at a fuel cell testing system (Dalian Yuke Innovation Technology Ltd.) and a single cell with serpentine flow fields. The cathode catalyst ink was prepared by ultrasonically mixing the HEA catalyst, a solvent (75 wt% n-propanol and 25 wt% H2O) and 20% Nafion ionomer solution at an ionomer to carbon ratio of 0.8 for 40 min. Then the uniform ink was sprayed (Xiniu Ultrasonic Coating System 260E) on the one side of Nafion membrane (GORE 12 µm) to form the cathode catalyst layer with an active area of 5 cm2. Extra care must be taken to avoid wrinkling of the membrane during catalyst coating. The anode catalyst layer was prepared by coating the commercial TKK 20% (TEC10E20E) Pt/C catalyst onto the other side of the Nafion membrane using the same method as above. The catalyst-coated membrane (CCM) containing cathode and anode catalyst layers, two gas diffusion layers (GDLs) including microporous layers (SGL 29BC) and two polyimide gaskets (used to better control the compression) were hot-pressed to obtain a MEA at 120oC for 5 min under a pressure of 4 MPa. Both the anode and cathode Pt loadings of the prepared MEAs have been measured by ICP-OES and repeated for at least 3 times to ensure an accurate evaluation.
H 2 -O 2 single fuel cell tests: PtCuNiCoN was used as the cathode catalyst at 0.05 mgPt cm− 2 and Pt/C was used for the anode at 0.05 mgPt cm− 2. The temperature of cells was maintained at 80oC. Pure hydrogen (99.999%) and O2 (99.999%) were with 100% relative humidity (RH) and the back pressure of both the anode and cathode was 150 kPaabs, outlet, while keeping the flow rate of H2/O2 at 200/500 mL min− 1. Before recording the polarization curves, MEAs were activated via keeping the cell voltage at 0.5 V until the stable current value was obtained. Pt mass activity (MA) was determined at 0.9 ViR−correct and H2/O2 flow rates of 200/500 mL min− 1.
H 2 -air single fuel cell tests: PtCuNiCoN cathodes with Pt loadings of 0.07 and 0.05 mg cm− 2 were selected, while Pt/C of 0.03, 0.01 mg cm− 2 was accordingly used for the anode. The cell temperature was kept at 80oC at 100% RH, the anode/cathode flow rates were fixed at 500/1000 mL min− 1, and the backpressure was 250 kPaabs,outlet. Prior to the performance testing, an activation step by maintaining the cell voltage at 0.5 V was conducted for at least three hours. We also measured the polarization curves using constant current density holds, which showed very similar performance as obtained via constant voltage holds. The H2-air polarization curves of TKK 20% Pt/C as cathode catalysts with Pt loadings of 0.07 and 0.05 mg cm− 2, respectively, were also measured using the same method.
Accelerated stress test (AST): In order to evaluate the durability of the catalysts, the AST was performed using 30,000 cycles of 0.6 to 0.95 V square-wave voltage with a duration of 3 s for each voltage, following DOE fuel cell test protocol. The fuel cell test condition in AST was set at 80oC, 150 kPaabs, outlet, and 100% RH with a H2/N2 flow rate of 200/75 mL min− 1 for the anode/cathode. After the completion of voltage cycles, the H2-O2 current-voltage polarization curve was recorded again to measure the change of MA. Chronoamperometric measurements were conducted on MEAs fabricated using catalyst loadings of 0.03 mg/cm2 TKK 20wt% Pt/C (anode), 0.07 mg/cm2 PtCuNiCoN (cathode) and 12 µm Nafion membranes, under 150 kPaabs, outlet H2-air conditions.
All fuel cell impedance spectra were acquired with an electrochemical workstation (Gamry instrument, reference 3000). High-frequency resistance (HFR) was in-situ measured during the polarization curve test. The perturbation amplitude for the AC impedance was 5% of the direct current, and the frequency was 10 kHz, which would not disturb the electrochemical reaction of fuel cell. The EIS was tested from 20 kHz to 1 Hz, and 1 to 15 A under the conditions of H2 (500 mL min− 1) and air (1000 mL min− 1) at 80ºC and 100% RH.
Calculation Methods
Spin-polarized density functional theory (DFT) calculations were performed using LASP (www.lasphub.com) program contained VASP 6.2.1 packages59 with projected augmented wave (PAW) pseudo-potentials60,61. The exchange-correlation energy was treated based on the generalized gradient approximation (GGA) by using Perdew-Burke-Ernzerhof (PBE) function62. The plane-wave cutoff energy was set to 450 eV. The DFT-D3(BJ) method of Grimme63,64 was employed to describe long-range vdW interactions. In this work, it was uniformly replaced in Pt bulk according to the experimental element ratio (Pt: Cu: Ni: Co = 25:8:13:2) to optimize the structure for PtCuNiCo high entropy alloy (HEA) bulk. The 11% N-doped was distributed in lower energy tetrahedron sites depended on the comparison of tetrahedron, octahedron and replace sites. The surface models for PtCuNiCo and PtCuNiCoN HEAs were cut out 5 layers from the bulk, which ensure top layer was the lowest energy surface with 16 atoms. For Co, Ni, Cu or Pt single-atom model, the most stable structure65,66 was selected. The Monkhorst-Pack scheme with a k-point separate on length of 0.05 Å−1 was utilized for sampling the first Brillion zone67. All atoms were fully relaxed in bulk optimizations, and 2 layers were relaxed while the rest layers were fixed in surface optimizations. The Quasi-Newton l-BFGS method was used for geometry relaxation until the maximal force on each degree of freedom less than 0.01 eV/Å.
To obtain the free energy profile, we have calculated Gibbs free energies for all the states based on the ZPE-corrected DFT total energy, which was set as the enthalpy at 0 K. It can be calculated with
where EDFT is the total energy getting from DFT optimization, EZPE is the zero-point vibrational energy using the vibrational frequency calculations were performed via the finite-difference approach68, CV is the heat capacity, T is the kelvin temperature, and S is the entropy. The free energy for the gas-phase molecule was computed from thermodynamics by utilizing the standard thermodynamics data at the standard state.
57. Ravel, B. & Newville, M. ATHENA, ARTEMIS, HEPHAESTUS: data analysis for X-ray absorption spectroscopy using IFEFFIT. Journal of Synchrotron Radiation 12, 537-541 (2005).
58. Zabinsky, S. I., Rehr, J. J., Ankudinov, A., Albers, R. C. & Eller, M. J. Multiple-scattering calculations of X-ray-absorption spectra. Phys. Rev. B 52, 2995-3009 (1995).
59. Kresse, G. & Furthmuller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Computational Materials Science 6, 15-50 (1996).
60. Blochl, P. E. Projector augmented-wave method. Phys. Rev. B 50, 17953-17979 (1994).
61. Kresse, G. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. Rev. B 59, 1758-1775 (1999).
62. Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 77, 3865-3868 (1996).
63. Grimme, S., Antony, J., Ehrlich, S. & Krieg, H. A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. J. Chem. Phys. 132, 154104 (2010).
64. Grimme, S., Ehrlich, S. & Goerigk, L. Effect of the damping function in dispersion corrected density functional theory. J. Comput. Chem. 32, 1456-1465 (2011).
65. Venegas, R. et al. Comparison of the catalytic activity for O2 reduction of Fe and Co MN4 adsorbed on graphite electrodes and on carbon nanotubes. Physical Chemistry Chemical Physics 19, 20441-20450 (2017).
66. Chen, C., Li, Y. F., Cheng, D., He, H. & Zhou, K. B. Graphite nanoarrays-confined Fe and Co single-atoms within graphene sponges as bifunctional oxygen electrocatalyst for ultralong lasting zinc-air battery. ACS Appl. Mater. Interfaces 12, 40415-40425 (2020).
67. Monkhorst, H. J. & Pack, J. D. Special points for Brillouin-zone integrations. Phys. Rev. B 13, 5188-5192 (1976).
68. Scott, A. P. & Radom, L. Harmonic vibrational frequencies: an evaluation of Hartree-Fock, Moller-Plesset, quadratic configuration interaction, density functional theory, and semiempirical scale factors. J. Phys. Chem. 100, 16502-16513 (1996).
69. Bu, L. et al. PtPb/PtNi Intermetallic core/atomic layer shell octahedra for efficient oxygen reduction electrocatalysis. J. Am. Chem. Soc. 139, 9576-9582 (2017).
70. Chen, S. et al. Pd–Pt tesseracts for the oxygen reduction reaction. J. Am. Chem. Soc. 143, 496-503 (2021).
71. Chen, T. et al. PtFeCoNiCu high-entropy solid solution alloy as highly efficient electrocatalyst for the oxygen reduction reaction. Iscience 26, 105890 (2023).
72. Cheng, Q. et al. High-loaded sub-6 nm Pt1Co1 intermetallic compounds with highly efficient performance expression in PEMFCs. Energ. Environ. Sci. 15, 278-286 (2022).
73. Gao, L. et al. Unconventional p–d hybridization interaction in PtGa ultrathin nanowires boosts oxygen reduction electrocatalysis. J. Am. Chem. Soc. 141, 18083-18090 (2019).
74. Guo, W. et al. A closely packed Pt1.5Ni1-x/Ni–N–C hybrid for relay catalysis towards oxygen reduction. Energ. Environ. Sci. 16, 148-156 (2023).
75. He, D. et al. Amorphous nickel boride membrane on a platinum–nickel alloy surface for enhanced oxygen reduction reaction. Nat. Commun. 7, 12362 (2016).
76. He, D. S. et al. Ultrathin icosahedral Pt-enriched nanocage with excellent oxygen reduction reaction activity. J. Am. Chem. Soc. 138, 1494-1497 (2016).
77. Huang, H. et al. Achieving remarkable activity and durability toward oxygen reduction reaction based on ultrathin Rh-doped Pt nanowires. J. Am. Chem. Soc. 139, 8152-8159 (2017).
78. Huang, J. et al. Experimental Sabatier plot for predictive design of active and stable Pt-alloy oxygen reduction reaction catalysts. Nat. Catal. 5, 513-523 (2022).
79. Huang, X. et al. High-performance transition metal-doped Pt3Ni octahedra for oxygen reduction reaction. Science 348, 1230-1234 (2015).
80. Jiang, K. et al. Ordered PdCu-based nanoparticles as bifunctional oxygen-reduction and ethanol-oxidation electrocatalysts. Angew. Chem. Int. Ed. 55, 9030-9035 (2016).
81. Li, C. et al. PtCoNi ternary intermetallic compounds anchored on Co, Ni and N co-doped mesoporous carbon: Synergetic effect between PtCoNi nanoparticles and doped mesoporous carbon promotes the catalytic activity. J. Energy Chem. 78, 340-349 (2023).
82. Li, J. et al. Fe stabilization by intermetallic L10-FePt and Pt catalysis enhancement in L10-FePt/Pt nanoparticles for efficient oxygen reduction reaction in fuel cells. J. Am. Chem. Soc. 140, 2926-2932 (2018).
83. Li, K. et al. One-nanometer-thick PtNiRh trimetallic nanowires with enhanced oxygen reduction electrocatalysis in acid media: integrating multiple advantages into one catalyst. J. Am. Chem. Soc. 140, 16159-16167 (2018).
84. Li, S. et al. Nanoporous high-entropy alloys with low Pt loadings for high-performance electrochemical oxygen reduction. J. Catal. 383, 164-171 (2020).
85. Liang, J. et al. Biaxial strains mediated oxygen reduction electrocatalysis on fenton reaction resistant L10-PtZn fuel cell cathode. Adv. Energy Mater. 10, 2000179 (2020).
86. Mao, Z. et al. Interstitial B-doping in Pt lattice to upgrade oxygen electroreduction performance. ACS Catal. 12, 8848-8856 (2022).
87. Sievers, G. W. et al. Self-supported Pt-CoO networks combining high specific activity with high surface area for oxygen reduction. Nat. Mater. 20, 208-213 (2021).
88. Tu, W. et al. Tungsten as "adhesive" in Pt2CuW0.25 ternary alloy for highly durable oxygen reduction electrocatalysis. Adv. Funct. Mater. 30, 1908230 (2020).
89. Wu, Z.-P. et al. Alloying–realloying enabled high durability for Pt–Pd-3d-transition metal nanoparticle fuel cell catalysts. Nat. Commun. 12, 859 (2021).
90. Xiong, Y. et al. N-doping induced tensile-strained Pt nanoparticles ensuring an excellent durability of the oxygen reduction reaction. J. Catal. 382, 247-255 (2020).
91. Zhang, N. et al. Superior bifunctional liquid fuel oxidation and oxygen reduction electrocatalysis enabled by PtNiPd core–shell nanowires. Adv. Mater. 29, 1603774 (2017).
92. Zaman, S. et al. Scalable molten salt synthesis of platinum alloys planted in metal–nitrogen–graphene for efficient oxygen reduction. Angew. Chem. Int. Ed. 61, e202115835 (2022).
93. Dionigi, F. et al. Controlling near-surface Ni composition in octahedral PtNi(Mo) nanoparticles by Mo doping for a highly active oxygen reduction reaction catalyst. Nano Lett. 19, 6876-6885 (2019).
94. Liu, X. et al. Inducing covalent atomic interaction in intermetallic Pt alloy nanocatalysts for high-performance fuel cells. Angew. Chem. Int. Ed. 62, e202302134 (2023).
95. Yin, P. et al. Machine-learning-accelerated design of high-performance platinum intermetallic nanoparticle fuel cell catalysts. Nat. Commun. 15, 415 (2024).
96. Gao, L. et al. Identifying the distinct roles of dual dopants in stabilizing the platinum-nickel nanowire catalyst for durable fuel cell. Nat. Commun. 15, 508 (2024).
97. Liang, J. et al. Metal bond strength regulation enables large-scale synthesis of intermetallic nanocrystals for practical fuel cells. Nat. Mater. (2024).
98. Zeng, Y. et al. Regulating catalytic properties and thermal stability of Pt and PtCo intermetallic fuel-cell catalysts via strong coupling effects between single-metal site-rich carbon and Pt. J. Am. Chem. Soc. 145, 17643-17655 (2023).