Source author record

G. Gräfener

G. Gräfener 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

27works
6topics
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

27 published item(s)

preprint2020arXiv

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

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

preprint2020arXiv

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

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

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.

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.

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

Light-travel-time diagnostics in early Supernova spectra: substantial mass loss of the IIb progenitor of SN 2013cu through a superwind

The progenitors of type-IIb supernovae (SNe) are believed to have lost their H-rich envelopes almost completely in the direct pre-SN phase. Recently the first 'flash spectrum' of a SN IIb (SN2013cu) has been presented, taken early enough to study its immediate circumstellar medium (CSM). Similar to a previous study by Groh (2014) we analyse the structure and chemical composition of the optically-thick CSM using non-LTE model atmospheres. For the first time we take light-travel time (LTT) effects on the spectrum formation into account, which affect the shapes and strengths of the observable emission lines, as well as the inferred SN luminosity. Based on the new CSM parameters we estimate a lower limit of ~0.3Msun for the CSM mass, which is a factor 10-100 higher than previous estimates. The spectral fit implies a CSM in the form of a homogeneous and spherically symmetric superwind whose mass-loss rate exceeds common expectations by up to two orders of magnitude. The derived chemical composition is in agreement with a progenitor that has just left, or is just about to leave the Red-Supergiant (RSG) stage, confirming the standard picture for the origin of SNe IIb. Due to its extreme mass loss the SN progenitor will likely appear as extreme RSG, Luminous Blue Variable (LBV), or Yellow Hypergiant (YHG). The direct detection of a superwind, and the high inferred CSM mass suggest that stellar wind mass loss may play an important role in the formation of SNe IIb.

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

Narrow He II emission in star-forming galaxies at low metallicity. Stellar wind emission from a population of Very Massive Stars

In a recent study star-forming galaxies with HeII emission between redshifts 2 and 4.6 have been found to occur in two modes, distinguished by the width of their HeII emission lines. Broad HeII emission has been attributed to stellar emission from a population of evolved Wolf-Rayet (WR) stars while narrow HeII emission has been attributed to nebular emission excited by a population of very hot PopIII stars formed in pockets of pristine gas at moderate redshifts. In this work we propose an alternative scenario for the origin of the narrow HeII emission, namely very massive stars (VMS) at low metallicity (Z) which form strong but slow WR-type stellar winds due to their proximity to the Eddington limit. We estimate the expected HeII line fluxes and equivalent widths based on wind models for VMS and population synthesis models, and compare the results with recent observations of star-forming galaxies at moderate redshifts. The observed HeII line strengths and equivalent widths are in line with what is expected for a population of VMS in one or more young super-clusters located within these galaxies. In our scenario the two observed modes of HeII emission originate from massive stellar populations in distinct evolutionary stages at low Z. If this interpretation is correct there is no need to postulate the existence of PopIII stars at moderate redshifts to explain the observed narrow HeII emission. An interesting possibility is the existence of self-enriched VMS with similar WR-type spectra at extremely low Z. Stellar HeII emission from such very early generations of VMS may be detectable in future studies of star-forming galaxies at high redshifts with the James Webb Space Telescope. The fact that the HeII emission of VMS is largely neglected in current population synthesis models will generally affect the interpretation of the integrated spectra of young stellar populations.

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

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

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 VLT-FLAMES Tarantula Survey XVI. The optical+NIR extinction laws in 30 Doradus and the photometric determination of the effective temperatures of OB stars

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

preprint2014arXiv

The Wolf-Rayet stars in the Large Magellanic Cloud: A comprehensive analysis of the WN class

Aims: Following our comprehensive studies of the WR stars in the Milky Way, we now present spectroscopic analyses of almost all known WN stars in the LMC. Methods: For the quantitative analysis of the wind-dominated emission-line spectra, we employ the Potsdam Wolf-Rayet (PoWR) model atmosphere code. By fitting synthetic spectra to the observed spectral energy distribution and the available spectra (ultraviolet and optical), we obtain the physical properties of 107 stars. Results: We present the fundamental stellar and wind parameters for an almost complete sample of WN stars in the LMC. Among those stars that are putatively single, two different groups can be clearly distinguished. While 12% of our sample are more luminous than 10^6 Lsun and contain a significant amount of hydrogen, 88% of the WN stars, with little or no hydrogen, populate the luminosity range between log (L/Lsun) = 5.3...5.8. Conclusions: While the few extremely luminous stars (log (L/Lsun) > 6), if indeed single stars, descended directly from the main sequence at very high initial masses, the bulk of WN stars have gone through the red-supergiant phase. According to their luminosities in the range of log (L/Lsun) = 5.3...5.8, these stars originate from initial masses between 20 and 40 Msun. This mass range is similar to the one found in the Galaxy, i.e. the expected metallicity dependence of the evolution is not seen. Current stellar evolution tracks, even when accounting for rotationally induced mixing, still partly fail to reproduce the observed ranges of luminosities and initial masses. Moreover, stellar radii are generally larger and effective temperatures correspondingly lower than predicted from stellar evolution models, probably due to subphotospheric inflation.

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

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

Stellar mass-loss near the Eddington limit. Tracing the sub-photospheric layers of classical Wolf-Rayet stars

