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.

Context: Tracing unstable isotopes produced in supernova nucleosynthesis provides a direct diagnostic of supernova explosion physics. Theoretical models predict an extensive variety of scenarios, which can be constrained through observations of the abundant isotopes $^{56}$Ni and $^{44}$Ti. Direct evidence of the latter was previously found only in two core-collapse supernova events, and appears to be absent in thermonuclear supernovae.Aims: We aim to to constrain the supernova progenitor types of Cas A, SN 1987A, Vela Jr., G1.9+0.3, SN1572, and SN1604 through their $^{44}$Ti ejecta masses and explosion kinematics. Methods: We analyzed INTEGRAL/SPI observations of the candidate sources utilizing an empirically motivated high-precision background model. We analyzed the three dominant spectroscopically resolved de-excitation lines at 68, 78, and 1157\,keV emitted in the decay chain of $^{44}$Ti. The fluxes allow the determination of the production yields of $^{44}$Ti. Remnant kinematics were obtained from the Doppler characteristics of the lines. Results: We find a significant signal for Cas A in all three lines with a combined significance of 5.4$\sigma$. The fluxes are $(3.3 \pm 0.9) \times 10^{-5}$ ph cm$^{-2}$ s$^{-1}$, and $(4.2 \pm 1.0) \times 10^{-5}$ ph cm$^{-2}$ s$^{-1}$ for the $^{44}$Ti and $^{44}$Sc decay, respectively. We obtain higher fluxes for $^{44}$Ti with our analysis of Cas A than were obtained in previous analyses. We discuss potential differences. Conclusions: We obtain a high $^{44}$Ti ejecta mass for Cas A that is in disagreement with ejecta yields from symmetric 2D models. Upper limits for the other core-collapse supernovae are in agreement with model predictions and previous studies. The upper limits we find for the three thermonuclear supernovae consistently exclude the double detonation and pure helium deflagration models as progenitors.

The signature of positron annihilation, namely the 511 keV $\gamma$-ray line, was first detected coming from the direction of the Galactic center in the 1970's, but the source of Galactic positrons still remains a puzzle. The measured flux of the annihilation corresponds to an intense steady source of positron production, with an annihilation rate on the order of $\sim10^{43}$~e$^{+}$/s. The 511 keV emission is the strongest persistent Galactic $\gamma$-ray line signal and it shows a concentration towards the Galactic center region. An additional low-surface brightness component is aligned with the Galactic disk; however, the morphology of the latter is not well constrained. The Compton Spectrometer and Imager (COSI) is a balloon-borne soft $\gamma$-ray (0.2--5 MeV) telescope designed to perform wide-field imaging and high-resolution spectroscopy. One of its major goals is to further our understanding of Galactic positrons. COSI had a 46-day balloon flight in May--July 2016 from Wanaka, New Zealand, and here we report on the detection and spectral and spatial analyses of the 511 keV emission from those observations. To isolate the Galactic positron annihilation emission from instrumental background, we have developed a technique to separate celestial signals utilizing the COMPTEL Data Space. With this method, we find a 7.2$\sigma$ detection of the 511 keV line. We find that the spatial distribution is not consistent with a single point source, and it appears to be broader than what has been previously reported.

