Source author record

J. S. Vink

J. S. Vink 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

59works
5topics
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

59 published item(s)

preprint2022arXiv

Wind properties of Milky Way and SMC massive stars: empirical Z dependence from CMFGEN models

Detailed knowledge about stellar winds and evolution at different metallicities is crucial for understanding stellar populations and feedback in the Local Group of galaxies and beyond. Despite efforts in the literature, we still lack a comprehensive, empirical view of the dependence of wind properties on metallicity ($Z$). Here, we investigate the winds of O and B stars in the Milky Way (MW) and Small Magellanic Cloud (SMC). We gathered a sample of 96 stars analyzed by means of the NLTE code CMFGEN. We explored their wind strengths and terminal velocities to address the $Z$ dependence, over a large luminosity range. The empirical wind-luminosity relation (WLR) obtained updates and extends previous results in the literature. It reveals a luminosity and $Z$ dependence, in agreement with the radiatively driven wind theory. For bright objects ($\log L/L_\odot \gtrsim 5.4$), we infer that $\dot{M} \sim Z^{0.5-0.8}$. However, this dependence seems to get weaker or vanish at lower luminosities. The analysis of the terminal velocities suggests a shallow $Z^n$ dependence, with $n \sim 0.1-0.2$, but it should be confirmed with a larger sample and more accurate $V_{\infty}$ determinations. Recent results on SMC stars based on the PoWR code support our inferred WLR. On the other hand, recent bow-shocks measurements stand mostly above our derived WLR. Theoretical calculations of the WLR are not precise, specially at low $L$, where the results scatter. Deviations between our results and recent predictions are identified to be due to the weak wind problem and the extreme terminal velocities predicted by the models. The Z dependence suggested by our analysis deserves further investigations, given its astrophysical implications.

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

COBRaS: The e-MERLIN 21 cm Legacy survey of Cygnus OB2

The role of massive stars is central to an understanding of galactic ecology. It is important to establish the details of how massive stars provide radiative, chemical, and mechanical feedback in galaxies. Central to these issues is an understanding of the evolution of massive stars, and the critical role of mass loss via strongly structured winds and stellar binarity. Ultimately, massive stellar clusters shape the structure and energetics of galaxies. We aim to conduct high-resolution, deep field mapping at 21cm of the core of the massive Cygnus OB2 association and to characterise the properties of the massive stars and colliding winds at this waveband. We used seven stations of the e-MERLIN radio facility, with its upgraded bandwidth and enhanced sensitivity to conduct a 21cm census of Cygnus OB2. Based on 42 hours of observations, seven overlapping pointings were employed over multiple epochs during 2014 resulting in 1 sigma sensitivities down to ~21microJy and a resolution of ~180mas. A total of 61 sources are detected at 21cm over a ~0.48deg x 0.48deg region centred on the heart of the Cyg OB2 association. Of these 61 sources, 33 are detected for the first time. We detect a number of previously identified sources including four massive stellar binary systems, two YSOs, and several known X-ray and radio sources. We also detect the LBV candidate (possible binary system) and blue hypergiant (BHG) star of Cyg OB2 #12. The 21cm observations secured in the COBRaS Legacy project provide data to constrain conditions in the outer wind regions of massive stars; determine the non-thermal properties of massive interacting binaries; examine evidence for transient sources, including those associated with young stellar objects; and provide unidentified sources that merit follow-up observations. The 21cm data are of lasting value and will serve in combination with other key surveys of Cyg OB2.

preprint2020arXiv

Optical/near-infrared observations of the Fried Egg Nebula: Multiple shell ejections on a 100 yr timescale from a massive yellow hypergiant

Context. The fate of a massive star during the latest stages of its evolution is highly dependent on its mass-loss rate/geometry and therefore knowing the geometry of the circumstellar material close to the star and its surroundings is crucial. Aims. We aim to study the nature (i.e. geometry, rates) of mass-loss episodes. In this context, yellow hypergiants are great targets. Methods. We analyse a large set of optical/near-infrared data, in spectroscopic and photometric (X-shooter/VLT), spectropolarimetric (ISIS/WHT), and interferometric GRAVITY-AMBER/VLTI) modes, toward the yellow hypergiant IRAS 17163-3907. We present the first model-independent reconstructed images of IRAS 17163-3907 at these wavelengths at milli-arcsecond scales. Lastly, we apply a 2D radiative transfer model to fit the dereddened photometry and the radial profiles of published VISIR images at 8.59 μm, 11.85 μm and 12.81 μm simultaneously, adopting the revised Gaia distance (DR2). Results. The interferometric observables around 2 μm show that the Brγ emission is more extended and asymmetric than the Na i and the continuum emission. In addition to the two known shells surrounding IRAS 17163-3907 we report on the existence of a third hot inner shell with a maximum dynamical age of only 30 yr. Conclusions. The interpretation of the presence of Na i emission at closer distances to the star compared to Brγ has been a challenge in various studies. We argue that the presence of a pseudophotosphere is not needed, but it is rather an optical depth effect. The three observed distinct mass-loss episodes are characterised by different mass-loss rates and can inform the theories on mass-loss mechanisms, which is a topic still under debate. We discuss these in the context of photospheric pulsations and wind bi-stability mechanisms.

