# Dynamics of Pd Subsurface Hydride Formation and Their Impact on the Selectivity Control for Selective Butadiene Hydrogenation Reaction

Esther Asedegbega-Nieto, Ana Iglesias-Juez, Marco Di Michiel, Marcos Fernandez-Garcia, Inmaculada Rodriguez-Ramos, Antonio Guerrero-Ruiz

► To cite this version:

Esther Asedegbega-Nieto, Ana Iglesias-Juez, Marco Di Michiel, Marcos Fernandez-Garcia, Inmaculada Rodriguez-Ramos, et al. Dynamics of Pd Subsurface Hydride Formation and Their Impact on the Selectivity Control for Selective Butadiene Hydrogenation Reaction. Nanomaterials, 2023, 13 (6), pp.1099. 10.3390/nano13061099. hal-05077986

HAL Id: hal-05077986

https://hal.science/hal-05077986v1

Submitted on 22 May 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# Dynamics of Pd Subsurface Hydride Formation and Their Impact on the Selectivity Control for Selective Butadiene Hydrogenation Reaction

Esther Asedegbega-Nieto <fnref n="1" />,* , Ana Iglesias-Juez <fnref n="2" />, Marco Di Michiel <fnref n="3" />, Marcos Fernandez-García <fnref n="2" /> ID,
Inmaculada Rodriguez-Ramos <fnref n="2" /> ID and Antonio Guerrero-Ruiz <fnref n="1" />,* ID

**Abstract:** Structure-sensitive catalyzed reactions can be influenced by a number of parameters. So far, it has been established that the formation of Pd-C species is responsible for the behavior of Pd nanoparticles employed as catalysts in a butadiene partial hydrogenation reaction. In this study, we introduce some experimental evidence indicating that subsurface Pd hydride species are governing the reactivity of this reaction. In particular, we detect that the extent of formation/decomposition of PdHx species is very sensitive to the Pd nanoparticle aggregate dimensions, and this finally controls the selectivity in this process. The main and direct methodology applied to determine this reaction mechanism step is time-resolved high-energy X-ray diffraction (HEXRD).

**Keywords:** palladium catalysts; subsurface hydride; butadiene; partial hydrogenation; HEXRD

<figure><img src="image_3.png" /></figure>

check for
updates

**Citation:** Asedegbega-Nieto, E.; Iglesias-Juez, A.; Di Michiel, M.; Fernandez-Garcia, M.; Rodriguez-Ramos, I.; Guerrero-Ruiz, A. Dynamics of Pd Subsurface Hydride Formation and Their Impact on the Selectivity Control for Selective Butadiene Hydrogenation Reaction. *Nanomaterials* **2023**, *13*, 1099. https://doi.org/10.3390/nano13061099

Academic Editor: Genqiang Zhang

Received: 20 February 2023

Revised: 10 March 2023

Accepted: 16 March 2023

Published: 19 March 2023

<figure><img src="image_4.png" /></figure>

**Copyright:** © 2023 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article
distributed under the terms and
conditions of the Creative Commons
Attribution (CC BY) license (https://
creativecommons.org/licenses/by/
4.0/).

## 1. Introduction

When studying heterogeneous catalysts, we encounter the concept of “structure sensitive reactions”. This term was first introduced by M. Boudart in the 1970s [1]. Briefly, for supported metal nanoparticles, some catalytic parameters, either specific activities for a given surface center (for example, turnover frequency) and/or the obtained selectivity values, vary drastically with the particle size.

The exact definition of catalytic surface sites is unclear due to the fact that, for instance, sometimes subsurface species from the bulk crystallites can participate in the reaction mechanism. For example, it has been reported that monoatomic dispersed Pd atoms over Cu surfaces constitute an excellent catalyst for the hydrogenation of acetylene to ethylene [2]. On the other hand, in the definition of a catalytic surface site, a discrete number (greater than one) of surface atoms can participate. These are the so-called ensemble structures or ensemble size effect, where the probability to find surface sites with the right size in terms of adjacent active centers governs the reaction rate [3]. In the surface site, a special arrangement of the metallic sites is considered to be involved in the definition of catalytic center. For example, it is assumed that B5-type surface sites, which consist of the arrangement of a layer of three atoms and another directly above of two atoms at the monoatomic step on a Ru (0001) terrace, are responsible for the activity of ruthenium catalysts [4]. Furthermore, J. A. Somorjai carried out studies where the relationship between reaction kinetics and atomic surface structure was questioned. This author deduced that structure sensitivity and insensitivity could be explained by the dynamic restructuring of metal surface induced by the presence of adsorbates [5]. More recently, a review dealing with the mechanistic and atomic-level insights into semihydrogenation reactions of alkynes and alkadienes was published [6]. In this sense, these authors classify under three categories, at the atomic

<footnote n="1">¹ Dpto. Química Inorgánica y Técnica, Facultad de Ciencias, UNED, Av. de Esparta s/n,
28232 Las Rozas, Madrid, Spain</footnote>

<footnote n="2">² Instituto de Catálisis y Petroleoquímica, CSIC, c/Marie Curie No. 2, Cantoblanco, 28049 Madrid, Spain</footnote>

