BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Jorick S. Vink

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

ResearcherUnclaimed source record
astro-ph.SRastro-ph.GAastro-ph.HEastro-ph.COastro-ph.IMastro-phgr-qc

Catalog footprint

What is connected

53works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

53 published item(s)

preprint2025arXiv

Metal-poor single Wolf-Rayet stars: The interplay of optically thick winds and rotation

The Small Magellanic Cloud (SMC) hosts 12 known Wolf-Rayet (WR) stars, seven of which are apparently single. Their formation is a challenge for current stellar evolution models because line-driven winds are generally assumed to be quenched at a metallicity of Z < 0.004. Here, we present a set of mesa models of single stars with zero-age main sequence masses of 20 - 80 Msun considering different initial rotation speeds (Ω = 0 - 0.7 Ω_c), metallicities (Z = 0.002 - 0.0045), and wind mass-loss models (optically thin and thick winds). We show that if we account for optically thick winds, fast rotating (Ω = 0.6 Ω_c) single metal-poor O-type stars (with M > 20 Msun) shed their envelope and become WR stars even at the low metallicity of the SMC. The luminosity, effective temperature, evolutionary timescale, surface abundance, and rotational velocity of our simulated WR stars are compatible to the WRs observed in the SMC. We speculate that this scenario can also alleviate the excess of giant stars across the Humphreys-Davidson limit. Our results have key implications for black hole masses, (pair instability) supernova explosions, and other observable signatures.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Potential Subpopulations and Assembling Tendency of the Merging Black Holes

The origins of coalescing binary black holes (BBHs) detected by the advanced LIGO/Virgo are still under debate, and clues may be present in the joint mass-spin distribution of these merger events. Here we construct phenomenological models containing two sub-populations to investigate the BBH population detected in gravitational wave observations. We find that our models can explain the GWTC-3 data rather well, and several constraints to our model are required by the data: first, the maximum mass for the component with a stellar-origin, $m_{\rm max}$, is $39.1^{+2.4}_{-2.7}M_{\odot}$ at 90\% credibility; second, about $15\%$ of the mergers happen in dynamical environments, in which $7-16\%$ of events are hierarchical mergers, and these BHs have an average spin magnitude significantly larger than the first-generation mergers, with ${\rm d}μ_{\rm a} > 0.4 $ at $99\%$ credibility; third, the dynamical component BHs tend to pair with each other with larger total mass and higher mass ratio. An independent analysis focusing on spins is also carried out, and we find that the spin amplitude of component BHs can be divided into two groups according to a division mass $m_{\rm d} = 46.1^{+5.6}_{-5.1}M_{\odot}$. These constraints can be naturally explained by current formation channels, and our results suggest that some of the observed events were likely from AGN disks.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

An X-ray quiet black hole born with a negligible kick in a massive binary within the Large Magellanic Cloud

Stellar-mass black holes are the final remnants of stars born with more than 15 solar masses. Billions are expected to reside in the Local Group, yet only few are known, mostly detected through X-rays emitted as they accrete material from a companion star. Here, we report on VFTS 243: a massive X-ray faint binary in the Large Magellanic Cloud. With an orbital period of 10.4-d, it comprises an O-type star of 25 solar masses and an unseen companion of at least nine solar masses. Our spectral analysis excludes a non-degenerate companion at a 5-sigma confidence level. The minimum companion mass implies that it is a black hole. No other X-ray quiet black hole is unambiguously known outside our Galaxy. The (near-)circular orbit and kinematics of VFTS 243 imply that the collapse of the progenitor into a black hole was associated with little or no ejected material or black-hole kick. Identifying such unique binaries substantially impacts the predicted rates of gravitational-wave detections and properties of core-collapse supernovae across the Cosmos.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Mass-loss implementation and temperature evolution of very massive stars

Very massive stars (VMS) dominate the physics of young clusters due to their ionising radiation and extreme stellar winds. It is these winds that determine their lifepaths until expiration. Observations in the Arches cluster show that VMS all have similar temperatures. The VLT-Flames Tarantula survey analysed VMS in the 30 Dor region of the LMC also finding a narrow range of temperatures, albeit at higher values - likely a metallicity effect. Using MESA, we study the main-sequence evolution of VMS with a new mass-loss recipe that switches from optically-thin O-star winds to optically-thick Wolf-Rayet type winds through the model-independent transition mass-loss rate of Vink & Gräfener. We examine the temperature evolution of VMS with mass loss that scales with the luminosity-over-mass (L/M) ratio and the Eddington parameter ($Γ_{\rm e}$), assessing the relevance of the surface hydrogen (H) abundance which sets the number of free electrons. We present grids of VMS models at Galactic and LMC metallicity and compare our temperature predictions with empirical results. Models with a steep $Γ_{\rm e}$-dependence evolve horizontally in the Hertzsprung-Russel (HR) diagram at nearly constant luminosities, requiring a delicate and unlikely balance between envelope inflation and enhanced mass loss over the entire VMS mass range. By contrast, models with a steep L/M-dependent mass loss are shown to evolve vertically in the HR-diagram at nearly constant Teff, naturally reproducing the narrow range of observed temperatures, as well as the correct trend with metallicity. This distinct behavior of a steeply dropping luminosity is a self-regulatory mechanism that keeps temperatures constant during evolution in the HR-diagram.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Resolving the core of R136 in the optical

The sharpest optical images of the R136 cluster in the Large Magellanic Cloud are presented, allowing for the first time to resolve members of the central core, including R136a1, the most massive star known. These data were taken using the Gemini speckle imager Zorro in medium-band filters with effective wavelengths similar to BVRI achieving angular resolutions between 30-40 mas. All stars previously known in the literature, having $V<$16 mag within the central $2''\times2''$ were recovered. Visual companions ($\geq$40 mas; 2000 au) were detected for the WN5h stars R136 a1, and a3. Photometry of the visual companion of a1 suggests it is of mid O spectral type. Based on new photometric luminosities using the resolved Zorro imaging, the masses of the individual WN5h stars are estimated to be between 150-200 $M_{\odot}$, lowering significantly the present-day masses of some of the most massive stars known. These mass estimates are critical anchor points for establishing the stellar upper-mass function.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Rubin Observatory LSST Transients and Variable Stars Roadmap

