Heterophase fcc-2H-fcc gold nanorods

The crystal phase-based heterostructures of noble metal nanomaterials are of great research interest for various applications, such as plasmonics and catalysis. However, the synthesis of unusual crystal phases of noble metals still remains a great challenge, making the construction of heterophase noble metal nanostructures difficult. Here, we report a one-pot wet-chemical synthesis of well-defined heterophase fcc-2H-fcc gold nanorods (fcc: face-centred cubic; 2H: hexagonal close-packed with stacking sequence of “AB”) at mild conditions. Single particle-level experiments and theoretical investigations reveal that the heterophase gold nanorods demonstrate a distinct optical property compared to that of the conventional fcc gold nanorods. Moreover, the heterophase gold nanorods possess superior electrocatalytic activity for the carbon dioxide reduction reaction over their fcc counterparts under ambient conditions. First-principles calculations suggest that the boosted catalytic performance stems from the energetically favourable adsorption of reaction intermediates, endowed by the unique heterophase characteristic of gold nanorods.

C rystal phase control is one of the effective strategies of phase engineering of nanomaterials (PEN) 1 to modulate the intrinsic physicochemical properties of materials via delicately tuning the order of atom arrangements [2][3][4][5][6][7] . Although noble metal nanomaterials have been widely investigated for decades, the previous studies almost focus on the size-, dimension-, shape-, architecture-, facet-, lattice strain-and composition-control of thermodynamically stable crystal phase [8][9][10][11][12][13][14][15][16][17][18][19] . Up to now, it remains difficult to synthesize noble metal nanomaterials with unusual crystal phase compared to their conventional thermodynamically stable phase. Taking gold (Au) as a typical example, although the colloidal synthesis of Au nanoparticles (NPs) dates back to the 1850s 20 , nearly all Au nanostructures synthesized by wet-chemical methods adopt the thermodynamically stable face-centred cubic (fcc) phase.
Recently, Au nanomaterials with the unusual hexagonal closepacked (hcp, the 2H type) and 4H phases have been successfully synthesized in solutions 21,22 . However, the high-yield synthesis of 2H and 2H-based Au nanomaterials is yet to be realized, restraining the exploration of their intriguing properties and potential applications.
Owing to their unique optical property, distinct electronic structure, good biocompatibility and extraordinary stability, noble metal nanomaterials possess diverse applications in plasmonics 23 , catalysis [14][15][16]18,24,25 , biosensing 26 and in vivo therapy 27 . In particularly, noble metal nanomaterials can be used as electrocatalysts toward carbon dioxide (CO 2 ) reduction reaction (CO 2 RR) that converts CO 2 into valuable chemicals, such as carbon monoxide (CO), ethylene, methane and formic acid [28][29][30][31][32] . Among these products, CO is an important product that can be readily used in downstream Fischer-Tropsch (FT) synthesis for various chemicals and fuels 33 . To date, although great efforts, including size-, composition-, morphology-, surface-and defectcontrol of noble metal catalysts have been devoted to enhancing CO production from CO 2 RR 34-38 , they still suffer from the low mass activity. Note that the crystal phase of almost all the reported noble metal nanocatalysts for CO 2 RR is limited to their thermodynamically stable crystal phase. On the other hand, crystal phase-based heterostructures have recently attracted increasing research interests because the synergistic effect between two different crystal phases (e.g. phase junction) of inorganic nanomaterials could significantly improve their physicochemical properties, such as optical, electronic and catalytic properties 4,5,[39][40][41] . Unfortunately, it is still a big challenge to synthesize well-defined heterophase noble metal nanomaterials.
Here, we report the direct synthesis of Au nanorods (NRs) with a well-defined fcc-2H-fcc heterophase under mild conditions. By combining experimental and theoretical studies at the single particle level, we identify the different optical properties in fcc, 2H and heterophase fcc-2H-fcc Au NRs. As a proof-of-concept application, we find that the heterophase fcc-2H-fcc Au NRs exhibit much greater electrocatalytic activity compared to the fcc Au NRs and Au NPs in the CO 2 RR. Theoretical calculations reveal that the energetically favourable adsorption of reaction intermediates accounts for the enhanced catalytic performance of heterophase Au NRs.