preprint2020arXiv

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

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

preprint2020arXiv

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

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

preprint2020arXiv

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

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

preprint2019arXiv

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

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

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

Linear spectropolarimetry across the optical spectrum of Herbig Ae/Be stars

We present the results of spectropolarimetric observations of 12 Herbig Ae/Be objects. Our data have the largest spectropolarimetric wavelength coverage, 4560 Å to 9480 Å, published to date. A change in linear polarisation across the Hα line, is detected in all objects. Such a line effect reveals the fact that stellar photons are scattered off free electrons that are not distributed in a spherically symmetric volume, suggesting the presence of small disks around these accreting objects. Thanks to the large wavelength coverage, we can report that Hα is the spectral line in the optical wavelength range that is most sensitive to revealing deviations from spherical symmetry, and the one most likely to show a line effect across the polarisation in such cases. Few other spectral lines display changes in polarisation across the line. In addition, Hα is the only line which shows an effect across its absorption component in some sources. We present a scenario explaining this finding and demonstrate that the detection of the line effect strongly relies on the number of photons scattered into our line of sight. We highlight the special case of R Mon, which is the only object in our sample to show many lines with a polarisation effect, which is much stronger than in all other objects. Given that the object and its nebulosity is spatially resolved, we argue that this is due to scattering of the stellar and emission spectrum off circumstellar dust.

preprint2016arXiv

Photoevaporation and close encounters: how the environment around Cygnus OB2 affects the evolution of protoplanetary disks

In our Galaxy, star formation occurs in a variety of environments, with a large fraction of stars formed in clusters hosting massive stars. OB stars have an important feedback on the evolution of protoplanetary disks around nearby young stars and likely on the process of planet formation occurring in them. The nearby massive association Cygnus OB2 is an outstanding laboratory to study this feedback. It is the closest massive association to our Sun, and hosts hundreds of massive stars and thousands of low mass members. In this paper, we analyze the spatial variation of the disk fraction in Cygnus OB2 and we study its correlation with the local values of Far and Extreme ultraviolet radiation fields and the local stellar surface density. We present definitive evidence that disks are more rapidly dissipated in the regions of the association characterized by intense local UV field and large stellar density. In particular, the FUV radiation dominates disks dissipation timescales in the proximity (i.e. within 0.5 pc) of the O stars. In the rest of the association, EUV photons potentially induce a significant mass loss from the irradiated disks across the entire association, but the efficiency of this process is reduced at increasing distances from the massive stars due to absorption by the intervening intracluster material. We find that disk dissipation due to close stellar encounters is negligible in Cygnus OB2, and likely to have affected 1% or fewer of the stellar population. Disk dissipation is instead dominated by photoevaporation. We also compare our results to what has been found in other young clusters with different massive populations, concluding that massive associations like Cygnus OB2 are potentially hostile to protoplanetary disks, but that the environments where disks can safely evolve in planetary systems are likely quite common in our Galaxy.

preprint2016arXiv

The deep OB star population in Carina from the VST Photometric H$α$ Survey (VPHAS+)

Massive OB stars are critical to the ecology of galaxies, and yet our knowledge of OB stars in the Milky Way, fainter than $V \sim 12$, remains patchy. Data from the VST Photometric H$α$ Survey (VPHAS+) permit the construction of the first deep catalogues of blue excess-selected OB stars, without neglecting the stellar field. A total of 14900 candidates with 2MASS cross-matches are blue-selected from a 42 square-degree region in the Galactic Plane, capturing the Carina Arm over the Galactic longitude range $282^{\circ} \lesssim \ell \lesssim 293^{\circ}$. Spectral energy distribution fitting is performed on these candidates' combined VPHAS+ $u,g,r,i$ and 2MASS $J,H,K$ magnitudes. This delivers: effective temperature constraints, statistically separating O from early-B stars; high-quality extinction parameters, $A_0$ and $R_V$ (random errors typically $< 0.1$). The high-confidence O-B2 candidates number 5915 and a further 5170 fit to later B spectral type. Spectroscopy of 276 of the former confirms 97% of them. The fraction of emission line stars among all candidate B stars is 7--8% . Greyer ($R_V > 3.5$) extinction laws are ubiquitous in the region, over the distance range 2.5--3 kpc to $\sim$10~kpc. Near prominent massive clusters, $R_V$ tends to rise, with particularly large and chaotic excursions to $R_V \sim 5$ seen in the Carina Nebula. The data reveal a hitherto unnoticed association of 108 O-B2 stars around the O5If$+$ star LSS 2063 ($\ell = 289.77^{\circ}$, $b = -1.22^{\circ}$). Treating the OB star scale-height as a constant within the thin disk, we find an orderly mean relation between extinction ($A_0$) and distance in the Galactic longitude range, $287.6^{\circ} < \ell < 293.5^{\circ}$, and infer the subtle onset of thin-disk warping. A halo around NGC 3603, roughly a degree in diameter, of $\sim$500 O-B2 stars with $4 < A_0 (\rm{mag}) < 7$ is noted.

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

