Source author record

R. Diehl

R. Diehl appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

24works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

24 published item(s)

preprint2020arXiv

Properties of gamma-ray decay lines in 3D core-collapse supernova models, with application to SN 1987A and Cas A

Comparison of theoretical line profiles to observations provides important tests for supernova explosion models. We study the shapes of radioactive decay lines predicted by current 3D core-collapse explosion simulations, and compare these to observations of SN 1987A and Cas A. Both the widths and shifts of decay lines vary by several thousand kilometers per second depending on viewing angle. The line profiles can be complex with multiple peaks. By combining observational constraints from 56Co decay lines, 44Ti decay lines, and Fe IR lines, we delineate a picture of the morphology of the explosive burning ashes in SN 1987A. For M_ZAMS=15-20 Msun progenitors exploding with ~1.5 *10^51 erg, ejecta structures suitable to reproduce the observations involve a bulk asymmetry of the 56Ni of at least ~400 km/s and a bulk velocity of at least ~1500 km/s. By adding constraints to reproduce the UVOIR bolometric light curve of SN 1987A up to 600d, an ejecta mass around 14 Msun is favoured. We also investigate whether observed decay lines can constrain the neutron star (NS) kick velocity. The model grid provides a constraint V_NS > V_redshift, and applying this to SN 1987A gives a NS kick of at least 500 km/s. For Cas A, our single model provides a satisfactory fit to the NuSTAR observations and reinforces the result that current neutrino-driven core-collapse SN models can achieve enough bulk asymmetry in the explosive burning material. Finally, we investigate the internal gamma-ray field and energy deposition, and compare the 3D models to 1D approximations.

preprint2019arXiv

Gamma-ray Emission of 60Fe and 26Al Radioactivities in our Galaxy

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 $γ$-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.

preprint2016arXiv

INTEGRAL upper limits on gamma-ray emission associated with the gravitational wave event GW150914

Using observations of the INTErnational Gamma-Ray Astrophysics Laboratory (INTEGRAL), we put upper limits on the gamma-ray and hard X-ray prompt emission associated with the gravitational wave event GW150914, discovered by the LIGO/Virgo collaboration. The omni-directional view of the INTEGRAL/SPI-ACS has allowed us to constrain the fraction of energy emitted in the hard X-ray electromagnetic component for the full high-probability sky region of LIGO trigger. Our upper limits on the hard X-ray fluence at the time of the event range from $F_γ=2 \times 10^{-8}$ erg cm$^{-2}$ to $F_γ=10^{-6}$ erg cm$^{-2}$ in the 75 keV - 2 MeV energy range for typical spectral models. Our results constrain the ratio of the energy promptly released in gamma-rays in the direction of the observer to the gravitational wave energy E$_γ/$E$_\mathrm{GW}<10^{-6}$. We discuss the implication of gamma-ray limits on the characteristics of the gravitational wave source, based on the available predictions for prompt electromagnetic emission.

preprint2016arXiv

Measurements of the Soft Gamma-ray Emission from SN2014J with Suzaku

The hard X-ray detector (HXD) onboard {\it Suzaku} measured soft $γ$-rays from the Type Ia supernova SN2014J at $77\pm2$ days after the explosion. Although the confidence level of the signal is about 90\% (i.e., $2 σ$), the $3 σ$ upper limit has been derived at $< 2.2 \times10^{-4}$ ph s$^{-1}$ cm$^{-2}$ in the 170 -- 250 keV band as the first independent measurement of soft $γ$-rays with an instrument other than {\it INTEGRAL}. For this analysis, we have examined the reproducibility of the NXB model of HXD/GSO using blank sky data. We find that the residual count rate in the 90 -- 500 keV band is distributed around an average of 0.19\% with a standard deviation of 0.42\% relative to the NXB rate. The averaged residual signals are consistent with that expected from the cosmic X-ray background. The flux of SN2014J derived from {\it Suzaku} measurements taken in one snapshot at $t=77\pm2$ days after the explosion is consistent with the {\it INTEGRAL} values averaged over the period between $t=$50 and 100 days and also with explosion models of single or double degenerate scenarios. Being sensitive to the total ejecta mass surrounding the radioactive material, the ratio between continuum and line flux in the soft gamma-ray regime might distinguish different progenitor models. The {\it Suzaku} data have been examined with this relation at $t=77\pm2$ days, but could not distinguish models between single and double degenerate-progenitors. We disfavor explosion models with larger $^{56}$Ni masses than 1 $M_\odot$, from our $1 σ$ error on the 170-250 keV X-ray flux of $(1.2\pm0.7) \times10^{-4}$ ph s$^{-1}$ cm$^{-2}$.

