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.

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.

Hα emission fluxes and lithium abundances of low-mass stars in the young open cluster IC 4665

As part of a long term effort to understand pre-main sequence Li burning, we have obtained high resolution spectroscopic observations of 14 late type stars (G0-M1) in the young open cluster IC 4665. Most of the stars have Hα filled-in and Li I absorption, as expected for their young age. From the equivalent widths of Hα emission excess (obtained using the spectral subtraction technique) and the the Li i λ6708 feature, we have derived Hα emission fluxes and photospheric Li abundances. The mean Li abundance of IC 4665 solar-type stars is log N(Li) = 3.1; the same as in other young clusters (α Per, Pleiades) and T Tauri stars. Our results support the conclusions from previous works that PMS Li depletion is very small for masses ∼ 1 M_⨀ . Among the IC 4665 late-G and early K-type stars, there is a spread in Li abundances of about one order of magnitude. The Li-poor IC 4665 members have low Hα excess and vsini≤10. Hence, the Li-activity-rotation connection which has been clearly established in the Pleiades also seems to hold in IC 4665. One M-type IC 4665 star that we have observed does not show Li, implying a very efficient Li depletion as observed in α Per stars of the same spectral type. The level of chromospheric activity and Li depletion among the low-mass stars of IC 4665 is similar to that in the Pleiades. In fact, we note that the Li abundance distributions in several young clusters (α Per, Pleiades, IC 2391, IC 4665) and in post T Tauri stars are strikingly similar. This result suggests that Hα emission and Li abundance not well correlated with age for low-mass stars between 20 and 100 Myr old. We argue that a finer age indicator, the ''LL-clock'', would be the luminosity at which the transition between efficient Li depletion and preservation takes place for fully convective objects. The LL-clock could allow in the near future to derive the relative ages of young open clusters, and clarify the study of PMS evolution of cool stars.

The VLT-FLAMES survey of massive stars: stellar parameters and rotational velocities in NGC 3293, NGC 4755 and NGC 6611

An analysis is presented of VLT-FLAMES spectroscopy for three Galactic clusters, NGC3293, NGC4755 and NGC6611. Non-LTE model atmosphere calculations have been used to estimate effective temperatures (from either the helium spectrum or the silicon ionization equilibrium) and gravities (from the hydrogen spectrum). Projected rotational velocities have been deduced from the helium spectrum (for fast and moderate rotators) or the metal line spectrum (for slow rotators). The origin of the low gravity estimates for apparently near main sequence objects is discussed and is related to the stellar rotational velocity. The atmospheric parameters have been used to estimate cluster distances (which are generally in good agreement with previous determinations) and these have been used to estimate stellar luminosities and evolutionary masses. The observed Hertzsprung-Russell diagrams are compared with theoretical predictions and some discrepancies including differences in the main sequence luminosities are discussed. Cluster ages have been deduced and evidence for non-coeval star formation is found for all three of the clusters. Projected rotational velocities for targets in the older clusters, NGC3293 and NGC4755, have been found to be systematically larger than those for the field, confirming recent results in other similar age clusters. The distribution of projected rotational velocities are consistent with a Gaussian distribution of intrinsic rotational velocities. For the relatively unevolved targets in the older clusters, NGC3293 and NGC4755, the peak of the velocity distribution would be 250 km s(-1) with a full-width-half-maximum of approximately 180 km s(-1). For NGC6611, the sample size is relatively small but implies a lower mean rotational velocity. This may be evidence for the spin-down effect due to angular momentum loss through stellar winds, although our results are consistent with those found for very young high mass stars. For all three clusters we deduce present day mass functions with Gamma-values in the range of -1.5 to -1.8, which are similar to other young stellar clusters in the Milky Way.

Low luminosity TypeII supernovae-II. Pointing towards moderate mass precursors

