Source author record

N. Langer

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

119works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

119 published item(s)

preprint2025arXiv

Populations of evolved massive binary stars in the Small Magellanic Cloud I: Predictions from detailed evolution models

Context. The majority of massive stars are born with a close binary companion. How this affects their evolution and fate is still largely uncertain, especially at low metallicity. Aims. We derive synthetic populations of massive post-interaction binary products and compare them with corresponding observed populations in the Small Magellanic Cloud (SMC). Methods. We analyse 53298 detailed binary evolutionary models computed with MESA. Our models include the physics of rotation, mass and angular momentum transfer, magnetic internal angular momentum transport, and tidal spin-orbit coupling. They cover initial primary masses of 5-100Msun, initial mass ratios of 0.3-0.95, and all initial periods for which interaction is expected. They are evolved through the first mass transfer and the donor star death, a possible ensuing Be/X-ray binary phase, and they end when the mass gainer leaves the main sequence. Results.In our fiducial synthetic population, 8% of the OB stars in the SMC are post-mass transfer systems, and 7% are merger products. In many of our models, the mass gainers are spun up and form Oe/Be stars. While our model underpredicts the number of Be/X-ray binaries in the SMC, it reproduces the main features of their orbital period distribution and the observed number of SMC binary WR stars. We expect $\sim$50 OB+BH binaries below and $\sim$170 above 20d orbital period. The latter might produce merging double BHs. However, their progenitors, the predicted long-period WR+OB binaries, are not observed. Conclusions. While the comparison with the observed SMC stars supports many physics assumptions in our high-mass binary models, a better match of the large number of observed OBe stars and Be/X-ray binaries likely requires a lower merger rate and/or a higher mass transfer efficiency during the first mass transfer. The fate of the initially wide O star binaries remains uncertain.

preprint2025arXiv

Populations of evolved massive binary stars in the Small Magellanic Cloud II: Predictions from rapid binary evolution

Massive star evolution plays a crucial role in astrophysics but bares large uncertainties. This problem becomes more severe by the majority of massive stars being born in close binary systems, whose evolution is affected by the interaction of their components. We want to constrain major uncertainties in massive binary star evolution, in particular the efficiency and the stability of the first mass transfer phase. We use the rapid population synthesis code ComBinE to generate synthetic populations of post-interaction binaries, assuming constant mass-transfer efficiency. We employ a new merger criterion that adjusts self-consistently to any prescribed mass-transfer efficiency. We tailor our synthetic populations to be comparable to the expected binary populations in the Small Magellanic Cloud (SMC). We find that the observed populations of evolved massive binaries can not be reproduced with a single mass-transfer efficiency. Instead, a rather high efficiency (>50%) is needed to reproduce the number of Be stars and Be/X-ray binaries in the SMC, while a low efficiency (~10%) leads to a better agreement with the observed number of Wolf-Rayet stars. We construct a corresponding mass-dependent mass-transfer efficiency recipe to produce our fiducial synthetic SMC post-interaction binary population. It reproduces the observed number and properties of the Be/X-ray and WR-binaries rather well, and is not in stark disagreement with the observed OBe star population. It further predicts two large, yet unobserved populations of OB+BH binaries, that is ~100 OB+BH systems with rather small orbital periods (<20 days) and ~40 longer period OBe+BH systems.

preprint2022arXiv

Hunting for red supergiant binaries: UVIT photometry of the SMC

We present UVIT/Astrosat UV photometry of the RSG population of the Small Cloud galaxy (SMC). As RSGs are extremely faint in the far-UV, these observations directly probe potential companion stars. From a sample of 861 SMC RSGs, we find 88 have detections at far-UV wavelengths: a clear signature of binarity. Stellar parameters are determined for both components, which allows us to study - for the first time - the mass-ratio (q) distribution of RSG binary systems. We find a flat mass-ratio distribution best describes the observations up to M{RSG}~15 Msun. We account for our main observing bias (i.e. the limiting magnitude of the UVIT survey) to determine the intrinsic RSG binary fraction of 18.8+/-1.5 %, for mass-ratios in the range 0.3 < q < 1.0 and orbital periods approximately in the range 3 < log P[ days] < 8.

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

Red supergiant stars in binary systems. I. Identification and characterisation in the Small Magellanic Cloud from the UVIT ultraviolet imaging survey

We aim to identify and characterise binary systems containing red supergiant (RSG) stars in the Small Magellanic Cloud (SMC) using a newly available ultraviolet (UV) point source catalogue obtained using the Ultraviolet Imaging Telescope (UVIT) on board AstroSat. We select a sample of 560 SMC RSGs based on photometric and spectroscopic observations at optical wavelengths and cross-match this with the far-UV point source catalogue using the UVIT F172M filter, finding 88 matches down to m$_{F172M}$=20.3 ABmag, which we interpret as hot companions to the RSGs. Stellar parameters (luminosities, effective temperatures and masses) for both components in all 88 binary systems are determined and we find mass distributions in the ranges 6.1 to 22.3 Solar masses for RSGs and 3.7 to 15.6 Solar masses for their companions. The most massive RSG binary system in the SMC has a combined mass of 32 $\pm$4 M$_\odot$, with a mass ratio (q) of 0.92. By simulating observing biases, we find an intrinsic multipliciy fraction of 18.8 $\pm$ 1.5% for mass ratios in the range 0.3 < q < 1.0 and orbital periods approximately in the range 3 < log P [days] < 8. By comparing our results with those of a similar mass on the main-sequence, we determine the fraction of single stars to be ~20% and argue that the orbital period distribution declines rapidly beyond log P ~ 3.5. We study the mass-ratio distribution of RSG binary systems and find that a uniform distribution best describes the data below 14 M$_\odot$. Above 15 M$_\odot$, we find a lack of high mass-ratio systems.

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

The VLT-FLAMES survey of massive stars. NGC2004#115: A triple system hosting a possible short period B+BH binary

The star NGC2004#115 in the LMC, originally classified as an (SB1) Be spectroscopic binary, bears some morphological resemblance to the Galactic systems LB-1 and HR 6819, both of which are proposed as either Be+black hole (BH) or Be+stripped He-star systems. Two data-sets (ESO/VLT and SALT) of multi-epoch optical spectra of NGC 2004#115, separated by baseline of $\sim$20 years, lead us conclude it is a triple system hosting an inner binary with a period of 2.92 d, eccentricity $\sim$0.0 and mass function $\sim$0.07 $M_\odot$. The inner binary harbours a B-type star (the primary) with projected rotational velocity of 10km/s, and luminosity $\log L/L_\odot$=3.87, contributing $\sim$60% of the V-band light to the system. The secondary is not detected, while the tertiary, which contributes 40% of the light, is tentatively identified as a less luminous B-type star with high projected rotational velocity. No ellipsoidal light variability is detected, with stringent limits being set by MACHO and Gaia data. Assuming the primary to be a main sequence star yields a mass of 8.6$ M_\odot$, while the additional assumption of synchronous rotation constrains the inclination to be almost pole-on with i~9 degrees, implying the secondary is a BH with a mass of $\sim$25 $M_\odot$. A low mass stripped star with similar luminosity is ruled out as a potential solution as its mass implies a Roche radius that is substantially smaller than the stellar radius. The outer period likely exceeds 120 days and, while the disk-like emission is variable (it is almost absent in the SALT dataset), it may be associated with the inner binary rather than the rapidly rotating tertiary. XMM-Newton provides an upper limit of 5x$10^{33}$ ergs/s on the X-ray flux, consistent with, though not constraining of, the system hosting a quiescent B+BH binary. A number of caveats to this scenario are discussed in the paper.

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

A dearth of young and bright massive stars in the Small Magellanic Cloud

Massive star evolution at low metallicity is closely connected to many fields in high-redshift astrophysics, but poorly understood. The Small Magellanic Cloud (SMC) is a unique laboratory to study this because of its metallicity of 0.2 Zsol, its proximity, and because it is currently forming stars. We used a spectral type catalog in combination with GAIA magnitudes to calculate temperatures and luminosities of bright SMC stars. By comparing these with literature studies, we tested the validity of our method, and using GAIA data, we estimated the completeness of stars in the catalog as a function of luminosity. This allowed us to obtain a nearly complete view of the most luminous stars in the SMC. When then compared with stellar evolution predictions. We also calculated the extinction distribution, the ionizing photon production rate, and the star formation rate. Our results imply that the SMS hosts only 30 very luminous main-sequence stars (M > 40 Msol; L > 10^5 Lsol), which are far fewer than expected from the number of stars in the luminosity range 3*10^4 < L/Lsol < 3*10^5 and from the typically quoted star formation rate in the SMC. Even more striking, we find that for masses above M > 20 Msol, stars in the first half of their hydrogen-burning phase are almost absent. This mirrors a qualitatively similar peculiarity that is known for the Milky Way and Large Magellanic Cloud. This amounts to a lack of hydrogen-burning counterparts of helium-burning stars, which is more pronounced for higher luminosities. We argue that a declining star formation rate or a steep initial mass function are unlikely to be the sole explanations for the dearth of young bright stars. Instead, many of these stars might be embedded in their birth clouds, although observational evidence for this is weak. We discuss implications for cosmic reionization and the top end of the initial mass function.

preprint2020arXiv

Galactic interstellar sulfur isotopes: A radial $^{32}$S$/$$^{34}$S gradient?

We present observations of $^{12}$C$^{32}$S, $^{12}$C$^{34}$S, $^{13}$C$^{32}$S and $^{12}$C$^{33}$S J=2$-$1 lines toward a large sample of massive star forming regions by using the Arizona Radio Observatory 12-m telescope and the IRAM\,30-m. Taking new measurements of the carbon $^{12}$C/$^{13}$C ratio, the $^{32}$S$/$$^{34}$S isotope ratio was determined from the integrated $^{13}$C$^{32}$S/$^{12}$C$^{34}$S line intensity ratios for our sample. Our analysis shows a $^{32}$S$/$$^{34}$S gradient from the inner Galaxy out to a galactocentric distance of 12\,kpc. An unweighted least-squares fit to our data yields $^{32}$S$/$$^{34}$S = (1.56 $\pm$ 0.17)$\rm D_{\rm GC}$ + (6.75 $\pm$ 1.22) with a correlation coefficient of 0.77. Errors represent 1$σ$ standard deviations. Testing this result by (a) excluding the Galactic center region, (b) excluding all sources with C$^{34}$S opacities $>$ 0.25, (c) combining our data and old data from previous study, and (d) using different sets of carbon isotope ratios leads to the conclusion that the observed $^{32}$S$/$$^{34}$S gradient is not an artefact but persists irrespective of the choice of the sample and carbon isotope data. A gradient with rising $^{32}$S$/$$^{34}$S values as a function of galactocentric radius implies that the solar system ratio should be larger than that of the local interstellar medium. With the new carbon isotope ratios we obtain indeed a local $^{32}$S$/$$^{34}$S isotope ratio about 10$\%$ below the solar system one, as expected in case of decreasing $^{32}$S$/$$^{34}$S ratios with time and increased amounts of stellar processing. However, taking older carbon isotope ratios based on a lesser amount of data, such a decrease is not seen. No systematic variation of $^{34}$S$/$$^{33}$S ratios along galactocentric distance was found.

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.

preprint2019arXiv

gamma Cas stars: Normal Be stars with disks impacted by the wind of a helium-star companion?

$γ$ Cas stars are a $\sim$1% minority among classical Be stars with hard but only moderately strong continuous thermal X-ray flux and mostly very early-B spectral type. The X-ray flux has been suggested to originate from matter accelerated via magnetic disk-star interaction, by a rapidly rotating neutron star (NS) companion via the propeller effect, or by accretion onto a white dwarf (WD) companion. In view of the growing number of identified $γ$ Cas stars and the only imperfect matches between these suggestions and the observations, alternative models should be pursued. Two of the three best-observed $γ$ Cas stars, $γ$ Cas itself and $π$ Aqr, have a low-mass companion with low optical flux; interferometry of BZ Cru is inconclusive. Binary-evolution models are examined for their ability to produce such systems. The OB+He-star stage of post-mass transfer binaries, which is otherwise observationally unaccounted, can potentially reproduce many observed properties of $γ$ Cas stars. The interaction of the fast wind of helium stars with the disk and/or with the wind of Be stars may give rise to the production of hard X-rays. While not modelling this process, it is shown that the energy budget is favourable, and that the wind velocities may lead to hard X-rays as observed in $γ$ Cas stars. Furthermore, their observed number appears to be consistent with the evolutionary models. Within the Be+He-star binary model, the Be stars in $γ$ Cas stars are conventional classical Be stars. They are encompassed by O-star+Wolf-Rayet systems towards higher mass, where no stable Be decretion disks exist, and by Be+sdO systems at lower mass where the sdO winds may be too weak to cause the $γ$ Cas phenomenon. In decreasing order of the helium-star mass, the descendants could be Be+black-hole, Be+NS or Be+WD binaries.

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

