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.

Massive stars exploding in a He-rich circumstellar medium - VI. Observations of two distant Type Ibn supernova candidates discovered by La Silla-QUEST

We present optical observations of the peculiar stripped-envelope supernovae (SNe) LSQ12btw and LSQ13ccw discovered by the La Silla-QUEST survey. LSQ12btw reaches an absolute peak magnitude of M-g = -19.3 +/- 0.2, and shows an asymmetric light curve. Stringent pre-discovery limits constrain its rise time to maximum light to less than 4 d, with a slower post-peak luminosity decline, similar to that experienced by the prototypical SN Ibn 2006jc. LSQ13ccw is somewhat different: while it also exhibits a very fast rise to maximum, it reaches a fainter absolute peak magnitude (M-g =-18.4 +/- 0.2), and experiences an extremely rapid post-peak decline similar to that observed in the peculiar SN Ib 2002bj. A stringent pre-discovery limit and an early marginal detection of LSQ13ccw allow us to determine the explosion time with an uncertainty of +/- 1 d. The spectra of LSQ12btw show the typical narrow He I emission lines characterizing Type Ibn SNe, suggesting that the SN ejecta are interacting with He-rich circumstellar material. The He I lines in the spectra of LSQ13ccw exhibit weak narrow emissions superposed on broad components. An unresolved H alpha line is also detected, suggesting a tentative Type Ibn/IIn classification. As for other SNe Ibn, we argue that LSQ12btw and LSQ13ccw likely result from the explosions of Wolf-Rayet stars that experienced instability phases prior to core collapse. We inspect the host galaxies of SNe Ibn, and we show that all of them but one are hosted in spiral galaxies, likely in environments spanning a wide metallicity range.

Early-type stars observed in the ESO UVES Paranal observatory project - IV. studies of CN, CH+ and CH in the interstellar medium

High spectral resolution (~80 000) and signal-to-noise observations from the Ultraviolet and Visual Echelle Spectrograph Paranal Observatory Project (UVES-POP) are used to study the interstellarmolecular lines CN (3874 Å), CH⁺ (3957, 4232 Å) and CH (3886, 4300 Å) towards 74 O- and B-type stellar sightlines. Additionally, archive data are presented for 140 ELODIE early-type stellar sightlines at R = 42 000, plus 25 FEROS at R = 48 000 and 3 UVES at R > 50 000, mainly in the CH⁺ (4232 Å) and CH (3886, 4300 Å) transitions. Detection rates are ~45 per cent for CN and ~67 per cent for the other lines in the POP sample, and ~10-15 per cent for CH⁺ and CH lines in the additional sample. CH and CH⁺ are well correlated between log[N(CH) cm^-2]~12-14, implying that these clouds are CH⁺-like CH and not CN-like CH. CH is also very well correlated with Na I D in the range log[N(Na I cm^-2]) ~12.2-14.2. A few sightlines show tentative velocity shifts of ~2 km s^-1 between CH and CH⁺, which appear to be caused by differences in component strength in blends, and hence do not provide firm evidence for shocks. Finally, we describe a search for ¹³CH⁺ in a sightline towards HD 76341. No ¹³CH⁺ is detected, placing a limit on the ¹³CH⁺ to ¹²CH⁺ ratio of ~0.01. If a formal fit is attempted, the equivalent width ratio in the two isotopes is a factor ~90 but with large errors. 

Early-type stars observed in the ESO UVES Paranal Observatory Project - V. Time-variable interstellar absorption

The structure and properties of the diffuse interstellar medium (ISM) on small scales, sub-au to 1 pc, are poorly understood. We compare interstellar absorption-lines, observed towards a selection of O- and B-type stars at two or more epochs, to search for variations over time caused by the transverse motion of each star combined with changes in the structure in the foreground ISM. Two sets of data were used: 83 VLT- UVES spectra with approximately 6 yr between epochs and 21 McDonald observatory 2.7m telescope echelle spectra with 6 - 20 yr between epochs, over a range of scales from 0 - 360 au. The interstellar absorption-lines observed at the two epochs were subtracted and searched for any residuals due to changes in the foreground ISM. Of the 104 sightlines investigated with typically five or more components in Na I D, possible temporal variation was identified in five UVES spectra (six components), in Ca II, Ca I and/or Na I absorption-lines. The variations detected range from 7\% to a factor of 3.6 in column density. No variation was found in any other interstellar species. Most sightlines show no variation, with 3{\sigma} upper limits to changes of the order 0.1 - 0.3 dex in Ca II and Na I. These variations observed imply that fine-scale structure is present in the ISM, but at the resolution available in this study, is not very common at visible wavelengths. A determination of the electron densities and lower limits to the total number density of a sample of the sightlines implies that there is no striking difference between these parameters in sightlines with, and sightlines without, varying components.