We present new data for five underluminous Type II-plateau supernovae (SNe IIP), namely SN 1999gn, SN 2002gd, SN 2003Z, SN 2004eg and SN 2006ov. This new sample of lowluminosity SNe IIP (LL SNe IIP) is analysed together with similar objects studied in the past. All of them show a flat light-curve plateau lasting about 100 d, an underluminous late-time exponential tail, intrinsic colours that are unusually red, and spectra showing prominent and narrow P Cygni lines. A velocity of the ejected material below 10³ km s^-1 is inferred from measurements at the end of the plateau. The 56Ni masses ejected in the explosion are very small (≤10^-2 M⊙). We investigate the correlations among ⁵⁶Ni mass, expansion velocity of the ejecta and absolute magnitude in the middle of the plateau, confirming the main findings of Hamuy, according to which events showing brighter plateau and larger expansion velocities are expected to produce more 56Ni. We propose that these faint objects represent the LL tail of a continuous distribution in parameters space of SNe IIP. The physical properties of the progenitors at the explosion are estimated through the hydrodynamical modelling of the observables for two representative events of this class, namely SN 2005cs and SN 2008in. We find that the majority of LL SNe IIP, and quite possibly all, originate in the core collapse of intermediate-mass stars, in the mass range 10-15 M⊙. 

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.

Imprints of fast-rotating massive stars in the Galactic Bulge

The first stars that formed after the Big Bang were probably massive(1), and they provided the Universe with the first elements heavier than helium (`metals`), which were incorporated into low-mass stars that have survived to the present(2,3). Eight stars in the oldest globular cluster in the Galaxy, NGC 6522, were found to have surface abundances consistent with the gas from which they formed being enriched by massive stars(4) (that is, with higher alpha-element/Fe and Eu/Fe ratios than those of the Sun). However, the same stars have anomalously high abundances of Ba and La with respect to Fe(4), which usually arises through nucleosynthesis in low-mass stars(5) (via the slow-neutron-capture process, or s-process). Recent theory suggests that metal-poor fast-rotating massive stars are able to boost the s-process yields by up to four orders of magnitude(6), which might provide a solution to this contradiction. Here we report a reanalysis of the earlier spectra, which reveals that Y and Sr are also over-abundant with respect to Fe, showing a large scatter similar to that observed in extremely metal-poor stars(7), whereas C abundances are not enhanced. This pattern is best explained as originating in metal-poor fast-rotating massive stars, which might point to a common property of the first stellar generations and even of the `first stars`.

The Gemini NICI planet-finding campaign: the frequency of giant planets around young B and A stars

We have carried out high contrast imaging of 70 young, nearby B and A stars to search for brown dwarf and planetary companions as part of the Gemini NICI Planet-Finding Campaign. Our survey represents the largest, deepest survey for planets around high-mass stars (≈1.5-2.5 M ☉) conducted to date and includes the planet hosts β Pic and Fomalhaut. We obtained follow-up astrometry of all candidate companions within 400 AU projected separation for stars in uncrowded fields and identified new low-mass companions to HD 1160 and HIP 79797. We have found that the previously known young brown dwarf companion to HIP 79797 is itself a tight (3 AU) binary, composed of brown dwarfs with masses 58$^{+21}_{-20}$ M Jup and 55$^{+20}_{-19}$ M Jup, making this system one of the rare substellar binaries in orbit around a star. Considering the contrast limits of our NICI data and the fact that we did not detect any planets, we use high-fidelity Monte Carlo simulations to show that fewer than 20% of 2 M ☉ stars can have giant planets greater than 4 M Jup between 59 and 460 AU at 95% confidence, and fewer than 10% of these stars can have a planet more massive than 10 M Jup between 38 and 650 AU. Overall, we find that large-separation giant planets are not common around B and A stars: fewer than 10% of B and A stars can have an analog to the HR 8799 b (7 M Jup, 68 AU) planet at 95% confidence. We also describe a new Bayesian technique for determining the ages of field B and A stars from photometry and theoretical isochrones. Our method produces more plausible ages for high-mass stars than previous age-dating techniques, which tend to underestimate stellar ages and their uncertainties.

