The VLT-FLAMES Tarantula Survey XXI. Stellar spin rates of O-type spectroscopic binaries

Context: The initial distribution of spin rates of massive stars is a fingerprint of their elusive formation process. It also sets a key initial condition for stellar evolution and is thus an important ingredient in stellar population synthesis. So far, most studies have focused on single stars. Most O stars are, however, found in multiple systems. 

Aims: By establishing the spin-rate distribution of a sizeable sample of O-type spectroscopic binaries and by comparing the distributions of binary subpopulations with one another and with that of presumed-single stars in the same region, we aim to constrain the initial spin distribution of O stars in binaries, and to identify signatures of the physical mechanisms that affect the evolution of the spin rates of massive stars. 

Methods: We use ground-based optical spectroscopy obtained in the framework of the VLT-FLAMES Tarantula Survey (VFTS) to establish the projected equatorial rotational velocities (ν_esini) for components of 114 spectroscopic binaries in 30 Doradus. The ν_esini values are derived from the full width at half maximum (FWHM) of a set of spectral lines, using a FWHM vs. ν_esini calibration that we derive based on previous line analysis methods applied to single O-type stars in the VFTS sample. 

Results: The overall ν_esini distribution of the primary stars resembles that of single O-type stars in the VFTS, featuring a low-velocity peak (at ν_esini<200 kms^-1) and a shoulder at intermediate velocities (200 <ν_esini<300 kms^-1). The distributions of binaries and single stars, however, differ in two ways. First, the main peak at ν_esini ~ 100kms^-1 is broader and slightly shifted towards higher spin rates in the binary distribution than that of the presumed-single stars. This shift is mostly due to short-period binaries (P_orb~<10 d). Second, the ν_esini distribution of primaries lacks a significant population of stars spinning faster than 300 kms^-1, while such a population is clearly present in the single-star sample. The ν_esini distribution of binaries with amplitudes of radial velocity variation in the range of 20 to 200 kms^-1 (mostly binaries with P_orb ~ 10-1000 d and/or with q<0.5) is similar to that of single O stars below ν_esini_~<170kms^-1. 

Conclusions: Our results are compatible with the assumption that binary components formed with the same spin distribution as single stars, and that this distribution contains few or no fast-spinning stars. The higher average spin rate of stars in short-period binaries may either be explained by spin-up through tides in such tight binary systems, or by spin-down of a fraction of the presumed-single stars and long-period binaries through magnetic braking (or by a combination of both mechanisms). Most primaries and secondaries of SB2 systems with P_orb~<10 d appear to have similar rotational velocities. This is in agreement with tidal locking in close binaries where the components have similar radii. The lack of very rapidly spinning stars among binary systems supports the idea that most stars with ν_esini_~> 300kms^-1 in the single-star sample are actually spun-up post-binary interaction products. Finally, the overall similarities (low-velocity peak and intermediate-velocity shoulder) of the spin distribution of binary and single stars argue for a massive star formation process in which the initial spin is set independently of whether stars are formed as single stars or as components of a binary system.

The VLT-FLAMES Tarantula Survey: XIX. B-type supergiants: Atmospheric parameters and nitrogen abundances to investigate the role of binarity and the width of the main sequence

Context: Model atmosphere analyses have been previously undertaken for both Galactic and extragalactic B-type supergiants. By contrast, little attention has been given to a comparison of the properties of single supergiants and those that are members of multiple systems. 

Aims: Atmospheric parameters and nitrogen abundances have been estimated for all the B-type supergiants identified in the VLT-FLAMES Tarantula survey. These include both single targets and binary candidates. The results have been analysed to investigate the role of binarity in the evolutionary history of supergiants. 

Methods: tlusty non-local thermodynamic equilibrium (LTE) model atmosphere calculations have been used to determine atmospheric parameters and nitrogen abundances for 34 single and 18 binary supergiants. Effective temperatures were deduced using the silicon balance technique, complemented by the helium ionisation in the hotter spectra. Surface gravities were estimated using Balmer line profiles and microturbulent velocities deduced using the silicon spectrum. Nitrogen abundances or upper limits were estimated from the Nii spectrum. The effects of a flux contribution from an unseen secondary were considered for the binary sample. Results. We present the first systematic study of the incidence of binarity for a sample of B-type supergiants across the theoretical terminal age main sequence (TAMS). To account for the distribution of effective temperatures of the B-type supergiants it may be necessary to extend the TAMS to lower temperatures. This is also consistent with the derived distribution of mass discrepancies, projected rotational velocities and nitrogen abundances, provided that stars cooler than this temperature are post-red supergiant objects. For all the supergiants in the Tarantula and in a previous FLAMES survey, the majority have small projected rotational velocities. The distribution peaks at about 50 km s^-1 with 65% in the range 30 km s^-1 ≤ ν_e sin i ≤ 60 km s^-1. About ten per cent have larger ve sin i (≥100 km s^-1), but surprisingly these show little or no nitrogen enhancement. All the cooler supergiants have low projected rotational velocities of ≤70 km s^-1 and high nitrogen abundance estimates, implying that either bi-stability braking or evolution on a blue loop may be important. Additionally, there is a lack of cooler binaries, possibly reflecting the small sample sizes. Single-star evolutionary models, which include rotation, can account for all of the nitrogen enhancement in both the single and binary samples. The detailed distribution of nitrogen abundances in the single and binary samples may be different, possibly reflecting differences in their evolutionary history. 

