Hydrogenation of 2-methyl-3-butyn-2-ol over a Pd/ZnO catalyst: kinetic model and selectivity study

S. Vernuccio *a, R. Goy b, Ph. Rudolf von Rohr a, J. Medlock b and W. Bonrath b
aInstitute of Process Engineering, ETH Zurich, Sonneggstrasse 3, 8092 Zurich, Switzerland. E-mail: svernuccio@ipe.mavt.ethz.ch; Fax: +44 446321325; Tel: +44 446322499
bDSM Nutritional Products, Research and Development, Basel, Switzerland

Received 4th May 2016 , Accepted 30th June 2016

First published on 7th July 2016


Abstract

The three-phase hydrogenation of 2-methyl-3-butyn-2-ol has been studied over a Pd/ZnO catalyst. A Langmuir–Hinshelwood mechanism was applied assuming noncompetitive adsorption between hydrogen and organic molecules on the catalyst active sites. All experimental runs used for the modeling have been obtained in the intrinsic kinetic regime in order to exclude any mass transfer limitation. An optimization procedure allowed the estimation of the kinetic and adsorption parameters governing the process. The results revealed that the proposed model accurately describes the behavior of the system in the typical operating ranges of industrial reactors. The performance of the catalyst in terms of selectivity to 2-methyl-3-buten-2-ol and initial activity is found to be higher compared with that of a commercial Lindlar catalyst under the same operating conditions. The mathematical model, successfully validated, is able to accurately predict the selectivity of the process.


Introduction

Catalytic hydrogenations of acetylenic alcohols to alkenols are carried out in industry to produce fine chemicals and intermediate products.1 One example is the hydrogenation of 2-methyl-3-butyn-2-ol (MBY) to 2-methyl-3-buten-2-ol (MBE), a key product in the manufacture of vitamins and aroma compounds.2,3

The reaction pathway of MBY hydrogenation is presented in Scheme 1. The selectivity of this process to MBE is affected by the over-hydrogenation reaction to 2-methyl-2-butanol (MBA) and by the formation of dimeric C10 products.


image file: c6re00093b-s1.tif
Scheme 1 Reaction network of MBY hydrogenation.

The industrial hydrogenation of MBY is conventionally conducted in slurry reactors with Lindlar catalyst, a 5 wt% Pd/CaCO3 powdered catalyst modified by lead.4

High selectivities to MBE are required due to the difficulties in separation of the mixture MBE/MBA.2 A typical approach to enhance the selectivity of the reaction to MBE is the addition of nitrogen-based modifiers (e.g. quinoline, pyridine, ammonia).5 These electron donor compounds inhibit alkene surface interactions due to their preferential adsorption on the surface of the catalyst.6 However, this incorporation requires further process steps, such as separation, resulting in additional costs.

Alternative approaches have mostly been directed on catalyst modifications. Bimetallic catalysts, for example, are known to improve olefin selectivity during alkyne hydrogenations.7,8

Zinc, in particular, is known to act as a promoter in the palladium-based catalyst because of its ability to form an intermetallic Pd–Zn phase affecting the adsorption strength of alkynes and alkenes.5,9

The performance of Pd-based catalysts during alkynes hydrogenations is also strongly influenced by the nature of the support. ZnO in particular is often considered a good alternative to CaCO3 with respect to reaction selectivity.10 High selectivities towards the alkene have been reported for monometallic Pd-based catalysts over ZnO.5,9,11 The electron donating effect of Zn results in an higher electron density of Pd. The selectivity of the process is enhanced due to a decreased alkene adsorption.12 This makes these catalysts a valuable alternative to the lead-containing Lindlar catalyst, well suited for industrial applications.

Semagina et al. studied the reaction kinetics of the water-assisted MBY hydrogenation over a Pd/ZnO catalyst coated on sintered metal fibers.5 In that work, the reaction kinetics were modeled at 308 K and 0.5 MPa, assuming a typical Langmuir–Hinshelwood kinetic mechanism with competitive adsorption of hydrogen and organic species on the catalyst surface. A similar mechanism was presented by Crespo-Quesada et al. to describe the solvent-free hydrogenation of MBY over a Pd/ZnO structured catalyst at 348 K and 0.8 MPa.9

