Search NASA⌕ Search

SEARCH · Search NASA

Results for “GAMMA 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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Interpolation of computed gamma-ray detector response functions

Gamma-ray spectra measured by traditional detectors contain features that result from a combination of the effects of detector materials/geometry, the incident gamma-ray energy, and the angle of entry. The features, such as the full-energy photopeak, Compton continuum, annihilation peak, and escape peaks, are governed by simple relationships depending on incident energy and have been known for a long time. Monte Carlo computer simulations of gamma rays interacting with a detector will show these features, and with a resolution function applied, the results should look similar to real measurements. The traditional approach to creating a detector response function requires many separate simulations of monoenergetic gamma rays striking the detector. This paper presents a new approach to developing computed detector response functions. The new approach involves a much smaller number of monoenergetic gamma-ray simulations and uses interpolation to quickly generate the responses of gamma rays that were not simulated. During the interpolation process, the underlying physics equations are used to accurately compute the response of a given energy gamma ray from the small set of simulations. Such work enables accelerated generation of synthetic radiation detector data.

Detector response↗
DOE OSTI2025

Low-energy enhancement of the magnetic dipole radiation in odd-mass lanthanides

We compute the magnetic dipole (M1) $\gamma$-ray strength functions ($\gamma$SF) for the odd-mass lanthanides $^{\textrm{143-151}}$Nd and $^{\textrm{147-153}}$Sm using the shell-model Monte Carlo method in combination with the static-path approximation and the maximum-entropy method. In particular, we quantify the statistical uncertainties in the calculated M1 $\gamma$SFs and show that they are under control for the excitation energies relevant to the experiments despite a Monte Carlo sign problem that originates in the projection onto an odd number of neutrons. We identify a low-energy enhancement (LEE) in the M1 $\gamma$SFs of these odd-mass lanthanides, which was recently observed experimentally in some of them. We also find a scissors mode resonance (SR) in the strongly deformed isotopes. We observe that the decrease in the LEE strength with neutron number along an isotopic chain is compensated for by an increase in the SR strength in the deformed nuclei. Furthermore, we compare our results with recent experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

New determination of the millisecond pulsar gamma-ray luminosity function and implications for the Galactic Center gamma-ray excess

It has been suggested that the Galactic Center gamma-ray excess (GCE) could be produced by a large number of centrally located millisecond pulsars. The fact that no such pulsar population has been detected implies that these sources must be very faint and very numerous. Here, in this study, we use the contents of Fermi’s recently released Third Pulsar Catalog (3PC) to measure the luminosity function of the millisecond pulsars in the Milky Way’s disk. We find that this source population exhibits a luminosity function with a mean γ -ray luminosity of ⟨ L γ ⟩ ∼ 6 × 10 32 erg / s (integrated above 0.1 GeV). If the GCE were generated by millisecond pulsars with the same luminosity function, we find that ∼ 20 such sources from the inner Galaxy population should have already been detected by Fermi and included in the 3PC. Given the lack of such observed sources, we exclude the hypothesis that the GCE is generated by pulsars with the same luminosity function as those in the Galactic disk with a significance of 3.4 σ . We conclude that either less than 39% of the GCE is generated by pulsars, or that the millisecond pulsars in the inner Galaxy are at least 5 times less luminous on average than those found in the Galactic disk.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Maximum Likelihood Estimation: Some Basics

The maximum likelihood estimation is a general estimation procedure. It is often compared to estimation procedures like the ordinary least squares regression or generalized method of moments, to name a few. We discuss some basics about the maximum likelihood estimation, its advantages and disadvantages, and provide an example application to a gamma distribution function.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Maximum Likelihood Estimation: Some Basics

The maximum likelihood estimation is a general estimation procedure. It is often compared to estimation procedures like the ordinary least squares regression or generalized method of moments, to name a few. We discuss some basics about the maximum likelihood estimation, its advantages and disadvantages, and provide an example application to a gamma distribution function.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Magnetic dipole γ-ray strength functions in the crossover from spherical to deformed neodymium isotopes

