Search NASA⌕ Search

SEARCH · Search NASA

Results for “gamma ray theory”

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 OSTI2024

Optimal celestial bodies for dark matter detection

A wide variety of celestial bodies have been considered as dark matter detectors. Which stands the best chance of delivering the discovery of dark matter? Which is the most powerful dark matter detector? We investigate a range of objects, including the Sun, Earth, Jupiter, Brown Dwarfs, White Dwarfs, Neutron Stars, Stellar populations, and Exoplanets. We quantify how different objects are optimal dark matter detectors in different regimes by deconstructing some of the in-built assumptions in these search sensitivities, including observation potential and particle model assumptions. We find new constraints and future sensitivities across a range of dark matter annihilation final states. We quantify mediator properties leading to detectable celestial-body energy injection or Standard Model fluxes, and show how different objects can be expected to deliver corroborating signals. We discuss different search strategies, their opportunities and limitations, and the interplay of regimes where different celestial objects are optimal dark matter detectors. Deconstructing the assumptions of these searches leads us to point out a new search using the Galactic center stellar population that can provide greater sensitivity to the dark matter-nucleon scattering cross section than the Sun, despite being significantly further away in our Galaxy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Searching for superheavy decaying particles with ultra-high-energy neutrino observatories

If there exist unstable but long-lived relics of the early universe, their decays could produce detectable fluxes of gamma rays and neutrinos. In this paper, we point out that the decays of superheavy particles, m χ ≳ 10 10 GeV, would produce an enhanced flux of ultra-high-energy neutrinos through the processes of muon and pion pair production in the resulting electromagnetic cascades. These processes transfer energy from electromagnetic decay products into neutrinos, relaxing the constraints that can be derived from gamma-ray observations, and increasing the sensitivity of high-energy neutrino telescopes to superheavy particle decays. Taking this into account, we derive new constraints on long-lived superheavy relics from the IceCube Neutrino Observatory, and from the Fermi Gamma-Ray Space Telescope. We find that IceCube-Gen2, and other next generation neutrino telescopes, will provide unprecedented sensitivity to the decays of superheavy dark matter particles and other long-lived relics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Neutrino and gamma-ray signatures of inelastic dark matter annihilating outside neutron stars

We present a new inelastic dark matter search: neutron stars in dark matter-rich environments capture inelastic dark matter which, for interstate mass splittings between about 45 - 285 MeV, will annihilate away before becoming fully trapped inside the object. This means a sizable fraction of the dark matter particles can annihilate while being outsidethe neutron star, producing neutron star-focused gamma-rays and neutrinos. We analyze this effect for the first time and target the neutron star population in the Galactic Center, where the large dark matter and neutron star content makes this signal most significant. Depending on the assumed neutron star and dark matter distributions, we set constraints on the dark matter-nucleon inelastic cross-section using existing H.E.S.S. observations. We also forecast the sensitivity of upcoming gamma-ray and neutrino telescopes to this signal, which can reach inelastic cross-sections as low as ∼ 2 × 10 -47 cm 2 .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Neutrino and pair creation in reconnection-powered coronae of accreting black holes

A ubiquitous feature of accreting black hole systems is their hard X-ray emission which is thought to be produced through Comptonization of soft photons by electrons and positrons in the vicinity of the black hole, in a region with optical depth of order unity. The origin and composition of this Comptonizing region, known as the corona, is a matter open for debate. In this paper we investigate the role of relativistic protons accelerated in black-hole magnetospheric current sheets for the pair enrichment and neutrino emission of AGN coronae. Our model has two free parameters, namely the proton plasma magnetization σ p , which controls the peak energy of the neutrino spectrum, and the Eddington ratio λ X,Edd (defined as the ratio between X-ray luminosity L X and Eddington luminosity L Edd ), which controls the amount of energy transferred to secondary particles. For sources with λ X,Edd ≳ λ Edd,crit (where λ Edd,crit ∼ 10 -1 for σ p = 10 5 or ∼ 10 -2 for σ p = 10 7 ), proton-photon interactions and γγ annihilation produce enough secondary pairs to achieve Thomson optical depths τ T ∼ 0.1-10. In the opposite case of λ X,Edd ≲ λ Edd,crit , the coronal pairs cannot originate only from hadronic interactions. Additionally, we find that the neutrino luminosity scales as L 2 X /L Edd for λ X,Edd ≲ λ Edd,crit , while it is proportional to L X for higher λ X,Edd values. We apply our model to four Seyfert galaxies, including NGC 1068, and discuss our results in light of recent IceCube observations.

active galactic nuclei↗
DOE OSTI2024

Extended Interacting Boson Model with Configuration Mixing Description For 98 Mo

