Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reaction”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

NCERC-CL TREAT Free Field Characterization Foil Reaction Rate Results

This report summarizes the analysis and preliminary MCNP reaction rate modeling of free field characterization measurements performed in the Transient Reaction Test Facility (TREAT) at Idaho National Laboratory (INL) between 12/18/2023 and 2/22/2024. A series of seven irradiations were performed using the Big-BUSTER (Broad Use Specimen Transient Experiment Rig) core configuration. Foil sets including 19.5% LEU-Zr alloy wire, S, Au, Fe, Ni, Co, Ti, and Zr were deployed in each irradiation. Identical foil sets were supplied to Los Alamos National Laboratory (LANL), Lawrence Livermore National Laboratory (LLNL), and INL. All foils were supplied by INL. Additional details on the experiment can be found in [2]. The foil sets were shipped to the National Criticality Experiments Research Center Counting Laboratory (NCERC-CL) branch at LANL (NCERC-CL NISC) and were received on 3/25/2024. Additional details on receipt and counting are reported in [3]. All reactor dosimetry measurements were performed in adherence to the ASTM standards applicable to reactor dosimetry. Measured reaction rates from Au, Co, Fe, Ni, and Ti are reported. Select reaction rate ratios and simulated reaction rates and ratios are discussed. Reaction rates for all measured reaction products were modeled with MCNP 6.3.1 and a TREAT input deck supplied by Edward Lum. Neutron emission estimates calculated using the measured and simulated reaction rates show excellent agreement between the three fast-threshold n-p reactions with the percent differences being < 9% between the three reactions. The capture reactions had poorer, but still reasonable, agreement with < 20% percent differences between the neutron emission estimates of the three capture reactions. The percent difference between the fast-threshold and capture reactions is ∼ 100% for each reaction. This supports the TREAT model is not accurately modeling the neutron spectrum and additional measurements are required to characterize TREAT and match measurements to simulation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Thermodynamically consistent forward and reverse degrees of rate control in reversible reactions

The net rate of a composite reaction is the difference between the forward and reverse reaction rates, which are kinetically distinct despite sharing elementary reaction steps and therefore have different rate-controlling transition states and species. Thus, degrees of rate control defined to identify rate-controlling transition states and species for the net rate confound contributions from the forward and reverse reactions. Herein, the forward and reverse degrees of rate control are defined to independently quantify the extent that species and transition states control the forward and reverse rates in reversible reactions. These degrees of rate control are defined as the relative change in the forward and reverse reaction rates per k B T decrease in the standard-state molecular free energies of transition states and species, and they are related to experimentally measurable quantities such as steady-state fractional coverages, reaction orders, and energies and entropies of activation of the forward and reverse reaction rates. The forward/reverse degrees of rate control represent stoichiometric coefficients for species and transition states in an equilibrium between the apparent initial states and the apparent transition states of the apparent rate-controlling steps of forward and reverse overall reactions. At equilibrium, the apparent transition states for the forward reaction and reverse reaction converge, and thus the forward and reverse rate-controlling steps combine to form a single apparent rate-controlling step at equilibrium. Here, this apparent rate-controlling step is comprised of an apparent initial state, transition state, and final state, where the apparent final state of the forward reaction is the apparent initial state of the reverse reaction. The apparent rate-controlling step behaves identically to an elementary step reaction at equilibrium with a pseudo-mass-action rate function given by the transition-state-theory (TST) form rate function (Foley and Bhan, 2020) with a stoichiometric number equal to the affinity-averaged stoichiometric number, $\bar{σ}$ The ratio of the forward and reverse TST-form rate functions is identical to the overall thermodynamic equilibrium relation, consistent with the principles of microscopic reversibility and detailed balance at equilibrium.

42 ENGINEERING↗
DOE OSTI2023

Description of the multinucleon transfer mechanism for Ca 48 + Pu 244 and Kr 86 + Pt 198 reactions in a quantal transport approach