WR 72: a born-again planetary nebula with hydrogen-poor knots

We report the discovery of a handful of optical hydrogen-poor knots in the central part of an extended infrared nebula centred on the [WO1] star WR 72, obtained by spectroscopic and imaging observations with the Southern African Large Telescope (SALT). Wide-field Infrared Survey Explorer (WISE) images of the nebula show that it is composed of an extended almost circular halo (of $\approx6$ arcmin or $\approx2.4$ pc in diameter) and an elongated and apparently bipolar inner shell (of a factor of six smaller size), within which the knots are concentrated. Our findings indicate that WR 72 is a new member of the rare group of hydrogen-poor planetary nebulae, which may be explained through a very late thermal pulse of a post-AGB star, or by a merger of two white dwarfs.

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

Asteroseismic test of rotational mixing in low-mass white dwarfs

We exploit the recent discovery of pulsations in mixed-atmosphere (He/H), extremely low-mass white dwarf precursors (ELM proto-WDs) to test the proposition that rotational mixing is a fundamental process in the formation and evolution of low-mass helium core white dwarfs. Rotational mixing has been shown to be a mechanism able to compete efficiently against gravitational settling, thus accounting naturally for the presence of He, as well as traces of metals such as Mg and Ca, typically found in the atmospheres of ELM proto-WDs. Here we investigate whether rotational mixing can maintain a sufficient amount of He in the deeper driving region of the star, such that it can fuel, through HeII-HeIII ionization, the observed pulsations in this type of stars. Using state-of-the-art evolutionary models computed with MESA, we show that rotational mixing can indeed explain qualitatively the very existence and general properties of the known pulsating, mixed-atmosphere ELM proto-WDs. Moreover, such objects are very likely to pulsate again during their final WD cooling phase.

preprint2016arXiv

B field in OB stars (BOB): The outstandingly strong magnetic field in the evolved He-strong star CPD-62 2124

The origin and evolution of magnetism in OB stars is far from being well understood. With approximately 70 magnetic OB stars known, any new object with unusual characteristics may turn out to be a key piece of the puzzle. We report the detection of an exceptionally strong magnetic field in the He-strong B2IV star CPD-62 2124. Spectropolarimetric FORS2 and HARPSpol observations were analysed by two independent teams and procedures, concluding on a strong longitudinal magnetic field of approximately 5.2 kG. The quantitative characterisation of the stellar atmosphere yields an effective temperature of 23650$\pm$250 K, a surface gravity of 3.95$\pm$0.10 dex and a surface helium fraction of 0.35$\pm$0.02 by number. The metal composition is in agreement with the cosmic abundance standard, except for Mg, Si and S, which are slightly non-solar. The strong and broad ($\sim$300 km/s) disc-like emission displayed by the H$α$ line suggests a centrifugal magnetosphere supported by the strong magnetic field. Our results imply that CPD-62 2124 is an early B-type star hosting one of the strongest magnetic fields discovered to date, and one of the most evolved He-strong stars known, with a fractional main-sequence lifetime of approximately 0.6.

preprint2016arXiv

B fields in OB stars (BOB): Detection of a magnetic field in the He-strong star CPD-57° 3509

We report the detection of a magnetic field in the helium-strong star CPD-57 3509 (B2 IV), a member of the Galactic open cluster NGC3293, and characterise the star's atmospheric and fundamental parameters. Spectropolarimetric observations with FORS2 and HARPSpol are analysed using two independent approaches to quantify the magnetic field strength. A high-S/N FLAMES/GIRAFFE spectrum is analysed using a hybrid non-LTE model atmosphere technique. Comparison with stellar evolution models constrains the fundamental parameters of the star. We obtain a firm detection of a surface averaged longitudinal magnetic field with a maximum amplitude of about 1 kG. Assuming a dipolar configuration of the magnetic field, this implies a dipolar field strength larger than 3.3 kG. Moreover, the large amplitude and fast variation (within about 1 day) of the longitudinal magnetic field implies that CPD-57 3509 is spinning very fast despite its apparently slow projected rotational velocity. The star should be able to support a centrifugal magnetosphere, yet the spectrum shows no sign of magnetically confined material; in particular, emission in Hα is not observed. Apparently, the wind is either not strong enough for enough material to accumulate in the magnetosphere to become observable or, alternatively, some leakage process leads to loss of material from the magnetosphere. The quantitative spectroscopic analysis of the star yields an effective temperature and a logarithmic surface gravity of 23750+-250 K and 4.05+-0.10, respectively, and a surface helium fraction of 0.28+-0.02 by number. The surface abundances of C, N, O, Ne, S, and Ar are compatible with the cosmic abundance standard, whereas Mg, Al, Si, and Fe are depleted by about a factor of 2. This abundance pattern can be understood as the consequence of a fractionated stellar wind. CPD-57 3509 is one of the most evolved He-strong stars known.

preprint2016arXiv

Common envelope ejection in massive binary stars - Implications for the progenitors of GW150914 and GW151226

The recently detected gravitational wave signals (GW150914 and GW151226) of the merger event of a pair of relatively massive stellar-mass black holes (BHs) calls for an investigation of the formation of such progenitor systems in general. We analyse the common envelope (CE) stage of the "traditional" formation channel in binaries where the first-formed compact object undergoes an in-spiral inside the envelope of its evolved companion star and ejects the envelope in that process. We calculate envelope binding energies of donor stars with initial masses between 4 and 115 Msun for metallicities of Z=Zsun/2 and Z=Zsun/50, and derive minimum masses of in-spiralling objects needed to eject these envelopes. We find that CE evolution, besides from producing WD-WD and NS-NS binaries, may, in principle, also produce massive BH-BH systems with individual BH component masses up to ~50-60 Msun, in particular for donor stars evolved to giants. However, the physics of envelope ejection of massive stars remains uncertain. We discuss the applicability of the energy-budget formalism, the location of the bifurcation point, the recombination energy and the accretion energy during in-spiral as possible energy sources, and also comment on the effect of inflated helium cores. Massive stars in a wide range of metallicities and with initial masses up to at least 115 Msun may possibly shed their envelopes and survive CE evolution, depending on their initial orbital parameters, similarly to the situation for intermediate mass and low-mass stars with degenerate cores. We conclude that based on stellar structure calculations, and in the view of the usual simple energy budget analysis, events like GW150914 and GW151226 could possibly be produced from the CE channel. Calculations of post-CE orbital separations, however, and thus the estimated LIGO detection rates, remain highly uncertain. [Abridged]

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

Imprints of the ejecta-companion interaction in Type Ia supernovae: main sequence, subgiant, and red giant companions

We study supernova ejecta-companion interactions in a sample of realistic semidetached binary systems representative of Type Ia supernova progenitor binaries in a single-degenerate scenario. We model the interaction process with the help of a high-resolution hydrodynamic code assuming cylindrical symmetry. We find that the ejecta hole has a half-opening angle of 40--50$^\circ$ with the density by a factor of 2-4 lower, in good agreement with the previous studies. Quantitative differences from the past results in the amounts and kinematics of the stripped companion material and levels of contamination of the companion with the ejecta material can be explained by different model assumptions and effects due to numerical diffusion.We analyse and, for the first time, provide simulation-based estimates of the amounts and of the thermal characteristics of the shock-heated material responsible for producing a prompt, soft X-ray emission. Besides the shocked ejecta material, considered in the original model by Kasen, we also account for the stripped, shock-heated envelope material of stellar companions, which we predict partially contributes to the prompt emission. The amount of the energy deposited in the envelope is comparable to the energy stored in the ejecta. The total energy budget available for the prompt emission is by a factor of about 2-4 smaller than originally predicted by Kasen. Although the shocked envelope has a higher characteristic temperature than the shocked ejecta, the temperature estimates of the shocked material are in good agreement with the Kasen's model. The hottest shocked plasma is produced in the subgiant companion case.

preprint2016arXiv

Models for the evolution of close binaries with He-Star and WD components towards Type Ia supernova explosions

Type Ia supernovae (SNe Ia) have been an important tool for astronomy for quite some time; however, the nature of their progenitors remains somewhat mysterious. Recent theoretical studies indicated the possibility of producing thermonuclear detonations of carbon-oxygen white dwarfs (CO WDs) at masses less than the Chandrasekhar mass through accretion of helium-rich matter, which would, depending on mass accretion rate, mass, and initial temperature of the WD, spectrally resemble either a normal SN Ia or a peculiar one. This study aims to further resolve the state of binary systems comprised of a sub-Chandrasekhar-mass CO WD and a helium star at the point where an accretion-induced detonation occurs and constrains the part of the initial parameter space where this kind of phenomenon is possible. Preexisting data obtained through simulations of single, constantly accreting CO WDs is used as an indicator for the behavior of new binary models in which the WD is treated as a point mass and which include the non-degenerate partner as a fully resolved stellar model. We parameterize the ignition of the accumulated helium layer, changes in the WD temperature, and changes in the CO core mass depending on the mass transfer rate. The initial conditions allowing for detonation do not form a single contiguous area in the parameter space, whose shape is profoundly influenced by the behavior of the donor star. Mass loss due to Nova outbursts acts in favor of detonation. According to our criteria, about 10% of the detonations in this study can be expected to show spectra consistent with ordinary SNe Ia; the rest exhibit peculiar features.

preprint2016arXiv

Rotating Stars and the Formation of Bipolar Planetary Nebulae II: Tidal Spin-up

We present new binary stellar evolution models that include the effects of tidal forces, rotation, and magnetic torques with the goal of testing Planetary Nebulae (PNe) shaping via binary interaction. We explore whether tidal interaction with a companion can spin up the AGB envelope. To do so we have selected binary systems with main sequence masses of 2.5 \Mo and of 0.8 \Mo and evolve them allowing initial separations of 5, 6, 7, and 8 AU. The binary stellar evolution models have been computed all the way to the PNe formation phase or until Roche lobe overflow (RLOF) is reached, whatever happens first. We show that with initial separations of 7 and 8 AU, the binary avoids entering into RLOF, and the AGB star reaches moderate rotational velocities at the surface ($\sim 3.5 $ and $\sim 2 $ \kms respectively) during the inter-pulse phases, but after the thermal pulses it drops to a final rotational velocity of only $\sim 0.03 $ \kms. For the closest binary separations explored, 5 and 6 AU, the AGB star reaches rotational velocities of $\sim 6 $ and $\sim 4 $ \kms respectively when the RLOF is initiated. We conclude that the detached binary models that avoid entering the RLOF phase during the AGB will not shape bipolar PNe, since the acquired angular momentum is lost via the wind during the last two thermal pulses. This study rules out tidal spin-up in non-contact binaries as a sufficient condition to form bipolar PNe.

preprint2016arXiv

The metallicity dependence of envelope inflation in massive stars

Recently it has been found that models of massive stars reach the Eddington limit in their interior, which leads to dilute extended envelopes. We perform a comparative study of the envelope properties of massive stars at different metallicities, with the aim to establish the impact of the stellar metallicity on the effect of envelope inflation. We analyse published grids of core-hydrogen burning massive star models computed with metallicities appropriate for massive stars in the Milky Way, the LMC and the SMC, the very metal poor dwarf galaxy I Zwicky 18, and for metal-free chemical composition. Stellar models of all the investigated metallicities reach and exceed the Eddington limit in their interior, aided by the opacity peaks of iron, helium and hydrogen, and consequently develop inflated envelopes. Envelope inflation leads to a redward bending of the zero-age main sequence and a broadening of the main sequence band in the upper part of the Hertzsprung-Russell diagram. We derive the limiting L/M-values as function of the stellar surface temperature above which inflation occurs, and find them to be larger for lower metallicity. While Galactic models show inflation above ~29 Msun, the corresponding mass limit for Population III stars is ~150 Msun. While the masses of the inflated envelopes are generally small, we find that they can reach 1-100 Msun in models with effective temperatures below ~8000 K, with higher masses reached by models of lower metallicity. Envelope inflation is expected to occur in sufficiently massive stars at all metallicities, and is expected to lead to rapidly growing pulsations, high macroturbulent velocities, and might well be related to the unexplained variability observed in Luminous Blue Variables like S Doradus and Eta Carina.

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

Asymmetric supernova remnants generated by Galactic, massive runaway stars

