J. Maíz Apellániz

The VLT-FLAMES Tarantula Survey VIII. Multiplicity properties of the O-type star population

Context. The Tarantula Nebula in the Large Magellanic Cloud is our closest view of a starburst region and is the ideal environment to investigate important questions regarding the formation, evolution and final fate of the most massive stars. Aims. We analyze the multiplicity properties of the massive O-type star population observed through multi-epoch spectroscopy in the framework of the VLT-FLAMES Tarantula Survey. With 360 O-type stars, this is the largest homogeneous sample of massive stars analyzed to date. Methods. We use multi-epoch spectroscopy and variability analysis to identify spectroscopic binaries. We also use a Monte-Carlo method to correct for observational biases. By modeling simultaneously the observed binary fraction, the distributions of the amplitudes of the radial velocity variations and the distribution of the time scales of these variations, we constrain the intrinsic current binary fraction and period and mass-ratio distributions. Results. We observe a spectroscopic binary fraction of 0.35 ± 0.03, which corresponds to the fraction of objects displaying statistically significant radial velocity variations with an amplitude of at least 20 km s-1. We compute the intrinsic binary fraction to be 0.51 ± 0.04. We adopt power-laws to describe the intrinsic period and mass-ratio distributions: f(log 10P/d) ~ (log 10P/d)π (with log 10P/d in the range 0.15−3.5) and f(q) ~ qκ with 0.1 ≤ q = M2/M1 ≤ 1.0. The power-law indexes that best reproduce the observed quantities are π = −0.45 ± 0.30 and κ = −1.0 ± 0.4. The period distribution that we obtain thus favours shorter period systems compared to an Opik law (π = 0). The mass ratio distribution is slightly skewed towards low mass ratio systems but remains incompatible with a random sampling of a classical mass function (κ = −2.35). The binary fraction seems mostly uniform across the field of view and independent of the spectral types and luminosity classes. The binary fraction in the outer region of the field of view (r > 7.8′, i.e. ≈117 pc) and among the O9.7 I/II objects are however significantly lower than expected from statistical fluctuations. The observed and intrinsic binary fractions are also lower for the faintest objects in our sample (Ks > 15.5 mag), which results from observational effects and the fact that our O star sample is not magnitude-limited but is defined by a spectral-type cutoff. We also conclude that magnitude-limited investigations are biased towards larger binary fractions. Conclusions. Using the multiplicity properties of the O stars in the Tarantula region and simple evolutionary considerations, we estimate that over 50% of the current O star population will exchange mass with its companion within a binary system. This shows that binary interaction is greatly affecting the evolution and fate of massive stars, and must be taken into account to correctly interpret unresolved populations of massive stars.

A magnetic confinement versus rotation classification of massive-star magnetospheres

Building on results from the Magnetism in Massive Stars (MiMeS) project, this paper shows how a two-parameter classification of massive-star magnetospheres in terms of the magnetic wind confinement (which sets the Alfv´ en radius RA) and stellar rotation (which sets the Kepler co-rotation radius RK) provides a useful organization of both observational signatures and theoretical predictions. We compile the first comprehensive study of inferred and observed values for relevant stellar and magnetic parameters of 64 confirmed magnetic OB stars with Teff 16 kK. Using these parameters, we locate the stars in the magnetic confinement–rotation diagram, a log–log plot of RK versus RA. This diagram can be subdivided into regimes of centrifugal magnetospheres (CM), with RA > RK ,v ersusdynamical magnetospheres (DM), with RK > RA. We show how key observational diagnostics, like the presence and characteristics of Hα emission, depend on a star’s position within the diagram, as well as other parameters, especially the expected wind mass-loss rates. In particular, we identify two distinct populations of magnetic stars with Hα emission: namely, slowly rotating O-type stars with narrow emission consistent with a DM, and more rapidly rotating B-type stars with broader emission associated with a CM. For O-type stars, the high mass-loss rates are sufficient to accumulate enough material for line emission even within the relatively short free-fall time-scale associated with a DM: this high mass-loss rate also leads to a rapid magnetic spindown of the stellar rotation. For the B-type stars, the longer confinement of a CM is required to accumulate sufficient emitting material from their relatively weak winds, which also lead to much longer spindown time-scales. Finally, we discuss how other observational diagnostics, e.g. variability of UV wind lines or X-ray emission, relate to the inferred magnetic properties of these stars, and summarize prospects for future developments in our understanding of massive-star magnetospheres.