Multinucleon transfer (MNT) reactions involving heavy projectile and target combinations stand as a promising method for synthesizing new neutron-rich exotic nuclei, which may not be possible using hot or cold fusion reactions or fragmentation. Exploring the mechanisms behind MNT reactions is essential and it requires a comprehensive theoretical framework that can explain the physical observables in these reactions. This work aims to show that the quantal diffusion approach based on the stochastic mean-field (SMF) theory is capable of explaining the reaction dynamics observed in MNT reactions. Primary product mass distributions in 48 Ca + 244 Pu reaction at E c.m. = 203.2 MeV and 86 Kr + 198 Pt reaction at E c.m. = 324.2 MeV are calculated and compared with the available experimental data. In this work, we utilize the time-dependent Hartree-Fock (TDHF) calculations to analyze the mean-field reaction dynamics computationally in the reactions 48 Ca + 244 Pu and 86 Kr + 198 Pt for a broad range of initial angular momenta. Quantal transport description based on the SMF approach is used to calculate quantal diffusion coefficients and mass variances in 48 Ca + 244 Pu and 86 Kr + 198 Pt systems. The primary products arising from quasifission reactions are described by joint probability distribution in the SMF approach and those arising from fusion-fission are estimated by using the statistical deexcitation code gemini + +. Mean values of charge and mass numbers, scattering angles of the primary reaction products, and the total kinetic energies after the collision are calculated within the TDHF framework for a broad range of initial angular momenta. Throughout all the collisions, drift toward the mass symmetry and large mass dispersion associated with this drift are observed. Here, the calculated primary fragment and mass distributions using the SMF approach successfully explain experimental observations for the 48 Ca + 244 Pu and 86 Kr + 198 Pt systems. The primary mass distributions, mean values of binary products, and mass dispersions are determined and results are compared with the available experimental data. The observed agreement between the experimental data and SMF results highlights the effectiveness of the quantal diffusion mechanism based on the SMF approach, which does not include any adjustable parameters other than standard parameters of Skyrme energy density functional.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Trojan horse method as an indirect approach to study resonant reactions in nuclear astrophysics

The primary goal of the Trojan horse method (THM) is to analyze resonant rearrangement reactions when the density of the resonance levels is low and statistical models cannot be applied. The main difficulty of the analysis is related with the facts that in the final state the THM reaction involves three particles and that the intermediate particle, which is transferred from the Trojan horse particle to the target nucleus to form a resonance state, is virtual. Another difficulty is associated with the Coulomb interaction between the particles, especially, taking into account that the goal of the THM is to study resonant rearrangement reactions at very low energies important for nuclear astrophysics. The exact theory of such reactions with three charged particles is very complicated and is not available. This is why different approximations are used to analyze THM reactions. In this review paper we describe a new approach based on a few-body formalism that provides a solid basis for deriving the THM reaction amplitude taking into account rescattering of the particles in the initial, intermediate and final states of the THM reaction. Since the THM uses a two-step reaction in which the first step is the transfer reaction populating a resonance state, we address the theory of the transfer reactions. The theory is based on the surface-integral approach and R-matrix formalism. We also discuss application of the THM to resonant reactions populating both resonances located on the second energy sheet and subthreshold resonances, which are subthreshold bound states located at negative energies close to thresholds. We consider the application of the THM to determine the astrophysical factors of resonant radiative-capture reactions at energies so low that direct measurements can hardly be performed due to the negligibly small penetrability factor in the entry channel of the reaction. Here, we elucidated the main ideas of the THM and outline necessary conditions to perform the THM experiments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Pickaxe: a Python library for the prediction of novel metabolic reactions

