Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “potential energy surface”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

On the performance of composite schemes in determining equilibrium molecular structures

Determination of equilibrium molecular structures is an essential ingredient in predicting spectroscopic parameters that help in identifying molecular carriers of microwave transitions. Here, in this study, the performance of two different ab initio composite approaches for obtaining equilibrium structures, “energy scheme” and “geometry scheme,” is explored and compared to semi-experimental equilibrium structures. This study is performed for a set of 11 molecules which includes diatomics, linear triatomics, and a few non-linear molecules. The ab initio calculations were performed using three tiers of composite chemical recipes. The current results show that as the overall rigor of calculation is increased, the semi-experimental and the ab initio numbers agree to within 0.0003 Å for all molecules in the test set. The composite approach based on correcting the potential energy surface (energy scheme) and the one based on correcting the geometry directly (geometry scheme) show excellent agreement with each other. This work represents a step toward development of efficient and highly accurate procedures for computing ab initio equilibrium structures.

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

Ab initio exploration of low–lying electronic states of linear and bent MNX + (M = Ca, Sr, Ba, Ra; X = O, S, Se, Te, Po) and their origins

High-level multireference and coupled cluster quantum calculations were employed to analyze low-lying electronic states of linear-MNX + and side-bonded-M[NX] + (M = Ca, Sr, Ba, Ra; X = O, S, Se, Te, Po) species. Their full potential energy curves (PECs), dissociation energies (D e s), geometric parameters, excitation energies (T e s), and harmonic vibrational frequencies (ω e s) are reported. The first three chemically bound electronic states of MNX + and M[NX] + are 3 Σ – , 1 Δ, 1 Σ + and 3 A", 1 A', 1 A", respectively. The 3 Σ – , 1 Δ, 1 Σ + of MNX + originate from the M + ( 2 D) + NX( 2 Π) fragments, whereas the 3 A", 1 A', 1A" states of M[NX] + dissociate to M + ( 2 S) + NX( 2 Π) as a result of avoided crossings. The MNX + and M[NX] + are real minima on the potential energy surface and their interconversions are possible. The M 2+ NX – /M 2+ [NX] – ionic structure is an accurate representation for their low-lying electronic states. The D e s of MNX + species were found to depend on the dipole moment (μ) of the corresponding NX ligands and a linear relationship between these two parameters was observed.

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

Reaction mechanisms of a cyclic ether intermediate: Ethyloxirane

Oxiranes are a class of cyclic ethers formed in abundance during low-temperature combustion of hydrocarbons and biofuels, either via chain-propagating steps that occur from unimolecular decomposition of β -hydroperoxyalkyl radicals ( β -˙QOOH) or from reactions of HOÓ with alkenes. Ethyloxirane is one of four alkyl-substituted cyclic ether isomers produced as an intermediate from n -butane oxidation. While rate coefficients for β -˙QOOH → ethyloxirane + ÓH are reported extensively, subsequent reaction mechanisms of the cyclic ether are not. As a result, chemical kinetics mechanisms commonly adopt simplified chemistry to describe ethyloxirane consumption by convoluting several elementary reactions into a single step, which may introduce mechanism truncation error—uncertainty derived from missing or incomplete chemistry. Here we provide fundamental insight on reaction mechanisms of ethyloxirane in support of ongoing efforts to minimize mechanism truncation error. Reaction mechanisms are inferred from the detection of products during chlorine atom-initiated oxidation experiments using multiplexed photoionization mass spectrometry conducted at 10 Torr and temperatures of 650 K and 800 K. To complement the experiments, calculations of stationary point energies were conducted using the ccCA-PS3 composite method on ˙R + O 2 potential energy surfaces for the four ethyloxiranyl radical isomers, which produced barrier heights for 24 reaction pathways. In addition to products from ˙QOOH → cyclic ether + ÓH and ˙R + O 2 → conjugate alkene + HOÓ, both of which were significant pathways and are prototypical to alkane oxidation, other species were identified from ring-opening of both ethyloxiranyl and ˙QOOH radicals. The latter occurs when the unpaired electron is localized on the ether group, causing the initial ˙QOOH structure to ring-open and form a resonance-stabilized ketohydroperoxide-type radical. The present work provides the first analysis of ethyloxirane oxidation chemistry, which reveals that consumption pathways are complex and may require an expansion of submechanisms to increase the fidelity of chemical kinetics mechanisms.

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

