Source author record

H. Sana

H. Sana appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.SR astro-ph.GA astro-ph.IM astro-ph.CO astro-ph.HE astro-ph

Catalog footprint

What is connected

106works

6topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

A spectroscopic multiplicity survey of Galactic Wolf-Rayet stars: II. The northern WNE sequence

Most massive stars reside in multiple systems that will interact over the course of their lifetime. Classical Wolf-Rayet (WR) stars represent the final end stages of stellar evolution at the upper-mass end. As part of a homogeneous, magnitude-limited ($V\leq12$) spectroscopic survey of northern Galactic WR stars, this paper aims to establish the observed and intrinsic multiplicity properties of the early-type nitrogen-rich WR population (WNE). We obtained high-resolution spectroscopic time series of the complete magnitude-limited sample of 16 WNE stars observable with the 1.2 m Mercator telescope at La Palma, typically providing a time base of about two to eight years. We measured relative radial velocities (RVs) using cross-correlation and used RV variations to flag binary candidates. Adopting a peak-to-peak RV variability threshold of 50 km/s as a criterion found an observed multiplicity fraction of 0.44$\pm$0.12. Using an updated Monte Carlo method with a Bayesian framework, we calculated the three-dimensional likelihood for the intrinsic binary fraction, the maximum period, and the power-law index for the period distribution for the WNE population. We also re-derived multiplicity parameters for the Galactic WC population. We found an intrinsic multiplicity fraction of $0.56\substack{+0.20 \\ -0.15}$ for the parent WNE population. For the Galactic WC population, we re-derive an intrinsic multiplicity fraction of $0.96\substack{+0.04 \\ -0.22}$. The derived multiplicity parameters for the WNE population are quite similar to those derived for main-sequence O binaries but differ from those of the WC population. The significant shift in the WC period distribution towards longer periods is too large to be explained via expansion of the orbit due to stellar winds, and we discuss possible implications of our results.

preprint2022arXiv

Detecting Stripped Stars While Searching for Quiescent Black Holes

While the number of stellar-mass black holes detected in X-rays or as gravitational wave sources is steadily increasing, the known population remains orders of magnitude smaller than predicted by stellar evolution theory. A significant fraction of stellar-mass black holes is expected to hide in X-ray-quiet binaries where they are paired with a "normal" star. Although a handful of such quiescent black hole candidates have been proposed, the majority have been challenged by follow-up investigations. A confusion that emerged recently concerns binary systems that appear to contain a normal B-type star with an unseen companion, believed to be a black hole. On closer inspection, some of these seemingly normal B-type stars instead turn out to be stars stripped of most of their mass through an interaction with their binary companion, which in at least two cases is a rapidly rotating star rather than a compact object. These contaminants in the search for quiescent black holes are themselves extremely interesting objects as they represent a rare phase of binary evolution, and should be given special attention when searching for binaries hosting black holes in large spectroscopic studies.

preprint2022arXiv

Identifying quiescent compact objects in massive Galactic single-lined spectroscopic binaries

Aims. To search for these rare objects, we study 32 Galactic O-type stars that were reported as SB1s in the literature. In our sample we include Cyg X-1, which is known to host an accreting stellar-mass BH, and HD 74194, a supergiant fast X-ray transient, in order to validate our methodology. The final goal is to characterise the nature of the unseen companions to determine if they are MS stars, stripped helium stars, triples, or compact objects such as neutron stars or stellar-mass BHs. Methods. After measuring radial velocities and deriving orbital solutions for all the systems in our sample, we performed spectral disentangling to extract putative signatures of faint secondary companions from the composite spectra. We derived stellar parameters for the visible stars and estimated the mass ranges of the secondary stars using the binary mass function. Variability observed in the photometric TESS light curves was also searched for indications of the presence of putative companions, degenerate or not. Results. In 17 of the 32 systems reported as SB1s, we extract secondary signatures, down to mass ratios of ~0.15. For the 17 newly detected double-lined spectroscopic binaries (SB2s), we derive physical properties of the individual components and discuss why they have not been detected as such before. Among the remaining systems, we identify nine systems with possible NS or low-mass MS companions. For Cyg X-1 and HD 130298, we are not able to extract any signatures for the companions, and the minimum masses of their companions are estimated to be about 7Msun. Our simulations show that secondaries with such a mass should be detectable from our dataset, no matter their nature: MS stars, stripped helium stars or even triples. While this is expected for Cyg X-1, confirming our methodology, our simulations also strongly suggest that HD 130298 could be another candidate to host a stellar-mass BH.

preprint2022arXiv

Modeling overcontact binaries, I. The effect of tidal deformation

In the realm of massive stars, strong binary interaction is commonplace. One extreme case are overcontact systems, which is expected to be part of the evolution of all stars evolving towards a merger, and is hypothesized to play a role in the formation of binary black holes. However, important simplifications are made to model the evolution of overcontact binaries, namely that tidal or rotational deformation is frequently ignored. Yet, both observation and theory show that overcontact stars are heavily tidally deformed, leaving a potentially important effect on the outer layers unaccounted for in models. In this work we develop the methodology to represented tidally deformed stars. Using numerical methods, we compute the structure correction factors to the 1D spherical stellar structure equations due to the binary Roche potential, and compare them to existing results and the structure corrections of single rotating stars. We implement the new structure correction factors into the stellar evolution code MESA and explore several case studies. We compare the differences between our simulations when no rotation is included, when we treat rotation using single star corrections (i.e. only accounting for centrifugal deformation) or when we use tidal deformation. We find that ignoring rotation in deformed detached eclipsing binaries can produce a radius discrepancy of up to 5%. The difference between tidal and single star centrifugal distortion models is more benign at 1%, showing that single rotating star models are a suitable approximation of tidally deformed stars in a binary system. In overcontact configurations, we find a similar 5% variation in surface properties as a result of tidal distortion with respect to non-rotating models, showing that it is inappropriate to model binary stars that fill their Roche lobe significantly, as non-rotating.

preprint2022arXiv

On the origin of close massive binaries in the M17 star-forming region

Spectroscopic multiplicity surveys of O stars in young clusters and OB associations have revealed that a large portion ($\sim$ 70%) of these massive stars (M$_{i}$ $\gt$ 15 $M_{\odot}$) belong to close and short-period binaries (physical separation d $\lt$few au). Despite the recent and significant progress, the formation mechanisms leading to such close massive multiple systems remain to be elucidated. As a result, young massive close binaries (or higher-order multiple systems) are unique laboratories to figure out the pairing mechanism of high-mass stars. We present the first VLTI/GRAVITY observations of six young O stars in the M17 star-forming region ($\lesssim$ 1 Myr) and two additional foreground stars. From the interferometric model fitting of visibility amplitudes and closure phases, we search for companions and measure their positions and flux ratios. Combining the resulting magnitude difference with atmosphere models and evolutionary tracks, we further constrain the masses of the individual components. All of the six high-mass stars are in multiple systems, leading to a multiplicity fraction (MF) of 100%, yielding a 68% confidence interval of 94-100%. We detect a total number of 9 companions with separations up to 120 au. Including previously identified spectroscopic companions, the companion fraction of the young O stars in our sample reaches 2.3$\pm$0.6. The derived masses span a wide range from 2.5 to 50 $M_{\odot}$, with a great tendency towards high-mass companions. While based on a modest sample, our results clearly indicate that the origin of the high degree of multiplicity is rooted in their star formation mechanism. No clear evidence for one of the competing concepts of massive star formation (core accretion or competitive accretion) could be found. However, our results are compatible with migration as a scenario for the formation of close massive binaries.

preprint2022arXiv

Properties of the Be-type stars in 30 Doradus

The evolutionary status of Be-type stars remains unclear, with both single-star and binary pathways having been proposed. Here, VFTS spectroscopy of 73 Be-type stars, in the spectral-type range, B0--B3, is analysed to estimate projected rotational velocities, radial velocities and stellar parameters. They are found to be rotating faster than the corresponding VFTS B-type sample but simulations imply that their projected rotational velocities are inconsistent with them all rotating at near critical velocities. The de-convolution of the projected rotational velocities estimates leads to a mean rotational velocity estimate of 320-350 km/s, approximately 100 km/s larger than that for the corresponding B-type sample. There is a dearth of targets with rotational velocities less than 0.4 of the critical velocity, with a broad distribution reaching up to critical rotation. Our best estimate for the mean or median of the rotational velocitiy is 0.68 of the critical velocity. Rapidly-rotating B-type stars are more numerous than their Be-type counterparts, whilst the observed frequency of Be-type stars identified as binary systems is significantly lower than that for normal B-type stars, consistent with their respective radial-velocity dispersions. The semi-amplitudes for the Be-type binaries are also smaller. Similar results are found for a SMC Be-type sample centred on NGC346 with no significant differences being found between the two samples. These results are compared with the predictions of single and binary stellar evolutionary models for Be-type stars. Assuming that a single mechanism dominated the production of classical Be-type stars, our comparison would favour a binary evolutionary history.

preprint2021arXiv

A relation between the radial velocity dispersion of young clusters and their age: Evidence for hardening as the formation scenario of massive close binaries

The majority of massive stars ($>8$ $\rm{M_{\odot}}$) in OB associations are found in close binary systems. Nonetheless, the formation mechanism of these close massive binaries is not understood yet. Using literature data, we measured the radial-velocity dispersion ($σ_\mathrm{RV}$) as a proxy for the close binary fraction in ten OB associations in the Galaxy and the Large Magellanic Cloud, spanning an age range from 1 to 6 Myrs. We find a positive trend of this dispersion with the cluster's age, which is consistent with binary hardening. Assuming a universal binary fraction of $f_\mathrm{bin}$ = 0.7, we converted the $σ_\mathrm{RV}$ behavior to an evolution of the minimum orbital period $P_\mathrm{cutoff}$ from $\sim$9.5 years at 1 Myr to $\sim$1.4 days for the oldest clusters in our sample at $\sim$6 Myr. Our results suggest that binaries are formed at larger separations, and they harden in around 1 to 2 Myrs to produce the period distribution observed in few million year-old OB binaries. Such an inward migration may either be driven by an interaction with a remnant accretion disk or with other young stellar objects present in the system. Our findings constitute the first empirical evidence in favor of migration as a scenario for the formation of massive close binaries.

preprint2021arXiv

Analytic, dust-independent mass-loss rates for red supergiant winds initiated by turbulent pressure

Context. Red supergiants are observed to undergo vigorous mass-loss. However, to date, no theoretical model has succeeded in explaining the origins of these objects' winds. This strongly limits our understanding of red supergiant evolution and Type II-P and II-L supernova progenitor properties. Aims. We examine the role that vigorous atmospheric turbulence may play in initiating and determining the mass-loss rates of red supergiant stars. Methods. We analytically and numerically solve the equations of conservation of mass and momentum, which we later couple to an atmospheric temperature structure, to obtain theoretically motivated mass-loss rates. We then compare these to state-of-the-art empirical mass-loss rate scaling formulae as well as observationally inferred mass-loss rates of red supergiants. Results. We find that the pressure due to the characteristic turbulent velocities inferred for red supergiants is sufficient to explain the mass-loss rates of these objects in the absence of the normally employed opacity from circumstellar dust. Motivated by this initial success, we provide a first theoretical and fully analytic mass-loss rate prescription for red supergiants. We conclude by highlighting some intriguing possible implications of these rates for future studies of stellar evolution, especially in light of the lack of a direct dependence on metallicity.

preprint2021arXiv

Binary-object spectral-synthesis in 3D (BOSS-3D) -- Modelling H-alpha emission in the enigmatic multiple system LB-1

Context: To quantitatively decode the information stored within an observed spectrum, detailed modelling of the physical state and accurate radiative transfer solution schemes are required. In the analysis of stellar spectra, the numerical model often needs to account for binary companions and 3D structures in the stellar envelopes. The enigmatic binary (or multiple) system LB-1 constitutes a perfect example of such a complex multi-D problem. Aims: To improve our understanding of the LB-1 system, we directly modelled the phase-dependent H-alpha line profiles of this system. To this end, we developed a multi-purpose binary-object spectral-synthesis code in 3D (BOSS-3D). Methods: BOSS-3D calculates synthetic line profiles for a given state of the circumstellar material. The standard pz-geometry commonly used for single stars is extended by defining individual coordinate systems for each involved object and by accounting for the appropriate coordinate transformations. The code is then applied to the LB-1 system, considering two main hypotheses, a binary containing a stripped star and Be star, or a B star and a black hole with a disc. Results: Comparing these two scenarios, neither model can reproduce the detailed phase-dependent shape of the H-alpha line profiles. A satisfactory match with the observations, however, is obtained by invoking a disc around the primary object in addition to the Be-star disc or the black-hole accretion disc. Conclusions: The developed code can be used to model synthetic line profiles for a wide variety of binary systems, ranging from transit spectra of planetary atmospheres, to post-asymptotic giant branch binaries including circumstellar and circumbinary discs and massive-star binaries with stellar winds and disc systems. For the LB-1 system, our modelling provides strong evidence that each object in the system contains a disc-like structure.

preprint2021arXiv

Detailed models of interacting short-period massive binary stars

About a quarter of massive binary stars undergo mass transfer while both stars burn hydrogen at their cores, first on the thermal and then on the nuclear timescale. The nuclear timescale mass transfer leads to observable counterparts: the semi-detached so-called massive Algol binaries. However, comprehensive model predictions for these systems are sparse. We study them using a large grid of ~10,000 detailed binary evolution models calculated with the stellar evolution code MESA, covering initial donor masses between 10-40 M$_{\odot}$ and initial orbital periods above 1.4 d, at a metallicity suitable for the Large Magellanic Cloud (LMC). Our models imply ~30, or ~3% of the ~1,000 core hydrogen burning O-star binaries in the LMC to be currently in the semi-detached phase. Our donor models are up to 25-times more luminous than single stars of identical mass and effective temperature, which agrees with the observed Algols. A comparison of our models with the observed orbital periods and mass ratios implies rather conservative mass transfer in some systems, while very inefficient one in others. This is generally well reproduced by our spin-dependent mass transfer algorithm, except for the lowest considered masses. The observations reflect the slow increase of the surface nitrogen enrichment of the donors during the semi-detached phase all the way to CNO equilibrium. We also investigate the properties of our models after core hydrogen depletion of the donor star, when these models correspond to Wolf-Rayet/helium+OB star binaries. A dedicated spectroscopic survey of massive Algol systems may allow to derive the dependence of the efficiency of thermal timescale mass transfer on the binary parameters, as well as the efficiency of semiconvective mixing in the stellar interior. This would be a crucial step towards reliable binary models up to the formation of supernovae and compact objects.

preprint2021arXiv

Probing the low-mass end of the companion mass function for O-type stars