Abstract Background Biochemical reaction prediction tools leverage enzymatic promiscuity rules to generate reaction networks containing novel compounds and reactions. The resulting reaction networks can be used for multiple applications such as designing novel biosynthetic pathways and annotating untargeted metabolomics data. It is vital for these tools to provide a robust, user-friendly method to generate networks for a given application. However, existing tools lack the flexibility to easily generate networks that are tailor-fit for a user’s application due to lack of exhaustive reaction rules, restriction to pre-computed networks, and difficulty in using the software due to lack of documentation. Results Here we present Pickaxe, an open-source, flexible software that provides a user-friendly method to generate novel reaction networks. This software iteratively applies reaction rules to a set of metabolites to generate novel reactions. Users can select rules from the prepackaged JN1224min ruleset, derived from MetaCyc, or define their own custom rules. Additionally, filters are provided which allow for the pruning of a network on-the-fly based on compound and reaction properties. The filters include chemical similarity to target molecules, metabolomics, thermodynamics, and reaction feasibility filters. Example applications are given to highlight the capabilities of Pickaxe: the expansion of common biological databases with novel reactions, the generation of industrially useful chemicals from a yeast metabolome database, and the annotation of untargeted metabolomics peaks from an E. coli dataset. Conclusion Pickaxe predicts novel metabolic reactions and compounds, which can be used for a variety of applications. This software is open-source and available as part of the MINE Database python package ( https://pypi.org/project/minedatabase/ ) or on GitHub ( https://github.com/tyo-nu/MINE-Database ). Documentation and examples can be found on Read the Docs ( https://mine-database.readthedocs.io/en/latest/ ). Through its documentation, pre-packaged features, and customizable nature, Pickaxe allows users to generate novel reaction networks tailored to their application.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Curating a comprehensive set of enzymatic reaction rules for efficient novel biosynthetic pathway design

Enzyme substrate promiscuity has significant implications for metabolic engineering. The ability to predict the space of possible enzymatic side reactions is crucial for elucidating underground metabolic networks in microorganisms, as well as harnessing novel biosynthetic capabilities of enzymes to produce desired chemicals. Reaction rule-based cheminformatics platforms have been implemented to computationally enumerate possible promiscuous reactions, relying on existing knowledge of enzymatic transformations to inform novel reactions. However, past versions of curated reaction rules have been limited by a lack of comprehensiveness in representing all possible transformations, as well as the need to prune rules to enhance computational efficiency in pathway expansion. To this end, we curated a set of 1224 most generalized reaction rules, automatically abstracted from atom-mapped MetaCyc reactions and verified to uniquely cover all common enzymatic transformations. Additionally, we developed a framework to systematically identify and correct redundancies and errors in the curation process, resulting in a minimal, yet comprehensive, rule set. These reaction rules were capable of reproducing more than 85% of all reactions in the KEGG and BRENDA databases, for which a large fraction of reactions is not present in MetaCyc. Our rules exceed all previously published rule sets for which reproduction was possible in this coverage analysis, which allows for the exploration of a larger space of known enzymatic transformations. By leveraging the entire knowledge of possible metabolic reactions through generalized enzymatic reaction rules, we are able to better utilize underground metabolic pathways and accelerate novel biosynthetic pathway design to enable bioproduction towards a wider range of new molecules.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

State-to-State Molecular Reactions in the Ultracold Regime (Final Scientific Report)

We have achieved the proposed goal to experimentally probe the AB + CD and AB + C types of reactions with state-to-state resolution, which we also compared to advanced theoretical calculations to help elucidate the role of quantum mechanics in the processes of bond breakage and formation. Our approach uses reactants that are prepared at ultracold temperatures (< 1µK) such that the quantum effects of translational motion are an important factor. Specific example reactions, including the potassium-rubidium metathesis reaction KRb + KRb → K 2 + Rb 2 as well as the atom exchange reaction Rb + KRb → Rb2 + K, are chosen because the technology of quantum internal and motional state control of these types of molecules is particularly advanced. The results for the entire funding period are fruitful. For the majority of this grant, we have constructed a one-of-the-kind quantum degenerate gas apparatus that integrates ion detection and velocity map imaging capabilities, allowing us to explore the KRb + KRb → K 2 + Rb 2 bimolecular reaction in detail. Specifically, we first verified such a reaction indeed proceed at ultracold temperatures by direct detection of reaction products. We then mapped out the complete product state distribution, which was compared to a state-counting model based on statistical theory. Our results show an overall agreement with the statistical state counting model, but also reveal several deviating state-pairs. An exact quantum calculation for molecule-molecule collisions, that is needed to understand these deviations, is however beyond the current state-of-the-art. Beside scrutinizing the reaction products, we also directly observe the reaction intermediate complex, which was quite a surprise to us. The intermediate complexes are long-lived and can interact with the inferred light that we use to trap the ultracold gas. After molecule-molecule collisions, we then explored the more theoretically tractable Rb + KRb reaction, which is endothermic. Surprisingly, we observed an exceedingly long-lived KRb$^*_2$ collisional complexes, with our experimentally measured complex lifetime deviating from conventional theoretical calculations by five orders of magnitude. This discrepancy has motivated many explorations of possible underlying causes, though no model yet captures this phenomenon completely. In the final year and the work that continues today, we extend upon these atom-molecule collision experiments to explore the origin of the long-lived KRb$^*_2$ complex lifetime and develop means to control the outcome of the reaction complex. The 5-year funded work advanced our understanding of chemical reactions at the lowest possible temperatures and at the same time opened up many new questions that are beyond our initial imaginations.

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

Theoretical studies of chemical reactions related to the formation and growth of polycyclic aromatic hydrocarbons (PAH) and molecular properties of their key intermediates (Final Progress Report)

The formation mechanisms of polycyclic aromatic hydrocarbons, (PAHs) – organic molecules carrying fused benzene rings – are of great interest to scientists and engineers due to their importance in combustion chemistry and astrochemistry. On Earth, PAHs are largely produced in incomplete combustion of fossil fuel and are considered as critical precursors to unwanted soot particles leading to combustion inefficiency and causing air pollution along with detrimental health effects. Simple PAH molecules initially formed in the gas phase, are further involved in a build-up process in combustion flames leading to larger PAH, bowl-shaped nanostructures, fullerenes, and solid-phase species including carbonaceous dust, graphene particles, and soot. In deep space, PAH and their derivatives are potential key intermediates and nucleation sites leading eventually to carbonaceous nanoparticles (“interstellar grains”). Therefore, the understanding of the key processes in the synthesis of PAHs along with their precursors and their degradation mechanisms in combustion systems and in interstellar, circumstellar, and planetary atmospheric environments will provide critical insights into how complex aromatic structures, carbonaceous nanoparticles, and fullerenes are formed and destroyed. Achieving this understanding is an important step in the development of the efficient combustion processes and of the ecofriendly devices with reduced environmental pollution as well as technological strategies for the production of hydrogen and solid carbon through thermal or plasma-assisted pyrolysis of natural gas and biomass. Also, the understanding of the key processes of PAH and soot growth will help in our comprehension of chemical evolution in the universe. Detailed information on the mechanisms and reliable rate constants of the key elementary chemical reactions involved in PAH formation and destruction processes and in inception of soot particles is often missing, with the main deficiencies being the absence of temperature- and pressure-dependent rate constants for the broad range of conditions occurring in various terrestrial and interstellar processes and the lack of data on the reaction products and their branching ratios. Complementary to experimental studies, these gaps in knowledge can be filled by using quantum chemical calculations of reaction potential energy surfaces providing us with accurate energies of reaction products, intermediates, and transition states, revealing the reaction mechanism, and giving the molecular properties required to compute rate constants for relevant reaction steps and product branching ratios using the RRKM-Master Equation (ME) method. Molecular dynamics (MD) simulations can be used in cases when a reaction rate cannot be properly described by statistical theories. During the terminal renewal project period we employed these ab initio/RRKM-ME and MD approaches to complete our studies on several key reactions relevant to the formation/growth of PAH and inception of soot particles including (1) the reaction mechanism and kinetics of the resonance stabilized fulvenallenyl radical with propargyl and C 3 H 4 isomers; (2) the reaction mechanism and kinetics for the C + indene and C 2 + styrene reactions producing naphthyl or azulenyl radicals in low-temperature environments; (3) the MD study of non-equilibrium dimerization of acepyrene and coronene and its radical. The information derived from our theoretical calculations contributed to a better fundamental understanding of the reaction mechanisms and provide missing critical kinetic data to improve combustion models of hydrocarbon fuels and astrochemical models of the growth of carbonaceous molecules and particles in cold molecular clouds, circumstellar envelopes, and planetary atmospheres.

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

Effects of surface species and homogeneous reactions on rates and selectivity in ethane oxidation on oxide catalysts

Selective alkane oxidations on metal oxide catalysts involve complex mechanisms with multiple reactions in series and parallel, different types of reduced and oxidized surface species, and potential contributions from gas-phase reactions. We report kinetics and thermodynamics of elementary steps involved in C 2 H 6 -O 2 reactions on SiO 2 -supported small vanadium oxide domains are determined using density functional theory. These surface reactions together with gas-phase mechanisms are incorporated in kinetic simulations to determine how surface and gaseous reactions interact and contribute to rates and selectivity. The results show that gas-phase reactions within pore volumes in contact with the catalyst contribute significantly to C 2 H 6 activation rates, even at conditions where gas-phase reactions in empty volumes without catalyst are negligible. The majority of C 2 H 6 activations occur on the surface, via H abstraction by vanadium oxo species present at terminal lattice oxygens. The gas-phase activations via H-abstraction by OH radicals also exhibit significant contributions. The reduced centers formed by reactions at vanadium oxo species are re-oxidized rapidly and, therefore, are present in very small concentrations at reaction conditions. The re-oxidation steps lead to the formation of HO 2 radicals and surface peroxo species that are also rapidly consumed and are present in small concentrations. The peroxo species preferentially convert C 2 H 4 to its epoxide product and influence selectivity even at low concentrations. The gas-phase reactions decrease the concentrations of peroxo species and improve selectivity slightly. The effects of reaction conditions and catalyst site density provide further insights into how factors beyond conversions at lattice oxygens influence rates and selectivity in alkane oxidation reactions of significant industrial importance.

42 ENGINEERING↗
DOE OSTI2022

Dissociation-induced depletion of high-energy reactant molecules as a mechanism for pressure-dependent rate constants for bimolecular reactions

In 1922, Lindemann proposed the now-well-known mechanism for pressure-dependent rate constants for unimolecular reactions: reactant molecules with sufficiently high energies dissociate more quickly than collisions can reestablish the Boltzmann distribution of the internal energies of the molecule during its dissociation at low pressures – yielding pressure-dependent rate constants for unimolecular reactions due to the preferential depletion of the high energy states capable of dissociation. In the last century, incredible progress has been made in achieving a far greater understanding of and quantitative predictions for unimolecular and association reactions. In the modern era, pressure-dependent phenomenological rate constants are now nearly universally used to describe the rates of unimolecular and associative reactions in phenomenological kinetic modeling. However, there is a second, more indirect, implication of Lindemann’s mechanism that relates to how these dissociation-induced non-equilibrium distributions impact bimolecular reactions, including non-associative bimolecular reactions – which are generally not considered to have pressure-dependent rate constants. Yet, as we show herein, the same high energy states depleted due to dissociation would otherwise react most rapidly in high-activation-energy bimolecular reactions – yielding a mechanism for pressure-dependent rate constants for bimolecular reactions (including non-associative reactions). Here, we present results from a case study for CH 2 O dissociation, isomerization, and bimolecular reaction with O2 to explore this question. Results from our master equation calculations indicate that the effect of dissociation-induced non-equilibrium distributions on bimolecular reactions can be substantial – even when chemical timescales are well separated from internal energy relaxational timescales (i.e. when the traditional rate constant description would be thought to apply). This effect is found to be more pronounced – and more complex – for bimolecular reactions involving molecular entities whose chemical timescales are merged with the internal energy relaxational timescales. Lastly, we present some ideas for discussion regarding what should be considered as “chemical species” in phenomenological kinetic models.

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

Revealing the role of redox reaction selectivity and mass transfer in current–voltage predictions for ensembles of photocatalysts

Photocatalysts are conceptually simple reaction units where nanoscale semiconductors integrated with catalysts drive a pair of redox reactions on illumination. However, the proximity of reaction sites performing cathodic and anodic reactions poses dire challenges to realize large light-to-fuel conversion efficiencies. In this study, a powerful, yet straightforward, equivalent-circuit detail-balance modeling framework is developed and applied to evaluate the performance of photocatalytic systems featuring multiple light absorbers. Specifically, low bandgap iridium-doped strontium titanate is modeled as a Z-scheme photocatalyst to achieve desirable hydrogen evolution and iron-based redox shuttle oxidation reactions. Our model has unique capabilities to simulate competing redox reactions and address mass-transfer limitations. In a significant departure from state-of-the-art circuit models, our study develops tools to perform load-line analyses by incorporating a net electrochemical load curve that includes both desired and competing redox reactions. Consequently, reaction selectivity is predicted from equivalent circuit models for photocatalytic and photoelectrochemical systems. Our investigation into ensembles comprised of multiple, semi-transparent light absorbers reveals their potential to outperform a single, optically thick light absorber, particularly when operated under mass-transfer-limited conditions. However, this outcome hinges on minimizing mass-transfer rates of select redox species to prevent undesired reactions of hydrogen oxidation and/or redox shuttle reduction. Our findings demonstrate that reaction selectivity can be achieved by tuning asymmetry in redox species mass-transfer even with perfectly symmetric electrocatalytic charge-transfer coefficients. The influences of various kinetic, mass-transfer, and thermodynamic parameters are explored to offer crucial insights for synthesis of the next-generation of photocatalysts and selective coatings, and reactor designs.

25 ENERGY STORAGE↗
DOE OSTI2027

A multi-region approach for the analysis of porous materials and gas-solid reactions using USAXS-SAXS-WAXS: application to CaO carbonation

The reaction rate of gas-solid non-catalytic reactions is typically investigated using reactant conversion data over time and ex-situ measurements of the porous solid reactant textural properties; these data enable the experimental estimation of the initial intrinsic reaction rate, whereas the evolution of the textural properties and the reaction rate over time are evaluated theoretically using reaction models. In this work, a different methodology is presented, based on: a) in-situ time-resolved USAXS-SAXS-WAXS measurements of the solid sample, and b) a multi-region modeling approach; this methodology allows for the estimation of the textural properties and the intrinsic reaction rate at any time during the reaction. Referring to a reaction in which a solid product is obtained from a solid reactant and a gaseous reactant, the porous particle is described as consisting of several distinct regions with different microstructural properties, and both a gas-solid reaction model and a SAXS model are derived and applied to the carbonation of porous CaO. The proposed reaction model highlights the role of the inaccessible reactant and accurately predicts the solid reactant conversion versus time profile. The SAXS model accurately predicts the measured linear trends of the Porod invariant and of the pre-factor of the power law scattering profile versus CaO mass fraction; in the absence of chemical reactions, the proposed equations extend classical SAXS theory to porous materials containing macropores and nonporous solid phases in addition to standard nanoscale inhomogeneous regions.

