Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reaction rates”

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 55 records · Page 3

Impact of Thermonuclear Reaction Rate Uncertainties on the Identification of Presolar Grains from Classical Novae

Approximately 30%–40% of classical novae generate dust between 20 and 100 days following the eruption. However, there has yet to be a definitive identification of presolar stardust grains originating from classical novae. While multiple studies have suggested a nova origin for specific grains, aligning simultaneously all measured isotopic ratios of a specific grain with those predicted from simulations remains challenging. Using Monte Carlo simulations, this work investigates how uncertainties in thermonuclear reaction rates influence the isotopic ratios predicted in simulations of classical novae, specifically impacting the identification of presolar grains. In particular, we address two questions: (i) What is the impact of uncertainties in reaction rates on the range of isotopic ratios predicted by classical nova simulations? (ii) Which reaction rate uncertainties most significantly influence the predicted abundance ratios in presolar grains? Our results show that current reaction rate uncertainties affect the isotopic ratios of 12 C/ 13 C, 14 N/ 15 N, 16 O/ 17 O, 16 O/ 18 O, 24 Mg/ 25 Mg, 24 Mg/ 26 Mg, 26 Al/ 27 Al, and 28 Si/ 29 Si by less than 20% in either carbon–oxygen or oxygen–neon (ONe) novae, especially when considering the mixing of matter throughout the entire envelope. However, the isotopic ratios of 28 Si/ 30 Si, 32 S/ 33 S, and 32 S/ 34 S in ONe novae are exceptions: their variability greatly exceeds a factor of 2 due to the uncertainties in the reaction rates of 30 P(p,γ) 31 S, 33 S(p,γ) 34 Cl, and 34 S(p,γ) 35 Cl, respectively. These results highlight the significant influence of specific reaction rates on the predicted abundance ratios and underscore the necessity for accurate nuclear measurements to reduce these uncertainties.

Classical novae↗

Uncertainties in the production of iron-group nuclides in core-collapse supernovae from Monte Carlo variations of reaction rates

Core-collapse supernovae, occurring at the end of massive star evolution, produce heavy elements, including those in the iron peak. Although the explosion mechanism is not yet fully understood, theoretical models can reproduce optical observations and observed elemental abundances. However, many nuclear reaction rates involved in explosive nucleosynthesis have large uncertainties, impacting the reliability of abundance predictions. To address this, we have previously developed a Monte Carlo-based nucleosynthesis code that accounts for reaction rate uncertainties and has been applied to nucleosynthesis processes beyond iron. Our framework is also well suited for studying explosive nucleosynthesis in supernovae. In this paper, we investigate 1D explosion models using the ‘PUSH method’ , focusing on progenitors with varying metallicities and initial masses around $M_{\rm ZAMS} = 16\, {\rm M}_{\odot }$. Detailed post-process nucleosynthesis calculations and Monte Carlo analyses are used to explore the effects of reaction rate uncertainties and to identify key reaction rates in explosive nucleosynthesis. We find that many reactions have little impact on the production of iron-group nuclei, as these elements are primarily synthesized in the nuclear statistical equilibrium. However, we identify a few ‘key reactions’ that significantly influence the production of radioactive nuclei, which may affect astrophysical observables. In particular, for the production of ${}^{44}{\rm Ti}$, we confirm that several traditionally studied nuclear reactions have a strong impact. However, determining a single reaction rate is insufficient to draw a definitive conclusion.

79 ASTRONOMY AND ASTROPHYSICS↗

Reevaluation of the O(+)(2P) reaction rate coefficients derived from Atmosphere Explorer C observations

O(+)(2P) is an important species for studies of the ionosphere and thermosphere: its emission at 7320 A can be used as a diagnostic of the thermospheric atomic oxygen density. Unfortunately, there are no laboratory measurements of the O and N2 reaction rates which are needed to determine the major sinks of (O+)(2p). We have recalculated the O and N2 reaction rates for O(+) (2P) using recent improvements in the solar EUV flux, cross sections, and photoelectron fluxes. For the standard solar EUV flux, the new N2 reaction rate of 3.4 +/- 1.5 x 10 exp -10 cu cm/s is close to the value obtained by Rusch et al. (1977), but the new O reaction rate of 4.0 +/- 1.9 x 10 exp -10 cu cm/sec is about 8 times larger. These new reaction rates are derived using neutral densities, electron density, and solar EUV fluxes measured by Atmosphere Explorer C in 1974 during solar minimum. The new theoretical emission rates are in good agreement with the data for the two orbits studied by Rusch et al.