The Vera C. Rubin Legacy Survey of Space and Time holds the potential to revolutionize time domain astrophysics, reaching completely unexplored areas of the Universe and mapping variability time scales from minutes to a decade. To prepare to maximize the potential of the Rubin LSST data for the exploration of the transient and variable Universe, one of the four pillars of Rubin LSST science, the Transient and Variable Stars Science Collaboration, one of the eight Rubin LSST Science Collaborations, has identified research areas of interest and requirements, and paths to enable them. While our roadmap is ever-evolving, this document represents a snapshot of our plans and preparatory work in the final years and months leading up to the survey's first light.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The R136 star cluster dissected with Hubble Space Telescope/STIS. III. The most massive stars and their clumped winds

Context: The star cluster R136 inside the LMC hosts a rich population of massive stars, including the most massive stars known. The strong stellar winds of these very luminous stars impact their evolution and the surrounding environment. We currently lack detailed knowledge of the wind structure that is needed to quantify this impact. Aims: To observationally constrain the stellar and wind properties of the massive stars in R136, in particular the parameters related to wind clumping. Methods: We simultaneously analyse optical and UV spectroscopy of 53 O-type and 3 WNh-stars using the FASTWIND model atmosphere code and a genetic algorithm. The models account for optically thick clumps and effects related to porosity and velocity-porosity, as well as a non-void interclump medium. Results: We obtain stellar parameters, surface abundances, mass-loss rates, terminal velocities and clumping characteristics and compare these to theoretical predictions and evolutionary models. The clumping properties include the density of the interclump medium and the velocity-porosity of the wind. For the first time, these characteristics are systematically measured for a wide range of effective temperatures and luminosities. Conclusions: We confirm a cluster age of 1.0-2.5 Myr and derive an initial stellar mass of $\geq 250 {\rm M}_\odot$ for the most massive star in our sample, R136a1. The winds of our sample stars are highly clumped, with an average clumping factor of $f_{\rm cl}=29\pm15$. We find tentative trends in the wind-structure parameters as a function of mass-loss rate, suggesting that the winds of stars with higher mass-loss rates are less clumped. We compare several theoretical predictions to the observed mass-loss rates and terminal velocities and find that none satisfactorily reproduces both quantities. The prescription of Krtička & Kubát (2018) matches best the observed mass-loss rates.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Theory and Diagnostics of Hot Star Mass Loss

Massive stars have strong stellar winds that direct their evolution through the upper Hertzsprung-Russell diagram and determine the black hole mass function. Secondly, wind strength dictates the atmospheric structure that sets the ionising flux. Thirdly, the wind directly intervenes with the stellar envelope structure, which is decisive for both single star and binary evolution, affecting predictions for gravitational wave events. Key findings of current hot-star research include: * The traditional line-driven wind theory is being updated with Monte Carlo and co-moving frame computations, revealing a rich multi-variate behaviour of the mass-loss rate dM/dt in terms of M, L, Eddington Gamma, Teff , and chemical composition Z. Concerning the latter, dM/dt is shown to depend on the iron (Fe) opacity, making Wolf-Rayet populations, and gravitational wave events dependent on host galaxy Z. * On top of smooth mass-loss behaviour, there are several transitions in the Hertzsprung-Russell diagram, involving bi-stability jumps around Fe recombination temperatures, leading to quasi-stationary episodic, and not necessarily eruptive, Luminous Blue Variable and pre-SN mass loss. * Moreover, there are kinks. At 100 Solar Masses, a high Gamma mass-loss transition implies that hydrogen-rich very massive stars have higher mass-loss rates than commonly considered. At the other end of the mass spectrum, low-mass stripped helium stars no longer appear as Wolf-Rayet stars, but as optically-thin stars. These stripped stars, in addition to very massive stars, are two newly identified sources of ionising radiation that could play a key role in local star formation as well as at high-redshift.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Ultraviolet Spectropolarimetry with Polstar: Clumping and Mass-loss Rate Corrections

The most massive stars are thought to lose a significant fraction of their mass in a steady wind during the main-sequence and blue supergiant phases. This in turn sets the stage for their further evolution and eventual supernova, with consequences for ISM energization and chemical enrichment. Understanding these processes requires accurate observational constraints on the mass-loss rates of the most luminous stars, which can also be used to test theories of stellar wind generation. In the past, mass-loss rates have been characterized via collisional emission processes such as H$α$ and free-free radio emission, but these so-called "density squared" diagnostics require correction in the presence of widespread clumping. Recent observational and theoretical evidence points to the likelihood of a ubiquitously high level of such clumping in hot-star winds, but quantifying its effects requires a deeper understanding of the complex dynamics of radiatively driven winds. Furthermore, large-scale structures arising from surface anisotropies and propagating throughout the wind can further complicate the picture by introducing further density enhancements, affecting mass-loss diagnostics. Time series spectroscopy of UV resonance lines with high resolution and high signal-to-noise are required to better understand this complex dynamics, and help correct "density squared" diagnostics of mass-loss rates. The proposed Polstar mission easily provides the necessary resolution at the sound-speed scale of 20 km s$^{-1}$, on three dozen bright targets with signal-to noise an order of magnitude above that of the celebrated IUE MEGA campaign, via continuous observations that track structures advecting through the wind in real time.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A theoretical investigation of the Humphreys-Davidson limit at high and low metallicity

Current massive star evolution grids are not able to simultaneously reproduce the empirical upper luminosity limit of red supergiants, the Humphrey-Davidson (HD) limit at high and low metallicity. In this study, we provide a better understanding of what drives massive star evolution to blue and red supergiant phases, with the ultimate aim of reproducing the HD limit at varied metallicities (Z). For solar, LMC, and SMC Z, we develop eight grids of MESA models for the mass range 20-60M to probe the effect of semiconvection and overshooting. We compare rotating and non-rotating models with efficient (alpha_semi = 100) and inefficient semi-convection (alpha_semi = 0.1), with high and low core overshooting (alpha_ov of 0.1 or 0.5). The red and blue supergiant evolutionary phases are investigated by comparing the fraction of core He-burning lifetimes spent in each phase. We find that the extension of the convective core by overshooting alpha_ov = 0.5 has an effect on the post-MS evolution which can disable semiconvection leading to more RSGs, but a lack of BSGs. We therefore implement alpha_ov = 0.1 which switches on semiconvective mixing, though for standard alpha_semi = 1, would result in an HD limit which is higher than observed at low Z. Therefore, we need to implement very efficient semiconvection of alpha_semi = 100 which reproduces the HD limit at log L ~ 5.5 for the Magellanic Clouds while simultaneously reproducing the Galactic HD limit of log L ~ 5.8 naturally. The effect of semiconvection is not active at high Z due to the depletion of the envelope structure by strong mass loss such that semiconvective regions could not form. Z-dependent mass loss plays an indirect, yet decisive role in setting the HD limit as a function of Z. For a combination of efficient semiconvection and low overshooting with standard Z-dependent mass loss, we find a natural HD limit at all metallicities.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Star Formation in the Ultraviolet

