Thermal Decomposition of 2-Pentanol: A Shock Tube Study and

Sep 22, 2016 - A single pulse shock tube was used to study the thermal decomposition of 2-pentanol in the temperatures between 1110 and 1325 K. Three ...
1 downloads 0 Views 1MB Size
Subscriber access provided by UNIV OF CALIFORNIA SAN DIEGO LIBRARIES

Article

Thermal Decomposition of 2-Pentanol: A Shock Tube Study and RRKM Calculations Parandaman Arathala, and Balla Rajakumar J. Phys. Chem. A, Just Accepted Manuscript • DOI: 10.1021/acs.jpca.6b06386 • Publication Date (Web): 22 Sep 2016 Downloaded from http://pubs.acs.org on September 23, 2016

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry A is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 42

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Thermal Decomposition of 2-Pentanol: A Shock Tube Study and RRKM Calculations A. Parandaman and B. Rajakumar* Department of chemistry, Indian Institute of Technology Madras, Chennai-600036, India. *Address for correspondence: [email protected] Telephone number for the corresponding author(s) +91-44-2257-4231

Abstract Single pulse shock tube was used to study the thermal decomposition of 2-pentanol in the temperatures between 1110 and 1325 K. Three major decomposition products are methane, ethylene and propylene. The minor products detected in lower concentration are ethane, acetylene, acetaldehyde, 1-pentene and 2-pentene. The rate coefficient for the overall  decomposition of 2-pentanol was found to be  (1110-1325 K) = (4.01 ± 0.51) × 109 exp (-

(36.2 ± 4.7)/RT) s-1, where the activation energies are given in kcal mol-1. To simulate reactant and product distribution over the experimentally studied temperatures between 1110 and 1325 K, a reaction scheme was constructed with 34 species and 39 reactions. In addition to this, the temperature and pressure dependent rate coefficients were computed for various unimolecular dissociation pathways using RRKM theory. The high pressure limit rate coefficient for overall  

decomposition of 2-pentanol was obtained to be  (500-2500 K) = (9.67 ± 1.11) × 1014 exp (-(67.7 ± 2.9)/RT) s-1. The calculated high pressure rate coefficients and experimentally measured rate constants are in good agreement with each other. The reaction is majorly governed by the unimolecular elimination of water.

1 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

1. Introduction The impact of fossil fuels on environment, increase in energy requirements and fuel cost have motivated serious search for alternative fuels in recent times. The suitable alternative fuels such as biofuels have been recommended in place of fossil fuels. Bioalcohol, biodiesel and biogas are known as biofuels, which are derived from biomass. It can be used in a pure form or in blend with fossil fuels.1-4 The interest in biofuels leads to consideration of alcohols as a fuel. The first generation biofuels, methanol and ethanol are extensively used in transportation sector all over the world. Approximately 90% of ethanol is used to produce biofuel all over the world,5,6 but ethanol have various drawbacks such as low energy density, high hygroscopicity and corrosiveness towards pipelines and engines.7,8 The next generation biofuels, butanol and pentanol were considered as potential biofuels,9-11 they exhibit better combustion properties including lower hygroscopicity and corrosivity, greater energy density and can be readily mixed with hydrocarbon fuels in fueling systems.12,13 The main disadvantage of biofuels is the formation of carbonyl compounds in the decomposition, which will further act as pollutants.6,14 Hence, the detailed mechanisms of the decomposition of these compounds are required in various conditions to get the best possible understanding of the combustion. The combustion of methanol and ethanol was studied by using various experimental techniques.15 The next generation alcohol is butanol, which is widely studied by combustion community over the last half a decade. Many experimental results on pyrolysis and oxidation of butanols were reported under various experimental conditions.9,16-18 Various computational methods were used to determine the rate coefficients and to understand various reaction 2 ACS Paragon Plus Environment

Page 2 of 42

Page 3 of 42

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

pathways involved in combustion of butanol.19,20 As a consequence, several research groups have developed chemical kinetic models for combustion of butanol in different conditions.16,18,21 Studies on the combustion of pentanol are very limited. The oxidation mechanism of isopentanol was studied at lower temperatures between 550 and 750 K by Weltz et al.22 using time-resolved tunable synchrotron photo ionization mass spectrometry. Togbe et al.6 studied decomposition of n-pentanol and isopentanol using jet stirred reactor (JSR) and reported the complete kinetic models for n-pentanol and isopentanol. Zhao et al.23 studied thermal decomposition of 1-pentanol and its isomers using computational methods. There are no experimental and computational studies on thermal decomposition of 2pentanol are reported so far. Thermal decomposition of 2-pentanol, was carried out using shock tube for the first time in the temperatures between 1110 and 1325 K. In this study, the most suitable mechanism for the decomposition was proposed by following the complete degradation of reactant and formation of products. To understand the experimental results further, unimolecular decomposition reactions of 2-pentanol and their rate coefficients were explored with computational methods coupled with RRKM theory. The calculated and measured rate coefficients were compared with each other. The possible mechanism for overall decomposition was proposed and simulated. This present study reports complete decomposition mechanism of 2-pentanol in the studied temperature range.

2. Experimental Section Single Pulse Shock Tube (SPST) was used to carry out experiments in the present study. The detailed description of the SPST set-up used in this work were discussed elsewhere.24,25 The samples containing 6 Torr of 2-pentanol and 60 Torr of diluted (1.7% in argon) 1,1,1-

3 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

