The death of massive stars - I. Observational constraints on the progenitors of Type II-P supernovae

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 m(min) = 8.5(-1.5)(+1) M-circle dot and the maximum mass for II-P progenitors is m(max) = 16.5 +/- 1.5 M-circle dot, assuming a Salpeter initial mass function holds for the progenitor population (in the range Gamma = -1.35(-0.7)(+0.3)). 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-circle dot and the minimum mass to produce a Wolf-Rayet star in single star evolution (between solar and LMC metallicity) is similarly 25-30 M-circle dot. The reason we have not detected any high-mass red supergiant progenitors above 17 M-circle dot is unclear, but we estimate that it is statistically significant at 2.4 sigma 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-circle dot 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 Ni-56 masses ejected in the SNe to the progenitor mass estimates and find that low-luminosity SNe with low Ni-56 production are most likely to arise from explosions of low-mass progenitors near the mass threshold that can produce a core-collapse.

Observational Constraints on the Catastrophic Disruption Rate of Small Main Belt Asteroids

We have calculated 90% confidence limits on the steady-state rate of catastrophic disruptions of main belt asteroids in terms of the absolute magnitude at which one catastrophic disruption occurs per year  as a function of the post-disruption increase in brightness (Δm) and subsequent brightness decay rate (τ  ). The confidence limits were calculated using the brightest unknown main belt asteroid (V=18.5) detected with the Pan-STARRS1 (Pan-STARRS1) telescope. We measured the Pan-STARRS1’s catastrophic disruption detection efficiency over a 453-day interval using the Pan-STARRS moving object processing system (MOPS) and a simple model for the catastrophic disruption event’s photometric behavior in a small aperture centered on the catastrophic disruption event. We then calculated the  contours in the ranges from  and  encompassing measured values from known cratering and disruption events and our model’s predictions. Our simplistic catastrophic disruption model suggests that  and  which would imply that H_0≳28—strongly inconsistent withH_{0,B2005=23.26±0.02} predicted by Bottke et al. (Bottke, W.F., Durda, D.D., Nesvorný, D., Jedicke, R., Morbidelli, A., Vokrouhlický, D., Levison, H.F. [2005]. Icarus, 179, 63–94.) using purely collisional models. However, if we assume that H_0=H0,B2005 our results constrain , inconsistent with our simplistic impact-generated catastrophic disruption model. We postulate that the solution to the discrepancy is that >99% of main belt catastrophic disruptions in the size range to which this study was sensitive (∼100 m) are not impact-generated, but are instead due to fainter rotational breakups, of which the recent discoveries of disrupted asteroids P/2013 P5 and P/2013 R3 are probable examples. We estimate that current and upcoming asteroid surveys may discover up to 10 catastrophic disruptions/year brighter than V=18.5.

Observational Constraints on the Progenitors of Core-Collapse Supernovae: The Case for Missing High-Mass Stars

Over the last 15 years, the supernova community has endeavoured to directly identify progenitor stars for core-collapse supernovae discovered in nearby galaxies. These precursors are often visible as resolved stars in high-resolution images from space-and ground-based telescopes. The discovery rate of progenitor stars is limited by the local supernova rate and the availability and depth of archive images of galaxies, with 18 detections of precursor objects and 27 upper limits. This review compiles these results (from 1999 to 2013) in a distance-limited sample and discusses the implications of the findings. The vast majority of the detections of progenitor stars are of type II-P, II-L, or IIb with one type Ib progenitor system detected and many more upper limits for progenitors of Ibc supernovae (14 in all). The data for these 45 supernovae progenitors illustrate a remarkable deficit of high-luminosity stars above an apparent limit of log L/L-circle dot similar or equal to 5.1 dex. For a typical Salpeter initial mass function, one would expect to have found 13 high-luminosity and high-mass progenitors by now. There is, possibly, only one object in this time-and volume-limited sample that is unambiguously high-mass (the progenitor of SN2009ip) although the nature of that supernovae is still debated. The possible biases due to the influence of circumstellar dust, the luminosity analysis, and sample selection methods are reviewed. It does not appear likely that these can explain the missing high-mass progenitor stars. This review concludes that the community's work to date shows that the observed populations of supernovae in the local Universe are not, on the whole, produced by high-mass (M greater than or similar to 18 M-circle dot) stars. Theoretical explosions of model stars also predict that black hole formation and failed supernovae tend to occur above an initial mass of M similar or equal to 18 M-circle dot. The models also suggest there is no simple single mass division for neutron star or black-hole formation and that there are islands of explodability for stars in the 8-120 M-circle dot range. The observational constraints are quite consistent with the bulk of stars above M similar or equal to 18 M-circle dot collapsing to form black holes with no visible supernovae.