Almost all massive stars are part of a binary system. Given the wide range of separations at which these companions are found, several observational techniques have been adopted to characterize them, but contrasts greater than 4 in the H-band have never been reached between 0".1 and 1". We used VLT/SPHERE to observe simultaneously with the IRDIS and IFS sub-systems 18 O-type stars within 6 kpc and ages between 1-5 Myrs to probe the existence of stellar companions in the angular separation range from 0".15 to 6" down to very low mass ratios. The IFS YJH- band observations have allowed us to probe the presence of sub-solar companions in a 1".7x1".7 field-of-view down to magnitude limits of deltaH=10 at 0".4. In the wider 12"x12" IRDIS field-of-view, we reached contrasts of deltaK=12 at 1", enabling us to look for even fainter companions. This paper presents five newly discovered intermediate (<1") separation companions, three of which are smaller than 0.2M_sun. If confirmed by future analyses of proper motions, these new detections represent the lowest-mass companions ever found around O-type stars. Assuming that all sources detected within 1" are physically bound, the observed fraction of companions for O-type stars between 0".15 and 0".9 is 0.39+/-0.15, whereas it increases to 1.6+/-0.3 in the separation range from 0".9 to 6". These findings clearly support the notion that massive stars form almost exclusively in multiple systems, and that larger AO-assisted coronagraphic surveys are crucial in placing constraints on the multiplicity properties of massive star companions in regions of the parameter space that have previously gone unexplored, and demonstrate that the companion mass function is populated down to the lowest stellar masses.

preprint2021arXiv

The CubeSpec space mission. I. Asteroseismology of massive stars from time series optical spectroscopy: science requirements and target list prioritisation

Context. There is currently a niche for providing high-cadence, high resolution, time-series optical spectroscopy from space, which can be filled by using a low-cost cubesat mission. The Belgian-led ESA CubeSpec mission is specifically designed to provide space-based, low-cost spectroscopy with specific capabilities that can be optimised for a particular science need. Approved as an ESA in-orbit demonstrator, the CubeSpec satellite's primary science objective will focus on obtaining high-cadence, high resolution optical spectroscopic data to facilitate asteroseismology of pulsating massive stars. Aims. In this first paper, we aim to search for pulsating massive stars suitable for the CubeSpec mission, specifically $β$ Cep stars, which typically require time series spectroscopy to identify the geometry of their pulsation modes. Methods. Based on the science requirements needed to enable asteroseismology of massive stars with the capabilities of CubeSpec's spectrograph, we combine a literature study for pulsation with the analysis of recent high-cadence time series TESS photometry to classify the variability for stars brighter than V < 4 mag and between O9 and B3 in spectral type. Results. Among the 90 stars that meet our magnitude and spectral type requirements, we identify 23 promising $β$ Cep stars with high-amplitude (non-)radial pulsation modes with frequencies below 7 d$^{-1}$. Using further constraints on projected rotational velocities, pulsation amplitudes and number of pulsation modes, we devise a prioritised target list for the CubeSpec mission according to its science requirements and the potential of the targets for asteroseismology. The full target catalogue further provide a modern TESS-based review of line profile and photometric variability properties among bright O9-B3 stars.

preprint2021arXiv

The inner circumstellar dust of the red supergiant Antares as seen with VLT/SPHERE/ZIMPOL

The processes by which red supergiants lose mass are not fully understood thus-far and their mass-loss rates lack theoretical constraints. The ambient surroundings of the nearby M0.5 Iab star Antares offers an ideal environment to obtain detailed empirical information on the outflow properties at its onset, and hence indirectly, on the mode(s) of mass loss. We present and analyse optical VLT/SPHERE/ZIMPOL polarimetric imaging with angular resolution down to 23 milli-arcsec, sufficient to spatially resolve both the stellar disk and its direct surroundings. We detect a conspicuous feature in polarised intensity that we identify as a clump containing dust, which we characterise through 3D radiative transfer modelling. The clump is positioned behind the plane of the sky, therefore has been released from the backside of the star, and its inner edge is only 0.3 stellar radii above the surface. The current dust mass in the clump is $1.3^{+0.2}_{-1.0} \times 10^{-8}$ M$_{\odot}$, though its proximity to the star implies that dust nucleation is probably still ongoing. The ejection of clumps of gas and dust makes a non-negligible contribution to the total mass lost from the star which could possibly be linked to localised surface activity such as convective motions or non-radial pulsations.

preprint2021arXiv

The observed multiplicity properties of B-type stars in the Galactic young open cluster NGC 6231

It is well known that massive O-stars are frequently (if not always) found in binary or higher-order multiple systems, but this fact has been less robustly investigated for the lower mass range of the massive stars, represented by the B-type stars. Obtaining the binary fraction and orbital parameter distributions of B-type stars is crucial to understand the impact of multiplicity on the archetypal progenitor of core-collapse supernovae as well as to properly investigate formation channels for gravitational wave progenitors. This work aims to characterise the multiplicity of the B-star population of the young open cluster NGC 6231 through multi-epoch optical spectroscopy of 80 B-type stars. We analyse 31 FLAMES/GIRAFFE observations of 80 B-type stars, monitoring their radial velocities (RVs) and performing a least-squares spectral analysis (Lomb-Scargle) to search for periodicity in those stars with statistically significant variability in their RVs. We constrain an observed spectroscopic binary fraction of $33\pm5$% for the B-type stars of NGC 6231, with a first order bias-correction giving a true spectroscopic binary fraction of $52\pm8$%. Out of 27 B-type binary candidates, we obtained orbital solutions for 20 systems: 15 single-lined (SB1) and 5 double-lined spectroscopic binaries (SB2s). We present these orbital solutions and the orbital parameter distributions associated with them. Our results indicate that Galactic B-type stars are less frequently found in binary systems than their more massive O-type counterparts, but their orbital properties generally resemble those of B- and O-type stars in both the Galaxy and Large Magellanic Cloud.

preprint2021arXiv

Uncovering astrometric black hole binaries with massive main-sequence companions with Gaia

The hunt for compact objects is on. Rarely seen massive binaries with a compact object are a crucial phase in the evolution towards compact object mergers. In Gaia data release 3 (DR3), the first Gaia astrometric orbital solutions for binary sources will become available. We investigate how many black holes (BH) with massive main-sequence dwarf companions (OB+BH binaries) are expected to be detected as binaries in DR3 and at the end of the nominal 5-yr mission (DR4). We estimate the fraction of identifiable OB+BH binaries and discuss the distributions of the masses of both components and the orbital periods. We study the impact of different BH-formation scenarios. Using tailored models for the massive star population, which assume a direct collapse and no kick upon BH formation (the fiducial case), we estimate the fraction of OB+BH systems that Gaia will detect as binaries. A distance distribution according to that of the second Alma Luminous Star catalogue (ALSII) is assumed. We investigate how many of the systems detected as binaries are identifiable as OB+BH binaries, using a method based on astrometric data. In the fiducial case we conservatively estimate that 77% of the OB+BH binaries in ALSII will be detected as binaries in DR3, of which 89% are identifiable as OB+BH binaries. This leads to a total of around 190 OB+BH binaries, a 20-fold increase in the known sample of OB+BH binaries, covering an uncharted parameter space of long-period binaries. The size and properties of the identifiable OB+BH population will contain crucial observational constraints to improve our understanding of BH formation. In DR4, the detected fraction will increase to 85%, of which 82% will be identifiable. Hence, an additional ~5 systems could be identified, which are expected to have either very short or long periods. The fractions become smaller for different BH-formation scenarios. (truncated)

preprint2020arXiv

Investigating the lack of main-sequence companions to massive Be stars

About 20% of all B-type stars are classical Be stars. The Be phenomenon is strongly correlated with rapid rotation, the origin of which remains unclear. It may be rooted in single- or binary-star evolution. In the framework of the binary channel, the initially more massive star transfers mass and angular momentum to the original secondary, which becomes a Be star. The system then evolves into a Be binary with a post-main-sequence companion, which may later be disrupted in a supernova event. Hence, if the binary channel dominates the formation of Be stars, one may expect a strong lack of close Be binaries with main sequence (MS) companions. Through an extensive, star-by-star review of the literature of a magnitude-limited sample of Galactic early-type Be stars, we investigate whether Be binaries with MS companions are known to exist. Our sample is constructed from the BeSS database and cross-matched with all available literature on the individual stars. Out of an initial list of 505 Be stars, we compile a final sample of 287 Galactic Be stars earlier than B1.5 with V<=12 mag. Out of those, 13 objects were reported as Be binaries with known post-MS companions and 11 as binaries with unknown, uncertain or debated companions. We find no confirmed reports of Be binaries with MS companions. For the remaining 263 targets, no significant reports of multiplicity exist in the literature, implying that they are either Be binaries with faint companions, or truly single. The clear lack of reported MS companions to Be stars, which stands in contrast to the high number of detected B+B MS binaries, strongly supports the hypothesis that early-type Be stars are binary interaction products that spun up after mass and angular momentum transfer from a companion star. Taken at face value, our results may suggest that a large majority of the early-type Be stars have formed through binary mass-transfer.

preprint2020arXiv

Is HR 6819 a triple system containing a black hole? -- An alternative explanation

HR 6819 was recently proposed to be a triple system consisting of an inner B-type giant + black hole binary with an orbital period of 40d and an outer Be tertiary. This interpretation is mainly based on two inferences: that the emission attributed to the outer Be star is stationary, and that the inner star, which is used as mass calibrator for the black hole, is a B-type giant. We re-investigate the properties of HR 6819 by spectral disentangling and an atmosphere analysis of the disentangled spectra to search for a possibly simpler alternative explanation for HR 6819. Disentangling implies that the Be component is not a static tertiary, but rather a component of the binary in the 40-d orbit. The inferred radial velocity amplitudes imply an extreme mass ratio of M_2/M_1 = 15 +/- 3. We infer spectroscopic masses of 0.4$^{+0.3}_{-0.1}$ Msun and 6$^{+5}_{-3}$ Msun for the primary and secondary, which agree well with the dynamical masses for an inclination of i = 32 deg. This indicates that the primary might be a stripped star rather than a B-type giant. Evolutionary modelling suggests that a possible progenitor system would be a tight (P_i ~ 2d) B+B binary system that experienced conservative mass transfer. While the observed nitrogen enrichment of the primary conforms with the predictions of the evolutionary models, we find no indications for the predicted He enrichment. We suggest that HR 6819 is a binary system consisting of a stripped B-type primary and a rapidly-rotating Be star that formed from a previous mass-transfer event. In the framework of this interpretation, HR 6819 does not contain a black hole. Interferometry can distinguish between these two scenarios by providing an independent measurement of the separation between the visible components.

preprint2020arXiv

Multiplicity of the red supergiant population in the young massive cluster NGC 330

The multiplicity properties of massive stars are one of the important outstanding issues in stellar evolution. Quantifying the binary statistics of all evolutionary phases is essential to paint a complete picture of how and when massive stars interact with their companions, and to determine the consequences of these interactions. We investigate the multiplicity of an almost complete census of red supergiant stars (RSGs) in NGC 330, a young massive cluster in the SMC. Using a combination of multi-epoch HARPS and MUSE spectroscopy, we estimate radial velocities and assess the kinematic and multiplicity properties of 15 RSGs in NGC 330. Radial velocities are estimated to better than +/-100 m/s for the HARPS data. The line-of-sight velocity dispersion for the cluster is estimated as 3.20 +0.69-0.52 km/s. When virial equilibrium is assumed, the dynamical mass of the cluster is log (M{dyn} /M{sun}) = 5.20+/-0.17, in good agreement with previous upper limits. We detect significant radial velocity variability in our multi-epoch observations and distinguish between variations caused by atmospheric activity and those caused by binarity. The binary fraction of NGC 330 RSGs is estimated by comparisons with simulated observations of systems with a range of input binary fractions. In this way, we account for observational biases and estimate the intrinsic binary fraction for RSGs in NGC 330 as f{RSG} = 0.3+/-0.1 for orbital periods in the range 2.3< log P [days] <4.3, with q>0.1. Using the distribution of the luminosities of the RSG population, we estimate the age of NGC 330 to be 45+/-5 Myr and estimate a red straggler fraction of 50%. We estimate the binary fraction of RSGs in NGC 330 and conclude that it appears to be lower than that of main-sequence massive stars, which is expected because interactions between an RSG and a companion are assumed to effectively strip the RSG envelope.

preprint2020arXiv

Properties of OB star-black hole systems derived from detailed binary evolution models

The recent gravitational wave measurements have demonstrated the existence of stellar mass black hole binaries. It is essential for our understanding of massive star evolution to identify the contribution of binary evolution to the formation of double black holes. A promising way to progress is investigating the progenitors of double black hole systems and comparing predictions with local massive star samples such as the population in 30 Doradus in the Large Magellanic Cloud (LMC). Methods. To this purpose, we analyse a large grid of detailed binary evolution models at LMC metallicity with initial primary masses between 10 and 40 Msun, and identify which model systems potentially evolve into a binary consisting of a black hole and a massive main sequence star. We then derive the observable properties of such systems, as well as peculiarities of the OB star component. We find that about 3% of the LMC late O and early B stars in binaries are expected to possess a black hole companion, when assuming stars with a final helium core mass above 6.6 M to form black holes. While the vast majority of them may be X-ray quiet, our models suggest that these may be identified in spectroscopic binaries, either by large amplitude radial velocity variations ( > 50 km s ) and simultaneous nitrogen surface enrichment, or through a moderate radial velocity ( > 10 km/s ) and simultaneously rapid rotation of the OB star. The predicted mass ratios are such that main sequence companions could be excluded in most cases. A comparison to the observed OB+WR binaries in the LMC, Be/X-ray binaries, and known massive BH binaries supports our conclusion. We expect spectroscopic observations to be able to test key assumptions in our models, with important implications for massive star evolution in general, and for the formation of double-black hole mergers in particular.

preprint2020arXiv

The 'hidden' companion in LB-1 unveiled by spectral disentangling

The intriguing binary LS V +22 25 (LB-1) has drawn much attention following claims of it being a single-lined spectroscopic binary with a 79-day orbit comprising a B-type star and a ~70 Msun black hole. Recent analyses have implied that the visible primary star is a stripped He-rich star. However, the nature of the secondary, which was proposed to be a black hole, a neutron star, or a main sequence star, remains unknown. Based on 26 newly acquired spectroscopic observations, we perform an orbital analysis and spectral disentangling of LB-1 to elucidate the nature of the system. Our analysis reveals that LB-1 contains two components of comparable brightness in the optical. The narrow-lined primary, which we estimate to contribute ~55% in the optical, has spectral properties that suggest that it is a stripped star: it has a small spectroscopic mass (~1 Msun) for a B-type star and it is He- and N-rich. The "hidden" secondary, which contributes about 45% of the optical flux, is a rapidly rotating (vsini ~ 300 km/s) B3 V star with a decretion disk -- a Be star. Hence, LB-1 does not contain a compact object. Instead, it is a rare Be binary system consisting of a stripped star (the former mass donor) and a Be star rotating at near its critical velocity (the former mass accretor). This system is a clear example that binary interactions play a decisive role in the production of rapid stellar rotators and Be stars.

preprint2020arXiv

The high-energy emission from HD~93129A near periastron

We conducted an observational campaign towards one of the most massive and luminous colliding wind binaries in the Galaxy, HD~93129A, close to its periastron passage in 2018. During this time the source was predicted to be in its maximum of high-energy emission. Here we present our data analysis from the X-ray satellites \textit{Chandra} and \textit{NuSTAR} and the $γ$-ray satellite \textit{AGILE}. High-energy emission coincident with HD~93129A was detected in the X-ray band up to $\sim$18~keV, whereas in the $γ$-ray band only upper limits were obtained. We interpret the derived fluxes using a non-thermal radiative model for the wind-collision region. We establish a conservative upper limit for the fraction of the wind kinetic power that is converted into relativistic electron acceleration, $f_\mathrm{NT,e} < 0.02$. In addition, we set a lower limit for the magnetic field in the wind-collision region as $B_\mathrm{WCR} > 0.3$~G. We also argue a putative interpretation of the emission from which we estimate $f_\mathrm{NT,e} \approx 0.006$ and $B_\mathrm{WCR} \approx 0.5$~G. We conclude that multi-wavelength, dedicated observing campaigns during carefully selected epochs are a powerful tool for characterising the relativistic particle content and magnetic field intensity in colliding wind binaries.

preprint2020arXiv

The VLT-FLAMES Tarantula Survey XXXII. Low-luminosity late O-type stars -- classification, main physical parameters, and silicon abundances

Analysis of late O-type stars observed in the Large Magellanic Cloud (LMC) by the VLT-FLAMES Tarantula Survey (VFTS) revealed a discrepancy between the physical properties estimated from model-atmosphere analysis and those expected from their morphological classifications. Here we revisit the analysis of 32 of these puzzling objects using new hydrogen-helium-silicon FASTWIND models and a different fitting approach to re-evaluate their physical properties. Our new analysis confirms that these stars indeed have properties that are typical of late O-type dwarfs. We also present the first estimates of silicon abundances for O-type stars in the 30 Dor clusters NGC 2060 and NGC 2070, with a weighted mean abundance for our sample of 7.05 +/- 0.03. Our values are about 0.20 dex lower than those previously derived for B-type stars in the LMC clusters N 11 and NGC 2004 using TLUSTY models. Various possibilities (e.g. differences in the analysis methods, effects of microturbulence, and real differences between stars in different clusters) were considered to account for these results. We also used our grid of FASTWIND models to reassess the impact of using the Galactic classification criteria for late O-type stars in the LMC by scrutinising their sensitivity to different stellar properties. At the cool edge of the O star regime the HeII 4686/HeI 4713 ratio used to assign luminosity class for Galactic stars can mimic giants or bright giants in the LMC, even for objects with high gravities (log_g > 4.0 dex). We argue that this line ratio is not a reliable luminosity diagnostic for late O-type stars in the LMC, and that the SiIV 4989/HeI4026 ratio is more robust for these types.

preprint2020arXiv

The Wolf-Rayet binaries of the nitrogen sequence in the Large Magellanic Cloud: spectroscopy, orbital analysis, formation, and evolution

Massive Wolf-Rayet (WR) stars dominate the radiative and mechanical energy budget of galaxies and probe a critical phase in the evolution of massive stars prior to core-collapse. It is not known whether core He-burning WR stars (classical WR, cWR) form predominantly through wind-stripping (w-WR) or binary stripping (b-WR). With spectroscopy of WR binaries so-far largely avoided due to its complexity, our study focuses on the 44 WR binaries / binary candidates of the Large Magellanic Cloud (LMC, metallicity Z~0.5 Zsun), identified on the basis of radial velocity variations, composite spectra, or high X-ray luminosities. Relying on a diverse spectroscopic database, we aim to derive the physical and orbital parameters of our targets, confronting evolution models of evolved massive stars at sub-solar metallicity, and constraining the impact of binary interaction in forming them. Spectroscopy is performed using the Potsdam Wolf-Rayet (PoWR) code and cross-correlation techniques. Disentanglement is performed using the code Spectangular or the shift-and-add algorithm. Evolutionary status is interpreted using the Binary Population and Spectral Synthesis (BPASS) code, exploring binary interaction and chemically-homogeneous evolution. No obvious dichotomy in the locations of apparently-single and binary WN stars on the Hertzsprung-Russell diagram is apparent. According to commonly used stellar evolution models (BPASS, Geneva), most apparently-single WN stars could not have formed as single stars, implying that they were stripped by an undetected companion. Otherwise, it must follow that pre-WR mass-loss/mixing (e.g., during the red supergiant phase) are strongly underestimated in standard stellar evolution models.

preprint2020arXiv

VLT/X-shooter spectroscopy of massive young stellar objects in the 30 Doradus region of the Large Magellanic Cloud

The process of massive star ($M\geq8~M_\odot$) formation is still poorly understood. Observations of massive young stellar objects (MYSOs) are challenging due to their rarity, short formation timescale, large distances, and high circumstellar extinction. Here, we present the results of a spectroscopic analysis of a population of MYSOs in the Large Magellanic Cloud (LMC). We took advantage of the spectral resolution and wavelength coverage of X-shooter (300-2500 nm), mounted on the European Southern Observatory Very Large Telescope, to detect characteristic spectral features in a dozen MYSO candidates near 30 Doradus, the largest starburst region in the Local Group hosting the most massive stars known. The X-shooter spectra are strongly contaminated by nebular emission. We used a scaling method to subtract the nebular contamination from our objects. We detect H$α,β$, [O I] 630.0 nm, Ca II infrared triplet, [Fe II] 1643.5 nm, fluorescent Fe II 1687.8 nm, H$_2$ 2121.8 nm, Br$γ$, and CO bandhead emission in the spectra of multiple candidates. This leads to the spectroscopic confirmation of 10 candidates as bona fide MYSOs. We compare our observations with photometric observations from the literature and find all MYSOs to have a strong near-infrared excess. We compute lower limits to the brightness and luminosity of the MYSO candidates, confirming the near-infrared excess and the massive nature of the objects. No clear correlation is seen between the Br$γ$ luminosity and metallicity. Combining our sample with other LMC samples results in a combined detection rate of disk features such as fluorescent Fe II and CO bandheads which is consistent with the Galactic rate (40\%). Most of our MYSOs show outflow features.

preprint2020arXiv

Why binary interaction does not necessarily dominate the formation of Wolf-Rayet stars at low metallicity

Classical Wolf-Rayet (WR) stars are massive, hydrogen depleted, post main-sequence stars that exhibit emission-line dominated spectra. For a given metallicity Z, stars exceeding a certain initial mass M_single(Z) can reach the WR phase through intrinsic mass-loss (single-star channel). Stars of lower masses can reach the WR phase via binary mass transfer (binary channel). It is commonly assumed that the binary channel dominates the formation of WR stars in environments with low Z such as the SMC and LMC. However, their reported WR binary fractions of 30-40% are comparable to that of the Galaxy. Here, we explain this apparent contradiction by considering the minimum initial mass M_spec(Z) needed for the stripped product to appear as a WR star. We calibrate M_spec(Z) using the lowest-luminosity WR stars in the Clouds and the Galaxy. A range of M_single(Z) values are explored using various evolution codes. We estimate the additional contribution of the binary channel by considering the interval [M_spec(Z), M_single(Z)], which characterises the initial-mass range in which binaries can form additional WR stars. Results: The WR-phenomenon ceases below luminosities of logL = 4.9, 5.25, and 5.6 [Lsun] in the Galaxy, LMC, and SMC, which translates to He-star masses of 7.5, 11, 17 Msun and initial masses of M_spec = 18, 23, 37 Msun. Stripped stars with lower initial masses in the respective galaxies would tend to not appear as WR stars. M_single lies in the range 20-30, 30-60, and > 40 Msun for the Galaxy, LMC, and SMC. We find that that the additional contribution of the binary channel is a non-trivial function of Z that cannot be conclusively claimed to be monotonically increasing with decreasing Z. Hence, one should not a-priori expect that binary interactions become increasingly important in forming WR stars at low Z, or that the WR binary fraction grows with decreasing Z.

preprint2019arXiv

The Tarantula Massive Binary Monitoring: III. Atmosphere analysis of double-lined spectroscopic systems

Accurate stellar parameters of individual objects in binary systems are essential to constrain the effects of binarity on stellar evolution. These parameters serve as a prerequisite to probing existing and future theoretical evolutionary models. We aim to derive the atmospheric parameters of the 31 SB2s in the TMBM sample. This sample, composed of detached, semi-detached and contact systems with at least one of the components classified as an O star, is an excellent test-bed to study how binarity can impact our knowledge of the evolution of massive stars. 32 epochs of FLAMES/GIRAFFE spectra are analysed using spectral disentangling to construct the individual spectra of 62 components. We apply the CMFGEN atmosphere code to determine their stellar parameters and their He, C and N surface abundances. From these properties, we show that the effects of tides on chemical mixing are limited. Components on longer-period orbits show higher nitrogen enrichment at their surface than those on shorter-period orbits, in contrast to expectations of rotational or tidal mixing, implying that other mechanisms play a role in this process. Components filling their Roche lobe are mass donors. They exhibit higher nitrogen content at their surface and rotate more slowly than their companions. By accreting new material, their companions spin faster and are rejuvenated. Their locations in the N-vsini diagram tend to show that binary products are good candidates to populate the two groups of stars (slowly rotating, nitrogen-enriched and rapidly rotating non-enriched) that cannot be reproduced through single-star population synthesis. This sample is the largest sample of binaries to be studied in such a homogeneous way. The study of these objects gives us strong observational constraints to test theoretical binary evolutionary tracks.

preprint2019arXiv

The Tarantula Massive Binary Monitoring. IV. Double-lined photometric binaries

A high fraction of massive stars are found to be binaries but only a few of them are reported as photometrically variable. By studying the populations of SB2 in the 30 Doradus region, we found a subset of them that have photometry from the OGLE project and that display variations in their light curves related to orbital motions. The goal of this study is to determine the dynamical masses and radii of the 26 binary components to investigate the mass-discrepancy problem and to provide an empirical mass-luminosity relation for the LMC. We use the PHOEBE programme to perform a systematic analysis of the OGLE V and I light curves obtained for 13 binary systems in 30 Dor. We adopt Teff, and orbital parameters derived previously to obtain the inclinations of the systems and the parameters of the individual components. Three systems display eclipses in their light curves, while the others only display ellipsoidal variations. We classify two systems as over-contact, five as semi-detached, and four as detached. The two remaining systems have uncertain configurations due to large uncertainties on their inclinations. The fact that systems display ellipsoidal variations has a significant impact on the inclination errors. From the dynamical masses, luminosities, and radii, we provide LMC-based empirical mass-luminosity and mass-radius relations, and we compare them to other relations given for the Galaxy, the LMC, and the SMC. These relations differ for different mass ranges, but do not seem to depend on the metallicity regimes. We also compare the dynamical, spectroscopic, and evolutionary masses of the stars in our sample. While the dynamical and spectroscopic masses agree with each other, the evolutionary masses are systematically higher, at least for stars in semi-detached systems. This suggests that the mass discrepancy can be partly explained by past or ongoing interactions between the stars.

preprint2019arXiv

The young massive SMC cluster NGC 330 seen by MUSE. I. Observations and stellar content

A majority of massive stars are part of binary systems, a large fraction of which will inevitably interact during their lives. Binary-interaction products (BiPs), i.e. stars affected by such interaction, are expected to be commonly present in stellar populations. BiPs are thus a crucial ingredient in the understanding of stellar evolution. We aim to identify and characterize a statistically significant sample of BiPs by studying clusters of 10-40 Myr, an age at which binary population models predict the abundance of BiPs to be highest. One example of such a cluster is NGC 330 in the Small Magellanic Cloud. Using MUSE WFM-AO observations of NGC 330, we resolve the dense cluster core for the first time and are able to extract spectra of its entire massive star population. We develop an automated spectral classification scheme based on the equivalent widths of spectral lines in the red part of the spectrum. We characterize the massive star content of the core of NGC 330 which contains more than 200 B stars, 2 O stars, 6 A-type supergiants and 11 red supergiants. We find a lower limit on the Be star fraction of $32 \pm 3\%$ in the whole sample. It increases to at least $46 \pm 10\%$ when only considering stars brighter than V=17 mag. We estimate an age of the cluster core between 35 and 40 Myr and a total cluster mass of $88^{+17}_{-18} \times 10^3 M_{\odot}$. We find that the population in the cluster core is different than the population in the outskirts: while the stellar content in the core appears to be older than the stars in the outskirts, the Be star fraction and the observed binary fraction are significantly higher. Furthermore, we detect several BiP candidates that will be subject of future studies.

preprint2019arXiv

The young stellar content of the giant HII regions M8, G333.6-0.2, and NGC6357 with VLT/KMOS

Context: The identification and characterisation of populations of young massive stars in (giant) HII regions provides important constraints on i) the formation process of massive stars and their early feedback on the environment, and ii) the initial conditions for population synthesis models predicting the evolution of ensembles of stars. Aims: We identify and characterise the stellar populations of the following young giant HII regions: M8, G333.6-0.2, and NGC6357. Methods: We acquired H- and K-band spectra of around 200 stars using The K-band KMOS on the ESO Very Large Telescope. The targets for M8 and NGC6357 were selected from the MYStIX project, which combines X-ray observations with near-infrared and mid-infrared data. For G333.6-0.2, the sample selection is based on the near-infrared colours combined with X-ray data. We introduce an automatic spectral classification method in order to obtain temperatures and luminosities for the observed stars. We analyse the stellar populations using their photometric, astrometric, and spectroscopic properties and compared the position of the stars in the Hertzprung-Russell diagram with stellar evolution models to constrain their ages and mass ranges. Results: We confirm the presence of candidate ionising sources in the three regions and report new ones, including the first spectroscopically identified O stars in G333.6-0.2. In M8 and NGC6357, two populations are identified: i) OB main-sequence stars ($M > 5~\rm{M_{\odot}}$) and ii) pre-main sequence stars ($M\approx0.5-5~\rm{M_{\odot}}$). The ages of the clusters are $\sim$1-3~Myr, $< 3$~Myr, and $\sim$0.5-3~Myr for M8, G333.6-0.2, and NGC6357, respectively. We show that MYStIX selected targets have $>$ 90\% probability of being members of the HII region, whereas a selection based on near infrared (NIR) colours leads to a membership probability of only $\sim$70\%.

