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

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

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

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

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

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

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

Phase-field modeling of solid-state metathesis reactions with the charge neutrality constraint

In this work, we present a phase-field model that captures the evolution of ionic concentrations and phase fractions during solid-state metathesis (SSM) reactions where diffusion limits the rate of transformation. The evolution of the mole fraction of each ion is obtained via governing equations that describe the reduction of a free energy, which includes an energy landscape with local minima located at compositions corresponding to stable products. We utilized two Lagrange multipliers to impose constraints of electroneutrality as well as on the sum of mole fractions, which were then eliminated to derive set of two partial differential equations that describe the dynamics of the mole fraction evolution. From these governing equations, the expressions for effective mobilities for the cations and the anions were obtained. We first study the effect of mobilities of ions on the reaction kinetics, using a simple model considering the ions with an identical absolute value of charge numbers. The simulation results show that the overall characteristic mobility, defined as the sum of the two effective ionic mobilities, provides an excellent measure of the rate at which reaction progresses and that the ratio of the effective mobilities of the anions and the cations signifies the manner by which the reaction progresses. We then generalize the model to consider ions with different charge numbers and tuned the mobility of ions based on their diffusion coefficients reported in the literature and experimental data from a thin-film experiment for the synthesis of FeS 2 to demonstrate the capability of the model to predict the phase evolution during SSM reactions. In particular, the simulation predicts nonplanar phase evolution, which is recently observed in thin-film reactions for the synthesis of FeS 2 via transmission electron microscopy. The approach can serve as a basis for models for phase transformations in other multiphase ionic mixtures, such as in all-solid-state batteries and in ionic liquids.

36 MATERIALS SCIENCE↗
DOE OSTI2023

The P( 4 S) + NH( 3 Σ – ) and N( 4 S) + PH( 3 Σ – )reactions as sources of interstellar phosphorus nitride

Phosphorus nitride (PN) is believed to be one of the major reservoirs of phosphorus in the interstellar medium (ISM). For this reason, understanding which reactions produce PN in space and predicting their rate coefficients is important for modelling the relative abundances of P-bearing species and clarifying the role of phosphorus in astrochemistry. Here, in this work, we explore the potential energy surfaces of the $\textrm{P}(^4\textrm{S}) + \textrm{NH}(^3\Sigma^-)$ and $\textrm{N}(^4\textrm{S}) + \textrm{PH}(^3\Sigma^-)$ reactions and the formation of $\textrm{H}(^2\textrm{S}) + \textrm{PN}(^1\Sigma^+)$ through high accuracy ab initio calculations and the variable reaction coordinate transition state theory (VRC-TST). We found that both reactions proceed without an activation barrier and with similar rate coefficients that can be described by a modified Arrhenius equation ( $k(T)=\alpha\!\left( T/300 \right)^{\beta} \exp\!{(\!-\!\gamma/T)})$ with $\alpha=0.93\times 10^{-10}\rm cm^3\,s^{-1}$ , $\beta=-0.18$ and $\gamma=0.24\, \rm K$ for the $\textrm{P} + \textrm{NH} \longrightarrow \textrm{H} + \textrm{PN}$ reaction and $\alpha=0.88\times 10^{-10}\rm cm^3\,s^{-1}$ , $\beta=-0.18$ and $\gamma=1.01\, \rm K$ for the $\textrm{N} + \textrm{PH} \longrightarrow \textrm{H} + \textrm{PN}$ one. Both reactions are expected to be relevant for modelling PN abundances even in the cold environments of the ISM. Given the abundance of hydrogen in space, we have also predicted rate coefficients for the destruction of PN via H + PN collisions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

A Combined Crossed Molecular Beam and Theoretical Investigation of the Elementary Reaction of Tricarbon (C 3 (X 1 Σ g + )) with Diacetylene (C 4 H 2 (X 1 Σ g + )): Gas Phase Formation of the Heptatriynylidyne Radical ( l -C 7 H(X 2 Π))