We calculate the magnetic dipole $\gamma$-ray strength functions in a chain of even-mass neodymium isotopes $^{144-152}$Nd in the framework of the configuration-interaction (CI) shell model. We infer the strength function by applying the maximum entropy method (MEM) to the exact imaginary-time response function calculated with the shell-model Monte Carlo (SMMC) method. The success of the MEM depends on the choice of a good strength function as a prior distribution. We investigate two choices for the prior strength function: the static path approximation (SPA) and the quasiparticle random-phase approximation (QRPA). We find that the QRPA is a better approximation at low temperatures (i.e., near the ground state), while the SPA is a better choice at finite temperatures. We identify a low-energy enhancement (LEE) in the MEM deexcitation $M1$ strength functions of the even-mass neodymium isotopes and compare with recent experimental results for the total deexcitation $\gamma$-ray strength functions. The LEE is already seen in the SPA strength function but not in the QRPA strength function, indicating the importance of large-amplitude static fluctuations around the mean field in reproducing the LEE. Our method is currently the only one which can reproduce LEE in heavy open-shell nuclei where conventional CI shell model calculations are prohibited. With the onset of deformation as number of neutrons increases along the chain of neodymium isotopes, we observe that some of the LEE strength transfers to a low-energy excitation, which we interpret as a finite-temperature ``scissors'' mode. Here, we also observe a finite-temperature spin-flip mode.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Measurement of the 70 Ge( d , n ) 71 As cross section from 2.1 MeV to 10 MeV

This work reports experimental 70 Ge(d,n) 71 As cross sections producing Arsenic-71 (t 1/2 = 65.3 h, 28% β + ), a potentially useful diagnostic radionuclide. Target stacks containing two nat Ge foils, a nat Ni monitor foil, and an Al degrader were irradiated with 5.5–10 MeV deuterons using an FN Tandem Van de Graaff accelerator. 71 As and 61 Cu, 56 Co, and 58 Co (from nat Ni(d,x)) were quantified in the foils by high purity germanium gamma spectrometry. Furthermore, the measured cross sections are larger than previously reported and suitable for producing preclinical quantities of diagnostic radiopharmaceuticals.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2025

Switching kinetics in ferroelectric zinc magnesium oxide thin films

The switching kinetics of RF-magnetron reactively sputtered ~200 nm thick Zn 1-x Mg x O (ZMO) ferroelectric thin films with x = 0.41 and x = 0.27 prepared on Pt/Ti/SiO 2 /Si substrates were studied at applied fields near the coercive fields, ranging from 3.4 to 5.1 MV cm -1 , and at temperatures ranging from room temperature to 100 °C. Polarization reversal in ZMO followss the Kolmogorov-Avrami-Ishibashi kinetics model for nucleation and growth at applied fields near the coercive field, and obeys an individual column switching (ICS) model at higher applied fields required to switch most of the spontaneous polarization. Switching in the high-applied field regime can be described using either the Gaussian or inverse gamma distribution functions depending on the mole fraction of magnesium in the film. The switching current transients of the high-Mg content ZMO film (x = 0.41) are always bi-modal, whereas the low Mg composition film (x = 0.27) is described by the Gaussian distribution initially following wake-up, and becomes bimodal with continued cycling. Rayleigh-like behavior of the dielectric constant revealed a simultaneous increase of the irreversible and decrease of the reversible contributions to the dielectric constant, which was ascribed to an increase in the density of mobile domain walls with cycling and resulted in faster switching.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2024

QED corrections to the thermal neutrino interaction rate