The isotopes $^{60}$Fe and $^{26}$Al originate from massive stars and their supernovae, reflecting ongoing nucleosynthesis in the Galaxy. We studied the gamma-ray emission from these isotopes at characteristic energies 1173, 1332, and 1809 keV with over 15 years of SPI data, finding a line flux in $^{60}$Fe combined lines of $(0.31\pm 0.06) \times 10^{-3}$ ph cm$^{-2}$ s$^{-1}$ and the $^{26}$Al line flux of $(16.8\pm 0.7) \times 10^{-4}$ ph cm$^{-2}$ s$^{-1}$ above the background and continuum emission for the whole sky. Based on the exponential-disk grid maps, we characterise the emission extent of $^{26}$Al to find scale parameters $R_0 =7.0^{+1.5}_{-1.0}$ kpc and $z_0=0.8^{+0.3}_{-0.2}$ kpc, however the $^{60}$Fe lines are too weak to spatially constrain the emission. Based on a point source model test across the Galactic plane, the $^{60}$Fe emission would not be consistent with a single strong point source in the Galactic center or somewhere else, providing a hint for a diffuse nature. We carried out comparisons of emission morphology maps using different candidate-source tracers for both $^{26}$Al and $^{60}$Fe emissions, and suggests that the $^{60}$Fe emission is more likely to be concentrated towards the Galactic plane. We determine the $^{60}$Fe/$^{26}$Al $\gamma$-ray flux ratio at $(18.4\pm4.2)\,\%$ , when using a parameterized spatial morphology model. Across the range of plausible morphologies, it appears possible that $^{26}$Al and $^{60}$Fe are distributed differently in the Galaxy. Using the best fitting maps for each of the elements, we constrain flux ratios in the range 0.2--0.4. We discuss its implications for massive star models and their nucleosynthesis.

The Compton Spectrometer and Imager (COSI) is a balloon-borne gamma-ray (0.2-5 MeV) telescope designed to study astrophysical sources. COSI employs a compact Compton telescope design utilizing 12 high-purity germanium double-sided strip detectors and is inherently sensitive to polarization. In 2016, COSI was launched from Wanaka, New Zealand and completed a successful 46-day flight on NASA's new Super Pressure Balloon. In order to perform imaging, spectral, and polarization analysis of the sources observed during the 2016 flight, we compute the detector response from well-benchmarked simulations. As required for accurate simulations of the instrument, we have built a comprehensive mass model of the instrument and developed a detailed detector effects engine which applies the intrinsic detector performance to Monte Carlo simulations. The simulated detector effects include energy, position, and timing resolution, thresholds, dead strips, charge sharing, charge loss, crosstalk, dead time, and detector trigger conditions. After including these effects, the simulations closely resemble the measurements, the standard analysis pipeline used for measurements can also be applied to the simulations, and the responses computed from the simulations are accurate. We have computed the systematic error that we must apply to measured fluxes at certain energies, which is 6.3% on average. Here we describe the detector effects engine and the benchmarking tests performed with calibrations.

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 because 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 Sun's 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 Sun's 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 Galaxy's 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.

In this Astro2020 APC White Paper, we describe a Small Explorer (SMEX) mission concept called the Compton Spectrometer and Imager. COSI is a Compton telescope that covers the bandpass often referred to as the "MeV Gap" because it is the least explored region of the whole electromagnetic spectrum. COSI provides a significant improvement in sensitivity along with high-resolution spectroscopy, enabling studies of 511 keV electron-positron annihilation emission and measurements of several radioactive elements that trace the Galactic history of supernovae. COSI also measures polarization of gamma-ray bursts (GRBs), accreting black holes, and pulsars as well as detecting and localizing multimessenger sources. In the following, we describe the COSI science, the instrument, and its capabilities. We highlight many Astro2020 science WPs that describe the COSI science in depth.

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 issue 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 $4\pm6$ km/s/deg. The 511 keV line width is constant as a function of longitude at $2.43\pm0.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 Galaxy's 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.

We utilize the joint capabilities of IBIS and SPI to perform a state-resolved 20--2000\,keV analysis of the microquasar Cygnus X-1. In both LHS and HSS, the spectral analysis reveals the presence of a high-energy tail above 400\,keV in addition to a standard, Compton-like continuum in the 20--400\,keV range. We study the polarisation properties of the hard X-ray radiation, as well as of the polarisation of this high-energy tail. Polarisation is detected in the LHS in agreement with previous work. We find potential variations of the LHS polarisation with time and no polarisation detection in the HSS.

The 511 keV line from positron annihilation in the Galaxy was the first $\gamma$-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 $\gamma$-ray astronomy. With a large flux of $\sim$10$^{-3}$ $\gamma$/cm$^{2}$/s, after 15 years in operation INTEGRAL/SPI has detected the 511 keV line at $>50\sigma$ 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 $\beta^+$-decay of nucleosynthesis products, such as $^{26}$Al, 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.

