Pyrolysis of furan in a microreactor

Search for transition states with external forces

It is much more difficult to find on the potential energy surface (PES) a transition state (TS) than a local minimum (LM). We propose a new methodology which makes this task much easier. Applying external forces to nuclei in a molecule we can locate on PES not just one particular TS but a number of consecutive transition states. With our approach it is possible to move over PES from one transition state to another without involving any local minima. The latter can be located in a separate step through the reaction path calculations performed for every transition state found before. Preliminary results for the 2-fluorofuran molecule illustrate the usefulness of the proposed method.

Gas-phase synthesis of benzene via the propargyl radical self-reaction

Polycyclic aromatic hydrocarbons (PAHs) have been invoked in fundamental molecular mass growth processes in our galaxy. We provide compelling evidence of the formation of the very first ringed aromatic and building block of PAHs-benzene-via the self-recombination of two resonantly stabilized propargyl (C₃H₃) radicals in dilute environments using isomer-selective synchrotron-based mass spectrometry coupled to theoretical calculations. Along with benzene, three other structural isomers (1,5-hexadiyne, fulvene, and 2-ethynyl-1,3-butadiene) and o-benzyne are detected, and their branching ratios are quantified experimentally and verified with the aid of computational fluid dynamics and kinetic simulations. These results uncover molecular growth pathways not only in interstellar, circumstellar, and solar systems environments but also in combustion systems, which help us gain a better understanding of the hydrocarbon chemistry of our universe.

Formation of polycyclic aromatic hydrocarbons, benzofuran, and dibenzofuran in fuel-rich oxidation of toluene using a flow reactor

Recently, polycyclic aromatic hydrocarbons (PAHs) and oxygenated PAHs (OPAHs) have been attracting considerable attention owing to their high toxicity. Understanding their formation mechanism during combustion processes is important to control their emission. However, there are few studies that have quantitatively investigated OPAH formation in the fuel-rich oxidation of hydrocarbons, despite the availability of several studies on PAH formation. In this study, benzofuran and dibenzofuran as OPAHs were quantified in the fuel-rich oxidation of toluene using a flow reactor at atmospheric pressure in a temperature range of 1050-1350 K at equivalence ratios from 3.0 to 12.0 and residence times from 0.2 to 1.5 s. In addition to benzofuran and dibenzofuran, 4 types of monocyclic aromatic hydrocarbons and 19 types of PAHs were also evaluated. The experimental data obtained in this study were compared with those of the ethylene oxidation performed in our previous study. The existing kinetic model for PAH growth was modified based on several theoretical studies to predict the behavior of OPAHs with furan structures. The modified model showed significant improvements in the prediction of benzofuran and dibenzofuran formation. Based on the rate of production and sensitivity analysis using the modified model, the dominant reaction pathways of benzofuran and dibenzofuran were investigated.

Isomerization and Decomposition of 2-Methylfuran with External Forces

The primary goal of this project was to evaluate the performance of the Standard and Enforced Geometry Optimization (SEGO) method which we have recently developed. The SEGO method has been designed for an automatic location of multiple minima on the molecular Potential Energy Surface (PES), and its usefulness has been demonstrated so far for three molecules only. In this project we applied the SEGO method to explore the 2-methylfuran (2MF) PES. Our choice was not accidental: this molecule recently gained a great deal of interest as a potential candidate for biofuel, and therefore its pyrolysis is extensively studied. To understand pyrolysis of 2MF a detailed knowledge about its PES is needed. In these studies we explored the 2MF PES and located a surprisingly large number of local minima corresponding to 2MF isomers and decomposition products. Some of the 2MF isomers and fragments have amazing structures which most likely were previously unknown. All structures presented in this paper were found in an automatic manner as a result of the enforced chemical reactions driven by the SEGO method. Thus, in these studies we have not only proven that the SEGO method is a very efficient tool for exploring molecular potential energy surfaces, but we also obtained very detailed knowledge about the 2MF PES.

Reactivity of the Indenyl Radical (C9 H7 ) with Acetylene (C2 H2 ) and Vinylacetylene (C4 H4 )