Massive stars exploding in a He-rich circumstellar medium - V. Observations of the slow-evolving SN Ibn OGLE-2012-SN-006

We present optical observations of the peculiar Type Ibn supernova (SN Ibn) OGLE-2012-SN-006, discovered and monitored by the Optical Gravitational Lensing Experiment-IV survey, and spectroscopically followed by Public ESO Spectroscopic Survey of Transient Objects (PESSTO) at late phases. Stringent pre-discovery limits constrain the explosion epoch with fair precision to JD = 245 6203.8 +/- 4.0. The rise time to the I-band light-curve maximum is about two weeks. The object reaches the peak absolute magnitude M-I = -19.65 +/- 0.19 on JD = 245 6218.1 +/- 1.8. After maximum, the light curve declines for about 25 d with a rate of 4 mag (100 d)(-1). The symmetric I-band peak resembles that of canonical Type Ib/c supernovae (SNe), whereas SNe Ibn usually exhibit asymmetric and narrower early-time light curves. Since 25 d past maximum, the light curve flattens with a decline rate slower than that of the Co-56-Fe-56 decay, although at very late phases it steepens to approach that rate. However, other observables suggest that the match with the Co-56 decay rate is a mere coincidence, and the radioactive decay is not the main mechanism powering the light curve of OGLE-2012-SN-006. An early-time spectrum is dominated by a blue continuum, with only a marginal evidence for the presence of He I lines marking this SN type. This spectrum shows broad absorptions bluewards than 5000 angstrom, likely O II lines, which are similar to spectral features observed in superluminous SNe at early epochs. The object has been spectroscopically monitored by PESSTO from 90 to 180 d after peak, and these spectra show the typical features observed in a number of SN 2006jc-like events, including a blue spectral energy distribution and prominent and narrow (v(FWHM) approximate to 1900 km s(-1)) He I emission lines. This suggests that the ejecta are interacting with He-rich circumstellar material. The detection of broad (10(4) km s(-1)) O I and Ca II features likely produced in the SN ejecta (including the [OI] lambda lambda 6300,6364 doublet in the latest spectra) lends support to the interpretation of OGLE-2012-SN-006 as a core-collapse event.

Massive stars exploding in a He-rich circumstellar medium - IV. Transitional Type Ibn supernovae

We present ultraviolet, optical and near-infrared data of the Type Ibn supernovae (SNe) 2010al and 2011hw. SN 2010al reaches an absolute magnitude at peak of M-R = -18.86 +/- 0.21. Its early light curve shows similarities with normal SNe Ib, with a rise to maximum slower than most SNe Ibn. The spectra are dominated by a blue continuum at early stages, with narrow P-Cygni He I lines indicating the presence of a slow-moving, He-rich circumstellar medium. At later epochs, the spectra well match those of the prototypical SN Ibn 2006jc, although the broader lines suggest that a significant amount of He was still present in the stellar envelope at the time of the explosion. SN 2011hw is somewhat different. It was discovered after the first maximum, but the light curve shows a double peak. The absolute magnitude at discovery is similar to that of the second peak (M-R = -18.59 +/- 0.25), and slightly fainter than the average of SNe Ibn. Though the spectra of SN 2011hw are similar to those of SN 2006jc, coronal lines and narrow Balmer lines are clearly detected. This indicates substantial interaction of the SN ejecta with He-rich, but not H-free, circumstellar material. The spectra of SN 2011hw suggest that it is a transitional SN Ibn/IIn event similar to SN 2005la. While for SN 2010al the spectrophotometric evolution favours a H-deprived Wolf-Rayet progenitor (of WN-type), we agree with the conclusion of Smith et al. that the precursor of SN 2011hw was likely in transition from a luminous blue variable to an early Wolf-Rayet (Ofpe/WN9) stage.

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.

Chemistry and Distribution of Daughter Species in the Circumstellar Envelopes of O-Rich AGB Stars

