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

We present the results of a 10.5-yr, volume-limited (28-Mpc) search for supernova (SN) progenitor stars. In doing so we compile all SNe discovered within this volume (132, of which 27 per cent are Type Ia) and determine the relative rates of each subtype from literature studies. The core-collapse SNe break down into 59 per cent II-P and 29 per cent Ib/c, with the remainder being IIb (5 per cent), IIn (4 per cent) and II-L (3 per cent). There have been 20 II-P SNe with high-quality optical or near-infrared pre-explosion images that allow a meaningful search for the progenitor stars. In five cases they are clearly red supergiants, one case is unconstrained, two fall on compact coeval star clusters and the other twelve have no progenitor detected. We review and update all the available data for the host galaxies and SN environments (distance, metallicity and extinction) and determine masses and upper mass estimates for these 20 progenitor stars using the STARS stellar evolutionary code and a single consistent homogeneous method. A maximum likelihood calculation suggests that the minimum stellar mass for a Type II-P to form is m(min) = 8.5(-1.5)(+1) M-circle dot and the maximum mass for II-P progenitors is m(max) = 16.5 +/- 1.5 M-circle dot, assuming a Salpeter initial mass function holds for the progenitor population (in the range Gamma = -1.35(-0.7)(+0.3)). The minimum mass is consistent with current estimates for the upper limit to white dwarf progenitor masses, but the maximum mass does not appear consistent with massive star populations in Local Group galaxies. Red supergiants in the Local Group have masses up to 25 M-circle dot and the minimum mass to produce a Wolf-Rayet star in single star evolution (between solar and LMC metallicity) is similarly 25-30 M-circle dot. The reason we have not detected any high-mass red supergiant progenitors above 17 M-circle dot is unclear, but we estimate that it is statistically significant at 2.4 sigma confidence. Two simple reasons for this could be that we have systematically underestimated the progenitor masses due to dust extinction or that stars between 17-25 M-circle dot produce other kinds of SNe which are not II-P. We discuss these possibilities and find that neither provides a satisfactory solution. We term this discrepancy the 'red supergiant problem' and speculate that these stars could have core masses high enough to form black holes and SNe which are too faint to have been detected. We compare the Ni-56 masses ejected in the SNe to the progenitor mass estimates and find that low-luminosity SNe with low Ni-56 production are most likely to arise from explosions of low-mass progenitors near the mass threshold that can produce a core-collapse.

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 Tarantula Survey XI. A census of the hot luminous stars and their feedback in 30 Doradus

Context. The VLT-FLAMES Tarantula Survey has an extensive view of the copious number of massive stars in the 30 Doradus (30 Dor) star forming region of the Large Magellanic Cloud. These stars play a crucial role in our understanding of the stellar feedback in more distant, unresolved star forming regions. Aims. The first comprehensive census of hot luminous stars in 30 Dor is compiled within a 10 arcmin (150 pc) radius of its central cluster, R136. We investigate the stellar content and spectroscopic completeness of the early type stars. Estimates were made for both the integrated ionising luminosity and stellar wind luminosity. These values were used to re-assess the star formation rate (SFR) of the region and determine the ionising photon escape fraction. Methods. Stars were selected photometrically and combined with the latest spectral classifications. Spectral types were estimated for stars lacking spectroscopy and corrections were made for binary systems, where possible. Stellar calibrations were applied to obtain their physical parameters and wind properties. Their integrated properties were then compared to global observations from ultraviolet (UV) to far-infrared (FIR) imaging as well as the population synthesis code, Starburst99. Results. Our census identified 1145 candidate hot luminous stars within 150 pc of R136 of which >700 were considered to be genuine early type stars and contribute to feedback. We assess the survey to be spectroscopically complete to 85% in the outer regions (>5 pc) but only 35% complete in the region of the R136 cluster, giving a total of 500 hot luminous stars in the census which had spectroscopy. Only 31 were found to be Wolf-Rayet (W-R) or Of/WN stars, but their contribution to the integrated ionising luminosity and wind luminosity was ~ 40% and ~ 50%, respectively. Similarly, stars with M > 100 M (mostly H-rich WN stars) also showed high contributions to the global feedback, ~ 25% in both cases. Such massive stars are not accounted for by the current Starburst99 code, which was found to underestimate the integrated ionising luminosity of R136 by a factor ~ 2 and the wind luminosity by a factor ~ 9. The census inferred a SFR for 30 Dor of 0.073 ± 0.04 M yr . This was generally higher than that obtained from some popular SFR calibrations but still showed good consistency with the far-UV luminosity tracer as well as the combined Hα and mid-infrared tracer, but only after correcting for Hα extinction. The global ionising output was also found to exceed that measured from the associated gas and dust, suggesting that ~6 % of the ionising photons escape the region. Conclusions. When studying the most luminous star forming regions, it is essential to include their most massive stars if one is to determine a reliable energy budget. Photon leakage becomes more likely after including their large contributions to the ionising output. If 30 Dor is typical of other massive star forming regions, estimates of the SFR will be underpredicted if this escape fraction is not accounted for.