preprint2015arXiv

Correlated optical, X-ray, and gamma-ray flaring activity seen with INTEGRAL during the 2015 outburst of V404 Cygni

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 $γ$-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.

preprint2014arXiv

5.9 keV Mn K-shell X-ray luminosity from the decay of 55Fe in Type Ia supernova models

We show that the X-ray line flux of the Mn Kalpha line at 5.9 keV from the decay of 55Fe is a promising diagnostic to distinguish between Type Ia supernova (SN Ia) explosion models. Using radiation transport calculations, we compute the line flux for two 3D explosion models: a near-Chandrasekhar mass delayed detonation and a violent merger of two white dwarfs. Both models are based on solar metallicity zero-age main sequence progenitors. Due to explosive nuclear burning at higher density, the delayed-detonation model synthesises 3.5 times more radioactive 55Fe than the merger model. As a result, we find that the peak Mn Kalpha line flux of the delayed-detonation model exceeds that of the merger model by a factor of 4.5. Since in both models the 5.9 keV X-ray flux peaks five to six years after the explosion, a single measurement of the X-ray line emission at this time can place a constraint on the explosion physics that is complementary to those derived from earlier phase optical spectra or light curves. We perform detector simulations of current and future X-ray telescopes to investigate the possibilities of detecting the X-ray line at 5.9 keV. For the delayed-detonation scenario, a line detection is feasible with Chandra up to 3 Mpc for an exposure time of 10^6 s. We find that it should be possible with currently existing X-ray instruments (with exposure times 5x10^5 s) to detect both of our models at sufficiently high S/N to distinguish between them for hypothetical events within the Local Group. The prospects for detection will be better with future missions. For example, the proposed Athena/X-IFU instrument could detect our delayed-detonation model out to a distance of 5 Mpc. This would make it possible to study future events occurring during its operational life at distances comparable to those of the recent supernovae SN 2011fe (6.4 Mpc) and SN 2014J (3.5 Mpc).

preprint2013arXiv

Gamma-ray diagnostics of Type Ia supernovae: Predictions of observables from three-dimensional modeling

Besides the fact that the gamma-ray emission due to radioactive decays is responsible for powering the light curves of Type Ia supernovae (SNe Ia), gamma rays themselves are of particular interest as a diagnostic tool because they provide a direct way to obtain deeper insights into the nucleosynthesis and the kinematics of these explosion events. Focusing on two of the most broadly discussed SN Ia progenitor scenarios - a delayed detonation in a Chandrasekhar-mass white dwarf (WD) and a violent merger of two WDs - we use three-dimensional explosion models and perform radiative transfer simulations to obtain synthetic gamma-ray spectra. Both chosen models produce the same mass of 56Ni and have similar optical properties that are in reasonable agreement with the recently observed supernova SN 2011fe. In contrast to the optical regime, the gamma-ray emission of our two chosen models proves to be rather different. The almost direct connection of the emission of gamma rays to fundamental physical processes occuring in SNe Ia permits additional constraints concerning several explosion model properties that are not easily accessible within other wavelength ranges. Proposed future MeV missions such as GRIPS will resolve all spectral details only for nearby SNe Ia, but hardness ratio and light curve measurements still allow for a distinction of the two different models at 10 and 16 Mpc for an exposure time of 10^6 s, respectively. The possibility to detect the strongest line features up to the Virgo distance will offer the opportunity to build up a first sample of SN Ia detections in the gamma-ray energy range and underlines the importance of future space observatories for MeV gamma rays.

preprint2013arXiv

Observation of SN2011fe with INTEGRAL. I. Pre--maximum phase

