Search NASA⌕ Search

SEARCH · Search NASA

Results for “PRESSURE DEPENDENCE”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Automated theoretical chemical kinetics: Predicting the kinetics for the initial stages of pyrolysis

Large scale implementation of high level computational theoretical chemical kinetics offers the prospect for dramatically improving the fidelity of combustion chemical modeling. To facilitate such efforts, we have developed a suite of codes, collectively referred to as AutoMech, that allow for the automatic prediction of the kinetics for large sets of reactions via ab initio transition-state-theory based master-equation calculations. The primary input is simply the mechanism, a dictionary relating chemically identifiable species descriptors (e.g., SMILES or InChIs) to species labels in the mechanism, and a specification of the electronic structure and transition state theory models to be implemented. Here we illustrate the current utility of AutoMech through a study of the initial stages of pyrolysis for 3 sets of fuels: sequences of alkanes, alcohols, and aldehydes. For simplicity, the analysis focuses on abstractions from the fuel by H, CH 3 , and OH, and the decomposition of the resulting radicals. Altogether, there are a total of 166 input channels in these sets (more than 363 forward reactions when expanded to the full set of elementary reactions). The code successfully produces high quality rate estimates (with apparent uncertainties less than a factor of two in limited comparisons with experiment) for > 95% of these. For the radical decomposition reactions, the analysis includes predictions for the pressure dependence of the kinetics. This wide-ranging exploration illustrates (i) the effect of different levels of prediction on the expected accuracy, (ii) the branching between abstractions at different sites for different abstractors, (iii) the dependence of the rates on the chemical structure, and (iv) the variation in radical stabilities across chemical families. These results, as well as the demonstrated feasibility of the methodology, should find further utility in the development of accurate rate expressions for arbitrary fuels. (c) 2020 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

Elliott, Sarah N.↗
DOE OSTI2020

Formation of phenanthrenyl radicals via the reaction of acenaphthyl with acetylene

Ab initio G3(MP2,CC)/B3LYP/6–31G** calculations of the potential energy surface for the reaction of 3-acenaphthyl radical with acetylene were combined with Rice-Ramsperger-Kassel-Marcus Master Equation calculations of the temperature- and pressure-dependent reaction rate constants and product branching ratios with the goal to examine the transformation of an edge five-membered ring into a six-membered one. The reaction mechanism has been shown to contain two stages: (1) the formation of 3-ethynylacenaphthylene + H in the 3-acenaphthyl + C 2 H 2 reaction via acetylene addition to the radical site and an H loss, (2) the formation of the 1-phenanthrenyl radical in the 3-ethynylacenaphthylene + H reaction via H addition to the five-membered ring followed by the ring flip. The results of kinetic modeling show this pathway being faster than the pathway proposed earlier featuring the CCH chemisorption on the six-membered ring without involvement of the five-membered ring. A kinetic test demonstrated the dominance of the former pathway for the formation of pyrene from naphthalene. Here, it is expected that the new mechanism plays a major role in the transformation of five-membered rings on edges of PAH structures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Probing the antiknock effect of anisole through an ignition, speciation and modeling study of its blends with isooctane

In order to unravel the reaction pathways relevant to anisole co-oxidation within a fuel blend, a detailed study of isooctane/anisole blends was performed with the ULille RCM. Ignition delays as well as mole fraction profiles were measured during a two-stage ignition delay using sampling and GC techniques. These results are used to validate a kinetic model developed from ab initio calculations for the most relevant rate constants which included H-atom abstraction reactions from anisole, and reactions on the potential energy surfaces of methoxyphenyl + O 2 and anisyl + O 2 . Pressure dependent rate constants were computed for the methoxyphenyl + O 2 and anisyl + O 2 reactive systems using master equation code analysis. The new kinetic model shows good agreement with the experimental data. Dual brute-force sensitivity analysis was performed, on both first- and second-stages of ignition, allowing the identification of the most important reactions in the prediction of both ignition delays. It was observed that while pure anisole does not show NTC behavior, a 60/40 isooctane/anisole blend displays such behavior, as well as two-stage ignition. This suggests anisole addition may not be as beneficial to knock resistance as expected from its high octane number. The kinetic modeling results demonstrate the importance of H-abstraction reactions both from the methoxy group and from the aryl ring in ortho-position and the addition of the resultant radicals to O2, mostly leading to the formation of polar or non-aromatic products.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

