SEARCH · Search NASA

Results for “Bimolecular reactions”

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.

SourceResearch connectionFrom yearThrough year

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

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Bimolecular Chemistry in the Ultracold Regime

Advances in atomic, molecular, and optical physics techniques allowed the cooling of simple molecules down to the ultracold regime ([Formula: see text]1 mK) and opened opportunities to study chemical reactions with unprecedented levels of control. This review covers recent developments in studying bimolecular chemistry at ultralow temperatures. We begin with a brief overview of methods for producing, manipulating, and detecting ultracold molecules. We then survey experimental works that exploit the controllability of ultracold molecules to probe and modify their long-range interactions. Further combining the use of physical chemistry techniques such as mass spectrometry and ion imaging significantly improved the detection of ultracold reactions and enabled explorations of their dynamics in the short range. We discuss a series of studies on the reaction KRb + KRb → K 2 + Rb 2 initiated below 1 μK, including the direct observation of a long-lived complex, the demonstration of product rotational state control via conserved nuclear spins, and a test of the statistical model using the complete quantum state distribution of the products.

Chemistry↗

DOE OSTI2026

ALD with Alternative Coreactants: Which Work, Which Do Not, and Why

In this paper, we present a combined experimental and theoretical study that systematically evaluates a series of alcohols (i.e., primary, secondary, and tertiary) as coreactants for trimethylaluminum (TMA)-based atomic layer deposition (ALD) of Al 2 O 3 . We employed in situ quartz crystal microbalance techniques and ex situ X-ray photoelectron spectroscopy to probe growth across a range of temperatures from T = 120 to 285 °C. Dispersion-inclusive hybrid density functional theory was employed to identify potential reaction pathways and compute corresponding Gibbs free-energy barriers and rate coefficients. Experimentally on an Al 2 O 3 surface, sustained thin-film growth, where the rate is constant over many cycles, was only observed with tertiary alcohols (tert-butanol and 2-methyl-2-butanol) at elevated temperatures (e.g., T = 285 °C); for tert-butanol specifically, sustained growth also occurred at T = 240 °C, with no sustained growth at or below T = 210 °C. Growth using primary and secondary alcohols decayed to a negligible value after only 2–3 cycles on an Al 2 O 3 surface at all tested temperatures. Intentional trace additions of H 2 O to anhydrous alcohols restored sustained growth but at reduced deposition rates relative to pure H 2 O. Theoretical analysis supports a bimolecular alkoxy β-H elimination mechanism as the primary reaction pathway: initiated by the formation of a bound alkoxy species, the subsequent β-H elimination step is rate-determining, with a barrier that systematically decreases as the degree of the alcohol coreactant is increased from primary to secondary to tertiary. Computed rate coefficients (k) also indicate a strong temperature dependence (with k increasing by ∼4–5 orders of magnitude from T = 120 to 285 °C) and are consistent with the experimental observation that only tertiary alcohols lead to steady thin-film growth at elevated temperatures. This integrated approach establishes why only certain alcohols can sustain Al 2 O 3 ALD and delivers a predictive framework for identifying effective alcohol coreactants.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Styrene Thermal Decomposition and Its Reaction with Acetylene under Shock Tube Pyrolysis Conditions: an Experimental and Kinetic Modeling Study

Styrene is an important compound for polymer production and a key intermediate in gas-phase kinetics of polycyclic aromatic hydrocarbons (PAHs). For the first time, the pyrolysis of styrene with and without the presence of acetylene is investigated in a single-pulse shock tube coupled to gas chromatography and mass spectrometry. For each reaction system, quantitative speciation profiles are probed within the temperature range of 1100-1730 K, nominal pressure of 20 bar, and reaction duration of similar to 4 ms. A kinetic model is built to simulate the results. The model explains how styrene is consumed under high-pressure pyrolytic conditions, how the secondary chemistry of intermediate products affect subsequent PAH formation, and how acetylene addition alters the reaction pathways. Throughout the temperature range, styrene breakdown is dominated by the bimolecular interaction between styrene and hydrogen atom, which produces benzene+vinyl or phenyl and ethylene through the stabilization of 2-phenylethyl and its subsequent dissociation. As a result, large amounts of phenyl accumulate, which react with styrene to form C14H12 species while simultaneously releasing H atoms through addition/elimination reactions. The reactivity of fuel consumption is preserved by the regeneration of H atoms as chain carriers. Several C14H10 compounds are formed as a result of the following breakdown of the C14H12 isomers, particularly stilbene, 1,1-diphenyl ethylene, and 9-methyl-9H-fluorene. The presence of acetylene as a co-reactant with styrene allows the Hydrogen-Abstraction-Acetylene-Addition (HACA) pathway to proceed from phenyl radical to enhance the production of phenylacetylene at very low temperatures and acenaphthylene. This hinders the formation of C14H12 isomers, exclusive products from pure styrene dissociation by competing with the styrene+phenyl routes.