<footnote n="3">³ ESRF—The European Synchrotron, 71 Avenue des Martyrs, 38000 Grenoble, France</footnote>

<footnote>* Correspondence: easedegbega@ccia.uned.es (E.A.-N.); aguerrero@ccia.uned.es (A.G.-R.)</footnote>
level, the regulations for the active sites: site isolation, local environment regulation, and oxygen vacancy and interfacial sites.

Particularly, the industrially relevant studies of the selective hydrogenation of triple bonds/or conjugated double bonds in hydrocarbons are often catalyzed by Pd nanocrystal-lites. Although this noble metal is known to show the highest activities and selectivities in these reactions [7,8], to date, the mechanism behind its behavior has not been established. In general, this group of partial hydrogenation reactions has been classified as structure sensitive [9], and in earlier studies [10], the differences in catalytic selectivity have been associated with the strength of chemisorption of reactant or product molecules. According to a study carried out by D. Teschner et al. [11], the formation of the Pd/C surface phase was directly linked to the nature of the reactive molecule. These authors who studied the hydrogenation of alkynes stated that the formed Pd/C species alter the surface chemistry by decreasing the amount of activated subsurface hydrogen, thereby controlling the reactivity of the Pd catalyst. In another study, involving the hydrogenation of 1,3-butadiene on Pd (110)-oriented surfaces, it was observed that the presence of both weakly π-bonded 1,3-butadiene and subsurface hydrogen atoms is required to achieve catalytic activity [12]. Furthermore, Wu et al. [13] also attributes enhancements in catalytic performance to subsurface chemistry. In this cited work, focused on the partial hydrogenation of acetylene over Pd catalysts, the authors claim that the main reason for the catalytic performance of the active site can be tuned based on the distinctive electronic and geometric structures of Pd derived as a result of the presence of the corresponding subsurface heteroatoms.

In this contribution, for the first time, we demonstrate how the participation of hydride species formed inside Pd crystallites can control the activity and selectivity and establish the most favorable conditions for optimum results. We show how the nature and response characteristics of PdHx species are metallic particle size-dependent. Further, we present an operando study of the reversible transformation of Pd to PdHx during the hydrogenation reaction of 1,3-butadiene.

## 2. Materials and Methods

### 2.1. *Synthesis of Materials*

All catalysts studied in this work are composed of carbon-based materials as support and Pd nanoparticles as active phase. The former is either commercial high-surface-area graphite (G) purchased from Timcal (SBET = 296 m²/g) or lab-prepared graphenic materials. The graphenic materials were synthesized starting from commercial natural graphite obtained from Alfa Aesar, which is oxidized to produce the corresponding graphite oxide via the modified Brodie method, as has been detailed in our previous work [14]. This can be briefly described as follows: 200 mL of fuming nitric acid was introduced in a round-bottom flask. Thereafter, 10 g of natural graphite (grain size of 100 mesh) was added. The mixture is kept at 0 °C and under agitation. Thus, 80 g of KClO₃ is slowly added over 1 h, after which the reaction is maintained at the abovementioned temperature for 21 h. Once completed, the obtained sample is filtered and washed until neutral pH is observed and then dried overnight in a vacuum furnace at 60 °C.

The graphite oxide (GO) prepared in the above paragraph is employed for the synthesis of two graphenic materials: GOE and GONE. GOE stands for exfoliated graphene oxide and is obtained by introducing the graphite oxide in a furnace and heating until 500 °C under inert atmosphere. As for GONE, which denotes N-doped exfoliated graphene oxide, this support is prepared by a solid-state synthesis method employing urea as the N precursor. For this, urea and GO are intimately mixed and ground, introduced in a furnace, and heated until 450 °C under nitrogen flow. Further details are specified in a previous work [15]. The three mentioned supports (G, GOE, and GONE) were impregnated with an ethanol–aqueous (1:1) solution of tetraamine palladium chloride. Quantities and concentrations were adjusted in order to obtain metal loadings of 1 and 2%. The prepared catalysts were dried overnight at 110 °C.
## 2.2. Characterization of Samples by TEM

Transmission Electronic Microscopy (TEM) was used to estimate the average particle size of Pd crystallites employing the formula $d_{Pd} = \sum n_i d_i / \sum n_i$. For this purpose, measurements were performed on a JEOL JEM-2100 field-emission gun electron microscope operating at 200 kV.

## 2.3. Catalytic Tests