After the death of a runaway massive star, its supernova shock wave interacts with the bow shocks produced by its defunct progenitor, and may lose energy, momentum, and its spherical symmetry before expanding into the local interstellar medium (ISM). We investigate whether the initial mass and space velocity of these progenitors can be associated with asymmetric supernova remnants. We run hydrodynamical models of supernovae exploding in the pre-shaped medium of moving Galactic core-collapse progenitors. We find that bow shocks that accumulate more than about 1.5 Mo generate asymmetric remnants. The shock wave first collides with these bow shocks 160-750 yr after the supernova, and the collision lasts until 830-4900 yr. The shock wave is then located 1.35-5 pc from the center of the explosion, and it expands freely into the ISM, whereas in the opposite direction it is channelled into the region of undisturbed wind material. This applies to an initially 20 Mo progenitor moving with velocity 20 km/s and to our initially 40 Mo progenitor. These remnants generate mixing of ISM gas, stellar wind and supernova ejecta that is particularly important upstream from the center of the explosion. Their lightcurves are dominated by emission from optically-thin cooling and by X-ray emission of the shocked ISM gas. We find that these remnants are likely to be observed in the [OIII] lambda 5007 spectral line emission or in the soft energy-band of X-rays. Finally, we discuss our results in the context of observed Galactic supernova remnants such as 3C391 and the Cygnus Loop.

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): FORS2 spectropolarimetric follow-up of the two rare rigidly rotating magnetosphere stars HD23478 and HD345439

Massive B-type stars with strong magnetic fields and fast rotation are very rare and provide a mystery for theories of both star formation and magnetic field evolution. Only two such stars, called sigma Ori E analogs, were previously known. Recently, a team involved in APOGEE, one of the Sloan Digital Sky Survey III programs, announced the discovery of two additional rigidly rotating magnetosphere stars, HD23478 and HD345439. The presence of magnetic fields in these newly discovered sigma Ori E analogs was not investigated in the past. In the framework of our ESO Large Programme, and one normal ESO programme, we carried out low-resolution FORS2 spectropolarimetric observations of HD23478 and HD345439. From the measurements using hydrogen lines, we discover a rather strong longitudinal magnetic field of the order of up to 1.5kG in HD23478, and up to 1.3kG using the entire spectrum. The analysis of HD345439 using four subsequent spectropolarimetric subexposures does not reveal the presence of a magnetic field at a significance level of 3sigma. On the other hand, the inspection of individual subexposures indicates that HD345439 may host a strong magnetic field, rapidly varying over 88 minutes. A hint at the fast rotation of HD345439 is also given by the behaviour of several metallic and He I lines in the low-resolution FORS2 spectra, showing profile variations already on such a short time scale.

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

Blue supergiants as descendants of magnetic main sequence stars