An elucidation of the underlying formation pathways to acyclic hydrocarbons such as polyynes (C n H 2 ), cumulenes (C n H 2 ), and linear resonantly stabilized linear radicals (l-C n H) is indispensable to understand the hydrocarbon chemistry in extreme low and high temperature environments. In this study, we exploited the crossed molecular beam technique to investigate the reaction of tricarbon C 3 (X 1 Σ g + ) with diacetylene (butadiyne; HCCCCH; X 1 Σ g + ) at a collision energy of 47 ± 1 kJ mol ⁻1 . The experimental data were merged with ab initio calculations of the singlet C 7 H 2 potential energy surface (PES) revealing that the reaction is initiated via the formation of an initial van der Waals reactant complex in the entrance channel. Subsequent rearrangements lead to various carbene-type and cyclic intermediates via ring-opening, ring-closure, and hydrogen migration processes eventually forming acyclic C 7 H 2 isomers prior to their barrierless unimolecular decomposition to the most stable linear isomer, heptatriynylidyne (C 7 H, X 2 Π) in an overall endoergic reaction (+57 kJ mol ⁻1 ). The reaction exhibits strong similarities to the tricarbon – acetylene (C 3 – C 2 H 2 ). Furthermore, the significant energy threshold suggests that the tricarbon reaction with (poly)acetylenes forming resonantly stabilized linear radicals are open in high-temperature environments such as combustion flames and circumstellar envelopes of carbon stars and planetary nebulae as their descendants; however, these reactions are closed in low-temperature environments as in cold molecular clouds and hydrocarbon-rich atmospheres of planets and their moons such as in Titan.

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

Deciphering a Reaction Network for the Switchable Production of Tetrahydroquinoline or Quinoline with MOF-Supported Pd Tandem Catalysts

A mechanistic study of heterogeneous tandem catalytic systems is crucial for understanding and improving catalyst activity and selectivity but remains challenging. Here, we demonstrate that a thorough mechanistic study of a multistep reaction can guide us to the controllable selective synthesis of phenyltetrahydroquinoline or phenylquinoline with easily accessible precursors. The one-pot production can be achieved, catalyzed by a well-defined, bifunctional metal–organic framework-supported Pd nanoparticles, with only water as the side product. Our mechanistic study identifies six transient intermediates and ten transformation steps from the operando magic angle spinning nuclear magnetic resonance study under 27.6 bar H 2 . In particular, reactive intermediate 2-phenyl-3,4-dihydroquinoline cannot be observed with conventional chromatographic techniques but is found to reach the maximal concentration of 0.11 mol L –1 under the operando condition. The most probable reaction network is further deduced based on the kinetic information of reaction species, obtained from both operando and ex situ reaction studies. This deep understanding of the complex reaction network enables the kinetic control of the conversions of key intermediate, 2-phenyl-3,4-dihydroquinoline, with the addition of a homogeneous co-catalyst, allowing the selective production of tetrahydroquinoline or quinoline on demand. Finally, the demonstrated methods in this work open up new avenues toward efficient modulation of reactions with a complex network to achieve desired selectivities.

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

Oxygen Atom Transfer Reactions of Colloidal Metal Oxide Nanoparticles

Redox transformations at metal oxide (MO x )/solution interfaces are broadly important, and oxygen atom transfer (OAT) is one of the simplest and most fundamental examples of such reactivity. OAT is a two-electron transfer process, well-known in gas/solid reactions and catalysis. However, OAT is rarely directly observed at oxide/water interfaces, whose redox reactions are typically proposed to occur in one-electron steps. Reported here are stoichiometric OAT reactions of organic molecules with aqueous colloidal titanium dioxide and iridium oxide nanoparticles (TiO 2 and IrO x NPs). Me 2 SO (DMSO) oxidizes reduced TiO 2 NPs with the formation of Me 2 S, and IrO x NPs transfer O atoms to a water-soluble phosphine and a thioether. The reaction stoichiometries were established and the chemical mechanisms were probed using typical solution spectroscopic techniques, exploiting the high surface areas and transparency of the colloids. Furthermore, these OAT reactions, including a catalytic example, utilize the ability of the individual NPs to accumulate many electrons and/or holes. Observing OAT reactions of two different materials, in opposite directions, is a step toward harnessing oxide nanoparticles for valuable multi-electron and multi-hole transformations.

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

Polaritonic Bright and Dark States Collectively Affect the Reactivity of a Hydrolysis Reaction

Vibrational strong coupling (VSC) has emerged as a means for modifying chemical reactivity. Despite the intriguing discoveries and progresses in the field, the precise mechanisms that govern polaritonic chemistry still deserve further interrogation. Herein, we use the hydrolysis of ammonia borane in D 2 O as an exemplary reaction and systematically investigate the influence of VSC on its reactivity. Experimental evidence of the coexistence of a resonant effect and reaction acceleration is observed in this system. In particular, we find that when the O-D stretching mode of D 2 O is strongly coupled to a cavity mode, reaction acceleration is observed. The reaction rate acceleration factor, mu, is consistently observed to be dependent on the coupling conditions between the vibrational and cavity modes, and reaches a minimum at zero mode detuning, suggesting that a resonant effect is likely in play. In addition, we find that mu decreases with an increasing Rabi splitting. Based on these experimental findings, we propose that the overall influence of VSC on this reaction is likely determined collectively by the polaritonic bright and dark states. In conclusion, these findings could help shed new light on the intricate effects of VSC on ground-state reaction landscapes.