trifluoroethane (CH3CF3) were added and diluted with argon up to the desired pressure (typically in the range of 200-600 Torr) depending on the required temperature in the test section. The temperatures behind the reflected shock waves were measured using chemical thermometric method.26,27 CH3CF3 was used as an internal standard to determine reaction temperature. It decomposes into CH2=CF2 + HF via a first order unimolecular reaction, whose rate constant of k = 8.14 × 1013exp (-71154/RT) s-1. Where the activation energy is given in cal mol-1 was reported earlier.28 The detailed description of measuring the reaction temperature using CH3CF3 as internal standard was given in our previous work.25 The experimental procedure was described in our recently reported articles24,25 and similar procedure was followed in the current experimentation. The post shock samples were withdrawn from the end of the driven section of the shock tube and were analyzed using gas chromatograph (Agilent Technologies 6890N) with a Flame Ionization Detector (FID). A 2 m long packed Hayesep-A column was used for the analysis and the oven temperature was programmed between 75 and 1500C. The method described in our previous work25 was followed here in measuring the concentrations of all the components of the post shocked mixture. The post shock mixture was also analyzed qualitatively using FTIR spectrometer (Bruker's VERTEX 70) as described earlier.24,25

3. Materials and chemicals The purities of the chemicals were as follows: Helium, HP certified to >99.995%, (Praxair); Argon (Ar), certified to >99.995%; 2-Pentanol, >99% (Sigma Aldrich); 1-pentene, >99.5% (Sigma Aldrich), 2-pentene, >99%, (Sigma Aldrich); acetaldehyde 99% (Sigma Aldrich), CH3CF3 (99%) and CH2CF2(99%) were purchased from SynQuest Laboratories are 4 ACS Paragon Plus Environment

Page 4 of 42

Page 5 of 42

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

used in these experiments. Methane (99.5%), propylene (99.5%), ethane (99.5%) and ethylene (99.5%) are purchased from Praxair and were used without further purification in the present work.

4. Computational 4.1. Quantum chemical calculations Geometries of reactant (2-pentanol), transition states (TSs) and products of possible unimolecular dissociation channels of 2-pentanol were optimized at M06-2X level of theory with basis set29,30 6-31+G(d,p) which are internally available in Gaussian 09 program suite.31 The basis set 6-31+G(d,p) includes polarization functions on over all atoms and diffuse functions on non-hydrogen atoms. Various previous studies25,32,33 used 6-31+G(d,p) basis set for predicting reliable geometries and vibrational frequency calculations of the stationary points for the similar molecules. Therefore, 6-31+G(d,p) was used as a basis set in the present calculations. The optimized reactant and products were identified with zero imaginary frequency and transition states identified with one imaginary frequency. The Intrinsic Reaction Coordinate calculations34 was carried out to obtain the minimum energy path (MEP) at M06-2X functional with 631+G(d,p) basis set to confirm that the designated transition states connect the reactant and products. The normal modes of all the stationary points are viewed in Gauss View.35 The pressure and temperature dependentrate coefficients were calculated for unimolecular dissociation reactions of 2-pentanol with RRKM theory using KiSThelP program.36 Lennard-Jones (L-J) potential was used to model the interaction between reactant and bath gas Ar. The Lennard Jones parameters σ and ε/k  for 2-pentanol and Ar were calculated using the following equations37

5 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60



σ = 2.44T ⁄P   ⁄ = 0.77!" Where TC and PC are critical temperature and critical pressure respectively taken from NIST database.38 kB is the Boltzmann constant. The L-J parameters for 2-pentanol σ = 6.046 Å and ε/kb = 431.2 K were used in present study. The L-J parameters for Ar bath gas are σ = 3.47 Å and ε/kb = 113.5 K, which were taken from KiSThelP program.36

5. Results and discussion 5.1 Experimental Section In the present study, a total of 38 experiments were performed to determine the distribution of reaction products over the temperature range of 1110-1325 K. The normalized concentrations of reactant, all other products of each experiment and the experimental conditions are given in Table 1. The decomposition products of 2-pentanol and CH3CF3 were identified to be methane, ethane, ethylene, acetylene, 1,1-difluoroethylene, propylene, acetaldehyde, 1pentene and 2-pentene. 1,1-difluoroethylene (CH2CF2) is the decomposition product of CH3CF3, which was used as internal standard. Remaining products were formed solely from the decomposition of 2-pentanol. To confirm this, only 2-pentanol diluted with argon was decomposed using shock tube in the studied temperature range and all these products were identified. Typical chromatogram of the post shocked mixture obtained where the test mixture was heated to 1276 K is shown in Figure 1. The FTIR spectrum of the post shock mixture of 2-pentanol decomposed at 1276 K is shown in Figure 2. The IR spectrum analysis indicates the presence of 2-pentanol, methane, ethylene, acetylene and CO. The medium intense bands in the region 3100-3400 cm-1 and at 975 6 ACS Paragon Plus Environment

Page 6 of 42

Page 7 of 42

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

cm-1 are assigned to the existence of 2-pentanol. The medium intensity bands in the region 20502220 cm-1 represents the presence of CO. The presence of acetylene was indicated by the bands at 3250-3300 cm-1 and high intense band at 750 cm-1. The distribution of the normalized concentrations of reactant and all other products in the studied temperature range are shown in Figure 3. Concentrations of products methane, ethylene, ethane, acetylene and acetaldehyde are found to be increasing with the increase in temperature. The concentrations of propylene was found to be increasing up to 1250 K and then started decreasing. The other products 1-pentene and 2-pentene are found to be increasing up to 1200 K and then decreasing. The first order rate coefficients for the total decomposition of 2-pentanol was determined by using the following equation, 1 (2 − )*'%+',&-  = − &' % (2 − )*'%+',&-. Where [2-pentanol]t and [2-pentanol]0 are experimentally quantified concentrations of 2pentanol at the end of the reaction time ‘t’ and initial concentration respectively. The Arrhenius plot of the experimentally measured rate coefficients is shown in Figure 4. The rate coefficients data were fit by linear least square method. The obtained Arrhenius expression for the overall  decomposition of 2-pentanol in the studied temperature range is  (1110-1325 K) = (4.01 ±

0.51) × 109exp (-(36.2 ± 4.7)/RT) s-1. The rate coefficients obtained using RRKM theory, are also appended to Figure 4 for comparison. The experimentally measured and computed rate coefficients shows very good agreement with each other in the studied temperatures between 1110 and 1325 K, vide infra. The detailed error analysis was performed to quantify the total error in the experimentally measured rate coefficients in the studied temperature range. The procedure for the estimation of the errors was reported in our recent work.25 The total uncertainties from all sources in the 7 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

experimental measurements were calculated using a root-sum-squared method. The estimated overall (2σ) uncertainty in the present experiments is ± 21%.

5.2 Computational 5.2.1 Electronic structures and energetics The geometries of all the stationary points optimized at M06-2X/6-31+G(d,p) level of theory are shown in Figure 5. The energy level diagram consists of reactant, transition states and products are labeled for the reactions R1-R7 is given in Figure S-I of the Supporting Information. The energy barriers for unimolecular decomposition of 2-pentanol estimated at same level of theory are given in Table 2. The structural parameters and vibrational frequencies of all the stationary points involved in the reactions R1-R7 obtained at the M06-2X/6-31+G(d,p) level of theory are given in the Supporting Information, SI (Tables S-I to S-II). Figure S-II in the Supporting Information reports R1-R7 reaction channels with corresponding structures identified for the unimolecular thermal decomposition of 2-pentanol. 2-pentanol undergoes seven unimolecular pathways which are labeled as R1-R7. These pathways yield stable reaction products through transition states. The molecular elimination channel R1 goes through four membered ring transition state (TS1) to form 1-pentene (P1) and H2O (P2) with barrier energy of 64.6 kcal mol-1. Similarly, reaction channel R2 goes through four membered ring transition state (TS2) to form 2-pentene (P3) and H2O (P4) but with the barrier of 66.9 kcal mol-1. Reaction channel R3 yields ethane (P5) and acetone (P6) by the transfer of –OH group hydrogen to methylene (-CH2-) group followed by C-C bond dissociation through transition state (TS3) with the energy barrier of 101 kcal mol-1. Reaction channel R4 forms ethylene (P7) and 2-hydroxy propane (P8) through transition state TS4 by intramolecular 8 ACS Paragon Plus Environment

Page 8 of 42

Page 9 of 42

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

transfer of hydrogen atom from –CH3 group to –CH2 group followed by CH2-CH2 bond dissociation with a barrier of 113 kcal mol-1. Reaction R5 forms butyraldehyde (P9) and methane (P10) through transition state TS5 by migration of hydrogen atom from –OH group to –CH3 group followed by CH(OH)-CH3 bond dissociation with barrier of 84 kcal mol-1. The reaction R6 produces methane (P11), ethylene (P12) and acetaldehyde (P13) by transfer of H atom from –OH group to methyl group followed by CH3-CH2 bond dissociation with a barrier of 94 kcal mol-1. The reaction R7 produces propane (P14) and acetaldehyde (P15) via migration of hydrogen atom from –OH group to –CH2 group followed by CH2-CH(OH) bond dissociation with a barrier of 85 kcal mol-1. It is clear that the reaction pathways R3 and R4 are very high barrier reactions and are less possible in the current experimental conditions. However, they are included in our computations for the completeness. It is obvious from these seven channels that, the reaction channel R1 is energetically more favorable when compared with all the other possible reaction channels.

5.2.2. Kinetic analysis The pressure and temperature dependent rate coefficients for the seven unimolecular reaction channels (R1-R7) were calculated using RRKM theory in the temperatures between 500 and 2500 K with KiSThelP program. Among all the channels, formation of 1-pentene (R1) reaction pathway is the major one. The rate coefficients were calculated by varying the pressure from low to high pressures i.e. 0.01 atm, 0.1 atm, 1 atm, 10 atm and 100 atm for all the seven pathways. The computed high pressure limit rate coefficients for all the reaction channels R1-R7 are plotted over the temperature between 500 and 2500 K are shown in Figure 6. The pressure and temperature dependent Arrhenius plot over the temperatures between 1300 and 2500 K and

9 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 42

pressure range of 0.01-100 atm for the major channel R1 is shown in Figure 7, as this channel is the major contributor. As seen from Figure 7, the rate coefficient is pressure independent at temperatures below 1400 K, but at temperatures above 1400 K the rate coefficient is found to be increasing with increase in pressure. The individual rate coefficients of the unimolecular reaction channels were added to obtain the overall rate coefficient for the title reaction. The comparison of reported rate coefficients with the current experimentally measured rate coefficients in the studied temperature range is shown in Figure 4. The linear least squares method was used to fit the computed rate coefficients between 500 and 2500 K. The obtained temperature dependent Arrhenius expression  

is  (500-2500 K) = (9.67 ± 1.11) × 1014exp [-(67.7 ± 2.9)/RT] s-1, where the activation energies are given in kcal mol-1. The rate coefficients calculated using RRKM theory and experimentally measured rate coefficients in the temperature between 1110 and 1325 K are in very good agreement with each other.

6. Kinetic simulations The reaction scheme given in Table 3 was proposed for understanding the reaction mechanism and also to model the distribution of the concentrations of reactant and products. The proposed reaction scheme contains 34 species and 39 elementary reactions. The rate coefficients are expressed as k = A exp(-Ea/RT) or k = A Tnexp(-Ea/RT) and are given in Table 3. The experimentally determined temperatures and reaction times were used in kinetic simulations. The plot of the concentrations obtained in the simulation using proposed reaction mechanism and experimentally measured concentrations are shown in Figure 3a-i. The simulated concentrations

