SEARCH · Search NASA

Results for “Reproducibility”

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.

SourceResearch connectionFrom yearThrough year

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

At least 937 records · Page 52

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Windowed multipole representation of R -matrix cross sections

Nuclear cross sections are basic inputs to any nuclear computation. Campaigns of experiments are fitted with the parametric R-matrix model of quantum nuclear interactions, and the resulting cross sections are documented—both pointwise and as resonance parameters (with uncertainties)—in standard evaluated nuclear data libraries (ENDF, JEFF, BROND, JENDL, CENDL, TENDL): these constitute our common knowledge of fundamental low-energy nuclear cross sections. In the past decade, a collaborative effort has been deployed to establish a new nuclear cross-section library format—the Windowed Multipole Library—with the goal of considerably reducing the computational cost of cross-section calculations in nuclear transport simulations. This work lays the theoretical foundations underpinning these efforts. From general R-matrix scattering theory, we derive the windowed multipole representation of nuclear cross sections. Though physically and mathematically equivalent to R-matrix cross sections, the windowed multipole representation is particularly well suited for subsequent temperature treatment of angle-integrated cross sections, in particular Doppler broadening, which is the averaging of cross sections over the thermal motion of the target atoms. Doppler broadening is of critical importance in neutron transport applications, as it ensures the stability of many nuclear reactors (negative thermal reactivity). Yet, Doppler broadening of nuclear cross sections has been a considerable bottleneck for nuclear transport computations, often requiring memory-costly pretabulations. We show that the windowed multipole representation can perform accurate Doppler broadening analytically (up to the first reaction threshold), from which we derive cross-section temperature derivatives to any order—all computable on the fly (without precalculations stored in memory). Furthermore, we here establish a way of converting the R-matrix resonance parameters uncertainty (covariance matrices) into windowed multipole parameters uncertainty. We show that generating stochastic nuclear cross sections by sampling from the resulting windowed multipole covariance matrix can reproduce the cross-section uncertainty in the original nuclear data file. The windowed multipole representation is therefore a novel nuclear physics formalism able to generate Doppler broadened stochastic nuclear cross sections on the fly, unlocking breakthrough computational gains for nuclear computations. Through this foundational paper, we hope to make the windowed multipole representation accessible, reproducible, and usable for the nuclear physics community, as well as provide the theoretical basis for future research on expanding its capabilities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Differential two-particle number and momentum correlations with the AMPT, UrQMD, and EPOS models in Pb-Pb collisions at s N N = 2.76 TeV

We report studies of charge-independent and charge-dependent two-particle differential number correlation functions, R 2 (Δη,Δφ), and transverse momentum correlation functions, P 2 (Δη,Δφ), of charged particles produced in Pb–Pb collisions at the LHC center-of-mass energy √sNN= 2.76 TeV with the UrQMD, AMPT, and EPOS models. Model calculations for R2 and P2 correlation functions are presented for inclusive charged hadrons (h ± ) in selected transverse momentum ranges and with full azimuthal coverage in the pseudorapidity range |η|<1.0. We compare these calculations for the strength, shape, and particularly the width of the correlation functions with recent measurements of these observables by the ALICE collaboration. Our analysis indicates that comparative studies of R 2 and P 2 correlation functions provide valuable insight towards the understanding of particle production in Pb–Pb collisions. We find, in particular, that these models produce quantitatively different magnitudes and shapes for these correlation functions and none reproduce the results reported by the ALICE collaboration. Accounting for quantum number conservation in models, particularly charge conservation, is mandatory to reproduce the detailed measurements of number and transverse momentum correlation functions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Applications of reduced-basis methods to the nuclear single-particle spectrum

Reduced-basis methods provide a powerful framework for building efficient and accurate emulators. Although widely applied in many fields to simplify complex models, reduced-basis methods have only been recently introduced into nuclear physics. In this Letter we build an emulator to study the single-particle structure of atomic nuclei. By scaling a suitable mean-field Hamiltonian, a “universal” reduced basis is constructed capable of accurately and efficiently reproduce the entire single-particle spectrum of a variety of nuclei. Indeed, the reduced-basis model reproduces both ground- and excited-state energies as well as the associated wave functions with remarkable accuracy. Here our results bode well for more demanding applications that use Bayesian optimization to calibrate nuclear energy density functionals.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Secondary neutron production from thick target deuteron breakup