Rabi splitting↗
DOE OSTI2021

Understanding Reaction Networks through Controlled Approach to Equilibrium Experiments Using Transient Methods

We report a combined experimental/theoretical approach to studying heterogeneous gas/solid catalytic processes using low-pressure pulse response experiments achieving a controlled approach to equilibrium that combined with quantum mechanics (QM)-based computational analysis provides information needed to reconstruct the role of the different surface reaction steps. We demonstrate this approach using model catalysts for ammonia synthesis/decomposition. Polycrystalline iron and cobalt are studied via low-pressure TAP (temporal analysis of products) pulse response, with the results interpreted through reaction free energies calculated using QM on Fe-BCC(110), Fe-BCC(111), and Co-FCC(111) facets. In TAP experiments, simultaneous pulsing of ammonia and deuterium creates a condition where the participation of reactants and products can be distinguished in both forward and reverse reaction steps. This establishes a balance between competitive reactions for D* surface species that is used to observe the influence of steps leading to nitrogen formation as the nitrogen product remains far from equilibrium. Here, the approach to equilibrium is further controlled by introducing delay timing between NH 3 and D 2 which allows time for surface reactions to evolve before being driven in the reverse direction from the gas phase. The resulting isotopic product distributions for NH 2 D, NHD 2 , and HD at different temperatures and delay times and NH 3 /D 2 pulsing order reveal the role of the N 2 formation barrier in controlling the surface concentration of NH x * species, as well as providing information on the surface lifetimes of key reaction intermediates. Conclusions derived for monometallic materials are used to interpret experimental results on a more complex and active CoFe bimetallic catalyst.

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

Mode Specificity in the OH + HO 2 → H 2 O + O 2 Reaction: Enhancement of Reactivity by Exciting a Spectator Mode

A reaction usually involves a few active modes while the other modes are largely preserved throughout the reaction as spectators. Excitation of an active mode is expected to promote the reaction, but depositing energy in a spectator mode typically has no effect, because of the differing ability for energy flow to the reaction coordinate. Here, we report a surprising case of mode specificity in a key radical–radical reaction OH + HO 2 → H 2 O + O 2 , where such canonical expectations fail to hold. Despite its spectator nature, the vibrational excitation of the OH reactant is shown at low collision energies to enhance the reactivity significantly. This unique effect can be attributed to the increased attraction with HO 2 due to the larger dipole of the stretched OH. At low collision energies, the stronger attraction increases the chance of capturing the reactants to form a hydrogen-bonded complex, thus of passing through the submerged barrier. The novel mechanism differs from the conventional vibrational enhancement via coupling to the reaction coordinate at the transition state, enriching our understanding of mode specificity in chemistry.

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

Effects of diffusion barriers on reaction wave stability in Co/Al reactive multilayers

Bimetallic, reactive multilayers are uniformly structured materials composed of alternating sputter-deposited layers that may be ignited to produce self-propagating mixing and formation reactions. These nanolaminates are most commonly used as rapid-release heat sources. The specific chemical composition at each metal/metal interface determines the rate of mass transport in a mixing and formation reaction. The inclusion of engineered diffusion barriers at each interface will not only inhibit solid-state mixing but also may impede the self-propagating reactions by introducing instabilities to wavefront morphology. Here, this work examines the effect of adding diffusion barriers on the propagation of reaction waves in Co/Al multilayers. The Co/Al system has been shown to exhibit a reaction propagation instability that is dependent on the bilayer thickness, which allows for the occurrence of unstable modes in otherwise stable designs from the inclusion of diffusion barriers. Based on the known stability criteria in the Co/Al multilayer system, the way in which the inclusion of diffusion barriers changes a multilayer's heat of reaction, thermal conductivity, and material mixing mechanisms can be determined. These factors, in aggregate, lead to changes in the wavefront velocity and stability.

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

A guided ion beam investigation of UO 2 + thermodynamics and f orbital participation: Reactions of U + + CO 2 , UO + + O 2 , and UO + + CO

