GaMMA: Gamma-Ray Modeling & Analysis

• GaMMA is a comprehensive framework that integrates experimental, analytical, and computational approaches to generate, detect, and analyze gamma-ray and spectral data.
• Its supervised NMF method in the gamma_flow package enables real-time spectral decomposition with high classification accuracy (>90%) and reconstruction fidelity (>98%).
• The framework supports a range of applications from high-energy astrophysics and collider physics to spectroscopy and materials science through modular, interpretable pipelines.

GaMMA refers to a set of frameworks, instruments, and computational methods that address the generation, detection, and analysis of gamma-ray or spectral data across multiple disciplines, ranging from high-energy particle physics and astrophysics to analytical spectroscopy and computational modeling. The various Gamma- or GaMMA-derived projects are distinguished by their specialized aims—either experimental (gamma-beam facilities, telescopes, colliders), analytical (matrix-based spectral analysis), or computational (hyperrelativistic fluid simulations). Each instantiation leverages the unique physical properties of gamma photons or the distinct characteristics of one-dimensional spectral data to enable scientific inference or technological application.

1. Supervised Factorization for Spectral Analysis: GaMMA (gamma_flow)

The GaMMA framework, as implemented in the gamma_flow open-source package, provides real-time automated analysis of one-dimensional spectral data—originally for gamma-ray spectra, though the methodology generalizes to other spectroscopic domains (IR, Raman, mass spectrometry, UV-Vis, stellar spectra). The core is a supervised non-negative matrix factorization (NMF) model for dimensionality reduction and interpretable decomposition:

• Formalism: Let $\mathbf{X}\in\mathbb{R}^{n\times m}_{\geq 0}$ represent $n$ spectra (each an $m$-dimensional channel count vector), and let $\mathbf{L}\in\mathbb{R}^{m\times k}_{\geq 0}$ denote the fixed basis whose $k$ columns are the mean spectra for each known training label. The score matrix $\mathbf{S}\in\mathbb{R}^{n\times k}_{\geq 0}$ encodes non-negative contributions of each basis spectrum to each observation, estimated via a row-wise non-negative least-squares (NNLS) optimization:

$$\min_{\mathbf{S}\geq 0}\ \frac{1}{2} \|\mathbf{X} - \mathbf{S}\mathbf{L}^\mathsf{T}\|^2_F + \frac{\lambda}{2}\|\mathbf{S}\|^2_F$$

• Workflow: The typical pipeline is:
  1. Preprocessing: Energy calibration, rebinning to a common channel grid, data aggregation.
  2. Decomposition: Compute $\mathbf{L}$ by averaging spectra per label in the training set. For each test spectrum $\mathbf{x}_i$, solve NNLS to infer $\mathbf{s}_i$.
  3. Denoising: Project input onto physical basis: $$\widehat{\mathbf{x}}_i = \mathbf{s}_i \mathbf{L}^\mathsf{T}$$.
  4. Classification: Assign label by highest normalized score or multi-label rule with threshold $\tau$.
  5. Outlier Detection: Use cosine similarity between input and reconstruction, $\rho_i$, with threshold $\eta$.

Table: Key Performance Metrics for GaMMA (gamma_flow) (Rädle et al., 12 Nov 2025)

Metric	Value / Method	Note
Classification accuracy	>90% (single-label), F1 > 0.92 (multi-label)	Demonstrated on test set of 5 isotopes+background
Denoising quality	$\rho > 0.98$, explained variance >98%	High-fidelity reconstruction
Outlier detection	Precision/Recall defined via standard formula	Cosine similarity decision boundary
Inference speed	$\sim$1 ms/spectrum (CPU), $m=10^2-10^3$	No GPU/massive memory required
Generalizability	Any 1D spectra, not just gamma-ray	Requires label-averaged basis per application

Applying supervised NMF with a fixed, physical basis yields interpretable latent axes (each corresponding to a label such as an isotope), offering a transparent alternative to black-box neural network models while maintaining operational efficiency and adaptability (Rädle et al., 12 Nov 2025).

2. Gamma-Ray Astrophysics Missions and Telescopes

Several missions and experimental platforms labeled as “GAMMA” or “GAMMA-LIGHT/400” are designed to probe fundamental questions in astrophysics, dark matter, and cosmic-ray propagation:

