Catalysis-tunable Heusler alloys in selective hydrogenation of alkynes: A new potential for old materials

See allHide authors and affiliations

Science Advances  19 Oct 2018:
Vol. 4, no. 10, eaat6063
DOI: 10.1126/sciadv.aat6063


Heusler alloys (X2YZ) are well-established intermetallic compound materials in various fields because their function can be precisely adjusted by elemental substitution (e.g., X2YZ1−xZx). Although intermetallic compound catalysts started attracting attention recently, catalysis researchers are not familiar with Heusler alloys. We report their potential as novel catalysts focusing on the selective hydrogenation of alkynes. We found that Co2MnGe and Co2FeGe alloys have great alkene selectivity. Mutual substitution of Mn and Fe (Co2MnxFe1−xGe) enhanced the reaction rate without changing selectivity. The substitution of Ga for Ge decreased the selectivity but increased the reaction rate monotonically with Ga composition. Elucidation of these mechanisms revealed that the fine tuning of catalytic properties is possible in Heusler alloys by separately using ligand and ensemble effects of elemental substitution.


Heusler alloys (X2YZ) are intermetallic compounds that were discovered by Heusler in 1903 (1). The crystal structure of the L21 Heusler alloy is illustrated in Fig. 1 (2). The Cu2MnAl Heusler alloy and other Heusler alloys have attracted attention because they exhibit ferromagnetism even without ferromagnetic elements. They are now particularly attractive in various research fields, including spintronics (3), thermoelectrics (4), and ferromagnetic shape memory alloys (5). Their features include having many possible sets of elements, as shown in Fig. 1 [the details are reported by Yin et al. (6)], and allowing partial substitution of components (e.g., X2YZ1–xZx), which is used to optimize electronic structures. The first-mentioned feature portends the discovery of novel catalysts, and the second one would enable the fine tuning of catalytic properties because catalysis is governed by the electronic structure (the ligand effect) and by the atomic arrangement at the surface (the ensemble effect).

Fig. 1 Concept of this study.

Effects of elemental substitution are schematically illustrated with the Heusler alloy’s crystal structure (L21-type X2YZ) (2). Periodic table represents typical elements consisting of X, Y, and Z (color painted). Other elements such as lanthanoids and those in groups 1 and 2 are also possible components of Heusler alloys [see (6)].

Intermetallic compounds are attracting much attention as novel catalysts because they form new electronic structures that differ from those of their components and because they have specific surface atomic arrangements (7, 8). To the best of our knowledge, there was no paper regarding Heusler alloys dealing with catalysis until our study. Previously, to represent their potential for new catalysts, we tested the potential of various Heusler alloys for hydrogenation of propyne and oxidation of carbon monoxide (9). Since then, we have found that Co2MnGe and Co2FeGe are very useful for selective hydrogenation of alkynes and now report that here.

Selective hydrogenation is an important process for removing alkyne impurities from alkene feedstock because the impurities poison the catalyst during alkene polymerization (10, 11). Palladium-based catalysts have been used in industry for selective hydrogenation of ethyne (acetylene) in ethene (ethylene) and propyne (methylacetylene) in propene (propylene) (1012). Replacement of Pd with a nonprecious metal is desirable in terms of cost performance and resource availability, and these catalysts have been developed in the laboratory using intermetallic compounds such as NiZn, NiZn3 (13), Al13Fe4 (14), Ni3Ga, and Ni3Sn2 (15). In such binary intermetallics, however, elemental substitution is restricted while keeping the original crystal structure. In contrast, this substitution with a wide composition range of various elements can be done in Heusler alloy catalysts. This would be a unique advantage enabling precise adjustment of the catalytic properties to the target reaction. Here, we demonstrate control of the catalytic properties of Heusler alloys by elemental substitution (Co2MnxFe1−xGayGe1−y) in the selective hydrogenation of alkynes (propyne and ethyne).


Hydrogenation of alkyne in the presence of alkene

