No Arabic abstract
The physics of core-collapse (CC) supernovae (SNe) and how the explosions depend on progenitor properties are central questions in astronomy. For only a handful of SNe, the progenitor star has been identified in pre-explosion images. Supernova remnants (SNRs), which are observed long after the original SN event, provide a unique opportunity to increase the number of progenitor measurements. Here, we systematically examine the stellar populations in the vicinities of 23 known SNRs in the Small Magellanic Cloud (SMC) using the star formation history (SFH) maps of Harris & Zaritsky (2004). We combine the results with constraints on the SNR metal abundances and environment from X-ray and optical observations. We find that 22 SNRs in the SMC have local SFHs and properties consistent with a CC explosion, several of which are likely to have been high-mass progenitors. This result supports recent theoretical findings that high-mass progenitors can produce successful explosions. We estimate the mass distribution of the CC progenitors and find that this distribution is similar to a Salpeter IMF (within the uncertainties), while this result is shallower than the mass distribution found in M31 and M33 by Jennings et al. (2014) and D{i}az-Rodr{i}guez et al. (2018) using a similar approach. Additionally, we find that a number of the SMC SNRs exhibit a burst of star formation between 50-200 Myr ago. As these sources are likely CC, this signature may be indicative of massive stars undergoing delayed CC as a consequence of binary interaction, rapid rotation, or low metallicity. In addition, the lack of Type Ia SNRs in the SMC is possibly a result of the short visibility times of these sources as they may fall below the sensitivity limits of current radio observations.
We model the present day, observable, normal radio pulsar population of the Small Magellanic Cloud (SMC). The pulsars are generated with SeBa, a binary population synthesis code that evolves binaries and the constituent stellar objects up to remnant formation and beyond. We define radio pulsars by selecting neutron stars that satisfy a selection of criteria defined by Galactic pulsars, and apply the detection thresholds of previous and future SMC pulsar surveys.The number of synthesised and recovered pulsars are exceptionally sensitive to the assumed star formation history and applied radio luminosity model, but is not affected extensively by the assumed common envelope model, metallicity, and neutron star kick velocity distribution. We estimate that the SMC formed (1.6$pm$0.3)$times 10^4$ normal pulsars during the last 100 Myrs. We study which pulsars could have been observed by the Parkes multibeam survey of the SMC, by applying the surveys specific selection effects, and recover 4.0$pm$0.8 synthetic pulsars.This is in agreement with their five observed pulsars. We also apply a proposed MeerKAT configuration for the upcoming SMC survey, and predict that the MeerKAT survey will detect 17.2$pm$2.5 pulsars.
We study the spatial correlations between the H$alpha$ emission and different types of massive stars in two local galaxies, the Large Magellanic Cloud (LMC) and Messier 33. We compare these to correlations derived for core-collapse supernovae (CCSNe) in the literature to connect CCSNe of different types with the initial masses of their progenitors and to test the validity of progenitor mass estimates which use the pixel statistics method. We obtain samples of evolved massive stars in both galaxies from catalogues with good spatial coverage and/or completeness, and combine them with coordinates of main-sequence stars in the LMC from the SIMBAD database. We calculate the spatial correlation of stars of different classes and spectral types with H$alpha$ emission. We also investigate the effects of distance, noise and positional errors on the pixel statistics method. A higher correlation with H$alpha$ emission is found to correspond to a shorter stellar lifespan, and we conclude that the method can be used as an indicator of the ages, and therefore initial masses, of SN progenitors. We find that the spatial distributions of type II-P SNe and red supergiants of appropriate initial mass ($gtrsim$9 $M_{odot}$) are consistent with each other. We also find the distributions of type Ic SNe and WN stars with initial masses $gtrsim$20 $M_{odot}$ consistent, while supergiants with initial masses around 15 $M_{odot}$ are a better match for type IIb and II-L SNe. The type Ib distribution corresponds to the same stellar types as type II-P, which suggests an origin in interacting binaries. On the other hand, we find that luminous blue variable stars show a much stronger correlation with H$alpha$ emission than do type IIn SNe.
What cosmic ray ionisation rate is required such that a non-ideal magnetohydrodynamics (MHD) simulation of a collapsing molecular cloud will follow the same evolutionary path as an ideal MHD simulation or as a purely hydrodynamics simulation? To investigate this question, we perform three-dimensional smoothed particle non-ideal magnetohydrodynamics simulations of the gravitational collapse of rotating, one solar mass, magnetised molecular cloud cores, that include Ohmic resistivity, ambipolar diffusion, and the Hall effect. We assume a uniform grain size of $a_text{g} = 0.1mu$m, and our free parameter is the cosmic ray ionisation rate, $zeta_text{cr}$. We evolve our models, where possible, until they have produced a first hydrostatic core. Models with $zeta_text{cr}gtrsim10^{-13}$ s$^{-1}$ are indistinguishable from ideal MHD models and the evolution of the model with $zeta_text{cr}=10^{-14}$ s$^{-1}$ matches the evolution of the ideal MHD model within one per cent when considering maximum density, magnetic energy, and maximum magnetic field strength as a function of time; these results are independent of $a_text{g}$. Models with very low ionisation rates ($zeta_text{cr}lesssim10^{-24}$ s$^{-1}$) are required to approach hydrodynamical collapse, and even lower ionisation rates may be required for larger $a_text{g}$. Thus, it is possible to reproduce ideal MHD and purely hydrodynamical collapses using non-ideal MHD given an appropriate cosmic ray ionisation rate. However, realistic cosmic ray ionisation rates approach neither limit, thus non-ideal MHD cannot be neglected in star formation simulations.
Most massive stars exchange mass with a companion, leading to evolution which is altered drastically from that expected of stars in isolation. Such systems are the result of unusual binary evolution pathways and, as such, may be used to place stringent constraints on the physics of these interactions. We use the R4 systems B[e] supergiant, which has been postulated to be the product of a binary stellar merger, to guide our understanding of such outcomes by comparing observations of R4 to the results of simulations of mergers performed with the 3d hydrodynamics code FLASH. Our approach tailors the simulation initial conditions to the observed properties of R4 and implements realistic stellar profiles generated by the 1d stellar evolution code MESA onto the 3d grid, resolving the merger inspiral to within $0.02, R_{odot}$. We then map the merger remnant into MESA to track its evolution on the HR diagram over a period of $10^4$ years. This generates models for a B[e] supergiant with stellar properties, age, and nebula structure in qualitative agreement with that of the R4 system. Our calculations provide concrete evidence to support the idea that R4 was originally a member of a triple system in which the inner binary merged after its most massive member evolved off the main sequence, producing a new object that is of similar mass yet significantly more luminous than the A supergiant companion. The potential applications of the code framework presented in this paper are wide ranging and can be used to generate models of a variety of merger stellar remnants.
We infer the progenitor mass distribution for 22 historic core-collapse supernovae (CCSNe) using a Bayesian hierarchical model. For this inference, we use the local star formation histories to estimate the age for each supernova (SN). These star formation histories often show multiple bursts of star formation; our model assumes that one burst is associated with the SN progenitor and the others are random bursts of star formation. The primary inference is the progenitor age distribution. Due to the limited number of historic SNe and highly uncertain star formation at young ages, we restrict our inference to the slope of the age distribution and the maximum age for CCSNe. Using single-star evolutionary models, we transform the progenitor age distribution into a progenitor mass distribution. Under these assumptions, the minimum mass for CCSNe is ${M_textrm{min}}~=~8.60^{+0.37}_{-0.41} M_odot$ and the slope of the progenitor mass distribution is $alpha = -2.61^{+1.05}_{-1.18}$. The power-law slope for the progenitor mass distribution is consistent with the standard Salpeter initial mass function ($alpha = -2.35$). These values are consistent with previous estimates using precursor imaging and the age-dating technique, further confirming that using stellar populations around SN and supernova remnants is a reliable way to infer the progenitor masses.