From dynamics to kinetics for the dissociation of •QOOH radicals derived from isopentane

Thermal dissociation rates for hydroperoxyalkyl (•QOOH) radicals are challenging to measure due to their inherent instability. Recent dynamics experiments have demonstrated that in situ synthesis of •QOOH radicals in a supersonic expansion in combination with IR action spectroscopy provides an effective route to accurately map their energy-resolved dissociation rates. Collaborative theoretical analyses provide a route for converting from those microcanonical rates to accurate temperature- and pressure-dependent thermal rate constants. Here, we present the results of this conversion for the dissociation of two isopentyl β-QOOH isomers of relevance to isopentane oxidation. The results are compared with existing literature values, as well as with results obtained in a similar fashion for related •QOOH systems. Extrapolations from benchmark calculations of barrier properties for smaller systems to accurate predictions for larger systems and multidimensional hindered rotor treatments are both shown to be important components of the theoretical analysis. For completeness, we extend the theoretical analysis to a treatment of the full RO 2 /QOOH system for isopentyl (1,1-dimethyl-propyl) radical.

Ab initio kinetics↗
DOE OSTI2020

Synthesis and characterization of Pt(Cu 0.67 Sn 0.33 )

Pt(Cu 0.67 Sn 0.33 ) has recently been found in a natural sample. In order to be able to characterize this new ternary compound, we synthesized it from the elements. Samples were characterized by X-ray powder diffraction, differential scanning calorimetry, thermal relaxation calorimetry, and scanning electron microscopy studies. Density functional theory-based model calculations complemented the experimental studies. Pt(Cu 0.67 Sn 0.33 )was already formed at a relatively low temperature of 773 K. Rietveld refinement of Pt(Cu 0.67 Sn 0.33 ) has been carried out in CuAu-type or L1 0 -type structure, space group P4/mmm, with Pt on 0,0,0 and disordered Cu and Sn on 1/2, 1/2, 1/2 and Z = 1. The lattice parameters are a = 2.823(1) Å, c = 3.64(1) Å, and V = 29.00(4) Å which are in good agreement with values obtained earlier on the natural sample and with the results of DFT calculations. The vibrational entropy for Pt(Cu 0.67 Sn 0.33 ) is $S_{298.15}^{vib}$ = 79.9(7) J mol -1 K -1 . The pressure dependence up to 36(2) GPa of the unit-cell volume and the lattice parameters and unit-cell volume have been obtained by synchrotron based powder diffraction using a diamond anvil cell. A fit of a 3rd-order Birch–Murnaghan equation of state to the Pt(Cu 0.67 Sn 0.33 )) (p,V)-data results in a bulk modulus of B 0 = 215(27) GPa and B' = 5(2).

36 MATERIALS SCIENCE↗
DOE OSTI2011

Influence of nitrogen background pressure on structure of niobium nitride films grown by pulsed laser deposition

Depositions of niobium nitride thin films on Nb using pulsed laser deposition (PLD) with different nitrogen background pressures (10.7 to 66.7 Pa) have been performed. The effect of nitrogen pressure on NbN formation in this process was examined. The deposited films were characterized by X-ray diffraction (XRD), scanning electron microscope (SEM), atomic force microscope (AFM), and energy dispersive X-ray (EDX) analysis. Hexagonal {beta}-Nb{sub 2}N and cubic {delta}-NbN phases resulted when growth was performed in low nitrogen background pressures. With an increase in nitrogen pressure, NbN films grew in single hexagonal {beta}-Nb{sub 2}N phase. The formation of the hexagonal texture during the film growth was studied. The c/a ratio of the hexagonal {beta}-Nb{sub 2}N unit cell parameter increases with increasing nitrogen pressure. Furthermore, the N:Nb ratio has a strong influence on the lattice parameter of the {delta}-NbN, where the highest value was achieved for this ratio was 1.19. It was found that increasing nitrogen background pressure leads to change in the phase structure of the NbN film. With increasing nitrogen pressure, the film structure changes from hexagonal to a mixed phase and then back to a hexagonal phase.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Stochastic modal velocity field in rough-wall turbulence