Time-resolved high-energy X-ray diffraction (and mass spectrometry (MS) were the techniques employed, at the ID15A beamline of the European Synchrotron (ESRF), to monitor the whole process, which consisted of two main stages. Firstly, temperature-programmed reduction (TPR) is carried out on the catalyst in order to study the reducibility and formation of Pd species in the process of heating under a hydrogen atmosphere until 450 °C and cooling down to room temperature. This stage is necessary prior to the catalytic test as Pd has to be in its reduced state. The reductive gas is composed of 15% hydrogen in helium and a total flow rate of 20 mL/min. Secondly, the catalytic test consisting of the partial hydrogenation of butadiene was performed. For this, samples were pelleted in a size range of 0.075 to 0.150 mm and introduced in a 2 mm ID quartz capillary tube, which served as a fixed-bed continuous flow reactor. Glass wool stoppers are inserted at the top and bottom of the catalyst bed. Total flow rate is the same as for the TPR process, although in this case, butadiene (BD) is introduced in the reaction feed, maintaining a ratio to hydrogen of 1:5.

As for the HEXRD acquisition parameters, two sets of conditions were employed. On one hand, a wavelength of 0.177 Å (70 keV) was used and samples were placed at a distance of 65 cm from the Perkin Elmer area detector where data acquisition time was 30 s. The second set of experiments, named fast, on the other hand, were carried out at a distance of 105 cm from the detector, and the registered wavelength was 0.237 Å (52.33 keV) with exposure time reduced to 1 s. Both sets of experiments include the initial TPR and its consequent cooling to room temperature, and the reaction study, which consists of the introduction of the reaction mixture for 30 min at a constant temperature (switch 1) followed by a switch (switch 2) to reductive mixture (hydrogen diluted in helium) and thereafter increasing the temperature for the next reaction test at a higher temperature. This process is repeated for a number of different temperatures. The exhaust products from the reactor were analyzed using a European Spectrometry ecoSyst-P Man-Portable mass spectrometer with capillary inlet and heated inlet tubes. This serves to estimate the catalytic conversion of the reaction as it registers the characteristic mass fragments of the reactant, butadiene ($m/z$ = 39), as well as the possible products, butene (partially hydrogenated, $m/z$ = 41) and butane (completely hydrogenated, $m/z$ = 43).

# 3. Results and Discussion

As expected, upon TPR, a diffraction peak at $q \sim 2.8$ Å⁻¹ identified as Pd (111) reflection demonstrates the presence of reduced Pd in all studied samples. Similar values have been reported elsewhere for this noble metal [16,17]. Figure 1 depicts the behavior for sample 1PdGOE. All other diffractogram sets are available in the Supplementary Materials.

The evolution of PdHx species (diffraction peak at $q \sim 2.7$ Å⁻¹) formed during the cooling step, under a hydrogen atmosphere, is quite different in each case and is intimately linked with the Pd particle sizes. In summary, samples of smallest nanoparticle sizes (1PdG) apparently had no PdHx phase formation (see Figure S1a,b), which is in agreement with the literature, where it has been said that below 2.5 nm, no hydride is formed [18]. In the case of the catalyst with largest particle sizes (1PdGOE), there was an almost complete transformation of Pd into PdHx on cooling the catalyst (Figure 1). This observation is consistent with a recent report, where large quantities of hydrogen absorbed on Pd and formation of palladium hydride are observed, at room temperature and under ambient pressure [19]. Surprisingly enough, the formed PdHx was consumed when put in contact with the reaction mixture (total flow 20 mL/min composed of H₂ and butadiene (ratio
5:1) diluted in helium) at room temperature, and an evolution to metallic Pd occurred during the reaction (see Figure S2). This could explain the difference in conversion and selectivity values observed when comparing sample 1PdG (where no PdHx was formed) and 1PDGOE, as reported in Table 1.

<figure><img src="image_5.png" /><figcaption>**Figure 1.** Formation of metallic Pd during TPR process for 1PdGOE (insert) and formation of PdHx upon cooling. Black initial ($t_0$) and red final ($t_f$) patterns.</figcaption></figure>

<table><caption>**Table 1.** Particle size of the Pd nanoparticles and steady-state conversions and selectivities at 50 °C reaction temperature.</caption><thead><tr><td>Samples</td><td>Particle Size nm (TEM)</td><td>Conversion (%)</td><td>Selectivity to Butenes (%)</td></tr></thead><tbody><tr><td>1PdG</td><td>3.2</td><td>89</td><td>27</td></tr><tr><td>2PdG</td><td>4.9</td><td>89</td><td>27</td></tr><tr><td>1PdGONE</td><td>8.9</td><td>15</td><td>88</td></tr><tr><td>1PdGOE</td><td>9.7</td><td>8.9</td><td>92</td></tr></tbody></table>

In general, for Pd supported on graphenic or graphitic materials, an increase in the average Pd crystallite sizes, ranging from 3 to 10 nm, provokes a higher selectivity to butenes (from 27% to 92%). Hence, selectivity towards butenes is intimately linked with PdHx species present in this sample (1PdGOE). This tendency would be in agreement with studies carried out by Bauer et al. [20] where smaller Pd nanoparticles (1–2 nm) offered complete hydrogenation of propyne to propane at high conversion values, while high selectivity towards the partially hydrogenated product propene was favored with larger nanoparticles (about 15 nm). In this referenced study, the presence of subsurface C seemed to play a key role in the product distribution orientation. Yet, at that stage, transformation of PdHx to Pd was not considered or discussed as we are introducing in the present research. Using Pd supported on activated carbon catalysts, working in the hydrogenation of acetylene at 100 °C, both carbide and hydride were detected [21]. In another study by the same authors, XANES was said to discriminate between both phases [22]. Furthermore, Liu et al. [23] followed the formation of subsurface carbon on carbon-supported Pd catalysts (employing XAS and in situ XRD) at different temperatures for the partial hydrogenation of acetylene. It was concluded that the presence of this subsurface carbon improved the selectivity of the catalyst towards the unsaturated product, ethylene. The reason behind this would be that the formation of this species (Pd-carbide)
hinders the formation of Pd hydride species and weakens the ethylene adsorption that gives rise to further hydrogenation into the saturated hydrocarbon. This was particularly relevant for reaction temperatures above 100 °C and especially at 250 °C, where the highest selectivity values towards ethylene were obtained. Hence, this behavior was said to be enhanced on increasing the reaction temperature where the penetration of carbon into Pd was more likely to occur. Nevertheless, at lower temperatures, such as 50 °C, as is in our case, the formation of the Pd-carbide phase is reduced, and the formation of the Pd-hydride is promoted. In summary, under our experimental conditions, the formation of Pd carbide cannot be excluded based only on the HEXRD results. Nevertheless, we could discard it for two main reasons. On one hand, the presence of this phase would yield methane, which was not detected among the products formed. Additionally, taking into consideration the low reaction temperature and the stoichiometric conversion of the peak into a metallic contribution during the hydrogenation process, the presence of PdHx is confirmed.

In order to have a deeper insight into the kinetic transformation of Pd to PdHx (necessary for the optimum performance of catalysts in this reaction) and vice versa, further experiments were performed. In this second set of experiments, some adjustments were made. On the one hand, and so as not to exclude the effect of support, graphite was still studied, although in this case, with higher Pd loading (2%) in order to provide larger particle sizes and discriminate whether formation of PdHx is dependent only on particle size and independent of the support. Acquisition time was reduced to 1 s to increase time resolution and experiments were performed with incident X-ray irradiation wavelength of 0.237 Å and a larger sample detector distance in order to improve the peak resolution (as detailed in the Materials and Methods section).

Therefore, three samples were studied: 2PdG, 1PdGOE, and 1PdGONE. All three introduced different and surprising discoveries on the formation of the Pd species. For 1PdGOE (Figure 1), both Pd and PdHx are formed during TPR and cooling to room temperature, respectively. Relevantly, for the 1PdGONE sample, the shape of the diffraction peak at 2.7 Å⁻¹ suggests the presence of more than one contribution of PdHx, as can be viewed in Figure 2 (see also Supplementary Materials where both (1PdGOE and 1PdGONE) are compared (Figure S10)). These two contributions correspond to two well-defined species: PdHx α (x ≤ 0.015) and PdHx β (0.5 ≤ x ≤ 0.75) [24]. This was most visible for 1PdGONE (Figure 2) but not for 2PdG, where only one contribution of Pd and PdHx species was observed (see also Supplementary Materials). It has been reported that the miscibility gap between both phases (α and β) is narrowed with decreasing particle size [24]. Hence, for 2PdG, there is just one phase, possibly α-hydride, because it is well established that α-PdHx is preferentially formed over the smaller Pd crystallites (see comparative average particle sizes for these two samples, 1PdGONE and 2PdG, in Table 1). As can also be seen in Figure 2, the temporal evolution of the two Pd hydride peaks is different, indicating a difference in stability as a function of temperature.

The effect of hydrogen concentration, pressure, and temperature was revealed on the formation of Pd-hydride species (alpha and beta) [25,26]. These parameters govern the adsorption (physical and chemical) and diffusion processes involved in the formation of Pd-H species. The behavior of Pd thin films was compared with that of Pd NPs. As expected, lattice contraction is higher in NP with respect to bulk. This would justify the slightly lower contraction observed for 2PdG in comparison to the two others (1PdGOE and 1PdGONE), as it has been described elsewhere that lattice contractions increase with decreasing particle size for systems such as Au [27]. At the same time, in previous referenced work [25], on decreasing the temperature, a transition from alpha to beta is observed. Thus, we can also observe this in the cooling step of our TPR experiment for 1PDGONE (see Figure S7).
<figure><img src="image_6.png" /><figcaption>**Figure 2.** Transformation of Pd to PdHx during cooling to RT for 1PdGONE. Black initial ($t_0$), red final ($t_f$).</figcaption></figure>

The transformation of PdHx to Pd during the hydrogenation reaction of butadiene is illustrated in Figure 3. When introducing the reaction mixture at a constant temperature of 30 °C, the Pd peak grows at the expense of the PdHx peak. Nonetheless, Figure 3 reveals that the transformation is not complete. On the other hand, PdHx is recovered after switching from the reaction mixture back to H₂. The rate of transformation was also increased with increasing reaction temperature (30–70 °C) (see Figures S5, S8, and S12). This reversible transformation was estimated to be about 80% (achieved after 4 min into the reaction) in the case of 1PdGOE in the whole range of reaction temperatures. Figure 4 and the Supplementary Materials offer a more detailed scheme of this transformation for all the samples.

<figure><img src="image_7.png" /><figcaption>**Figure 3.** Transformation of PdHx to Pd during butadiene hydrogenation at 30 °C for 1PdGOE. Figure based on data detailed in Supplementary Materials section.</figcaption></figure>
<figure><img src="image_8.png" /><figcaption>**Figure 4.** Final steady-state HEXRD patterns (selected Pd(111) region) after H₂ exposure at 50 °C, subsequent switching to reaction mixture (Bd + H₂), and back switching to H₂ again for Pd samples: (a) 1PdGOE, (b) 1PdGONE, and (c) 2PdG.</figcaption></figure>

As for the 1PdGONE sample, subsequent switching to a reaction mixture (Bd + H₂) and switching back to H₂ again also leads to reversible transformation of PdHx into Pd, although to a lesser extent (Figure S8). In this case, the transformation achieved was about 66% after 4 min in the reaction (Figure 4).

As can be observed in Figure 4 (view also Supplementary Materials Figures S8 and S10), on switching from H₂ atmosphere to reaction mixture, both phases of PdHx are not equally reduced. The β phase is preferentially and rapidly consumed. At this point, it can be deduced that the hydrogen content of the hydride is a determining factor for the optimal performance of Pd nanoparticles. The results presented in Figure 4 prove that the participation of the PdHx α phase is lower for 1PdGOE than for 1PdGONE and its concentration is inversely proportional to the selectivity towards butenes (Table 1).

Once more, the catalyst 2PdG exhibited the most extreme behavior. Firstly, at the lowest reaction temperature, 30 °C, PdHx transformation into Pd is negligible. At an intermediate temperature, 50 °C, this transformation is much slower than for the other two catalysts under the same conditions, reaching only about 20% throughout the reaction (Figure 4). At a higher temperature, 70 °C, although the transformation rate was still lower than with the other catalysts, a 100% reversible transformation after 20 min was reached (Figure S12). These results are once more related to the particle size of the Pd catalysts, which, as have been summarized in Table 1, follow the order: 1PdGOE > 1PdGONE > 2PdG.

At the same time, the evolution of the reaction products was followed by a mass spectrometer, so both butadiene conversion and selectivity towards butane and butenes were simultaneously analyzed with the Pd evolution. The results (Table 1) prove that the characteristics of the Pd phase (see Supplementary Materials) and its evolution when switching reactants from H₂ to the reaction C₄H₆/H₂ mixture depend on the carbon support, which defines Pd nanoparticle size and, thus, the exposed PdHx phase and the active reaction temperature.

1PdGOE exhibits low conversion values, which increase slightly with increasing temperature. Referring to selectivity, the highest values towards butenes were observed with this sample, leaving butane production lower than 10% in all circumstances. In the case of 1PdGONE, the conversion was also low, although higher than that of 1PdGOE. It was also selective towards butenes, while butane production increased with increasing reaction
temperatures. 2PdG gave high conversion values and low selectivity towards butenes, proving that butane formation is favored when small Pd nanoparticles are obtained, as no beta hydride phase is formed [28].

Therefore, time-resolved analysis of the combined in situ HEXRD and MS data shows that the reversible formation decomposition of the PdHx, which is dependent on particle size, is intimately linked with the conversion and the selectivity for a hydrogenation reaction. We conclude that the β-hydride phase is the species responsible for the production of butenes, while also suppressing the butane formation.

Taking into account the importance of the formed PdHx and its reversible transformation, it was necessary to take a deeper look into the subject. For this, a kinetic study at different reaction temperatures of the PdHx transformation during the hydrogenation reaction was performed. This was carried out employing the Avrami model, which is very suitable for phase transformations [29] (see Supplementary Materials for data and analysis description). Table 2 summarizes the obtained transformation rate constants, k, associated with the apparent activation energy, while *n* is related to the critical rate-limiting step [30]. The value of k was highest for the sample of largest particle size 1PdGOE and lowest for that of the smallest, 2PdG. The exponential growth index, *n*, is particularly high for 2PdG. It is accepted that values of *n* > 4 are found in homogeneous transient nucleation, while heterogeneous nucleation gives exponent values <4 [31]. As transformation progresses, *n* becomes lower. The degree of transformation of 2PdG at 50 °C was only about 20%, hence, its high *n*. Values of about 1.5 can be ascribed to diffusion-controlled growth in pre-existing nuclei [32].

<table><caption>**Table 2.** Kinetic parameters of the metal to hydride conversion.</caption><thead><tr><td>Samples</td><td>Particle Size nm (TEM)</td><td>Ln k</td><td>n</td></tr></thead><tbody><tr><td>2PdG</td><td>4.9</td><td>-25.2</td><td>8.6</td></tr><tr><td>1PdGONE</td><td>8.9</td><td>-3.5</td><td>1.9</td></tr><tr><td>1PdGOE</td><td>9.7</td><td>-2.7</td><td>1.6</td></tr></tbody></table>

## 4. Conclusions

Depending on the Pd particle size, the reversible transformation of PdHx to Pd is more or less rapid. For instance, in the case of 1PdGOE, the transformation of the Pd hydride (having high β PdHx proportion) occurs with the highest velocity in comparison with the others. At the same time, kinetic studies prove that the hydride transformation of 2PdG proceeds at a slower rate. Relevantly, butadiene hydrogenation reaction behavior (activity and selectivity) as well as reversible Pd to PdHx transformation are intimately linked with the Pd nanoparticle size. Decreases in Pd particle size, increases in specific catalytic activity, and decreases in butene selectivity follow the same pattern: 1PdGOE > 1PdGONE > 2PdG. Larger Pd nanoparticles favor H solubility, reaching H-richer hydride species such as beta-hydride. Higher H content accelerated the transformation of PdHx into Pd during the reaction and, consecutively, the partial hydrogenation of butadiene to butenes. On the other hand, Pd and alfa-PdH lead to the butane production.

In this study, by simultaneously exploring the formation of subsurface Pd-hydride and the catalytic properties (selectivity and activity in the butadiene hydrogenation), it is demonstrated that the participation of hydride species formed inside Pd crystallites is controlling the dynamics of the hydrogen activation and insertion (reversible transformation of Pd to PdHx) in the reaction products.

**Supplementary Materials:** The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/nano13061099/s1, Figure S1. HEXRD of Pd (111) during temperature programmed reduction (TPR) (a), upon cooling (b), during reaction (temporal evolution when switching from H₂ to reaction mixture) at different temperatures 30 (c), 40 (d) and 60 °C (e) and conversion/selectivity (f) for sample 1PdG. Figure S2. HEXRD of Pd (111) during temperature programmed reduction (TPR) (a), upon cooling (b), temporal evolution on switching from H₂ to reaction
mixture during reaction at different temperatures 30 (c), and 60 °C (d) and conversion/selectivity results (e) for sample 1PdGOE. Figure S3. HEXRD of Pd (111) during temperature programmed reduction (TPR) (a), upon cooling (b), temporal evolution on switching from H₂ to reaction mixture during reaction at different temperatures 30 (c), 40 (d) and 60 °C (e) and conversion/selectivity (f) for sample 1PdGONE. Figure S4. HEXRD of Pd (111) for sample 1PdGOE during temperature programmed reduction (TPR) (a) and upon cooling (b). Figure S5. HEXRD of Pd (111) for sample 1PdGOE during reaction at different temperatures: 30 °C: (a–d); 40 °C: (e–h); 45 °C: (i–l); 50 °C: (m–p); and 60 °C: (q–t). Switch (1): Hydrogen to butadiene reaction mixture. Switch (2): butadiene reaction mixture to hydrogen. Figure S5. (a,c,e,g,l,k,m,o,q,s) correspond to the acquisition during steady state before switching to H₂ or to reaction mixture. Figure S6. Conversion/selectivity for sample 1PdGOE fast. Figure S7. HEXRD of Pd (111) for sample 1PdGONE during temperature programmed reduction (TPR) (a) and upon cooling (b). Figure S8. HEXRD of Pd (111) for sample 1PdGONE during reaction at different temperatures: 30 °C: (a–d); 40 °C: (e–h); and 50°C: (I–l). Switch (1): Hydrogen to butadiene reaction mixture. Switch (2): butadiene reaction mixture to Hydrogen. Figure S8. (a,c,e,g,l,k) correspond to the acquisition during steady state before switching to H₂ or to reaction mixture. Figure S9. Conversion/selectivity for sample 1PdGONE fast. Figure S10. Deconvolution of PdHx peak for 1PdGOE and for 1PdGONE samples under H₂ atmosphere (a,b) and under reaction mixture conditions (c,d). Figure S11. HEXRD of Pd (111) for sample 2PdG during temperature programmed reduction (TPR) (a) and upon cooling (b). Figure S12. HEXRD of Pd (111) for sample 2PdG during reaction at different temperatures: 30 °C: (a–d); 50 °C: (e–h); and 70 °C: (i–l). Switch (1): Hydrogen to butadiene reaction mixture. Switch (2): butadiene reaction mixture to Hydrogen. Figure S12. (a,c,e,g,l,k) correspond to the acquisition during steady state before switching to H₂ or to reaction mixture. Figure S13. Conversion/selectivity for sample 2PdG fast. Figure S14. TEM images (a–c) and particle size histogram (d) of 1PdGOE. Figure S15. TEM images (a–c) and particle size histogram (d) of 1PdGONE. Figure S16. TEM images (a–c) and particle size histogram (d) of 2PdG. Figure S17. TEM images (a–c) and particle size histogram (d) of 1PdG. Figure S18. Kinetics of transformation of PdHx to Pd: 2PdG (a,b), 1PdGONE (c,d), 1PdGOE (e,f).

**Author Contributions:** Conceptualization, A.G.-R., I.R.-R. and M.F.-G.; experimental methodology, E.A.-N., M.D.M., A.I.-J., A.G.-R. and M.F.-G.; data analysis: E.A.-N., M.D.M., A.I.-J. and M.F.-G.; investigation, all authors; data curation, A.I.-J., M.D.M. and M.F.-G.; writing—original draft preparation, E.A.-N. and A.I.-J.; writing—review and editing, E.A.-N., I.R.-R. and A.G.-R.; supervision, A.G.-R., M.F.-G. and I.R.-R. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research received no external funding apart from that mentioned in the acknowledgments.

**Data Availability Statement:** Not applicable.

**Acknowledgments:** The authors acknowledge financial support from Spanish Ministry MINECO (projects PID2020-119160RB-C21 and -C22). They are also grateful to the ESRF for granting beamtime at ID15.

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Boudart, M. *Advances in Catalysis*; Herman, P., Weisz, P.B., Eley, D.D., Eds.; Catalysis by Supported Metals; Academic Press: Cambridge, MA, USA, 1969; Volume 20, pp. 153–166.

2. Kyriakou, G.; Boucher, M.B.; Jewell, A.D.; Lewis, E.A.; Lawton, T.J.; Baber, A.E.; Tierney, H.L.; Flytzani-Stephanopoulos, M.; Sykes, E.C. Isolated metal atom geometries as a strategy for selective heterogeneous hydrogenations. *Science* **2012**, *335*, 1209–1212. [CrossRef]

3. Uzio, D.; Berhault, G. Factors Governing the Catalytic Reactivity of Metallic Nanoparticles. *Catal. Rev. 2010*, *52*, 106–131. [CrossRef]

4. Bielawa, H.; Hinrichsen, O.; Birkner, A.; Muhler, M. The Ammonia-Synthesis Catalyst of the Next Generation: Barium-Promoted Oxide-Supported Ruthenium. *Angew. Chem. Int. Ed.* **2001**, *6*, 1061–1063. [CrossRef]

5. Somorjai, G.A. The structure sensitivity and insensitivity of catalytic reactions in light of the adsorbate induced dynamic restructuring of surfaces. *Catal. Lett.* **1990**, *7*, 169. [CrossRef]

6. Li, Y.; Yan, K.; Cao, Y.; Ge, X.; Zhou, X.; Yuan, W.; Chen, D.; Duan, X. Mechanistic and Atomic-Level Insights into Semihydrogenation Catalysis to Light Olefins. *ACS Catal.* **2022**, *12*, 12138–12161. [CrossRef]

7. Bernsmeier, D.; Chuenchom, L.; Paul, B.; Rümmler, S.; Smarsly, B.; Kraehnert, R. Highly Active Binder-Free Catalytic Coatings for Heterogeneous Catalysis and Electrocatalysis: Pd on Mesoporous Carbon and Its Application in Butadiene Hydrogenation and Hydrogen Evolution. *ACS Catal.* **2016**, *6*, 8255–8263. [CrossRef]
8. Chesnokov, V.V.; Podyacheva, O.Y.; Richards, R.M. Influence of carbon nanomaterials on the properties of Pd/C catalysts in selective hydrogenation of acetylene. *Mater. Res. Bull.* **2017**, *88*, 78–84. [CrossRef]

9. Silvestre-Albero, J.; Rupprecter, G.; Freund, H.-J. From Pd nanoparticles to single crystals: 1,3-butadiene hydrogenation on well-defined model catalysts. *Chem. Comm.* **2006**, *80–82*. [CrossRef] [PubMed]

10. Guerrero, A.; Reading, M.; Grillet, Y.; Rouquerol, J.; Boitiaux, J.P.; Cosyns, J. Influence of dispersion on the energies of adsorption: H2, CO, propylene and propyne on supported Pd or Pt. *Z. Phys. D At. Mol. Clusters* **1989**, *12*, 583–586. [CrossRef]

11. Teschner, D.; Revay, Z.; Borsodi, J.; Havecker, M.; Knop-Gericke, A.; Schlögl, R.; Milroy, D.; Jackson, S.D.; Torres, D.; Sautet, P. Understanding Palladium Hydrogenation Catalysts: When the Nature of the Reactive Molecule Controls the Nature of the Catalyst Active Phase. *Angew. Chem. Int. Ed.* **2008**, *47*, 9274–9278. [CrossRef]

12. Katano, S.; Kato, H.S.; Kawai, M.; Domen, K. Self-Activated Catalyst Layer for Partial Hydrogenation of 1,3-Butadiene on a Hydrogen-Precovered Pd(110) Surface. *J. Phys. Chem. C* **2009**, *113*, 14872–14878. [CrossRef]

13. Wu, Y.; Zhao, W.; Wang, Y.; Wang, B.; Fan, M.; Zhang, R. Enhancing Catalytic performance through Subsurface Chemistry: The Case of C₂H₂ Semihydrogenation over Pd Catalysts. *ACS Appl. Mater. Interfaces* **2022**, *14*, 56743–56757. [CrossRef]

14. Morales, M.V.; Asedegbega-Nieto, E.; Bachiller-Baeza, B.; Guerrero-Ruiz, A. Bioethanol dehydrogenation over copper supported on functionalized graphene materials and a high surface area graphite. *Carbon* **2016**, *102*, 426–436. [CrossRef]

15. Castillejos, E.; Bachiller-Baeza, B.; Asedegbega-Nieto EGuerrero-Ruiz, A.; Rodriguez-Ramos, I. Selective 1,3-Butadiene hydrogenation by gold nanoparticles deposited&precipitated onto nano-carbon materials. *RSC Adv.* **2015**, *5*, 81583–81598.

16. Kumara, L.S.R.; Sakata, O.; Kobayashi, H.; Song, C.; Kohara, S.; Ina, T.; Yoshimoto, T.; Yoshioka, S.; Matsumura, S.; Kitagawa, H. Hydrogen storage and stability properties of Pd–Pt solid-solution nanoparticles revealed via atomic and electronic structure. *Sci. Rep.* **2017**, *7*, 14606. [CrossRef]

17. Gross, O.; Neuber, N.; Kuball, A.; Bochtler, B.; Hechler, S.; Frey, M.; Busch, R. Signatures of structural differences in Pt–P- and Pd–P-based bulk glass-forming liquids. *Commun. Phys.* **2019**, *2*, 83. [CrossRef]

18. Langhammer, C.; Zhdanov, V.P.; Zoric, I.; Kasemo, B. Size-dependent hysteresis in the formation and decomposition of hydride in metal nanoparticles. *Chem. Phys. Lett.* **2010**, *488*, 62–66. [CrossRef]

19. Suzana, A.F.; Wu, L.; Assefa, T.A.; Williams, B.P.; Harder, R.; Cha, W.; Kuo, C.-H.; Tsung, C.-K.; Robinson, I.K. Structure of a seeded palladium nanoparticle and its dynamics during the hydride phase transformation. *Commun. Chem.* **2021**, *4*, 64. [CrossRef]

20. Bauer, M.; Schoch, R.; Shao, L.; Zhang, B.; Knop-Gericke, A.; Willinger, M.; Schlögl, R.; Teschner, D. Structure–Activity Studies on Highly Active Palladium Hydrogenation Catalysts by X-ray Absorption Spectroscopy. *J. Phys. Chem. C* **2012**, *116*, 22375–22385. [CrossRef]

21. Bugaev, A.L.; Guda, A.A.; Lazzarini, A.; Lomachenko, K.A.; Groppo, E.; Pellegrini, R.; Piovano, A.; Emerich, H.; Soldatov, A.V.; Bugaev, L.A.; et al. In situ formation of hydrides and carbides in palladium catalyst: When XANES is better than EXAFS and XRD. *Catal. Today* **2017**, *283*, 119–126. [CrossRef]

22. Bugaev, A.L.; Usoltsev, O.A.; Lazzarini, A.; Lomachenko, K.A.; Guda, A.A.; Pellegrini, R.; Carosso, M.; Vitillo, J.G.; Groppo, E.; van Bokhoven, J.A.; et al. Time-resolved operando studies of carbon supported Pd nanoparticles under hydrogenation reactions by X-ray diffraction and absorption. *Faraday Discuss.* **2018**, *208*, 187–205. [CrossRef] [PubMed]

23. Liu, Y.; Fu, F.; McCue, A.; Jones, W.; Rao, D.; Feng, J.; He, Y.; Li, D. Adsorbate-Induced Structural Evolution of Pd Catalyst for Selective Hydrogenation of Acetylene. *ACS Catal.* **2020**, *10*, 15048–15059. [CrossRef]

24. Vogel, W.; He, W.; Huang, Q.-H.; Zou, Z.; Zhang, X.-G. Palladium nanoparticles “breathe” hydrogen; a surgical view with X-ray diffraction. *Int. J. Hydrogen Energy* **2010**, *35*, 8609–8620. [CrossRef]

25. Khanuja, M.; Mehta, B.R.; Agar, P.; Kulriya, P.K.; Avasthi, D.K. Hydrogen induced lattice expansion and crystallinity degradation in palladium nanoparticles: Effect of hydrogen concentration, pressure, and temperature. *J. Appl. Phys.* **2009**, *106*, 093515. [CrossRef]

26. Tew, M.W.; Miller, J.T.; van Bokhoven, J.A. Particle Size Effect of Hydride Formation and Surface Hydrogen Adsorption of Nanosized Palladium Catalysts: L₃ Edge vs K Edge X-ray Absorption Spectroscopy. *J. Phys. Chem. C* **2009**, *113*, 15140–15147. [CrossRef]

27. Zhang, P.; Sham, T.K. X-Ray Studies of the Structure and Electronic Behavior of Alkanethiolate-Capped Gold Nanoparticles: The Interplay of Size and Surface Effects. *Phys. Rev. Lett.* **2003**, *90*, 245502. [CrossRef] [PubMed]

28. Borodziński, A.; Bond, G.C. Selective Hydrogenation of Ethyne in Ethene-Rich Streams on Palladium Catalysts, Part 2: Steady-State Kinetics and Effects of Palladium Particle Size, Carbon Monoxide, and Promoters. *Catal. Rev. 2008*, *50*, 379–469. [CrossRef]

29. Woldt, E. The relationship between isothermal and non-isothermal description of Johnson-Mehl-Avrami-Kolmogorov kinetics. *J. Phys. Chem. Solids* **1992**, *53*, 521–527. [CrossRef]

30. Li, H.; Gai, K.; He, L.; Zhang, C.; Cui, H.; Li, M. Non-isothermal phase-transformation kinetics model for evaluating the austenization of 55CrMo steel based on Johnson-Mehl-Avrami equation. *Mater. Des.* **2016**, *92*, 731–741. [CrossRef]

31. Sinha, I.; Mandal, R.K. Avrami exponent under transient and heterogeneous nucleation transformation conditions. *J. Non-Cryst. Solids* **2011**, *357*, 919–925. [CrossRef]

32. Bruna, P.; Crespo, D.; Gonzalez-Cinca, R.; Pineda, E. On the validity of Avrami formalism in primary crystallization. *J. Appl. Phys.* **2006**, *100*, 054907. [CrossRef]

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.