Kinetic Modeling↗

DOE OSTI2022

Catalytic Upgrading of Ethanol to C8+ Distillate Range Ethers via Guerbet Coupling and Etherification

Guerbet coupling followed by etherification is a two step process for conversion of ethanol into distillate range molecules. This involves first alcohol oligomerization, followed by acid-catalyzed dehydration. As we will show in this presentation, ethanol can be converted to higher linear and α-branched alcohols by C-C coupling reactions (ethanol oligomerization) as well as to high molecular weight esters by C-O coupling reactions with Cu/MgxAlOy (CuHT) catalysts. Alcohols may be converted to ethers via bimolecular dehydration in a subsequent step. We have identified more than 160 number of species in the products including alcohols, esters, aldehydes, ketonces, and olefins. Alcohols range from C4 to C10. Both alcohols and esters follow a Schulz-Flory chain growth model. We show the relationship between the catalyst properties (BET surface area, acid and base site count and Cu loading and synthesis method) and the performance in the reactions. We also show that physical mixtures of CuHT and HT can have similar product selectivity of low loading CuHT catalysts. The selectivity towards diesel fuel precursor compounds (hereafter ‘DFPC’) increased with conversion until reaching a plateau at high ethanol conversion (~70%). Alcohol selectivity follows a Schultz-Flory distribution at all studied conversions, and that adsorbed ethanol-derived species may undergo surface oligomerization into 1-butanol and higher alcohols before desorbing in a chain-growth mechanism. Zeolite catalysts convert the C4+ alcohols into C8+ ethers in both batch and continuous flow reactors. Selectivities of up to 80% to C8+ ethers at around 70% conversion using a single pass continuous flow system are able to be achieved. The final product feedstocks obtained from etherification have been used in technoeconomic (TEA) and lifecycle analysis (LCA), which indicate a reduction in greenhouse gas (GHG) emissions of 50% relative to conventional diesel. These final feedstocks are undergoing engine testing to elucidate the physiochemical properties and will be compared to diesel #2 ASTM standards as a reference and will serve as a basis to improve the final bio-diesel blendstocks from etherification.

Huber, George↗

DOE OSTI2023

Catalytic upgrading of ethanol to C8+ distillate range ethers via Guerbet coupling and etherification

In this presentation we describe the catalytic approaches for conversion of ethanol into diesel fuel ethers. This involves first alcohol oligomerization, followed by acid-catalyzed dehydration. As we will show in this presentation, ethanol can be converted to higher linear and α-branched alcohols by C-C coupling reactions (ethanol oligomerization) as well as to high molecular weight esters by C-O coupling reactions with Cu/MgxAlOy (CuHT) catalysts. Alcohols may be converted to ethers via bimolecular dehydration in a subsequent step. We have identified more than 160 number of species in the products including alcohols, esters, aldehydes, ketones, and olefins. Alcohols range from C4 to C10. Both alcohols and esters follow a Schulz-Flory chain growth model. We show the relationship between the catalyst properties (BET surface area, acid and base site count and Cu loading and synthesis method) and the performance in the reactions. We also show that physical mixtures of CuHT and HT can have similar product selectivity of low loading CuHT catalysts. The selectivity towards diesel fuel precursor compounds (hereafter ‘DFPC’) increased with conversion until reaching a plateau at high ethanol conversion (~70%). Alcohol selectivity follows a Schultz-Flory distribution at all studied conversions, and that adsorbed ethanol-derived species may undergo surface oligomerization into 1-butanol and higher alcohols before desorbing in a chain-growth mechanism. Zeolite catalysts convert the C4+ alcohols into C8+ ethers in both batch and continuous flow reactors. Selectivities of up to 80% to C8+ ethers at around 70% conversion using a single pass continuous flow system are achieved. The final product feedstocks obtained from etherification have been used in technoeconomic (TEA) and lifecycle analysis (LCA), which indicate a reduction in greenhouse gas (GHG) emissions of 50% relative to conventional diesel and diesel fuel prices that are lower than biodiesel. These final feedstocks are undergoing engine testing to elucidate the physiochemical properties and compare well to diesel #2 ASTM standards.