• GAMMA-LIGHT covers the 10 MeV–100 GeV energy range, bridging the observational gap left by prior instruments (COMPTEL, AGILE, Fermi-LAT). It features:
  • A high-resolution silicon tracker (41 trays of microstrip detectors) for sub-degree point-spread function ($\theta_{68\%} \sim 0.1^\circ$ at 1 GeV).
  • A CsI calorimeter and anticoincidence systems for full gamma-ray event reconstruction.
  • Energy resolution parameterized as $$\Delta E/E \sim 0.15 (E/100\,\rm MeV)^{-0.3} \oplus 0.05$$; effective area up to 2000 cm² at 1 GeV.
  • Sensitivities enable decisive studies of diffuse galactic emission, SNR pion decay, dark matter annihilation/decay, and transients like GRBs (Morselli et al., 2014).
• GAMMA-400 is a next-generation telescope with lateral and top-down detection, unique for its ability to register gamma-ray bursts (GRBs) from lateral directions using a 16 $X_0$ CsI(Tl) calorimeter (CC2) and lateral detectors. Salient metrics:
  • Lateral effective area $\sim0.13$ m² per side ($\sim6$ steradian total FoV).
  • Energy resolution $\sim$10–15% (10–100 MeV), $\sim$2% (100 GeV).
  • Angular resolution $<1^\circ$ (100 MeV), improving to $<0.01^\circ$ (100 GeV).
  • Simulations predict up to 320 GRB detections/yr (lateral mode), facilitating high-statistics prompt and afterglow studies in the 10–100 MeV regime (Leonov et al., 2021).

Such missions generate critical high-energy datasets necessary for resolving the physical origin of cosmic gamma-ray emission and enabling cross-correlation with other spectral bands.

3. Gamma-Gamma and Gamma Factories: Photon Colliders and High-Flux Sources

“Gamma factories” and “gamma-gamma colliders” designate accelerator-based sources for the production of high-intensity, energy-tunable, quasi-monochromatic gamma-ray beams, often with the capability of producing secondary polarized particle beams:

• Gamma Factory (CERN): Utilizes resonant laser excitation of partially stripped ions (PSI) in the LHC to generate gamma rays with photon energies in the 1–400 MeV range. The upscattering process:

$$E_\gamma(\theta) = \frac{4\gamma_L^2 h\nu_L}{1 + a_0^2 + \gamma_L^2\theta^2}$$

where $\gamma_L$ is the ion Lorentz factor, $h\nu_L$ is incident photon energy, and $a_0$ the normalized laser strength. With $N_\mathrm{ions} \sim 10^9$ and $P_\mathrm{RF} \sim 100$ kW, photon fluxes up to $10^{17}$ $\gamma$/s are anticipated. Polished secondary beams—polarized $e^+$, $\mu^+$, cold neutrons—are key deliverables for future collider and neutrino-factory concepts.

• SAPPHiRE γγ Higgs Factory: Proposes a pair of 10 GeV recirculating linacs to create 80 GeV electrons, which collide with high-power lasers near the IP, Compton upscattering photons to achieve $E_{CM}(\gamma\gamma) \sim 125$ GeV suitable for resonant Higgs production. The collider achieves:
  • Peak luminosity $$L_{\gamma\gamma} \sim 0.36 \times 10^{34}\,\mathrm{cm}^{-2}\mathrm{s}^{-1}$$
  • $\sim$20,000 Higgs events/year (integrated luminosity $\sim$120 fb$^{-1}$)
  • Energy spectrum width $\Delta E_{\gamma\gamma}/E_{\gamma\gamma} \sim$ 10%
  • Statistical precision: $\sim$2% ($h\to b\bar b$), $\sim$5% ($h\to WW^*$), $\sim$8% ($h\to \gamma\gamma$), 100 MeV mass scan (Bogacz et al., 2012, Krasny, 2015)

These facilities exploit either Compton-backscattering (electron+laser) or resonant scattering from PSI, offering unprecedented beam intensity and tunability for fundamental physics and applied research.

4. DAΦNE–GAΜMΑ: Storage-Ring Compton Gamma-Ray Source

