Search NASA⌕ Search

SEARCH · Search NASA

Results for “ab initio theory”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Graph theory approach to determine configurations of multidentate and high coverage adsorbates for heterogeneous catalysis

Abstract Heterogeneous catalysts constitute a crucial component of many industrial processes, and to gain an understanding of the atomic-scale features of such catalysts, ab initio density functional theory is widely employed. Recently, growing computational power has permitted the extension of such studies to complex reaction networks involving either high adsorbate coverages or multidentate adsorbates, which bind to the surface through multiple atoms. Describing all possible adsorbate configurations for such systems, however, is often not possible based on chemical intuition alone. To systematically treat such complexities, we present a generalized Python-based graph theory approach to convert atomic scale models into undirected graph representations. These representations, when combined with workflows such as evolutionary algorithms, can systematically generate high coverage adsorbate models and classify unique minimum energy multidentate adsorbate configurations for surfaces of low symmetry, including multi-elemental alloy surfaces, steps, and kinks. Two case studies are presented which demonstrate these capabilities; first, an analysis of a coverage-dependent phase diagram of absorbate NO on the Pt 3 Sn(111) terrace surface, and second, an investigation of adsorption energies, together with identifying unique minimum energy configurations, for the reaction intermediate propyne (CHCCH 3 *) adsorbed on a PdIn(021) step surface. The evolutionary algorithm approach reproduces high coverage configurations of NO on Pt 3 Sn(111) using only 15% of the number of simulations required for a brute force approach. Furthermore, the screening of potentially hundreds of multidentate adsorbates is shown to be possible without human intervention. The strategy presented is quite general and can be applied to a spectrum of complex atomic systems.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Shake loss intensities in x-ray photoelectron spectroscopy: Theory, experiment, and atomic composition accuracy for MgO and related compounds

The relative intensities of XPS core levels, scaled by their photoionization cross sections, are regularly used to determine sample atomic composition. Cross sections, however, give the intensity to all possible final states for the core ionizations, not just to the main peak. This includes all intrinsic satellite structure (shake states and, for open shell systems, the different ionic multiplets). In practice, for solids, this is usually experimentally impossible to determine accurately because such a satellite structure sits on the inelastically scattered electron background and cannot be easily separated. Therefore, usually, only the intensity of the main peak is used. This limits the ultimate possible accuracy of XPS composition determination. The purpose of the present paper is to examine the contributions that a theoretical analysis of losses of intensity can make to improve quantitation. For an MgO single crystal, we show that the correct stoichiometry of 1:1 can be recovered using the theoretical analysis of the experimental MgO peak ratio intensities. For materials with a sufficient bandgap for the XPS main peaks to be separated from the scattered background, the intensity of main peaks can often be accurately determined. Thus, if one uses theory to calculate that fraction of the total intensity lost from a main peak into all its satellite structure, the intensity of just main peaks could then be used to more accurately determine relative atom % composition. This work tests this approach using a single crystal MgO (50% Mg, 50% O) standard. Ab initio electronic structure theory of representative MgO clusters is used to determine Hartree–Fock wave functions for the ground state and final ionized states corresponding to the main Mg 2p and O1s XPS peaks of the oxide. The sudden approximation, SA, is used to determine the fractional losses from these main peaks to shake satellites, which is found to be greater for O1s than Mg2p. This results in predicted “apparent composition” for stoichiometric MgO of 55.2% Mg, 44.8% O instead of the true 50% Mg, 50% O. Equivalent theory for CaO results in a predicted apparent Ca value of 53.4%. Experimentally, using Mg2s or 2p intensity ratio to O1s, we find values between 52.2% and 56.0% Mg using two crystals and four different instrument electron pass energies. The average value of the measurements is 54.5% Mg when corrected for the presence of an adventitious carbon overlayer and slight surface hydroxide. Though this agreement with theory may be somewhat fortuitous, given the potential experimental errors, which are fully discussed, it is similar to that in our earlier study on LiF. We also present preliminary experimental data on Mg(OH) 2 and MgSO 4 , which show a similar trend of apparently higher than 50% Mg, but we have no theory values. We are not yet able to experimentally test for validation of the difference between apparent composition for MgO (55.2% Mg) and CaO (53.4% Ca), owing to significant carbonate formation at the surface of the single crystal CaO. Here, an important conclusion is that the theoretical determination of shake losses, obtained with ab initio wavefunctions and the SA, is likely to be a useful way to calibrate the accuracy and reliability of compositions obtained from XPS intensities and merits further study.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Spiers Memorial Lecture: Theory of unimolecular reactions