The Type IIb SN 2011dh: Two years of observations and modelling of the lightcurves

We present optical and near-infrared (NIR) photometry and spectroscopy as well as modelling of the lightcurves of the Type IIb supernova (SN) 2011dh. Our extensive dataset, for which we present the observations obtained after day 100, spans two years, and complemented with Spitzer mid-infrared (MIR) data, we use it to build an optical-to-MIR bolometric lightcurve between days 3 and 732. To model the bolometric lightcurve before day 400 we use a grid of hydrodynamical SN models, which allows us to determine the errors in the derived quantities, and a bolometric correction determined with steady-state non-local thermodynamic equilibrium (NLTE) modelling. Using this method we find a helium core mass of 3.1^+0.7<inf>-0.4</inf> M<inf>⊙</inf> for SN 2011dh, consistent within error bars with previous results obtained using the bolometric lightcurve before day 80. We compute bolometric and broad-band lightcurves between days 100 and 500 from spectral steady-state NLTE models, presented and discussed in a companion paper. The preferred 12 M<inf>⊙</inf> (initial mass) model, previously found to agree well with the observed spectra, shows a good overall agreement with the observed lightcurves, although some discrepancies exist. Time-dependent NLTE modelling shows that after day ∼600 a steady-state assumption is no longer valid. The radioactive energy deposition in this phase is likely dominated by the positrons emitted in the decay of ⁵⁶Co, but seems insufficient to reproduce the lightcurves, and what energy source is dominating the emitted flux is unclear. We find an excess in the K and the MIR bands developing between days 100 and 250, during which an increase in the optical decline rate is also observed. A local origin of the excess is suggested by the depth of the He I 20 581 Å absorption. Steady-state NLTE models with a modest dust opacity in the core (τ = 0.44), turned on during this period, reproduce the observed behaviour, but an additional excess in the Spitzer 4.5 μm band remains. Carbon-monoxide (CO) first-overtone band emission is detected at day 206, and possibly at day 89, and assuming the additional excess to bedominated by CO fundamental band emission, we find fundamental to first-overtone band ratios considerably higher than observed in SN 1987A. The profiles of the [O i] 6300 Å and Mg i] 4571 Å lines show a remarkable similarity, suggesting that these lines originate from a common nuclear burning zone (O/Ne/Mg), and using small scale fluctuations in the line profiles we estimate a filling factor of ≲ 0.07 for the emitting material. This paper concludes our extensive observational and modelling work on SN 2011dh. The results from hydrodynamical modelling, steady-state NLTE modelling, and stellar evolutionary progenitor analysis are all consistent, and suggest an initial mass of ∼12 M<inf>⊙</inf> for the progenitor.

A FOSSIL BULGE GLOBULAR CLUSTER REVEALED BY VERY LARGE TELESCOPE MULTI-CONJUGATE ADAPTIVE OPTICS