Protasova et al. tested ZnO nanowires coated with Pd nanoparticles in the hydrogenation of MBY at 323 K and 0.5 MPa. High MBE selectivities of 95% were observed at 99% of conversion and the results were modeled assuming a typical Eley–Rideal kinetic mechanism.11

No attempt were made in these previous studies to describe the behavior of the system under a wide range of operating conditions.

In a previous experimental study Vernuccio et al. presented a kinetic investigation of the selective hydrogenation of MBY over a commercial Lindlar catalyst based on the noncompetitive adsorption between hydrogen and organics species on the catalyst surface.13 In the aforementioned study a general kinetic model is proposed to predict the evolution of the system under varying experimental conditions. The present work aims to extend this model to simulate MBY hydrogenation reaction over a Pd/ZnO powdered catalyst under varying temperatures (313–353 K), pressures (0.3–1.0 MPa) and catalyst loadings (0.25–0.75 wt%). All the reactions involved in the network are completely characterized from a kinetic point of view with the estimation of all the kinetic and adsorption parameters (pre-exponential factors, activation energies and enthalpies of adsorption). The performance of this catalyst in terms of selectivity are compared with the ones achievable with a commercial Lindlar catalyst under the same operating conditions and the same palladium loading.

Experimental

Materials

2-Methyl-3-butyn-2-ol (purum, >98%) was supplied by DSM Nutritional Products. 2-Methyl-3-buten-2-ol (purum, ≥98%) and 2-methyl-2-butanol (purum, ≥99%) were purchased from Sigma-Aldrich. The physical properties of MBY, measured by DSM Nutritional Products at Kaiseraugust/Switzerland are listed in the ESI. A Pb-modified Lindlar catalyst (5 wt% Pd/CaCO3), purchased from Sigma-Aldrich was used for comparison. Pure hydrogen (99.995%) and nitrogen (99.995%), for inertization purposes, were obtained from Pangas. All the reagents were used as received.

Pd/ZnO catalyst preparation and characterization

EOS MaragingSteel MS1 metal powder (20.0 g) was heated in an oven under air at 723 K for 3 h. The resulting metal powder was initially coated with an Al2O3/ZnO base layer. Separately a solution of Pd nanoparticles was prepared and deposited on the metal powder containing the Al2O3/ZnO layer.

A solution of Al(NO3)3·9H2O (200.0 g) in water (700 mL) was heated to 368 K and ZnO powder (43.4 g) was added slowly. The mixture was stirred until completely dissolved. After cooling to room temperature, the solution was filtered. The oxidized metal powder (20.0 g) was added to 50 mL of the Al/Zn solution and the mixture was stirred at room temperature for 15 min, filtered through a 0.45 μm membrane filter and dried under vacuum at 313 K for a minimum of 2 h. This process was repeated two additional times.

Sodium molybdate dihydrate (162 mg) and anhydrous palladium(II)chloride (108 mg) were added in 30 mL of deionized water. The solution was evaporated slowly under stirring until a solid residue was formed. Afterwards, 30 mL of deionized water were added to the solid residue under stirring. This evaporation-dissolving cycle was repeated two times in order to completely dissolve palladium salt. Finally, 50 mL of hot water was added to the solid residue. The deep brown solution was cooled down to room temperature, filtered and water added to adjust the volume to 60 mL. Hydrogen gas was slowly bubbled through the solution at room temperature for 60 min.

The obtained solution of Pd(0) was added to a suspension of the metal powder – Al2O3/ZnO in water (60 mL). The mixture was stirred for 30 min and then filtered. The catalyst was dried under vacuum at 313 K overnight. The catalyst was thermally activated in a tube furnace at 573 K for 4 h (temperature ramp 10 K min−1) under H2 – Ar flow (1[thin space (1/6-em)]:[thin space (1/6-em)]9; total flow rate: 450 mL min−1).

The microstructure of the catalyst was observed by FIB/SEM (Zeiss NVision 40) equipped with EDX apparatus. Samples were opened at the region of interest through milling with gallium-ions at 30 kV/6 nA perpendicular to the sample surface. The cross section was polished at 30 kV/1.5 nA and imaged by detection of secondary electrons (SE) using an in-lens detector. EDX measurements were performed at 9 kV.