About 10$\%$ of the massive main sequence stars have recently been found to host a strong, large scale magnetic field. Both, the origin and the evolutionary consequences of these fields are largely unknown. We argue that these fields may be sufficiently strong in the deep interior of the stars to suppress convection near the outer edge of their convective core. We performed parametrised stellar evolution calculations and assumed a reduced size of the convective core for stars in the mass range 16 M$_{\odot}$ to 28 M$_{\odot}$ from the zero age main sequence until core carbon depletion. We find that such models avoid the coolest part of the main sequence band, which is usually filled by evolutionary models that include convective core overshooting. Furthermore, our `magnetic' models populate the blue supergiant region during core helium burning, i.e., the post-main sequence gap left by ordinary single star models, and some of them end their life in a position near that of the progenitor of Supernova 1987A in the HR diagram. Further effects include a strongly reduced luminosity during the red supergiant stage, and downward shift of the limiting initial mass for white dwarf and neutron star formation.

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

Evolution of Mass Functions of Coeval Stars through Wind Mass Loss and Binary Interactions

Accurate determinations of stellar mass functions and ages of stellar populations are crucial to much of astrophysics. We analyse the evolution of stellar mass functions of coeval main sequence stars including all relevant aspects of single- and binary-star evolution. We show that the slope of the upper part of the mass function in a stellar cluster can be quite different to the slope of the initial mass function. Wind mass loss from massive stars leads to an accumulation of stars which is visible as a peak at the high mass end of mass functions, thereby flattening the mass function slope. Mass accretion and mergers in close binary systems create a tail of rejuvenated binary products. These blue straggler stars extend the single star mass function by up to a factor of two in mass and can appear up to ten times younger than their parent stellar cluster. Cluster ages derived from their most massive stars that are close to the turn-off may thus be significantly biased. To overcome such difficulties, we propose the use of the binary tail of stellar mass functions as an unambiguous clock to derive the cluster age because the location of the onset of the binary tail identifies the cluster turn-off mass. It is indicated by a pronounced jump in the mass function of old stellar populations and by the wind mass loss peak in young stellar populations. We further characterise the binary induced blue straggler population in star clusters in terms of their frequency, binary fraction and apparent age.

preprint2015arXiv

Low-metallicity massive single stars with rotation. Evolutionary models applicable to I Zwicky 18

Massive rotating single stars with an initial metal composition appropriate for the dwarf galaxy I Zw 18 ([Fe/H]=$-$1.7) are modelled during hydrogen burning for initial masses of 9-300 M$_{\odot}$ and rotational velocities of 0-900 km s$^{-1}$. Internal mixing processes in these models were calibrated based on an observed sample of OB-type stars in the Magellanic Clouds. Even moderately fast rotators, which may be abundant at this metallicity, are found to undergo efficient mixing induced by rotation resulting in quasi chemically-homogeneous evolution. These homogeneously-evolving models reach effective temperatures of up to 90 kK during core hydrogen burning. This, together with their moderate mass-loss rates, make them Transparent Wind Ultraviolet INtense stars (TWUIN star), and their expected numbers might explain the observed HeII ionizing photon flux in I Zw 18 and other low-metallicity HeII galaxies. Our slowly rotating stars above $\sim$80 M$_{\odot}$ evolve into late B- to M-type supergiants during core hydrogen burning, with visual magnitudes up to 19$^{\mathrm{m}}$ at the distance of I Zw 18. Both types of stars, TWUIN stars and luminous late-type supergiants, are only predicted at low metallicity. Massive star evolution at low metallicity is shown to differ qualitatively from that in metal-rich environments. Our grid can be used to interpret observations of local star-forming dwarf galaxies and high-redshift galaxies, as well as the metal-poor components of our Milky Way and its globular clusters.

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

The blue supergiant MN18 and its bipolar circumstellar nebula

We report the results of spectrophotometric observations of the massive star MN18 revealed via discovery of a bipolar nebula around it with the Spitzer Space Telescope. Using the optical spectrum obtained with the Southern African Large Telescope, we classify this star as B1 Ia. The evolved status of MN18 is supported by the detection of nitrogen overabundance in the nebula, which implies that it is composed of processed material ejected by the star. We analysed the spectrum of MN18 by using the code CMFGEN, obtaining a stellar effective temperature of \approx 21 kK. The star is highly reddened, E(B-V)\approx 2 mag. Adopting an absolute visual magnitude of M_V=-6.8\pm0.5 (typical of B1 supergiants), MN18 has a luminosity of log L/Lsun \approx 5.42\pm0.30, a mass-loss rate of \approx (2.8-4.5)\times10^{-7} Msun/yr, and resides at a distance of \approx 5.6^{+1.5} _{-1.2} kpc. We discuss the origin of the nebula around MN18 and compare it with similar nebulae produced by other blue supergiants in the Galaxy (Sher 25, HD 168625, [SBW2007] 1) and the Large Magellanic Cloud (Sk-69 202). The nitrogen abundances in these nebulae imply that blue supergiants can produce them from the main sequence stage up to the pre-supernova stage. We also present a K-band spectrum of the candidate luminous blue variable MN56 (encircled by a ring-like nebula) and report the discovery of an OB star at \approx 17 arcsec from MN18. The possible membership of MN18 and the OB star of the star cluster Lynga 3 is discussed.

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

A VLT/FLAMES survey for massive binaries in Westerlund 1. IV. Wd1-5 - binary product and a pre-supernova companion for the magnetar CXOU J1647-45?

The first soft gamma-ray repeater was discovered over three decades ago, and subsequently identified as a magnetar. However there is currently no consenus on the formation channel(s) of these objects. The presence of a magnetar in the starburst cluster Westerlund 1 implies a >40Msun progenitor, favouring formation in a binary that was disrupted at supernova. To test this hypothesis we searched for the putative pre-SN companion via a radial velocity survey to identify high-velocity runaways, with subsequent atmospheric analysis of the resultant candidate, Wd1-5. Wd1-5 appears overluminous for its spectroscopic mass and we find evidence of He- and N-enrichment, O-depletion, and critically C-enrichment, a combination of properties that is difficult to explain for a single star. We infer a pre-SN history for Wd1-5 which supposes an initial close binary comprising two stars of comparable (~41Msun+35Msun) masses. Efficient mass transfer from the initially more massive component leads to the mass-gainer evolving more rapidly, initiating luminous blue variable/common envelope evolution. Reverse, wind-driven mass transfer during its subsequent WC Wolf-Rayet phase leads to the carbon pollution of Wd1-5, before a type Ibc supernova disrupts the binary system. Under the assumption of a physical association between Wd1-5 and J1647-45, the secondary is identified as the magnetar progenitor; its common envelope evolutionary phase prevents spin-down of its core prior to SN and the seed magnetic field for the magnetar forms either in this phase or during the earlier episode of mass transfer in which it was spun-up. Our results suggest that binarity is a key ingredient in the formation of at least a subset of magnetars by preventing spin-down via core-coupling and potentially generating a seed magnetic field.

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

B fields in OB stars (BOB): The discovery of a magnetic field in a multiple system in the Trifid Nebula, one of the youngest star forming regions

Recent magnetic field surveys in O- and B-type stars revealed that about 10% of the core-hydrogen-burning massive stars host large-scale magnetic fields. The physical origin of these fields is highly debated. To identify and model the physical processes responsible for the generation of magnetic fields in massive stars, it is important to establish whether magnetic massive stars are found in very young star-forming regions or whether they are formed in close interacting binary systems. In the framework of our ESO Large Program, we carried out low-resolution spectropolarimetric observations with FORS2 in 2013 April of the three most massive central stars in the Trifid nebula, HD164492A, HD164492C, and HD164492D. These observations indicated a strong longitudinal magnetic field of about 500-600G in the poorly studied component HD164492C. To confirm this detection, we used HARPS in spectropolarimetric mode on two consecutive nights in 2013 June. Our HARPS observations confirmed the longitudinal magnetic field in HD164492C. Furthermore, the HARPS observations revealed that HD164492C cannot be considered as a single star as it possesses one or two companions. The spectral appearance indicates that the primary is most likely of spectral type B1-B1.5V. Since in both observing nights most spectral lines appear blended, it is currently unclear which components are magnetic. Long-term monitoring using high-resolution spectropolarimetry is necessary to separate the contribution of each component to the magnetic signal. Given the location of the system HD164492C in one of the youngest star formation regions, this system can be considered as a Rosetta Stone for our understanding of the origin of magnetic fields in massive stars.

preprint2014arXiv

Discovery of a new Galactic bona fide luminous blue variable with Spitzer

We report the discovery of a circular mid-infrared shell around the emission-line star Wray 16-137 using archival data of the Spitzer Space Telescope. Follow-up optical spectroscopy of Wray 16-137 with the Southern African Large Telescope revealed a rich emission spectrum typical of the classical luminous blue variables (LBVs) like P Cygni. Subsequent spectroscopic and photometric observations showed drastic changes in the spectrum and brightness during the last three years, meaning that Wray 16-137 currently undergoes an S Dor-like outburst. Namely, we found that the star has brightened by \approx 1 mag in the V and I_c bands, while its spectrum became dominated by Fe ii lines. Taken together, our observations unambiguously show that Wray 16-137 is a new member of the family of Galactic bona fide LBVs.

preprint2014arXiv

Discovery of a new Wolf-Rayet star and a candidate star cluster in the Large Magellanic Cloud with Spitzer

We report the first-ever discovery of a Wolf-Rayet (WR) star in the Large Magellanic Cloud via detection of a circular shell with the Spitzer Space Telescope. Follow-up observations with Gemini-South resolved the central star of the shell into two components separated from each other by approx 2 arcsec (or approx 0.5 pc in projection). One of these components turns out to be a WN3 star with H and He lines both in emission and absorption (we named it BAT99 3a using the numbering system based on extending the Breysacher et al. catalogue). Spectroscopy of the second component showed that it is a B0 V star. Subsequent spectroscopic observations of BAT99 3a with the du Pont 2.5-m telescope and the Southern African Large Telescope revealed that it is a close, eccentric binary system, and that the absorption lines are associated with an O companion star. We analyzed the spectrum of the binary system using the non-LTE Potsdam Wolf-Rayet (PoWR) code, confirming that the WR component is a very hot (approx 90 kK) WN star. For this star, we derived a luminosity of log L/Lsun =5.45 and a mass-loss rate of 10^{-5.8} Msun/yr, and found that the stellar wind composition is dominated by helium with 20 per cent of hydrogen. Spectroscopy of the shell revealed an He iii region centred on BAT99 3a and having the same angular radius (approx 15 arcsec) as the shell. We thereby add a new example to a rare class of high-excitation nebulae photoionized by WR stars. Analysis of the nebular spectrum showed that the shell is composed of unprocessed material, implying that the shell was swept-up from the local interstellar medium. We discuss the physical relationship between the newly identified massive stars and their possible membership of a previously unrecognized star cluster.

preprint2014arXiv

Instability of Magnetic Equilibria in Barotropic Stars

In stably stratified stars, numerical magneto-hydrodynamics simulations have shown that arbitrary initial magnetic fields evolve into stable equilibrium configurations, usually containing nearly axisymmetric, linked poloidal and toroidal fields that stabilize each other. In this work, we test the hypothesis that stable stratification is a requirement for the existence of such stable equilibria. For this purpose, we follow numerically the evolution of magnetic fields in barotropic (and thus neutrally stable) stars, starting from two different types of initial conditions, namely random disordered magnetic fields, as well as linked poloidal-toroidal configurations resembling the previously found equilibria. With many trials, we always find a decay of the magnetic field over a few Alfvén times, never a stable equilibrium. This strongly suggests that there are no stable equilibria in barotropic stars, thus clearly invalidating the assumption of barotropic equations of state often imposed on the search of magnetic equilibria. It also supports the hypothesis that, as dissipative processes erode the stable stratification, they might destabilize previously stable magnetic field configurations, leading to their decay.

preprint2014arXiv

Models of the circumstellar medium of evolving, massive runaway stars moving through the Galactic plane

At least 5 per cent of the massive stars are moving supersonically through the interstellar medium (ISM) and are expected to produce a stellar wind bow shock. We explore how the mass loss and space velocity of massive runaway stars affect the morphology of their bow shocks. We run two-dimensional axisymmetric hydrodynamical simulations following the evolution of the circumstellar medium of these stars in the Galactic plane from the main sequence to the red supergiant phase. We find that thermal conduction is an important process governing the shape, size and structure of the bow shocks around hot stars, and that they have an optical luminosity mainly produced by forbidden lines, e.g. [OIII]. The Ha emission of the bow shocks around hot stars originates from near their contact discontinuity. The H$α$ emission of bow shocks around cool stars originates from their forward shock, and is too faint to be observed for the bow shocks that we simulate. The emission of optically-thin radiation mainly comes from the shocked ISM material. All bow shock models are brighter in the infrared, i.e. the infrared is the most appropriate waveband to search for bow shocks. Our study suggests that the infrared emission comes from near the contact discontinuity for bow shocks of hot stars and from the inner region of shocked wind for bow shocks around cool stars. We predict that, in the Galactic plane, the brightest, i.e. the most easily detectable bow shocks are produced by high-mass stars moving with small space velocities.

preprint2014arXiv

Observational properties of low redshift pair instability supernovae

So called superluminous supernovae have been recently discovered in the local Universe. It appears possible that some of them originate from stellar explosions induced by the pair instability mechanism. Recent stellar evolution models also predict pair instability supernovae (PISNe) from very massive stars at fairly high metallicities (i.e. Z~0.004). We provide supernova (SN) models and synthetic light curves for two progenitor models, a 150 Msun red-supergiant and a 250 Msun yellow-supergiant at a metallicity of Z=0.001, for which the evolution from the main sequence to collapse, and the initiation of the PISN itself, has been previously computed in a realistic and self-consistent way. We are using the radiation hydrodynamics code STELLA to describe the SN evolution of both models over a time frame of about 500 days. We describe the shock-breakout phases of both SNe which are characterized by a higher luminosity, a longer duration and a lower effective temperature than those of ordinary SNeIIP. We derive the bolometric as well as the U, B, V, R and I light curves of our PISN models, which show a long-lasting plateau phase with maxima at Mbol=-19.3 mag and -21.3 mag for our lower and higher mass model, respectively. We also describe the photospheric composition and velocity as function of time. We conclude that the light curve of the explosion of our initially 150 Msun star resembles those of relatively bright SNeIIP, whereas its photospheric velocity at early times is smaller. Its 56Ni mass of 0.04 Msun also falls well into the range found in ordinary core collapse SNe. The light curve and photospheric velocity of our 250 Msun models has a striking resemblance with that of the superluminous SN2007bi, strengthening its interpretation as PISN. We conclude that PISNe may occur more frequently in the local universe than previously assumed.

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

Single Rotating Stars and the Formation of Bipolar Planetary Nebula

We have computed new stellar evolution models that include the effects of rotation and magnetic torques under different hypothesis. The goal is to test if a single star can sustain in the envelope the rotational velocities needed for the magneto hydrodynamical (MHD) simulations to shape bipolar Planetary Nebulae (PNe) when the high mass-loss rates take place. Stellar evolution models with main sequence masses of 2.5 and 5 Mo, and initial rotational velocities of 250 km/s have been followed all the way to the PNe formation phase. We find that stellar cores have to be spun down using magnetic torques in order to reproduce the rotation rates observed for white dwarfs. During the asymptotic giant branch phase and beyond, the magnetic braking of the core has a practically null effect in increasing the rotational velocity of the envelope since the stellar angular momentum is removed efficiently by the wind. We have, as well, tested best possible case scenarios in rather non-physical contexts to give enough angular momentum to the envelope. We find that we cannot get the envelope of a single star rotating at the speeds needed by the MHD simulations to form bipolar PNe. We conclude that single stellar rotators are unlikely to be the progenitors of bipolar PNe under the current MHD model paradigm.

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 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 spectroscopic Hertzsprung-Russell diagram

The Hertzsprung-Russell diagram is an essential diagnostic diagram for stellar structure and evolution, which has now been in use for more than 100 years. Our spectroscopic Hertzsprung-Russell (sHR) diagram shows the inverse of the flux-mean gravity versus the effective temperature. Observed stars whose spectra have been quantitatively analyzed can be entered in this diagram without the knowledge of the stellar distance or absolute brightness. Observed stars can be as conveniently compared to stellar evolution calculations in the sHR diagram as in the Hertzsprung-Russell diagram. However, at the same time, our ordinate is proportional to the stellar mass-to-luminosity ratio, which can thus be directly determined. For intermediate- and low-mass star evolution at constant mass, we show that the shape of an evolutionary track in the sHR diagram is identical to that in the Hertzsprung-Russell diagram. We also demonstrate that for hot stars, their stellar Eddington factor can be directly read off the sHR diagram. For stars near their Eddington limit, we argue that a version of the sHR diagram may be useful where the gravity is exchanged by the effective gravity. We discuss the advantages and limitations of the sHR diagram, and show that it can be fruitfully applied to Galactic stars, but also to stars with known distance, e.g., in the LMC or in galaxies beyond the Local Group.

preprint2014arXiv

The spectroscopic Hertzsprung-Russell diagram of Galactic massive stars

The distribution of stars in the Hertzsprung-Russell diagram narrates their evolutionary history and directly assesses their properties. Placing stars in this diagram however requires the knowledge of their distances and interstellar extinctions, which are often poorly known for Galactic stars. The spectroscopic Hertzsprung-Russell diagram (sHRD) tells similar evolutionary tales, but is independent of distance and extinction measurements. Based on spectroscopically derived effective temperatures and gravities of almost 600 stars, we derive for the first time the observational distribution of Galactic massive stars in the sHRD. While biases and statistical limitations in the data prevent detailed quantitative conclusions at this time, we see several clear qualitative trends. By comparing the observational sHRD with different state-of-the-art stellar evolutionary predictions, we conclude that convective core overshooting may be mass-dependent and, at high mass ($\geq 15\,M_\odot$), stronger than previously thought. Furthermore, we find evidence for an empirical upper limit in the sHRD for stars with $T_{\rm{eff}}$ between 10000 and 32000 K and, a strikingly large number of objects below this line. This over-density may be due to inflation expanding envelopes in massive main-sequence stars near the Eddington limit.

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.

preprint2014arXiv

Two spotted and magnetic early B-type stars in the young open cluster NGC2264 discovered by MOST and ESPaDOnS

Star clusters are known as superb tools for understanding stellar evolution. In a quest for understanding the physical origin of magnetism and chemical peculiarity in about 7% of the massive main-sequence stars, we analysed two of the ten brightest members of the ~10 Myr old Galactic open cluster NGC 2264, the early B-dwarfs HD47887 and HD47777. We find accurate rotation periods of 1.95 and 2.64 days, respectively, from MOST photometry. We obtained ESPaDOnS spectropolarimetric observations, through which we determined stellar parameters, detailed chemical surface abundances, projected rotational velocities, and the inclination angles of the rotation axis. Because we found only small (<5 km/s) radial velocity variations, most likely caused by spots, we can rule out that HD47887 and HD47777 are close binaries. Finally, using the least-squares deconvolution technique, we found that both stars possess a large-scale magnetic field with an average longitudinal field strength of about 400 G. From a simultaneous fit of the stellar parameters we determine the evolutionary masses of HD47887 and HD47777 to be 9.4+/-0.7 M0 and 7.6+/-0.5 M0. Interestingly, HD47777 shows a remarkable helium underabundance, typical of helium-weak chemically peculiar stars, while the abundances of HD47887 are normal, which might imply that diffusion is operating in the lower mass star but not in the slightly more massive one. Furthermore, we argue that the rather slow rotation, as well as the lack of nitrogen enrichment in both stars, can be consistent with both the fossil and the binary hypothesis for the origin of the magnetic field. However, the presence of two magnetic and apparently single stars near the top of the cluster mass-function may speak in favour of the latter.

preprint2013arXiv

Exploring the origin of magnetic fields in massive stars: II. New magnetic field measurements in cluster and field stars

Theories on the origin of magnetic fields in massive stars remain poorly developed, because the properties of their magnetic field as function of stellar parameters could not yet be investigated. To investigate whether magnetic fields in massive stars are ubiquitous or appear only in stars with a specific spectral classification, certain ages, or in a special environment, we acquired 67 new spectropolarimetric observations for 30 massive stars. Among the observed sample, roughly one third of the stars are probable members of clusters at different ages, whereas the remaining stars are field stars not known to belong to any cluster or association. Spectropolarimetric observations were obtained during four different nights using the low-resolution spectropolarimetric mode of FORS2 (FOcal Reducer low dispersion Spectrograph) mounted on the 8-m Antu telescope of the VLT. Furthermore, we present a number of follow-up observations carried out with the high-resolution spectropolarimeters SOFIN mounted at the Nordic Optical Telescope (NOT) and HARPS mounted at the ESO 3.6m between 2008 and 2011. To assess the membership in open clusters and associations, we used astrometric catalogues with the highest quality kinematic and photometric data currently available. The presence of a magnetic field is confirmed in nine stars previously observed with FORS1/2: HD36879, HD47839, CPD-282561, CPD-472963, HD93843, HD148937, HD149757, HD328856, and HD164794. New magnetic field detections at a significance level of at least 3sigma were achieved in five stars: HD92206c, HD93521, HD93632, CPD-468221, and HD157857. Among the stars with a detected magnetic field, five stars belong to open clusters with high membership probability. According to previous kinematic studies, five magnetic O-type stars in our sample are candidate runaway stars.

preprint2013arXiv

IRC-10414: a bow-shock-producing red supergiant star

Most runaway OB stars, like the majority of massive stars residing in their parent clusters, go through the red supergiant (RSG) phase during their lifetimes. Nonetheless, although many dozens of massive runaways were found to be associated with bow shocks, only two RSG bow-shock-producing stars, Betelgeuse and μCep, are known to date. In this paper, we report the discovery of an arc-like nebula around the late M-type star IRC-10414 using the SuperCOSMOS H-alpha Survey. Our spectroscopic follow-up of IRC-10414 with the Southern African Large Telescope (SALT) showed that it is a M7 supergiant, which supports previous claims on the RSG nature of this star based on observations of its maser emission. This was reinforced by our new radio- and (sub)millimeter-wavelength molecular line observations made with the Atacama Pathfinder Experiment (APEX) 12 meter telescope and the Effelsberg 100 m radio telescope, respectively. The SALT spectrum of the nebula indicates that its emission is the result of shock excitation. This finding along with the arc-like shape of the nebula and an estimate of the space velocity of IRC-10414 (\approx 70\pm20 km/s) imply the bow shock interpretation for the nebula. Thus, IRC-10414 represents the third case of a bow-shock-producing RSG and the first one with a bow shock visible at optical wavelengths. We discuss the smooth appearance of the bow shocks around IRC-10414 and Betelgeuse and propose that one of the necessary conditions for stability of bow shocks generated by RSGs is the ionization of the stellar wind. Possible ionisation sources of the wind of IRC-10414 are proposed and discussed.

preprint2013arXiv

Magnetic Fields in Stars: Origin and Impact

Various types of magnetic fields occur in stars: small scale fields, large scale fields, and internal toroidal fields. While the latter may be ubiquitous in stars due to differential rotation, small scale fields (spots) may be associated with envelop convection in all low and high mass stars. The stable large scale fields found in only about 10 per cent of intermediate mass and massive stars may be understood as a consequence of dynamical binary interaction, e.g., the merging of two stars in a binary. We relate these ideas to magnetic fields in white dwarfs and neutron stars, and to their role in core-collapse and thermonuclear supernova explosions.

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

On the stability of bow shocks generated by red supergiants: the case of IRC-10414

In this Letter, we explore the hypothesis that the smooth appearance of bow shocks around some red supergiants (RSGs) might be caused by the ionization of their winds by external sources of radiation. Our numerical simulations of the bow shock generated by IRC-10414 (the first-ever RSG with an optically detected bow shock) show that the ionization of the wind results in its acceleration by a factor of two, which reduces the difference between the wind and space velocities of the star and makes the contact discontinuity of the bow shock stable for a range of stellar space velocities and mass-loss rates. Our best fit model reproduces the overall shape and surface brightness of the observed bow shock and suggests that the space velocity and mass-loss rate of IRC-10414 are $\approx$50 ${\rm km} \, {\rm s}^{-1}$ and $\approx$$10^{-6}$ $M_\odot \, {\rm yr}^{-1}$, respectively, and that the number density of the local ISM is $\approx$3 ${\rm cm}^{-3}$. It also shows that the bow shock emission comes mainly from the shocked stellar wind. This naturally explains the enhanced nitrogen abundance in the line-emitting material, derived from the spectroscopy of the bow shock. We found that photoionized bow shocks are $\approx$15$-$50 times brighter in optical line emission than their neutral counterparts, from which we conclude that the bow shock of IRC-10414 must be photoionized.

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

Search for Stable Magnetohydrodynamic Equilibria in Barotropic Stars

It is now believed that magnetohydrodynamic equilibria can exist in stably stratified stars due to the seminal works of Braithwaite & Spruit (2004) and Braithwaite & Nordlund (2006). What is still not known is whether magnetohydrodynamic equilibria can exist in a barotropic star, in which stable stratification is not present. It has been conjectured by Reisenegger (2009) that there will likely not exist any magnetohydrodynamical equilibria in barotropic stars. We aim to test this claim by presenting preliminary MHD simulations of barotropic stars using the three dimensional stagger code of Nordlund & Galsgaard (1995).

preprint2013arXiv

Spectroscopic and physical parameters of Galactic O-type stars. II. Observational constraints on projected rotational and extra broadening velocities as a function of fundamental parameters and stellar evolution

Rotation is of key importance for the evolution of hot massive stars, however, the rotational velocities of these stars are difficult to determine. Based on our own data for 31 Galactic O stars and incorporating similar data for 86 OB supergiants from the literature, we aim at investigating the properties of rotational and extra line-broadening as a function of stellar parameters and at testing model predictions about the evolution of stellar rotation. Fundamental stellar parameters were determined by means of the code FASTWIND. Projected rotational and extra broadening velocities originate from a combined Ft + GOF method. Model calculations published previously were used to estimate the initial evolutionary masses. The sample O stars with Minit > 50 Msun rotate with less that 26% of their break-up velocity, and they also lack objects with v sin i < 50 km/s. For the stars with Minit > 35 Msun on the hotter side of the bi-stability jump, the observed and predicted rotational rates agree quite well; for those on the cooler side of the jump, the measured velocities are systematically higher than the predicted ones. In general, the derived extra broadening velocities decrease toward cooler Teff, whilst for later evolutionary phases they appear, at the same v sin i, higher for high-mass stars than for low-mass ones. None of the sample stars shows extra broadening velocities higher than 110 km/s. For the majority of the more massive stars, extra broadening either dominates or is in strong competition with rotation. Conclusions: For OB stars of solar metallicity, extra broadening is important and has to be accounted for in the analysis. When appearing at or close to the zero-age main sequence, most of the single and more massive stars rotate slower than previously thought. Model predictions for the evolution of rotation in hot massive stars may need to be updated.

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

TYC 3159-6-1: a runaway blue supergiant

We report the results of optical spectroscopy of a candidate evolved massive star in the Cygnus X region, TYC 3159-6-1, revealed via detection of its curious circumstellar nebula in archival data of the Spitzer Space Telecope. We classify TYC 3159-6-1 as an O9.5-O9.7 Ib star and derive its fundamental parameters by using the stellar atmosphere code FASTWIND. The He and CNO abundances in the photosphere of TYC 3159-6-1 are consistent with the solar abundances, suggesting that the star only recently evolved off the main sequence. Proper motion and radial velocity measurements for TYC 3159-6-1 show that it is a runaway star. We propose that Dolidze 7 is its parent cluster. We discuss the origin of the nebula around TYC 3159-6-1 and suggest that it might be produced in several successive episodes of enhanced mass-loss rate (outbursts) caused by rotation of the star near the critical, Ω-limit.

preprint2013arXiv

Ultra-stripped Type Ic supernovae from close binary evolution

Recent discoveries of weak and fast optical transients raise the question of their origin. We investigate the minimum ejecta mass associated with core-collapse supernovae (SNe) of Type Ic. We show that mass transfer from a helium star to a compact companion can produce an ultra-stripped core which undergoes iron core collapse and leads to an extremely fast and faint SN Ic. In this Letter, a detailed example is presented in which the pre-SN stellar mass is barely above the Chandrasekhar limit, resulting in the ejection of only ~0.05-0.20 M_sun of material and the formation of a low-mass neutron star. We compute synthetic light curves of this case and demonstrate that SN 2005ek could be explained by our model. We estimate that the fraction of such ultra-stripped to all SNe could be as high as 0.001-0.01. Finally, we argue that the second explosion in some double neutron star systems (for example, the double pulsar PSR J0737-3039B) was likely associated with an ultra-stripped SN Ic.

preprint2012arXiv

3D Simulations of Betelgeuse's Bow Shock

Betelgeuse, the bright, cool red supergiant in Orion, is moving supersonically relative to the local interstellar medium. The star emits a powerful stellar wind which collides with this medium, forming a cometary structure, a bow shock, pointing in the direction of motion. We present the first 3D hydrodynamic simulations of the formation and evolution of Betelgeuse's bow shock. The models include realistic low temperature cooling and cover a range of plausible interstellar medium densities and stellar velocities between 0.3 - 1.9 cm-3 and 28 - 73 km/s. We show that the flow dynamics and morphology of the bow shock differ substantially due to the preferential growth of Rayleigh-Taylor or Kelvin-Helmholtz instabilities in the models. The former dominate the models with slow stellar velocities resulting in a clumpy bow shock sub-structure, whereas the latter produce a smoother, more layered sub-structure in the fast models. If the mass in the bow shock shell is low, as seems to be implied by the AKARI luminosities (~0.003 Msun), then Betelgeuse's bow shock is very young and is unlikely to have reached a steady state. The circular nature of the bow shock shell is consistent with this conclusion. Thus, our results suggest that Betelgeuse entered the red supergiant phase only recently.

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

Discovery of two new Galactic candidate luminous blue variables with WISE

We report the discovery of two new Galactic candidate luminous blue variable (cLBV) stars via detection of circular shells (typical of known confirmed and cLBVs) and follow-up spectroscopy of their central stars. The shells were detected at 22 um in the archival data of the Mid-Infrared All Sky Survey carried out with the Wide-field Infrared Survey Explorer (WISE). Follow-up optical spectroscopy of the central stars of the shells conducted with the renewed Southern African Large Telescope (SALT) showed that their spectra are very similar to those of the well-known LBVs P Cygni and AG Car, and the recently discovered cLBV MN112, which implies the LBV classification for these stars as well. The LBV classification of both stars is supported by detection of their significant photometric variability: one of them brightened in the R- and I-bands by 0.68\pm0.10 mag and 0.61\pm0.04 mag, respectively, during the last 13-18 years, while the second one (known as Hen 3-1383) varies its B,V,R,I and K_s brightnesses by \simeq 0.5-0.9 mag on time-scales from 10 days to decades. We also found significant changes in the spectrum of Hen 3-1383 on a timescale of \simeq 3 months, which provides additional support for the LBV classification of this star. Further spectrophotometric monitoring of both stars is required to firmly prove their LBV status. We discuss a connection between the location of massive stars in the field and their fast rotation, and suggest that the LBV activity of the newly discovered cLBVs might be directly related to their possible runaway status.

preprint2012arXiv

Evolution of massive population III stars with rotation and magnetic fields

[Abridged] We present a new grid of massive population III star models including the effects of rotation on the stellar structure and chemical mixing, and magnetic torques for the transport of angular momentum. Based on the grid, we also present a phase diagram for the expected final fates of rotating massive Pop III stars. Our non-rotating models become redder than the previous models in the literature, given the larger overshooting parameter adopted in this study. In particular, convective dredge-up of the helium core material into the hydrogen envelope is observed in our non-rotating very massive star models (>~200 Msun), which is potentially important for the chemical yields. On the other hand, the stars become bluer and more luminous with a higher rotational velocity. With the Spruit-Tayler dynamo, our models with a sufficiently high initial rotational velocity can reach the critical rotation earlier and lose more mass as a result, compared to the previous models without magnetic fields. The most dramatic effect of rotation is found with the so-called chemically homogeneous evolution (CHE), which is observed for a limited mass and rotational velocity range. CHE has several important consequences: 1) Both primary nitrogen and ionizing photons are abundantly produced. 2) Conditions for gamma-ray burst progenitors are fulfilled for an initial mass range of 13 - 84 Msun. 3) Pair instability supernovae of type Ibc are expected for 84 -190 Msun and 4) Both a pulsational pair instability supernova and a GRB may occur from the same progenitor of about 56 - 84 Msun, which might significantly influence the consequent GRB afterglow. We find that CHE does not occur for very massive stars (> 190 Msun), in which case the hydrogen envelope expands to the red-supergiant phase and the final angular momentum is too low to make any explosive event powered by rotation.

preprint2012arXiv

Mass-loss predictions for evolved very metal-poor massive stars

(shortened) The first couple of stellar generations may have been massive, of order 100 Msun, and to have played a dominant role in galaxy formation and the chemical enrichment of the early Universe. Some fraction of these objects may have died as pair-instability supernovae or gamma-ray bursts. The winds if these stars may have played an important role in determining these outcomes. As the winds are driven by radiation pressure on spectral lines, their strengths are expected to vary with metallicity. Until now, most mass-loss predictions for metal-poor O-type stars have assumed a scaled-down solar-abundance pattern. However, Population III evolutionary tracks show significant surface enrichment through rotational mixing of CNO-processed material, because even metal-poor stars switch to CNO-burning early on. We address the question of whether the CNO surface enhanced self-enrichment in the first few generations of stars could impact their mass-loss properties. For this, we employ Monte Carlo simulations to establish the local line-force and solve for the momentum equation of the stellar outflow, testing whether an outflow can actually be established by assessing the net acceleration at the sonic point of the flow. Stellar evolution models of rotating metal-poor stars are used to specify the surface chemical composition, focussing on the phases of early enrichment. We find that the mass-loss rates of CNO enhanced metal-poor stars are higher than those of non-enriched stars, but they are much lower than those rates where the CNO abundance is included in the total abundance Z. We present a heuristic formula that provides mass-loss estimates for CNO-dominated winds in relation to scaled-down solar abundances.

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

Nitrogen chronology of massive main sequence stars

Rotational mixing in massive main sequence stars is predicted to monotonically increase their surface nitrogen abundance with time. We use this effect to design a method for constraining the age and the inclination angle of massive main sequence stars, given their observed luminosity, effective temperature, projected rotational velocity and surface nitrogen abundance. This method relies on stellar evolution models for different metallicities, masses and rotation rates. We use the population synthesis code STARMAKER to show the range of applicability of our method. We apply this method to 79 early B-type main sequence stars near the LMC clusters NGC 2004 and N 11 and the SMC clusters NGC 330 and NGC 346. From all stars within the sample, 17 were found to be suitable for an age analysis. For ten of them, which are rapidly rotating stars without a strong nitrogen enhancement, it has been previously concluded that they did not evolve as rotationally mixed single stars. This is confirmed by our analysis, which flags the age of these objects as highly discrepant with their isochrone ages. For the other seven stars, their nitrogen and isochrone ages are found to agree within error bars, what validates our method. Constraints on the inclination angle have been derived for the other 62 stars,with the implication that the nitrogen abundances of the SMC stars, for which mostly only upper limits are known, fall on average significantly below those limits. Nitrogen chronology is found to be a new useful tool for testing stellar evolution and to constrain fundamental properties of massive main sequence stars. A web version of this tool is provided.

preprint2012arXiv

Pre-Supernova Evolution of Massive Single and Binary Stars

Massive stars are essential to understand a variety of branches of astronomy including galaxy and star cluster evolution, nucleosynthesis and supernovae, pulsars and black holes. It has become evident that massive star evolution is very diverse, being sensitive to metallicity, binarity, rotation, and possibly magnetic fields. While the problem to obtain a good statistical observational database is alleviated by current large spectroscopic surveys, it remains a challenge to model these diverse paths of massive stars towards their violent end stage. We show that the main sequence stage offers the best opportunity to gauge the relevance of the various possible evolutionary scenarios. This also allows to sketch the post-main sequence evolution of massive stars, for which observations of Wolf-Rayet stars give essential clues. Recent supernova discoveries due to the current boost in transient searches allow tentative mappings of progenitor models with supernova types, including pair instability supernovae and gamma-ray bursts.

preprint2012arXiv

Recycling Pulsars: spins, masses and ages

Although the first millisecond pulsars (MSPs) were discovered 30 years ago we still do not understand all details of their formation process. Here, we present new results from Tauris, Langer & Kramer (2012) on the recycling scenario leading to radio MSPs with helium or carbon-oxygen white dwarf companions via evolution of low- and intermediate mass X-ray binaries (LMXBs, IMXBs). We discuss the location of the spin-up line in the (P,Pdot)-diagram and estimate the amount of accreted mass needed to obtain a given spin period and compare with observations. Finally, we constrain the true ages of observed recycled pulsars via calculated isochrones in the (P,Pdot)-diagram.

preprint2012arXiv

Rotating Wolf-Rayet stars in a post RSG/LBV phase. An evolutionary channel towards long-duration GRBs?

(shortened) We investigate the properties of Galactic WR stars and their environment to identify evolutionary channels that may lead to the formation of LGRBs. To this purpose we compile available information on the spectropolarimetric properties, the presence of CS ejecta, and the CS velocities in the environment of Galactic WR stars. We use linear line-depolarization as an indicator of rotation, nebular morphology as an indicator of stellar ejecta, and velocity patterns in UV absorption features as an indicator of increased velocities in the CS environment. We find that the ~23% WR stars with "possible ejecta nebulae" dominate the population of WR stars with spectropolarimetric signatures of rotation, while WR stars without such nebulae only rarely show indications of rotation. The corresponding objects are most likely in an early stage after a preceding RSG or LBV phase, and have not yet lost their angular momenta due to the strong mass-loss in the WR phase. From their photometric periods we estimate rotation parameters in the range omega = 0.04...0.25, corresponding to moderate rotation speeds of 36...120 km/s. These values are very uncertain, but comply with the specific surface angular momentum requirement for LGRB progenitors. Our results indicate that, in the Galaxy, mainly "young" WR stars shortly after a RSG/LBV phase show spectropolarimetric signatures of rotation. Their rotation rates are thus likely enhanced with respect to the majority of Galactic WR stars. According to their estimated specific surface angular momenta, a subgroup of stars exploding in this phase may represent an evolutionary channel towards LGRBs at high metallicities, comparable to the Galaxy.

preprint2012arXiv

The evolution and masses of the neutron star and donor star in the high mass X-ray binary OAO 1657-415

We report near-infrared radial velocity measurements of the recently identified donor star in the high mass X-ray binary system OAO 1657-415 obtained in the H band using ISAAC on the VLT. Cross-correlation methods were employed to construct a radial velocity curve with a semi-amplitude of $22.1 \pm 3.5$ km s$^{-1}$. Combined with other measured parameters of this system this provides a dynamically determined neutron star mass of $1.42 \pm 0.26 M$_{\odot}$ and a mass of $14.3 \pm 0.8 M$_{\odot}$ for the Ofpe/WN9 highly evolved donor star. OAO 1657-415 is an eclipsing High Mass X-ray binary pulsar with the largest eccentricity and orbital period of any within its class. Of the ten known eclipsing X-ray binary pulsars OAO 1657-415 becomes the ninth with a dynamically determined neutron star mass solution and only the second in an eccentric system. Furthermore, the donor star in OAO 1657-415 is much more highly evolved than the majority of the supergiant donors in other High Mass X-ray binaries (HMXBs), joining a small but growing list of HMXBs donors with extensive hydrogen depleted atmospheres. Considering the evolutionary development of OAO 1657-415, we have estimated the binding energy of the envelope of the mass donor and find that there is insufficient energy for the removal of the donor's envelope via spiral-in, ruling out a Common Envelope evolutionary scenario. With its non-zero eccentricity and relatively large orbital period the identification of a definitive evolutionary pathway for OAO1657-415 remains problematic, we conclude by proposing two scenarios which may account for OAO1657-415 current orbital configuration.

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

Zeta Oph and the weak-wind problem

Mass-loss rate, $\dot{M}$, is one of the key parameters affecting evolution and observational manifestations of massive stars, and their impact on the ambient medium. Despite its importance, there is a factor of ~100 discrepancy between empirical and theoretical $\dot{M}$ of late-type O dwarfs, the so-called weak-wind problem. In this Letter, we propose a simple novel method to constrain $\dot{M}$ of runaway massive stars through observation of their bow shocks and Stromgren spheres, which might be of decisive importance for resolving the weak-wind problem. Using this method, we found that $\dot{M}$ of the well-known runaway O9.5 V star zeta Oph is more than an order of magnitude higher than that derived from ultraviolet (UV) line-fitting (Marcolino et al. 2009) and is by a factor of 6 to 7 lower than those based on the theoretical recipe by Vink et al. (2000) and the Halpha line (Mokiem et al. 2005). The discrepancy between $\dot{M}$ derived by our method and that based on UV lines would be even more severe if the stellar wind is clumpy. At the same time, our estimate of $\dot{M}$ agrees with that predicted by the moving reversing layer theory by Lucy (2010).

preprint2011arXiv

A VLT/FLAMES survey for massive binaries in Westerlund 1. III. The WC9d binary W239 and implications for massive stellar evolution

There is growing evidence that a treatment of binarity amongst OB stars is essential for a full theory of stellar evolution. However the binary properties of massive stars - frequency, mass ratio and orbital separation - are still poorly constrained. In order to address this shortcoming we have undertaken a multiepoch spectroscopic study of the stellar population of the young massive cluster Westerlund 1. In this paper we present an investigation into the nature of the dusty Wolf-Rayet star and candidate binary W239. We determined an orbital period of only 5.05 days - making it one of the most compact WC binaries yet identified - while we interpret a significant near-IR flare between 2004-6 as evidence for a third massive stellar component in the system in a wide eccentric orbit. The presence of a near-IR excess characteristic of hot dust at every epoch is consistent with the expectation that the subset of persistent dust forming WC stars are short period binaries. In conjunction with previous results, the simultaneous presence of both short period Wolf-Rayet binaries and cool hypergiants within Wd 1 provides compelling evidence for a bifurcation in the post-Main Sequence evolution of massive stars due to binarity. Short period O+OB binaries will evolve directly to the Wolf-Rayet phase, either due to an episode of binary mediated mass loss or via chemically homogenous evolution. Conversely, long period binaries and single stars will instead undergo a red loop across the HR diagram via a cool hypergiant phase. Future analysis of the full spectroscopic dataset for Wd 1 will constrain the proportion of massive stars experiencing each pathway; hence quantifying the importance of binarity in massive stellar evolution up to and beyond supernova and the resultant production of relativistic remnants and X-ray binaries. [ABRIDGED]

preprint2011arXiv

Exploring the origin of magnetic fields in massive stars: A survey of O-type stars in clusters and in the field

To investigate statistically whether magnetic fields in massive stars are ubiquitous or appear in stars with specific spectral classification, certain ages, or in a special environment, we acquired 41 new spectropolarimetric observations for 36 stars. Among the observed sample roughly half of the stars are probable members of clusters at different ages, whereas the remaining stars are field stars not known to belong to any cluster or association. Spectropolarimetric observations were obtained during three different nights using the low-resolution spectropolarimetric mode of FORS2 (FOcal Reducer low dispersion Spectrograph) mounted on the 8-m Antu telescope of the VLT. To assess the membership in open clusters and associations, we used astrometric catalogues with the best currently available kinematic and photometric data. A field at a significance level of 3sigma was detected in ten O-type stars. Importantly, the largest longitudinal magnetic fields were measured in two Of?p stars: <B_z>=-381+-122G for CPD-282561 and <B_z>=-297+-62G for HD148937, previously detected by us as magnetic. The obtained observations of HD148937 on three different nights indicate that the magnetic field is slightly variable. Our new measurements support our previous conclusion that large-scale organized magnetic fields with polar field strengths in excess of 1kG are not widespread among O-type stars. Among the stars with a detected magnetic field, only one star, HD156154, belongs to an open cluster at high membership probability. According to previous kinematic studies, four magnetic O-type stars in the sample are well-known candidate runaway stars.

preprint2011arXiv

GRIPS - Gamma-Ray Imaging, Polarimetry and Spectroscopy

We propose to perform a continuously scanning all-sky survey from 200 keV to 80 MeV achieving a sensitivity which is better by a factor of 40 or more compared to the previous missions in this energy range. The Gamma-Ray Imaging, Polarimetry and Spectroscopy (GRIPS) mission addresses fundamental questions in ESA's Cosmic Vision plan. Among the major themes of the strategic plan, GRIPS has its focus on the evolving, violent Universe, exploring a unique energy window. We propose to investigate $γ$-ray bursts and blazars, the mechanisms behind supernova explosions, nucleosynthesis and spallation, the enigmatic origin of positrons in our Galaxy, and the nature of radiation processes and particle acceleration in extreme cosmic sources including pulsars and magnetars. The natural energy scale for these non-thermal processes is of the order of MeV. Although they can be partially and indirectly studied using other methods, only the proposed GRIPS measurements will provide direct access to their primary photons. GRIPS will be a driver for the study of transient sources in the era of neutrino and gravitational wave observatories such as IceCUBE and LISA, establishing a new type of diagnostics in relativistic and nuclear astrophysics. This will support extrapolations to investigate star formation, galaxy evolution, and black hole formation at high redshifts.

preprint2011arXiv

Magnetic field studies of massive main sequence stars

We report on the status of our spectropolarimetric observations of massive stars. During the last years, we have discovered magnetic fields in many objects of the upper main sequence, including Be stars, beta Cephei and Slowly Pulsating B stars, and a dozen O stars. Since the effects of those magnetic fields have been found to be substantial by recent models, we are looking into their impact on stellar rotation, pulsation, stellar winds, and chemical abundances. Accurate studies of the age, environment, and kinematic characteristics of the magnetic stars are also promising to give us new insight into the origin of the magnetic fields. Furthermore, longer time series of magnetic field measurements allow us to observe the temporal variability of the magnetic field and to deduce the stellar rotation period and the magnetic field geometry. Studies of the magnetic field in massive stars are indispensable to understand the conditions controlling the presence of those fields and their implications on the stellar physical parameters and evolution.

preprint2011arXiv

Predictions for mass-loss rates and terminal wind velocities of massive O-type stars

Mass loss forms an important aspect of the evolution of massive stars, as well as for the enrichment of the surrounding ISM. Our goal is to predict accurate mass-loss rates and terminal wind velocities. These quantities can be compared to empirical values, thereby testing radiation-driven wind models. One specific issue is that of the "weak-wind problem", where empirically derived mass-loss rates fall orders of magnitude short of predicted values. We employ an established Monte Carlo model and a recently suggested new line acceleration formalism to solve the wind dynamics consistently. We provide a new grid of mass-loss rates and terminal wind velocities of O stars, and compare the values to empirical results. Our models fail to provide mass-loss rates for main-sequence stars below a luminosity of log(L/Lsun) = 5.2, where we run into a fundamental limit. At luminosities below this critical value there is insufficient momentum transferred in the region below the sonic point to kick-start the acceleration. This problem occurs at the location of the onset of the weak-wind problem. For O dwarfs, the boundary between being able to start a wind, and failing to do so, is at spectral type O6/O6.5. The direct cause of this failure is a combination of the lower luminosity and a lack of Fe V lines at the wind base. This might indicate that another mechanism is required to provide the necessary driving to initiate the wind. For stars more luminous than log(L/Lsun) = 5.2, our new mass-loss rates are in excellent agreement with the mass-loss prescription by Vink et al. 2000. This implies that the main assumption entering the method of the Vink et al. prescriptions - i.e. that the momentum equation is not explicitly solved for - does not compromise the reliability of the Vink et al. results for this part of parameter space (Abridged).

preprint2011arXiv

Resolved photometry of extragalactic young massive star clusters

We present colour-magnitude diagrams (CMDs) for a sample of seven young massive clusters in the galaxies NGC 1313, NGC 1569, NGC 1705, NGC 5236 and NGC 7793. The clusters have ages in the range 5-50 million years and masses of 10^5 -10^6 Msun. Although crowding prevents us from obtaining photometry in the central regions of the clusters, we are still able to measure up to 30-100 supergiant stars in each of the richest clusters, along with the brighter main sequence stars. The resulting CMDs and luminosity functions are compared with photometry of artificially generated clusters, designed to reproduce the photometric errors and completeness as realistically as possible. In agreement with previous studies, our CMDs show no clear gap between the H-burning main sequence and the He-burning supergiant stars, contrary to predictions by common stellar isochrones. In general, the isochrones also fail to match the observed number ratios of red-to-blue supergiant stars, although the difficulty of separating blue supergiants from the main sequence complicates this comparison. In several cases we observe a large spread (1-2 mag) in the luminosities of the supergiant stars that cannot be accounted for by observational errors. This spread can be reproduced by including an age spread of 10-30 million years in the models. However, age spreads cannot fully account for the observed morphology of the CMDs and other processes, such as the evolution of interacting binary stars, may also play a role.

preprint2011arXiv

The Eddington factor as the key to understand the winds of the most massive stars. Evidence for a Gamma-dependence of Wolf-Rayet type mass loss

The most massive stars are thought to be hydrogen-rich Wolf-Rayet stars of late spectral subtype (WNh stars). In previous theoretical studies the enhanced mass loss of these stars has been attributed to their proximity to the Eddington limit. Here we investigate observed trends in the mass-loss properties of such young, very massive stars. We derive theoretical mass-luminosity relations for very massive stars, based on a large grid of stellar structure models. Using these relations, we estimate Eddington factors for a sample of stars, under different assumptions of their evolutionary status. We evaluate the resulting mass-loss relations, and compare them with theoretical predictions. We find observational evidence that the mass loss in the WR regime is dominated by the Eddington parameter Gamma_e, which has important consequences for the way we understand Wolf-Rayet stars and their mass loss. In addition, we derive wind masses that support the picture that the WNh stars in young stellar clusters are very massive, hydrogen-burning stars. Our findings suggest that the proximity to the Eddington limit is the physical reason for the onset of Wolf-Rayet type mass loss. This means that, e.g. in stellar evolution models, the Wolf-Rayet stage should be identified by large Eddington parameters, instead of a helium-enriched surface composition. The latter is most likely only a consequence of strong mass loss, in combination with internal mixing. For very massive stars, the enhanced Gamma-dependent mass loss is responsible for the formation of late WNh subtypes with high hydrogen surface abundances, partly close to solar. Because mass loss dominates the evolution of very massive stars, we expect a strong impact of this effect on their end products, in particular on the potential formation of black holes, and Gamma-Ray Bursts, as well as the observed upper mass limit of stars.

preprint2011arXiv

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

NGC346-013 is a peculiar double-lined eclipsing binary in the Small Magellanic Cloud discovered by the VLT-FLAMES survey of massive stars. Spectra obtained with VLT-FLAMES are used to construct a radial velocity curve and photometry obtained with the Faulkes Telescope South is then used to derive orbital parameters, while spectra of the secondary are compared with synthetic spectra from TLUSTY model atmospheres. The orbital period is found to be 4.20381(12) days, with masses of 19.1+/-1.0 and 11.9+/-0.6 Msun. The primary is a rapidly rotating late-O dwarf while the secondary, an early-B giant, displays near-synchronous rotation and has filled its Roche lobe, implying that it was originally the more massive component with recent mass transfer `spinning up' the primary to near-critical rotation. Comparison with synthetic spectra finds temperatures of 34.5kK and 24.5kK for the primary and secondary respectively, with the nitrogen abundance of the secondary enhanced compared to baseline values for the SMC, consistent with the predictions of models of interacting binaries. NGC346-013 likely evolved via non-conservative mass transfer in a system with initial masses ~22+15Msun, with the well-constrained orbital solution and atmospheric parameters making it an excellent candidate for tailored modelling with binary evolution codes. This system will form a cornerstone in constraining the physics of thermal timescale mass transfer, and the associated mass transfer efficiency, in massive close binary systems.

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.