The coded-mask spectrometer-telescope SPI on board the INTEGRAL observatory records photons in the energy range between 20 and 8000 keV. A robust and versatile method to model the dominating instrumental background (BG) radiation is difficult to establish for such a telescope in the rapidly changing space environment. From long-term monitoring of SPI's Germanium detectors, we built up a spectral parameter data base, which characterises the instrument response as well as the BG behaviour. We aim to build a self-consistent and broadly applicable BG model for typical science cases of INTEGRAL/SPI, based on this data base. The general analysis method for SPI relies on distinguishing between illumination patterns on the 19-element Germanium detector array from BG and sky in a maximum likelihood framework. We illustrate how the complete set of measurements, even including the exposures of the sources of interest, can be used to define a BG model. We apply our method to different science cases, including point-like, diffuse, continuum, and line emission, and evaluate the adequacy in each case. From likelihood values and the number of fitted parameters, we determine how strong the impact of the unknown BG variability is. We find that the number of fitted parameters, i.e. how often the BG has to be re-normalised, depends on the emission type (diffuse with many observations over a large sky region, or point-like with concentrated exposure around one source), and the spectral energy range and bandwidth. A unique time scale, valid for all analysis issues, is not applicable for INTEGRAL/SPI, but must and can be inferred from the chosen data set. We conclude that our BG modelling method is usable in a large variety of INTEGRAL/SPI science cases, and provides nearly systematics-free and robust results.

Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions.

We exploited observational constraints on stars, gas and nucleosynthesis ashes for the closest region of recent massive-star formation, Scorpius-Centaurus OB2, and combined them with 3D hydrodynamical simulations, in order to address physics and history for the case of the Scorpius-Centaurus superbubble. We used published cold gas observations through PLANCK survey data processing, HERSCHEL and APEX, continuum and molecular line observations. We analysed the Galactic All Sky Survey (GASS) to investigate shell structures in atomic hydrogen, and used HIPPARCOS and Gaia data in combination with interstellar absorption against stars to obtain new constraints for the distance to the Hi features. Hot gas is traced in soft X-rays via the ROSAT all sky survey. Nucleosynthesis ejecta from massive stars were traced with new INTEGRAL spectrometer observations via 26Al radioactivity. We also performed 3D hydrodynamical simulations for the Sco-Cen superbubble. To investigate the impact of massive star feedback on extended clouds, we simulate the interaction of a turbulent cloud with the hot, pressurised gas in a superbubble. The hot gas fills the tenuous regions of the cloud and compresses the denser parts. Stars formed in these dense clumps would have distinct spatial and kinematic distributions. The combined results from observations and simulations are consistent with a scenario where dense gas was initially distributed in a band elongated in the direction now occupied by the OB association. Superbubbles powered by massive stars would then repeatedly break out of the elongated parent cloud, surround and squash the denser parts of the gas sheet and thus induce more star formation. The expected spatial and kinematic distribution of stars is consistent with observations of Sco-Cen. The scenario might apply to many similar regions in the Galaxy and also to AGN-related superbubbles.

e-ASTROGAM (enhanced ASTROGAM) is a breakthrough Observatory space mission, with a detector composed by a Silicon tracker, a calorimeter, and an anticoincidence system, dedicated to the study of the non-thermal Universe in the photon energy range from 0.3 MeV to 3 GeV - the lower energy limit can be pushed to energies as low as 150 keV for the tracker, and to 30 keV for calorimetric detection. The mission is based on an advanced space-proven detector technology, with unprecedented sensitivity, angular and energy resolution, combined with polarimetric capability. Thanks to its performance in the MeV-GeV domain, substantially improving its predecessors, e-ASTROGAM will open a new window on the non-thermal Universe, making pioneering observations of the most powerful Galactic and extragalactic sources, elucidating the nature of their relativistic outflows and their effects on the surroundings. With a line sensitivity in the MeV energy range one to two orders of magnitude better than previous generation instruments, e-ASTROGAM will determine the origin of key isotopes fundamental for the understanding of supernova explosion and the chemical evolution of our Galaxy. The mission will provide unique data of significant interest to a broad astronomical community, complementary to powerful observatories such as LIGO-Virgo-GEO600-KAGRA, SKA, ALMA, E-ELT, TMT, LSST, JWST, Athena, CTA, IceCube, KM3NeT, and LISA.