10 ACS Paragon Plus Environment

Page 11 of 42

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

of all the products and the reactants are in good agreement with experimentally measured concentrations. 2-pentanol decomposes via three possible channels, which are unimolecular. These reactions include water elimination, intramolecular hydrogen transfer and simple bond scission reactions (C-C, C-O, C-H and O-H). The rate coefficients for reactions R1-R7 were taken from the RRKM calculations and all other reactions rate coefficients were taken from past research works from various groups.39-52 The H2O elimination always has the lowest energy barrier than the other small molecular elimination reactions in thermal decomposition of 2-pentanol. The energy barrier for water elimination reactions is 10 kcal mol-1 lesser than the other molecular elimination reaction (methane, ethylene etc.). CH3-CH2-CH2-CH(OH)-CH3 → CH3-CH2-CH2-CH=CH2 + H2O k1 = 1.58 × 1014exp (-67.1/RT) s-1

(R1)

CH3-CH2-CH2-CH(OH)-CH3 → CH3-CH2-CH=CH-CH3 + H2O k2 = 3.01 × 1014exp (-69.5/RT) s-1

(R2)

In the decomposition of 2-pentanol, four homolytic C-C bond scissions and one C-O bond scission are quite possible, particularly in the studied temperature range. The C-C bond and C-O bond energies are much lower than those of C-H and O-H bonds.18,53,54 Therefore, the former type of bond scission processes would not compete with the latter ones and hence, they were not included in the proposed reaction scheme. CH3-CH2-CH2-CH(OH)-CH3 → CH3 + CH2-CH2-CH(OH)-CH3 k8 = 5.12 × 1026(T)-2.92exp (-90.5/RT) s-1 11 ACS Paragon Plus Environment

(R8)

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 42

CH3-CH2-CH2-CH(OH)-CH3 → CH3-CH2 + CH2-CH(OH)-CH3 k9 = 8.54 × 1027 (T)-3.23exp (-91.7/RT) s-1

(R9)

CH3-CH2-CH2-CH(OH)-CH3 → CH3-CH2-CH2 + CH(OH)-CH3 k10 = 3.01 × 1025 (T)-2.57exp (-89.9/RT) s-1

(R10)

CH3-CH2-CH2-CH(OH)-CH3 → CH3-CH2-CH2-CH-CH3 + OH k11 = 2.82 × 1029 (T)-3.60exp (-101.9/RT) s-1

(R11)

CH3-CH2-CH2-CH(OH)-CH3 → CH3-CH2-CH2-CH(OH) + CH3 k12 = 5.13 × 1026 (T)-2.91exp (-90.5/RT) s-1

(R12)

6.1 Methane In the studied experimental temperature range, the maximum concentration of methane was found to beabout 24%. The methyl radicals formed in the decomposition of 2-pentanol play a major role in the formation of methane through reactions R8, R12 and R24. The formed CH3 radicals react with H atom (R16), abstract H atom from ethylene (R19) and water (R26) and form methane. CH3 + H → CH4

k16 = 6.93 × 1013(T)0.18cm3mol-1s-1

(R16)

C2H4 + CH3 → C2H3 + CH4

k19 = 4.16 × 1012exp (-11.1/RT) cm3mol-1s-1

(R19)

CH3 + H2O → CH4 + OH

k26 = 0.64 × 101(T)3.31exp (-12.6/RT) cm3mol-1s-1

(R26)

The reactions R5 and R6 are added to the reaction scheme for the direct formation of methane from 2-pentanol, whose rate coefficients were calculated using RRKM theory. In reaction R6, 12 ACS Paragon Plus Environment

Page 13 of 42

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

The transfer of H atom from –OH group to –CH3 group followed by CH3-CH2 bond dissociation by forming the products CH4 and CH3-CH(O)-CH2-CH2. The formed unstable biradical product CH3-CH(O)-CH2-CH2 will immediately undergo decomposition at CH(O)-CH bond by forming stable products acetaldehyde and ethylene. CH3-CH2-CH2-CH(OH)-CH3 → CH4 + CH3-CH2-CH2-CHO k5 = 1.92 × 1014exp (-86.2/RT) s-1

(R5)

CH3-CH2-CH2-CH(OH)-CH3 → CH4 + C2H4 + CH3-CHO k6 = 1.60 × 1014exp (-96.1/RT) s-1

(R6)

Methane is also formed from acetaldehyde via direct molecular elimination of CO (R17). The contribution of this reaction is very less, as maximum concentration of acetaldehyde formed in the reaction is 7%. In addition, the activation energy for this reaction is very high (85 kcal mol-1). CH3-C(O)-H → CH4 + CO

k17 = 1.00 × 1015exp (-85.0/RT) s-1

(R17)

On the other hand, the formed methane react with H atom (R37) and methyl radicals (R38) to regenerate methyl and ethyl radicals respectively. CH4 + H → H2 + CH3

k37 = 1.77 × 1014 exp (-13.8/RT) cm3mol-1s-1

(R37)

CH4 + CH3 → H2 + CH3-CH2

k38 = 1.00 × 1013exp (-23.1/RT) cm3mol-1s-1

(R38)

The rate coefficients of all the elementary reactions were varied by a factor of 10 in order to analyze the sensitivity towards the formation and decomposition of methane. The sensitivity of each reaction in terms of percentage change in the concentrations with respect to the reaction numbers are given in Table 4. The reactions R8 (C4-C5 bond scission in the reactant), R12 (C4C5 bond scission in the reactant) and R24 (C4-C5 bond scission in 1-pentene, which is formed 13 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 42

via 1,2- unimolecular elimination of H2O from 2-pentanol through the reaction R1) are very sensitive to the formation of CH4. Therefore, methyl radicals plays major role for the formation of methane. Reactions R1 and R2 (2,3- unimolecular elimination of H2O from 2-pentanol) also shows major impact on the formation of methane.

6.2 Ethylene The decomposition of 2-pentanol leads to the formation of ethylene, which is one of the major product. Maximum concentrations of C2H4 were found to be 39% which is formed via reactions R4 and R6. The rate coefficients for these reactions were taken from the present RRKM calculations. R6 is a low barrier (96.1 kcal mol-1) reaction than that of R4 (115.9 kcal mol-1) and hence it contributes more to the formation of ethylene. CH3-CH2-CH2-CH(OH)-CH3 → C2H4 + CH3-CH(OH)-CH3 k4 = 6.01 × 1014exp (-115.9/RT) s-1

(R4)

CH3-CH2-CH2-CH(OH)-CH3 → CH4 + C2H4 + CH3-CHO k6 = 1.60 × 1014exp (-96.1/RT) s-1

(R6)

Initially, it was assumed that ethylene could be the product from reactions of methyl and other radicals formed in the decomposition of the reactant via following reactions CH3-CH2 → C2H4 + H

k14 = 3.06 × 1010(T)0.95exp (-36.9/RT) s-1

(R14)

2 CH3 → C2H4 + H2

k18 = 9.90 × 1015exp (-32.9/RT) cm3mol-1s-1

(R18)

C2H3 + H → C2H4

k20 = 3.88 × 1013(T)0.20cm3mol-1s-1

(R20)

14 ACS Paragon Plus Environment

Page 15 of 42

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

C2H2 + H2 → C2H4

k39 = 3.01 × 1011exp(-38.9/RT) cm3mol-1s-1

(R39)

All the above reactions involve CH3, C2H5, C2H3 radicals. Among all the above reactions involving CH3, C2H5, C2H3 radicals, reaction R14 and R18 are expected to contribute more as the concentrations of methyl and ethyl radicals are higher when compared to that of C2H3 radicals. The sensitivity analysis of the proposed reaction scheme is carried out for the formation of C2H4 and the results are given in Table 4. Reactions R1(1,2- unimolecular elimination of H2O from 2pentanol to form 1-pentene which further decompose via C-C bond scission to form C2H5 radicals in R25), R9 (C3-C4 bond scission in the reactant), R14 (decomposition of ethyl radical to form ethylene), R23 (ethane reacts with H atom to form C2H5 radicals) and R25 show significant sensitivity towards the formation of C2H4.

6.3 Propylene It was observed that the propylene concentration was increasing up to 1250 K and further increase in the temperature leads to decrement in concentration of propylene. The main possible reaction for the formation of propylene is the decomposition of propyl radical through reaction R15. The other minor reaction channels for the formation of propylene are R30 and R33. CH3-CH2-CH2 → C3H6 + H