With the launch of JWST and the upcoming installation of extremely large telescopes, the first galaxies in our Universe will finally be revealed. Their light will be dominated by massive stars, which peak in in the ultra-violet (UV) part of the electromagnetic spectrum. Star formation is the key driver of the evolution of our Universe. At young ages, within 10 Million years, both high and low mass stars generate complex UV emission processes which are poorly understood yet are vital for interpreting high red-shift line emission. For these reasons, the Hubble Space Telescope (HST) will devote 1000 orbits to obtaining a UV Legacy Library of Young Stars as Essential Standards (ULLYSES). The purpose of this Overview is to outline the basic physical principles driving UV emission processes from local (within 100 parsecs of) star formation, ranging from huge star-forming complexes containing hundreds of massive and very-massive stars (VMS), such as 30 Doradus (the Tarantula Nebula) in the neighbouring Magellanic Clouds (only 50 kpc away), to galaxies near and far, out to the epoch of Cosmic Reionization.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Constraining the progenitor evolution of GW 150914

One of the largest surprises from the LIGO results regarding the first gravitational wave detection (GW 150914) was the fact the black holes (BHs) were "heavy", of order 30 - 40 solar masses. The most promising explanation for this obesity is that the BH-BH merger occurred at low metallicity (Z): when the iron (Fe) contents is lower this is expected to result in weaker mass loss during the Wolf-Rayet (WR) phase. We therefore critically evaluate the claims for the reasons of heavy BHs as a function of Z in the literature. Furthermore, weaker stellar winds might lead to more rapid stellar rotation, allowing WR and BH progenitor evolution in a chemically homogeneous manner. However, there is as yet no empirical evidence for more rapid rotation amongst WR stars in the low Z environment of the Magellanic Clouds. Due to the intrinsic challenge of determining WR rotation rates from emission lines, the most promising avenue to constrain rotation-rate distributions amongst various WR subgroups is through the utilisation of their emission lines in polarised light. We thus provide an overview of linear spectro-polarimetry observations of both single and binary WRs in the Galaxy, as well as the Large and Small Magellanic Clouds, at 50% and 20% of solar Z, respectively. Initial results suggest that the route of chemically homogeneous evolution (CHE) through stellar rotation is challenging, whilst the alternative of a post-LBV or common envelope evolution is more likely.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Massive star evolution revealed in the Mass-Luminosity plane

Massive star evolution is dominated by key physical processes such as mass loss, convection and rotation, yet these effects are poorly constrained, even on the main sequence. We utilise a detached, eclipsing binary HD166734 as a testbed for single star evolution to calibrate new MESA stellar evolution grids. We introduce a novel method of comparing theoretical models with observations in the 'Mass-Luminosity Plane', as an equivalent to the HRD (see Higgins & Vink, 2018). We reproduce stellar parameters and abundances of HD166734 with enhanced overshooting (alpha=0.5), mass loss and rotational mixing. When comparing the constraints of our testbed to the systematic grid of models we find that a higher value of alpha = 0.5 (rather than 0.1) results in a solution which is more likely to evolve to a neutron star than a black hole, due to a lower value of the compactness parameter.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Polarimetry as a tool to study multi-dimensional winds and disks

I start with a discussion of spherical winds and small-scale clumping, before continuing with various theories that have been proposed to predict how mass loss depends on stellar rotation -- both in terms of wind strength, as well as the latitudinal dependence of the wind. This very issue is crucial for our general understanding of angular momentum evolution in massive stars, and the B[e] phenomenon in particular. I then discuss the tool of linear polarimetry that allows us to probe the difference between polar and equatorial mass loss, allowing us to test B[e] and related disk formation theories.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Two bi-stability jumps in theoretical wind models for massive stars and the implications for Luminous Blue Variable supernovae

Luminous Blue Variables have been suggested to be the direct progenitors of supernova types IIb and IIn, with enhanced mass loss prior to explosion. However, the mechanism of this mass loss is not yet known. Here, we investigate the qualitative behaviour of theoretical stellar wind mass-loss as a function of Teff across two bi-stability jumps in blue supergiant regime and also in proximity to the Eddington limit, relevant for LBVs. To investigate the physical ingredients that play a role in the radiative acceleration we calculate blue supergiant wind models with the CMFGEN non-LTE model atmosphere code over an effective temperature range between 30000 and 8800 K. Although our aim is not to provide new mass-loss rates for BA supergiants, we study and confirm the existence of two bi-stability jumps in mass-loss rates predicted by Vink, de Koter, & Lamers (1999). However, they are found to occur at somewhat lower Teff (20000 and 9000 K, respectively) than found previously, which would imply that stars may evolve towards lower Teff before strong mass-loss is induced by the bi-stability jumps. When the combined effects of the second bi-stability jump and the proximity to Eddington limit are accounted for, we find a dramatic increase in the mass-loss rate by up to a factor of 30. Further investigation of both bi-stability jumps is expected to lead to a better understanding of discrepancies between empirical modelling and theoretical mass-loss rates reported in the literature, and to provide key inputs for the evolution of both normal AB supergiants and LBVs, as well as their subsequent supernova type II explosions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Linear line spectropolarimetry as a new window to measure 2D and 3D wind geometries

Various theories have been proposed to predict how mass loss depends on the stellar rotation rate, both in terms of its strength, as well as its latitudinal dependence, crucial for our understanding of angular momentum evolution. Here we discuss the tool of linear spectropolarimetry that can probe the difference between mass loss from the pole versus the equator. Our results involve several groups of O stars and Wolf-Rayet stars, involving Oe stars, Of?p stars, Onfp stars, as well as the best candidate gamma-ray burst progenitors identified to date.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Linear line spectropolarimetry of Herbig Ae/Be stars

Accretion is the prime mode of star formation, but the exact mode has not yet been identified in the Herbig Ae/Be mass range. We provide evidence that the the maximum variation in mass-accretion rate is reached on a rotational timescale, which suggests that rotational modulation is the key to understanding mass accretion. We show how spectropolarimetry is uniquely capable of resolving the innermost (within 0.1 AU) regions between the star and the disk, allowing us to map the 3D geometry of the accreting gas, and test theories of angular momentum evolution. We present Monte Carlo line-emission simulations showing how one would observe changes in the polarisation properties on rotational timescales, as accretion columns come and go into our line of sight.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The True origin of Wolf-Rayet stars