Classical T Tauri stars with VPHAS$+$: I: H$α$ and $u$-band accretion rates in the Lagoon Nebula M8

We estimate the accretion rates of 235 Classical T Tauri star (CTTS) candidates in the Lagoon Nebula using $ugri$H$α$ photometry from the VPHAS+ survey. Our sample consists of stars displaying H$α$-excess, the intensity of which is used to derive accretion rates. For a subset of 87 stars, the intensity of the $u$-band excess is also used to estimate accretion rates. We find the mean variation in accretion rates measured using H$α$ and $u$-band intensities to be $\sim$ 0.17 dex, agreeing with previous estimates (0.04-0.4 dex) but for a much larger sample. The spatial distribution of CTTS align with the location of protostars and molecular gas suggesting that they retain an imprint of the natal gas fragmentation process. Strong accretors are concentrated spatially, while weak accretors are more distributed. Our results do not support the sequential star forming processes suggested in the literature.

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

New OB star candidates in the Carina Arm around Westerlund 2 from VPHAS+

O and early B stars are at the apex of galactic ecology, but in the Milky Way, only a minority of them may yet have been identified. We present the results of a pilot study to select and parametrise OB star candidates in the Southern Galactic plane, down to a limiting magnitude of $g=20$. A 2 square-degree field capturing the Carina Arm around the young massive star cluster, Westerlund 2, is examined. The confirmed OB stars in this cluster are used to validate our identification method, based on selection from the $(u-g, g-r)$ diagram for the region. Our Markov Chain Monte Carlo fitting method combines VPHAS+ $u, g, r, i$ with published $J, H, K$ photometry in order to derive posterior probability distributions of the stellar parameters $\log(\rm T_{\rm eff})$ and distance modulus, together with the reddening parameters $A_0$ and $R_V$. The stellar parameters are sufficient to confirm OB status while the reddening parameters are determined to a precision of $σ(A_0)\sim0.09$ and $σ(R_V)\sim0.08$. There are 489 objects that fit well as new OB candidates, earlier than $\sim$B2. This total includes 74 probable massive O stars, 5 likely blue supergiants and 32 reddened subdwarfs. This increases the number of previously known and candidate OB stars in the region by nearly a factor of 10. Most of the new objects are likely to be at distances between 3 and 6 kpc. We have confirmed the results of previous studies that, at these longer distances, these sight lines require non-standard reddening laws with $3.5<R_V<4$.

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 XVIII. Classifications and radial velocities of the B-type stars

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

preprint2015arXiv

The VLT-FLAMES Tarantula Survey 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

A 3D extinction map of the Northern Galactic Plane based on IPHAS photometry

We present a three dimensional map of extinction in the Northern Galactic Plane derived using photometry from the IPHAS survey. The map has fine angular ($\sim 10$ arcmin) and distance (100 pc) sampling allied to a significant depth ($\gtrsim 5$ kpc). We construct the map using a method based on a hierarchical Bayesian model as previously described by Sale (2012). In addition to mean extinction, we also measure differential extinction, which arises from the fractal nature of the ISM, and show that it will be the dominant source of uncertainty in estimates of extinction to some arbitrary position. The method applied also furnishes us with photometric estimates of the distance, extinction, effective temperature, surface gravity, and mass for $\sim 38$ million stars. Both the extinction map and the catalogue of stellar parameters are made publicly available via http://www.iphas.org/extinction .

preprint2014arXiv

Pre-main sequence accretion in the low metallicity Galactic star-forming region Sh 2-284