CO2 solid sorbents↗
DOE OSTI2023

Concepts Relevant for the Kinetic Analysis of Reversible Reaction Systems

The net rate of a reversible chemical reaction is the difference between unidirectional rates of traversal along forward and reverse reaction paths. In a multistep reaction sequence, the forward and reverse trajectories, in general, are not the microscopic reverse of one another; rather, each unidirectional route is comprised of distinct rate-controlling steps, intermediates, and transition states. Consequently, traditional descriptors of rate (e.g., reaction orders) do not reflect intrinsic kinetic information but instead conflate unidirectional contributions determined by (i) the microscopic occurrence of forward/reverse reactions (i.e., unidirectional kinetics) and (ii) the reversibility of reaction (i.e., nonequilibrium thermodynamics). This review aims to provide a comprehensive resource of analytical and conceptual tools which deconvolute the contributions of reaction kinetics and thermodynamics to disambiguate unidirectional reaction trajectories and precisely identify rate- and reversibility-controlling molecular species and steps in reversible reaction systems. The extrication of mechanistic and kinetic information from bidirectional reactions is accomplished through equation-based formalisms (e.g., De Donder relations) grounded in principles of thermodynamics and interpreted in the context of theories of chemical kinetics developed in the past 25 years. Here, the aggregate of mathematical formalisms detailed herein is general to thermochemical and electrochemical reactions and encapsulates a diverse body of scientific literature encompassing chemical physics, thermodynamics, chemical kinetics, catalysis, and kinetic modeling.