Chang, T.↗

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↗

The Influence of Particle Charge on Heterogeneous Reaction Rate Coefficients

The effects of particle charge on heterogeneous reaction rates are presented. Many atmospheric particles, whether liquid or solid are charged. This surface charge causes a redistribution of charge within a liquid particle and as a consequence a perturbation in the gaseous uptake coefficient. The amount of perturbation is proportional to the external potential and the square of the ratio of debye length in the liquid to the particle radius. Previous modeling has shown how surface charge affects the uptake coefficient of charged aerosols. This effect is now included in the heterogeneous reaction rate of an aerosol ensemble. Extension of this analysis to ice particles will be discussed and examples presented.

Aikin, A. C.↗

Monte Carlo analysis of uncertainty propagation in a stratospheric model. 2: Uncertainties due to reaction rates

A concise stratospheric model was used in a Monte-Carlo analysis of the propagation of reaction rate uncertainties through the calculation of an ozone perturbation due to the addition of chlorine. Two thousand Monte-Carlo cases were run with 55 reaction rates being varied. Excellent convergence was obtained in the output distributions because the model is sensitive to the uncertainties in only about 10 reactions. For a 1 ppby chlorine perturbation added to a 1.5 ppby chlorine background, the resultant 1 sigma uncertainty on the ozone perturbation is a factor of 1.69 on the high side and 1.80 on the low side. The corresponding 2 sigma factors are 2.86 and 3.23. Results are also given for the uncertainties, due to reaction rates, in the ambient concentrations of stratospheric species.

Stolarski, R. S.↗

Experimental Determination of 𝛼 Widths of 21 Ne Levels in the Region of Astrophysical Interest: New 17 O +𝛼 Reaction Rates and Impact on the Weak 𝑠 Process

The efficiency of the weak 𝑠 process in low-metallicity rotating massive stars depends strongly on the rates of the competing 17 O ⁡(𝛼,𝑛)⁢ 20 Ne and 17 O ⁡(𝛼,𝛾)⁢ 21 Ne reactions that determine the potency of the 16 O neutron poison. Their reaction rates are poorly known in the astrophysical energy range of interest for core helium burning in massive stars because of the lack of spectroscopic information (partial widths, spin parities) for the relevant states in the compound nucleus 21 Ne . In this Letter, we report on the first experimental determination of the 𝛼-particle spectroscopic factors and partial widths of these states using the 17 O ⁡( 7 Li, 𝑡)⁢ 21 Ne 𝛼-transfer reaction. With these the 17 O ⁡(𝛼, 𝑛)⁢ 20 Ne and 17 O ⁡(𝛼, 𝛾)⁢ 21 Ne reaction rates were evaluated with uncertainties reduced by a factor more than 3 with respect to previous evaluations and the present 17 O ⁡(𝛼, 𝑛)⁢ 20 Ne reaction rate is more than 20 times larger. Finally, the present (𝛼, 𝑛)/(𝛼, 𝛾) rate ratio favors neutron recycling and suggests an enhancement of the weak 𝑠 process in the Zr-Nd region by more than 1.5 dex in metal-poor rotating massive stars.

direct reactions↗

Supervised learning and the finite-temperature string method for computing committor functions and reaction rates

