Massive stars exploding in a He-rich circumstellar medium - II. The transitional case of SN 2005la

We present photometric and spectroscopic data of the peculiar SN 2005la, an object which shows an optical light curve with some luminosity fluctuations and spectra with comparably strong narrow hydrogen and helium lines, probably of circumstellar nature. The increasing full width at half-maximum velocity of these lines is indicative of an acceleration of the circumstellar material. SN 2005la exhibits hybrid properties, sharing some similarities with both Type IIn supernovae and 2006jc-like (Type Ibn) events. We propose that the progenitor of SN 2005la was a very young Wolf-Rayet (WN-type) star which experienced mass ejection episodes shortly before core collapse.

The VLT-FLAMES survey of massive stars: wind properties and evolution of hot massive stars in the Large Magellanic Cloud

We have studied the optical spectra of a sample of 28 O- and early B-type stars in the Large Magellanic Cloud, 22 of which are associated with the young star forming region N11. Our observations sample the central associations of LH9 and LH10, and the surrounding regions. Stellar parameters are determined using an automated fitting method ( Mokiem et al. 2005), which combines the stellar atmosphere code fastwind ( Puls et al. 2005) with the genetic algorithm based optimisation routine PIKAIA ( Charbonneau 1995). We derive an age of 7.0 +/- 1.0 and 3.0 +/- 1.0 Myr for LH9 and LH10, respectively. The age difference and relative distance of the associations are consistent with a sequential star formation scenario in which stellar activity in LH9 triggered the formation of LH10. Our sample contains four stars of spectral type O2. From helium and hydrogen line fitting we find the hottest three of these stars to be similar to 49- 54 kK ( compared to similar to 45- 46 kK for O3 stars). Detailed determination of the helium mass fraction reveals that the masses of helium enriched dwarfs and giants derived in our spectroscopic analysis are systematically lower than those implied by non-rotating evolutionary tracks. We interpret this as evidence for efficient rotationally enhanced mixing leading to the surfacing of primary helium and to an increase of the stellar luminosity. This result is consistent with findings for SMC stars by Mokiem et al. ( 2006). For bright giants and supergiants no such mass discrepancy is found; these stars therefore appear to follow tracks of modestly or non-rotating objects. The set of programme stars was sufficiently large to establish the mass loss rates of OB stars in this Z similar to 1/2 Z(circle dot) environment sufficiently accurate to allow for a quantitative comparison with similar objects in the Galaxy and the SMC. The mass loss properties are found to be intermediate to massive stars in the Galaxy and SMC. Comparing the derived modified wind momenta D-mom as a function of luminosity with predictions for LMC metallicities by Vink et al. ( 2001) yields good agreement in the entire luminosity range that was investigated, i.e. 5.0

The VLT-FLAMES survey of massive stars: observations centered on the Magellanic Cloud clusters NGC 330, NGC 346, NGC 2004, and the N11 region

We present new observations of 470 stars using the Fibre Large Array Multi-Element Spectrograph ( FLAMES) instrument in fields centered on the clusters NGC330 and NGC346 in the Small Magellanic Cloud (SMC), and NGC2004 and the N11 region in the Large Magellanic Cloud (LMC). A further 14 stars were observed in the N11 and NGC330 fields using the Ultraviolet and Visual Echelle Spectrograph (UVES) for a separate programme. Spectral classifications and stellar radial velocities are given for each target, with careful attention to checks for binarity. In particular, we have investigated previously unexplored regions around the central LH9/LH10 complex of N11, finding similar to 25 new O-type stars from our spectroscopy. We have observed a relatively large number of Be-type stars that display permitted Fe II emission lines. These are primarily not in the cluster cores and appear to be associated with classical Be-type stars, rather than pre main-sequence objects. The presence of the Fe II emission, as compared to the equivalent width of Ha, is not obviously dependent on metallicity. We have also explored the relative fraction of Be- to normal B-type stars in the field-regions near to NGC330 and NGC2004, finding no strong evidence of a trend with metallicity when compared to Galactic results. A consequence of service observations is that we have reasonable time-sampling in three of our FLAMES fields. We find lower limits to the binary fraction of O- and early B-type stars of 23 to 36%. One of our targets (NGC346-013) is especially interesting with a massive, apparently hotter, less luminous secondary component.

The VLT-FLAMES survey of massive stars: NGC 346-013 as a test case for massive close binary evolution

Context. NGC 346-013 is a peculiar double-lined eclipsing binary in the Small Magellanic Cloud (SMC) discovered by the VLT-FLAMES survey of massive stars.

High-contrast imaging search for planets and brown dwarfs around the most massive stars in the solar neighborhood

