Search NASA⌕ Search

SEARCH · Search NASA

Results for “truncation uncertainty”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Factorial growth at low orders in perturbative QCD: control over truncation uncertainties

A method, known as “minimal renormalon subtraction”, relates the factorial growth of a perturbative series (in QCD) to the power p of a power correction Λ p /Q p . (Λ is the QCD scale, Q some hard scale.) Here, the derivation is simplified and generalized to any p , more than one such correction, and cases with anomalous dimensions. Strikingly, the well-known factorial growth is seen to emerge already at low or medium orders, as a consequence of constraints on the Q dependence from the renormalization group. The effectiveness of the method is studied with the gluonic energy between a static quark and static antiquark (the “static energy”). Truncation uncertainties are found to be under control after next-to-leading order, despite the small exponent of the power correction ( p = 1) and associated rapid growth seen in the first four coefficients of the perturbative series.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Uncertainties in ab initio nuclear structure calculations with chiral interactions

We present theoretical ground state energies and their uncertainties for p -shell nuclei obtained from chiral effective field theory internucleon interactions as a function of chiral order, fitted to two- and three-body data only. We apply a Similary Renormalization Group transformation to improve the numerical convergence of the many-body calculations, and discuss both the numerical uncertainties arising from basis truncations and those from omitted induced many-body forces, as well as chiral truncation uncertainties. With complete Next-to-Next-to-Leading (N 2 LO) order two- and three-body interactions, we find significant overbinding for the ground states in the upper p -shell, but using higher-order two-body potentials, in combination with N 2 LO three-body forces, our predictions agree with experiment throughout the p -shell to within our combined estimated uncertainties. The uncertainties due to chiral order truncation are noticeably larger than the numerical uncertainties, but they are expected to become comparable to the numerical uncertainties at complete N 3 LO.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Ab initio nucleon-nucleus elastic scattering with chiral effective field theory uncertainties

Effective interactions for nucleon-nucleus elastic scattering from first principles require the use of the same nucleon-nucleon interaction in the structure and reaction calculations and a consistent treatment of the relevant operators at each order. Systematic investigations of the effect of truncation uncertainties of chiral nucleon-nucleon ( N N ) forces have been carried out for scattering observables in the two- and three-nucleons system as well as for bound-state properties of light nuclei. Here we extend this type of study to proton and neutron elastic scattering for 16 O and 12 C. Using the frameworks of the spectator expansion of multiple scattering theory as well as the no-core shell model, we employ one specific chiral interaction from the LENPIC collaboration and consistently calculate the leading-order effective nucleon-nucleus interaction up to the third chiral order, from which we extract elastic scattering observables. Then we apply pointwise as well as correlated uncertainty quantification for the estimation of the chiral truncation error. We calculate and analyze proton elastic scattering observables for 16 O and 12 C as well as neutron elastic scattering observables for 12 C between 65 and 185 MeV projectile kinetic energy. We find in all cases qualitatively similar results for the chiral truncation uncertainties as in few-body systems and assess them using similar diagnostic tools. The order-by-order convergence of the scattering observables for 16 O and 12 C is very reasonable around 100 MeV, while for higher energies the chiral expansion parameter becomes too large for convergence. Comparing proton and neutron differential cross sections reveals that their description is comparable up to around 100 MeV. Here, we also find a nearly perfect correlation between the differential cross section for neutron scattering and the N N Wolfenstein amplitudes for small momentum transfers. The diagnostic tools for studying order-by-order convergence in observables in few-body systems can be employed for observables in nucleon-nucleus scattering with only minor modifications provided the momentum scale in the problems is not too large. We also find that the chiral N N interaction on which our study is based gives a very good description of differential cross sections for 16 O and 12 C as low as 65-MeV projectile energy, particularly in the forward direction. In addition, the very forward direction of the neutron differential cross section mirrors the behavior of the N N interaction amazingly well.

6 ≤ A ≤ 19↗
DOE OSTI2022

Bayesian analysis of nuclear polarizability corrections to the Lamb shift of muonic H-atoms and He-ions