Stars and brown dwarfs in the σ Orionis cluster (Research Note) III. OSIRIS/GTC low-resolution spectroscopy of variable sources

Context. Although many studies have been performed so far, there are still dozens of low-mass stars and brown dwarfs in the young σ Orionis open cluster without detailed spectroscopic characterisation. Aims. We look for unknown strong accretors and disc hosts that were undetected in previous surveys. Methods. We collected low-resolution spectroscopy (R ~ 700) of ten low-mass stars and brown dwarfs in σ Orionis with OSIRIS at the Gran Telescopio Canarias under very poor weather conditions. These objects display variability in the optical, infrared, Hα, and/or X-rays on time scales of hours to years. We complemented our spectra with optical and near-/mid-infrared photometry. Results. For seven targets, we detected lithium in absorption, identified Hα, the calcium doublet, and forbidden lines in emission, and/or determined spectral types for the first time. We characterise in detail a faint, T Tauri-like brown dwarf with an 18 h-period variability in the optical and a large Hα equivalent width of –125  ±  15 Å, as well as two M1-type, X-ray-flaring, low-mass stars, one with a warm disc and forbidden emission lines, the other with a previously unknown cold disc with a large inner hole. Conclusions. New unrevealed strong accretors and disc hosts, even below the substellar limit, await discovery among the list of known σ Orionis stars and brown dwarfs that are variable in the optical and have no detailed spectroscopic characterisation yet.

Molecular line and continuum studies of the early stages of star formation

New stars form in dense interstellar clouds of gas and dust called molecular clouds. The actual sites where the process of star formation takes place are the dense clumps and cores deeply embedded in molecular clouds. The details of the star formation process are complex and not completely understood. Thus, determining the physical and chemical properties of molecular cloud cores is necessary for a better understanding of how stars are formed. Some of the main features of the origin of low-mass stars, like the Sun, are already relatively well-known, though many details of the process are still under debate. The mechanism through which high-mass stars form, on the other hand, is poorly understood. Although it is likely that the formation of high-mass stars shares many properties similar to those of low-mass stars, the very first steps of the evolutionary sequence are unclear. Observational studies of star formation are carried out particularly at infrared, submillimetre, millimetre, and radio wavelengths. Much of our knowledge about the early stages of star formation in our Milky Way galaxy is obtained through molecular spectral line and dust continuum observations. The continuum emission of cold dust is one of the best tracers of the column density of molecular hydrogen, the main constituent of molecular clouds. Consequently, dust continuum observations provide a powerful tool to map large portions across molecular clouds, and to identify the dense star-forming sites within them. Molecular line observations, on the other hand, provide information on the gas kinematics and temperature. Together, these two observational tools provide an efficient way to study the dense interstellar gas and the associated dust that form new stars. The properties of highly obscured young stars can be further examined through radio continuum observations at centimetre wavelengths. For example, radio continuum emission carries useful information on conditions in the protostar+disk interaction region where protostellar jets are launched. In this PhD thesis, we study the physical and chemical properties of dense clumps and cores in both low- and high-mass star-forming regions. The sources are mainly studied in a statistical sense, but also in more detail. In this way, we are able to examine the general characteristics of the early stages of star formation, cloud properties on large scales (such as fragmentation), and some of the initial conditions of the collapse process that leads to the formation of a star. The studies presented in this thesis are mainly based on molecular line and dust continuum observations. These are combined with archival observations at infrared wavelengths in order to study the protostellar content of the cloud cores. In addition, centimetre radio continuum emission from young stellar objects (YSOs; i.e., protostars and pre-main sequence stars) is studied in this thesis to determine their evolutionary stages. The main results of this thesis are as follows: i) filamentary and sheet-like molecular cloud structures, such as infrared dark clouds (IRDCs), are likely to be caused by supersonic turbulence but their fragmentation at the scale of cores could be due to gravo-thermal instability; ii) the core evolution in the Orion B9 star-forming region appears to be dynamic and the role played by slow ambipolar diffusion in the formation and collapse of the cores may not be significant; iii) the study of the R CrA star-forming region suggests that the centimetre radio emission properties of a YSO are likely to change with its evolutionary stage; iv) the IRDC G304.74+01.32 contains candidate high-mass starless cores which may represent the very first steps of high-mass star and star cluster formation; v) SiO outflow signatures are seen in several high-mass star-forming regions which suggest that high-mass stars form in a similar way as their low-mass counterparts, i.e., via disk accretion. The results presented in this thesis provide constraints on the initial conditions and early stages of both low- and high-mass star formation. In particular, this thesis presents several observational results on the early stages of clustered star formation, which is the dominant mode of star formation in our Galaxy.