Thick target deuteron breakup is a variable-energy accelerator-based source of high-energy neutrons, with applications in fundamental and applied nuclear science and engineering. However, the breakup mechanism remains poorly understood, and data on neutron yields from thick target breakup remain relatively scarce. In this work, the double-differential neutron yields from deuteron breakup have been measured on a thick beryllium target at ϵ d = 33 and 40 MeV, using both time-of-flight and activation techniques. We have also introduced a simple hybrid model for the double-differential deuteron breakup cross section, applicable in the ϵ d = 10-100 MeV energy range on light (Z ≤ 6) targets. This model features four empirical parameters that have been fit to reproduce experimental breakup measurements on beryllium targets, using the method of least squares. It was shown that these parameters extrapolate well to higher energies, and to other low-Z target materials. Here, we also include optimization of the parameters that modify the Kalbach systematics for compound and preequilibrium reactions, in order to better reproduce the experimental data for beryllium targets at large angles.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Density dependence of the symmetry energy in the post–PREX-CREX era

The recently published CREX results suggest a rather peculiar picture for the density dependence of the symmetry energy. Whereas PREX favors a large neutron skin thickness in 208 Pb, thereby suggesting a stiff equation of state, CREX suggests instead a much softer equation of state. Here, this discrepancy has caused a large spur in the theoretical community, since no model has been able to simultaneously reproduce within 1σ the PREX and CREX results. Motivated by a novel correlation between a CREX observable and a combination of bulk symmetry energy parameters, we calibrate three new covariant energy density functionals that reproduce binding energies and charge radii of spherical nuclei—and also accommodate the constraints imposed by PREX and CREX. Given that these models suggest a stiff equation of state at high densities, predictions for neutron star properties are also discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Single-channel and single-energy partial-wave analysis with continuity improved through minimal phase constraints

Single-energy partial-wave analysis has often been applied as a way to fit data with minimal model dependence. However, remaining unconstrained, partial waves at neighboring energies will vary discontinuously because the overall amplitude phase cannot be determined through single-channel measurements. This problem can be mitigated through the use of a constraining penalty function based on an associated energy-dependent fit. However, the weight given to this constraint results in a biased fit to the data. In this paper, for the first time, we explore a constraining function which does not influence the fit to data. The constraint comes from the overall phase found in multichannel fits which, in the present study, are the Bonn-Gatchina and Jülich-Bonn multichannel analyses. The data are well reproduced and weighting of the penalty function does not influence the result. The method is applied to K⁢Λ photoproduction data and all observables can be maximally well reproduced. While the employed multichannel analyses display very different multipole amplitudes, we show that the major difference between two sets of multipoles can be related to the different overall phases.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Magnetic dipole γ-ray strength functions in the crossover from spherical to deformed neodymium isotopes

We calculate the magnetic dipole $\gamma$-ray strength functions in a chain of even-mass neodymium isotopes $^{144-152}$Nd in the framework of the configuration-interaction (CI) shell model. We infer the strength function by applying the maximum entropy method (MEM) to the exact imaginary-time response function calculated with the shell-model Monte Carlo (SMMC) method. The success of the MEM depends on the choice of a good strength function as a prior distribution. We investigate two choices for the prior strength function: the static path approximation (SPA) and the quasiparticle random-phase approximation (QRPA). We find that the QRPA is a better approximation at low temperatures (i.e., near the ground state), while the SPA is a better choice at finite temperatures. We identify a low-energy enhancement (LEE) in the MEM deexcitation $M1$ strength functions of the even-mass neodymium isotopes and compare with recent experimental results for the total deexcitation $\gamma$-ray strength functions. The LEE is already seen in the SPA strength function but not in the QRPA strength function, indicating the importance of large-amplitude static fluctuations around the mean field in reproducing the LEE. Our method is currently the only one which can reproduce LEE in heavy open-shell nuclei where conventional CI shell model calculations are prohibited. With the onset of deformation as number of neutrons increases along the chain of neodymium isotopes, we observe that some of the LEE strength transfers to a low-energy excitation, which we interpret as a finite-temperature ``scissors'' mode. Here, we also observe a finite-temperature spin-flip mode.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Negative components in the representation $\frac{𝑑^{2}⁡(𝜎⁢𝐸)}{𝑑⁡𝐸^2}$ of heavy-ion fusion cross sections