An extended version of the interacting boson model (IBM), encompassing two-particle–two-hole (2p-2h) excitations originating from below a valence shell and one-particle–one-hole (1p-1h) excitations extending to the upper next-nearest major shell, is developed to investigate cross-shell excitation effects in the transitional nucleus 98 Mo. This extended IBM framework is utilized to reproduce experimentally observed properties such as low-lying positive-parity level energies, 𝐵⁡(𝐸⁢2) values, and the electric quadrupole moment of the 2$^{+}_{1}$ state, 𝜌 2 ⁡(𝐸⁢0, 0$^{+}_{2}$ → 0$^{+}_{g}$), as well as the strength functions of the isoscalar giant monopole resonance (ISGMR) and the isoscalar giant quadrupole resonance (ISGQR) in 98 Mo. It is shown that the extended consistent-𝑄 IBM, incorporating configuration mixing of both 2p-2h excitations from below the valence shell and 1p-1h excitations to the upper next-nearest major shell, describes the low-lying structural properties and the ISGMR and ISGQR strength distributions of 98 Mo quite well, with the ISGMR and the ISGQR strength distributions comparable to the results of the self-consistent quasiparticle random-phase approximation with Skyrme interactions.

90 ≤ A ≤ 149↗
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 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

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 OSTI2024

Radioisotope Identification with List-Mode Gamma Ray Data: A rigorous assessment on the value of temporal information applied to radioisotope identification.

This work explores the potential of utilizing temporal data from gamma-ray detectors, known as list-mode data, to enhance radioisotope identification. Traditional identification methods, which rely on full gamma-ray spectrum analysis, often require long dwell times and struggle with “confuser” sources, or spectra with similarly spaced spectral peaks. We hypothesize that by leveraging the probabilistic nature of nuclear decay and the time-encoded information from decay sequences and interactions with surrounding materials, we can improve classification accuracy over static spectral analysis. This research rigorously examines the temporal content of list-mode data through exploratory data analysis via correlation discovery and information theory. We further propose a basic classification model that can utilize spectral or temporal data (or both) to determine if the incorporation of temporal information can improve radioisotope identification. The findings suggest that the temporal information present in list-mode gamma-ray data has merit and should be further investigated.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
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

Updating the MCNP neutron guide extension for supermirror absorption, gamma production, and neutron polarization

As supermirror technology has developed, there is a trend towards higher m-values for supermirrors which propagate more neutrons with higher divergence to experimental areas. These supermirrors have thousands of alternating layers of typically nickel and titanium with varying thicknesses deposited on the surface. Neutrons of grazing incidence are reflected from these layers with a probability dependent on the angle of incidence and the wavelength of the neutron. Those neutrons that are not reflected may be captured in the supermirror layers and produce high energy gamma rays that require shielding. The supermirror extension for MCNP has been updated to estimate the capture of neutrons in the supermirror layers of a neutron guide through the banking of photons in a virtual supermirror layer and the reduction of the weight of the neutrons transmitted into the substrate. A new functionality has been added to simulate neutron polarization specifically in the context of a supermirror polarizer for neutron beamlines. In conclusion, the underlying theory, code implementation in MCNP, and examples of gamma production and a V-cavity polarizer will be discussed.

Gamma production↗
DOE OSTI2025

A Continuous Galactic Line Source of Axions: The Remarkable Case of 23 Na

Here, we argue that $^{23}$Na is a potentially significant source of galactic axions. For temperatures $\gtrsim 7 \times 10^8$K -- characteristic of carbon burning in the massive progenitors of supernovae and ONeMg white dwarfs -- the 440 keV first excited state of $^{23}$Na is thermally populated, with its repeated decays pumping stellar energy into escaping axions. Odd-A nuclear abundances are typically very low in high-temperature stellar environments (or absent entirely due to burn-up). $^{23}$Na is an exception: $\approx 0.1 M_\odot$ of the isotope is synthesized during carbon burning then maintained at $\approx 10^9$K for times ranging up to $6 \times 10^4$y. Using MESA simulations, a galactic model, and sampling over progenitor masses, locations, and evolutionary stages, we find a continuous flux at earth of $\langle ϕ_a \rangle \approx 22$/cm$^2$s for $g^\mathrm{eff}_{aNN} = 10^{-9}$. Some fraction of these axions convert to photons as they propagate through the galactic magnetic field, producing a distinctive 440 keV line $γ$ ray detectable by all-sky detectors like the Compton Spectrometer and Imager (COSI). Assuming a 1$μ$G galactic magnetic field and a sufficiently light axion mass, we find that COSI will be able to probe $| g_{aNN}^\mathrm{eff} g_{a γγ} | \gtrsim1.8 \times 10^{-22}$ GeV$^{-1}$ at $3σ$ after two years of surveying.

Dark matter↗
DOE OSTI2025

Magnetosphere Evolution and Precursor-driven Electromagnetic Signals in Merging Binary Neutron Stars