SN2011fe was detected by the Palomar Transient Factory on August 24th 2011 in M101 a few hours after the explosion. From the early optical spectra it was immediately realized that it was a Type Ia supernova thus making this event the brightest one discovered in the last twenty years. The distance of the event offered the rare opportunity to perform a detailed observation with the instruments on board of INTEGRAL to detect the gamma-ray emission expected from the decay chains of $^{56}$Ni. The observations were performed in two runs, one before and around the optical maximum, aimed to detect the early emission from the decay of $^{56}$Ni and another after this maximum aimed to detect the emission of $^{56}$Co. The observations performed with the instruments on board of INTEGRAL (SPI, IBIS/ISGRI, JEMX and OMC) have been analyzed and compared with the existing models of gamma-ray emission from such kind of supernovae. In this paper, the analysis of the gamma-ray emission has been restricted to the first epoch. Both, SPI and IBIS/ISGRI, only provide upper-limits to the expected emission due to the decay of $^{56}$Ni. These upper-limits on the gamma-ray flux are of 7.1 $\times$ 10$^{-5}$ ph/s/cm$^2$ for the 158 keV line and of 2.3 $\times$ 10$^{-4}$ ph/s/cm$^2$ for the 812 keV line. These bounds allow to reject at the $2σ$ level explosions involving a massive white dwarf, $\sim 1$ M$\odot$ in the sub--Chandrasekhar scenario and specifically all models that would have substantial amounts of radioactive $^{56}$Ni in the outer layers of the exploding star responsible of the SN2011fe event. The optical light curve obtained with the OMC camera also suggests that SN2011fe was the outcome of the explosion, possibly a delayed detonation although other models are possible, of a CO white dwarf that synthesized $\sim 0.55$ M$_\odot$ of $^{56}$Ni. For this specific model.

preprint2013arXiv

Observations of SN2011fe with INTEGRAL

SN2011fe was detected by the Palomar Transient Factory on August 24th 2011 in M101 few hours after the explosion. From the early spectra it was immediately realized that it was a Type Ia supernova thus making this event the brightest one discovered in the last twenty years. In this paper the observations performed with the instruments on board of INTEGRAL (SPI, IBIS/ISGRI, JEM-X and OMC) before and after the maximum of the optical light as well as the interpretation in terms of the existing models of $γ$--ray emission from such kind of supernovae are reported. All INTEGRAL high-energy have only been able to provide upper limits to the expected emission due to the decay of $^{56}$Ni. These bounds allow to reject explosions involving a massive white dwarf in the sub--Chandrasekhar scenario. On the other hand, the optical light curve obtained with the OMC camera suggests that the event was produced by a delayed detonation of a CO white dwarf that produced $\sim 0.5$ M$\odot$ of $^{56}$Ni. In this particular case, INTEGRAL would have only been able to detect the early $γ$--ray emission if the supernova had occurred at a distance of 2 -3 Mpc, although the brightest event could be visible up to distances larger by a factor two.

preprint2012arXiv

Energetic feedback and $^{26}$Al from massive stars and their supernovae in the Carina region

We study the populations of massive stars in the Carina region and their energetic feedback and ejection of $^{26}$Al. We did a census of the stellar populations in young stellar clusters within a few degrees of the Carina Nebula. For each star we estimated the mass, based on the spectral type and the host cluster age. We used population synthesis to calculate the energetic feedback and ejection of $^{26}$Al from the winds of the massive stars and their supernova explosions. We used 7 years of INTEGRAL observations to measure the $^{26}$Al signal from the region. The INTEGRAL $^{26}$Al signal is not significant with a best-fit value of about 1.5e-5 ph/cm^2/s, approximately half of the published Compton Gamma Ray Observatory (CGRO) result, but in agreement with the latest CGRO estimates. Our analysis of the stellar populations in the young clusters leads to an expected signal of half the observed value, but the results are consistent within 2 sigma. We find that the fraction of $^{26}$Al ejected in Wolf-Rayet winds is high, and the observed signal is unlikely to be caused by $^{26}$Al ejected in supernovae alone, indicating a strong wind ejection of $^{26}$Al. Due to the lack of prominent O stars, regions with ages $\gtrsim$10 Myr are often neglected in studies of OB associations. We find that in the Carina region such clusters contribute significantly to the stellar mass and the energetics of the region.