Towards the end of their evolution hot massive stars develop strong stellar winds and appear as emission line stars, such as WR stars or LBVs. The quantitative description of the mass loss in these important pre-SN phases is hampered by unknowns such as clumping and porosity due to an in-homogeneous wind structure, and by an incomplete theoretical understanding of optically thick stellar winds. In this work we investigate the conditions in deep atmospheric layers of WR stars to find out whether these comply with the theory of optically thick winds, and whether we find indications of clumping in these layers. We use a new semi-empirical method to determine sonic-point optical depths, densities, and temperatures for a large sample of WR stars of the carbon (WC) and oxygen (WO) sequence. Based on an artificial model sequence we investigate the reliability of our method and its sensitivity to uncertainties in stellar parameters. We find that the WR stars in our sample obey an approximate relation with P_rad/P_gas~80 at the sonic point. This 'wind condition' is ubiquitous for radiatively driven, optically thick winds, and sets constraints on possible wind/envelope solutions affecting radii, mass-loss rates, and clumping properties. Our results suggest that the presence of an optically thick wind may force many stars near the Eddington limit to develop clumped, radially extended sub-surface zones. The clumping in these zones is most likely sustained by the non-linear strange-mode instability, and may be the origin of the observed wind clumping. The properties of typical late-type WC stars comply with this model. Solutions without sub-surface clumping and inflation are also possible but demand for compact stars with comparatively low mass-loss rates. These objects may resemble the small group of WO stars with their exceptionally hot stellar temperatures and highly ionized winds.

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.

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 VLT-FLAMES Tarantula Survey IV: Candidates for isolated high-mass star formation in 30 Doradus

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

preprint2011arXiv

Stellar envelope inflation near the Eddington limit. Implications for the radii of Wolf-Rayet stars and luminous blue variables

(shortened) It has been proposed that the envelopes of luminous stars may be subject to substantial radius inflation. The inflation effect has been discussed in relation to the radius problem of WR stars, but has yet failed to explain the large observed radii of Galactic WR stars. We wish to obtain a physical perspective of the inflation effect, and study the consequences for the radii of WR stars, and LBVs. For WR stars the observed radii are up to an order of magnitude larger than predicted by theory, whilst S Doradus-type LBVs are subject to humongous radius variations, which remain as yet ill-explained. We use a dual approach to investigate the envelope inflation, based on numerical models for stars near the Eddington limit, and a new analytic formalism to describe the effect. An additional new aspect is that we take the effect of density inhomogeneities (clumping) within the outer stellar envelopes into account. Due to the effect of clumping we are able to bring the observed WR radii in agreement with theory. Based on our new formalism, we find that the radial inflation is a function of a dimensionless parameter W, which largely depends on the topology of the Fe-opacity peak, i.e., on material properties. For W>1, we discover an instability limit, for which the stellar envelope becomes gravitationally unbound, i.e. there no longer exists a static solution. Within this framework we are also able to explain the S Doradus-type instabilities for LBVs like AG Car, with a possible triggering due to changes in stellar rotation. The stellar effective temperatures in the upper HR diagram are potentially strongly affected by the inflation effect. This may have particularly strong effects on the evolved massive LBV and WR stars just prior to their final collapse, as the progenitors of SNe Ibc, SNe II, and long GRBs.

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

preprint2010arXiv

The central star of the planetary nebula PB 8: a Wolf-Rayet-type wind of an unusual WN/WC chemical composition

A considerable fraction of the central stars of planetary nebulae (CSPNe) are hydrogen-deficient. As a rule, these CSPNe exhibit a chemical composition of helium, carbon, and oxygen with the majority showing Wolf-Rayet-like emission line spectra. These stars are classified as CSPNe of a spectral type [WC]. We perform a spectral analysis of CSPN PB 8 with the Potsdam Wolf-Rayet (PoWR) models for expanding atmospheres. The source PB 8 displays wind-broadened emission lines from strong mass loss. Most strikingly, we find that its surface composition is hydrogen-deficient, but not carbon-rich. With mass fractions of 55% helium, 40% hydrogen, 1.3% carbon, 2% nitrogen, and 1.3% oxygen, it differs greatly from the 30-50% of carbon which are typically seen in [WC]-type central stars. The atmospheric mixture in PB 8 has an analogy in the WN/WC transition type among the massive Wolf-Rayet stars. Therefore we suggest to introduce a new spectral type [WN/WC] for CSPNe, with PB 8 as its first member. The central star of PB 8 has a relatively low temperature of T=52kK, as expected for central stars in their early evolutionary stages. Its surrounding nebula is less than 3000 years old, i.e. relatively young. Existing calculations for the post-AGB evolution can produce hydrogen-deficient stars of the [WC] type, but do not predict the composition found in PB 8. We discuss various scenarios that might explain the origin of this unique object.

preprint2003arXiv

Comprehensive modelling of the planetary nebula LMC-SMP 61 and its [WC]-type central star

We present a comprehensive study of the Magellanic Cloud planetary nebula SMP 61 and of its nucleus, a Wolf-Rayet type star classified [WC 5-6]. We have performed a detailed spectral analysis of the central star, using the Potsdam code for expanding atmospheres in non-LTE. The fluxes from the model stellar atmosphere were used to compute photoionization models of the nebula. All the available observations, within their error bars, were used to constrain these models. We find that the ionizing fluxes predicted by the stellar model are basically consistent with the fluxes needed by the photoionization model to reproduce the nebular emission, within the error margins. However, there are indications that the stellar model overestimates the number and hardness of Lyman continuum photons. The photoionization models imply a clumped density structure of the nebular material. The observed CIII] 1909/CII 4267 line ratio implies the existence of carbon-rich clumps in the nebula. Such clumps are likely produced by stellar wind ejecta, possibly mixed with the nebular material. We discuss our results with regard to the stellar and nebular post-AGB evolution.