One hundred years ago, at an earlier Faraday Discussion meeting, Lindemann presented a mechanism that provides the foundation for contemplating the pressure dependence of unimolecular reactions. Since that time, our ability to model and predict the kinetics of such reactions has grown in leaps and bounds through the synergy of ever more sophisticated experimental studies with increasingly high level theoretical analyses. Herein, this review begins with a brief historical overview of the progress from the Lindemann mechanism to the master equation, which now provides the cornerstone for most theoretical analyses. The current status of ab initio transition state theory based implementations of the master equation is then reviewed, beginning with a discussion of the energy resolved chemical conversion rates, followed by a review of the pressure dependence of the thermal kinetics. The latter discussion focusses on the collisional energy and angular momentum transfer rates as well as the complexities of non-thermal effects and of multiple-well multiple-channel reactions. The synergy with recent state-of-the-art experiments is used to motivate these discussions. Attempts to automate the calculations are also briefly reviewed. Throughout the discussion, we provide our perspective on current understanding and continuing challenges and opportunities. One key challenge relates to the coupling of sequences of reactions, which frequently leads to deviations from the classic presumption of thermal reactants.

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

Construction of meta-GGA functionals through restoration of exact constraint adherence to regularized SCAN functionals

The strongly constrained and appropriately normed (SCAN) meta-GGA exchange–correlation functional is constructed as a chemical environment-determined interpolation between two separate energy densities: one describes single-orbital electron densities accurately and another describes slowly varying densities accurately. To conserve constraints known for the exact exchange–correlation functional, the derivatives of this interpolation vanish in the slowly varying limit. While theoretically convenient, this choice introduces numerical challenges that degrade the functional’s efficiency. We have recently reported a modification to the SCAN meta-GGA, termed restored-regularized-SCAN (r 2 SCAN) , that introduces two regularizations into SCAN, which improve its numerical performance at the expense of not recovering the fourth order term of the slowly varying density gradient expansion for exchange. Here, we show the derivation of a progression of density functional approximations [regularized SCAN (rSCAN), r++SCAN, r 2 SCAN, and r 4 SCAN] with increasing adherence to exact conditions while maintaining a smooth interpolation. The greater smoothness of r 2 SCAN seems to lead to better general accuracy than the additional exact constraint of SCAN or r 4 SCAN does.

, correlation energy↗
DOE OSTI2026

Cross-slip and easy-glide CRSS in titanium: Theoretical predictions and in-situ TEM measurements

This study investigates the mechanics of prismatic and first-order pyramidal $\langle$a$\rangle$ slip in titanium (Ti), elucidating the physics of easy-glide and cross-slip through a combination of theory and experiments. Screw-character prismatic (Pr) dislocations in Ti are of particular interest because their complex cores can be stable or unstable, leading to activation by either cross-slip or planar glide. To investigate these mechanisms, site-specific micro tensile samples were prepared using focused ion beam (FIB) milling and mounted on a push-to-pull (PTP) device for in-situ transmission electron microscopy (TEM) tensile testing. The in-situ experiments provide direct observations of the onset of dislocation motion and the precise determination of the critical resolved shear stress (CRSS) for the activated mechanisms, and its evolution with load cycling. A comprehensive theory has been developed to predict the CRSS values for easy glide, cross-slip, and multiplication of dislocations. Predicted critical stresses for pyramidal (π)-to-Pr and reverse cross-slip agree closely with the experimental measurements. The latter cross-slip stress is a factor of two higher than that of unobstructed planar slip. The model accounts for overlapping dislocation cores and employs a Wigner-Seitz based cell to evaluate misfit energies. By combining ab initio density functional theory (DFT) with anisotropic elasticity, the framework identifies minimum energy pathways for dislocation glide, which can be intermittent and zig-zag. A simplified expression utilizing (π) and Pr Schmid factor ratios is proposed for critical stress corresponding to (π)-to-Pr cross-slip transition. The results are strongly dependent on crystal orientation, underscoring non-Schmid behavior. Overall, this study explores key critical stress parameters essential for informing higher-scale simulations of plasticity in Ti.

Critical stress↗
DOE OSTI2024

Breaking the Million-Electron and 1 EFLOP/s Barriers: Biomolecular-Scale Ab Initio Molecular Dynamics Using MP2 Potentials