preprint2012arXiv

Galactic annihilation emission from nucleosynthesis positrons

The Galaxy hosts a widespread population of low-energy positrons revealed by successive generations of gamma-ray telescopes through a bright annihilation emission from the bulge region, with a fainter contribution from the inner disk. The exact origin of these particles remains currently unknown. We estimate the contribution to the annihilation signal of positrons generated in the decay of radioactive 26Al, 56Ni and 44Ti. We adapted the GALPROP propagation code to simulate the transport and annihilation of radioactivity positrons in a model of our Galaxy. Using plausible source spatial distributions, we explored several possible propagation scenarios to account for the large uncertainties on the transport of ~1MeV positrons in the interstellar medium. We then compared the predicted intensity distributions to the INTEGRAL/SPI observations. We obtain similar intensity distributions with small bulge-to-disk ratios, even for extreme large-scale transport prescriptions. At least half of the positrons annihilate close to their sources, even when they are allowed to travel far away. In the high-diffusion, ballistic case, up to 40% of them escape the Galaxy. In proportion, this affects bulge positrons more than disk positrons because they are injected further off the plane in a tenuous medium, while disk positrons are mostly injected in the dense molecular ring. The predicted intensity distributions are fully consistent with the observed longitudinally-extended disk-like emission, but the transport scenario cannot be strongly constrained by the current data. Nucleosynthesis positrons alone cannot account for the observed annihilation emission in the frame of our model. An additional component is needed to explain the strong bulge contribution, and the latter is very likely concentrated in the central regions if positrons have initial energies in the 100keV-1MeV range.

preprint2012arXiv

Prospect of Studying Hard X- and Gamma-Rays from Type Ia Supernovae

We perform multi-dimensional, time-dependent radiation transfer simulations for hard X-ray and gamma-ray emissions, following radioactive decays of 56Ni and 56Co, for two-dimensional delayed detonation models of Type Ia supernovae (SNe Ia). The synthetic spectra and light curves are compared with the sensitivities of current and future observatories for an exposure time of 10^6 seconds. The non-detection of the gamma-ray signal from SN 2011fe at 6.4 Mpc by SPI on board INTEGRAL places an upper limit for the mass of 56Ni of \lesssim 1.0 Msun, independently from observations in any other wavelengths. Signals from the newly formed radioactive species have not been convincingly measured yet from any SN Ia, but the future X-ray and gamma-ray missions are expected to deepen the observable horizon to provide the high energy emission data for a significant SN Ia sample. We predict that the hard X-ray detectors on board NuStar (launched in 2012) or ASTRO-H (scheduled for launch in 2014) will reach to SNe Ia at \sim15 Mpc, i.e., one SN every few years. Furthermore, according to the present results, the soft gamma-ray detector on board ASTRO-H will be able to detect the 158 keV line emission up to \sim25 Mpc, i.e., a few SNe Ia per year. Proposed next generation gamma-ray missions, e.g., GRIPS, could reach to SNe Ia at \sim20 - 35 Mpc by MeV observations. Those would provide new diagnostics and strong constraints on explosion models, detecting rather directly the main energy source of supernova light.

preprint2011arXiv

Fermi/GBM observations of the ultra-long GRB 091024: A burst with an optical flash

In this paper we examine gamma-ray and optical data of GRB 091024, a gamma-ray burst (GRB) with an extremely long duration of T90~1020 s, as observed with the Fermi Gamma-Ray Burst Monitor (GBM). We present spectral analysis of all three distinct emission episodes using data from Fermi/GBM. Because of the long nature of this event, many ground-based optical telescopes slewed to its location within a few minutes and thus were able to observe the GRB during its active period. We compare the optical and gamma-ray light curves. Furthermore, we estimate a lower limit on the bulk Lorentz factor from the variability and spectrum of the GBM light curve and compare it with that obtained from the peak time of the forward shock of the optical afterglow. From the spectral analysis we note that, despite its unusually long duration, this burst is similar to other long GRBs, i.e. there is spectral evolution (both the peak energy and the spectral index vary with time) and spectral lags are measured. We find that the optical light curve is highly anti-correlated to the prompt gamma-ray emission, with the optical emission reaching the maximum during an epoch of quiescence in the prompt emission. We interpret this behavior as the reverse shock (optical flash), expected in the internal-external shock model of GRB emission but observed only in a handful of GRBs so far. The lower limit on the initial Lorentz factor deduced from the variability time scale ($Γ_{min}=195_{-110}^+{90}$)is consistent within the error to the one obtained using the peak time of the forward shock ($Γ_0=120$) and is also consistent with Lorentz factors of other long GRBs.