preprint2018arXiv

Response to comment on "An excess of massive stars in the local 30 Doradus starburst"

Farr and Mandel reanalyse our data, finding initial-mass-function slopes for high mass stars in 30 Doradus that agree with our results. However, their reanalysis appears to underpredict the observed number of massive stars. Their technique results in more precise slopes than in our work, strengthening our conclusion that there is an excess of massive stars above $30\,\mathrm{M}_\odot$ in 30 Doradus.

preprint2016arXiv

A new prescription for the mass-loss rates of WC and WO stars

We present a new empirical prescription for the mass-loss rates of carbon and oxygen sequence Wolf-Rayet stars as a function of their luminosity, surface chemical composition, and initial metallicity. The new prescription is based on results of detailed spectral analyses of WC and WO stars, and improves the often applied Nugis & Lamers (2000) relation. We find that the mass-loss rates of WC and WO stars (with $X=0$ and $Y < 0.98$) can be expressed as $\log{\dot{M}} = -9.20 + 0.85\log{(L/L_{\odot})} + 0.44\log{Y} + 0.25\log{(Z_{\mathrm{Fe}}/Z_{\mathrm{Fe}, \odot})}$. This relation is based on mass-loss determinations that assume a volume-filling factor of 0.1, but the prescription can easily be scaled to account for other volume-filling factors. The residual of the fit is $σ= 0.06$ dex. We investigated whether the relation can also describe the mass loss of hydrogen-free WN stars and showed that it can when an adjustement of the metallicty dependence ($\log{\dot{M}} \propto 1.3\log{(Z_{\mathrm{Fe}}/Z_{\mathrm{Fe}, \odot})}$) is applied. Compared to Nugis & Lamers (2000), $\dot{M}$ is less sensitive to the luminosity and the surface abundance, implying a stronger mass loss of massive stars in their late stages of evolution. The modest metallicity dependence implies that if WC or WO stars are formed in metal deficient environments, their mass-loss rates are higher than currently anticipated. These effects may result in a larger number of type Ic supernovae and less black holes to be formed, and may favour the production of superluminous type Ic supernovae through interaction with C and O rich circumstellar material or the dense stellar wind.

preprint2016arXiv

Evidence of magnetic field decay in massive main-sequence stars

A significant fraction of massive main-sequence stars show strong, large-scale magnetic fields. The origin of these fields, their lifetimes, and their role in shaping the characteristics and evolution of massive stars are currently not well understood. We compile a catalogue of 389 massive main-sequence stars, 61 of which are magnetic, and derive their fundamental parameters and ages. The two samples contain stars brighter than magnitude 9 in the V band and range in mass between 5 and 100 Msun. We find that the fractional main-sequence age distribution of all considered stars follows what is expected for a magnitude limited sample, while that of magnetic stars shows a clear decrease towards the end of the main sequence. This dearth of old magnetic stars is independent of the choice of adopted stellar evolution tracks, and appears to become more prominent when considering only the most massive stars. We show that the decreasing trend in the distribution is significantly stronger than expected from magnetic flux conservation. We also find that binary rejuvenation and magnetic suppression of core convection are unlikely to be responsible for the observed lack of older magnetic massive stars, and conclude that its most probable cause is the decay of the magnetic field, over a time span longer than the stellar lifetime for the lowest considered masses, and shorter for the highest masses. We then investigate the spin-down ages of the slowly rotating magnetic massive stars and find them to exceed the stellar ages by far in many cases. The high fraction of very slowly rotating magnetic stars thus provides an independent argument for a decay of the magnetic fields.

preprint2016arXiv

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

We introduce a HST/STIS stellar census of R136a, the central ionizing star cluster of 30 Doradus. We present low resolution far-ultraviolet STIS/MAMA spectroscopy of R136 using 17 contiguous 52x0.2 arcsec slits which together provide complete coverage of the central 0.85 parsec (3.4 arcsec). We provide spectral types of 90% 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 CIV 1548-51, including 16 O2-3 stars. For the first time we spectroscopically classify all Weigelt & 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.3_-0.7 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 HeII 1640 emission line. This emission is totally dominated by very massive stars with initial masses above ~100 Msun. The presence of strong HeII 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 100 Msun. We include montages of ultraviolet spectroscopy for LMC O stars in the Appendix. Future studies in this series will focus on optical STIS/CCD medium resolution observations.

preprint2015arXiv

B fields in OB stars (BOB): Detection of a strong magnetic field in the O9.7 V star HD54879

The number of magnetic stars detected among massive stars is small; nevertheless, the role played by the magnetic field in stellar evolution cannot be disregarded. Links between line profile variability, enhancements/depletions of surface chemical abundances, and magnetic fields have been identified for low-mass B-stars, but for the O-type domain this is almost unexplored. Based on FORS2 and HARPS spectropolarimetric data, we present the first detection of a magnetic field in HD54879, a single slowly rotating O9.7 V star. Using two independent and different techniques we obtained the firm detection of a surface average longitudinal magnetic field with a maximum amplitude of about 600 G, in modulus. A quantitative spectroscopic analysis of the star with the stellar atmosphere code FASTWIND results in an effective temperature and a surface gravity of 33000$\pm1000$ K and 4.0$\pm0.1$ dex. The abundances of carbon, nitrogen, oxygen, silicon, and magnesium are found to be slightly lower than solar, but compatible within the errors. We investigate line-profile variability in HD54879 by complementing our spectra with spectroscopic data from other recent OB-star surveys. The photospheric lines remain constant in shape between 2009 and 2014, although H$α$ shows a variable emission. The H$α$ emission is too strong for a standard O9.7 V and is probably linked to the magnetic field and the presence of circumstellar material. Its normal chemical composition and the absence of photospheric line profile variations make HD54879 the most strongly magnetic, non-variable single O-star detected to date.

preprint2015arXiv

B fields in OB stars (BOB): low-resolution FORS2 spectropolarimetry of the first sample of 50 massive stars

Within the context of the collaboration "B fields in OB stars (BOB)", we used the FORS2 low-resolution spectropolarimeter to search for a magnetic field in 50 massive stars, including two reference magnetic massive stars. Because of the many controversies of magnetic field detections obtained with the FORS instruments, we derived the magnetic field values with two completely independent reduction and analysis pipelines. We compare and discuss the results obtained from the two pipelines. We obtained a general good agreement, indicating that most of the discrepancies on magnetic field detections reported in the literature are caused by the interpretation of the significance of the results (i.e., 3-4 sigma detections considered as genuine, or not), instead of by significant differences in the derived magnetic field values. By combining our results with past FORS1 measurements of HD46328, we improve the estimate of the stellar rotation period, obtaining P = 2.17950+/-0.00009 days. For HD125823, our FORS2 measurements do not fit the available magnetic field model, based on magnetic field values obtained 30 years ago. We repeatedly detect a magnetic field for the O9.7V star HD54879, the HD164492C massive binary, and the He-rich star CPD -57 3509. We obtain a magnetic field detection rate of 6+/-4%, while by considering only the apparently slow rotators we derive a detection rate of 8+/-5%, both comparable with what was previously reported by other similar surveys. We are left with the intriguing result that, although the large majority of magnetic massive stars is rotating slowly, our detection rate is not a strong function of the stellar rotational velocity.

preprint2015arXiv

Discovery of the massive overcontact binary VFTS 352: Evidence for enhanced internal mixing

The contact phase expected to precede the coalescence of two massive stars is poorly characterized due to the paucity of observational constraints. Here we report on the discovery of VFTS 352, an O-type binary in the 30 Doradus region, as the most massive and earliest spectral type overcontact system known to date. We derived the 3D geometry of the system, its orbital period $P_{\rm orb}=1.1241452(4)$ d, components' effective temperatures -- $T_1=42\,540\pm280$ K and $T_2=41\,120\pm290$ K -- and dynamical masses -- $M_1=28.63\pm0.30 M_{\odot}$ and $M_2 = 28.85\pm0.30 M_{\odot}$. Compared to single-star evolutionary models, the VFTS 352 components are too hot for their dynamical masses by about 2700 and 1100 K, respectively. These results can be explained naturally as a result of enhanced mixing, theoretically predicted to occur in very short-period tidally-locked systems. The VFTS 352 components are two of the best candidates identified so far to undergo this so-called chemically homogeneous evolution. The future of VFTS 352 is uncertain. If the two stars merge, a very rapidly rotating star will be produced. Instead, if the stars continue to evolve homogeneously and keep shrinking within their Roche Lobes, coalescence can be avoided. In this case, tides may counteract the spin down by winds such that the VFTS 352 components may, at the end of their life, fulfill the requirements for long gamma-ray burst (GRB) progenitors in the collapsar scenario. Independently of whether the VFTS 352 components become GRB progenitors, this scenario makes VFTS 352 interesting as a progenitor of a black hole binary, hence as a potential gravitational wave source through black hole-black hole merger.

