Reticulum II (Ret II) is a satellite galaxy of the Milky Way and presents a prime target to investigate the nature of dark matter (DM) because of its high mass-to-light ratio. We evaluate a dedicated INTEGRAL observation campaign data set to obtain $
gamma$-ray fluxes from Ret II and compare those with expectations from DM. Ret II is not detected in the $gamma$-ray band 25--8000 keV, and we derive a flux limit of $lesssim 10^{-8},mathrm{erg,cm^{-2},s^{-1}}$. The previously reported 511 keV line is not seen, and we find a flux limit of $lesssim 1.7 times 10^{-4},mathrm{ph,cm^{-2},s^{-1}}$. We construct spectral models for primordial black hole (PBH) evaporation and annihilation/decay of particle DM, and subsequent annihilation of positrons produced in these processes. We exclude that the totality of DM in Ret II is made of a monochromatic distribution of PBHs of masses $lesssim 8 times 10^{15},mathrm{g}$. Our limits on the velocity-averaged DM annihilation cross section into $e^+e^-$ are $langle sigma v rangle lesssim 5 times 10^{-28} left(m_{rm DM} / mathrm{MeV} right)^{2.5},mathrm{cm^3,s^{-1}}$. We conclude that analysing isolated targets in the MeV $gamma$-ray band can set strong bounds on DM properties without multi-year data sets of the entire Milky Way, and encourage follow-up observations of Ret II and other dwarf galaxies.
We use 15 years of $gamma$-ray data from INTEGRAL/SPI in a refined investigation of the morphology of the Galactic bulge positron annihilation signal. Our spatial analysis confirms that the signal traces the old stellar population in the bulge and re
veals for the first time that it traces the boxy bulge and nuclear stellar bulge. Using a 3D smoothing kernel, we find that the signal is smeared out over a characteristic length scale of $150 pm 50,$pc, suggesting either annihilation in situ at astrophysical sources kicked at formation or positron propagation away from sources. The former is disfavoured by its requiring kick velocities different between the Galactic nucleus ($gtrsim 50,mathrm{km,s^{-1}}$) and wider bulge ($lesssim 15,mathrm{km,s^{-1}}$) source. Positron propagation prior to annihilation can explain the overall phenomenology of the 511 keV signal for positrons injection energies $lesssim 1.4,$MeV, suggesting a nucleosynthesis origin.
Classical novae are among the most frequent transient events in the Milky Way, and key agents of ongoing nucleosynthesis. Despite their large numbers, they have never been observed in soft $gamma$-ray emission. Measurements of their $gamma$-ray signa
tures would provide both, insights on explosion mechanism as well as nucleosynthesis products. Our goal is to constrain the ejecta masses of $mathrm{^7Be}$ and $mathrm{^{22}Na}$ from classical novae through their $gamma$-ray line emissions at 478 and 1275 keV. We extract posterior distributions on the line fluxes from archival data of the INTEGRAL/SPI spectrometer telescope. We then use a Bayesian hierarchical model to link individual objects and diffuse emission and infer ejecta masses from the whole population of classical novae in the Galaxy. Individual novae are too dim to be detectable in soft $gamma$-rays, and the upper bounds on their flux and ejecta mass uncertainties cover several orders of magnitude. Within the framework of our hierarchical model, we can, nevertheless, infer tight upper bounds on the $mathrm{^{22}Na}$ ejecta masses, given all uncertainties from individual objects as well as diffuse emission, of $<2.0 times 10^{-7},mathrm{M_{odot}}$ (99.85th percentile). In the context of ONe nucleosynthesis, the $mathrm{^{22}Na}$ bounds are consistent with theoretical expectations, and exclude that most ONe novae happen on white dwarfs with masses around $1.35,mathrm{M_{odot}}$. The upper bounds from $mathrm{^{7}Be}$ are uninformative. From the combined ejecta mass estimate of $mathrm{^{22}Na}$ and its $beta^+$-decay, we infer a positron production rate of $<5.5 times 10^{42},mathrm{e^+,s^{-1}}$, which would make at most 10% of the total annihilation rate in the Milky Way.
The soft MeV gamma-ray sky, from a few hundred keV up to several MeV, is one of the least explored regions of the electromagnetic spectrum. The most promising technology to access this energy range is a telescope that uses Compton scattering to detec
t the gamma rays. Going from the measured data to all-sky images ready for scientific interpretation, however, requires a well-understood detector setup and a multi-step data-analysis pipeline. We have developed these capabilities for the Compton Spectrometer and Imager (COSI). Starting with a deep understanding of the many intricacies of the Compton measurement process and the Compton data space, we developed the tools to perform simulations that match well with instrument calibrations and to reconstruct the gamma-ray path in the detector. Together with our work to create an adequate model of the measured background while in flight, we are able to perform spectral and polarization analysis, and create images of the gamma-ray sky. This will enable future telescopes to achieve a deeper understanding of the astrophysical processes that shape the gamma-ray sky from the sites of star formation (26-Al map), to the history of core-collapse supernovae (e.g. 60-Fe map) and the distributions of positron annihilation (511-keV map) in our Galaxy.
The balloon-borne Compton Spectrometer and Imager (COSI) had a successful 46-day flight in 2016. The instrument is sensitive to photons in the energy range $0.2$-$5$ MeV. Compton telescopes have the advantage of a unique imaging response and provide
the possibility of strong background suppression. With its high-purity germanium detectors, COSI can precisely map $gamma$-ray line emission. The strongest persistent and diffuse $gamma$-ray line signal is the 511 keV emission line from the annihilation of electrons with positrons from the direction of the Galactic centre. While many sources have been proposed to explain the amount of positrons, $dot{N}_{mathrm{e^+}} sim 10^{50},mathrm{e^+,yr^{-1}}$, the true contributions remain unsolved. In this study, we aim at imaging the 511 keV sky with COSI and pursue a full-forward modelling approach, using a simulated and binned imaging response. For the strong instrumental background, we describe an empirical approach to take the balloon environment into account. We perform two alternative methods to describe the signal: Richardson-Lucy deconvolution, an iterative method towards the maximum likelihood solution, and model fitting with pre-defined emission templates. Consistently with both methods, we find a 511 keV bulge signal with a flux between $0.9$ and $3.1 times 10^{-3},mathrm{ph,cm^{-2},s^{-1}}$, confirming earlier measurements, and also indications of more extended emission. The upper limit we find for the 511 keV disk, $< 4.3 times 10^{-3},mathrm{ph,cm^{-2},s^{-1}}$, is consistent with previous detections. For large-scale emission with weak gradients, coded aperture mask instruments suffer from their inability to distinguish isotropic emission from instrumental background, while Compton-telescopes provide a clear imaging response, independent of the true emission.
We illustrate a method for estimating the vertical position of the Sun above the Galactic plane by $gamma$-ray observations. Photons of $gamma$-ray wavelengths are particularly well suited for geometrical and kinematic studies of the Milky Way becaus
e they are not subject to extinction by interstellar gas or dust. Here, we use the radioactive decay line of $mathrm{^{26}Al}$ at $1.809,mathrm{MeV}$ to perform maximum likelihood fits to data from the spectrometer SPI on board the INTEGRAL satellite as a proof-of-concept study. Our simple analytic 3D emissivity models are line-of-sight integrated, and varied as a function of the Suns vertical position, given a known distance to the Galactic centre. We find a vertical position of the Sun of $z_0 = 15 pm 17,mathrm{pc}$ above the Galactic plane, consistent with previous studies, finding $z_0$ in a range between $5$ and $29,mathrm{pc}$. Even though the sensitivity of current MeV instruments is several orders of magnitude below that of telescopes for other wavelengths, this result reveals once more the disregarded capability of soft $gamma$-ray telescopes. We further investigate possible biases in estimating the vertical extent of $gamma$-ray emission if the Suns position is set incorrectly, and find that the larger the true extent, the less is it affected by the observer position. In the case of $mathrm{^{26}Al}$ with an exponential scale height of $150,mathrm{pc}$ ($700,mathrm{pc}$) in the inner (full) Galaxy, this may lead to misestimates of up to $25,%$.
Context. The diffuse gamma-ray emission of $^{26}{rm Al}$ at 1.8 MeV reflects ongoing nucleosynthesis in the Milky Way, and traces massive-star feedback in the interstellar medium due to its 1 Myr radioactive lifetime. Interstellar-medium morphology
and dynamics are investigated in astrophysics through 3D hydrodynamic simulations in fine detail, as only few suitable astronomical probes are available. Aims. We compare a galactic-scale hydrodynamic simulation of the Galaxys interstellar medium, including feedback and nucleosynthesis, with gamma-ray data on $^{26}{rm Al}$ emission in the Milky Way extracting constraints that are only weakly dependent on the particular realisation of the simulation or Galaxy structure. Methods. Due to constraints and biases in both the simulations and the gamma-ray observations, such comparisons are not straightforward. For a direct comparison, we perform maximum likelihood fits of simulated sky maps as well as observation-based maximum entropy maps to measurements with INTEGRAL/SPI. To study general morphological properties, we compare the scale heights of $^{26}{rm Al}$ emission produced by the simulation to INTEGRAL/SPI measurements.} Results. The direct comparison shows that the simulation describes the observed inner Galaxy well, but differs significantly from the observed full-sky emission morphology. Comparing the scale height distribution, we see similarities for small scale height features and a mismatch at larger scale heights. We attribute this to the prominent foreground emission sites that are not captured by the simulation.
The annihilation of cosmic positrons ($e^+$) with electrons in the interstellar medium (ISM) results in the strongest persistent gamma-ray line signal in the sky. For 50 years, this 511 keV emission has puzzled observers and theoreticians. A key issu
e for understanding $e^+$-astrophysics is found in cosmic-ray propagation, especially at low kinetic energies (< 10 MeV). We want to shed light on how $e^+$s propagate and the resulting morphology of the emission. We approach this positron puzzle by inferring kinematic information of the 511 keV line in the inner radian of the Galaxy. This constrains propagation scenarios and source populations. By dissecting the 511 keV emission as measured with INTEGRAL/SPI, we derive spectra for individual regions in the sky. The centroid energies are converted into Doppler-shifts, representing the line-of-sight velocity along different longitudes. This results in a longitude-velocity diagram of $e^+$-annihilation. We also determine Doppler-broadenings to study annihilation conditions as they vary across the Galaxy. We find line-of-sight velocities in the 511 keV line that are consistent with zero, as well as with galactic rotation from CO measurements, and measurements of radioactive Al-26. The velocity gradient in the inner 60 deg is determined to be $4pm6$ km/s/deg. The 511 keV line width is constant as a function of longitude at $2.43pm0.14$ keV. The positronium fraction is found to be 1.0 along the galactic plane. The weak signals in the disk leave open the question whether $e^+$-annihilation is associated with the high velocities seen in Al-26 or rather with ordinarily rotating components of the Galaxys ISM. We confirm previous results that $e^+$s are slowed down to the 10 eV energy scale before annihilation, and constrain bulk Doppler-broadening contributions to <1.25 keV. Consequently, the true annihilation conditions remain unclear.
INTEGRAL observed the nova V5668 Sgr around the time of its optical maximum on March 21, 2015. Studies at UV wavelengths showed spectral lines of freshly produced Be-7. This could be measurable also in gamma-rays at 478 keV from the decay to Li-7. No
vae are also expected to synthesise Na-22 which decays to Ne-22, emitting a 1275 keV photon. About one week before the optical maximum, a strong gamma-ray flash on time-scales of hours is expected from short-lived radioactive nuclei, such as N-13 and F-18. These beta-plus-unstable nuclei should yield emission up to 511 keV, but which has never been observed. The spectrometer SPI aboard INTEGRAL pointed towards V5668 by chance. We use these observations to search for possible gamma-ray emission of decaying Be-7, and to directly measure the synthesised mass during explosive burning. We also aim to constrain possible burst-like emission days to weeks before the optical maximum using the SPI anticoincidence shield (ACS). We extract spectral and temporal information to determine the fluxes of gamma-ray lines at 478 keV, 511 keV, and 1275 keV. A measured flux value directly converts into abundances produced by the nova. The SPI-ACS rates are analysed for burst-like emission using a nova model light-curve. For the obtained nova flash candidate events, we discuss possible origins. No significant excess for the expected gamma-ray lines is found. Our upper limits on the synthesised Be-7 and Na-22 mass depend on the uncertainties of the distance to the nova: The Be-7 mass is constrained to less than $4.8times10^{-9},(d/kpc)^2$, and Na-22 to less than $2.4times10^{-8},(d/kpc)^2$ solar masses. For the Be-7 mass estimate from UV studies, the distance to V5668 Sgr must be larger than 1.2 kpc. During three weeks before the optical maximum, we find 23 burst-like events in the ACS rate, of which six could possibly be associated with V5668 Sgr.