The Wolf-Rayet (WR) phenomenon is widespread in astronomy. It involves classical WRs, very massive stars (VMS), WR central stars of planetary nebula CSPN [WRs], and supernovae (SNe). But what is the root cause for a certain type of object to turn into an emission-line star? In this contribution, I discuss the basic aspects of radiation-driven winds that might reveal the ultimate difference between WR stars and canonical O-type stars. I discuss the aspects of (i) self-enrichment via CNO elements, (ii) high effective temperatures Teff, (iii) an increase in the helium abundance Y, and finally (iv) the Eddington factor Gamma. Over the last couple of years, we have made a breakthrough in our understanding of Gamma-dependent mass loss, which will have far-reaching consequences for the evolution and fate of the most massive stars in the Universe. Finally, I discuss the prospects for studies of the WR phenomenon in the highest redshift Ly-alpha and He II emitting galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Internal entrainment and the origin of jet-related broad-band emission in Centaurus A

The dimensions of Fanaroff-Riley class I jets and the stellar densities at galactic centres imply that there will be numerous interactions between the jet and stellar winds. These may give rise to the observed diffuse and 'knotty' structure of the jets in the X-ray, and can also mass load the jets. We performed modelling of internal entrainment from stars intercepted by Centaurus A's jet, using stellar evolution- and wind codes. From photometry and a code-synthesised population of 12 Gyr (Z = 0.004), 3 Gyr (Z = 0.008) and 0 - 60 Myr (Z = 0.02) stars, appropriate for the parent elliptical NGC 5128, the total number of stars in the jet is ~ 8 x 10^8. Our model is energetically capable of producing the observed X-ray emission, even without young stars. We also reproduce the radio through X-ray spectrum of the jet, albeit in a downstream region with distinctly fewer young stars, and recover the mean X-ray spectral index. We derive an internal entrainment rate of ~ 2.3 x 10^-3 Msun yr^-1 which implies substantial jet deceleration. Our absolute nucleosynthetic yields for the AGB stellar population in the jet show the highest amounts for 4He, 16O, 12C, 14N and 20Ne. If some of the events at >= 55 EeV detected by the Pierre Auger Observatory originate from internal entrainment in Centaurus A, we predict that their composition will be largely intermediate-mass nuclei with 16O, 12C and 14N the key isotopes.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Luminous Blue Variables and superluminous supernovae from binary mergers

