Detecting the progenitors of core collapse supernovae

The masses and the evolutionary states of the progenitors of core-collapse supernovae are not well constrained by direct observations. Stellar evolution theory generally predicts that massive stars with initial masses less than about 30M_sol should undergo core-collapse when they are cool M-type supergiants. However the only two detections of a SN progenitor before explosion are SN1987A and SN1993J, and neither of these was an M-type supergiant. Attempting to identify the progenitors of supernovae is a difficult task, as precisely predicting the time of explosion of a massive star is impossible for obvious reasons. There are several different types of supernovae which have different spectral and photometric evolution, and how exactly these are related to the evolutionary states of the progenitor stars is not currently known. I will describe a novel project which may allow the direct identification of core-collapse supernovae progenitors on pre-explosion images of resolved, nearby galaxies. This project is now possible with the excellent image archives maintained by several facilities and will be enhanced by the new initiatives to create Virtual Observatories, the earliest of which ASTROVIRTEL is already producing results.

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.

Machine learning for transient discovery in Pan-STARRS1 difference imaging

Efficient identification and follow-up of astronomical transients is hindered by the need for humans to manually select promising candidates from data streams that contain many false positives. These artefacts arise in the difference images that are produced by most major ground-based time-domain surveys with large format CCD cameras. This dependence on humans to reject bogus detections is unsustainable for next generation all-sky surveys and significant effort is now being invested to solve the problem computationally. In this paper, we explore a simple machine learning approach to real-bogus classification by constructing a training set from the image data of similar to 32 000 real astrophysical transients and bogus detections from the Pan-STARRS1 Medium Deep Survey. We derive our feature representation from the pixel intensity values of a 20 x 20 pixel stamp around the centre of the candidates. This differs from previous work in that it works directly on the pixels rather than catalogued domain knowledge for feature design or selection. Three machine learning algorithms are trained (artificial neural networks, support vector machines and random forests) and their performances are tested on a held-out subset of 25 per cent of the training data. We find the best results from the random forest classifier and demonstrate that by accepting a false positive rate of 1 per cent, the classifier initially suggests a missed detection rate of around 10 per cent. However, we also find that a combination of bright star variability, nuclear transients and uncertainty in human labelling means that our best estimate of the missed detection rate is approximately 6 per cent.

Energy levels, radiative rates and electron impact excitation rates for transitions in AlX

Energy levels, radiative rates and lifetimes are calculated among the lowest 98 levels of the n ≤4 configurations of Be-like Al X. The GRASP (General-purpose Relativistic Atomic Structure Package) is adopted and data are provided for all E1, E2, M1 and M2 transitions. Similar data are also obtained with the FAC (Flexible Atomic Code) to assess the accuracy of the calculations. Based on comparisons between calculations with the two codes as well as with available measurements, our listed energy levels are assessed to be accurate to better than 0.3 per cent. However, the accuracy for radiative rates and lifetimes is estimated to be about 20 per cent. Collision strengths are also calculated for which the DARC (Dirac Atomic R-matrix Code) is used. A wide energy range (up to 380 Ryd) is considered and resonances resolved in a fine energy mesh in the thresholds region. The collision strengths are subsequently averaged over a Maxwellian velocity distribution to determine effective collision strengths up to a temperature of 1.6 × 10⁷ K. Our results are compared with the previous (limited) atomic data and significant differences (up to a factor of 4) are noted for several transitions, particularly those which are not allowed in jj coupling. 

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. 

Energy levels, radiative rates and electron impact excitation rates for transitions in Si II