There has been a long-standing discussion in the literature as to whether core accretion or disk instability is the dominant mode of planet formation. Over the last decade, several lines of evidence have been presented showing that core accretion is most likely the dominant mechanism for the close-in population of planets probed by radial velocity and transits. However, this does not by itself prove that core accretion is the dominant mode for the total planet population, since disk instability might conceivably produce and retain large numbers of planets in the far-out regions of the disk. If this is a relevant scenario, then the outer massive disks of B-stars should be among the best places for massive planets and brown dwarfs to form and reside. In this study, we present high-contrast imaging of 18 nearby massive stars of which 15 are in the B2-A0 spectral-type range and provide excellent sensitivity to wide companions. By comparing our sensitivities to model predictions of disk instability based on physical criteria for fragmentation and cooling, and using Monte Carlo simulations for orbital distributions, we find that ~85% of such companions should have been detected in our images on average. Given this high degree of completeness, stringent statistical limits can be set from the null-detection result, even with the limited sample size. We find that

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.

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.

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

GEMINI near-infrared spectroscopic observations of young massive stars embedded in molecular clouds

K-band spectra of young stellar candidates in four Southern hemisphere clusters have been obtained with the Gemini Near-Infrared Spectrograph in Gemini South. The clusters are associated with IRAS sources that have colours characteristic of ultracompact H II regions. Spectral types were obtained by comparison of the observed spectra with those of a near-infrared (NIR) library; the results include the spectral classification of nine massive stars and seven objects confirmed as background late-type stars. Two of the studied sources have K-band spectra compatible with those characteristic of very hot stars, as inferred from the presence of C IV, N III and N V emission lines at 2.078, 2.116 and 2.100 mu m, respectively. One of them, I16177_IRS1, has a K-band spectrum similar to that of Cyg OB2 7, an O3If* supergiant star. The nebular K-band spectrum of the associated Ultra-Compact (UC) H II region shows the s-process [Kr III] and [Se IV] high excitation emission lines, previously identified only in planetary nebula. One young stellar object was found in each cluster, associated with either the main IRAS source or a nearby resolved Midecourse Space eXperiment (MSX) component, confirming the results obtained from previous NIR photometric surveys. The distances to the stars were derived from their spectral types and previously determined JHK magnitudes; they agree well with the values obtained from the kinematic method, except in the case of IRAS 15408-5356, for which the spectroscopic distance is about a factor of 2 smaller than the kinematic value.

Collective scattering of massive stars in higher order gravity

In the weak field approximation of higher order gravity theory a gravitational potential is described by a Newtonian plus a Yukawa-like term. This new term is used to explain some aspects of galactic dynamics, without considering dark matter. Its presence modifies the scattering probability of a massive intruder star and relaxation time of the stellar system.

3D MHD simulations of subsurface convection in OB stars

During their main sequence evolution, massive stars can develop convective regions very close to their surface. These regions are caused by an opacity peak associated with iron ionization. Cantiello et al. (2009) found a possible connection between the presence of sub-photospheric convective motions and small scale stochastic velocities in the photosphere of early-type stars. This supports a physical mechanism where microturbulence is caused by waves that are triggered by subsurface convection zones. They further suggest that clumping in the inner parts of the winds of OB stars could be related to subsurface convection, and that the convective layers may also be responsible for stochastic excitation of non-radial pulsations. Furthermore, magnetic fields produced in the iron convection zone could appear at the surface of such massive stars. Therefore subsurface convection could be responsible for the occurrence of observable phenomena such as line profile variability and discrete absorption components. These phenomena have been observed for decades, but still evade a clear theoretical explanation. Here we present preliminary results from 3D MHD simulations of such subsurface convection.

Dynamical stability of nascent neutron stars

This thesis presents a study of the dynamical stability of nascent neutron stars resulting from the accretion induced collapse of rapidly rotating white dwarfs.

Chapter 2 and part of Chapter 3 study the equilibrium models for these neutron stars. They are constructed by assuming that the neutron stars have the same masses, angular momenta, and specific angular momentum distributions as the pre-collapse white dwarfs. If the pre-collapse white dwarf is rapidly rotating, the collapsed object will contain a high density central core of size about 20 km, surrounded by a massive accretion torus extending to hundreds of kilometers from the rotation axis. The ratio of the rotational kinetic energy to gravitational binding energy, β, of these neutron stars is all found to be less than 0.27.

Chapter 3 studies the dynamical stability of these neutron stars by numerically evolving the linearized hydrodynamical equations. A dynamical bar-mode instability is observed when the β of the star is greater than the critical value β_d ≈ 0.25. It is expected that the unstable mode will persist until a substantial amount of angular momentum is carried away by gravitational radiation. The detectability of these sources is studied and it is estimated that LIGO II is unlikely to detect them unless the event rate is greater than 10^-6/year/galaxy.

All the calculations on the structure and stability of the neutron stars in Chapters 2 and 3 are carried out using Newtonian hydrodynamics and gravity. Chapter 4 studies the relativistic effects on the structure of these neutron stars. New techniques are developed and used to construct neutron star models to the first post-Newtonian (1PN) order. The structures of the 1PN models are qualitatively similar to the corresponding Newtonian models, but the values of β are somewhat smaller. The maximum β for these 1PN neutron stars is found to be 0.24, which is 8% smaller than the Newtonian result (0.26). However, relativistic effects will also change the critical value β_d. A detailed post-Newtonian stability analysis has yet to be carried out to study the relativistic effects on the dynamical stability of these neutron stars.