Using the Binary Population and Spectral Synthesis code, bpass, we have calculated the rates, time-scales and mass distributions for binary black hole (BH) mergers as a function of metallicity. We ...consider these in the context of the recently reported first Laser Interferometer Gravitational-Wave Observatory (LIGO) event detection. We find that the event has a very low probability of arising from a stellar population with initial metallicity mass fraction above Z = 0.010 (Z ≳ 0.5 Z⊙). Binary BH merger events with the reported masses are most likely in populations below 0.008 (Z ≲ 0.4 Z⊙). Events of this kind can occur at all stellar population ages from 3 Myr up to the age of the Universe, but constitute only 0.1–0.4 per cent of binary BH mergers between metallicities of Z = 0.001 and 0.008. However at metallicity Z = 10−4, 26 per cent of binary BH mergers would be expected to have the reported masses. At this metallicity, the progenitor merger times can be close to ≈10 Gyr and rotationally mixed stars evolving through quasi-homogeneous evolution, due to mass transfer in a binary, dominate the rate. The masses inferred for the BHs in the binary progenitor of GW 150914 are amongst the most massive expected at anything but the lowest metallicities in our models. We discuss the implications of our analysis for the electromagnetic follow-up of future LIGO event detections.
Literature data are collated for 38 stripped-envelope core-collapse supernovae (SE SNe; i.e. SNe IIb, Ib, Ic and Ic-BL) that have good light-curve coverage in more than one optical band. Using ...bolometric corrections derived in previous work, the bolometric light curve of each SN is recovered and template bolometric light curves provided. Peak light distributions and decay rates are investigated; SNe subtypes are not cleanly distinguished in this parameter space, although some grouping of types does occur and there is a suggestion of a Phillips-like relation for most SNe Ic-BL. The bolometric light curves are modelled with a simple analytical prescription and compared to results from more detailed modelling. Distributions of the explosion parameters show the extreme nature of SNe Ic-BL in terms of their 56Ni mass and the kinetic energy, however ejected masses are similar to other subtypes. SNe Ib and Ic have very similar distributions of explosion parameters, indicating a similarity in progenitors. SNe IIb are the most homogeneous subtype and have the lowest average values for 56Ni mass, ejected mass, and kinetic energy. Ejecta masses for each subtype and SE SNe as a whole are inconsistent with those expected from very massive stars. The majority of the ejecta mass distribution is well described by more moderately massive progenitors in binaries, indicating these are the dominant progenitor channel for SE SNe.
The relationship between stellar populations and the ionizing flux with which they irradiate their surroundings has profound implications for the evolution of the intergalactic medium (IGM). We ...quantify the ionizing flux arising from synthetic stellar populations which incorporate the evolution of interacting binary stars. We determine that these show ionizing flux boosted by 60 per cent at 0.05 ≤ Z ≤ 0.3 Z⊙ and a more modest 10–20 per cent at near-solar metallicities relative to star-forming populations in which stars evolve in isolation. The relation of ionizing flux to observables such as 1500 Å continuum and ultraviolet spectral slope is sensitive to attributes of the stellar population including age, star formation history and initial mass function (IMF). For a galaxy forming 1 M⊙ yr−1, observed at >100 Myr after the onset of star formation, we predict a production rate of photons capable of ionizing hydrogen, N
ion = 1.4 × 1053 s−1 at Z = Z⊙ and 3.5 × 1053 s−1 at 0.1 Z⊙, assuming a Salpeter-like IMF. We evaluate the impact of these issues on the ionization of the IGM, finding that the known galaxy populations can maintain the ionization state of the Universe back to z ∼ 9, assuming that their luminosity functions continue to M
UV = −10, and that constraints on the IGM at z ∼ 2–5 can be satisfied with modest Lyman-continuum photon escape fractions of 4–24 per cent depending on assumed metallicity.
We present the results of a 10.5-yr, volume-limited (28-Mpc) search for supernova (SN) progenitor stars. In doing so we compile all SNe discovered within this volume (132, of which 27 per cent are ...Type Ia) and determine the relative rates of each subtype from literature studies. The core-collapse SNe break down into 59 per cent II-P and 29 per cent Ib/c, with the remainder being IIb (5 per cent), IIn (4 per cent) and II-L (3 per cent). There have been 20 II-P SNe with high-quality optical or near-infrared pre-explosion images that allow a meaningful search for the progenitor stars. In five cases they are clearly red supergiants, one case is unconstrained, two fall on compact coeval star clusters and the other twelve have no progenitor detected. We review and update all the available data for the host galaxies and SN environments (distance, metallicity and extinction) and determine masses and upper mass estimates for these 20 progenitor stars using the stars stellar evolutionary code and a single consistent homogeneous method. A maximum likelihood calculation suggests that the minimum stellar mass for a Type II-P to form is mmin= 8.5+1−1.5 M⊙ and the maximum mass for II-P progenitors is mmax= 16.5 ± 1.5 M⊙, assuming a Salpeter initial mass function holds for the progenitor population (in the range Γ=−1.35+0.3−0.7). The minimum mass is consistent with current estimates for the upper limit to white dwarf progenitor masses, but the maximum mass does not appear consistent with massive star populations in Local Group galaxies. Red supergiants in the Local Group have masses up to 25 M⊙ and the minimum mass to produce a Wolf–Rayet star in single star evolution (between solar and LMC metallicity) is similarly 25–30 M⊙. The reason we have not detected any high-mass red supergiant progenitors above 17 M⊙ is unclear, but we estimate that it is statistically significant at 2.4σ confidence. Two simple reasons for this could be that we have systematically underestimated the progenitor masses due to dust extinction or that stars between 17–25 M⊙ produce other kinds of SNe which are not II-P. We discuss these possibilities and find that neither provides a satisfactory solution. We term this discrepancy the ‘red supergiant problem’ and speculate that these stars could have core masses high enough to form black holes and SNe which are too faint to have been detected. We compare the 56Ni masses ejected in the SNe to the progenitor mass estimates and find that low-luminosity SNe with low 56Ni production are most likely to arise from explosions of low-mass progenitors near the mass threshold that can produce a core-collapse.
The progenitors of many Type II core-collapse supernovae (SNe) have now been identified directly on pre-discovery imaging. Here, we present an extensive search for the progenitors of Type Ibc SNe in ...all available pre-discovery imaging since 1998. There are 12 Type Ibc SNe with no detections of progenitors in either deep ground-based or Hubble Space Telescope archival imaging. The deepest absolute BVR magnitude limits are between −4 and − 5 mag. We compare these limits with the observed Wolf-Rayet population in the Large Magellanic Cloud and estimate a 16 per cent probability that we have failed to detect such a progenitor by chance. Alternatively, the progenitors evolve significantly before core-collapse or we have underestimated the extinction towards the progenitors. Reviewing the relative rates and ejecta mass estimates from light-curve modelling of Ibc SNe, we find both incompatible with Wolf-Rayet stars with initial masses >25 M being the only progenitors. We present binary evolution models that fit these observational constraints. Stars in binaries with initial masses 20 M lose their hydrogen envelopes in binary interactions to become low-mass helium stars. They retain a low-mass hydrogen envelope until 104 yr before core-collapse; hence, it is not surprising that Galactic analogues have been difficult to identify.
We propose a simple model to explain the velocity of young neutron stars. We attempt to confirm a relationship between the amount of mass ejected in the formation of the neutron star and the ‘kick’ ...velocity imparted to the compact remnant resulting from the process. We assume that the velocity is given by v
kick = α (M
ejecta/M
remnant) + β . To test this simple relationship, we use the bpass (Binary Population and Spectral Synthesis) code to create stellar population models from both single and binary star evolutionary pathways. We then use our Remnant Ejecta and Progenitor Explosion Relationship (reaper) code to apply different α and β values, and three different ‘kick’ orientations then record the resulting velocity probability distributions. We find that while a single star population provides a poor fit to the observational data, the binary population provides an excellent fit. Values of α = 70 km s−1 and β = 110 km s−1 reproduce the Hobbs et al. observed two-dimensional velocities, and α = 70 km s−1 and β = 120 km s−1 reproduce their inferred three-dimensional velocity distribution for nearby single neutron stars with ages less than 3 Myr. After testing isotropic, spin-axis aligned and orthogonal to spin-axis ‘kick’ orientations, we find no statistical preference for a ‘kick’ orientation. While ejecta mass cannot be the only factor that determines the velocity of supernova compact remnants, we suggest that it is a significant contributor and that the ejecta-based ‘kick’ should replace the Maxwell–Boltzmann velocity distribution currently used in many population synthesis codes.
We report and discuss post-explosion observations of supernova iPTF13bvn. We find that the brightness of the supernova (SN) at +740 d is below the level of the pre-explosion source and thus confirm ...that the progenitor has exploded. We estimate that the late-time brightness is still dominated by the SN, which constrains the magnitude and thus mass of a possible companion star to below approximately 10 M⊙. In turn, this implies that the progenitor's initial mass is constrained to a narrow range between 10 and 12 M⊙. The progenitor of iPTF13bvn would have been a helium giant rather than a Wolf–Rayet star. In addition, we suggest that sufficiently deep observations acquired in 2016 would now stand a chance to directly observe the companion star.
Re-evaluating old stellar populations Stanway, E R; Eldridge, J J
Monthly Notices of the Royal Astronomical Society,
09/2018, Volume:
479, Issue:
1
Journal Article
ABSTRACT We compare the impacts of uncertainties in both binary population synthesis models and the cosmic star formation history on the predicted rates of gravitational wave (GW) compact binary ...merger events. These uncertainties cause the predicted rates of GW events to vary by up to an order of magnitude. Varying the volume-averaged star formation rate density history of the Universe causes the weakest change to our predictions, while varying the metallicity evolution has the strongest effect. Double neutron star merger rates are more sensitive to assumed neutron star kick velocity than the cosmic star formation history. Varying certain parameters affects merger rates in different ways depending on the mass of the merging compact objects; thus some of the degeneracy may be broken by looking at all the event rates rather than restricting ourselves to one class of mergers.