preprint2011arXiv

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 $γ$-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.

preprint2011arXiv

Quasi-Periodic Pulsations in Solar Flares: new clues from the Fermi Gamma-Ray Burst Monitor

In the last four decades it has been observed that solar flares show quasi-periodic pulsations (QPPs) from the lowest, i.e. radio, to the highest, i.e. gamma-ray, part of the electromagnetic spectrum. To this day, it is still unclear which mechanism creates such QPPs. In this paper, we analyze four bright solar flares which show compelling signatures of quasi-periodic behavior and were observed with the Gamma-Ray Burst Monitor (\gbm) onboard the Fermi satellite. Because GBM covers over 3 decades in energy (8 keV to 40 MeV) it can be a key instrument to understand the physical processes which drive solar flares. We tested for periodicity in the time series of the solar flares observed by GBM by applying a classical periodogram analysis. However, contrary to previous authors, we did not detrend the raw light curve before creating the power spectral density spectrum (PSD). To assess the significance of the frequencies we made use of a method which is commonly applied for X-ray binaries and Seyfert galaxies. This technique takes into account the underlying continuum of the PSD which for all of these sources has a P(f) ~ f^{-α} dependence and is typically labeled red-noise. We checked the reliability of this technique by applying it to a solar flare which was observed by the Reuven Ramaty High-Energy Solar Spectroscopic Imager (RHESSI) which contains, besides any potential periodicity from the Sun, a 4 s rotational period due to the rotation of the spacecraft around its axis. While we do not find an intrinsic solar quasi-periodic pulsation we do reproduce the instrumental periodicity. Moreover, with the method adopted here, we do not detect significant QPPs in the four bright solar flares observed by GBM. We stress that for the purpose of such kind of analyses it is of uttermost importance to appropriately account for the red-noise component in the PSD of these astrophysical sources.

preprint2011arXiv

Rest-frame properties of 32 gamma-ray bursts observed by the Fermi Gamma-Ray Burst Monitor

Aims: In this paper we study the main spectral and temporal properties of gamma-ray bursts (GRBs) observed by Fermi/GBM. We investigate these key properties of GRBs in the rest-frame of the progenitor and test for possible intra-parameter correlations to better understand the intrinsic nature of these events. Methods: Our sample comprises 32 GRBs with measured redshift that were observed by GBM until August 2010. 28 of them belong to the long-duration population and 4 events were classified as short/hard bursts. For all of these events we derive, where possible, the intrinsic peak energy in the $νF_ν$ spectrum (\eprest), the duration in the rest-frame, defined as the time in which 90% of the burst fluence was observed (\tninetyrest) and the isotropic equivalent bolometric energy (\eiso). Results: The distribution of \eprest has mean and median values of 1.1 MeV and 750 keV, respectively. A log-normal fit to the sample of long bursts peaks at ~800 keV. No high-\ep population is found but the distribution is biased against low \ep values. We find the lowest possible \ep that GBM can recover to be ~ 15 keV. The \tninetyrest distribution of long GRBs peaks at ~10 s. The distribution of \eiso has mean and median values of $8.9\times 10^{52}$ erg and $8.2 \times 10^{52}$ erg, respectively. We confirm the tight correlation between \eprest and \eiso (Amati relation) and the one between \eprest and the 1-s peak luminosity ($L_p$) (Yonetoku relation). Additionally, we observe a parameter reconstruction effect, i.e. the low-energy power law index $α$ gets softer when \ep is located at the lower end of the detector energy range. Moreover, we do not find any significant cosmic evolution of neither \eprest nor \tninetyrest.

preprint2010arXiv

DUAL Gamma-Ray Mission