preprint2011arXiv

Wind modelling of very massive stars up to 300 solar masses

Some studies have claimed a universal stellar upper-mass limit of 150 Msun. A factor that is often overlooked is that there might be a difference between the current and initial masses of the most massive stars, as a result of mass loss. We present Monte Carlo mass-loss predictions for very massive stars in the range 40-300 Msun, with large luminosities and Eddington factors Gamma. Using our new dynamical approach, we find an upturn in the mass-loss vs. Gamma dependence, at the point where the winds become optically thick. This coincides with the location where wind efficiency numbers surpass the single-scattering limit of Eta = 1, reaching values up to Eta = 2.5. Our modelling suggests a transition from common O-type winds to Wolf-Rayet characteristics at the point where the winds become optically thick. This transitional behaviour is also revealed with respect to the wind acceleration parameter beta, which starts at values below 1 for the optically thin O-stars, and naturally reaches values as high as 1.5-2 for the optically thick Wolf-Rayet models. An additional finding concerns the transition in spectral morphology of the Of and WN characteristic He II line at 4686 Angstrom. When we express our mass-loss predictions as a function of the electron scattering Gamma_e (=L/M) only, we obtain a mass-loss Gamma dependence that is consistent with a previously reported power-law Mdot propto Gamma^5 (Vink 2006) that was based on our semi-empirical modelling approach. When we express Mdot in terms of both Gamma and stellar mass, we find Mdot propto M^0.8 Gamma^4.8 for our high Gamma models. Finally, we confirm that the Gamma-effect on the mass-loss predictions is much stronger than that of an increased helium abundance, calling for a fundamental revision in the way mass loss is incorporated in evolutionary models of the most massive stars.