A guided ion beam tandem mass spectrometer was employed to study the reactions of U + + CO 2 , UO + + O 2 , and the reverse of the former, UO + + CO. Reaction cross sections as a function of kinetic energy over about a three order of magnitude range were studied for all systems. The reaction of U + + CO 2 proceeds to form UO + + CO with an efficiency of 118% ± 24% as well as generating UO 2 + + C and UCO + + O. The reaction of UO + + O 2 forms UO 2 + in an exothermic, barrierless process and also results in the collision-induced dissociation of UO + to yield U + . In the UO + + CO reaction, the formation of UO 2 + in an endothermic process is the dominant reaction, but minor products of UCO + + O and U + + (O + CO) are also observed. Analysis of the kinetic energy dependences observed provides the bond energies, D 0 (U + –O) = 7.98 ± 0.22 and 8.05 ± 0.14 eV, D 0 (U + –CO) = 0.73 ± 0.13 eV, and D 0 (OU + –O) = 7.56 ± 0.12 eV. The values obtained for D 0 (U + –O) and D 0 (OU + –O) agree well with the previously reported literature values. To our knowledge, this is the first experimental measurement of D 0 (U + –CO). Furthermore, an analysis of the oxide bond energies shows that participation of 5f orbitals leads to a substantial increase in the thermodynamic stability of UO 2 + relative to ThO 2 + and especially transition metal dioxide cations.

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

