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.
The proton induced NIELs for representative spacecraft materials are presented for the energy range between the displacement thresholds of the material to 1000 MeV. All interaction mechanisms (Coulomb and nuclear elastic/inelastic) are fully accounted in the present NIEL calculations.
To meet the challenge of future deep space programs, an accurate and efficient engineering code for analyzing the shielding requirements against high-energy galactic heavy radiation is needed. To address this need, a new Green's function code capable of simulating high charge and energy ions with either laboratory or space boundary conditions is currently under development. The computational model consists of combinations of physical perturbation expansions based on the scales of atomic interaction, multiple scattering, and nuclear reactive processes with use of the Neumann-asymptotic expansions with non-perturbative corrections. The code contains energy loss due to straggling, nuclear attenuation, nuclear fragmentation with energy dispersion and downshifts. Previous reports show that the new code accurately models the transport of ion beams through a single slab of material. Current research efforts are focused on enabling the code to handle multiple layers of material and the present paper reports on progress made towards that end.
We report on our study of the design of a neutrino detector, shielding and veto array needed to operate a neutrino detector in space close to the Sun. This study also took into account the expected rates of Galactic gamma and cosmic rays in addition to the particles from the Sun.These preliminary studies show that we can devise a detector such that a small signal of neutrino interactions can be extracted from a large random number of events from the background sources using a double timing method from the conversion electron produced in the neutrino interaction and a secondary delayed signal from the nuclear excited state produced from the initial neutrino interaction; in our case the conversion of Ga 69 or 71 into Ge 69 or 71, but this method could apply to other nuclei with large neutrino cross sections such as Ir 115. Although these types of events need to be above 0.405 megaelectronvolt (MeV) neutrino energy and are only 66 percent of all conversion neutrino interactions on Gallium, this is a small price to pay for an increase of 10,000 by going close to the Sun to enhance the neutrino rate over the background combatorical fake-signal events. The conclusion of this Phase-1 study is very positive in that we can get the backgrounds less than 20 percent fake signals, and in addition to this we have devised another shielding method that makes the Galactic gamma-ray rate a hundred fold less which will make further improvements over these initial estimates. Although these studies are very encouraging it suggests that the next step is a NIAC Phase-II to actually build a test device,measuring basic principles such as light attention within the scintillator with high dopants and to take data in the lab with a cosmic-ray test stand and triggered X-ray source for comparison with simulated expected performance of the detector. This would be the perfect lead into a future proposal beyond a NIAC (NASA Innovative Advanced Concepts) Phase-II for a test flight of a small one-pint detector in orbit of the detector concept beyond Earth outside of the radiation belts.
This paper summarizes my work during the Summer 2025 SULI internship, which focused on two main projects and broader scientific development. The first project involved improving the particle identification (PID) of protons, deuterons, and tritons using PIDA distributions and template fitting, with the goal of modeling nuclear final-state interactions (FSI) and testing the robustness of the method against systematic uncertainties. These techniques pave the way for future application to LArTPC data from the ICARUS detector. The second project centered on the optimization and data-taking of the TinyTPC detector, a compact LArTPC used for high-resolution low-energy measurements. I adjusted gain and threshold parameters, performed hardware validation tests, and developed analysis strategies to extract meaningful physics from collected data. Throughout the summer, I also enhanced my scientific communication and mentorship skills through presentations, collaborative analysis, and peer guidance.
The corrosion and oxidation of actinide metals, leading to the formation of metal-oxide surface layers with the catalytic evolution of hydrogen, impacts the management of nuclear materials. In this work, the interaction of hydrogen with actinide dioxide (AnO 2 , An = U, Np, or Pu) (011) surfaces by Hubbard corrected density functional theory (PBEsol+U) has been studied, including spin–orbit interactions and non-collinear 3k anti-ferromagnetic behavior. The actinide dioxides crystalize in the fluorite-type structure, and although the (111) surface dominates the crystal morphology, the (011) surface energetics may lead to more significant interaction with hydrogen. The dissociative adsorption of hydrogen on the UO 2 (0.44 eV), NpO 2 (-0.47 eV), and PuO 2 (-1.71 eV) (011) surfaces has been calculated. It is found that hydrogen dissociates on the PuO 2 (011) surface; however, UO 2 (011) and NpO 2 (011) surfaces are relatively inert. Recombination of hydrogen ions is likely to occur on the UO 2 (011) and NpO 2 (011) surfaces, whereas hydroxide formation is shown to occur on the PuO 2 (011) surface, which distorts the surface structure.
The Li ge Intranuclear Cascade (INCL) model is a nuclear-physics model that simulates hadron (baryon, anti-baryon and meson) reactions on nuclei, for incident energies ranging from a few tens of MeV to 10-20 GeV. The INCL model has been well validated by hadron scattering data. In my work, I implement an interface in GENIE to use the INCL nuclear model in the simulations of both the initial state of the target nucleus and the Final State Interaction in neutrino-nucleus interaction. It has a consistent treatment of nuclear models in both neutrino interaction and hadron rescattering. A full event record including neutrino vertex and each vertex of hadron rescattering has been accomplished. Several processes, e.g. cluster production, Delta transportation and de-excitation will also be included as benefits of the implementation of the INCL model in GENIE. I will show some initial simulation results showcasing the new GENIE features and discuss plans for making them available for use in experimental analyses.
The Liège Intranuclear Cascade (INCL) model is a nuclear-physics model that simulates hadron (baryon, anti-baryon and meson) reactions on nuclei, for incident energies ranging from a few tens of MeV to 10-20 GeV. The INCL model has been well validated by hadron scattering data. In my work, I implement an interface in GENIE to use the INCL nuclear model in the simulations of both the initial state of the target nucleus and the Final State Interaction in neutrino-nucleus interaction. It has a consistent treatment of nuclear models in both neutrino interaction and hadron rescattering. A full event record including neutrino vertex and each vertex of hadron rescattering has been accomplished. Several processes, e.g. cluster production, Delta transportation and de-excitation will also be included as benefits of the implementation of the INCL model in GENIE. I will show some initial simulation results showcasing the new GENIE features and discuss plans for making them available for use in experimental analyses.
Ab initio calculations have an essential role in our fundamental understanding of quantum many-body systems across many subfields, from strongly correlated fermions to quantum chemistry and from atomic and molecular systems to nuclear physics. One of the primary challenges is to perform accurate calculations for systems where the interactions may be complicated and difficult for the chosen computational method to handle. Here we address the problem by introducing an approach called wavefunction matching. Wavefunction matching transforms the interaction between particles so that the wavefunctions up to some finite range match that of an easily computable interaction. This allows for calculations of systems that would otherwise be impossible owing to problems such as Monte Carlo sign cancellations. We apply the method to lattice Monte Carlo simulations of light nuclei, medium-mass nuclei, neutron matter and nuclear matter. We use high-fidelity chiral effective field theory interactions and find good agreement with empirical data. These results are accompanied by insights on the nuclear interactions that may help to resolve long-standing challenges in accurately reproducing nuclear binding energies, charge radii and nuclear-matter saturation in ab initio calculations.
Two-neutrino double-β decay of certain nuclear isotopes is one of the rarest Standard Model processes observed in nature. Its neutrinoless counterpart is an exotic lepton-number nonconserving process that is widely searched to determine if the neutrinos are Majorana fermions. In order to connect the rate of these processes to the Standard Model and beyond the Standard Model interactions, it is essential that the corresponding nuclear matrix elements are constrained reliably from theory. Lattice quantum chromodynamics (LQCD) and low-energy effective field theories (EFTs) are expected to play an essential role in constraining the matrix element of the two-nucleon subprocess, which could in turn provide the input into ab initio nuclear-structure calculations in larger isotopes. Focusing on the two-neutrino process $nn\rightarrow pp(ee\bar{v}_{e}\bar{v}_e)$, the amplitude is constructed in this work in pionless EFT at next-to-leading order, demonstrating the emergence of a renormalization-scale independent amplitude and the absence of any new low-energy constant at this order beyond those present in the single-weak process. Most importantly, it is shown how a LQCD four-point correlation function in Euclidean and finite-volume spacetime can be used to constrain the Minkowski infinite-volume amplitude in the EFT. The same formalism is provided for the related single-weak process, which is an input to the double-β decay formalism. The LQCD-EFT matching procedure outlined for the double-weak amplitude paves the road toward constraining the two-nucleon matrix element entering the neutrinoless double-β decay amplitude with a light Majorana neutrino.
Nuclear charge radii are sensitive probes of different aspects of the nucleon–nucleon interaction and the bulk properties of nuclear matter, providing a stringent test and challenge for nuclear theory. Experimental evidence suggested a new magic neutron number at N = 32 in the calcium region, whereas the unexpectedly large increases in the charge radii open new questions about the evolution of nuclear size in neutron-rich systems. By combining the collinear resonance ionization spectroscopy method with β-decay detection, we were able to extend charge radii measurements of potassium isotopes beyond N = 32. Here we provide a charge radius measurement of 52 K. It does not show a signature of magic behaviour at N = 32 in potassium. The results are interpreted with two state-of-the-art nuclear theories. The coupled cluster theory reproduces the odd–even variations in charge radii but not the notable increase beyond N = 28. This rise is well captured by Fayans nuclear density functional theory, which, however, overestimates the odd–even staggering effect in charge radii. These findings highlight our limited understanding of the nuclear size of neutron-rich systems, and expose problems that are present in some of the best current models of nuclear theory.
Monoenergetic muon neutrinos at 236 MeV are readily produced in intense medium-energy proton facilities (≳2–3GeV) when a positive kaon decays at rest (KDAR; K + →μ + ν μ ). These neutrinos provide a unique opportunity to both study the neutrino interaction and probe the nucleus with a monoenergetic weak-interaction-only tool. We present cross-section calculations for quasielastic scattering of these 236-MeV neutrinos off 12 C and 40 Ar, paying special attention to low-energy aspects of the scattering process. Our model takes the description of the nucleus in a mean-field approach as the starting point, where we solve Hartree-Fock equations using a Skyrme type nucleon-nucleon interaction. Thereby, we introduce long-range nuclear correlations by means of a continuum random-phase approximation (CRPA) framework where we solve the CRPA equations using a Green's function method. The model successfully describes (e,e') data on 12 C and 40 Ca in the kinematic region that overlaps with the KDAR ν μ phase space. In addition to these results, we present future prospects for precision KDAR cross-section measurements and applications of our calculations in current and future experiments that will utilize these neutrinos.
Not Available
There is a renewed attention to nuclear fusion as a commercial source of carbon-free energy, however there are many scientific needs that must be addressed to enable the future success of fusion as an economical energy option. Among these is the proper description of the impact of radiation produced in the fusion vessel chamber and all other components of the reactor. In this work we will focus on the nuclear data needs to describe the interaction between primary and secondary neutron radiation and the medium- and high-mass nuclei commonly present in structural (such as stainless steel) and superconducting (e.g., electromagnets) materials.
The U.S. Nuclear Regulatory Commission (NRC) has interacted with vendors pursuing the commercialization of micro-reactors (i.e., reactors capable of producing about 1 MW(th) to 20 MW(th) of energy from nuclear fission). It is envisioned that micro-reactors could be assembled and fueled in a factory and shipped to a site. Many of the sites are expected to be remote locations requiring off-grid power or in some cases military bases. The objective of this effort is to explore the technical issues and the approach required to reach a finding of "reasonable assurance of public health and safety" for this new and different class of reactors. The analysis performed here leverages available micro-reactor design and testing data available from national laboratory experience as well as commercial design information to explore technical issues. Some factors considered include source term, accidents that would need to be analyzed, and the extent of the probabilistic risk assessment (PRA). The technical evaluation was performed within the framework of the Licensing Modernization Project (LMP) to identify licensing basis events, classification of structures, systems and components, and defense-in-depth needed to provide regulatory certainty. With this framework and technical evaluation in mind, the scope and content of a micro-reactor licensing application is discussed.
The properties of excitons, or correlated electron–hole pairs, are of paramount importance to optoelectronic applications of materials. A central component of exciton physics is the electron–hole interaction, which is commonly treated as screened solely by electrons within a material. However, nuclear motion can screen this Coulomb interaction as well, with several recent studies developing model approaches for approximating the phonon screening of excitonic properties. While these model approaches tend to improve agreement with experiment, they rely on several approximations that restrict their applicability to a wide range of materials, and thus far they have neglected the effect of finite temperatures. Here, we develop a fully first-principles, parameter-free approach to compute the temperature-dependent effects of phonon screening within the ab initio GW -Bethe–Salpeter equation framework. We recover previously proposed models of phonon screening as well-defined limits of our general framework, and discuss their validity by comparing them against our first-principles results. We develop an efficient computational workflow and apply it to a diverse set of semiconductors, specifically AlN, CdS, GaN, MgO, and SrTiO 3 . We demonstrate under different physical scenarios how excitons may be screened by multiple polar optical or acoustic phonons, how their binding energies can exhibit strong temperature dependence, and the ultrafast timescales on which they dissociate into free electron–hole pairs.
The complex interactions of emerging and disruptive technologies (EDTs) could significantly impact nuclear decision-making, particular in an escalating regional conventional conflict. Such conflicts may present governments with a range of nuclear decisions: whether to introduce a nuclear dimension to a crisis, whether to cross the nuclear threshold through limited nuclear use, how to respond to a limited nuclear attack, whether to expand the scope and intensity of initial limited attacks, and whether to escalate to an all-out nuclear war. At each decision point, EDTs create potential risks as well as rewards. EDTs are likely to influence the context for nuclear decision-making and the choices between different courses of action. EDTs could impact the context of nuclear decision-making by improving or degrading situational assessment, the ability to deliberate, and the ability to manage one’s nuclear forces. EDTs could influence the choice between nuclear restraint or escalation by affecting the perceived strategic benefits, escalatory risks, and operational requirements associated with different courses of action. Even though particular combinations of EDTs could precipitate nuclear use in some scenarios, they could encourage restraint in others. The impact and relevance of the same combinations of EDTs might be different at various nuclear decision points. The availability of specific combinations of EDTs at different stages of a conflict would also vary because of the attrition and one-time-use nature of some capabilities. In later stages of a conflict, the decision maker’s confidence in different combinations of EDTs would depend on their previous experience in using them. While the interactions of EDTs are likely to bring additional complexity to a nuclear decision-making process, EDTs are also not the only source of complexity. Broader strategic, military, operational, legal, moral, and emotional factors are also likely to play an important role. These factors may dominate decision-making in a range of potential cases.
Advanced control rooms (ACR's) will utilize human-system interface (HSI) technologies that may have significant implications for plant safety in that they will affect the operator's overall role and means of interacting with the system. The Nuclear Regulatory Commission (NRC) reviews the human factors engineering (HFE) aspects of HSI's to ensure that they are designed to good HFE principles and support performance and reliability in order to protect public health and safety. However, the only available NRC guidance was developed more than ten years ago, and does not adequately address the human performance issues and technology changes associated with ACR's. Accordingly, a new approach to ACR safety reviews was developed based upon the concept of 'convergent validity'. This approach to ACR safety reviews is described.
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. . .