Search NASA⌕ Search

SEARCH · Search NASA

Results for “CONTINUUM RADIATION”

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

Three-dimensional Skyrme Hartree-Fock-Bogoliubov solver in coordinate-space representation

The coordinate-space representation of the Hartree-Fock-Bogoliubov theory is the method of choice to study weakly bound nuclei whose properties are affected by the quasiparticle continuum space. To describe such systems, we developed a three-dimensional Skyrme-Hartree-Fock-Bogoliubov solver HFBFFT based on the existing, highly optimized and parallelized Skyrme-Hartree-Fock code Sky3D. The code does not impose any self-consistent spatial symmetries such as mirror inversions or parity. The underlying equations are solved in HFBFFT directly in the canonical basis using the fast Fourier transform. To remedy the problems with pairing collapse, we implemented the soft energy cutoff and pairing annealing. The convergence of HFB solutions was improved by a sub-iteration method. The Hermiticity violation of differential operators brought by Fourier-transform-based differentiation has also been solved. Furthermore, the accuracy and performance of HFBFFT were tested by benchmarking it against other HFB codes, both spherical and deformed, for a set of nuclei, both well-bound and weakly-bound.

3D coordinate-space representation↗
DOE OSTI2024

Further investigation on the fusion of 6 Li with 209 Bi target at near-barrier energies

The complete and incomplete fusion cross sections for 6 Li+ 209 Bi were measured using the in-beam γ-ray method around the Coulomb barrier. The cross sections of (deuteron captured) incomplete fusion (ICF) products were re-quantified experimentally for this reaction system. Here, the results reveal that the ICF cross section is equivalent to that of complete fusion (CF) above the Coulomb barrier and dominant near or below the Coulomb barrier. A theoretical calculation based on the continuum discretized coupled channel (CDCC) method was performed for the aforementioned CF and ICF cross sections; the result is consistent with the experiments. The universal fusion function (UFF) was also compared with the measured CF cross section for different barrier parameters, demonstrating that the CF suppression factor is significantly influenced by the choice of potential, which can reflect both dynamic and static effects of breakup on the fusion process.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Neutron-unbound states in 31 Ne

Background: The Island of Inversion near the N = 20 shell gap is home to nuclei with reordered single-particle energy levels compared to the spherical shell model. Studies of 31 Ne have revealed that its ground state has a halo component, characterized by a valence neutron orbiting a deformed 30 Ne core. This lightly-bound nucleus with a separation energy of only S n = 170 keV is expected to have excited states that are neutron unbound. Purpose: The purpose of this experiment was to investigate the low-lying excited states in 31 Ne that decay by the emission of a single neutron. Methods: An 89 MeV/u 33 Mg beam impinged on a segmented Be reaction target. Neutron unbound states in 31 Ne were populated via a two-proton knockout reaction. The 30 Ne fragment and associated neutron from the decay of 31 Ne were detected by the MoNA-LISA-Sweeper experimental setup at the National Superconducting Cyclotron Laboratory. Invariant mass spectroscopy was used to reconstruct the two-body decay energy ( 30 Ne + n). Results: The two-body decay energy spectrum exhibits two features: a low-lying peak at 0.30 ± 0.17 MeV and a broad enhancement at 1.50 ± 0.33 MeV, each fit with an energy-dependent asymmetric Breit-Wigner line shape representing a resonance in the continuum. Accompanying shell model calculations using the FSU interaction within NuShellX, combined with cross-section cal culations using the eikonal reaction theory, indicate that these peaks in the decay energy spectrum are caused by multiple resonant states in 31 Ne. Conclusions: In this study, excited states in 31 Ne were observed for the first time. Transitions from calculated shell model final states in 31 Ne to bound states in 30 Ne are in good agreement with the measured decay energy spectrum. Cross-section calculations for the two-proton knockout populating 31 Ne states as well as spectroscopic factors pertaining to the decay of 31 Ne into 30 Ne are used to examine the results within the context of the shell model expectations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Superfluid condensate fraction and pairing wave function of the unitary Fermi gas