We present optical spectra of pre-main sequence (PMS) candidates around the H$α$ region taken with the Southern African Large Telescope, SALT, in the low metallicity ($Z$) Galactic region Sh 2-284, which includes the open cluster Dolidze 25 with an atypical low metallicity of $Z$ $\sim$ 1/5 $Z_{\odot}$. It has been suggested on the basis of both theory and observations that PMS mass-accretion rates, $\dot M_{\rm{acc}}$, are a function of $Z$. We present the first sample of spectroscopic estimates of mass-accretion rates for PMS stars in any low-$Z$ star-forming region. Our data-set was enlarged with literature data of H$α$ emission in intermediate-resolution R-band spectroscopy. Our total sample includes 24 objects spanning a mass range between 1 - 2 $M_{\odot}$ and with a median age of approximately 3.5 Myr. The vast majority (21 out of 24) show evidence for a circumstellar disk on the basis of 2MASS and Spitzer infrared photometry. We find $\dot M_{\rm{acc}}$ in the 1 - 2 $M_{\odot}$ interval to depend quasi-quadratically on stellar mass, with $\dot M_{\rm{acc}}$ $\propto$ $M_{\ast}^{2.4\,\pm\,0.35}$, and inversely with stellar age $\dot M_{\rm{acc}}$ $\propto$ $t_{\ast}^{-\,0.7\,\pm\,0.4}$. Furthermore, we compare our spectroscopic $\dot M_{\rm{acc}}$ measurements with solar $Z$ Galactic PMS stars in the same mass range, but, surprisingly find no evidence for a systematic change in $\dot M_{\rm{acc}}$ with $Z$. We show that literature accretion-rate studies are influenced by detection limits, and we suggest that $\dot M_{\rm{acc}}$ may be controlled by factors other than $Z_{\ast}$, $M_{\ast}$, and age.

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 Second Data Release of the INT Photometric H-Alpha Survey of the Northern Galactic Plane (IPHAS DR2)

The INT/WFC Photometric H-Alpha Survey of the Northern Galactic Plane (IPHAS) is a 1800 square degrees imaging survey covering Galactic latitudes |b| < 5 deg and longitudes l = 30 to 215 deg in the r, i and H-alpha filters using the Wide Field Camera (WFC) on the 2.5-metre Isaac Newton Telescope (INT) in La Palma. We present the first quality-controlled and globally-calibrated source catalogue derived from the survey, providing single-epoch photometry for 219 million unique sources across 92% of the footprint. The observations were carried out between 2003 and 2012 at a median seeing of 1.1 arcsec (sampled at 0.33 arcsec/pixel) and to a mean 5σ-depth of 21.2 (r), 20.0 (i) and 20.3 (H-alpha) in the Vega magnitude system. We explain the data reduction and quality control procedures, describe and test the global re-calibration, and detail the construction of the new catalogue. We show that the new calibration is accurate to 0.03 mag (rms) and recommend a series of quality criteria to select the most reliable data from the catalogue. Finally, we demonstrate the ability of the catalogue's unique (r-Halpha, r-i) diagram to (1) characterise stellar populations and extinction regimes towards different Galactic sightlines and (2) select H-alpha emission-line objects. IPHAS is the first survey to offer comprehensive CCD photometry of point sources across the Galactic Plane at visible wavelengths, providing the much-needed counterpart to recent infrared surveys.

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 VST Photometric Halpha Survey of the Southern Galactic Plane and Bulge (VPHAS+)

The VST Photometric Halpha Survey of the Southern Galactic Plane and Bulge (VPHAS+) is surveying the southern Milky Way in u, g, r, i and Halpha at 1 arcsec angular resolution. Its footprint spans the Galactic latitude range -5 < b < +5 at all longitudes south of the celestial equator. Extensions around the Galactic Centre to Galactic latitudes +/-10 bring in much of the Galactic Bulge. This ESO public survey, begun on 28th December 2011, reaches down to 20th magnitude (10-sigma) and will provide single-epoch digital optical photometry for around 300 million stars. The observing strategy and data pipelining is described, and an appraisal of the segmented narrowband Halpha filter in use is presented. Using model atmospheres and library spectra, we compute main-sequence (u - g), (g - r), (r - i) and (r - Halpha) stellar colours in the Vega system. We report on a preliminary validation of the photometry using test data obtained from two pointings overlapping the Sloan Digital Sky Survey. An example of the (u - g, g - r) and (r - Halpha, r - i) diagrams for a full VPHAS+ survey field is given. Attention is drawn to the opportunities for studies of compact nebulae and nebular morphologies that arise from the image quality being achieved. The value of the u band as the means to identify planetary-nebula central stars is demonstrated by the discovery of the central star of NGC 2899 in survey data. Thanks to its excellent imaging performance, the VST/OmegaCam combination used by this survey is a perfect vehicle for automated searches for reddened early-type stars, and will allow the discovery and analysis of compact binaries, white dwarfs and transient sources.

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 - XI. A census of the hot luminous stars and their feedback in 30 Doradus

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