Huber, George↗

DOE OSTI2023

Sulfur ( 3 P) Reaction with Conjugated Dienes Gives Cyclization to Thiophenes under Single Collision Conditions

Here, we combine crossed-beam velocity map imaging with high-level ab initio/transition state theory modeling of the reaction of S( 3 P) with 1,3-butadiene and isoprene under single collision conditions. For the butadiene reaction, we detect both H and H 2 loss from the initial adduct, and from reaction with isoprene, we see both H loss and methyl loss. Theoretical calculations confirm these arise following intersystem crossing to the singlet surface forming long-lived intermediates. For the butadiene reaction, these lose H 2 to form thiophene as the dominant channel, H to form the detected 2H-thiophenyl radical, or ethene, giving thioketene. For isoprene, additional reaction products are suggested by theory, including the observed H and methyl loss radicals, but also methyl thiophene, thioformaldehyde, and thioketene. The results for S( 3 P) + 1,3-butadiene, showing direct cyclization to the aromatic product and yielding few bimolecular product channels, are in striking contrast to those for the analogous O( 3 P) reaction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Renormalization of excitonic properties by polar phonons

Here we employ quasiparticle path integral molecular dynamics to study how the excitonic properties of model semiconductors are altered by electron–phonon coupling. We describe ways within a path integral representation of the system to evaluate the renormalized mass, binding energy, and radiative recombination rate of excitons in the presence of a fluctuating lattice. To illustrate this approach, we consider Fröhlich-type electron–phonon interactions and employ an imaginary time influence functional to incorporate phonon-induced effects nonperturbatively. The effective mass and binding energies are compared with perturbative and variational approaches, which provide qualitatively consistent trends. We evaluate electron-hole recombination rates as mediated through both trap-assisted and bimolecular processes, developing a consistent statistical mechanical approach valid in the reaction limited regime. These calculations demonstrate how phonons screen electron–hole interactions, generically reducing exciton binding energies and increasing their radiative lifetimes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Gas-Phase Formation of 1,3,5,7-Cyclooctatetraene (C 8 H 8 ) through Ring Expansion via the Aromatic 1,3,5-Cyclooctatrien-7-yl Radical (C 8 H 9 • ) Transient

Gas-phase 1,3,5,7-cyclooctatetraene (C 8 H 8 ) and triplet aromatic 1,3,5,7-cyclooctatetraene (C 8 H 8 ) were formed for the first time through the bimolecular methylidyne radicals (CH) - 1,3,5-cycloheptatriene (C 7 H 8 ) reactions under single collision conditions on the doublet surface. The reaction involves methylidyne radical addition to the olefinic pi electron system of 1,3,5-cycloheptatriene followed by isomerization and ring expansion to aromatic 1,3,5-cyclooctatrien-7-yl radical (C 8 H 9 •). The chemically activated doublet radical intermediate undergoes unimolecular decomposition to 1,3,5,7-cyclooctatetraene. Substituted 1,3,5,7-cyclooctatetraene molecules can be prepared in the gas phase with hydrogen atom(s) in the 1,3,5-cycloheptatriene reactant being replaced by organic side groups. Furthermore, these findings are also of potential interest to organometallic chemists by expanding the synthesis of exotic transition metal complexes incorporating substituted 1,3,5,7-cyclooctatetraene dianion (C 8 H 8 2- ) ligands and to untangle the unimolecular decomposition of chemically activated and substituted 1,3,5-cyclooctatrien-7-yl radical eventually gaining a fundamental insight of their bonding chemistry, electronic structures, and stabilities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Radical–Radical Reactions in Molecular Weight Growth: The Phenyl + Propargyl Reaction