Negative components are observed in the representation $\frac{𝑑^{2}⁡(𝜎⁢𝐸)}{𝑑⁡𝐸^2}$ in heavy-ion fusions related to recent studies of noticeable oscillations in the high-energy region of the fusion excitation functions. A modified multiGaussian model has been developed, which includes contributions for both positive and negative components in the $\frac{𝑑^{2}⁡(𝜎⁢𝐸)}{𝑑⁡𝐸^2}$ representation. The model reproduces the heavy-ion fusion excitation functions very well, including oscillations. The property of these negative components in the $\frac{𝑑^{2}⁡(𝜎⁢𝐸)}{𝑑⁡𝐸^2}$ spectrum and the reaction mechanism of fusion at energies above the Coulomb barrier are investigated. Finally, the model also reproduces well light heavy-ion fusion reactions between 12 C and 16 O at high energy; results imply that the compound nucleus channel effect may have an important contribution at the high energy region for all fusion systems.

low & intermediate energy heavy-ion reactions↗

DOE OSTI2017

Anisotropic fluid dynamics for Gubser flow

Investigating a variety of closing schemes to the infinite hierarchy of momentum moments of the exactly solvable Boltzmann equation for systems undergoing Gubser flow, we study the precision with which the resulting hydrodynamic equations reproduce the exact evolution of hydrodynamic moments of the distribution function. We find that anisotropic hydrodynamics, obtained by expanding the distribution function around a dynamically evolving locally anisotropic background whose evolution is matched to exactly reproduce the macroscopic pressure anisotropy caused by the different longitudinal and transverse expansion rates in Gubser ow, provides the most accurate macroscopic description of the microscopic kinetic evolution. This confirms a similar earlier finding for Bjorken ow [Molnar, Niemi and Rischke, Phys. Rev. D 94, 125003 (2016)]. Also, we explain the physics behind this optimal matching procedure and show that one can efficiently correct for a nonoptimized matching choice by adding a residual shear stress to the energy-momentum tensor whose evolution is again determined by the Boltzmann equation. Further insights to guide the optimal choice of a macroscopic anisotropic hydrodynamic framework for strongly-coupled systems that do not admit a microscopic kinetic description are reported.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2017

Neutron transfer in the C 13 + Au 197 reaction from gold isotope residuals

Residual gold nuclei were produced in this paper via neutron transfer at multiple energies using a 130-MeV 13 C beam incident on a stacked-foil target consisting of alternating layers of 197 Au and 27 Al. Production cross sections, over an energy range of 56 to 129 MeV, for seven gold isotopes and two gold isomers were determined through activation analysis. By using the Wilczynski binary transfer model with a modified version of the recoil formula and a standard evaporation model, we were able to reproduce the isotopic production cross sections at high beam energy, with some disagreement at lower beam energies. Finally, this limiting angular momentum model does not predict the transfer of sufficient angular momentum to reproduce the observed isomeric populations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2018

Δ isobars and nuclear saturation

In this paper, we construct a nuclear interaction in chiral effective field theory with explicit inclusion of the $\mathrm{{\Delta}}$-isobar $\mathrm{{\Delta}}(1232)$ degree of freedom at all orders up to next-to-next-to-leading order (NNLO). We use pion-nucleon (${\pi}N$) low-energy constants (LECs) from a Roy-Steiner analysis of ${\pi}N$ scattering data, optimize the LECs in the contact potentials up to NNLO to reproduce low-energy nucleon-nucleon scattering phase shifts, and constrain the three-nucleon interaction at NNLO to reproduce the binding energy and point-proton radius of $^{4}\mathrm{He}$. For heavier nuclei we use the coupled-cluster method to compute binding energies, radii, and neutron skins. We find that radii and binding energies are much improved for interactions with explicit inclusion of $\mathrm{{\Delta}}(1232)$, while $\mathrm{{\Delta}}$-less interactions produce nuclei that are not bound with respect to breakup into ${\alpha}$ particles. Finally, the saturation of nuclear matter is significantly improved, and its symmetry energy is consistent with empirical estimates.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2018