preprint2011arXiv

Wind Models for Very Massive Stars in the Local Universe

Some studies have claimed the existence of a stellar upper-mass limit of 150 Msun. A factor that is often overlooked concerns the issue that there might be a significant difference between the present-day and the initial mass of the most massive stars - as a result of mass loss. The upper-mass limit may be substantially higher, possibly exceeding 200 Msun. The issue of the upper mass-limit will however remain uncertain as long as there is only limited quantitative knowledge of mass loss in close proximity to the Eddington (= Gamma) limit. For this reason, we present mass-loss predictions from Monte Carlo radiative transfer models for very massive stars up to 300 Msun. Using our new dynamical approach, we find an upturn or "kink" in the mass-loss versus Gamma dependence, at the point where our model winds become optically thick. These are the first mass-loss predictions where the transition from optically thin O-star winds to optically thick Wolf-Rayet winds has been resolved.

preprint2010arXiv

A VLT/FLAMES survey for massive binaries in Westerlund 1: II. Dynamical constraints on magnetar progenitor masses from the eclipsing binary W13

Westerlund 1 is a young, massive Galactic starburst cluster that contains a rich coeval population of Wolf-Rayet stars, hot- and cool-phase transitional supergiants, and a magnetar. We use spectroscopic and photometric observations of the eclipsing double-lined binary W13 to derive dynamical masses for the two components, in order to determine limits for the progenitor masses of the magnetar CXOU J164710.2-455216 and the population of evolved stars in Wd1. W13 has an orbital period of 9.2709+/-0.0015 days and near-contact configuration. The shallow photometric eclipse rules out an inclination greater than 65 degrees, leading to lower limits for the masses of the emission-line optical primary and supergiant optical secondary of 21.4+/-2.6Msun and 32.8+/-4.0Msun respectively, rising to 23.2 +3.3/-3.0Msun and 35.4 +5.0/-4.6 Msun for our best-fit inclination 62 +3/-4 degrees. Comparison with theoretical models of Wolf-Rayet binary evolution suggest the emission-line object had an initial mass in excess of 35Msun, with the most likely model featuring highly non-conservative late-Case-A/Case-B mass transfer and an initial mass in excess of 40Msun. This confirms the high magnetar progenitor mass inferred from its membership in Wd1, and represents the first dynamical constraint on the progenitor mass of any magnetar. The red supergiants in Wd1 must have similar progenitor masses to W13 and are therefore amongst the most massive stars to undergo a red supergiant phase, representing a challenge for population models that suggest stars in this mass range end their redwards evolution as yellow hypergiants. [ABRIDGED]