preprint2015arXiv

HST astrometry in the 30 Doradus region: measuring proper motions of individual stars in the Large Magellanic Cloud

We present measurements of positions and relative proper motions in the 30 Doradus region of the Large Magellanic Cloud (LMC). We detail the construction of a single-epoch astrometric reference frame, based on specially-designed observations obtained with the two main imaging instruments ACS/WFC and WFC3/UVIS onboard the Hubble Space Telescope (HST). Internal comparisons indicate a sub milli-arc-second (mas) precision in the positions and the presence of semi-periodic systematics with a mean amplitude of ~0.8 mas. We combined these observations with numerous archival images taken with WFPC2 and spanning 17 years. The precision of the resulting proper motions for well-measured stars around the massive cluster R 136 can be as good as ~20 microarcsec/yr, although the true accuracy of proper motions is generally lower due to the residual systematic errors. The observed proper-motion dispersion for our highest-quality measurements is ~0.1 mas/yr. Our catalog of positions and proper motions contains 86,590 stars down to V~25 and over a total area of ~70 square arcmin. We examined the proper motions of 105 relatively bright stars and identified a total of 6 candidate runaway stars. We are able to tentatively confirm the runaway status of star VFTS 285, consistent with the findings from line-of-sight velocities, and to show that this star has likely been ejected from R 136. This study demonstrates that with HST it is now possible to reliably measure proper motions of individual stars in the nearest dwarf galaxies such as the LMC.

preprint2015arXiv

Hubble Tarantula Treasury Project. III. Photometric Catalog and Resulting Constraints on the Progression of Star Formation in the 30 Doradus Region

We present and describe the astro-photometric catalog of more than 800,000 sources found in the Hubble Tarantula Treasury Project (HTTP). HTTP is a Hubble Space Telescope (HST) Treasury program designed to image the entire 30 Doradus region down to the sub-solar (~0.5 solar masses) mass regime using the Wide Field Camera 3 (WFC3) and the Advanced Camera for Surveys (ACS). We observed 30 Doradus in the near ultraviolet (F275W, F336W), optical (F555W, F658N, F775W), and near infrared (F110W, F160W) wavelengths. The stellar photometry was measured using point-spread function (PSF) fitting across all the bands simultaneously. The relative astrometric accuracy of the catalog is 0.4 mas. The astro-photometric catalog, results from artificial star experiments and the mosaics for all the filters are available for download. Color-magnitude diagrams are presented showing the spatial distributions and ages of stars within 30 Dor as well as in the surrounding fields. HTTP provides the first rich and statistically significant sample of intermediate and low mass pre-main sequence candidates and allows us to trace how star formation has been developing through the region. The depth and high spatial resolution of our analysis highlight the dual role of stellar feedback in quenching and triggering star formation on the giant HII region scale. Our results are consistent with stellar sub-clustering in a partially filled gaseous nebula that is offset towards our side of the Large Magellanic Cloud.

preprint2015arXiv

Masses of the components of SB2 binaries observed with Gaia. II. Masses derived from PIONIER interferometric observations for Gaia validation

In anticipation of the Gaia astrometric mission, a sample of spectroscopic binaries is being observed since 2010 with the Sophie spectrograph at the Haute--Provence Observatory. Our aim is to derive the orbital elements of double-lined spectroscopic binaries (SB2s) with an accuracy sufficient to finally obtain the masses of the components with relative errors as small as 1 % when combined with Gaia astrometric measurements. In order to validate the masses derived from Gaia, interferometric observations are obtained for three SB2s in our sample with F-K components: HIP 14157, HIP 20601 and HIP 117186. The masses of the six stellar components are derived. Due to its edge-on orientation, HIP 14157 is probably an eclipsing binary. We note that almost all the derived masses are a few percent larger than the expectations from the standard spectral-type-mass calibration and mass-luminosity relation. Our calculation also leads to accurate parallaxes for the three binaries, and the Hipparcos parallaxes are confirmed.

preprint2015arXiv

Massive stars on the verge of exploding: the properties of oxygen sequence Wolf-Rayet stars

Context. Oxygen sequence Wolf-Rayet (WO) stars represent a very rare stage in the evolution of massive stars. Their spectra show strong emission lines of helium-burning products, in particular highly ionized carbon and oxygen. The properties of WO stars can be used to provide unique constraints on the (post-)helium burning evolution of massive stars, as well as their remaining lifetime and the expected properties of their supernovae. Aims. We aim to homogeneously analyse the currently known presumed-single WO stars to obtain the key stellar and outflow properties and to constrain their evolutionary state. Methods. We use the line-blanketed non-local thermal equilibrium atmosphere code cmfgen to model the X-Shooter spectra of the WO stars and deduce the atmospheric parameters. We calculate dedicated evolutionary models to determine the evolutionary state of the stars. Results. The WO stars have extremely high temperatures that range from 150 kK to 210 kK, and very low surface helium mass fractions that range from 44% down to 14%. Their properties can be reproduced by evolutionary models with helium zero-age main sequence masses of $M_{\mathrm{He, ini}} = 15-25 M_{\odot}$ that exhibit a fairly strong (on the order of a few times $10^{-5} M_{\odot} \mathrm{yr}^{-1}$), homogeneous ($f_\mathrm{c} > 0.3$) stellar wind. Conclusions. WO stars represent the final evolutionary stage of stars with estimated initial masses of $M_{\mathrm{ini}} = 40-60 M_{\odot}$. They are post core-helium burning and predicted to explode as type Ic supernovae within a few thousand years.

preprint2015arXiv

Molecfit: A general tool for telluric absorption correction II. Quantitative evaluation on ESO-VLT X-Shooter spectra

Context: Absorption by molecules in the Earth's atmosphere strongly affects ground-based astronomical observations. The resulting absorption line strength and shape depend on the highly variable physical state of the atmosphere, i.e. pressure, temperature, and mixing ratio of the different molecules involved. Usually, supplementary observations of so-called telluric standard stars (TSS) are needed to correct for this effect, which is expensive in terms of telescope time. We have developed the software package molecfit to provide synthetic transmission spectra based on parameters obtained by fitting narrow ranges of the observed spectra of scientific objects. These spectra are calculated by means of the radiative transfer code LBLRTM and an atmospheric model. In this way, the telluric absorption correction for suitable objects can be performed without any additional calibration observations of TSS. Aims: We evaluate the quality of the telluric absorption correction using molecfit with a set of archival ESO-VLT X-Shooter visible and near-infrared spectra. Methods: Thanks to the wavelength coverage from the U to the K band, X-Shooter is well suited to investigate the quality of the telluric absorption correction with respect to the observing conditions, the instrumental set-up, input parameters of the code, the signal-to-noise of the input spectrum, and the atmospheric profiles. These investigations are based on two figures of merit, I_off and I_res, that describe the systematic offsets and the remaining small-scale residuals of the corrections. We also compare the quality of the telluric absorption correction achieved with moelcfit to the classical method based on a telluric standard star. (Abridged)

preprint2015arXiv

Molecfit: A general tool for telluric absorption correction. I. Method and application to ESO instruments

Context: The interaction of the light from astronomical objects with the constituents of the Earth's atmosphere leads to the formation of telluric absorption lines in ground-based collected spectra. Correcting for these lines, mostly affecting the red and infrared region of the spectrum, usually relies on observations of specific stars obtained close in time and airmass to the science targets, therefore using precious observing time. Aims: We present molecfit, a tool for correcting for telluric absorption lines based on synthetic modelling of the Earth's atmospheric transmission. Molecfit is versatile and can be used with data obtained with various ground-based telescopes and instruments. Methods: Molecfit combines a publicly available radiative transfer code, a molecular line database, atmospheric profiles, and various kernels to model the instrument line spread function. The atmospheric profiles are created by merging a standard atmospheric profile representative of a given observatory's climate, of local meteorological data, and of dynamically retrieved altitude profiles for temperature, pressure, and humidity. We discuss the various ingredients of the method, its applicability, and its limitations. We also show examples of telluric line correction on spectra obtained with a suite of ESO Very Large Telescope (VLT) instruments. Results: Compared to previous similar tools, molecfit takes the best results for temperature, pressure, and humidity in the atmosphere above the observatory into account. As a result, the standard deviation of the residuals after correction of unsaturated telluric lines is frequently better than 2% of the continuum. Conclusion: Molecfit is able to accurately model and correct for telluric lines over a broad range of wavelengths and spectral resolutions. (Abridged)

preprint2015arXiv

The evolution of rotating very massive stars with LMC composition

We present a dense model grid with tailored input chemical composition appropriate for the Large Magellanic Cloud. We use a one-dimensional hydrodynamic stellar evolution code, which accounts for rotation, transport of angular momentum by magnetic fields, and stellar wind mass loss to compute our detailed models. We calculate stellar evolution models with initial masses of 70-500 Msun and with initial surface rotational velocities of 0-550 km/s, covering the core-hydrogen burning phase of evolution. We find our rapid rotators to be strongly influenced by rotationally induced mixing of helium, with quasi-chemically homogeneous evolution occurring for the fastest rotating models. Above 160 Msun, homogeneous evolution is also established through mass loss, producing pure helium stars at core hydrogen exhaustion independent of the initial rotation rate. Surface nitrogen enrichment is also found for slower rotators, even for stars that lose only a small fraction of their initial mass. For models above 150 MZAMS, and for models in the whole considered mass range later on, we find a considerable envelope inflation due to the proximity of these models to their Eddington limit. This leads to a maximum zero-age main sequence surface temperature of 56000 K, at 180 Msun, and to an evolution of stars in the mass range 50-100 Msun to the regime of luminous blue variables in the HR diagram with high internal Eddington factors. Inflation also leads to decreasing surface temperatures during the chemically homogeneous evolution of stars above 180 Msun. The cool surface temperatures due to the envelope inflation in our models lead to an enhanced mass loss, which prevents stars at LMC metallicity from evolving into pair-instability supernovae. The corresponding spin-down will also prevent very massive LMC stars to produce long-duration gamma-ray bursts, which might, however, originate from lower masses.

preprint2015arXiv

The mass of the very massive binary WR21a

We present multi-epoch spectroscopic observations of the massive binary system WR21a, which include the January 2011 periastron passage. Our spectra reveal multiple SB2 lines and facilitate an accurate determination of the orbit and the spectral types of the components. We obtain minimum masses of $64.4\pm4.8 \ M_{\odot}$ and $36.3\pm1.7 \ M_{\odot}$ for the two components of WR21a. Using disentangled spectra of the individual components, we derive spectral types of O3/WN5ha and O3Vz~((f*)) for the primary and secondary, respectively. Using the spectral type of the secondary as an indication for its mass, we estimate an orbital inclination of $i=58.8\pm2.5^{\mathrm{o}}$ and absolute masses of $103.6\pm10.2 \ M_{\odot}$ and $58.3\pm3.7 \ M_{\odot}$, in agreement with the luminosity of the system. The spectral types of the WR21a components indicate that the stars are very young (1$-$2 Myr), similar to the age of the nearby Westerlund 2 cluster. We use evolutionary tracks to determine the mass-luminosity relation for the total system mass. We find that for a distance of 8 kpc and an age of 1.5 Myr, the derived absolute masses are in good agreement with those from evolutionary predictions.

preprint2015arXiv

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

TLUSTY non-LTE model atmosphere calculations have been used to determine atmospheric parameters and nitrogen (N) abundances for 34 single and 18 binary B-type supergiants (BSGs). The effects of flux contribution from an unseen secondary were considered for the binary sample. We present the first systematic study of the incidence of binarity for a sample of BSGs across the theoretical terminal age main sequence (TAMS). To account for the distribution of effective temperatures of the BSGs it may be necessary to extend the TAMS to lower temperatures. This is consistent with the derived distribution of mass discrepancies, projected rotational velocities (vsini) and N abundances, provided that stars cooler than this temperature are post RSG objects. For the BSGs in the Tarantula and previous FLAMES surveys, most have small vsini. About 10% have larger vsini (>100 km/s) but surprisingly these show little or no N enhancement. All the cooler BSGs have low vsini of <70km/s and high N abundance estimates, implying that either bi-stability braking or evolution on a blue loop may be important. A lack of cool binaries, possibly reflects the small sample size. Single star evolutionary models, which include rotation, can account for the N enhancement in both the single and binary samples. The detailed distribution of N abundances in the single and binary samples may be different, possibly reflecting differences in their evolutionary history. The first comparative study of single and binary BSGs has revealed that the main sequence may be significantly wider than previously assumed, extending to Teff=20000K. 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 evolution of BSGs. Full abstract in paper

preprint2015arXiv

The VLT-FLAMES Tarantula Survey XVIII. Classifications and radial velocities of the B-type stars

We present spectral classifications for 438 B-type stars observed as part of the VLT-FLAMES Tarantula Survey (VFTS) in the 30 Doradus region of the Large Magellanic Cloud. Radial velocities are provided for 307 apparently single stars, and for 99 targets with radial-velocity variations which are consistent with them being spectroscopic binaries. We investigate the spatial distribution of the radial velocities across the 30 Dor region, and use the results to identify candidate runaway stars. Excluding potential runaways and members of two older clusters in the survey region (SL 639 and Hodge 301), we determine a systemic velocity for 30 Dor of 271.6 +/- 12.2 km/s from 273 presumed single stars. Employing a 3-sigma criterion we identify nine candidate runaway stars (2.9% of the single stars with radial-velocity estimates). The projected rotational velocities of the candidate runaways appear to be significantly different to those of the full B-type sample, with a strong preference for either large (>345 km/s) or small (<65 km/s) rotational velocities. Of the candidate runaways, VFTS 358 (classified B0.5: V) has the largest differential radial velocity (-106.9 +/- 16.2 km/s), and a preliminary atmospheric analysis finds a significantly enriched nitrogen abundance of 12+log(N/H) > ~8.5. Combined with a large rotational velocity (vsini = 345 +/- 22 km/s), this is suggestive of past binary interaction for this star.

preprint2015arXiv

The VLT-FLAMES Tarantula survey XX. The nature of the X-ray bright emission line star VFTS 399

The stellar population of the 30 Doradus star-forming region in the Large Magellanic Cloud contains a subset of apparently single, rapidly rotating O-type stars. The physical processes leading to the formation of this cohort are currently uncertain. One member of this group, the late O-type star VFTS 399, is found to be unexpectedly X-ray bright for its bolometric luminosity - in this study we aim to determine its physical nature and the cause of this behaviour. We find VFTS 399 to be an aperiodic photometric variable with an apparent near-IR excess. Its optical spectrum demonstrates complex emission profiles in the lower Balmer series and select HeI lines - taken together these suggest an OeBe classification. The highly variable X-ray luminosity is too great to be produced by a single star, while the hard, non-thermal nature suggests the presence of an accreting relativistic companion. Finally, the detection of periodic modulation of the X-ray lightcurve is most naturally explained under the assumption that the accretor is a neutron star. VFTS 399 appears to be the first high-mass X-ray binary identified within 30 Dor. Comparison of the current properties of VFTS 399 to binary-evolution models suggests a progenitor mass in excess of 25Msun for the putative neutron star, which may host a magnetic field comparable in strength to those of magnetars. VFTS 399 is now the second member of the cohort of rapidly rotating `single' O-type stars in 30 Dor to show evidence of binary interaction resulting in spin-up, suggesting that this may be a viable evolutionary pathway for the formation of a subset of this stellar population.

preprint2015arXiv

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

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. By establishing the spin-rate distribution of a sizeable sample of O-type spectroscopic binaries and by comparing the distributions of binary sub-populations with one another as well as 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 massive stars spin rates. We use ground-based optical spectroscopy obtained in the framework of the VLT-FLAMES Tarantula Survey (VFTS) to establish the projected equatorial rotational velocities (\vrot) for components of 114 spectroscopic binaries in 30 Doradus. The \vrot\ values are derived from the full-width at half-maximum (FWHM) of a set of spectral lines, using a FWHM vs. \vrot\ calibration that we derive based on previous line analysis methods applied to single O-type stars in the VFTS sample. The overall \vrot\ distribution of the primary stars resembles that of single O-type stars in the VFTS, featuring a low-velocity peak (at $\vrot < 200$ kms) and a shoulder at intermediate velocities ($200 < \vrot < 300$ kms). The distributions of binaries and single stars however differ in two ways. First, the main peak at $\vrot \sim$100 kms is broader and slightly shifted toward higher spin rates in the binary distribution compared to that of the presumed-single stars. Second, the \vrot distribution of primaries lacks a significant population of stars spinning faster than 300 kms while such a population is clearly present in the single star sample.

preprint2015arXiv

The VLT-FLAMES Tarantula Survey XXII. Multiplicity properties of the B-type stars

We investigate the multiplicity properties of 408 B-type stars observed in the 30 Doradus region of the Large Magellanic Cloud with multi-epoch spectroscopy from the VLT-FLAMES Tarantula Survey (VFTS). We use a cross-correlation method to estimate relative radial velocities from the helium and metal absorption lines for each of our targets. Objects with significant radial-velocity variations (and with an amplitude larger than 16 km/s) are classified as spectroscopic binaries. We find an observed spectroscopic binary fraction (defined by periods of <10^3.5 d and mass ratios >0.1) for the B-type stars, f_B(obs) = 0.25 +/- 0.02, which appears constant across the field of view, except for the two older clusters (Hodge 301 and SL 639). These two clusters have significantly lower fractions of 0.08 +/- 0.08 and 0.10 +/- 0.09, respectively. Using synthetic populations and a model of our observed epochs and their potential biases, we constrain the intrinsic multiplicity properties of the dwarf and giant (i.e. relatively unevolved) B-type stars in 30 Dor. We obtain a present-day binary fraction f_B(true) = 0.58 +/- 0.11, with a flat period distribution. Within the uncertainties, the multiplicity properties of the B-type stars agree with those for the O stars in 30 Dor from the VFTS.

preprint2015arXiv

The VLT-FLAMES Tarantula Survey. XXIII: two massive double-lined binaries in 30~Doradus

We investigate the characteristics of two newly discovered short-period, double-lined, massive binary systems, VFTS 450 (O9.7$\;$II--Ib$\,$+$\,$O7::) and VFTS 652 (B1$\;$Ib$\,+\,$O9:$\;$III:). We perform model-atmosphere analyses to characterise the photospheric properties of both members of each binary (denoting the `primary' as the spectroscopically more conspicuous component). Radial velocities and optical photometry are used to estimate the binary-system parameters. We estimate $T_{\rm eff}=27$ kK, $\log{(g)}=2.9$ (cgs) for the VFTS 450 primary spectrum (34kK, 3.6: for the secondary spectrum); and $T_{\rm eff} = 22$kK, $\log{(g)}=2.8$ for the VFTS 652 primary spectrum (35kK, 3.7: for the secondary spectrum). Both primaries show surface nitrogen enrichments (of more than 1 dex for VFTS 652), and probable moderate oxygen depletions relative to reference LMC abundances. We determine orbital periods of 6.89d and 8.59d for VFTS 450 and VFTS 652, respectively, and argue that the primaries must be close to filling their Roche lobes. Supposing this to be the case, we estimate component masses in the range $\sim$20--50M$_\odot$. The secondary spectra are associated with the more massive components, suggesting that both systems are high-mass analogues of classical Algol systems, undergoing case-A mass transfer. Difficulties in reconciling the spectroscopic analyses with the light-curves and with evolutionary considerations suggest that the secondary spectra are contaminated by (or arise in) accretion disks.