k15 = 4.14 × 1012(T)0.17exp(-35.6/RT) s-1

CH3-CH2-CH2-CH-CH3→ C3H6 + CH3-CH2 k30 = 1.22 × 1012(T)0.64exp(-29.4/RT) s-1 CH2-CH(OH)-CH3 → CH3-CH=CH2 + OH

k33 = 3.70 × 1012exp(-26.8/RT) s-1

15 ACS Paragon Plus Environment

(R15) (R30) (R33)

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 42

The concentration of propylene decreases above 1250 K due to its own dissociation (R34) to form CH3 and C2H3 radicals and reaction with H atom to form n-C3H7 and iso-C3H7 radicals through reactions R35 and R36 respectively. C3H6 → CH3 + C2H3

k34 = 8.00 × 1014exp(-88.0/RT) s-1

(R34)

C3H6 + H → n-C3H7

k35 = 4.57 × 1012(T)0.51exp(-2.6/RT) cm3mol-1s-1

(R35)

C3H6 + H → iso-C3H7

k36 = 7.75 × 1012(T)0.51exp(-1.2/RT) cm3mol-1s-1

(R36)

The sensitivity analysis in Table 4 shows that reaction R10 (C2-C3 bond scission in the reactant to form C3H7 radicals) will have significant sensitivity towards the formation of propylene.

6.4 Minor products The detectedminor reaction products in the decomposition of 2-pentanol were ethane, acetylene, acetaldehyde, 1-pentene and 2-pentene. On comparing it with main reaction products, the concentrations of these individual minor products are found to be less than 10%. Another insignificant product that was identified was CO, which was qualitatively analyzed by FTIR spectroscopy. The concentrations of products ethane, acetylene and acetaldehyde were observed to be increasing in the studied temperature range. The concentrations of the other two products 1pentene and 2-pentene seem to be increasing up to 1200 K and beyond this the concentrations were found to be decreasing with the increase in temperature. The major channel for the production of C2H6 is the molecular elimination of ethane through reaction R3 and reaction between ethyl radical with H atom and water molecule through reactions R13 and R15 respectively. CH3-CH2-CH2-CH(OH)-CH3 → C2H6 + CH3-C(O)-CH3 16 ACS Paragon Plus Environment

Page 17 of 42

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

k3 = 1.90 × 1014exp (-103.4/RT) s-1

(R3)

CH3-CH2 + H → C2H6

k13 = 5.44 × 1013(T)0.16cm3mol-1s-1

(R13)

C2H5 + H2O → C2H6 + OH

k52 = 3.39 × 106(T)1.44exp (-20.3/RT) cm3mol-1s-1

(R52)

The reaction channels through which acetylene can be formed using the present scheme are R22 and R31. The maximum concentration of acetylene in present experimental study was measured to be 6%. C2H3 + H → C2H2 + H2

k22 = 1.21 × 1013cm3mol-1s-1

(R22)

C3H5 → C2H2 + CH3

k31 = 1.26 × 1013exp (-33.7/RT) s-1

(R31)

Acetaldehyde concentration in the present experiments was measured to be around 7%. It was formed through reactions R7 and R21. CH3-CH2-CH2-CH(OH)-CH3 → C3H8 + CH3-CHO k7 = 9.44 × 1013exp (-86.4/RT) s-1 (R7) CH(OH)-CH3 → CH3-CHO + H

k21 = 5.00 × 1013exp (-21.9/RT) s-1

(R21)

1-pentene and 2-pentenewere formed through unimolecular elimination of water from 2pentanol through reactions R1 and R2. These two reaction products were observed to get decomposed beyond 1200 K and the reactions were added in the reaction scheme. The rate coefficients for C-C bond scission reactions of 1-pentene and 2-pentene are not available in the literature. Therefore, the rate coefficients for the C-C bond scissions were taken from the analogous systems such as 1-butene and 2-butene.47 CH3-CH2-CH2-CH=CH2 → CH3+ CH2-CH2-CH=CH2

17 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

k24 = 1.10 × 1016exp (-72.9/RT) s-1

Page 18 of 42

(R24)

CH3-CH2-CH2-CH=CH2 → CH3-CH2+ CH2-CH=CH2 k25 = 1.10 × 1016exp (-70.7/RT) s-1

(R25)

CH3-CH2-CH=CH-CH3 → CH3 + CH2-CH=CH-CH3 k28 = 3.16 × 1017exp (-99.4/RT) s-1

(R28)

CH3-CH2-CH=CH-CH3 → CH3 + CH3-CH2-CH=CH k27 = 3.16 × 1017exp (-99.4/RT) s-1

(R27)

CH3-CH2-CH=CH-CH3 → CH3-CH2 + CH3-CH=CH k29 = 1.00 × 1016exp (-75.7/RT) s-1

(R29)

Sensitivity analysis was carried out for the formation and decomposition of minor products and the results are given in Table 4. The reactions R1, R2, R9 and R32 shows more sensitivity towards the formation of ethane. Reaction R29 has shown more sensitivity towards the formation of C2H2. Reaction R10 is more sensitive for the formation of acetaldehyde. Reactions R1 and R2 have shown significant sensitivity towards the formation of 1-pentene and 2-pentene.

7. Conclusions Thermal decomposition of 2-pentanol was studied in the temperatures between 1110 and 1325 K. The main reaction products identified in the experiments were methane, ethylene and propylene; and the minor reaction products detected in lower concentrations were ethane, acetylene, acetaldehyde, 1-pentene and 2-pentene. The rate coefficient for the overall

18 ACS Paragon Plus Environment

Page 19 of 42

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

 decomposition of 2-pentanol was found to be  (1110-1325 K) = (4.01 ± 0.51) × 109exp (-

(36.2 ± 4.7)/RT) s-1. The major unimolecular decomposition pathways of 2-pentanol were investigated using computational methods. The RRKM theory was used to calculate the rate coefficients for unimolecular decomposition reactions. Among the unimolecular decomposition pathways, elimination of H2O was found to be the major channel giving 1-pentene. The obtained  

temperature dependent Arrhenius expression in studied temperature range is  (500-2500 K) = (9.67 ± 1.11) × 1014exp (-(67.7 kcal mol-1 ± 2.9)/RT) s-1. A reaction scheme for the decomposition of 2-pentanol was constructed with 34 species and 39 reactions. The simulated concentrations of all the products and the reactants shows good agreement with experimentally measured concentrations using proposed reaction scheme and the results can be used for further improvement in the accuracy and completeness of the 2-pentanol kinetic mechanism. The formation of methane, ethane and other lower hydrocarbons could be explained only when C-C bond scission was taken into consideration in the kinetic mechanism. The decomposition of the test molecule by C-C bond scission at various positions was competing with H2O elimination. The formation of carbon monoxide in the decomposition of 2-pentanol should be taken into consideration, before its use as a biofuel additive.

Supporting Information Table S-I. Optimized geometries of the reactant, transition states and products at M06-2X level of theory with 6-31+G( d, p) basis set. Table S-II. Vibrational frequencies (cm-1) of the reactant, transition states and products at M062X level of theory with 6-31+G (d, p) basis set.

19 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 42

Figure S-I. Potential energy surfaces of reaction pathways R1-R7 computed using the M062X/6-31+G(d,p) level of theory. The energies are given in the units of kcal mol-1. Figure S-II. The considered reaction pathways (R1-R7) for the unimolecular decomposition of 2-pentanol (The rate coefficients of these reactions were calculated using RRKM theory were used in the kinetic simulations).

Acknowledgment B.R. thanks Defence Research and Development Organization (DRDO), Government of India for funding. A.P. is very grateful to CSIR for providing a research fellowship. A.P. also thank Mr. G. Sudhakar and Mr. G. Srinivasulu for their help in doing the experiments.