The mechanism for hydrocarbon ring growth in sooting environments is still the subject of considerable debate. The reaction of phenyl radical (C 6 H 5 ) with propargyl radical (H 2 CCCH) provides an important prototype for radical–radical ring-growth pathways. We studied this reaction experimentally over the temperature range of 300–1000 K and pressure range of 4–10 Torr using time-resolved multiplexed photoionization mass spectrometry. We detect both the C 9 H 8 and C 9 H 7 + H product channels and report experimental isomer-resolved product branching fractions for the C 9 H 8 product. We compare these experiments to theoretical kinetics predictions from a recently published study augmented by new calculations. Here, these ab initio transition state theory-based master equation calculations employ high-quality potential energy surfaces, conventional transition state theory for the tight transition states, and direct CASPT2-based variable reaction coordinate transition state theory (VRC-TST) for the barrierless channels. At 300 K only the direct adducts from radical–radical addition are observed, with good agreement between experimental and theoretical branching fractions, supporting the VRC-TST calculations of the barrierless entrance channel. As the temperature is increased to 1000 K we observe two additional isomers, including indene, a two-ring polycyclic aromatic hydrocarbon, and a small amount of bimolecular products C 9 H 7 + H. Our calculated branching fractions for the phenyl + propargyl reaction predict significantly less indene than observed experimentally. We present further calculations and experimental evidence that the most likely cause of this discrepancy is the contribution of H atom reactions, both H + indenyl (C 9 H 7 ) recombination to indene and H-assisted isomerization that converts less stable C 9 H 8 isomers into indene. Especially at low pressures typical of laboratory investigations, H-atom-assisted isomerization needs to be considered. Regardless, the experimental observation of indene demonstrates that the title reaction leads, either directly or indirectly, to the formation of the second ring in polycyclic aromatic hydrocarbons.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Divergent Bimetallic Mechanisms in Copper(II)-Mediated C–C, N–N, and O–O Oxidative Coupling Reactions

Copper-catalyzed aerobic oxidative coupling of diaryl imines provides a route for conversion of ammonia to hydrazine. Here the present study uses experimental and density functional theory computational methods to investigate the mechanism of N–N bond formation, and the data support a mechanism involving bimolecular coupling of Cu-coordinated iminyl radicals. Computational analysis is extended to Cu II -mediated C–C, N–N, and O–O coupling reactions involved in the formation of cyanogen (NC–CN) from HCN, 1,3-butadiyne from ethyne (i.e., Glaser coupling), hydrazine from ammonia, and hydrogen peroxide from water. The results reveal two different mechanistic pathways. Heteroatom ligands with an uncoordinated lone pair (iminyl, NH 2 , OH) undergo charge transfer to Cu II , generating ligand-centered radicals that undergo facile bimolecular radical–radical coupling. Ligands lacking a lone pair (CN and CCH) form bridged binuclear diamond-core structures that undergo C–C coupling. This mechanistic bifurcation is rationalized by analysis of spin densities in key intermediates and transition states, as well as multiconfigurational calculations. Radical–radical coupling is especially favorable for N–N coupling owing to energetically favorable charge transfer in the intermediate and thermodynamically favorable product formation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Theoretical study of the reaction mechanism and kinetics of the phenyl + propargyl association

Potential energy surface for the phenyl + propargyl radical recombination reaction has been studied at the CCSD(T)-F12/cc-pVTZ-f12//B3LYP/6-311G** level of theory for the closed-shell singlet species and at the triplet–singlet gap CASPT2/cc-pVTZ-CCSD(T)-F12/cc-pVTZ-f12//CASSCF/cc-pVTZ level of theory for the diradical species. High-pressure limit rate constants for the barrierless channels were evaluated with variable reaction coordinate transition state theory (VRC-TST). Rice–Ramsperger–Kassel–Marcus Master Equation (RRKM-ME) calculations have been performed to assess temperature- and pressure-dependent phenomenological rate constants and product branching ratios. Here, the entrance channels of the radical association reaction produce 3-phenyl-1-propyne and phenylallene which can further dissociate/isomerize into a variety of unimolecular and bimolecular products. Theoretical evidence is presented that, at combustion relevant conditions, the phenyl + propargyl recombination provides a feasible mechanism for the addition of a second five-member ring to the first six-member aromatic ring producing the prototype two-ring species indene and indenyl. Rate expressions for all important reaction channels in a broad range of temperatures and pressures have been generated for kinetic modeling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Atomic-scale mechanism of carbon nucleation from a deep crustal fluid by replica exchange reactive molecular dynamics simulation