The extraction of nuclear charge radii from spectroscopy experiments in muonic atoms is currently limited by the large uncertainties associated with the theoretical evaluation of the nuclear polarizability effects. To facilitate calculations, these polarizability corrections are conventionally expressed as an expansion in a dimensionless parameter η, which has been argued in previous literature to hold an approximate value of 0.33 in light-nuclear systems. In this work, we check this claim by doing a Bayesian analysis of the nuclear-polarizability corrections to the Lamb shift in µ 2 H and µ 3 H atoms and in µ 3 He + and µ 4 He + ions at various orders in the η-expansion. As presented here, our analysis supports the claim that η $\ll$ 1 in these systems and finds truncation uncertainties that are similar to the previous estimate, the only exception being the truncation uncertainties in the µ 3 He + ion, which are found to be larger.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Order-by-order uncertainties of nucleon-nucleon Wolfenstein amplitudes in chiral effective field theory

Quantum mechanical invariance principles dictate the most general operator structure that can be present in the nucleon-nucleon (NN) interaction. Five independent operators appear in the on-shell NN amplitude together with five corresponding coefficient functions. The usual choice for these coefficient functions is known as the NN Wolfenstein amplitudes. We analyze the order-by-order convergence of each of the five NN Wolfenstein amplitudes predicted by a semilocal coordinate space potential implementation of chiral effective field theory (𝜒⁢EFT). We do this at laboratory kinetic energies between 25 and 200 MeV for both neutron-proton and proton-proton scattering. Our analysis uses the Gaussian-process methods developed by the BUQEYE Collaboration to describe the contributions of each 𝜒⁢EFT order, and so yields truncation uncertainties for each Wolfenstein amplitude that are correlated across scattering angles. We combine information on the size of different orders in the EFT to infer the 𝜒⁢EFT breakdown scale for each amplitude, finding, on average, Λ 𝑏 between 750 and 800 MeV. Furthermore, with this choice of Λ 𝑏 , the EFT truncation uncertainties cover both higher-order results and empirical Wolfenstein amplitudes well for all orders other than the leading order.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Emulators for Scarce and Noisy Data: Application to Auxiliary-Field Diffusion Monte Carlo for Neutron Matter

Understanding the equation of state (EOS) of pure neutron matter is necessary for interpreting multimessenger observations of neutron stars. Reliable data analyses of these observations require well-quantified uncertainties for the EOS input, ideally propagating uncertainties from nuclear interactions directly to the EOS. This, however, requires calculations of the EOS for a prohibitively larger number of nuclear Hamiltonians, solving the nuclear many-body problem for each one. Quantum Monte Carlo methods, such as auxiliary-field diffusion Monte Carlo (AFDMC), provide precise and accurate results for the neutron matter EOS, but they are very computationally expensive, making them unsuitable for the fast evaluations necessary for uncertainty propagation. Here, we employ parametric matrix models to develop fast emulators for AFDMC calculations of neutron matter and use them to directly propagate uncertainties of coupling constants in the Hamiltonian to the EOS. As these uncertainties include estimates of the effective field theory truncation uncertainty, this approach provides robust uncertainty estimates for use in astrophysical data analyses. In conclusion, this Letter will enable novel applications such as using astrophysical observations to put constraints on coupling constants for nuclear interactions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Analyzing rotational bands in odd-mass nuclei using effective field theory and Bayesian methods

We recently developed an Effective Field Theory (EFT) for rotational bands in odd-mass nuclei. Here we use EFT expressions to perform a Bayesian analysis of data on the rotational energy levels of 99 Tc, 155,157 Gd, 159 Dy, 167,169 Er, 167,169 Tm, 183 W, 235 U and 239 Pu. The error model in our Bayesian analysis includes both experimental and EFT truncation uncertainties. It also accounts for the fact that low-energy constants (LECs) at even and odd orders are expected to have different sizes. We use Markov Chain Monte Carlo (MCMC) sampling to explore the joint posterior of the EFT and error-model parameters and show both the LECs and the breakdown scale can be reliably determined. We extract the LECs up to fourth order in the EFT and find that, provided we correctly account for EFT truncation errors in our likelihood, results for lower-order LECs are stable as we go to higher orders. LEC results are also stable with respect to the addition of higher-energy data. We extract the expansion parameter for all the nuclei listed above and find a clear correlation between the extracted and the expected value of the inverse breakdown scale, W, based on the single-particle and vibrational energy scales. However, the W that actually determines the convergence of the EFT expansion is markedly smaller than would be naively expected based on those scales.