The VLT-FLAMES Tarantula Survey - XII. Rotational velocities of the single O-type stars

Context. The 30 Doradus (30 Dor) region of the Large Magellanic Cloud, also known as the Tarantula nebula, is the nearest starburst region. It contains the richest population of massive stars in the Local Group, and it is thus the best possible laboratory to investigate open questions on the formation and evolution of massive stars. Aims. Using ground-based multi-object optical spectroscopy obtained in the framework of the VLT-FLAMES Tarantula Survey (VFTS), we aim to establish the (projected) rotational velocity distribution for a sample of 216 presumably single O-type stars in 30 Dor. The sample is large enough to obtain statistically significant information and to search for variations among subpopulations – in terms of spectral type, luminosity class, and spatial location – in the field of view. Methods. We measured projected rotational velocities, ν_esini, by means of a Fourier transform method and a profile fitting method applied to a set of isolated spectral lines. We also used an iterative deconvolution procedure to infer the probability density, P(ν_e), of the equatorial rotational velocity, ν_e. Results. The distribution of νesini shows a two-component structure: a peak around 80 kms^(-1) and a high-velocity tail extending up to ~600 kms^(-1). This structure is also present in the inferred distribution P(ν_e) with around 80% of the sample having 0 < ν_e ≤ 300 kms^(-1) and the other 20% distributed in the high-velocity region. The presence of the low-velocity peak is consistent with what has been found in other studies for late O- and early B-type stars. Conclusions. Most of the stars in our sample rotate with a rate less than 20% of their break-up velocity. For the bulk of the sample, mass loss in a stellar wind and/or envelope expansion is not efficient enough to significantly spin down these stars within the first few Myr of evolution. If massive-star formation results in stars rotating at birth with a large portion of their break-up velocities, an alternative braking mechanism, possibly magnetic fields, is thus required to explain the present-day rotational properties of the O-type stars in 30 Dor. The presence of a sizeable population of fast rotators is compatible with recent population synthesis computations that investigate the influence of binary evolution on the rotation rate of massive stars. Even though we have excluded stars that show significant radial velocity variations, our sample may have remained contaminated by post-interaction binary products. That the high-velocity tail may be populated primarily (and perhaps exclusively) by post-binary interaction products has important implications for the evolutionary origin of systems that produce gamma-ray bursts.

The VLT-FLAMES Tarantula Survey. III. A very massive star in apparent isolation from the massive cluster R136

VFTS 682 is located in an active star-forming region, at a projected distance of 29 pc from the young massive cluster R136 in the Tarantula Nebula of the Large Magellanic Cloud. It was previously reported as a candidate young stellar object, and more recently spectroscopically revealed as a hydrogen-rich Wolf-Rayet (WN5h) star. Our aim is to obtain the stellar properties, such as its intrinsic luminosity, and to investigate the origin of VFTS 682. To this purpose, we model optical spectra from the VLT-FLAMES Tarantula Survey with the non-LTE stellar atmosphere code cmfgen, as well as the spectral energy distribution from complementary optical and infrared photometry. We find the extinction properties to be highly peculiar (RV ∼ 4.7), and obtain a surprisingly high luminosity log(L/L� ) = 6.5 ± 0.2, corresponding to a present-day mass of ∼150 M� . The high effective temperature of 52.2 ± 2.5 kK might be explained by chemically homogeneous evolution – suggested to be the key process in the path towards long gamma-ray bursts. Lightcurves of the object show variability at the 10% level on a timescale of years. Such changes are unprecedented for classical WolfRayet stars, and are more reminiscent of Luminous Blue Variables. Finally, we discuss two possibilities for the origin of VFTS 682: (i) the star either formed in situ, which would have profound implications for the formation mechanism of massive stars, or (ii) VFTS 682 is a slow runaway star that originated from the dense cluster R136, which would make it the most massive runaway known to date.

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 Minit > 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-1. 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+55-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 VLT-FLAMES Tarantula Survey. XIV. The O-type stellar content of 30 Doradus