We detail new force-free simulations to investigate magnetosphere evolution and precursor electromagnetic (EM) signals from binary neutron stars. Our simulations fully follow a representative inspiral motion, capturing the intricate magnetospheric dynamics and their impact on EM outflows. We explore a range of stellar magnetic moment orientations and relative strengths, finding that the magnetospheres and Poynting flux evolution are strongly configuration dependent. The Poynting flux exhibits pulsations at twice the orbital frequency, 2Ω, and is highly anisotropic, following a power-law dependence on orbital frequency. The index ranges from 1 to 6, shaped by the intricate dynamics of the magnetospheres. Furthermore, we present the first computation of (1) the EM forces acting on the star surfaces, revealing the presence of torques that, for highly magnetized stars, could influence the orbital dynamics or break the crust; (2) the high-energy emission signals from these systems by adopting the established isolated pulsar theory. Assuming curvature radiation in the radiation-reaction limit, we find that photons could reach TeV–PeV energies in the last ∼ms for magnetic field strengths 10 10 –10 15 G. However, our analysis of single photon magnetic pair production suggests that these photons are unlikely to escape, with the MeV band emerging as a promising observational window for precursor high-energy emission. In this framework, we construct high-energy emission skymaps and light curves, exploring observational implications. Finally, we propose potential precursor radio emission and delayed afterglow echoes from magnetized outflows, which may contribute to late-time rebrightening in short gamma-ray bursts or to orphan afterglows.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

New Tests of Low-Scale Quantum Gravity with Cosmic-Ray Collisions

Cosmic ray collisions at high center of mass energy could enable graviton and black hole production as expected in theories of low-scale quantum gravity, such as extra-dimensions, many species, or some versions of string theory. Here we propose three novel phenomenological tests of these theories. We first consider the collision of cosmic rays with ambient protons, electrons and photons in Active Galactic Nuclei (AGN), finding that high-energy neutrino data from the blazar TXS 0506+056 places a constraint on the fundamental scale of gravity of $M_f \gtrsim 0.3$ TeV, and future high-energy neutrino data could raise this bound to $M_f \gtrsim 200$ TeV. We then point out that collisions of pairs of cosmic rays could occur at a sizable rate in AGN where the accelerated cosmic rays are not collimated, or on supermassive black hole binaries. This consideration could potentially let us test unprecedented large fundamental scales of $M_f \gtrsim 2$ PeV. We further compute the corresponding thermal neutrino emission arising from the Hawking evaporation of black holes produced in cosmic ray collisions, finding a spectrum that clearly differs from that expected in meson decays. Finally, we speculate with an scenario which would produce high-energy neutrino and gamma-ray emission from regions in the sky where no multi-wavelength counterparts would be expected, via graviton propagation from a different brane, which then decays in our Universe.

Ettengruber, Manuel [IPhT, Saclay]↗
DOE OSTI2025

VARI3D & PERSENT: Perturbation and Sensitivity Analysis

The nodal diffusion method is one of the most widely used approaches in modern reactor analysis. In the nodal diffusion method, a coarse multi-group set of “homogenized” parameters is constructed such that the complex geometry of a reactor core along with the energy dependence of neutron and gamma ray cross sections in a nuclear reactor are conserved in the simpler geometry. The homogenization is typically done on a fuel assembly level as is the case in the DIF3D code developed at Argonne National Laboratory. The nodal methodology is used primarily to predict fuel cycle behavior of nuclear systems of which there is a substantial amount of validation in the literature. Another use of the nodal method is to obtain reactivity coefficients and kinetics parameters for use in a safety analysis of a given nuclear reactor. While there are many ways to obtain reactivity worth and kinetics parameters, the work presented in this manuscript is unique as it provides the user with the ability to compute reactivity worths, kinetics parameters, and cross section sensitivities with a Cartesian and hexagonal geometry based transport code. This manuscript serves as a single manual for two separate codes: VARI3D and PERSENT. The VARI3D code (VARIational 3D) is based upon the classic finite difference diffusion theory solver available in DIF3D. The PERSENT code (PERturbation and SENitivity for Transport) is based upon the variational nodal method employed in DIF3D termed VARIANT. The VARIANT solver was added to DIF3D in 1995 and has seen continued development and use for the last 18 years. Because VARI3D primarily uses deprecated coding practices, rather than incorporating the perturbation and sensitivity treatments for transport within VARI3D, a new coding development was built using modern Fortran coding. The primary purpose of this manual is to describe the theory behind PERSENT (and by convenience, that of VARI3D) and discuss the input and output of PERSENT along with giving potential users an idea of how to use it. While this manuscript does describe the input and output of VARI3D, the PERSENT code is intended to be the replacement capability of VARI3D as PERSENT can generate nearly identical (if not superior) diffusion theory results. In this manuscript, the relevant aspects of generalized perturbation theory and exact perturbation theory that apply to both VARI3D and PERSENT are covered. The input and output of VARI3D is displayed by excerpting several of the example problems. Similarly, the input and output of PERSENT is displayed along with tips on how best to use the code. Note that the input and output of the inhomogeneous solver wrapped around DIF3D (DIF3D_IFS) is also discussed as it is needed to carry out some of the sensitivities in PERSENT such as reaction rate ratios. This manuscript describes several perturbation and sensitivity problems, and the results computed using PERSENT. From these sections, potential users should find that PERSENT provides not only the typical tables of numbers desired in perturbation and sensitivity analysis work, but also can visually plot the result for a more thorough understanding of the space and energy distribution (Section 5). Overall, PERSENT is observed to produce accurate reactivity worths and sensitivities for the displayed set of test problems and clearly demonstrates the need to have a transport-based sensitivity capability as evident from the thousands of percent errors observed in the 21-group hexagonal fast reactor problem (covered in Section 7). The uncertainty calculation capability is described in Section 3 and demonstrated in Section 7.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Measurement of the Isolated Nuclear Two-Photon Decay in Ge 72