Conclusions: The first comparative study of single and binary B-type supergiants has revealed that the main sequence may be significantly wider than previously assumed, extending to T_eff = 20 000 K. Some marginal differences in single and binary atmospheric parameters and abundances have been identified, possibly implying non-standard evolution for some of the sample. This sample as a whole has implications for several aspects of our understanding of the evolutionary status of blue supergiants.

Broad balmer wings in BA Hyper/Supergiants distorted by diffuse interstellar bands: five examples in the 30 Doradus Region from the VLT-Flames Tarantula Survey

Extremely broad emission wings at Hβ and Hα have been found in VLT-FLAMES Tarantula Survey data for five very luminous BA supergiants in or near 30 Doradus in the Large Magellanic Cloud. The profiles of both lines are extremely asymmetrical, which we have found to be caused by very broad diffuse interstellar bands (DIBs) in the longward wing of Hβ and the shortward wing of Hα. These DIBs are well known to interstellar but not to many stellar specialists, so that the asymmetries may be mistaken for intrinsic features. The broad emission wings are generally ascribed to electron scattering, although we note difficulties for that interpretation in some objects. Such profiles are known in some Galactic hyper/supergiants and are also seen in both active and quiescent Luminous Blue Variables (LBVs). No prior or current LBV activity is known in these 30 Dor stars, although a generic relationship to LBVs is not excluded; subject to further observational and theoretical investigation, it is possible that these very luminous supergiants are approaching the LBV stage for the first time. Their locations in the HRD and presumed evolutionary tracks are consistent with that possibility. The available evidence for spectroscopic variations of these objects is reviewed, while recent photometric monitoring does not reveal variability. A search for circumstellar nebulae has been conducted, with an indeterminate result for one of them.

Rotational velocities of single and binary O-type stars in the Tarantula Nebula

Rotation is a key parameter in the evolution of massive stars, affecting their evolution, chemical yields, ionizing photon budget, and final fate. We determined the projected rotational velocity, υ e sin i, of ~330 O-type objects, i.e. ~210 spectroscopic single stars and ~110 primaries in binary systems, in the Tarantula nebula or 30 Doradus (30 Dor) region. The observations were taken using VLT/FLAMES and constitute the largest homogeneous dataset of multi-epoch spectroscopy of O-type stars currently available. The most distinctive feature of the υ e sin i distributions of the presumed-single stars and primaries in 30 Dor is a low-velocity peak at around 100 km s^-1. Stellar winds are not expected to have spun-down the bulk of the stars significantly since their arrival on the main sequence and therefore the peak in the single star sample is likely to represent the outcome of the formation process. Whereas the spin distribution of presumed-single stars shows a well developed tail of stars rotating more rapidly than 300 km s^-1, the sample of primaries does not feature such a high-velocity tail. The tail of the presumed-single star distribution is attributed for the most part - and could potentially be completely due - to spun-up binary products that appear as single stars or that have merged. This would be consistent with the lack of such post-interaction products in the binary sample, that is expected to be dominated by pre-interaction systems. The peak in this distribution is broader and is shifted toward somewhat higher spin rates compared to the distribution of presumed-single stars. Systems displaying large radial velocity variations, typical for short period systems, appear mostly responsible for these differences.

The VLT-FLAMES Tarantula Survey. XV. VFTS 822: A candidate Herbig B[e] star at low metallicity

We report the discovery of the B[e] star VFTS 822 in the 30 Doradus star-forming region of the Large Magellanic Cloud, classified by optical spectroscopy from the VLT-FLAMES Tarantula Survey and complementary infrared photometry. VFTS 822 is a relatively low-luminosity (log L = 4.04 ± 0.25 L_·) B8[e] star. In this Letter, we evaluate the evolutionary status of VFTS 822 and discuss its candidacy as a Herbig B[e] star. If the object is indeed in the pre-main sequence phase, it would present an exciting opportunity to spectroscopically measure mass accretion rates at low metallicity, to probe the effect of metallicity on accretion rates. 

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, 3e sin i, 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(3e), of the equatorial rotational velocity, 3e. Results. The distribution of 3e sin i shows a two-component structure: a peak around 80 km s1 and a high-velocity tail extending up to 600 km s^-1 This structure is also present in the inferred distribution P(3e) with around 80% of the sample having 0 <3e ≤ 300 km s^-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 highvelocity 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. © 2013 Author(s).

The VLT-FLAMES tarantula survey:X. Evidence for a bimodal distribution of rotational velocities for the single early B-type stars