Detailed spectral classifications are presented for 352 O-B0 stars in the VLT-FLAMES Tarantula Survey ESO Large Programme, of which 213 O-type are judged of sufficiently high quality for further morphological analysis. Among them, six subcategories of special interest are distinguished. (1) Several new examples of the earliest spectral types O2-O3 have been found, while a previously known example has been determined to belong to the nitrogen-rich ON2 class. (2) A group of extremely rapidly rotating main-sequence objects has been isolated, including the largest vsini values known, the spatial and radial-velocity distributions of which suggest ejection from the two principal ionizing clusters NGC 2070 and NGC 2060. (3) Several new examples of the evolved, rapidly rotating Onfp class show similar evidence, although at least some of them are spectroscopic binaries. (4) No fewer than 48 members of the Vz category, hypothesized to be on or near the zero-age main sequence, are found in this sample; in contrast to the rapid rotators, they are strongly concentrated to the ionizing clusters and a newly recognized region of current and recent star formation to the north, supporting their interpretation as very young objects, as do their relatively faint absolute magnitudes. (5) A surprisingly large fraction of the main-sequence spectra belong to the recently recognized V((fc)) class, with C iii emission lines of similar strength to the usual N iii in V((f)) spectra, although a comparable number of the latter are also present, as well as six objects with very high-quality data but no trace of either emission feature, presenting new challenges to physical interpretations. (6) Two mid-O Vz and three late-O giant/supergiant spectra with morphologically enhanced nitrogen lines have been detected. Absolute visual magnitudes have been derived for each star with individual extinction laws, and composite Hertzsprung-Russell diagrams provide evidence of the multiple generations present in this field. Spectroscopic binaries, resolved visual multiples, and possible associations with X-ray sources are noted. Astrophysical and dynamical analyses of this unique dataset underway will provide new insights into the evolution of massive stars and starburst clusters.

The R136 star cluster dissected with Hubble Space Telescope/STIS. I. Far-ultraviolet spectroscopic census and the origin of He ii λ1640 in young star clusters

We introduce a Hubble Space Telescope (HST)/Space Telescope Imaging Spectrograph (STIS) stellar census of R136a, the central ionizing star cluster of 30 Doradus. We present low resolution far-ultraviolet STIS spectroscopy of R136 using 17 contiguous 52 arcsec x 0.2 arcsec slits which together provide complete coverage of the central 0.85 parsec (3.4 arcsec). We provide spectral types of 90 per cent of the 57 sources brighter than m(F555W) = 16.0 mag within a radius of 0.5 parsec of R136a1, plus 8 additional nearby sources including R136b (O4 If/WN8). We measure wind velocities for 52 early-type stars from C IV lambda lambda 1548-51, including 16 O2-3 stars. For the first time, we spectroscopically classify all Weigelt and Baier members of R136a, which comprise three WN5 stars (a1-a3), two O supergiants (a5-a6) and three early O dwarfs (a4, a7, a8). A complete Hertzsprung-Russell diagram for the most massive O stars in R136 is provided, from which we obtain a cluster age of 1.5(-0.7)(+0.3) Myr. In addition, we discuss the integrated ultraviolet spectrum of R136, and highlight the central role played by the most luminous stars in producing the prominent He II lambda 1640 emission line. This emission is totally dominated by very massive stars with initial masses above similar to 100M(circle dot). The presence of strong He II lambda 1640 emission in the integrated light of very young star clusters (e.g. A1 in NGC 3125) favours an initial mass function extending well beyond a conventional upper limit of 100M(circle dot). We include montages of ultraviolet spectroscopy for Large Magellanic Cloud O stars in the appendix. Future studies in this series will focus on optical STIS medium resolution observations.