Context. Thanks to the advent of Herschel and ALMA, new high-quality observations of molecules present in the circumstellar envelopes of asymptotic giant branch (AGB) stars are being reported that reveal large differences from the existing chemical models. New molecular data and more comprehensive models of the chemistry in circumstellar envelopes are now available.
Aims: The aims are to determine and study the important formation and destruction pathways in the envelopes of O-rich AGB stars and to provide more reliable predictions of abundances, column densities, and radial distributions for potentially detectable species with physical conditions applicable to the envelope surrounding IK Tau.
Methods: We use a large gas-phase chemical model of an AGB envelope including the effects of CO and N2 self-shielding in a spherical geometry and a newly compiled list of inner-circumstellar envelope parent species derived from detailed modeling and observations. We trace the dominant chemistry in the expanding envelope and investigate the chemistry as a probe for the physics of the AGB phase by studying variations of abundances with mass-loss rates and expansion velocities.
Results: We find a pattern of daughter molecules forming from the photodissociation products of parent species with contributions from ion-neutral abstraction and dissociative recombination. The chemistry in the outer zones differs from that in traditional PDRs in that photoionization of daughter species plays a significant role. With the proper treatment of self-shielding, the N → N2 and C+→ CO transitions are shifted outward by factors of 7 and 2, respectively, compared with earlier models. An upper limit on the abundance of CH4 as a parent species of (≲2.5 × 10-6 with respect to H2) is found for IK Tau, and several potentially observable molecules with relatively simple chemical links to other parent species are determined. The assumed stellar mass-loss rate, in particular, has an impact on the calculated abundances of cations and the peak-abundance radius of both cations and neutrals: as the mass-loss rate increases, the peak abundance of cations generally decreases and the peak-abundance radius of all species moves outwards. The effects of varying the envelope expansion velocity and cosmic-ray ionization rate are not as significant.

The VLT-FLAMES Tarantula Survey XXII. Multiplicity properties of the B-type stars

We investigate the multiplicity properties of 408 B-type stars observed in the 30 Doradus region of the Large Magellanic Cloud with multi-epoch spectroscopy from the VLT-FLAMES Tarantula Survey (VFTS). We use a cross-correlation method to estimate relative radial velocities from the helium and metal absorption lines for each of our targets. Objects with significant radial-velocity variations (and with an amplitude larger than 16 km s^-1) are classified as spectroscopic binaries. We find an observed spectroscopic binary fraction (defined by periods of <10^3.5 d and mass ratios >0.1) for the B-type stars, f<inf>B</inf>(obs) = 0.25 ± 0.02, which appears constant across the field of view, except for the two older clusters (Hodge 301 and SL 639). These two clusters have significantly lower binary fractions of 0.08 ± 0.08 and 0.10 ± 0.09, respectively. Using synthetic populations and a model of our observed epochs and their potential biases, we constrain the intrinsic multiplicity properties of the dwarf and giant (i.e. relatively unevolved) B-type stars in 30 Dor. We obtain a present-day binary fraction f<inf>B</inf>(true) = 0.58 ± 0.11, with a flat period distribution. Within the uncertainties, the multiplicity properties of the B-type stars agree with those for the O stars in 30 Dor from the VFTS.

Coronal properties of planet-bearing stars

Context. Do extrasolar planets affect the activity of their host stars? Indications for chromospheric activity enhancement have been found for a handful of targets, but in the X-ray regime, conclusive observational evidence is still missing. Aims: We want to establish a sound observational basis to confirm or reject major effects of Star-Planet Interactions (SPI) in stellar X-ray emissions. Methods: We therefore conduct a statistical analysis of stellar X-ray activity of all known planet-bearing stars within 30 pc distance for dependencies on planetary parameters such as mass and semimajor axis. Results: In our sample, there are no significant correlations of X-ray luminosity or the activity indicator L_X/L_bol with planetary parameters which cannot be explained by selection effects. Conclusions: Coronal SPI seems to be a phenomenon which might only manifest itself as a strong effect for a few individual targets, but not to have a major effect on planet-bearing stars in general.

Indications for an influence of hot Jupiters on the rotation and activity of their host stars

Context. The magnetic activity of planet-hosting stars is an importantfactor for estimating the atmospheric stability of close-in exoplanetsand the age of their host stars. It has long been speculated thatclose-in exoplanets can influence the stellar activity level. However,testing for tidal or magnetic interaction effects in samples ofplanet-hosting stars is difficult because stellar activity hindersexoplanet detection, so that stellar samples with detected exoplanetsshow a bias toward low activity for small exoplanets.

Aims: Weaim to test whether exoplanets in close orbits influence the stellarrotation and magnetic activity of their host stars.

Methods: Wedeveloped a novel approach to test for systematic activity-enhancementsin planet-hosting stars. We use wide (several 100 AU) binary systems inwhich one of the stellar components is known to have an exoplanet, whilethe second stellar component does not have a detected planet andtherefore acts as a negative control. We use the stellar coronal X-rayemission as an observational proxy for magnetic activity and analyzeobservations performed with Chandra and XMM-Newton.

