New approach to intranuclear cascades with quantum Monte Carlo configurations
Not Available
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
SEARCH · Engineering Papers
Results for “nuclear configuration interaction”
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.
Cross-shell excited configurations in the structure of 34 Si
The cross-shell excited states of 34 Si have been investigated via β decays of the 4 − ground state and the 1 + isomeric state of 34 Al. Since the valence protons and valence neutrons occupy different major shells in the ground state as well as the intruder 1 + isomeric state of 34 Al, intruder levels of 34 Si are populated via allowed β decays. Spin assignments to such intruder levels of 34 Si were established through γ-γ angular correlation analysis for the negative-parity states. The configurations of such intruder states play crucial roles in our understanding of the N = 20 shell gap evolution. A configuration interaction model derived from the FSU Hamiltonian was utilized in order to interpret the intruder states in 34 Si. Furthermore, shell model interaction derived from a more fundamental theory with the valence space in medium similarity renormalization group method was also employed to interpret the structure of 34 Si.
20 ≤ A ≤ 38↗
The DUNE-DAQ Application Framework
The deep underground neutrino experiment (DUNE) is a next-generation neutrino experiment that will probe the properties of these elusive particles with unparalleled precision. It will also act as an observatory for neutrino bursts caused by nearby supernovae, in the event that one occurs, while the experiment is in operation. Given these goals, the DUNE trigger and DAQ system must be able to maintain extremely high uptime and provide a path for full readout of the detectors for very long times (up to 100 s). To achieve these ends, we have designed the DUNE DAQ system around a flexible “application framework,” which provides a modular interface for specific tasks while handling the interconnections between them. The application framework collects modules into applications, which can then be interacted with as units by the control, configuration, and monitoring systems. One of the key features of the framework is its communication abstraction layer, which allows for modules to interact with both internal queues and external network connections with a single transport-agnostic interface. We will report on the architecture and features of the framework.
46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
A modular apparatus for use in high-precision measurements of parity violation in polarized eV neutron transmission
We describe a modular apparatus for use in parity-violation measurements in epithermal neutron–nucleus resonances with high instantaneous neutron fluxes at the Manuel Lujan Jr. Neutron Scattering Center at Los Alamos National Laboratory. This apparatus is designed to conduct high-precision measurements of the parity-odd transmission asymmetry of longitudinally polarized neutrons through targets containing nuclei with p-wave neutron–nucleus resonances in the 0.1–10 eV energy regime and to accommodate a future search for time reversal violation in polarized neutron transmission through polarized nuclear targets. The apparatus consists of an adjustable neutron and gamma collimation system, a 3 He- 4 He ion chamber neutron flux monitor, two identical cryostats for target cooling, an adiabatic eV-neutron spin flipper, a near-unit efficiency 6 Li- 7 Li scintillation detector operated in current mode, a flexible CAEN data acquisition system, and a neutron spin filter based on spin-exchange optical pumping of 3 He gas. We describe the features of the apparatus design devoted to the suppression of systematic errors in parity-odd asymmetry measurements. We describe the configuration of the apparatus used to conduct a precision measurement of parity violation at the 0.7 eV p-wave resonance in 139 La which employs two identical 139 La targets, one to polarize the beam on the p-wave resonance using the weak interaction and one to analyze the polarization.
46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
Engineered Barrier System R&D and International Collaborations – LANL (FY20): Spent Fuel and Waste Disposition
The United States Department of Energy’s Spent Fuel and Waste Disposition (SFWD) program is investigating the design and safety function of generic nuclear geologic repositories in a variety of geologic settings (salt, argillite, and crystalline rock). Different configurations and loadings of spent nuclear fuel and waste within disposal canisters are also being investigated, some of which have the potential to generate repository temperatures higher than previously considered (i.e., temperatures >100ºC) by foreign and domestic concepts. This report expands on engineered barrier material stability in a high temperature crystalline rock repository through high temperature hydrothermal experiments. Experiments were designed to develop engineered barrier system (EBS) concepts in a hightemperature crystalline environment in 1) bentonite-Grimsel Granodiorite interactions, 2) bentonite-cement reactions, and 3) interaction between waste canister materials and bentonite. Experiment results are applied to understanding long-term repository performance in terms of radionuclide isolation. One hydrothermal experiment was completed in the rocking autoclaves at LANL in FY-20: IEBS-6 (Grimsel Granodiorite + Wyoming bentonite + cured ordinary Portland cement + Grimsel Granodiorite synthetic groundwater, 250ºC/150 bar, 8 weeks). Several other experiments were planned but were delayed due to the COVID-19 pause in laboratory work. In addition, some characterization of the reaction products of experiments conducted in FY-20 was prevented by COVID-19. Mineral phase chemistry was not measured via electron microprobe analyses for IEBS-6. Quantitative X-ray diffraction results from IEBS-6 and HBT-1 were not completed. The missing analyses and discussion of the results will be included in next year’s report. New characterization that was completed includes scanning electron microscopy of reaction products from IEBS-6, quantitative X-ray diffraction results are presented from IEBS-1 through IBES-5, preliminary scanning electron microscope images and chemical analyses for IEBS-6, X-ray diffraction of the clay fraction from all experiments, measurement/imaging of mineral growth on the surface of steel coupons. Major observations pertaining to bentonite stability in a Grimsel Granodiorite environment include the stability of Na-montmorillonite at 250°C, the formation of trace CSH phases, and the formation of bentonite colloids on experiment cooling. The addition of a cured chip of Portland cement to the bentonite-Grimsel system results in slightly higher pH values and the formation of diverse secondary mineral phases that were not observed in the previous experiments (e.g., analcime, garronite, CSH phases). The new characterization efforts related to the interaction of stainless-steel coupons and bentonite clay focused on thickness and mineralogy of phases that formed at the steel surface. In Wyoming bentonite + Grimsel Granodiorite systems, newly formed minerals at the bentonitesteel coupon interface included alteration of the outermost steel edge to Fe,Cr-oxide phases, followed by Fe-rich phyllosilicates (Fe-saponite, chlorite) and interbedded Fe,Cr,Ni-sulfide phases (pentlandite). Hydrothermal experiments were completed to assess uranium-carbonate complexation at conditions relevant to high-temperature disposal. Autoclave solubility experiments were conducted at 150 to 250ºC with a range of carbonate and uranium concentrations. The experiment results were characterized via situ UV-Visible spectroscopy and synchrotron-based in situ XAS techniques. Results show a significant decrease in the stability of uranyl-carbonate complexes at temperatures above 100ºC along with the precipitation of uranium oxides. Further, at T > 200ºC, results show that uranyl-hydroxyl complexes control solubility of uranium instead of the previously predicted uranyl-carbonate species. These results are significant for understanding the mobility of uranium in the EBS, which will likely contain carbonate-rich fluids. International research efforts focused on three main areas: 1) participation in international conferences, 2) building collaborations with foreign repository programs, and 3) the initiation of an experimental program to complement the full-scale HotBENT test at the Grimsel test site. This experiment included Wyoming bentonite + low carbon steel + Grimsel Granodiorite synthetic groundwater and was run at the planned maximum temperature of the HotBENT test (200°C). Complete characterization of reaction products was hindered by disruptions to laboratory work but will be reported in the next FY. The experimental results obtained in FY-20 continue to document the wide-ranging effects of bulk composition and pressure-temperature conditions in the mineralogical and geochemical evolution of a high-temperature repository environment. Concepts developed will be used to inform models of long-term material stability in a generic crystalline rock-hosted repository.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Two Delayed Critical 15-Inch-Diameter Interacting Enriched (93.14) Uranium Metal Cylinders without Moderator and Reflector
This report documents very accurately the configuration and the materials for experiments with two unmoderated, unreflected, interacting, coaxial, highly enriched, 15-in.-diameter, uranium metal cylinders performed at the Oak Ridge Critical Experiments Facility (ORCEF) in May to Aug 1963 and described in logbook E-19 and E -20 associated with experiments in the East cell of ORCEF. Measurements were also performed in April and May of 1965 and described in logbooks E-22 and E-23 The information is sufficiently accurate that it can be used as the basis for preparation of benchmarks for International Criticality Safety Benchmark Program (ICSBEP) at Idaho National laboratory. The thickness of the cylinders was varied and the spacing between them was adjusted to achieve a delayed critical configuration. The average enrichment of the uranium metal was 94.14 wt. % 235 U. The heights of the 15-in.-diameter, equal-height cylinders varied from 1-5/8 to 3.0-inches. All interacting cylinders were assembled coaxially with their flat faces parallel and their combined masses varied between 182 and 325 kilograms of HEU metal. The data from these 12 experiments described would be acceptable for use as criticality safety benchmark experiments for the ICSBEP and EURATOM’s Nuclear Energy Agency nuclear criticality safety benchmark program, once the uncertainty analysis is completed. Based on previous ICSBEP benchmarks with this enriched uranium metal at ORCEF, it is expected that the uncertainties in k eff could be as low as ± 0.0002.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Three Delayed Critical 15-Inch-Diameter Interacting Enriched (93.14) Uranium Metal Cylinders Without Moderator and Reflector
This report documents very accurately the configuration and the materials for experiments with three unmoderated, unreflected, interacting, coaxial, highly enriched, 15-in.-diameter, uranium metal cylinders performed at the Oak Ridge Critical Experiments Facility (ORCEF) in August and September 1963 and described in logbook E 20 associated with experiments in the East cell of ORCEF. The information is sufficiently accurate that it can be used as the basis for preparation of benchmarks for International Criticality Safety Benchmark Program (ICSBEP) at Idaho National laboratory. The thickness of the cylinders was varied and the spacing between them was adjusted to achieve a delayed critical configuration. The average enrichment of the uranium metal was 94.14 wt. % 235U. The heights of the 15-in.-diameter, equal-height cylinders varied from 1-1/8to 2.00 inches. All interacting cylinders were assembled coaxially with their flat faces parallel and their combined masses varied between 182 and 325 kilograms of HEU metal. The data from these six experiments described are judged to be acceptable for use as criticality safety benchmark experiments for the ICSBEP and EURATOM’s Nuclear Energy Agency nuclear criticality safety benchmark program, once the uncertainty analysis is completed. Based on previous ICSBEP benchmarks with this enriched uranium metal at ORCEF, it is expected that the uncertainties in measured keff could be as low as ± 0.0002.
36 MATERIALS SCIENCE↗
Generation of a controllable TNSA deuteron beam using deuterated metal targets
This study investigated the feasibility of using deuterated titanium targets for the generation of a target normal sheath acceleration (TNSA) deuteron beam. Commercial 25-μm-thick titanium foil was cut into 500 x 500-μm 2 squares and subsequently deuterated using different approaches. The spectra and total yields of all emitted ions were measured using a Thomson parabola ion spectrometer. It was found that a 24-h exposure to 1 atm of D 2 gas at 400°C is the most-efficient deuteration method, producing yields in the mid-10 11 deuterons per shot. The deuteron energy spectra changed from an exponential shape at low ion yields to exponentially modified Gaussians at higher yields. Simulations suggest that this effect is due to increased Coulomb interactions between the ions, which suppress the low energies. Separate campaigns will utilize the presented approach to produce a TNSA triton beam. Furthermore, two such targets will be used in a pitcher/catcher configuration to study the tritium–tritium reaction T(t, 2n)α.
46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
Ab initio many-fermion structure calculations on a quantum computer
To overcome the limitations of existing algorithms for solving self-bound quantum many-body problems—such as those encountered in nuclear and particle physics—that access only a restricted subset of energy levels and provide limited structural information, we introduce and demonstrate a novel quantum-classical approach capable of resolving the complete bound-state spectrum. This method also provides the total angular momentum 𝐽 associated with each eigenstate. Here, our approach is based on expressing the Hamiltonian in second-quantized form within a novel input model combined with a scan scheme, enabling broad applicability to configuration-interaction calculations across diverse fields. We apply this hybrid method to compute, for the first time, the bound-state spectrum together with corresponding 𝐽 values of 20 O using a realistic strong-interaction Hamiltonian. Our approach applies to hadron spectra and 𝐽 values solved in the relativistic basis light-front quantization approach.
Du, Weijie [Chinese Academy of Sciences (CAS), Lan↗
Coupled-cluster theory for strong entanglement in nuclei
Atomic nuclei can exhibit shape coexistence and multireference physics that enters in their ground states, and accurately capturing the ensuing correlations and entanglement is challenging. Here, we address this problem by applying single-reference coupled-cluster theory based on spherical and deformed reference states and the tailored coupled-cluster method. The latter combines configuration interactions to capture static correlations with coupled-cluster theory for dynamic correlations. We compute the atomic nuclei 12 C, 28 Si, and 56 Ni and find that the tailored coupled-cluster method and the single-reference approach based on a deformed Hartree-Fock state yield the most accurate results.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Investigating the crust of neutron stars with neural-network quantum states
An accurate description of low-density nuclear matter is crucial for explaining the physics of neutron star crusts. In the density range between approximately 0.01 fm −3 and 0.1 fm −3 , matter transitions from neutron-rich nuclei to various higher-density pasta shapes, before ultimately reaching a uniform liquid. In this work, we introduce a variational Monte Carlo method based on a neural Pfaffian-Jastrow quantum state, which allows us to model the transition from the liquid phase to neutron-rich nuclei microscopically. At low densities, nuclear clusters dynamically emerge from the microscopic interactions among protons and neutrons, which we model based on pionless effective field theory. Our variational Monte Carlo approach represents a significant improvement over the state-of-the-art auxiliary-field diffusion Monte Carlo method, which is severely hindered by the fermion-sign problem in this low-density regime and cannot capture the onset of clusters. In addition to computing the energy per particle of symmetric nuclear matter and pure neutron matter, we analyze an intermediate isospin-asymmetry configuration to elucidate the formation of nuclear clusters. We also provide evidence that the presence of such nuclear clusters influences the amount of protons in the crust compared to protons in beta-equilibrated, neutrino-transparent matter.
Nuclear astrophysics↗
A Multireference Approach to Electron and Electron–Nuclear Dynamics in Nanomaterials (Final Report)
Many important chemical and physical phenomena involve dynamics on large number of electronic states. Thus, there is a critical need to develop methods to simulate dynamics in dense manifolds of states. Towards this end, we have: a) developed the multiple cloning in dense manifolds of states (MCDMS) method, which is capable of accurately modeling the quantum mechanical coherence between populations on a large number of electronic states, b) implemented MCDMS into the free, open-source PySpawn software package, c) developed graphics processing unit-accelerated algorithms modeling electron dynamics in light fields via Floquet time-dependent configuration interaction (F-TDCI), and d) critically compared different orbital bases in order to achieve an accurate and efficient F-TDCI expansion. This grant ended in August 2020, when our group moved to from Michigan State University to Stony Brook University, where this project continues under grant number DE-SC0021643.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Radiation Dose Modeling for Niowave’s Accelerator Driven Uranium Target Assembly 3
Molybdenum-99 is a high-value radionuclide commonly used for medical purposes within the United States. The National Nuclear Security Administration (NNSA) seeks to reliably produce the radioisotope 99 Mo without the use of highly enriched uranium. NNSA’s Office of Material Management and Minimization (M3) provides funding and government laboratory expertise to private companies to expedite the production process domestically and currently funds designs that use low-enriched uranium or other 99 Mo production pathways. Several production designs are being explored across the industry, including uranium fission and photonuclear conversion of 100 Mo targets. Niowave Inc. seeks to produce 99 Mo via a high-energy electron accelerator that strikes a lead-bismuth eutectic target that ultimately produces a consistent neutron flux. The neutron flux then interacts in a subcritical reactor core configuration to produce fission in low-enriched or natural uranium targets. These fissionable targets are then processed to extract 99 Mo. The purpose of this work is to estimate the neutron and photon dose response across Niowave’s proposed facility for worker safety during operation. Owing to the size of the proposed Niowave facility and necessary shielding, unbiased Monte Carlo radiation transport is impractical, and variance reduction methods are required. This work focuses on the weight window variance reduction method to produce high confidence dose response results within a Monte Carlo radiation transport code. Specifically, an adjoint-informed weight window methodology was created to improve the dose response estimates for accelerator-driven subcritical reactor designs. This adjoint-informed methodology was implemented for Niowave’s proposed design and improved dose results at far-field locations across the facility. Acceptable dose rate contours for the proposed facility were generated across the facility and are presented in this work.
07 ISOTOPE AND RADIATION SOURCES↗
Understanding baryon stopping at the BNL Relativistic Heavy Ion Collider top energies
The nucleon exhibits a rich internal structure governed by quantum chromodynamics (QCD), where its electric charge arises from valence quarks, while its spin and mass emerge from complex interactions among valence quarks, sea (anti)quarks, and gluons. At the advent of QCD, an alternative hypothesis emerged suggesting, at high energies, the transport of a nucleon's baryon number could be traced by a nonperturbative configuration of gluon fields connecting its three valence quarks, forming a 𝑌-shaped topology known as the gluon junction. Recent measurements by the STAR experiment are compatible with this scenario. In light of these measurements, this study aims to explore the mechanisms of baryon transport in high-energy nuclear collisions using the pythia-8 framework, which incorporates a state-of-the-art hadronization model with advanced color flow (CF) and color reconnection (CR) mechanisms that mimic signatures of a baryon junction. Within this model setup, we investigate (i) the rapidity slope of the net-baryon distributions in photon-included processes (𝛾 + 𝑝) and (ii) baryon over charge transport in the isobaric (Ru + Ru and Zr + Zr) collisions. Our study highlights the importance of the CF and CR mechanisms in pythia-8, which play a crucial role in baryon transport. The results show that the CF and CR schemes significantly affect the isobaric baryon-to-charge ratio, leading to different predictions for baryon stopping and underscoring the need to account for CF and CR effects in comparisons with experimental measurements.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Radiation-hard miniature optical engine with high bandwidth
Future silicon trackers will be operated in an intense radiation environment and require large volumes of data to be transmitted off detector. In addition, the optical modules must be of low mass in order to limit multiple scattering and nuclear interactions that would degrade the overall performance of the detector. Here, we present a miniature optical engine that satisfies these constraints. The optical engine consists of an ASIC driving a VCSEL (Vertical Cavity Surface Emitting Laser) array in an optical package. Two ASICs are designed to operate a 12-channel VCSEL array at 1.28 or 5.12 Gb/s per channel, which yields a total data rate of up to 60 Gb/s. The core transistors are fabricated in a 65 nm CMOS process which enhance the radiation-hardness. Each channel contains equalizer (CTLE) and clock-data recovery circuits (CDR) so that the ASIC can restore the highly distorted electrical signal after propagating through several meters of cables of small diameter. The equalizer, CDR, and VCSEL driver are configured via a digital I 2 C chip interface with triple redundant memory to mitigate single event upset (SEU) effects. The bias and modulation currents are controlled by a digital-to-analog converter (DAC). We present the design of the circuit together with the results of the simulations and preliminary measurements.
47 OTHER INSTRUMENTATION↗
A New 1D Model for Thermal Mixing and Stratification in Advanced Reactor Transients
Thermal mixing and stratification in large pools and enclosures play a critical role in the safety and performance of pool-type nuclear reactors, particularly during transient scenarios involving significant temperature differences between incoming and bulk coolant. Accurate modeling of these phenomena is essential for predicting system behavior and supporting passive safety features such as natural circulation. Here, this paper presents a new 1D model for thermal mixing and stratification, developed and implemented in the SAM code. The model represents a large pool as 1D coolant jet channels and zero-dimensional bulk pool volumes, enabling the simulation of a wide range of flow configurations, including hot and cold jet interactions, stratified layers, and the influence of complex geometries such as ceilings, free surfaces, and internal obstacles. Heat exchange between jet and pool regions is governed by closure relations calibrated against 3D computational fluid dynamics (CFD) simulations. The model improves upon earlier approaches by incorporating time-dependent jet characteristics and capturing the associated delay effects more accurately. Code-to-code comparisons and validation against experimental data from the Thermal Stratification Test Facility demonstrate the model’s accuracy and flexibility. This work offers two key contributions: (1) an efficient and robust method for simulating thermal mixing and stratification at the system level, eliminating the need for external coupling between system analysis codes and CFD, and (2) a significant enhancement of SAM’s capabilities to analyze thermal stratification phenomena in advanced reactor systems.
SAM↗
Ground and excited state gradients with end-to-end differentiable semiempirical quantum chemistry
Accurate and efficient gradients of molecular energy with respect to nuclear degrees of freedom are essential for geometry optimization and molecular dynamics, including simulations that go beyond the Born–Oppenheimer regime. A common approach involves deriving analytical formulas for new electronic structure methods, which is often conceptually difficult and requires tedious coding. Here, we implement analytical, semi-numerical, and automatic differentiation (AD)-based gradient pathways for semiempirical Hamiltonian models in the PYSEQM software package, leveraging both graphics processing unit (GPU) and central processing unit (CPU) architectures. We further extend these capabilities to excited states calculated using the configuration interaction singles and time-dependent Hartree–Fock ansätze. We benchmark wall time, peak memory usage, and accuracy across three molecular families of varying chemical complexity, including systems of up to a thousand atoms. For ground-state simulations, analytical and AD gradients achieve near-identical GPU runtimes, while semi-numerical gradients are slower on GPU but remain competitive on CPU. For excited states, both analytical and custom AD approaches using implicit differentiation show similar performance and low memory requirements, whereas gradients with full AD are memory-limited. AD gradients match analytical ones in accuracy across all tested systems, aided by a quaternion-based diatomic frame rotation for two-center quantities that ensures smooth energy surfaces. Overall, automatic differentiation emerges as a practical alternative to analytical gradients in semiempirical quantum chemistry, offering high accuracy while allowing seamless integration in AI-driven workflows and popular packages, such as PyTorch and JAX. Our results provide actionable guidance for selecting optimal gradient strategies in large-scale ground- and excited-state molecular dynamics simulations.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Relativistic configuration-interaction and many-body-perturbation-theory calculations of U i hyperfine constants
Neutral uranium (U i) is a very difficult atom for theoretical calculations due to a large number (six) of valence electrons, strong valence-valence and valence-core correlations, high density of states, and relativistic effects. Configuration-interaction many-body perturbation theory (CI-MBPT) can efficiently treat valence-core correlations and relativistic effects, but because the formalism was developed for a Dirac-Hartree-Fock (DHF) starting potential that does not contain valence electrons, quite large CI space is needed to compensate for the + 6 charge of such a potential. Much more efficient is the relativistic configuration-interaction (RCI) approach, which uses a relatively accurate starting DHF potential that includes some valence electrons to make the valence-electron Hamiltonian diagonally dominated for some states. Here we report calculations of U i hyperfine constants of several low-energy states using the RCI method with the starting potential that includes four f valence electrons. With this starting potential, it is possible to use the single-configuration approximation or small basis sets to obtain quite accurate results for hyperfine-structure constants. In fact, by scaling the nuclear magnetic moment, the agreement for five levels was within 5% and a new magnetic moment can be recommended, 0.43(2). When two states are mixed, it is difficult to predict the exact mixing fractions, so in this case we proposed a method in which the mixing fraction is found by optimization of one parameter in the MBPT correction to obtain correct g factors. The resulting hyperfine constants are significantly improved. Finally, the CI-MBPT approach was also tested, with some success for low-energy states, with the limited basis set. The methods investigated here can be further developed to include more extensive data sets to improve accuracy and can be applied to other atoms and for calculations of other properties, for example, relevant to fundamental symmetry tests.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