42 ENGINEERING↗
DOE OSTI2020

Reactants, products, and transition states of elementary chemical reactions based on quantum chemistry

Reaction times, activation energies, branching ratios, yields, and many other quantitative attributes are important for precise organic syntheses and generating detailed reaction mechanisms. Often, it would be useful to be able to classify proposed reactions as fast or slow. However, quantitative chemical reaction data, especially for atom-mapped reactions, are difficult to find in existing databases. Therefore, we used automated potential energy surface exploration to generate 12,000 organic reactions involving H, C, N, and O atoms calculated at the ω B97X-D3/def2-TZVP quantum chemistry level. We report the results of geometry optimizations and frequency calculations for reactants, products, and transition states of all reactions. Additionally, we extracted atom-mapped reaction SMILES, activation energies, and enthalpies of reaction. We believe that this data will accelerate progress in automated methods for organic synthesis and reaction mechanism generation—for example, by enabling the development of novel machine learning models for quantitative reaction prediction.

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

Probing the Potential Energy Profile of the I + (H 2 O) 3 → HI + (H 2 O) 2 OH Forward and Reverse Reactions: High Level CCSD(T) Studies with Spin-Orbit Coupling Included

Three different pathways for the atomic iodine plus water trimer reaction I + (H 2 O) 3 → HI + (H 2 O) 2 OH were preliminarily examined by the DFT-MPW1K method. Related to previous predictions for the F/Cl/Br + (H 2 O) 3 reactions, three pathways for the I + (H 2 O) 3 reaction are linked in terms of geometry and energetics. To legitimize the results, the “gold standard” CCSD(T) method was employed to investigate the lowest-lying pathway with the correlation-consistent polarized valence basis set up to cc-pVQZ(-PP). According to the CCSD(T)/cc-pVQZ(-PP)//CCSD(T)/cc-pVTZ(-PP) results, the I + (H 2 O) 3 → HI + (H 2 O) 2 OH reaction is predicted to be endothermic by 47.0 kcal mol -1 . The submerged transition state is predicted to lie 43.7 kcal mol -1 above the separated reactants. The I···(H 2 O) 3 entrance complex lies below the separated reactants by 4.1 kcal mol -1 , and spin-orbit coupling has a significant impact on this dissociation energy. The HI···(H 2 O) 2 OH exit complex is bound by 4.3 kcal mol -1 in relation to the separated products. Compared with simpler I + (H 2 O) 2 and I + H 2 O reactions, the I + (H 2 O) 3 reaction is energetically between them in general. It is speculated that the reaction between the iodine atom and the larger water clusters may be energetically analogous to the I + (H 2 O) 3 reaction. The iodine reaction I + (H 2 O) 3 is connected with the analogous valence isoelectronic bromine/chlorine reactions Br/Cl + (H 2 O) 3 but much different from the F + (H 2 O) 3 reaction. Significant difference with other halogen systems, especially for barrier heights, are seen for the iodine systems.

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