preprint2010arXiv

Binaries are the best single stars

Stellar models of massive single stars are still plagued by major uncertainties. Testing and calibrating against observations is essential for their reliability. For this purpose one preferably uses observed stars that have never experienced strong binary interaction, i.e. "true single stars". However, the binary fraction among massive stars is high and identifying "true single stars" is not straight forward. Binary interaction affects systems in such a way that the initially less massive star becomes, or appears to be, single. For example, mass transfer results in a widening of the orbit and a decrease of the luminosity of the donor star, which makes it very hard to detect. After a merger or disruption of the system by the supernova explosion, no companion will be present. The only unambiguous identification of "true single stars" is possible in detached binaries, which contain two main-sequence stars. For these systems we can exclude the occurrence of mass transfer since their birth. A further advantage is that binaries can often provide us with direct measurements of the fundamental stellar parameters. Therefore, we argue these binaries are worth the effort needed to observe and analyze them. They may provide the most stringent test cases for single stellar models.

preprint2010arXiv

Chemically homogeneous evolution in massive binaries

Rotation can have severe consequences for the evolution of massive stars. It is now considered as one of the main parameters, alongside mass and metallicity that determine the final fate of single stars. In massive, fast rotating stars mixing processes induced by rotation may be so efficient that helium produced in the center is mixed throughout the envelope. Such stars evolve almost chemically homogeneously. At low metallicity they remain blue and compact, while they gradually evolve into Wolf-Rayet stars and possibly into progenitors of long gamma-ray bursts. In binaries this type of evolution may occur because of (I) tides in very close binaries, as a result of (II) spin up by mass transfer, as result of (III) a merger of the two stars and (IV) when one of the components in the binary was born with a very high initial rotation rate. As these stars stay compact, the evolutionary channels are very different from what classical binary evolutionary models predict. In this contribution we discuss examples of nearly chemically homogeneous evolution in very close tidally-locked binaries. Even in such very close massive binaries, the stars may remain compact and avoid mass transfer, while Roche lobe overflow and a merger would be inevitable in the classical picture. This type of evolution may provide an alternative path to form tight Wolf-Rayet binaries and massive black hole binaries.

preprint2010arXiv

Fast rotating stars resulting from binary evolution will often appear to be single

Rapidly rotating stars are readily produced in binary systems. An accreting star in a binary system can be spun up by mass accretion and quickly approach the break-up limit. Mergers between two stars in a binary are expected to result in massive, fast rotating stars. These rapid rotators may appear as Be or Oe stars or at low metallicity they may be progenitors of long gamma-ray bursts. Given the high frequency of massive stars in close binaries it seems likely that a large fraction of rapidly rotating stars result from binary interaction. It is not straightforward to distinguish a a fast rotator that was born as a rapidly rotating single star from a fast rotator that resulted from some kind of binary interaction. Rapidly rotating stars resulting from binary interaction will often appear to be single because the companion tends to be a low mass, low luminosity star in a wide orbit. Alternatively, they became single stars after a merger or disruption of the binary system during the supernova explosion of the primary. The absence of evidence for a companion does not guarantee that the system did not experience binary interaction in the past. If binary interaction is one of the main causes of high stellar rotation rates, the binary fraction is expected to be smaller among fast rotators. How this prediction depend on uncertainties in the physics of the binary interactions requires further investigation.

preprint2010arXiv

Light elements in massive single and binary stars

We highlight the role of the light elements (Li, Be, B) in the evolution of massive single and binary stars, which is largely restricted to a diagnostic value, and foremost so for the element boron. However, we show that the boron surface abundance in massive early type stars contains key information about their foregoing evolution which is not obtainable otherwise. In particular, it allows to constrain internal mixing processes and potential previous mass transfer event for binary stars (even if the companion has disappeared). It may also help solving the mystery of the slowly rotating nitrogen-rich massive main sequence stars.

preprint2010arXiv

The nature of B supergiants: clues from a steep drop in rotation rates at 22000 K. The possibility of Bi-stability braking

The location of B supergiants in the Hertzsprung-Russell diagram (HRD) represents a long-standing problem in massive star evolution. Here we propose their nature may be revealed utilising their rotational properties, and we highlight a steep drop in massive star rotation rates at an effective temperature of 22000 K. We discuss two potential explanations for it. On the one hand, the feature might be due to the end of the main sequence, which could potentially constrain the core overshooting parameter. On the other hand, the feature might be the result of enhanced mass loss at the predicted location of the bi-stability jump. We term this effect "bi-stability breaking" and discuss its potential consequences for the evolution of massive stars.