Stochastically generated instantaneous velocity profiles are used to reproduce the outer region of rough-wall turbulent boundary layers in a range of Reynolds numbers extending from the wind tunnel to field conditions. Each profile consists in a sequence of steps, defined by the modal velocities and representing uniform momentum zones (UMZs), separated by velocity jumps representing the internal shear layers. Height-dependent UMZ is described by a minimal set of attributes: thickness, mid-height elevation, and streamwise (modal) and vertical velocities. These are informed by experimental observations and reproducing the statistical behaviour of rough-wall turbulence and attached eddy scaling, consistent with the corresponding experimental datasets. Sets of independently generated profiles are reorganized in the streamwise direction to form a spatially consistent modal velocity field, starting from any randomly selected profile. The operation allows one to stretch or compress the velocity field in space, increases the size of the domain and adjusts the size of the largest emerging structures to the Reynolds number of the simulated flow. By imposing the autocorrelation function of the modal velocity field to be anchored on the experimental measurements, we obtain a physically based spatial resolution, which is employed in the computation of the velocity spectrum, and second-order structure functions. The results reproduce the Kolmogorov inertial range extending from the UMZ and their attached-eddy vertical organization to the very-large-scale motions (VLSMs) introduced with the reordering process. The dynamic role of VLSM is confirmed in the –u'w' co-spectra and in their vertical derivative, representing a scale-dependent pressure gradient contribution.

42 ENGINEERING↗
DOE OSTI2022

Predicting Oxygen Off-Stoichiometry and Hydrogen Incorporation in Complex Perovskite Oxides

Chemically and structurally complex solid compounds, including those with significant off-stoichiometry, are rapidly extending new material functionality across a variety of applications. Accelerated development of these compounds requires accurate predictions of material defect properties including effective defect formation energies and equilibrium defect concentrations. Traditional first-principles approaches typically examine dilute defect concentrations and relatively ordered atomic structures to identify the lowest energy defect sites. These approaches are rarely suitable for describing the disorder present in these systems and its influence on defect formation, which can lead to unphysically large predictions for defect concentrations. Here, we demonstrate a new method to accurately predict the temperature and pressure dependence of oxygen vacancy concentrations and proton interstitial concentrations in complex oxides. This method extends standard dilute defect calculations to incorporate atomic and magnetic disorder, employs the ensemble descriptions of defect sites resulting in improved predictions of defect formation energies, and accounts for effects beyond the dilute defect limit. To demonstrate our method, we show that the predicted defect concentrations in perovskites used as ceramic fuel cell cathodes, including Ba0.5Sr0.5Fe0.8Zn0.2O3-d, Ba0.5Sr0.5Co0.8Fe0.2O3-d, and BaCo1-x-y-zFexZryYzO3-d, are in good agreement with experimental values, thereby opening the door for predictive design of complex oxides by these applications.

DFT↗
DOE OSTI2016

Investigation of Iso-octane Ignition and Validation of a Multizone Modeling Method in an Ignition Quality Tester