Inertia-induced mixing and reaction maximization in laminar porous media flows

Solute transport and biogeochemical reactions in porous and fractured media flows are controlled by mixing, as are subsurface engineering operations such as contaminant remediation, geothermal energy production, and carbon sequestration. Porous media flows are generally regarded as slow, so the effects of fluid inertia on mixing and reaction are typically ignored. Here, we demonstrate through microfluidic experiments and numerical simulations of mixing-induced reaction that inertial recirculating flows readily emerge in laminar porous media flows and dramatically alter mixing and reaction dynamics. An optimal Reynolds number that maximizes the reaction rate is observed for individual pore throats of different sizes. This reaction maximization is attributed to the effects of recirculation flows on reactant availability, mixing, and reaction completion, which depend on the topology of recirculation relative to the boundary of the reactants or mixing interface. Recirculation enhances mixing and reactant availability, but a further increase in flow velocity reduces the residence time in recirculation, leading to a decrease in reaction rate. The reaction maximization is also confirmed in a flow channel with grain inclusions and randomized porous media. Interestingly, the domain-wide reaction rate shows a dramatic increase with increasing Re in the randomized porous media case. This is because fluid inertia induces complex three-dimensional flows in randomized porous media, which significantly increases transverse spreading and mixing. This study shows how inertial flows control reaction dynamics at the pore scale and beyond, thus having major implications for a wide range of environmental systems.