References (1) Azad, A. K.; Rasul, M. G.; Khan, M. M. K.; Sharma, S. C.; Hazrat, M. A. Prospect of biofuels as an alternative transport fuel in Australia. Renew Sust Energy Rev 2015, 43, 331-351. (2) Varuvel, E. G.; Mrad, N.; Tazerout, M.; Aloui, F. Experimental analysis of biofuel as an alternative fuel for diesel engines. Appl Energy 2012, 94, 224-231. (3) Kohse-Hoinghaus, K.; Osswald, P.; Cool, T. A.; Kasper, T.; Hansen, N.; Qi, F.; Westbrook, C. K.; Westmoreland, P. R. Biofuel combustion chemistry: From ethanol to biodiesel. Angew. Chem.2010, 49, 3572-3597. (4) Christian, W.; Alexander, F.; Feline, B.; Dawid, B.; Heiko, D.; Thorsten, S.; Eckhard, B.; Trends and challenges in the microbial production of lingo cellulosic bioalcohol fuels. Appl. Microbiol. Biotechnol. 2010, 87, 1303-1315. (5) Tang, C.; Wei, L.; Man, X.; Zhang, J.; Huang, Z. High temperature ignition delay times of C5 primary alcohols. Combust. Flame 2012, 160, 520-529. (6) Togbe, C.; Halter, F.; Foucher, F.; Mounaim-Rousselle, C.; Dagaut, P. Experimental and detailed kinetic modeling study of 1-pentanol oxidation in a JSR and combustion in a bomb. Proc. Combust. Inst. 2011, 33, 367-374.

20 ACS Paragon Plus Environment

Page 21 of 42

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(7) Black, G.; Curran, H. J.; Pichon, S.; Simmie, J. M.; Zhukov, V. Bio-butanol: combustion properties and detailed chemical kinetic model. Combust. Flame 2010, 157, 363-373. (8) Tsujimura, T.; Pitz, W. J.; Gillespie, F.; Curran, H. J.; Weber, B. W.; Zhang Y.; Sung C. J. Development of isopentanol reaction mechanism reproducing autoignition character at high and low temperatures. Energy Fuels 2012, 26, 4871-4886. (9) Rosado-Reyes, C. M.; Tsang, W. Shock tube studies on the decomposition of 2-butanol. J. Phys. Chem. A 2012, 116, 9599-9606. (10) Li, L.; Wang, J.; Wang, Z.; Liu, H. Combustion and emissions of compression ignition in a direct injection diesel engine fueled with pentanol. Energy 2015, 80, 575-581. (11) Nielsen, D. R.; Amarasiriwardena, G. S.; Prather, K. L. J. Predicting the adsorption of second generation biofuels by polymeric resins with applications for in situ product recovery (ISPR). Bioresour. Technol. 2010, 101, 2762-2769. (12) Johnson, M. V.; Goldsborough, S. S.; Serinyel, Z.; O’Toole, P.; Larkin, E.; O’Malley, G.; Curran, H. J. A shock tube study of n- and iso-propanol ignition. Energy Fuels 2009, 23, 58865898. (13) Demirbas, A. Progress and recent trends in biofuels Prog. Energy Combust. Sci. 2007, 33, 1-18. (14) Yang, M.; Karra, S. B.; Senkan, S. M. Equilibrium analysis of combustion/incineration. Hazard. Waste. Hazard. Mater., 1988, 4, 55-68. (15) Cooke, D. F.; Dodson, M. G.; Williams, A. A shock-tube study of the ignition of methanol and ethanol with oxygen. Combust. Flame 1971, 16, 233-236. (16) Hansen, N.; Harper, M. R.; Green, W. H. High-temperature oxidation chemistry of nbutanol-experiments in low-pressure premixed flames and detailed kinetic modeling. Phys. Chem. Chem. Phys. 2011, 13, 20262-20274. (17) Veloo, P. S.; Egolfopoulos, F. N. Flame propagation of butanol isomers/air mixtures. Proc. Combust. Inst. 2010, 33, 987-993. (18) Sarathy, S. M.; Thomson, M. J.; Togbe, C.; Dagaut, P.; Halter, F.; Mounaim-Rousselle, C. An experimental and kinetic modeling study of n-butanol combustion. Combust. Flame 2009, 156, 852-864.

21 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 42

(19) El-Nahas, A. M.; Mangood, A. H.; Takeuchi, H.; Taketsugu, T. Thermal decomposition of 2-butanol as a potential non fossil fuel: A computational study. J. Phys. Chem. A 2011, 115, 2837-2846. (20) Zhou, C. W.; Simmie, J. M.; Curran, H. J. Rate constants for hydrogen-abstraction by image from n-butanol. Combust. Flame 2011, 158, 726-731. (21) Black, G.; Curran, H. J.; Pichon, S.; Simmie, J. M.; Zhukov, V. Bio-butanol: combustion properties and detailed chemical kinetic model. Combust. Flame 2010, 157, 363-373. (22) Welz, O.; Zádor, J.; Savee, J. D.; Ng, M. Y.; Meloni, G.; Fernandes, R. X.; Sheps, L.; Simmons, B. A.; Lee, T. S.; Osborn, D. L.; Taatjes, C. A. Low-temperature combustion chemistry of biofuels: pathways in the initial low-temperature (550 K–750 K) oxidation chemistry of isopentanol. Phys. Chem. Chem. Phys. 2012, 14, 3112-3127. (23) Zhao, L.; Ye, L.; Zhang, F.; Zhang, L.Thermal decomposition of 1‑pentanol and Its isomers: A theoretical study. J. Phys. Chem. A 2012, 116, 9238-9244. (24) Sudhakar, G.; Rajakumar, B. Thermal decomposition of 1-chloropropane behind the reflected shock waves in the temperature range of 1015-1220 K: Single pulse shock tube and computational studies. J. Chem. Sci. 2014, 126, 897-909. (25) Parandaman, A.; Balaganesh, M.; Rajakumar, B. Experimental and theoretical study on thermal decomposition of methyl butanoate behind reflected shock waves. RSC Adv., 2015, 5, 86536–86550. (26) Lifshitz, A.; Tamburu, C.; Suslensky, A. Decomposition and isomerization of 1,2benzisoxazole: Single-pulse shock-tube experiments, quantum chemical and transition-state theory calculations. J. Phys. Chem. A 2006, 110, 11677-11683. (27) Manion, J. A.; Awan, I. A. Evaluated kinetics of terminal and non-terminal addition of hydrogen atoms to 1-alkene: a shock tube study of H + 1-butene. J. Phys. Chem. A, 2015, 119, 429-441. (28) Akira, M.; Kenji, Y.; Hiroumi, S. Thermal decomposition of 1,1,1-Trifluoroethane revisited. J. Phys. Chem. A 2014, 118, 11688-11695. (29) Zhao, Y.; Truhlar, D. G. Density functionals with broad applicability in chemistry. Acc. Chem. Res., 2008, 41, 157-161. (30) Zhao, Y.; Truhlar, D. G. The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and 22 ACS Paragon Plus Environment

Page 23 of 42

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

transition elements: two new functionals and systematic testing of four M06-class functionals and 12 other functional. Theor. Chem. Acc., 2008, 120, 215-241. (31) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. A. Gaussian 09, Revision B.01, Gaussian, Inc., Wallingford, CT. 2010. (32) Li, Q. S.; Zhang,Y.; Zhang, S. Dual Level Direct ab Initio and Density-Functional Theory Dynamics Study on the Unimolecular Decomposition of CH3OCH2 Radical J. Phys. Chem. A 2004, 108, 2014-2019. (33) Mishra, B. K.; Chakrabartty, A. K.; Bhattacharjee, D.; Deka, R. C. Theoretical investigation on unimolecular decomposition of malonic acid: a potential sink for ketene. RSC Adv., 2014, 4, 38034–38039. (34) Gonzalez, C.; Schlegel, H. B. An improved algorithm for reaction path following. J. Chem. Phys. 1989, 90, 2154-2161. (35) Dennington, I. I. R.; Keith, T.; Millam, J.; Eppinnett, K.; Hovell, W. L.; Gilliland, R. Gauss View, Version 3.09; Semichem, Inc.: Shawnee Mission, KS.2003. (36) Canneaux, S.; Bohr, F.; Hénon, E. KiSThelP: a program to predict thermodynamic properties and rate constants from quantum chemistry results, J. Comp. Chem.2014, 35, 82-93. (37) Welty, J. R.; Wicks, C. E.; Wilson, R. E.; Rorrer, G. L. Fundamentals of Momentum, Heat and Mass Transfer, 4th ed.; JohnWiley and Sons Ltd.: New York, 2001. (38) Gude, M.; Teja, A. S. Vapor-liquid critical properties of elements and compounds. 4. aliphatic alkanols, J. Chem. Eng. Data,1995, 40, 1025-1036. (39) Harding, L. B.; Georgievskii, Y.; Klippenstein, S. J. Predictive theory for hydrogen atomHydrocarbon radical association kinetics. J. Phys. Chem. A 2005, 109, 4646-4656. (40) Curran, H. J. Rate constant estimation for C-1 to C-4 alkyl and alkoxyl radical decomposition. Int. J. Chem. Kinet. 2006, 38, 250-275. (41) Yasunaga, K.; Kubo, S.; Hoshikawa, H.; Kamesawa, T.; Hidaka, Y. Shock-tube and modeling study of acetaldehyde pyrolysis and oxidation. Int. J. Chem. Kinet. 2008, 40, 73-102. (42) Kern, R. D.; Singh, H. J.; Wu, C. H. Thermal decomposition of 1,2 butadiene. Int. J. Chem. Kinet.1988, 20, 731-747.

