$^{44}\rm Ti$ ejecta in young supernova remnants

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.


Introduction
Supernova explosions play a crucial role in the chemical and kinematic evolution of the Universe. Self-consistent detailed models for the explosion mechanism and the ensuing kinematics of the ejected material are still lacking. Despite the frequent occurrence of supernova explosions (1.9 ± 1.1 per century in the Galaxy ), observational constraints on the explosion mechanisms are still sparse, due to the large variety of progenitor systems and large parameter space of their models. For modeling core-collapse supernova explosions, reviving the stalled shock and triggering an explosion presents a major challenge (Janka 2012;Burrows et al. 2018), as it has long been understood that the prompt explosion mechanism following core bounce cannot explode the star. Energy deposition by neutrinos in a gain region close to the stalled shock is considered as the driving force of the shock revival. Heating-induced creation of hydrodynamic effects, observed in 2D and 3D models, such as neutrino-driven convection and the standing accretion shock instability (SASI) (Blondin & Shaw 2007;Marek & Janka 2009;Bruenn et al. 2013;Hix et al. 2016;Bruenn et al. 2016) enhance the explodability, however, different implementation schemes favor the dominance of either neutrino-driven convection or the SASI (Pan et al. 2016;Summa et al. 2016;Skinner et al. 2016). Inclusion of additional, microphysical effects, for example, strangeness corrections (Melson et al. 2015) and rota-tion from the progenitor Summa et al. 2018;Takiwaki et al. 2016;Iwakami et al. 2014) affect the explodability. Model calculations of supernova explosions performed by various groups lead to successful explosions in a mass range 8 M/M 25. Asymmetries evolved from hydrodynamic effects become frozen in the explosion and are reflected in the kinematics and mass distribution of the ejecta (Nomoto et al. 1995;Buras et al. 2006;Fryer & Kusenko 2006;Takiwaki et al. 2012;Wongwathanarat et al. 2015;Orlando et al. 2016;Wongwathanarat et al. 2017). The variety of observed type Ia supernova luminosity and temporal behavior of light curves in the first few hundred days (Phillips et al. 1999) cannot be reproduced by a single progenitor type (Wang & Han 2012;Hillebrandt et al. 2013). Mergers of binary white dwarfs (van Kerkwijk et al. 2010;Pakmor et al. 2010;Ruiter et al. 2012;Pakmor et al. 2013;Kashyap et al. 2018) and mass accretion on single white dwarf stars (single degenerate scenario) are expected to lead to a central thermonuclear runaway, disrupting the white dwarf star in the process. On single degenerates, both stable mass accretion on white dwarfs towards the Chandresekhar mass limit (Nomoto et al. 1984;Parthasarathy et al. 2007;Maeda et al. 2010;Woosley & Kasen 2011;Hachisu et al. 2011;Seitenzahl et al. 2013;Chen et al. 2014;Fink et al. 2014) and surface helium detonation (Livne 1990;Fink et al. 2010 & Nomoto 2020) can increase the central density sufficiently to ignite nuclear fusion. The high opacity of the ejected material does not allow direct observations of the first stages of the explosions. One of the most promising methods of deducing physical constraints for the explosion mechanisms is to directly observe the decay of radioactive isotopes produced during explosive nucleosynthesis. The best candidates for observing supernova interior physics through their nucleosynthesis imprints are the isotopes 56 Ni and 44 Ti due to their high abundances (≈ 10 −2 M for core-collapse supernovae, 0.5 M for SNe Ia in 56 Ni, 10 −6 − 10 −4 M in 44 Ti, see also below). Both isotopes have, for gamma-ray observations, ideal radioactive decay times: 6 d and 77 d for 56 Ni and its daughter nucleus 56 Co, respectively, allow only a very short postexplosion observation window in which the majority of emitted gamma-rays are expected to be absorbed by the dense stellar material. With a half life of 58.9 ± 0.3 yr (Ahmad et al. 2006), a considerable amount of 44 Ti is still decaying even after centuries when the supernova remnant has long become optically thin to X-and gamma-ray emission. In contrast to the decay of synthesized 26 Al and 60 Fe with half lives of 700 kyr and 2 Myr, producing a diffuse emission throughout the Milky Way (Plüschke et al. 2001;Bouchet et al. 2015;Siegert & Diehl 2017), the emission signature of 44 Ti decay is expected to be the one of individual point sources. Spatially co-produced in core-collapse supernovae, 44 Ti and 56 Ni are mainly synthesized during alpha-rich freeze out (Woosley et al. 1973) deep in the central region of the supernova, where nucleosynthesis is strongly dependent on the thermodynamic conditions of the inner ejecta (Magkotsios et al. 2010). Due to the hydrodynamic instabilities required for an effective explosion, mixing and asymmetric expansion of the burning volume invalidates the 1D model characteristics of a mass cut, an idealized radius (Woosley & Weaver 1995), separating the gravitationally bound material from the ejecta. Since the nuclear burning occurs close to such a mass cut, the final amount of ejected 44 Ti is subject to uncertainty and ranges between 10 −5 − 10 −4 M (Timmes et al. 1996;Limongi & Chieffi 2018). Type Ia supernovae typically produce ≈ 10 −6 M 44 Ti in centrally ignited pure deflagration models, and up to a few times 10 −5 M in delayed detonation models (Maeda et al. 2010;Seitenzahl et al. 2013;Fink et al. 2014). Double detonation and surface He deflagration of sub-Chandrasekhar mass white dwarfs produce yields up to 10 −3 M of 44 Ti Woosley & Kasen 2011;Moll & Woosley 2013). The subclass of helium surface explosions synthesizes 56 Ni very little, potentially producing subluminous type Ia supernova events. Simulations (Waldman et al. 2011) and observations (Perets et al. 2010) suggest special configurations of binary white dwarfs, which, when exploding as peculiar type Ia supernovae, can produce of the order of 10 −2 M of 44 Ti. Observations of late optical spectra allow for the determination of mass ratios of synthesized elements, constraining the burning mechanisms (Eriksen et al. 2009;Pakmor et al. 2010;Jerkstrand et al. 2015;Maguire et al. 2018;Mori et al. 2018). However, the total ejected mass remains a free parameter. Nucleosynthesis yields can be estimated from the bolometric light curves of explosion, however, inferring the ejected mass of radioactive material is highly model dependent . Comparing optical and infrared spectra to the late-time light curve, an ejected 44 Ti mass of (1.5 ± 0.5) × 10 −4 M is obtained for SN 1987A (Jerkstrand et al. 2015) in agreement with NuSTAR findings of (1.5 ± 0.3) × 10 −4 M (Boggs et al. 2015) based on the 68 and 78 keV line. This is, however, in disagreement with the results obtained from a multicomponent long-time light-curve model of (0.6 ± 0.2) × 10 −4 M (Seitenzahl et al. 2013) and direct detection of 44 Ti decay with INTEGRAL/IBIS at 68 and 78 keV of (3.1 ± 0.8) × 10 −4 M (Grebenev et al. 2012) . Direct observational evidence for the production of 44 Ti can be obtained through the decay chain of 44 Ti → 44 Sc → 44 Ca, where the prominent decay lines arise at energies of 68 keV and 78 keV for the decay of 44 Ti, and at 1157 keV for the 44 Sc decay. These lines have a probability of 93.0 %, 96.4 %, and 99.9 % per decay, respectively (Chen et al. 2011). In addition, a fluorescence photon is emitted from shell transitions in 44 Sc with a probability of 16.7 % at 4.1 keV. With a significant difference of half life times of 60 yr in the first decay and 4 h (Audi et al. 2003) in the subsequent decay, the activity of all three decay channels can be safely assumed to be identical after correcting for the branching ratios. In this work, we utilized INTEGRAL/SPI data to search for emission of 44 Ti in all three decay lines simultaneously in the young close by supernova remnants Cassiopeia A, Tycho, Kepler, G1.9+0.3, Vela Junior, and the extragalactic (but very young) SN 1987A. We aim to constrain both 44 Ti ejecta yields and explosion kinematics for remnants with ages of less than a few centuries where a substantial amount of 44 Ti may still be present. We re-evaluated previous analyses of these remnants concerning the decay of 44 Ti (see Section 2) in the 68 and 78 keV lines and improved on them by including decay signature of the 44 Sc daughter nucleus at 1157 keV, which can only be seen in SPI. The 1157 keV line was studied for the remnant Cassiopeia A and Vela Jr. This paper is structured as follows: Section 2 provides an overview of the six target remnants. Section 3 describes SPI data and spectral analysis, followed by Section 4, which includes results for the six remnants. Finally, in Section 5, we discuss our results and give a summary of astrophysical implications.