preprint2013arXiv

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

Aims. Using ground based multi-object optical spectroscopy obtained in the framework of the VLT-FLAMES Tarantula Survey (VFTS), we aim to establish the (projected) rotational velocity distribution for a sample of 216 presumably single O-type stars in 30 Doradus (30 Dor). Methods. We measured projected rotational velocities, \vrot, by means of a Fourier transform method and a profile fitting method applied on a set of isolated spectral lines. We also used an iterative deconvolution procedure to infer the probability density, $\rm{P(\veq)}$, of the equatorial rotational velocity, \veq. Results. The distribution of \vrot\ shows a two-component structure: a peak around 80 \kms\ and a high-velocity tail extending up to $\sim$600 \kms. This structure is also present in the inferred distribution $\rm{P(\veq)}$ with around 80% of the sample having 0 $<$ \veq\, $\leq\, 300$ \kms\ and the other 20% distributed in the high-velocity region. Conclusions. Most of the stars in our sample rotate with a rate less than 20%\ of their break-up velocity. For the bulk of the sample, mass-loss in a stellar wind and/or envelope expansion is not efficient enough to significantly spin down these stars within the first few Myr of evolution. The presence of a sizeable population of fast rotators is compatible with recent population synthesis computations that investigate the influence of binary evolution on the rotation rate of massive stars. Despite the fact that we have excluded stars that show significant radial velocity variations, our sample may have remained contaminated by post-interaction binary products. The fact that the high-velocity tail may be preferentially (and perhaps even exclusively), populated by post-binary interaction products, has important implications for the evolutionary origin of systems that produce gamma-ray bursts.

preprint2013arXiv

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

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

preprint2012arXiv

Energetic feedback and $^{26}$Al from massive stars and their supernovae in the Carina region

We study the populations of massive stars in the Carina region and their energetic feedback and ejection of $^{26}$Al. We did a census of the stellar populations in young stellar clusters within a few degrees of the Carina Nebula. For each star we estimated the mass, based on the spectral type and the host cluster age. We used population synthesis to calculate the energetic feedback and ejection of $^{26}$Al from the winds of the massive stars and their supernova explosions. We used 7 years of INTEGRAL observations to measure the $^{26}$Al signal from the region. The INTEGRAL $^{26}$Al signal is not significant with a best-fit value of about 1.5e-5 ph/cm^2/s, approximately half of the published Compton Gamma Ray Observatory (CGRO) result, but in agreement with the latest CGRO estimates. Our analysis of the stellar populations in the young clusters leads to an expected signal of half the observed value, but the results are consistent within 2 sigma. We find that the fraction of $^{26}$Al ejected in Wolf-Rayet winds is high, and the observed signal is unlikely to be caused by $^{26}$Al ejected in supernovae alone, indicating a strong wind ejection of $^{26}$Al. Due to the lack of prominent O stars, regions with ages $\gtrsim$10 Myr are often neglected in studies of OB associations. We find that in the Carina region such clusters contribute significantly to the stellar mass and the energetics of the region.

preprint2012arXiv

LAMP: The Long-term Accretion Monitoring Program of T Tauri stars in Chamaeleon-I

We present the results of a variability study of accreting young stellar objects in the Chameleon I star-forming region which is based on ~300 high resolution optical spectra from the multi-object fibre spectrograph FLAMES/GIRAFFE at the ESO/VLT. Twenty five objects with spectral types from G2-M5.75 were observed 12 times over the course of 15 months. Using the emission lines Ha (6562.81 A) and Ca II (8662.1 A) as accretion indicators we found 10 accreting and 15 non-accreting objects. We derived accretion rates for all accretors in the sample using the Ha equivalent width, Ha 10% width and the CaII equivalent width. The mean amplitude of variations in derived accretion rate from Ha equivalent width was ~ 0.37 dex, from Ca II equivalent width ~0.83 dex and from Ha 10% width ~1.11 dex. Based on the large amplitude of variations in accretion rates derived from the Ha 10% width with respect to the other diagnostics, we do not consider it to be a reliable accretion rate estimator. Taking the variations in Ha equivalent width and CaII equivalent width accretion rates to be closer to the true value, they suggest that the spread which has been found around the accretion rate to stellar mass relation is not due to the variability of individual objects on time-scales of weeks to ~1 year. From these variations we can also infer that the accretion rates are stable within < 0.37 dex over time-scales of less than 15 months. A major portion of the accretion variability was found to occur on less than the shortest time-scales in our observations, 8-25 days, which is comparable with the rotation periods of these young stellar objects. This could be an indication that what we are probing is spatial structure in the accretion flows, and also suggests that observations on time-scales of ~a couple of weeks are sufficient to limit the total extent of accretion rate variations in typical young stars.