23 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 24 of 42

(43) Baulch, D. L.; Cobos, C. J.; Cox, R. A.; Esser, C.; Frank, P.; Just, T.; Kerr, J. A.; Pilling, M. J.; Troe, J.; Walker, R. W et al. Evaluated kinetic data for combustion modelling. J. Phys. Chem. Ref. Data, 1992, 21, 411-429. (44) Natarajan, K.; Bhaskaran, K. A. Experimental and analytical investigation of high temperature ignition of ethanol. No. AD-P-000310/3. Indian Inst. of Tech., Madras. Dept. of Mechanical Engineering. 1981. (45) Cao, J. R.; Back, M. H. Kinetics of the reaction H + C2H6 → H2 + C2H5 in the temperature reaction of 1000 K. Can. J. Chem. 1984, 62, 86-91. (46) Ma, S.; Liu, R. Theoretical studies on the reaction path and dynamics of the reaction CH3+H2O → CH4 + OH. Sci. China Ser. B 1996, 39, 37-44. (47) Dean, A. M. Predictions of pressure and temperature effects upon radical addition and recombination reactions. J. Phys. Chem. 1985, 89, 4600-4608. (48) Tsang, W.; Hampson, R. F. Chemical kinetic data base for combustion chemistry. Part I. Methane and related compounds. J. Phys. Chem. Ref. Data, 1986, 15, 1087-1279. (49) Dunlop, J. R.; Tully, F. P. Catalytic dehydration of alcohols by OH. 2-propanol: an intermediate case. J. Phys. Chem. 1993, 97, 6457-6464. (50) Hidaka, Y.; Nakamura, T.; Tanaka, H.; Jinno, A.; Kawano, H. Shock tube and modeling study of propene pyrolysis. Int. J. Chem. Kinet. 1992, 24, 761-780. (51) Sutherland, J. W.; Su, M. C.; Michael, J. V. Rate constants for H + CH4, CH3 + H2, and CH4 dissociation at high temperature. Int J. Chem. Kinet. 2001, 33, 669-684. (52) Tabayashi, K.; Bauer, S. H. The early stages of pyrolysis and oxidation of methane. Combust. Flame 1979, 34, 63-83. (53) Park, J.; Zhu, R. S.; Lin, M. C. Thermal decomposition of ethanol. I. Ab Initio molecular orbital/Rice–Ramsperger–Kassel–Marcus prediction of rate constant and product branching ratios. J. Chem. Phys. 2002, 117, 3224-3230. (54) Bui, B. H.; Zhu, R. S.; Lin, M. C. Thermal decomposition of iso-propanol: First-principles prediction of total and product-branching rate constants. J. Chem. Phys. 2002, 117, 1118811195.

24 ACS Paragon Plus Environment

Page 25 of 42

F 1276 K

200

Current (pA)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

150 A

B G J

100

50

C

0

20

K

I

E D

H

40

60

80

100

120

Retention time (min) Figure 1. Gas chromatogram showing the products of the post shock mixture of the experiment carried out at 1276 K. The peaks labeled in the chromatogram are A: methane, B: ethane, C: ethylene, D: 1,1-difluoroethylene, E: acetylene, F: 1,1,1-trifluoroethane G: propylene, H: acetaldehyde, I: 1-pentene, J: 2-pentene and K: 2-pentanol. The peaks shown in black color corresponds to1,1-difluoroethylene and 1,1,1-trifluoroethane.

25 ACS Paragon Plus Environment

The Journal of Physical Chemistry

40

H 1276 K

-3

30 Absorbance×10

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 26 of 42

B 20

D E C A, D

10

C

A

F

G, B

G

0 500

1000

1500

2000

2500

3000

3500

4000

-1

wave number (cm )

Figure 2. FTIR spectrum showing the products of the post shock mixture of the experiment carried out at 1276 K. The peaks labeled in the spectrum are A: 2-pentanol, B: methane, C: ethane, D: ethylene, E: acetylene, F: acetaldehyde, G: CO and H: 1-pentene.

26 ACS Paragon Plus Environment

Page 27 of 42

0.25

1.0

(b)

(a) 0.20 [CH4]t/[C5H12O]0

[C5H12O]t/[C5H12O]0

0.8

0.6

0.4

0.15

0.10

0.2

0.05

0.0

0.00

1100

1150

1200 1250 1300 Temperature (K)

1350

1100

1150

1200 1250 Temperature (K)

1300

1350

0.4

(c)

(d)

-3

0.3

[C2H6]t/[C5H12O]0x10

[C2H4]t/[C5H12O]0

80

60

0.2

40

0.1

20

0.0

0

1100

1150

1200 1250 1300 Temperature (K)

1100

1350

1150 1200 1250 Temperature (K)

1300

1350

0.14

(e)

(f) 0.12

50 [C3H6]t/[C5H12O]0

-3

60

[C2H2]t/[C5H12O]0x10

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

40 30 20 10

0.10 0.08 0.06 0.04 0.02

0

0.00

1100

1150

1200 1250 1300 Temperature (K)

1350

1100

1150

27 ACS Paragon Plus Environment

1200 1250 Temperature (K)

1300

1350

The Journal of Physical Chemistry

30

(h)

(g) [1-C5H10]t/[C5H12O]0x10

-3

-3

80 [C2H4O]t/[C5H12O]0x10

60

40

20

0

25 20 15 10 5 0

1100

1150

1200 1250 1300 Temperature (K)

1350

1100

1150

1200 1250 Temperature (K)

1300

1350

0.12

(i) 0.10 [2-C5H10]t/[C5H12O]0

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 28 of 42

0.08 0.06 0.04 0.02 0.00 1100

1150

1200 1250 1300 Temperature (K)

1350

Figure 3. Comparison between the experimental and simulated concentrations of (a) 2-pentanol, (b) methane, (c) ethylene, (d) ethane, (e) acetylene, (f) propylene, (g) acetaldehyde, (h) 1-pentene and (i) 2-pentene. Filled circles are experimental concentrations and the unfilled circles are simulated concentrations.

28 ACS Paragon Plus Environment

Page 29 of 42

3000 2000 1500

1000

T (K) 800 700

600

500

40 20

This work (Experiment) This work (RRKM)

0

ln(k)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

-20 10 9

-40

8 7

-60

6 5

-80

4 0.75

0.80

0.85

0.90

-100 0.5

1.0

1.5 1000/T (K ) -1

2.0

Figure 4. Arrhenius plot for the overall decomposition of 2-pentanol in the temperature range of  1110-1325 K. The obtained temperature dependent rate coefficient is  (1110-1325 K) =

(4.01 ± 0.51) × 109 exp (-(36.2 kcal mol-1 ± 4.7)/RT) s-1. The insert is the zoom of the data obtained in the present experiments.

29 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

30 ACS Paragon Plus Environment

Page 30 of 42

Page 31 of 42

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 5. Geometries of the reactant, transition states, and products optimized at the M06-2X/631+G(d,p) level of theory. Blue color represents hydrogen, silver color represents carbon and red color represents oxygen in the structures. The bond lengths (Å) given on the structures are obtained at the M06-2X/6-31+G(d,p) level of theory.

31 ACS Paragon Plus Environment

The Journal of Physical Chemistry

T [K] 2500

1500

1000

800

700

40

600

500

R1 R2 R3 R4 R5 R6 R7

20 0 ln(k)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 32 of 42

-20 -40 -60 -80 -100 0.5

1.0

1.5

2.0

-1

1000/T [K ]

Figure 6. Calculated high pressure limit rate constants of reaction pathways R1-R7 over the temperature range of 500-2500 K.

32 ACS Paragon Plus Environment

Page 33 of 42

20

18

16

14 ln(k)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

12

TS1 0.01 atm 0.1 atm 1 atm 10 atm 100 atm

10

8

6 1400

1600