The nuclear two-photon or double-gamma (2 γ ) decay is a second-order electromagnetic process whereby a nucleus in an excited state emits two gamma rays simultaneously. To be able to directly measure the 2⁢ γ decay rate in the low-energy regime below the electron-positron pair-creation threshold, we combined the isochronous mode of a storage ring with Schottky resonant cavities. The newly developed technique can be applied to isomers with excitation energies down to ~100 keV and half-lives as short as ~10 ms. The half-life for the 2⁢ γ decay of the first-excited 0 + state in bare 72 Ge ions was determined to be 23.9(6) ms, which strongly deviates from expectations.

59 ≤ A ≤ 89↗
DOE OSTI2026

Burn mode signatures enabled by high dynamic range fusion reaction history

The evolution of the fusion burn of a compressed inertial confinement fusion (ICF) implosion gives information about the evolution of the temperature, mass, and volume of the hot spot. Currently, the fusion reaction history has been measured with about a decade of dynamic range, giving information about just the peak of the fusion burn. There are proposals for extending the dynamic range to 1000×, measuring the rising edge of the burn earlier in time. Using fusion hot spot theory and a set of ICF simulations, we identify and categorize different stages of fusion burn and what signatures could be measured. For ice layered implosions, the details, conditions, and dynamics of the fusion burn propagation could be observed. For double shell and OMEGA scale implosions, the evolution of the final compression and rebounding shock before ignition could be observed. Higher dynamic range expands investigation into the logarithmic derivative, α, as a signature of various mechanisms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Programmable Phase Selection between Altermagnetic and Noncentrosymmetric Polymorphs of MnTe on InP via Molecular Beam Epitaxy

This dataset contains DFT input and output files supporting the theoretical modeling in the associated publication (ACS Appl. Mater. Interfaces 2026, 18, 15654-15664). The calculations model the interfacial energetics of two MnTe polymorphs — NiAs-MnTe (hexagonal, alpha phase) and ZnS-MnTe (cubic, gamma phase) — on InP(111) substrates with two surface terminations: In-terminated InP(111)A and P-terminated InP(111)B. This gives four interface configurations: NiAs on In-terminated (experimentally observed), NiAs on P-terminated (computed for comparison), ZnS on In-terminated (computed for comparison), and ZnS on P-terminated (experimentally observed). The dataset is organized into four calculation types, each covering all four polymorph/termination combinations: (i) Slabs: Pristine MnTe/InP heterostructure slabs used to compute total energies and interface energy densities (Eint) for all four configurations, as reported in Fig. 6 of the main text. (ii) Disorder: Same slab geometries with a P_Te + Te_P antisite defect pair introduced near the interface, used to assess chemical intermixing effects on interface stability (Fig. S8, SI). (iii) Strain: Pristine slab calculations with in-plane lattice parameters strained by -1% and +1% relative to the InP lattice constant, used to evaluate strain-dependent interface energetics (Fig. S9, SI). (iv) Charge_Density: Single-point calculations on the full heterostructure, the isolated InP slab, and the isolated MnTe slab at fixed geometry, used to compute differential charge density plots showing interfacial charge accumulation and depletion as a function of surface termination (Fig. S10, SI). Each calculation folder contains INCAR, KPOINTS, POSCAR, CONTCAR, OUTCAR, and POTCAR_info.txt (PAW potential information, excluding the full POTCAR due to VASP licensing restrictions). The calculations were performed using VASP 6.4.3 with PBE exchange-correlation, PAW potentials, a Hubbard correction of Ueff = 5 eV on Mn d-states, and A-type AFM spin initialization.

36 MATERIALS SCIENCE↗