Star formation in giant extragalactic H II regions

Star formation properties in Giant Extragalactic H II Regions (GEHRs) are investigated using optical photometry and evolutionary population synthesis models. Photometric data in $BVR$ bands and in the emission line of H-alpha are obtained by CCD imaging at Vainu Bappu Observatory, Kavalur. Aperture photometry is performed for 180 GEHRs in galaxies NGC 1365, 1566, 2366, 2903, 2997, 3351, 4303, 4449, 4656 and 5253. Thirty six of these GEHRs having published spectroscopic data are studied for star formation properties. The population synthesis model is constructed based on Maeder's stellar evolutionary and Kurucz stellar atmosphere models, to synthesize observational quantities of embedded clusters in GEHRs. The observed H-alpha luminosity is a measure of the number of massive stars while the contribution to BVR bands is from intermediate mass (5-15 solar mass) stars when the cluster is young and from evolving supergiants when the cluster is old (age >/= 6~Myr). Differential reddening between gas and embedded stars is essential to constrain the dereddened cluster colors within the range of youngest clusters. Obscuring dust closely associated with gas, which is distributed in filaments and clumps, as in the case of 30 Doradus, is the most likely configuration giving rise to net reduction of extinction towards stars. The fraction of the stellar photons escaping the nebula unattenuated is estimated to be 50%. GEHRs are rarely found to be simple systems containing stars from single generation. In the present sample such regions in addition to being older than 3~Myr, have their Lyman continuum luminosity reduced by as much as 60%, compared to the observed $B$ band luminosity for a normal IMF. The missing ionizing photons may be escaping the nebula, leading to the ionization of extra-H II region ionized medium. Co-existence of young (age stars producing ionizing photons) and old populations (~10~Myr; Red Supergiants) is found to be common in GEHRs. The emission and continuum knots are seen spatially separated (40-100 pc) on CCD images in NGC 2997, 4303 and 4449 and may be direct evidences for the co-existence of young and old populations in giant star forming complexes. Triggering of star formation from earlier bursts is the most likely cause of new generation of stars, and may be a common phenomenon in GEHRs. Spatial separation between the young and old stars (~30 pc) had been earlier reported in 30 Doradus. Thus GEHRs in nearby galaxies share many of the properties shown by 30 Dor, the nearest GEHR. (SECTION: Dissertation Summaries)

Grain growth in protoplanetary disks

The majority of young, low-mass stars are surrounded by optically thick accretion disks. These circumstellar disks provide large reservoirs of gas and dust that will eventually be transformed into planetary systems. Theory and observations suggest that the earliest stage toward planet formation in a protoplanetary disk is the growth of particles, from sub-micron-sized grains to centimeter- sized pebbles. Theory indicates that small interstellar grains are well coupled into the gas and are incorporated to the disk during the proto-stellar collapse. These dust particles settle toward the disk mid-plane and simultaneously grow through collisional coagulation in a very short timescale. Observationally, grain growth can be inferred by measuring the spectral energy distribution at long wavelengths, which traces the continuum dust emission spectrum and hence the dust opacity. Several observational studies have indicated that the dust component in protoplanetary disks has evolved as compared to interstellar medium dust particles, suggesting at least 4 orders of magnitude in particle- size growth. However, the limited angular resolution and poor sensitivity of previous observations has not allowed for further exploration of this astrophysical process.