Figure 2 shows the alkyne conversions and alkene selectivities in the hydrogenations of propyne (C3H4) in the presence of propene (C3H6), and ethyne (C2H2) in the presence of ethene (C2H4) for powder samples sieved to a particle size of 20 to 63 μm. Ordinary catalysts such as pure metals usually exhibit high alkene selectivity when alkyne conversion is low, and selectivity decreases as conversion increases (13, 16, 17) as it does for Co2FeGa, as shown in Fig. 2F. This behavior is attributed to strongly adsorbed alkyne molecules inhibiting the adsorption of alkene molecules, while a decrease in inhibiters (alkyne) at a high conversion enables alkene adsorption.

Fig. 2 Results of alkyne hydrogenation in the presence of alkene.

Alkyne conversions (circles) and alkene selectivities (triangles) for (A) Co2FeGe, (B) Co2Mn0.5Fe0.5Ge, (C) Co2MnGe, (D) Co2FeGa0.25Ge0.75, (E) Co2FeGa0.5Ge0.5, and (F) Co2FeGa. Closed and open symbols represent results for hydrogenations of C3H4 in C3H6 and C2H2 in C2H4, respectively. Reactants were [0.1% alkynes/10% alkenes/40% H2/He balance]. Selectivities in low-conversion region are not shown due to large errors.

In contrast, Co2MnxFe1−xGe showed high alkene selectivities even at a 100% alkyne conversion (Fig. 2, A to C). This is surprising because the hydrogen/alkyne ratios in the feeds were 400 (0.1% alkyne/10% alkene/40% H2/He balance), which is much higher than the reported ratios (10 to 20) (10, 1315). In addition, Co2MnxFe1−xGe poorly catalyzed the hydrogenation of C3H6 in the absence of C3H4 (0.1% C3H6/40% H2/He balance), which indicates that it has almost no ability to hydrogenate alkenes (fig. S3). Since alkene selectivity usually decreases with the hydrogen/alkyne ratio (10, 18), operating conditions in industrial plants must be precisely controlled to avoid hydrogenating the alkene feedstock and to maximize the removal rate of alkyne (10, 11). Therefore, Co2MnxFe1−xGe is expected to be useful under any operating conditions.

The Mn substitution for Fe resulted in certain changes in the conversions and selectivities (Fig. 2, A to C). The Ga substitution for Ge increased the conversions at lower temperatures and decreased the selectivities at higher temperatures (Fig. 2, D to F).

Hydrogenation of C3H4 in the absence of C3H6

To determine the effects of elemental substitution, C3H4 hydrogenation in the absence of C3H6 was investigated, excluding the complexities due to the presence of alkene molecules and the oligomerization that often occurs in C2H2 hydrogenation (10, 11, 18). Figure 3 shows the effects of Mn substitution for Fe (Co2MnxFe1−xGe; Fig. 3A) and Ga substitution for Ge (Co2FeGayGe1−y; Fig. 3B) on the reaction rates of C3H4 normalized to the values for Co2FeGe and on C3H6 selectivities with a C3H4 conversion of 100% at 200°C (overall data are shown in fig. S1). Similar to the results shown in Fig. 2, Co2MnxFe1−xGe showed high C3H6 selectivity, and Ga substitution increased the rate but decreased the selectivity. The Mn substitution did not change the selectivity but increased the rate in the partial substitution range (0.25 ≤ x ≤ 0.75). These tendencies were also observed for Co2MnxFe1−xGa and Co2MnGayGe1−y (fig. S2), indicating that the Mn-Fe substitution and the Ga-Ge substitution act independently of each other.

Fig. 3 Effects of elemental substitution in C3H4 hydrogenation in the absence of C3H6.

Changes in reaction rates of C3H4 per surface area and selectivities for C3H6 due to (A) Mn substitution (Co2MnxFe1−xGe) and (B) Ga substitution (Co2FeGayGe1−y). The reaction rates were normalized by those of Co2FeGe (x = y = 0). In (A), the rates at 100°C (diamonds), 75°C (squares), and 50°C (circles) are shown. Blue triangles represent the selectivities for 100% C3H4 conversion at 200°C. The reactant was [0.1% C3H4/40% H2/He balance].