The globular cluster HP 1 is projected on the bulge, very close to the Galactic center. The Multi-Conjugate Adaptive Optics Demonstrator on the Very Large Telescope allowed us to acquire high-resolution deep images that, combined with first epoch New Technology Telescope data, enabled us to derive accurate proper motions. The cluster and bulge fields` stellar contents were disentangled through this process and produced an unprecedented definition in color-magnitude diagrams of this cluster. The metallicity of [Fe/H] approximate to -1.0 from previous spectroscopic analysis is confirmed, which together with an extended blue horizontal branch imply an age older than the halo average. Orbit reconstruction results suggest that HP 1 is spatially confined within the bulge.

Multi-Conjugate Adaptive Optics VLT imaging of the distant old open cluster FSR 1415

We employ the recently installed near-infrared Multi-Conjugate Adaptive Optics demonstrator (MAD) to determine the basic properties of a newly identified, old and distant, Galactic open cluster (FSR 1415). The MAD facility remarkably approaches the diffraction limit, reaching a resolution of 0.07 arcsec (in K), that is also uniform in a field of similar to 1.8 arcmin in diameter. The MAD facility provides photometry that is 50 per cent complete at K similar to 19. This corresponds to about 2.5 mag below the cluster main-sequence turn-off. This high-quality data set allows us to derive an accurate heliocentric distance of 8.6 kpc, a metallicity close to solar and an age of similar to 2.5 Gyr. On the other hand, the deepness of the data allows us to reconstruct (completeness-corrected) mass functions (MFs) indicating a relatively massive cluster, with a flat core MF. The Very Large Telescope/MAD capabilities will therefore provide fundamental data for identifying/analysing other faint and distant open clusters in the Galaxy III and IV quadrants.

Transiting exoplanets from the CoRoT space mission XX. CoRoT-20b: A very high density, high eccentricity transiting giant planet

We report the discovery by the CoRoT space mission of a new giant planet, CoRoT-20b. The planet has a mass of 4.24 +/- 0.23 M-Jup and a radius of 0.84 +/- 0.04 R-Jup. With a mean density of 8.87 +/- 1.10 g cm(-3), it is among the most compact planets known so far. Evolutionary models for the planet suggest a mass of heavy elements of the order of 800 M-circle plus if embedded in a central core, requiring a revision either of the planet formation models or both planet evolution and structure models. We note however that smaller amounts of heavy elements are expected by more realistic models in which they are mixed throughout the envelope. The planet orbits a G-type star with an orbital period of 9.24 days and an eccentricity of 0.56. The star's projected rotational velocity is v sin i = 4.5 +/- 1.0 km s(-1), corresponding to a spin period of 11.5 +/- 3.1 days if its axis of rotation is perpendicular to the orbital plane. In the framework of Darwinian theories and neglecting stellar magnetic breaking, we calculate the tidal evolution of the system and show that CoRoT-20b is presently one of the very few Darwin-stable planets that is evolving toward a triple synchronous state with equality of the orbital, planetary and stellar spin periods.

From planetesimals to planets: volatile molecules

Context. Solar and extrasolar planets are the subject of numerous studies aiming to determine their chemical composition and internal structure. In the case of extrasolar planets, the composition is important as it partly governs their potential habitability. Moreover, observational determination of chemical composition of planetary atmospheres are becoming available, especially for transiting planets. Aims. The present works aims at determining the chemical composition of planets formed in stellar systems of solar chemical composition. The main objective of this work is to provide valuable theoretical data for models of planet formation and evolution, and future interpretation of chemical composition of solar and extrasolar planets. Methods. We have developed a model that computes the composition of ices in planets in different stellar systems with the use of models of ice and planetary formation. Results. We provide the chemical composition, ice/rock mass ratio and C:O molar ratio for planets in stellar systems of solar chemical composition. From an initial homogeneous composition of the nebula, we produce a wide variety of planetary chemical compositions as a function of the mass of the disk and distance to the star. The volatile species incorporated in planets are mainly composed of H2O, CO, CO2, CH3OH, and NH3. Icy or ocean planets have systematically higher values of molecular abundances compared to giant and rocky planets. Gas giant planets are depleted in highly volatile molecules such as CH4, CO, and N2 compared to icy or ocean planets. The ice/rock mass ratio in icy or ocean and gas giant planets is, respectively, equal at maximum to 1.01 ± 0.33 and 0.8 ± 0.5, and is different from the usual assumptions made in planet formation models, which suggested this ratio to be 2–3. The C:O molar ratio in the atmosphere of gas giant planets is depleted by at least 30% compared to solar value.

Helioseismology: Probing the interior of a star

Helioseismology offers, for the first time, an opportunity to probe in detail the deep interior of a star (our Sun). The results will have a profound impact on our understanding not only of the solar interior, but also neutrino physics, stellar evolution theory, and stellar population studies in astrophysics.

Isotopes in cosmochemistry: recipe for a Solar System

Extreme isotopic variations among extraterrestrial materials provide great insights into the origin and evolution of the Solar System. In this tutorial review, we summarize how the measurement of isotope ratios can expand our knowledge of the processes that took place before and during the formation of our Solar System and its subsequent early evolution. The continuous improvement of mass spectrometers with high precision and increased spatial resolution, including secondary ion mass spectrometry (SIMS), thermal ionization mass spectrometry (TIMS) and multi collector-inductively coupled plasma-mass spectrometry (MC-ICP-MS), along with the ever growing amounts of available extraterrestrial samples have significantly increased the temporal and spatial constraints on the sequence of events that took place since and before the formation of the first Solar System condensates (i.e., Ca-Al-rich inclusions). Grains sampling distinct stellar environments with a wide range of isotopic compositions were admixed to, but possibly not fully homogenized in, the Sun's parent molecular cloud or the nascent Solar System. Before, during and after accretion of the nebula, as well as the formation and subsequent evolution of planetesimals and planets, chemical and physical fractionation processes irrevocably changed the chemical and isotopic compositions of all Solar System bodies. Since the formation of the first Solar System minerals and rocks 4.568 Gyr ago, short-and long-lived radioactive decay and cosmic ray interaction also contributed to the modification of the isotopic framework of the Solar System, and permit to trace the formation and evolution of directly accessible and inferred planetary and stellar isotopic reservoirs.

COST AND PERFORMANCE TRADEOFFS OF ALTERNATIVE SOLAR-DRIVEN S-CO2 BRAYTON CYCLE CONFIGURATIONS

This paper evaluates cost and performance tradeoffs of alternative supercritical carbon dioxide (s-CO2) closed-loop Brayton cycle configurations with a concentrated solar heat source. Alternative s-CO2 power cycle configurations include simple, recompression, cascaded, and partial cooling cycles. Results show that the simple closed-loop Brayton cycle yielded the lowest power-block component costs while allowing variable temperature differentials across the s-CO2 heating source, depending on the level of recuperation. Lower temperature differentials led to higher sensible storage costs, but cycle configurations with lower temperature differentials (higher recuperation) yielded higher cycle efficiencies and lower solar collector and receiver costs. The cycles with higher efficiencies (simple recuperated, recompression, and partial cooling) yielded the lowest overall solar and power-block component costs for a prescribed power output.

Spin-orbit angle measurements for six southern transiting planets. New insights into the dynamical origins of hot Jupiters

Context. Several competing scenarios for planetary-system formation and evolution seek to explain how hot Jupiters came to be so close to their parent stars. Most planetary parameters evolve with time, making it hard to distinguish between models. The obliquity of an orbit with respect to the stellar rotation axis is thought to be more stable than other parameters such as eccentricity. Most planets, to date, appear aligned with the stellar rotation axis; the few misaligned planets so far detected are massive (> 2 MJ). Aims: Our goal is to measure the degree of alignment between planetary orbits and stellar spin axes, to search for potential correlations with eccentricity or other planetary parameters and to measure long term radial velocity variability indicating the presence of other bodies in the system. Methods: For transiting planets, the Rossiter-McLaughlin effect allows the measurement of the sky-projected angle ß between the stellar rotation axis and a planet's orbital axis. Using the HARPS spectrograph, we observed the Rossiter-McLaughlin effect for six transiting hot Jupiters found by the WASP consortium. We combine these with long term radial velocity measurements obtained with CORALIE. We used a combined analysis of photometry and radial velocities, fitting model parameters with the Markov Chain Monte Carlo method. After obtaining ß we attempt to statistically determine the distribution of the real spin-orbit angle ?. Results: We found that three of our targets have ß above 90°: WASP-2b: ß = 153°+11-15, WASP-15b: ß = 139.6°+5.2-4.3 and WASP-17b: ß = 148.5°+5.1-4.2; the other three (WASP-4b, WASP-5b and WASP-18b) have angles compatible with 0°. We find no dependence between the misaligned angle and planet mass nor with any other planetary parameter. All six orbits are close to circular, with only one firm detection of eccentricity e = 0.00848+0.00085-0.00095 in WASP-18b. No long-term radial acceleration was detected for any of the targets. Combining all previous 20 measurements of ß and our six and transforming them into a distribution of ? we find that between about 45 and 85% of hot Jupiters have ? > 30°. Conclusions: Most hot Jupiters are misaligned, with a large variety of spin-orbit angles. We find observations and predictions using the Kozai mechanism match well. If these observational facts are confirmed in the future, we may then conclude that most hot Jupiters are formed from a dynamical and tidal origin without the necessity to use type I or II migration. At present, standard disc migration cannot explain the observations without invoking at least another additional process.

Rotating massive main-sequence stars II. Simulating a population of LMC early B-type stars as a test of rotational mixing

Context. Rotational mixing in massive stars is a widely applied concept, with far-reaching consequences for stellar evolution, nucleosynthesis, and stellar explosions.

Log gf values for astrophysically important transitions Fe II

Aims. In a recent measurement, Meléndez & Barbuy (2009, A&A, 497, 611) report accurate log gf values for 142 important astrophysical lines with wavelengths in the range 4000 Å to 8000 Å. Their results include both solar and laboratory measurements. In this paper, we describe a theoretical study of these lines. Methods. The CIV3 structure codes, combined with our "fine-tuning" extrapolation process, are used to undertake a large-scale CI calculation involving the lowest 262 fine-structure levels belonging to the 3d4s, 3d, 3d4s, 3d4p, and 3d4s4p configurations. Results. We find that many of the 142 transitions are very weak intercombination lines. Other transitions are weak because the dominant configurations in the two levels differ by two orbitals. Conclusions. The comparison between our log gf values and the experimental values generally shows good agreement for most of these transitions, with our theoretical values agreeing slightly more closely with the solar than with the laboratory measurements. A detailed analysis of the small number of transitions for which the agreement between theory and experiment is not as good shows that such disagreements largely arise from severe cancellation due to CI mixing.

Intensity Enhancement of O VI Ultraviolet Emission Lines in Solar Spectra due to Opacity

Opacity is a property of many plasmas. It is normally expected that if an emission line in a plasma becomes optically thick, then its intensity ratio to that of another transition that remains optically thin should decrease. However, radiative transfer calculations undertaken both by ourselves and others predict that under certain conditions the intensity ratio of an optically thick to an optically thin line can show an increase over the optically thin value, indicating an enhancement in the former. These conditions include the geometry of the emitting plasma and its orientation to the observer. A similar effect can take place between lines of differing optical depths. While previous observational studies have focused on stellar point sources, here we investigate the spatially resolved solar atmosphere using measurements of the I(1032 Å)/I(1038 Å) intensity ratio of O VI in several regions obtained with the Solar Ultraviolet Measurements of Emitted Radiation instrument on board the Solar and Heliospheric Observatory satellite. We find several I(1032 Å)/I(1038 Å) ratios observed on the disk to be significantly larger than the optically thin value of 2.0, providing the first detection (to our knowledge) of intensity enhancement in the ratio arising from opacity effects in the solar atmosphere. The agreement between observation and theory is excellent and confirms that the O VI emission originates from a slab-like geometry in the solar atmosphere, rather than from cylindrical structures.