Here we present a mechanistic model of carbon nucleation and growth from a fluid at elevated temperature (T) and pressure conditions, typical of those found in the shallow Earth’s lithosphere. Our model uses a replica exchange reactive molecular dynamics framework in which molecular configurations are swapped between adjacent T replica at regular intervals according to underlying statistical mechanics. This framework allows predicting complex molecular structures and thermodynamics while remaining computationally efficient. Here we simulate the reactivity of an unstable mixture of CO 2 and CH 4 at 1000 K and 1 GPa. We find that the path to thermodynamic equilibrium is initially entropy-driven, producing a diversity of short-lived species, including various alcohols with intermediate carbon oxidation states. Cyclic and polycyclic radicals that are sometimes resonance-stabilized form next and set the stage for carbon nucleation. The carbon exsolution process releases abundant water, is exothermic and starts with the nucleation of a large aggregate of hydrogenated graphene flakes from covalently bonded polycyclic units. The carbon backbone of this nucleus subsequently grows into a hydrogen-depleted fullerene-like structure, before evolving toward a partially bilayered graphene layer. Overall, our results show that the mechanism of graphitic C formation is certainly not bimolecular, and that it may involve a combination of key condensation and radical chain reactions. This will help understand the isotopic, and reactive characteristics of carbon-bearing fluids during their upward transit through the Earth’s mantle and crust. Moreover, the mechanistic insights outlined here present intriguing similarities with the process of soot and interstellar dust formation, which suggests that the widespread distribution of abiotic polyaromatic and graphitic material on Earth and beyond may reflect the prevalence of a fundamental chemical pathway.

58 GEOSCIENCES↗

DOE OSTI2026

Data for "Photoinduced Chemomimetic Biocatalysis for Enantioselective Intermolecular Radical Conjugate Addition"

Exploiting nature’s catalysts for non-natural transformations that are inaccessible to chemocatalysis is highly desirable but challenging. On the one hand, the widespread nicotinamide-dependent oxidoreductases have not been utilized for single-electron-transfer-induced bimolecular cross-couplings; on the other, the addition of catalytic asymmetric radical conjugate to terminal alkenes remains a challenge owing to strong racemic background reaction and unselective termination of prochiral radical species. Here we report a chemomimetic biocatalysitic approach for construction of alpha-carbonyl stereocentres via an unnatural intermolecular conjugate addition of N-(acyloxy)phthalimides-derived radicals with acceptor-substituted terminal alkenes, by combination of visible-light excitation and nicotinamide-dependent ketoreductases (KREDs). Based on protein crystal structure, we engineered KREDs via a semi-rational mutagenesis strategy to improve reaction outcomes with a small and high-quality variants library. Mechanistic investigations combining wet experiments, crystallographic studies and computational simulations demonstrate that the repurposed biocatalyst can suppress racemic background reaction and unselected side reactions, yielding enantioselectivity that is challenging to achieve by chemocatalysis.

Catalysis↗

DOE OSTI2020

Development of a high temperature pyrolysis mechanism for cyclopentanone, a potential biofuel derived from biomass

Cyclopentanone (CPO) is a promising biofuel for spark-ignition engines due to its ring strain and high auto-ignition resistance. The Department of Energy has identified CPO as a potential biofuel candidate because of its unique chemical properties. Understanding CPO decomposition is crucial to building a high-temperature combustion model. In this study, we present a comprehensive kinetic model for high-temperature pyrolysis of CPO with verified results from high-pressure shock tube (HPST) measurements. The time histories of carbon monoxide (CO), ethylene (C2H4), and CPO absorbances over the temperature range of 1156-1416 K and 8.53-10.06 atm were measured. A corresponding kinetic model was generated using the Reaction Mechanism Generator (RMG) with the parameters of dominant pathways from either previous studies or quantum calculation within current work. The final model containing 821 species and 79,859 reactions exhibited a good agreement with the experimental results. In this study, the C2H4/CO ratio was used as an essential factor to validate the model and to prove that radical-involved bimolecular pathways were as significant as unimolecular decomposition of CPO. The rate of production (ROP) analysis showed that H radicals play major roles in both primary reactions and product pyrolysis. The insights from this work can be used to generate a better CPO combustion model and help evaluate CPO as a biofuel. Such understanding can also be helpful when proposing pyrolysis or combustion models for other biofuels.

cyclopentanone, reaction mechanism, pyrolysis, sho↗

DOE OSTI2021

Molecularly tunable heterogeneous catalysts by edge functionalization of graphitic carbons