Energies for the lowest 56 levels, belonging to the 3s² 3p, 3s 3p², 3p³, 3s² 3d, 3s 3p 3d, 3s² 4ℓ and 3s² 5ℓ configurations of Si II, are calculated using the General-purpose Relativistic Atomic Structure Package (GRASP) code. Analogous calculations have also been performed (for up to 175 levels) using the FlexibleAtomicCode (FAC). Furthermore, radiative rates are calculated for all E1, E2, M1 and M2 transitions. Extensive comparisons are made with available theoretical and experimental energy levels, and the accuracy of the present results is assessed to be better than 0.1Ryd. Similarly, the accuracy for radiative rates (and subsequently lifetimes) is estimated to be better than 20 per cent for most of the (strong) transitions. Electron impact excitation collision strengths are also calculated, with the Dirac Atomic R-matrix Code (DARC), over a wide energy range up to 13 Ryd. Finally, to determine effective collision strengths, resonances are resolved in a fine energy mesh in the thresholds region. These collision strengths are averaged over a Maxwellian velocity distribution and results listed over a wide range of temperatures, up to 10^5.5 K. Our data are compared with earlier R-matrix calculations and differences noted, up to a factor of 2, for several transitions. Although scope remains for improvement, the accuracy for our results of collision strengths and effective collision strengths is assessed to be about 20 per cent for a majority of transitions. 

Transit timing variations in WASP-10b induced by stellar activity

The hot-JupiterWASP-10bwas reported by Maciejewski et al. to showtransit timing variations (TTVs) with an amplitude of ~3.5 min. These authors proposed that the observed TTVs were caused by a 0.1M_Jup perturbing companion with an orbital period of ~5.23 d, and hence, close to the outer 5:3 mean-motion resonance with WASP-10b. To test this scenario, we present eight new transit light curves of WASP-10b obtained with the Faulkes Telescope North and the Liverpool Telescope. The new light curves, together with 22 previously published ones, were modelled with a Markov Chain Monte Carlo transit fitting code. Transit depth differences reported forWASP-10b are thought to be due to starspot-induced brightness modulation of the host star. Assuming the star is brighter at the activity minimum, we favour a small planetary radius. We find Rp = 1.039^+0.043 _-0.049R_Jup in agreement with Johnson et al. and Maciejewski et al. Recent studies find no evidence for a significant eccentricity in this system. We present consistent system parameters for a circular orbit and refine the orbital ephemeris ofWASP-10b. Our homogeneously derived transit times do not support the previous claimed TTV signal, which was strongly dependent on two previously published transits that have been incorrectly normalized. Nevertheless, a linear ephemeris is not a statistically good fit to the transit times of WASP-10b. We show that the observed transit time variations are due to spot occultation features or systematics. We discuss and exemplify the effects of occultation spot features in the measured transit times and show that despite spot occultation during egress and ingress being difficult to distinguish in the transit light curves, they have a significant effect in the measured transit times. We conclude that if we account for spot features, the transit times of WASP-10b are consistent with a linear ephemeris with the exception of one transit (epoch 143) which is a partial transit. Therefore, there is currently no evidence for the existence of a companion to WASP-10b. Our results support the lack of TTVs of hot-Jupiters reported for the Kepler sample. 

Expanded iron uta spectra - Probing the thermal stability limits in AGN clouds

The Fe unresolved transition arrays (UTAs) produce prominent features in the 15-17 Å wavelength range in the spectra of active galactic nuclei (AGNs). Here, we present new calculations of the energies and oscillator strengths of inner-shell lines from Fe XIV, Fe XV, and Fe XVI. These are crucial ions since they are dominant at inflection points in the gas thermal stability curve, and UTA excitation followed by autoionization is an important ionization mechanism for these species. We incorporate these, and data reported in previous papers, into the plasma simulation code Cloudy. This updated physics is subsequently employed to reconsider the thermally stable phases in absorbing media in AGNs. We show how the absorption profile of the Fe XIV UTA depends on density, due to the changing populations of levels within the ground configuration. © 2013. The American Astronomical Society. All rights reserved.

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.

Energy levels, radiative rates and electron impact excitation rates for transitions in C III

