Modeling Aerial Gamma-Ray Backgrounds using Non-negative Matrix Factorization

Airborne gamma-ray surveys are useful for many applications, ranging from geology and mining to public health and nuclear security. In all these contexts, the ability to decompose a measured spectrum into a linear combination of background source terms can provide useful insights into the data and lead to improvements over techniques that use spectral energy windows. Multiple methods for the linear decomposition of spectra exist but are subject to various drawbacks, such as allowing negative photon fluxes or requiring detailed Monte Carlo modeling. We propose using Non-negative Matrix Factorization (NMF) as a data-driven approach to spectral decomposition. Using aerial surveys that include flights over water, we demonstrate that the mathematical approach of NMF finds physically relevant structure in aerial gamma-ray background, namely that measured spectra can be expressed as the sum of nearby terrestrial emission, distant terrestrial emission, and radon and cosmic emission. These NMF background components are compared to the background components obtained using Noise-Adjusted Singular Value Decomposition (NASVD), which contain negative photon fluxes and thus do not represent emission spectra in as straightforward a way. Finally, we comment on potential areas of research that are enabled by NMF decompositions, such as new approaches to spectral anomaly detection and data fusion.Comment: 14 pages, 12 figures, accepted for publication in IEEE Transactions on Nuclear Scienc

Geant4-related R&D for new particle transport methods

A R&D project has been launched in 2009 to address fundamental methods in radiation transport simulation and revisit Geant4 kernel design to cope with new experimental requirements. The project focuses on simulation at different scales in the same experimental environment: this set of problems requires new methods across the current boundaries of condensed-random-walk and discrete transport schemes. An exploration is also foreseen about exploiting and extending already existing Geant4 features to apply Monte Carlo and deterministic transport methods in the same simulation environment. An overview of this new R&D associated with Geant4 is presented, together with the first developments in progress.Comment: 4 pages, to appear in proceedings of the Nuclear Science Symposium and Medical Imaging Conference 2009, Orland

Component characterization and commissioning of a gamma-PET prototype detector system

Hybrid imaging systems, comprising PET and Compton camera modules, have recently gained in interest, due to their capability to simultaneously detect positron annihilation photons and γ-rays from single-photon emitting sources as also used in SPECT. A unique feature of such systems, however, is the capability to also be operated in a so called γ-PET mode. Here, specific β+- emitting radioisotopes (such as 44Sc, 1°C or 14O) are used to detect triple-coincidences between two annihilation γ-rays (in PET imaging) and a third, prompt photon (in Compton imaging), that is emitted by the deexcitation of the decay’s daughter nucleus. Consequently, an intersection between the line-of-response (LOR) and the Compton cone can be determined, which (in principle) allows to localize the photons’ emission vertices on a single decay basis. In practice, however, a few tens of events are required to localize a point source, which still results in a considerable sensitivity improvement compared to conventional PET imaging. For a proof-of-principle study, we used a pixelated GAGG crystal array (16 × 16 crystals; 1.45 × 1.45 × 6 mm3 crystal volume; 25 μm SPAD SiPMs as readout) as Compton camera scatterer and PET detectors, and a three-layered LYSO crystal array (1.2 × 1.2 × 6.66 mm3 crystal volume; 50 μm SPAD SiPMs as readout) as Compton camera absorber. We characterized the individual detector components with regard to their energy resolution and the capability to identify the various scintillator array’s individual crystals. Our first γ-PET prototype was tested in PET-only and Compton-only imaging mode, in which spatial resolutions of 3.2–3.5 mm FWHM (PET-only mode) and 14.4–19.3 mm FWHM (Compton-only mode at 1,274 keV) were achieved, respectively, using a22Na point source and 10 iterations of an ML-EM reconstruction algorithm. By using triple-coincidences in a γ-PET mode (event-wise intersection of the LOR and the Compton cone), we could demonstrate the capability of the prototype to perform a full 3D point source reconstruction using only 77 events

Positron Annihilation in the Galaxy

The 511 keV line from positron annihilation in the Galaxy was the first γ-ray line detected to originate from outside our solar system. Going into the fifth decade since the discovery, the source of positrons is still unconfirmed and remains one of the enduring mysteries in γ-ray astronomy. With a large flux of ∼10−3 γ/cm2/s, after 15 years in operation INTEGRAL/SPI has detected the 511 keV line at >50σ and has performed high-resolution spectral studies which conclude that Galactic positrons predominantly annihilate at low energies in warm phases of the interstellar medium. The results from imaging are less certain, but show a spatial distribution with a strong concentration in the center of the Galaxy. The observed emission from the Galactic disk has low surface brightness and the scale height is poorly constrained, therefore, the shear number of annihilating positrons in our Galaxy is still not well know. Positrons produced in β+-decay of nucleosynthesis products, such as 26Al, can account for some of the annihilation emission in the disk, but the observed spatial distribution, in particular the excess in the Galactic bulge, remains difficult to explain. Additionally, one of the largest uncertainties in these studies is the unknown distance that positrons propagate before annihilation. In this paper, we will summarize the current knowledge base of Galactic positrons, and discuss how next-generation instruments could finally provide the answers.Non peer reviewedFinal Accepted Versio

Locating the most energetic electrons in Cassiopeia A