α decay of high- K isomers in Ds 270 and Hs 266 in a superfluid tunneling model

Here, we use the superfluid tunneling model (STM) to calculate the half-lives of ground-state α decays of even-even superheavy nuclei (SHN) with Z ≥ 100 . The experimental data are reproduced to accuracies comparable to other contemporary models of α decay of SHN. We apply the STM to the case of the α decaying high-K isomers identified in the decay chains of 270 Ds . By accounting for the α -decay Q values, Q α , the angular momentum difference between initial and final states, L, and a reduction in the pairing gap, Δ, we are able to reproduce the observed α decay of the isomers, including the unusual competition between L ≈ 10 and L ≈ 0 α branches seen for the K isomer in 270 Ds ( Z = 110 ).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2018

Quantum Monte Carlo formalism for dynamical pions and nucleons

In most simulations of nonrelativistic nuclear systems, the wave functions found solving the many-body Schrodinger equations describe the quantum-mechanical amplitudes of the nucleonic degrees of freedom. In those simulations the pionic contributions are encoded in nuclear potentials and electroweak currents, and they determine the low-momentum behavior. In this work we present an alternative quantum Monte Carlo formalism in which both relativistic pions and nonrelativistic nucleons are explicitly included in the quantum-mechanical states of the system. We report the renormalization of the nucleon mass as a function of the momentum cutoff, a Euclidean time density correlation function that deals with the short-time nucleon diffusion, and the pion cloud density and momentum distributions. In the two-nucleon sector we show that the interaction of two static nucleons at large distances reduces to the one-pion exchange potential, and we fit the low-energy constants of the contact interactions to reproduce the binding energy of the deuteron and two neutrons in finite volumes. We show that the method can be readily applied to light-nuclei.

Madeira, Lucas↗

DOE OSTI2018

Microscopic optical potentials for calcium isotopes

We construct nucleonic microscopic optical potentials by combining the Green's function approach with the coupled-cluster method for 40 Ca and 48 Ca. For the computation of the ground states of 40 Ca and 48 Ca, we use the coupled-cluster method in the singles-and-doubles approximation, while for the A = ±1 nuclei we use particle-attached and particle-removed equation-of-motion method truncated at two-particle–one-hole and one-particle–two-hole excitations, respectively. Our calculations are based on the chiral nucleon-nucleon and three-nucleon interaction NNLO sat , which reproduces the charge radii of 40 Ca and 48 Ca, and the chiral nucleon-nucleon interaction NNLO opt . In all cases considered here, we observe that the overall form of the neutron scattering cross section is reproduced for both interactions, but the imaginary part of the potential, which reflects the loss of flux in the elastic channel, is negligible. Furthermore, the latter points to neglected many-body correlations that would appear beyond the coupled-cluster truncation level considered in this work. We show that, by artificially increasing the parameter η in the Green's function, practical results can be further improved.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2019

Assessing saturation physics explanations of collectivity in small collision systems with the ip-jazma model

Experimental measurements in relativistic collisions of small systems from $p + p$ to $p/d/ ^3He + A$ at the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC) reveal particle emission patterns that are strikingly similar to those observed in $A + A$ collisions of large nuclei. One explanation of these patterns is the formation of small droplets of quark-gluon plasma (QGP) followed by hydrodynamic evolution. A geometry engineering program was proposed to further investigate these emission patterns, and the experimental data from that program in $p + Au, d + Au, ^3 He + Au$ collisions for elliptic and triangular anisotropy coefficients $v_2$ and $v_3$ follow the pattern predicted by hydrodynamic calculations [C. Aidala et al. (PHENIX Collaboration), Nat. Phys. 15, 214 (2019)]. One alternative approach, referred to as initial-state correlations, suggests that for small systems the patterns observed in the final-state hadrons are encoded at the earliest moments of the collision and therefore require no final-state parton scattering or hydrodynamic evolution. Recently, new calculations using only initial-state correlations, in the dilute-dense approximation of gluon saturation physics, reported striking agreement with the $v_2$ patterns observed in $p/d/ ^3He + Au$ data at RHIC [M. Mace, V. V. Skokov, P. Tribedy, and R. Venugopalan, Phys. Rev. Lett. 121, 052301 (2018)]. The results reported by Mace, Skokov, Tribedy and Venugopalan (MSTV) are counterintuitive and thus we aim here to reproduce some of the basic features of these calculations. In this first investigation, we provide a description of our publicly available model, ip-jazma, and investigate its implications for saturation scales, multiplicity distributions, and eccentricities, reserving for later work the analysis of momentum spectra and azimuthal anisotropies. We find that our implementation of the saturation physics model reproduces the results of the MSTV calculation of the multiplicity distribution in $d + Au$ collisions at RHIC. However, additional aspects of studies, together with existing data, call into question some of the essential elements reported by MSTV. Resolution of these issues will require further developments of ip-jazma, in order to determine if it can replicate the qualitative agreement with the $v_2$ reported by MSTV. Both the work reported here and future studies will establish which features in the experimental data are uniquely attributable to the color glass condensate description.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Investigating the XENON1T low-energy electronic recoil excess using NEST