The reactions of the indenyl radicals with acetylene (C₂ H₂ ) and vinylacetylene (C₄ H₄ ) is studied in a hot chemical reactor coupled to synchrotron based vacuum ultraviolet ionization mass spectrometry. These experimental results are combined with theory to reveal that the resonantly stabilized and thermodynamically most stable 1-indenyl radical (C₉ H₇ ^. ) is always formed in the pyrolysis of 1-, 2-, 6-, and 7-bromoindenes at 1500 K. The 1-indenyl radical reacts with acetylene yielding 1-ethynylindene plus atomic hydrogen, rather than adding a second acetylene molecule and leading to ring closure and formation of fluorene as observed in other reaction mechanisms such as the hydrogen abstraction acetylene addition or hydrogen abstraction vinylacetylene addition pathways. While this reaction mechanism is analogous to the bimolecular reaction between the phenyl radical (C₆ H₅ ^. ) and acetylene forming phenylacetylene (C₆ H₅ CCH), the 1-indenyl+acetylene→1-ethynylindene+hydrogen reaction is highly endoergic (114 kJ mol^-1 ) and slow, contrary to the exoergic (-38 kJ mol^-1 ) and faster phenyl+acetylene→phenylacetylene+hydrogen reaction. In a similar manner, no ring closure leading to fluorene formation was observed in the reaction of 1-indenyl radical with vinylacetylene. These experimental results are explained through rate constant calculations based on theoretically derived potential energy surfaces.

An Experimental and Theoretical Investigation of 1-Butanol Pyrolysis

Bioalcohols are a promising family of biofuels. Among them, 1-butanol has a strong potential as a substitute for petrol. In this manuscript, we report on a theoretical and experimental characterization of 1-butanol thermal decomposition, a very important process in the 1-butanol combustion at high temperatures. Advantage has been taken of a flash pyrolysis experimental set-up with mass spectrometric detection, in which the brief residence time of the pyrolyzing mixture inside a short, resistively heated SiC tube allows the identification of the primary products of the decomposing species, limiting secondary processes. Dedicated electronic structure calculations of the relevant potential energy surface have also been performed and RRKM estimates of the rate coefficients and product branching ratios up to 2,000 K are provided. Both electronic structure and RRKM calculations are in line with previous determinations. According to the present study, the H₂O elimination channel leading to 1-butene is more important than previously believed. In addition to that, we provide experimental evidence that butanal formation by H₂ elimination is not a primary decomposition route. Finally, we have experimental evidence of a small yield of the CH₃ elimination channel.

Computational insight into excited states of the ring-opening radicals from the pyrolysis of furan biofuels

The low-lying valence excited states and Rydberg states of the radical species from the ring-opening reactions in pyrolysis of furan biofuels have been determined by extensive density functional theory and sophisticated wave function theory calculations. The radicals 1-C₄ H₅ O-2, 2-furylCH₂ , and 4-C₆ H₇ O with the delocalized π-type single electron are predicted to be most stable among the reactive species here for furan, 2-methyfuran, and 2,5-dimethylfuran, respectively. Predicted vertical transition energies by TD-CAM-B3LYP show good agreement with those by CASPT2. Some among the electronic excitations to low-lying states can take place in the visible light region, and they may be involved in the combustion process. Further surface hopping dynamics simulations on the excited states of the most stable ring-opening radical 1-C₄ H₅ O-2 of furan as an example reveal that 89.9% sampling trajectories at the initial excited state of 2² A"(π¹ π*² ) decay to the 1² A'(n¹ π*² ) state within an average of 384 fs, and then 81.2% trajectories at the 1² A' state go to the ground state within an average of 114 fs. At the end of the simulation for 1000 fs, 18.8% trajectories still stay on the excited states of 2² A" and 1² A', suggesting that the reactive radicals in the ground state are mainly responsible for the combustion chemistry of furan biofuels. © 2018 Wiley Periodicals, Inc.

How to add a five-membered ring to polycyclic aromatic hydrocarbons (PAHs) – molecular mass growth of the 2-naphthyl radical (C10H7) to benzindenes (C13H10) as a case study

The reaction of aryl radicals with allene/methylacetylene leads to five-membered ring addition in PAH growth processes.

Furan Fragmentation in the Gas Phase: New Insights from Statistical and Molecular Dynamics Calculations