The morphology of the supported Pd-nanoparticles was studied by HAADF-STEM conducted on an aberration-corrected dedicated STEM microscope HD-2700CS (Hitachi). The microscope was operated at an acceleration potential of 200 kV.

Temperature programmed oxidation/reduction (TPO/TPR) were performed in a Micrometrics AutoChem 2920 II instrument using 5 vol% of O2/H2 in Ar. The data were collected at a ramp of 8 K min−1.

Experimental setup

In the typical experiment for the kinetic modeling, 200.0 g of MBY were charged with the desired quantity of catalyst in a 400 mL semi-batch reactor (Premex Reactor AG, Lengnau, Switzerland). Isothermal conditions were kept by an electrical heating jacket and a cooling system using cold water. The reactor was equipped with a system of four equidistant baffles (width = 12 mm, thickness = 2.5 mm) and a four-blades gas-entrainment stirrer (width = 12 mm, diameter = 38 mm).

Pure hydrogen was supplied only after stabilization of the temperature at the required value under nitrogen atmosphere. The pressure in the reactor was maintained constant during the course of the experiments by supplying hydrogen from an external reservoir. The reaction mixture was stirred at 20 Hz.

Liquid samples were periodically withdrawn from the reactor and analyzed by gas-chromatography (Bruker GC-450 with flame ionization detector). The analysis was conducted by means of a VF-Wax ms column (25 m × 0.25 mm, coating thickness = 0.25 μm) operated at the following conditions: FID temperature, 493 K; injector temperature, 523 K; oven temperature, 353 K to 473 K with a temperature program of 20 K min−1.

Results and discussion

Catalyst characterization

Some of the physical characteristics of the Pd/ZnO catalyst are listed in Table 1.
Table 1 Catalyst characteristics
Property Value
Cumulative particle size distribution μm <73 for 100%
<43 for 73%
<36 for 58%
<21 for 20%
<10 for 3%
<4 for 1.5%
<2 for 1%
d 50 μm 32
BET surface area m2 g−1 1.5
Pd loading wt% 0.22
Density, ρS kg m−3 ∼8000


Analysis of the catalyst surface using BET techniques gives a surface of ∼1.5 m2 g−1, in addition to a total specific pore volume of ∼2 × 10−3 cm3 g−1 and a mean pore diameter of 6.0 nm. The extremely low pore volume is beneficial in order to avoid internal diffusion limitations. TPO and TPR heating ramps of the catalyst are depicted in the ESI.

SEM microphotographs of a catalyst particle cross section with two different magnifications are showed in Fig. 1.


image file: c6re00093b-f1.tif
Fig. 1 SEM microphotographs of a catalyst particle with diameter of about 50 μm. Detector: in-lens, voltage: 9 kV. (a) Cross section of the catalyst particle. (b) Magnified image of the selected area (red box) in (a).

The metal particle appears to be fully coated with an homogenous base layer. The coating layer (about 700 nm thick) remains intact after Pd deposition and thermal activation. EDX analysis reveals that Pd is dispersed only on the external surface of the catalyst while Al and Zn are present in the base layer (see ESI).

The surface of the catalyst appears to be decorated by large metal clusters (sizes 100–200 nm) containing Pd according to EDX measurement.

HAADF-STEM images provided in Fig. 2 reveal that the supported Pd-nanoparticles exhibit different diameters up to ca. 10 nm. Nevertheless they are often agglomerated in certain regions. EDX analysis of two different areas of Fig. 2(a) are provided in the ESI.


image file: c6re00093b-f2.tif
Fig. 2 a, b HAADF-STEM images of the catalyst surface with Pd nanoparticles dispersed on Al2O3/ZnO.

Mass transport resistances

For the purpose of kinetic study of the hydrogenation of MBY, it is necessary to characterize quantitatively the gas–liquid (G–L) and liquid–solid (L–S) mass transport resistances. The following criteria, proposed by Ramachandran and Chaudari14 were used:
 
image file: c6re00093b-t1.tif(1)
 
image file: c6re00093b-t2.tif(2)
where r0 is the observed initial reaction rate of MBY hydrogenation.

The gas absorption method15 allowed the experimental estimation of the G–L mass transfer coefficient in the absence of the catalyst.

The theoretical calculation of the L–S mass transfer coefficient kSaS is extensively described in a previous work.11 The catalyst particles were approximated as spherical, according to the particle size distribution by weight basis (see Table 1).

