-
Galactic Population Synthesis of Radioactive Nucleosynthesis Ejecta
Authors:
Thomas Siegert,
Moritz M. M. Pleintinger,
Roland Diehl,
Martin G. H. Krause,
Jochen Greiner,
Christoph Weinberger
Abstract:
Diffuse gamma-ray line emission traces freshly produced radioisotopes in the interstellar gas, providing a unique perspective on the entire Galactic cycle of matter from nucleosynthesis in massive stars to their ejection and mixing in the interstellar medium. We aim at constructing a model of nucleosynthesis ejecta on galactic scale which is specifically tailored to complement the physically most…
▽ More
Diffuse gamma-ray line emission traces freshly produced radioisotopes in the interstellar gas, providing a unique perspective on the entire Galactic cycle of matter from nucleosynthesis in massive stars to their ejection and mixing in the interstellar medium. We aim at constructing a model of nucleosynthesis ejecta on galactic scale which is specifically tailored to complement the physically most important and empirically accessible features of gamma-ray measurements in the MeV range, in particular for decay gamma-rays such as $^{26}$Al, $^{60}$Fe or $^{44}$Ti. Based on properties of massive star groups, we developed a Population Synthesis Code which can instantiate galaxy models quickly and based on many different parameter configurations, such as the star formation rate, density profiles, or stellar evolution models. As a result, we obtain model maps of nucleosynthesis ejecta in the Galaxy which incorporate the population synthesis calculations of individual massive star groups. Based on a variety of stellar evolution models, supernova explodabilities, and density distributions, we find that the measured $^{26}$Al distribution from INTEGRAL/SPI can be explained by a Galaxy-wide population synthesis model with a star formation rate of $4$-$8\,\mathrm{M_{\odot}\,yr^{-1}}$ and a spiral-arm dominated density profile with a scale height of at least 700 pc. Our model requires that most massive stars indeed undergo a supernova explosion. This corresponds to a supernova rate in the Milky Way of $1.8$-$2.8$ per century, with quasi-persistent $^{26}$Al and $^{60}$Fe masses of $1.2$-$2.4\,\mathrm{M_{\odot}}$ and $1$-$6\,\mathrm{M_{\odot}}$, respectively. Comparing the simulated morphologies to SPI data suggests that a frequent merging of superbubbles may take place in the Galaxy, and that an unknown but strong foreground emission at 1.8 MeV could be present.
△ Less
Submitted 24 January, 2023;
originally announced January 2023.
-
$^{26}$Al gamma rays from the Galaxy with INTEGRAL/SPI
Authors:
Moritz M. M. Pleintinger,
Roland Diehl,
Thomas Siegert,
Jochen Greiner,
Martin G. H. Krause
Abstract:
The presence of radioactive $^{26}$Al at 1.8 MeV reflects ongoing nucleosynthesis in the Milky Way. Diffuse emission from its decay can be measured with gamma-ray telescopes in space. The intensity, line shape, and spatial distribution of the $^{26}$Al emission allow a study of these nucleosynthesis sources. The line parameters trace massive-star feedback in the interstellar medium due to its 1~My…
▽ More
The presence of radioactive $^{26}$Al at 1.8 MeV reflects ongoing nucleosynthesis in the Milky Way. Diffuse emission from its decay can be measured with gamma-ray telescopes in space. The intensity, line shape, and spatial distribution of the $^{26}$Al emission allow a study of these nucleosynthesis sources. The line parameters trace massive-star feedback in the interstellar medium due to its 1~My lifetime. We aim to deepen previous studies of the $^{26}$Al emission in the Milky Way, using all gamma-ray data including single and double events as collected with SPI on INTEGRAL from 2003 until 2020. We apply improved spectral response and background as evaluated from tracing spectral details over the entire mission. The exposure for Galactic $^{26}$Al emission is enhanced using all event types measured within SPI. We re-determine the intensity of Galactic $^{26}$Al emission across the entire sky, through maximum likelihood fits of simulated and model-built sky distributions to SPI spectra for single and for double detector hits. We find an all-sky flux of (1.84$\pm$0.03$)\times$10$^{-3}$~ph~cm$^{-2}$s$^{-1}$ in the 1.809~MeV line from $^{26}$Al, determined as fitted to sky distributions from previous observations with COMPTEL. Significant emission from higher latitudes indicate an origin from nearby massive-star groups and superbubbles, also supported by a bottom-up population synthesis model. The line centroid is found at (1809.83$\pm$0.04~keV, and line broadening from source kinematics integrated over the sky is (0.62$\pm0.3$)~keV (FWHM).
△ Less
Submitted 21 December, 2022;
originally announced December 2022.
-
The Radioactive Nuclei $^{\textbf{26}}$Al and $^{\textbf{60}}$Fe in the Cosmos and in the Solar System
Authors:
Roland Diehl,
Maria Lugaro,
Alexander Heger,
Andre Sieverding,
Xiaodong Tang,
KuoAng Li,
Ertao Li,
Carolyn L. Doherty,
Martin G. H. Krause,
Anton Wallner,
Nikos Prantzos,
Hannah E. Brinkman,
Jaqueline W. den Hartogh,
Benjamin Wehmeyer,
Andrés Yagüe López,
Moritz M. M. Pleintinger,
Projjval Banerjee,
Wei Wang
Abstract:
The cosmic evolution of the chemical elements from the Big Bang to the present time is driven by nuclear fusion reactions inside stars and stellar explosions. A cycle of matter recurrently re-processes metal-enriched stellar ejecta into the next generation of stars. The study of cosmic nucleosynthesis and of this matter cycle requires the understanding of the physics of nuclear reactions, of the c…
▽ More
The cosmic evolution of the chemical elements from the Big Bang to the present time is driven by nuclear fusion reactions inside stars and stellar explosions. A cycle of matter recurrently re-processes metal-enriched stellar ejecta into the next generation of stars. The study of cosmic nucleosynthesis and of this matter cycle requires the understanding of the physics of nuclear reactions, of the conditions at which the nuclear reactions are activated inside the stars and stellar explosions, of the stellar ejection mechanisms through winds and explosions, and of the transport of the ejecta towards the next cycle, from hot plasma to cold, star-forming gas. Due to the long timescales of stellar evolution, and because of the infrequent occurrence of stellar explosions, observational studies are challenging. Due to their radioactive lifetime of million years, the 26Al and 60Fe isotopes are suitable to characterise simultaneously the processes of nuclear fusion reactions and of interstellar transport. We describe and discuss the nuclear reactions involved in the production and destruction of 26Al and 60Fe, the key characteristics of the stellar sites of their nucleosynthesis and their interstellar journey after ejection from the nucleosynthesis sites. We connect the theoretical astrophysical aspects to the variety of astronomical messengers, from stardust and cosmic-ray composition measurements, through observation of gamma rays produced by radioactivity, to material deposited in deep-sea ocean crusts and to the inferred composition of the first solids that have formed in the Solar System. We show that considering measurements of the isotopic ratio of 26Al to 60Fe eliminate some of the unknowns when interpreting astronomical results, and discuss the lessons learned from these two isotopes on cosmic chemical evolution.
△ Less
Submitted 5 October, 2021; v1 submitted 17 September, 2021;
originally announced September 2021.
-
Steady-state nucleosynthesis throughout the Galaxy
Authors:
Roland Diehl,
Martin G. H. Krause,
Karsten Kretschmer,
Michael Lang,
Moritz M. M. Pleintinger,
Thomas Siegert,
Wei Wang,
Laurent Bouchet,
Pierrick Martin
Abstract:
Measurement and astrophysical interpretation of characteristic gamma-ray lines from nucleosynthesis was one of the prominent science goals of the INTEGRAL mission and in particular its spectrometer SPI. Emission from 26Al and from 60Fe decay lines originates from accumulated ejecta of nucleosynthesis sources, and appears diffuse in nature. 26Al and 60Fe are believed to originate mostly from massiv…
▽ More
Measurement and astrophysical interpretation of characteristic gamma-ray lines from nucleosynthesis was one of the prominent science goals of the INTEGRAL mission and in particular its spectrometer SPI. Emission from 26Al and from 60Fe decay lines originates from accumulated ejecta of nucleosynthesis sources, and appears diffuse in nature. 26Al and 60Fe are believed to originate mostly from massive star clusters. Gamma-ray observations open an interesting window to trace the fate and flow of nucleosynthesis ejecta, after they have left the immediate sources and their birth sites, and on their path to mix with ambient interstellar gas. The INTEGRAL 26Al emission image confirms earlier findings of clumpiness and an extent along the entire plane of the Galaxy, supporting its origin from massive-star groups. INTEGRAL spectroscopy resolved the line and found Doppler broadenings and systematic shifts from large-scale galactic rotation. But an excess velocity of ~200 km/s suggests that 26Al decays preferentially within large superbubbles that extend in forward directions between spiral arms. The detection of 26Al line emission from nearby Orion and the Eridanus superbubble supports this interpretation. Positrons from beta+ decays of 26Al and other nucleosynthesis ejecta have been found to not explain the morphology of positron annihilation gamma-rays at 511 keV that have been measured by INTEGRAL. The 60Fe signal measured by INTEGRAL is diffuse but too weak for an imaging interpretation, an origin from point-like/concentrated sources is excluded. The 60Fe/26Al ratio is constrained to a range 0.2-0.4. Beyond improving precision of these results, diffuse nucleosynthesis contributions from novae (through 22Na radioactivity) and from past neutron star mergers in our Galaxy (from r-process radioactivity) are exciting new prospects for the remaining mission extensions.
△ Less
Submitted 12 November, 2020;
originally announced November 2020.
-
$^{44}\rm Ti$ ejecta in young supernova remnants
Authors:
Christoph Weinberger,
Roland Diehl,
Moritz M. M. Pleintinger,
Thomas Siegert,
Jochen Greiner
Abstract:
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 t…
▽ More
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$σ$. 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.
△ Less
Submitted 27 April, 2020;
originally announced April 2020.
-
Gamma-ray Emission of 60Fe and 26Al Radioactivities in our Galaxy
Authors:
W. Wang,
T. Siegert,
Z. G. Dai,
R. Diehl,
J. Greiner,
A. Heger,
M. Krause,
M. Lang,
M. M. M. Pleintinger,
X. L. Zhang
Abstract:
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 l…
▽ More
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.
△ Less
Submitted 17 December, 2019;
originally announced December 2019.
-
Comparing simulated $^{26}$Al maps to gamma-ray measurements
Authors:
Moritz M. M. Pleintinger,
Thomas Siegert,
Roland Diehl,
Yusuke Fujimoto,
Jochen Greiner,
Martin G. H. Krause,
Mark R. Krumholz
Abstract:
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 availa…
▽ More
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.
△ Less
Submitted 14 October, 2019;
originally announced October 2019.
-
Constraints on positron annihilation kinematics in the inner Galaxy
Authors:
Thomas Siegert,
Roland M. Crocker,
Roland Diehl,
Martin G. H. Krause,
Fiona H. Panther,
Moritz M. M. Pleintinger,
Christoph Weinberger
Abstract:
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…
▽ More
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.
△ Less
Submitted 2 June, 2019;
originally announced June 2019.
-
Background modelling for $γ$-ray spectroscopy with INTEGRAL/SPI
Authors:
Thomas Siegert,
Roland Diehl,
Christoph Weinberger,
Moritz M. M. Pleintinger,
Jochen Greiner,
Xiaoling Zhang
Abstract:
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 p…
▽ More
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.
△ Less
Submitted 4 March, 2019;
originally announced March 2019.
-
Gamma-Ray Observations of Nova Sgr 2015 No. 2 with INTEGRAL
Authors:
Thomas Siegert,
Alain Coc,
Laura Delgado,
Roland Diehl,
Jochen Greiner,
Margarita Hernanz,
Pierre Jean,
Jordi Jose,
Paolo Molaro,
Moritz M. M. Pleintinger,
Volodymyr Savchenko,
Sumner Starrfield,
Vincent Tatischeff,
Christoph Weinberger
Abstract:
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…
▽ More
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.
△ Less
Submitted 19 March, 2018;
originally announced March 2018.