An ignition quality tester was used to characterize the autoignition delay times of iso-octane. The experimental data were characterized between temperatures of 653 and 996 K, pressures of 1.0 and 1.5 MPa, and global equivalence ratios of 0.7 and 1.05. A clear negative temperature coefficient behavior was seen at both pressures in the experimental data. These data were used to characterize the effectiveness of three modeling methods: a single-zone homogeneous batch reactor, a multizone engine model, and a three-dimensional computational fluid dynamics (CFD) model. A detailed 874 species iso-octane ignition mechanism (Mehl, M.; Curran, H. J.; Pitz, W. J.; Westbrook, C. K.Chemical kinetic modeling of component mixtures relevant to gasoline. Proceedings of the European Combustion Meeting; Vienna, Austria, April 14-17, 2009) was reduced to 89 species for use in these models, and the predictions of the reduced mechanism were consistent with ignition delay times predicted by the detailed chemical mechanism across a broad range of temperatures, pressures, and equivalence ratios. The CFD model was also run without chemistry to characterize the extent of mixing of fuel and air in the chamber. The calculations predicted that the main part of the combustion chamber was fairly well-mixed at longer times (> ~30 ms), suggesting that the simpler models might be applicable in this quasi-homogeneous region. The multizone predictions, where the combustion chamber was divided into 20 zones of temperature and equivalence ratio, were quite close to the coupled CFD-kinetics results, but the calculation time was ~11 times faster than the coupled CFD-kinetics model. Although the coupled CFD-kinetics model captured the observed negative temperature coefficient behavior and pressure dependence, discrepancies remain between the predictions and the observed ignition time delays, suggesting improvements are still needed in the kinetic mechanism and/or the CFD model. This approach suggests a combined modeling approach, wherein the CFD calculations (without chemistry) can be used to examine the sensitivity of various model inputs to in-cylinder temperature and equivalence ratios. In conclusion, these values can be used as inputs to the multizone model to examine the impact on ignition delay. Additionally, the speed of the multizone model also makes it feasible to quickly test more detailed kinetic mechanisms for comparison to experimental data and sensitivity analysis.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

Identifying the Gate-Opening Mechanism in the Flexible Metal–Organic Framework UTSA-300