The accurate simulation of complex biochemical phenomena has historically been hampered by the computational requirements of high-fidelity molecular-modeling techniques. Quantum mechanical methods, such as ab initio wave-function (WF) theory, deliver the desired accuracy, but have impractical scaling for modeling biosystems with thousands of atoms. Combining molecular fragmentation with MP2 perturbation theory, this study presents an innovative approach that enables biomolecular-scale ab initio molecular dynamics (AIMD) simulations at WF theory level. Leveraging the resolution-of-the-identity approximation for Hartree-Fock and MP2 gradients, our approach eliminates computationally intensive four-center integrals and their gradients, while achieving near-peak performance on modern GPU architectures. The introduction of asynchronous time steps minimizes time step latency, overlapping computational phases and effectively mitigating load imbalances. Utilizing up to 9,400 nodes of Frontier and achieving 59% (1006.7 PFLOP/s) of its double-precision floating-point peak, our method enables us to break the million-electron and 1EFLOP/s barriers for AIMD simulations with quantum accuracy.

Kurzak, Jakub↗
DOE OSTI2021

Computational study of the electrostatic potential and charges of multivalent ionic liquid molecules

The electrostatic interactions in ionic liquids (ILs) can be difficult to quantify, but this information is very valuable for improving intermolecular potentials for performing molecular-level simulations. Although polarizable forcefields can provide improved accuracy, the computational efficiency of fixed charge molecular models is often preferred for larger systems. Here, we propose a multi-scale screening approach for analyzing the charge scaling behavior and other electrostatic properties of ILs, which is more attainable than previous analysis methods that rely on the IL crystalline structures. Based on isolated molecules and small clusters, we study the electrostatic potential properties and charges of multivalent IL molecules using density-functional theory calculations and ab initio molecular dynamics simulations. The charges of the IL molecules are calculated and compared using several different charge estimation techniques. Although there are some differences among the estimated partial charges on the atoms, the total anion and cation charges are very similar. In conclusion, the charge scaling factors are substantial, and they are only slightly affected by the IL cluster size model.

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

Combined density functional theory/kinetic Monte Carlo investigation of surface morphology during cycling of Li-Cu electrodes

In this work, we employed a combined density functional theory and kinetic Monte Carlo approach in order to investigate the nanoscale behavior of lithium stripping and plating on lithium-copper nanoalloys. Density functional theory calculations were used to develop an atom-by-atom understanding of the nucleation behavior of lithium on copper surfaces with different surface morphologies. We found that nanometric pits resulted in significantly more favorable lithium deposition when compared to pristine or lightly-defective Cu surfaces, which is in agreement with previous experimental studies. We then used density functional theory and constrained ab initio molecular dynamics to develop interaction and reaction parameters for a larger-scale kinetic Monte Carlo model in order to provide insight into the stripping and plating behavior of nanostructured Li-Cu nanoalloy anodes. Substantial losses in capacity and strong indications of morphological change were observed for a wide range of Cu:Li ratios over 100 cycles. We observed a complex relationship between the nanoparticle composition and the properties of the formed solid-electrolyte interphase.

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

Accelerated Modeling of Lithium Diffusion in Solid State Electrolytes using Artificial Neural Networks

Abstract Previous efforts to understand structure‐function relationships in high ionic conductivity materials for solid state batteries have predominantly relied on density functional theory (DFT‐) based ab initio molecular dynamics (MD). Such simulations, however, are computationally demanding and cannot be reasonably applied to large systems containing more than a hundred atoms. Here, an artificial neural network (ANN) is trained to accelerate the calculation of high accuracy atomic forces and energies used during such MD simulations. After carefully training a robust ANN for four and five element systems, nearly identical lithium ion diffusivities are obtained for Li 10 GeP 2 S 12 (LGPS) when benchmarking the ANN‐MD results with DFT‐MD. Applying the ANN‐MD approach, the effect of chlorine doping on the lithium diffusivity is calculated in an LGPS‐like structure and it is found that a dopant concentration of 1.3% maximizes ionic conductivity. The optimal concentration balances the competing consequences of effective atomic radii and dielectric constants on lithium diffusion and agrees with the experimental composition. Performing simulations at the resolution necessary to model experimentally relevant and optimal concentrations would be infeasible with traditional DFT‐MD. Systems that require a large number of simulated atoms can be studied more efficiently while maintaining high accuracy with the proposed ANN‐MD framework.

Rao, Karun K.↗
DOE OSTI2023

The α-WC(0001) Surface as a Hydrogen Sponge: A First Principle Study of H 2 Dissociation and Formation of Low and High Coverages