The parameters α1 and α2 calculated for the 14 experimental runs are listed in the ESI. The reported values indicate the absence of any external mass-transfer resistance.

The absence of microporosity, confirmed by BET analysis, justifies the lack of internal liquid–solid diffusion limitations.

The initial reaction rate was found to vary linearly with the catalyst loading in the investigated range of temperature as depicted in Fig. 3. This trend was considered an additional verification that the reaction is under kinetic control.


image file: c6re00093b-f3.tif
Fig. 3 Effect of catalyst loading on the initial rate of MBY hydrogenation. Conditions: 0.7 MPa.

The solubility of hydrogen in MBY was calculated from the following relationship:13

 
image file: c6re00093b-t3.tif(3)

Kinetic model

During a previous study Vernuccio et al.13 developed a general kinetic model to describe the hydrogenation of MBY over a commercial Pd-based catalyst. In that work a Langmuir–Hinshelwood model with noncompetitive adsorption of hydrogen and organics on the catalyst surface and dissociative adsorption of hydrogen was assumed. Similar hypothesis have been previously validated to develop the kinetic expressions of the reaction rates of the hydrogenation of 1-butyne16 and n-butenes.17

From that revision, the kinetic mechanism depicted in Table 2 was proposed to describe the hydrogenation of MBY. Similar reaction steps were assumed for the over hydrogenation to MBA (r2) and the formation of dimers (r3). The hypothesis of noncompetitive adsorption results in the presence of different active sites for hydrogen (○) and organics (●) adsorption.

Table 2 Kinetic mechanism for MBY hydrogenation13
image file: c6re00093b-t4.tif (e1) image file: c6re00093b-t5.tif (e4)
image file: c6re00093b-t6.tif (e2) image file: c6re00093b-t7.tif (e5)
image file: c6re00093b-t8.tif (e3) image file: c6re00093b-t9.tif (e6)


All adsorption/desorption steps were assumed to be reversible while the last hydrogen addition was considered the rate determining step.

The application of the aforementioned mechanism leads to the following simplified kinetic expressions under the assumption of weak adsorption of hydrogen and alkane:

 
image file: c6re00093b-t10.tif(4)
 
image file: c6re00093b-t11.tif(5)
 
image file: c6re00093b-t12.tif(6)
where ki are the apparent kinetic constants for the 3 hydrogenation reactions.

The hydrogen concentration in the liquid bulk CH was considered equal to the hydrogen solubility calculated with eqn (3).

The variation of the concentrations of the species involved in the system can be expressed in the kinetic regime by the following set of ordinary differential equations:

 
image file: c6re00093b-t13.tif(7)
 
image file: c6re00093b-t14.tif(8)
 
image file: c6re00093b-t15.tif(9)
 
image file: c6re00093b-t16.tif(10)

The temperature dependence of the kinetic constants ki and of the equilibrium adsorption constants Ki was assumed to obey an Arrhenius-type equation, according to the relations:

 
image file: c6re00093b-t17.tif(11)
 
image file: c6re00093b-t18.tif(12)

Conversion of MBY and selectivity of the reaction to MBE were defined as:

 
image file: c6re00093b-t19.tif(13)
 
image file: c6re00093b-t20.tif(14)

During this study, 14 runs on MBY hydrogenation have been carried out in a variety of experimental conditions listed in Table 3. The first 12 of these runs have been used to estimate the parameters governing the process (pre-exponential factors, activation energies and enthalpies of adsorption) while the 2 remaining have been employed for a validation of the model.

Table 3 Overview of the operating conditions adopted during the hydrogenation experiments on pure MBY with the overall percentage standard deviations obtained with the model
Run T K p MPa Catalyst loading wt% Overall σ%
1 333 0.7 0.50 2.5
2 353 0.7 0.50 6.2
3 333 0.7 0.75 4.7
4 353 0.3 0.75 3.9
5 333 1.0 0.50 6.5
6 333 0.4 0.50 5.7
7 313 0.7 0.75 5.8
8 353 0.7 0.25 4.8
9 313 0.7 0.50 6.8
10 353 0.4 0.50 4.5
11 353 0.9 0.50 3.1
12 333 0.7 0.25 2.1
13 343 0.7 0.50 5.6
14 333 0.5 0.75 5.0