preprint2009arXiv

Massive binaries and the enrichment of the interstellar medium in globular clusters

Abundance anomalies observed in globular cluster stars indicate pollution with material processed by hydrogen burning. Two main sources have been suggested: asymptotic giant branch stars and massive stars rotating near the break-up limit. We discuss the potential of massive binaries as an interesting alternative source of processed material. We discuss observational evidence for mass shedding from interacting binaries. In contrast to the fast, radiatively driven winds of massive stars, this material is typically ejected with low velocity. We expect that it remains inside the potential well of a globular cluster and becomes available for the formation or pollution of a second generation of stars. We estimate that the amount of processed low-velocity material that can be ejected by massive binaries is larger than the contribution of two previously suggested sources combined.

preprint2009arXiv

Massive binaries as the source of abundance anomalies in globular clusters

Abundance anomalies observed in globular cluster stars indicate pollution with material processed by hydrogen burning. Two main sources have been suggested: asymptotic giant branch (AGB) stars and massive stars rotating near the break-up limit (spin stars). We propose massive binaries as an alternative source. We compute the evolution of a 20 Msun star in a close binary and find that it sheds about 10 Msun. The ejecta are enriched in He, N, Na, and Al and depleted in C and O, similar to the abundance patterns observed in gobular cluster stars. In contrast to the fast, radiatively driven winds of massive stars, this material is typically ejected with low velocity. We expect that it remains inside the potential well of a globular cluster and becomes available for the formation or pollution of a second generation of stars. We estimate that the amount of processed low-velocity material ejected by massive binaries is greater than the contribution of AGB stars and spin stars combined, assuming that the majority of massive stars in a proto-globular cluster interact. If we take the possible contribution of intermediate mass stars in binaries into account and assume that the ejecta are diluted with an equal amount of unprocessed material, we find that this scenario can provide enough material to form a second generation of low-mass stars, which is as numerous as the first generation of low-mass stars. In this scenario there is no need to make commonly adopted assumptions, such as preferential loss of the first generation of stars, external pollution of the cluster, or an anomalous initial mass function. [Abridged]

preprint2009arXiv

Sub-surface convection zones in hot massive stars and their observable consequences

We study the convection zones in the outer envelope of hot massive stars which are caused by opacity peaks associated with iron and helium ionization. We determine the occurrence and properties of these convection zones as function of the stellar parameters. We then confront our results with observations of OB stars. A stellar evolution code is used to compute a grid of massive star models at different metallicities. In these models, the mixing length theory is used to characterize the envelope convection zones. We find the iron convection zone (FeCZ) to be more prominent for lower surface gravity, higher luminosity and higher initial metallicity. It is absent for luminosities below about $10^{3.2}\Lsun$, $10^{3.9}\Lsun$, and $10^{4.2}$\Lsun$ for the Galaxy, LMC and SMC, respectively. We map the strength of the FeCZ on the Hertzsprung-Russell diagram for three metallicities, and compare this with the occurrence of observational phenomena in O stars: microturbulence, non-radial pulsations, wind clumping, and line profile variability. The confirmation of all three trends for the FeCZ as function of stellar parameters by empirical microturbulent velocities argues for a physical connection between sub-photospheric convective motions and small scale stochastic velocities in the photosphere of O- and B-type stars. We further suggest that clumping in the inner parts of the winds of OB stars could be caused by the same mechanism, and that magnetic fields produced in the FeCZ could appear at the surface of OB stars as diagnosed by discrete absorption components in ultraviolet absorption lines.

preprint2009arXiv

Supernova progenitor stars in the initial range of 23 to 33 solar masses and their relation with the SNR Cas A

Multi wavelength observations of Cassiopeia A (Cas A) have provided us with a strong evidence for the presence of circumstellar material surrounding the progenitor star. It has been suggested that its progenitor was a massive star with a strong mass loss. But, despite the large amount of observational data from optical, IR, radio and x-ray observations, the identity of Cas A progenitor is still elusive. In this work, we compute stellar and circumstellar numerical models to look for the progenitor of Cas A. The models will be compared with the observational constraints. We have computed stellar evolution models to get time-dependent wind parameters and surface abundances. We have chosen a set of probable progenitor stars, with initial masses of 23, 28, 29, 30 and 33 \Mo, with mass loss. The derived mass loss rates and wind terminal velocities are used to simulate the evolution of the circumstellar medium. Our stellar set gives distinct SN progenitors: RSG, luminous blue super giants (LBSGs) and Wolf-Rayet (WR) stars. The 23 \Mo star explodes as a RSG in a $\rm ρ\sim r^{-2}$ dense, free streaming wind surrounded by a thin, compressed, RSG shell. The 28 $\Mo$ star explodes as a LBSG, and the SN blast wave interacts with a low density, free streaming wind surrounded by an unstable and massive ''RSG+LBSG'' shell. Finally, the 30 and 33 \Mo stars explode as WR stars surrounded by fast, WR winds that terminate in highly fragmented ''WR+RSG shell''. We have compared the surface chemical abundances of our stellar models with the observational abundances in Cas A. The abundance analysis shows that the progenitor was a star with an initial mass of the order of 30 \Mo, while the hydrodynamical analysis favors progenitors of initial masses around 23.

preprint2009arXiv

The effect of rotation on the stability of nuclear burning in accreting neutron stars

Hydrogen and/or helium accreted by a neutron star from a binary companion may undergo thermonuclear fusion. At different mass accretion rates different burning regimes are discerned. Theoretical models predict helium fusion to proceed as a thermonuclear runaway for accretion rates below the Eddington limit and as stable burning above this limit. Observations, however, place the boundary close to 10% of the Eddington limit. We study the effect of rotationally induced transport processes on the stability of helium burning. For the first time detailed calculations of thin helium shell burning on neutron stars are performed using a hydrodynamic stellar evolution code including rotation and rotationally induced magnetic fields. We find that in most cases the instabilities from the magnetic field provide the dominant contribution to the chemical mixing, while Eddington-Sweet circulations become important at high rotation rates. As helium is diffused to greater depths, the stability of the burning is increased, such that the critical accretion rate for stable helium burning is found to be lower. Combined with a higher heat flux from the crust, as suggested by recent studies, turbulent mixing could explain the observed critical accretion rate. Furthermore, close to this boundary we find oscillatory burning, which previous studies have linked to mHz QPOs. In models where we continuously lower the heat flux from the crust, the period of the oscillations increases by up to several tens of percents, similar to the observed frequency drift, suggesting that this drift could be caused by the cooling of deeper layers.

preprint2009arXiv

The hydrodynamics of the supernova remnant Cas A: The influence of the progenitor evolution on the velocity structure and clumping

There are large differences in the proposed progenitor models for the Cas A SNR. One of these differences is the presence or absence of a Wolf-Rayet (WR) phase of the progenitor star. The mass loss history of the progenitor star strongly affects the shape of the Supernova remnant (SNR). In this paper we investigate whether the progenitor star of Cas A had a WR phase or not and how long it may have lasted. We performed two-dimensional multi-species hydrodynamical simulations of the CSM around the progenitor star for several WR life times, each followed by the interaction of the supernova ejecta with the CSM. We then looked at the influence of the length of the WR phase and compared the results of the simulations with the observations of Cas A. The difference in the structure of the CSM, for models with different WR life times, has a strong impact on the resulting SNR. With an increasing WR life time the reverse shock velocity of the SNR decreases and the range of observed velocities in the shocked material increases. Furthermore, if a WR phase occurs, the remainders of the WR shell will be visible in the resulting SNR. Comparing our results with the observations suggests that the progenitor star of Cas A did not have a WR phase. We also find that the quasi-stationary flocculi (QSF) in Cas A are not consistent with the clumps from a WR shell that have been shocked and accelerated by the interaction with the SN ejecta. We can also conclude that for a SN explosion taking place in a CSM that is shaped by the wind during a short < 15000 yr WR phase, the clumps from the WR shell will be visible inside the SNR.

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.

preprint2006arXiv

Single star progenitors of long gamma-ray bursts I: Model grids and redshift dependent GRB rate

We present grids of massive star evolution models at four different metallicities (Z=0.004, 0.002, 0.001, 0.00001). The effects of rotation on the stellar structure and the transport of angular momentum and chemical elements through the Spruit-Tayler dynamo and rotationally induced instabilities are considered. After discussing uncertainties involved with the adopted physics, we elaborate the final fate of massive stars as a function of initial mass and spin rate, at each considered metallicity. In particular, we investigate for which initial conditions long gamma-ray bursts (GRBs) are expected to be produced in the frame of the collapsar model. Then, using an empirical spin distribution of young massive metal-poor stars and a specified metallicity-dependent history of star-formation, we compute the expected GRB rate as function of metallicity and redshift based on our stellar evolution models. The GRB production in our models is limited to metallicities of Z \lsim 0.004, with the consequence that about 50 % of all GRBs are predicted to be found at redshifts above z = 4, with most supernovae occurring at redshifts below z\simeq 2.2. The average GRB/SN ratio predicted by our model is about 1/200 globally, and 1/1250 at low redshift. Future strategies for testing the considered GRB progenitor scenario are briefly discussed.

preprint2004arXiv

The first binary star evolution model producing a Chandrasekhar mass white dwarf

Today, Type Ia supernovae are essential tools for cosmology, and recognized as major contributors to the chemical evolution of galaxies. The construction of detailed supernova progenitor models, however, was so far prevented by various physical and numerical difficulties in simulating binary systems with an accreting white dwarf component, e.g., unstable helium shell burning which may cause significant expansion and mass loss. Here, we present the first binary evolution calculation which models both stellar components and the binary interaction simultaneously, and where the white dwarf mass grows up to the Chandrasekhar limit by mass accretion. Our model starts with a 1.6 Msun helium star and a 1.0 Msun CO white dwarf in a 0.124 day orbit. Thermally unstable mass transfer starts when the CO core of the helium star reaches 0.53 Msun, with mass transfer rates of 1...8 times 10^{-6} Msun/yr. The white dwarf burns the accreted helium steadily until the white dwarf mass has reached ~ 1.3 Msun and weak thermal pulses follow until carbon ignites in the center when the white dwarf reaches 1.37 Msun. Although the supernova production rate through this channel is not well known, and this channel can not be the only one as its progenitor life time is rather short (~ 10^7 - 10^8 yr), our results indicate that helium star plus white dwarf systems form a reliable route for producing Type Ia supernovae.

preprint2001arXiv

Formation of High Mass X-ray Black Hole Binaries

The discrepancy in the past years of many more black-hole soft X-ray transients (SXTs), of which a dozen have now been identified, had challenged accepted wisdom in black hole evolution. The sheer number of inferred SXTs requires that many, if not most, stars of ZAMS masses $20-35\msun$ end up as black holes (Portegies Zwart et al. 1997; Ergma & van den Heuvel 1998). In this paper we show that this can be understood by challenging the accepted wisdom that the result of helium core burning in a massive star is independent of whether the core is covered by a hydrogen envelope, or "naked" while it burns. The latter case occurs in binaries when the envelope of the more massive star is transferred to the companion by Roche Lobe overflow while in either main sequence or red giant stage. For solar metallicity, whereas the helium cores which burn while naked essentially never go into high-mass black holes, those that burn while clothed do so, beginning at ZAMS mass $\sim 20\msun$, the precise mass depending on the $^{12}C(α,γ)^{16}$O rate as we outline. In this way the SXTs can be evolved, provided that the H envelope of the massive star is removed only following the He core burning. We also show that in order to evolve a black hole of mass $\gsim 10\msun$ such as observed in Cyg X-1, even employing extremely massive progenitors of ZAMS mass $\gsim 60\msun$ for the black hole, the core must be covered by hydrogen during a substantial fraction of the core burning. In other words, the progenitor must be a WNL star. We evolve Cyg X-1 in an analogous way to which the SXTs are evolved, the difference being that the companion in Cyg X-1 is more massive than those in the SXTs, so that Cyg X-1 shines continuously.

preprint1997arXiv

B[e] supergiants: What is their evolutionary status?

In this paper, we investigate the evolutionary status of B[e]~stars from the point of view of stellar evolution theory. We try to answer to the question of how massive hot supergiants --- i.e. evolved stars --- can be capable of producing a circumstellar disk. We find and discuss three possibilities: very massive evolved main sequence stars close to critical rotation due to their proximity to their Eddington-limit, blue supergiants which have just left the red supergiant branch, and single star merger remnants of a close binary system. While the latter process seems to be required to understand the properties of the spectroscopic binary R4 in the LMC, the other two scenarios may be capable of explaining the distribution of the B[e] stars in the HR~diagram. The three scenarios make different predictions about the duration of the B[e]~phase, the time integrated disk mass and the stellar properties during the B[e]~phase, which may ultimately allow to distinguish them observationally.