Tungsten carbide (WC) displays a Pt-like behavior in catalysis, applied in hydrogenation processes. Numerous theoretical studies have modeled the behavior and use of adsorbed hydrogen without obtaining a general picture, missing basic links between H 2 dissociation and generation of high surface coverage (Θ H >0.5 ML). Here, in this study, the capability of C- and W-terminations of the α-WC(0001) surface is analyzed to dissociate several H 2 molecules to produce coverages, Θ H , ranging from low to very high values (0.13<Θ H <2.00 ML). Density functional theory and an ab initio atomistic thermodynamic were used to achieve the conditions for H 2 dissociation. The WC–C surface has higher capacity to dissociate H 2 molecules than WC–W. However, both surfaces can reach full surface coverage, Θ H =1 ML, at mild ambient conditions, T=300 K and P=1 atm, and even up to 500 K at low and high pressures. The H-adatoms on WC–W are more labile than on WC–C. The binding of adsorbates is hindered at high Θ H , implying a need to modulate Θ H according to the application. The results give the basis to understand the capabilities of WC-based catalysts in hydrogenation-related reactions, with the advantage of WC being a hydrogen reservoir at mild practical catalytic conditions.

03 NATURAL GAS↗
DOE OSTI2021

Effect of Ligand Polarity on the Internal Dipoles and Ferroelectric Distortion in BaTiO 3 Nanocubes

Here, surface adsorbates and surrounding matrix species have been demonstrated to affect the properties of nanoscale ferroelectrics and nanoscale ferroelectric composites; potentially counteracting performance losses that can occur in small particle sizes. In this work, the effects of nonpolar oleic acid (OA) and polar tetrafluoroborate (BF 4 - ) ligand capping on the surface of various sizes of BaTiO 3 nanocubes have been investigated with combined neutron diffraction and neutron pair distribution function (PDF), density functional theory (DFT), and ab initio molecular dynamics (AIMD) methods. The low real space PDF region provides an unobstructed view of rhombohedral (split short and long) Ti-O distances in BaTiO 3 nanocubes, mimicking the well-established order-disorder local structure found in bulk BaTiO 3 . Interestingly, the intermediate-range order in nanocubes is found to be orthorhombic, rather than tetragonal. It is concluded that polar ligands adsorbed at BaTiO 3 surfaces stabilize the correlation length scale of local rhombohedral distortions in ferroelectric nanoparticles relative to nonpolar ligands.

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

Barriers to carriers: faults and recombination in non-stoichiometric perovskite scintillators

Abstract Tuning the efficiency and speed of charge carrier recombination in inorganic scintillators can potentially improve their performance in diverse applications. Recent work suggests that this maybe be achieved via a two-phase scintillator AB that naturally phase separates into A -rich and B -rich domains. In addition, a favorable electronic structure and band-edge alignment such that the charge carriers are confined or are thermodynamically driven to preferentially accumulate in one of the two domains, might lead to an improved radiative recombination rate. Here, we use density functional theory computations and ab initio molecular dynamics (AIMD), including non-adiabatic molecular dynamics (NAMD) simulations, to examine an alternative phase structure and its potential impact on recombination. Using a model perovskite SrTiO $$_3$$ 3 system with one-, two- and three-dimensional Ruddlesden–Popper (RP) phases, we demonstrate that RP faults induce band structure changes in the material that can act as barriers to carrier transport. Our AIMD/NAMD simulations indicate competing effects of a lower mean free path (potentially enhancing the desired radiative recombination and overall scintillating efficiency) and faster non-radiative recombination (undesired) due to enhanced electron–phonon coupling in the faulted system. Full exploitation of such a rational design approach would require tuning of the effective scintillation efficiency by varying the perovskite chemistry using appropriate arrangements of RP faults in the bulk material. Finally, other effects, such as the tendency of point defects to segregate at the interface, that might affect the overall performance, are briefly discussed. We expect the basic results found here to apply to other nanostructured scintillators. Graphical Abstract

36 MATERIALS SCIENCE↗
DOE OSTI2020

Predicting mesoscale spectral thermal conductivity using advanced deterministic phonon transport techniques

