Context. Massive stars play a vital role in the Universe, however, their evolution even on the main-sequence is not yet well understood. Aims. Because of the steep mass-luminosity relation, massive ...main-sequence stars become extremely luminous. This brings their envelopes very close to the Eddington limit. We analyse stellar evolutionary models in which the Eddington limit is reached and exceeded, explore the rich diversity of physical phenomena that take place in their envelopes, and investigate their observational consequences. Methods. We use published grids of detailed stellar models, computed with a state-of-the-art, one-dimensional hydrodynamic stellar evolution code using LMC composition, to investigate the envelope properties of core hydrogen burning massive stars. Results. We find that the Eddington limit is almost never reached at the stellar surface, even for stars up to 500 M⊙. When we define an appropriate Eddington limit locally in the stellar envelope, we can show that most stars more massive than ~40 M⊙ actually exceed this limit, in particular, in the partial ionisation zones of iron, helium, or hydrogen. While most models adjust their structure such that the local Eddington limit is exceeded at most by a few per cent, our most extreme models do so by a factor of more than seven. We find that the local violation of the Eddington limit has severe consequences for the envelope structure, as it leads to envelope inflation, convection, density inversions, and, possibly to, pulsations. We find that all models with luminosities higher than 4 × 105L⊙, i.e. stars above ~40 M⊙ show inflation, with a radius increase of up to a factor of about 40. We find that the hot edge of the S Dor variability region coincides with a line beyond which our models are inflated by more than a factor of two, indicating a possible connection between S Dor variability and inflation. Furthermore, our coolest models show highly inflated envelopes with masses of up to several solar masses, and appear to be candidates for producing major luminous blue variable eruptions. Conclusions. Our models show that the Eddington limit is expected to be reached in all stars above ~40 M⊙ in the LMC, even in lower mass stars in the Galaxy, or in close binaries or rapid rotators. While our results do not support the idea of a direct super-Eddington wind driven by continuum photons, the consequences of the Eddington limit in the form of inflation, pulsations and possibly eruptions may well give rise to a significant enhancement of the time averaged mass-loss rate.
ABSTRACT
We present Herschel PACS spectroscopy of the O iii 88.4 μm fine-structure line in the nearby WC8+O binary system γ Vel to determine its oxygen abundance. The critical density of this line ...corresponds to several 105R* such that it is spatially extended in PACS observations at the 336 pc distance to γ Vel. Two approaches are used, the first involving a detailed stellar atmosphere analysis of γ Vel using cmfgen, extending to Ne ∼ 100 cm−3 in order to fully sample the line formation region of O iii 88.4 $\mu$m. The second approach involves the analytical model introduced by Barlow et al. and revised by Dessart et al., additionally exploiting ISO LWS spectroscopy of O iii 51.8 $\mu$m. We obtain higher luminosities for the WR and O components of γ Vel with respect to De Marco et al., log L/L⊙ = 5.31 and 5.56, respectively, primarily due to the revised (higher) interferometric distance. We obtain an oxygen mass fraction of XO = 1.0 ± 0.3 per cent for an outer wind volume filling factor of f = 0.5 ± 0.25, favouring either standard or slightly reduced Kunz et al. rates for the 12C(α, γ)16O reaction from comparison with BPASS binary population synthesis models. We also revisit neon and sulphur abundances in the outer wind of γ Vel from ISO SWS spectroscopy of S iv 10.5 $\mu$m, Ne ii 12.8 $\mu$m, and Ne iii 15.5 $\mu$m. The inferred neon abundance XNe = 2.0$_{-0.6}^{+0.4}$ per cent is in excellent agreement with BPASS predictions, while the sulphur abundance of XS = 0.04 ± 0.01 per cent agrees with the solar abundance, as expected for unprocessed elements.
The VLT-FLAMES Tarantula Survey Schneider, F. R. N.; Ramírez-Agudelo, O. H.; Tramper, F. ...
Astronomy and astrophysics (Berlin),
10/2018, Letnik:
618
Journal Article
Recenzirano
Odprti dostop
The 30 Doradus (30 Dor) nebula in the Large Magellanic Cloud (LMC) is the brightest HII region in the Local Group and a prototype starburst similar to those found in high redshift galaxies. It is ...thus a stepping stone to understand the complex formation processes of stars in starburst regions across the Universe. Here, we have studied the formation history of massive stars in 30 Dor using masses and ages derived for 452 mainly OB stars from the spectroscopic VLT-FLAMES Tarantula Survey (VFTS). We find that stars of all ages and masses are scattered throughout 30 Dor. This is remarkable because it implies that massive stars either moved large distances or formed independently over the whole field of view in relative isolation. We find that both channels contribute to the 30 Dor massive star population. Massive star formation rapidly accelerated about 8 Myr ago, first forming stars in the field before giving birth to the stellar populations in NGC 2060 and NGC 2070. The R136 star cluster in NGC 2070 formed last and, since then, about 1 Myr ago, star formation seems to be diminished with some continuing in the surroundings of R136. Massive stars within a projected distance of 8 pc of R136 are not coeval but show an age range of up to 6 Myr. Our mass distributions are well populated up to 200 M⊙. The inferred IMF is shallower than a Salpeter-like IMF and appears to be the same across 30 Dor. By comparing our sample of stars to stellar models in the Hertzsprung–Russell diagram, we find evidence for missing physics in the models above log L/L⊙ = 6 that is likely connected to enhanced wind mass loss for stars approaching the Eddington limit. Our work highlights the key information about the formation, evolution and final fates of massive stars encapsulated in the stellar content of 30 Dor, and sets a new benchmark for theories of massive star formation in giant molecular clouds.
The VLT-FLAMES Tarantula Survey Ramirez-Agudelo, O H; Sana, H; de Koter, A ...
Astronomy and astrophysics (Berlin),
4/2017, Letnik:
600
Journal Article
Recenzirano
Context. The Tarantula region in the Large Magellanic Cloud (LMC) contains the richest population of spatially resolved massive O-type stars known so far. This unmatched sample offers an opportunity ...to test models describing their main-sequence evolution and mass-loss properties. Aims. Using ground-based optical spectroscopy obtained in the framework of the VLT-FLAMES Tarantula Survey (VFTS), we aim to determine stellar, photospheric and wind properties of 72 presumably single O-type giants, bright giants and supergiants and to confront them with predictions of stellar evolution and of line-driven mass-loss theories. Methods. We apply an automated method for quantitative spectroscopic analysis of O stars combining the non-LTE stellar atmosphere model fastwind with the genetic fitting algorithm pikaia to determine the following stellar properties: effective temperature, surface gravity, mass-loss rate, helium abundance, and projected rotational velocity. The latter has been constrained without taking into account the contribution from macro-turbulent motions to the line broadening. Results. We present empirical effective temperature versus spectral subtype calibrations at LMC-metallicity for giants and supergiants. The calibration for giants shows a +1kK offset compared to similar Galactic calibrations; a shift of the same magnitude has been reported for dwarfs. The supergiant calibrations, though only based on a handful of stars, do not seem to indicate such an offset. The presence of a strong upturn at spectral type O3 and earlier can also not be confirmed by our data. In the spectroscopic and classical Hertzsprung-Russell diagrams, our sample O stars are found to occupy the region predicted to be the core hydrogen-burning phase by state-of-the-art models. For stars initially more massive than approximately 60M sub(middot in circle), the giant phase already appears relatively early on in the evolution; the supergiant phase develops later. Bright giants, however, are not systematically positioned between giants and supergiants at M sub(init)> or = 25M sub(middot in circle). At masses below 60M sub(middot in circle), the dwarf phase clearly precedes the giant and supergiant phases; however this behavior seems to break down at M sub(init)< or = 18M sub(middot in circle). Here, stars classified as late OIII and II stars occupy the region where O9.5-9.7 V stars are expected, but where few such late OV stars are actually seen. Though we can not exclude that these stars represent a physically distinct group, this behavior may reflect an intricacy in the luminosity classification at late O spectral subtype. Indeed, on the basis of a secondary classification criterion, the relative strength of Siiv to Hei absorption lines, these stars would have been assigned a luminosity class IV or V. Except for five stars, the helium abundance of our sample stars is in agreement with the initial LMC composition. This outcome is independent of their projected spin rates. The aforementioned five stars present moderate projected rotational velocities (i.e., nu sub(e) sini< 200kms super(-1)) and hence do not agree with current predictions of rotational mixing in main-sequence stars. They may potentially reveal other physics not included in the models such as binary-interaction effects. Adopting theoretical results for the wind velocity law, we find modified wind momenta for LMC stars that are ~0.3 dex higher than earlier results. For stars brighter than 10 super(5)L sub(middot in circle), that is, in the regime of strong stellar winds, the measured (unclumped) mass-loss rates could be considered to be in agreement with line-driven wind predictions if the clump volume filling factors were f sub(V)~ 1/8 to 1/6.
Mass loss by stellar wind is a key agent in the evolution and spectroscopic appearance of massive main sequence and post-main sequence stars. In Wolf–Rayet stars the winds can be so dense and so ...optically thick that the photosphere appears in the highly supersonic part of the outflow, veiling the underlying subsonic part of the star, and leaving the initial acceleration of the wind inaccessible to observations. Here we investigate the conditions and the structure of the subsonic part of the outflow of Galactic Wolf–Rayet stars, in particular of the WNE subclass; our focus is on the conditions at the sonic point of their winds. We compute 1D hydrodynamic stellar structure models for massive helium stars adopting outer boundaries at the sonic point. We find that the outflows of our models are accelerated to supersonic velocities by the radiative force from opacity bumps either at temperatures of the order of 200 kK by the iron opacity bump or of the order of 50 kK by the helium-II opacity bump. For a given mass-loss rate, the diffusion approximation for radiative energy transport allows us to define the temperature gradient based purely on the local thermodynamic conditions. For a given mass-loss rate, this implies that the conditions in the subsonic part of the outflow are independent from the detailed physical conditions in the supersonic part. Stellar atmosphere calculations can therefore adopt our hydrodynamic models as ab initio input for the subsonic structure. The close proximity to the Eddington limit at the sonic point allows us to construct a sonic HR diagram, relating the sonic point temperature to the luminosity-to-mass ratio and the stellar mass-loss rate, thereby constraining the sonic point conditions, the subsonic structure, and the stellar wind mass-loss rates of WNE stars from observations. The minimum stellar wind mass-loss rate necessary to have the flow accelerated to supersonic velocities by the iron opacity bump is derived. A comparison of the observed parameters of Galactic WNE stars to this minimum mass-loss rate indicates that these stars have their winds launched to supersonic velocities by the radiation pressure arising from the iron opacity bump. Conversely, stellar models which do not show transonic flows from the iron opacity bump form low-density extended envelopes. We derive an analytic criterion for the appearance of envelope inflation and of a density inversion in the outer sub-photospheric layers.
ABSTRACT
Upcoming large-scale spectroscopic surveys with e.g. WEAVE (William herschel telescope Enhanced Area Velocity Explorer) and 4MOST (4-metre Multi-Object Spectroscopic Telescope) will provide ...thousands of spectra of massive stars, which need to be analysed in an efficient and homogeneous way. Usually, studies of massive stars are limited to samples of a few hundred objects, which pushes current spectroscopic analysis tools to their limits because visual inspection is necessary to verify the spectroscopic fit. Often uncertainties are only estimated rather than derived and prior information cannot be incorporated without a Bayesian approach. In addition, uncertainties of stellar atmospheres and radiative transfer codes are not considered as a result of simplified, inaccurate, or incomplete/missing physics or, in short, idealized physical models. Here, we address the question of ‘How to compare an idealized model of complex objects to real data?’ with an empirical Bayesian approach and maximum a posteriori approximations. We focus on application to large-scale optical spectroscopic studies of complex astrophysical objects like stars. More specifically, we test and verify our methodology on samples of OB stars in 30 Doradus region of the Large Magellanic Clouds using a grid of fastwind model atmospheres. Our spectroscopic model de-idealization analysis pipeline takes advantage of the statistics that large samples provide by determining the model error to account for the idealized stellar atmosphere models, which are included into the error budget. The pipeline performs well over a wide parameter space and derives robust stellar parameters with representative uncertainties.
The VLT-FLAMES Tarantula Survey Apellaniz, J Maiz; Evans, C J; Barba, R H ...
Astronomy and astrophysics (Berlin),
04/2014, Letnik:
564
Journal Article
Recenzirano
Context. The commonly used extinction laws of Cardelli et al. (1989, ApJ, 345, 245) have limitations that, among other issues, hamper the determination of the effective temperatures of O and early B ...stars from optical and near-infrared (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 and 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 and 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.
We report the results of spectrophotometric observations of the massive star MN18 revealed via discovery of a bipolar nebula around it with the Spitzer Space Telescope. Using the optical spectrum ...obtained with the Southern African Large Telescope, we classify this star as B1 Ia. The evolved status of MN18 is supported by the detection of nitrogen overabundance in the nebula, which implies that it is composed of processed material ejected by the star. We analysed the spectrum of MN18 by using the code cmfgen, obtaining a stellar effective temperature of ≈21 kK. The star is highly reddened, E(B − V) ≈ 2 mag. Adopting an absolute visual magnitude of M
V = −6.8 ± 0.5 (typical of B1 supergiants), MN18 has a luminosity of log L/L⊙ ≈ 5.42 ± 0.30, a mass-loss rate of ≈(2.8-4.5) × 10− 7 M⊙ yr− 1, and resides at a distance of ≈5.6
$^{+1.5} _{-1.2}$
kpc. We discuss the origin of the nebula around MN18 and compare it with similar nebulae produced by other blue supergiants in the Galaxy
(Sher 25, HD 168625, SBW2007 1) and the Large Magellanic Cloud (Sk−69
$\deg$
202). The nitrogen abundances in these nebulae imply that blue supergiants can produce them from the main-sequence stage up to the pre-supernova stage. We also present a K-band spectrum of the candidate luminous blue variable MN56 (encircled by a ring-like nebula) and report the discovery of an OB star at ≈17 arcsec from MN18. The possible membership of MN18 and the OB star of the star cluster Lynga 3 is discussed.
The VLT-FLAMES Tarantula Survey Doran, E I; Crowther, P A; de Koter, A ...
Astronomy and astrophysics (Berlin),
10/2013, Letnik:
558
Journal Article
Recenzirano
The VLT-FLAMES Tarantula Survey has an extensive view of the copious number of massive stars in the 30 Doridus (30 Dor) star forming region of the Large Magellanic Cloud. The first comprehensive ...census of hot luminous stars in 30 Dor is compiled within a l0 arcmin (l50pc) radius of its central cluster, R136. We investigate the stellar content and spectroscopic completeness of the early type stars. Estimates were made for both the integrated ionising luminosity and stellar wind luminosity. These values were used to re-assess the star formation rate (SFR) of the region and determine the ionising photon escape fraction. When studying the most luminous star forming regions, it is essential to include their most massive stars if one is to determine a reliable energy budget. Photon leakage becomes more likely after including their large contributions to the ionising output. If 30 Dor is typical of other massive star forming regions, estimates of the SFR will be underpredicted if this escape fraction is not accounted for.
Context. The envelopes of stars near the Eddington limit are prone to various instabilities. A high Eddington factor in connection with the iron opacity peak leads to convective instability, and a ...corresponding envelope inflation may induce pulsational instability. Here, we investigate the occurrence and consequences of both instabilities in models of Wolf-Rayet stars. Aims. We determine the convective velocities in the sub-surface convective zones to estimate the amplitude of the turbulent velocity at the base of the wind that potentially leads to the formation of small-scale wind structures, as observed in several Wolf-Rayet stars. We also investigate the effect of stellar wind mass loss on the pulsations of our stellar models. Methods. We approximated solar metallicity Wolf-Rayet stars in the range 2−17 M⊙ by models of mass-losing helium stars, computed with the Bonn stellar evolution code. We characterized the properties of convection in the envelope of these stars adopting the standard mixing length theory. Results. Our results show the occurrence of sub-surface convective regions in all studied models. Small (≈1 km s-1) surface velocity amplitudes are predicted for models with masses below ≈10 M⊙. For models with M ≳ 10 M⊙, the surface velocity amplitudes are of the order of 10 km s-1. Moreover we find the occurrence of pulsations for stars in the mass range 9−14 M⊙, while mass loss appears to stabilize the more massive Wolf-Rayet stars. We confront our results with observationally derived line variabilities of 17 WN stars, of which we analysed eight here for the first time. The data suggest variability to occur for stars above 10 M⊙, which is increasing linearly with mass above this value, in agreement with our results. We further find our models in the mass range 9−14M⊙ to be unstable to radial pulsations, and predict local magnetic fields of the order of hundreds of gauss in Wolf-Rayet stars more massive than ≈10 M⊙. Conclusions. Our study relates the surface velocity fluctuations induced by sub-surface convection to the formation of clumping in the inner part of the wind. From this mechanism, we expect a stronger variability in more massive Wolf-Rayet stars, and a weaker variability in corresponding low metallicity Wolf-Rayet stars.