Relationship between electronic structure and activation energy

Considering the features of Heusler alloys, the substitution effects most likely originated from the electronic structure. The d-band theory combined with the Brønsted-Evans-Polanyi (BEP) relation is widely accepted to explain the relationship between catalytic activity and electronic structure (1921). This theory suggests a linear relationship between the activation energy (Ea) and the mean energy of the d-band (εd, so-called d-band center). Figure 4 (A1 and A2) shows the densities of states (DOSs) calculated using density functional theory; they reveal continuous band shifts with Mn and Ga compositions. This behavior was verified by hard x-ray photoelectron spectroscopy (HAXPES) of the valence bands, as shown in Fig. 4 (B1 and B2), which indicated that the electronic structures were properly controlled by elemental substitutions in actual samples. The high intensities in the high binding energy region were due to large photoionization cross sections of the s-bands under this measurement condition (22). Figure 4C shows the Ea obtained experimentally for the C3H4 hydrogenation and the εd estimated from the partial DOSs of the d-bands. Their changes with substitution were similar. The Ea versus εd plot in Fig. 4D shows a linear relationship, indicating that changes in the electronic structure contributed to changes in the catalytic properties.

Fig. 4 Electronic structures and activation energy for C3H4 hydrogenation in the absence of C3H6.

(A1 and A2) DOSs. f.u., formula unit. (B1 and B2) Valence band HAXPES spectra. (C) Changes in experimental Ea and calculated εd due to Mn and Ga substitutions. (D) Relationship between Ea and εd.


Validity of relationship between Ea and εd

To confirm the validity of the changes in Ea due to εd, we look for insights into the reaction mechanism. The reaction orders with respect to C3H4 and H2 were roughly estimated to be −0.2 to 0 and 0.8 to 1, respectively (fig. S5 and table S1). These values are usually obtained for pure transition metal catalysts (17, 23). The adsorption processes of C3H4 and H2 were considered to be almost noncompetitive because the C3H4 reaction order is typically −1 or −2 for competitive adsorption in the Langmuir-Hinshelwood (LH) model. In our case, H2 (g) ⇒ 2H (ad) or C3H4 (ad) + 2H (ad) ⇒ C3H6 (ad) is likely to be the rate-determining step corresponding to the reaction orders obtained, and following equations were derived under the assumption of the LH model with noncompetitive adsorption (see the Supplementary Materials) (23).

[H2 (g) ⇒ 2H (ad) is rate determining]Embedded Image(1)

[C3H4 (ad) + 2H (ad) => C3H6 (ad) is rate determining]Embedded Image(2)where r, k, PH2, and KH2,ad are the reaction rate, the rate constant, the partial pressure of H2, and the equilibrium constant for dissociative adsorption of H2, respectively. The k is expanded to the Arrhenius equation using the activation energy of the rate-determining step (Ea,H2,ad or Ea,s) and a frequency factor (A). The KH2,ad is expanded to two exponential parts containing the standard entropy (ΔS°H2,ad) and enthalpy (ΔH°H2,ad) of the H2 adsorption. Accordingly, the Ea value obtained experimentally is described by Ea = Ea,s + ΔH°H2,ad, when r follows Eq. 2.