EFT↗
DOE OSTI2021

Light nuclei with semilocal momentum-space regularized chiral interactions up to third order

In this work, we present a systematic investigation of few-nucleon systems and light nuclei using the current Low Energy Nuclear Physics International Collaboration interactions comprising semilocal momentum-space regularized two- and three-nucleon forces up to third chiral order (N 2 LO). Following our earlier study utilizing the coordinate-space regularized interactions, the two low-energy constants entering the three-body force are determined from the triton binding energy and the differential cross-section minimum in elastic nucleon-deuteron scattering. Predictions are made for selected observables in elastic nucleon-deuteron scattering and in the deuteron breakup reactions, for properties of the A = 3 and A = 4 nuclei, and for spectra of p-shell nuclei up to A = 16. A comprehensive error analysis is performed including an estimation of correlated truncation uncertainties for nuclear spectra. The obtained predictions are generally found to agree with experimental data within errors. Similarly to the coordinate-space regularized chiral interactions at the same order, a systematic overbinding of heavier nuclei is observed, which sets in for A~10 and increases with A.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Factorial growth in perturbation theory, power corrections: precise extraction of quark masses and $\alpha_\text{s}$

These proceedings summarize a newly found connection between the factorial growth of coefficients in perturbative QCD and power corrections to the perturbation series, discussed in refs. [1-4]. The improved convergence is shown for three quantities four which four terms in the series are available: the static energy, the quark pole mass, and the polarized Bjorken sum rule. Prospects for determinations of $\alpha_\text{s}$ with controlled truncation uncertainties are discussed, as was found earlier in quark-mass determinations [3,5].