Atomic-level understanding of the gate-opening phenomenon in flexible porous materials is an important step toward learning how to control, design, and engineer them for applications such as the separation of gases from complex mixtures. Here, we report such mechanistic insight through an in-depth study of the pressure-induced gate-opening phenomenon in our earlier reported metal–organic framework (MOF) Zn(dps) 2 (SiF 6 ) (dps = 4,4'-dipyridylsulfide), also called UTSA-300, using isotherm and calorimetry measurements, in situ infrared spectroscopy, and ab initio simulations. UTSA-300 is shown to selectively adsorb acetylene (C 2 H 2 ) over ethylene (C 2 H 4 ) and ethane (C 2 H 6 ) and undergoes an abrupt gate-opening phenomenon, making this framework a highly selective gas separator of this complex mixture. The selective adsorption is confirmed by pressure-dependent in situ infrared spectroscopy, which, for the first time, shows the presence of multiple C 2 H 2 species with varying strengths of bonding. A rare energetic feature at the gate-opening condition of the flexible MOF is observed in our differential heat energies, directly measured by calorimetry, showcasing the importance of this tool in adsorption property exploration of flexible frameworks and offering an energetic benchmark for further energy-based fundamental studies. Based on the agreement of this feature with ab initio-based adsorption energies of C 2 H 2 in the closed-pore structure UTSA-300a (“a” refers to the activated form), this feature is assigned to the weakening of the H-bond C–H···F formed between C 2 H 2 and fluorine of the MOF. Furthermore, our analysis identifies the weakening of this H-bond, the expansion of the closed-pore MOF upon successive C 2 H 2 coadsorption until its volume is close to that of the open-pore MOF, and the spontaneous gate opening to energetically favor C 2 H 2 adsorption in the open-pore structure as crucial steps in the gate-opening mechanism in this system.

36 MATERIALS SCIENCE↗
DOE OSTI2024

New Nonreactive Force Field for Accurate Molecular Dynamics Simulations of TATB at Extreme Conditions

Insensitive high explosives based on TATB (1,3,5-triamino-2,4,6-trinitrobenzene) are needed in applications when safety is of paramount importance, but the basic material properties that give rise to its insensitivity are not fully understood. Molecular dynamics modeling using empirical force fields (FFs) has been the main route to characterize many complicated dynamical properties of TATB single crystal, but these FFs have not been comprehensively tested at extreme conditions typical of detonation. We collect a benchmark data set of (quasi)static TATB physical properties as determined by experiments and electronic structure calculations and apply this data set to validate four existing TATB FFs along with a new TATB FF that we develop here and denote as the CEA-LLNL-Missouri (CLM) FF. Benchmark data include vibrational spectra, the TATB crystal temperature–pressure–volume equation of state and lattice parameters, properties of TATB crystal polymorphs and transitions to the gaseous and liquid states, dimer energy landscapes, the pressure-dependent elastic tensor, and the energy landscape for inelastic deformation via sliding of TATB crystal layers. As a general assessment, we find that the two existing nonreactive FFs are more accurate in describing TATB’s physical properties compared to the two variants of the ReaxFF reactive FF considered. The new CLM FF is found to consistently yield similar or better agreement with experiments and electronic structure theory than any of the existing FF models, and it presents a distinct improvement in accurately modeling TATB elasticity and equation of state. So this work is expected to help improve the accuracy of FF-based modeling of complicated dynamic responses that ultimately govern the safety and performance characteristics of this material.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Permutationally Invariant Polynomial Expansions with Unrestricted Complexity

A general strategy is presented for constructing and validating permutationally invariant polynomial (PIP) expansions for chemical systems of any stoichiometry. Demonstrations are made for three categories of gas-phase dynamics and kinetics: collisional energy-transfer trajectories for predicting pressure-dependent kinetics, three-body collisions for describing transient van der Waals adducts relevant to atmospheric chemistry, and nonthermal reactivity via quasiclassical trajectories. In total, 30 systems are considered with up to 15 atoms and 39 degrees of freedom. Permutational invariance is enforced in PIP expansions with as many as 13 million terms and 13 permutationally distinct atom types by taking advantage of petascale computational resources. The quality of the PIP expansions is demonstrated through the systematic convergence of in-sample and out-of-sample errors with respect to both the number of training data and the order of the expansion, and these errors are shown to predict errors in the dynamics for both reactive and nonreactive applications. Here, the parallelized code distributed as part of this work enables the automation of PIP generation for complex systems with multiple channels and flexible user-defined symmetry constraints and for automatically removing unphysical unconnected terms from the basis set expansions, all of which are required for simulating complex reactive systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Formation of Phenanthrene via Recombination of Indenyl and Cyclopentadienyl Radicals: A Theoretical Study

This work presents quantum chemical G3(MP2,CC)//B2PLYPD3/6-311G(d,p) calculations of the potential energy surface for the indenyl (C 9 H 7 ) + cyclopentadienyl (C 5 H 5 ) reaction followed by unimolecular decomposition of the C 14 H 11 radicals formed as the primary products, as well as the Rice–Ramsperger–Kassel–Marcus master equation (RRKM-ME) calculations to predict temperature- and pressure-dependent reaction rate constants and product branching ratios. The reaction begins with the barrierless recombination of indenyl and cyclopentadienyl forming a C 14 H 12 molecule with a new C–C bond connecting two five-membered rings, which subsequently dissociates to C 14 H 11 radicals by H losses. The primary products of the C 9 H 7 + C 5 H 5 → C 14 H 11 + H reaction can directly decompose by another H loss to benzofulvalene, and this pathway is most favorable in terms of the entropy factor and hence is preferable at higher temperatures. Otherwise, the initial C 14 H 11 isomers can undergo significant structural rearrangements before eliminating an H atom and producing phenanthrene, anthracene, or benzoazulenes, among which the formation of phenanthrene via the “spiran” pathway is clearly preferred. The calculated barriers along the computed favorable dissociation pathways are relatively low, in the ~30–40 kcal/mol range, making the C 14 H 11 radicals unstable at temperatures above 1000–1250 K at 1 atm. The results of RRKM-ME calculations show that, under typical combustion conditions, the decomposition of the C 14 H 11 radicals predominantly leads to benzofulvalene. However, the latter can be rapidly converted to phenanthrene via H-assisted isomerization with the rate constant for the benzofulvalene + H → phenanthrene + H reaction being close to 10 –11 cm 3 molecule –1 s –1 at 1000–1500 K and 1 atm. Here, the results provide further support for the hypothesis that recombination of two π radicals containing five-membered rings can lead to a growth of PAH with the formation of two fused six-membered rings, but the reaction mechanism may not be direct and is likely to involve two consecutive H atom losses leading to a fulvalene-like product, with subsequent H-assisted isomerization of the latter to a benzenoid PAH.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Theoretical Kinetics Predictions for Reactions on the NH 2 O Potential Energy Surface

Recent modeling studies of ammonia oxidation, which are motivated by the prospective role of ammonia as a zero-carbon fuel, have indicated significant discrepancies among the existing literature mechanisms. In this study, high-level theoretical kinetics predictions have been obtained for reactions on the NH 2 O potential energy surface, including the NH 2 + O, HNO + H, and NH + OH reactions. These reactions have previously been highlighted as important reactions in NH 3 oxidation with high sensitivity and high uncertainty. The potential energy surface is explored with coupled cluster calculations, including large basis sets and high-level corrections to yield high-accuracy (~0.2 kcal/mol 2σ uncertainty) estimates of the stationary point energies. Variational transition state theory is used to predict the microcanonical rate constants, which are then incorporated in master equation treatments of the temperature- and pressure-dependent kinetics. For radical–radical channels, the microcanonical rates are obtained from variable reaction coordinate transition state theory implementing directly evaluated multireference electronic energies. The analysis yields predictions for the total rate constants as well as the branching ratios. We find that the NO + H 2 channel contributes 10% of the total NH 2 + O flux at combustion temperatures, although this channel is not included in modern NH 3 oxidation mechanisms. Modeling is used to illustrate the ramifications of these rate predictions on the kinetics of NH 3 oxidation and NO x formation. Finally, the present results for NH 2 + O are important for predicting the chain branching and formation of NO in the oxidation of NH 3 and thermal DeNO x .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Association Kinetics for Perfluorinated n -Alkyl Radicals

Radical-radical reaction channels are important in the pyrolysis and oxidation chemistry of perfluoroalkyl substances (PFAS). In particular, unimolecular dissociation reactions within unbranched n-perfluoroalkyl chains, and their corresponding reverse barrierless association reactions, are expected to be significant contributors to the gas-phase thermal decomposition of families of species such as perfluorinated carboxylic acids and perfluorinated sulfonic acids. Unfortunately, experimental data for these reactions are scarce and uncertain. Furthermore, obtaining reliable theoretical predictions for such reactions is a laborious and computationally intensive task. Here, in this work, the chemical kinetics of the various association/decomposition reactions producing/decomposing the C 2 -C 4 series of unbranched n-perfluoroalkanes (C 2 F 6 , C 3 F 8 , and C 4 F 10 ) are examined using state-of-the-art ab initio transition-state-theory-based master-equation calculations. The variable-reaction-coordinate transition-state theory (VRC-TST) formalism is employed in computing the microcanonical and canonical rates for the association reactions. Reaction thermochemistry is obtained via composite quantum chemistry calculations and the laddering of error-canceling reaction schemes via a connectivity-based hierarchy approach employing ANL1/ANL0-style reference energies. Lennard-Jones collision model parameters for the considered systems were estimated by a direct dynamics approach, and collisional energy transfer parameters were obtained from analogies to systems of similar size and heavy-atom connectivity. A one-dimensional master equation approach was used to convert the microcanonical rate coefficients from the VRC-TST analysis into temperature- and pressure-dependent rate constants for the association reactions and the reverse dissociation reactions. The data are reported in standardized formats for usage in comprehensive chemical kinetic models for PFAS thermal destruction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

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 - CH 3 bond, approximately 190 kJ mol -1 (45 kcal mol -1 ); thus, the primary decomposition pathway is through bond scission resulting in CH 3 and 2-furanyloxy (O - C 4 H 3 O) radicals. Final products from the ring opening of the furanyloxy radical include 2 CO, HC$\equiv$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 CH 3 addition to the parent and reaction with 2-furanyloxy were also observed and include CH 2 =CH-CHO, CH 3 -CH=CH-CHO, CH 3 =CO-CH=CH 2 , 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-C 4 H 3 O), 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.

