Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “direct reaction theory”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 19 records

Neutron transfer reactions in halo effective field theory

Direct reaction experiments exhibit a powerful tool to probe the structure of neutron-rich nuclei like beryllium-11. We use halo effective field theory to calculate the cross section of the deuteron-induced neutron transfer reaction 10 Be(d, p) 11 Be. The effective theory contains dynamical fields for the beryllium-10 core, the neutron, and the proton. In contrast, the deuteron and the beryllium-11 halo nucleus are generated dynamically from contact interactions using experimental and ab initio input. Breakup contributions are then included by construction. The reaction amplitude is constructed up to next-to-leading order in an expansion in the ratio of the length scales characterizing the core and the halo. The Coulomb repulsion between core and proton is treated perturbatively. Lastly, we compare our results to cross-section data and other calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Statistical Hauser-Feshbach theory with width-fluctuation correction including direct reaction channels for neutron-induced reactions at low energies

Here, a model to calculate particle-induced reaction cross sections with statistical Hauser-Feshbach theory including direct reactions is given. The energy average of the scattering matrix from the coupled-channels optical model is diagonalized by the transformation proposed by Engelbrecht and Weidenmüller [C. A. Engelbrecht and H. A. Weidenmüller, Phys. Rev. C 8, 859 (1973)]. The ensemble average of S -matrix elements in the diagonalized channel space is approximated by a model of Moldauer [P. A. Moldauer, Phys. Rev. C 12, 744 (1975)] using the newly parametrized channel degree-of-freedom ν a to better describe the Gaussian orthogonal ensemble (GOE) reference calculations. The Moldauer approximation is confirmed by a Monte Carlo study using a randomly generated S matrix, as well as the GOE threefold integration formula. Lastly, the method proposed is applied to the 238U(n,n') cross-section calculation in the fast-energy range, showing an enhancement in the inelastic scattering cross sections.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Constraining transfer cross sections using Bayes' theorem

Background: Being able to rigorously quantify the uncertainties in reaction models is crucial to moving this field forward. Even though Bayesian methods are becoming increasingly popular in nuclear theory, they have yet to be implemented and applied in reaction theory problems. Purpose: The purpose of this work is to investigate, using Bayesian methods, the uncertainties in the optical potentials generated from fits to elastic-scattering data and the subsequent uncertainties in the transfer predictions. We also study the differences in two reaction models where the parameters are constrained in a similar manner, as well as the impact of reducing the experimental error bars on the data used to constrain the parameters. Method: We use Bayes' theorem combined with a Markov chain Monte Carlo to determine posterior distributions for the parameters of the optical model, constrained by neutron-, proton-, and/or deuteron-target elastic scattering. These potentials are then used to predict transfer cross sections within the adiabatic wave approximation or the distorted-wave Born approximation. Results: We study a number of reactions involving deuteron projectiles with energies in the range 10–25 MeV/nucleon on targets with mass A = 48–208. The case of 48 Ca(d,p) 49 Ca transfer to the ground state is described in detail. A comparative study of the effect of the size of experimental errors is also performed. Five transfer reactions are studied, and their results are then compiled in order to systematically identify trends.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Uncertainty quantification for optical model parameters

Although uncertainty quantification has been making its way into nuclear theory, these methods have yet to be explored in the context of reaction theory. For example, it is well known that different parameterizations of the optical potential can result in different cross sections, but these differences have not been systematically studied and quantified. The purpose of our work is to investigate the uncertainties in nuclear reactions that result from fitting a given model to elastic-scattering data, as well as to study how these uncertainties propagate to the inelastic and transfer channels. We use statistical methods to determine a best fit and create corresponding 95% confidence bands. A simple model of the process is fit to elastic-scattering data and used to predict either inelastic or transfer cross sections. In this initial work, we assume that our model is correct, and the only uncertainties come from the variation of the fit parameters. Here, we study a number of reactions involving neutron and deuteron projectiles with energies in the range of 5–25 MeV/u, on targets with mass A=12–208. We investigate the correlations between the parameters in the fit. The case of deuterons on 12 C is discussed in detail: the elastic-scattering fit and the prediction of 12 C(d,p) 13 C transfer angular distributions, using both uncorrelated and correlated χ 2 minimization functions. The general features for all cases are compiled in a systematic manner to identify trends. This work shows that, in many cases, the correlated χ 2 functions (in comparison to the uncorrelated χ 2 functions) provide a more natural parameterization of the process. These correlated functions do, however, produce broader confidence bands. Further optimization may require improvement in the models themselves and/or more information included in the fit.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Nearside-farside, local angular momentum and resummation theories: Useful tools for understanding the dynamics of complex-mode reactions

A valuable tool for understanding the dynamics of direct reactions is Nearside-Farside (NF) scattering theory. It makes a decomposition of the (resummed) partial wave series for the scattering amplitude, both for the differential cross section (DCS) and the Local Angular Momentum (LAM). This paper makes the first combined application of these techniques to complex-mode reactions. We ask if NF theory is a useful tool for their identification, in particular, can it distinguish complex-mode from direct-mode reactions? We also ask whether NF theory can identify NF interference oscillations in the full DCSs of complex-mode reactions. Our investigation exploits the fact that accurate quantum scattering matrix elements have recently become available for complex-mode reactions. We first apply NF theory to two simple models for the scattering amplitude of a complex-mode reaction: One involves a single Legendre polynomial; the other involves a single Legendre function of the first kind, whose form is suggested by complex angular momentum theory. We then study, at fixed translational energies, four state-to-state complex-mode reactions. They are: S({sup 1}D) + HD → SH + D, S({sup 1}D) + DH → SD + H, N({sup 2}D) +H{sub 2} → NH + H, and H{sup +} + D{sub 2} → HD + D{sup +}. We compare the NF results for the DCSs and LAMs with those for a state-to-state direct reaction, namely, F + H{sub 2} → FH + H. We demonstrate that NF theory is a valuable tool for identifying and analyzing the dynamics of complex-mode reactions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Variational Transition State Theory

This is the final report on a project involving the development and applications of variational transition state theory. This project involved the development of variational transition state theory for gas-phase reactions, including optimized multidimensional tunneling contributions and the application of this theory to gas-phase reactions with a special emphasis on developing reaction rate theory in directions that are important for applications to combustion. The development of variational transition state theory with optimized multidimensional tunneling as a useful computational tool for combustion kinetics involved eight objectives.

02 PETROLEUM↗

Advancing the Theory of Nuclear Reactions with Rare Isotopes: From the Laboratory to the Cosmos

The mission of the TORUS Topical Collaboration is to develop new methods that will advance nuclear reaction theory for unstable isotopes by using three-body techniques to improve direct-reaction calculations, and, by using a new partial-fusion theory, to integrate descriptions of direct and compound-nucleus reactions. Ohio University concentrates its efforts on the first part of the mission. Since direct measurements are often not feasible, indirect methods, e.g. (d,p) reactions, should be used. Those (d,p) reactions may be viewed as three-body reactions and described with Faddeev techniques. Faddeev equations in momentum space have a long tradition of utilizing separable interactions in order to arrive at sets of coupled integral equations in one variable. While there exist several separable representations for the nucleon-nucleon interaction, the optical potential between a neutron (proton) and a nucleus is not readily available in separable form. For this reason we first embarked in introducing a separable representation for complex phenomenological optical potentials of Woods-Saxon type.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

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

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

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ultradoping Boron on Si(100) via Solvothermal Chemistry

Ultradoping introduces unprecedented dopant levels into Si, which transforms its electronic behavior and enables its use as a next-generation electronic material. Commercialization of ultradoping is currently limited by gas-phase ultra-high vacuum requirements. Solvothermal chemistry is amenable to scale-up. However, an integral part of ultradoping is a direct chemical bond between dopants and Si, and solvothermal dopant-Si surface reactions are not well-developed. This work provides the first quantified demonstration of achieving ultradoping concentrations of boron (~1e14 cm 2 ) by using a solvothermal process. Surface characterizations indicate the catalyst cross-reacted, which led to multiple surface products and caused ambiguity in experimental confirmation of direct surface attachment. Density functional theory computations elucidate that the reaction results in direct B-Si surface bonds. Finally, this proof-of-principle work lays groundwork for emerging solvothermal ultradoping processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Implementing Chiral Three-Body Forces in Terms of Medium-Dependent Two-Body Forces

Three-nucleon (3N) forces are an indispensable ingredient for accurate few-body and many- body nuclear structure and reaction theory calculations. While the direct implementation of chiral 3N forces can be technically very challenging, a simpler approach is given by employing instead a medium-dependent NN interaction V med that reflects the physics of three-body forces at the two-body normal-ordered approximation. We review the derivation and construction of Vmed from the chiral 3N interaction at next-to-next-to-leading order (N2LO), consisting of a long-range 2π-exchange term, a mid-range 1π-exchange component and a short-range contact-term. Several applications of V med to the equation of state of cold nuclear and neutron matter, the nucleon single-particle potential in nuclear matter, and the nuclear quasiparticle interaction are discussed. We also explore differences in using local vs. nonlocal regulating functions on 3N forces and make direct comparisons to exact results at low order in perturbation theory expansions for the equation of state and single-particle potential. We end with a discussion and numerical calculation of the in-medium NN potential V med from the next-to-next-to-next-to-leading order (N3LO) chiral 3N force, which consists of a series of long-range and short-range terms.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The Trojan Horse method for nuclear astrophysics: Recent results for direct reactions

The Trojan Horse method is a powerful indirect technique to determine the astrophysical factor for binary rearrangement processes A+x→b+B at astrophysical energies by measuring the cross section for the Trojan Horse (TH) reaction A+a→B+b+s in quasi free kinematics. The Trojan Horse Method has been successfully applied to many reactions of astrophysical interest, both direct and resonant. In this paper, we will focus on direct sub-processes. The theory of the THM for direct binary reactions will be shortly presented based on a few-body approach that takes into account the off-energy-shell effects and initial and final state interactions. Examples of recent results will be presented to demonstrate how THM works experimentally.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Theoretical kinetics predictions for NH 2 + HO 2

Recent modeling studies of NH 3 oxidation, which are motivated by the prospective role of ammonia as a zero-carbon fuel, have indicated significant discrepancies between existing literature mechanisms. In this study high level theoretical kinetics predictions have been obtained for the reaction of NH 2 with HO 2 , which has previously been highlighted as an important reaction with high sensitivity and high uncertainty. The potential energy surface is explored with coupled cluster calculations including large basis sets and high-level corrections to yield high accuracy (~0.2 kcal/mol) estimates of the stationary point energies. Variational transition state theory is used to predict the microcanonical rate constants, which are then incorporated in master equation treatments of the temperature and pressure dependent kinetics. For the radical-radical channels, the microcanonical rates are obtained from variable reaction coordinate transition state theory implementing directly evaluated multireference electronic energies. The analysis yields predictions for the total rate constant as well as the branching to the NH 3 + O 2 , H 2 NO + OH, and HNO + H 2 O channels. Rate constants are also reported for the H 2 NO + OH reaction as they arise naturally from the analysis. The rate constant and branching fraction determined in this work for the NH 2 + HO 2 reaction deviate significantly from values used in most previous modeling studies. The fact that the main product channel is chain terminating, rather than propagating, has strong implications for modeling NH 3 ignition and oxidation, in particular at intermediate temperatures and elevated pressure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

TORUS: Theory of Reactions for Unstable iSotopes.Topical Collaboration for Nuclear Theory Project. Period: June 1, 2010 - May 31, 2015

The work of this collaboration during its existence is summarized. The mission of the TORUS Topical Collaboration was to develop new methods that advance nuclear reaction theory for unstable isotopes by using three-body techniques to improve direct reaction calculations. This multi-institution collaborative effort was and remains directly relevant to three areas of interest: the properties of nuclei far from stability, microscopic studies of nuclear input parameters for astrophysics, and microscopic nuclear reaction theory. The TORUS project focused on understanding the details of (d,p) reactions for neutron transfer to heavier nuclei. The bulk of the work fell into three areas: coupled channel theory, modeling (d,p) reactions with a Faddeev-AGS approach, and capture reactions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Theoretical Kinetics Predictions for Reactions on the NH 2 O Potential Energy Surface

Recent modeling studies of ammonia oxidation, which are motivated by the prospective role of ammonia as a zero-carbon fuel, have indicated significant discrepancies among the existing literature mechanisms. In this study, high-level theoretical kinetics predictions have been obtained for reactions on the NH 2 O potential energy surface, including the NH 2 + O, HNO + H, and NH + OH reactions. These reactions have previously been highlighted as important reactions in NH 3 oxidation with high sensitivity and high uncertainty. The potential energy surface is explored with coupled cluster calculations, including large basis sets and high-level corrections to yield high-accuracy (~0.2 kcal/mol 2σ uncertainty) estimates of the stationary point energies. Variational transition state theory is used to predict the microcanonical rate constants, which are then incorporated in master equation treatments of the temperature- and pressure-dependent kinetics. For radical–radical channels, the microcanonical rates are obtained from variable reaction coordinate transition state theory implementing directly evaluated multireference electronic energies. The analysis yields predictions for the total rate constants as well as the branching ratios. We find that the NO + H 2 channel contributes 10% of the total NH 2 + O flux at combustion temperatures, although this channel is not included in modern NH 3 oxidation mechanisms. Modeling is used to illustrate the ramifications of these rate predictions on the kinetics of NH 3 oxidation and NO x formation. Finally, the present results for NH 2 + O are important for predicting the chain branching and formation of NO in the oxidation of NH 3 and thermal DeNO x .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Direct Demonstration of Unified Brønsted-Evans-Polanyi Relationships for Proton-Coupled Electron Transfer Reactions on Transition Metal Surfaces

Brønsted-Evans-Polanyi (BEP) relationships, which relate elementary reaction barriers to reaction thermodynamics, have long been discussed in electrochemical science, but demonstration of their existence in heterogeneous electrocatalysis, across many different catalyst surfaces and voltages, is lacking. In this contribution, the BEP model is demonstrated to describe both multiple catalyst surfaces and variable voltages in electrochemical environments. For proton-coupled electron transfer (PCET) to surface nitrogen (N*) and nitric oxide (NO*), reaction energies and activation barriers are calculated using Density Functional Theory (DFT) on a parallel plate capacitor model at three different potentials and for nine different transition metal surfaces. Linear BEP relations that describe all potentials and catalyst surfaces are obtained for these elementary reactions, and the slopes of the correlations are shown to be directly related to the fractional coordinate of the transition states (FCTS) of the reactions. The results, which are explained in terms of Marcus Theory, prove a direct equivalence between unified BEP coefficients, describing both variable catalyst surfaces and voltages, and electrochemical symmetry factors and provide a straightforward means of estimating this quantify for PCET reactions on transition metal surfaces. The resulting relationships, in turn, could lead to predictions of electrocatalytic reactivity trends of enhanced accuracy and efficiency.

Electrochemistry↗

Communication: Direct comparison between theory and experiment for correlated angular and product-state distributions of the ground-state and stretching-excited O({sup 3}P) + CH{sub 4} reactions

Motivated by a recent experiment [H. Pan and K. Liu, J. Chem. Phys. 140, 191101 (2014)], we report a quasiclassical trajectory study of the O({sup 3}P) + CH{sub 4}(v{sub k} = 0, 1) → OH + CH{sub 3} [k = 1 and 3] reactions on an ab initio potential energy surface. The computed angular distributions and cross sections correlated to the OH(v = 0, 1) + CH{sub 3}(v = 0) coincident product states can be directly compared to experiment for O + CH{sub 4}(v{sub 3} = 0, 1). Both theory and experiment show that the ground-state reaction is backward scattered, whereas the angular distributions shift toward sideways and forward directions upon antisymmetric stretching (v{sub 3}) excitation of the reactant. Theory predicts similar behavior for the O + CH{sub 4}(v{sub 1} = 1) reaction. The simulations show that stretching excitation enhances the reaction up to about 15 kcal/mol collision energy, whereas the O + CH{sub 4}(v{sub k} = 1) reactions produce smaller cross sections for OH(v = 1) + CH{sub 3}(v = 0) than those of O + CH{sub 4}(v = 0) → OH(v = 0) + CH{sub 3}(v = 0). The former finding agrees with experiment and the latter awaits for confirmation. The computed cold OH rotational distributions of O + CH{sub 4}(v = 0) are in good agreement with experiment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A “Master Curve” Describing Reaction Inhomogeneity and Plating Onset during Fast-Charging of Graphite Electrodes

Lithium plating in porous graphite electrodes is a major limitation for fast charging. Theoretical evidence suggests plating during fast charging is largely due to inhomogeneous intercalation through the electrode thickness (caused by ionic/mass transfer limitations and inherent thermodynamic properties of the graphite material). Numerical and analytical predictions of plating onset during fast charging have been proposed but not compared directly with experiments. This work validates these model predictions against plating onset measurements via the “dOCV” method in graphite half-cells, for various electrode thicknesses and C-rates. Remarkably, it is shown that experimental and theoretical trends in plating onset (vs C-rate and electrode thickness) collapse to a single “master curve” using a nondimensional “reaction inhomogeneity” parameter. This observation supports the hypothesized reaction inhomogeneity mechanism for accelerating plating onset and provides practical guidance for electrode design. Furthermore, this work develops theory to quantify reaction inhomogeneity in situ directly from voltage V vs capacity Q data, using dQ / dV analysis. Here again, experiments and numerical predictions show good agreement, where peaks that correspond to LiC 6 correlate with plating onset. This work provides experimental validation of theoretical tools that can predict plating onset, aid electrode design, and give insight for plating onset mechanisms during fast charging.

25 ENERGY STORAGE↗

Indirect reactions and connection with R-matrix theory

It is often the case that nuclear reactions that are important for societal applications or basic science are difficult to measure directly in existing facilities for accelerated beams. The difficulty might be associated with an exceedingly small cross section, as the ones obtained at beam energies well below the Coulomb barrier, or with the impossibility of devising a short-lived target, as is the case for neutron-induced reactions on unstable isotopes. The fruitful line of experimental research addressing this issue with alternative (indirect) reactions has been developed in parallel to the theory needed to make the connection between the observed data and the desired cross sections of the (direct) reaction under study. Within this context, we present here an explicit connection between the indirect cross section and the R-matrix parameters which describe the energy-dependent cross section of the desired direct reaction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