The d-band theory and the BEP relation indicate that increase in εd decreases (negatively increases) the H2 adsorption energy and eventually decreases Ea,H2,ad. Thus, Eq. 1 agrees with the trend in Fig. 4D. In the case of Eq. 2, increase in εd strengthens adsorption of C and H atoms, which probably increases Ea,s as reported on the hydrogenation of C2Hx molecules (24). Thus, the increase in εd should decrease ΔH°H2,ad more than the increase in Ea,s to produce the trend in Fig. 4D. According to the literature (20, 25, 26), Ea,s = 114 + 13εd [kJ mol−1] and ΔH°H2,ad = −106 − 31εd [kJ mol−1] (unit of εd, eV) can be obtained, although the Ea,s versus εd relationship is rough because they were extrapolated from only four data points in two reports (20, 25). The slope of Ea,s + ΔH°H2,ad = 8 − 18εd agrees well with the slope of the experimental values in Fig. 4D. Therefore, both Eqs. 1 and 2 qualitatively explain the experimental results. However, it is hard to conclude which model is more reasonable, because there is uncertainty regarding the reaction mechanism (e.g., eq. S3) and the adsorption configurations of C3H4 and H2 for Heusler alloys differ from those reported for pure metals (discussed below).

Effects of Mn substitution

Although uncertainty remains, the reaction rate would basically follow Eq. 1 or Eq. 2. That is, the rate should monotonically increase with Mn composition (x) since the Ea value decreased with x, but the rate turned downward at x > 0.5, as shown in Fig. 3A. This behavior is attributed to the enthalpy-entropy compensation effect often observed in various thermodynamic phenomena including catalytic reactions (2729). Although the mechanism of this effect is still not fully understood, taking the H2 adsorption as an example, we can qualitatively understand that strengthening the binding between the catalyst and H atoms decreases (negatively increases) ΔH°H2,ad and also decreases ΔS°H2,ad since a stronger binding reduces the degree of freedom of adsorbed H atoms. Therefore, a decrease in ΔS°H2,ad due to the compensation effect might suppress an increase in the rate if the rate followed Eq. 2. According to the transition state theory (27), A ∝ exp(ΔSTS/R) and exp(−Ea/RT) ∝ exp(−ΔHTS/RT) are derived, where ΔSTS and ΔHTS are changes in entropy and enthalpy in the transition state, respectively. Thus, the compensation could happen even if the rate follows Eq. 1.

In many cases of this compensation, the changes in enthalpy and entropy show a linear relationship, dH = TdS, corresponding to full compensation (28). In Fig. 3A, the compensation for x = 1 was apparently almost full, while that for 0.25 ≤ x ≤ 0.75 was not full. This behavior is not strange given that a nonlinear relationship between dH and dS is sometimes observed (29). However, the reason for this behavior was not elucidated in this study due to various complexities, including the uncertainty of the reaction model and the compensation between Ea,s and ΔSTS in the case of Eq. 2. We suppose that the behavior is due to a random distribution of Mn and Fe atoms at Y sites (of X2YZ), which would cause the interactions between adsorbed species to be inhomogeneous, because the interaction between adsorbates affects kinetics and thermodynamics on the surface such that the enthalpy and entropy of adsorption depend on the adsorbate coverage (30, 31).

Effects of Ga substitution

In contrast to Mn substitution, Ga substitution increased the reaction rate monotonically with Ga composition (y) along with the decrease in Ea (Figs. 3B and 4C), which seems reasonable. However, why was there no compensation for the rate increase? Furthermore, why did the selectivity decrease significantly? We attribute these results to the same mechanism discussed below. In the hydrogenation of C3H6 in the absence of C3H4, the C3H6 conversion was very low for the Co2FeGe, while it monotonically increased with Ga substitution (fig. S3). This indicates that Ga substitution brought the ability to adsorb C3H6.

In selective hydrogenation of C2H2 using Pd-based catalysts, the high C2H4 selectivity is attributed to a small adsorption site surrounded by poisonous species that sterically hinder the adsorption of larger molecules, C2H4, but allow the adsorption of smaller molecules, C2H2 (10, 11). Figure 5A shows schematic illustrations of close-packed (110) surfaces for Co2FeGe, Co2FeGa0.5Ge0.5, and Co2FeGa. Assuming that Ge atoms do not have the ability to adsorb hydrocarbons while Ga atoms do, the size of sites available for adsorption is small on Co2FeGe(110), while it is enlarged by Ga substitution as we can see on Co2FeGa0.5Ge0.5(110). We believe that this is the mechanism of the high alkene selectivity for Co2MnxFe1−xGe and the reduction in selectivity due to Ga substitution. Under this assumption, the increase in the rate due to Ga substitution can also be explained by an increase in the number of sites available for adsorption.