Abstract Motivated by precision computations of neutrino decoupling at MeV temperatures, we show how QED corrections to the thermal neutrino interaction rate can be related to the electron-positron spectral function as well as an effective$$\overline{\nu }\nu \gamma $$vertex. The spectral function is needed both in a timelike and in a spacelike domain, and for both of its physical polarization states (transverse and longitudinal with respect to spatial momentum). Incorporating an NLO evaluation of this spectral function, an estimate of the$$\overline{\nu }\nu \gamma $$vertex, and HTL resummation of scatterings mediated by soft Bose-enhancedt-channel photons, we compute the interaction rate as a function of the neutrino momentum and flavour. Effects on the – (0. . .2)% level are found, noticeably smaller than a previous estimate of a related quantity.

Physics↗
DOE OSTI2026

Evaluation of simulated HPGe detector efficiencies in OpenMC compared to MCNP

In late 2023, the open-source radiation transport code OpenMC introduced a pulse-height tally (PHT) feature, enabling users to track the total energy deposited by individual photons in OpenMC cells. This function represents an important improvement of the OpenMC code because it allows users to simulate the response function of a gamma detector without relying on closed-source alternatives. Despite this, limited work has been published evaluating OpenMC as a radiation transport code for simulating gamma spectroscopy experiments. This study attempts to demonstrate the usefulness of OpenMC in this space by directly comparing its PHT output to MCNP (a trusted industry-standard Monte Carlo code). In the first half of this study, a radiation transport experiment is described in which a detector with a complex internal geometry is exposed to various gamma-emitting isotopes ( 133 Ba, 137 Cs, and 60 Co) over a range of distances. The setup of this experiment was modeled in MCNP with sufficient detail to capture the efficiency characteristics of a high purity germanium (HPGe) detector for the primary gammas of each isotope. After good agreement between the MCNP model and experiment was achieved, an identical model was produced in OpenMC to allow direct comparison between the full energy peak (FEP) values produced in OpenMC and MCNP. The results show strong agreement between OpenMC and MCNP across the full range of tested energies, with each model’s FEP values typically within 2% of each other and most FEP areas within 4% of experimental data. Worse agreement was seen between the Monte Carlo codes and experiment below 300 keV (an expected result). For the 662 keV line of 137 Cs, both codes were found in poor agreement with experiment and each other over the full range of distances tested (possibly indicating an error with the 137 Cs experimental data). Ignoring the anomalous results of the 662 keV line, all other data show good qualitative and quantitative agreement between MCNP and OpenMC. This result demonstrates the accuracy of OpenMC’s PHT feature for spectroscopic applications in which detector efficiency is a primary concern.

07 - ISOTOPES AND RADIATION SOURCES↗
DOE OSTI2025

Convective shells in the interior of Cepheid variable stars: Overshooting models based on hydrodynamic simulations

Context. Because Cepheid variable stars have long been used as a cosmic benchmark for scaling distances in our Galaxy and beyond, the accuracy of stellar evolution models for Cepheids have wide-reaching effects. However, our understanding of the dynamics in the interiors of these physically complex stars is limited. Aims. Our goal is to provide a detailed multi-dimensional picture of hydrodynamic convection and convective boundary mixing in the interior of Cepheids. Methods. Using the Modules for Experiments in Stellar Astrophysics (MESA), we studied the structure of intermediate-mass stars that cross the instability strip. Then, we performed two-dimensional hydrodynamic simulations of six stars with the fully compressible Multidimensional Stellar Implicit Code (MUSIC). Our simulations did not model the radial pulsations but focused on the interior structure of this family of stars. We developed and applied a new statistical analysis to examine convection and convective boundary mixing in the interior of these stellar simulations. Results. Based on a grid of MESA models, we demonstrated that a common structure for intermediate mass Cepheids includes an interior convective shell as well as a thin outer convective envelope. Using the extreme value theory approach to analyze our MUSIC simulation data, we found that overshooting above the convective shell fills the space between these convectively unstable layers. We developed a new statistical analysis that provides a clearer picture of how overshooting fills this layer; it also allowed us to formulate a detailed comparison between overshooting above and below the convective shell. Our analysis effectively decomposes the overshooting layer into two layers: a weak overshooting layer and a strong overshooting layer. Statistically, this is accomplished by decomposing the strongly non-Gaussian probability density function into a mixture of gamma distributions. Using our mixture model, we showed that the ratio of overshooting lengths above and below the convective shell depends directly on the radial extent of the convective shell as well as its depth in the star. We proposed a new form for the diffusion coefficient that addresses the need for overlapping overshooting layers between convective shells. We introduced the idea of a “super-mixing layer” where overshooting from both the convective shell and the convective envelope results in efficient mixing and could be viewed as merging the two adjacent convective zones.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Solid, structured composite neutron detectors with high dynamic range capability