For each experiment, approximately 10 concentration data points at different reaction times were recorded. The system of differential and algebraic equations representing the proposed mathematical model (eqn (3)–(12)) has been numerically solved according to proper initial conditions. A numerical optimization procedure was implemented to minimize an objective function expressed as:

 
image file: c6re00093b-t21.tif(15)
where yi,j,k and ci,j,k are the experimental and the calculated concentrations, m is the number of experimental points recorded during each run, n the number of detected species and q the number of the experiments involved in the optimization. The analysis was conducted using 456 data points to estimate 14 parameters.

Some examples of calculated and experimental concentrations for single experiments are shown in Fig. 4 and 5 (overall percentage standard deviations 4.7% and 6.2% respectively).


image file: c6re00093b-f4.tif
Fig. 4 Experimental points and calculated kinetic curves. Conditions: 333 K, 0.7 MPa, 0.75 wt% of catalyst.

image file: c6re00093b-f5.tif
Fig. 5 Experimental points and calculated kinetic curves. Conditions: 353 K, 0.7 MPa, 0.50 wt% of catalyst.

The total percentage standard deviation for all the performed experiments are listed in Table 3. In all cases these overall deviations are comparable with those related to the single components.

Fig. 6 and 7 show the comparison of simulated results with the experimental concentrations of MBY and MBE for different temperatures and pressures. The model appears to be consistent with the experimental results in the investigated range of operating conditions.


image file: c6re00093b-f6.tif
Fig. 6 Effect of temperature on MBY hydrogenation: MBY (a) and MBE (b) experimental and calculated concentrations. Conditions: 0.7 MPa, 0.50 wt% of catalyst.

image file: c6re00093b-f7.tif
Fig. 7 Effect of pressure on MBY hydrogenation: MBY (a) and MBE (b) experimental and calculated concentrations. Conditions: 333 K, 0.50 wt% of catalyst.

Table 4 presents the best-estimated values of the kinetic and adsorption parameters with their 95% confidence intervals, resulting from the optimization procedure.

Table 4 Kinetic parameters estimates with 95% confidence intervals
k 01 [mol molPd−1 s−1] (1.04 ± 0.05) × 108 E 1 [kJ mol−1] 26.3 ± 1.3
k 02 [mol molPd−1 s−1] (8.52 ± 0.5) × 109 E 2 [kJ mol−1] 37.5 ± 1.9
k 03 [mol molPd−1 s−1] (1.45 ± 0.2) × 108 E 3 [kJ mol−1] 28.9 ± 1.4
K 0Y [m3 mol−1] (3.76 ± 0.2) × 10−8 −ΔHadY [kJ mol−1] 26.7 ± 1.3
K 0E [m3 mol−1] (6.52 ± 0.3) × 10−9 −ΔHadE [kJ mol−1] 18.4 ± 0.9
K 0H [m3 mol−1] (1.78 ± 0.1) × 10−5 −ΔHadH [kJ mol−1] 4.9 ± 0.2
K [m3 mol−1] (3.67 ± 0.2) × 10−5 −ΔHadφ [kJ mol−1] 6.5 ± 0.3


The values of the pre-exponential factors k0i and K0i previously estimated for the Lindlar catalyst13 were used as initial data for the optimization procedure. The best-estimated pre-exponential factors do not differ significantly from the initial data with the only exception of k03. This discrepancy is related to the much lower concentration of dimers detected in the reaction system with respect to the hydrogenation conducted over the Lindlar catalyst. However, the estimation of this constant is less reliable compared to the others since the associated error is greater than 10%.

The activation energy of MBY hydrogenation E1 is found to be in good agreement with the value of 25.2 ± 1.6 kJ mol−1 reported by Crespo-Quesada et al. for a structured Pd/ZnO catalyst.9

The ratio between the kinetic coefficients of alkene and alkyne hydrogenation k2/k1 results to be higher than 1 in the investigated range of temperatures. According to the thermodynamic of the alkynes the over-hydrogenation reaction is favored from a kinetic point of view.18 The reason of the high selectivity of the process to MBE is attributed to the stronger adsorption of the alkyne compared to the alkene. However, the adsorption equilibrium constant of MBY is found to be 5.83 × 10−4 m3 mol−1 at 333 K, slightly lower than the ones previously published for Lindlar catalyst, resulting in a lower adsorption strength of the alkyne. On the other hand the ratio of the adsorption equilibrium constants of MBY and MBE KY/KE ranges between 97 and 140 depending on the investigated temperature. This value is 40% higher than the corresponding ones previously calculated for the Lindlar catalyst.13,19

