Search NASA⌕ Search

SEARCH · Search NASA

Results for “Gamma Distribution Function”

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.

At least 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Three-State Majority-vote Model on Scale-Free Networks and the Unitary Relation for Critical Exponents

Abstract We investigate the three-state majority-vote model for opinion dynamics on scale-free and regular networks. In this model, an individual selects an opinion equal to the opinion of the majority of its neighbors with probability 1 − q , and different to it with probability q . The parameter q is called the noise parameter of the model. We build a network of interactions where z neighbors are selected by each added site in the system, a preferential attachment network with degree distribution k − λ , where λ = 3 for a large number of nodes N . In this work, z is called the growth parameter. Using finite-size scaling analysis, we obtain that the critical exponents $$\beta /\bar{\nu }$$ β / ν ¯ and $$\gamma /\bar{\nu }$$ γ / ν ¯ associated with the magnetization and the susceptibility, respectively. Using Monte Carlo simulations, we calculate the critical noise parameter q c as a function of z for the scale-free networks and obtain the phase diagram of the model. We find that the critical exponents add up to unity when using a special volumetric scaling, regardless of the dimension of the network of interactions. We verify this result by obtaining the critical noise and the critical exponents for the two and three-state majority-vote model on cubic lattice networks.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

On the Prospects of Gamma-Ray Burst Detection in the TeV Band

A gamma-ray burst (GRB) jet running into an external medium is expected to generate luminous GeV–TeV emission lasting from minutes to several hours. The high-energy emission results from inverse Compton upscattering of prompt and afterglow photons by shock-heated thermal plasma. At its peak the high-energy radiation carries a significant fraction of the power dissipated at the forward shock. We discuss in detail the expected TeV luminosity, using a robust “minimal” emission model. Then, using the statistical properties of the GRB population (luminosity function, redshift distribution, afterglow energy), we simulate the expected detection rates of GRBs by current and upcoming atmospheric Cherenkov instruments. We find that GRBs exploding into a low-density interstellar medium must produce TeV emission that would have already been detected by the currently operating Cherenkov telescopes. The absence of detections is consistent with explosions into a dense wind of the GRB progenitor. If, as suggested by the recent analysis of Fermi LAT data, the typical environment of long GRBs is a Wolf–Rayet progenitor wind with the density parameter A∼10{sup 11} g cm{sup −1}, then 10%–20% of the bursts that trigger the space-borne detectors should also be detectable by the upcoming Cherenkov Telescope Array (CTA) under favorable observing conditions. Since absorption by the extragalactic background light limits the detectability above 0.1 TeV for all but the most nearby bursts (z≲1), the reduced energy threshold of CTA is the key improvement over current instruments, which should increase the number of detectable bursts by at least a factor of 3 compared with currently operating facilities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

$\gamma$-ray angular correlations, distributions and linear polarization in tracking arrays

With $\gamma $-ray tracking arrays, the angular correlation and, especially, the linear polarization of $\gamma $ rays can be measured with high precision. However, because of the complexity of these spectrometers, measuring angular distributions is more involved. The techniques and associated formalisms were developed and established in the 1950s. The experimental setups of that era were relatively simple when compared to the advanced $\gamma $-ray tracking arrays available today. This paper aims to review and refine these techniques, unify the nomenclatures and, specifically, adapt them for use with the newest tracking arrays. In addition, new methods for producing the required response functions for these arrays are introduced.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Methodology for inferring reactor core power distribution from an optical fiber based gamma thermometer array

Here, a 4-step data analytic methodology has been devised for the purpose of inferring the distribution of power in a reactor core, based on the response of an array of optical fiber based gamma thermometers (OFBGTs). This data analytic methodology is crucial for the development of a system of OFBGTs for the purpose of calibrating local power range monitors in boiling water reactors. Such a system would be an improvement to the present calibration system in boiling water reactors, in terms of safety, efficiency, and permanence. The first step of this methodology is to establish an energy balance method. In this method, one uses MCNP to determine response functions, which allow one to convert from gamma thermometer response to power. The gamma thermometers and the reactor core are segmented, such that each gamma thermometer segment provides an estimate of power for each reactor core segment. The estimates of power for the reactor core segments (hereafter referred to as fuel assembly segments) are calculated based on the measured response of the gamma thermometers. The second step of the methodology is to employ a weighting scheme to combine the estimates of the power of the various fuel assembly segments, based on the response of the various OFGBT segments and the incremental dose rates of the various fuel assembly segments to the OFBGT segments. The third step is to iteratively calculate the estimates of the power of the fuel assembly segments, until a convergence criterion is met, which indicates that the calculation has converged. The fourth and final step of the methodology is to estimate the uncertainty of the power for each fuel assembly segment. The mathematical basis for this step is not the focus of this paper. We have used a 3-D homogeneous reactor model to demonstrate the data analytic methodology; and have found that the data analytic methodology operates as intended.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