preprint2012arXiv

On the nature and detectability of Type Ib/c supernova progenitors

The progenitors of many Type II supernovae have been observationally identified but the search for Type Ibc supernova (SN Ibc) progenitors has thus far been unsuccessful, despite the expectation that they are luminous Wolf-Rayet (WR) stars. We investigate how the evolution of massive helium stars affects their visual appearances, and discuss the implications for the detectability of SN Ibc progenitors. Massive WR stars that rapidly lose their helium envelopes through stellar-wind mass-loss end their lives when their effective temperatures -- related to their hydrostatic surfaces -- exceed about 150kK.At their pre-supernova stage, their surface properties resemble those of hot Galactic WR stars of WO sub-type. These are visually faint with narrow-band visual magnitudes Mv = -1.5 ~ -2.5, despite their high bolometric luminosities (log L/Lsun = 5.6 ~ 5.7), compared to the bulk of Galactic WR stars (Mv < -4). In contrast, relatively low-mass helium stars that retain a thick helium envelope appear fairly bright in optical bands, depending on the final masses and the history of the envelope expansion during the late evolutionary stages. We conclude that SNe Ibc observations have so far not provided strong constraints on progenitor bolometric luminosities and masses, even with the deepest searches. We also argue that Ic progenitors are more challenging to identify than Ib progenitors in any optical images.

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.

preprint2012arXiv

The VLT-FLAMES Tarantula Survey V. The peculiar B[e]-like supergiant, VFTS698, in 30 Doradus

We present an analysis of a peculiar supergiant B-type star (VFTS698/Melnick 2/Parker 1797) in the 30 Doradus region of the Large Magellanic Cloud which exhibits characteristics similar to the broad class of B[e] stars. We analyse optical spectra from the VLT-FLAMES survey, together with archival optical and infrared photometry and X-ray imaging to characterise the system. We find radial velocity variations of around 400 km/s in the high excitation Si IV, N III and He II spectra, and photometric variability of ~0.6 mag with a period of 12.7 days. In addition, we detect long-term photometric variations of ~0.25 mag, which may be due to a longer-term variability with a period of ~400 days. We conclude that VFTS698 is likely an interacting binary comprising an early B-type star secondary orbiting a veiled, more massive companion. Spectral evidence suggests a mid-to-late B-type primary, but this may originate from an optically-thick accretion disc directly surrounding the primary.

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

VLTI/AMBER observations of the binary B[e] supergiant HD 327083

HD 327083 is a luminous B type star which exhibits emission lines and an infrared excess and is therefore classified as a supergiant B[e] star. In addition, the star is the primary of a close binary system. It is not clear whether the B[e] behaviour of HD 327083 is related to its binarity or its evolutionary state. Here we address this issue by studying its circumstellar environment with high spatial resolution. To this end, we have observed HD 327083 with the VLTI and AMBER in the medium resolution K-band setting. 13CO bandhead emission is detected, confirming HD 327083 is a post-main sequence object. The observations spatially resolve the source of the NIR continuum and the Br-gamma and CO line emission. In addition, differential phase measurements allow us to probe the origin of the observed Br-gamma emission with sub-mas precision. Using geometrical models, we find that the visibilities and closure phases suggest that the close binary system is surrounded by a circum-binary disk. We also find that in the case of the binary HD 327083, the relative sizes of the continuum and Br-gamma emitting regions are different to those of a single supergiant B[e] star where the standard dual outflow scenario is thought to apply. These findings are consistent with the hypothesis that the mass loss of HD 327083 is related to its binary nature.

preprint2011arXiv

On the alignment between the circumstellar disks and orbital planes of Herbig Ae/Be binary systems

The majority of the pre-main-sequence Herbig Ae/Be stars reside in binary systems. As these systems are young, their properties may contain an imprint of the star formation process at intermediate masses. However, these systems are generally spatially unresolved, making it difficult to probe their circumstellar environment to search for manifestations of their formation process, such as accretion disks. Here we investigate the formation mechanism of Herbig Ae/Be (HAe/Be) binary systems by studying the relative orientation of their binary orbits and circumstellar disks. We present linear spectropolarimetric observations of HAe/Be stars over the Halpha line, which are used to determine the orientation of their circumstellar disks. In conjunction with data from the literature, we obtain a sample of 20 binaries with known disk position angles (PAs). We compare our disk PA data to a model to investigate whether HAe/Be binary systems and their disks are co-planar. Moreover, in the light of a recent suggestion that some HAe/Be star spectropolarimetric signatures may not necessarily be related to circumstellar disks, we re-assess the relationship between spectropolarimetric signatures and disk PAs. We do this by comparing spectropolarimetric and high spatial resolution observations of young stellar objects. We find that spectropolarimetric observations of pre-main-sequence stars do trace circumstellar disks. This finding is significant above the 3 sigma level. Using this result, our data are entirely consistent with the situation in which HAe/Be binary systems and circumstellar disks are co-planar, while random orientations can be rejected at the 2.2 sigma level. We conclude that HAe/Be binary systems and the disks around the primaries are likely to be co-planar, which is consistent with the notion that these systems form via monolithic collapse and subsequent disk fragmentation (abridged).

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