09 BIOMASS FUELS↗
DOE OSTI2015

Ab Initio Unimolecular Reaction Kinetics of CH 2 C(＝)OCH 3 and CH 3 C(＝O)OCH 2 Radicals

The unimolecular dissociation and isomerization kinetics of the methyl ester combustion intermediates methoxycarbonylmethyl (CH 2 C(＝O)OCH 3 ) and acetyloxylmethyl (CH 3 C(＝O)OCH 2 ) are theoretically investigated using high-level ab initio methods and the Rice–Ramsperger–Kassel–Marcus (RRKM)/master equation (ME) theory. Potential energy surfaces (PESs) are obtained using coupled cluster singles and doubles with perturbative triples correction (CCSD(T)), multireference singles and doubles configuration interaction (MRSDCI) with the Davidson–Silver (DS) correction, and multireference averaged coupled pair functional (MRACPF2) theory. The transition states exhibit high T1 diagnostics in coupled cluster calculations, suggesting the need for a multireference correlated wave function treatment. MRSDCI+DS and MRACPF2 capture their multiconfigurational character well, yielding lower barrier heights than CCSD(T) for these reactions. Additionally, the rate coefficients are computed using the RRKM/ME theory over a 500–2500 K temperature range and at a pressure range of 0.01 atm to the high-pressure limit. The temperature- and pressure-dependent rate coefficients are given in modified Arrhenius expressions. The β-scission of CH 2 C(＝O)OCH 3 is predicted to have a much higher barrier than the corresponding isomerization reaction and the β-scission of CH 3 C(＝O)OCH 2 . Consequently, the rate coefficients for β-scission of CH 2 C(＝O)OCH 3 are the smallest among the three reactions and the isomerization followed by decomposition to CH 3 C(＝O) and HCHO is the dominant reaction pathway for CH 2 C(＝O)OCH 3 . Furthermore, both radicals CH 2 C(＝O)OCH 3 and CH 3 C(＝O)OCH 2 are predicted to mainly decompose to CH 3 C(＝O) + HCHO rather than to the bimolecular product CH 2 C(＝O) + CH 3 O. A newly developed MA combustion mechanism, using our theoretical rate coefficients for the MA-related reactions, predicts combustion properties in good agreement with available experimental data.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