In the Milky Way galaxy, positrons, which are responsible for the diffuse $511\,\mathrm{keV}$ gamma ray emission observed by space-based gamma ray observatories, are thought to annihilate predominantly through charge exchange interactions with neutral hydrogen. These charge exchange interactions can only take place if positrons have energies greater than $6.8\,\mathrm{eV}$, the minimum energy required to liberate the electron bound to the hydrogen atom and then form positronium, a short-lived bound state composed of a positron-electron pair. Here we demonstrate the importance of positron interactions with neutral alkali metals in the warm interstellar medium (ISM). Positrons may undergo charge exchange with these atoms at any energy. In particular, we show that including positron interactions with sodium at solar abundance in the warm ISM can significantly reduce the annihilation timescale of positrons with energies below $6.8\,\mathrm{eV}$ by at least an order of magnitude. We show that including these interactions in our understanding of positron annihilation in the Milky Way rules out the idea that the number of positrons in the Galactic ISM could be maintained in steady state by injection events occurring at a typical periodicity $>\mathrm{Myr}$.

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. Novae 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.8\times10^{-9}\,(d/kpc)^2$, and Na-22 to less than $2.4\times10^{-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.

Thanks to its high orbit and a set of complementary detectors providing continuous coverage of the whole sky, the INTEGRAL satellite has unique capabilities for the identification and study of the electromagnetic radiation associated to gravitational waves signals and, more generally, for multi-messenger astrophysics. Here we briefly review the results obtained during the first two observing runs of the advanced LIGO/Virgo interferometers.

INTEGRAL/SPI is a space-based coded mask telescope featuring a 19-element Germanium detector array for high-resolution gamma-ray spectroscopy, encapsulated in a scintillation detector assembly that provide a veto for background from charged particles. In space, cosmic rays irradiate spacecraft and instrument, which results in a large instrumental background from activation of those materials, and leads to deterioration of the charge collection properties of the Ge detectors. We use 13.5 years of INTEGRAL/SPI spectra for each detector and for each pointing of the satellite. Spectral fits provide us with details about separated background components. From the strongest background lines, we first determine how the spectral response changes with time. We then determine how the instrumental background components change in intensities and other characteristics, most-importantly their relative distribution among detectors. Spectral resolution of Ge detectors in space degrades with time, up to 15% within half a year, consistently for all detectors, and across the SPI energy range. Semi-annual annealing operations recover these losses, yet there is a small long-term degradation. The intensity of instrumental background varies anti-correlated to solar activity, in general. There are significant differences among different lines and with respect to continuum. Background lines are found to have a characteristic, well-defined and long-term consistent intensity ratio among detectors. We use this to categorise lines in groups of similar behaviour. The dataset of spectral-response and background parameters as fitted across the INTEGRAL mission allows studies of SPI spectral response and background behaviour in a broad perspective, and efficiently supports precision modelling of instrumental background.

The electron-positron annihilation gamma-ray signal at 511 keV in the Milky Way is investigated towards a possible dark matter interpretation. If all bulge positrons were created by dark matter particle annihilation, the satellite galaxies of the Milky Way, apparently being dominated by dark matter, should also show measurable 511 keV signals. Using INTEGRAL/SPI, we test for emission in 39 neighbouring dwarf satellite galaxies, and found a consistent trend against a dark matter scenario. One galaxy, Reticulum II, shows up as a strong source of annihilation emission, which we interpret as the presence of a microquasar, ejecting pair-plasma into the galaxy's interstellar medium.

During the first observing run of LIGO, two gravitational wave events and one lower-significance trigger (LVT151012) were reported by the LIGO/Virgo collaboration. At the time of LVT151012, the INTErnational Gamma-Ray Astrophysics Laboratory (INTEGRAL) was pointing at a region of the sky coincident with the high localization probability area of the event and thus permitted us to search for its electromagnetic counterpart (both prompt and afterglow emission). The imaging instruments on-board INTEGRAL (IBIS/ISGRI, IBIS/PICsIT, SPI, and the two JEM-X modules) have been exploited to attempt the detection of any electromagnetic emission associated with LVT151012 over 3 decades in energy (from 3 keV to 8 MeV). The omni-directional instruments on-board the satellite, i.e. the SPI-ACS and IBIS monitored the entire LVT151012 localization region at energies above 75 keV. We did not find any significant transient source that was spatially and/or temporally coincident with LVT151012, obtaining tight upper limits on the associated hard X-ray and $\gamma$-ray radiation. For typical spectral models, the upper limits on the fluence of the emission from any 1 s long-lasting counterpart of LVT151012 ranges from $F_{\gamma}=$3.5$\times$10$^{-8}$ erg cm$^{-2}$ (20 - 200 keV) to $F_{\gamma}$=7.1$\times$10$^{-7}$ erg cm$^{-2}$ (75 - 2000 keV), constraining the ratio of the isotropic equivalent energy released in the electromagnetic emission to the total energy of the gravitational waves: $E_{75-2000~keV}/E_{GW}<$4.4$\times$10$^{-5}$. Finally, we provide an exhaustive summary of the capabilities of all instruments on-board INTEGRAL to hunt for $\gamma$-ray counterparts of gravitational wave events, exploiting both serendipitous and pointed follow-up observations. This will serve as a reference for all future searches.

The measurement of gamma rays from the diffuse afterglow of radioactivity originating in massive-star nucleosynthesis is considered a laboratory for testing models, when specific stellar groups are investigated, at known distance and with well-constrained stellar population. Regions which have been exploited for such studies include Cygnus, Carina, Orion, and Scorpius-Centaurus. The Orion region hosts the Orion OB1 association and its subgroups at about 450~pc distance. We report the detection of $^{26}$Al gamma rays from this region with INTEGRAL/SPI.

The positron annihilation gamma-ray signal in the Milky Way (MW) shows a puzzling morphology: a very bright bulge and a very low surface-brightness disk. A coherent explanation of the positron origin, propagation through the Galaxy and subsequent annihilation in the interstellar medium has not yet been found. Tentative explanations involve positrons from radioactivity, X-ray binaries, and dark matter (DM). Dwarf satellite galaxies (DSGs) are believed to be DM-dominated and hence promising candidates in the search for 511 keV emission as a result of DM annihilation into electron-positron pairs. The goal of this study is to constrain possible 511 keV gamma-ray signals from 39 DSGs of the MW and to test the annihilating DM scenario. We use the spectrometer SPI on INTEGRAL to extract individual spectra for the studied objects. As the diffuse galactic emission dominates the signal, the large scale morphology of the MW has been modelled accordingly and was included in a maximum likelihood analysis. Alternatively, a distance-weighted stacked spectrum has been determined. Only Reticulum II (Ret II) shows a 3.1 sigma signal. Five other sources show tentative 2 sigma signals. The mass-to-511-keV-luminosity-ratio shows a marginal trend towards higher values for intrinsically brighter objects, opposite to the V band mass-to-light-ratio, which is generally used to uncover DM in DSGs. All derived flux values are above the level implied by a DM interpretation of the MW bulge signal. The signal from Ret II is unlikely to be related to a DM origin alone, otherwise, the MW bulge would be about 100 times brighter than what is seen. Ret II is exceptional considering the DSG sample, and rather points to enhanced recent star formation activity, if its origins are similar to processes in the MW. Understanding this emission may provide further clues regarding the origin of the annihilation emission in the MW.

Microquasars are stellar-mass black holes accreting matter from a companion star and ejecting plasma jets at almost the speed of light. They are analogues of quasars that contain supermassive black holes of $10^6$ to $10^{10}$ solar masses. Accretion in microquasars varies on much shorter timescales than in quasars and occasionally produces exceptionally bright X-ray flares. How the flares are produced is unclear, as is the mechanism for launching the relativistic jets and their composition. An emission line near 511 kiloelectronvolts has long been sought in the emission spectrum of microquasars as evidence for the expected electron-positron plasma. Transient high-energy spectral features have been reported in two objects, but their positron interpretation remains contentious. Here we report observations of $\gamma$-ray emission from the microquasar V404 Cygni during a recent period of strong flaring activity. The emission spectrum around 511 kiloelectronvolts shows clear signatures of variable positron annihilation, which implies a high rate of positron production. This supports the earlier conjecture that microquasars may be the main sources of the electron-positron plasma responsible for the bright diffuse emission of annihilation $\gamma$-rays in the bulge region of our Galaxy. Additionally, microquasars could be the origin of the observed megaelectronvolt continuum excess in the inner Galaxy.

The annihilation of positrons in the Galaxy's interstellar medium produces characteristic gamma-rays with a line at 511 keV. This emission has been observed with the spectrometer SPI on INTEGRAL, confirming a puzzling morphology with bright emission from an extended bulge-like region, and faint disk emission. Most plausible sources of positrons are believed to be distributed throughout the disk of the Galaxy. We aim to constrain characteristic spectral shapes for different spatial components in the disk and bulge with the high-resolution gamma-ray spectrometer SPI, based on a new instrumental background method and detailed multi-component sky model fitting. We confirm the detection of the main extended components of characteristic annihilation gamma-ray signatures at 58$\sigma$ significance in the line. The total Galactic line intensity amounts to $(2.7\pm0.3)\times10^{-3}~\mathrm{ph~cm^{-2}~s^{-1}}$ for our assumed spatial model. We derive spectra for the bulge and disk, and a central point-like and at the position of Sgr A*, and discuss spectral differences. The bulge shows a 511 keV line intensity of $(0.96\pm0.07)\times10^{-3}~\mathrm{ph~cm^{-2}~s^{-1}}$ and an ortho-positronium continuum equivalent to a positronium fraction of ($1.08\pm0.03$). The 2D-Gaussian representing the disk has an extent of $60^{+10}_{-5}$ deg in longitude and a latitudinal extent of $10.5^{+2.5}_{-1.5}$ deg; the line intensity is $(1.7\pm0.4)\times10^{-3}~\mathrm{ph~cm^{-2}~s^{-1}}$ with a marginal detection of ortho-positronium and an overall diffuse Galactic continuum of $(5.85\pm1.05)\times10^{-5}~\mathrm{ph~cm^{-2}~s^{-1}~keV^{-1}}$ at 511 keV. The disk shows no flux asymmetry between positive and negative longitudes, although spectral details differ. The flux ratio between bulge and disk is ($0.6\pm0.1$). The central source has an intensity of $(0.8\pm0.2)\times10^{-4}~\mathrm{ph~cm^{-2}~s^{-1}}$.

After 25 years of quiescence, the microquasar V404 Cyg entered a new period of activity in June 2015. This X-ray source is known to undergo extremely bright and variable outbursts seen at all wavelengths. It is therefore an object of prime interest to understand the accretion-ejection connections. These can, however, only be probed through simultaneous observations at several wavelengths. We made use of the INTEGRAL instruments to obtain long, almost uninterrupted observations from 2015 June 20$^{\mathrm{th}}$, 15:50 UTC to June 25$^{\mathrm{th}}$, 4:05 UTC, from the optical V-band, up to the soft $\gamma$-rays. V404 Cyg was extremely variable in all bands, with the detection of 18 flares with fluxes exceeding 6 Crab (20--40 keV) within 3 days. The flare recurrence can be as short as $\sim$ 20~min from peak to peak. A model-independent analysis shows that the $>$6 Crab flares have a hard spectrum. A simple 10--400 keV spectral analysis of the off-flare and flare periods shows that the variation in intensity is likely to be due to variations of a cut-off power law component only. The optical flares seem to be at least of two different types: one occurring in simultaneity with the X-ray flares, the other showing a delay greater than 10 min. The former could be associated with X-ray reprocessing by either an accretion disk or the companion star. We suggest that the latter are associated with plasma ejections that have also been seen in radio.

The 340-year old supernova remnant Cassiopeia A at 3.4 kpc distance is the best-studied young core-collapse supernova remnant. Nucleosynthesis yields in radioactive isotopes have been studied with different methods, in particular for production and ejection of $^{44}$Ti and $^{56}$Ni which originate from the innermost regions of the supernova. $^{44}$Ti was first discovered in this remnant, but is not seen consistently in other core-collapse sources. We analyse the observations accumulated with the SPI spectrometer on INTEGRAL, together with an improved instrumental background method, to achieve high spectroscopic resolution which enables interpretation towards a velocity constraint on $^{44}$Ti ejecta from the 1.157 MeV $\gamma$-ray line of the $^{44}$Sc decay. We observe both the hard X-ray line at 78 keV and the $\gamma$-ray line at 1157 keV from the $^{44}$Ti decay chain, at a combined significance of 3.8 $\sigma$. Measured fluxes are $(2.1\pm0.4)~10^{-5}~\mathrm{ph~cm^{-2}~s^{-1}}$ and $(3.5\pm1.2)~10^{-5}~\mathrm{ph~cm^{-2}~s^{-1}}$, which corresponds to $(1.5\pm0.4)~10^{-4}$ and $(2.4\pm0.9)~10^{-4}~\mathrm{M}_{\odot}$ of $^{44}$Ti, respectively. The measured Doppler broadening of the lines implies expansion velocities of $4300$ and $2200~\mathrm{km~s^{-1}}$, respectively. Combining our results with previous studies, we determine a more precise estimate of ejected $^{44}$Ti of $(1.37\pm0.19)~10^{-4}~\mathrm{M}_{\odot}$. The measurements of both lines are consistent with previous studies. The flux in the line originating from excited $^{44}$Ca is significantly higher than the flux determined in the lines from $^{44}$Sc. Cosmic ray acceleration within the supernova remnant may be responsible for an additional contribution to this line from nuclear de-excitation following energetic particle collisions in the remnant and swept-up material.

High energy resolution spectroscopy of the 1.8 MeV radioactive decay line of 26Al with the SPI instrument on board the INTEGRAL satellite has recently revealed that diffuse 26Al has large velocities in comparison to other components of the interstellar medium in the Milky Way. 26Al shows Galactic rotation in the same sense as the stars and other gas tracers, but reaches excess velocities up to 300 km/s. We investigate if this result can be understood in the context of superbubbles, taking into account the statistics of young star clusters and H I supershells, as well as the association of young star clusters with spiral arms. We derive energy output and 26Al mass of star clusters as a function of the cluster mass via population synthesis from stellar evolutionary tracks of massive stars. [...] We link this to the size distribution of HI supershells and assess the properties of likely 26Al-carrying superbubbles. 26Al is produced by star clusters of all masses above about 200 solar masses, roughly equally contributed over a logarithmic star cluster mass scale, and strongly linked to the injection of feedback energy. The observed superbubble size distribution cannot be related to the star cluster mass function in a straight forward manner. In order to avoid that the added volume of all superbubbles exceeds the volume of the Milky Way, individual superbubbles have to merge frequently. If any two superbubbles merge, or if 26Al is injected off-centre in a bigger HI supershell we expect the hot 26Al-carrying gas to obtain velocities of the order of the typical sound speed in superbubbles, about 300 km/s before decay. [...]

On 21 January 2014, SN2014J was discovered in M82 and found to be the closest type Ia supernova (SN Ia) in the last four decades. INTEGRAL observed SN2014J from the end of January until late June for a total exposure time of about 7 Ms. SNe Ia light curves are understood to be powered by the radioactive decay of iron peak elements of which $^{56}$Ni is dominantly synthesized during the thermonuclear disruption of a CO white dwarf (WD). The measurement of $\gamma$-ray lines from the decay chain $^{56}$Ni$\rightarrow$$^{56}$Co$\rightarrow$$^{56}$Fe provides unique information about the explosion in supernovae. Canonical models assume $^{56}$Ni buried deeply in the supernova cloud, absorbing most of the early $\gamma$-rays, and only the consecutive decay of $^{56}$Co should become directly observable through the overlaying material several weeks after the explosion when the supernova envelope dilutes as it expands. Surprisingly, with the spectrometer on INTEGRAL, SPI, we detected $^{56}$Ni $\gamma$-ray lines at 158 and 812 keV at early times with flux levels corresponding to roughly 10% of the total expected amount of $^{56}$Ni, and at relatively small velocities. This implies some mechanism to create a major amout of $^{56}$Ni at the outskirts, and at the same time to break the spherical symmetry of the supernova. One plausible explanation would be a belt accreted from a He companion star, exploding, and triggering the explosion of the white dwarf. The full set of observations of SN2014J show $^{56}$Co $\gamma$-ray lines at 847 and 1238 keV, and we determine for the first time a SN Ia $\gamma$-ray light curve. The irregular appearance of these $\gamma$-ray lines allows deeper insights about the explosion morphology from its temporal evolution and provides additional evidence for an asymmetric explosion, from our high-resolution spectroscopy and comparisons with recent models.

The measurement of gamma-ray lines from the decay chain of 56Ni provides unique information about the explosion in supernovae. The 56Ni freshly-produced in the supernova powers the optical light curve, as it emits gamma-rays upon its radioactive decay to 56Co and then 56Fe. Gamma-ray lines from 56Co decay are expected to become directly visible through the overlying white dwarf material several weeks after the explosion, as they progressively penetrate the overlying material of the supernova envelope, diluted as it expands. The lines are expected to be Doppler-shifted or broadened from the kinematics of the 56Ni ejecta. With the SPI spectrometer on INTEGRAL and using an improved instrumental background method, we detect the two main lines from 56Co decay at 847 and 1238 keV from SN2014J at 3.3 Mpc, significantly Doppler-broadened, and at intensities (3.65+/-1.21) 10^-4 and (2.27+/-0.69) 10^-4 ph cm^-2s^-1, respectively, at brightness maximum. We measure their rise towards a maximum after about 60-100 days and decline thereafter. The intensity ratio of the two lines is found consistent with expectations from 56Co decay (0.62+/-0.28 at brightness maximum, expected is 0.68). We find that the broad lines seen in the late, gamma-ray transparent phase are not representative for the early gamma-ray emission, and rather notice the emission spectrum to be complex and irregular until the supernova is fully transparent to gamma-rays, with progressive uncovering of the bulk of 56Ni. We infer that the explosion morphology is not spherically symmetric, both in the distribution of 56Ni and of the unburnt material which occults the 56Co emission. Comparing light curves from different plausible models, the resulting 56Ni mass is determined as 0.49+/-0.09 M_sol.

Type-Ia supernovae result from binary systems that include a carbon-oxygen white dwarf, and these thermonuclear explosions typically produce 0.5 M_solar of radioactive 56Ni. The 56Ni is commonly believed to be buried deeply in the expanding supernova cloud. Surprisingly, in SN2014J we detected the lines at 158 and 812 keV from 56Ni decay ({\tau}~8.8 days) earlier than the expected several-week time scale, only ~20 days after the explosion, and with flux levels corresponding to roughly 10% of the total expected amount of 56Ni. Some mechanism must break the spherical symmetry of the supernova, and at the same time create a major amount of 56Ni at the outskirts. A plausible explanation is that a belt of helium from the companion star is accreted by the white dwarf, where this material explodes and then triggers the supernova event.