Lattice QCD at LANL (2022 PR) [Slides]
Highlights of research using LANL Institutional Computing resources under allocation Y22_nme.
SEARCH · Engineering Papers
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.
Highlights of research using LANL Institutional Computing resources under allocation Y22_nme.
The competing 22 Ne(α, γ) 26 Mg and 22 Ne(α,n) 25 Mg reactions control the production of neutrons for the weak s process in massive and asymptotic giant branch (AGB) stars. In both systems, the ratio between the corresponding reaction rates strongly impacts the total neutron budget and strongly influences the final nucleosynthesis. A number of experimental studies have been performed over recent years which necessitate the reevaluation of the 26 Mg and 22 Ne(α, γ) 25 Mg reactions control the production 26 Mg and 22 Ne(α, n) 25 Mg reaction rates. Evaluations of the reaction rates following the collection of new nuclear data presently show differences of up to a factor of 500, resulting in considerable uncertainty in the resulting nucleosynthesis. To reevaluate the 22 Ne(α, γ) 26 Mg and 22 Ne(α,n) 25 Mg reaction rates using updated nuclear data from a number of sources including updating spin and parity assignments. With updated spin and parity assignments, the levels which can contribute to the reaction rates are identified. The reaction rates are computed using a Monte Carlo method which has been used for previous evaluations of the reaction rates in order to focus solely on the changes due to modified nuclear data. The evaluated 22 Ne(α, γ) 26 Mg reaction rate remains substantially similar to that of Longland et al. but, including recent results from Texas A&M, the 22 Ne(α, n) 25 Mg reaction rate is lower at a range of astrophysically important temperatures. Stellar models computed with newton and mesa predict decreased production of the weak branch s process due to the decreased efficiency of 22 Ne as a neutron source. Using the new reaction rates in the mesa model results in 96 Zr / 94 Zr and 135 Ba / 136 Ba ratios in much better agreement with the measured ratios from presolar SiC grains. The 22 Ne + α reaction rates 22 Ne(α, γ) 26 Mg and 22 Ne(α, n) 25 Mg have been recalculated based on more recent nuclear data. The 22 Ne(α, γ) 26 Mg reaction rate remains substantially unchanged since the previous evaluation but the 22 Ne(α, γ) 25 Mg reaction rate is substantially decreased due to updated nuclear data. This results in significant changes to the nucleosynthesis in the weak branch of the s process.
Unexplained variations of the decay rate parameter for weak interaction decays such as β ± decay, electron capture, as well as strong interaction α decay have been reported. Because these variations have been presented by a number of groups, at various locations, using various types of detectors, different isotopes, and over extended periods of time, some researchers have interpreted the source of these variations as not from ambient environmental factors such as temperature, pressure, and humidity but via an unexplained fundamental interaction. The state of decay rate parameter variations experiments is reviewed, and the reported results are placed into a common comparable context by defining a phenomenological cross section. After that we make decay parameter measurements as a function of time at the level of 10 -5 in the presences of an antineutrino flux with an on-off cycle time of ≈ 30 days. This level of precision requires a detailed understanding of both systematic and statistical errors; otherwise, systematic errors in the measurement may mimic the decay events of fundamental interactions. The experiment searched for variation of the $^{54}_{25}$Mn, e - capture decay rate parameter and $^{137}_{55}$Cs, β - decay rate parameter both to a level of precision of 1 part in ≈10 5 by comparing the difference between the decay rate in the presence of an antineutrino flux ≈ 3 x 10 12 $\bar{ν}$ cm -2 s -1 and no flux measurements. The experiment was located 6.53 m from the reactor core of the High Flux Isotope Reactor at the Oak Ridge National Laboratory. Two weak interaction decays, one via electron capture and the other via β - decay, were selected because the final state and the time reverse state each contain a neutrino and antineutrino, covering arguments that the antineutrino flux may interact differently or not at all in one of the cases. The $\gamma$ spectra from both decays were collected and analyzed independently. The measured variation in the decay rate parameters are found to be δλ/λ =(0.034 ± 1.38) x10 -5 for $^{54}_{25}$Mn and δλ/λ =(0.67±1.56) x10 -5 for $^{137}_{55}$Cs. These results are consistent with no measurable decay rate parameter variation due to an antineutrino flux, yielding a 68% confidence level upper limit sensitivity for $^{54}_{25}$Mn δλ/λ ≤ 1.31 x 10 -5 or σ ≤ 1.29 x10 -25 cm 2 in cross section and for $^{137}_{55}$Cs δλ/λ ≤ 2.23 x 10 -5 or σ ≤ 5.69 x10 -27 cm 2 . The cross-section upper limit obtained in these null or no observable effect measurements are ≈ 10 4 times more sensitive than past experiments reporting positive results in 54 Mn and 137 Cs. . .
Indirect methods have become the predominant approach in experimental nuclear astrophysics for studying several low-energy nuclear reactions occurring in stars, as direct measurements of many of these relevant reactions are rendered infeasible due to their low reaction probability. Such indirect methods, however, require theoretical input that in turn can have significant poorly quantified uncertainties, which can then be propagated to the reaction rates and have a large effect on our quantitative understanding of stellar evolution and nucleosynthesis processes. Here we present two such examples involving α-induced reactions, 13 C (α,n) 16 O and 12 C (α,γ) 16 O, for which the low-energy cross sections have been constrained with ( 6 Li,d) transfer data. In this Letter, we discuss how a first-principle calculation of 6 Li leads to a 21% reduction of the 12 C(α,γ) 16 O cross sections with respect to a previous estimation. This calculation further resolves the discrepancy between recent measurements of the 13 C (α,n) 16 O reaction and points to the need for improved theoretical formulations of nuclear reactions.
Sterile neutrinos are natural extensions to the standard model of particle physics and provide a possible portal to the dark sector. We report a new search for the existence of sub-MeV sterile neutrinos using the decay-momentum reconstruction technique in the decay of 7 Be. The experiment measures the total energy of the 7 Li daughter atom from the electron capture decay of 7 Be implanted into sensitive superconducting tunnel junction (STJ) quantum sensors. Here, this first experiment presents data from a single STJ operated at a low count rate for a net total of 28 days, and provides exclusion limits on sterile neutrinos in the mass range from 100 to 850 keV that improve upon previous work by up to an order of magnitude.
In transport calculations, it is well known how S{sub N} method is extremely inefficient in problems where the particle physics is dominated by streaming. The ray-effect eventually produced by the insufficient angular discretization, appears to be extremely persistent with respect to the refinement of the angular quadrature. The MP{sub N} method, that relies on continuous angular representation, offers a robust remedy to such an issue. MP{sub N} is based on the decomposition of the unit sphere into solid angles and on a piecewise continuous definition of interface fluxes, which are expanded in polynomials in each solid angle. This allows propagating more than one angular degree of freedom simultaneously while maintaining unaltered the block-diagonal pattern of the displacement plus removal operator. The method is therefore well suited for the flux resolution by means of a conventional sweep algorithm. Furthermore, unlike the S{sub N} method, MP{sub N} does not rely on discrete directions and, thus, on an angular quadrature formula, but rather constructs a set of linear equations solving for the angular moments of the flux for all discrete solid angles within the sweep. MP{sub N} shows an error convergence rate higher than S{sub N} at the expense of an increased size of the coefficient matrices, so of the computational cost. Although MP{sub N} is not free from ray-effect, the latter is effectively mitigated and less persistent with respect to the increase of the angular refinement order. (authors)
The thorium fuel cycle is emerging as an attractive alternative to conventional nuclear fuel cycles, as it does not require the enrichment of uranium for long-term sustainability. The operating principle of this fuel cycle is the irradiation of 232 Th to produce 233 U, which is fissile and sustains the fission chain reaction. 233 U poses unique challenges for nuclear safeguards, as it is associated with a uniquely extreme γ-ray environment from 232 U contamination, which limits the feasibility of the γ-ray-based assay, as well as more conservative accountability requirements than for 235 U set by the International Atomic Energy Agency. Consequently, instrumentation used for safeguarding 235 U in traditional fuel cycles may be inapplicable. It is essential that the nondestructive signatures of 233 U be characterized so that nuclear safeguards can be applied to thorium fuel-cycle facilities as they come online. In this work, a set of 233 U 3 O 8 plates, containing 984 g 233 U, was measured at the National Criticality Experiments Research Center. A high-pressure 4 He gaseous scintillation detector, which is insensitive to γ-rays, was used to perform a passive fast neutron spectral signature measurement of 233 U 3 O 8 , and was used in conjunction with a pulsed deuterium-tritium neutron generator to demonstrate the differential die-away signature of this material. Furthermore, an array of 3 He detectors was used in conjunction with the same neutron generator to measure the delayed neutron time profile of 233 U, which is unique to this nuclide. These measurements provide a benchmark for future nondestructive assay instrumentation development, and demonstrate a set of key neutron signatures to be leveraged for nuclear safeguards in the thorium fuel cycle.
The astrophysical 29 Si(p,γ) reaction is expected to play a key role in determining the final 29 Si yields ejected in nova explosions. Such yields are used to accurately identify the stellar origins of meteoritic stardust and recently, distinctive silicon isotopic ratios have been extracted from a number of presolar grains. Here, the light-ion 28 Si( 3 He,p) fusion-evaporation reaction was used to populate low-spin proton-unbound excited states in the nucleus 30 P that govern the rate of the astrophysical 29 Si(p,γ) reaction. In particular, γ decays were observed from resonances up to E r = 500keV, and key resonances at 217 and 315 keV have now been identified as 2 + and 2¯ levels, respectively. Here, the present paper provides the first estimate of the 217-keV resonance strength and indicates that the strength of the 315-keV resonance, which dominates the rate of the 29 Si(p,γ) reaction over the entire peak temperature range of oxygen-neon novae, is higher than previously expected. As such, the abundance of 29 Si ejected during nova explosions is likely to be less than that predicted by the most recent theoretical models.
Nonyrast, excited states in neutron-rich 186 W were populated via inelastic-scattering reactions using beams of 136 Xe nuclei accelerated to 725 and 800 MeV. Levels populated in the reactions were investigated via particle-γ coincidence techniques using the Gammasphere array of high-purity germanium detectors and the compact heavy-ion counter, CHICO2. The K π = 2 + (γ ), K π = 0 + and K π = 2 – (octupole) rotational side bands were extended to spins 14h¯, 12h¯, and 13h¯, respectively. A staggering pattern observed in the energies of levels in the K π = 2 + band was found to be consistent with a potential that gets softer to vibration in the γ degree of freedom with increasing spin. Furthermore, the odd-even staggering of states in the K π = 2 – band was found to exhibit a phase opposite to that seen in the γ band; an effect most probably associated with Coriolis coupling to other, unobserved octupole vibrational bands in 186 W.
Approximately half of all atomic nuclei heavier than iron are synthesized by the slow neutroncapture process. The weak component of this process is not well understood and the reaction rates of each isotope in the s-process path affect nucleosynthesis abundances downstream.
Background: Detailed spectroscopy of neutron-rich, heavy, deformed nuclei is of broad interest for nuclear astrophysics and nuclear structure. Nuclei in the r-process path and following freeze-out region impact the resulting r-process abundance distribution, and the structure of nuclei midshell in both proton and neutron number helps to understand the evolution of subshell gaps and large deformation in these nuclei. Purpose: We aim to improve the understanding of the nuclear structure of 160 Gd, specifically the K π = 4 + bands, as well as study the β decay of 160 Eu into 160 Gd. Methods: High-statistics decay spectroscopy of 160 Gd resulting from the β-decay of 160 Eu was collected using the GRIFFIN spectrometer at the TRIUMF-ISAC facility. Results: Two new excited states and ten new transitions were observed in 160 Gd. The β-decaying half-lives of the low- and high-spin isomers in 160 Eu were determined, and the low-spin state's half-life was measured to be t 1/2 = 26.0 (8) s, ≈ 16% shorter than previous measurements. Lifetimes of the two K π = 4 + bandheads in 160 Gd were measured for the first time, as well as γ – γ angular correlations and mixing ratios of intense transitions out of those bandheads. Conclusions: Lifetimes and mixing ratios suggest that the hexadecapole phonon model of the K π = 4 + bandheads in 160 Gd is preferred over a simple two-state strong mixing scenario, although further theoretical calculations are needed to fully understand these states. Additionally, the 1999.0-keV state in 160 Gd heavily populated in β decay is shown to have positive parity, which raises questions regarding the structure of the high-spin β-decaying state in 160 Eu.
The 25 Al(p,γ) 26 Si reaction plays a key role in nucleosynthesis pathways that influence the galactic abundance of 26 Al. A more precise determination of the proton strength of the lowest ℓ = 0 proton resonance in 26 Si is needed to improve reaction-rate calculations. This work measures the 25 Al(d,n) 26 Si proton-transfer reaction in inverse kinematics using a radioactive ion beam at RESOLUT, determining excitation energies and cross sections for the lowest ℓ = 0 resonance associated with the 3 + 3 state at 5.92(2) MeV. Coupled-reaction-channels calculations (FRESCO) are used to extract the ℓ = 0 spectroscopic factor and the corresponding proton width, yielding Γₚ = 2.19(45) eV and a (p,γ) resonance strength of 26(10) meV. This resonance is found to dominate the 25 Al(p,γ) 26 Si reaction rate above 0.2 GK.
Here, in situ synchrotron high-energy x-ray diffraction experiments and detailed transmission electron microscopy (TEM) characterization were conducted on as-fabricated and neutron-irradiated yttrium hydrides. The high-resolution synchrotron x-ray diffraction revealed minor α yttrium and major δ yttrium hydride phases in all specimens. Specimens were subject to heat treatments (heating-cooling cycles), and the intensity of α yttrium partially and completely disappeared in as-fabricated and neutron-irradiated specimens, respectively. The disappearance of α yttrium was unforeseen because hydrogen was expected to leave δ phase, causing an increase in α yttrium diffraction peak intensity. This observation indicated a surplus of hydrogen in the specimens where it was odd for hydride-forming early transition metal elements. The subsequent through-focus TEM characterization discovered nanometric cavities in both as-fabricated and neutron-irradiated yttrium hydride specimens for the first time. Two types of cavities were identified as fabrication-caused and irradiation-induced. The fabrication-caused cavities were associated with regions having linear deformation features, interfaces, and inclusions. The irradiation-induced cavities were observed as being formed isolated in the yttrium hydride phase. The presence of such nanometric cavities was considered as potential hydrogen storage pockets where the overall hydrogen storing capacity of yttrium hydride would be enhanced.
Understanding microstructural and strain evolutions induced by noble gas production in the nuclear fuel matrix or plasma-facing materials is crucial for designing next generation nuclear reactors, as they are responsible for volumetric swelling and catastrophic failure. In this paper we describe a multimodal approach combining synchrotron-based nanoscale X-ray imaging techniques with atomic-scale electron microscopy techniques for mapping chemical composition, morphology and lattice distortion in a single crystal W induced by Kr irradiation. We report that Kr-irradiated single crystal W undergoes surface deformation, forming Kr containing cavities. Furthermore, positive strain fields are observed in Kr-irradiated regions, which lead to compression of underlying W matrix.
There are programs for high-Z shell experiments at the National Ignition Facility (NIF). For shells made of actinide material, a quantitative fission diagnostic is needed in order to determine how much fission took place and whether the fission was sufficient to produce a non-negligible heat source in the burning capsule. Here, we present a viable coupled experimental and theoretical technique for making quantitative fission measurements possible. The proposed scheme involves using a well-characterized set of depleted uranium foils outside an NIF capsule to verify the conversion of xenon and krypton fission fragments collected at the Radiochemical Analysis of Gaseous Samples (RAGS) facility into total fission yield. We present the calculations needed for this conversion, including the decays in and out of fission fragment chains during the RAGS pump-down of the NIF chamber.
The thermonuclear rate of the 29 Si(p,γ) 30 P reaction impacts the 29 Si abundance in classical novae. A reliable reaction rate is essential for testing the nova paternity of presolar stardust grains. At present, the fact that no classical nova grains have been unambiguously identified in primitive meteorites among thousands of grains studied is puzzling, considering that classical novae are expected to be prolific producers of dust grains. We investigated the 29 Si + p reaction at center-of-mass energies of 200–420 keV, and present improved values for resonance energies, level excitation energies, resonance strengths, and branching ratios. One new resonance was found at a center-of-mass energy of 303 keV. For an expected resonance at 215 keV, an experimental upper limit could be determined for the strength. We evaluated the level structure near the proton threshold, and present new reaction rates based on all the available experimental information. Our new reaction rates have much reduced uncertainties compared with previous results at temperatures of T ≥ 140 MK, which are most important for classical nova nucleosynthesis. Furthermore, future experiments to improve the reaction rates at lower temperatures are discussed.
The background index (BI) is an important quantity to project and calculate the half-life sensitivity of neutrinoless double-𝛽 decay (0𝜈𝛽𝛽) experiments. An analysis framework is presented to calculate the BI using the specific activities, masses, and simulated efficiencies of an experiments components as distributions. This Bayesian framework includes a unified approach to combine specific activities from assay. Monte Carlo uncertainty propagation is used to build a BI distribution from the specific activity, mass, and efficiency distributions. This method is applied to the M AJORANA D EMONSTRATOR , which deployed arrays of high-purity Ge detectors enriched in 76 Ge to search for 0𝜈𝛽𝛽. The original assay-based projection is requantified in the new framework, using the as-built geometry of the Demonstrator and additional assay information. While 47% higher than the original projection, the resulting BI of [8.95±0.36]×10 −4 cts/(keVkgyr) from the 232 Th and 238 U decay chains does not account for the higher-than-expected BI observed by the D EMONSTRATOR . Finally, this method enables us to demonstrate the statistical incompatibility between the D EMONSTRATOR 's observed background and the assay results.
The Z = 50 shell closure, near N = 82, is unique in the sense that it is the only shell closure with the spin-orbit partner orbitals, π g 9 / 2 and π g 7 / 2 , enclosing the magic gap. The interaction of the proton hole/particle in the above-mentioned orbitals with neutrons in the νh 11/2 orbital is an important prerequisite to the understanding of the nuclear structure near N = 82 and the νπ interaction. To explore the structural similarity between the high-spin isomeric states in In ( Z = 49), Sn ( Z = 50), and Sb ( Z = 51) isotopes from a microscopic point of view. In addition, to understand the role of a proton hole or particle in the spin-orbit partner orbitals, π g 9 / 2 and π g 7 / 2 , respectively, with neutron holes in the ν h 11 / 2 orbital on these aforementioned isomers. The fusion and transfer induced fission reaction Be 9 ( U 238 , f ) with 6.2 MeV/u beam energy, using a unique setup consisting of AGATA, VAMOS++, and EXOGAM detectors, was used to populate through the fission process and study the neutron-rich In 119 , 121 isotopes. This setup enabled the prompt-delayed γ -ray spectroscopy of isotopes in the time range of 100 ns – 200 μ s . In the odd- A In 119 , 121 isotopes, indications of a short half-life 19/2 - isomeric state, in addition to the previously known 25/2 + isomeric state, were observed from the present data. Further, new prompt transitions above the 25/2 + isomer in In 121 were identified along with reevaluation of its half-life. The experimental data were compared with the theoretical results obtained in the framework of large-scale shell-model calculations in a restricted model space. Furthermore, the ( π g 9 / 2 ν h 11 / 2 ; I | H ^ | π g 9 / 2 ν h 11 / 2 ; I ) two-body matrix elements of residual interaction were modified to explain the excitation energies and the B ( E 2) transition probabilities in the neutron-rich In isotopes. The (i) decreasing trend of E (29/2 + ) - E (25/2 + ) in odd-In (with dominant configuration π g 9 / 2 - 1 ν h 11 / 2 - 2 and maximum aligned spin of 29/2 + ) and (ii) increasing trend of E (27/2 + ) - E (23/2 + ) in odd-Sb (with dominant configuration π g 7 / 2 + 1 ν h 11 / 2 - 2 and maximum aligned spin of 27/2 + ) with increasing neutron number could be understood as a consequence of hole-hole and particle-hole interactions, respectively.