The unitary Fermi gas is a many-body system of two-component fermions with zero-range interactions tuned to infinite scattering length. Despite much activity and interest in unitary Fermi gas and its universal properties, there have been great difficulties in performing accurate calculations of the superfluid condensate fraction and pairing wave function. In this paper, we present auxiliary-field lattice Monte Carlo simulations using a lattice interaction which accelerates the approach to the continuum limit, thereby allowing for robust calculations of these difficult observables. As a benchmark test, we compute the ground-state energy of 33 spin-up and 33 spin-down particles. As a fraction of the free Fermi gas energy $E_{\text{FG}}$, we find $E_0/E_{\text{FG}}$ = 0.369(2), 0.372(2), using two different definitions of the finite-system energy ratio, in agreement with the latest theoretical and experimental results. We then determine the condensate fraction by measuring off-diagonal long-range order in the two-body density matrix. We find that the fraction of condensed pairs is α = 0.43(2). Further, we also extract the pairing wave function and find the pair correlation length to be $ζ_pk_F$ = 1.8(3)ℏ, where $k_F$ is the Fermi momentum. Provided that the simulations can be performed without severe sign oscillations, the methods we present here can be applied to superfluid neutron matter as well as more exotic $\textit{P}$-wave and $\textit{D}$-wave superfluids.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

A Coupled Analysis of Atmospheric Mass Loss and Tidal Evolution in XUV Irradiated Exoplanets: The TRAPPIST-1 Case Study

Exoplanets residing close to their stars can experience evolution of both their physical structures and their orbits due to the influence of their host stars. In this work, we present a coupled analysis of dynamical tidal dissipation and atmospheric mass loss for exoplanets in X-ray and ultraviolet (XUV) irradiated environments. As our primary application, we use this model to study the TRAPPIST-1 system and place constraints on the interior structure and orbital evolution of the planets. We start by reporting on an ultraviolet continuum flux measurement (centered around ∼1900 Å) for the star TRAPPIST-1, based on 300 ks of Neil Gehrels Swift Observatory data, and which enables an estimate of the XUV-driven thermal escape arising from XUV photodissociation for each planet. We find that the X-ray flaring luminosity, measured from our X-ray detections, of TRAPPIST-1 is 5.6 × 10{sup −4} L {sub *}, while the full flux including non-flaring periods is 6.1 × 10{sup −5} L {sub *}, when L {sub *} is TRAPPIST-1's bolometric luminosity. We then construct a model that includes both atmospheric mass loss and tidal evolution and requires the planets to attain their present-day orbital elements during this coupled evolution. We use this model to constrain the ratio Q{sup ′}=3Q/2k{sub 2} for each planet. Finally, we use additional numerical models implemented with the Virtual Planet Simulator VPLanet to study ocean retention for these planets using our derived system parameters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Which nuclear data can be validated with LLNL pulsed-sphere experiments?

It is shown here that 14-MeV D+T LLNL pulsed-sphere experiments bring complementary information into the process of validating nuclear data compared to experiments that are traditionally used for this purpose—such as critical assemblies. To be more specific, the 14-MeV D+T LLNL pulsed-sphere neutron-leakage spectra enable to validate scattering and fission nuclear data up to 15 MeV (compared to approximately up to 5 MeV when using criticality experiments) and employ to this end simple compound targets containing only few isotopes. In this work, sensitivity profiles of the spectra to nuclear data are calculated in order to understand in detail which isotopes, observables, and energy ranges of nuclear data contribute significantly to their simulation. These profiles are presented for a few selected spheres containing 16 O, 12 C, 56 Fe, and 239 Pu. It is shown that the neutron-leakage spectra of spheres containing light elements are mostly sensitive to elastic- and inelastic-scattering cross sections on discrete levels and corresponding angular distributions. Spheres of structural materials are sensitive to elastic- and inelastic-scattering cross sections, including scattering on discrete levels and the continuum, and double-differential cross sections. Actinide spheres are also strongly sensitive to the fission observables, in particular to the total-fission neutron spectrum. Thin spheres (in which neutrons experience on average less than one scatter) are mostly sensitive to data near the elastic peak, in the energy range from 12–15 MeV, while thicker ones can be sensitive to data at lower incident-neutron energies due to multiple-scattering effects. This information is brought together with simulations of 71 pulsed-sphere neutron-leakage spectra using the ENDF/B-VII.1 and ENDF/B-VIII.0 nuclear-data libraries. This analysis highlights ENDF/B-VIII.0 data that could be further investigated for potential shortcomings ( 6 Li, 12 C, 16 O, 24-26 Mg, 27 Al, 48 Ti, 56 Fe, and 208 Pb) or are likely reliable ( 1,2 H, 7 Li, 9 Be, 14 N, 235,238 U, and 239 Pu) as indicated by validating with LLNL pulsed-sphere experiments.