A DOUBLE CLUSTER AT THE CORE OF 30 DORADUS

Based on an analysis of data obtained with the Wide Field Camera 3 on the Hubble Space Telescope we report the identification of two distinct stellar populations in the core of the giant H II region 30 Doradus in the Large Magellanic Cloud. The most compact and richest component coincides with the center of R136 and is ~1 Myr younger than a second more diffuse clump, located ~5.4 pc toward the northeast. We note that published spectral types of massive stars in these two clumps lend support to the proposed age difference. The morphology and age difference between the two sub-clusters suggests that an ongoing merger may be occurring within the core of 30 Doradus. This finding is consistent with the predictions of models of hierarchical fragmentation of turbulent giant molecular clouds, according to which star clusters would be the final products of merging smaller sub-structures.

Pismis 24-1 : The stellar upper mass limit preserved

Is there a stellar upper mass limit? Recent statistical work seems to indicate that there is and that it is in the vicinity of 150 M☉. In this paper we use HST and ground-based data to investigate the brightest members of the cluster Pismis 24, one of which (Pismis 24-1) was previously inferred to have a mass greater than 200 M☉, in apparent disagreement with that limit. We determine that Pismis 24-1 is composed of at least three objects, the resolved Pismis 24-1SW and the unresolved spectroscopic binary Pismis 24-1NE. The evolutionary zero-age masses of Pismis 24-1SW, the unresolved system Pismis 24-1NE, and the nearby star Pismis 24-17 are all ≈100 M☉, very large but under the stellar upper mass limit.

NGC 1624-2: a slowly rotating, X-ray luminous Of?cp star with an extraordinarily strong magnetic field

This paper presents a first observational investigation of t he faint Of?p star NGC 1624-2, yielding important new constraints on its spectral and physical characteristics, rotation, magnetic field strength, X-ray emission and magnetospheric pro perties. Modeling the spectrum and spectral energy distribution, we conclude that NGC 1624-2 is a main sequence star of mass M≃ 30 M⊙, and infer an effective temperature of 35± 2 kK and log g = 4.0± 0.2. Based on an extensive time series of optical spectral observation s we report significant variability of a large number of spectral lines, and infer a unique period of 157.99± 0.94 d which we interpret as the rotational period of the star. We report the detec tion of a very strong - 5.35± 0.5 kG - longitudinal magnetic field h Bzi , coupled with probable Zeeman splitting of Stokes I profiles of metal lines confirming a surface field modulus h Bi of 14± 1 kG, consistent with a surface dipole of polar strength ∼ 20 kG. This is the largest magnetic field ever detected in an O-type star, and the first report of Zeeman splitting of Stoke s I profiles in such an object. We also report the detection of reversed Stokes V profiles associated with weak, high-excitation emission lines of Oiii, which we propose may form in the close magnetosphere of the star. We analyze archival Chandra ACIS-I X-ray data, inferring a very hard spectrum with an X-ray effi ciency log Lx/Lbol =−6.4, a factor of 4 larger than the canonical value for O-type sta rs and comparable to that of the young magnetic O-type starθ 1 Ori C and other Of?p stars. Finally, we examine the probable magnetospheric properties of the star, reporting in particular very strong magnetic confinement of the stellar wind, with η∗≃ 1.5× 10 4 , and a very large Alfven radius, RAlf = 11.4 R∗.