Results
Synthesis and characterization. The heterophase fcc-2H-fcc Au NRs were synthesized by a one-pot approach in which the potassium tetrachloroaurate was reduced in a mixture of oleylamine and n-dodecylamine at 65°C (See "Methods" for details). Two crystal phases, i.e. fcc and 2H, of Au are schematically illustrated in Fig. 1a, b, respectively, in which the close-packed planes exhibit different characteristic stacking sequences, i.e. "ABC" for the fcc phase and "AB" for the 2H phase. The X-ray photoelectron spectroscopy (XPS) core level spectrum of the Au 4f doublet confirms the metallic state of Au in the fcc-2H-fcc Au NRs ( Supplementary Fig. 1). The transmission electron microscopy (TEM) and high-angle annular dark-field scanning TEM (HAADF-STEM) images reveal the size distribution of the obtained Au NRs (Fig. 1c, d), i.e. the aspect ratio is~2, with a length of 28.2 ± 5.7 nm and a width of 14.6 ± 3.2 nm (Supplementary Fig. 2). It is worth mentioning that the correct choice of Au precursor (Supplementary Fig. 3 and Supplementary Note 1), solvent (Supplementary Fig. 4 and Supplementary Note 2) and gas atmosphere (Supplementary Fig. 5 and Supplementary Note 3) in the reaction is important for the successful synthesis of heterophase fcc-2H-fcc Au NRs.
The crystal phase of as-prepared Au NRs was examined with aberration-corrected HAADF-STEM. Figure 1e shows a representative aberration-corrected high-resolution HAADF-STEM image of an Au NR, which exhibits distinct types of atom arrangements along its long axis. Selected-area fast Fourier transform (FFT) patterns from the two ends of Au NRs match well with the typical [101] f -zone axis diffraction pattern of fcc phase, exhibiting the diffraction spots of (020) f and (11 1) f planes (Fig. 1f, h). In contrast, the selected-area FFT pattern from the middle of Au NRs is consistent with the characteristic [110] hzone axis diffraction pattern of the 2H phase, demonstrating the diffraction spots of the (1 10) h and (002) h planes (Fig. 1g). As indicated by the dashed-dotted white lines in Fig. 1e, there are sharp interfaces between two different phases, and a well-defined sandwich structure of fcc-2H-fcc heterophase along the closepacked directions of [111] f /[001] h is formed. Figure 1i shows the integrated pixel intensities of 2H (002) h and fcc (11 1) f lattices from the selected areas in Fig. 1e. The average interlayer spacing of (002) h planes was calculated to be 2.35 Å (red curve), which is slightly (ca. 1.3%) greater than that of (11 1) f planes (2.32 Å, green and purple curves). Based on the aforementioned result, an atomic model of the as-prepared heterophase fcc-2H-fcc Au NRs is schematically illustrated in Fig. 1j.
To further corroborate the heterophase structure in the synthesized fcc-2H-fcc Au NRs, high-resolution TEM (HRTEM) images from the other two directions, i.e. Y direction and Z direction ( Supplementary Fig. 6), were collected. The corresponding FFT pattern of HRTEM image collected from the Y direction agrees well with the diffraction patterns of [1 10] Fig. 8), providing further evidence for the unique fcc-2H-fcc heterophase structure of Au NRs. Besides, the structure of fcc-2H-fcc Au NRs was also confirmed by the X-ray diffraction (XRD) pattern, which demonstrates two sets of diffraction peaks belonging to 2H and fcc phases ( Supplementary Fig. 9).
Optical property. Monochromated electron energy loss spectroscopy (EELS) was used to investigate the electromagnetic eigenmodes of fcc-2H-fcc Au NRs at the single particle level (Fig. 2a, see "Methods" for details). The monochromated EELS experiments were conducted by raster-scanning an electron probe with beam size of~1 nm 2 within a rectangular area containing an individual fcc-2H-fcc Au NR. Figure 2b shows the HRTEM image of a fcc-2H-fcc Au NR with length, width and aspect ratio of 29.8 nm, 14.6 nm and 2.0, respectively. The corresponding EELS spectra of a single fcc-2H-fcc Au NR were collected from different locations, as marked in Fig. 2b. Two distinct eigenmode resonances were observed with one at 2.48 eV and the other at 2.07 eV excited by placing the electron beam (E-beam) in the middle and end of fcc-2H-fcc Au NR, respectively (Fig. 2c). In addition, the monochromated EELS maps demonstrate that the high-energy eigenmode is distributed over the whole surface of fcc-2H-fcc Au NR, while the low-energy eigenmode is excited at the two ends (Fig. 2d, e).
In order to reveal the effect of crystal phase on the electromagnetic eigenmodes of Au NRs, density functional theory (DFT) computation ( Supplementary Fig. 10) and threedimensional electromagnetic full-wave EELS simulation were used to calculate the EELS spectra (see "Methods" for details). Figure 2f-h show the calculated EELS spectra for three different types of Au NRs, i.e. pure fcc Au, pure 2H Au, and heterophase fcc-2H-fcc Au, as schematically illustrated in the insets. To mimic the EELS experiment, we excite Au NRs with the E-beam at two different locations. A low-energy peak is dominant in the end (red curve) of the Au NR, while a high-energy peak is more prominent in its middle (blue curve), consistent with the experimental observation (Fig. 2c). In comparison with the pure fcc Au NR (Fig. 2f), the optical responses of pure 2H and heterophase fcc-2H-fcc Au NRs exhibit both position shifts and intensity variations for all the eigenmode resonances in the range of 1.5-4.0 eV (Fig. 2g, h). To uncover the nature of the eigenmodes associated with the two peaks of fcc-2H-fcc Au NRs, we performed a Helmholtz decomposition analysis to reveal the contributions of the transverse and longitudinal fields. The peaks observed in the EELS probability plot ( Supplementary  Fig. 11) are confirmed to be resonances on surface plasmon polaritons. The fields, excited by the E-beam located in the middle of fcc-2H-fcc Au NRs, are shown in Fig. 2i, j. The electric field is concentrated at two ends for the low-energy mode (Fig. 2i), while for the high-energy mode, it distributes over the entire particle (Fig. 2j), agreeing well with the experimental results (Fig. 2d, e).   Catalytic property. It is well recognized that the surface property of noble metal catalysts plays a key role in dictating their catalytic performance 42,43 . In electrochemical CO 2 RR, the size 31 , shape 44 , grain boundary 37,45 and surface ligand 46 of catalysts have been identified as key factors that can significantly alter their surface properties, and hence their catalytic performance. These studies provide deeper understanding of the CO 2 RR mechanism, which in turn promotes the development of CO 2 reduction technology.
Here, by using the synthesized unique heterophase fcc-2H-fcc Au NR as a new model catalyst, we explore the critical role of the unusual 2H phase as well as the interface between 2H and fcc phases in CO 2 RR. We investigated the CO 2 RR performance of fcc-2H-fcc Au NRs by using the H-type cell with three-electrode system under ambient conditions and analysed the gaseous products with an on-line gas chromatography (See "Methods" for details). For comparison, fcc Au NRs and fcc Au NPs were synthesized (Supplementary Figs. 12 and 13) and used for the CO 2 RR test under the same experimental conditions. All catalysts were thoroughly washed to remove their surface capping ligands ( Supplementary Fig. 14). As shown in Fig. 3a, fcc-2H-fcc Au NRs exhibited over 90% selectivity towards CO production from −0.5 to −0.8 V (vs. RHE) and reached the highest Faradaic efficiency (FE) of 98.2% at −0.6 V (vs. RHE). The competing hydrogen evolution reaction was almost entirely suppressed (Supplementary Fig. 15). In contrast, the highest FE for fcc Au NRs and fcc Au NPs were only 81.6% and 63.6%, respectively. In addition, the linear sweep voltammetry (LSV) measurements showed that fcc-2H-fcc Au NRs have lower onset potential and much higher current density compared to fcc Au NRs and fcc Au NPs (Fig. 3b). As a result, fcc-2H-fcc Au NRs showed greater CO partial current density (J CO ) in a wide range of potentials (Fig. 3c). At the applied potential of −0.6 V (vs. RHE), J CO reached −8.03 mA cm −2 , which was 2.2 and 3.2 times of that of fcc Au NR and fcc Au NP, respectively. These results suggest that the heterophase fcc-2H-fcc Au nanorods possess superior catalytic selectivity and activity in CO 2 RR, compared with their fcc counterparts. Notably, the mass activity at the potential of −0.6 V (vs. RHE) reached 30.8 A g −1 ( Supplementary Fig. 16), which is higher compared to the previously reported noble metal catalysts (Supplementary Table 1).
To gain insight into the origin of the superior CO 2 RR performance of fcc-2H-fcc Au NRs, we first measured the specific activity by normalizing the current density to the electrochemically active surface area (ECSA). Note that ECSAs of the three Au catalysts were estimated by underpotential deposition of copper on their surfaces ( Supplementary Fig. 17). As shown in Fig. 3d, fcc-2H-fcc Au NRs showed much higher specific activity than that of their fcc counterparts. In specific, at the applied potential of −0.6 V (vs. RHE), fcc-2H-fcc Au NR demonstrated a high specific activity of −0.64 mA cm −2 , which is 2.1 and 3.4 times that of fcc Au NR and fcc Au NP, respectively. These results indicate that the intrinsic catalytic activity of the heterophase fcc-2H-fcc Au NRs is much higher than that of their fcc counterparts.
As shown in Fig. 3e, based on the plot of the logarithm of geometric J CO versus overpotential, the Tafel slope was calculated by linear fitting at low overpotentials. As a result, the heterophase fcc-2H-fcc Au NR exhibited a much smaller Tafel slope of 64.9 mV dec −1 , compared to the fcc Au NR (84.6 mV dec −1 ) and fcc Au NP (89.4 mV dec −1 ). The values of these Tafel slopes revealed that all three Au catalysts shared a common reaction mechanism. It means that a fast one-electron transfer to the adsorbed CO 2 to generate the surface-adsorbed CO 2 •− occurs first, and the subsequent combination of CO 2 •− and a proton to form *COOH is the rate-limiting step (Fig. 3f) 30,46 . Therefore, the smaller Tafel slope of fcc-2H-fcc Au NR suggests a faster electron exchange rate (i.e. faster reaction kinetics).
First-principles calculations were conducted to further elucidate the reaction mechanism. Free energy diagrams for CO 2 RR towards the CO production were calculated for the three Au catalysts. The calculation models are shown in Supplementary  Fig. 18. Consistent with the Tafel analysis (Fig. 3e), the first-step reduction of CO 2 into *COOH is the rate-limiting step, as verified by the positive reaction free energy (ΔG *COOH ) for all studied surfaces (Fig. 3g). Note that the values of ΔG *COOH for 2H (110), 2H (110)/fcc (101) interface, 2H (1 10), fcc (100), fcc (101) and fcc (111) are 0.95, 1.02, 1.25, 1.06, 1.10 and 1.52 eV, respectively. In the next step, the produced *COOH reacts with a proton/electron pair to produce *CO and H 2 O, which was found to be exothermic. In general, 2H (110) surface and 2H (110)/fcc (101) interface possess lower ΔG *COOH compared to the representative fcc Au surfaces, indicating a lower kinetic barrier for CO 2 RR. Therefore, the excellent CO 2 RR catalytic performance of the heterophase fcc-2H-fcc Au NRs could be attributed to the presence of unusual 2H phase and the interface between 2H and fcc phases. A long-term durability test was carried out to evaluate the stability of fcc-2H-fcc Au NRs towards CO 2 RR. Notably, at the applied potential of −0.6 V (vs. RHE), the FE for CO production showed negligible decrease, and was maintained above 95% after an almost 24-h test (Fig. 3h). Meanwhile, the structure characterization of Au NRs after the durability test identified that the fcc-2H-fcc heterophase of Au NRs was well preserved, suggesting that the heterophase Au NRs possess outstanding electrochemical stability (Supplementary Fig. 19).