In accordance with published data, the formation of a Pd–Zn alloy can change the adsorption strength of alkynes and alkenes affecting the selectivity of the hydrogenation.7,20 The addition of Zn, in particular, is considered to be beneficial for the selectivity of the process towards the alkene.7 It affects the electronic structure of Pd resulting in a decrease of the heat of adsorption of alkynes and alkenes.8,21 As a result of the observed increase of the ratio KY/KE, the selectivity of the process to MBE is significantly enhanced. This tendency is depicted in Fig. 8 and 9 where the selectivity of the process to MBE over the Pd/ZnO is compared with the one obtained with the Lindlar catalyst under the same experimental conditions and the same palladium loading. In both cases the Pd/ZnO exhibits an higher selectivity compared with the Lindlar catalyst at various conversions.


image file: c6re00093b-f8.tif
Fig. 8 Selectivity to MBE vs. MBY conversion during the hydrogenation of MBY over Pd/ZnO (experimental and calculated) and Lindlar catalyst (experimental). Conditions: 333 K, 0.7 MPa, 1.1 × 10−3 wt% of Pd.

image file: c6re00093b-f9.tif
Fig. 9 Selectivity to MBE vs. MBY conversion during the hydrogenation of MBY over Pd/ZnO (experimental and calculated) and Lindlar catalyst (experimental). Conditions: 353 K, 0.9 MPa, 1.1 × 10−3 wt% of Pd.

The continuous lines represent the theoretical selectivity calculated using the best-estimated values of the parameters. It is worthy to note that the model delivers good predictions of the system behavior even at high conversions of MBY. Therefore, it allows estimation of the experimental selectivity of the process with good accuracy.

The reaction carried out at 333 K and 0.7 MPa using the commercial Lindlar catalyst showed a lower initial activity of 0.9 mol molPd−1 s−1 compared with the 3.0 mol molPd−1 s−1 of the Pd/ZnO catalyst. Similar trends were obtained during previous studies on the hydrogenation of MBY over Lindlar and Pd/ZnO catalysts.5,9 The higher activity may be ascribed to the presence of the Pd–Zn alloy in the active phase of the catalyst suppressing the overhydrogenation of MBE to MBA.

The last two experiments listed in Table 3 (not employed during the parametric optimization) were used to validate the proposed model. During this phase, the best-estimated values of the parameters have only been used to predict the concentration profiles of the species involved in the system without any further adjustment. A comparison between experimental and simulated data, achieved during this validation procedure, is shown in Fig. 10.


image file: c6re00093b-f10.tif
Fig. 10 Comparison between experimental and simulated kinetic curves. Conditions: 343 K, 0.7 MPa, 0.50 wt% of catalyst.

The good agreement between experimental data and kinetic curves, supported by the overall standard deviations in line with those of the other experiments, confirms the reliability of the model.

Catalyst reuse

Stability of the catalyst was evaluated performing hydrogenation experiments over 6 consecutive reaction cycles. Between each run the catalyst was filtered, washed with isopropanol and dried. MBY initial reaction rate and selectivity to MBE for 6 consecutive runs conducted at 353 K and 0.7 MPa are presented in Fig. 11.
image file: c6re00093b-f11.tif
Fig. 11 Selectivity to MBE at 97% of conversion and MBY initial reaction rate for reaction over Pd/ZnO catalyst over consecutive runs. Conditions: 353 K, 0.7 MPa, 0.50 wt% of catalyst.

Both initial reaction rate and selectivity result constant over about 35 hours of reaction. The high stability of the catalyst can be attributed to the strong interaction between Pd nanoparticles and ZnO leading to the formation of Pd–Zn alloy.9,22

Conclusions