We present a complete exploration of the different fragmentation mechanisms of furan (C₄H₄O) operating at low and high energies. Three different theoretical approaches are combined to determine the structure of all possible reaction intermediates, many of them not described in previous studies, and a large number of pathways involving three types of fundamental elementary mechanisms: isomerization, fragmentation, and H/H₂ loss processes (this last one was not yet explored). Our results are compared with the existing experimental and theoretical investigations for furan fragmentation. At low energies the first processes to appear are isomerization, which always implies the breaking of one C-O bond and one or several hydrogen transfers; at intermediate energies the fragmentation of the molecular skeleton becomes the most relevant mechanism; and H/H₂ loss is the dominant processes at high energy. However, the three mechanisms are active in very wide energy ranges and, therefore, at most energies there is a competition among them.

Ammonolysis of ketene as a potential source of acetamide in the troposphere: a quantum chemical investigation

Quantum chemical calculations at the CCSD(T)/CBS//MP2/aug-cc-pVTZ levels of theory have been carried out to investigate a potential new source of acetamide in Earth's atmosphere through the ammonolysis of the simplest ketene.

A vacuum ultraviolet photoionization study on high-temperature decomposition of JP-10 (exo-tetrahydrodicyclopentadiene)

High-temperature pyrolysis of JP-10 in flow reactors were performed both experimentally and theoretically. Dozens of products were detected and the decomposition pathways of JP-10 were discussed.

Hydrogen-Abstraction/Acetylene-Addition Exposed

Polycyclic aromatic hydrocarbons (PAHs) are omnipresent in the interstellar medium (ISM) and also in carbonaceous meteorites (CM) such as Murchison. However, the basic reaction routes leading to the formation of even the simplest PAH-naphthalene (C₁₀ H₈ )-via the hydrogen-abstraction/acetylene-addition (HACA) mechanism still remain ambiguous. Here, by revealing the uncharted fundamental chemistry of the styrenyl (C₈ H₇ ) and the ortho-vinylphenyl radicals (C₈ H₇ )-key transient species of the HACA mechanism-with acetylene (C₂ H₂ ), we provide the first solid experimental evidence on the facile formation of naphthalene in a simulated combustion environment validating the previously postulated HACA mechanism for these two radicals. This study highlights, at the molecular level spanning combustion and astrochemistry, the importance of the HACA mechanism to the formation of the prototype PAH naphthalene.

An optically accessible pyrolysis microreactor

We report an optically accessible pyrolysis micro-reactor suitable for in situ laser spectroscopic measurements. A radiative heating design allows for completely unobstructed views of the micro-reactor along two axes. The maximum temperature demonstrated here is only 1300 K (as opposed to 1700 K for the usual SiC micro-reactor) because of the melting point of fused silica, but alternative transparent materials will allow for higher temperatures. Laser induced fluorescence measurements on nitric oxide are presented as a proof of principle for spectroscopic characterization of pyrolysis conditions.

Absolute photoionization cross sections of furanic fuels: 2-ethylfuran, 2-acetylfuran and furfural

Absolute photoionization cross sections of the molecules 2-ethylfuran, 2-acetylfuran and furfural, including partial ionization cross sections for the dissociative ionized fragments, are measured for the first time. These measurements are important because they allow fuel quantification via photoionization mass spectrometry and the development of quantitative kinetic modeling for the complex combustion of potential fuels. The experiments are carried out using synchrotron photoionization mass spectrometry with an orthogonal time-of-flight spectrometer used for mass analysis at the Advanced Light Source of Lawrence Berkeley National Laboratory. The CBS-QB3 calculations of adiabatic ionization energies and appearance energies agree well with the experimental results. Several bond dissociation energies are also derived and presented.

Pyrolysis Pathways of the Furanic Ether 2-Methoxyfuran