Kronfeld, Andreas S. [Fermilab; TUM-IAS, Munich] (↗
DOE OSTI2023

Chiral uncertainties in ab initio elastic nucleon-nucleus scattering

The effective interaction between a nucleon and a nucleus is one of the most important ingredients for reaction theories. Theoretical formulations were introduced early by Feshbach and Watson, and efforts of deriving and computing those ‘optical potentials’ in a microscopic fashion have a long tradition. However, only recently the leading order term in the Watson multiple scattering approach could be calculated fully ab initio, meaning that the same nucleon-nucleon (NN) interaction enters both the structure as well as the reaction pieces on equal footing. This allows the uncertainties from the underlying chiral effective NN interaction to be systematically explored in nucleon-nucleus elastic scattering observables. In this contribution the main ingredients for arriving at the ab initio leading order of the effective nucleon-nucleus interaction in the Watson approach will be reviewed. Concentrating on one specific chiral NN interaction from the LENPIC collaboration and light nuclei with a 0 + ground state, the leading order nucleon-nucleus interaction is calculated using up to the third chiral order (N2LO) in the nucleon-nucleon potential, and elastic scattering observables are extracted. Then pointwise as well as correlated uncertainty quantification is used for the estimation of the chiral truncation error. Elastic scattering observables for 4 He, 12 C, and 16 O for between 65 and 200 MeV projectile energy will be analyzed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Reaction mechanisms of a cyclic ether intermediate: Ethyloxirane

Oxiranes are a class of cyclic ethers formed in abundance during low-temperature combustion of hydrocarbons and biofuels, either via chain-propagating steps that occur from unimolecular decomposition of β -hydroperoxyalkyl radicals ( β -˙QOOH) or from reactions of HOÓ with alkenes. Ethyloxirane is one of four alkyl-substituted cyclic ether isomers produced as an intermediate from n -butane oxidation. While rate coefficients for β -˙QOOH → ethyloxirane + ÓH are reported extensively, subsequent reaction mechanisms of the cyclic ether are not. As a result, chemical kinetics mechanisms commonly adopt simplified chemistry to describe ethyloxirane consumption by convoluting several elementary reactions into a single step, which may introduce mechanism truncation error—uncertainty derived from missing or incomplete chemistry. Here we provide fundamental insight on reaction mechanisms of ethyloxirane in support of ongoing efforts to minimize mechanism truncation error. Reaction mechanisms are inferred from the detection of products during chlorine atom-initiated oxidation experiments using multiplexed photoionization mass spectrometry conducted at 10 Torr and temperatures of 650 K and 800 K. To complement the experiments, calculations of stationary point energies were conducted using the ccCA-PS3 composite method on ˙R + O 2 potential energy surfaces for the four ethyloxiranyl radical isomers, which produced barrier heights for 24 reaction pathways. In addition to products from ˙QOOH → cyclic ether + ÓH and ˙R + O 2 → conjugate alkene + HOÓ, both of which were significant pathways and are prototypical to alkane oxidation, other species were identified from ring-opening of both ethyloxiranyl and ˙QOOH radicals. The latter occurs when the unpaired electron is localized on the ether group, causing the initial ˙QOOH structure to ring-open and form a resonance-stabilized ketohydroperoxide-type radical. The present work provides the first analysis of ethyloxirane oxidation chemistry, which reveals that consumption pathways are complex and may require an expansion of submechanisms to increase the fidelity of chemical kinetics mechanisms.

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

Isomer-Dependent Reaction Mechanisms of Cyclic Ether Intermediates: cis -23-Dimethyloxirane and trans -23-Dimethyloxirane.

Oxiranes are a class of cyclic ethers formed in abundance during low-temperature combustion of hydrocarbons and biofuels, either via chain-propagating steps that occur from unimolecular decomposition of β-hydroperoxyalkyl radicals (β-˙QOOH) or from reactions of HÒO with alkenes. The cis- and trans-isomers of 2,3-dimethyloxirane are intermediates of n-butane oxidation, and while rate coefficients for β-˙QOOH → 2,3-dimethyloxirane + OH are reported extensively, subsequent reaction mechanisms of the cyclic ethers are not. As a result, chemical kinetics mechanisms commonly adopt simplified chemistry to describe the consumption of 2,3-dimethyloxirane by convoluting several elementary reactions into a single step, which may introduce mechanism truncation error—uncertainty derived from missing or incomplete chemistry. The present research examines the isomerdependence of 2,3-dimethyloxirane reaction mechanisms in support of ongoing efforts to minimize mechanism truncation error. Reaction mechanisms are inferred via the detection of products from Cl-initiated oxidation of both cis-2,3-dimethyloxirane and trans-2,3-dimethyloxirane using multiplexed photoionization mass spectrometry (MPIMS). The experiments were conducted at 10 Torr and temperatures of 650 K and 800 K. To complement the experiments, the enthalpies of stationary points on the ˙R + O 2 surfaces were computed at the ccCA-PS3 level of theory. In total, 28 barrier heights were computed on the 2,3-dimethyloxiranylperoxy surfaces. Two notable aspects are low-lying pathways that form resonance-stabilized ketohydroperoxide-type radicals caused by ˙QOOH ring-opening when the unpaired electron is localized adjacent to the ether group, and cis-trans isomerization of ˙R and ˙QOOH radicals, via inversion, which enable reaction pathways otherwise restricted by stereochemistry. Several species were identified in the MPIMS experiments from ring opening of 2,3-dimethyloxiranyl radicals. Neither of the two conjugate alkene isomers prototypical of ˙R + O 2 reactions were detected. Products were also identified from decomposition of ketohydroperoxide-type radicals. The present work provides the first analysis of 2,3-dimethyloxirane oxidation chemistry and reveals that consumption pathways are complex and require the expansion of submechanisms in chemical kinetics mechanisms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2012

System Identification and Uncertainty Quantification Using Orthogonal Excitations and the Semi-span Super Sonic Transport (S4T) Model

Orthogonal harmonic multisine excitations were utilized in a wind tunnel test and in simulation of the SemiSpan Supersonic Transport model to assess aeroservoelastic characteristics. Fundamental issues associated with analyzing sinusoidal signals were examined, including spectral leakage, excitation truncation, and uncertainties on frequency response functions and mean-square coherence. Simulation allowed for evaluation of these issues relative to a truth model, while wind tunnel data introduced real-world implementation issues.

Heeg, Jennifer↗
DOE OSTI2024

Probing O 2 -dependence of cyclopentyl reactions via isomer-resolved speciation

Modeling chemical kinetics relevant to low-temperature combustion requires complete description of reactions involving critical species such as hydroperoxyalkyl radicals, Q̇OOH, which undergo competing unimolecular reactions and bimolecular reactions with O 2 . The balance of flux across the two pathways affects rates of chain-branching and depends on temperature, pressure, and oxygen concentration. Accordingly, the influence of [O 2 ] on product formation from alkyl + O 2 reactions and the subsequent fate of Q̇OOH and related products is central to the development of an accurate chemical kinetics mechanism. Furthermore, chemical reactions consuming Q̇OOH-mediated species are often simplified to such a degree that mechanism truncation error (uncertainty derived from incomplete reaction networks) becomes significant.

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

Effective field theory analysis of 3 He–α scattering data

We treat low-energy 3 He–α elastic scattering in an effective field theory (EFT) that exploits the separation of scales in this reaction. We compute the amplitude up to next-to-next-to-leading order, developing a hierarchy of the effective-range parameters (ERPs) that contribute at various orders. We use the resulting formalism to analyse data for recent measurements at center-of-mass energies of 0.38–3.12 MeV using the scattering of nuclei in inverse kinematics (SONIK) gas target at TRIUMF as well as older data in this energy regime. We employ a likelihood function that incorporates the theoretical uncertainty due to truncation of the EFT and use Markov chain Monte Carlo sampling to obtain the resulting posterior probability distribution. We find that the inclusion of a small amount of data on the analysing power A y is crucial to determine the sign of the p-wave splitting in such an analysis. Here, the combination of A y and SONIK data constrains all ERPs up to O(p 4 ) in both s- and p-waves quite well. The asymptotic normalisation coefficients and s-wave scattering length are consistent with a recent EFT analysis of the capture reaction 3 He(α, γ) 7 Be.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Rigorous constraints on three-nucleon forces in chiral effective field theory from fast and accurate calculations of few-body observables

We explore the constraints on the three-nucleon force (3NF) of chiral effective field theory (χ EFT) that are provided by bound-state observables in the A = 3 and A = 4 sectors. Our statistically rigorous analysis incorporates experimental error, computational method uncertainty, and the uncertainty due to truncation of the χ EFT expansion at next-to-next-to-leading order. A consistent solution for the 3 H binding energy, the 4 He binding energy and radius, and the 3 H β-decay rate can only be obtained if χ EFT truncation errors are included in the analysis. Here, the β-decay rate is the only one of these that yields a nondegenerate constraint on the 3NF low-energy constants, which makes it crucial for the parameter estimation. We use eigenvector continuation for fast and accurate emulation of no-core shell model calculations of the few-nucleon observables. This facilitates sampling of the posterior probability distribution, allowing us to also determine the distributions of the parameters that quantify the truncation error. We find a χ EFT expansion parameter of Q = 0.33 ± 0.06 for these observables.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1982

Synthesis of a modern/classical controller for large space structures

The paper is concerned with wide-sense stationary random vibrations of a representative large space structure subject to nonwhite excitations, and with synthesis of its controller amid modeling uncertainties due to truncated modes. Dynamics is analyzed by residue calculus in frequency domain. Steady state variance of a modal coordinate is shown to depend on, among other things, power spectral density of excitation at the modal frequency. Modal cost analysis is performed for an order reduction with a performance objective of attitude, attitude rate and energy control. Optimal and suboptimal output feedback controllers are compared. By treating deleted modes as an additive and a multiplicative perturbation in an idealized transfer function, the low frequency modes omitted according to modal cost analysis are found to cause instability. Finally, the size of the idealized model is shown to be governed by light damping and stringency in performance and stability specifications.

Hablani, H. B.↗
NASA NTRS1981

Model error estimation for large flexible spacecraft

The basic model obtained for spacecraft control system design applications is used to develop the control and estimation algorithms which constitute integral elements of the control system. The performance of this system is limited on account of modeling errors. The presence of such errors is inevitable in connection with truncated dynamics, parameter uncertainties, neglected nonlinearities, and external disturbances. In many cases, an approximate knowledge of the model errors can only be established by means of an estimation procedure. The present investigation has the objective to outline an approach to estimate errors inherent in a reduced-order model and to illustrate its application to modeling problems arising in spacecraft control system design. The considered procedure makes use of the principle of least-squares. In the discussed application this principle is employed to minimize a quadratic functional of the model errors.

Rodriguez, G.↗