A Pd/ZnO catalyst has been developed and tested in the hydrogenation of MBY under solvent free conditions. External mass transfer limitations have been studied in order to ensure the kinetic regime of the experimental runs. The intrinsic kinetics of the reaction were assumed to obey a Langmuir–Hinshelwood mechanism with noncompetitive adsorption between hydrogen and organics. The kinetic and the adsorption coefficients governing the process were estimated with relative confidence intervals.

The proposed kinetic model reliably estimates the concentration profiles of the species involved in the system in the investigated ranges of operating conditions (313–353 K, 0.3–1.0 MPa) and for varying catalyst loadings (0.25–0.75 wt%). The goodness of the optimization is confirmed by the overall standard deviation lower than 7.0% obtained for each experimental run. The model also allows the selectivity of the process to MBE to be estimated with high accuracy.

The validation of the model was conducted using the best-estimated values of the kinetic parameters, without additional adjustments, to simulate the concentration profiles of the species involved in the system.

The high selectivity of the process to MBE is found to be related to the thermodynamics of the alkyne hydrogenation. Zn acts as a promoter in a Pd-based catalyst resulting in a lower adsorption strength of alkynes and alkenes. As result, the Pd/ZnO catalyst exhibits higher selectivity and higher initial activity compared to the lead-containing Pd catalysts. For this reason, it can be considered as a valuable alternative for the industrial hydrogenation of MBY with respect to the classic Lindlar catalyst.

Nomenclature

Roman symbols

a Specific gas–liquid contact area, m−1
a S Specific liquid–solid contact area, m−1
image file: c6re00093b-t22.tif Hydrogen solubility in the liquid at equilibrium, mol m−3
C i Concentration of species i, mol m−3
E i Activation energy of reaction i, kJ mol−1
F Objective function, (mol m−3)2
k i Kinetic constant of reaction i, mol (molPd s)−1
K i Equilibrium adsorption constant, m3 mol−1
k L Gas–liquid mass transfer coefficient, m s−1
k L a Volumetric gas–liquid mass transfer coefficient, s−1
k S Liquid–solid mass transfer coefficient, m s−1
k S a S Volumetric liquid–solid mass transfer coefficient, s−1
k 0i Pre-exponential factor of the kinetic constant, mol (molPd s)−1
K 0i Pre-exponential factor of the equilibrium constant, m3 mol−1
n Pd Number of moles of palladium, mol
p Pressure, MPa
p H Hydrogen partial pressure, MPa
R Universal gas constant, kJ mol−1 K−1
r 0 Initial reaction rate, mol m−3 s−1
r i Reaction rate of species i, mol m−3 s−1
S E Selectivity of the process to MBE, −
T Temperature, K
V L Liquid volume, m3
X Y Conversion of MBY, −

Greek symbols

α 1, α2Parameters defined by eqn (1) and (2), −
σOverall percentage standard deviation, %
ΔEAbsorption activation energy, kJ mol−1
ΔHadiEnthalpy of adsorption, kJ mol−1
ρ S Density of catalyst particles, kg m−3

Abbreviations

DDimers
GGas
HHydrogen
LLiquid
MBA, A2-Methyl-2-butanol
MBE, E2-Methyl-3-buten-2-ol
MBY, Y2-Methyl-3-butyn-2-ol
RxRadical originating in step (e4)
SSolid
Y2Complex originating in step (e3)

Acknowledgements

The Scientific Center for Optical and Electron Microscopy of ETH Zurich and in particular Anne Greet Bittermann and Frank Krumeich are gratefully acknowledge for providing SEM/STEM facility access and for supporting the analysis. The authors acknowledge financial support from the Swiss Commission for Technology and Innovation (CTI) through grant Nr. 16706.1 PFIW-IW. The authors thank Robert Büchel for help with BET analysis and temperature programmed studies and Georges Siddiqui for valuable discussion concerning STEM analysis.