Chen, Michael A. (ORCID:0000000294614710)↗
DOE OSTI2020

Improved Performance of Cu(InGa)(SeS) 2 PV Modules Using the Reaction of Metal Precursors. Final Report

This project “Improved Performance of Cu(InGa)(SeS) 2 PV Modules using the Reaction of Metal Precursors” was a partnership led by the Institute of Energy Conversion (IEC) at the University of Delaware with Columbia University and the Molecular Foundry at the Lawrence Berkeley National Laboratory. The aim was to develop pathways to improve Cu(InGa)(SeS) 2 (CIGSS) thin film photovoltaic modules using processes compatible with low manufacturing cost. The CIGSS approach investigated was a two-step process including deposition of metal precursor films following by reaction in hydride gases utilizing IEC’s novel reactor. The process was similar to that under commercial development by the project’s industry partner Stion. When Stion went out of business mid-project the focus changed to a rapid thermal process considered more commercially viable. Approaches to improve the performance of solar cells using the reacted films focused on two material innovations. First, the overall Ga content was increased to increase the operating voltage, which is desirable for scale-up to commercial modules. Second, the processing and performance advantages arising from Ag alloying were investigated. Advanced characterization guided process and material development including control of relative composition gradients. Research on the formation of Cu-Ga-In metal precursors utilized sputtering deposition which is normally used in commercial applications. The work resulted in processes for deposition of precursor stacks with increased relative Ga content and effects of deposition parameters on morphology and phase composition were established. It was shown that the metal precursor films have comparable phase composition and morphology so subsequent reaction follows from the same starting point. The addition of Ag to the metal precursors gave more uniform morphology and improved adhesion of reacted films which enable higher reaction temperature for faster processing. A novel outcome was the discovery of a previously undocumented material phase in sputter-deposited and evaporated Ag-Cu-In-Ga thin films. Hydride gas reaction processes including time-temperature-concentration profiles were developed for different precursor compositions. This enables control of composition profiles to engineer through-film gradients for solar cell optimization with characterization and simulations used to correlate measured film composition profiles to measurements of devices. In particular, the gradient of sulfur at the front of the CIGSS film was found to be critical. The simulations guided process development leading to improved reproducibility of devices improved performance with higher Ga content and higher voltage. With Ag-alloyed precursors, the reaction pathways leading were determined. A significant finding was that Ag-alloying increases the reaction rate to completely convert precursor films to the final chalcopyrite which could enable reduced reaction time to benefit manufacturability. To maintain potential commercial viability, the process under investigation was refocused to a rapid thermal process that could potentially be incorporated into an in-line process for manufacturing. Precursors with different composition were capped with an extra selenium layer and reacted in hydrogen sulfide 5-15 minutes, compared to typically 2 hours in the previous multi-step batch process. Critical RTP parameters were identified to control the reaction. Further optimization would be needed for high efficiency solar cells but pathways to high quality devices with further optimization and improved heating uniformity were developed. The project also developed new optoelectronic characterization approaches with a focus on development and application of spatial- and time-resolved photoluminescence and a custom mapping photoluminescence microscope built. It was shown how critical electronic transport properties strongly depend on the chemical composition of the material and that a wide range of samples show inhomogeneity on a length scale larger than the grains in the films. Additionally, two-photon excitation capability was developed to distinguish bulk vs surface losses. The project advances the state-of-the -art for precursor reaction processes in several ways that could impact manufacturing. This includes validation of approaches to increase voltage and establishment of model-guided control to form optimal composition profiles. The application of process control approaches with knowledge of phase formation and reaction pathways can be critically valuable in designing a large-scale process.