The DAΦNE–GAΜMΑ project defines a storage-ring–based gamma factory using Compton backscattering of laser photons off a high-current ($I=1.5$ A), low-emittance ($\epsilon_x = 0.10$ mm·mrad) electron ring and a %%%%50$$\widehat{\mathbf{x}}_i = \mathbf{s}_i \mathbf{L}^\mathsf{T}$$51%%%% finesse Fabry–Pérot cavity ($\sim$37 kW stored power):

• Photon Energy: Tunable via electron energy (250–900 MeV) and laser wavelength (0.5–10 μm), yielding $E_{\gamma,\max}$ in 2–9 MeV range.
• Key Source Metrics (benchmark case):
  • $N_\gamma \sim 10^{12}$ ph/s, $\Delta E/E \sim 0.5\%$, spectral density $>5\times 10^4$ ph/s/eV.
  • Minimal perturbation to e-beam (energy spread/emittance) due to separation of collision and damping timescales.
  • Rapid energy tuning and compact layouts facilitate deployment in nuclear physics, radiology, and materials science.
• Comparison: Demonstrates flux and bandwidth competitive with or surpassing other sources (Duke HIGS, ELI-NP, Mega-Ray, IRIDE), at MHz repetition and continuous tunability (Alesini et al., 2014).

5. Computational Methods: GAMMA for Relativistic Blastwave Modeling

The GAMMA code represents a modern approach for simulating relativistic hydrodynamics and associated non-thermal emission, particularly in the context of gamma-ray burst afterglows:

• Scheme: Implements ALE (arbitrary Lagrangian–Eulerian) SRHD on a moving mesh—advecting along the dominant fluid direction to avoid mesh entanglement, maximizing local resolution at shocks, and enabling efficient large-scale blast wave evolution.
• Microphysics: Includes in situ shock detection, injection of power-law electron distributions ($N_e'(\gamma)\propto\gamma^{-p}$ for $\gamma_{\min} < \gamma < \gamma_{\max}$), local radiative cooling via synchrotron and inverse-Compton processes, and per-zone broadband synchrotron spectra.
• Astrophysical Impact: Demonstrates that the local treatment of synchrotron cooling results in a critical frequency shift (cooling break) by a factor $\sim40$ above predictions from global, spatially averaged models. The package provides validated, high-performance light-curve synthesis from early relativistic to late Newtonian regimes (Ayache et al., 2021).

6. Scientific and Industrial Applications

The ensemble of GaMMA-derived technologies and methods underpins a broad array of scientific and technological endeavors:

• High-energy astrophysics: Studies of cosmic-ray acceleration, non-thermal processes in supernova remnants, dark matter indirect detection, and gamma-ray burst phenomenology.
• Collider physics: Precision Higgs property measurements, electroweak quartic coupling studies, and searches for new physics via dedicated $\gamma\gamma$ and $e\gamma$ collisions.
• Nuclear and materials science: Exploitation of high-brilliance, quasi-monochromatic gamma-ray sources for isotope production, waste transmutation, spectroscopy, neutron radiography, industrial tomography, and non-destructive materials testing.
• Spectroscopy and analytical chemistry: Real-time, robust analysis of multi-component spectra for research and industrial quality control across diverse modalities.

In all manifestations, GaMMA solutions strive for high performance, physical interpretability, and flexibility—leveraging real-time computation, high photon flux, and/or modular hybrid numerical strategies as appropriate for the application domain.

---

Key References:

• Supervised NMF for spectral analysis: "GAMMA_FLOW: Guided Analysis of Multi-label spectra by MAtrix Factorization for Lightweight Operational Workflows" (Rädle et al., 12 Nov 2025)
• High-energy astrophysics missions: "GAMMA-LIGHT: High-Energy Astrophysics above 10 MeV" (Morselli et al., 2014); "Capabilities of the GAMMA-400 gamma-ray telescope..." (Leonov et al., 2021)
• Gamma factories and colliders: "The Gamma Factory proposal for CERN" (Krasny, 2015), "SAPPHiRE: a Small Gamma-Gamma Higgs Factory" (Bogacz et al., 2012)
• Storage-ring Compton sources: "Daφne gamma-rays factory" (Alesini et al., 2014)
• Computational GRB modeling: "GAMMA: a new method for modeling relativistic hydrodynamics and non-thermal emission on a moving mesh" (Ayache et al., 2021)