Young supernova remnants
The astrophysical parameters of the six most promising candidates to coherently observe the 44 Ti decay chain are listed in Tab. 1. Evidence for the signature of the 44 Ti decay has been claimed in the majority of these historic supernova explosions. Within SPI's narrow-line sensitivity verification of the 44 Ti signal and detection of the ensuing decay of 44 Sc is feasible within a limited age-distance volume including these remnants. We searched for signatures of the 44 Ti decay at 68 and 78 keV due to the de-excitation of 44 Sc*, and at 1157 keV due to the subsequent de-excitation of 44 Ca*. We utilized data from the INTE-GRAL mission from 2003 to 2019, meaning INTEGRAL revolutions 43 to 2047. For our study, we applied an average date of 01.01.2011 AD as an observation date for all calculations and included age uncertainties in our mass derivation.

Cassiopeia A
Cassiopeia A is one of the best studied supernova remnants in the Milky Way (Vink 2004). With an approximate age of 340 years, it is the youngest known Galactic supernova remnant, attributed to a core-collapse explosion. From the detection of hydrogen and weak helium lines in a supernova light echo, attributed to Cassiopeia A, the explosion is characterized as a type IIb supernova (Krause et al. 2008a). Explosions of this type are typically produced by 15-25 M stars (Young et al. 2006). The decay of 44 Ti has been consistently measured (Siegert et al. 2015;Grefenstette et al. 2014;Iyudin et al. 1994;Vink et al. 2001) Table 1. Astrophysical parameters for the analyzed objects including the distances of supernova remnants used in the mass estimate for 44 Ti. The age of the remnant is given for 01.01.2011 AD, which is the average observation date also used in Eq. 3. Age uncertainties of the time bin are included in the mass estimates. Distance to Vela Jr. is estimated from 44 Ti yields. Distance to G1.9+0.3 is estimated from absorption towards the Galactic center. Exposure is dead time corrected, effective INTEGRAL/SPI exposure. References: (1) (Alarie et al. 2014); (2) (Reed et al. 1995); (3) (Krause et al. 2008a); (4) (Pietrzyński et al. 2019); (5) (Gall et al. 2015); (6) (Aschenbach et al. 1999); (7) (Hayato et al. 2010); (8) (Krause et al. 2008b); (9) (Sankrit et al. 2016); (10) (Reynolds et al. 2007); (11) (Reynolds et al. 2008); (12) (Borkowski et al. 2013

SN 1987A
The explosion of SN 1987A occurred in the Large Magellanic Cloud (LMC) on February 24 1987, giving rise to a peculiar light curve containing a plateau phase (type II-P supernova). Detection of 56 Co lines (Matz et al. 1988;Tueller et al. 1990) for the first time directly has confirmed that supernova light is indeed powered by this isotope, produced in the inner regions at the time of core-collapse. Recent refinements of astronomical precision yield a distance of 49.6 ± 0.5 kpc (Pietrzyński et al. 2019) to the LMC. Direct proof for the decay of 44 Ti has been reported from measurements with IBIS/INTEGRAL and NuSTAR in the hard X-ray lines of 44 Sc. However, the two measurements show a discrepancy in the ejecta mass, with (3.1 ± 0.8) × 10 −4 M and (1.5 ± 0.3) × 10 −4 M , respectively (Grebenev et al. 2012;Boggs et al. 2015). NuSTAR constrains the ejecta kinematics to an expansion velocity of less than 4100 km s −1 , and a bulk Doppler shift suggests an asymmetric explosion.

Vela Jr.
COMPTEL discovered a significant gamma-ray emission in the energy range centered at the 1157 keV 44 Ca line. The signal was located in the direction of the Vela region with a flux of (3.8 ± 0.7) × 10 −5 ph cm −2 s −1 (Iyudin et al. 1998). This emission was attributed to the decay of 44 Ti in a previously unknown type II supernova remnant, which had been identified through detailed analysis of the X-ray emission of this region, and was called RX J0852.0-4622 or "Vela Jr." (Aschenbach 1998). An estimated age of ≈680 yr and a distance of ≈200 pc is derived (Aschenbach et al. 1999 G1.9+0.3 is presumably the youngest supernova remnant seen in the Galaxy so far. It has been identified by Reynolds et al. (2008) in the years 1985 and 2008 in the radio and X-ray regime. It is identified as a type Ia explosion. Radio observations suggest a distance of 8.5 kpc, placing the remnant in the Galactic center (Reynolds et al. 2008). From the apparent increase in size, an age of ≈ 100 yr is deduced, which suggests a very high expansion velocity of 14000 km s −1 for the shock front. Using Chandra data, a detection of a soft X-ray component at 4.1 keV from the fluorescence line of 44 Sc has been reported (Borkowski et al. 2010) inferring a 44 Ti ejecta mass of (1 − 7) × 10 −5 M . Extrapolating the detected flux in the 44 Sc fluorescence line to the hard X-ray lines provides a line flux estimate for the 68 keV 44 Ti line, which is below the upper limits determined from NuSTAR and IBIS/INTEGRAL instruments ((0.7 − 1.5) × 10 −5 ph cm −2 s −1 , 9 × 10 −6 ph cm −2 s −1 ; Zoglauer et al. 2015;Tsygankov et al. 2016, respectively).

Tycho
The Tycho supernova remnant is attributed to an explosion of type Ia, from measurements of the light echo and a comparison of model light curves with X-ray spectra (Badenes et al. 2006;Krause et al. 2008a). Tycho exploded in 1572 AD. Detection of synchrotron emission from a thin shell supports the idea of particle acceleration in young supernova remnants (e.g., Slane et al. 2014). X-ray line emission from intermediate to iron group elements in the interior of the supernova remnant has been found to be clumped (XMM Newton Miceli et al. 2015), and the hard Xray lines of the 44 Ti decay chain were also detected (Swift/BAT Troja et al. 2014). These measurements suggest the presence of titanium both in the shocked shell and the interior region of the remnants. However, upper limits obtained from NuSTAR measurements exclude the presence of 44 Ti within a 2' remnant radius at the Swift/BAT detection level over a large range of expansion velocities (Lopez et al. 2015). The distance to the remnant is somewhat uncertain, with estimates ranging from 1.7 to 5.1 kpc (Hayato et al. 2010;Slane et al. 2014;Albinson et al. 1986;Völk et al. 2008). We adopted a distance of 4.1 ± 1 kpc (Hayato et al. 2010).

Kepler
Johannes Kepler detected this supernova in 1604 AD. Also named G4.5+6.8, this is the youngest Galactic supernova for which an optical transient has been observed. This supernova occurred at a distance of 5.1 ± 0.8 kpc (Sankrit et al. 2016), and is located 460 − 600 pc above the Galactic plane. Due to the detection of strong iron lines in the ejecta, it is attributed to a Fig. 1. Intensity distribution of a celestial source located at the center of the plot (blue dot) (Siegert et al. 2019). The orientation of the SPI instrument is changed in steps of two degrees to visit the locations marked on the map with black dots. The instrument remains centered on each celestial position for ≈2000s, which is called "one pointing". The source located at the blue position is folded through the instrumental response function to determine the intensity distribution in the detector plane.
type Ia explosion (Reynolds et al. 2007), also supported by the absence of a central compact object. 44 Ti decay lines have not been found, with upper limits of 1.8 × 10 −5 ph cm −2 s −1 in the 1157 keV line (Dupraz et al. 1997) determined from COMP-TEL data and 6.3 × 10 −6 ph cm −2 s −1 determined from INTE-GRAL/IBIS for the 68 and 78 keV line (Tsygankov et al. 2016).

Instrument and analysis method
ESA's gamma-ray space observatory INTEGRAL (Winkler et al. 2003) carries two main instruments on board, the imager IBIS and the spectrometer SPI. The SPI camera ) is a germanium detector array consisting of 19 hexagonally shaped detectors, optimized for high-resolution spectroscopy in the energy range between 18 keV and 8 MeV, with a spectral resolution of ≈ 2.3 keV (full width at half maximum, FWHM) at 1 MeV. SPI electronics records 16384 energy channels in the range of 18 keV − 2 MeV, which is called SPI's "low-energy range". We analyzed SPI data for the signatures of the 44 Ti decay lines, which have centroid energies in the laboratory at 68, 78 keV and 1157 keV. Our analysis is focused on the energy bands 30 − 100 keV and 1090 − 1210 keV in order to constrain potential underlying continuum emission at lower energies and account for potential large line broadening at higher energies. SPI data after initial energy calibration and pre-processing comprises spectra in 0.5 keV bins for each of the 19 detectors, accumulated over exposures of typically 2000 s, called pointings. The orientation of the satellite is shifted by ≈ 2 • between each consecutive exposure in a rectangular-shaped dithering pattern consisting of 5×5 sets of coordinates around the observation target. We include data in our analysis in which the celestial objects of interest are within the partially coded field of view of SPI of 34 • × 34 • . In general, celestial photons entering the aperture of SPI are partially blocked by a coded mask placed 171 cm above the detector plane, imprinting shadow grams on the camera. For an idealized source at long integration times, this creates rela-tive detector intensity distributions (detector patterns), since the absolute number of measured photons per detector is governed by the visibility of the source through the mask (Fig. 1). Sources within a field of view of 16 • ×16 • are fully coded by the tungsten mask, with decreasing coding fraction towards the coding limit at 34 • × 34 • . The main challenge of SPI data analysis consists of extracting the sparsely populated celestial detector pattern above a large, underlying, instrumental background. The latter is introduced by interaction and activation of satellite and instrument material by cosmic ray bombardment. In our spectroscopic analysis method, we compare the combination of the celestial detector patterns and the background detector patterns to the time series of measured patterns for the 19 detectors by fitting (time dependent) scaling parameters for both contributors. The celestial detector intensity distribution is calculated by applying the energy-dependent image response function (IRF) specific to SPI's tungsten mask and the re-orientations during the dithering exposures. The model of the background consists of two separate components, one for continuum emission and one for nuclear de-excitation lines at specific energies. Both components are determined over a broad energy range and from multiple years of data. Degradation effects of detectors and time dependent variance of background level are taken into account by modeling the background per detector on an orbital timescale. The mathematical description of our modeling method is given by which means that the data and model per energy bin k are represented by the sum of the celestial components i of the total number of N l celestial sources convolved through the image response function R per detector j and the sum of all background components N b of detector j. No prior knowledge concerning the energy spectrum of the celestial sources is assumed in our analysis. In general, the model is fit to the data by minimizing the Cash Statistic (Cash 1979), adjusting the scaling parameters θ i in Eq. (1), where different timescales for the scaling of the components are allowed. We use the spimodfit analysis tool (Strong et al. 2005;Halloin 2009), which applies a Levenberg-Marquardt algorithm to determine the maximum likelihood solution for all intensity parameters θ i . The software is based on the ISDC software spiros (Dubath et al. 2005), however optimized for high spectroscopic resolution of low signal to noise sources. Unless otherwise stated, uncertainties are given as 1σ. We use the Pearson χ 2 as an absolute goodness-of-fit criterion. We note that the chosen absolute goodness-of-fit criterion (χ 2 ! = 1.0) itself carries an uncertainty (Andrae et al. 2010).The smallest possible timescale we utilized is a single pointing. To minimize the contamination of our data set by known periods of increased background, we excluded orbit phases below 0.10 and above 0.88, during which the satellite passes through the Van Allen radiation belts.

Background modeling
SPI instrumental background mainly originates from the bombardment of the satellite by cosmic ray particles. Interaction of these cosmic rays can induce nuclear reactions in the satellite materials. Subsequent decays from excited nuclear levels lead to the emission of nuclear de-excitation lines, which fall into the energy range of the SPI detectors. Among others, bremsstrahlung is a second dominant contribution to the background, forming an underlying continuum. To determine the temporal and spectral behavior of the background, we used the knowledge gained from 17 yr of integrated mission data. Long-term temporal variation is introduced by the degradation of the lattice structure of the Germanium detectors and the absolute production rate of cosmic rays, which is anticorrelated with the solar cycle. Short-term variations are introduced by solar flares.
The detailed spectral shapes of the continuum emission and nuclear de-excitation lines are determined separately. Since a physically based model is difficult to construct and calibrate at the required precision, we used an empirical description of the background. This is based on previous attempts to model the highly variable instrumental background in SPI, as, for example, in Knödlseder et al. (2004) andJean et al. (2003, and supersedes the standard on-off methods as presented in Dubath et al. (2005). All background components are determined as a linear superposition of an underlying continuum normalized to a central pivot energy, superimposed by emission lines. The line shapes are represented by Gaussian functions convolved with a degradation function, which accounts for the degradation of the germanium charge collection efficiency (Kretschmer 2011). We determine the spectral shape on a 3 d (one orbit) period separately for each detector to trace the time-dependent degradation of the detectors. This timescale is chosen as the best compromise between accumulating sufficient statistics and appropriate determination of temporal variations of the spectral shape. Secondary contributions to the background are smeared out in our coded-mask analysis by accumulating data over multiple pointings (Siegert et al. 2019). The consistency of this high-resolution, time-dependent background modeling approach is demonstrated, for example, in Siegert et al. (2016); Siegert (2017), and Diehl et al. (2018).

Spectral analysis
With the modest spatial resolution of SPI of ≈ 3 • , supernova remnants cannot be resolved in separate clumps of ejecta. To enhance sensitivity for the relatively low-intensity total celestial signal in a line, a model for the line shape has to be adopted. We describe the emission produced by radioactive decay with Gaussian shaped lines, plus a power-law-shaped continuum.
where F 0 is the measured line flux, E 0 is the energy of the Doppler shifted-line centroid, and σ is line width. We interpret any broadening of the line, which would be additional to the detector resolution, as Doppler broadening caused by the expansion velocity of the ejecta. We determine the line parameters separately for each line when possible. No changes in the kinematics of the ejecta are expected, allowing for a combined three-line fit, assuming identical Doppler parameters for all three lines. We further allow for an underlying celestial continuum accounting for bremsstrahlung processes, with normalization parameter A 0 and power-law index α. The mass of the ejected 44 Ti per line is determined by where F L is the flux of the specific line L = (68, 78, 1157) keV, d is the distance to the source, N = 44 is the number of nuclei in 44 Ti, u is the atomic mass number, τ = 86.6 yr is the decay constant of 44 Ti, and t is the age of the supernova remnant. All line fluxes F 0 are normalized with the branching ratio b L 1 of the specific line L. For comparison, all flux values are stated as the normalized flux F L = F 0 /b L . In a multiline fit, the branching ratios determine the relative line intensities, and the normalized flux F L is fit as the parameter of interest. Velocities corresponding to the Doppler broadening of the lines are calculated for the ejecta from the FWHM of the line, assuming a uniformly expanding sphere. While this model might not adequately describe asymmetries as seen in the supernova remnant Cassiopeia A (Grefenstette et al. 2017), it provides a reasonable first-order approach for determining fluxes from the remaining unresolved sources. Confidence intervals for our results of the spectral fits are estimated from the 68th percentile interval of a Metropolis-Hastings algorithm, minimizing the Pearson χ 2 as test statistics. Upper limits are given at 2σ (i.e., ∆χ 2 = 4, for one degree of freedom, dof). We derived our upper limits by varying only the integrated flux of the respective line, assuming values for Doppler broadening and shift, and utilizing the best fit values for the underlying continuum.

Cassiopeia A
We used all available data for Cassiopeia A up to 2019 AD, containing a total exposure of 11.2 Ms for our analysis. Figure 2 shows the spectrum of Cassiopeia A in the energy ranges of interest. The average reduced χ 2 per fit energy bin is 1.001 (χ 2 /do f = 92,772/92,658). We adopted a uniform power law underlying the line emission across the entire energy range between 30 keV and 1200 keV with a fit power-law index of α = −2.6 ± 0.4. The line signal with the highest significance for a single Gaussian shaped line is found for the 78 keV line, with a significance of 3.6σ. The strong background lines of germanium lead to relatively large flux variations in the energy range between ≈ 50 − 65 keV, to such an extent that the 68 keV decay line is only marginally detectable. In order to validate our findings and to avoid spurious detection, we searched for an emission that could mimic 44 Ti decay from any celestial point source in a 20 • × 20 • area centered at Cassiopeia A. Possible emission was determined at locations in a spherical square-shaped grid, where each point is separated by 2 • from the adjacent point, yielding a total of 121 grid points. Figure 3 shows the map of source significances for a 78+1157 keV signal. Detection at 4.9σ is only found at the location consistent with Cassiopeia A, while other test points do not show significant signals. We additionally show (Fig. 3) the 1σ uncertainty band determined form these 121 source locations, which, if interpreted as statistical signal fluctuations of a zero signal, bracket the spectrum of Cassiopeia A. A significant excess above the uncertainty band is clearly visible at the 78 keV line, while the majority of the 68 keV line lies within the statistical backgrounduncertainty band. The 78 keV line, if represented by a Gaussian-shaped line, represents an integrated flux of (3.3 ± 0.9) × 10 −5 ph cm −2 s −1 and a 3.6σ detection level. With Eq. (3), this flux corresponds to an ejected 44 Ti mass of (2.1 ± 0.6) × 10 −4 M , for a remnant age of 330 yr and a distance of 3300±100 pc. The uncertainty is mainly due to the uncertainty in the determined flux, however, the uncertainty in the distance estimate is also incorporated in the result. The line is centered at (77.7 ± 0.5) keV, which is slightly red-shifted with respect to the laboratory-determined decay energy of 78.3 keV (Firestone et al. 2003). This Doppler shift translates into a bulk motion of (2400 ± 1500) km s −1 away from the observer. The full width at half maximum (FWHM) of the line is (2.3 ± 0.8) keV. This is broadened with respect to the instrumental resolution of 1.6 keV FWHM at 78 keV. We interpret this broadening of the line as Doppler broadening due to the expansion of the supernova remnant; this translates the FWHM of the line into an expansion velocity of (5500 ± 2700) km s −1 .
Representing the 1157 keV line with a simple Gaussian on top of the uniform power law, our best fit values suggest a significantly higher flux of (9.5 ± 3.0) × 10 −5 ph cm −2 s −1 . This flux can be overestimated due to the assumption of a single underlying continuum fit across the broad energy range. The continuum flux from our best fit power law (largely determined at energies below 100 keV) in the energy region between 1125 -1175 keV is 1.0±7.0×10 −7 ph cm −2 s −1 , consistent with zero. We estimated the potential offset in the fit continuum flux density by allowing a separate, constant offset in the energy range 1090 -1210 keV, which accounts for a flux of 3.1 ± 1.5 × 10 −5 ph cm −2 s −1 in the energy region between 1125 -1175 keV. Therefore, the line flux at 1157 keV can be reduced to (6.4±3.4)×10 −5 ph cm −2 s −1 . The Gaussian is centered at (1151 ± 7.9) keV, which is red shifted but compatible with the laboratory-determined energy of 1157 keV. The line width is (40.0 ± 6.7) keV FWHM. This corresponds to (8900 ± 1500) km s −1 expansion velocity. We find a combined signal with a significance of 4.9σ by simultaneously fitting two lines with identical Doppler shift, Doppler broadening, and integrated flux. The overall Doppler shift of the lines corresponds to a bulk ejecta velocity of (1800±800) km s −1 . The Doppler broadening for both lines translates to (6400 ± 1900) km s −1 expansion velocity, in agreement with values determined for the 78 keV line alone. Due to the higher relative spectral resolving power of SPI at higher energies, the expansion velocity can be better constrained including the 1157 keV line. The combined fit contains a flux of (4.2 ± 1.0) × 10 −5 ph cm −2 s −1 per line. This higher flux corresponds to a 44 Ti mass of (2.6 ± 0.6) × 10 −4 M . Table 2 contains measured line parameters and derived physical quantities. Even though the uncertainties are high, we also include the 68 keV line in our analysis. However, we adopted the kinematic values determined from the combined 78 and 1157 keV line fit for this line, as the strong fluctuations induced by the strong background lines might lead to an artificial broadening of the line. The 68 keV line is then observed with a single-line significance of 2.2σ. When linked to a common origin, the total significance for the three fit lines is then increased to 5.4σ.

SN 1987A
SPI was pointed towards the LMC including SN 1987A for a total of 7 Ms. With SPI's angular resolution of ≈ 2.7 • , we cannot distinguish between SN 1987A and other potential or known sources of high energy emission. In particular, the pulsar PSR B0540-69 and the high-mass X-ray binary LMC-X1 are located less than 1 • apart from SN 1987A. We obtain an average reduced χ 2 of 1.000 (χ 2 /do f = 51,716/51,731) per energy bin. In Fig. A.4, we show the spectrum obtained from our data for a source located at the position of SN 1987A. We find no significant flux excess in either energy region that could be attributed to the decay chain emission of 44 Ti. Tueller et al. (1990) determined an expansion velocity of 3100 km s −1 from the line profiles of measured radioactive 56 Co. Assuming co-moving 44 Ti ejecta, we adopted this value to determine 2σ upper limits on the flux. In our analysis, we searched for the combined signal of all three lines simultaneously, for which we determined a value of 1.8 × 10 −5 ph cm −2 s −1 per line, corresponding to an upper mass limit of ejected 44 Ti of 6.9 × 10 −4 M for a distance of 49.6 kpc, and a remnant age of 24 yr. No systematically increased flux is observed in the 1157 keV line with respect to the lines at 68 and 78 keV.

Vela Junior
We modeled the supernova remnant as a source of extended emission with a 2D Gaussian emission profile and a diameter of 0.6 • for the width of the remnant. This means that the 2 • diameter contains ≈ 90% of the expected 44 Ti signal. The obtained fit is satisfactory with a reduced χ 2 of 0.998 (χ 2 /do f = 65,423/65,535) per energy bin. We find no signal for the decay of 44 Ti (Fig. A.3). We determined a 2σ upper limit of 2.1 × 10 −5 ph cm −2 s −1 for the combined signal of all three lines, assuming no bulk motion and an expansion velocity of 3000 km s −1 . This corresponds to an upper limit for the ejected 44 Ti mass of 3.3 × 10 −5 M for the remnant age and distance of 690 yr and 200 pc, respectively. Considering the updated age and distance estimates of 2.4 -5.1 kyr and 700 ± 200 pc (Allen et al. 2014), the ejecta mass limits determined from our results significantly increase to a value ≤ 2.2 × 10 −1 M for the lower age limit of 2.4 kyr.

Tycho's supernova remnant
SPI was pointed towards the region containing Tycho's supernova remnant for a total of 10 Ms. Figures A.1   79,982/80,196) per energy bin. We find no significant excess for the emission in the three lines of the 44 Ti decay chain. To determine our upper limits, we adopted an expansion velocity of 5000 km s −1 . This value is in agreement with the expansion velocities found in the central ejecta (Sato & Hughes (2017). We determined a 2σ upper limit of 1.4 × 10 −5 ph cm −2 s −1 for each line in the 44 Ti decay, corresponding to an 44 Ti ejecta-mass limit of 4.8 × 10 −4 M for a distance of 4.1 kpc and a remnant age of 438 yr.

G1.9+0.3
Due to its location close to the Galactic center, several hard Xray sources (Bird et al. 2016  safe assumption that potential flux excess in the 68 and 78 keV regime can be attributed to the emission from G1.9+0.3. The average reduced χ 2 is 1.011 (χ 2 /do f = 318,356/315,005). We interpret this as due to the possible presence of unresolved sources in SPI's field of view in the Galactic central region. The resulting spectra (Figs. A.1 and A.2) show an underlying continuum from the spatially coincident sources, upon which we search for the imprints of the three decay lines. We find no significant excess in both energy ranges, determining a corresponding upper limit of 1.0 × 10 −5 ph cm −2 s −1 for an assumed expansion velocity of 5000 km s −1 for 44 Ti containing ejecta. This translates into a 44 Ti yield of 0.3 × 10 −4 M for a remnant age of 120 yr and a distance of 8.5 kpc. Even though Article number, page 7 of 15 A&A proofs: manuscript no. 44Ti_In_Young_SN_Remnants_arXiv  this velocity is lower than the expansion of the remnant's blast wave, we believe our assumption is plausible, as the distribution of the ejecta containing radioactive 44 Ti is uncertain anyway and may consist of clumps as seen for Cassiopeia A. We determined velocity-dependent limits, which depend on the expected line width for expansion velocities between 0 and 15000 km s −1 . We obtain limits in the range from 0.7 to 1.5 × 10 −5 ph cm −2 s −1 , assuming the same Doppler velocities for all three lines.

Kepler's supernova remnant
In our analysis, we do not find emission from Kepler's supernova remnant in the two energy bands (Figs. A.1 and A.2). The (2σ) upper limit is 1.1 × 10 −5 ph cm −2 s −1 . For a remnant age of 406 yr and distance of 5.1 kpc, the flux limit corresponds to a 44 Ti ejecta mass limit of 4.0 × 10 −4 M .

Cassiopeia A
For six analyzed supernova remnants, we find significant detection only for Cassiopeia A, with an integrated flux of (4.2 ± 1.0) × 10 −5 ph cm −2 s −1 corresponding to an 44 Ti ejecta mass of (2.6 ± 0.6) × 10 −4 M . Conventional models of core-collapse supernova explosions (Timmes et al. 1996;Magkotsios et al. 2010), including models specific to the progenitor evolution of Cassiopeia A (Young et al. 2006), suggest 44 Ti ejecta of less than 1.0×10 −4 M , significantly lower than the amount that we determine for Cassiopeia A. This underproduction of 44 Ti in models is also supported by measurements from other instruments, consistently showing higher ejecta mass for Cassiopeia A ((2.4 ± 0.9) × 10 −4 M Siegert et al. 2015) and ((1.5 ± 0.2) × 10 −4 M Grefenstette et al. 2017) than the values from models of ≤ 1.0 × 10 −4 M obtained for a 30 M star (Timmes et al. 1996;Limongi & Chieffi 2018). Figure 4 shows the expected yield of 44 Ti for various supernova scenarios. The green shaded area represents standard, mostly pistondriven 2 explosion models, in which nucleosynthesis is calculated by post-processing from the modeled thermodynamic evolution of the remnant. Harris et al. (2017) point out that including nucleosynthesis networks into the simulation, rather than post processing yields, can change the production of alpha nuclei, especially at intermediate mass range A = 36 -52, by an order of magnitude. The high 44 Ti mass seen in Cassiopeia A shows that a more detailed treatment of explosive nucleosynthesis appears necessary. The 3D model of Wongwathanarat et al. (2017), representing the special case of Cassiopeia A, suggests higher 44 Ti masses, especially considering the clumpy and asymmetric distribution in the supernova remnant. The measured expansion velocity (6400 ± 1900) km s −1 is compatible with models including Rayleigh-Taylor instabilities that lead to large-scale mixing of the inner ejecta with overlying stellar shells in type II-b supernova models (Nomoto et al. 1995). We determined a bulk motion of (2200 ± 1300) km s −1 . This suggests that the bulk of the ejecta is receding from the observer. Both the kinematics and the ejected mass of 44 Ti support the interpretation that Cassiopeia A is an asymmetric supernova explosion. While we determined kinematic constraints from a spectral analysis, other evidence for an asymmetric explosion is provided from the spatially resolved analysis of the remnant with the NuSTAR telescope. Grefenstette et al. (2017) have found that the majority of the 44 Ticontaining ejecta is expelled in a large solid angle, where the bulk of the ejecta moves away from the observer. Our measurements and the resulting velocity spread, determined from the Doppler broadening, is consistent with NuSTAR findings, which suggest a clumped nature of 44 Ti -containing ejecta (Grefenstette et al. 2017). Despite concurring kinematic constraints, we determine a higher integrated flux in the 44 Ti decay lines. Due to the different angular resolutions of the NuSTAR telescope (18" FWHM Harrison et al. 2013) and the SPI spectrometer (2.7 • FWHM Vedrenne et al. 2003), different spherical surface areas for the integration of the flux are considered in both analyses. Grefenstette et al. (2017) considered the flux of an integrated emission from a region of 120" radius centered on Cassiopeia A, containing all spatial points in which emission from 44 Ti decay is detected in their analysis, and they gave upper limits on regions outside the 120" radius. X-ray measurements both suggest a forward shock radius of 153" (Gotthelf et al. 2001) and the presence of iron at radii between 110" -170" (Willingale et al. 2002). Co-moving 44 Ti can be present at large radii extending as far outwards as the observed iron distribution. Flux from these unresolved regions contributes to the total flux of the 44 Ti emission. Within the SPI's angular resolution, the entire surface area of the remnant is included, constituting the increased integrated flux we measured in our analysis. The increased line-of-sight beam width in SPI of 2.7 • FWHM also includes areas outside of the supernova remnant. This also means that unresolved or previously unknown sources can contribute to the total flux in SPI measurements. For the first time, we also identified a very broad decay signature in the high-energy decay line at 1157 keV, which also reveals the kinematic evolution of the supernova remnant. We determine an expansion velocity of (8900 ± 1500) km s −1 and a line that is not significantly redshifted with (1600 ± 2000) km s −1 . This line contains an integrated flux of (9.5±3.0)×10 −5 ph cm −2 s −1 ((6.4±3.4)×10 −5 ph cm −2 s −1 ). A systematic offset of the highenergy line has been also observed with COMPTEL (Iyudin 1999) and with earlier SPI data (Siegert et al. 2015), however with less exposure on Cassiopeia A. In contrast to our updated values, Siegert et al. (2015) provided a more constraining expansion velocity and an overall different kinematic behavior of 44 Sc decay in comparison to the 78 keV line for a line centered at (1158.0 ± 3.6) keV. As shown by Grefenstette et al. (2017) 44 Ti is ejected in clumps in Cassiopeia A. Each clump translates into a separate peak in the energy range between 1130 -1180 keV, which blends into a broadened line. We described the entire emission by one Gaussian, which captures the overall expansion of the entire remnant. The line measured by Siegert et al. (2015) in the narrow energy window around 1157 keV only captures parts, or one separate ejecta clump, of the total emission in the 1157 keV regime. We find an increased flux in the 1157 keV line, that cannot be explained by systematic effects alone. We speculate that the flux included in this line could be enhanced for different reasons: 1) Excitation of the nuclear transition in 44 Ca, in addition to the decay of 44 Sc. Interaction with ambient material can lead to excitation of the nucleus. An excitation of the stable 44 Ca nucleus by cosmic rays in the shock region of the supernova envelope might thus contribute to the flux in the 1157 keV line. This mechanism would only influence the flux of the 44 Ca * line, as the half life of 44 Sc is too short for efficient cosmic-ray-induced excitation. More analysis, in particular of other candidate nuclear de-excitation lines, is required to support this hypothesis. The most promising approach would be the detection of the de-excitation lines at 4.4 MeV and 6.1 MeV, which are the most prominent de-excitation lines caused by cosmic ray interaction in the shock front (Summa et al. 2011). A first search for these lines shows that flux values as high as those suggested by Summa et al. (2011) can be excluded. 2) McKinnon et al. (2016) pointed out that two thirds of dust in the Milky Way-like galaxies can be produced by type II supernova events. The presence of dust grains composed of ejecta material in the vicinity of the supernova remnant or in the line of sight towards the supernova remnant can alter the observed flux ratios beyond the branching ratios. Attenuation coefficients for 68 and 78 keV photons are higher than for the 1157 keV line (Iyudin et al. 2019) for common dust grain compositions. Including correction for branching ratios, we derived the following flux ratio: This suggests that 16 -80% of the emission in the 78 keV line could be absorbed by dust, located between INTEGRAL and Cassiopeia A.
3) The assumption of a Gaussian-shaped line for the total emission does not correctly represent the ejecta kinematics as found in Cassiopeia A. This can artificially lead to an increased flux in the 1157 keV line.

SN 1987A
We determine an upper limit of 1.8 × 10 −5 ph cm −2 s −1 per decay line, assuming all lines share identical Doppler characteristics. We attribute this flux to the decay of 44 Ti, corresponding to an upper ejecta mass limit of 6.9 × 10 −4 M for a remnant age of 24 yr and a distance of 49.6 kpc. We find no evidence for a systematically increased flux in the high energy line at 1157 keV. Assuming that the 44 Ti ejecta are contained in the central region of the supernova, the expansion velocity of the 44 Ti ejecta should be lower than 1800 km s −1 (McCray 2017). In addition to an upper limit determined from the combined three lines, we give upper limits only for the 1157 keV line for expansion velocities corresponding to the interior of the supernova core (v exp ≤ 1800 km s −1 ). For this velocity range, we determine flux limits that range between (1.7 − 3.2) × 10 −5 ph cm −2 s −1 . The flux limit we derived for the 44 Ti decay chain combined fit is compatible with direct detection of 44 Ti, as found by IBIS/INTEGRAL and NuSTAR (Grebenev et al. 2012;Boggs et al. 2015). Both analyses suggest narrow line broadening, which is compatible with slowly expanding 44 Ti ejecta. Given the NuSTAR and IBIS/INTEGRAL fluxes, a significant offset in flux of the 1157 keV line, either from less efficient absorption at higher energies or from an additional excitation, as seen in Cassiopeia A, should be detectable, albeit with small significance, within the SPI's sensitivity.

Vela Junior
Vela Jr. still poses a mystery 20 yr after the serendipitous detection of gamma-ray emission in the 1157 keV 44 Ca line by COMPTEL (Iyudin et al. 1998). Our upper limit of 2.1 × 10 −5 ph cm −2 s −1 assumes an extended source for the combined signal in all three decay lines, excluding a signal at the level reported by Iyudin et al. (1998) from COMPTEL data. Emission at the COMPTEL level was also excluded from the nondetection of the scandium fluorescence line Slane et al. (2001). More recent studies with IBIS/INTEGRAL (Tsygankov et al. 2016) found no excess in the energy bands of the 68 and 78 keV lines, with an upper limit for the 44 Ti flux of 1.8 × 10 −5 ph cm −2 s −1 , which excludes the COMPTEL detection. Tsygankov et al. (2016) point out that they considered the remnant as a point-like source neglecting the apparent 2 • diameter of the remnant. Measurements of the radial displacement in the northern rim of Vela Jr. suggest that the remnant age is 2.4 − 5.1 kyr (Allen et al. 2014) at a distance of 0.5 − 1.0 kpc, which is in contrast with the 0.7 kyr age and 200 pc distance estimates discussed by Iyudin et al. (1998). Our results substantiate the higher age and larger distance.

Tycho's supernova remnant
Our upper limit of 1.4 × 10 −5 ph cm −2 s −1 for the Tycho supernova remnant is in agreement with the upper limits determined with NuSTAR and INTEGRAL/IBIS (Wang & Li 2014;Lopez et al. 2015) of ≥ 10 −5 ph cm −2 s −1 for moderate expansion velocities, with 44 Ti spatially distributed over the entire remnant. However, detection has been claimed from observations with Swift/BAT (Troja et al. 2014) at a flux of 1.3 × 10 −5 ph cm −2 s −1 and 1.4 × 10 −5 ph cm −2 s −1 for the 68 and 78 keV lines, respectively. For an assumed distance of 4.1 kpc, the upper limits correspond to an ejecta mass of the order of 10 −4 M of 44 Ti. Lopez et al. (2015) provided an ejecta-mass upper limit of 2.4×10 −4 M of 44 Ti for a distance of 2.3 kpc. In all cases, we can exclude the double detonation and Ca-rich models as explosion scenario for this type Ia explosion remnant. The results are, however, in agreement with delayed detonation models. This model may be favored as it best reproduces the measured X-ray spectra from Tycho (Badenes et al. 2006).

G1.9+0.3
We find no significant excess of 44 Ti line emission at the position of the supernova remnant G1.9+0.3. We determine an upper flux limit of 1.0 × 10 −5 ph cm −2 s −1 . This translates into a mass limit of 0.3 × 10 −4 M at a distance of 8.5 kpc and a remnant age of 120 yr, which excludes both double detonation and Ca-rich models for this candidate type Ia explosion. Predictions for classical delayed detonation models (Maeda et al. 2010) suggest less than 10 −5 M of 44 Ti, which would be in agreement with our results. We compare our results with the signal detected by Borkowski et al. (2010). Using Chandra data, they reported a line at 4.1 keV, which they attribute to a fluorescence transition in 44 Sc following the electron capture decay on 44 Ti. Their inferred mass of (1 − 7) × 10 −5 M of ejected 44 Ti translates into a flux of (0.3 − 1.9) × 10 −5 ph cm −2 s −1 in the decay lines at 68, 78, and 1157 keV. Our upper limits (0.7 to 1.5 × 10 −5 ph cm −2 s −1 ) are in conflict with the extrapolated fluxes ((0.3 − 1.9) × 10 −5 ph cm −2 s −1 ) expected from the 44 Sc fluorescence line for expansion velocity below ∼15000 km s −1 .
The discrepancy between the 4.1 keV fluorescence line and the decay emission at hard X-ray energies suggests that the fluorescence line is not necessarily produced from the decay of 44 Ti alone. Since the line is produced by the emission of a Kα photon, which is independent of the scandium isotope, other stable isotopes (e.g., 45 Sc ) co-produced in the supernova explosion can contribute to the fluorescence emission at 4.1 keV. Zoglauer et al. (2015) determined a 2σ upper limit for the flux in the 68 keV line of 1.5 × 10 −5 ph cm −2 s −1 for a non-shifted line with a 4 keV width (1σ), using measurements of the NuSTAR telescope. Results from the IBIS telescope (Tsygankov et al. 2016) yield a 3σ upper limit of 9 × 10 −6 ph cm −2 s −1 .
The derived mass limit disagrees with the double detonation and Ca-rich scenario; however, due to the uncertain distance to the remnant, ranging between 4.4 kpc and 5.9 kpc, the double detonation scenario cannot be explicitly excluded. Our result is in agreement with the nondetection of 44 Ti emission with the COMPTEL and IBIS instrument (Iyudin 1999;Dupraz et al. 1997;Tsygankov et al. 2016).

Implications for supernova models
The detection of 44 Ti in only one Galactic supernova remnant provides a striking and significant conclusion regarding supernovae in the Galaxy. With an average core-collapse supernova rate of ≈ 1 − 3 century −1 van den Bergh & Tammann 1991), and current understanding of nucleosynthesis in supernovae, five supernova remnants with a 44 Ti -decay line flux of more than 10 −5 ph cm −2 s −1 are expected to be visible in the Galaxy. Detection of a single remnant at the position of Cassiopeia A is unlikely, with a probability of less than 2.7% Dufour & Kaspi 2013). With the high 44 Ti mass measured in Cassiopeia A (and also SN 1987A), it is possible that the 44 Ca content in the Galaxy is produced by a few, rare, 44 Ti producing supernovae. It is possible that both Cassiopeia A and SN 1987A are the prototypes for asymmetric explosions producing a high ejecta 44 Ti mass, whereas the majority of core-collapse supernovae explode in a more symmetric scenario, producing less 44 Ti ejecta.
We use the solar abundance value [ 44 Ca/ 56 Fe] = 1.2 × 10 −3 (Anders & Grevesse 1989) of the 44 Ca to 56 Fe ratio, which are the end products of the 44 Ti and 56 Ni decay chain respectively, as another criterion to judge supernova model types, assuming that each of those supernova types would be the sole source of 44 Ca and 56 Fe as found in the Sun. Figure. 5 shows the 44 Ti to 56 Ni ratio of the candidate sources in our analysis (dot symbols), together with modeled values for several supernova types (star symbols), and the solar 44 Ca to 56 Fe ratio (red line).
For Vela Jr., we adopted a 56 Ni ejecta-mass estimate from the explosion model of a 25 M (Maeda & Nomoto 2003). We note that the values for Vela Jr. may deviate from the solar ratio due to the assumed distance and age (200 pc; 690 yr), which underestimates the updated values determined by Allen et al. (2014). The ejected 56 Ni mass for Cassiopeia A is inferred from near-infrared spectral analysis (Eriksen et al. 2009). For SN 1987A, we used the 56 Ni ejecta-mass estimate based on the early bolometric light curve by Woosley et al. (1987). For type Ia supernova remnant candidates, 56 Ni masses are obtained from Badenes et al. (2006, Tycho), Patnaude et al. (2012, Kepler), and Borkowski et al. (2013, G1.9+0.3), who inferred ejecta masses by comparing measured X-ray spectra with long-term simulated remnant models.
We infer from Fig. 5 that the measured 44 Ti-to-56 Ni ratios of Cassiopeia A and SN 1987A are plausibly consistent with asymmetric core-collapse supernova models being responsible for the solar abundance ratio. In contrast, the majority of type Ia supernova scenarios do not plausibly reproduce the measured solar 44 Ca to 56 Fe ratio. Type Ia supernovae appear to consistently produce a 56 Ni ejecta mass between 10 −1 − 10 0 M (Dhawan et al. 2016;Wang et al. 2008;Stritzinger et al. 2006). In model calculations, nucleosynthesis yields of intermediate mass elements (such as 44 Ti) are highly dependent on physical conditions during nuclear burning. In centrally ignited, pure-deflagration, and delayed-detonation scenarios, nucleosynthesis occurs mainly in a high-density regime, producing only a small amount of 44 Ti, while the majority of nuclear fuel is completely burned to iron group elements (Seitenzahl et al. 2013;Fink et al. 2014). On the other hand, nucleosynthesis in double-detonation supernova scenarios also occurs during the initial burning of the surface helium layer at lower densities, allowing for the production of large amounts of 44 Ti on the surface Sim et al. 2012). Pure-deflagration and delayed-detonation type Ia supernovae (green and red star symbols, Fig.5) therefore could not be major contributors to solar 44 Ca, while double detonation supernovae (blue star symbols in Fig.5

Summary
In this work, we searched for the signature of gamma rays produced in the decay chain of 44 Ti in the six young nearby supernova remnants Cassiopeia A, SN 1987A, Vela Jr., Tycho's supernova, Kepler's supernova, and G1.9+0.3. In Tab. 3, we list the mass estimates we derived from SPI/INTEGRAL data acquired over the entire mission duration of 17 yr. Table 3. Values for the masses and fluxes of the six young supernova remnants. Fluxes are determined for the 44 Sc and 44 Ti decay separately and also with a combined fit for the the most stringent constraints in our analysis. We also include results determined with the NuSTAR and Integral/IBIS telescope, separated by the dashed line. References: (1) (Grefenstette et al. 2017); (2) (Tsygankov et al. 2016); (3) (Boggs et al. 2015); (4) (two lines combined Grebenev et al. 2012); (5) (Lopez et al. 2015); (6) (Wang & Li 2014); and (7) (Zoglauer et al. 2015 (Eriksen et al. 2009;Woosley et al. 1987;Badenes et al. 2006;Patnaude et al. 2012;Borkowski et al. 2013), respectively. Vela Jr. 56 Ni ejecta mass is modeled from a 25 M star (Maeda & Nomoto 2003). We only detect emission in the supernova remnant Cassiopeia A. Inferred masses of more than 2−3×10 −4 M ejected 44 Ti exceed theoretical predictions. Upper limits determined for Vela Jr. exclude the detection of the 44 Sc decay line found with COMPTEL by Iyudin (1999). The detection of 44 Ti hard X-ray lines in SN 1987A (Grebenev et al. 2012;Boggs et al. 2015) cannot be confirmed. We exclude models predicting high yields of 44 Ti such as the double-detonation models and Ca-rich explosion (Waldman et al. 2011;Perets et al. 2010) (Iyudin 1999) are close to the Galactic north pole with latitude 60 • <b < 80 • and longitude −320 • < l < -270 • containing a total of 11.5 Ms exposure for the entire INTEGRAL mission time. This region contains only a few sources, making it ideal for testing the analysis approach in an "empty" celestial region. The sources can be outside the fully coded field of view for various observation times, reducing the exposure on specific locations. The image response function was not calculated for angles between the pointing direction and source location of more than 25 degree offset, and we assumed no source contribution in these cases. Figures B.1 and B.2 show the spectra for the galaxy cluster NGC4388 and "empty" space. With our background modeling approach, the strong background lines in the region between 50 and 68 keV, 90 keV, as well as 1115-1125 keV are adequately suppressed. The photon index of the continuum emission (γ = 1.740±0.074) of NGC4388 (l = 279.1 • ; b = 74.3 • ) is consistent in the spectral index with the measurements of 1.72 ± 0.05 (Beckmann et al. 2004). We found no spurious 44 Ti signatures in either location.
In both regions, fluctuations around the expected continuum and baseline are present, however, these fluctuations are on scales smaller than the instrumental energy resolution at the respective energies and compatible with statistical fluctuations. We allow for a constant offset to account for possible diffuse emission sources located within the 2.7 • angular resolution of SPI. We do not find significant flux excess for any of the 44 Ti decay lines. The 2σ upper-limit flux determined for an expansion velocity of 3000 km s −1 is 1.8 × 10 −5 ph cm −2 s −1 corresponding to 7.0 × 10 −4 M synthesized 44 Ti (green line). Red dashed line corresponds to the IBIS flux (Grebenev et al. 2012) as it would be seen in SPI data.