RadSim: Math, Utility & RTK

This package includes three parts, (1) gov.llnl.math, (2)gov.llnl.utility and (3)gov.llnl.rtk, which are utilized in the development of Radiation Detector Simulator (RadSim) project. RadSim is being developed to provide the capability to: (1) simulate radiation source emissions, (2) interpolate results from radiation transport tools into a common format to prepare incident flux, and (3) model radiation detector response to rapidly produce synthetic radiation measurement templates. RadSim is targeted for open-source release, which will enable researchers and industry partners to model gamma-ray detectors response to simulated flux from the transport tool of their choice. The techniques and implementation will be entirely transparent, which will allow for improvements and boutique modifications by future researchers beyond the lifespan of this specific project. The first tool of the package, gov.llnl.math, includes classes and functions to define and perform basic math operations. Some of the example features available in the package include defining statistical distributions and performing algebra and matrix operations, all of which are already accessible on publicly available software packages such as MATLAB and ROOT. The second package gov.llnl.utility includes tools commonly used to enable optimization and readability of various data structures such as Java lists and external xml files. Lastly, the gov.llnl.rtk package includes classes and functions to implement methods commonly used in radiation physics, such as data structures to represent and characterize photon spectra and tools to apply well-defined and published methods to calibrate a given spectra.

Cheung, Hoi Sing↗
DOE OSTI2024

Expected Gamma-Ray Burst Detection Rates and Redshift Distributions for the BlackCAT CubeSat Mission

We report the results of an extensive set of simulations exploring the sensitivity of the BlackCAT CubeSat to long-duration gamma-ray bursts (GRBs). BlackCAT is a NASA APRA-funded CubeSat mission for the detection and real-time subarcminute localization of high-redshift (z ≳ 3.5) GRBs. Thanks to their luminous and long-lived afterglow emissions, GRBs are uniquely valuable probes of high-redshift star-forming galaxies and the intergalactic medium. In addition, each detected GRB with a known redshift serves to localize a region of high-redshift star formation in three dimensions, enabling deep follow-on searches for host galaxies and associated local and large-scale structures. We explore two distinct models for the GRB redshift distribution and luminosity function, both consistent with Swift observations. We find that, for either model, BlackCAT is expected to detect a mean of 42 bursts per year on orbit, with 6.7% to 10% of these at z > 3.5. BlackCAT bursts will be localized to an r 90 ≲ 55'' precision and reported to the community within seconds. Due to the mission orbit and pointing scheme, bursts will be located in the night sky and well placed for deep multiwavelength follow-up observations. BlackCAT is on schedule to achieve launch readiness in 2025.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Investigation of Liquid Cloud Microphysical Properties of Deep Convective Systems: 2. Parameterization of Raindrop Size Distribution and its Application for Convective Rain Estimation

The liquid cloud microphysical properties for stratiform rain (SR) have been investigated in Part 1 of this series. Since convective rain (CR) has characteristics in raindrop size distribution (DSD) and precipitation that are distinct from SR, we investigate the CR properties in this study by using 20 hr of CR samples collected by 17 Automatic Parsivel Units disdrometers during the Midlatitude Continental Convective Clouds Experiment over the Atmospheric Radiation Measurement Southern Great Plains site. A full spectrum of DSD is constructed based on a total of 23 size channels (0.321 to 9.785 mm), and both Gamma and Exponential fitted functions are applied to extract the DSD shape parameters. Compared to SR properties, CR has distinct features including broader size range and narrower exponential slope parameter (λ E ). These results indicate that the assumption of constant exponential intercept parameter ( N 0E ) is inappropriate for CR rain rate estimates. Additionally, the subsetting scheme of the CR DSD spectra also has a strong impact on the Gamma/Exponential functions. Therefore, a new CR DSD parameterization scheme is developed by choosing the appropriate spectra subset with the constraint of rain rate and using the constant λ E instead of constant N 0E . With the input of the radar reflectivity at the lowest observed height, the newly calculated CR rain rates match well with the collocated surface rain gauge measurements (127 Mesonet stations and 17 Automatic Parsivel Units), while the rain rates calculated using traditional Z-R relationship are 3–4 times larger, indicating that constant λ E is a better assumption for CR DSD.

58 GEOSCIENCES↗
DOE OSTI2024

SCALE 6.3 Validation: Radiation Shielding