14-MeV D+T LLNL pulsed-sphere neutron-leakage spec↗
DOE OSTI2021

Improvements to the macroscopic-microscopic approach of nuclear fission

The well-established macroscopic-microscopic (mac-mic) description of nuclear fission enables the prediction of fission-fragment yields for a broad range of fissioning systems. In this work, we present several key enhancements to this approach. We improve upon the microscopic sector of nuclear potential-energy surfaces by magnifying the resolution of the Lipkin-Nogami equations and strengthening the Strutinsky procedure, thus reducing spurious effects from the continuum. We further present a novel deterministic method for calculating fission dynamics under the assumption of strongly damped nucleonic motion. Our technique directly determines the evolution of the scissioned shape distribution according to the number of random-walk steps rather than the statistical accumulation of fission events. We show that our new technique is equivalent to the Metropolis random walk pioneered over the past decade by Randrup and colleagues. It further improves upon it because we remove the need for altering the nuclear landscape via a biased potential. With our final improvement, we calculate fission fragments mass and charge distributions using particle number projection, which affords the simultaneous calculation of both mass- and charge-yield distributions. Fission fragments are thus calculated from the quantum-mechanical A -body states of the potential-energy surface rather than from the collective mass asymmetry variable α g of the finite-range liquid-drop model used in past work. We highlight the success of our enhancements by predicting the odd-even staggering and the charge polarization for the neutron-induced fission of 233 U and 235 U .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Surrogate-constructed scalable-circuits adaptive variational quantum eigensolver in the Schwinger model

Inspired by recent advancements in simulating periodic systems on quantum computers, we develop an approach to further advance the simulation of these systems, named (SC) 2 -ADAPT-VQE. Our approach extends the scalable-circuits ADAPT-VQE framework, which builds an ansatz from a pool of coordinate-invariant operators defined for arbitrarily large, though not arbitrarily small, volumes. Our method uses a classically tractable “surrogate constructed” method to remove irrelevant operators from the pool, reducing the minimum size for which the scalable circuits are defined. Bringing together the scalable circuits and the surrogate constructed approaches forms the core of the (SC) 2 methodology. Our approach allows for a wider set of classical computations on small volumes, which can be used for a more robust extrapolation protocol. While developed in the context of lattice models, the surrogate construction portion is applicable to a wide variety of problems where information about the relative importance of operators in the pool is available. As an example, we use it to compute the properties of the Schwinger model—quantum electrodynamics for a single, massive fermion in 1 +1 dimensions—and show that our method can be used to accurately extrapolate to the continuum limit.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Enabling event-by-event precision in γ-ray cascades for neutron-induced reactions

Neutron-induced γ-ray spectra provide key inputs for modern active interrogation applications. A precise modeling of the nuclear reaction and subsequent emission of γ rays is challenging and often impossible due to limitations on evaluated data file formats and nuclear transport simulation codes. We present a framework that addresses these challenges by combining experimental data and reaction-model calculation outputs into an extended candidate version of the Generalized Nuclear Data Structure (GNDS) file, the successor format for the legacy Evaluated Nuclear Data File (ENDF-6). This proposed GNDS hierarchical format contains all the necessary ingredients for inline γ-ray cascade reproduction with event-by-event precision, including continuum–continuum and continuum–discrete transitions following neutron-capture and inelastic neutron scattering reactions. Cascade-event generation based on our approach demonstrates improved energy conservation on an event-by-event basis and permits the use of γ-γ coincidences in applications. This work offers, for the first time, a method to generate neutron-capture and inelastic neutron-scattering γ-ray cascades where energy conservation, correlations, and experimental primaries are fully accounted for.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Reexamination of isoscalar giant resonances in 12 C and 93 Nb through 6 Li scattering