preprint2014arXiv

B fields in OB stars (BOB): on the detection of weak magnetic fields in the two early B-type stars beta CMa and epsilon CMa

Within the context of the "B fields in OB stars (BOB)" collaboration, we used the HARPSpol spectropolarimeter to observe the early B-type stars beta CMa (HD44743; B1 II/III) and epsilon CMa (HD52089; B1.5 II). For both stars, we consistently detected the signature of a weak (<30 G in absolute value) longitudinal magnetic field. We determined the physical parameters of both stars and characterise their X-ray spectrum. For beta CMa, our mode identification analysis led to determining a rotation period of 13.6+/-1.2 days and of an inclination angle of the rotation axis of 57.6+/-1.7 degrees, with respect to the line of sight. On the basis of these measurements and assuming a dipolar field geometry, we derived a best fitting obliquity of ~22 degrees and a dipolar magnetic field strength (Bd) of ~100 G (60<Bd<230 G within 1 sigma), below what is typically found for other magnetic massive stars. For epsilon CMa we could only determine a lower limit on the dipolar magnetic field strength of 13 G. For this star, we determine that the rotation period ranges between 1.3 and 24 days. Both stars are expected to have a dynamical magnetosphere. We also conclude that both stars are most likely core hydrogen burning and that they have spent more than 2/3 of their main sequence lifetime. A histogram of the distribution of the dipolar magnetic field strength for the magnetic massive stars known to date does not show the magnetic field "desert" observed instead for intermediate-mass stars. The biases involved in the detection of (weak) magnetic fields in massive stars with the currently available instrumentation and techniques imply that weak fields might be more common than currently observed. Our results show that, if present, even relatively weak magnetic fields are detectable in massive stars and that more observational effort is probably still needed to properly access the magnetic field incidence.

preprint2014arXiv

HD 152246 - a new high-mass triple system and its basic properties

Analyses of multi-epoch, high-resolution (R ~ 50.000) optical spectra of the O-type star HD 152246 (O9 IV according to the most recent classification), complemented by a limited number of earlier published radial velocities, led to the finding that the object is a hierarchical triple system, where a close inner pair (Ba-Bb) with a slightly eccentric orbit (e = 0.11) and a period of 6.0049 days revolves in a 470-day highly eccentric orbit (e = 0.865) with another massive and brighter component A. The mass ratio of the inner system must be low since we were unable to find any traces of the secondary spectrum. The mass ratio A/(Ba+Bb) is 0.89. The outer system has recently been resolved using long-baseline interferometry on three occasions. The interferometry confirms the spectroscopic results and specifies elements of the system. Our orbital solutions, including the combined radial-velocity and interferometric solution indicate an orbital inclination of the outer orbit of 112° and stellar masses of 20.4 and 22.8 solar masses. We also disentangled the spectra of components A and Ba and compare them to synthetic spectra from two independent programmes, TLUSTY and FASTWIND. In either case, the fit was not satisfactory and we postpone a better determination of the system properties for a future study, after obtaining observations during the periastron passage of the outer orbit (the nearest chance being March 2015). For the moment, we can only conclude that component A is an O9 IV star with v*sin(i) = 210 +\- 10 km/s and effective temperature of 33000 +\- 500 K, while component Ba is an O9 V object with v*sin(i) = 65 +/- 3 km/s and T_eff = 33600 +\- 600 K.

preprint2014arXiv

Molecfit: A Package for Telluric Absorption Correction

Correcting for the sky signature usually requires supplementary calibration data which are very expensive in terms of telescope time. In addition, the scheduling flexibility is restricted as these data have to be taken usually directly before/after the science observations due to the high variability of the telluric absorption which depends on the state and the chemical composition of the atmosphere at the time of observations. Therefore, a tool for sky correction, which does not require this supplementary calibration data, saves a significant amount of valuable telescope time and increases its efficiency. We developed a software package aimed at performing telluric feature corrections on the basis of synthetic absorption spectra.

preprint2014arXiv

Relating jet structure to photometric variability: the Herbig Ae star HD 163296

Herbig Ae/Be stars are intermediate-mass pre-main sequence stars surrounded by circumstellar dust disks. Some are observed to produce jets, whose appearance as a sequence of shock fronts (knots) suggests a past episodic outflow variability. This "jet fossil record" can be used to reconstruct the outflow history. We present the first optical to near-infrared (NIR) VLT/X-shooter spectra of the jet from the Herbig Ae star HD 163296. We determine physical conditions in the knots, as well as their kinematic "launch epochs". Knots are formed simultaneously on either side of the disk, with a regular interval of ~16 yr. The velocity dispersion versus jet velocity and the energy input are comparable in both lobes. However, the mass loss rate, velocity, and shock conditions are asymmetric. We find Mjet/Macc ~ 0.01-0.1, consistent with magneto-centrifugal jet launching models. No evidence for dust is found in the high-velocity jet, suggesting it is launched within the sublimation radius (<0.5 au). The jet inclination measured from proper motions and radial velocities confirms it is perpendicular to the disk. A tentative relation is found between the structure of the jet and the photometric variability of the source. Episodes of NIR brightening were previously detected and attributed to a dusty disk wind. We report for the first time significant optical fadings lasting from a few days up to a year, coinciding with the NIR brightenings. These are likely caused by dust lifted high above the disk plane; this supports the disk wind scenario. The disk wind is launched at a larger radius than the high-velocity atomic jet, although their outflow variability may have a common origin. No significant relation between outflow and accretion variability could be established. Our findings confirm that this source undergoes periodic ejection events, which may be coupled with dust ejections above the disk plane.