Gamma-ray astronomy presents an extraordinary scientific potential for the study of the most powerful sources and the most violent events in the Universe. In order to take full advantage of this potential, the next generation of instrumentation for this domain will have to achieve an improvement in sensitivity over present technologies of at least an order of magnitude. The DUAL mission concept takes up this challenge in two complementary ways: a very long observation of the entire sky, combined with a large collection area for simultaneous observations of Type Ia SNe. While the Wide-Field Compton Telescope (WCT) accumulates data from the full gamma-ray sky (0.1-10 MeV) over the entire mission lifetime, the Laue-Lens Telescope (LLT) focuses on 56Co emission from SNe Ia (0.8-0.9 MeV), collecting gamma-rays from its large area crystal lens onto the WCT. Two separated spacecraft flying in formation will maintain the DUAL payloads at the lens' focal distance.

preprint2010arXiv

Probing the evolving massive star population in Orion with kinematic and radioactive tracers

We assemble a census of the most massive stars in Orion, then use stellar isochrones to estimate their masses and ages, and use these results to establish the stellar content of Orion's individual OB associations. From this, our new population synthesis code is utilized to derive the history of the emission of UV radiation and kinetic energy of the material ejected by the massive stars, and also follow the ejection of the long-lived radioactive isotopes 26Al and 60Fe. In order to estimate the precision of our method, we compare and contrast three distinct representations of the massive stars. We compare the expected outputs with observations of 26Al gamma-ray signal and the extent of the Eridanus cavity. We find an integrated kinetic energy emitted by the massive stars of 1.8(+1.5-0.4)times 10^52 erg. This number is consistent with the energy thought to be required to create the Eridanus superbubble. We also find good agreement between our model and the observed 26Al signal, estimating a mass of 5.8(+2.7-2.5) times 10^-4 Msol of 26Al in the Orion region. Our population synthesis approach is demonstrated for the Orion region to reproduce three different kinds of observable outputs from massive stars in a consistent manner: Kinetic energy as manifested in ISM excavation, ionization as manifested in free-free emission, and nucleosynthesis ejecta as manifested in radioactivity gamma-rays. The good match between our model and the observables does not argue for considerable modifications of mass loss. If clumping effects turn out to be strong, other processes would need to be identified to compensate for their impact on massive-star outputs. Our population synthesis analysis jointly treats kinematic output and the return of radioactive isotopes, which proves a powerful extension of the methodology that constrains feedback from massive stars.

preprint2010arXiv

The 511 keV emission from positron annihilation in the Galaxy

The first gamma-ray line originating from outside the solar system that was ever detected is the 511 keV emission from positron annihilation in the Galaxy. Despite 30 years of intense theoretical and observational investigation, the main sources of positrons have not been identified up to now. Observations in the 1990's with OSSE/CGRO showed that the emission is strongly concentrated towards the Galactic bulge. In the 2000's, the SPI instrument aboard ESA's INTEGRAL gamma-ray observatory allowed scientists to measure that emission across the entire Galaxy, revealing that the bulge/disk luminosity ratio is larger than observed in any other wavelength. This mapping prompted a number of novel explanations, including rather "exotic ones (e.g. dark matter annihilation). However, conventional astrophysical sources, like type Ia supernovae, microquasars or X-ray binaries, are still plausible candidates for a large fraction of the observed total 511 keV emission of the bulge. A closer study of the subject reveals new layers of complexity, since positrons may propagate far away from their production sites, making it difficult to infer the underlying source distribution from the observed map of 511 keV emission. However, contrary to the rather well understood propagation of high energy (>GeV) particles of Galactic cosmic rays, understanding the propagation of low energy (~MeV) positrons in the turbulent, magnetized interstellar medium, still remains a formidable challenge. We review the spectral and imaging properties of the observed 511 keV emission and we critically discuss candidate positron sources and models of positron propagation in the Galaxy.

preprint2009arXiv

Using population synthesis of massive stars to study the interstellar medium near OB associations