As part of my thesis, I embarked in an observational program to search for evidence of radial variations in the dust properties across a protoplanetary disk, which may be indicative of grain growth. By making use of high angular resolution observations obtained with CARMA, VLA, and SMA, I searched for radial variations in the dust opacity inside protoplanetary disks. These observations span more than an order of magnitude in wavelength (from sub-millimeter to centimeter wavelengths) and attain spatial resolutions down to 20 AU. I characterized the radial distribution of the circumstellar material and constrained radial variations of the dust opacity spectral index, which may originate from particle growth in these circumstellar disks. Furthermore, I compared these observational constraints with simple physical models of grain evolution that include collisional coagulation, fragmentation, and the interaction of these grains with the gaseous disk (the radial drift problem). For the parameters explored, these observational constraints are in agreement with a population of grains limited in size by radial drift. Finally, I also discuss future endeavors with forthcoming ALMA observations.

Stellar jitter from variable gravitational redshift: implications for radial velocity confirmation of habitable exoplanets

A variation of gravitational redshift, arising from stellar radius fluctuations, will introduce astrophysical noise into radial velocity measurements by shifting the centroid of the observed spectral lines. Shifting the centroid does not necessarily introduce line asymmetries. This is fundamentally different from other types of stellar jitter so far identified, which do result from line asymmetries. Furthermore, only a very small change in stellar radius, ˜0.01 per cent, is necessary to generate a gravitational redshift variation large enough to mask or mimic an Earth-twin. We explore possible mechanisms for stellar radius fluctuations in low-mass stars. Convective inhibition due to varying magnetic field strengths and the Wilson depression of starspots are both found to induce substantial gravitational redshift variations. Finally, we investigate a possible method for monitoring/correcting this newly identified potential source of jitter and comment on its impact for future exoplanet searches.

The AstraLux large M-dwarf multiplicity survey

We present the results of an extensive high-resolution imaging survey of M-dwarf multiplicity using the Lucky Imaging technique. The survey made use of the AstraLux Norte camera at the Calar Alto 2.2m telescope and the AstraLux Sur camera at the ESO New Technology Telescope in order to cover nearly the full sky. In total, 761 stars were observed (701M-type and 60 late K-type), among which 182 new and 37 previously known companions were detected in 205 systems. Most of the targets have been observed during two or more epochs, and could be confirmed as physical companions through common proper motion, often with orbital motion being confirmed in addition. After accounting for various bias effects, we find a total M-dwarf multiplicity fraction of 27% ± 3% within the AstraLux detection range of 008-6? (semimajor axes of ~3-227 AU at a median distance of 30pc). We examine various statistical multiplicity properties within the sample, such as the trend of multiplicity fraction with stellar mass and the semimajor axis distribution. The results indicate that M-dwarfs are largely consistent with constituting an intermediate step in a continuous distribution from higher-mass stars down to brown dwarfs. Along with other observational results in the literature, this provides further indications that stars and brown dwarfs may share a common formation mechanism, rather than being distinct populations. © 2012. The American Astronomical Society. All rights reserved.

Possible binary progenitors for the Type Ib supernova iPTF13bvn

Cao et al. reported a possible progenitor detection for the Type Ib supernovae iPTF13bvn for the first time. We find that the progenitor is in fact brighter than the magnitudes previously reported by approximately 0.7-0.2 mag with a larger error in the bluer filters. We compare our new magnitudes to our large set of binary evolution models and find that many binary models with initial masses in the range of 10-20M(circle dot) match this new photometry and other constraints suggested from analysing the supernova. In addition, these lower mass stars retain more helium at the end of the model evolution indicating that they are likely to be observed as Type Ib supernovae rather than their more massive, Wolf-Rayet counter parts. We are able to rule out typical Wolf-Rayet models as the progenitor because their ejecta masses are too high and they do not fit the observed SED unless they have a massive companion which is the observed source at the supernova location. Therefore only late-time observations of the location will truly confirm if the progenitor was a helium giant and not a Wolf-Rayet star.