Universality of weak interactions and rare pion decays
Explore the source record for details and available documents.
Atomic, Nuclear and Particle Physics↗
SEARCH · Engineering Papers
Results for “nuclear interactions”
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.
Measurement of $\pi^+$ -- Argon Absorption and Charge Exchange Interactions Using ProtoDUNE-SP
Explore the source record for details and available documents.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Quantum Computing for Neutrino-Nucleus Scattering
Neutrino-nucleus cross section uncertainties are expected to be a dominant systematic in future accelerator neutrino experiments. The cross sections are determined by the linear response of the nucleus to the weak interactions of the neutrino, and are dominated by energy and distance scales of the order of the separation between nucleons in the nucleus. These response functions are potentially an important early physics application of quantum computers. Here we present an analysis of the resources required and their expected scaling for scattering cross section calculations. The current estimates of Trotter steps needed to achieve an energy resolution of 10 MeV and the number of CNOT gates for analyzing Ar40 highlights the need for significant improvements in algorithms. We also examine simple small-scale neutrino-nucleus models on modern quantum hardware. In this paper, we use variational methods to obtain the ground state of a three nucleon system (the triton) and then implement the relevant time evolution. To tame the errors in present-day NISQ devices, we explore the use of different error-mitigation techniques to increase the fidelity of the calculations.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Novel machine-learning method for spin classification of neutron resonances
The performance of nuclear reactors and other nuclear systems depends on a precise understanding of the neutron interaction cross sections for materials used in these systems. These cross sections exhibit resonant structure whose shape is determined in part by the angular-momentum quantum numbers of the resonances. The correct assignment of the quantum numbers of neutron resonances is, therefore, paramount. In this project, we apply machine learning to automate the quantum number assignments using only the resonances' energies and widths and not relying on detailed transmission or capture measurements. The classifier used for quantum number assignment is trained using stochastically generated resonance sequences whose distributions mimic those of real data. Here we explore the use of several physics-motivated features for training our classifier. These features amount to out-of-distribution tests of a given resonance's widths and resonance-pair spacings. We pay special attention to situations where either capture widths cannot be trusted for classification purposes or where there is insufficient information to classify resonances by the total spin J. We demonstrate the efficacy of our classification approach using simulated and actual 52 Cr resonance data.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Expansion of Machine-Learning Method for Classifying Neutron Resonances
The understanding of astrophysics processes and the performance of nuclear reactors and other nuclear systems depend on a precise description of the neutron interaction cross sections for materials and nuclei present in these environments. At low neutron energies, these cross sections exhibit resonance structure represented by sharp enhancements when the neutron energy is sufficiently close to excited levels in a compound nucleus. Such resonances can be characterized by their quantum numbers relative to angular momenta, which are often deduced in an ad hoc and irreproducible manner from the shape of the cross sections. The correct assignment of the quantum numbers of neutron resonances is therefore of paramount importance. To address this we have developed a machine-learning method to automate the identification and correction of these spin assignments. The algorithm is trained from simulated data, generated from statistical properties of resonance data for a given nucleus, to mimic the errors found in real data. In this project we describe five independent approaches to further develop and expand the applicability of the machine-learning spin classifier: i) Feature impact; ii) Integration with the Atlas; iii) Training optimization; iv) Spacings systematics; and v) Validation with polarized data. The premises, methods, results, and future perspectives are discussed.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Machine learning applied to classifying neutron resonances
The performance of nuclear reactors and other nuclear systems depends on a precise understanding of the neutron interaction cross sections for materials used in these systems. These cross sections exhibit resonance structure whose shape is determined in part by the angular momentum quantum numbers of the resonances. The correct assignment of the quantum numbers of neutron resonances is therefore of paramount importance. In this project, we apply a machine learning technique, namely decision trees, to automate the quantum number assignments. The tree is trained from simulated data generated to mimic the errors found in real data. We explore the use of several physics-motivated features for training our trees, including the nearest neighbor spacing distribution, cumulative level distribution, and channel width distributions. Initial results using random matrix theory motivated fits which demonstrated that we can determine resonance spin groups somewhat reliably. If we use these fits as features in our trees, we can train them to spot outliers corresponding to misassigned resonances. We found that with the large number of features used in this project that the decision tree tended to over t training data resulting in poor performance with respect to the test data. By reducing the number of features, we can achieve nearly perfect assignment of quantum numbers with our training data.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Neutron-neutron scattering length from the 6 He(p,pα)nn reaction
We propose a novel method to measure the neutron-neutron scattering length using the 6 He(p,pα)nn reaction in inverse kinematics at high energies. The method is based on the final-state interaction (FSI) between the neutrons after the sudden knockout of the α particle. We show that the details of the neutron-neutron relative-energy distribution allow for a precise extraction of the s-wave scattering length. Here, we present the state of the art in regard to the theory of this distribution. The distribution is calculated in two steps. First, we calculate the ground-state wave function of 6 He as a αnn three-body system. For this purpose we use Halo effective field theory, which also provides uncertainty estimates for the results. We compare our results at this stage to model calculations done with the computer code face. In a second step we determine the effects of the nn FSI using the nn t-matrix. We compare these FSI results to approximate FSI approaches based on standard FSI enhancement factors. While the final distribution is sensitive to the nn scattering length, it depends only weakly on the effective range. Throughout we emphasize the impact of theoretical uncertainties on the neutron-neutron relative-energy distribution, and discuss the extent to which those uncertainties limit the extraction of the neutron-neutron scattering length from the reaction 6 He(p,pα)nn.
6 ≤ A ≤ 19↗
Accelerating eigenvalue computation for nuclear structure calculations via perturbative corrections
Subspace projection methods utilizing perturbative corrections have been proposed for computing the lowest few eigenvalues and corresponding eigenvectors of large Hamiltonian matrices. In this paper, we build upon these methods and introduce the term Subspace Projection with Perturbative Corrections (SPPC) method to refer to this approach. We tailor the SPPC for nuclear many-body Hamiltonians represented in a truncated configuration interaction subspace, i.e., the no-core shell model (NCSM). We use the hierarchical structure of the NCSM Hamiltonian to partition the Hamiltonian as the sum of two matrices. The first matrix corresponds to the Hamiltonian represented in a small configuration space, whereas the second is viewed as the perturbation to the first matrix. Eigenvalues and eigenvectors of the first matrix can be computed efficiently. Because of the split, perturbative corrections to the eigenvectors of the first matrix can be obtained efficiently from the solutions of a sequence of linear systems of equations defined in the small configuration space. These correction vectors can be combined with the approximate eigenvectors of the first matrix to construct a subspace from which more accurate approximations of the desired eigenpairs can be obtained. We show by numerical examples that the SPPC method can be more efficient than conventional iterative methods for solving large-scale eigenvalue problems such as the Lanczos, block Lanczos and the locally optimal block preconditioned conjugate gradient (LOBPCG) method. The method can also be combined with other methods to avoid convergence stagnation.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Gas-phase stability of large lanthanide:diglycolamide clusters evaluated using collision-induced dissociation
Efficient utilization of long-term deep geologic storage repositories is critical for the large-scale deployment of nuclear energy. As the minor actinides (americium and curium) are the largest contributor to the decay heat of used nuclear fuel after a few hundred years, their removal from used fuel prior to disposal can significantly increase the amount of waste that can be stored in a given volume. The development of ligands that can effectively separate the minor actinides from other components of used nuclear fuel are crucial to efficient utilization of geological storage repositories. N,N,N’,N’-tetraoctyl diglycolamide (TODGA) is a promising extractant for the separation of lanthanides and the minor actinides from other components of used nuclear fuel. While the use of TODGA in f-element separations has been investigated in process-based formulations, gas-phase metal ion cluster experiments enable the study of covalent interactions in reprocessing systems absent from solvent effects. Exploring fundamental differences in lanthanide-ligand covalent interactions can impact the development and implementation of actinide-lanthanide separation systems in nuclear fuels reprocessing. In this study, lanthanide-TODGA clusters were synthesized in the gas-phase and identified using mass spectrometry fragmentation experiments. Large europium and samarium clusters were identified that contained up to 8 and 10 bound TODGA ligands, respectively; this was surprising due to the size and multidentate binding that is normal for TODGA. In addition, while smaller clusters showed evidence of sequential ligand fragmentation, larger clusters displayed the loss of neutral TODGA with applied collision voltage. Interestingly, the voltage required for this removal decreased as more TODGA ligands were bound to the metal, suggesting that the metal coordination sphere was becoming more saturated and TODGA ligands were more weakly bound as the clusters got larger.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Pursuing the Ultimate Power of Xenon Dark Matter Detectors (Annual Progress Report)
Liquid xenon-based experiments have been leading direct searches for dark matter – a cornerstone of modern cosmology and particle physics. Despite rapid improvement of experimental sensitivities in the past two decades, no definitive dark matter interactions have been observed. This project aims to expand the physics reach of existing and future xenon dark matter experiments, especially for low-mass dark matter interactions that would fall below the energy thresholds of current detectors. The main approach is to thoroughly characterize and to suppress the low-energy electron background observed in dual-phase xenon Time Projection Chambers (TPCs), which has so far prevented these detectors from achieving lower energy thresholds. Per recent discussion with and approval from the program manager, we have added a new task of experimentally measuring the Migdal effect to this project. The Migdal effect predicts that ultra-low energy dark matter interactions may produce detectable electron recoil signals in liquid xenon at the keV level in addition to much lower energy nuclear recoils. If this effect is experimentally measured, it will drastically improve xenon detectors’ sensitivity to subGeV dark matter interactions. This new task shares strong synergy with the original project goal of pursuing the ultimate power of xenon dark matter experiments.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Ab initio leading order effective potentials for elastic nucleon-nucleus scattering
Calculating microscopic effective interactions (optical potentials) for elastic nucleon-nucleus scattering has already in the past led to a large body of work. For first-order calculations a nucleon-nucleon ( NN ) interaction and a one-body density of the nucleus were taken as input to rigorous calculations of microscopic full-folding calculations. Based on the spectator expansion of the multiple scattering series we employ a chiral next-to-next-to-leading order (NNLO) nucleon-nucleon interaction on the same footing in the structure as well as in the reaction calculation to obtai n an in leading-order consistent effective potential for nucleon-nucleus elastic scattering, which includes the spin of the struck target nucleon. The first order effective folding potential is computed by first deriving a nonlocal scalar density as well as a spin-projected momentum distribution. Those are then integrated with the off-shell Wolfenstein amplitudes $A$, $C$, and $M $. The resulting nonlocal potential serves as input to a momentum-space Lippmann-Schwinger equation, whose solutions are summed to obtain the nucleon-nucleus scattering observables. As a result, we calculate elastic scattering observables for $^4$He, $^6$He, $^8$He, $^{12}$C, and $^{16}$O in the energy r egime between 100 and 200~MeV projectile kinetic energy, and compare to available data. We also explore the extension down to about 70 MeV, and study the effect of ignoring the spin of the struck nucleon in the In our calculations we contrast elastic scattering off closed-shell and open-shell nuclei. We find that fo r closed-shell nuclei the approximation of ignoring the spin of the struck target nucleon is excellent. We only see effects of the spin of the struck target nucleon when considering $^6$He and $^8$He, which are nuclei with a $N/Z$ ratio larger than 1.
6 ≤ A ≤ 19↗
Two-nucleon 𝑆-wave interactions at the SU(3) flavor-symmetric point with 𝑚 𝑢𝑑 ≈ 𝑚$^{phys}_𝑠$: A first lattice QCD calculation with the stochastic Laplacian Heaviside method
We report on the first application of the stochastic Laplacian Heaviside method for computing multiparticle interactions with lattice QCD to the two-nucleon system. Like the Laplacian Heaviside method, this method allows for the construction of interpolating operators which can be used to construct a set of positive-definite two-nucleon correlation functions, unlike nearly all other applications of lattice QCD to two nucleons in the literature. It also allows for a variational analysis in which optimal linear combinations of the interpolating operators are formed that couple predominantly to the eigenstates of the system. Utilizing such methods has become of paramount importance to help resolve the discrepancy in the literature on whether two nucleons in either isospin channel form a bound state at pion masses heavier than physical, with the discrepancy persisting even in the SU(3)-flavor-symmetric point with all quark masses near the physical strange quark mass. This is the first in a series of papers aimed at resolving this discrepancy. In the present work, we employ the stochastic Laplacian Heaviside method without a hexaquark operator in the basis at a lattice spacing of 𝑎 ≈0.086 fm, lattice volume of 𝐿 = 48𝑎 ≈ 4.1 fm and pion mass 𝑚 𝜋 ≈ 714 MeV. With this setup, the observed spectrum of two-nucleon energy levels strongly disfavors the presence of a bound state in either the deuteron or dineutron channel.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Bayesian analysis of nucleon-nucleon scattering data in pionless effective field theory
We perform Bayesian model calibration of two-nucleon (NN) low-energy constants (LECs) appearing in an NN interaction based on pionless effective field theory (πEFT). The calibration is carried out for potentials constructed using naive dimensional analysis in NN relative momenta (p) up to next-to-leading order [NLO, O(p 2 )] and next-to-next-to-next-to-leading order [N3LO, O(p 4 )]. We consider two classes of pionless πEFT potential: one that acts in all partial waves and another that is dominated by s-wave physics. The two classes produce broadly similar results for calibrations to NN data up to E lab = 5 MeV. Our analysis accounts for the correlated uncertainties that arise from the truncation of the pionless πEFT. We simultaneously estimate both the πEFT LECs and the parameters that quantify the truncation error. This permits the first quantitative estimates of the pionless πEFT breakdown scale, Λ b : the 95% intervals are Λ b ∈[50.11,63.03] MeV at NLO and Λ b ∈[72.27,88.54] MeV at N3LO. Furthermore, invoking naive dimensional analysis for the NN potential, therefore, does not lead to consistent results across orders in pionless πEFT. This exemplifies the possible use of Bayesian tools to identify inconsistencies in a proposed EFT power counting.
Bayesian methods↗
High-Precision Spectroscopy of 20 O Benchmarking Ab Initio Calculations in Light Nuclei
The excited states of unstable 20 O were investigated via γ-ray spectroscopy following the 19 O(d, p) 20 O reaction at 8 A MeV. By exploiting the Doppler shift attenuation method, the lifetimes of the $2^+_2$ and $3^+_1$ states were firmly established. From the γ-ray branching and E2/M1 mixing ratios for transitions deexciting the $2^+_2$ and $3^+_1$ states, the B(E2) and B(M1) were determined. Various chiral effective field theory Hamiltonians, describing the nuclear properties beyond ground states, along with a standard USDB interaction, were compared with the experimentally obtained data. Such a comparison for a large set of γ-ray transition probabilities with the valence space in medium similarity renormalization group ab initio calculations was performed for the first time in a nucleus far from stability. It was shown that the ab initio approaches using chiral effective field theory forces are challenged by detailed high-precision spectroscopic properties of nuclei. Finally, the reduced transition probabilities were found to be a very constraining test of the performance of the ab initio models.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Fallout Cloud Regimes
The U.S. Department of Defense (DOD), Department of Energy (DOE), and other organizations maintain operational nuclear explosion and atmospheric dispersion models to provide critical guidance on the expected effects of an accidental or deliberate explosion of a nuclear weapon (in this paper simply referred to as “device”). To be effective, these models must represent, as accurately as possible, the complex interactions of the blast, fire, and residual radiological hazards with the environment and population. One hundred atmospheric nuclear tests that form the basis for many models were conducted at the Nevada Test Site (NTS) (now referred to as the Nevada Nuclear Security Site, NNSS) in a dry desert environment. Other environments should be studied, but have less data available and are beyond the scope of the work presented in this paper. The debris clouds produced by the NTS tests, frequently called “mushroom clouds,” are familiar, with common structural elements such as a buoyant cap connected to a skirt of raised dust at the desert surface by a thin, dirt-filled stem. The film scanning project at LLNL has investigated historical film records of nuclear weapons tests. Here, we summarize findings showing that the mushroom cloud behavior for historic U.S. tests conducted in Nevada, has similar characteristics based on the distance of the device from the ground surface or Height of Burst (HOB), scaled by the energy release, or yield, of the device. This scaled height is referred to as the scaled-height-of-burst (SHOB). The findings discussed below show that mushroom clouds look and behave similarly when detonated at the same SHOB. The amount of residual radiation that is produced by a nuclear detonation is proportional to the yield. But, the amount of that residual radiation that actually becomes local fallout is strongly dependent on the SHOB and the type of surface over which the detonation occurs. In order to develop a more comprehensive model that predicts the fraction of the residual radiation that becomes local fallout, it is convenient to define a series of regimes based on SHOB values in which all detonations that occur within a given regime can be modeled using the same algorithms. The purpose of this paper is to provide a framework for defining different regimes, and, in a qualitative way, a basic understanding of the fundamental characteristics of each of these regimes.
45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
Grain boundary self-diffusion and point defect interactions in α -U via molecular dynamics
Though metallic U-Zr fuel has been used in nuclear reactors since the 1960s, many of its fundamental and thermodynamic properties are still unknown. The a-U phase, which has a highly anisotropic crystal structure and physical properties, is present in U-Zr fuel. The character and behavior of a-U grain boundaries will strongly impact fuel thermophysical performance under irradiation. Here, we study the interaction of point defects with grain boundaries, diffusion along grain boundaries, and the predicted diffusional creep behavior of a-U via molecular dynamics. We calculate the segregation energy of vacancies and interstitials to grain boundaries and quantify the biased sink strength of the grain boundaries, and observe that this sink strength is not strongly dependent on the grain boundary orientation. We also find that grain boundary diffusivity is strongly dependent on the grain boundary energy and grain boundary orientation. The presence of point defects within the grain boundary can induce diffusion in grain boundaries with low formation energies and can enhance diffusion in high-energy grain boundaries. We also find that diffusional creep of a-U at prototypical metallic fuel operation conditions is extremely high and could help explain observed metallic fuel swelling behaviors.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Identification of preferred multimodal ligand‐binding regions on IgG1 F C using nuclear magnetic resonance and molecular dynamics simulations
Abstract In this study, the binding of multimodal chromatographic ligands to the IgG1 F C domain were studied using nuclear magnetic resonance and molecular dynamics simulations. Nuclear magnetic resonance experiments carried out with chromatographic ligands and a perdeuterated 15 N‐labeled F C domain indicated that while single‐mode ion exchange ligands interacted very weakly throughout the F C surface, multimodal ligands containing negatively charged and aromatic moieties interacted with specific clusters of residues with relatively high affinity, forming distinct binding regions on the F C . The multimodal ligand‐binding sites on the F C were concentrated in the hinge region and near the interface of the C H 2 and C H 3 domains. Furthermore, the multimodal binding sites were primarily composed of positively charged, polar, and aliphatic residues in these regions, with histidine residues exhibiting some of the strongest binding affinities with the multimodal ligand. Interestingly, comparison of protein surface property data with ligand interaction sites indicated that the patch analysis on F C corroborated molecular‐level binding information obtained from the nuclear magnetic resonance experiments. Finally, molecular dynamics simulation results were shown to be qualitatively consistent with the nuclear magnetic resonance results and to provide further insights into the binding mechanisms. An important contribution to multimodal ligand‐F C binding in these preferred regions was shown to be electrostatic interactions and π–π stacking of surface‐exposed histidines with the ligands. This combined biophysical and simulation approach has provided a deeper molecular‐level understanding of multimodal ligand–F C interactions and sets the stage for future analyses of even more complex biotherapeutics.
Gudhka, Ronak B.↗
Verification of numerical models for seismic fluid-structure interaction analysis of internal components in liquid-filled advanced reactors
Earthquake shaking of a liquid-filled advanced reactor induces fluid-structure interaction (FSI) between the reactor vessel (tank), its internal components, and the contained liquid. Seismic design, qualification, and risk assessment of such reactors must consider fluid-structure responses, their geometries and support conditions, and three-directional seismic inputs, all of which require numerical simulations. Numerical models used for analysis of safety-related nuclear equipment must be verified and validated. Here, this paper verifies numerical models of submerged components using analytical solutions. Seismic FSI analysis of the numerical models is performed using the Arbitrary Lagrangian-Eulerian (ALE) and Incompressible Computational Fluid Dynamics (ICFD) solvers in LS-DYNA. Prior analytical solutions are reworked in this paper, and calculation errors are identified and corrected. The solutions address frequencies of two concentric cylindrical pipes filled with liquid: the inner pipe is a submerged component, and the outer pipe contains the liquid. A unitless frequency coefficient is defined and presented for pipes of different materials, filled with different liquids, and with a range of dimensions for application to advanced reactors. The numerical models are verified here by comparing the lateral frequencies of submerged components with those calculated using the corrected analytical solutions. Recommendations for verification of numerical models of internal components in advanced reactors for seismic FSI analysis are provided. Although the unitless frequency coefficients and verification procedures are developed for application to advanced nuclear reactors, they are broadly applicable to FSI analysis of submerged components in liquid-filled vessels such as storage tanks, boilers, and steam generators.
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