Xiaofeng Tang*a,
Xiaoxiao Lina,
Yupeng Zhuab,
Xiangkun Wubc,
Zuoying Wenab,
Lidong Zhangd,
Fuyi Liud,
Xuejun Gua and
Weijun Zhang*ae
aLaboratory of Atmospheric Physico-Chemistry, Anhui Institute of Optics and Fine Mechanics, Chinese Academy of Sciences, Hefei, 230031 Anhui, China. E-mail: tangxf@aiofm.ac.cn; wjzhang@aiofm.ac.cn
bGraduate School, University of Science and Technology of China, Hefei, 230026 Anhui, China
cHefei National Laboratory for Physical Sciences at the Microscale, Department of Chemical Physics, University of Science and Technology of China, Hefei 230026, China
dNational Synchrotron Radiation Laboratory, University of Science and Technology of China, Hefei 230029, China
eSchool of Environmental Science and Optoelectronic Technology, University of Science and Technology of China, Hefei, 230026 Anhui, China
First published on 31st May 2017
We present here the coupling of a flash pyrolysis micro-reactor with threshold photoelectron photoion coincidence (TPEPICO) spectroscopy at Hefei synchrotron to investigate thermal decomposition of n-butane. Primary products in the pyrolysis were determined with TPEPICO time-of-flight mass spectra and mass-selected threshold photoelectron spectra (TPES). Vibrational structures in the TPES of the CH3, C2H5 and C3H6 products were observed and their ionization energies were measured and compared satisfactorily to literature data. Moreover, the potential energy surface involved in the decomposition was calculated by using high-level quantum chemistry and several transition states were recognized. Then the detailed mechanisms of the pyrolysis of n-butane have been discussed and revealed.
The objective of the present study, n-butane is the prototype of hydrocarbons and an essential model fuel to study combustion behaviors. The thermal decomposition reactions of n-butane are important initiation steps in the high-temperature oxidation and strongly influence on the ignition time,16 which have attracted a lot of attention in the past decades. A series of products from both the primary and the secondary pyrolysis reactions of n-butane have been identified and analyzed by using the methods of GC, Raman spectroscopy, electron impact ionization mass spectrometry and PIMS.6,7,17–20 For example, for the major stable products of CH4, C2H4, C2H6 and C3H6, it was believed that the thermal decomposition reactions responsible for the production were mostly proceeded through free-radical chain reactions.6 But, these stable products may also be formed from the direct decomposition of n-butane and in recent quantum chemical calculations several transition states were determined on the potential energy surfaces of the dissociation of n-butane and can correlate to the production of the stable products.21
Recently, we have developed a molecular beam flash pyrolysis micro-reactor22–24 in combination with a threshold PEPICO (TPEPICO) spectrometer10 at Hefei synchrotron to explore VUV spectroscopy of organic radicals25 and thermal decomposition of hydrocarbons. Benefited from its high temperature and short contact time, the flash pyrolysis micro-reactor allows for the thermal decomposition with high efficiency and minimizes secondary reactions.26–28 The TPEPICO scheme provides a strategy to obtain pure spectra of products without contamination from other byproducts. As a representative example, thermal decomposition of n-butane has been selected and investigated by using the flash pyrolysis micro-reactor and the TPEPICO spectrometer. The primary pyrolysis products were identified and determined from TPEPICO time-of-flight (TOF) mass spectra and mass-selected TPES. In addition, the potential energy surface involved in the decomposition was also calculated by using high-level quantum chemistry and several transition states were recognized. With the results of experiments and theoretical calculations, a clear picture of the dissociation mechanisms has been revealed.
Synchrotron photons emitted from an undulator of 800 MeV electron ring at NSRL were dispersed by a 6 m monochromator equipped with three gratings (370, 740 and 1250 grooves per mm) covering the photon energy range of 7.5–124 eV. To maximize photon flux, presently the slits of the monochromator were expanded with widths of 200 μm and the photon energy resolution is about 20 meV (full width at half maximum, FWHM) at hν = 10 eV with the 370 grooves per mm grating. A gas filter filled with argon was installed to suppress the high-order harmonic photons emitted from the undulator. The absolute photon energy of the monochromator was calibrated on-line with the absorption resonant lines of Ar in the gas filter. The photons were softly focused onto a spot with ∼1 mm (V) × 2 mm (H) size in the photoionization region. A photodiode just located behind the photoionization region was used to measure the photon flux for normalizing signals in the photon energy scans.
The TPEPICO spectrometer10 is mainly composed of two TOF analyzers with double velocity map imaging29 designs for ions and electrons respectively operated in coincidence. A mask with a small central hole and a concentric ring was installed and located at the front of electron micro-channel plates (MCPs) with a chevron configuration. The small central hole and the concentric ring were used to collect near zero kinetic energy photoelectrons (threshold photoelectrons) and hot electrons, respectively, which were velocity focused as shown by the black and green trajectories in Fig. 1.30 After amplified, the electron signals were transferred to the start connector of a multiple-event time digitizers (P7888-2, FAST Comtec, Germany) and the ion signals to the stop connector. Then the electron signals provide start to record ion events in coincidence. The contribution of hot electrons that happen to the center of the electron detector was corrected by using the method of subtraction.30 TPES, especially mass-selected TPES corresponding to each species can be acquired in the coincidence mode by scanning synchrotron photon energy.
The flash pyrolysis micro-reactor25 is similar to the design of Chen et al.26 and primarily consists an orifice nozzle (50 μm diameter), a ceramic insulation plate and a silicon carbide tube (SiC, 1 mm inner diameter). The jet of n-butane seeded in argon with a stagnation pressure of 2 atm in the nozzle was adiabatically expanded into the SiC tube in vacuum and a continuous molecular beam with a subsonic speed was formed. The SiC tube was electrically connected to a DC power supply via a vacuum feedthrough port and the pyrolysis power can be controlled outside. The temperature inside the SiC tube was measured off-line with a thermocouple and can be approached at ∼1600 K. The high temperature and the short contact time (∼several tens of microseconds) due to the short heating length of the SiC tube (∼15 mm) and the fast speed of molecular beam allow for the thermal decomposition with high efficiency. Moreover, the flash pyrolysis micro-reactor was mounted on a manually-controlled XYZ manipulator in the source chamber of the TPEPICO spectrometer and its position can be tuned real-time to optimize the signals. The source and ionization chambers of the TPEPICO spectrometer were separated by a skimmer (1 mm diameter) and their backing pressures were 4 × 10−2 Pa and 5 × 10−5 Pa respectively with the molecular beam on.
Fig. 2 TPEPICO TOF mass spectra recorded at fixed photon energy of hν = 10.60 eV with (a) no pyrolysis, (b) pyrolysis power P = 29 W, and (c) pyrolysis power P = 36 W. |
m/z | Species | IE (eV) | Methods | Ref. |
---|---|---|---|---|
a PES: photoelectron spectroscopy; PFI-ZEKE: pulsed-field-ionization zero-kinetic-energy photoelectron spectroscopy. | ||||
58 | n-C4H10 | 10.53 ± 0.02 | Thermochemistry | 31 |
57 | 1-C4H9 | 8.02 ± 0.04 | PES | 31 |
43 | 1-C3H7 | 8.09 ± 0.01 | PES | 33 |
42 | C3H6 | 9.7452 ± 0.0005 | PFI-ZEKE | 34 |
30 | C2H6 | 11.5 ± 0.1 | TPES | 35 |
29 | C2H5 | 8.117 ± 0.008 | PIMS | 36 |
28 | C2H4 | 10.51 ± 0.02 | PES | 37 |
16 | CH4 | 12.618 ± 0.004 | PFI-ZEKE | 38 |
15 | CH3 | 9.8387 ± 0.0018 | PFI-ZEKE | 39 |
In the mass spectra the widths of the CH3+, C2H4+, C2H5+, and C3H6+ peaks are very narrow, similar to the peak of C4H10+ parent ion, indicating that they were also from photoionization of their neutral molecules. Their intensities increase with pyrolysis power. The peaks of C3H7+ are much wider than those of the other ions and the C3H7+ ions may be produced from different processes. Photoionization and dissociative photoionization of n-butane have been studied well previously by the method of photoionization mass spectrometry and C3H7+ is the major fragment ion with the largest intensity in the dissociation.32 The appearance energy of C3H7+ in dissociative photoionization of n-butane was measured at 11.18 eV (ref. 32) which is above the present photon energy of hν = 10.60 eV. But, in the present experimental condition to arrive the photoionization region n-butane needed to go past the high-temperature SiC tube and the increased temperature can shift the onset of the dissociative photoionization. This effect of hot band has also been observed in the TPES of CH3, as shown in Fig. 4(a). Some kinetic energies should have been released in the dissociative photoionization and contribute to the wide width of the C3H7+ peak in the mass spectra.
The pyrolysis of n-butane may also produce the neutral C3H7 radical through the direct fission of the C1–C2 bond of n-butane in the SiC tube and the above observed CH3 can be its accompanied product. The signal of the C3H7 product is very weak, as demonstrated by the weak CH3+ signal in the mass spectra of Fig. 2, and may be buried in the C3H7+ wide peak. So we tried to change synchrotron photon energy and look for the evidence of the C3H7 product in TPEPICO TOF mass spectra. For example, with the fixed pyrolysis power of P = 36 W, TPEPICO TOF mass spectra were measured at fixed photon energies of hν = 9.50, 9.90 and 10.20 eV, all of which are far below the appearance energy of C3H7+ in dissociative photoionization of n-butane, and are displayed in Fig. 3. Except the absence of the C4H10+ parent ion, the peaks of CH3+, C2H4+, C2H5+, C3H6+ and C3H7+ ions can be determined in Fig. 3. All of them are very narrow in the mass spectra. Especially, in contrary to the case of Fig. 2, presently the peaks of C3H7+ in the mass spectra do not take a wide width again indicating that they were only from photoionization of the propyl radicals.
Fig. 3 TPEPICO TOF mass spectra recorded at different photon energies with pyrolysis power of P = 36 W. (a) hν = 9.50 eV, (b) hν = 9.90 eV and (c) hν = 10.20 eV. |
Fig. 4 Mass-selected threshold photoelectron spectra (TPES) of (a) CH3, (b) C2H5 and (c) C3H6 products from pyrolysis of n-butane. |
The TPES of ethyl radical is exhibited with a broad band in Fig. 4(b). The electron signal appears at ∼8.2 eV and its intensity approaches the maximal at hν = 8.7 eV, contributing to the adiabatic and vertical ionization energies of ethyl radical, respectively. For propene product, as shown in Fig. 4(c), two vivid vibrational peaks can be observed in the TPES. The adiabatic ionization energy of propene is measured at 9.75 eV and corresponds to the first peak. The other peak with a spacing of 165 meV (1330 cm−1) to the first peak is ascribed to the CC stretching vibrational population of C3H6+.34,41 We had also tried to measure the TPES of the propyl radical from pyrolysis of n-butane. But, as we discussed in the above part of 3.1, the signals of C3H7+ were from two processes, photoionization of the propyl radical and dissociative photoionization of n-butane, and presently it is very difficult for us to discriminate them due to the restriction of the detectors. This problem can be overcome by employing position sensitive detectors (PSDs) to collect electrons and ions in coincidences and then the pure TPES of the propyl radical could be acquired with an ion-kinetic-energy-filtered method.42,43
C4H10 + Δ → C4H9 + H | (1) |
C4H9 → CH3 + C3H6 | (2) |
C4H9 → C2H5 + C2H4 | (3) |
CH3 + C4H10 → CH4 + C4H9 | (4) |
C2H5 + C4H10 → C2H6 + C4H9 | (5) |
CH3, C2H5 and C4H9 free radicals are the intermediates of the reactions. But, the above stable products can be formed directly via the primary decomposition reactions of n-butane too, not involving the secondary reactions.
C4H10 + Δ → CH4 + C3H6 | (6) |
C4H10 + Δ → C2H6 + C2H4 | (7) |
In addition, the primary decomposition reactions of n-butane can also produce CH3, C2H5 and C3H7 radical products.
C4H10 + Δ → CH3 + C3H7 | (8) |
C4H10 + Δ → C2H5 + C2H5 | (9) |
The successful identification of the reaction intermediates in experiment, especially the reactive radicals, is the key to reveal and demonstrate its decomposition mechanisms. The thermal decomposition of n-butane via the reactions (8) and (9) had been approved in shock tube by Hanson et al.4 In previous flow reactor experiment CH3 and C2H5 radical products, with no C3H7 and C4H9 radicals, were also detected by synchrotron VUV PIMS.19
As shown in Fig. 2 and 3, in the present flash pyrolysis study of n-butane the products including the stable species of C2H4 and C3H6 and the radicals of CH3, C2H5 and C3H7 have been observed in the TPEPICO TOF mass spectra. C3H7 radical products are directly and successfully identified in the mass spectra for the first time. No C4H9 radical products, whose ionization energy is located at 8.02 eV,31 was observed in our experiments. Due to the high temperature and the short contact time in the SiC tube, the thermal decomposition of n-butane can be taken place with high efficiency. So we believe that all the products observed in our experiments were mainly formed from the reactions of (6)–(9).
To elucidate the present experimental results, the decomposition of n-butane was also investigated by using quantum chemical calculations. The method of Becke three-parameter exchange functional with the Lee, Yang, and Parr correlation functional (B3LYP) with the 6-311+G(d,p) basis set was employed for molecular geometry optimizations and harmonic vibrational frequency analysis.45 The intrinsic reaction coordinate (IRC) calculations were also carried out to probe transition states between the reactants and products and to characterize the reaction profiles and elucidate the decomposition channels.46,47 Then the high-level method of coupled cluster method with single, double, and non-iterative triple excitations (CCSD(T))48 with the aug-cc-pVTZ basis set was adopted to refine electronic energies based on the optimized B3LYP/6-311+G(d,p) structures. The corrections of zero point vibrational energies (ZPE) had been taken into account with the method of B3LYP/6-311+G(d,p). All the calculations were performed with the Gaussian 09 program.49
The potential energy surface together with the optimized molecular structures and their relative energies to the ground state of n-butane is presented in Fig. 5. It is shown that n-butane has two conformers with trans and gauche structures. The energy of the trans conformer is a litter smaller than that of the gauche and has a more stable configuration. The isomerization between the trans and gauche conformers can be proceeded via a transition state (TS1) with a small height barrier of 0.14 eV. The barrier is so small that the products in the high temperature pyrolysis can be formed from the both conformers.
The pyrolysis pathways observed in the present experiments were calculated and then the decomposition mechanisms of n-butane can be inferred. As shown in Fig. 5, CH3 and its accompanying C3H7 radical products corresponding to the reaction (8) were produced from the direct fission of the C1–C2 bond of n-butane and no barrier was found on the decomposition route. The C1–C2 bond energy of n-butane is calculated to be 3.76 eV. Similarly, two C2H5 radical products were formed directly through the C2–C3 bond fission of n-butane and the bond energy is 3.70 eV. But, the dissociation mechanisms to produce the C2H4 or C3H6 products observed in experiments are a litter complicated than those of the CH3 and C2H5 products. The calculations show that the process of H atom migration should be involved in the decomposition of n-butane to produce CH4 and C3H6 species. As shown in Fig. 5, this speculation is attested and the two transition states (TS2 and TS3) with a complex structure of CH4–CHCH2CH3 were also identified along the decomposition routes. For the productions of C2H4 and C2H6 species, two four-member-ring transition states of TS4 and TS5 are determined along the routes. The energy barriers for this production channel from the trans and gauche conformers of n-butane are calculated at 4.85 and 4.82 eV, respectively.
The potential energy surface related to the thermal decomposition of n-butane has been calculated with the high-level theoretical method of CCSD(T)/aug-cc-pVTZ//B3LYP/6-311+G(d,p) + ZPE and the detailed dissociation mechanisms have been discussed and revealed. All the products observed in our experiments were mainly formed from the primary decomposition reactions of n-butane. The trans and gauche conformers of n-butane should be involved in the pyrolysis. The productions of the CH3 radical and its accompanying C3H7 products, and two C2H5 radical products, are simple and just through the direct bond fission of n-butane. But, the dissociation mechanisms to produce C3H6 or C2H4 fragments are a litter more complicated than those of the CH3 and C2H5 products. Transition states with four-member-ring structures were determined in the production of C2H4 and C2H6. In addition, to explain the pyrolysis to the CH4 and C3H6 products, a mechanism of H atom migration should be taken into account.
This journal is © The Royal Society of Chemistry 2017 |