Neutron detectors are essential across disciplines such as fundamental science, nuclear security, safeguards, and civilian applications. While 3 He-filled gas proportional counters have long been revered for their efficacy in detecting thermal neutrons and praised for their efficiency, neutron/gamma discrimination, and stability, the scarcity of 3He has spurred a search for alternatives. Here, we explore a solid structured scintillating particle composite (SPC) consisting of 6 Li-containing scintillating glass particles within an acrylic matrix as a neutron detector for high dynamic range applications. We show for the first time that an SPC neutron detector can boast an intrinsic detection efficiency of 0.261% for pure 252 Cf fission neutrons and an overall neutron detection efficiency of (0.546 ± 0.003)% at the Neutron Free-in-Air facility while being able to function in an intense gamma-ray environment. We also show that the SPC neutron detector supports fast neutron capture times and enables a dual-readout scheme that extends the detector dynamic range to high incident neutron fluxes. A scalable fabrication process allows for tailoring the SPC detector properties to the requirements of specific applications. Good agreement is found between the experimental results taken with a National Institute of Standards and Technology traceable 252 Cf source and the coupled MCNP6 and optical-ray-tracing simulations.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Origin of the Low-Energy Enhancement of the γ-Ray Strength Function

The low-energy enhancement (LEE) in γ-ray strength functions has been experimentally identified in a large number of nuclei during the past two decades; however, the origin of the enhancement is not fully understood. Building on previous theoretical work, we investigate the LEE and its relation to the scissors mode (SM) with an independent theoretical approach. We apply a novel angular-momentum-projected shell-model method that explicitly endows degrees of freedom to describe the scissors motion. Taking the recently measured γ-ray strength functions in Neodymium isotopes as examples, we find that the LEE arises from a quasi-free scissors motion appearing only in weakly-deformed nuclei, which can be viewed as an approximate free-rotation of neutrons with respect to protons. This leads us to propose a new type of collective motion, scissors rotation, to contrast the scissors vibration widely known in well-deformed nuclei. The observed LEE is naturally interpreted as the first evidence for this collective excitation mode.

Collective models↗
DOE OSTI2026

Extraction of neutron-capture cross sections on 92 Zr using the charge-exchange Oslo method

The 93 Nb (𝑡, 3 He ) reaction at 115 MeV/nucleon was studied to demonstrate that nuclear level densities and 𝛾-ray strength functions can be extracted from charge-exchange reactions at intermediate energies using the Oslo technique. The matrix of excitation energy in 93 Zr, reconstructed from the (𝑡, 3 He ) reaction, versus the energy of 𝛾 rays emitted by the excited 93 Zr nuclei, was obtained in an experiment with the S800 Spectrograph operated in coincidence with the GRETINA 𝛾 -ray detector. The extracted level density and 𝛾 -ray strength function obtained by applying the Oslo method to this matrix were used to estimate the 92 Zr⁢(𝑛,𝛾) ⁢93 Zr cross section by combining the new results with other experimental data and theoretical calculations for 𝐸⁢1 and 𝑀⁢1 strength functions at higher energies. Good agreement with direct measurements of the 92 Zr⁢(𝑛,𝛾)⁢ 93 Zr cross section was found. The contribution from the upbend in the extracted 𝛾-ray strength function was important to achieve the consistency, as the neutron-capture cross section without this contribution is significantly below the direct measurements otherwise. Since charge-exchange reactions at intermediate energies have long been used for extracting Gamow-Teller strengths, the successful demonstration of the charge-exchange Oslo method enables experiments in which (𝑛,𝛾) cross sections and Gamow-Teller strengths can be measured simultaneously, which is of benefit for astrophysical studies.