Fig. 5 Effects of Ga substitution for Ge of Co2FeGe.

(A) Schematic illustration of (110) surfaces for Co2FeGe, Co2FeGa0.5Ge0.5, and Co2FeGa. (B) Local partial DOSs (LPDOSs) of Ge at Co2FeGe(110) (blue) and Ga at Co2FeGa(110) (red) surfaces.

This hypothesis is supported by theoretical calculation. A series of theoretical studies by Krajčí and Hafner revealed that typical element atoms such as Al, Ga, and Zn contribute to the adsorption of hydrocarbons on the surface of intermetallic compounds with transition metals (3234). Taking the surface of Al13Co4 as an example, we see that C2H2 molecules adsorb on a bridge site between two Al atoms rather than on a site containing Co atoms (32, 33). In terms of the electronic structure, the hydrocarbon adsorption ability of typical elements can be related to the hybridization of the p-bands of the typical elements with the d-bands of the transition metals (32). Figure 5B shows the p-bands of Ge in Co2FeGe and of Ga in Co2FeGa for (110) surfaces. We can see that the Ga p-band is located at a shallower level from the Fermi level than the Ge p-band. The Ga atoms have a higher hydrocarbon adsorption ability than Ge atoms, similar to the d-band theory. Therefore, the assumption of differing adsorption abilities between Ga and Ge is reasonable.


Here, we have presented the advantage of using Heusler alloys as new catalysts and have demonstrated the control of catalytic properties by elemental substitution in the selective hydrogenation of alkyne. The Co2MnxFe1−xGe catalysts showed high alkene selectivity even with 100% alkyne conversion and even with a very high H2/alkyne ratio of 400; they poorly catalyzed the C3H6 hydrogenation in the absence of C3H4. These results indicate that the Co2MnxFe1−xGe catalysts can be useful under any operating condition.

The elemental substitution resulted in a change in Ea along with εd, which indicates that the electronic structure was accurately controlled and that it contributed to the catalytic properties. The significant increase in reaction rate and decrease in selectivity with Ga substitution for Ge were attributed to the difference in hydrocarbon adsorption ability between Ga and Ge atoms in Co2FeGayGe1−y. This seems to be an ensemble effect in a conventional manner, in contrast to the effect of Mn-Fe substitution changing only the electronic structure, the so-called ligand effect. This means that Heusler alloys can separately use the ensemble and ligand effects depending on the purpose. Although Mn substitution for Fe did not greatly increase the reaction rate with the decrease in Ea probably due to the compensation effect, this ligand effect could be effective for other reactions in which a correspondence between an activity and an electronic structure-derived parameter (e.g., Ea, εd, and adsorption energy) has been experimentally observed, such as methanation from CO and an oxygen reduction reaction (21, 35).

In addition, Heusler alloys can be used as a well-defined platform for revealing the catalysis of intermetallic compounds because catalysis with various electronic structures and various surface elements can be investigated using elemental substitution under the same crystal structure. Last, since many Heusler alloys are free of precious metals, as shown in Fig. 1, the future development of new practical catalysts is promising.


Sample preparation and characterization

Heusler alloy ingots were prepared from pure metallic sources (purity, >99.9%) by arc melting followed by annealing during which the condition was typically 1000°C for 72 hours followed by 600°C for 72 hours and then 500°C for 72 hours under an Ar atmosphere. Powder catalysts were obtained by crushing the ingots using a mortar and pestle, followed by sieving to a particle size of 20 to 63 μm.