We present a review, demonstration and simulation of phonon transport for the purposes of predicting materials performance at the mesoscale. We focus primarily on the development and implementation of a unified methodology to enable predictive heat transport. We report on the current state of the art as it pertains to deterministic phonon transport methodologies, discussing various topics concerning phonons. In application, we focus on the Self-Adjoint Angular Flux (SAAF) formulation of the Boltzmann transport equation for phonons, and develop the spatial, angular, and material property discretization required to accurately simulate the predictive physics of heat transport in dielectrics. We discuss thermal interfacial resistance and present our formulation of the diffuse mismatch model for simulating phonon interactions at internal boundaries. We have recently developed a deterministic, spectral phonon transport method for predicting effective thermal conductivity ($\kappa_{\textrm{eff}}$), using Bose-Einstein source terms coupled through an average material temperature. In this method, we introduce a closure term to the phonon transport system which acts as a redistribution function for the total energy of the system, and serves as a guide for the amount of non-equilibrium behavior occurring in the system. This method predicts thermal conductivity and equilibrium temperature distributions in homogeneous and heterogeneous materials using data generated by ab initio density functional theory methods. We employ polarization, density of states and full dispersion spectra to resolve thermal conductivity with numerous angular and spatial discretizations. Our implementation utilizes a Richardson iteration on a modified version of the phonon scattering source. The equations associated with this method are solved via a modification of traditional source iteration. We compare the performance of source iteration applied to an existing uncoupled, traditional SAAF method to our new method and comment on the iterative performance of each. We observe ballistic and diffusive phonon scattering as acoustic thickness of the domain changes, and are able to make comparisons between the accuracy and efficiency of both methods.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Interactive multiscale modeling to bridge atomic properties and electrochemical performance in Li-CO 2 battery design

Li-CO 2 batteries are promising energy storage systems due to their high theoretical energy density and CO 2 fixation capability, relying on reversible Li 2 CO 3 /C formation during discharge/charge cycles. Here, we present a multiscale modeling framework integrating Density Functional Theory (DFT), Ab-Initio Molecular Dynamics (AIMD), classical Molecular Dynamics (MD), and Finite Element Analysis (FEA) to investigate atomic and cell-level properties. The considered Li-CO 2 battery consists of a lithium metal anode, an ionic liquid electrolyte, and a carbon cloth cathode with Sb 0.67 Bi 1.33 Te 3 catalyst. DFT and AIMD determined the electrical conductivities of Sb 0.67 Bi 1.33 Te 3 and Li 2 CO 3 using the Kubo–Greenwood formalism and studied the CO 2 reduction mechanism on the cathode catalyst. MD simulations calculated the CO 2 diffusion coefficient, Li + transference number, ionic conductivity, and Li + solvation structure. The FEA model, parameterized with atomistic simulation data, reproduced the available experimental voltage–capacity profile at 1 mA/cm 2 and revealed spatio-temporal variations in Li 2 CO 3 /C deposition, porosity, and CO 2 concentration dependence on discharge rates in the cathode. Accordingly, Li 2 CO 3 can form large and thin film deposits, leading to dispersed and local porosity changes at 0.1 mA/cm 2 and 1 mA/cm 2 , respectively. The capacity decreases exponentially from 81,570 mAh/g at 0.1 mA/cm 2 to 6200 mAh/g at 1 mA/cm 2 , due to pore clogging from excessive discharge product deposition that limits CO 2 transport to the cathode interior. Therefore, the performance of Li-CO 2 batteries can be improved by enhancing CO 2 transport, regulating Li 2 CO 3 deposition, and optimizing cathode architecture.

Battery performance↗
DOE OSTI2020

First-principle investigation on catalytic hydrogenation of benzaldehyde over Pt-group metals

We report that understanding the hydrogenation of organic compounds in the aqueous phase has always been fundamentally important for improving carbon neutral pathways to fuels and value-added chemicals. In this study, we investigated both thermodynamic and kinetic profiles of benzaldehyde hydrogenation over the Pd(111) and Pt(111) metal surfaces using density functional theory (DFT) and ab initio molecular dynamic (AIMD) simulations. The adsorption of H 2 shows the mixed preference of H adsorption sites on the Pt(111), while the fcc adsorption site is dominant for H on the Pd(111). When benzaldehyde is added to the systems, we observe a strong reduction of benzaldehyde on charged Pd (111) surface compared with that on neutral surface. In contrast, charged state of the Pt(111) surface does not change their interaction. Subsequent hydrogenation reaction of benzaldehyde over Pd(111), proceeding via Langmuir-Hinshelwood mechanism, is affected by two major factors: the presence of H 2 O solvent and surface charge. The presence of H 2 O solvent greatly reduces the activation energy of C-H and O-H bond formation during the hydrogenation process. Furthermore, the hydrogenation step via C-H bond formation is preferred thermodynamically and kinetically over O-H bond formation during thermocatalytic hydrogenation, while the opposite trend holds true during electrocatalytic hydrogenation.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