A central object in the computational studies of rare events is the committor function. Though costly to compute, the committor function encodes complete mechanistic information of the processes involving rare events, including reaction rates and transition-state ensembles. Under the framework of transition path theory, Rotskoff et al. [Proceedings of the 2nd Mathematical and Scientific Machine Learning Conference, Proceedings of Machine Learning Research (PLMR, 2022), Vol. 145, pp. 757–780] proposes an algorithm where a feedback loop couples a neural network that models the committor function with importance sampling, mainly umbrella sampling, which collects data needed for adaptive training. Here, in this work, we show additional modifications are needed to improve the accuracy of the algorithm. The first modification adds elements of supervised learning, which allows the neural network to improve its prediction by fitting to sample-mean estimates of committor values obtained from short molecular dynamics trajectories. The second modification replaces the committor-based umbrella sampling with the finite-temperature string (FTS) method, which enables homogeneous sampling in regions where transition pathways are located. We test our modifications on low-dimensional systems with non-convex potential energy where reference solutions can be found via analytical or finite element methods, and show how combining supervised learning and the FTS method yields accurate computation of committor functions and reaction rates. We also provide an error analysis for algorithms that use the FTS method, using which reaction rates can be accurately estimated during training with a small number of samples. The methods are then applied to a molecular system in which no reference solution is known, where accurate computations of committor functions and reaction rates can still be obtained.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Advancing a Watershed Hydro-Biogeochemical Theory: Linking Water Travel Time and Reaction Rates Under Changing Climate

The major goal of the project is to answer the question how and to what degree does a warming climate alter hydrological functioning and (biogeochemical and chemical weathering) reaction rates at the watershed scale? In particular, how does hydrological functioning (water storage, flow paths, and travel time) change responding to a warming climate? What are the corresponding alterations in biogeochemical and chemical weathering rates at the watershed scale? What is the general watershed-scale reaction rate law that causally links reaction rates to metrics of hydrological functioning that are determined by external hydroclimatic conditions and internal watershed structure characteristics? To answer these questions, we have used field measurements and reactive transport modeling. We have found that: 1) stream chemistry has altered tremendously in warmer years; 2) In this water-limited environment, hydrology is the dominant driver for soil respiration and carbon lateral export at the watershed scale; 3) deeper carbon respiration plays an essential role in regulating the patterns of stream carbon in low and high flow conditions, and deeper groundwater will become increasingly important as warming intensifies. This indicates in mountain streams, a warming climate will modify not only the timing and magnitude of streamflow as expected, but also the chemistry and quality of stream water.

54 ENVIRONMENTAL SCIENCES↗

Fast estimation of reaction rates in spherical and non-spherical porous catalysts

Here we present a methodology for modeling multi-step reaction rates in porous catalyst particles for use in CFD-DEM and two fluid models. Single-step effectiveness factors based on a Thiele modulus, while useful, cannot accurately capture the cascading reaction systems common in high temperature vapor-phase chemical reactors like fluidized catalytic cracking units and catalytic biomass fast pyrolysis systems. Instead, multi-step effectiveness vectors derived from steady-state solutions to the governing reaction-diffusion equations are needed. Solutions for various catalyst shapes are presented, including spheres, cylinders, and prisms. Computational challenges inherent in repeated evaluation of reaction rates with diffusion limitations are discussed, and an efficient implementation based on pre-computed lookup tables is proposed and demonstrated on a simulation of a fluidized bed reactor. Open-source code is provided for the compilation of reaction rate tables for use in ODE, DEM, and two-fluid models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Thermonuclear reaction rate of 30 Si($p, γ$) 31 P

Silicon synthesis in high-temperature hydrogen burning environments presents one possible avenue for the study of abundance anomalies in globular clusters. This was suggested in a prior study, which found that the large uncertainties associated with the Si 30 ( p , γ ) P 31 reaction rate preclude a firm understanding of the stellar conditions that give rise to the Mg-K anticorrelation observed in the globular cluster NGC 2419. In an effort to improve the reaction rate, we present new strength measurements of the E r lab = 435 keV and E r lab = 501 keV resonances in Si 30 ( p , γ ) P 31 . For the former, which was previously unobserved, we obtain a resonance strength of ω γ = ( 1.14 ± 0.25 ) × 10 – 4 eV. For the latter, we obtain a value of ω γ = ( 1.88 ± 0.14 ) × 10 – 1 eV, which has a smaller uncertainty compared to previously measured strengths. Based on these results, the thermonuclear reaction rate has been re-evaluated. The impact of the new measurements is to lower the reaction rate by a factor of ≈ 10 at temperatures important to the study of NGC 2419. The rate uncertainty at these temperatures has also been reduced significantly.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Reaction rate constant of HO2+O3 measured by detecting HO2 from photofragment fluorescence