14 SOLAR ENERGY↗
DOE OSTI2025

Molecular Mass Growth Processes to Polycyclic Aromatic Hydrocarbons through Radical–Radical Reactions Exploiting Photoionization Reflectron Time-of-Flight Mass Spectrometry

Polycyclic aromatic hydrocarbons (PAHs) represent critical building blocks in molecular mass growth processes to carbonaceous nanoparticles, referred to as interstellar and circumstellar grains along with soot particles in astrophysical environments and combustion systems, respectively. Recent advancements on elucidating elementary steps to PAHs have utilized reactions of aromatic radicals, resonantly stabilized free radicals, and aliphatic radicals with closed shell hydrocarbons. However, the role of radical–radical reactions (RRRs) leading to PAHs has remained largely unexplored on the molecular level due to preceding experimental challenges in producing sufficiently high number densities of radical reactants for isomer-selective detection of products from bimolecular and termolecular reactions. This Account offers the latest developments in our knowledge on the mechanisms and pathways to PAHs via RRRs probed in a chemical microreactor at temperatures as high as 1600 K. Product preservation in a molecular beam coupled with synchrotron vacuum ultraviolet photoionization reflectron time-of-flight mass spectrometry and photoelectron photoion coincidence spectroscopy enabled isomer-selective detection of PAHs of up to three rings by their photoionization efficiency curves, which were fit with a linear combination of reference curves for identification. Experiments were combined with computational fluid dynamics modeling of the physicochemical processes in the microreactor, as well as high-level electronic structure calculations to reveal the reaction pathways of each system. Six distinct reaction mechanisms were discovered in this work: propargyl addition─benzannulation (PABA), methyl addition─ring expansion (MARE), cyclopentadienyl addition─naphthylization (CPAN), fulvenallenyl addition─cyclization─aromatization (FACA), benzyl addition─aromatization (BAA), and phenyl addition─pentacyclization (PAP). By systematically varying the number of carbon atoms in the radical reactants, molecular mass growth processes involving reactions between radicals with odd numbers of carbon atoms access aromatics carrying one, two, or three six-membered rings, whereas reactions between even- and odd-carbon-numbered radicals produce aromatics combining five- and six-membered rings. Our investigations reveal unconventional cycloadditions on excited state triplet surfaces, additions of radicals to low spin density carbon-centered radicals, spiroaromatic and fulvene-type intermediates, and highly strained bicyclic reaction intermediates, challenging current perceptions of PAH molecular mass growth processes. All of the listed mechanisms, except for FACA, feature endoergic reactions or barriers which lie above the separated reactants and therefore might be central to circumstellar environments of carbon-rich stars and planetary nebulae as their descendants, but they play no role in the gas phase of cold molecular clouds where temperatures as low as 10 K dominate. Altogether, this work provides detailed reaction mechanisms of PAH growth processes, advancing our knowledge of the chemistry of carbonaceous matter in the universe.

Addition reactions↗