Mechanistic study of direct coupling of CO 2 and C 2 H 4 over atomically dispersed metal at graphene edges

Direct coupling of CO 2 and ethylene (hereinafter DCCE) to acrylic acid is valuable for valorizing CO 2 to manufacture acrylate-derived products. However, previous studies in DCCE have been limited on molecular catalysts with challenges in improving catalytic performance. In this work, we employed density functional theory calculations and ab initio molecular dynamics simulations to investigate the heterogeneous catalysis of DCCE over atomically dispersed metal centers at nitrogen-doped zigzag edge of graphene. Based on competitive adsorption and structural stability, Mo, Cr, V, Ru, and Ni active sites are chosen to explore the reaction kinetics. Here, we find that the activation barriers are determined by the charge redistribution at transition states, which explains the trend of activity for the C-C coupling and the hydrogen transfer, two key steps in DCCE. Furthermore, we show that the intramolecular hydrogen transfer (rate-limiting step) is hindered due to the lack of local coordinate at the active sites. We thus propose to use co-adsorbed water as a “proton-exchanger” following a water-assisted route, and show that the activation barriers are reduced over all metal centers. Particularly, water promotes the hydrogen transfer over metals with strong CO 2 -ethylene co-activation and facile C-C coupling kinetics, which could be considered promising for DCCE. In both mechanisms, the stability of metallactone intermediate can be used to predict the catalytic activity. It is anticipated that the insights from this work can provide guidelines for mimicking well-defined multifunctional active sites in molecular catalysts to design heterogeneous catalysts for such C-C coupling, which advances catalytic utilization of CO 2 .

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

Experimental and modeling study of the autoignition behavior of a saturated heterocycle: Pyrrolidine

Experiments are conducted in both rapid compression machine (RCM) and shock tube (ST) to better quantify autoignition behavior (e.g., ignition delay, heat release) and understand heteroatomic effects in heterocyclic compounds, which are important reference components for the combustion of biomass-derived liquid fuels. These tests focus on the nitrogen-containing, five-membered saturated ring, pyrrolidine, at diluted conditions covering pressures of 20 and 50 bar, temperatures of 720–1450 K and a range of stoichiometries (ϕ = 0.5–2). A chemical kinetic model is developed and coupled to an existing combustion kinetics framework describing key nitrogen containing intermediates (e.g. pyrrole, ammonia and NOx). H-abstraction reactions by OH, H, CH 3 and HO 2 , are determined using ab-initio transition state theory methods, while analogies to cyclopentane are adopted for many other reactions, such as ring-opening. Here, the autoignition measurements reveal the lack of negative temperature coefficient (NTC) behavior and low-temperature chemistry for pyrrolidine, as opposed to its saturated hydrocarbon analogue, cyclopentane. Interestingly, at the lowest temperatures (T < 750 K), the reactivity of cyclopentane is greater than pyrrolidine, while at higher temperatures, pyrrolidine becomes more reactive. Agreement between the experimental measurements and the model is good, and it is found that H-abstraction reactions by HO 2 and ensuing chemistry play key roles in controlling the reactivity of this cyclic amine. Most of the fluxes, i.e., >70 %, are predicted to move through 1- or 2-pyrroline (C 4 H 7 N) and then the cyclic C 4 H 6 N radical, at both lower and higher temperatures, to form either CH 2 CHCHCHNH via ring-opening or pyrrole via β-scission. It appears that the ring opens more easily at lower temperature whereas the C–H β-scission dominates at higher temperature and lower pressure, such that the reaction of the fuel radical intermediate carrying an unpaired electron on the nitrogen atom with HO 2 is the next most notable in promoting oxidation. When comparing pyrrolidine and cyclopentane, which exhibits distinct pathways in different temperature regimes, the pyrrolidine pathways and sensitivity analysis align more closely to the high temperature case of cyclopentane where the important role of HO 2 radicals is seen to provide chain branching through HO 2 reaction with the fuel, accompanied by H 2 O 2 formation and decomposition to OH. The formation of 5-membered diene rings and ring opening reactions are also found to be highly relevant. Of particular note, it is found that there is little influence of small molecule nitrogen-chemistry, e.g., NH 2 , HCN, NO/NO 2 on the reactivity of the pyrrolidine mixtures investigated here where no recirculated combustion gases are included.

Autoignition↗