The structural properties were evaluated by x-ray diffraction (XRD) using powders with a size of <20 μm after removing the strain and defects by annealing at 600°C for 1 hour under an H2 atmosphere. The XRD was performed using the Bragg-Brentano geometry with Cu Kα radiation (Rigaku, Ultima IV diffractometer). The XRD patterns shown in fig. S6 indicate the formation of a (nearly) single phase L21 structure for all samples, although negligible unknown peaks were detected in several samples. The degree of atomic ordering was evaluated using Webster’s model with factors S and α. The S factor corresponds to a long-range order parameter in a binary alloy. S = 1 ensures that all Co atoms occupy the X sites in X2YZ, that is, the ordering between the Co and YZ elements. S = 0 means that all atoms randomly occupy the X, Y, and Z sites. α describes the “disordering” between the Y and Z elements. α = 0.5 means no Y-Z order, while α = 0 means no Y-Z “disorder.” Thus, satisfaction of both S = 1 and α = 0 means perfect L21 ordering. These factors were estimated by usingEmbedded Image(3)Embedded Image(4)where I200, I111, and I400 are the integrated intensities of the 200 and 111 superlattice peaks and the 400 fundamental peak, respectively, and Ifund is the fundamental peak of 220, 400, 422, or 440 diffraction. The numerators are experimental values, while the denominators are theoretical values in the perfect-order cases calculated from atomic scattering factors including anomalous scattering terms, the multiplicity factor, the Lorentz-polarization factor, and the Debye-Waller factors. The obtained S and α factors are listed in table S1. Sufficiently high S and low α were obtained, indicating high ordering into the L21 structure. Kojima et al. (9) for more information about the estimation of S and α.

The surface area of the catalyst after the C3H4 hydrogenation in the absence of C3H6 (table S1) was estimated using the Brunauer-Emmett-Teller method with Kr adsorption (MicrotracBEL, BELSORP-max volumetric adsorption instrument) to calculate the reaction rates shown in Fig. 3 and fig. S2.

Catalytic measurements and evaluation

The catalytic properties were evaluated for the reactions using [0.1% C3H4/40% H2/He balance], [0.1% C3H6/40% H2/He balance], [0.1% C3H4/10% C3H6/40% H2/He balance], and [0.1% C2H2/10% C2H4/40% H2/He balance] reactants. The measurements were conducted in standard flow reactors with gas chromatographs: (i) an Agilent 490 Micro GC equipped with a thermal conductivity detector (TCD) and a PoraPLOT Q column for the C3H4 hydrogenation and (ii) a Shimadzu GC-8A equipped with a TCD and a Shincarbon-ST column for the C2H2 hydrogenation. The catalyst (400 mg) was supported on quartz wool in a quartz tube with an internal diameter of 4 mm and surrounded by an electric furnace. After the catalyst was heated under an H2 gas flow at 600°C for 1 hour to remove the surface oxides, the reactant mixture was introduced at 30 ml min−1 (standard temperature and pressure) at ambient temperature and pressure through mass flow controllers into the catalyst channel (space velocity, about 20,000 hours−1). The catalyst was then left standing for 1 hour, and then, analysis of the gaseous species and heating were started. The unreacted reactants and products were analyzed 30 min after heating ended at every 25°C interval from 25° to 250°C.

The alkyne conversion and the alkene selectivity were estimated by usingEmbedded Image(5)Embedded Image(6)Embedded Image(7)where Cfeed, Cunreact, Calkene, and Calkane were the concentrations of the feed alkyne, the unreacted alkyne, the produced alkene, and the produced alkane, respectively. Clost in Eq. 6 is the concentration of the carbon species lost due probably to oligomerization (Clost = CfeedCunreactCalkeneCalkane). For the alkyne hydrogenation in the presence of alkene, the alkene selectivity was estimated using Eq. 7 under the assumption that all the reacted alkynes were converted into alkenes and that Clost = 0 because the very large amount of the feed alkene made it impossible to estimate small amounts of produced alkene and the lost carbon.