Emission lifetimes of a fluorescent dye under shock compression

The emission lifetimes of rhodamine 6G (R6G), were measured under shock compression to 9.1 GPa, with the dual intent of better understanding molecular photophysics in extreme environments and assessing the usefulness of fluorescence lifetime microscopy to measure spatially-dependent pressure distributions in shocked microstructured media. R6G was studied as free dye dissolved in poly-methyl methacrylate (PMMA), or dye encapsulated in silica microparticles suspended in PMMA. Thin layers of these materials in impedance-matched geometries were subjected to planar single-stage shocks created by laser-driven flyer plates. A synchronized femtosecond laser excited the dye at selected times relative to flyer plate arrival and the emission lifetimes were measured with a streak camera. Lifetimes decreased when shocks arrived. The lifetime decrease was attributed to a shock-induced enhancement of R6G nonradiative relaxation. At least part of the relaxation involved shock-enhanced intersystem crossing. For free dye in PMMA, the lifetime decrease during the shock was shown to be a linear function of shock pressure from 0-9 GPa, with a slope of -0.22 ns·GPa -1 . Furthermore, the linear relationship makes it simple to convert lifetimes into pressures. Lifetime measurements in shocked microenvironments may be better than emission intensity measurements, since lifetimes are sensitive to the surrounding environment, but insensitive to intensity variations associated with the motion and optical properties of a dynamically changing structure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