Synthesis and characterization of PdPtRuCuNi/CNFsThe PdPtRuCuNi HEA nanoparticles were synthesized with acetylacetone salts as metal precursors (Supplementary Fig. 1)39. The atomic ratio of Pd, Pt, Ru, Cu, and Ni is nearly 21:13:20:24:22 (summarized in Supplementary Table 1), as measured by inductively coupled plasma optical emission spectroscopy (ICP-OES). The HEA particles are supported uniformly on CNFs and have an average size (dave) of approximately 17.5 ± 12.5 nm (Fig. 2a, b and supplementary Fig. 2a). The high-resolution transmission electron microscopy (HRTEM) suggests that the PdPtRuCuNi HEA particles are nanocrystalline with an interplanar spacing of about 2.20 Å (Fig. 2c and supplementary Fig. 2b). High-angle annular dark-field scanning transmission electron microscopy (HAADF-STEM) and energy-dispersive X-ray spectroscopy (EDS) mapping reveal that all elements are homogeneously distributed within the PdPtRuCuNi HEA nanoparticles (Fig. 2d and Supplementary Fig. 2c). The PdPtRuCuNi/CNFs exhibit two broad X-ray diffraction (XRD) peaks at 2θ of 41.4 ° and 47.4 °, corresponding to the (111) and (200) crystalline planes of PdPtRuCuNi HEA (Fig. 2e and Supplementary Table 3), indicating a single-phase alloy with a fcc structure. X-ray photoelectron spectroscopy (XPS) analysis shows that the metals exist mainly in their zero-valence states (Fig. 2f-j). A small number of metal oxides are present, likely originating from surface oxidation. Compared to their individual metals, the binding energies of Pd 3d, Pt 4 f, Ru 3p, and Cu 2p shift to the lower energies (Fig. 2f–i), while Ni 2p exhibits a slight shift to higher energies (Fig. 2j). This suggests that alloying these metals result in a strong electronic interaction within the PdPtRuCuNi HEA, which could change the surface electronic structure and potentially enhance its catalytic performance40. These results have also been demonstrated by X-ray absorption near edge structure (XANES) spectra (Supplementary Fig. 4a-c). The Pt, Cu, and Ni absorption energies for PdPtRuCuNi/CNFs are similar to those of metal foil, indicating the dominance of metallic-stated elements. The last half is slightly changed in shape and intensity, indicating the interaction among each element. The Fourier transforms of these EXAFS are displayed in Supplementary Fig. 4d-f. The distances of Me−Me in PdPtRuCuNi HEA are overall shorter than that of metal foils, and the line intensity of PdPtRuCuNi HEA is much lower than that of the metal foil, indicating the change in bond lengths and coordination numbers. The bond length (R) and coordination numbers are summarized in Supplementary Table 4.Fig. 2: Characterization of PdPtRuCuNi/CNFs.a, b The morphologies. c HRTEM images. d EDS mapping. e XRD pattern. f-j The high-resolution XPS spectra of PdPtRuCuNi HEA and monometal catalyst. f Pd 3d, g Pt 4 f, h Ru 3p, i Cu 2p, and j Ni 2p.Catalytic hydrogenation of alkynyl and phenyl groupsThe solid DEB molecule was selected as a representative reactant to evaluate the catalytic activity of PdPtRuCuNi/CNFs. They were mixed uniformly with PdPtRuCuNi/CNFs and Pd/CNFs to form loose composites, providing more accessible pathways for H2 diffusion (Supplementary Fig. 5 and 6). The hydrogen uptake tests were performed at different temperatures under ≤1 bar H2, recording the time and uptake capacity automatically. Both the PdPtRuCuNi/CNFs-DEB and Pd/CNFs-DEB composites exhibit increasing hydrogen uptake with reaction time at 25 °C (Fig. 3d and Supplementary Fig. 7a). Notably, the final hydrogen uptake capacities for PdPtRuCuNi/CNFs-DEB and Pd/CNFs composites are 762 cm3/g and 198 cm3/g, significantly surpassing those of the catalysts themselves without DEB ( ~ 2 cm3/g for PdPtRuCuNi/CNFs, and ~30 cm3/g for Pd/CNFs, Fig. 3a). This indicates that the hydrogen uptake is primarily attributed to the chemisorption (hydrogenation) of DEB molecules, affecting significantly by the element type. In addition, the commercial Pd/C and Ru/C are used as contrast catalysts, the hydrogen uptakes of Pd/C-DEB and Ru/C-DEB composites are 186 cm3/g and 0 cm3/g (Fig. 3b and supplementary Fig. 7), respectively, which are significantly lower than PdPtRuCuNi/CNFs, and Ru/C with no catalytic activity for DEB in mild condition.Fig. 3: The hydrogenation performances of DEB catalyzed by PdPtRuCuNi/CNFs and Pd/CNFs at 25 °C under ≤ 1 bar H2.a The final hydrogen uptake capacities of PdPtRuCuNi/CNFs-DEB and Pd/CNFs-DEB composites compared with catalysts themselves. b The final hydrogen uptake capacities of commercial Pd/C-DEB, Ru/C-DEB compared with PdPtRuCuNi/CNFs-DEB. c The hydrogen uptake curves with time at different temperature. d The curves of hydrogen uptake and reaction time of PdPtRuCuNi/CNFs-DEB and Pd/CNFs-DEB. e 13C NMR spectra of DEB before and after hydrogenation. f The conversion of C-C unsaturated bonds (alkynyl and phenyl groups) in the hydrogenated DEB molecules.We use the slopes of the hydrogen uptake curve to denote the reaction rate of DEB, and the linearly fitted rates of different stages are marked in Fig. 3c, d. Firstly, the hydrogen uptake curves of PdPtRuCuNi/CNFs-DEB were evaluated at various temperatures (25, 50, and 100 °C, Fig. 3c and supplementary Fig. 8). The results show that the similar hydrogen uptake values are obtained at 25 and 50 °C, while some decrease of hydrogen uptake at 10 °C because increasing the temperature promote the chemical reaction equilibrium moving in the opposite direction due to the exothermic hydrogenation reaction. Moreover, the reaction rates are elevated from 7.12 cm3/(g·h) to 25.6 cm3/(g·h) as the temperature from 25 °C to 100 °C in the initial hydrogenation process. However, the reaction rate falls sharply as the times increase for 50 °C and 100 °C, and the curves have a significant inflection point, which means the reaction orders and kinetics have changed. So PdPtRuCuNi has a high activity at 25 °C. Next, the Pd/CNFs-DEB as the control group to further elevate the reaction rate, Pd/CNFs-DEB reacts rapidly with H2, with a rate of about 12.85 cm3/(g·h), and then the reaction rate decays with time. In contrast, PdPtRuCuNi/CNFs-DEB displays a nearly constant reaction rate of ~7.12 cm3/(g·h) until reaching its maximum hydrogen uptake. This suggests that the PdPtRuCuNi HEA catalyst demonstrates excellent catalytic performance in DEB hydrogenation and operates through a different catalytic mechanism with Pd/CNFs. This difference is likely caused by the synergistic effects of multiple elements in HEA during the solventless hydrogenation process.The hydrogenated products are analyzed using NMR spectra. As shown in Fig. 3e, f and Supplementary Fig. 7c, the Pd/CNFs catalyst only hydrogenates the alkynyl group in the DEB molecules (88% conversion of all the alkynyl groups), forming 1,4-bis(phenylethyl)benzene (8H-DEB). For the PdPtRuCuNi/CNFs, they display outstanding catalytic activity towards both alkynyl and phenyl groups, achieving complete hydrogenation and forming 1,4-bis(ethylcyclohexane)cyclohexane (26H-DEB). This transformation brings almost 100% conversion of the C-C unsaturated bonds in DEB catalyzed by PdPtRuCuNi/CNFs, while achieving ~27% conversion by Pd/CNFs, approximately three times increased, which agrees well with the hydrogen uptake capacities. Moreover, the PdPtRuCuNi/CNFs-DEB composite evolves from powders into melts after hydrogenation, indirectly indicating the deep hydrogenation of DEB molecules (Supplementary Fig. 6). Overall, the utilization of a novel PdPtRuCuNi HEA nanocatalyst has led to a breakthrough in achieving the simultaneous and complete hydrogenation of alkynyl and phenyl groups under ≤1 bar hydrogen pressure, marking a notable advancement in this field.Synergistic catalytic effect of PdPtRuCuNiTo delve deeper into the synergistic effect of the PdPtRuCuNi catalyst and clarify the roles of each element in the hydrogenation of DEB, we have synthesized a series of quaternary nanocatalysts with similar loadings (supplementary table 5) by removing one element at a time (i.e., PtRuCuNi/CNFs, PdRuCuNi/CNFs, PdPtCuNi/CNFs, PdPtRuNi/CNFs, and PdPtRuCu/CNFs), and their phase structures, chemical compositions, and morphologies have been characterized in detail (Supplementary Fig. 9-19). The catalytic performances of these catalysts for the alkynyl and phenyl groups in the DEB molecules have been investigated and the reaction rates and hydrogenated products are analyzed following the same procedures stated above (Fig. 4 and Supplementary Fig. 20). The reduction in reaction rate occurs in descending order as the Pd, Pt, Ru, Cu, and Ni are sequentially removed, with the largest decrease occurring when Pd is removed and the smallest when Ni is removed (Fig. 4a and Supplementary Fig. 20a). It is noteworthy that the reaction rates of these quaternary catalysts are all lower than that of the PdPtRuCuNi, indicating that all five elements work synergetically to regulate the reaction kinetics of DEB. Importantly, when any one of the Pd, Pt, or Ru elements are absent, the composites display not only decreased reaction rates but also a reduction in final hydrogen uptake capacity and conversion of DEB (Fig. 4b and Supplementary Fig. 20b). In contrast, when these three essential elements—Pd, Pt, and Ru—coexist, the composites exhibit robust hydrogen uptake performance, achieving nearly 100% DEB conversion. These findings have also been demonstrated by PdPtRu/CNFs catalyst (Supplementary Fig. 21 and 22), indicating that the Pd, Pt, and Ru elements play a vital role in the hydrogenation of DEB.Fig. 4: The synergistic catalytic effect of PdPtRuCuNi HEA.a The reaction rates of the DEB-contained composites catalyzed by quaternary and quinary catalysts. b The overall conversion of C-C unsaturated bonds in DEB and the relative conversion of alkynyl (pink) and phenyl (blue) groups. c 13C NMR spectra of DEB before and after hydrogenation. d The yields of different hydrogenated DEB products.A detailed comparison from NMR analysis shows that the difference in the total conversion of these quaternary catalysts mainly stems from the hydrogenation of phenyl groups (Fig. 4c and Supplementary Fig. 20c). Specifically, under the catalysis of PtRuCuNi and PdRuCuNi, the primary products are 8H-DEB, with only alkynyl groups hydrogenated. When using PdPtCuNi as the catalyst, the products transform into a mixture of 1-ethylcyclohexane-4-phenylethyl benzene and 1,4-bis(ethylcyclohexane)benzene, termed 14H- and 20H-DEB, with one or two sides of phenyl groups further hydrogenation. When Pd, Pt, and Ru coexist, the product becomes 26H-DEB, with alkynyl and phenyl groups undergoing complete hydrogenation (Supplementary Fig. 20d). Thus, we deduce that alkynyl groups are most readily hydrogenated, then one or two sides of phenyl groups, and finally the middle phenyl group. This hierarchy in the hydrogenation process of DEB can be attributed to the different activation energies associated with phenyl groups at different sites.Understanding the hydrogenation process over PdPtRuCuNi HEAWe proceed to gain insights into the synergistic effect of elements within the PdPtRuCuNi HEA in the context of hydrogenation. As the activation of both H2 and aromatic molecules are critical for hydrogenation, we utilize density functional theory (DFT) calculations to unveil their interactions with the close-packed (111) surfaces of PdPtRuCuNi HEA and Pd, Pt, Ru, Cu, Ni single metals. To save computation costs, we employ phenylacetylene (C6H5–C ≡ CH), possessing phenyl and alkynyl groups, as a model molecule of DEB (abbreviated as PhA hereafter). In the later section, we also verify that the conversion of PhA is ~98.5%, close to 100%. As to the favored adsorption sites, we consider both hollow and bridge sites for H2 and H species (Ha). For PhA, previous studies, verified by our tests, indicate that the optimal adsorption configuration involves the phenyl group sitting above a 3-fold hollow site (πph), with the alkynyl group forming di-μ bonding41. Furthermore, we’ve constructed six random slabs for the PdPtRuCuNi HEA and explored all possible adsorption sites on the slabs to evaluate the impact of local atomic arrangement on its catalytic activity (Methods).We start with the adsorption of H2 molecules. The adsorption energies (denoted as ∆Eads-H2) and the optimized H-H distance for 480 adsorption sites (192 hollow and 288 bridge sites, see Methods) on the PdPtRuCuNi slabs are shown in Supplementary Fig. 26. All H2 molecules are activated on the PdPtRuCuNi surface, with H-H bonds elongated (dH-H > 0.74 Å). At ~47.5% of the adsorption sites, the H2 molecule dissociates, with a final H-H distance > 1.48 Å (twice the equilibrium H-H bond length) and Ha as the product. This dissociation also pulls the Ha closer to the HEA surface. In comparison, the H2 dissociates completely on the Pd catalyst’s surface (Supplementary Table 6). Thus, we infer that the slight reduction in the reaction rate of PdPtRuCuNi/CNFs-DEB compared with Pd/CNFs-DEB at the initial stage in Fig. 3a might be attributed to an insufficient Ha supply under low H2 pressure. Supplementary Fig. 27 illustrates the spatial variation of ∆Eads-H2 on the six random slabs, where site-to-site heterogeneity is observed, suggesting a large variety on a single HEA surface. We selected an active site near Pd from Slab 6 and derived the H2 dissociation energy via the climbing image nudged elastic band method (CI-NEB)42. This site displays the lowest H2 dissociation energy compared to any other single metal, despite all the metals except Cu having relatively low dissociation energies (<0.1 eV) (Supplementary Fig. 28a).The migration of Ha from the catalyst to the reactant is crucial for hydrogenation, significantly affecting the reaction rate and ultimate hydrogenation conversion. We calculate the binding energies of the dissociated Ha with the adsorption sites on Slab 6 of PdPtRuCuNi (Fig. 5a). Compared with single metals, the PdPtRuCuNi HEA display a broader distribution of binding energies. About 86.5% and 40.6% of the sites possess lower binding energies than those on the hollow and bridge sites of the pure Pd surface, due to the weak binding of Pt and Cu with Ha. It indicates that most Ha bind weakly with PdPtRuCuNi and are easy to transfer from the catalyst to the reactant. Our temperature-programmed desorption (TPD) data validates these DFT calculation results, as evidenced by the broader desorption peak of H2 and a ~ 7.8 °C decrease in the maximum desorption peak of PdPtRuCuNi/CNFs compared to that of Pd/CNFs (Supplementary Fig. 28b).Fig. 5: DFT calculation of H2 and phenylacetylene over the PdPtRuCuNi (111) surface.a Binding energies of dissociated H species at the hollow and bridge sites of a representative PdPtRuCuNi slab (upper panel) and at the hollow (with subscript “h”) and bridge (with subscript “b”) sites of Pd, Pt, Ru, Cu, Ni slabs (lower panel). b The absorption energy of PhA over the (111) surface of PdPtRuCuNi (above) and Pd, Pt, Ru, Cu, Ni (below) catalysts. c The first-step hydrogenation reaction barriers of the alkynyl and phenyl groups in PhA over PdPtRuCuNi (111). d The first-step hydrogenation reaction barriers of the phenyl groups in phenylethane over PdPtRuCuNi (111), Pd (111), and Ru (111). e Reaction pathways for the hydrogenation of DEB over PdPtRuCuNi HEA nanocatalysts.Next, we examine the interaction of PhA with the HEA surface. Figure 5b displays the adsorption energy distribution of PhA (∆Eads-PhA) on the HEA and monometallic slabs. In the monometallic case, the ∆Eads-PhA on Pd, Pt, Ru, Cu, and Ni surfaces are −3.13, −2.66, −4.99, −0.88, and −3.78 eV, respectively. After alloying into PdPtRuCuNi, the ∆Eads-PhA exhibits a widened distribution from −1.47 to −4.41 eV, reminiscent of a Gaussian shape with its peak near the Pd catalyst. Supplementary Fig. 29 gives the spatial distribution of ∆Eads-PhA on the slabs as in Supplementary Fig. 27. This diversity can effectively regulate the adsorption energies (or behaviors) of PhA and its intermediates, creating an undulating landscape for them to adsorb and desorb29. Supplementary Fig. 30a presents the increase in C-C bond distances within the alkynyl and phenyl groups after activation by Pd and PdPtRuCuNi HEA catalysts. We see the C-C distances in the two groups increase, with the elongation within the alkynyl group being more pronounced. The alkynyl group is also dragged closer to the PdPtRuCuNi HEA surface compared to the phenyl group (Supplementary Fig. 30b). It is reasonable to consider that the alkynyl group is more prone to hydrogenation than the phenyl group. To validate this conjecture, we calculated the first-step hydrogenation reaction barriers for the alkynyl and phenyl groups in PhA over the PdPtRuCuNi surface (Fig. 5c). Only 0.13 eV is needed for the alkynyl group to add the first Ha, while the phenyl group demands 0.43 eV, corroborating the experimental observations in Fig. 3d and Fig. 4b. Consequently, the phenyl hydrogenation is a thorny problem for achieving complete hydrogenation. Since the alkynyl group is more prone to hydrogenation and all catalysts can hydrogenate the alkynyl group, we proceed to utilize the product after alkynyl group hydrogenation, namely phenylethane (C6H5–CH2-CH3), as a model molecule to investigate the subsequent hydrogenation of the phenyl group. Similarly, its first-step hydrogenation reaction barriers on PdPtRuCuNi, Pd, and Ru surfaces are calculated (Fig. 5d). The carbon in the para-position (Supplementary Fig. 31) has the lowest hydrogenation reaction barrier. We find that the first-step hydrogenation reaction barriers of phenylethane on the Pd and Ru are 0.58 and 0.56 eV, while on PdPtRuCuNi, 0.35 eV. These findings further indicate that PdPtRuCuNi can be more favorable for hydrogenating the phenyl group than the Pd and Ru catalysts, explaining why PdPtRuCuNi can achieve complete hydrogenation of DEB.Based on the above experiment and calculation results, we propose the following reaction mechanism (Fig. 5e): Initially, H2 is activated by the PdPtRuCuNi/CNFs HEA catalyst, dissociating into active Ha. Due to the low binding energy and undulating energy landscape on the surfaces, Ha can migrate to other sites through diffusion or hydrogen spillover (i). Meanwhile, the DEB molecules are adsorbed and activated on the PdPtRuCuNi HEA, especially around the active sites containing Pd, Pt, and Ru atoms (ii). Then, the Ha meets the activated DEB molecules, and the alkynyl groups are preferentially hydrogenated (iii). Next, one or two sides of phenyl groups in DEB undergo hydrogenation, followed by the middle one, until all phenyl groups are hydrogenated (iv). Finally, the hydrogenated products detach from the PdPtRuCuNi HEA surface. In addition to the reduced reaction barriers for hydrogenation, the existence of diverse and multi-functional sites within the HEA surface enables that the H2 activation, reactant absorption, and hydrogenation steps can occur separately at different sites, making the hydrogenation more efficient.Universal hydrogenation of diverse aromatic derivativesTo assess the versatility of our PdPtRuCuNi HEA catalyst in the hydrogenation of aromatics, we conduct experiments with different aromatic substrates containing phenyl groups. These investigations were conducted at 25 °C in a solvent-free state, with 10 bar H2 pressure to speed up the reaction. The resulting products were analyzed by NMR, and the conversion was derived (Supplementary Figs. 32–36). Remarkably, nearly 100% reactivity is achieved for toluene, biphenyl, and diphenylacetylene, 98.6% for styrene, and 98.5% for PhA (Fig. 1b). All C-C unsaturated bonds, including ene/alkyne and phenyl groups in different liquid or solid reactants, can be fully hydrogenated at ambient temperature. Notably, these hydrogenations occur with no additional solvents, promoters, or energy inputs. The final hydrogenated products are saturated cyclanes, regardless of the presence of substituted groups around the benzene ring (alkyl/ene/alkyne or phenyl), the location of phenyl groups (terminal or internal), or whether the catalysis occurs in a liquid or solid phase. This performance significantly expands the applications of HEA nanocatalysts in the catalytic hydrogenation of various aromatic derivatives.To highlight the performance of our PdPtRuCuNi HEA catalyst in the hydrogenation of aromatics, we first use toluene as an example and summarize the hydrogenation conditions of toluene catalyzed by different catalysts in Fig. 6a. Compared with Pd, Pt, Ru, and Rh-contained catalysts17,18,19,20,21,43,44,45,46, we achieve a complete conversion and a 100% yield of methylcyclohexane under mild conditions (25 °C, 10 bar H2), without any solvent or promoter. Although similar reaction conditions (30 °C, 10 bar H2) have been reported for benzene hydrogenation over Ru/α-Al2O3 or Ru-B/MIL-53 (B = Al, Cr, and AlCr)43,45, these reactions require alcohol solvents such as isopropanol or ethanol (Supplementary Table 7). Moreover, after being used four times, our PdPtRuCuNi HEA catalyst still maintains 100% conversion for toluene, underscoring its stability (Supplementary Fig. 32c and d). Figure 6b and Supplementary Table 7 further give the temperature-pressure conditions for hydrogenating benzene, toluene, styrene, PhA, biphenyl, phenol, nitrobenzene, and DEB. We see that our reaction temperature is the lowest to achieve complete hydrogenation. As stated above, solid organic hydrogen getters are usually employed to eliminate the hydrogen in confined spaces, ensuring hydrogen safety. Figure 6c and Supplementary Table 8 compare the hydrogen uptake capacity of DEB-containing composites using Pd and Pd-based multifunctional catalysts. Our PdPtRu-contained nanocatalysts exhibit the highest catalytic activity for DEB than other Pd catalysts22,23,24,47,48.Fig. 6: Catalytic hydrogenation of aromatic substrates over PdPtRuCuNi HEA catalyst at 25 °C.a Catalytic hydrogenation of toluene with different catalysts with or without solvent. b The comparison of the hydrogenation conditions of the present study with the reported results in the literature (Supplementary Table 7). c The comparison of hydrogen uptake capacities of DEB-contained composites catalyzed by Pd and PdPtRuCuNi HEA catalysts (Supplementary Table 8). d The extensive application of aromatic derivatives hydrogenated at ambient conditions.In summary, our HEA catalyst has a universal ability for ambient hydrogenation of diverse aromatic derivatives, offering significant potential applications in petrochemical refining, pharmaceuticals, chemical industries, and hydrogen safety and storage (Fig. 6d). Importantly, hydrogenation at ambient conditions and low pressure requires no specific equipment or additional energy, allowing for minimal energy consumption and cost-effectiveness, offering advantages for fine chemical products and organic synthesis.