The apparent Ea for the C3H4 hydrogenation in the absence of C3H6 was evaluated by continuous cyclic measurement with analysis performed at every 1°C during the heating and cooling cycle at a rate of 0.5°C min−1 between 50° and 100°C. Appropriate amounts of the catalysts and the flow rates were used so that the conversion was appropriate for estimating the Ea value. The Ea value was estimated by averaging the values for both the heating and cooling processes after the change in the slope of the Arrhenius plot had sufficiently settled, as shown in fig. S4. The values and the ranges of the temperature and the conversion used for the estimation are shown in table S1. The error bar corresponds to the original values before the averaging.

The reaction orders with respect to C3H4 and H2 for the C3H4 hydrogenation in the absence of C3H6 were evaluated at 100°C from the reaction rates for four different concentrations in the feed. The reactants were [n% C3H4/40% H2/He balance (n = 0.05, 0.08, 0.13, and 0.2)] and [0.1% C3H4/m% H2/He balance (m = 20, 30, 50, and 80)]. Appropriate amounts of the catalysts and the flow rates were used so that the conversion was appropriate for estimating the reaction orders. Data were obtained 30 min after changing the concentration. Figure S5 shows the double logarithmic plots of the reaction rate versus concentration. The reaction orders were measured twice; their averages are shown in table S1.

Calculation of electronic structures

The electronic structures were calculated using the full-potential linearized augmented plane wave method (36). The generalized gradient approximation developed by Perdew et al. (37) was used for the exchange-correlation potential. The plane wave cutoff was RKmax = 7.0, where R is the smallest atomic sphere radius and Kmax is the magnitude of the largest K vector. For the atomic sphere radius, we used 2.22 atomic unit (a.u.) for the 3d transition metals and 2.09 a.u. for Ga and Ge. To calculate the electronic structure of each Heusler alloy listed in table S2, we used the model structure with the space group given in the table. The optimized lattice constants are also listed in the table. The lattice vectors for Co2Mn0.5Fe0.5Ge and Co2FeGe0.5Ga0.5 are a′ = (a + b)/2, b′ = (−a + b)/2, and c′ = c, where a, b, and c are lattice vectors for the cubic lattice. We used a 20 by 20 by 20 K-mesh for the Brillouin zone integration for ternary alloys such as Co2FeGe. For quaternary alloys, the K-meshes were chosen so that the size of the K-mesh corresponds to that for the ternary alloys.

The most stable surface of Heusler alloys seems to be a close-packed (110) surface, which has been indicated by calculation for seven types of surfaces, as shown in table S3. To simulate a thin (110) film of Co2FeGe and Co2FeGa used for Fig. 5B, we considered a slab in a supercell specified by three primitive vectors: a′ = (−a + b)/2 = a(−0.5,0.5,0), b′ = c = a(0,0,1), and c′ = 7(a + b) = a(7,7,0), where a denotes the length of the side of the conventional unit cell of a cubic lattice. The film used was constructed on 11 layers with a total thickness of 5| c′|/14 and an empty (vacuum) region with a thickness of 9| c′|/14. We used a = 5.741 and 5.720 Å, which gave a minimum total energy for the bulk ferromagnetic Co2FeGe and Co2FeGa, respectively (see table S2). For both films, we estimated the equilibrium positions of the atoms in the direction of c′. The K-mesh adopted was 14 by 10 by 2. See (38) for more information.

The d-band center (εd) is defined as (20)Embedded Image(8)where D(ε) is the PDOS of the d-band at an energy, ε. We estimated the εd values from the sum of D(ε) for the transition metal components with the integration range from −13.000 to 6.457 eV, which was enough for a relative comparison of εd. The εd values were estimated from the bulk electronic structures because the calculation of the surface electronic structures for Co2Mn0.75Fe0.25Ge, Co2Mn0.25Fe0.75Ge, Co2FeGa0.25Ge0.75, and Co2FeGa0.75Ge0.25 made it difficult to assume appropriate supercells and slabs and would have very high calculation costs. The difference in DOS between the surface and the bulk originates from symmetry breaking at the surface, which indicates that materials with the same crystal structure and similar component elements make similar changes in the DOS from the bulk to the surface (39). Thus, the bulk εd values could be used for relative comparisons among our samples.