The death of massive stars - II. Observational constraints on the progenitors of Type Ibc supernovae

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.

Lucky Imaging of transiting planet host stars with LuckyCam

We obtained high-resolution, high-contrast optical imaging in the Sloan Digital Sky Survey i′ band with the LuckyCam camera mounted on the 2.56 m Nordic Optical Telescope, to search for faint stellar companions to 16 stars harbouring transiting exoplanets. The Lucky imaging technique uses very short exposures to obtain near diffraction-limited images yielding sub-arcsecond sensitivity, allowing us to search for faint stellar companions within the seeing disc of the primary planet host. Here, we report the detection of two candidate stellar companions to the planet host TrES-1 at separations <6.5 arcsec and we confirm stellar companions to CoRoT-2, CoRoT-3, TrES-2, TrES-4 and HAT-P-7 already known in the literature. We do not confirm the candidate companions to HAT-P-8 found via Lucky imaging by Bergfors et al., however, most probably because HAT-P-8 was observed in poor seeing conditions. Our detection sensitivity limits allow us to place constraints on the spectral types and masses of the putative bound companions to the planet host stars in our sample. If bound, the stellar companions identified in this work would provide stringent observational constraints to models of planet formation and evolution. In addition, these companions could affect the derived physical properties of the exoplanets in these systems.

Observations of ubiquitous compressive waves in the Sun's chromosphere

The details of the mechanism(s) responsible for the observed heating and dynamics of the solar atmosphere still remain a mystery. Magnetohydrodynamic waves are thought to have a vital role in this process. Although it has been shown that incompressible waves are ubiquitous in off-limb solar atmospheric observations, their energy cannot be readily dissipated. Here we provide, for the first time, on-disk observation and identification of concurrent magnetohydrodynamic wave modes, both compressible and incompressible, in the solar chromosphere. The observed ubiquity and estimated energy flux associated with the detected magnetohydrodynamic waves suggest the chromosphere is a vast reservoir of wave energy with the potential to meet chromospheric and coronal heating requirements. We are also able to propose an upper bound on the flux of the observed wave energy that is able to reach the corona based on observational constraints, which has important implications for the suggested mechanism(s) for quiescent coronal heating. © 2012 Macmillan Publishers Limited. All rights reserved.

Glycolaldehyde Formation via the Dimerisation of the Formyl Radical

Glycolaldehyde, the simplest monosaccharide sugar, has recently been detected in low- and high-mass star-forming cores. Following our previous investigation into glycolaldehyde formation, we now consider a further mechanism for the formation of glycolaldehyde that involves the dimerization of the formyl radical, HCO. Quantum mechanical investigation of the HCO dimerization process upon an ice surface is predicted to be barrierless and therefore fast. In an astrophysical context, we show that this mechanism can be very efficient in star-forming cores. It is limited by the availability of the formyl radical, but models suggest that only very small amounts of CO are required to be converted to HCO to meet the observational constraints.

SN 2009ib: a Type II-P supernova with an unusually long plateau