Inelastic Li 6 scattering at 100 MeV/u on C 12 and Nb 93 has been measured with the high-resolution magnetic spectrometer Grand Raiden. The magnetic-rigidity settings of the spectrometer covered excitation energies from 10 to 40 MeV and scattering angles in the range 0 ° < θ lab. < 2 ° . The isoscalar giant monopole resonance was selectively excited in the present data. Measurements free of instrumental background and the very favorable resonance-to-continuum ratio of Li 6 scattering allowed for precise determination of the E 0 strengths in C 12 and Nb 93 . Additionally, it was found that the monopole strength in C 12 exhausts 52 ± 3 (stat.) ± 8 (sys.) % of the energy-weighted sum rule (EWSR), which is considerably higher than results from previous α -scattering experiments. The monopole strength in Nb 93 exhausts 92 ± 4 (stat.) ± 10 (sys.) % of the EWSR, and it is consistent with measurements of nuclei with mass number of A ≈ 90 . Such comparison indicates that the isoscalar giant monopole resonance distributions in these nuclei are very similar, and no influence due to nuclear structure was observed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Gauge-symmetrization method for energy-momentum tensors in high-order electromagnetic field theories

For electromagnetic field theories, canonical energy-momentum conservation laws can be derived from the underpinning spacetime translation symmetry according to the Noether procedure. However, the canonical energy-momentum tensors (EMTs) are neither symmetric nor gauge-symmetric (gauge invariant). The Belinfante-Rosenfeld (BR) method is a well-known procedure to symmetrize the EMTs, which also renders them gauge symmetric for first-order field theories. High-order electromagnetic field theories appear in the study of gyrokinetic systems for magnetized plasmas and the Podolsky system for the radiation reaction of classical charged particles. For these high-order field theories, gauge-symmetric EMTs are not necessarily symmetric and vice versa. In the present study, we develop a new gauge-symmetrization method for EMTs in high-order electromagnetic field theories. The Noether procedure is carried out using the Faraday tensor $F_{μν}$, instead of the 4-potential $A_{μ}$, to derive a canonical EMT $T^{μν}_{N}$. We show that the gauge-dependent part of $T^{μν}_{N}$ can be removed using the displacement-potential tensor $F^{σμν}$ ≡ $D^{σμ}A^{ν}/4π$, where $D^{σμ}$ is the antisymmetric electric displacement tensor. This method gauge-symmetrizes the EMT without necessarily making it symmetric, which is adequate for applications not involving general relativity. For first-order electromagnetic field theories, such as the standard Maxwell system, $F^{σμν}$ reduces to the familiar BR superpotential $S^{σμν}$, and the method developed can be used as a simpler procedure to calculate $S^{σμν}$ without employing the angular momentum tensor in 4D spacetime. When the electromagnetic system is coupled to classical charged particles, the gauge-symmetrization method for EMTs is shown to be effective as well.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Evaluation of Compton suppression for enhancing trace element identification in neutron activation analysis of reference materials

Neutron activation analysis (NAA) is a powerful technique for identifying and quantifying trace elements in materials. However, challenges such as high dead times, spectral interferences, and high Compton continuum often arise. This study employs a Compton suppression system (CSS) to enhance NAA sensitivity by reducing the Compton continuum, thereby improving the peak-to-Compton ratio. National Institute of Standards and Technology certified reference materials 1632d, 1633c, and Canadian National Research Council TORT-1 were irradiated in a thermal and epithermal neutron flux under various irradiation, decay and counting times and analyzed using high-resolution gamma-ray spectroscopy with and without Compton suppression. The reduction factor was calculated to evaluate the effectiveness of the CSS, demonstrating significant background reduction and improved detection limits for trace elements. Additional experiments with a 137Cs point source demonstrated the impact of detector-source geometry on system performance, showing a decrease in the reduction factor as the source was moved further from the NaI detector. An optimum distance between the source and HPGe detector was observed, yielding the highest peak-to-Compton ratio. The results highlight the CSS's ability to minimize spectral interference and enhance elemental identification.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

9 Li($\textit{d,p}$) reaction as a specific probe of 10 Li, the paradigm of parity-inverted nuclei around the $\textit{N}$ = 6 closed shell