Hard x-ray photoelectron spectroscopy

HAXPES measurements at room temperature were performed to evaluate the electronic structures experimentally at the BL15XU undulator beamline (40) of SPring-8. The sample ingots were cut into about 1-mm by 1-mm by 10-mm rods, fixed on a copper holder using a silver paste, and then placed in a preparation chamber (base pressure, 2.0 × 10−7 Pa) connected with a measurement chamber (base pressure, <1.0 × 10−7 Pa). After evacuation of the preparation chamber, the rods were fractured in situ to obtain a clean surface and then moved into the measurement chamber. The excitation photon energy and total energy resolution were set to 5.95 keV and 240 meV, respectively. The binding energy was referred to the Fermi level of an evaporated gold film. Details of the experimental setup are described elsewhere (40). The electric field vector of the incident beam was parallel to the photoelectron trajectory to the analyzer. Because of this geometry, a high-intensity structure originating from the s-bands at binding energies >4 eV is evident in Fig. 4 (B1 and B2) (22). We also performed the measurement for Co2FeGa0.5Ge0.5 powders before and after the alkyne hydrogenation in the presence of alkene, which indicated no surface segregation and no carbon deposition (coking or oil production) within the sensitivity, as shown in fig. S7.


Supplementary material for this article is available at

Supplementary Text

Fig. S1. Results of C3H4 hydrogenation in the absence of C3H6.

Fig. S2. Details of elemental substitution effects in C3H4 hydrogenation in the absence of C3H6.

Fig. S3. Results of C3H6 hydrogenation.

Fig. S4. An example Arrhenius plot for Co2FeGa in C3H4 hydrogenation in the absence of C3H6.

Fig. S5. C3H4 reaction rate at 100°C versus feed concentrations in C3H4 hydrogenation in the absence of C3H6.

Fig. S6. XRD patterns for Heusler alloy samples.

Fig. S7. Core-level HAXPES spectra for Co2FeGa0.5Ge0.5 powders before and after alkyne hydrogenation in the presence of alkene.

Table S1. List of Heusler alloy samples and related information.

Table S2. Information on calculation of electronic structures.

Table S3. List for calculation on surface stability.

This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.


Acknowledgments: We are grateful to K. Nozawa of Kagoshima University for helpful discussion regarding the relation between the electronic structure and catalysis and to H. Iwamoto and R. Tsukuda of Tohoku University for help in the HAXPES measurement performed at BL15XU at SPring-8 under the approval of NIMS Synchrotron X-ray Station (proposal no. 2017B4905). Funding: A part of this work was supported by the Hattori Hokokai Foundation, the Iwatani Naoji Foundation, the Noguchi Institute, and NIMS microstructural characterization platform as a program of “Nanotechnology Platform” (project no. 12024046) of the Ministry of Education, Culture, Sports, Science and Technology (MEXT), Japan. We also thank the support from Dynamic Alliance for Open Innovation Bridging Human, Environment and Materials and the Frontier Research Institute for Interdisciplinary Sciences of Tohoku University. Author contributions: T.K. designed this research, performed all experiments and their analyses, and wrote the manuscript. S.K. and A.-P.T. supported and discussed the experiments and the analyses. S.F. performed the calculation of electronic structure and wrote the method of calculation in the manuscript. S.U. performed the HAXPES measurement and its analysis together with T.K. All authors discussed the results and supported the manuscript writing. Competing interests: T.K., S.K., and A.-P.T. are authors on a patent application related to this work filed by Tohoku University (JP patent application no. 2017–220616, filed on 16 November 2017). All other authors declare that they have no competing interests. Data and materials availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from the authors.
View Abstract

Stay Connected to Science Advances

Navigate This Article