We present optical and near-infrared photometry and spectroscopy of SN 2009ib, a Type II-P supernova in NGC 1559. This object has moderate brightness, similar to those of the intermediate-luminosity SNe 2008in and 2009N. Its plateau phase is unusually long, lasting for about 130 d after explosion. The spectra are similar to those of the subluminous SN 2002gd, with moderate expansion velocities. We estimate the Ni-56 mass produced as 0.046 +/- A 0.015 M-aS (TM). We determine the distance to SN 2009ib using both the expanding photosphere method (EPM) and the standard candle method. We also apply EPM to SN 1986L, a Type II-P SN that exploded in the same galaxy. Combining the results of different methods, we conclude the distance to NGC 1559 as D = 19.8 +/- A 3.0 Mpc. We examine archival, pre-explosion images of the field taken with the Hubble Space Telescope, and find a faint source at the position of the SN, which has a yellow colour [(V - I)(0) = 0.85 mag]. Assuming it is a single star, we estimate its initial mass as M-ZAMS = 20 M-aS (TM). We also examine the possibility, that instead of the yellow source the progenitor of SN 2009ib is a red supergiant star too faint to be detected. In this case, we estimate the upper limit for the initial zero-age main sequence (ZAMS) mass of the progenitor to be similar to 14-17 M-aS (TM). In addition, we infer the physical properties of the progenitor at the explosion via hydrodynamical modelling of the observables, and estimate the total energy as similar to 0.55 x 10(51) erg, the pre-explosion radius as similar to 400 R-aS (TM), and the ejected envelope mass as similar to 15 M-aS (TM), which implies that the mass of the progenitor before explosion was similar to 16.5-17 M-aS (TM).

A simple approach to the supernova progenitor–explosion connection

We present a new approach to understand the landscape of supernova explosion energies, ejected nickel masses, and neutron star birth masses. In contrast to other recent parametric approaches, our model predicts the properties of neutrino-driven explosions based on the pre-collapse stellar structure without the need for hydrodynamic simulations. The model is based on physically motivated scaling laws and simple differential equations describing the shock propagation, the contraction of the neutron star, the neutrino emission, the heating conditions, and the explosion energetics. Using model parameters compatible with multi-D simulations and a fine grid of thousands of supernova progenitors, we obtain a variegated landscape of neutron star and black hole formation similar to other parametrized approaches and find good agreement with semi-empirical measures for the ‘explodability’ of massive stars. Our predicted explosion properties largely conform to observed correlations between the nickel mass and explosion energy. Accounting for the coexistence of outflows and downflows during the explosion phase, we naturally obtain a positive correlation between explosion energy and ejecta mass. These correlations are relatively robust against parameter variations, but our results suggest that there is considerable leeway in parametric models to widen or narrow the mass ranges for black hole and neutron star formation and to scale explosion energies up or down. Our model is currently limited to an all-or-nothing treatment of fallback and there remain some minor discrepancies between model predictions and observational constraints.

Time-resolved optical spectroscopy of the pulsating DA white dwarf HS0507+0434B - New constraints on mode identification and pulsation properties

We present a detailed analysis of time-resolved optical spectra of the ZZ Ceti white dwarf, HS 0507+0434B. Using the wavelength dependence of observed mode amplitudes, we deduce the spherical degree, l, of the modes, most of which have l = 1. The presence of a large number of combination frequencies (linear sums or differences of the real modes) enabled us not only to test theoretical predictions but also to indirectly infer spherical and azimuthal degrees of real modes that had no observed splittings. In addition to the above, we measure line-of-sight velocities from our spectra. We find only marginal evidence for periodic modulation associated with the pulsation modes: at the frequency of the strongest mode in the lightcurve, we measure an amplitude of 2.6 +/- 1.0 kms(-1), which has a probability of 2% of being due to chance; for the other modes, we find lower values. Our velocity amplitudes and upper limits are smaller by a factor of two compared to the amplitudes found in ZZ Psc. We find that this is consistent with expectations based on the position of HS 0507+0434B in the instability strip. Combining all the available information from data such as ours is a first step towards constraining atmospheric properties in a convectionally unstable environment from an observational perspective.