In this work, we show, within the framework of renormalized nuclear field theory and of the induced reaction surrogate formalism, that the highly debated 10 Li structure, studied in a recent high statistics 9 Li ($\textit{d,p}$) 10 Li one-neutron transfer experiment, is consistent with, or better, requires, the presence of a virtual 1/2 + state of similar single-particle strength than that of the 1/2 – resonance at 0.45 ± 0.03 MeV. Based on continuum spectroscopy self-energy techniques, we find that the physical mechanism responsible for parity inversion in $^{10}_3$Li is the same as that at the basis of the similar phenomenon observed in $^{11}_4$Be and as that needed in 11 Li to have an important $\textit{s}$-wave ground-state component. In particular the strong dynamical coupling between the $s_{1/2}$ and the $d_{5/2}$ states, mediated by the quadrupole vibration of the core 9 Li. A phenomenon which also affects the strength distribution of the $d_{5/2}$ state, in particular, in the energy range of 3–4.5 MeV. Furthermore, this mechanism is also consistent with the (normal) sequence of the $1_{p_{1/2}}$ and $2_{s_{1/2}}$ levels in the $\textit{N}$ = 7 isotones $^{12}_5$B and $^{13}_6$C. The main aim of the present Rapid Communication is that of treating structure and reactions on equal footing and in a common language. In other words, the calculation of the 9 Li ($\textit{d,p}$) 10 Li reaction as a single conceptual step from individual single-particle motion and collective vibrations to absolute double differential cross sections of renormalized virtual and resonant final states, which can be directly compared with experiment.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Classical and quantum computing of shear viscosity for ( 2 + 1 ) D SU(2) gauge theory

We perform a nonperturbative calculation of the shear viscosity for ( 2 + 1 )-dimensional SU(2) gauge theory by using the lattice Hamiltonian formulation. The retarded Green’s function of the stress-energy tensor is calculated from real time evolution via exact diagonalization of the lattice Hamiltonian with a local Hilbert space truncation, and the shear viscosity is obtained via the Kubo formula. When taking the continuum limit, we account for the renormalization group flow of the coupling but no additional operator renormalization. We find the ratio of the shear viscosity and the entropy density η s is consistent with a well-known holographic result 1 4 π at several temperatures on a 4 × 4 honeycomb lattice with the local electric representation truncated at j max = 1 2 . We also find the ratio of the spectral function and frequency ρ x y ( ω ) ω exhibits a peak structure when the frequency is small. Both the exact diagonalization method and simple matrix product state classical simulation method beyond j max = 1 2 on bigger lattices require exponentially growing resources. So we develop a quantum computing method to calculate the retarded Green’s function and analyze various systematics of the calculation including j max truncation and finite size effects, Trotter errors and the thermal state preparation efficiency. Our thermal state preparation method still requires resources that grow exponentially with the lattice size, but with a very small prefactor at high temperature. We test our quantum circuit on both the Quantinuum emulator and the IBM simulator for a small lattice and obtain results consistent with the classical computing ones. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Hubble Space Telescope UV and Hα Measurements of the Accretion Excess Emission from the Young Giant Planet PDS 70 b

Recent discoveries of young exoplanets within their natal disks offer exciting opportunities to study ongoing planet formation. In particular, a planet’s mass accretion rate can be constrained by observing the accretion-induced excess emission. So far, planetary accretion is only probed by the Hα line, which is then converted to a total accretion luminosity using correlations derived for stars. However, the majority of the accretion luminosity is expected to emerge from hydrogen continuum emission, and is best measured in the ultraviolet (UV). In this paper, we present HST/WFC3/UVIS F336W (UV) and F656N (Hα) high-contrast imaging observations of PDS 70. Applying a suite of novel observational techniques, we detect the planet PDS 70 b with signal-to-noise ratios of 5.3 and 7.8 in the F336W and F656N bands, respectively. This is the first time that an exoplanet has been directly imaged in the UV. Our observed Hα flux of PDS 70 b is higher by 3.5σ than the most recent published result. However, the light curve retrieved from our observations does not support greater than 30% variability in the planet’s Hα emission in six epochs over a five month timescale. We estimate a mass accretion rate of 1.4±0.2×10{sup −8}M{sub Jup} yr{sup −1}. Hα accounts for 36% of the total accretion luminosity. Such a high proportion of energy released in line emission suggests efficient production of Hα emission in planetary accretion, and motivates using the Hα band for searches of accreting planets. These results demonstrate HST/WFC3/UVIS’s excellent high-contrast imaging performance and highlight its potential for planet formation studies.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

New narrow resonances observed in the unbound nucleus F 15