1800

2000 T [K]

2200

2400

2600

Figure 7. Arrhenius plot for the formation of 1-pentene via four membered ring transition state (TS1) through reaction channel R1 in the temperature range of 1300-2500 K and pressure range of 0.1-100 atm with Ar as diluent.

33 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48

Page 34 of 42

Table 1. Experimental conditions and distribution of normalized concentrations of reactant and reaction products in the decomposition of 2-pentanol.

T5

P5

(K)

(atm)

1

1325

15.9

2

1272

3

Reaction

[CH4]t/

[C2H4]t/

[C2H6]t/

[C2H2]t/

[C3H6]t/

[C2H4O]t/

[1-C5H10]t/

[2-C5H10]t/

[2-C5H12O]t/

[C5H12O]0

[C5H12O]0

[C5H12O]0

[C5H12O]0

[C5H12O]0

[C5H12O]0

[C5H12O]0

[2-C5H12O]0

[2-C5H12O]0

620

0.2254

0.3517

0.0775

0.0529

0.0797

0.0658

0.0017

0.0170

0.1284

18.1

600

0.1849

0.3177

0.0794

0.0227

0.0977

0.0499

0.0103

0.0641

0.1734

1268

18.6

650

0.1983

0.3416

0.0702

0.0208

0.1089

0.0262

0.0112

0.0528

0.1700

4

1183

16.0

557

0.0624

0.1204

0.0204

0.0022

0.0631

0.0042

0.0157

0.1143

0.5973

5

1215

16.3

600

0.0991

0.1938

0.0366

0.0044

0.0896

0.0120

0.0140

0.0988

0.4518

6

1179

16.0

518

0.0476

0.0900

0.0169

0.0021

0.0503

0.0048

0.0156

0.0910

0.6815

7

1176

16.1

544

0.0609

0.1163

0.0196

0.0013

0.0617

0.0075

0.0184

0.0934

0.6209

8

1187

16.1

462

0.0709

0.1404

0.0235

0.0022

0.0689

0.0174

0.0128

0.0977

0.5661

9

1238

16.6

594

0.1280

0.2595

0.0526

0.0092

0.0997

0.0533

0.0105

0.0679

0.3194

10

1245

16.2

515

0.1248

0.2567

0.0501

0.0080

0.0988

0.0517

0.0079

0.0652

0.3369

11

1221

16.5

483

0.0867

0.1783

0.0372

0.0024

0.0792

0.0397

0.0178

0.0679

0.4907

12

1208

16.7

477

0.0797

0.1646

0.0269

0.0028

0.0763

0.0307

0.0189

0.0713

0.5289

13

1158

15.7

542

0.0305

0.0627

0.0067

0.0007

0.0359

0.0068

0.0126

0.0414

0.8026

14

1174

16.8

454

0.0486

0.0974

0.0131

0.0014

0.0498

0.0148

0.0311

0.0601

0.6837

15

1110

16.3

471

0.0191

0.0378

0.0045

0.0006

0.0214

0.0030

0.0086

0.0384

0.8667

16

1271

16.8

590

0.1638

0.3226

0.0647

0.0191

0.1012

0.0579

0.0109

0.0406

0.2193

S.No

time (µs)

34 ACS Paragon Plus Environment

Page 35 of 42

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48

The Journal of Physical Chemistry

17

1213

15.3

507

0.0801

0.1742

0.0275

0.0038

0.0809

0.0329

0.0160

0.0506

0.5340

18

1236

15.3

598

0.1198

0.2564

0.0487

0.0080

0.0991

0.0504

0.0183

0.0508

0.3486

19

1240

15.3

600

0.1138

0.2425

0.0455

0.0084

0.0944

0.0485

0.0151

0.0428

0.3889

20

1241

15.4

570

0.1326

0.2728

0.0496

0.0105

0.1027

0.0497

0.0120

0.0427

0.3273

21

1248

15.2

602

0.1264

0.2641

0.0470

0.0076

0.0999

0.0512

0.0149

0.0682

0.3208

22

1180

17.4

464

0.0480

0.1134

0.0119

0.0014

0.0583

0.0156

0.0124

0.0469

0.6922

23

1177

16.7

528

0.0512

0.1121

0.0162

0.0013

0.0572

0.0170

0.0212

0.0836

0.6402

24

1152

16.1

476

0.0303

0.0655

0.0070

0.0011

0.0376

0.0052

0.0098

0.0363

0.8072

25

1190

13.7

475

0.0602

0.1464

0.0227

0.0017

0.0747

0.0272

0.0180

0.0487

0.6005

26

1294

17.3

580

0.2008

0.3710

0.0725

0.0322

0.0968

0.0411

0.0045

0.0243

0.1567

27

1231

12.4

521

0.0933

0.1765

0.0337

0.0113

0.0535

0.0221

0.0058

0.0317

0.5720

28

1165

16.4

512

0.0301

0.0734

0.0088

0.0009

0.0422

0.0076

0.0130

0.0285

0.7954

29

1315

17.2

590

0.2016

0.3772

0.0751

0.0400

0.0860

0.0378

0.0046

0.0193

0.1586

30

1266

16.0

612

0.1461

0.3116

0.0577

0.0101

0.1015

0.0504

0.0164

0.0490

0.2573

31

1284

16.9

590

0.1731

0.3523

0.0745

0.0228

0.0950

0.0612

0.0076

0.0258

0.1877

32

1225

14.6

550

0.0998

0.2428

0.0373

0.0069

0.1006

0.0501

0.0171

0.0718

0.3736

33

1140

15.5

521

0.0145

0.0317

0.0042

0.0001

0.0236

0.0055

0.0116

0.1010

0.8077

34

1127

15.6

415

0.0198

0.0368

0.0035

0.0005

0.0226

0.0060

0.0105

0.0408

0.8595

35

1276

15.9

626

0.1781

0.3522

0.0667

0.0196

0.1036

0.0543

0.0088

0.0353

0.1813

36

1310

15.6

616

0.1926

0.3677

0.0768

0.0333

0.0929

0.0635

0.0077

0.0181

0.1473

37

1253

14.3

600

0.1475

0.3278

0.0629

0.0115

0.1124

0.0561

0.0141

0.0340

0.2337

38

1324

15.4

590

0.1963

0.3740

0.0791

0.0357

0.0894

0.0704

0.0076

0.0177

0.1299

35 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Table 2. Energy barriers (kcal mol-1) for thermal decomposition of 2-pentanol reaction via pathways R1-R7 calculated at M06-2X/6-31+G(d,p).

TSs

M06-2X/6-31+G(d,p)

TS1

64.6

TS2

66.9

TS3

101.1

TS4

112.8

TS5

84.1

TS6

93.6

TS7

84.7

36 ACS Paragon Plus Environment

Page 36 of 42

Page 37 of 42

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48

The Journal of Physical Chemistry

Table 3. Proposed reaction scheme for the decomposition of 2-pentanol with 39 reaction species and 37 elementary reactions.a

S.No Reaction

A

n

Ea

Reference

R1

CH3-CH2-CH2-CH(OH)-CH3 → CH3-CH2-CH2-CH=CH2 + H2O

1.58 × 1014

0.00

67.1

This work

R2

CH3-CH2-CH2-CH(OH)-CH3 → CH3-CH2-CH=CH-CH3 + H2O

3.01 × 1014

0.00

69.5

This work

R3

CH3-CH2-CH2-CH(OH)-CH3 → C2H6 + CH3-C(O)-CH3

1.90 × 1014

0.00

103.4

This work

14

R4

CH3-CH2-CH2-CH(OH)-CH3 → C2H4 + CH3-CH(OH)-CH3

6.01 × 10

0.00

115.9

This work

R5

CH3-CH2-CH2-CH(OH)-CH3 → CH4 + CH3-CH2-CH2-CHO

1.92 × 1014

0.00

86.2

This work

R6

CH3-CH2-CH2-CH(OH)-CH3 → CH4 + C2H4 + CH3-CHO

1.60 × 1014

0.00

96.1

This work

R7

CH3-CH2-CH2-CH(OH)-CH3 → C3H8 + CH3-CHO

9.44 × 1013

0.00

86.4

This work

R8

CH3-CH2-CH2-CH(OH)-CH3 → CH3 + CH2-CH2-CH(OH)-CH3

5.12 × 1026

-2.92

90.5

23

R9

CH3-CH2-CH2-CH(OH)-CH3 → CH3-CH2 + CH2-CH(OH)-CH3