Results: Wefind that in two systems for which strong tidal interaction can beexpected the planet-hosting primary displays a much higher magneticactivity level than the planet-free secondary. In three systems forwhich weaker tidal interaction can be expected the activity levels ofthe two stellar components agree with each other.

Conclusions:Our observations indicate that the presence of Hot Jupiters may inhibitthe spin-down of host stars with thick outer convective layers. Possiblecauses for this effect include a transfer of angular momentum from theplanetary orbit to the stellar rotation through tidal interaction, ordifferences during the early evolution of the system, where the hoststar may decouple from the protoplanetary disk early because of a gapopened by the forming Hot Jupiter.

Magnetospheric Structure and Atmospheric Joule Heating of Habitable Planets Orbiting M-dwarf Stars

We study the magnetospheric structure and the ionospheric Joule Heating of planets orbiting M-dwarf stars in the habitable zone using a set of magnetohydrodynamic models. The stellar wind solution is used to drive a model for the planetary magnetosphere, which is coupled with a model for the planetary ionosphere. Our simulations reveal that the space environment around close-in habitable planets is extreme, and the stellar wind plasma conditions change from sub- to super-Alfvénic along the planetary orbit. As a result, the magnetospheric structure changes dramatically with a bow shock forming in the super-Alfvénic sectors, while no bow shock forms in the sub-Alfvénic sectors. The planets reside most of the time in the sub-Alfvénic sectors with poor atmospheric protection. A significant amount of Joule Heating is provided at the top of the atmosphere as a result of the intense stellar wind. For the steady-state solution, the heating is about 0.1%-3% of the total incoming stellar irradiation, and it is enhanced by 50% for the time-dependent case. The significant Joule Heating obtained here should be considered in models for the atmospheres of habitable planets in terms of the thickness of the atmosphere, the top-side temperature and density, the boundary conditions for the atmospheric pressure, and particle radiation and transport. Here we assume constant ionospheric Pedersen conductance similar to that of the Earth. The conductance could be greater due to the intense EUV radiation leading to smaller heating rates. We plan to quantify the ionospheric conductance in future study.

Explaining the Coexistence of Large-scale and Small-scale Magnetic Fields in Fully Convective Stars

Despite the lack of a shear-rich tachocline region, low-mass fully convective (FC) stars are capable of generating strong magnetic fields, indicating that a dynamo mechanism fundamentally different from the solar dynamo is at work in these objects. We present a self-consistent three-dimensional model of magnetic field generation in low-mass FC stars. The model utilizes the anelastic magnetohydrodynamic equations to simulate compressible convection in a rotating sphere. A distributed dynamo working in the model spontaneously produces a dipole-dominated surface magnetic field of the observed strength. The interaction of this field with the turbulent convection in outer layers shreds it, producing small-scale fields that carry most of the magnetic flux. The Zeeman–Doppler-Imaging technique applied to synthetic spectropolarimetric data based on our model recovers most of the large-scale field. Our model simultaneously reproduces the morphology and magnitude of the large-scale field as well as the magnitude of the small-scale field observed on low-mass FC stars.

An age-activity calibration for old low-mass stars

Low-mass stars are highly interesting targets: we are able to detect planets in their habitable zones, and upcoming searches for biomarkers in exoplanet atmospheres will focus on low-mass star systems due to their ubiquity and proximity. We aim to develop an age-activity calibration for old low-mass stars, using wide binary systems consisting of an M or K dwarf and a white dwarf. The age of the system is determined by the WD cooling time plus its progenitor lifetime, yielding reliable ages in the regime >1 Gyr. For an exploratory sample of 7 systems where we have already derived ages, we propose to perform Chandra ACIS-S observations to determine the X-ray luminosities of the M dwarfs and correlate their stellar activity with age. We ask for a total observing time of 110 ks.

Do exoplanets spin up their host stars?

We have collected initial evidence that tidal interaction between a late-type star and its close-in, massive planet can lead to a spin-up of the host star. We propose to explore this further by studying a small sample of proper motion pairs in which one of the stars is orbited by a Hot Jupiter. We will determine if the activity-estimated age appears to be strongly different for the two stars, which would indicate a tidal spin up of the Hot Jupiter host star. We propose to observe 4 such systems with Chandra/ACIS-S, and to perform a similar observation of one additional system with large angular separation using XMM-Newton/EPIC. The total proposed exposure times are 141 ks (Chandra) and 38 ks (XMM).