The structure of the unbound 15 F nucleus is investigated using the inverse kinematics resonant scattering of a radioactive 14 O beam impinging on a CH 2 target. The analysis of 1 H( 14 O,p) 14 O and 1 H( 14 O,2p) 13 N reactions allowed the confirmation of the previously observed narrow 1/2 - resonance, near the two-proton decay threshold, and the identification of two new narrow 5/2 - and 3/2 - resonances. The newly observed levels decay by 1p emission to the ground of 14 O, and by sequential 2p emission to the ground state (g.s.) of 13 N via the 1 - resonance of 14 O. Gamow shell model (GSM) analysis of the experimental data suggests that the wave functions of the 5/2 - and 3/2 - resonances may be collectivized by the continuum coupling to nearby 2p- and 1p- decay channels. Finally, the observed excitation function 1 H( 14 O, p) 14 O and resonance spectrum in 15 F are well reproduced in the unified framework of the GSM.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Probing the Wind Component of Radio Emission in Luminous High-redshift Quasars

We discuss a probe of the contribution of wind-related shocks to the radio emission in otherwise radio-quiet quasars. Given (1) the nonlinear correlation between UV and X-ray luminosity in quasars, (2) that such a correlation leads to higher likelihood of radiation-line-driven winds in more luminous quasars, and (3) that luminous quasars are more abundant at high redshift, deep radio observations of high-redshift quasars are needed to probe potential contributions from accretion disk winds. We target a sample of 50 z ≃ 1.65 color-selected quasars that span the range of expected accretion disk wind properties as traced by broad C iv emission. 3 GHz observations with the Very Large Array to an rms of ≈10 μJy beam -1 probe to star formation rates of ~400 M ⊙ yr -1 , leading to 22 detections. Supplementing these pointed observations are survey data of 388 sources from the LOFAR Two-meter Sky Survey Data Release 1 that reach comparable depth (for a typical radio spectral index), where 123 sources are detected. These combined observations reveal a radio detection fraction that is a nonlinear function of C iv emission-line properties and suggest that the data may require multiple origins of radio emission in radio-quiet quasars. We find evidence for radio emission from weak jets or coronae in radio-quiet quasars with low Eddington ratios, with either (or both) star formation and accretion disk winds playing an important role in optically luminous quasars and correlated with increasing Eddington ratio. Additional pointed radio observations are needed to fully establish the nature of radio emission in radio-quiet quasars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Quantifying Uncertainties In Neutron-Alpha Scattering With Chiral Nucleon-Nucleon And Three-Nucleon Forces

We report that modern ab initio theory combined with high-quality nucleon-nucleon (NN) and three-nucleon (3N) interactions from chiral effective field theory (EFT) can provide a predictive description of low-energy light-nuclei reactions relevant for astrophysics and fusion-energy applications. However, the high cost of computations has so far impeded a complete analysis of the uncertainty budget of such calculations. Starting from NN potentials up to fifth order (N 4 LO) combined with leading-order 3N forces, we study how the order-by-order convergence of the chiral expansion and confidence intervals for the 3N contact and contact-plus-one-pion-exchange low-energy constants (c E and c D ) contribute to the overall uncertainty budget of many-body calculations of neutron- 4 He (n-α) elastic scattering. We compute structure and reaction observables for three-, four-, and five-nucleon systems within the ab initio frameworks of the no-core shell model and no-core shell model with continuum. Using a small set of design runs, we construct a Gaussian process model (GPM) that acts as a statistical emulator for the theory. With this, we gain insight into how uncertainties in the 3N low-energy constants propagate throughout the calculation and determine the Bayesian posterior distribution of these parameters with Markov-Chain Monte Carlo. We find rapidly converging n-α phase shifts with respect to the chiral order. With the adopted leading-order 3N force, calculations based on the NN interaction at N 4 LO of Entem, Machleidt, and Nosyk are unable to reproduce the experimental phase shifts in the 3/2 - channel within the estimated chiral truncation errors. Closer agreement with empirical data is found when using an older parametrization of the NN interaction at order N 3 LO, and the position and width of the P-wave resonances can be used to reduce the uncertainty of the 3N low-energy constants. The present results point to a lack of spin-orbit strength when the newer parametrization of the chiral NN force up to fifth order is combined with the leading-order 3N force. The inclusion of higher-order 3N-force terms may be required to recover the missing strength. GPMs can act as fast and accurate emulators of ab initio many-body calculations of low-energy scattering and reactions of light nuclei, opening the way to a robust quantification of theoretical uncertainties grounded in the description of the underlying chiral Hamiltonian.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