The search for dark matter, the missing mass of the Universe, is one of the most active fields of study within particle physics. The XENON1T experiment recently observed a 3.5⁢𝜎 excess potentially consistent with dark matter, or with solar axions. Here, we will use the Noble Element Simulation Technique (NEST) software to simulate the XENON1T detector, reproducing the excess. We utilize different detector efficiency and energy reconstruction models, but they primarily impact sub-keV energies and cannot explain the XENON1T excess. However, using NEST, we can reproduce their excess in multiple, unique ways, most easily via the addition of 31 ± 11 37 Ar decays. Furthermore, this results in new, modified background models, reducing the significance of the excess to ≤2.2⁢𝜎 at least using non-Profile Likelihood Ratio (PLR) methods. This is independent confirmation that the excess is a real effect, but potentially explicable by known physics. Many cross-checks of our 37 Ar hypothesis are presented.

particle dark matter↗

DOE OSTI2022

Dark Energy Survey Year 3 results: Cosmology from combined galaxy clustering and lensing validation on cosmological simulations

In this work, we present a validation of the Dark Energy Survey Year 3 (DES Y3) 3×2-point analysis choices by testing them on Buzzard2.0, a new suite of cosmological simulations that is tailored for the testing and validation of combined galaxy clustering and weak-lensing analyses. We show that the buzzard2.0 simulations accurately reproduce many important aspects of the DES Y3 data, including photometric redshift and magnitude distributions, and the relevant set of two-point clustering and weak-lensing statistics. We then show that our model for the 3×2-point data vector is accurate enough to recover the true cosmology in simulated surveys assuming the true redshift distributions for our source and lens samples, demonstrating robustness to uncertainties in the modeling of the nonlinear matter power spectrum, nonlinear galaxy bias, and higher-order lensing corrections. Additionally, we demonstrate for the first time that our photometric redshift calibration methodology, including information from photometry, spectroscopy, clustering cross-correlations, and galaxy–galaxy lensing ratios, is accurate enough to recover the true cosmology in simulated surveys in the presence of realistic photometric redshift uncertainties.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Unified description of DGLAP, CSS, and BFKL evolution: TMD factorization bridging large and small x

This paper introduces a transverse-momentum dependent (TMD) factorization scheme designed to unify both large and small Bjorken-x regimes. We compute the next-to-leading order quantum chromodynamics (QCD) corrections to the gluon TMD operator for an unpolarized hadron within this proposed scheme. This leads to the emergence of a new TMD evolution, incorporating those in transverse momentum, rapidity, and Bjorken-x. When matched to the collinear factorization scheme, our factorization scheme faithfully reproduces the well-established Dokshitzer-Gribov-Lipatov-Altarelli-Parisi and Collins-Soper-Sterman (CSS) evolutions. Conversely, matching with high-energy factorization not only yields the Balitsky-Fadin-Kuraev-Lipatov evolution but also reveals distinctive signatures of CSS logarithms. The development of this novel TMD factorization scheme, capable of seamlessly reconciling disparate Bjorken-x regimes and faithfully reproducing established QCD evolution equations, has the potential to significantly advance our comprehension of high-energy processes and three-dimensional parton structures of hadrons. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