12 C(p,p') 12 C Reaction ( E p =19.5–30 MeV) for Active Interrogation of Special Nuclear Material

Passive detection of special nuclear material (SNM) is challenging due to its inherently low rate of spontaneous emission of penetrating radiation, the relative ease of shielding, and the fluctuating and frequently overwhelming background. Active interrogation, the use of external radiation to increase the emission rate of characteristic radiation from SNM, has long been considered to be a promising method to overcome those challenges. Current active-interrogation systems that incorporate radiography tend to use bremsstrahlung beams, which can deliver high radiation doses. Low-energy ion-driven nuclear reactions that produce multiple monoenergetic photons may be used as an alternative. The 12 C ( p , p ' ) 12 C reaction is one such reaction that could produce large yields of highly penetrating 4.4- and 15.1-MeV gamma rays. This reaction does not directly produce neutrons below the approximately 19.7 MeV threshold, and the 15.1-MeV gamma-ray line is well matched to the photofission cross section of 235 U and 238 U . In this article, we report the measurements of thick-target gamma-ray yields at 4.4 and 15.1 MeV from the 12 C ( p , p ' ) 12 C reaction at proton energies of 19.5, 25, and 30 MeV. Measurements are made with two 3 -in. EJ-309 cylindrical liquid scintillation detectors and thermoluminescent dosimeters placed at 0 ° and 90 ° , with an additional 1.5 -in. Na I ( Tl ) cylindrical scintillation detector at 0 ° . We estimate the highest yields of the 4.4- and 15.1-MeV gamma rays of 1.65 × 10 10 and 4.47 × 10 8 sr - 1 μ C - 1 at a proton energy of 30 MeV, respectively. The yields in all experimental configurations are greater than in a comparable deuteron-driven reaction that produces the same gamma-ray energies— 11 B ( d , n γ ) 12 C . However, a significant increase of the neutron radiation dose accompanies the proton energy increase from 19.5 to 30 MeV.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2025

Deuteron-induced reactions on natural Zr from threshold to 50 MeV: production of 86g Y

Two stacks of thin Zr foils were irradiated with 30 and 50 MeV deuterons, respectively, using the Lawrence Berkeley National Laboratory 88-Inch Cyclotron, and 19 excitation functions for nat Zr(d,x) reactions were measured over a beam energy range of 6.3–47.64 MeV, where the independent cross sections for nat Zr(d,x) 88 Nb and nat Zr(d,x) 86m,g Y were measured for the first time. The well-characterized nat Fe(d,x) 56 Co, nat Ni(d,x) 56 Co, nat Ni(d,x) 58 Co, nat Ni(d,x) 61 Cu, nat Ti(d,x) 46 Sc and nat Ti(d,x) 48 V monitor reactions were used to determine the deuteron beam current throughout the stacks. All cross sections were determined using High Purity Germanium (HPGe) detector γ-ray spectroscopy. A variance minimization technique was employed to simultaneously constrain the deuteron beam currents with multiple monitor reactions, thus reducing systematic uncertainties. An additional 16 channels are reported for reactions on the nickel, titanium, and iron monitor foils, leading to a total of 35 excitation functions, with seven reaction channels reported for the first time in this work. The measured excitation functions are compared to calculations provided by the reaction modeling codes TALYS – 2.0, ALICE – 2020, CoH – 3.5.3 and EMPIRE – 3.2.3, as well as the TENDL – 2023 data library. The degree of agreement between theory and experiments is discussed. The possible production of the important PET radionuclide 86g Y via the nat Zr(d,x) route was critically examined. The physical yields for nat Zr(d,x) 86 Y and other yttrium isotopes produced were calculated and compared to other production pathways. Due to high-level of associated radionuclide impurities, this route cannot deliver 86g Y suitable for medical applications.

86gY↗
DOE OSTI2020

HMX detonation calculation using multi-reaction chain

Detonation of HMX has been simulated using a multi-reaction chain and species equation-of-states. A reacting flow code capable of handling arbitrary number of reactions and species has been used. Computational results are compared to the results produced by a pressure-based single reaction model (AJEX) and a temperature-based single-reaction model. Results show that the multi-reaction approach is capable of producing detonation calculations comparable to the AJEX model calculations, demonstrating that the multi-reaction approach possesses the properties required for the detonation modeling, a prerequisite of deflagration-to-detonation transition modeling.

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

How Well Can Quantum Embedding Method Predict the Reaction Profiles for Hydrogenation of Small Li Clusters?

Quantum computing leverages the principles of quantum mechanics in novel ways to tackle complex chemistry problems that cannot be accurately addressed using traditional quantum chemistry methods. However, the high computational cost and available number of physical qubits with high fidelity limit its application to small chemical systems. This work employed a quantum-classical framework which features a quantum active space-embedding approach to perform simulations of chemical reactions that require up to 14 qubits. This framework was applied to prototypical example metal hydrogenation reactions: the coupling between hydrogen and Li 2 , Li 3 , and Li 4 clusters. Particular attention was paid to the computation of barriers and reaction energies. The predicted reaction profiles compare well with advanced classical quantum chemistry methods, demonstrating the potential of the quantum embedding algorithm to map out reaction profiles of realistic gas-phase chemical reactions to ascertain qualitative energetic trends. Additionally, the predicted potential energy curves provide a benchmark to compare against both current and future quantum embedding approaches.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Thermonuclear Reaction Rates and Primordial Nucleosynthesis

Assuming the best numerical value for the cosmic baryonic density and the existence of three neutrino flavors, standard Big Bang nucleosynthesis is a parameter-free model. It is important to assess if the observed primordial abundances can be reproduced by simulations. Numerous studies have shown that the simulations overpredict the primordial 7 Li abundance by a factor of ≈3 compared to the observations. The discrepancy may be caused by unknown systematics in 7 Li observations, poorly understood depletion of lithium in stars, errors in thermonuclear rates that take part in the lithium and beryllium synthesis, or physics beyond the standard model. Here, we focus on the likelihood of a nuclear physics solution. The status of the key nuclear reaction rates is summarized. Big Bang nucleosynthesis simulations are performed with the most recent reaction rates, and the uncertainties of the predicted abundances are established using a Monte Carlo technique. Correlations between abundances and reaction rates are investigated based on the metric of mutual information. The rates of four reactions impact the primordial 7 Li abundance: 3 He(α,γ) 7 Be, d(p,γ) 3 He, 7 Be(d,p)2α, and 7 Be(n,p) 7 Li. We employ a genetic algorithm to search for simultaneous rate changes in these four reactions that may account for all observed primordial abundances. Furthermore, when the search is performed for reaction rate ranges that are much wider than recently reported uncertainties, no acceptable solutions are found. Based on the currently available evidence, we conclude that it is highly unlikely for the cosmological lithium problem to have a nuclear physics solution.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

The PFLOTRAN Reaction Sandbox

Abstract. As modern reactive transport simulators evolve to accommodate the demands of a user community, researchers need a platform for prototyping new biogeochemical processes, many of which are niche and specific to laboratory or field experiments. The PFLOTRAN Reaction Sandbox leverages modern, object-oriented Fortran in an attempt to provide such an environment within an existing reactive transport simulator. This work describes the PFLOTRAN Reaction Sandbox concept and implementation through several illustrative examples. Reaction Sandbox Biodegradation Hill customizes the existing microbially mediated biodegradation reaction formulation within PFLOTRAN to better match empirical data. Reaction Sandbox Simple provides an isolated environment for testing numerous preconfigured kinetic rate expressions and developing user intuition. Reaction Sandbox Example serves as a template for creating new sandboxes within PFLOTRAN.

58 GEOSCIENCES↗