We present deep ($>$2.4 Ms) observations of the Cassiopeia A supernova remnant with {\it NuSTAR}, which operates in the 3--79 keV bandpass and is the first instrument capable of spatially resolving the remnant above 15 keV. We find that the emission is not entirely dominated by the forward shock nor by a smooth "bright ring" at the reverse shock. Instead we find that the $>$15 keV emission is dominated by knots near the center of the remnant and dimmer filaments near the remnant's outer rim. These regions are fit with unbroken power-laws in the 15--50 keV bandpass, though the central knots have a steeper ($\Gamma \sim -3.35$) spectrum than the outer filaments ($\Gamma \sim -3.06$). We argue this difference implies that the central knots are located in the 3-D interior of the remnant rather than at the outer rim of the remnant and seen in the center due to projection effects. The morphology of $>$15 keV emission does not follow that of the radio emission nor that of the low energy ($<$12 keV) X-rays, leaving the origin of the $>$15 keV emission as an open mystery. Even at the forward shock front we find less steepening of the spectrum than expected from an exponentially cut off electron distribution with a single cutoff energy. Finally, we find that the GeV emission is not associated with the bright features in the {\it NuSTAR} band while the TeV emission may be, suggesting that both hadronic and leptonic emission mechanisms may be at work.Comment: 12 pages, 11 figures, accepted for publication in Ap

GRIPS - Gamma-Ray Imaging, Polarimetry and Spectroscopy

We propose to perform a continuously scanning all-sky survey from 200 keV to 80 MeV achieving a sensitivity which is better by a factor of 40 or more compared to the previous missions in this energy range. The Gamma-Ray Imaging, Polarimetry and Spectroscopy (GRIPS) mission addresses fundamental questions in ESA's Cosmic Vision plan. Among the major themes of the strategic plan, GRIPS has its focus on the evolving, violent Universe, exploring a unique energy window. We propose to investigate $\gamma$-ray bursts and blazars, the mechanisms behind supernova explosions, nucleosynthesis and spallation, the enigmatic origin of positrons in our Galaxy, and the nature of radiation processes and particle acceleration in extreme cosmic sources including pulsars and magnetars. The natural energy scale for these non-thermal processes is of the order of MeV. Although they can be partially and indirectly studied using other methods, only the proposed GRIPS measurements will provide direct access to their primary photons. GRIPS will be a driver for the study of transient sources in the era of neutrino and gravitational wave observatories such as IceCUBE and LISA, establishing a new type of diagnostics in relativistic and nuclear astrophysics. This will support extrapolations to investigate star formation, galaxy evolution, and black hole formation at high redshifts.Comment: to appear in Exp. Astron., special vol. on M3-Call of ESA's Cosmic Vision 2010; 25 p., 25 figs; see also www.grips-mission.e

An exceptionally bright flare from SGR1806-20 and the origins of short-duration gamma-ray bursts

Soft-gamma-ray repeaters (SGRs) are galactic X-ray stars that emit numerous short-duration (about 0.1 s) bursts of hard X-rays during sporadic active periods. They are thought to be magnetars: strongly magnetized neutron stars with emissions powered by the dissipation of magnetic energy. Here we report the detection of a long (380 s) giant flare from SGR 1806-20, which was much more luminous than any previous transient event observed in our Galaxy. (In the first 0.2 s, the flare released as much energy as the Sun radiates in a quarter of a million years.) Its power can be explained by a catastrophic instability involving global crust failure and magnetic reconnection on a magnetar, with possible large-scale untwisting of magnetic field lines outside the star. From a great distance this event would appear to be a short-duration, hard-spectrum cosmic gamma-ray burst. At least a significant fraction of the mysterious short-duration gamma-ray bursts therefore may come from extragalactic magnetars.Comment: 21 pages, 5 figures. Published in Natur

NuSTAR observations of the bullet cluster: constraints on inverse compton emission

The search for diffuse non-thermal inverse Compton (IC) emission from galaxy clusters at hard X-ray energies has been undertaken with many instruments, with most detections being either of low significance or controversial. Because all prior telescopes sensitive at E &gt; 10 keV do not focus light and have degree-scale fields of view, their backgrounds are both high and difficult to characterize. The associated uncertainties result in lower sensitivity to IC emission and a greater chance of false detection. In this work, we present 266 ks NuSTAR observations of the Bullet cluster, which is detected in the energy range 3-30 keV. NuSTAR's unprecedented hard X-ray focusing capability largely eliminates confusion between diffuse IC and point sources; however, at the highest energies, the background still dominates and must be well understood. To this end, we have developed a complete background model constructed of physically inspired components constrained by extragalactic survey field observations, the specific parameters of which are derived locally from data in non-source regions of target observations. Applying the background model to the Bullet cluster data, we find that the spectrum is well-but not perfectly-described as an isothermal plasma with kT = 14.2 ± 0.2 keV. To slightly improve the fit, a second temperature component is added, which appears to account for lower temperature emission from the cool core, pushing the primary component to kT ~ 15.3 keV. We see no convincing need to invoke an IC component to describe the spectrum of the Bullet cluster, and instead argue that it is dominated at all energies by emission from purely thermal gas. The conservatively derived 90% upper limit on the IC flux of 1.1 × 10-12 erg s-1 cm-2 (50-100 keV), implying a lower limit on B ≳ 0.2 μG, is barely consistent with detected fluxes previously reported. In addition to discussing the possible origin of this discrepancy, we remark on the potential implications of this analysis for the prospects for detecting IC in galaxy clusters in the future