Notes and references

  1. F. Roessler, Chimia, 2003, 57, 791–798 CrossRef CAS.
  2. W. Bonrath, J. Medlock, J. Schütz, B. Wüstenberg and T. Netscher, in Hydrogenation, ed. I. Karamé, Intech Open, Rijeka, Croatia, 2012 Search PubMed.
  3. 3-buten-2-ol, 2-methyl-, CAS No. 115-18-4, in Chemicals Screening Information Data Set (SIDS) for High Volume Chemicals, UNEP Publications, Ispra, Italy, 1995 Search PubMed.
  4. H. Lindlar, Helv. Chim. Acta, 1952, 35, 446 CrossRef CAS.
  5. N. Semagina, M. Grasemann, N. Xanthopoulos, A. Renken and L. Kiwi-Minsker, J. Catal., 2007, 251, 213–222 CrossRef CAS.
  6. M. Garcia-Mota, J. Gomez-Diaz, G. Novell-Leruth, C. Vargas-Fuentes, L. Bellarosa, B. Bridier, J. Perez-Ramirez and N. Lopez, Theor. Chem. Acc., 2011, 128, 663–673 CrossRef CAS.
  7. A. Yarulin, I. Yuranov, F. Cardenas-Lizana, P. Abdulkin and L. Kiwi-Minsker, J. Phys. Chem. C, 2013, 117, 13424–13434 CAS.
  8. I. S. Mashkovsky, G. N. Baeva, A. Y. Stakheev, M. N. Vargaftik, N. Y. Kozitsyna and I. I. Moiseev, Mendeleev Commun., 2014, 24, 355–357 CrossRef CAS.
  9. M. Crespo-Quesada, M. Grasemann, N. Semagina, A. Renken and L. Kiwi-Minsker, Catal. Today, 2009, 147, 247–254 CrossRef CAS.
  10. T. O. Omarkulov, Z. Mukataev, S. V. Goncharova and D. V. Sokol'skii, Russ. J. Phys. Chem., 1987, 59, 137 Search PubMed.
  11. L. N. Protasova, E. V. Rebrov, K. L. Choy, S. Y. Pung, V. Engels, M. Cabaj, A. E. H. Wheatley and J. C. Schouten, Catal. Sci. Technol., 2011, 1, 768–777 Search PubMed.
  12. E. V. Ramos-Fernandez, A. F. P. Ferreira, A. Sepulveda-Escribano, F. Kapteijn and F. Rodriguez-Reinoso, J. Catal., 2008, 258, 52 CrossRef CAS.
  13. S. Vernuccio, P. Rudolf von Rohr and J. Medlock, Ind. Eng. Chem. Res., 2015, 54, 11543–11551 CrossRef CAS.
  14. P. A. Ramachandran and R. V. Chaudari, in Three phase catalytic reactors, Gordon and Breach Science Publishers, New York, 1983 Search PubMed.
  15. E. Dietrich, C. Mathieu, H. Delmas and J. Jenk, Chem. Eng. Sci., 1992, 47, 3597–3604 CrossRef CAS.
  16. J. A. Alves, S. P. Bressa, O. M. Martinez and G. F. Barreto, Chem. Eng. J., 2007, 125, 131–138 CrossRef CAS.
  17. S. P. Bressa, O. M. Martinez and G. F. Barreto, Ind. Eng. Chem. Res., 2003, 42, 2081–2092 CrossRef CAS.
  18. A. Molnar, A. Sarkany and M. Varga, J. Mol. Catal. A: Chem., 2001, 173, 185–221 CrossRef CAS.
  19. A. Bruehwiler, N. Semagina, M. Grassmann, A. Renken, L. Kiwi-Minsker, A. Saaler, H. Lehmann, W. Bonrath and F. Roessler, Ind. Eng. Chem. Res., 2008, 47, 6862–6869 CrossRef CAS.
  20. A. Sarkany, Z. Zsoldos, B. Furlong, J. W. Hightower and L. Guczi, J. Catal., 1993, 141, 566–582 CrossRef CAS.
  21. F. Studt, F. Abild-Pedersen, Th. Bligaard, R. Z. Sørensen, C. H. Christensen and J. K. Nørskov, Science, 2008, 320, 1320 CrossRef CAS PubMed.
  22. A. Yarulin, I. Yuranov, F. Cardenas-Lizana, D. T. L. Alexander and L. Kiwi-Minsker, Appl. Catal., A, 2014, 478, 186–193 CrossRef CAS.

Footnote

Electronic supplementary information (ESI) available: Physical properties of MBY, EDX spectra relative to SEM and STEM images, TPR and TPO profiles, summary of G–L and L–S mass transfer resistances. See DOI: 10.1039/c6re00093b

This journal is © The Royal Society of Chemistry 2016
Click here to see how this site uses Cookies. View our privacy policy here.