We developed a new population synthesis code for groups of massive stars, where we model the emission of different forms of energy and matter from the stars of the association. In particular, the ejection of the two radioactive isotopes 26Al and 60Fe is followed, as well as the emission of hydrogen ionizing photons, and the kinetic energy of the stellar winds and supernova explosions. We investigate various alternative astrophysical inputs and the resulting output sensitivities, especially effects due to the inclusion of rotation in stellar models. As the aim of the code is the application to relatively small populations of massive stars, special care is taken to address their statistical properties. Our code incorporates both analytical statistical methods applicable to small populations, as well as extensive Monte Carlo simulations. We find that the inclusion of rotation in the stellar models has a large impact on the interactions between OB associations and their surrounding interstellar medium. The emission of 26Al in the stellar winds is strongly enhanced, compared to non-rotating models with the same mass-loss prescription. This compensates the recent reductions in the estimates of mass-loss rates of massive stars due to the effects of clumping. Despite the lower mass-loss rates, the power of the winds is actually enhanced for rotating stellar models. The supernova power (kinetic energy of their ejecta) is decreased due to longer lifetimes of rotating stars, and therefore the wind power dominates over supernova power for the first 6 Myr after a burst of star-formation. For populations typical of nearby star-forming regions, the statistical uncertainties are large and clearly non-Gaussian.

preprint2005arXiv

SPI/INTEGRAL observation of the Galactic central radian: contribution of discrete sources and implication for the diffuse emission

The INTEGRAL observatory has been performing a deep survey of the Galactic central radian since 2003, with the goal of both extracting a catalog of sources and gaining insight into the Galactic diffuse emission. This paper concentrates on the estimation of the total point sources emission contribution. It is now clear that unresolved point sources contribute to the observed diffuse emission; the increasing sensitivity of instruments with time has lead to a steady decrease in estimates of this ``diffuse emission''. We have analysed the first year data obtained with the spectrometer and imager SPI on board INTEGRAL. First, a catalog of 63 hard X-ray sources detected, time-averaged, during our 2003 Galactic plane survey, is derived. Second, after extracting the spectra of the sources detected by SPI, their combined contribution is compared to the total (resolved and unresolved) emission from the Galactic ridge. The data analysis is complex: it requires us to split the total emission into several components, as discrete sources and diffuse emission are superimposed in SPI data. The main result is that point source emission dominates in the hard X-ray/soft $γ$-ray domain, and contributes around 90 % of the total emission around 100 keV, while above 250 keV, diffuse electron-positron annihilation, through its three-photon positronium continuum with a positronium fraction $\sim$ 0.97 and the 511 keV electron-positron line, dominates over the sources.

preprint2002arXiv

COMPTEL upper limits for the 56Co gamma-ray emission from SN1998bu

Supernova 1998bu in the galaxy M96 was observed by COMPTEL for a total of 88 days starting 17 days after the explosion. We searched for a signal in the 847 keV and 1238 keV lines of radioactive 56Co from this type Ia supernova. Using several different analysis methods, we did not detect SN1998bu. Our measurements should have been sensitive enough to detect 60Co gamma-rays as predicted from supernova models. Our 2-sigma flux limit is 2.3 10^{-5} photons cm^{-2} s^{-1}; this would correspond to 0.35 solar mass of ejected 56Ni, if SN1998bu were at a distance of 11.3 Mpc and transparent to MeV gamma rays for the period of our measurements. We discuss our measurements in the context of common supernova models, and conclude disfavoring a supernova event with large mixing and major parts of the freshly-generated radioactivity in outer layers.

preprint2001arXiv

Radioactivities in Population Studies: 26Al and 60Fe from OB Associations

The observation of the interstellar 1.809 MeV decay-line of radioactive 26Al by the imaging gamma-ray telescope COMPTEL have let to the conclusion, that massive stars and their subsequent core-collapse supernovae are the dominant sources of the interstellar 26Al abundance. Massive stars are known to affect the surrounding interstellar medium by their energetic stellar winds and by the emission of ionising radiation. We present a population synthesis model allowing the correlated investigation of the gamma-ray emission characteristics with integrated matter, kinetic energy and extreme ultra-violet radiation emission of associations of massive stars. We study the time evolution of the various observables. In addition, we discuss systematic as well as statistical uncertainties affecting the model. Beside uncertainties in the input stellar physics such as stellar rotation, mass loss rates or internal mixing modifications due to a unknown binary component may lead to significant uncertainties.