preprint2014arXiv

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, $v_e\sin i$, of $\sim$330 O-type objects, i.e. $\sim$210 spectroscopic single stars and $\sim$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 $v_e\sin i$ distributions of the presumed-single stars and primaries in 30 Dor is a low-velocity peak at around 100\,$\rm{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\,$\rm{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.

preprint2014arXiv

Southern Massive Stars at High Angular Resolution: Observational Campaign and Companion Detection

Multiplicity is one of the most fundamental observable properties of massive O-type stars and offers a promising way to discriminate between massive star formation theories. Nevertheless, companions at separations between 1 and 100 mas remain mostly unknown due to intrinsic observational limitations. [...] The Southern MAssive Stars at High angular resolution survey (SMASH+) was designed to fill this gap by providing the first systematic interferometric survey of Galactic massive stars. We observed 117 O-type stars with VLTI/PIONIER and 162 O-type stars with NACO/SAM, respectively probing the separation ranges 1-45 and 30-250mas and brightness contrasts of Delta H < 4 and Delta H < 5. Taking advantage of NACO's field-of-view, we further uniformly searched for visual companions in an 8''-radius down to Delta H = 8. This paper describes the observations and data analysis, reports the discovery of almost 200 new companions in the separation range from 1mas to 8'' and presents the catalog of detections, including the first resolved measurements of over a dozen known long-period spectroscopic binaries. Excluding known runaway stars for which no companions are detected, 96 objects in our main sample (DEC < 0 deg; H<7.5) were observed both with PIONIER and NACO/SAM. The fraction of these stars with at least one resolved companion within 200mas is 0.53. Accounting for known but unresolved spectroscopic or eclipsing companions, the multiplicity fraction at separation < 8'' increases to f_m = 0.91 +/- 0.03. The fraction of luminosity class V stars that have a bound companion reaches 100% at 30mas while their average number of physically connected companions within 8'' is f_c = 2.2 +/- 0.3. This demonstrates that massive stars form nearly exclusively in multiple systems. Additionally, the nine non-thermal (NT) radio emitters observed by SMASH+ are all resolved [...]

preprint2014arXiv

The B Fields in OB Stars (BOB) Survey

The B fields in OB stars (BOB) survey is an ESO large programme collecting spectropolarimetric observations for a large number of early-type stars in order to study the occurrence rate, properties, and ultimately the origin of magnetic fields in massive stars. As of July 2014, a total of 98 objects were observed over 20 nights with FORS2 and HARPSpol. Our preliminary results indicate that the fraction of magnetic OB stars with an organised, detectable field is low. This conclusion, now independently reached by two different surveys, has profound implications for any theoretical model attempting to explain the field formation in these objects. We discuss in this contribution some important issues addressed by our observations (e.g., the lower bound of the field strength) and the discovery of some remarkable objects.

preprint2014arXiv

The BinaMIcS project: understanding the origin of magnetic fields in massive stars through close binary systems

It is now well established that a fraction of the massive (M>8 Msun) star population hosts strong, organised magnetic fields, most likely of fossil origin. The details of the generation and evolution of these fields are still poorly understood. The BinaMIcS project takes an important step towards the understanding of the interplay between binarity and magnetism during the stellar formation and evolution, and in particular the genesis of fossil fields, by studying the magnetic properties of close binary systems. The components of such systems are most likely formed together, at the same time and in the same environment, and can therefore help us to disentangle the role of initial conditions on the magnetic properties of the massive stars from other competing effects such as age or rotation. We present here the main scientific objectives of the BinaMIcS project, as well as preliminary results from the first year of observations from the associated ESPaDOnS and Narval spectropolarimetric surveys.

preprint2014arXiv

The properties of single WO stars

The enigmatic oxygen sequence Wolf-Rayet (WO) stars represent a very late stage in massive star evolution, although their exact nature is still under debate. The spectra of most of the WO stars have never been analysed through detailed modelling with a non-local thermodynamic equilibrium expanding atmosphere code. Here we present preliminary results of the first homogeneous analysis of the (apparently) single WOs.

preprint2014arXiv

The properties of ten O-type stars in the low-metallicity galaxies IC 1613, WLM and NGC 3109

Massive stars likely played an important role in the reionization of the Universe, and the formation of the first black holes. Massive stars in low-metallicity environments in the local Universe are reminiscent of their high redshift counterparts. In a previous paper, we reported on indications that the stellar winds of low-metallicity O stars may be stronger than predicted, which would challenge the current paradigm of massive star evolution. In this paper, we aim to extend our initial sample of six O stars in low-metallicity environments by four. We aim to derive their stellar and wind parameters, and compare these to radiation-driven wind theory and stellar evolution models. We have obtained intermediate-resolution VLT/X-Shooter spectra of our sample of stars. We derive the stellar parameters by fitting synthetic fastwind line profiles to the VLT/X-Shooter spectra using a genetic fitting algoritm. We compare our parameters to evolutionary tracks and obtain evolutionary masses and ages. We also investigate the effective temperature versus spectral type calibration for SMC and lower metallicities. Finally, we reassess the wind momentum versus luminosity diagram. The derived parameters of our target stars indicate stellar masses that reach values of up to 50 $M_{\odot}$. The wind strengths of our stars are, on average, stronger than predicted from radiation-driven wind theory and reminiscent of stars with an LMC metallicity. We discuss indications that the iron content of the host galaxies is higher than originally thought and is instead SMC-like. We find that the discrepancy with theory is lessened, but remains significant for this higher metallicity. This may imply that our current understanding of the wind properties of massive stars, both in the local universe as well as at cosmic distances, remains incomplete.

preprint2014arXiv

The Three body system delta Circini

Delta Cir is known as an O7.5 III eclipsing and spectroscopic binary with an eccentric orbit. Penny et al. discovered the presence of a third component in the IUE spectra. The eclipsing binary and the third body revolve around a common centre of gravity with a period of 1644 days in an eccentric orbit with the semimajor axis of 10 AU. We demonstrate the presence of the apsidal-line rotation with a period of about 141 years, which is considerably longer than its theoretically predicted value, based on the published radii of the binary components derived from the Hipparchos Hp light curve. However, our new solution of the same light curve resulted in smaller radii and a better agreement between the observed and predicted period of the apsidal- line advance. There are indications that the third body is a binary. The object was resolved by VLTI with the PIONIER combiner; in June 2012 the separation was 3.78 mas, magnitude difference in the H region 1.75 mag. This result means that (assuming the distance 770 pc) the inclination of the long orbit is 87.7 degrees.

preprint2014arXiv

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, of which 213 O-type are of sufficient 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. (2) A group of extremely rapidly rotating main-sequence objects has been isolated, including the largest $v\sin i$ values known, the spatial and radial-velocity distributions of which suggest ejection from the two principal ionizing clusters. (3) Several new examples of the evolved, rapidly rotating Onfp class show similar evidence. (4) No fewer than 48 members of the Vz category, hypothesized to be on or near the ZAMS, are found in this sample; in contrast to the rapid rotators, they are strongly concentrated to the ionizing clusters, 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; there are also six objects with very high-quality data but no trace of either mission feature, presenting new challenges to physical interpretations. (6) Five spectra with morphologically enhanced nitrogen lines have been detected. Absolute visual magnitudes have been derived for each star with individual extinction laws, and composite HRDs provide evidence of the multiple generations present in this field. Associations with X-ray sources are noted. Further analyses of this unique dataset underway will provide new insights into the evolution of massive stars and starburst clusters.

preprint2014arXiv

The VLT-FLAMES Tarantula Survey XVI. The optical+NIR extinction laws in 30 Doradus and the photometric determination of the effective temperatures of OB stars

Context: The commonly used extinction laws of Cardelli et al. (1989) have limitations that, among other issues, hamper the determination of the effective temperatures of O and early B stars from optical+NIR photometry. Aims: We aim to develop a new family of extinction laws for 30 Doradus, check their general applicability within that region and elsewhere, and apply them to test the feasibility of using optical+NIR photometry to determine the effective temperature of OB stars. Methods: We use spectroscopy and NIR photometry from the VLT-FLAMES Tarantula Survey and optical photometry from HST/WFC3 of 30 Doradus and we analyze them with the software code CHORIZOS using different assumptions such as the family of extinction laws. Results: We derive a new family of optical+NIR extinction laws for 30 Doradus and confirm its applicability to extinguished Galactic O-type systems. We conclude that by using the new extinction laws it is possible to measure the effective temperatures of OB stars with moderate uncertainties and only a small bias, at least up to E(4405-5495) ~ 1.5 mag.

preprint2014arXiv

The VLT-FLAMES Tarantula Survey XVII. Physical and wind properties of massive stars at the top of the main sequence

The evolution and fate of very massive stars (VMS) is tightly connected to their mass-loss properties. Their initial and final masses differ significantly as a result of mass loss. VMS have strong stellar winds and extremely high ionising fluxes, which are thought to be critical sources of both mechanical and radiative feedback in giant Hii regions. However, how VMS mass-loss properties change during stellar evolution is poorly understood. In the framework of the VLT-Flames Tarantula Survey (VFTS), we explore the mass-loss transition region from optically thin O to denser WNh star winds, thereby testing theoretical predictions. To this purpose we select 62 O, Of, Of/WN, and WNh stars, an unprecedented sample of stars with the highest masses and luminosities known. We perform a spectral analysis of optical VFTS as well as near-infrared VLT/SINFONI data using the non-LTE radiative transfer code CMFGEN to obtain stellar and wind parameters. For the first time, we observationally resolve the transition between optically thin O and optically thick WNh star winds. Our results suggest the existence of a kink between both mass-loss regimes, in agreement with recent MC simulations. For the optically thick regime, we confirm the steep dependence on the Eddington factor from previous theoretical and observational studies. The transition occurs on the MS near a luminosity of 10^6.1Lsun, or a mass of 80...90Msun. Above this limit, we find that - even when accounting for moderate wind clumping (with f = 0.1) - wind mass-loss rates are enhanced with respect to standard prescriptions currently adopted in stellar evolution calculations. We also show that this results in substantial helium surface enrichment. Based on our spectroscopic analyses, we are able to provide the most accurate ionising fluxes for VMS known to date, confirming the pivotal role of VMS in ionising and shaping their environments.

preprint2013arXiv

On the nature of the WO3 star DR1 in IC 1613

We present the results of a quantitative spectroscopic analysis of the oxygen-sequence Wolf- Rayet star DR1 in the low-metallicity galaxy IC 1613. Our models suggest that the strong oxygen emission lines are the result of the high temperature of this WO3 star and do not necessarily reflect a more advanced evolutionary stage than WC stars.

preprint2013arXiv

On the nature of WO stars: a quantitative analysis of the WO3 star DR1 in IC 1613

Context. Oxygen sequence Wolf-Rayet (WO) stars are thought to represent the final evolutionary stage of the most massive stars. The characteristic strong O vi emission possibly originates from an enhanced oxygen abundance in the stellar wind. Alternatively, the O vi emission can be caused by the high temperature of these stars, in which case the WO stars are the high-temperature extension of the more common carbon sequence Wolf-Rayet (WC) stars. Aims. By constraining the physical properties and evolutionary status of DR1, a WO star in the low-metallicity Local Group dwarf galaxy IC 1613 and one of only two objects of its class known in a SMC-like metallicity environment, we aim to investigate the nature of WO stars and their evolutionary connection with WC stars. Methods. We use the non-LTE atmosphere code cmfgen to model the observed spectrum of DR1 and to derive its stellar and wind parameters. We compare our values with other studies of WC and WO stars, as well as with the predictions of evolutionary models. We also model the surrounding nebula using the photo-ionization code cloudy. Results. The oxygen and carbon abundances that we obtain are comparable to values found for WC stars. The temperature and luminosity are, however, higher than those of WC stars. DR1 is embedded in the hottest known H ii region in the Local Group. The nebular properties can be consistently reproduced by cloudy models adopting DR1 as central ionizing source. Conclusions. Comparison of the abundances and temperature of DR1 with core helium-burning models show that DR1 is currently well into the second half of helium burning. If the properties of DR1 are representative for the WO class, it would imply that WO stars are the high-temperature and high-luminosity extension of the WC stars, and do not necessarily represent a later evolutionary stage.

preprint2013arXiv

R144 revealed as a double-lined spectroscopic binary

R144 is a WN6h star in the 30 Doradus region. It is suspected to be a binary because of its high luminosity and its strong X-ray flux, but no periodicity could be established so far. Here, we present new Xshooter multi-epoch spectroscopy of R144 obtained at the ESO Very Large Telescope (VLT). We detect variability in position and/or shape of all the spectral lines. We measure radial velocity variations with an amplitude larger than 250 km/s in NIV and NV lines. Furthermore, the NIII and NV line Doppler shifts are anti-correlated and the NIV lines show a double-peaked profile on six of our seven epochs. We thus conclude that R144 is a double-lined spectroscopic binary. Possible orbital periods range from 2 to 6 months, although a period up to one year is allowed if the orbit is highly eccentric. We estimate the spectral types of the components to be WN5-6h and WN6-7h, respectively. The high luminosity of the system (log Lbol/Lsun ~ 6.8) suggests a present-day total mass content in the range of about 200 to 300 Msun, depending on the evolutionary stage of the components. This makes R144 the most massive binary identified so far, with a total mass content at birth possibly as large as 400 Msun. We briefly discuss the presence of such a massive object 60 pc away from the R136 cluster core in the context of star formation and stellar dynamics.

preprint2013arXiv

RCW36: characterizing the outcome of massive star formation

Massive stars play a dominant role in the process of clustered star formation, with their feedback into the molecular cloud through ionizing radiation, stellar winds and outflows. The formation process of massive stars is poorly constrained because of their scarcity, the short formation timescale and obscuration. By obtaining a census of the newly formed stellar population, the star formation history of the young cluster and the role of the massive stars within it can be unraveled. We aim to reconstruct the formation history of the young stellar population of the massive star-forming region RCW 36. We study several dozens of individual objects, both photometrically and spectroscopically, look for signs of multiple generations of young stars and investigate the role of the massive stars in this process. We obtain a census of the physical parameters and evolutionary status of the young stellar population. Using a combination of near-infrared photometry and spectroscopy we estimate ages and masses of individual objects. We identify the population of embedded young stellar objects (YSO) by their infrared colors and emission line spectra. RCW 36 harbors a stellar population of massive and intermediate-mass stars located around the center of the cluster. Class 0/I and II sources are found throughout the cluster. The central population has a median age of 1.1 +/- 0.6 Myr. Of the stars which could be classified, the most massive ones are situated in the center of the cluster. The central cluster is surrounded by filamentary cloud structures; within these, some embedded and accreting YSOs are found. Our age determination is consistent with the filamentary structures having been shaped by the ionizing radiation and stellar winds of the central massive stars. The formation of a new generation of stars is ongoing, as demonstrated by the presence of embedded protostellar clumps, and two exposed jets.

preprint2013arXiv

Rotational properties of the O-type star population in the Tarantula region

The 30 Doradus (30\,Dor) region in the Large Magellanic Cloud (also known as the Tarantula Nebula) is the nearest massive starburst region, containing the richest sample of massive stars in the Local Group. It is the best possible laboratory to investigate aspects of the formation and evolution of massive stars. Here, we focus on rotation which is a key parameter in the evolution of these objects. We establish the projected rotational velocity, $v_{e}\sin i$, distribution of an unprecedented sample of 216 radial velocity constant ($\rm{ΔRV\, \leq\, 20 \,km s^{-1}}$) O-type stars in 30\,Dor observed in the framework of the VLT-FLAMES Tarantula Survey (VFTS). The distribution of $v_{e}\sin i$ shows a two-component structure: a peak around 80 $\rm{km s^{-1}}$ and a high-velocity tail extending up to $\sim$600 $\rm{km s^{-1}}$. Around 75% of the sample has 0 $\leq\, v_{e}\sin i \leq$ 200 $\rm{km s^{-1}}$ with the other 25% distributed in the high-velocity tail. The presence of the low-velocity peak is consistent with that found in other studies of late-O and early-B stars. The high-velocity tail is compatible with expectations from binary interaction synthesis models and may be predominantly populated by post-binary interaction, spun-up, objects and mergers. This may have important implications for the nature of progenitors of long-duration gamma ray bursts.

preprint2013arXiv

The incidence of stellar mergers and mass gainers among massive stars

Because the majority of massive stars are born as members of close binary systems, populations of massive main-sequence stars contain stellar mergers and products of binary mass transfer. We simulate populations of massive stars accounting for all major binary evolution effects based on the most recent binary parameter statistics and extensively evaluate the effect of model uncertainties. Assuming constant star formation, we find that $8^{+9}_{-4}\,\%$ of a sample of early type stars to be the product of a merger resulting from a close binary system. In total we find that $30^{+10}_{-15}\,\%$ of massive main-sequence stars are the product of binary interaction. We show that the commonly adapted approach to minimize the effects of binaries on an observed sample by excluding systems detected as binaries through radial velocity campaigns can be counterproductive. Systems with significant radial velocity variations are mostly pre-interaction systems. Excluding them substantially enhances the relative incidence of mergers and binary products in the non radial velocity variable sample. This poses a challenge for testing single stellar evolutionary models. It also raises the question of whether certain peculiar classes of stars, such as magnetic O-stars, are the result of binary interaction and it emphasizes the need to further study the effect of binarity on the diagnostics that are used to derive the fundamental properties (star-formation history, initial mass function, mass to light ratio) of stellar populations nearby and at high redshift.

preprint2013arXiv

The VLT-FLAMES Tarantula Survey - XI. A census of the hot luminous stars and their feedback in 30 Doradus

We compile the first comprehensive census of hot luminous stars in the 30 Doradus (30 Dor) star forming region of the LMC. The census extends to a radius of 10arcmin (150pc) from the central cluster, R136. Stars were selected photometrically and combined with the latest spectral types. 1145 candidate hot luminous stars were identified of which >700 were considered genuine early type stars that contribute to feedback. We assess the spectroscopic completeness to be 85% in outer regions (>5pc) but fall to 35% in the vicinity of R136, giving a total of 500 hot luminous stars with spectroscopy. Stellar calibrations and models were used to obtain their physical parameters before integrated values were compared to global observations and the population synthesis code, Starburst99. The 31 W-R and Of/WN stars made large contributions to the total ionising and wind luminosities of ~40% and ~50%, respectively. Stars with Minit>100Msun 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 underestimated the ionising and wind luminosities of R136 by factors of ~2 and ~9, respectively. The census inferred a SFR of 0.073+/-0.04Msun/yr for 30 Dor, typically higher than results from popular SFR calibrations. However, it remained consistent with a far-UV luminosity tracer and a combined Halpha and mid-infrared tracer, but only after correcting for Halpha extinction. The global ionising output exceeded measurements from the associated gas and dust, suggesting ~6(+55/-6)% of ionising photons escape the region. When studying very luminous star forming regions, it is therefore essential to include the most massive stars to ensure a reliable energy budget. If 30 Dor is typical of other large star forming regions, estimates of the SFR will be underpredicted if this escape fraction is not accounted for.(abridged)

preprint2013arXiv

The VLT-FLAMES Tarantula Survey XII. Rotational velocities of the single O-type 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 Doradus (30 Dor). Methods. We measured projected rotational velocities, \vrot, by means of a Fourier transform method and a profile fitting method applied on a set of isolated spectral lines. We also used an iterative deconvolution procedure to infer the probability density, $\rm{P(\veq)}$, of the equatorial rotational velocity, \veq. Results. The distribution of \vrot\ shows a two-component structure: a peak around 80 \kms\ and a high-velocity tail extending up to $\sim$600 \kms. This structure is also present in the inferred distribution $\rm{P(\veq)}$ with around 80% of the sample having 0 $<$ \veq\, $\leq\, 300$ \kms\ and the other 20% distributed in the high-velocity region. 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. 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. Despite the fact that we have excluded stars that show significant radial velocity variations, our sample may have remained contaminated by post-interaction binary products. The fact that the high-velocity tail may be preferentially (and perhaps even exclusively), populated by post-binary interaction products, has important implications for the evolutionary origin of systems that produce gamma-ray bursts.

preprint2013arXiv

The VLT-FLAMES Tarantula Survey. XIII: On the nature of O Vz stars in 30 Doradus

AIMS. We test the hypothesis of O Vz stars (characterized by having HeII4686 stronger in absorption than other He lines in their blue-violet spectra) being at a younger evolutionary stage than are normal O-type dwarfs. METHODS. We have performed a quantitative spectroscopic analysis of a sample of 38 O Vz and 46 O V stars, identified by the VLT-FLAMES Tarantula Survey in the 30 Doradus region of the Large Magellanic Cloud (LMC). We obtained the stellar and wind parameters of both samples using the FASTWIND stellar atmosphere code and the IACOB-GBAT grid-based automatic tool. In the framework of a differential study, we compared the physical and evolutionary properties of both samples, regarding Teff, logg, logQ and logL. We also investigated the predictions of the FASTWIND code about the O Vz phenomenon. RESULTS. We find a differential distribution of objects in terms of effective temperature, with O Vz stars dominant at intermediate values. The O Vz stars in 30 Doradus tend to be younger and less luminous, and they have weaker winds than the O V stars, but we also find examples with ages of 2-4 Myr and with luminosities and winds that are similar to those of normal O dwarfs. Moreover, the O Vz stars do not appear to have higher gravities than the O V stars. In addition to effective temperature and wind strength, our FASTWIND predictions indicate how important it is to take other stellar parameters (gravity and projected rotational velocity) into account for correctly interpreting the O Vz phenomenon. CONCLUSIONS. In general, the O Vz stars appear to be on or very close to the ZAMS, but there are some examples where the Vz classification does not necessarily imply extreme youth. In particular, the presence of O Vz stars in our sample at more evolved phases than expected is likely a consequence of modest O-star winds owing to the low-metallicity environment of the LMC.

preprint2013arXiv

Three-dimensional orbits of the triple-O stellar system HD 150136

Context. HD 150136 is a triple hierarchical system and a non-thermal radio emitter. It is formed by an O3-3.5 V + O5.5-6 V close binary and a more distant O6.5-7 V tertiary. So far, only the inner orbital properties have been reliably constrained. Aims. To quantitatively understand the non-thermal emission process, accurate knowledge of the physical and orbital properties of the object is crucial. Here, we aim to investigate the orbital properties of the wide system and to constrain the inclinations of the inner and outer binaries, and with these the absolute masses of the system components. Methods. We used the PIONIER combiner at the Very Large Telescope Interferometer to obtain the very first interferometric measurements of HD 150136. We combined the interferometric observations with new and existing high resolution spectroscopic data to derive the orbital solution of the outer companion in the three-dimensional space. Results. The wide system is clearly resolved by PIONIER, with a projected separation on the plane of the sky of about 9 milli-arcsec. The best-fit orbital period, eccentricity, and inclination are 8.2 yr, 0.73 and 108 degr. We constrain the masses of the three stars of the system to 63 +/- 10, 40 +/- 6, and 33 +/- 12 Msun for the O3-3.5 V, O5.5-6 V and O6.5-7 V components. Conclusions. The dynamical masses agree within errors with the evolutionary masses of the components. Future interferometric and spectroscopic monitoring of HD 150136 should allow one to reduce the uncertainties to a few per cent only and to accurately constrain the distance to the system. This makes HD 150136 an ideal system to quantitatively test evolutionary models of high-mass stars as well as the physics of non-thermal processes occurring in O-type systems.

preprint2012arXiv

9 Sgr: uncovering an O-type spectroscopic binary with an 8.6 year period

The O-type object 9 Sgr is a well-known synchrotron radio emitter. This feature is usually attributed to colliding-wind binary systems, but 9 Sgr was long considered a single star. We have conducted a long-term spectroscopic monitoring of this star to investigate its multiplicity and search for evidence for wind-wind interactions. Radial velocities are determined and analysed using various period search methods. Spectral disentangling is applied to separate the spectra of the components of the binary system. We derive the first ever orbital solution of 9 Sgr. The system is found to consist of an O3.5 V((f+)) primary and an O5- 5.5 V((f)) secondary moving around each other on a highly eccentric (e = 0.7), 8.6 year orbit. The spectra reveal no variable emission lines that could be formed in the wind interaction zone in agreement with the expected properties of the interaction in such a wide system. Our results provide further support to the paradigm of synchrotron radio emission from early-type stars being a manifestation of interacting winds in a binary system.

preprint2012arXiv

Binary interaction dominates the evolution of massive stars

The presence of a nearby companion alters the evolution of massive stars in binary systems, leading to phenomena such as stellar mergers, X-ray binaries and gamma-ray bursts. Unambiguous constraints on the fraction of massive stars affected by binary interaction were lacking. We simultaneously measured all relevant binary characteristics in a sample of Galactic massive O stars and quantified the frequency and nature of binary interactions. Over seventy per cent of all massive stars will exchange mass with a companion, leading to a binary merger in one third of the cases. These numbers greatly exceed previous estimates and imply that binary interaction dominates the evolution of massive stars, with implications for populations of massive stars and their supernovae.

preprint2012arXiv

Evidence for a physically bound third component in HD 150136

Context. HD150136 is one of the nearest systems harbouring an O3 star. Although this system was for a long time considered as binary, more recent investigations have suggested the possible existence of a third component. Aims. We present a detailed analysis of HD 150136 to confirm the triple nature of this system. In addition, we investigate the physical properties of the individual components of this system. Methods. We analysed high-resolution, high signal-to-noise data collected through multi-epoch runs spread over ten years. We applied a disentangling program to refine the radial velocities and to obtain the individual spectra of each star. With the radial velocities, we computed the orbital solution of the inner system, and we describe the main properties of the orbit of the outer star such as the preliminary mass ratio, the eccentricity, and the orbital-period range. With the individual spectra, we determined the stellar parameters of each star by means of the CMFGEN atmosphere code. Results. We offer clear evidence that HD 150136 is a triple system composed of an O3V((f\ast))-3.5V((f+)), an O5.5-6V((f)), and an O6.5-7V((f)) star. The three stars are between 0-3 Myr old. We derive dynamical masses of about 64, 40, and 35 Msun for the primary, the secondary and the third components by assuming an inclination of 49°. It currently corresponds to one of the most massive systems in our galaxy. The third star moves with a period in the range of 2950 to 5500 d on an outer orbit with an eccentricity of at least 0.3. This discovery makes HD 150136 the first confirmed triple system with an O3 primary star. However, because of the long orbital period, our dataset is not sufficient to constrain the orbital solution of the tertiary component with high accuracy.

preprint2012arXiv

Multiplicity of massive O stars and evolutionary implications

Nearby companions alter the evolution of massive stars in binary systems. Using a sample of Galactic massive stars in nearby young clusters, we simultaneously measure all intrinsic binary characteristics relevant to quantify the frequency and nature of binary interactions. We find a large intrinsic binary fraction, a strong preference for short orbital periods and a flat distribution for the mass-ratios. Our results do not support the presence of a significant peak of equal-mass `twin' binaries. As a result of the measured distributions, we find that over seventy per cent of all massive stars exchange mass with a companion. Such a rate greatly exceeds previous estimates and implies that the majority of massive stars have their evolution strongly affected by interaction with a nearby companion.

preprint2012arXiv

The long period eccentric orbit of the particle accelerator HD167971 revealed by long baseline interferometry

Using optical long baseline interferometry, we resolved for the first time the two wide components of HD167971, a candidate hierarchical triple system known to efficiently accelerate particles. Our multi-epoch VLTI observations provide direct evidence for a gravitational link between the O8 supergiant and the close eclipsing O + O binary. The separation varies from 8 to 15 mas over the three-year baseline of our observations, suggesting that the components evolve on a wide and very eccentric orbit (most probably e>0.5). These results provide evidence that the wide orbit revealed by our study is not coplanar with the orbit of the inner eclipsing binary. From our measurements of the near-infrared luminosity ratio, we constrain the spectral classification of the components in the close binary to be O6-O7, and confirm that these stars are likely main-sequence objects. Our results are discussed in the context of the bright non-thermal radio emission already reported for this system, and we provide arguments in favour of a maximum radio emission coincident with periastron passage. HD167971 turns out to be an efficient O-type particle accelerator that constitutes a valuable target for future high angular resolution radio imaging using VLBI facilities.

preprint2012arXiv

The outflow history of two Herbig-Haro jets in RCW 36: HH 1042 and HH 1043

Jets around low- and intermediate-mass young stellar objects (YSOs) contain a fossil record of the recent accretion and outflow activity of their parent star-forming systems. We aim to understand whether the accretion/ejection process is similar across the entire stellar mass range of the parent YSOs. To this end we have obtained VLT/X-shooter spectra of HH 1042 and HH 1043, two newly discovered jets in the massive star-forming region RCW 36. HH 1042 is associated with the intermediate-mass YSO 08576nr292. Over 90 emission lines are detected in the spectra. High-velocity (up to 220 km/s) blue- and redshifted emission from a bipolar flow is observed in typical shock tracers. Low-velocity emission from the background cloud is detected in nebular tracers, including lines from high ionization species. We applied combined optical and infrared spectral diagnostic tools in order to derive the physical conditions (density, temperature, and ionization) in the jets. The measured mass outflow rates are Mjet ~ 10^-7 Msun/yr. We measure a high accretion rate for HH 1042 (Macc ~ 10^-6 Msun/yr) and Mjet/Macc ~ 0.1, comparable to low-mass sources and consistent with models for magneto-centrifugal jet launching. The knotted structure and velocity spread in both jets are interpreted as fossil signatures of a variable outflow rate. The mean velocities in both lobes of the jets are comparable, but the variations in Mjet and velocity in the two lobes are not symmetric, suggesting that the launching mechanism on either side of the accretion disk is not synchronized. For HH 1042, we have constructed an interpretative physical model with a stochastic or periodic outflow rate and a description of a ballistic flow as its constituents. The knotted structure and velocity spread can be reproduced qualitatively with the model, indicating that the outflow velocity varies on timescales on the order of 100 yr.

preprint2012arXiv

The rotation rates of massive stars: the role of binary interaction through tides, mass transfer and mergers

Rotation is thought to be a major factor in the evolution of massive stars, especially at low metallicity, with consequences for their chemical yields, ionizing flux and final fate. Determining the natal rotation-rate distribution of stars is of high priority given its importance as a constraint on theories of massive star formation and as input for models of stellar populations in the local Universe and at high redshift. Recently, it has become clear that the majority of massive stars interact with a binary companion before they die. We investigate how this affects the distribution of rotation rates. For this purpose, we simulate a massive binary-star population typical for our Galaxy assuming continuous star formation. We find that, because of binary interaction, 20^+5_-10% of all massive main-sequence stars have projected rotational velocities in excess of 200km/s. We evaluate the effect of uncertain input distributions and physical processes and conclude that the main uncertainties are the mass transfer efficiency and the possible effect of magnetic braking, especially if magnetic fields are generated or amplified during mass accretion and stellar mergers. The fraction of rapid rotators we derive is similar to that observed. If indeed mass transfer and mergers are the main cause for rapid rotation in massive stars, little room remains for rapidly rotating stars that are born single. This implies that spin down during star formation is even more efficient than previously thought. In addition, this raises questions about the interpretation of the surface abundances of rapidly rotating stars as evidence for rotational mixing. Furthermore, our results allow for the possibility that all early-type Be stars result from binary interactions and suggest that evidence for rotation in explosions, such as long gamma-ray bursts, points to a binary origin.

preprint2012arXiv

The VLT-FLAMES Tarantula Survey IV: Candidates for isolated high-mass star formation in 30 Doradus

Whether massive stars can occasionally form in relative isolation or if they require a large cluster of lower-mass stars around them is a key test in the differentiation of star formation theories as well as how the initial mass function of stars is sampled. Previous attempts to find O-type stars that formed in isolation were hindered by the possibility that such stars are merely runaways from clusters, i.e., their current isolation does not reflect their birth conditions. We introduce a new method to find O-type stars that are not affected by such a degeneracy. Using the VLT-FLAMES Tarantula Survey and additional high resolution imaging we have identified stars that satisfy the following constraints: 1) they are O-type stars that are not detected to be part of a binary system based on RV time series analysis; 2) they are designated spectral type O7 or earlier ; 3) their velocities are within 1σof the mean of OB-type stars in the 30 Doradus region, i.e. they are not runaways along our line-of-sight; 4) the projected surface density of stars does not increase within 3 pc towards the O-star (no evidence for clusters); 5) their sight lines are associated with gaseous and/or dusty filaments in the ISM, and 6) if a second candidate is found in the direction of the same filament with which the target is associated, both are required to have similar velocities. With these criteria, we have identified 15 stars in the 30 Doradus region, which are strong candidates for being high-mass stars that have formed in isolation. Additionally, we employed extensive MC stellar cluster simulations to confirm that our results rule out the presence of clusters around the candidates. Eleven of these are classified as Vz stars, possibly associated with the zero-age main sequence. We include a newly discovered W-R star as a candidate, although it does not meet all of the above criteria.

preprint2012arXiv

The VLT-FLAMES Tarantula Survey VI: Evidence for rotation of the young massive cluster R136

Although it has important ramifications for both the formation of star clusters and their subsequent dynamical evolution, rotation remains a largely unexplored characteristic of young star clusters (few Myr). Using multi-epoch spectroscopic data of the inner regions of 30 Doradus in the Large Magellanic Cloud (LMC) obtained as part of the VLT-FLAMES Tarantula Survey, we search for rotation of the young massive cluster R136. From the radial velocities of 36 apparently single O-type stars within a projected radius of 10 pc from the centre of the cluster, we find evidence, at the 95% confidence level, for rotation of the cluster as a whole. We use a maximum likelihood method to fit simple rotation curves to our data and find a typical rotational velocity of ~3 km/s. When compared to the low velocity dispersion of R136, our result suggests that star clusters may form with at least ~20% of the kinetic energy in rotation.

preprint2012arXiv

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

Aims. We analyze the multiplicity properties of the massive O-type star population. 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. Results. We observe a spectroscopic binary fraction of 0.35\pm0.03, which corresponds to the fraction of objects displaying statistically significant radial velocity variations with an amplitude of at least 20km/s. We compute the intrinsic binary fraction to be 0.51\pm0.04. We adopt power-laws to describe the intrinsic period and mass-ratio distributions: f_P ~ (log P)^π (with 0.15 < log P < 3.5) and f_q ~ q^κ with 0.1 < q=M_2/M_1 < 1.0. The power-law indexes that best reproduce the observed quantities are π= -0.45 +/- 0.30 and κ= -1.0\pm0.4. The obtained period distribution thus favours shorter period systems compared to an Oepik law. The mass ratio distribution is slightly skewed towards low mass ratio systems but remains incompatible with a random sampling of a classical mass function. 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. approx117 pc) and among the O9.7 I/II objects are however significantly lower than expected from statistical fluctuations. Conclusions. Using simple evolutionary considerations, we estimate that over 50% of the current O star population in 30 Dor 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.

preprint2012arXiv

The VLT-FLAMES Tarantula Survey. VII. A low velocity dispersion for the young massive cluster R136

Detailed studies of resolved young massive star clusters are necessary to determine their dynamical state and evaluate the importance of gas expulsion and early cluster evolution. In an effort to gain insight into the dynamical state of the young massive cluster R136 and obtain the first measurement of its velocity dispersion, we analyse multi-epoch spectroscopic data of the inner regions of 30 Doradus in the Large Magellanic Cloud (LMC) obtained as part of the VLT-FLAMES Tarantula Survey. Following a quantitative assessment of the variability, we use the radial velocities of non-variable sources to place an upper limit of 6 km/s on the line-of-sight velocity dispersion of stars within a projected distance of 5 pc from the centre of the cluster. After accounting for the contributions of undetected binaries and measurement errors through Monte Carlo simulations, we conclude that the true velocity dispersion is likely between 4 and 5 km/s given a range of standard assumptions about the binary distribution. This result is consistent with what is expected if the cluster is in virial equilibrium, suggesting that gas expulsion has not altered its dynamics. We find that the velocity dispersion would be ~25 km/s if binaries were not identified and rejected, confirming the importance of the multi-epoch strategy and the risk of interpreting velocity dispersion measurements of unresolved extragalactic young massive clusters.

preprint2011arXiv

On the Mass-Loss Rates of Massive Stars in the Low-Metallicity Galaxies IC 1613, WLM and NGC 3109

We present a spectroscopic analysis of VLT/X-Shooter observations of six O-type stars in the low-metallicity (Z ~ 1/7 Z\odot) galaxies IC 1613, WLM and NGC 3109. The stellar and wind parameters of these sources allow us, for the first time, to probe the mass-loss versus metallicity dependence of stellar winds below that of the Small Magellanic Cloud (at Z ~ 1/5Z\odot) by means of a modified wind momentum versus luminosity diagram. The wind strengths that we obtain for the objects in WLM and NGC 3109 are unexpectedly high and do not agree with theoretical predictions. The objects in IC 1613 tend towards a higher than expected mass-loss rate, but remain consistent with predictions within their error bars. We discuss potential systematic uncertainties in the mass-loss determinations to explain our results. However, if further scrutinization of these findings point towards an intrinsic cause for this unexpected sub-SMC mass-loss behavior, implications would include a higher than anticipated number of Wolf-Rayet stars and Ib/Ic supernovae in low-metallicity environments, but a reduced number of long-duration gamma-ray bursts produced through a single-star evolutionary channel.

preprint2011arXiv

Sparse Aperture Masking of Massive Stars

We present the earliest results of our NACO/VLT sparse aperture masking (SAM) campaign to search for binarity in a sample of 60 O-type stars. We detect Delta Ks < 5 mag companions for 20-25% of our targets with separations in the range 30-100 mas (typically, 40 - 200 A.U.). Most of these companions were unknown, shedding thus new light on the multiplicity properties of massive stars in a separation and brightness regime that has been difficult to explore so far. Adding detections from other techniques (spectroscopy, interferometry, speckle, lucky imaging, AO), the fraction of O stars with at least one companion is 85% (51/60 targets). This is the largest multiplicity fraction ever found.

preprint2011arXiv

The FLAMES Tarantula Survey

The Tarantula survey is an ESO Large Programme which has obtained multi-epochs spectroscopy of over 800 massive stars in the 30 Dor region in the Large Magelanic Cloud. Here we briefly describe the main drivers of the survey and the observational material derived.

preprint2011arXiv

The massive star binary fraction in young open clusters - III. IC 2944 and the Cen OB2 association

Using an extended set of multi-epoch high-resolution high signal-to-noise ratio optical spectra, we readdress the multiplicity properties of the O-type stars in IC 2944 and in the Cen OB2 association. We present new evidence of binarity for five objects and we confirm the multiple nature of another two. We derive the first orbital solutions for HD 100099, HD 101436 and HD 101190 and we provide additional support for HD 101205 being a quadruple system. The minimal spectroscopic binary fraction in our sample is f_min = 0.57. Using numerical simulations, we show that the detection rate of our observational campaign is close to 90 per cent, leaving thus little room for undetected spectroscopic binary systems. The statistical properties of the O-star population in IC 2944 are similar, within the uncertainties, to the results obtained in the earlier papers in this series despite the fact that sample size effects limit the significance of the comparison. Using newly derived spectroscopic parallaxes, we reassess the distance to IC 2944 and obtained 2.3+/-0.3 kpc, in agreement with previous studies. We also confirm that, as far as the O stars are concerned, the IC 2944 cluster is most likely a single entity.

preprint2011arXiv

The non-thermal radio emitter HD 93250 resolved by long baseline interferometry

As the brightest O-type X-ray source in the Carina nebula, HD 93250 (O4 III(fc)) is X-ray overluminous for its spectral type and has an unusually hard X-ray spectrum. Two different scenarios have been invoked to explain its X-ray properties: wind-wind interaction and magnetic wind confinement. Yet, HD 93250 shows absolutely constant radial velocities over time scales of years suggesting either a single star, a binary system seen pole-one or a very long period and/or highly eccentric system. Using the ESO Very Large Telescope Interferometer, we resolved HD 93250 as a close pair with similar components. We measured a near-infrared flux ratio of 0.8+/-0.1 and a separation of 1.5+/-0.2 x 10E-03 arcsec. At the distance of Carina, this corresponds to a projected physical distance of 3.5 A.U. While a quantitative investigation would require a full characterization of the orbit, the binary nature of HD 93250 allows us to qualitatively explain both its X-ray flux and hardness and its non-thermal radio emission in the framework of a colliding wind scenario. We also discuss various observational biases. We show that, due to line-blending of two similar spectral components, HD 93250 could have a period as short as one to several years despite the lack of measurable radial velocity variations.

preprint2011arXiv

The VLT-FLAMES Tarantula Survey I: Introduction and observational overview

The VLT-FLAMES Tarantula Survey (VFTS) is an ESO Large Programme that has obtained multi-epoch optical spectroscopy of over 800 massive stars in the 30 Doradus region of the Large Magellanic Cloud (LMC). Here we introduce our scientific motivations and give an overview of the survey targets, including optical and near-infrared photometry and comprehensive details of the data reduction. One of the principal objectives was to detect massive binary systems via variations in their radial velocities, thus shaping the multi-epoch observing strategy. Spectral classifications are given for the massive emission-line stars observed by the survey, including the discovery of a new Wolf-Rayet star (VFTS 682, classified as WN5h), 2' to the northeast of R136. To illustrate the diversity of objects encompassed by the survey, we investigate the spectral properties of sixteen targets identified by Gruendl & Chu from Spitzer photometry as candidate young stellar objects or stars with notable mid-infrared excesses. Detailed spectral classification and quantitative analysis of the O- and B-type stars in the VFTS sample, paying particular attention to the effects of rotational mixing and binarity, will be presented in a series of future articles to address fundamental questions in both stellar and cluster evolution.

preprint2011arXiv

The VLT-FLAMES Tarantula Survey II: R139 revealed as a massive binary system

We report the discovery that R139 in 30 Doradus is a massive spectroscopic binary system. Multi-epoch optical spectroscopy of R139 was obtained as part of the VLT-FLAMES Tarantula Survey, revealing a double-lined system. The two components are of similar spectral types; the primary exhibits strong C III 4650 emission and is classified as an O6.5 Iafc supergiant, while the secondary is an O6 Iaf supergiant. The radial-velocity variations indicate a highly eccentric orbit with a period of 153.9 days. Photometry obtained with the Faulkes Telescope South shows no evidence for significant variability within an 18 month period. The orbital solution yields lower mass limits for the components of M1sin^3 i = 78 \pm 8 Msun and M2sin^3 i = 66 \pm 7 Msun. As R139 appears to be the most massive binary system known to contain two evolved Of supergiants, it will provide an excellent test for atmospheric and evolutionary models.

preprint2011arXiv

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/Lsun) = 6.5 \pm 0.2, corresponding to a present-day mass of ~150Msun. The high effective temperature of 52.2 \pm 2.5kK 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 Wolf-Rayet 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.

preprint2011arXiv

The VLT-FLAMES Tarantula Survey: The fastest rotating O-type star and shortest period LMC pulsar - remnants of a supernova disrupted binary?

We present a spectroscopic analysis of an extremely rapidly rotating late O-type star, VFTS102, observed during a spectroscopic survey of 30 Doradus. VFTS102 has a projected rotational velocity larger than 500\kms\ and probably as large as 600\kms; as such it would appear to be the most rapidly rotating massive star currently identified. Its radial velocity differs by 40\kms\ from the mean for 30 Doradus, suggesting that it is a runaway. VFTS102 lies 12 pcs from the X-ray pulsar PSR J0537-6910 in the tail of its X-ray diffuse emission. We suggest that these objects originated from a binary system with the rotational and radial velocities of VFTS102 resulting from mass transfer from the progenitor of PSR J0537-691 and the supernova explosion respectively.

preprint2010arXiv

A MAD view of Trumpler 14

We present adaptive optics (AO) near-infrared observations of the core of the Tr 14 cluster in the Carina region obtained with the ESO multi-conjugate AO demonstrator, MAD. Our campaign yields AO-corrected observations with an image quality of about 0.2 arcsec across the 2 arcmin field of view, which is the widest AO mosaic ever obtained. We detected almost 2000 sources spanning a dynamic range of 10 mag. The pre-main sequence (PMS) locus in the colour-magnitude diagram is well reproduced by Palla & Stahler isochrones with an age of 3 to 5 1E+05 yr, confirming the very young age of the cluster. We derive a very high (deprojected) central density n0~4.5(+/-0.5) \times 10^4 pc^-3 and estimate the total mass of the cluster to be about ~4.3^{+3.3}_{-1.5} \times 10^3 Msun, although contamination of the field of view might have a significant impact on the derived mass. We show that the pairing process is largely dominated by chance alignment so that physical pairs are difficult to disentangle from spurious ones based on our single epoch observation. Yet, we identify 150 likely bound pairs, 30% of these with a separation smaller than 0.5 arcsec (~1300AU). We further show that at the 2-sigma level massive stars have more companions than lower-mass stars and that those companions are respectively brighter on average, thus more massive. Finally, we find some hints of mass segregation for stars heavier than about 10 Msun. If confirmed, the observed degree of mass segregation could be explained by dynamical evolution, despite the young age of the cluster.

preprint2010arXiv

A massive runaway star from 30 Doradus

We present the first ultraviolet (UV) and multi-epoch optical spectroscopy of 30 Dor 016, a massive O2-type star on the periphery of 30 Doradus in the Large Magellanic Cloud. The UV data were obtained with the Cosmic Origins Spectrograph on the Hubble Space Telescope as part of the Servicing Mission Observatory Verification program after Servicing Mission 4, and reveal #016 to have one of the fastest stellar winds known. From analysis of the CIV 1548-51 doublet we find a terminal velocity, v_infty=3450 +/- 50km/s. Optical spectroscopy is from the VLT-FLAMES Tarantula Survey, from which we rule out a massive companion (with 2d<P<1yr) to a confidence of 98%. The radial velocity of #016 is offset from the systemic value by -85km/s, suggesting that the star has traveled the 120pc from the core of 30 Doradus as a runaway, ejected via dynamical interactions.

preprint2010arXiv

Chemical evolution of the Galactic bulge as traced by microlensed dwarf and subgiant stars. II. Ages, metallicities, detailed elemental abundances, and connections to the Galactic thick disc

The Bulge is the least understood major stellar population of the Milky Way. Most of what we know about the formation and evolution of the Bulge comes from bright giant stars. The underlying assumption that giants represent all the stars, and accurately trace the chemical evolution of a stellar population, is under debate. In particular, recent observations of a few microlensed dwarf stars give a very different picture of the evolution of the Bulge from that given by the giant stars. [ABRIDGED] We perform a detailed elemental abundance analysis of dwarf stars in the Galactic bulge, based on high-resolution spectra that were obtained while the stars were optically magnified during gravitational microlensing events. [ABRIDGED] We present detailed elemental abundances and stellar ages for six new dwarf stars in the Galactic bulge. Combining these with previous events, here re-analysed with the same methods, we study a homogeneous sample of 15 stars, which constitute the largest sample to date of microlensed dwarf stars in the Galactic bulge. We find that the stars span the full range of metallicities from [Fe/H]=-0.72 to +0.54, and an average metallicity of <[Fe/H]>=-0.08+/-0.47, close to the average metallicity based on giant stars in the Bulge. Furthermore, the stars follow well-defined abundance trends, that for [Fe/H]<0 are very similar to those of the local Galactic thick disc. This suggests that the Bulge and the thick disc have had, at least partially, comparable chemical histories. At sub-solar metallicities we find the Bulge dwarf stars to have consistently old ages, while at super-solar metallicities we find a wide range of ages. Using the new age and abundance results from the microlensed dwarf stars we investigate possible formation scenarios for the Bulge.

preprint2010arXiv

Optical atmospheric extinction over Cerro Paranal

Aims: The present study was conducted to determine the optical extinction curve for Cerro Paranal under typical clear-sky observing conditions, with an accuracy of 0.01 mag/airmass. Methods: The extinction curve of Paranal was obtained through low-resolution spectroscopy of 8 spectrophotometric standard stars observed with FORS1 mounted at the 8.2 m Very Large Telescope, covering a spectral range 3300-8000 A. A total of 600 spectra were collected on more than 40 nights distributed over six months, from October 2008 to March 2009. The average extinction curve was derived using a global fit algorithm, which allowed us to simultaneously combine all the available data. The main atmospheric parameters were retrieved using the LBLRTM radiative transfer code, which was also utilised to study the impact of variability of the main molecular bands of O2, O3, and H2O, and to estimate their column densities. Results: In general, the extinction curve of Paranal appears to conform to those derived for other astronomical sites in the Atacama desert, like La Silla and Cerro Tololo. However, a systematic deficit with respect to the extinction curve derived for Cerro Tololo before the El Chichon eruption is detected below 4000 A. We attribute this downturn to a non standard aerosol composition, probably revealing the presence of volcanic pollutants above the Atacama desert. An analysis of all spectroscopic extinction curves obtained since 1974 shows that the aerosol composition has been evolving during the last 35 years. The persistence of traces of non meteorologic haze suggests the effect of volcanic eruptions, like those of El Chichon and Pinatubo, lasts several decades. The usage of the standard CTIO and La Silla extinction curves implemented in IRAF and MIDAS produce systematic over/under-estimates of the absolute flux.

preprint2010arXiv

The multiplicity of massive stars

Binaries are excellent astrophysical laboratories that provide us with direct measurements of fundamental stellar parameters. Compared to single isolated star, multiplicity induces new processes, offering the opportunity to confront our understanding of a broad range of physics under the extreme conditions found in, and close to, astrophysical objects. In this contribution, we will discuss the parameter space occupied by massive binaries, and the observational means to investigate it. We will review the multiplicity fraction of OB stars within each regime, and in different astrophysical environments. In particular we will compare the O star spectroscopic binary fraction in nearby open clusters and we will show that the current data are adequately described by an homogeneous fraction of f~0.44. We will also summarize our current understanding of the observed parameter distributions of O+OB spectroscopic binaries. We will show that the period distribution is overabundant in short period binaries and that it can be described by a bi-modal Oepik law with a break point around P~10d. The distribution of the mass-ratios shows no indication for a twin population of equal mass binaries and seems rather uniform in the range 0.2< q=M_2/M_1<1.0.

preprint2009arXiv

A new investigation of the binary HD 48099

With an orbital period of about 3.078 days, the double-lined spectroscopic binary HD 48099 is, still now, the only short-period O+O system known in the Mon OB2 association. Even though an orbital solution has already been derived for this system, few information are available about the individual stars. We present, in this paper, the results of a long-term spectroscopic campaign. We derive a new orbital solution and apply a disentangling method to recover the mean spectrum of each star. To improve our knowledge concerning both components, we determine their spectral classifications and their projected rotational velocities. We also constrain the main stellar parameters of both stars by using the CMFGEN atmosphere code and provide the wind properties for the primary star through the study of IUE spectra. This investigation reveals that HD 48099 is an O5.5 V((f))+O9 V binary with M_1 sin^3 i = 0.70 M_{\sun} and M_2 sin^3 i = 0.39 M_{\sun}, implying a rather low orbital inclination. This result, combined with both a large effective temperature and log g, suggests that the primary star (v sini ~ 91 km s^-1) is actually a fast rotator with a strongly clumped wind and a nitrogen abundance of about 8 times the solar value.

preprint2009arXiv

Elemental abundances in the Galactic bulge from microlensed dwarf stars

We present elemental abundances of 13 microlensed dwarf and subgiant stars in the Galactic bulge, which constitute the largest sample to date. We show that these stars span the full range of metallicity from Fe/H=-0.8 to +0.4, and that they follow well-defined abundance trends, coincident with those of the Galactic thick disc.

preprint2009arXiv

On the velocity dispersion of young star clusters: super-virial or binaries?

Many young extra-galactic clusters have a measured velocity dispersion that is too high for the mass derived from their age and total luminosity, which has led to the suggestion that they are not in virial equilibrium. Most of these clusters are confined to a narrow age range centred around 10 Myr because of observational constraints. At this age the cluster light is dominated by luminous evolved stars, such as red supergiants, with initial masses of ~13-22 Msun for which (primordial) binarity is high. In this study we investigate to what extent the observed excess velocity dispersion is the result of the orbital motions of binaries. We demonstrate that estimates for the dynamical mass of young star clusters, derived from the observed velocity dispersion, exceed the photometric mass by up-to a factor of 10 and are consistent with a constant offset in the square of the velocity dispersion. This can be reproduced by models of virialised star clusters hosting a massive star population of which ~25 is in binaries, with typical mass ratios of ~0.6 and periods of ~1000 days. We conclude that binaries play a pivotal role in deriving the dynamical masses of young (~10 Myr) moderately massive and compact (<1e5 Msun; > 1 pc) star clusters.

preprint2009arXiv

The massive star binary fraction in young open clusters - II. NGC 6611 (Eagle Nebula)

Based on a set of over 100 medium- to high-resolution optical spectra collected from 2003 to 2009, we investigate the properties of the O-type star population in NGC6611 in the core of the Eagle Nebula (M16). Using a much more extended data set than previously available, we revise the spectral classification and multiplicity status of the nine O-type stars in our sample. We confirm two suspected binaries and derive the first SB2 orbital solutions for two systems. We further report that two other objects are displaying a composite spectrum, suggesting possible long-period binaries. Our analysis is supported by a set of Monte-Carlo simulations, allowing us to estimate the detection biases of our campaign and showing that the latter do not affect our conclusions. The absolute minimal binary fraction in our sample is f_min=0.44 but could be as high as 0.67 if all the binary candidates are confirmed. As in NGC6231 (see Paper I), up to 75% of the O star population in NGC6611 are found in an O+OB system, thus implicitly excluding random pairing from a classical IMF as a process to describe the companion association in massive binaries. No statistical difference could be further identified in the binary fraction, mass-ratio and period distributions between NGC6231 and NGC6611, despite the difference in age and environment of the two clusters.

preprint2009arXiv

The nature of the line profile variability in the spectrum of the massive binary HD 152219

HD 152219 is a massive binary system with O9.5 III + B1-2 V/III components and a short orbital period of 4.2 d. Its primary component further displays clear line profile variability (LPV). The primary component being located within the pulsational instability domain predicted for high-luminosity stars, we previously suggested that the observed LPV could be associated with non-radial pulsations. The aim of the present work is to determine the nature of the observed LPV in the spectrum of the primary component of HD 152219.During a 4-night FEROS monitoring campaign, we collected a new set of 134 high signal-to-noise spectra. These new observations were then used to re-investigate the variability of different line profiles in the spectrum of HD 152219. Based on the present analysis, we discard the non-radial pulsations and point out the Rossiter-McLaughlin effect as the cause of the LPV in HD 152219. The upper limit on the amplitude of possible weak pulsations is set at a few parts per thousand of the continuum level.

preprint2009arXiv

The VLT-FLAMES Tarantula Survey

The Tarantula Survey is an ambitious ESO Large Programme that has obtained multi-epoch spectroscopy of over 1,000 massive stars in the 30 Doradus region of the Large Magellanic Cloud. Here we introduce the scientific motivations of the survey and give an overview of the observational sample. Ultimately, quantitative analysis of every star, paying particular attention to the effects of rotational mixing and binarity, will be used to address fundamental questions in both stellar and cluster evolution.

preprint2007arXiv

Constraining the fundamental parameters of the O-type binary CPD-41degr7733

Using a set of high-resolution spectra, we studied the physical and orbital properties of the O-type binary CPD-41 7733, located in the core of \ngc. We report the unambiguous detection of the secondary spectral signature and we derive the first SB2 orbital solution of the system. The period is 5.6815 +/- 0.0015 d and the orbit has no significant eccentricity. CPD-41 7733 probably consists of stars of spectral types O8.5 and B3. As for other objects in the cluster, we observe discrepant luminosity classifications while using spectroscopic or brightness criteria. Still, the present analysis suggests that both components display physical parameters close to those of typical O8.5 and B3 dwarfs. We also analyze the X-ray light curves and spectra obtained during six 30 ks XMM-Newton pointings spread over the 5.7 d period. We find no significant variability between the different pointings, nor within the individual observations. The CPD-41 7733 X-ray spectrum is well reproduced by a three-temperature thermal mekal model with temperatures of 0.3, 0.8 and 2.4 keV. No X-ray overluminosity, resulting e.g. from a possible wind interaction, is observed. The emission of CPD-41 7733 is thus very representative of typical O-type star X-ray emission.

H. Sana

What is connected

Connect this record

See the researcher in context

Building this map preview

106 published item(s)

A spectroscopic multiplicity survey of Galactic Wolf-Rayet stars: II. The northern WNE sequence

Detecting Stripped Stars While Searching for Quiescent Black Holes

Identifying quiescent compact objects in massive Galactic single-lined spectroscopic binaries

Modeling overcontact binaries, I. The effect of tidal deformation

On the origin of close massive binaries in the M17 star-forming region

Properties of the Be-type stars in 30 Doradus

A relation between the radial velocity dispersion of young clusters and their age: Evidence for hardening as the formation scenario of massive close binaries

Analytic, dust-independent mass-loss rates for red supergiant winds initiated by turbulent pressure

Binary-object spectral-synthesis in 3D (BOSS-3D) -- Modelling H-alpha emission in the enigmatic multiple system LB-1

Detailed models of interacting short-period massive binary stars

Probing the low-mass end of the companion mass function for O-type stars

The CubeSpec space mission. I. Asteroseismology of massive stars from time series optical spectroscopy: science requirements and target list prioritisation

The inner circumstellar dust of the red supergiant Antares as seen with VLT/SPHERE/ZIMPOL

The observed multiplicity properties of B-type stars in the Galactic young open cluster NGC 6231

Uncovering astrometric black hole binaries with massive main-sequence companions with Gaia

Investigating the lack of main-sequence companions to massive Be stars

Is HR 6819 a triple system containing a black hole? -- An alternative explanation

Multiplicity of the red supergiant population in the young massive cluster NGC 330

Properties of OB star-black hole systems derived from detailed binary evolution models

The 'hidden' companion in LB-1 unveiled by spectral disentangling

The high-energy emission from HD~93129A near periastron

The VLT-FLAMES Tarantula Survey XXXII. Low-luminosity late O-type stars -- classification, main physical parameters, and silicon abundances

The Wolf-Rayet binaries of the nitrogen sequence in the Large Magellanic Cloud: spectroscopy, orbital analysis, formation, and evolution

VLT/X-shooter spectroscopy of massive young stellar objects in the 30 Doradus region of the Large Magellanic Cloud

Why binary interaction does not necessarily dominate the formation of Wolf-Rayet stars at low metallicity

The Tarantula Massive Binary Monitoring: III. Atmosphere analysis of double-lined spectroscopic systems

The Tarantula Massive Binary Monitoring. IV. Double-lined photometric binaries

The young massive SMC cluster NGC 330 seen by MUSE. I. Observations and stellar content

The young stellar content of the giant HII regions M8, G333.6-0.2, and NGC6357 with VLT/KMOS

Response to comment on "An excess of massive stars in the local 30 Doradus starburst"

A new prescription for the mass-loss rates of WC and WO stars

Evidence of magnetic field decay in massive main-sequence stars

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

B fields in OB stars (BOB): Detection of a strong magnetic field in the O9.7 V star HD54879

B fields in OB stars (BOB): low-resolution FORS2 spectropolarimetry of the first sample of 50 massive stars

Discovery of the massive overcontact binary VFTS 352: Evidence for enhanced internal mixing

HST astrometry in the 30 Doradus region: measuring proper motions of individual stars in the Large Magellanic Cloud

Hubble Tarantula Treasury Project. III. Photometric Catalog and Resulting Constraints on the Progression of Star Formation in the 30 Doradus Region

Masses of the components of SB2 binaries observed with Gaia. II. Masses derived from PIONIER interferometric observations for Gaia validation

Massive stars on the verge of exploding: the properties of oxygen sequence Wolf-Rayet stars

Molecfit: A general tool for telluric absorption correction II. Quantitative evaluation on ESO-VLT X-Shooter spectra

Molecfit: A general tool for telluric absorption correction. I. Method and application to ESO instruments

The evolution of rotating very massive stars with LMC composition

The mass of the very massive binary WR21a

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

The VLT-FLAMES Tarantula Survey XVIII. Classifications and radial velocities of the B-type stars

The VLT-FLAMES Tarantula survey XX. The nature of the X-ray bright emission line star VFTS 399

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

The VLT-FLAMES Tarantula Survey XXII. Multiplicity properties of the B-type stars

The VLT-FLAMES Tarantula Survey. XXIII: two massive double-lined binaries in 30~Doradus

B fields in OB stars (BOB): on the detection of weak magnetic fields in the two early B-type stars beta CMa and epsilon CMa

HD 152246 - a new high-mass triple system and its basic properties

Molecfit: A Package for Telluric Absorption Correction

Relating jet structure to photometric variability: the Herbig Ae star HD 163296

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

Southern Massive Stars at High Angular Resolution: Observational Campaign and Companion Detection

The B Fields in OB Stars (BOB) Survey

The BinaMIcS project: understanding the origin of magnetic fields in massive stars through close binary systems

The properties of single WO stars

The properties of ten O-type stars in the low-metallicity galaxies IC 1613, WLM and NGC 3109

The Three body system delta Circini

The VLT-FLAMES Tarantula Survey XIV. The O-Type Stellar Content of 30 Doradus

The VLT-FLAMES Tarantula Survey XVI. The optical+NIR extinction laws in 30 Doradus and the photometric determination of the effective temperatures of OB stars

The VLT-FLAMES Tarantula Survey XVII. Physical and wind properties of massive stars at the top of the main sequence

On the nature of the WO3 star DR1 in IC 1613

On the nature of WO stars: a quantitative analysis of the WO3 star DR1 in IC 1613

R144 revealed as a double-lined spectroscopic binary

RCW36: characterizing the outcome of massive star formation

Rotational properties of the O-type star population in the Tarantula region

The incidence of stellar mergers and mass gainers among massive stars

The VLT-FLAMES Tarantula Survey - XI. A census of the hot luminous stars and their feedback in 30 Doradus

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

The VLT-FLAMES Tarantula Survey. XIII: On the nature of O Vz stars in 30 Doradus

Three-dimensional orbits of the triple-O stellar system HD 150136