Discussion
In summary, we have realized a one-pot synthesis of well-defined heterophase fcc-2H-fcc Au NRs. This method overcomes the limitation of the classical seeded growth approach, which usually generates Au NRs with the thermodynamically stable fcc phase [47][48][49] . Compared with the common fcc Au NRs, the heterophase fcc-2H-fcc Au NRs demonstrate novel optical features like the energy shift and intensity modulation of surface plasmon polariton resonances. Furthermore, the electrocatalytic activity of heterophase fcc-2H-fcc Au NRs toward CO 2 RR outperforms their fcc counterparts due to the energetically more favourable adsorption of reaction intermediates on the 2H-phase surface and the heterophase 2H/fcc interface, which was confirmed by firstprinciples calculations. It is believed that the controlled synthesis of well-defined heterophase nanostructures not only provides an effective strategy to fine-tune the physicochemical properties and catalytic performances of noble metal nanomaterials, but also opens up an avenue for fundamentally understanding the crystallization behaviour of noble metals at the nanometer scale.
Synthesis of heterophase fcc-2H-fcc Au NRs. In a typical synthesis, 960 mg of ndodeylamine were first dissolved into 0.8 mL of oleylamine in a glass vial to form a homogeneous solution, denoted as Solution A. Then 3.0 mg of potassium tetrachloroaurate were added into the glass vial containing Solution A, followed by purging oxygen gas into the solution for about 5 min. After the glass vial was sealed, the mixture was sonicated for about 2 min. The resultant homogenous growth solution was then heated in an oil bath at 65°C for 17 h. After the reaction was completed, 2 mL of the mixture of cyclohexane and ethanol (v/v = 5/1) were added into the resulted solution, followed by sonication for about 1 min. The colloidal product was collected by centrifugation at 7,000 r.p.m. for 3 min, thoroughly washed with a solvent mixture of cyclohexane/ethanol (v/v = 5/1), and then re-dispersed into cyclohexane.
Synthesis of fcc Au NRs. The fcc Au NRs were synthesized according to a previously reported seeded growth method with a slight modification 50  Synthesis of fcc Au NPs. In a typical experiment, 6.12 mg of HAuCl 4 , 330 µL of oleylamine and 5.31 mL of hexane were successively added in a glass vial. Then the glass vial was carefully sealed and gently shaken to dissolve HAuCl 4 in the mixture of hexane and oleylamine. The resultant growth solution in the glass vial was heated at 58°C for 16 h. After the reaction was completed, the colloidal product was collected by centrifugation at 14,500 r.p.m. for 5 min, thoroughly washed with hexane, and then re-dispersed into cyclohexane.
Eigenmode characterization of individual heterophase fcc-2H-fcc Au NR. HRTEM images of individual fcc-2H-fcc Au NR were taken with a FEI Titan TEM with Schottky electron source, operated at 200 kV. The same instrument in STEM mode was used at 80 kV for the eigenmode characterization. A Wien-type monochromator dispersed the electron beam in energy, and a narrow energyselecting slit formed a monochrome electron beam with typical full-width at half maximum values of~80 meV. Using a convergence semi-angle of 13 mrad, ã 1 nm 2 electron probe was formed and scanned in rectangular areas of pixels (for mapping), or placed at specific locations to acquire individual spectra using a Gatan Tridiem ER EELS detector. All spectra were normalized at their maximum, i.e., the top of zero-loss peak (ZLP). The spectral background signal was corrected by fitting and subtracting a high-quality background spectrum that was measured elsewhere on the same sample without particles nearby.
DFT computation. The dielectric permittivities of fcc Au and 2H Au were computed using standard DFT approaches, i.e., independent particle approximation and the dipole-transition limit as implemented in density functional theory 51,52 WIEN2k package (see www.wien2k.at for details). In order to choose the density functional, we have benchmarked the computed dielectric permittivity for bulk fcc Au against measured permittivity in the low-energy region up to 4 eV. It turned out that the modified Becke-Johnson (MBJ) 53 exchange correlated potential developed by Tran and Blaha 54 in their latest parametrization 55 gives the closest match. This exchange correlation potential was subsequently applied to compute the response for Au in the 2H phase.
EELS simulation. The electron beam in EELS simulations was modelled as a monoenergetic flow of relativistic particles moving with constant velocity. This model is validated 56 by the average energy loss being negligible compared to the initial energy of electrons (80 keV), as well as by a small change of the beam momentum during its interaction with a thin sample. The electric E(ω, r) and magnetic H(ω, r) fields excited by the electron beam at frequency ω were numerically computed in COMSOL Multiphysics following the macroscopic Maxwell's equations, where ε 0 and μ 0 are electric and magnetic constants, respectively, and J(ω, r) is the current density of the electron beam. The space-dependent dielectric permittivity ε ω; r ð Þ was composed from the DFT-derived complex permittivities depending on the nanorod composition. Finally, the energy loss probability was classically computed though the ohmic loss 57 as follows, where the integration is performed over a volume of the electron beam (V).
Preparation of working electrodes. Typically, 6 mg of Vulcan XC-72R carbon black were mixed with 1 mL of catalyst solution, and the obtained solution was sonicated for 60 min. The solid was collected by centrifuging at 14,000 r.p.m. for 10 min, and re-dispersed into 1 mL of ethanol. Subsequently, 20 µL of Nafion solution were added into the solution containing catalyst and carbon black, and then the mixed solution was sonicated for another 15 min to form a well-dispersed catalyst ink. After that, 100 µL of the resultant catalyst ink were dropped onto both sides of a carbon paper (1.0 cm × 0.5 cm) with a working area of 1.0 cm 2 . The loading amount of Au catalysts was 266 µg. The obtained electrode was dried at room temperature.
Electrochemical CO 2 RR. CO 2 RR was performed in a gas-tight two-compartment H-type cell separated by an ion exchange membrane (Nafion 117). The electrolyte was CO 2 -saturated 0.5 M KHCO 3 . The H-type cell was connected to an electrochemical workstation (CHI 760E) using a three-electrode system. Platinum mesh and Ag/AgCl (with saturated KCl as the filling solution) were used as counter electrode and reference electrode, respectively. All electrode potentials were converted to the reversible hydrogen electrode (RHE) reference scale using the following equation: Before the test, the electrolyte in the cathodic compartment was bubbled with CO 2 gas for 30 min. During the test, the electrolyte in the cathodic compartment was stirred at 450 r.p.m. CO 2 gas was bubbled into the cathodic compartment at a rate of 20.0 standard cubic centimeters per minute (s.c.c.m.), and was routed into an on-line gas chromatograph (GC, Agilent 7890B). The GC analysis was set up to split the gas sample into three aliquots. One aliquot passed a thermal conductivity detector (TCD), and the other two aliquots were routed through two different flame ionization detectors (FIDs). FE calculation. The FE and partial current density of CO were calculated from GC chromatogram peak areas, where C CO is the volume concentration of CO based on the calibration of the GC. The equation is shown as follows.
where i total is the measured current, v is CO 2 flow rate, F is Faradaic constant, p 0 is pressure (1.013 bar), T is temperature (300 K), and R is the ideal gas constant (8.314 J mol K −1 ).
ECSA measurements. The ECSA of all the Au nanostructures were measured based on the charge associated with the stripping of an underpotential-deposited Cu monolayer. Briefly, after the working electrode was immersed into N 2 -saturated 0.5 M H 2 SO 4 solution, cyclic voltammograms (CVs) from 0.05 V to 0.50 V (vs. Ag/ AgCl) at a scan rate of 10 mV s −1 were acquired repeatedly until traces converged. Subsequently, the working electrode was immersed into a N 2 -saturated 0.5 M H 2 SO 4 solution containing 0.1 M CuSO 4 . Note that the electrolyte was continuously purged with N 2 . CVs were recorded from 0.05 to 0.50 V (vs. Ag/AgCl) at a scan rate of 10 mV s −1 . The anodic stripping waves were integrated. The factor used to convert the stripping charge to surface area was 0.42 mC cm −2 .
First-principles calculations. All the calculations were performed by using Vienna ab initio simulation package (VASP) 58,59 . The ion-electron interactions were described by the projector augmented wave (PAW) method 60 . The generalized gradient approximation in the revised Perdew-Burke-Ernzerhof (RPBE) functional and a cut-off energy of 400 eV for plane-wave basis set were adopted 61 . The convergence threshold was 10 −5 eV and 0.03 eV/Å for energy and force, respectively. To avoid the interaction between two periodic units, a vacuum space of 20 Å was used. The models were constructed from four layers of the corresponding Au facets, where the bottom two layers were fixed. The Brillouin zones were sampled by a Monkhorst-Pack k-point mesh with a 3 × 3 × 1 k-point grid. The free energy change (ΔG) of each elementary reaction was calculated as follows.
where E, E ZPE , T, S, C P and G sol are the electronic energy, zero-point energy, temperature, entropy, heat capacity and correction to the solvation effect on the free energy, respectively. All the corrections we used in this work were summarized in Supplementary Table 2.
Characterization. The TEM samples were prepared by directly dropping the suspension of different samples onto the full carbon-coated copper grids (200 mesh), followed by drying at ambient conditions. The TEM, HRTEM and HAADF-STEM images were taken on a JEOL JEM-2100F TEM operated at 200 kV. The aberration-corrected HAADF-STEM images were obtained by a JEOL ARM200F (JEOL, Tokyo, Japan) aberration-corrected transmission electron microscope operated at 200 kV with a cold field emission gun and double hexapole C s correctors (CEOS GmbH, Heidelberg, Germany). The XPS samples were prepared by directly dropping the suspension of different samples onto the Si substrates, followed by drying at ambient conditions. The XPS measurements were carried out on a VG ESCALAB 220i-XL instrument (base pressure <10 −5 mbar) equipped with a monochromatic Al Kα (1486.7 eV) X-ray source. The XRD characterization was conducted on an X-ray diffractometer (Shimadzu, XRD-6000) operated at 30 mA and 40 kV. The concentration of catalyst was estimated by the inductively coupled plasma-optical emission spectroscopy (ICP-OES, Dual-view Optima 5300 DV). The catalysts after the durability measurements were collected from the electrodes with sonication in ethanol and then used for the TEM characterization.

Data availability
The data that support the findings of this study are available from the corresponding author upon reasonable request.