Automated reaction kinetics and network exploration (Arkane): A statistical mechanics, thermodynamics, transition state theory, and master equation software

The open-source statistical mechanics software described here, Arkane–Automated Reaction Kinetics and Network Exploration–facilitates computations of thermodynamic properties of chemical species, high-pressure limit reaction rate coefficients, and pressure-dependent rate coefficient over multi-well molecular potential energy surfaces (PES) including the effects of collisional energy transfer on phenomenological kinetics. Arkane can use estimates to fill in information for molecules or reactions where quantum chemistry information is missing. The software solves the internal energy master equation for complex unimolecular reaction systems. Inputs to the software include converged electronic structure computations performed by the user using a variety of supported software packages (Gaussian, Molpro, Orca, TeraChem, Q-Chem, Psi4). The software outputs high-pressure limit rate coefficients and pressure-dependent phenomenological rate coefficients, as well as computed thermodynamic properties (enthalpy, entropy, and constant pressure heat capacity) with added energy corrections. Some of the key features of Arkane include treatment of 1D, 2D or ND hindered internal rotation modes, treatment of free internal rotation modes, quantum tunneling effect consideration, transition state theory (TST) and Rice-Ramsperger-Kassel-Marcus (RRKM) rate coefficient computations, master equation solution with four implemented methods, inverse-Laplace transform of high-pressure limit rate coefficients into the energy domain, energy corrections based on bond-additivity or isodesmic reactions, automated and efficient PES exploration, and PES sensitivity analysis. The present work describes the design of Arkane, how it should be used, and refers to the theory that it employs. Arkane is distributed via the RMG-Py software suite (https://github.com/ReactionMechanismGenerator/RMG-Py).

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

The role of precursor decomposition in the formation of samarium doped ceria nanoparticles via solid-state microwave synthesis

The impact on the final morphology of ceria (CeO 2 ) nanoparticles made from different precursors (commercial: cerium acetate/nitrate) and in house: cerium tri(methylsilyl)amide (Ce-TMSA)) via a microwave solid state reaction has been determined. In all instances, powder X-ray diffraction indicated that the cubic fluorite CeO 2 phase (PDF# 04–004-9150, with the space group Fm-3 m) had formed. Scanning electron microscopy (SEM) images revealed spherical nanoparticles were produced from the Ce-TMSA precursor. The commercial acetate and nitrate precursors produced particles with irregular morphology. The roles of the precursor decomposition and binding energy in the synthesis of the nanocrystals with various morphologies, as well as a possible growth mechanism, were evaluated based on experimental and computational data. The formation of spherical shaped nanoparticles was determined to be due to the preferential single-step decomposition of the Ce-TMSA as well as the low activation energy to overcome decomposition. Due to the complicated decomposition of the commercial precursors and high activation energy the resulting particles adopted an irregular morphology. Highly uniform samarium doped ceria (Sm x Ce 1-x O 2-δ ) nanospheres were also synthesized from Ce-TMSA and samarium tri(methylsilyl)amide (Sm-TMSA). The effects of reaction time and temperature, on the final morphology were observed through SEM. The rapid single-step decomposition of TMSA-based precursors as observed through thermogravimetric analysis (TGA) and confirmed through the calculation of potential energy surfaces and binding energies from density functional theory (DFT) calculations, indicated that nanoparticle formation follows LaMer’s classical nucleation theory.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Theoretical and kinetic modeling study of H 2 S pyrolysis

Hydrogen sulfide pyrolysis was investigated theoretically and through chemical kinetic modeling. Reactions on the SHH potential energy surface, primarily S + H 2 (+Ar) ⇌ H 2 S (+ Ar) (R1) and S + H 2 ⇌ SH + H (R6b) were characterized by ab initio calculations. Results for k 1 were in good agreement with experiment, but deviated strongly below 2000 K from values previously used in modeling. Collider efficiencies for H 2 S, S 2 , and N 2 compared to Ar were calculated for R1. Hydrogen sulfide decomposition experiments reported in literature were re-examined in terms of an updated detailed chemical kinetic model. Concentration profiles for the atomic S at high temperature in shock tubes supported the present value of k 1 and served to constrain the rate constants for reaction of S with SH and H 2 S. To explain results from batch and flow reactors, conducted at high H 2 S concentrations in the 900–1400 K range, a very fast rate constant was required for HSS + H ⇌ SH + SH. Under dilute conditions, the gas-phase chemistry was too slow to compete and the decomposition of H 2 S was controlled by loss on the reactor surface.

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

Theoretical and kinetic modeling study of hydrazine oxidation

The present work constitutes the first theoretical and kinetic modeling study of hydrazine oxidation, which may be important for burnout in ammonia-fueled combustion. The kinetics of the oxidation of N 2 H 4 , N 2 H 3 and tHNNH by molecular oxygen were investigated via a quantum chemistry/canonical transition state theory approach. Geometries and anharmonic frequencies were obtained with density functional theory, and energies from coupled cluster calculations (CCSD(T)) extrapolated to the infinite basis set limit, with corrections for core-valence electron correlation, scalar relativistic effects, and higher level correlation up to lambda coupled cluster, CCSDT(Q) Λ . The key reactions occurred on the N 2 H 4 O 2 potential energy surface, where the results indicated a fast reaction of N 2 H 3 with HO 2 via singlet adducts to yield tHNNH + H 2 O 2 and HNN(H)O + H 2 O, while reaction on the triplet surface proceeds via a bound complex followed by a tight, submerged barrier to yield N 2 H 4 + O 2 . The results were incorporated in a detailed reaction mechanism, which was used to interpret the shock tube results from Michel and Wagner (1965) on the effect of O 2 on hydrazine conversion at 1100–1400 K. The kinetic model captured qualitatively the observed behavior, but underestimated the reaction rate under oxidizing conditions. The hydrazine pyrolysis chemistry dominated conversion at reducing conditions and/or high temperature. At oxidizing conditions and intermediate temperatures (≲ 1400 K), reactions of N 2 -amines with HO 2 and O 2 were important for the oxidation rate.

Ab initio calculations↗
DOE OSTI2021

Theoretical kinetics predictions for NH 2 + HO 2

Recent modeling studies of NH 3 oxidation, which are motivated by the prospective role of ammonia as a zero-carbon fuel, have indicated significant discrepancies between existing literature mechanisms. In this study high level theoretical kinetics predictions have been obtained for the reaction of NH 2 with HO 2 , which has previously been highlighted as an important reaction 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) 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 the 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 constant as well as the branching to the NH 3 + O 2 , H 2 NO + OH, and HNO + H 2 O channels. Rate constants are also reported for the H 2 NO + OH reaction as they arise naturally from the analysis. The rate constant and branching fraction determined in this work for the NH 2 + HO 2 reaction deviate significantly from values used in most previous modeling studies. The fact that the main product channel is chain terminating, rather than propagating, has strong implications for modeling NH 3 ignition and oxidation, in particular at intermediate temperatures and elevated pressure.

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