90 ≤ A ≤ 149↗
DOE OSTI2025

First simultaneous measurement of the γ-ray and neutron emission probabilities in inverse kinematics at a heavy-ion storage ring

The probabilities for γ-ray and particle emission as a function of the excitation energy of a decaying nucleus are valuable observables for constraining the ingredients of the models that describe the deexcitation of nuclei near the particle emission threshold. These models are essential in nuclear astrophysics and applications. In this paper, we have for the first time simultaneously measured the γ-ray and neutron emission probabilities of 208 Pb. The measurement was performed in inverse kinematics at the Experimental Storage Ring (ESR) of the GSI/FAIR facility, where a 208 Pb beam interacted through the 208 Pb (p,p') reaction with a hydrogen gas jet target. Instead of detecting the γ rays and neutrons emitted by 208 Pb, we detected the heavy beamlike residues produced after γ and neutron emission. These heavy residues were fully separated by a dipole magnet of the ESR and were detected with outstanding efficiencies. The comparison of the measured probabilities with model calculations has allowed us to test and select different descriptions of the γ-ray strength function and the nuclear level density available in the literature.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Deviations from the Porter-Thomas Distribution due to Nonstatistical 𝛾 Decay below the 150 Nd Neutron Separation Threshold

We introduce a new method for the study of fluctuations of partial transition widths based on nuclear resonance fluorescence experiments with quasimonochromatic linearly polarized photon beams below particle separation thresholds. It is based on the average branching of decays of 𝐽=1 states of an even-even nucleus to the 2$^{+}_{1}$ state in comparison to the ground state. Between 5 and 7 MeV, a constant average branching ratio for 𝛾 decays from 1 − states of 0.490(16) is observed for the nuclide 150 Nd. Assuming 𝜒 2 -distributed partial transition widths, this average branching ratio is related to a degree of freedom of 𝜈 = 1.93⁢(12), rejecting the validity of the Porter-Thomas distribution, requiring 𝜈 = 1. The observed deviation can be explained by nonstatistical effects in the 𝛾-decay behavior with contributions in the range of 9.4(10)% up to 94(10)%.

150 ≤ A ≤ 189↗
DOE OSTI2025

Ab initio prediction of FeCr sigma (001)/FCC Fe(111) interfacial energy: Effect of interfacial doping of C, B, and N

Interfacial energy plays a crucial role for high-temperature material systems in determining their microstructure evolution and mechanical properties. Optimizing the interfacial energy through doping can improve the stability, performance, and overall functionality of these materials. Here, the present study delivers theoretical investigation on the sigma FeCr (001)//FCC Fe (111) interfacial energy using ab initio methods, focusing on the implications of interfacial doping with important interstitial elements carbon (C), boron (B), and nitrogen (N). The calculated interfacial energy without doping is 0.183J/m 2 . Upon doping, notable reductions in interfacial energy were observed. Doping with B decreased the interfacial energy to 0.046J/m 2 , with single C doping to 0.126J/m 2 , with two C doping to 0.098J/m 2 , and with a single N doping to 0.071J/m 2 . These findings not only shed light on the atomistic origin for the effect of stabilizing the microstructure of stainless steel at elevated temperatures by interstitial doping, but also present a systematic, thorough ab initio approach to predict interfacial doping properties and guide alloy design to enhance the performance and stability.

FeCr sigma phase↗