Aims. Projected rotational velocities (ve sin i) have been estimated for 334 targets in the VLT-FLAMES Tarantula Survey that do not manifest significant radial velocity variations and are not supergiants. They have spectral types from approximately O9.5 to B3. The estimates have been analysed to infer the underlying rotational velocity distribution, which is critical for understanding the evolution of massive stars. Methods. Projected rotational velocities were deduced from the Fourier transforms of spectral lines, with upper limits also being obtained from profile fitting. For the narrower lined stars, metal and non-diffuse helium lines were adopted, and for the broader lined stars, both non-diffuse and diffuse helium lines; the estimates obtained using the different sets of lines are in good agreement. The uncertainty in the mean estimates is typically 4% for most targets. The iterative deconvolution procedure of Lucy has been used to deduce the probability density distribution of the rotational velocities. Results. Projected rotational velocities range up to approximately 450 kms-1 and show a bi-modal structure. This is also present in the inferred rotational velocity distribution with 25% of the sample having 0 <ve <100 km s^-1 and the high velocity component having v_e ∼ 250 km s^-1. There is no evidence from the spatial and radial velocity distributions of the two components that they represent either field and cluster populations or different episodes of star formation. Be-type stars have also been identified. Conclusions. The bi-modal rotational velocity distribution in our sample resembles that found for late-B and early-A type stars.While magnetic braking appears to be a possible mechanism for producing the low-velocity component, we can not rule out alternative explanations. © ESO 2013.

Transiting Exoplanet Studies and Community Targets for JWST's Early Release Science Program

The James Webb Space Telescope (JWST) will likely revolutionize transiting exoplanet atmospheric science, due to a combination of its capability for continuous, long duration observations and its larger collecting area, spectral coverage, and spectral resolution compared to existing space-based facilities. However, it is unclear precisely how well JWST will perform and which of its myriad instruments and observing modes will be best suited for transiting exoplanet studies. In this article, we describe a prefatory JWST Early Release Science (ERS) Cycle 1 program that focuses on testing specific observing modes to quickly give the community the data and experience it needs to plan more efficient and successful transiting exoplanet characterization programs in later cycles. We propose a multi-pronged approach wherein one aspect of the program focuses on observing transits of a single target with all of the recommended observing modes to identify and understand potential systematics, compare transmission spectra at overlapping and neighboring wavelength regions, confirm throughputs, and determine overall performances. In our search for transiting exoplanets that are well suited to achieving these goals, we identify 12 objects (dubbed “community targets”) that meet our defined criteria. Currently, the most favorable target is WASP-62b because of its large predicted signal size, relatively bright host star, and location in JWST's continuous viewing zone. Since most of the community targets do not have well-characterized atmospheres, we recommend initiating preparatory observing programs to determine the presence of obscuring clouds/hazes within their atmospheres. Measurable spectroscopic features are needed to establish the optimal resolution and wavelength regions for exoplanet characterization. Other initiatives from our proposed ERS program include testing the instrument brightness limits and performing phase-curve observations. The latter are a unique challenge compared to transit observations because of their significantly longer durations. Using only a single mode, we propose to observe a full-orbit phase curve of one of the previously characterized, short-orbital-period planets to evaluate the facility-level aspects of long, uninterrupted time-series observations.

The large- and small-scale Ca IIK structure of the Milky Way from observations of Galactic and Magellanic sightlines

Aims. The large and small-scale (pc) structure of the Galactic interstellar medium can be investigated by utilising spectra of early-type stellar probes of known distances in the same region of the sky. This paper determines the variation in line strength of Ca ii at 3933.661 Å as a function of probe separation for a large sample of stars, including a number of sightlines in the Magellanic Clouds. 

Methods. FLAMES-GIRAFFE data taken with the Very Large Telescope towards early-type stars in 3 Galactic and 4 Magellanic open clusters in Ca ii are used to obtain the velocity, equivalent width, column density, and line width of interstellar Galactic calcium for a total of 657 stars, of which 443 are Magellanic Cloud sightlines. In each cluster there are between 43 and 111 stars observed. Additionally, FEROS and UVES Ca ii K and Na i D spectra of 21 Galactic and 154 Magellanic early-type stars are presented and combined with data from the literature to study the calcium column density - parallax relationship. 

Results. For the four Magellanic clusters studied with FLAMES, the strength of the Galactic interstellar Ca ii K equivalent width on transverse scales from ∼0.05-9 pc is found to vary by factors of ∼1.8-3.0, corresponding to column density variations of ∼0.3-0.5 dex in the optically-thin approximation. Using FLAMES, FEROS, and UVES archive spectra, the minimum and maximum reduced equivalent widths for Milky Way gas are found to lie in the range ∼35-125 mÅ and ∼30-160 mÅ for Ca ii K and Na i D, respectively. The range is consistent with a previously published simple model of the interstellar medium consisting of spherical cloudlets of filling factor ∼0.3, although other geometries are not ruled out. Finally, the derived functional form for parallax (π) and Ca ii column density (N_CaII) is found to be π(mas) = 1 / (2.39 × 10^-13 × N_CaII (cm^-2) + 0.11). Our derived parallax is ∼25 per cent lower than predicted by Megier et al. (2009, A&A, 507, 833) at a distance of ∼100 pc and ∼15 percent lower at a distance of ∼200 pc, reflecting inhomogeneity in the Ca ii distribution in the different sightlines studied.