Safe and reliable use of scientific and engineering computer codes requires validation for the types of applications in which they will be used. An example in the nuclear reactor engineering and licensing field is radiation transport employed in shielding analyses. The validity of computer codes for shielding applications is demonstrated in this report for SCALE version 6.3.0. Representative benchmarks corresponding to shielding analyses are selected for the validation study. Typical measurement results analyzed from these benchmarks include neutron fluxes, detector count rates, detector energy response functions, neutron and gamma dose rates, neutron activation rates and activities, neutron leakage fluxes, and skyshine dose rates. Thousands of points of comparison between measurement and calculation are presented in this work. Other than rare outliers typically explained by either a lack of information or large uncertainties in the experiment conditions, material, or dimensions, the Monaco with Automated Variance Reduction using Importance Calculations (MAVRIC) radiation transport computer code with built-in variance reduction methods distributed with the SCALE computer code system agrees well with the measurement results. In selected benchmarks, MAVRIC is also compared to Monte Carlo N- Particle® (MCNP® ) 1 calculations. Both computer codes generally agree well within the estimated uncertainties. With the release of SCALE 6.3.0, Shift was integrated as an alternative transport solver in MAVRIC, denoted MAVRIC-Shift. Although the traditional MAVRIC using Monaco was used primarily in this validation study, many results have also been generated using MAVRIC-Shift. Agreement between MAVRIC-Monaco and MAVRIC-Shift is generally very good. The benchmarks presented in this report were obtained from reliable sources such as the International Criticality Safety Benchmark Evaluation Project Handbook, the Shielding Integral Benchmark Archive & Database, and other shielding validation work found in the literature. Additional datapoints and benchmarks will be added to future versions of this report to expand the shielding validation suite.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2022

Quark and gluon helicity evolution at small x: revised and updated

We revisit the problem of small Bjorken-x evolution of the gluon and flavor-singlet quark helicity distributions in the shock wave (s-channel) formalism. Earlier works on the subject in the same framework resulted in an evolution equation for the gluon field-strength F 12 and quark “axial current” $\overline{\psi}\gamma$ + γ 5 ψ operators (sandwiched between the appropriate light-cone Wilson lines) in the double-logarithmic approximation (summing powers of α s ln 2 (1/x) with α s the strong coupling constant). In this work, we observe that an important mixing of the above operators with another gluon operator, ${}_D{}^{\leftarrow i} {}_D{}^{i}$, also sandwiched between the light-cone Wilson lines (with the repeated transverse index i = 1, 2 summed over), was missing in the previous works. This operator has the physical meaning of the sub-eikonal (covariant) phase: its contribution to helicity evolution is shown to be proportional to another sub-eikonal operator, ${}_D{}^{i} {}_D{}^{\leftarrow i}$, which is related to the Jaffe-Manohar polarized gluon distribution. In this work we include this new operator into small-x helicity evolution, and construct novel evolution equations mixing all three operators (${}_D{}^{i} {}_D{}^{\leftarrow i}$, F 12 , and $\overline{\psi}\gamma$ + γ 5 ψ), generalizing the results of . We also construct closed double-logarithmic evolution equations in the large-N c and large-N c &N f limits, with N c and N f the numbers of quark colors and flavors, respectively. Solving the large-N c equations numerically we obtain the following small-x asymptotics of the quark and gluon helicity distributions ΔΣ and ΔG, along with the g 1 structure function, $$ \Delta \Sigma \left(x,{Q}^2\right)\sim \Delta G\left(x,{Q}^2\right)\sim {g}_1\left(x,{Q}^2\right)\sim {\left(\frac{1}{x}\right)}^{3.66\sqrt{\frac{\alpha_s{N}_c}{2\pi }}} $$ in complete agreement with the earlier work by Bartels, Ermolaev and Ryskin.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

GADRAS-DRF Validation for Safeguards and Custom Peak Fit Enhancements

In previous years, SGTech funded enhancements to the isotopics routine in the software called Gamma Detector Response and Analysis Software-Detector Response Function (GADRAS-DRF), including the addition of peak fit customization capabilities. A project was also funded that focused on implementing a peak-based model fitting routine, allowing model fitting to be performed without dependence on export-controlled cross-sections. In FY25 significant improvements were made to the custom peak fitting interface, accompanied by several validation studies within GADRAS-DRF. These studies encompassed IsotopeID performance, distributed source analysis, isotopics validation, and activity estimation. Additionally, the peak-only model fitting option was validated using an HPGe measurement of a rotating drum with line sources.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Gamma-ray spectroscopy using angular distribution of Compton scattering

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.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2018

RAINIER: A simulation tool for distributions of excited nuclear states and cascade fluctuations

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.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Forecasting constraints on the high-z IGM thermal state from the Lyman-α forest flux autocorrelation function

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.

Wolfson, Molly↗
DOE OSTI2019

Low-lying octupole isovector excitation in Nd 144

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.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Application of Gaussian Bayes classifier to differentiate chlorine-based chemical agents

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.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2018

Neutron knockout from 68,70 Ni ground and isomeric states

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.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Gamma ray Imaging Based Methods to Improve the Accuracy of Uranium Holdup Quantification

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.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