preprint2010arXiv

Are the Stellar Winds in IC1613 stronger than expected?

In this poster we present the results of our analyses of three early massive stars in IC 1613, whose spectra have been observed with VIMOS and analyzed with CMFGEN and FASTWIND. One of the targets resulted a possible LBV and the other two are Of stars with unexpectedly strong winds. The Of stars seem to be strongly contaminated by CNO products. Our preliminary results may represent a challenge for the theory of stellar atmospheres, but they still have to be confirmed by the analysis of more objects and a more complete coverage of the parameter space.

preprint2010arXiv

Multiple Major Outbursts from a Restless Luminous Blue Variable in NGC 3432

We present new photometric and spectroscopic observations of an unusual luminous blue variable (LBV) in NGC 3432, covering three major outbursts in October 2008, April 2009 and November 2009. Previously, this star experienced an outburst also in 2000 (known as SN 2000ch). During outbursts the star reached an absolute magnitude between -12.1 and -12.8. Its spectrum showed H, He I and Fe II lines with P-Cygni profiles during and soon after the eruptive phases, while only intermediate-width lines in pure emission (including He II 4686A were visible during quiescence. The fast-evolving light curve soon after the outbursts, the quasi-modulated light curve, the peak magnitude and the overall spectral properties are consistent with multiple episodes of variability of an extremely active LBV. However, the widths of the spectral lines indicate unusually high wind velocities (1500-2800 km/s), similar to those observed in Wolf-Rayet stars. Although modulated light curves are typical of LBVs during the S-Dor variability phase, the luminous maxima and the high frequency of outbursts are unexpected in S-Dor variables. Such extreme variability may be associated with repeated ejection episodes during a giant eruption of an LBV. Alternatively, it may be indicative of a high level of instability shortly preceding the core-collapse or due to interaction with a massive, binary companion. In this context, the variable in NGC 3432 shares some similarities with the famous stellar system HD 5980 in the Small Magellanic Cloud, which includes an erupting LBV and an early Wolf-Rayet star.

preprint2010arXiv

Probing the evolving massive star population in Orion with kinematic and radioactive tracers

We assemble a census of the most massive stars in Orion, then use stellar isochrones to estimate their masses and ages, and use these results to establish the stellar content of Orion's individual OB associations. From this, our new population synthesis code is utilized to derive the history of the emission of UV radiation and kinetic energy of the material ejected by the massive stars, and also follow the ejection of the long-lived radioactive isotopes 26Al and 60Fe. In order to estimate the precision of our method, we compare and contrast three distinct representations of the massive stars. We compare the expected outputs with observations of 26Al gamma-ray signal and the extent of the Eridanus cavity. We find an integrated kinetic energy emitted by the massive stars of 1.8(+1.5-0.4)times 10^52 erg. This number is consistent with the energy thought to be required to create the Eridanus superbubble. We also find good agreement between our model and the observed 26Al signal, estimating a mass of 5.8(+2.7-2.5) times 10^-4 Msol of 26Al in the Orion region. Our population synthesis approach is demonstrated for the Orion region to reproduce three different kinds of observable outputs from massive stars in a consistent manner: Kinetic energy as manifested in ISM excavation, ionization as manifested in free-free emission, and nucleosynthesis ejecta as manifested in radioactivity gamma-rays. The good match between our model and the observables does not argue for considerable modifications of mass loss. If clumping effects turn out to be strong, other processes would need to be identified to compensate for their impact on massive-star outputs. Our population synthesis analysis jointly treats kinematic output and the return of radioactive isotopes, which proves a powerful extension of the methodology that constrains feedback from massive stars.

preprint2010arXiv

The MiMeS Project: Overview and Current Status

The Magnetism in Massive Stars (MiMeS) Project is a consensus collaboration among many of the foremost international researchers of the physics of hot, massive stars, with the basic aim of understanding the origin, evolution and impact of magnetic fields in these objects. At the time of writing, MiMeS Large Programs have acquired over 950 high-resolution polarised spectra of about 150 individual stars with spectral types from B5-O4, discovering new magnetic fields in a dozen hot, massive stars. The quality of this spectral and magnetic matériel is very high, and the Collaboration is keen to connect with colleagues capable of exploiting the data in new or unforeseen ways. In this paper we review the structure of the MiMeS observing programs and report the status of observations, data modeling and development of related theory.

preprint2010arXiv

The nature of V39: an LBV candidate or LBV impostor in the very low metallicity galaxy IC 1613?

[abridged] Context: Very few examples of luminous blue variable (LBV) stars or LBV candidates (LBVc) are known, particularly at metallicities below the SMC. The LBV phase is crucial for the evolution of massive stars, and its behavior with metallicity is poorly known. V39 in IC 1613 is a well-known photometric variable, with B-band changes larger than 1mag. over its period. The star, previously proposed to be a projection of a Galactic W Virginis and an IC 1613 red supergiant, shows features that render it a possible LBVc. Method: We investigate mid-resolution blue and red VLT-VIMOS spectra of V39, covering a time span of 40 days, and perform a quantitative analysis of the combined spectrum using the model atmosphere code CMFGEN. Results: We identify strong Balmer and FeII P-Cygni profiles, and a hybrid spectrum resembling a B-A supergiant in the blue and a G-star in the red. No significant Vrad variations are detected, and the spectral changes are small over the photometric period. Our analysis places V39 in the low-luminosity part of the LBV and LBVc region, but it is also consistent with a sgB[e] star. Conclusions: The radial velocity indicates that V39 belongs to IC 1613. The lack of Vrad changes and spectroscopic variations excludes binary scenarios. The features observed are not consistent with a W Virginis star, and this possibility is also discarded. We propose that the star is a B-A LBVc or sgB[e] star surrounded by a thick disk precessing around it. If confirmed, V39 would be the lowest metallicity resolved LBV candidate known to date. Alternatively, it could represent a new transient phase of massive star evolution, an LBV impostor.

preprint2009arXiv

IPHAS and the symbiotic stars. II. New discoveries and a sample of the most common mimics

In a previous paper [arXiv:0712.2391], we presented the selection criteria needed to search for symbiotic stars in IPHAS, the INT Halpha survey of the Northern Galactic plane. IPHAS gives us the opportunity to make a systematic, complete search for symbiotic stars in a magnitude-limited volume. Follow-up spectroscopy at different telescopes worldwide of a sample of sixty two symbiotic star candidates is presented. Seven out of nineteen S-type candidates observed spectroscopically are confirmed to be genuine symbiotic stars. The spectral type of their red giant components, as well as reddening and distance, were computed by modelling the spectra. Only one new D-type symbiotic system, out of forty-three candidates observed, was found. This was as expected (see discussion in our paper on the selection criteria). The object shows evidence for a high density outflow expanding at a speed larger than 65 km/s. Most of the other candidates are lightly reddened classical T Tauri stars and more highly reddened young stellar objects that may be either more massive young stars of HAeBe type or classical Be stars. In addition, a few notable objects have been found, such as three new Wolf-Rayet stars and two relatively high-luminosity evolved massive stars. We also found a helium-rich source, possibly a dense ejecta hiding a WR star, which is surrounded by a large ionized nebula.

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

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.

preprint2009arXiv

Using population synthesis of massive stars to study the interstellar medium near OB associations

We developed a new population synthesis code for groups of massive stars, where we model the emission of different forms of energy and matter from the stars of the association. In particular, the ejection of the two radioactive isotopes 26Al and 60Fe is followed, as well as the emission of hydrogen ionizing photons, and the kinetic energy of the stellar winds and supernova explosions. We investigate various alternative astrophysical inputs and the resulting output sensitivities, especially effects due to the inclusion of rotation in stellar models. As the aim of the code is the application to relatively small populations of massive stars, special care is taken to address their statistical properties. Our code incorporates both analytical statistical methods applicable to small populations, as well as extensive Monte Carlo simulations. We find that the inclusion of rotation in the stellar models has a large impact on the interactions between OB associations and their surrounding interstellar medium. The emission of 26Al in the stellar winds is strongly enhanced, compared to non-rotating models with the same mass-loss prescription. This compensates the recent reductions in the estimates of mass-loss rates of massive stars due to the effects of clumping. Despite the lower mass-loss rates, the power of the winds is actually enhanced for rotating stellar models. The supernova power (kinetic energy of their ejecta) is decreased due to longer lifetimes of rotating stars, and therefore the wind power dominates over supernova power for the first 6 Myr after a burst of star-formation. For populations typical of nearby star-forming regions, the statistical uncertainties are large and clearly non-Gaussian.