Molecular weight growth by the phenyl + cyclopentadienyl reaction: Well-skipping, ring-opening, and dissociation

Radical-radical reactions contribute to molecular weight growth of polycyclic aromatic hydrocarbons (PAHs), but their role is not well understood. The phenyl (Ph, C 6 H 5 ) + cyclopentadienyl (C 5 H 5 ) reaction has not been studied but can serve as a surrogate reaction for many reactions of larger aryl-σ radicals with resonance-stabilized π-radicals. The adducts of these reactions easily lose an H atom to produce a new π-radical. Here, through a combination of dilute flash pyrolysis experiments, theory, and simulation, we find that for the title reaction, substantial Ph—C 5 H 4 + H is formed, dominantly by a well-skipping pathway, at conditions around 30 Torr and 1400 K. This radical is thermally stable at moderate temperatures but decomposes at higher temperatures following ring opening of the 5-member ring. Products corresponding to the loss of H, C 2 H 2 , and C 3 H 3 are observed experimentally. Calculation of the C 11 H 9 potential energy surface affirms that these are expected decomposition products and reveals multiple possible pathways.

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

The role of energy transfer and competing bimolecular reactions in characterizing the unimolecular dissociations of allylic radicals

Resonantly stabilized radicals (RSR's) such as allylic radicals typically have bond energies exceeding 45 kcal/mol and are therefore reasonably long-lived up to quite high temperatures. Consequently, kinetics of the unimolecular dissociations and bimolecular reactions of these long-lived RSR's assume significance in high temperature chemically reacting systems such as in flames. Our recent analyses (J. Cho et al., Proc. Comb. Inst., 2022) of the uncertainty propagated by the energy transfer parameters for 1-methylallyl (1MA) dissociation to the laminar flame speed, indicates an intricate coupling between the kinetics of 1MA dissociation (chain propagation) and its reaction with H-atoms (chain-termination). Here, this work extends such analyses to other C 3 – C 5 allylic radicals: allyl, 2-methylallyl, 1,1-dimethylallyl, 1,2-dimethylallyl, and 1,3-dimethylallyl. Theoretical kinetics for the C 3 – C 4 allylic radicals were taken from prior literature studies (J.A. Miller et al., J. Phys. Chem., 2008; J. Cho et al., Proc. Comb. Inst., 2022; R.S. Tranter et al., Proc. Comb. Inst., 2017). Potential energy surfaces were newly computed for the dissociations of the larger C5 allylic radicals and the reactions of these radicals with H-atoms. These PES's were then used to calculate rate constants for the relevant unimolecular and bimolecular reactions. Energy transfer has a crucial role to play in these pressure-dependent reactions and therefore the present kinetics calculations relied on a-priori predictions of collisional energy transfer parameters to best characterize the branching between chain propagating unimolecular dissociations of these allylic radicals and chain terminating bimolecular reactions of these RSR's with H-atoms. Prompt dissociations are now established to be a universal feature of free radical dissociations. Here we have also characterized the propensity for prompt dissociations in these more-stable allylic RSR's and their role in combustion simulations.

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

A theoretical kinetic study of ĊH 3 + ṄH 2 : From electronic structure to NH 3 /CH 4 combustion modelling implications

Carbon–nitrogen interaction reactions play an important role in governing the reactivity of ammonia blended fuels. However, there remains uncertainties regarding their detailed reaction pathways and rate constants, hampering the development of high-fidelity chemical kinetic models. In this study, the kinetics of ĊH 3 + ṄH 2 , a key C–N interaction reaction in ammonia/methane blend combustion have been investigated. The potential energy surface has been explored using the high-level ANL0F method, yielding highly accurate stationary point energies that agree with ATcT values within 0.1 kcal mol –1 . Variable reaction coordinate transition state theory is used to treat the barrierless association and decomposition reaction channels, based on directly sampled radical-radical interaction energies at the CASPT2-F12(2e,2o)/cc-pVTZ-F12 level of theory. The minimum transitional mode numbers of states obtained are then coupled with the RRKM/master equation to calculate temperature- and pressure-dependent rate constants. Our a priori calculations capture available experimental measurements from the literature very well. The calculated rate constants have been incorporated into an NH 3 /CH 4 chemical kinetic model currently under development at the University of Galway. The effect of the updated kinetic data for ĊH 3 + ṄH 2 on model predicted NH 3 /CH 4 fuel reactivity is elucidated.

ab initio↗
DOE OSTI2021

Saddle point sampling using scaled normal coordinates

The sampling of saddle points (SPs) on a potential energy surface (PES) is critical for describing the dynamics and transport properties of materials. In this paper we propose a saddle point search (SPS) method that enables efficient sampling of practically meaningful SPs around a local minimum by converting the PES to the scaled normal coordinates (SNCs), together with the dimer method. We find that the pathway leading to a SP using SNCs is highly efficient and is independent of the system size. This results in the significantly increased SPS efficiency and the enhanced possibility of a complete catalog of the SPs. In addition, we perform SPSs using the SNCs for the diffusion of a vacancy and a dumbbell in body-centered cubic Fe and demonstrate (1) the required number of system force evaluations for a SPS decreases by at least an order of magnitude, (2) the effects of the system size on the number of force evaluations are greatly decreased, and (3) the effects of the number of atoms involved in the SPS on the probability of finding practically meaningful SPs is mostly eliminated. We also discuss the computational cost of introducing the SNCs. It should be highlighted that the SNCs could also be applied to other minimum-mode following methods, demonstrating the general versatility of the proposed method.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Saddle point search with dynamic active volume

Sampling potential energy surface (PES) is critical for many problems in materials science, chemistry, physics, and biology and requires highly efficient saddle point searches (SPS). In the study, we introduce the concept of dynamic active volume (DAV) in addition to the active volume in self-evolving atomistic kinetic Monte Carlo (SEAKMC). The DAV method has further reduced the dimensionality of the PES at the elevation stage of a SPS. At the subsequent converging stage, the dynamic boundary in DAV is lifted to allow the system to converge to the right location of a saddle point. Coupled with the dimer method, the DAV method not only significantly reduces the time cost for a given search attempt, but also dramatically increases the probability of finding relevant saddle points for PES sampling. Further, a Python software package within the framework SEAKMC (SEAKMC_py) with the DAV method has been developed.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Pilgrim : A thermal rate constant calculator and a chemical kinetics simulator

Pilgrim is a program written in Python and designed to use direct dynamics in the calculation of thermal rate constants of chemical reactions by the variational transition state theory (VTST), based on electronic structure calculations for the potential energy surface. Pilgrim can also simulate reaction mechanisms using kinetic Monte Carlo (KMC). For reaction processes with many elementary steps, the rate constant of each of these steps can be calculated by means of conventional transition state theory (TST) or by using VTST. In the current version, Pilgrim can evaluate thermal rates using the canonical version of reaction-path VTST, which requires the calculation of the minimum energy path (MEP) associated with each elementary step or transition structure. Multi-dimensional quantum effects can be incorporated through the small-curvature tunneling (SCT) approximation. These methodologies are available both for reactions involving a single structure of the reactants and the transition state and also for reactions involving flexible molecules with multiple conformations of the reactant and/or of the transition state. For systems with many conformers, the program can evaluate each of the elementary reaction rate constants by multipath canonical VTST or multi-structural VTST. Moreover, the reactant can be unimolecular or bimolecular. Torsional anharmonicity can be incorporated through either the MSTor or the Q2DTor programs. Dual-level calculations are also available in Pilgrim: automatic high-level single-point energies can be used to correct the energy of reactants, transition states, products, and MEP points using the interpolated single-point energies (ISPE) algorithm. When the rate constants of all the chemical processes of interest are known, by means of their calculation using Pilgrim or alternatively through analytical fits to the rate constants as functions of temperature, it is possible to simulate a multistep mechanism under specified laboratory conditions using KMC. Finally, this algorithm allows performing a kinetic simulation to monitor the evolution of each chemical species with time and obtain the product yields.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Anion photoelectron spectroscopic and relativistic coupled-cluster studies of uranyl dichloride anion, UO 2 Cl 2 –

A joint relativistic coupled-cluster and experimental photoelectron (PE) spectroscopic study of the uranyl dichloride anion, UO 2 Cl 2 – , is reported. Sophisticated electronic-structure calculations predict the photodetachment of UO 2 Cl 2 – to involve a U 5f electron and to be followed by significant geometry relaxation. Therefore, the adiabatic electron affinity (EA a ) of the uranyl dichloride neutral molecule, UO 2 Cl 2 , and the vertical detachment energy (VDE) of its anion, UO 2 Cl 2 – , provide valuable information about its uranium 5f orbital energies. The EA a value was computed to be 3.15 eV. The VDE value was calculated to be 3.55 eV by augmenting the computed EA a with a shift derived from a Franck–Condon simulation using coupled-cluster potential energy surfaces. Here, the VDE, which corresponds to the highest intensity peak in the PE spectrum, was measured to be 3.69 ± 0.20 eV, in good agreement with the computed value. The origin transition in the PE spectrum, whose electron binding energy corresponds to the EA a , was assigned to the feature at 3.2 ± 0.20 eV, consistent with the computed EA a .

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
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 OSTI2020

Insights on keto-hydroperoxide formation from O 2 addition to the beta-tetrahydrofuran radical

Cyclic ethers provide an interesting case study of low-temperature oxidation chemistry, especially in relevance to biofuels. A recent experimental study revealed new questions regarding the low-temperature oxidation mechanism of tetrahydrofuran (THF) concerning the formation of keto-hydroperoxides. In particular, keto-hydroperoxides originating from the THF-β radical were not captured accurately by current literature models, motivating this work to calculate the energetics of the first and second O 2 -addition pathways for THF radicals. Electronic structure calculations at the CCSD(T)/cc-pV∞Z//M06-2X/cc-pVTZ level of theory were used to generate potential energy surfaces for the α-C and β-C THF radicals and subsequent pathways to the formation of the keto-hydroperoxide isomers. These are the first theoretical calculations of the second O 2 -addition radical pathways for the THF-β radical. Results from the theorical work provided further insight into the low-temperature oxidation of THF. This included identifying the pathways most likely to form the keto-hydroperoxide isomers observed in prior experimental work; and detecting that the shortcomings in prior models are likely due to uncertainties in R + THF abstraction reaction rates. Here, these conclusions will motivate future work for accurate THF kinetic model development.

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↗