We report energy levels, radiative rates (A-values) and lifetimes for the astrophysically important Be-like ion C III. For the calculations, 166 levels belonging to the n ≤ 5 configurations are considered and the GRASP (General-purpose Relativistic Atomic Structure Package) is adopted. Einstein A-coefficients are provided for all E1, E2, M1 and M2 transitions, while lifetimes are compared with available measurements as well as theoretical results, and no large discrepancies noted. Our energy levels are assessed to be accurate to better than 1 per cent for a majority of levels, and A-values to better than 20 per cent for most transitions. Collision strengths are also calculated, for which the Dirac Atomic R-matrix Code (DARC) is used. A wide energy range, up to 21 Ryd, is considered and resonances resolved in a fine energy mesh in the thresholds region. The collision strengths are subsequently averaged over a Maxwellian velocity distribution to determine effective collision strengths up to a temperature of 8.0 ×10[5]K, sufficient for most astrophysical applications. Our data are compared with the recent R-matrix calculations of Fernández-Menchero et al., and significant differences (up to over an order of magnitude) are noted for several transitions over the complete temperature range of the results.

Electron impact excitation of Be-like ions: a comparison of DARC and ICFT results

Emission lines of Be-like ions are frequently observed in astrophysical plasmas, and many are useful for density and temperature diagnostics. However, accurate atomic data for energy levels, radiative rates (A-values) and effective electron excitation collision strengths ($\Upsilon$) are required for reliable plasma modelling. In general it is reasonably straightforward to calculate energy levels and A- values to a high level of accuracy. By contrast, considerable effort is required to calculate $\Upsilon$, and hence it is not always possible to assess the accuracy of available data. Recently, two independent calculations (adopting the $R$-matrix method) but with different approaches (DARC and ICFT) have appeared for a range of Be-like ions. Therefore, in this work we compare the two sets of $\Upsilon$, highlight the large discrepancies for a significant number of transitions and suggest possible reasons for these.

Energy levels, radiative rates and electron impact excitation rates for transitions in Fe XIV

Energies and lifetimes are reported for the lowest 136 levels of Fe XIV, belonging to the (1s(2)2s(2)2p(6)) 3s(2)3p, 3s(3)p(2), 3s(2)3d, 3p(3), 3s(3)p(3)d, 3p(2)3d, 3s(3)d(2), 3p(3)d(2) and 3s(2)4l configurations. Additionally, radiative rates, oscillator strengths and line strengths are calculated for all electric dipole (E1), magnetic dipole (M1), electric quadrupole (E2) and magnetic quadrupole (M2) transitions. Theoretical lifetimes determined from these radiative rates for most levels show satisfactory agreement with earlier calculations, a swell as with measurements. Electron impact excitation collision strengths are also calculated with the Dirac atomic R-matrix code (DARC) over a wide energy range up to 260 Ryd. Furthermore, resonances have been resolved in a fine energy mesh to determine effective collision strengths, obtained after integrating the collision strengths over a Maxwellian distribution of electron velocities. Results are listed for all 9180 transitions among the 136 levels over a wide range of electron temperatures, up to 10(7.1) K. Comparisons are made with available results in the literature, and the accuracy of the present data is assessed.

Large- and small-scale structure of the intermediate- and high-velocity clouds towards the LMC and SMC.

We employ Ca II K and Na I D interstellar absorption-line spectroscopy of early-type stars in the Large and Small Magellanic Clouds (LMC, SMC) to investigate the large- and small-scale structure in foreground intermediate- and high-velocity clouds (I/HVCs). Data include FLAMES-GIRAFFE Ca II K observations of 403 stars in four open clusters, plus FEROS or UVES spectra of 156 stars in the LMC and SMC. The FLAMES observations are amongst the most extensive probes to date of Ca II structures on ∼20 arcsec scales in Magellanic I/HVCs. From the FLAMES data within a 0 ∘.∘.∘.5 field of view, the Ca II K equivalent width in the I/HVC components towards three clusters varies by factors of ≥10. There are no detections of molecular gas in absorption at intermediate or high velocities, although molecular absorption is present at LMC and Galactic velocities towards some sightlines. The FEROS/UVES data show Ca II K I/HVC absorption in ∼60 per cent of sightlines. The range in the Ca II/Na I ratio in I/HVCs is from –0.45 to +1.5 dex, similar to previous measurements for I/HVCs. In 10 sightlines we find Ca II/O I ratios in I/HVC gas ranging from 0.2 to 1.5 dex below the solar value, indicating either dust or ionization effects. In nine sightlines I/HVC gas is detected in both H I and Ca II at similar velocities, implying that the two elements form part of the same structure.

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).