Supporting Information. First principles kinetic study on the effect of zeolite framework on 1-butanol dehydration
|
|
- Lesley Eaton
- 5 years ago
- Views:
Transcription
1 Supporting Information First principles kinetic study on the effect of zeolite framework on 1-butanol dehydration Mathew John, Konstantinos Alexopoulos, Marie-Françoise Reyniers* and Guy B. Marin Laboratory for Chemical Technology, Ghent University, Technologiepark 914, B-9052 Gent, Belgium * Corresponding author. Tel.: Fax: MarieFrancoise.Reyniers@ugent.be. S1
2 Table of Contents S1. Zeolite structure and unit cell parameter of H FAU, H ZSM 5, H ZSM 22 and H- FER optimized using the periodic DFT D method are shown in Figure S1 and Table S1. S2. Statistical thermodynamic calculation. S3. Thermodynamic correction for gas-phase species S4. Comparison with experimental observations: Comparison of simulated and experimental TOFs for butanol dehydration in H-ZSM-5 is shown in Table S2. S5. Comparison of immobile and mobile adsorbate approach: overview of mobile adsorbate approach and frequencies corresponding to translational and rotational modes shown in Table S3 and Table S4, respectively. Comparison of immobile and mobile adsorbate approach shown in Table S5 and Figure S2. S6. Ab initio based thermodynamic and kinetic parameters: Standard reaction enthalpy ( H r o ), reaction entropy ( S r o ), Arrhenius activation energies and pre-exponential factors for each of elementary step over H-FAU, H-ZSM-5, H-ZSM-22 and H-FER is listed in supporting information Tables S6, S7, S8 and S9, respectively S7. Adsorption of butanol/ether/water on zeolites: The adsorbed butanol monomer (M1), butanol dimer (D1), protonated ether (DBE*) and co-adsorbed butanol-water (C2) in H-FAU/ZSM-5/H-ZSM-22/H-FER zeolites is shown in Figure S3. S8. Transition state structures, TS-1 to TS-11 for all four zeolites are shown in Figure S4- S14, Comparison of geometric sizes of transition states for syn-elimination, antielimination and SN2 substitution is shown in Figure S15. S9. Standard free energy diagrams for direct dehydration (path A), ether formation (path B) and ether decomposition (path C) for H-FAU, H-ZSM-5, H-ZSM-22 and H- FER are shown in Figures S16-S27. The standard Gibbs free energy barriers for key activated steps corresponding to each path are listed in Table S10. Standard activation enthalpy ( H 0 ) and entropy ( S 0 ) for key activated steps corresponding to each path are shown in Figure S28 S10. Surface coverages: effect of 1-butanol conversion and feed butanol partial pressure on the surface coverage are shown in Figure S29 and Figure 30. Summary on the effect of reaction condition on the surface coverage is shown in Figure 31 S11. Analytical expression for the TOF of 1-butanol dehydration in zeolites. S12. References S2
3 S1. Zeolites Structure and unit cell parameters: Table S1. Unit cell parameters optimized using the periodic DFT D method for H FAU, H ZSM 5, H ZSM 22 and H FER. H FAU H ZSM 5 H ZSM 22 H FER a / pm b / pm c / pm α / β / γ / Figure S1: Framework structure and location of Brønsted acid site for H-FAU, H-ZSM-5, H-ZSM-22 and H-FER S3
4 S2. Statistical thermodynamic calculation The molecular partition function q for gas-phase species is calculated from statistical thermodynamics 1 as follows: q( V, T) q ( V, T) q ( T) q ( T) = (S.1) trans vib rot q vib, q trans and q rot correspond to vibrational, translational and rotational molecular partition function, respectively. Translational partition function q trans 3 2 2πmkBT ( V, T ) = V 2 h (S.2) Rotational partition function for non-linear molecule q rot 2 πi AI BIC 8π kbt ( T ) = 2 σ h 3 2 (S.3) Vibrational partition function considering the zero-point vibrational energy (ZPVE) to be the zero reference point q vib ( T ) = 3n 6 i 1 hν i 1 exp kbt (S.4) where, n is the number of atoms in the molecule, h is the Planck constant, k B is the Boltzmann factor, T is the temperature,ν i represent the vibrational frequencies obtained by a full hessian calculation with VASP, m is the molecular mass. For gas phase molecules, the volume V is calculated using an ideal gas equation. I A, I B, and I C are the principal moments of inertia, σ is symmetry number. S4
5 On the other hand, the molecular partition function for immobile surface species (q immobile ) is calculated from statistical thermodynamics as follows: q( T) q ( T) q ( T) = (S.5) immobile = vib q vib ( T ) = 3n i 1 hν i 1 exp kbt (S.6) where, n is the number of atoms considered in the partial hessian calculation, h is the Planck constant, k B is the Boltzmann factor, T is the temperature and ν i represent the vibrational frequencies obtained by a partial hessian calculation with VASP. The total ensemble partition function Q(N,V,T) is given as: [ q( V, T )] N Q( N, V, T ) = (S.7) N! where, N is the number of particles and q is the molecular partition function. Statistical thermodynamic calculation for internal energy (U), enthalpy (H), entropy (S) and Gibbs free energy (G) 1 U E DFT D2 + E ZPVE 2 lnq + kbt T = (S.8) H U+ PV= U+ RT N, V = (S.9) lnq S= kb lnq+ kbt (S.10) T N, V G = H TS (S.11) where, E DFT-D2 is energy from the DFT-D2 calculation, Q is total partition function, E ZPVE is the zero point mod es 1 vibrational energy calculated as the sum of all the molecular vibrational energies at 0K ( EZPVE = hν i ). 2 i S5
6 S3. Thermodynamic correction for gas-phase species The error in the calculated standard Gibbs free reaction energy can be ascribed to errors associated with the DFT method and the harmonic oscillator approximation when calculating the gas-phase reaction enthalpy and entropy. NIST database DFT-D2 Temperature: 298 K, Pressure: kpa H 0 S 0 G 0 H 0 S 0 G 0 kj/mol kj/mol/k kj/mol kj/mol J/mol/K kj/mol 1-butanol (g) 1-butene (g) + H 2 O (g) x 1-butanol (g) di-1-butyl ether (g) + H 2 O (g) di-1-butyl ether (g) 1-butene (g) +1-butanol (g) The following approach was used for the thermodynamic correction for the investigated gasphase reactions. The standard enthalpy of formation (H 0 ), standard entropy of formation (S 0 ) and heat capacity (Cp(T)) at different temperatures were obtained for each of the gas phase species (water, 1-butanol, 1-butene and di-butyl ether) from the NIST website 2. The temperature dependence of heat capacity were fitted to a third order polynomial which was then used to obtain enthalpy H (T) and entropy S (T) values at specific temperatures Cp (T) = a + bt + c't 2 + dt 3 (S.12) S T = + S.13 H(T) = H 0 + (S.14) S6
7 These values of H (T) and entropy S (T) were used for calculating the change in Gibbs free energy G (T) for each of the three reaction paths Path A (mechanism # m 1 -m 5 ) 1-butanol (g) 1-butene (g) + H 2 O (g) Path B (mechanism # m 6 -m 8 ) 2 x 1-butanol (g) di-1-butyl ether (g) + H 2 O (g) Path C (mechanism # m 9 -m 10 ) di-1-butyl ether (g) 1-butene (g) +1-butanol (g) The difference between the standard Gibbs free energy G 0 (T) DFT-D2 + harmonic oscillator approximation and the one calculated on the basis of NIST values 2 are used to calculate error associated with each of the paths (ε GpathA, ε GpathB and ε GpathC).. NIST database DFT-D2 Error dg 0 (500K) kj/mol dg 0 (500K) kj/mol dg 0 error (500K) 1-butanol (g) 1-butene (g) + H 2 O (g) x 1-butanol (g) di-1-butyl ether (g) + H 2 O (g) di-1-butyl ether (g) 1-butene (g) + 1-butanol (g) In order to eliminate this error along each path, the error is distributed amongst the adsorption/ desorption steps in each of the mechanisms. As seen from the following mechanisms, we need to solve a set of algebraic equations (equations S15-S24) in order to calculate errors associated with adsorption and desorption steps. As the number of independent equations is less than number of unknowns (error terms). Simplifying assumption were made to reduce number of unknowns. S7
8 The error for all steps involving butanol adsorption/desorption step are assumed to same. 1.e. error in elementary reaction involving adsorption of butanol on Zeolite forming butanol monomer (M1) is same as that of error in formation of dimer (D1) from (M1) A similar assumption was made for the adsorption/ desorption steps of other species. Although, G error may be different for different adsorbing species, but error associated with the adsorption/desorption of the same species is the same. For the activated steps that involved desorption of product, as most of the transition states the error associated with the gas-phase thermodynamic correction is assumed to be all present in the reverse step, with no change in the activation energy barrier and pre-exponential factor for the forward reaction. The distribution of error amongst the elementary steps involving adsorption/ desorption at 450 K is shown in Table S3-S6 of section S4. Comparison of the free energy diagram (for H-ZSM-5 at 500 K) for thermodynamic corrected (dotted lines) and the uncorrected (full line) case for different mechanisms are shown below : S8
9 S9
10 S10
11 S11
12 S12
13 S13
14 S4. Comparison with experimental observations. To validate the results of our ab initio based microkinetic model, a comparison is made with the experimental observation of Makarova et al. 3 for 1-butanol dehydration in H-ZSM-5. The microkinetic model was used to calculate the turnover frequency (TOF) for the production of dibutyl ether (DBE) and butene at their experimental conditions (400K, butanol feed mole fraction 0.7 %, total pressure of 1 bar and site time of 36.6 mol H+ s mol -1 ). Table S2 provides a comparison between current theoretical results with thermodynamic correction, our previous theoretical result without thermodynamic correction 4 and experimental TOFs for production of dibutyl ether (DBE) and butene. The inclusion of thermodynamic correction for gas-phase species (to account for errors associated with DFT and harmonic approximation) leads to significant improvement in results which are in line with the experimental results. Table S2. Comparison between current theoretical results with thermodynamic correction, our previous theoretical result without thermodynamic correction and experimental TOFs for production of dibutyl ether (DBE) and butene at 400K, butanol feed mole fraction 0.7%, total pressure of 1 bar and site time of 36.6 mol H+ s mol -1. This work Experimental 3 Ref 4 1 TOF Butene (mol /mol H+ /s) TOF DBE (mol /mol H+ /s) Conversion (mol %) 2.1 ~ S14
15 S5. Comparison of immobile and mobile adsorbate approach for the physisorbed 1- butene. As seen from Table S3, 2D free translational and 1D free rotational modes were identified for physisorbed 1-butene in H-FAU, H-ZSM-5 and H-ZSM-22 by visual analysis of the vibrational modes 5, while only 2D free translational modes could be identified for physisorbed 1-butene in H-FER (with free rotation being restricted due to the location of the acid site). The harmonic frequencies corresponding to these translational and rotational modes are shown in Table S4. In the mobile adsorbate approach, these frequencies are removed from the calculation of the vibrational partition and are replaced by free translational and/or rotational contributions, as described by De Moor et al. 5 Table S3: Overview of the number of free translations and rotations considered in mobile approach for physisorbed 1-butene in different zeolites free translation molecular surface area (pm pm) free rotation number zeolite number H-FAU H-ZSM H-ZSM H-FER Table S4: List of harmonic frequencies (cm -1 ) corresponding to free translational and rotational modes of physisorbed 1-butene in different zeolites translational mode rotational mode H-FAU i H-ZSM H-ZSM H-FER A comparison of the standard physisorption entropy obtained using the immobile and mobile approach is presented in Table S5. As expected, the calculated entropy losses are significantly lower if mobile physisorption is assumed. S15
16 Table S5: Immobile versus mobile physisorption entropies (J mol -1 K -1 ) of 1-butene in different zeolites zeolite Immobile approach Mobile approach H-FAU H-ZSM H-ZSM H-FER Nevertheless, as seen from Figure S2, no significant difference was observed for 1-butanol conversion, while a marginal increase in 1-butene yield is observed at higher site time (conversion) when using mobile adsorbate approach instead of the immobile adsorbate approach. Figure S2: Comparison of immobile and mobile adsorbate approach for the physisorbed 1- butene. (Full and dotted lines correspond to result obtained using immobile and mobile adsorbate approach, respectively) S16
17 S6. Ab initio derived Thermodynamic and kinetic parameters Table S6:. Standard reaction enthalpy (kj/mol), reaction entropy (J/mol/K), activation energy (kj/mol), preexponential factor (s -1 ), forward reaction rate coefficient k f (s -1 ) at 450K, thermodynamic correction term on standard Gibbs free energy at 450 K (kj/mol) and equilibrium coefficient at 450K (10-2 kpa -1, 10 2 kpa or dimensionless for adsorption, desorption and surface transformation, respectively) for the elementary steps (numbered as indicated in Scheme 1 ) for H-FAU. Elementary steps H r o S r o E a(f) A f G 0 r ε G(T) k f (T) K eq (T) corr (R1) 1-BuOH(g) + * M E+03 (R2) M1 W + 1-Butene(g) E E E+00 (R3) W H2O(g) + * E-02 (R4) M1 C E E E-06 (R5) C1 W + 1-Butene(g) E+06 (R6) M1 M E-08 (R7) M2 1-Butene*+ H2O(g) E E E+06 (R8) 1-Butene* 1-Butene(g) + * E+00 (R9) M2 Butoxy + H2O(g) E E E+05 (R10) Butoxy 1-Butene* E E E+01 (R11) M1 + BuOH(g) D E+02 (R12) D1 D E-14 (R13) D2 C2+1-Butene(g) E E E+12 (R14) C2 M1 + H2O(g) E+01 (R15) D2 DBE* + H2O(g) E E E+14 (R16) DBE* DBE(g) + * E-07 (R17) Butoxy + BuOH(g) C E-03 (R18) C3 DBE* (Sn2) E E E+08 (R19) C3 DBE* (Sn1) E E E+08 (R20) DBE* C E E E-10 (R21) C4 1-Butene*+ BuOH(g) E+04 (R22) DBE* DBE E-09 (R23) DBE2 1-Butene*+ BuOH(g) E E E+03 S17
18 Table S7: Standard reaction enthalpy (kj/mol), reaction entropy (J/mol/K), activation energy (kj/mol), preexponential factor (s -1 ), forward reaction rate coefficient k f (s -1 ) at 450K, thermodynamic correction term on standard Gibbs free energy at 450 K (kj/mol) and equilibrium coefficient at 450K (10-2 kpa -1, 10 2 kpa or dimensionless for adsorption, desorption and surface transformation, respectively) for the elementary steps (numbered as indicated in Scheme 1 ) for H-ZSM-5. (adapted from ref 4 including thermodynamic correction ) Elementary steps H r o S r o E a(f) A f G 0 r ε G(T) k f (T) K eq (T) corr (R1) 1-BuOH(g) + * M E+07 (R2) M1 W + 1-Butene(g) E E E-02 (R3) W H2O(g) + * E-02 (R4) M1 C E E E-05 (R5) C1 W + 1-Butene(g) E+02 (R6) M1 M E-10 (R7) M2 1-Butene*+ H2O(g) E E E+07 (R8) 1-Butene* 1-Butene(g) + * E-02 (R9) M2 Butoxy + H2O(g) E E E+06 (R10) Butoxy 1-Butene* E E E+01 (R11) M1 + BuOH(g) D E+05 (R12) D1 D E-06 (R13) D2 C2+1-Butene(g) E E E+00 (R14) C2 M1 + H2O(g) E+02 (R15) D2 DBE* + H2O(g) E E E+06 (R16) DBE* DBE(g) + * E-13 (R17) Butoxy + BuOH(g) C E+02 (R18) C3 DBE* (Sn2) E E E+08 (R19) C3 DBE* (Sn1) E E E+08 (R20) DBE* C E E E-10 (R21) C4 1-Butene*+ BuOH(g) E+00 (R22) DBE* DBE E-07 (R23) DBE2 1-Butene*+ BuOH(g) E E E-03 S18
19 Table S8:. Standard reaction enthalpy (kj/mol), reaction entropy (J/mol/K), activation energy (kj/mol), preexponential factor (s -1 ), forward reaction rate coefficient k f (s -1 ) at 450K, thermodynamic correction term on standard Gibbs free energy at 450 K (kj/mol) and equilibrium coefficient at 450K (10-2 kpa -1, 10 2 kpa or dimensionless for adsorption, desorption and surface transformation, respectively) for the elementary steps (numbered as indicated in Scheme 1 ) for H-ZSM-22. Elementary steps H r o S r o E a(f) A f G 0 r ε G(T) k f (T) K eq (T) corr (R1) 1-BuOH(g) + * M E+07 (R2) M1 W + 1-Butene(g) E E E-04 (R3) W H2O(g) + * E-01 (R4) M1 C E E E-09 (R5) C1 W + 1-Butene(g) E+04 (R6) M1 M E-06 (R7) M2 1-Butene*+ H2O(g) E E E+02 (R8) 1-Butene* 1-Butene(g) + * E-02 (R9) M2 Butoxy + H2O(g) E E E+02 (R10) Butoxy 1-Butene* E E E-01 (R11) M1 + BuOH(g) D E+04 (R12) D1 D E-11 (R13) D2 C2+1-Butene(g) E E E+01 (R14) C2 M1 + H2O(g) E-02 (R15) D2 DBE* + H2O(g) E E E+11 (R16) DBE* DBE(g) + * E-13 (R17) Butoxy + BuOH(g) C E-01 (R18) C3 DBE* (Sn2) E E E+09 (R19) C3 DBE* (Sn1) E E E+09 (R20) DBE* C E E E-08 (R21) C4 1-Butene*+ BuOH(g) E-01 (R22) DBE* DBE E-09 (R23) DBE2 1-Butene*+ BuOH(g) E E E-01 S19
20 Table S9: Standard reaction enthalpy (kj/mol), reaction entropy (J/mol/K), activation energy (kj/mol), preexponential factor (s -1 ), forward reaction rate coefficient k f (s -1 ) at 450K, thermodynamic correction term on standard Gibbs free energy at 450 K (kj/mol) and equilibrium coefficient at 450K (10-2 kpa -1, 10 2 kpa or dimensionless for adsorption, desorption and surface transformation, respectively) for the elementary steps (numbered as indicated in Scheme 1 ) for H-FER. Elementary steps H r o S r o E a(f) A f G 0 r ε G(T) k f (T) K eq (T) corr (R1) 1-BuOH(g) + * M E+00 (R2) M1 W + 1-Butene(g) E E E+02 (R3) W H2O(g) + * E-01 (R4) M1 C E E E-02 (R5) C1 W + 1-Butene(g) E+04 (R6) M1 M E-05 (R7) M2 1-Butene*+ H2O(g) E E E+02 (R8) 1-Butene* 1-Butene(g) + * E+04 (R9) M2 Butoxy + H2O(g) E E E+04 (R10) Butoxy 1-Butene* E E E-02 (R11) M1 + BuOH(g) D E+06 (R12) D1 D E-02 (R13) D2 C2+1-Butene(g) E E E-01 (R14) C2 M1 + H2O(g) E-05 (R15) D2 DBE* + H2O(g) E E E-01 (R16) DBE* DBE(g) + * E-05 (R17) Butoxy + BuOH(g) C E+00 (R18) C3 DBE* (Sn2) E E E+05 (R19) C3 DBE* (Sn1) E E E+05 (R20) DBE* C E E E-04 (R21) C4 1-Butene*+ BuOH(g) E-05 (R22) DBE* DBE E-02 (R23) DBE2 1-Butene*+ BuOH(g) E E E-06 S20
21 S7. Adsorption of butanol/ether/water on Zeolites Figure S3: The adsorbed butanol monomer (M1 ), butanol dimer (D1), protonated ether (DBE*) and co-adsorbed butanol-water (C2) in H-FAU/ZSM-5/H-ZSM-22/H-FER zeolites. S21
22 S8. Transition state structures (TS): Figure S4. Transition state (TS-1) for conversion of 1-butanol monomer (M1) to 1-butene and adsorbed water (W) via E1-elimination (step 2 of mechanism m 1 ) over H-FAU/H-ZSM- 5/H-ZSM-22/H-FER. Color code: silicon - light blue, aluminum - pink, oxygen - red, hydrogen - white, carbon - grey, hydrogen bonds (distance < 250 pm) blue lines, bond breaking/forming - black lines. S22
23 Figure S5. Transition state (TS-2) for conversion of M1 (1-butanol monomer) to C1 (butanol and water on zeolite ) via syn-elimination (step 4 of mechanism m 2 ) over H-FAU/H-ZSM- 5/H-ZSM-22/H-FER.Color code: silicon - light blue, aluminum - pink, oxygen - red, hydrogen - white, carbon - grey, hydrogen bonds (distance < 250 pm) blue lines, bond breaking/forming - black lines. S23
24 Figure S6. Transition state (TS-3) for conversion of reoriented 1-butanol monomer (M2) to 1- butene and H 2 O via anti-elimination (step 7 of mechanism m 3 ) over H-FAU/H-ZSM-5/H- ZSM-22/H-FER.Color code: silicon - light blue, aluminum - pink, oxygen - red, hydrogen - white, carbon - grey, hydrogen bonds (distance < 250 pm) blue lines, bond breaking/forming - black lines. S24
25 Figure S7: Transition state (TS-4) for conversion of reoriented 1-butanol monomer (M2) to 1-butoxide and H 2 O via SN2 type substitution ( step 9, common to mechanism m 4/ m 7 and m 8 ) over H-FAU/H-ZSM-5/H-ZSM-22/H-FER.Color code: silicon - light blue, aluminum - pink, oxygen - red, hydrogen - white, carbon - grey, hydrogen bonds (distance < 250 pm) blue lines, bond breaking/forming - black lines. S25
26 Figure S8: Transition state (TS-5) for deprotonation of 1-butoxide to 1-butene (step 10 of mechanism m 4, over H-FAU/H-ZSM-5/H-ZSM-22/H-FER. Color code: silicon - light blue, aluminum - pink, oxygen - red, hydrogen - white, carbon - grey, hydrogen bonds (distance < 250 pm) blue lines, bond breaking/forming - black lines. S26
27 Figure S9. Transition state (TS-6) for conversion of reoriented 1-butanol dimer (D2) to 1butene,co-adsorbed butanol and H 2 O via syn-elimination ( step 13 of mechanism m 5 ) over H-FAU/H-ZSM-5/H-ZSM-22/H-FER.Color code: silicon - light blue, aluminum - pink, oxygen - red, hydrogen - white, carbon - grey, hydrogen bonds (distance < 250 pm) blue lines, bond breaking/forming - black lines. S27
28 Figure S10. Transition state (TS-7) for conversion of reoriented 1-butanol dimer (D2) to di-1- butyl ether (DBE*) and H 2 O via SN2 type substitution ( step 15 of mechanism m 6 ) over H- FAU/H-ZSM-5/H-ZSM-22/H-FER.Color code: silicon - light blue, aluminum - pink, oxygen - red, hydrogen - white, carbon - grey, hydrogen bonds (distance < 250 pm) blue lines, bond breaking/forming - black lines. S28
29 Figure S11. Transition state (TS-8) for conversion of butanol and 1-butoxide (C3) to di-1- butyl ether (DBE*) via SN2 type substitution ( step 18 of mechanism m 7 ) over H-FAU/H- ZSM-5/H-ZSM-22/H-FER.Color code: silicon - light blue, aluminum - pink, oxygen - red, hydrogen - white, carbon - grey, hydrogen bonds (distance < 250 pm) blue lines, bond breaking/forming - black lines. S29
30 Figure S12. Transition state (TS-9) for conversion of butanol and 1-butoxide (C3) to di-1- butyl ether (DBE*) via SN1 type substitution (step 19 of mechanism m 8 ) over H-FAU/H- ZSM-5/H-ZSM-22/H-FER.Color code: silicon - light blue, aluminum - pink, oxygen - red, hydrogen - white, carbon - grey, hydrogen bonds (distance < 250 pm) blue lines, bond breaking/forming - black lines. S30
31 Path C ( TS-10) Figure S13: Transition state (TS-10) for conversion of adsorbed di-1-butyl ether (DBE*) to co-adsorbed 1-butene and butanol (C4) via syn-elimination (step 20 of mechanism m 9) over H- FAU/H-ZSM-5/H-ZSM-22/H-FER. Color code: silicon - light blue, aluminum - pink, oxygen - red, hydrogen - white, carbon - grey, hydrogen bonds (distance < 250 pm) blue lines, bond breaking/forming - black lines. S31
32 Figure S14: Transition state (TS-11) for conversion of re-oriented adsorbed di-1-butyl ether (DBE2) to adsorbed 1-butene and butanol (g) via anti-elimination (step 23 of mechanism m 10) over H-FAU/H-ZSM-5/H-ZSM-22/H-FER. Color code: silicon - light blue, aluminum - pink, oxygen - red, hydrogen - white, carbon - grey, hydrogen bonds (distance < 250 pm) blue lines, bond breaking/forming - black lines. S32
33 Figure S15 : Comparison of transition state geometric size for syn,/anti-elimination and SN2 reaction (all dimensions in pm) S33
34 S9. Standard free energy diagram for butanol dehydration on zeolites Energy Diagram for H-FAU zeolite at 450K Figure S16: Standard Gibbs free energy for 1-butanol dehydration to 1-butene (path A) in H-FAU at 450K Figure S17: Standard Gibbs free energy for 1-butanol dehydration to di-1-butyl ether (path B) in H-FAU at 450K Figure S18: Standard Gibbs free energy for ether decomposition reaction (path C) in H-FAU at 450K S34
35 Energy Diagram for H-ZSM-5 zeolite at 450K Figure S19: Standard Gibbs free energy for 1-butanol dehydration to 1-butene (path A) in H-ZSM-5 at 450K Figure S20: Standard Gibbs free energy for 1-butanol dehydration to di-1-butyl ether (path B) in H-ZSM-5 at 450K Figure S21: Standard Gibbs free energy for ether decomposition reaction (path C) in H-ZSM-5 at 450K S35
36 Energy Diagram for H-ZSM-22 zeolite at 450K Figure S22: Standard Gibbs free energy for 1-butanol dehydration to 1-butene (path A) in H-ZSM-22 at 450K Figure S23: Standard Gibbs free energy for 1-butanol dehydration to di-1-butyl ether (path B) in H-ZSM-22 at 450K Figure S24: Standard Gibbs free energy for ether decomposition reaction (path C) in H-ZSM-22 at 450K S36
37 Energy Diagram for H-FER zeolite at 450K Figure S25: Standard Gibbs free energy for 1-butanol dehydration to 1-butene (path A) in H-FER at 450K Figure S26: Standard Gibbs free energy for 1-butanol dehydration to di-1-butyl ether (path B) in H-FER at 450K Figure S27: Standard Gibbs free energy for ether decomposition reaction (path C) in H-FER at 450K S37
38 Table S10: The standard Gibbs free energy barriers for direct dehydration of butanol ( G 0 A = G 0 TS3-M1, ether formation ( G 0 B = G 0 TS7-D1 and ether decomposition ( G 0 C = G 0 TS10-DBE ) over H-FAU, H-ZSM-5, H-ZSM- 22 and H-FER. (see Figure 4 for the standard Gibbs free energy profile). G 0 M1ads G 0 A = G TS3 G M1 G 0 D1ads G 0 B = G TS7 G D1 G 0 C = G TS10 G DBE* (kj/mol) (kj/mol) (kj/mol) (kj/mol) (kj/mol) H-FAU H-ZSM H-ZSM H-FER Figure S28 : Standard activation enthalpy ( H 0 ) and entropy ( S 0 ) of direct dehydration of butanol (butanol monomer (M1) to TS-3, mechanism m 3 of path A), ether formation (butanol dimer (D1) to TS-7, mechanism m 6 of Path B) and ether decomposition path (adsorbed di butyl ether DBE to TS-7, mechanism m 9 of path C) over H-FAU, H-ZSM-5, H- ZSM-22 and H-FER. S38
39 S10. Surface coverages Figure S29: Effect of 1-butanol conversion on the surface coverage of adsorbed 1-butanol monomer M1 ( ), dimer D1 ( ), di-1-butyl ether DBE* ( ), water W ( ) and free sites θ* ( ) for different zeolites at reaction temperature of 450K and butanol partial pressure of 10 kpa (these results correspond to plug-flow reactor simulations (Eqs. 3-5) with the ab-initio based microkinetic model. See Tables S3, S4, S5 and S6, for H-ZSM-5, H-FAU, H-ZSM-22 and H-FER, respectively, for the model parameters) S39
40 Figure S30 : Effect of 1-butanol partial pressure on surface coverage of adsorbed 1-butanol monomer M1 ( ), dimer D1 ( ), di-1-butyl ether DBE* ( ), water W ( ) and free sites θ* ( ) for different zeolites at reaction temperature of 450K and 1-butanol conversion of 10%. (these results correspond to plug-flow reactor simulations (Eqs. 3-5) with the ab-initio based microkinetic model. See Tables S3, S4, S5 and S6, for H-ZSM-5, H-FAU, H-ZSM-22 and H- FER, respectively, for the model parameters) S40
41 Figure 31: Summary on the effect of reaction temperature and butanol partial pressure on the surface coverage of adsorbed 1-butanol monomer M1 ( ), dimer D1 ( ), di-1-butyl ether DBE* ( ) at constant conversion of 10 %.(these results correspond to plug-flow reactor simulations (Eqs. 3-5) with the ab-initio based microkinetic model. See Tables S3, S4, S5 and S6, for H-ZSM-5, H-FAU, H-ZSM-22 and H-FER, respectively, for the model parameters) S41
42 S11. Analytical expression for the TOF of 1-butanol dehydration in zeolites Derivation of simplified analytical equation for path A, path B and path C. Site balance: θ +θ +θ +θ +θ +θ +θ +θ +θ +θ +θ +θ =1 (S25) where, θ is the fraction coverage for i th surface species. A detailed microkinetic simulation for all the four zeolites indicates that M1, D1 and DBE* are the key surface species and depending upon the reaction condition one of them becomes the most abundant reaction intermediates (MARI). Moreover, under the simulated range of conditions (temperature K, 1-butanol partial pressure 10-3 to 100 kpa and conversion less than 30%) the summation of surface converges of M1, D1, DBE* and θ ~1. Accordingly, the site balance equation simplifies to: θ +θ +θ +θ =1 (S26) Based on quasi steady state assumption and site balance we can derive the values for surface coverages of each species as θ = K 1 P BuOH θ θ = K 1 K 11 (P BuOH ) 2 θ (S27) (S28) θ = K (P BuOH ) 2 θ (S29) Equation S26 becomes : θ +K 1 P BuOH θ + K 1 K 11 (P BuOH ) 2 θ + K (P BuOH ) 2 θ =1 (S30) Which gives θ = (S31) TOFs for each path are obtained by summing up the contributions from all the mechanisms associated with each of the paths. TOF, is obtained by summing up the contributions from each of the 5 mechanisms (m 1 - m 5 ) and is given as : TOF = TOF +TOF +TOF +TOF +TOF (S32) TOF, is obtained by summing up the contributions from each of the 3 mechanisms (m 6 - m 7 ) and is given as : TOF = TOF +TOF +TOF (S33) S42
43 As the mechanism m 6 remain dominant for most of the reaction conditions the equation becomes: TOF = TOF (S34) TOF, is obtained by summing up the contributions from each of the 3 mechanisms (m 6 - m 8 ) and is given as : TOF = TOF +TOF +TOF (S35) As the mechanism m 6 remain dominant for most of the reaction conditions the equation becomes: TOF = TOF (S36) TOF, is obtained by summing up the contributions from mechanism m 9 and m 10 is given as : TOF = TOF +TOF (S37) A simplified form of equation can be derived for conversion x= 0, as the contribution from backward (reverse ) reaction is zero TOF =0 =k θ = K k P θ TOF =0 =k θ = K k P θ TOF =0 =k θ = K K k P θ TOF =0 =γk θ =γ K K k P θ (S38) (S39) (S40) (S41) where γ is the fraction of butoxide that gets converted into butene and is given as the ratio of the TOF for formation of butene and ether from butoxide (γ= ). TOF =0 =k θ = K K K k P θ θ (S42) Thus the TOF at zero conversion is given as : TOF at X= 0 = K k +K k +K K k +K K k γ P +K K K k P 1+ K P + K K P S43 Similar approach is used for formulating TOF expression for path B and path C. S43
44 These simplified analytical expression can serve as a tool for easy comparison of theoretical TOFs with that obtained from experimental measurements of initial reaction rates. At extremely low P BuOH and high temperature (K 1 and K 11 decrease with increase in temperature), free sites are much more abundant than M1 and D1, i.e. K K P + K 1 P BuOH < 1, and a first order dependence of TOF path A on P BuOH is expected (based on equation S.43) : TOF = K k +K k +K K k +K K k γ P S.44 A region of positive order pressure dependence for path A is seen for H-FAU and H-FER (see Figure 3c ) at very low butanol partial pressure (P BuOH below 10-2 kpa). At moderate and high P BuOH, K K P and K 1 P BuOH > > 1, leading to: TOF = k +k +K k + K k γ +K K k P 1+ K P S.45 Here, a negative reaction order is seen with increasing butanol partial pressure (ascribed to the P independent term in the numerator corresponding to monomolecular path A mechanisms) followed by a regime of constant TOF = K k. A negative order pressure dependence is seen for all four zeolites at moderate to high pressure (see Figure 3c) and is consistent with the experimental observation of Chiang and Bhan 6 for ethanol dehydration within zeolites and of Macht et al. for 2-butanol dehydration on POM 7. Within the range of our simulation, a zero order pressure dependence is clearly seen for H-ZSM-5 above 10 kpa butanol partial pressure (see Figure 3c). Such a profile has been reported by Zhi et al. for 1-propanol dehydration in H-ZSM-5 8. S44
45 For DBE formation (path B), the initial TOF for formation of DBE is essentially determined by the dimer mediated mechanism m 6. Accordingly, the initial (at conversion, X= 0) TOF for path B is given as: TOF at X= 0 = K K K k P 1+K P +K K P S.46 When (K K P 1+K P ), a positive order pressure dependence is possible and is seen in the low pressure regime for H-FAU and H-FER (P BuOH below 10 2 kpa) (see Figure 3c). Typically, DBE formation is favored at high pressure (K K P > >1+ K P ) and shows a zero order dependence with TOF = K k. A zero order dependence regime for path B is seen for H-ZSM-5 and H-ZSM-22 at butanol partial pressures below 1kPa and at moderate butanol partial pressures (between 1-50 kpa) for H- FER (see Figure 3c). A zero order is also reported for zeolite catalyzed conversion of methanol to dimethyl ether 9. For DBE decomposition (path C), the TOF for path C (TOF is obtained by summing up TOF contributions from each of the 2 mechanisms m 9 and m 10 and is given as: TOF = K k +K k P 1+ K P + K K P + K P S.47 Here K is the equilibrium coefficient for elementary step 22 (DBE* DBE2) and k and k are forward reaction rate coefficients for elementary steps 20 (m 9 via syn-elimination) and 23 (m 10 via anti-elimination), respectively. A simplified form of equation (S.22) is obtained under conditions when DBE* is MARI (i.e. 1+ K P + K K P < K P, see Figure S31 for H-ZSM-5 and H- ZSM-22 in a butanol partial pressure range of 10-2 to 1 kpa). Under these conditions a zero order dependence is seen with TOF = k +K k, which is observed in Figure 3c S45
46 for H-ZSM-5 and H-ZSM-22. A more complex dependence on 1-butanol partial pressure is seen (see Figure 3c), when DBE is not MARI (see Figure S30 for surface coverages). When the TOFs of DBE formation and DBE decomposition are comparable (see Figure 3c, TOFs of path B and C for H-FER), the TOF for DBE decomposition shows the same pressure dependence as that of DBE formation path. At high butanol partial pressures (above 50 kpa), when K P + K K P 1+ K P ), a negative order dependence on 1-butanol partial pressure is seen. S46
47 S10. References: (1) Cramer, C. J. Essentials of computational chemistry : theories and models; 2nd ed.; Wiley: Chichester Hoboken, NJ, (2) NIST Chemistry WebBook, NIST Standard Reference Database Number 69; National Institute of Standards and Technology, (3) Makarova, M. A.; Paukshtis, E. A.; Thomas, J. M.; Williams, C.; Zamaraev, K. I. J Catal 1994, 149, (4) John, M.; Alexopoulos, K.; Reyniers, M. F.; Marin, G. B. J Catal 2015, 330, (5) De Moor, B. A.; Reyniers, M. F.; Marin, G. B. Phys Chem Chem Phys 2009, 11, (6) Chiang, H.; Bhan, A. J Catal 2010, 271, (7) Macht, J.; Janik, M. J.; Neurock, M.; Iglesia, E. J Am Chem Soc 2008, 130, (8) Zhi, Y.; Shi, H.; Mu, L.; Liu, Y.; Mei, D.; Camaioni, D. M.; Lercher, J. A. J Am Chem Soc 2015, 137, (9) Jones, A. J.; Iglesia, E. Angew Chem Int Edit 2014, 53, S47
Supporting Information. DFT Study of Methane Synthesis from Syngas on Ce Doped Ni(111) Surface
Supporting Information DFT Study of Methane Synthesis from Syngas on Ce Doped Ni(111) Surface Kai Li, 1 Cong Yin, 2 Yi Zheng, 3 Feng He, 1 Ying Wang, 1 Menggai Jiao, 1 Hao Tang, 2,* Zhijian Wu 1,* 1 State
More informationFirst-principles based design of Pt- and Pd-based catalysts for benzene hydrogenation
1 1 First-principles based design of Pt- and Pd-based catalysts for benzene hydrogenation Maarten K. Sabbe, Gonzalo Canduela, Marie- Françoise Reyniers, Guy B. Marin Introduction: benzene hydrogenation
More informationExpress the transition state equilibrium constant in terms of the partition functions of the transition state and the
Module 7 : Theories of Reaction Rates Lecture 33 : Transition State Theory Objectives After studying this Lecture you will be able to do the following. Distinguish between collision theory and transition
More informationFoundations of Chemical Kinetics. Lecture 12: Transition-state theory: The thermodynamic formalism
Foundations of Chemical Kinetics Lecture 12: Transition-state theory: The thermodynamic formalism Marc R. Roussel Department of Chemistry and Biochemistry Breaking it down We can break down an elementary
More informationThermodynamics and Kinetics
Thermodynamics and Kinetics C. Paolucci University of Notre Dame Department of Chemical & Biomolecular Engineering What is the energy we calculated? You used GAMESS to calculate the internal (ground state)
More informationKinetic and Theoretical Insights into the Mechanism of Alkanol Dehydration on Solid Brønsted Acid Catalysts
pubs.acs.org/jpcc Kinetic and Theoretical Insights into the Mechanism of Alkanol Dehydration on Solid Brønsted Acid Catalysts William Knaeble and Enrique Iglesia* Department of Chemical and Biomolecular
More informationUnusual Entropy of Adsorbed Methane on Zeolite Templated Carbon. Supporting Information. Part 2: Statistical Mechanical Model
Unusual Entropy of Adsorbed Methane on Zeolite Templated Carbon Supporting Information Part 2: Statistical Mechanical Model Nicholas P. Stadie*, Maxwell Murialdo, Channing C. Ahn, and Brent Fultz W. M.
More informationDocument Version Publisher s PDF, also known as Version of Record (includes final page, issue and volume numbers)
On the temperature dependence of the Arrhenius activation energy for hydroisomerization catalyzed by Pt/Mordenite Runstraat, van de, A.; van Grondelle, J.; van Santen, R.A. Published in: Journal of Catalysis
More informationReaction Dynamics (2) Can we predict the rate of reactions?
Reaction Dynamics (2) Can we predict the rate of reactions? Reactions in Liquid Solutions Solvent is NOT a reactant Reactive encounters in solution A reaction occurs if 1. The reactant molecules (A, B)
More informationTopic 4 Thermodynamics
Topic 4 Thermodynamics Thermodynamics We need thermodynamic data to: Determine the heat release in a combustion process (need enthalpies and heat capacities) Calculate the equilibrium constant for a reaction
More informationCatalytic Consequences of Spatial Constraints and Acid Site Location for Monomolecular Alkane Activation on Zeolites
Catalytic Consequences of Spatial Constraints and Acid Site Location for Monomolecular Alkane Activation on Zeolites Rajamani Gounder and Enrique Iglesia Supporting Information S.1. 27 Al MAS NMR Spectroscopy
More informationRole of products and intermediates in bioethanol conversion to hydrocarbons on H-ZSM-5: A time-resolved study
Role of products and intermediates in bioethanol conversion to hydrocarbons on H-ZSM-5: A time-resolved study Rakesh Batchu, Vladimir V. Galvita, Konstantinos Alexopoulos, Kristof Van der Borght, Hilde
More informationSupporting Information: Ethene Oligomerization in Ni-containing Zeolites: Theoretical Discrimination of Reaction. Mechanisms
Supporting Information: Ethene Oligomerization in Ni-containing Zeolites: Theoretical Discrimination of Reaction Mechanisms Rasmus Y. Brogaard and Unni Olsbye Department of Chemistry, Centre for Materials
More informationTransition Theory Abbreviated Derivation [ A - B - C] # E o. Reaction Coordinate. [ ] # æ Æ
Transition Theory Abbreviated Derivation A + BC æ Æ AB + C [ A - B - C] # E A BC D E o AB, C Reaction Coordinate A + BC æ æ Æ æ A - B - C [ ] # æ Æ æ A - B + C The rate of reaction is the frequency of
More informationSupporting Online Material (1)
Supporting Online Material The density functional theory (DFT) calculations were carried out using the dacapo code (http://www.fysik.dtu.dk/campos), and the RPBE (1) generalized gradient correction (GGA)
More informationGibbs Free Energy. Evaluating spontaneity
Gibbs Free Energy Evaluating spontaneity Predicting Spontaneity An increase in entropy; Changing from a more structured to less structured physical state: Solid to liquid Liquid to gas Increase in temperature
More informationGherman Group Meeting. Thermodynamics and Kinetics and Applications. June 25, 2009
Gherman Group Meeting Thermodynamics and Kinetics and Applications June 25, 2009 Outline Calculating H f, S, G f Components which contribute to H f, S, G f Calculating ΔH, ΔS, ΔG Calculating rate constants
More informationDefinitions and Concepts
2 Definitions and Concepts It is important that precise and unambiguous terms be used when dealing with rates of reaction and reaction modeling of a chemical system. Many of the definitions provided here
More informationCH361/361H Week 7 Lecture Gas Chromatography and Thermodynamics
CH36/36H Week 7 Lecture Gas Chromatography and Thermodynamics Experiment II Overview: Preparation & Equilibration of Isomeric 2,3-Dimethylbutenes 2 3 4 H + cat. 2 3 4 K eq =? ΔG =? ΔH =? ΔS =? 2,3-dimethyl-2-butene
More informationMechanistic Consequences of Composition in Acid Catalysis by Polyoxometalate Keggin Clusters
Published on Web 7/1/28 Mechanistic Consequences of Composition in Acid Catalysis by Polyoxometalate Keggin Clusters Josef Macht, Michael J. Janik, Matthew Neurock, and Enrique Iglesia*, Department of
More informationInsights into Different Products of Nitrosobenzene and Nitrobenzene. Hydrogenation on Pd(111) under the Realistic Reaction Condition
Insights into Different Products of Nitrosobenzene and Nitrobenzene Hydrogenation on Pd(111) under the Realistic Reaction Condition Lidong Zhang a, Zheng-Jiang Shao a, Xiao-Ming Cao a,*, and P. Hu a,b,*
More informationEfficient Synthesis of Ethanol from CH 4 and Syngas on
Efficient Synthesis of Ethanol from CH 4 and Syngas on a Cu-Co/TiO 2 Catalyst Using a Stepwise Reactor Zhi-Jun Zuo 1, Fen Peng 1,2, Wei Huang 1,* 1 Key Laboratory of Coal Science and Technology of Ministry
More informationHONOUR SCHOOL OF NATURAL SCIENCE. Final Examination GENERAL PHYSICAL CHEMISTRY I. Answer FIVE out of nine questions
HONOUR SCHOOL OF NATURAL SCIENCE Final Examination GENERAL PHYSICAL CHEMISTRY I Monday, 12 th June 2000, 9.30 a.m. - 12.30 p.m. Answer FIVE out of nine questions The numbers in square brackets indicate
More informationKinetic modeling in heterogeneous catalysis.
Kinetic modeling in heterogeneous catalysis. C. Franklin Goldsmith Modern Methods in Heterogeneous Catalysis Research Friday, January 27 th, 2012 Microkinetic modeling is a valuable tool in catalysis.
More informationCrystallographic Dependence of CO Activation on Cobalt Catalysts: HCP versus FCC
Crystallographic Dependence of CO Activation on Cobalt Catalysts: HCP versus FCC Jin-Xun Liu, Hai-Yan Su, Da-Peng Sun, Bing-Yan Zhang, and Wei-Xue Li* State Key Laboratory of Catalysis, Dalian Institute
More informationMechanisms of H- and OH-assisted CO activation as well as C-C coupling on the flat Co(0001) surface Revisited
Electronic Supplementary Material (ESI) for Catalysis Science & Technology. This journal is The Royal Society of Chemistry 2016 Mechanisms of H- and OH-assisted CO activation as well as C-C coupling on
More informationStatistical Thermodynamics. Lecture 8: Theory of Chemical Equilibria(I)
Statistical Thermodynamics Lecture 8: Theory of Chemical Equilibria(I) Chemical Equilibria A major goal in chemistry is to predict the equilibria of chemical reactions, including the relative amounts of
More informationModelling of Reaction Mechanisms KJE-3102
Modelling of Reaction Mechanisms KJE-3102 Kathrin H. Hopmann Kathrin.hopmann@uit.no Outline Potential energy surfaces Transition state optimization Reaction coordinates Imaginary frequencies Verification
More informationConsequences of Confinement in Zeolite Acid Catalysis. Rajamani Pachayappan Gounder. A dissertation submitted in partial satisfaction of the
Consequences of Confinement in Zeolite Acid Catalysis By Rajamani Pachayappan Gounder A dissertation submitted in partial satisfaction of the requirements for the degree of Doctor of Philosophy in Chemical
More informationSupplementary Information for:
Supplementary Information for: Towards Active and Stable Oxygen Reduction Cathode: A Density Functional Theory Survey on Pt 2 M skin alloys Guang-Feng Wei and Zhi-Pan Liu* Shanghai Key Laboratory of lecular
More informationStatistical thermodynamics for MD and MC simulations
Statistical thermodynamics for MD and MC simulations knowing 2 atoms and wishing to know 10 23 of them Marcus Elstner and Tomáš Kubař 22 June 2016 Introduction Thermodynamic properties of molecular systems
More informationCHAPTER 21: Reaction Dynamics
CHAPTER 21: Reaction Dynamics I. Microscopic Theories of the Rate Constant. A. The Reaction Profile (Potential Energy diagram): Highly schematic and generalized. A---B-C B. Collision Theory of Bimolecular
More informationSome properties of the Helmholtz free energy
Some properties of the Helmholtz free energy Energy slope is T U(S, ) From the properties of U vs S, it is clear that the Helmholtz free energy is always algebraically less than the internal energy U.
More informationChapter 12 What Is Temperature? What Is Heat Capacity?
Chapter 2 What Is Temperature? What Is Heat Capacity?. Heat capacity peaks from phase transitions. The peak heat capacity in the figure below shows that helium gas adsorbed on graphite undergoes a transition
More informationUnderstanding Chemical Reactions through Computer Modeling. Tyler R. Josephson University of Delaware 4/21/16
Understanding Chemical Reactions through Computer Modeling Tyler R. Josephson University of Delaware 4/21/16 A little about me B.S. in Chem E from U of M, 2011 Currently, Ph.D. student at University of
More information5.62 Physical Chemistry II Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 5.62 Physical Chemistry II Spring 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 5.62 Spring 2007 Lecture
More informationIntro/Review of Quantum
Intro/Review of Quantum QM-1 So you might be thinking I thought I could avoid Quantum Mechanics?!? Well we will focus on thermodynamics and kinetics, but we will consider this topic with reference to the
More information3. RATE LAW AND STOICHIOMETRY
Page 1 of 39 3. RATE LAW AND STOICHIOMETRY Professional Reference Shelf R3.2 Abbreviated Lecture Notes Full Lecture Notes I. Overview II. Introduction A. The Transition State B. Procedure to Calculate
More informationTowards More Complete Descriptors of Reactivity in Catalysis by Solid Acids
Supporting Information Towards More Complete Descriptors of Reactivity in Catalysis by Solid Acids Prashant Deshlahra 1 1, 2, * and Enrique Iglesia 1 Department of Chemical Engineering, University of California
More information[ A] 2. [ A] 2 = 2k dt. [ A] o
Chemistry 360 Dr Jean M Standard Problem Set 3 Solutions The reaction 2A P follows second-order kinetics The rate constant for the reaction is k350 0 4 Lmol s Determine the time required for the concentration
More informationTransition State Enthalpy and Entropy Effects on Reactivity. and Selectivity in Hydrogenolysis of n-alkanes
Transition State Enthalpy and Entropy Effects on Reactivity and Selectivity in Hydrogenolysis of n-alkanes David W. Flaherty, Enrique Iglesia * Department of Chemical Engineering, University of California
More informationIntro/Review of Quantum
Intro/Review of Quantum QM-1 So you might be thinking I thought I could avoid Quantum Mechanics?!? Well we will focus on thermodynamics and kinetics, but we will consider this topic with reference to the
More informationBrønsted/Lewis Acid Synergistically Promote the Initial C-C Bond Formation in MTO Reaction
Electronic Supplementary Material (ESI) for Chemical Science. This journal is The Royal Society of Chemistry 2018 Supporting Information for Brønsted/Lewis Acid Synergistically Promote the Initial C-C
More informationCH1101 Physical Chemistry Tutorial 1. Prof. Mike Lyons.
CH111 Physical Chemistry Tutorial 1. Prof. Mike Lyons. CH111 Section A Annual 1 Internal Energy Units: Joules J Internal Energy (U) : total kinetic & potential energy of system. e.g. Gas in container with
More information) 3 2. The nuclear and electronic partition functions are both neglected, and the total partition function q can be written as below:
Supporting Information Stable Pd-Doped Ceria Structures for CH 4 Activation and CO Oxidation Ya-Qiong Su, Ivo A.W. Filot, Jin-Xun Liu, and Emiel J.M. Hensen* Laboratory of Inorganic Materials Chemistry,
More informationMolecular Interactions F14NMI. Lecture 4: worked answers to practice questions
Molecular Interactions F14NMI Lecture 4: worked answers to practice questions http://comp.chem.nottingham.ac.uk/teaching/f14nmi jonathan.hirst@nottingham.ac.uk (1) (a) Describe the Monte Carlo algorithm
More informationDirect Synthesis of H 2 O 2 on AgPt Octahedra: The Importance of Ag-Pt Coordination for High H 2 O 2 Selectivity
Supporting Information Direct Synthesis of H 2 O 2 on AgPt Octahedra: The Importance of Ag-Pt Coordination for High H 2 O 2 Selectivity Neil M. Wilson, 1 Yung-Tin Pan, 1 Yu-Tsun Shao, 2 Jian-Min Zuo, 2
More informationThermal and Statistical Physics Department Exam Last updated November 4, L π
Thermal and Statistical Physics Department Exam Last updated November 4, 013 1. a. Define the chemical potential µ. Show that two systems are in diffusive equilibrium if µ 1 =µ. You may start with F =
More information1+e θvib/t +e 2θvib/T +
7Mar218 Chemistry 21b Spectroscopy & Statistical Thermodynamics Lecture # 26 Vibrational Partition Functions of Diatomic Polyatomic Molecules Our starting point is again the approximation that we can treat
More informationSupplementary Information. Rational Screening Low-Cost Counter Electrodes for Dye-Sensitized Solar Cells
Supplementary Information Rational Screening Low-Cost Counter Electrodes for Dye-Sensitized Solar Cells Yu Hou, Dong Wang, Xiao Hua Yang, Wen Qi Fang, Bo Zhang, Hai Feng Wang, Guan Zhong Lu, P. Hu, Hui
More informationCatalyst structure and C-O activation during FTS: new ideas from computational catalysis. Mark Saeys
Catalyst structure and C-O activation during FTS: new ideas from computational catalysis Mark Saeys Laboratory for Chemical Technology, Ghent University Laboratory for Chemical http://www.lct.ugent.be
More informationSupporting Information. Water-Gas Shift Activity of Atomically Dispersed Cationic Platinum versus Metallic Platinum Clusters on Titania Supports
Supporting Information Water-Gas Shift Activity of Atomically Dispersed Cationic Platinum versus Metallic Platinum Clusters on Titania Supports Salai Cheettu Ammal and Andreas Heyden * Department of Chemical
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION Using first principles to predict bimetallic catalysts for the ammonia decomposition reaction Danielle A. Hansgen, Dionisios G. Vlachos, Jingguang G. Chen SUPPLEMENTARY INFORMATION.
More informationCHEM 4641 Fall questions worth a total of 32 points. Show your work, except on multiple-choice questions. 1 V α=
Physical Chemistry I Final Exam Name: KEY CHEM 4641 Fall 017 15 questions worth a total of 3 points. Show your work, except on multiple-choice questions. 1 V 1 V α= κt = V T P V P T Gas constant R = 8.314
More informationChemistry 123: Physical and Organic Chemistry Topic 2: Thermochemistry
Recall the equation. w = -PΔV = -(1.20 atm)(1.02 L)( = -1.24 10 2 J -101 J 1 L atm Where did the conversion factor come from? Compare two versions of the gas constant and calculate. 8.3145 J/mol K 0.082057
More informationFIRST PUBLIC EXAMINATION. Trinity Term Preliminary Examination in Chemistry SUBJECT 3: PHYSICAL CHEMISTRY. Time allowed: 2 ½ hours
FIRST PUBLIC EXAMINATION Trinity Term 004 Preliminary Examination in Chemistry SUBJECT 3: PHYSICAL CHEMISTRY Wednesday, June 9 th 004, 9.30 a.m. to 1 noon Time allowed: ½ hours Candidates should answer
More informationSupporting Information. Heterostructures of MXene and N-doped graphene as highly. active bifunctional electrocatalysts
Electronic Supplementary Material (ESI) for Nanoscale. This journal is The Royal Society of Chemistry 2018 Supporting Information Heterostructures of MXene and N-doped graphene as highly active bifunctional
More informationAdvanced Physical Chemistry CHAPTER 18 ELEMENTARY CHEMICAL KINETICS
Experimental Kinetics and Gas Phase Reactions Advanced Physical Chemistry CHAPTER 18 ELEMENTARY CHEMICAL KINETICS Professor Angelo R. Rossi http://homepages.uconn.edu/rossi Department of Chemistry, Room
More informationwhere R = universal gas constant R = PV/nT R = atm L mol R = atm dm 3 mol 1 K 1 R = J mol 1 K 1 (SI unit)
Ideal Gas Law PV = nrt where R = universal gas constant R = PV/nT R = 0.0821 atm L mol 1 K 1 R = 0.0821 atm dm 3 mol 1 K 1 R = 8.314 J mol 1 K 1 (SI unit) Standard molar volume = 22.4 L mol 1 at 0 C and
More informationAdsorption of gases on solids (focus on physisorption)
Adsorption of gases on solids (focus on physisorption) Adsorption Solid surfaces show strong affinity towards gas molecules that it comes in contact with and some amt of them are trapped on the surface
More information= 16! = 16! W A = 3 = 3 N = = W B 3!3!10! = ΔS = nrln V. = ln ( 3 ) V 1 = 27.4 J.
Answer key: Q1A Both configurations are equally likely because the particles are non-interacting (i.e., the energy does not favor one configuration over another). For A M = 16 N = 6 W A = 16! 0.9 101 =
More informationJournal of Catalysis
Journal of Catalysis 278 (211) 78 93 Contents lists available at ScienceDirect Journal of Catalysis journal homepage: www.elsevier.com/locate/jcat Catalytic consequences of acid strength in the conversion
More informationChemical Kinetics. Topic 7
Chemical Kinetics Topic 7 Corrosion of Titanic wrec Casón shipwrec 2Fe(s) + 3/2O 2 (g) + H 2 O --> Fe 2 O 3.H 2 O(s) 2Na(s) + 2H 2 O --> 2NaOH(aq) + H 2 (g) Two examples of the time needed for a chemical
More informationSupplementary Information. The role of copper particle size in low pressure methanol synthesis via CO 2 hydrogenation over Cu/ZnO catalysts
Electronic Supplementary Material (ESI) for Catalysis Science & Technology. This journal is The Royal Society of Chemistry 2014 Supplementary Information The role of copper particle size in low pressure
More informationCBE 450 Chemical Reactor Fundamentals Fall, 2009 Homework Assignment #1
CBE 45 Chemical Reactor Fundamentals Fall, 9 omework Assignment #1 1. Review of Mass Balances of Reacting Systems atural gas (assume pure methane) is being burned in air (assume.79 mol% and.1 mol% O ).
More information= k 2 [CH 3 *][CH 3 CHO] (1.1)
Answers to Exercises Last update: Tuesday 29 th September, 205. Comments and suggestions can be sent to i.a.w.filot@tue.nl Exercise d[ch 4 ] = k 2 [CH 3 *][CH 3 CHO].) The target is to express short-lived
More informationMicroscopic Treatment of the Equilibrium Constant. Lecture
Microscopic Treatment of the Equilibrium Constant Lecture The chemical potential The chemical potential can be expressed in terms of the partition function: μ = RT ln Q j j N j To see this we first expand
More informationAn Introduction to Chemical Kinetics
An Introduction to Chemical Kinetics Michel Soustelle WWILEY Table of Contents Preface xvii PART 1. BASIC CONCEPTS OF CHEMICAL KINETICS 1 Chapter 1. Chemical Reaction and Kinetic Quantities 3 1.1. The
More information09 Intro to Mass Dependent Fractionation
09 Intro to Mass Dependent Fractionation Reading: White #26 Valley and Cole, Chapter 1 Guide Questions: 1) In a general sense why do heavier isotopes of an element behave differently from lighter isotopes?
More informationKINETICS PROBLEMS. (Alberty R A, Physical Chemistry, 7 th ed., Wiley, NY, 1987) (Answers on the last page)
KINETICS PROBLEMS (Alberty R A, Physical Chemistry, 7 th ed., Wiley, NY, 1987) (Answers on the last page) 20.1 The half-life of a first-order chemical reaction A --> B is 10 min. What percent of A remains
More informationTutorial 1 (not important for 2015)
Tutorial 1 (not important for 2015) 1 st Law of thermodynamics and other basic concepts Do No. 5 (05-03-2015) 1. One mole of an ideal gas is allowed to expand against a piston which supports 41 atm pressures.
More informationLecture 20. Chemical Potential
Lecture 20 Chemical Potential Reading: Lecture 20, today: Chapter 10, sections A and B Lecture 21, Wednesday: Chapter 10: 10 17 end 3/21/16 1 Pop Question 7 Boltzmann Distribution Two systems with lowest
More informationTOPIC 6: Chemical kinetics
TOPIC 6: Chemical kinetics Reaction rates Reaction rate laws Integrated reaction rate laws Reaction mechanism Kinetic theories Arrhenius law Catalysis Enzimatic catalysis Fuente: Cedre http://loincognito.-iles.wordpress.com/202/04/titanic-
More informationWashington University in St. Louis Chemistry Tournament Sample Problems for Individual Round #3: Kinetics, Electrochemistry, and Thermodynamics
Individual Exam #3: Kinetics, Electrochemistry, and Thermodynamics 1) A galvanic cell is to be constructed using the following half reactions under acidic conditions: ε 1.49 V ε 1.66 V Given that the potentials
More informationSelectivity in the initial C-H bond cleavage of n-butane on PdO(101)
Supporting Information for Selectivity in the initial C-H bond cleavage of n-butane on PdO(101) Can Hakanoglu (a), Feng Zhang (a), Abbin Antony (a), Aravind Asthagiri (b) and Jason F. Weaver (a) * (a)
More informationChemical reactions as network of rare events: Kinetic MonteCarlo
Chemical reactions as network of rare events: Kinetic MonteCarlo Extending the scale Length (m) 1 10 3 Potential Energy Surface: {Ri} 10 6 (3N+1) dimensional 10 9 E Thermodynamics: p, T, V, N continuum
More informationPHYSICAL CHEMISTRY CHEM330
PHYSICAL CHEMISTRY CHEM330 Duration: 3 hours Total Marks: 100 Internal Examiner: External Examiner: Professor B S Martincigh Professor J C Swarts University of the Free State INSTRUCTIONS: 1. Answer five
More informationReaction rate. reaction rate describes change in concentration of reactants and products with time -> r = dc j
Reaction rate ChE 400 - Reactive Process Engineering reaction rate describes change in concentration of reactants and products with time -> r = dc j /dt r is proportional to the reactant concentrations
More informationYou MUST sign the honor pledge:
CHEM 3411 MWF 9:00AM Fall 2010 Physical Chemistry I Exam #2, Version B (Dated: October 15, 2010) Name: GT-ID: NOTE: Partial Credit will be awarded! However, full credit will be awarded only if the correct
More informationPHYS 352 Homework 2 Solutions
PHYS 352 Homework 2 Solutions Aaron Mowitz (, 2, and 3) and Nachi Stern (4 and 5) Problem The purpose of doing a Legendre transform is to change a function of one or more variables into a function of variables
More informationLecture 5. Solid surface: Adsorption and Catalysis
Lecture 5 Solid surface: Adsorption and Catalysis Adsorbtion G = H T S DG ads should be negative (spontaneous process) DS ads is negative (reduced freedom) DH should be negative for adsorption processes
More informationToluene total oxidation over CuO-CeO 2 /Al 2 O 3 catalyst: nature and role of oxygen species
Toluene total oxidation over CuO-CeO 2 /Al 2 O 3 catalyst: nature and role of oxygen species Vladimir V. Galvita, Unmesh Menon, Konstantinos Alexopoulos, Guy B. Marin Laboratory for Chemical Technology,
More informationGAS-SURFACE INTERACTIONS
Page 1 of 16 GAS-SURFACE INTERACTIONS In modern surface science, important technological processes often involve the adsorption of molecules in gaseous form onto a surface. We can treat this adsorption
More informationDFT modeling of novel materials for hydrogen storage
DFT modeling of novel materials for hydrogen storage Tejs Vegge 1, J Voss 1,2, Q Shi 1, HS Jacobsen 1, JS Hummelshøj 1,2, AS Pedersen 1, JK Nørskov 2 1 Materials Research Department, Risø National Laboratory,
More informationCO 2 abatement by two-dimensional MXene carbides
for Journal of Materials Chemistry A. This journal is The Royal Society of Chemistry 2018 Electronic Supplementary Material (ESI) CO 2 abatement by two-dimensional MXene carbides Ángel Morales-García,
More informationThe kinetic isotope effect of Fischer-Tropsch synthesis on a cobalt stepped surface
Eindhoven University of Technology MASTER The kinetic isotope effect of Fischer-Tropsch synthesis on a cobalt stepped surface Vermeer, J.S. Award date: 2018 Link to publication Disclaimer This document
More information5. Systems in contact with a thermal bath
5. Systems in contact with a thermal bath So far, isolated systems (micro-canonical methods) 5.1 Constant number of particles:kittel&kroemer Chap. 3 Boltzmann factor Partition function (canonical methods)
More informationProbe Reactions of Alcohols and Alkanes for Understanding Catalytic Properties of Microporous Materials and Alumina Oxide Solid Acid Catalysts
Probe Reactions of Alcohols and Alkanes for Understanding Catalytic Properties of Microporous Materials and Alumina Oxide Solid Acid Catalysts A Dissertation SUBMITTED TO THE FACULTY OF UNIVERSITY OF MINNESOTA
More informationUnderstanding Brønsted-Acid Catalyzed Monomolecular Reactions of Alkanes in Zeolite Pores by Combining Insights from Experiment and Theory
Understanding Brønsted-Acid Catalyzed Monomolecular Reactions of Alkanes in Zeolite Pores by Combining Insights from Experiment and Theory Jeroen Van der Mynsbrugge, 1,3 Amber Janda, 1 Li-Chiang Lin, 4
More informationMolecular Scaffolding Strategy with Synergistic Active Centers to Facilitate Electrocatalytic CO2 Reduction to Hydrocarbon/Alcohol
Supporting Information Molecular Scaffolding Strategy with Synergistic Active Centers to Facilitate Electrocatalytic CO2 Reduction to Hydrocarbon/Alcohol Yan Jiao 1,, Yao Zheng 1,, Ping Chen 1,2,, Mietek
More informationMolecular-Level Description of Reaction Equilibrium
Molecular-Level Description of Reaction Equilibrium CBE 450: Reactor Engineering Fundamentals Fall, 2009 Prof. David Keffer dkeffer@utk.edu Thermodynamic Properties Quantum Mechanics can be combined with
More information3.5. Kinetic Approach for Isotherms
We have arrived isotherm equations which describe adsorption from the two dimensional equation of state via the Gibbs equation (with a saturation limit usually associated with monolayer coverage). The
More informationHandout 11: Ideal gas, internal energy, work and heat. Ideal gas law
Handout : Ideal gas, internal energy, work and heat Ideal gas law For a gas at pressure p, volume V and absolute temperature T, ideal gas law states that pv = nrt, where n is the number of moles and R
More informationChapter 19 Chemical Thermodynamics Entropy and free energy
Chapter 19 Chemical Thermodynamics Entropy and free energy Learning goals and key skills: Understand the meaning of spontaneous process, reversible process, irreversible process, and isothermal process.
More informationCHAPTER 9 LECTURE NOTES
CHAPTER 9 LECTURE NOTES 9.1, 9.2: Rate of a reaction For a general reaction of the type A + 3B 2Y, the rates of consumption of A and B, and the rate of formation of Y are defined as follows: Rate of consumption
More informationFIRST PUBLIC EXAMINATION SUBJECT 3: PHYSICAL CHEMISTRY
CCHE 4273 FIRST PUBLIC EXAMINATION Trinity Term 2005 Preliminary Examination in Chemistry SUBJECT 3: PHYSICAL CHEMISTRY Wednesday June 8 th 2005, 9:30am Time allowed: 2 ½ hours Candidates should answer
More informationChapter 17.3 Entropy and Spontaneity Objectives Define entropy and examine its statistical nature Predict the sign of entropy changes for phase
Chapter 17.3 Entropy and Spontaneity Objectives Define entropy and examine its statistical nature Predict the sign of entropy changes for phase changes Apply the second law of thermodynamics to chemical
More information2 Reaction kinetics in gases
2 Reaction kinetics in gases October 8, 2014 In a reaction between two species, for example a fuel and an oxidizer, bonds are broken up and new are established in the collision between the species. In
More informationRole of Heteronuclear Interactions in Selective H 2 Formation from HCOOH Decomposition on Bimetallic Pd/M (M=Late Transition FCC Metals) Catalysts
Supporting Information Role of Heteronuclear Interactions in Selective H 2 Formation from HCOOH Decomposition on Bimetallic Pd/M (M=Late Transition FCC Metals) Catalysts Jinwon Cho 1, #, Sangheon Lee 1,4,#,
More informationThe Effect of Model Internal Flexibility Upon NEMD Simulations of Viscosity
Draft: September 29, 1999 The Effect of Model Internal Flexibility Upon NEMD Simulations of Viscosity N. G. Fuller 1 and R. L. Rowley 1,2 Abstract The influence of model flexibility upon simulated viscosity
More information