Evidence suggests that the direct progenitor stars of some core-collapse supernovae (CCSNe) are luminous blue variables (LBVs), perhaps including some `superluminous supernovae' (SLSNe). We examine models in which massive stars gain mass soon after the end of core hydrogen burning. These are mainly intended to represent mergers following a brief contact phase during early Case B mass transfer, but may also represent stars which gain mass in the Hertzsprung Gap or extremely late during the main-sequence phase for other reasons. The post-accretion stars spend their core helium-burning phase as blue supergiants (BSGs), and many examples are consistent with being LBVs at the time of core collapse. Other examples are yellow supergiants at explosion. We also investigate whether such post-accretion stars may explode successfully after core collapse. The final core properties of post-accretion models are broadly similar to those of single stars with the same initial mass as the pre-merger primary star. More surprisingly, when early Case B accretion does affect the final core properties, the effect appears likely to favour a successful SN explosion, i.e., to make the core properties more like those of a lower-mass single star. However, the detailed structures of these cores sometimes display qualitative differences to any single-star model we have calculated. The rate of appropriate binary mergers may match the rate of SNe with immediate LBV progenitors; for moderately optimistic assumptions we estimate that the progenitor birthrate is ~1% of the CCSN rate.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the H$α$ Behaviour of Blue Supergiants: Rise and Fall over the Bi-stability Jump

The evolutionary state of blue supergiants is still unknown. Stellar wind mass loss is one of the dominant processes determining the evolution of massive stars, and it may provide clues on the evolutionary properties of blue supergiants. As the H$α$ line is the most oft-used mass-loss tracer in the OB-star regime, we provide a detailed analysis of the H$α$ line for OB supergiant models over an $T_{\rm eff}$ range between 30000 and 12500K. We find a maximum in the H$α$ equivalent width at 22500 K - at the location of the bi-stability jump. The H$α$ line-profile behaviour is characterised by two branches of $T_{\rm eff}$: (i) a "hot" branch between 30000 and 22500 K, where H$α$ emission becomes stronger with decreasing $T_{\rm eff}$, and (ii) a "cool" branch between 22500 and 12500 K, where the line becomes weaker. Our models show that this non-monotonic H$α$ behaviour is related to the optical depth of Ly$α$, finding that at the "cool" branch the population of the 2nd level of hydrogen is enhanced in comparison to the 3rd level. This is expected to increase line absorption, leading to weaker H$α$ flux when $T_{\rm eff}$ drops from 22500 K downwards. We also show that for late B supergiants (at $T_{\rm eff}$ below ~15000 K), the differences in the H$α$ line between homogeneous and clumpy winds becomes insignificant. Moreover, we show that at the bi-stability jump H$α$ changes its character completely, from an optically thin to an optically thick line, implying that macro-clumping should play an important role at temperatures below the bi-stability jump. This would not only have consequences for the character of observed H$α$ line profiles, but also for the reported discrepancies between theoretical and empirical mass-loss rates.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the missing second generation AGB stars in NGC6752

[Abridged] In recent years the view of Galactic globular clusters as simple stellar populations has changed dramatically, as it is now thought that basically all GCs host multiple stellar populations, each with its own chemical abundance pattern and colour-magnitude diagram sequence. Recent spectroscopic observations of asymptotic giant branch stars in the GC NGC6752 have disclosed a low [Na/Fe] abundance for the whole sample, suggesting that they are all first-generation stars, and that all second-generation stars fail to reach the AGB in this cluster. A scenario proposed to explain these observations invokes strong mass loss in second-generation horizontal branch stars possibly induced by the metal enhancement associated to radiative levitation. This enhanced mass loss would prevent second generation stars from reaching the AGB phase, thus explaining at the same time the low value of the ratio between HB and AGB stars (the R_2 parameter) observed in NGC6752. We have critically discussed this scenario, finding that the required mass-loss rates are of the order of 10^{-9} Mo/yr, significantly higher than current theoretical and empirical constraints. By making use of synthetic HB simulations, we demonstrate that our modelling predicts correctly the R_2 parameter for NGC6752, without the need to invoke very efficient mass loss during the core He-burning stage. Our simulations for NGC6752 HB predict however the presence of a significant fraction - at the level of about 50% - second generation stars along the cluster AGB. We conclude that there is no simple explanation for the lack of second generation stars in the spectroscopically surveyed sample, although the interplay between mass loss (with low rates) and radiative levitation may play a role in explaining this puzzle.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Rotating massive O stars with non-spherical 2D winds

We present solutions for the velocity field and mass-loss rates for 2D axisymmetric outflows, as well as for the case of mass accretion through the use of the Lambert W-function. For the case of a rotating radiation-driven wind the velocity field is obtained analytically using a parameterised description of the line acceleration that only depends on radius r at any given latitude $θ$. The line acceleration g(r) is obtained from Monte-Carlo multi-line radiative transfer calculations. The critical/sonic point of our equation of motion varies with latitude $θ$. Furthermore, an approximate analytical solution for the supersonic flow of a rotating wind is derived, which is found to closely resemble the exact solution. For the simultaneous solution of the mass-loss rate and velocity field, we use the iterative method of our 1D method extended to the non-spherical 2D case. We apply the new theoretical expressions with our iterative method to the stellar wind from a differentially rotating 40 $M_{sun}$ O5-V main sequence star as well as to a 60 $M_{sun}$ O-giant star, and we compare our results to previous studies that are extensions of the Castor et al. (1975, ApJ, 195, 157) CAK formalism. Next, we account for the effects of oblateness and gravity darkening. Our numerical results predict an equatorial decrease of the mass-loss rate, which would imply that (surface-averaged) total mass-loss rates are lower than for the spherical 1D case, in contradiction to the Maeder & Meynet (2000, A&A, 361, 159) formalism that is oftentimes employed in stellar evolution calculations for rotating massive stars. To clarify the situation in nature we discuss observational tests to constrain the shapes of large-scale 2D stellar winds.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Temperaments of young stars: Rapid mass-accretion rate changes in T Tauri and Herbig Ae stars

Variability in emission lines is a characteristic feature in young stars and can be used as a tool to study the physics of the accretion process. Here we present a study of Hα variability in 15 T Tauri and Herbig Ae stars (K7-B2) over a wide range of time windows, from minutes, to hours, to days, and years. We assess the variability using linewidth measurements and the time series of line profiles. All objects show gradual, slow profile changes on time-scales of days. In addition, in three cases there is evidence for rapid variations in Hα with typical time-scales of 10 min, which occurs in 10% of the total covered observing time. The mean accretion-rate changes, inferred from the line fluxes,are 0.01-0.07 dex for time-scales of < 1 hour, 0.04-0.4 dex for time-scales of days, and 0.13-0.52 dex for time-scales of years. In Costigan et al. 2012 we derived an upper limit finding that the intermediate (days) variability dominated over longer (years) variability. Here our new results, based on much higher cadence observations, also provide a lower limit to accretion-rate variability on similar time-scales (days), thereby constraining the accretion rate variability physics in a much more definitive way. A plausible explanation for the gradual variations over days is an asymmetric accretion flow resulting in a rotational modulation of the accretion-related emission, although other interpretations are possible as well. In conjunction with our previous work, we find that the time-scales and the extent of the variability is similar for objects ranging in mass from ~ 0.1 to ~ 5 Msol. This confirms that a single mode of accretion is at work from T Tauri to Herbig Ae stars - across a wide range of stellar masses.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Very Massive Stars in the Local Universe

Recent studies suggest the existence of very massive stars (VMS) up to 300 solar masses in the local Universe. As this finding may represent a paradigm shift for the canonical stellar upper-mass limit of 150 solar masses, it is timely to evaluate the physics specific to VMS, which is currently missing. For this reason, we decided to construct a book entailing both a discussion of the accuracy of VMS masses (Martins), as well as the physics of VMS formation (Krumholz), mass loss (Vink), instabilities (Owocki), evolution (Hirschi), and fate (theory -- Woosley & Heger; observations -- Smith).

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Very Massive Stars (VMS) in the Local Universe

Recent studies have claimed the existence of very massive stars (VMS) up to 300 solar masses in the local Universe. As this finding may represent a paradigm shift for the canonical stellar upper-mass limit of 150 Msun, it is timely to discuss the status of the data, as well as the far-reaching implications of such objects. We held a Joint Discussion at the General Assembly in Beijing to discuss (i) the determination of the current masses of the most massive stars, (ii) the formation of VMS, (iii) their mass loss, and (iv) their evolution and final fate. The prime aim was to reach broad consensus between observers and theorists on how to identify and quantify the dominant physical processes.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Bayesian inference of T Tauri star properties using multi-wavelength survey photometry

There are many pertinent open issues in the area of star and planet formation. Large statistical samples of young stars across star-forming regions are needed to trigger a breakthrough in our understanding, but most optical studies are based on a wide variety of spectrographs and analysis methods, which introduces large biases. Here we show how graphical Bayesian networks can be employed to construct a hierarchical probabilistic model which allows pre-main sequence ages, masses, accretion rates, and extinctions to be estimated using two widely available photometric survey databases (IPHAS r/i/Halpha and 2MASS J-band magnitudes.) Because our approach does not rely on spectroscopy, it can easily be applied to homogeneously study the large number of clusters for which Gaia will yield membership lists. We explain how the analysis is carried out using the Markov Chain Monte Carlo (MCMC) method and provide Python source code. We then demonstrate its use on 587 known low-mass members of the star-forming region NGC 2264 (Cone Nebula), arriving at a median age of 3.0 Myr, an accretion fraction of 20+/-2% and a median accretion rate of 10^-8.4 Msol/yr. The Bayesian analysis formulated in this work delivers results which are in agreement with spectroscopic studies already in the literature, but achieves this with great efficiency by depending only on photometry. It is a significant step forward from previous photometric studies, because the probabilistic approach ensures that nuisance parameters, such as extinction and distance, are fully included in the analysis with a clear picture on any degeneracies.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Linear Spectropolarimetry and the Circumstellar Media of Young and Massive Stars

Linear spectropolarimetry is a powerful tool to probe circumstellar structures on spatial scales that cannot yet be achieved through direct imaging. In this review I discuss the role that emission-line polarimetry can play in constraining geometrical and physical properties of a wide range of circumstellar environments, varying from the accretion disks around pre-main sequence T Tauri and Herbig Ae/Be stars, to the issue of stellar wind clumping, and the aspherical outflows from the massive star progenitors of supernovae and long gamma-ray bursts at low metallicity.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Mass loss and fate of the most massive stars

The fate of massive stars up to 300 Msun is highly uncertain. Do these objects produce pair-instability explosions, or normal Type Ic supernovae? In order to address these questions, we need to know their mass-loss rates during their lives. Here we present mass-loss predictions for very massive stars (VMS) in the range of 60-300 Msun. We use a novel method that simultaneously predicts the wind terminal velocities (vinf) and mass-loss rate (dM/dt) as a function of the stellar parameters: (i) luminosity/mass Gamma, (ii) metallicity Z, and (iii) effective temperature Teff. Using our results, we evaluate the likely outcomes for the most massive stars.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The transition mass-loss rate: Calibrating the role of line-driven winds in massive star evolution

A debate has arisen regarding the importance of stationary versus eruptive mass loss for massive star evolution. The reason is that stellar winds have been found to be clumped, which results in the reduction of unclumped empirical mass-loss rates. Most stellar evolution models employ theoretical mass-loss rates which are already reduced by a moderate factor of ~2-3 compared to non-corrected empirical rates. A key question is whether these reduced rates are of the correct order of magnitude, or if they should be reduced even further, which would mean that the alternative of eruptive mass loss becomes necessary. Here we introduce the transition mass-loss rate (dM/dt)_trans between O and Wolf-Rayet (WR) stars. Its novelty is that it is model independent. All that is required is postulating the spectroscopic transition point in a given data-set, and determining the stellar luminosity, which is far less model dependent than the mass-loss rate. The transition mass-loss rate is subsequently used to calibrate stellar wind strength by its application to the Of/WNh stars in the Arches cluster. Good agreement is found with two alternative modelling/theoretical results, suggesting that the rates provided by current theoretical models are of the right order of magnitude in the ~50Msun mass range. Our results do not confirm the specific need for eruptive mass loss as Luminous Blue Variables, and current stellar evolution modelling for Galactic massive stars seems sound. Mass loss through alternative mechanisms might still become necessary at lower masses, and/or metallicities, and the quantification of alternative mass loss is desirable.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Chandra observations of Cygnus OB2

Cygnus OB2 is the nearest example of a massive star forming region, containing over 50 O-type stars and hundreds of B-type stars. We have analyzed two Chandra pointings in Cyg OB2, detecting ~1700 X-ray sources, of which ~1450 are thought to be members of the association. Optical and near-IR photometry has been obtained for ~90% of these sources from recent deep Galactic plane surveys. We have performed isochrone fits to the near-IR color-magnitude diagram, deriving ages of 3.5(+0.75,-1.0) and 5.25(+1.5,-1.0) Myrs for sources in the two fields, both with considerable spreads around the pre-MS isochrones. The presence of a second population in the region, somewhat older than the present-day O-type stars, has been suggested by other authors and fits with the ages derived here. The fraction of sources with inner circumstellar disks (as traced by the K-band excess) is found to be very low, but appropriate for a population of age ~5 Myrs. We measure the stellar mass functions and find a power-law slope of Gamma = -1.09 +/- 0.13, in good agreement with the global mean value estimated by Kroupa. A steepening of the mass function at high masses is observed and we suggest this is due to the presence of the previous generation of stars that have lost their most massive members. Finally, combining our mass function and an estimate of the radial density profile of the association suggests a total mass of Cyg OB2 of ~30,000 Msun, similar to that of many of our Galaxy's most massive star forming regions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Eta Carinae and the Luminous Blue Variables

We evaluate the place of Eta Carinae amongst the class of luminous blue variables (LBVs) and show that the LBV phenomenon is not restricted to extremely luminous objects like Eta Car, but extends luminosities as low as log(L/Lsun) = 5.4 - corresponding to initial masses ~25 Msun, and final masses as low as ~10-15 Msun. We present a census of S Doradus variability, and discuss basic LBV properties, their mass-loss behaviour, and whether at maximum light they form pseudo-photospheres. We argue that those objects that exhibit giant Eta Car-type eruptions are most likely related to the more common type of S Doradus variability. Alternative atmospheric models as well as sub-photospheric models for the instability are presented, but the true nature of the LBV phenomenon remains as yet elusive. We end with a discussion on the evolutionary status of LBVs - highlighting recent indications that some LBVs may be in a direct pre-supernova state, in contradiction to the standard paradigm for massive star evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

In pursuit of gamma-ray burst progenitors: the identification of a sub-population of rotating Wolf-Rayet stars

Long gamma-ray bursts involve the most powerful cosmic explosions since the Big Bang. Whilst it has been established that GRBs are related to the death throes of massive stars, the identification of their progenitors has proved challenging. Theory suggests that rotating Wolf-Rayet stars are the best candidates, but their strong stellar winds shroud their surfaces, preventing a direct measurement of their rotation. Fortunately, linear spectropolarimetry may be used to probe the flattening of their winds due to stellar spin. Spectropolarimetry surveys show that an 80% majority of WR stars have spherically symmetric winds and are thus rotating slowly, yet a small 20% minority display a spectropolarimetric signature indicative of rotation. Here we find a highly significant correlation between WR objects that carry the signature of stellar rotation and the subset of WR stars with ejecta nebulae that have only recently transitioned from a red sugergiant or luminous blue variable phase. As these youthful WR stars have yet to spin-down due to mass loss, they are the best candidate GRB progenitors identified to date. When we take recently published WR ejecta nebula numbers we find that five out of the six line-effect WR stars are surrounded by ejecta nebulae. The statistics imply that the null hypothesis of no correlation between line-effect WR stars and ejecta nebulae can be rejected at the 0.0004% level. Given that four line-effect and WR ejecta nebula have spectroscopically been confirmed to contain nucleosynthetic products, we argue that the correlation is both statistically significant and physically convincing. The implication is that we have identified a WR sub-population that fulfills the necessary criteria for making GRBs. Finally, we discuss the potential of identifying GRB progenitors via spectropolarimetry with extremely large telescopes.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Rotating Massive Main-Sequence Stars I: Grids of Evolutionary Models and Isochrones

We present a dense grid of evolutionary tracks and isochrones of rotating massive main-sequence stars. We provide three grids with different initial compositions tailored to compare with early OB stars in the Small and Large Magellanic Clouds and in the Galaxy. Each grid covers masses ranging from 5 to 60 Msun and initial rotation rates between 0 and about 600 km/s. To calibrate our models we used the results of the VLT-FLAMES Survey of Massive Stars. We determine the amount of convective overshooting by using the observed drop in rotation rates for stars with surface gravities log g < 3.2 to determine the width of the main sequence. We calibrate the efficiency of rotationally induced mixing using the nitrogen abundance determinations for B stars in the Large Magellanic cloud. We describe and provide evolutionary tracks and the evolution of the central and surface abundances. In particular, we discuss the occurrence of quasi-chemically homogeneous evolution, i.e. the severe effects of efficient mixing of the stellar interior found for the most massive fast rotators. We provide a detailed set of isochrones for rotating stars. Rotation as an initial parameter leads to a degeneracy between the age and the mass of massive main sequence stars if determined from its observed location in the Hertzsprung-Russell diagram. We show that the consideration of surface abundances can resolve this degeneracy.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

T Tauri candidates and accretion rates using IPHAS: method and application to IC1396

The INT Photometric H-Alpha Survey (IPHAS) is a 1800 deg^2 survey of the Northern Galactic Plane, reaching down to r'~21. We demonstrate how the survey can be used to (1) reliably select classical T Tauri star candidates and (2) constrain the mass accretion rates with an estimated relative uncertainty of 0.6 dex. IPHAS is a necessary addition to spectroscopic surveys because it allows large and uniform samples of accretion rates to be obtained with a precise handle on the selection effects. We apply the method on a region of 7 deg^2 towards the HII region IC1396 in Cepheus OB2 and identify 158 pre-main sequence candidates with masses between 0.2 and 2.0 Msun and accretion rates between 10^-9.2 and 10^-7.0 Msun yr^-1. We find a power-law dependency between the stellar mass and the accretion rates with a slope of α= 1.1 +/- 0.2, which is less steep than indicated by previous studies. We discuss the influence of method-dependent systematic effects on the determination of this relationship. The majority of our sample consists of faint, previously unknown, low-mass T Tauri candidates (56 per cent between 0.2 and 0.5 Msun). Many candidates are clustered in front of three bright-rimmed molecular clouds, which are being ionized by the massive star HD206267 (O6.5V). We discover a spatio-temporal gradient of increasing accretion rates, increasing Spitzer infrared excess, and younger ages away from the ionizing star, providing a strong indication that the formation of these clusters has been sequentially triggered by HD206267 during the last ~1 Myr.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The bi-stability jump as the origin for multiple P-Cygni absorption components in Luminous Blue Variables

Luminous Blue Variables (LBVs) oftentimes show double-troughed absorption in their strong Halpha lines, which are as yet not understood. Intriguingly, the feature has also been seen in the interacting supernova SN 2005gj, which was for this reason suggested to have an LBV progenitor. Our aims are to understand the double-troughed absorption feature in LBVs and investigate whether this phenomenon is related to wind variability. To this purpose, we perform time-dependent radiative transfer modeling using CMFGEN. We find that abrupt changes in the wind-terminal velocity - as expected from the bi-stability jump - are required to explain the double-troughed absorption profiles in LBVs. This strengthens scenarios that discuss the link between LBVs and SNe utilizing the progenitor's wind variability resulting from the bi-stability jump. We also discuss why the presence of double-troughed P-Cygni components may become an efficient tool to detect extra-galactic LBVs and how to analyze their mass-loss history on the basis of just one single epoch of spectral observations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The theory of stellar winds

We present a brief overview of the theory of stellar winds with a strong emphasis on the radiation-driven outflows from massive stars. The resulting implications for the evolution and fate of massive stars are also discussed. Furthermore, we relate the effects of mass loss to the angular momentum evolution, which is particularly relevant for the production of long and soft gamma-ray bursts. Mass-loss rates are not only a function of the metallicity, but are also found to depend on temperature, particularly in the region of the bi-stability jump at 21 000 Kelvin. We highlight the role of the bi-stability jump for Luminous Blue Variable (LBV) stars, and discuss suggestions that LBVs might be direct progenitors of supernovae. We emphasize that radiation-driven wind studies rely heavily on the input opacity data and linelists, and that these are thus of fundamental importance to both the mass-loss predictions themselves, as well as to our overall understanding of the lives and deaths of massive stars.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The VLT-FLAMES Tarantula Survey III: A very massive star in apparent isolation from the massive cluster R136

VFTS 682 is located in an active star-forming region, at a projected distance of 29 pc from the young massive cluster R136 in the Tarantula Nebula of the Large Magellanic Cloud. It was previously reported as a candidate young stellar object, and more recently spectroscopically revealed as a hydrogen-rich Wolf-Rayet (WN5h) star. Our aim is to obtain the stellar properties, such as its intrinsic luminosity, and to investigate the origin of VFTS 682. To this purpose, we model optical spectra from the VLT-FLAMES Tarantula Survey with the non-LTE stellar atmosphere code CMFGEN, as well as the spectral energy distribution from complementary optical and infrared photometry. We find the extinction properties to be highly peculiar (RV ~4.7), and obtain a surprisingly high luminosity log(L/Lsun) = 6.5 \pm 0.2, corresponding to a present-day mass of ~150Msun. The high effective temperature of 52.2 \pm 2.5kK might be explained by chemically homogeneous evolution - suggested to be the key process in the path towards long gamma-ray bursts. Lightcurves of the object show variability at the 10% level on a timescale of years. Such changes are unprecedented for classical Wolf-Rayet stars, and are more reminiscent of Luminous Blue Variables. Finally, we discuss two possibilities for the origin of VFTS 682: (i) the star either formed in situ, which would have profound implications for the formation mechanism of massive stars, or (ii) VFTS 682 is a slow runaway star that originated from the dense cluster R136, which would make it the most massive runaway known to date.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Wind modelling of very massive stars up to 300 solar masses

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Wind Models for Very Massive Stars in the Local Universe

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

On The Maximum Mass of Stellar Black Holes

We present the spectrum of compact object masses: neutron stars and black holes that originate from single stars in different environments. In particular, we calculate the dependence of maximum black hole mass on metallicity and on some specific wind mass loss rates (e.g., Hurley et al. and Vink et al.). Our calculations show that the highest mass black holes observed in the Galaxy M_bh = 15 Msun in the high metallicity environment (Z=Zsun=0.02) can be explained with stellar models and the wind mass loss rates adopted here. To reach this result we had to set Luminous Blue Variable mass loss rates at the level of about 0.0001 Msun/yr and to employ metallicity dependent Wolf-Rayet winds. With such winds, calibrated on Galactic black hole mass measurements, the maximum black hole mass obtained for moderate metallicity (Z=0.3 Zsun=0.006) is M_bh,max = 30 Msun. This is a rather striking finding as the mass of the most massive known stellar black hole is M_bh = 23-34 Msun and, in fact, it is located in a small star forming galaxy with moderate metallicity. We find that in the very low (globular cluster-like) metallicity environment the maximum black hole mass can be as high as M_bh,max = 80 Msun (Z=0.01 Zsun=0.0002). It is interesting to note that X-ray luminosity from Eddington limited accretion onto an 80 Msun black hole is of the order of about 10^40 erg/s and is comparable to luminosities of some known ULXs. We emphasize that our results were obtained for single stars only and that binary interactions may alter these maximum black hole masses (e.g., accretion from a close companion). This is strictly a proof-of-principle study which demonstrates that stellar models can naturally explain even the most massive known stellar black holes.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Massive Star Forming Region Cygnus OB2. II. Integrated Stellar Properties and the Star Formation History

Cygnus OB2 is the nearest example of a massive star forming region, containing over 50 O-type stars and hundreds of B-type stars. We have analysed the properties of young stars in two fields in Cyg OB2 using the recently published deep catalogue of Chandra X-ray point sources with complementary optical and near-IR photometry. Our sample is complete to 1 Msun (excluding A and B-type stars that do not emit X-rays), making this the deepest study of the stellar properties and star formation history in Cyg OB2 to date. From Siess et al. (2000) isochrone fits to the near-IR color-magnitude diagram, we derive ages of 3.5 (+0.75/-1.0) and 5.25 (+1.5/-1.0) Myrs for sources in the two fields, both with considerable spreads around the pre-MS isochrones. The presence of a stellar population somewhat older than the present-day O-type stars, also fits in with the low fraction of sources with inner circumstellar disks (as traced by the K-band excess) that we find to be very low, but appropriate for a population of age ~5 Myrs. We also find that the region lacks a population of highly embedded sources that is often observed in young star forming regions, suggesting star formation in the vicinity has declined. We measure the stellar mass functions in this limit and find a power-law slope of Gamma = -1.09 +/- 0.13, in good agreement with the global mean value estimated by Kroupa (2002). A steepening of the slope at higher masses is observed and suggested as due to the presence of the previous generation of stars that have lost their most massive members. Finally, combining our mass function and an estimate of the radial density profile of the association suggests a total mass of Cyg OB2 of ~3 x 10^4 Msun, similar to that of many of our Galaxy's most massive star forming regions.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

On the presence and absence of disks around O-type stars

(abridged) As the favoured progenitors of long-duration gamma-ray bursts, massive stars may represent our best signposts of individual objects in the early Universe, but special conditions seem required to make these bursters, which might originate from the progenitor's rapid rotation and associated asymmetry. To obtain empirical constraints on the interplay between stellar rotation and wind asymmetry, we perform linear Halpha spectropolarimetry on a sample of 18 spectroscopically peculiar massive O stars, including OVz, Of?p, Oe, and Onfp stars, supplemented by an earlier sample of 20 O supergiants. Despite their rapid rotation (with vsin(i) up to 400 km/s) most O-type stars are found to be spherically symmetric, but with notable exceptions amongst specific object classes. We divide the peculiar O stars into four distinct categories: Groups III and IV include the Oe stars and Onfp stars, which are on the high-end tail of the O star rotation distribution and have in the past been claimed to be embedded in disks. Here we report the detection of a classical depolarization ``line effect'' in the Oe star HD 45314, but the overall incidence of line effects amongst Oe stars is significantly lower (1 out of 6) than amongst Be stars. The chance that the Oe and Be datasets are drawn from the same parent population is negligible (with 95% confidence). This implies there is as yet no evidence for a disk hypothesis in Oe stars, providing relevant constraints on the physical mechanism that is responsible for the Be phenomenon. Finally, we find that 3 out of 4 of the group IV Onfp stars show evidence for complex polarization effects, which are likely related to rapid rotation, and we speculate on the evolutionary links to B[e] stars.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Modelling the clumping-induced polarimetric variability of hot star winds

Clumping in the winds of massive stars may significantly reduce empirical mass-loss rates, and which in turn may have a large impact on our understanding of massive star evolution. Here, we investigate wind-clumping through the linear polarization induced by light scattering off the clumps. Through the use of an analytic wind clumping model, we predict the time evolution of the linear polarimetry over a large parameter space. We concentrate on the Luminous Blue Variables, which display the greatest amount of polarimetric variability and for which we recently conducted a spectropolarimetric survey. Our model results indicate that the observed level of polarimetric variability can be reproduced for two regimes of parameter space: one of a small number of massive, optically-thick clumps; and one of a very large number of low-mass clumps. Although a systematic time-resolved monitoring campaign is required to distinguish between the two scenarios, we currently favour the latter, given the short timescale of the observed polarization variability. As the polarization is predicted to scale linearly with mass-loss rate, we anticipate that all hot stars with very large mass-loss rates should display polarimetric variability. This is consistent with recent findings that intrinsic polarization is more common in stars with strong H$α$ emission.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Probing the circumstellar structure of pre-main sequence stars

We present Halpha spectropolarimetry of a large sample of pre-main sequence (PMS) stars of low and intermediate mass, and argue that the technique is a powerful tool in studying the circumstellar geometry around these objects. For the intermediate mass (2 -- 15 Msun) Herbig Ae/Be stars we find that 16 out of 23 show a line effect, which immediately implies that flattening is common among these objects. Furthermore, we find a significant difference in Halpha spectropolarimetry behaviour between the Herbig Be and Ae groups. For the Herbig Be stars, the concept of an electron scattering disc is shown to be a useful concept to explain the depolarizations seen in this spectral range. At lower masses, more complex Halpha polarimetry behaviour starts to appear. The concept of a compact source of Halpha emission that is formed close to the stellar surface, for instance by hot spots due to magnetospheric accretion, is postulated as a working hypothesis to qualitatively explain the Halpha spectropolarimetry behaviour around Herbig Ae and lower mass (M < 2 Msun) T Tauri stars. The striking resemblance in spectropolarimetric behaviour between the T Tauri star RY Tau and the Herbig Ae stars suggests a common origin of the polarized line photons, and hints that low and higher mass pre-main sequence stars may have more in common than had hitherto been suspected.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

Hot Horizontal Branch Stars: Predictions for Mass Loss

We predict mass-loss rates for the late evolutionary phases of low-mass stars, with special emphasis on the consequences for the morphology of the Horizontal Branch (HB). We show that the computed rates, as predicted by the most plausible mechanism of radiation pressure on spectral lines, are too low to produce EHB/sdB stars. This invalidates the scenario recently outlined by Yong et al. (2000) to create these objects by mass loss on the HB. We argue, however, that mass loss plays a role in the distribution of rotational velocities of hot HB stars, and may -- together with the enhancement of heavy element abundances due to radiative levitation -- provide an explanation for the so-called ``low gravivity problem. The mass loss recipe derived for hot HB (and extreme HB, sdB, sdOB) stars may also be applied to post-HB (AGB-manque, UV-bright) stars over a range in effective temperatures between 12500 -- 40000 K.

0 citations0 reviewsNo code linkedOpen access