8.54 × 1027

-3.23

91.7

23

R10

CH3-CH2-CH2-CH(OH)-CH3 → CH3-CH2-CH2 + CH(OH)-CH3

3.01 × 1025

-2.57

89.9

23

R11

CH3-CH2-CH2-CH(OH)-CH3 → CH3-CH2-CH2-CH-CH3 + OH

2.82 × 1029

-3.60

101.9

23

R12

CH3-CH2-CH2-CH(OH)-CH3 → CH3-CH2-CH2-CH(OH) + CH3

5.13 × 1026

-2.91

90.5

23

R13

CH3-CH2 + H → C2H6

5.44 × 1013

0.16

0.0

39

R14

CH3-CH2 → C2H4 + H

10

3.06 × 10

0.95

36.9

40

R15

CH3-CH2-CH2 → C3H6 + H

4.14 × 1012

0.17

35.6

40

R16

CH3 + H → CH4

6.93 × 1013

0.18

0.0

39

37 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48

Page 38 of 42

R17

CH3-C(O)-H → CH4 + CO

1.00 × 1015

0.00

85.0

41

R18

2 CH3 → C2H4 + H2

9.90 × 1015

0.00

32.9

42

R19

C2H4 + CH3 → C2H3 + CH4

4.16 × 1012

0.00

11.1

43

R20

C2H3 + H → C2H4

3.88 × 1013

0.20

0.00

39

R21

CH(OH)-CH3 → CH3-C(O)-H + H

5.00 × 1013

0.00

21.9

44

R22

C2H3 + H → C2H2 + H2

1.21 × 1013

0.00

0.00

43

R23

C2H6 + H → C2H5 + H2

5.95 × 1014

0.00

12.9

45

R24

CH3-CH2-CH2-CH=CH2 → CH3+ CH2-CH2-CH=CH2

1.00 × 1016

0.00

72.9

47

R25

CH3-CH2-CH2-CH=CH2 → CH3-CH2 + CH2-CH=CH2

1.00 × 1016

0.00

70.7

47

R26

CH3 + H2O → CH4 + OH

0.64 × 101

3.31

12.6

46

17

R27

CH3-CH2-CH=CH-CH3 → CH3 + CH3-CH2-CH=CH

3.16 × 10

0.00

99.4

47

R28

CH3-CH2-CH=CH-CH3 → CH3 + CH2-CH=CH-CH3

3.16 × 1017

0.00

99.4

47

R29

CH3-CH2-CH=CH-CH3 → CH3-CH2 + CH=CH-CH3

1.00 × 1016

0.00

75.7

47

R30

CH3-CH2-CH2-CH-CH3→ C3H6 + CH3-CH2

1.22 × 1012

0.64

29.4

40

R31

C3H5 → C2H2 + CH3

1.26 × 1013

0.00

33.7

47

R32

C2H5 + H2O → C2H6 + OH

3.39 × 106

1.44

20.3

48

R33

CH2-CH(OH)-CH3 → CH3-CH=CH2 + OH

3.70 × 1012

0.00

26.8

49

R34

C3H6 → CH3 + C2H3

8.00 × 1014

0.00

88.0

50

R35

C3H6 + H → n-C3H7

4.57 × 1012

0.51

2.6

40

R36

C3H6 + H → iso-C3H7

7.75 × 1012

0.51

1.2

40

38 ACS Paragon Plus Environment

Page 39 of 42

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48

The Journal of Physical Chemistry

a

R37

CH4 + H → H2 + CH3

1.77 × 1014

0.00

13.8

51

R38

CH4 + CH3 → H2 + CH3-CH2

1.00 × 1013

0.00

23.1

52

R39

C2H2 + H2 → C2H4

3.01 × 1011

0.00

38.9

48

Rate expressions are given in the form of k = A exp(-Ea/RT) and k = ATnexp(-Ea/RT). The units of the rate coefficients are s-1 and

cm3mol-1s-1 for first and second order reactions respectively. The units for the activation barrier are kcal mol-1.

39 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48

Page 40 of 42

Table 4. The percentage changes in the concentrations of individual species, when the rate coefficients of the reactions in the proposed reaction scheme were varied by a factor of 10 in order to analyze the sensitivity towards the formation of all the products at 1268 K. No

Reaction

CH4

C2H4

C3H6

1-C5H10

CH3CHO

C2H6

C2H2

2-C5H10

R1

CH3-CH2-CH2-CH(OH)-CH3 → CH3-CH2-CH2-CH=CH2 + H2O

101

149

-58

15

-29

202

-19

-66

R2

CH3-CH2-CH2-CH(OH)-CH3 → CH3-CH2-CH=CH-CH3 + H2O

114

-89

-51

-14

-26

138

181

377

R5

CH3-CH2-CH2-CH(OH)-CH3 → CH4 + CH3-CH2-CH2-CHO

0

0

0

2

-1

-1

0

0

R6

CH3-CH2-CH2-CH(OH)-CH3 → CH4 + C2H4 + CH3-CHO

0

0

0

3

-1

-1

0

0

R8

CH3-CH2-CH2-CH(OH)-CH3 → CH3 + CH2-CH2-CH(OH)-CH3

178

-130

-22

-6

-15

10

-5

-31

R9

CH3-CH2-CH2-CH(OH)-CH3 → CH3-CH2 + CH2-CH(OH)-CH3

-30

306

354

-8

-12

35

-20

-31

R10

CH3-CH2-CH2-CH(OH)-CH3 → CH3-CH2-CH2 + CH(OH)-CH3

-49

19

306

-8

317

-41

-23

-25

R11

CH3-CH2-CH2-CH(OH)-CH3 → CH3-CH2-CH2-CH-CH3 + OH

-4

-10

-22

0

0

0

0

0

R12

CH3-CH2-CH2-CH(OH)-CH3 → CH3-CH2-CH2-CH(OH) + CH3

180

-130

-28

-6

-13

9

-6

-27

R14

CH3-CH2 → C2H4 + H

6

48

4

0

-4

-41

0

0

R15

CH3-CH2-CH2 → C3H6 + H

0

0

0

4

-1

-3

0

0

R19

C2H4 + CH3 → C2H3 + CH4

8

-18

1

0

0

0

5

1

R21

CH(OH)-CH3 → CH3-C(O)-H + H

0

0

0

3

0

-1

0

0

R22

C2H3 + H → C2H2 + H2

1

-23

-2

5

-2

4

19

0

R23

C2H6 + H → C2H5 + H2

3

38

-1

3

1

-35

-6

-1

R24

CH3-CH2-CH2-CH=CH2 → CH3+ CH2-CH2-CH=CH2

72

-133

3

-15

-2

-5

10

-1

R25

CH3-CH2-CH2-CH=CH2 → CH3-CH2 + CH2-CH=CH2

-34

78

0

-17

-2

-6

-9

0

40 ACS Paragon Plus Environment

Page 41 of 42

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48

The Journal of Physical Chemistry

R26

CH3 + H2O → CH4 + OH

6

16

-6

2

-4

-9

-4

5

R27

CH3-CH2-CH=CH-CH3 → CH3 + CH3-CH2-CH=CH

0

0

0

0

0

0

-3

2

R28

CH3-CH2-CH=CH-CH3 → CH3 + CH2-CH=CH-CH3

0

0

0

0

0

0

-3

2

R29

CH3-CH2-CH=CH-CH3 → CH3-CH2 + CH=CH-CH3

87

48

6

5

-5

37

69

-85

R31

C3H5 → C2H2 + CH3

0

0

0

0

0

0

-1

0

R32

C2H5 + H2O → C2H6 + OH

-11

-147

-2

5

-2

111

6

3

R34

C3H6 → CH3 + C2H3

0

0

0

0

0

0

0

0

R37

CH4 + H → H2 + CH3

-8

-38

2

1

0

18

0

0

R38

CH4 + CH3 → H2 + CH3-CH2

-43

31

-3

0

0

0

0

0

R39

C2H2 + H2 → C2H4

3

13

-4

0

0

0

-21

-1

41 ACS Paragon Plus Environment

The Journal of Physical Chemistry

Table of content (TOC)

3000 2000 1500

1000

T (K) 800 700

600

500

40 20

This work (Experiment) This work (RRKM)

0

ln(k)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 42 of 42

-20 10 9

-40

8 7

-60

6 5

-80

4 0.75

0.80

0.85

0.90

-100 0.5

1.0

1.5 1000/T (K ) -1

42 ACS Paragon Plus Environment

2.0