Substituted furans, including furanic ethers, derived from nonedible biomass have been proposed as second-generation biofuels. In order to use these molecules as fuels, it is important to understand how they break apart thermally. In this work, a series of experiments were conducted to study the unimolecular and low-pressure bimolecular decomposition mechanisms of the smallest furanic ether, 2-methoxyfuran. Electronic structure (CBS-QB3) calculations indicate this substituted furan has an unusually weak O-CH3 bond, approximately 190 kJ mol(-1) (45 kcal mol(-1)); thus, the primary decomposition pathway is through bond scission resulting in CH3 and 2-furanyloxy (O-C4H3O) radicals. Final products from the ring opening of the furanyloxy radical include 2 CO, HC≡CH, and H. The decomposition of methoxyfuran is studied over a range of concentrations (0.0025-0.1%) in helium or argon in a heated silicon carbide (SiC) microtubular flow reactor (0.66-1 mm i.d., 2.5-3.5 cm long) with reactor wall temperatures from 300 to 1300 K. Inlet pressures to the reactor are 150-1500 Torr, and the gas mixture emerges as a skimmed molecular beam at a pressure of approximately 10 μTorr. Products formed at early pyrolysis times (100 μs) are detected by 118.2 nm (10.487 eV) photoionization mass spectrometry (PIMS), tunable synchrotron VUV PIMS, and matrix infrared absorption spectroscopy. Secondary products resulting from H or CH3 addition to the parent and reaction with 2-furanyloxy were also observed and include CH2═CH-CHO, CH3-CH═CH-CHO, CH3-CO-CH═CH2, and furanones; under the conditions in the reactor, we estimate these reactions contribute to at most 1-3% of total methoxyfuran decomposition. This work also includes observation and characterization of an allylic lactone radical, 2-furanyloxy (O-C4H3O), with the assignment of several intense vibrational bands in an Ar matrix, an estimate of the ionization threshold, and photoionization efficiency. A pressure-dependent kinetic mechanism is also developed to model the decomposition behavior of methoxyfuran and provide pathways for the minor bimolecular reaction channels that are observed experimentally.

Chirped-Pulse millimeter-Wave spectroscopy for dynamics and kinetics studies of pyrolysis reactions

A Chirped-Pulse millimeter-Wave (CPmmW) spectrometer is applied to the study of chemical reaction products that result from pyrolysis in a Chen nozzle heated to 1000-1800 K. Millimeter-wave rotational spectroscopy unambiguously determines, for each polar reaction product, the species, the conformers, relative concentrations, conversion percentage from precursor to each product, and, in some cases, vibrational state population distributions. A chirped-pulse spectrometer can, within the frequency range of a single chirp, sample spectral regions of up to ∼10 GHz and simultaneously detect many reaction products. Here we introduce a modification to the CPmmW technique in which multiple chirps of different spectral content are applied to a molecular beam pulse that contains the pyrolysis reaction products. This technique allows for controlled allocation of its sensitivity to specific molecular transitions and effectively doubles the bandwidth of the spectrometer. As an example, the pyrolysis reaction of ethyl nitrite, CH3CH2ONO, is studied, and CH3CHO, H2CO, and HNO products are simultaneously observed and quantified, exploiting the multi-chirp CPmmW technique. Rotational and vibrational temperatures of some product molecules are determined. Subsequent to supersonic expansion from the heated nozzle, acetaldehyde molecules display a rotational temperature of 4 ± 1 K. Vibrational temperatures are found to be controlled by the collisional cooling in the expansion, and to be both species- and vibrational mode-dependent. Rotational transitions of vibrationally excited formaldehyde in levels ν4, 2ν4, 3ν4, ν2, ν3, and ν6 are observed and effective vibrational temperatures for modes 2, 3, 4, and 6 are determined and discussed.

A Signature of Roaming Dynamics in the Thermal Decomposition of Ethyl Nitrite: Chirped-Pulse Rotational Spectroscopy and Kinetic Modeling

Chirped-pulse (CP) Fourier transform rotational spectroscopy is uniquely suited for near-universal quantitative detection and structural characterization of mixtures that contain multiple molecular and radical species. In this work, we employ CP spectroscopy to measure product branching and extract information about the reaction mechanism, guided by kinetic modeling. Pyrolysis of ethyl nitrite, CH3CH2ONO, is studied in a Chen type flash pyrolysis reactor at temperatures of 1000-1800 K. The branching between HNO, CH2O, and CH3CHO products is measured and compared to the kinetic models generated by the Reaction Mechanism Generator software. We find that roaming CH3CH2ONO → CH3CHO + HNO plays an important role in the thermal decomposition of ethyl nitrite, with its rate, at 1000 K, comparable to that of the radical elimination channel CH3CH2ONO → CH3CH2O + NO. HNO is a signature of roaming in this system.