Disclosed are surface immobilized (electro)catalysts that may be prepared by a condensation reaction that generates an aromatic unit that is robust to acid and base and elevated temperatures. Among their many desirable characteristics, the catalysts are far less prone to the bimolecular deactivation pathways commonly observed for homogeneous catalysts, and may be used in solvents with a range of polarities and dielectric strengths. The catalysts are suitable for a wide array of thermal catalytic reactions (polymerization, oxidation, hydrogenation, cross-coupling etc.) and as anodes and/or cathodes in fuel cells, electrolyzers, and in batteries and supercapacitors.

Surendranath, Yogesh↗

DOE OSTI2020

Oxide-catalyzed self- and cross-condensation of cycloketones. Kinetically relevant steps that determine product distribution

Metal oxides with acid-base properties can be used to catalyze aldol condensation and potentially obtain valuable chemicals from biomass-derived oxygenates. Understanding the reaction mechanisms could enable the development of effective catalysts with controllable selectivity. In this work, we investigate the reaction parameters that affect the kinetics and product distribution of self- and cross-aldol condensation of cyclopentanone (CPO) and cyclohexanone (CHO) over CeO 2 and ZrO 2 catalysts by combining detailed analysis of experiments and DFT calculations. It is well-known that either the deprotonation of an α-C–H that generates an enolate or the C–C coupling that produces the dimers can be rate-limiting step. Indeed, the detailed kinetics analysis shows that the observed data can be fitted by unimolecular or bimolecular Langmuir-Hinshelwood models. Specifically, on the ZrO 2 catalyst, the unimolecular α-C–H activation is the rate limiting step for the CPO self-condensation, while the bimolecular C–C coupling limits the CHO self-condensation. By contrast, on CeO 2 , the bimolecular C–C coupling step seems to be rate-limiting for both CPO and CHO self-condensation reactions. The adsorption parameters derived from kinetics fitting indicate that when the two reactants are co-fed, the surface is preferentially covered by CHO. So, the condensation rate of [CPO]-activated products ([CPO]CPO + [CPO]CHO) decreases in the presence of CHO, due to the competitive adsorption. However, despite a higher surface coverage, CHO is not an effective electrophile because it imposes stronger steric constraints to the formation of the C–C bond than CPO. Therefore, for a given enolate, the product distribution shows that [CPO]CPO > [CPO]CHO and, similarly [CHO]CPO > [CHO]CHO. Likewise, the yields of [CPO]- and [CHO]-activated products in mixed CPO/CHO feeds are similar because the higher surface coverage and stability of the [CHO] enolate is counterbalanced by the inhibition of C–C coupling by steric constraints when CHO is involved. In summary, the overall rank of products observed experimentally can be explained in terms of a combination of the following effects: relative surface coverage of the reactants, steric constraints for the C–C coupling, stability of the enolate, and effectiveness of the electrophile to accept an electron.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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 OSTI2024

Hydrogen ejection from hydrocarbons: Characterization and relevance in soot formation and interstellar chemistry

Polycyclic aromatic hydrocarbons (PAHs) play a major role in the chemistry of combustion, pyrolysis, and the interstellar medium. Production (or activation) of radical PAHs and propagation of their resulting reactions require efficient dehydrogenation, but the preferred method of hydrogen loss is not well understood. Unimolecular hydrogen ejection (i.e., direct C─H bond fission) and bimolecular radical abstraction are two main candidate pathways. We performed a computational study to characterize the role of H ejection, particularly as a driver for radical-centric hydrocarbon-growth mechanisms and particle formation. Electronic structure calculations establish that C─H bond strengths span a broad range of energies, which can be weaker than 30 kcal/mol in some C 9 and C 13 PAH radicals. At T > 1200 K, calculated thermal rates for hydrogen ejection from weak C─H bonds at zigzag sites on PAH radicals are significantly larger than typical H-abstraction rates. These results are highly relevant in the context of chain reactions of radical species and soot inception under fuel-rich combustion conditions. Furthermore, calculated microcanonical rates that include the additional internal energy released by bond formation (e.g., ring closure to yield C 9 H 9 ) yield significantly higher rates than those associated with full thermalization. These microcanonical considerations are relevant to the astrochemical processes associated with hydrocarbon growth and processing in the low-density interstellar environment.

08 HYDROGEN↗