A room-temperature discharge-flow system investigation of the rate constant for the reaction 'HO2 + O3 yields OH + 2O2' has detected HO2 through the OH(A-X) fluorescence produced by photodissociative excitation of HO2 at 147 nm. A reaction rate constant of 1.9 + or - 0.3 x 10 to the -15th cu cm/molecule per sec is obtained from first-order decay of HO2 in excess O3; this agrees well with published data.

Manzanares, E. R.↗

Pycnonuclear reaction rates for binary ionic mixtures

Through a combination of compositional scaling arguments and examinations of Monte Carlo simulation results for the interparticle separations in binary-ionic mixture (BIM) solids, we have derived parameterized expressions for the BIM pycnonuclear rates as generalizations of those in one-component solids obtained previously by Salpeter and Van Horn and by Ogata et al. We have thereby discovered a catalyzing effect of the heavier elements, which enhances the rates of reactions among the lighter elements when the charge ratio exceeds a critical value of approximately 2.3.

Ichimaru, S.↗

Temperature-dependent reaction-rate expression for oxygen recombination at Shuttle entry conditions

A temperature-dependent oxygen surface reaction-rate coefficient has been determined from experimental STS-2 heating and wall temperature data at altitudes of 77.91 km, 74.98 km, and 71.29 km. The coefficient is presented in an Arrhenius form and is shown to be less temperature dependent than previous results. Finite-rate viscous-shock-layer heating rates based on this present expression have been compared with predicted heating rates using the previous rate coefficients and with experimental heating data obtained over an extensive range of STS-2 and STS-3 entry conditions. A substantial improvement is obtained in comparison of experimental data and predicted heating rates using the present oxygen reaction-rate expression.

Zoby, E. V.↗

Univariate Prediction of Hammett Parameters and Select Relative Reaction Rates Using Loewdin Atomic Charges

Loewdin charges from density functional theory calculations were used here to obtain general, univariate linear correlations for the prediction of experimental Hammett parameters and relative reaction rates. While previous studies have established that Hirshfeld and CM5 charges perform strongly as univariate predictors, the near-ubiquitous Loewdin charges have not yet been evaluated. To this end, we assess the predictive capability of Loewdin charges for three chemical systems. First, we show that Loewdin charges outperform Hirshfeld and CM5 charges for Hammett parameter prediction. Second, we see that Loewdin charges generally perform comparably to Hirshfeld charges for predicting the relative rates of olefin cleavage by photoexcited nitroarenes. The single case of poor correlation, between relative rates and the Loewdin charges on nitrogen sites, is ameliorated when considering the net charge on the NO 2 group. Third, we show that Loewdin, Hirshfeld, and CM5 charges all perform very well for generating correlations for relative reaction rates for C–H activation of 9-(4-X-phenyl)-9H-fluorene substrates by a transition metal catalyst. The equations generated throughout the study enable the prediction of Hammett parameters and relative reaction rates. Finally, these tools can accelerate synthetic and experimental studies by enabling the in silico prediction of uncharacterized chemical properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Study of the ${{^{25}}}$Mg(d,p)${{^{26}}}$Mg reaction to constrain the ${{^{25}}}$Al(p,${\gamma }$)${{^{26}}}$Si resonant reaction rates in nova burning conditions

The rate of the 25 Al( p,γ) 26 Si reaction is one of the few key remaining nuclear uncertainties required for predicting the production of the cosmic γ-ray emitter 26 Al in explosive burning in novae. This reaction rate is dominated by three key resonances (J π = 0 + , 1 + and 3 + ) in 26 Si. Only the 3 + resonance strength has been directly constrained by experiment. A high resolution measurement of the 25 Mg(d, p) reaction was used to determine spectroscopic factors for analog states in the mirror nucleus, 26 Mg. A first spectroscopic factor value is reported for the 0 + state at 6.256 MeV, and a strict upper limit is set on the value for the 1 + state at 5.691 MeV, that is incompatible with an earlier ( 4 He, 3 He) study. These results are used to estimate proton partial widths, and resonance strengths of analog states in 26 Si contributing to the 25 Al( p,γ) 26 Si reaction rate in nova burning conditions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