Uncertainty analysis of Response Functions and y-Backgrounds on Tion and t0 Measurements from Cherenkov Neutron Detectors at the National Ignition Facility (NIF) a)
Explore the source record for details and available documents.
SEARCH · Search NASA
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.
Explore the source record for details and available documents.
We present a novel method to perform gamma-ray spectroscopy suitable for extreme pileup situations, wherein traditional pulse-height spectroscopy is infeasible or impractical. One example of such application is spectral characterization of intense laser-based radiation sources that employ laser wakefield acceleration, in which a large number of photons can be produced in a well-collimated beam on a picosecond or shorter timescale. The method relies on the angularly-resolved calorimetric measurement of Compton scattering. The probability for a photon to undergo Compton scattering into a given solid angle is a function of its energy and is described by the Klein–Nishina formula. By first Compton scattering the incident beam and then measuring the total energies deposited in multiple detectors placed at various scattering angles, the incident beam energy spectrum can be reconstructed without resolving the individual events in detectors. Here, we compare the performance of an maximum likelihood expectation maximization algorithm with that of an artificial neural network to measure the energy spectrum of a beam of monoenergetic gamma rays from 137Cs decay. We demonstrate that the method can be used to reconstruct the characteristic inverse Compton scattering spectra.
In this paper, a new code has been developed named RAINIER that simulates the γ-ray decay of discrete and quasi-continuum nuclear levels for a user-specified range of energy, angular momentum, and parity including a realistic treatment of level spacing and transition width fluctuations. A similar program, DICEBOX, uses the Monte Carlo method to simulate level and width fluctuations but is restricted in its initial level population algorithm. On the other hand, modern reaction codes such as TALYS and EMPIRE populate a wide range of states in the residual nucleus prior to γ-ray decay, but do not go beyond the use of deterministic functions and therefore neglect cascade fluctuations. This combination of capabilities allows RAINIER to be used to determine quasi-continuum properties through comparison with experimental data. Finally, several examples are given that demonstrate how cascade fluctuations influence experimental high-resolution γ-ray spectra from reactions that populate a wide range of initial states.
ABSTRACT The autocorrelation function of the Lyman-$\alpha$ (Ly $\alpha$) forest flux from high-z quasars probes the small-scale structure of the intergalactic medium (IGM). The thermal state of the IGM, determined by the physics of reionization, sets the small-scale power observed in the Ly $\alpha$ forest. To explore the sensitivity of the autocorrelation function to the IGM’s thermal state, we compute the autocorrelation function from a cosmological hydrodynamical simulation with an instantaneous reionization model and 135 post-processed thermal states. Using mock data sets of 20 quasars, we forecast constraints on $T_0$ and $\gamma$, which characterize the post-processed IGM thermal state, at $5.4 \le z \le 6$. While this model simplifies the IGM’s thermal state, it serves as a key first step in assessing future observational prospects. We also perform an inference test on mocks and re-weight out posterior distributions to guarantee that they exhibit statistically correct behaviour. At $z = 5.4$, we find that an idealized data set constrains $T_0$ to 59 per cent and $\gamma$ to 16 per cent at the 1$\sigma$ equivalent confidence level. To explore more realistic, non-instantaneous reionization scenarios, we analyse four models combining temperature and ultraviolet background (UVB) fluctuations at $z = 5.8$. We find that mock data generated from a model with both temperature and UVB fluctuations can rule out a model with only temperature fluctuations at the $> 1\sigma$ level 73.9 per cent of the time.
The nature of low-lying 3 – levels in 144 Nd was investigated in the 143 Nd(n, γγ) cold neutron-capture reaction. The combination of the high neutron flux from the research reactor at the Institut Laue-Langevin and the high γ-ray detection efficiency of the EXILL setup allowed the recording of γγ coincidences. From the coincidence data precise branching ratios were extracted. Furthermore, the octagonal symmetry of the setup allowed angular-distribution measurements to determine multipole-mixing ratios. Additionally, in a second measurement the ultra-high resolution spectrometer GAMS6 was employed to conduct lifetime measurements using the gamma-ray induced Doppler-shift technique (GRID). The confirmed strong M1 component in the 3 – 3 → 3 – 1 decay strongly supports the assignment of the 3 – 3 level at 2779 keV as low-lying isovector octupole excitation. Here, the microscopic calculations within the quasiparticle phonon model confirm an isovector component in the wave function of the 3 – 3 level, firmly establishing this fundamental mode of nuclear excitation in near-spherical nuclei.
The Portable Isotopic Neutron Spectroscopy (PINS) is a commercialized system developed by Idaho National Laboratory (INL) to examine chemical warfare agents (CWA) non-destructively, utilizing Prompt Gamma Neutron Activation Analysis (PGNAA) techniques. The PINS system takes advantage of a high-resolution gamma-ray spectrum from a mechanically-cooled high-purity germanium (HPGe) detector. One of the difficult technical challenges is to discriminate the chlorine-based chemical agents. Especially, CN, CNB, CNS and CG have similar chemical compositions to make it hard to discriminate them with a higher confidence. Current identification algorithms for PINS systems with 252-Cf sources have been improved and updated continuously as more field data became available, and new algorithms was studied to complement the current algorithms by adopting the Gaussian Bayes classifier. These new algorithms were intended to be applied to a subset of chlorine-based chemical agents, and their main goal is discriminate CN, CNB, CNS and CG with their ratios of the chlorine neutron inelastic 1763keV peak to the chlorine thermal neutron capture 1959keV peak, which is referred to as the “Cl i/c” or “clic” ratio in this study. The Cl i/c ratios were assumed to follow Gaussian distributions with the means and the standard deviations unique to their corresponding chemical agents. The prior probabilities of these four chemical agents were optimized with a collection of field data to achieve the best performance in terms of precision or positive predictive value (PPV). Finally, their posterior probabilities as functions of the Cl i/c ratio were implemented in the current version of PINS analysis software in order to be tested with more field data.
Neutron-rich isotopes are an important source of new information on nuclear physics. Specifically, the spin-isospin components in the nucleon-nucleon (NN) interaction, e.g., the proton-neutron tensor force, are expected to modify shell structure in exotic nuclei. These potential changes in the intrinsic shell structure are of fundamental interest. The study of the excitation energy of states corresponding to specific configurations in even-even isotopes, together with the single-particle character of the first excited states of odd-A, neutron-rich Ni isotopes, probes the evolution of the neutron orbitals around the Fermi surface as a function of the neutron number a step forward in the understanding of the region and the nature of the NN interaction at large N/Z ratios. In an experiment carried out at the National Superconducting Cyclotron Laboratory, new spectroscopic information was obtained for 68 Ni and the distribution of single-particle strengths in 67,69 Ni was characterized by means of single-neutron knockout from 68,70 Ni secondary beams. The spectroscopic strengths, deduced from the measured partial cross sections to the individual states tagged by their de-exciting gamma rays, is used to identify and quantify configurations that involve neutron excitations across the N = 40 harmonic oscillator shell closure. Finally, the de-excitation rays were measured with the GRETINA tracking array.
Deposits in commercial uranium processing facilities of interest to nuclear safeguards are predominantly low enriched uranium (LEU). The shapes and sizes of holdup deposits can vary significantly, thus making it difficult to model the deposits accurately for calibration purposes. This presents a challenge to currently employed methods such as the Generalized Geometry Holdup (GGH) that rely on quantifying 235U mass by simplifying deposit shapes as a point, a line, or an area. In this work, gamma ray imaging using high energy resolution Germanium Gamma Imagers (GeGI™) is employed to determine the distribution of uranium inside the source containment. Two types of imaging methods are employed; coded-aperture imaging and Compton imaging. Gamma-ray emissions seen in different regions of the image are quantified using an inverse gamma-ray transport solver being developed at ORNL. Using measured data, the inverse method will be used to solve for unknown source parameters such as source matrix thickness, density, and attenuation due to container wall and shielding. The intrinsic efficiency of the detector is determined based on a library of built-in response functions and key dimensions of the detector (e.g., thickness and radius of the detector crystal). Measurements are performed in both imaging modalities using uranium sources of well-known masses and enrichments configured inside mocked-up holdup fixtures. Results from the two imaging modalities are intercompared. The paper presents preliminary results of the imaging measurements and the progress made in the computational methods to quantify uranium mass.
This study reports on the quantification of deuterium in ion-irradiated gamma-LiAlO2 pellets as a function of dose and temperature. The gamma-LiAlO2 pellets were sequentially irradiated with He+ and D2+ ions to the same fluences of 5E16, 1E17 and 2E17 He++D+/cm2 at 188 K. Additional irradiation was performed to 1E17, 2E17 and 4E17 He++D+/cm2 at 573 K. A set of the pellets irradiated at 188 K was shipped and stored at low temperatures from 80 to 132 K and characterized using time-of-flight secondary ion mass spectrometry at ~173 K. The deuterium depth profiles show a Gaussian-like distribution in the low-temperature pellets. The total deuterium retention is found to be directly proportional to the ion fluence. About 27 at.% of the implanted deuterium atoms were released from the pellet irradiated to 2E17 He++D+/cm2 at 188 K during storage at room temperature for ~1 month. Retention of the trapped or bound deuterium during ion irradiation at 573 K increases initially with ion fluence and tends to saturate at a high fluence. The amount of the released deuterium is observed to be quantitatively consistent with that of the released tritium from similar standard pellets during neutron irradiation at 573 K.
We study formation of the Vibrational Distribution Function (VDF) in a molecular gas at low pressure, when vibrational levels are excited by electron impact and deactivated in collisions with walls and show that this problem has a convenient analytical solution that can be used to obtain VDF and its dependence on external parameters. The VDF is determined by excitation of vibrational levels by an external source and deactivation in collisions with the wall. Deactivation in wall collisions is little known process. However, we found that the VDF is weakly dependent on the functional form of the actual form of probability ${\gamma }_{\upsilon ^{\prime} \to \upsilon }$ for a vibrational number $\upsilon ^{\prime} $ to transfer into a lower level $\upsilon $ at the wall. Because for a given excitation source of vibrational states, the problem is linear the solution for VDF involves solving linear matrix equation. The matrix equation can be easily solved if we approximate probability, in the form: ${\gamma }_{\upsilon ^{\prime} \to \upsilon }=(1/\upsilon ^{\prime} )\theta (\upsilon ^{\prime} -\upsilon ).$ In this case, the steady-state solution for VDF($\upsilon $) simply involves a sum of source rates for levels above $\upsilon ,$ with a factor of $1/(\upsilon +1).$ As an example of application, we study the vibrational kinetics in a hydrogen gas and verify the analytical solution by comparing with a full model.
A Gamma-Ray and Neutron Spectrometer (GRNS) instrument has been developed as part of the science payload for NASA’s Discovery Program Psyche mission to the M-class asteroid (16) Psyche. The GRNS instrument is designed to measure the elemental composition of Psyche with the goal to understand the origin of this mysterious, potentially metal-rich planetary body. The GRNS will measure the near-surface abundances for the elements Ni, Fe, Si, K, S, Al, and Ca, as well as the spatial distribution of Psyche’s metal-to-silicate fraction (or metal fraction). These measurements address three of the five Psyche mission science objectives: determine if Psyche is a core; determine whether small metal bodies incorporate light elements into the metal phase; and determine whether Psyche was formed under reducing conditions. The Gamma-Ray Spectrometer (GRS) uses a cryocooled, high-purity Ge (HPGe) sensor to detect cosmic-ray generated gamma rays in the 60 to 9000-keV energy range. The HPGe sensor is surrounded by a borated plastic anticoincidence shield that provides three functions: active background rejection from charged particle interactions in the HPGe sensor; fast neutron measurements; and direct measurements of the incident galactic cosmic ray flux. The Neutron Spectrometer (NS) uses three 3 He gas proportional sensors, each with different material wraps to measure thermal (<0.4 eV), low-energy epithermal (0.4 eV to 1 keV), and high-energy epithermal (up to 100 keV) neutrons. This paper provides an overview of the Psyche GRNS, including: its science and measurement objectives; the design of the instrument hardware, software, and operation; pre-launch performance measurements and its initial performance in space; and an overview of its data products and expected operation for different Psyche mission phases.
GRB 160709A is one of the few bright short gamma-ray bursts detected by both the Gamma-ray Burst Monitor and the Large Area Telescope on board the Fermi Gamma-ray Space Telescope. The γ-ray prompt emission of GRB 160709A is adequately fitted by combinations of three distinct components: (i) a nonthermal component described by a power law (PL) with a high-energy exponential cutoff, (ii) a thermal component modeled with a Planck function, and (iii) a second nonthermal component shaped by an additional PL crossing the whole γ-ray spectrum. While the thermal component dominates during ~0.12 s of the main emission episode of GRB 160709A with an unusually high temperature of ~340 keV, the nonthermal components dominate in the early and late time. The thermal component is consistent with the photospheric emission resulting in the following parameters: the size of the central engine, ${R}_{0}={3.8}_{-1.8}^{+5.9}$ × 10 8 cm, the size of the photosphere, ${R}_{\mathrm{ph}}={7.4}_{-1.2}^{+0.8}$ × 10 10 cm, and a bulk Lorentz factor, ${\rm{\Gamma }}={728}_{-93}^{+75}$, assuming a redshift of 1. The slope of the additional PL spectrum stays unchanged throughout the burst duration; however, its flux decreases continuously as a function of time. A standard external shock model has been tested for the additional PL component using the relation between the temporal and spectral indices (the closure relation). Each set of spectral and temporal indices from two energy bands (200 keV–40 MeV and 100 MeV–10 GeV) satisfies a distinct closure relation. From the closure relation test we derived the index for the electron spectral distribution, p = 2.5 ± 0.1. The interaction of the jet with the interstellar environment is preferred over the interaction with the wind medium.
We report the discovery of 1.97 ms period gamma-ray pulsations from the 75 minute orbital-period binary pulsar now named PSR J1653-0158. The associated Fermi Large Area Telescope gamma-ray source 4FGL J1653.6-0158 has long been expected to harbor a binary millisecond pulsar. Despite the pulsar-like gamma-ray spectrum and candidate optical/X-ray associations—whose periodic brightness modulations suggested an orbit—no radio pulsations had been found in many searches. The pulsar was discovered by directly searching the gamma-ray data using the GPU-accelerated Einstein@Home distributed volunteer computing system. The multidimensional parameter space was bounded by positional and orbital constraints obtained from the optical counterpart. More sensitive analyses of archival and new radio data using knowledge of the pulsar timing solution yield very stringent upper limits on radio emission. Any radio emission is thus either exceptionally weak, or eclipsed for a large fraction of the time. The pulsar has one of the three lowest inferred surface magnetic-field strengths of any known pulsar with B surf ≈ 4 × 10 7 G. The resulting mass function, combined with models of the companion star's optical light curve and spectra, suggests a pulsar mass gsim 2 M ⊙ . The companion is lightweight with mass ~0.01 M ⊙ , and the orbital period is the shortest known for any rotation-powered binary pulsar. This discovery demonstrates the Fermi Large Area Telescope's potential to discover extreme pulsars that would otherwise remain undetected.
Despite the widespread application of statistical downscaling tools, uncertainty remains regarding the role of model formulation in determining model skill for daily maximum and minimum temperature (Tmax and Tmin), and precipitation occurrence and intensity. Impacts of several key aspects of statistical transfer function form on model skill are evaluated using a framework resistant to model overspecification. We focus on: (a) model structure: simple (generalized linear models, GLMs) versus complex (artificial neural networks, ANNs) models. (b) Predictor selection: Fixed number of predictors chosen a priori versus stepwise selection of predictors and inclusion of grid point values versus predictors derived from application of principal components analysis (PCA) to spatial fields. We also examine the influence of domain size on model performance. For precipitation downscaling, we consider the role of the threshold used to characterize a wet day and apply three approaches (Poisson and Gamma distributions in GLM and ANN) to downscale wet-day precipitation amounts. While no downscaling formulation is optimal for all predictands and at 10 locations representing diverse U.S. climates, and due to the exclusion of variance inflation all of the downscaling formulations fail to reproduce the range of observed variability, models with larger suites of prospective predictors generally have higher skill. For temperature downscaling, ANNs generally outperform GLM, with greater improvements for T min than T max . Furthermore, use of PCA-derived predictors does not systematically improve model skill, but does improve skill for temperature extremes. Model skill for precipitation occurrence generally increases as the wet-day threshold increases and models using PCA-derived predictors tend to outperform those based on grid cell predictors. Each model for wet-day precipitation intensity overestimates annual total precipitation and underestimates the proportion derived from extreme precipitation events, but ANN-based models and those with larger predictor suites tend to have the smallest bias.
Quantifying the radioactive sources present in gamma spectra is an ever-present and growing national security mission and a time-consuming process for human analysts. While machine learning models exist that are trained to estimate radioisotope proportions in gamma spectra, few address the eventual need to provide explanatory outputs beyond the estimation task. In this work, we develop two machine learning models for a NaI detector measurements: one to perform the estimation task, and the other to characterize the first model’s ability to provide reasonable estimates. To ensure the first model exhibits a behavior that can be characterized by the second model, the first model is trained using a custom, semi-supervised loss function which constrains proportion estimates to be explainable in terms of a spectral reconstruction. The second auxiliary model is an out-of-distribution detection function (a type of meta-model) leveraging the proportion estimates of the first model to identify when a spectrum is sufficiently unique from the training domain and thus is out-of-scope for the model. In demonstrating the efficacy of this approach, we encourage the use of meta-models to better explain ML outputs used in radiation detection and increase trust.
Very hot plasmas with ion temperature exceeding 10 10 K can be formed in certain astrophysical environments. The distinct radiation signature of such plasmas is the γ-ray emission dominated by the prompt deexcitation nuclear lines and π 0 -decay γ rays. Using a large nuclear reaction network, we compute the time evolution of the chemical composition of such hot plasmas and their γ-ray line emissivity. At higher energies, we provide simple but accurate analytical presentations for the π 0 -meson production rate and the corresponding π 0 → 2γ emissivity derived for the Maxwellian distribution of protons. Furthermore, we discuss the impact of the possible deviation of the high energy tail of the particle distribution function from the “nominal” Maxwellian distribution on the plasma γ-ray emissivity.
The planned conversion of the High Flux Isotope Reactor (HFIR) at Oak Ridge National Laboratory from highly enriched uranium (HEU) to low-enriched uranium (LEU) fuel requires detailed evaluation of experiment-relevant parameters to ensure continued performance for materials testing and isotope production. Here, this study presents the first comprehensive assessment of displacements per atom (dpa) and heat deposition rates in target materials within the HFIR flux trap with both HEU and candidate LEU core configurations. Seven analyses were conducted to evaluate key performance metrics, including fast neutron flux distribution, cross section response functions, cross section data, and local dpa and heat deposition rates using mesh- and cell-based tallies. Simulations employed Shift, Monte Carlo N-Particle (MCNP), and the HIFR Controller (HFIRCON) tool suite for high-fidelity transport and depletion modeling. The LEU designs—using U 3 Si 2 -Al dispersion fuel and operating at 95 MW—were compared to the current 85 MW HEU configuration. Results show that while the candidate LEU cores exhibit higher dpa rates due to a harder spectrum and extended cycle lengths, they also demonstrate reduced heat deposition rates in irradiation experiments, primarily due to increased gamma self-shielding from higher 238 U content in the core. These findings confirm that LEU conversion can maintain HFIR’s materials irradiation capabilities but may require redesigning existing experimental hardware.
Here, this study explores the relationship between experimentally observed inter-event time-interval distributions (TIDs) and measured dead times for pulses generated by a HPGe gamma detector and processed by a digital spectrometer. The system utilizes a fast and slow channel, pile-up rejector, trapezoid filtering, and flash analog-to-digital converter. The experimentally derived TIDs were compared with theory for validation. The results demonstrate that the theoretical model reliably describes measured TIDs up to 40% dead time. However, significant distortion effects become increasingly pronounced at higher input rates. It appears that the deviation between the measured and calculated time-interval distributions take the shape of higher-order convolutions of the TID, which represent the time difference between counts for more than two successive events. In this work, theoretical functions for the TID are expanded to reproduce measurements up to 90%. This refinement in the interpretation and treatment of the measured TIDs provides improved accuracy and precision in the prediction of the true event rate and measured dead time in the counter. Even though it reflects occasional count loss that is specific to the set-up, it is expected that a similar adjustment may be applicable to other detection systems as well.