-
Connecting stellar and galactic scales: Energetic feedback from stellar wind bubbles to supernova remnants
Authors:
Yvonne A. Fichtner,
Jonathan Mackey,
Luca Grassitelli,
Emilio Romano-Díaz,
Cristiano Porciani
Abstract:
Energy and momentum feedback from stars is a key element of models for galaxy formation and interstellar medium dynamics, but resolving the relevant length scales to directly include this feedback remain out of reach of current-generation simulations. We aim to constrain the energy feedback by winds, photoionisation and supernovae (SNe) from massive stars. We measure the thermal and kinetic energy…
▽ More
Energy and momentum feedback from stars is a key element of models for galaxy formation and interstellar medium dynamics, but resolving the relevant length scales to directly include this feedback remain out of reach of current-generation simulations. We aim to constrain the energy feedback by winds, photoionisation and supernovae (SNe) from massive stars. We measure the thermal and kinetic energy imparted to the interstellar medium on various length scales, calculated from high-resolution 1D radiation-hydrodynamics simulations. Our grid of simulations covers a broad range of densities, metallicities, and state-of-the-art evolutionary models of single and binary stars. We find that a single star or binary system can carve a cavity of tens-of-pc size into the surrounding medium. During the pre-SN phase, post-main-sequence stellar winds and photoionisation dominate. While SN explosions dominate the total energy budget, the pre-SN feedback is of great importance by reducing the circumstellar gas density and delaying the onset of radiative losses in the SN remnant. Contrary to expectations, the metallicity dependence of the stellar wind has little effect on the cumulative energy imparted by feedback to the ISM; the only requirement is the existence of a sufficient level of pre-SN radiative and mechanical feedback. The ambient medium density determines how much and when feedback energy reaches to distance $\gtrsim 10-20$ pc and affects the division between kinetic and thermal feedback. Our results can be used as a sub-grid model for feedback in large-scale simulations of galaxies. The results reinforce that the uncertain mapping of stellar evolution sequences to SN explosion energy is very important to determining the overall feedback energy from a stellar population.
△ Less
Submitted 21 August, 2024; v1 submitted 16 February, 2024;
originally announced February 2024.
-
X-Shooting ULLYSES: massive stars at low metallicity. I. Project Description
Authors:
Jorick S. Vink,
A. Mehner,
P. A. Crowther,
A. Fullerton,
M. Garcia,
F. Martins,
N. Morrell,
L. M. Oskinova,
N. St-Louis,
A. ud-Doula,
A. A. C. Sander,
H. Sana,
J. -C. Bouret,
B. Kubatova,
P. Marchant,
L. P. Martins,
A. Wofford,
J. Th. van Loon,
O. Grace Telford,
Y. Gotberg,
D. M. Bowman,
C. Erba,
V. M. Kalari,
M. Abdul-Masih,
T. Alkousa
, et al. (56 additional authors not shown)
Abstract:
Observations of individual massive stars, super-luminous supernovae, gamma-ray bursts, and gravitational-wave events involving spectacular black-hole mergers, indicate that the low-metallicity Universe is fundamentally different from our own Galaxy. Many transient phenomena will remain enigmatic until we achieve a firm understanding of the physics and evolution of massive stars at low metallicity…
▽ More
Observations of individual massive stars, super-luminous supernovae, gamma-ray bursts, and gravitational-wave events involving spectacular black-hole mergers, indicate that the low-metallicity Universe is fundamentally different from our own Galaxy. Many transient phenomena will remain enigmatic until we achieve a firm understanding of the physics and evolution of massive stars at low metallicity (Z). The Hubble Space Telescope has devoted 500 orbits to observe 250 massive stars at low Z in the ultraviolet (UV) with the COS and STIS spectrographs under the ULLYSES program. The complementary ``X-Shooting ULLYSES'' (XShootU) project provides enhanced legacy value with high-quality optical and near-infrared spectra obtained with the wide-wavelength coverage X-shooter spectrograph at ESO's Very Large Telescope.
We present an overview of the XShootU project, showing that combining ULLYSES UV and XShootU optical spectra is critical for the uniform determination of stellar parameters such as effective temperature, surface gravity, luminosity, and abundances, as well as wind properties such as mass-loss rates in function of Z. As uncertainties in stellar and wind parameters percolate into many adjacent areas of Astrophysics, the data and modelling of the XShootU project is expected to be a game-changer for our physical understanding of massive stars at low Z.
To be able to confidently interpret James Webb Space Telescope (JWST) spectra of the first stellar generations, the individual spectra of low Z stars need to be understood, which is exactly where XShootU can deliver.
△ Less
Submitted 1 June, 2023; v1 submitted 10 May, 2023;
originally announced May 2023.
-
Inverse Compton cooling of thermal plasma in colliding-wind binaries
Authors:
Jonathan Mackey,
Thomas A. K. Jones,
Robert Brose,
Luca Grassitelli,
Brian Reville,
Arun Mathew
Abstract:
The inverse-Compton effect (IC) is a widely recognized cooling mechanism for both relativistic and thermal electrons in various astrophysical environments, including the intergalactic medium and X-ray emitting plasmas. Its effect on thermal electrons is however frequently overlooked in theoretical and numerical models of colliding-wind binaries (CWB). In this article, we provide a comprehensive in…
▽ More
The inverse-Compton effect (IC) is a widely recognized cooling mechanism for both relativistic and thermal electrons in various astrophysical environments, including the intergalactic medium and X-ray emitting plasmas. Its effect on thermal electrons is however frequently overlooked in theoretical and numerical models of colliding-wind binaries (CWB). In this article, we provide a comprehensive investigation of the impact of IC cooling in CWBs, presenting general results for when the photon fields of the stars dominate the cooling of the thermal plasma and when shocks at the stagnation point are expected to be radiative. Our analysis shows that IC cooling is the primary cooling process for the shocked-wind layer over a significant portion of the relevant parameter space, particularly in eccentric systems with large wind-momentum ratios, e.g., those containing a Wolf-Rayet and O-type star. Using the binary system WR 140 as a case study, we demonstrate that IC cooling leads to a strongly radiative shocked wind near periastron, which may otherwise remain adiabatic if only collisional cooling was considered. Our results are further supported by 2D and 3D simulations of wind-wind collisions. Specifically, 3D magnetohydrodynamic simulations of WR 140 show a significant decrease in hard-X-ray emission around periastron, in agreement with observations but in contrast to equivalent simulations that omit IC cooling. A novel method is proposed for constraining mass-loss rates of both stars in eccentric binaries where the wind-collision zone switches from adiabatic to radiative approaching periastron. IC scattering is an important cooling process in the thermal plasma of CWBs.
△ Less
Submitted 13 September, 2023; v1 submitted 31 January, 2023;
originally announced January 2023.
-
Bringing Stellar Evolution & Feedback Together: Summary of proposals from the Lorentz Center Workshop, 2022
Authors:
Sam Geen,
Poojan Agrawal,
Paul A. Crowther,
B. W. Keller,
Alex de Koter,
Zsolt Keszthelyi,
Freeke van de Voort,
Ahmad A. Ali,
Frank Backs,
Lars Bonne,
Vittoria Brugaletta,
Annelotte Derkink,
Sylvia Ekström,
Yvonne A. Fichtner,
Luca Grassitelli,
Ylva Götberg,
Erin R. Higgins,
Eva Laplace,
Kong You Liow,
Marta Lorenzo,
Anna F. McLeod,
Georges Meynet,
Megan Newsome,
G. André Oliva,
Varsha Ramachandran
, et al. (12 additional authors not shown)
Abstract:
Stars strongly impact their environment, and shape structures on all scales throughout the universe, in a process known as ``feedback''. Due to the complexity of both stellar evolution and the physics of larger astrophysical structures, there remain many unanswered questions about how feedback operates, and what we can learn about stars by studying their imprint on the wider universe. In this whit…
▽ More
Stars strongly impact their environment, and shape structures on all scales throughout the universe, in a process known as ``feedback''. Due to the complexity of both stellar evolution and the physics of larger astrophysical structures, there remain many unanswered questions about how feedback operates, and what we can learn about stars by studying their imprint on the wider universe. In this white paper, we summarize discussions from the Lorentz Center meeting `Bringing Stellar Evolution and Feedback Together' in April 2022, and identify key areas where further dialogue can bring about radical changes in how we view the relationship between stars and the universe they live in.
△ Less
Submitted 31 January, 2023;
originally announced January 2023.
-
Mechanical feedback from stellar winds with an application to galaxy formation at high redshift
Authors:
Yvonne A. Fichtner,
Luca Grassitelli,
Emilio Romano-Díaz,
Cristiano Porciani
Abstract:
We compute different sets of stellar evolutionary tracks in order to quantify the energy, mass, and metals yielded by massive main-sequence and post-main-sequence winds. Our aim is to investigate the impact of binary systems and of a metallicity-dependent distribution of initial rotational velocities on the feedback by stellar winds. We find significant changes compared to the commonly used non-ro…
▽ More
We compute different sets of stellar evolutionary tracks in order to quantify the energy, mass, and metals yielded by massive main-sequence and post-main-sequence winds. Our aim is to investigate the impact of binary systems and of a metallicity-dependent distribution of initial rotational velocities on the feedback by stellar winds. We find significant changes compared to the commonly used non-rotating, single-star scenario. The largest differences are noticeable at low metallicity, where the mechanical-energy budget is substantially increased. So as to establish the maximal (i.e. obtained by neglecting dissipation in the near circumstellar environment) influence of winds on the early stages of galaxy formation, we use our new feedback estimates to simulate the formation and evolution of a sub-$L_*$ galaxy at redshift 3 (hosted by a dark-matter halo with a mass of $1.8\times 10^{11}$ M$_\odot$) and compare the outcome with simulations in which only supernova (SN) feedback is considered. Accounting for the continuous energy injection by winds reduces the total stellar mass, the metal content, and the burstiness of the star-formation rate as well as of the outflowing gas mass. However, our numerical experiment suggests that the enhanced mechanical feedback from the winds of rotating and binary stars has a limited impact on the most relevant galactic properties compared to the non-rotating single-star scenario. Eventually, we look at the relative abundance between the metals entrained in winds and those ejected by SNe and find that it stays nearly constant within the simulated galaxy and its surrounding halo.
△ Less
Submitted 22 April, 2022; v1 submitted 18 January, 2022;
originally announced January 2022.
-
Modeling non-thermal emission from SN 1987A
Authors:
Robert Brose,
Jonathan Mackey,
Sean Kelly,
Nathan Grin,
Luca Grassitelli
Abstract:
The remnant of SN 1987A is the best-studied object of its kind. The rich data-set of its thermal and non-thermal emission across the electromagnetic spectrum poses a unique testbed for the elaboration of particle-acceleration theory.
We use 2D simulations of the progenitor's wind to obtain hydro-profiles for the medium around the supernova explosion. Various cones along prominent features of the…
▽ More
The remnant of SN 1987A is the best-studied object of its kind. The rich data-set of its thermal and non-thermal emission across the electromagnetic spectrum poses a unique testbed for the elaboration of particle-acceleration theory.
We use 2D simulations of the progenitor's wind to obtain hydro-profiles for the medium around the supernova explosion. Various cones along prominent features of the ambient medium are then used in our time-dependent acceleration code RATPaC to model the evolution of the emission of SN 1987A and compare it to observational data. We solve for the transport of cosmic rays and the hydrodynamical flow, in the test-particle limit.The simulation code relies on 1D profiles but the large expansion speed of the young remnant renders lateral transport unimportant.
We find that the increase in thermal X-ray emission predates the increase in the low-energy gamma-ray brightness by several years. The increase of the gamma-ray brightness at lower energies is followed by a smooth increase at the highest energies. The gamma-ray spectrum at the highest energies appears soft during the brightening but hardens as more material in the equatorial ring gets shocked. The X-ray and gamma-ray brightness remain almost constant once the SNR blast-wave passed the region of peak-density in the equatorial plane.
△ Less
Submitted 13 September, 2021;
originally announced September 2021.
-
On the Origin of Stochastic, Low-Frequency Photometric Variability in Massive Stars
Authors:
Matteo Cantiello,
Daniel Lecoanet,
Adam S. Jermyn,
Luca Grassitelli
Abstract:
High-precision photometric observations have revealed ubiquitous stochastic low-frequency photometric variability in early type stars. It has been suggested that this variability arises due to either subsurface convection or internal gravity waves launched by the convective core. Here we show that relevant properties of convection in subsurface convective layers correlate very well with the timesc…
▽ More
High-precision photometric observations have revealed ubiquitous stochastic low-frequency photometric variability in early type stars. It has been suggested that this variability arises due to either subsurface convection or internal gravity waves launched by the convective core. Here we show that relevant properties of convection in subsurface convective layers correlate very well with the timescale and amplitude of stochastic low-frequency photometric variability, as well as with the amplitude of macroturbulence. We suggest that low-frequency, stochastic photometric variability and surface turbulence in massive stars are caused by the the presence of subsurface convection. We show that an explanation for the observed surface photometric variability and macroturbulence relying on convective core driven internal gravity waves encounters a number of difficulties and seems unlikely to be able to explain the observed trends.
△ Less
Submitted 19 May, 2021; v1 submitted 10 February, 2021;
originally announced February 2021.
-
Wind-envelope interaction as the origin of the slow cyclic brightness variations of luminous blue variables
Authors:
Luca Grassitelli,
Norbert Langer,
Jonathan Mackey,
Goetz Graefener,
Nathan Grin,
Andreas Sander,
Jorick Vink
Abstract:
Luminous blue variables (LBVs) are hot, very luminous massive stars displaying large quasi-periodic variations in brightness, radius,and photospheric temperature, on timescales of years to decades. The physical origin of this variability, called S Doradus cycle after its prototype, has remained elusive. Here, we study the feedback of stellar wind mass-loss on the envelope structure in stars near t…
▽ More
Luminous blue variables (LBVs) are hot, very luminous massive stars displaying large quasi-periodic variations in brightness, radius,and photospheric temperature, on timescales of years to decades. The physical origin of this variability, called S Doradus cycle after its prototype, has remained elusive. Here, we study the feedback of stellar wind mass-loss on the envelope structure in stars near the Eddington limit. We perform a time-dependent hydrodynamic stellar evolutionary calculation, applying a stellar wind mass-loss prescription with a temperature-dependence inspired by the predicted systematic increase in mass-loss rates below 25 kK. We find that when the wind mass-loss rate crosses a well-defined threshold, a discontinuous change in the wind base conditions leads to a restructuring of the stellar envelope. The induced drastic radius and temperature changes, which occur on the thermal timescale of the inflated envelope, impose in turn mass-loss variations that reverse the initial changes, leading to a cycle that lacks a stationary equilibrium configuration. Our proof-of-concept model broadly reproduces the typical observational phenomenology of the S Doradus variability. We identify three key physical ingredients needed to trigger the instability: inflated envelopes in close proximity to the Eddington limit, a temperature range where decreasing opacities do not lead to an accelerating outflow, and a mass-loss rate that increases with decreasing temperature, crossing a critical threshold value within this temperature range. Our scenario and model provide testable predictions, and open the door for a consistent theoretical treatment of the LBV phase in stellar evolution, with consequences for their further evolution as single stars or in binary systems.
△ Less
Submitted 3 December, 2020; v1 submitted 30 November, 2020;
originally announced December 2020.
-
Turbulent Gas in Lensed Planck-selected Starbursts at redshifts 1-3.5
Authors:
Kevin C. Harrington,
Axel Weiss,
Min S. Yun,
Benjamin Magnelli,
C. E. Sharon,
T. K. D. Leung,
A. Vishwas,
Q. D. Wang,
E. F. Jimenez-Andrade,
D. T. Frayer,
D. Liu,
P. Garcia,
E. Romano-Diaz,
B. L. Frye,
S. Jarugula,
T. Badescu,
D. Berman,
H. Dannerbauer,
A. Diaz-Sanchez,
L. Grassitelli,
P. Kamieneski,
W. J. Kim,
A. Kirkpatrick,
J. D. Lowenthal,
H. Messias
, et al. (4 additional authors not shown)
Abstract:
Dusty star-forming galaxies at high redshift (1 < z < 3) represent the most intense star-forming regions in the Universe. Key aspects to these processes are the gas heating and cooling mechanisms. Although it is well known that these galaxies are gas-rich, little is known about the gas excitation conditions. Here we examine these processes in a sample of 24 strongly lensed star-forming galaxies id…
▽ More
Dusty star-forming galaxies at high redshift (1 < z < 3) represent the most intense star-forming regions in the Universe. Key aspects to these processes are the gas heating and cooling mechanisms. Although it is well known that these galaxies are gas-rich, little is known about the gas excitation conditions. Here we examine these processes in a sample of 24 strongly lensed star-forming galaxies identified by the \textit{Planck} satellite (LPs) at z ~ 1.1 - 3.5. We analyze 162 CO rotational transitions (ranging from Jupper = 1 - 12) and 37 atomic carbon fine-structure lines ([CI]) in order to characterize the physical conditions of the gas in sample of LPs. We simultaneously fit the CO and [CI] lines, and the dust continuum emission, using two different non-LTE, radiative transfer models. The first model represents a two component gas density, while the second assumes a turbulence driven log-normal gas density distribution. These LPs are among the most gas-rich, infrared (IR) luminous galaxies ever observed ($μ_{\rm L}$L$_{\rm IR(8-1000μm) } \sim 10^{13-14.6} $\Lsun; $< μ_{\rm L}$M$_{\rm ISM}> = 2.7 \pm 1.2 \times 10^{12}$ \Msun, with $μ_{\rm L} \sim 10-30$ the average lens magnification factor). Our results suggest that the turbulent ISM present in the LPs can be well-characterized by a high turbulent velocity dispersion ($<ΔV_{\rm turb}> \sim 100 $ \kms) and gas kinetic temperature to dust temperature ratios $<T_{\rm kin}$/$T_{\rm d}> \sim 2.5$, sustained on scales larger than a few kpc. We speculate that the average surface density of the molecular gas mass and IR luminosity $Σ_{\rm M_{\rm ISM}}$ $\sim 10^{3 - 4}$ \Msun pc$^{-2}$ and $Σ_{\rm L_{\rm IR}}$ $\sim 10^{11 - 12}$ \Lsun kpc$^{-2}$, arise from both stellar mechanical feedback and a steady momentum injection from the accretion of intergalactic gas.
△ Less
Submitted 30 October, 2020;
originally announced October 2020.
-
Properties of OB star-black hole systems derived from detailed binary evolution models
Authors:
N. Langer,
C. Schürmann,
K. Stoll,
P. Marchant,
D. J. Lennon,
L. Mahy,
S. E. de Mink,
M. Quast,
W. Riedel,
H. Sana,
P. Schneider,
A. Schootemeijer,
Chen Wang,
L. A. Almeida,
J. M. Bestenlehner,
J. Bodensteiner,
N. Castro,
S. Clark,
P. A. Crowther,
P. Dufton,
C. J. Evans,
L. Fossati,
G. Gräfener,
L. Grassitelli,
N. Grin
, et al. (16 additional authors not shown)
Abstract:
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…
▽ More
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.
△ Less
Submitted 10 April, 2020; v1 submitted 20 December, 2019;
originally announced December 2019.
-
The "Red Radio Ring": Ionised and Molecular Gas in a Starburst/Active Galactic Nucleus at $z \sim 2.55$
Authors:
K. C. Harrington,
A. Vishwas,
A. Weiss,
B. Magnelli,
L. Grassitelli,
M. Zajacek,
E. F. Jimenez-Andrade,
T. K. D. Leung,
F. Bertoldi,
E. Romano-Diaz,
D. T. Frayer,
P. Kamieneski,
D. Riechers,
G. J. Stacey,
M. S. Yun,
Q. D. Wang
Abstract:
We report the detection of the far-infrared (FIR) fine-structure line of singly ionised nitrogen, \Nplusa, within the peak epoch of galaxy assembly, from a strongly lensed galaxy, hereafter ``The Red Radio Ring''; the RRR, at z = 2.55. We combine new observations of the ground-state and mid-J transitions of CO (J$_{\rm up} =$ 1,5,8), and the FIR spectral energy distribution (SED), to explore the m…
▽ More
We report the detection of the far-infrared (FIR) fine-structure line of singly ionised nitrogen, \Nplusa, within the peak epoch of galaxy assembly, from a strongly lensed galaxy, hereafter ``The Red Radio Ring''; the RRR, at z = 2.55. We combine new observations of the ground-state and mid-J transitions of CO (J$_{\rm up} =$ 1,5,8), and the FIR spectral energy distribution (SED), to explore the multi-phase interstellar medium (ISM) properties of the RRR. All line profiles suggest that the HII regions, traced by \Nplusa, and the (diffuse and dense) molecular gas, traced by the CO, are co-spatial when averaged over kpc-sized regions. Using its mid-IR-to-millimetre (mm) SED, we derive a non-negligible dust attenuation of the \Nplusa line emission. Assuming a uniform dust screen approximation results a mean molecular gas column density $> 10^{24}$\, cm$^{-2}$, with a molecular gas-to-dust mass ratio of 100. It is clear that dust attenuation corrections should be accounted for when studying FIR fine-structure lines in such systems. The attenuation corrected ratio of $L_{\rm NII205} / L_{\rm IR(8-1000μm)} = 2.7 \times 10^{-4}$ is consistent with the dispersion of local and $z >$ 4 SFGs. We find that the lower-limit, \Nplusa -based star-formation rate (SFR) is less than the IR-derived SFR by a factor of four. Finally, the dust SED, CO line SED and $L_{\rm NII205}$ line-to-IR luminosity ratio of the RRR is consistent with a starburst-powered ISM.
△ Less
Submitted 23 June, 2019;
originally announced June 2019.
-
Subsonic structure and optically thick winds from Wolf--Rayet stars
Authors:
Luca Grassitelli,
Norbert Langer,
Nathan J. Grin,
Jonathan Mackey,
Joachim M. Bestenlehner,
Goetz Graefener
Abstract:
Wolf-Rayet star's 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. We investigate the conditions and the structure of the subsonic part of the outflow of Galactic WR stars, in particular of the WNE su…
▽ More
Wolf-Rayet star's 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. We investigate the conditions and the structure of the subsonic part of the outflow of Galactic WR stars, in particular of the WNE subclass; our focus is on the conditions at the sonic point. 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 200kK by the Fe opacity bump or of the order of 50kK by the HeII opacity bump. For a given mass-loss rate, the conditions in the subsonic part of the outflow are independent from the detailed physical conditions in the supersonic part. 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 L/M ratio and the stellar mass-loss rate, thereby constraining the sonic point conditions, the subsonic structure, and the stellar wind mass-loss rates from observations. The minimum mass-loss rate necessary to have the flow accelerated to supersonic velocities by the Fe opacity bump is derived. A comparison of the observed parameters of Galactic WNE stars to this minimum mass-loss rate indicates that their winds are launched to supersonic velocities by the radiation pressure arising from the Fe-bump. Conversely, models which do not show transonic flows from the Fe opacity bump form inflated envelopes. We derive an analytic criterion for the appearance of envelope inflation in the subphotospheric layers.
△ Less
Submitted 8 March, 2018;
originally announced March 2018.
-
On the optically-thick winds of Wolf-Rayet stars
Authors:
G. Gräfener,
S. P. Owocki,
L. Grassitelli,
N. Langer
Abstract:
(abridged) The strong winds of Wolf-Rayet (WR) stars are important for the mechanical and chemical feedback of the most massive stars and determine whether they end their lives as neutron stars or black holes. In this work we investigate theoretically the mass-loss properties of H-free WR stars of the nitrogen sequence (WN stars). We connect stellar structure models for He stars with wind models f…
▽ More
(abridged) The strong winds of Wolf-Rayet (WR) stars are important for the mechanical and chemical feedback of the most massive stars and determine whether they end their lives as neutron stars or black holes. In this work we investigate theoretically the mass-loss properties of H-free WR stars of the nitrogen sequence (WN stars). We connect stellar structure models for He stars with wind models for optically-thick winds and assess how both types of models can simultaneously fulfill their respective sonic-point conditions. Fixing the outer wind law and terminal wind velocity, we obtain unique solutions for the mass-loss rates of optically-thick, radiatively-driven winds of WR stars in the phase of core He-burning. The resulting mass-loss relations as a function of stellar parameters, agree well with previous empirical relations. Furthermore, we encounter stellar mass limits below which no continuous solutions exist. While these mass limits agree with observations of WR stars in the Galaxy, they are in conflict with observations in the LMC. While our results confirm in particular the slope of oft-used empirical mass-loss relations, they imply that only part of the observed WN population can be understood in the framework of the standard assumptions of a smooth transonic flow and compact stellar core. This means that alternative approaches, such as a clumped and inflated wind structure, or deviations from the diffusion limit at the sonic point may have to be invoked. Qualitatively, the existence of mass limits for the formation of WR-type winds may be relevant for the non-detection of low-mass WR stars in binary systems, which are believed to be progenitors of Type Ib/c supernovae. The sonic-point conditions derived in this work may provide a possibility to include optically-thick winds in stellar evolution models in a more physically motivated form than in current models.
△ Less
Submitted 12 October, 2017;
originally announced October 2017.
-
The metallicity dependence of envelope inflation in massive stars
Authors:
D. Sanyal,
N. Langer,
D. Szécsi,
S. -C Yoon,
L. Grassitelli
Abstract:
Recently it has been found that models of massive stars reach the Eddington limit in their interior, which leads to dilute extended envelopes. We perform a comparative study of the envelope properties of massive stars at different metallicities, with the aim to establish the impact of the stellar metallicity on the effect of envelope inflation. We analyse published grids of core-hydrogen burning m…
▽ More
Recently it has been found that models of massive stars reach the Eddington limit in their interior, which leads to dilute extended envelopes. We perform a comparative study of the envelope properties of massive stars at different metallicities, with the aim to establish the impact of the stellar metallicity on the effect of envelope inflation. We analyse published grids of core-hydrogen burning massive star models computed with metallicities appropriate for massive stars in the Milky Way, the LMC and the SMC, the very metal poor dwarf galaxy I Zwicky 18, and for metal-free chemical composition. Stellar models of all the investigated metallicities reach and exceed the Eddington limit in their interior, aided by the opacity peaks of iron, helium and hydrogen, and consequently develop inflated envelopes. Envelope inflation leads to a redward bending of the zero-age main sequence and a broadening of the main sequence band in the upper part of the Hertzsprung-Russell diagram. We derive the limiting L/M-values as function of the stellar surface temperature above which inflation occurs, and find them to be larger for lower metallicity. While Galactic models show inflation above ~29 Msun, the corresponding mass limit for Population III stars is ~150 Msun. While the masses of the inflated envelopes are generally small, we find that they can reach 1-100 Msun in models with effective temperatures below ~8000 K, with higher masses reached by models of lower metallicity. Envelope inflation is expected to occur in sufficiently massive stars at all metallicities, and is expected to lead to rapidly growing pulsations, high macroturbulent velocities, and might well be related to the unexplained variability observed in Luminous Blue Variables like S Doradus and Eta Carina.
△ Less
Submitted 22 November, 2016;
originally announced November 2016.
-
Asteroseismic test of rotational mixing in low-mass white dwarfs
Authors:
A. G. Istrate,
G. Fontaine,
A. Gianninas,
L. Grassitelli,
P. Marchant,
T. M. Tauris,
N. Langer
Abstract:
We exploit the recent discovery of pulsations in mixed-atmosphere (He/H), extremely low-mass white dwarf precursors (ELM proto-WDs) to test the proposition that rotational mixing is a fundamental process in the formation and evolution of low-mass helium core white dwarfs. Rotational mixing has been shown to be a mechanism able to compete efficiently against gravitational settling, thus accounting…
▽ More
We exploit the recent discovery of pulsations in mixed-atmosphere (He/H), extremely low-mass white dwarf precursors (ELM proto-WDs) to test the proposition that rotational mixing is a fundamental process in the formation and evolution of low-mass helium core white dwarfs. Rotational mixing has been shown to be a mechanism able to compete efficiently against gravitational settling, thus accounting naturally for the presence of He, as well as traces of metals such as Mg and Ca, typically found in the atmospheres of ELM proto-WDs. Here we investigate whether rotational mixing can maintain a sufficient amount of He in the deeper driving region of the star, such that it can fuel, through HeII-HeIII ionization, the observed pulsations in this type of stars. Using state-of-the-art evolutionary models computed with MESA, we show that rotational mixing can indeed explain qualitatively the very existence and general properties of the known pulsating, mixed-atmosphere ELM proto-WDs. Moreover, such objects are very likely to pulsate again during their final WD cooling phase.
△ Less
Submitted 26 October, 2016;
originally announced October 2016.
-
The IACOB project: III. New observational clues to understand macroturbulent broadening in massive O- and B-type stars
Authors:
S. Simón-Díaz,
M. Godart,
N. Castro,
A. Herrero,
C. Aerts,
J. Puls,
J. Telting,
L. Grassitelli
Abstract:
We aim to provide new empirical clues about macroturbulent spectral line broadening in O- and B-type stars to evaluate its physical origin.
We use high-resolution spectra of ~430 stars with spectral types in the range O4-B9 (all luminosity classes). We characterize the line-broadening of adequate diagnostic metal lines using a combined FT and GOF technique. We perform a quantitative spectroscopi…
▽ More
We aim to provide new empirical clues about macroturbulent spectral line broadening in O- and B-type stars to evaluate its physical origin.
We use high-resolution spectra of ~430 stars with spectral types in the range O4-B9 (all luminosity classes). We characterize the line-broadening of adequate diagnostic metal lines using a combined FT and GOF technique. We perform a quantitative spectroscopic analysis of the whole sample using automatic tools coupled with a huge grid of FASTWIND models. We also incorporate quantitative information about line asymmetries to our observational description of the characteristics of the line-profiles, and present a comparison of the shape and type of line-profile variability found in a small sample of O stars and B supergiants with still undefined pulsational properties and B main sequence stars with variable line-profiles.
We present a homogeneous and statistically significant overview of the (single snapshot) line-broadening properties of stars in the whole O and B star domain. We find empirical evidence of the existence of various types of non-rotational broadening agents acting in the realm of massive stars. Even though all of them could be quoted and quantified as a macroturbulent broadening from a practical point of view, their physical origin can be different. Contrarily to the early- to late-B dwarfs/giants, which present a mixture of cases in terms of line-profile shape and variability, the whole O-type and B supergiant domain (or, roughly speaking, stars with M_ZAMS > 15 M_sol) is fully dominated by stars with a remarkable non-rotational broadening component and very similar profiles (including type of variability). We provide some examples illustrating how this observational dataset can be used to evaluate scenarios aimed at explaining the existence of sources of non-rotational broadening in massive stars.
△ Less
Submitted 19 August, 2016;
originally announced August 2016.
-
Metallicity dependence of turbulent pressure and macroturbulence in stellar envelopes
Authors:
Luca Grassitelli,
Luca Fossati,
Norbert Langer,
Sergio Simon-Diaz,
Norberto Castro,
Debashis Sanyal
Abstract:
Macroturbulence, introduced as a fudge to reproduce the width and shape of stellar absorption lines, reflects gas motions in stellar atmospheres. While in cool stars, it is thought to be caused by convection zones immediately beneath the stellar surface, the origin of macroturbulence in hot stars is still under discussion. Recent works established a correlation between the turbulent-to-total press…
▽ More
Macroturbulence, introduced as a fudge to reproduce the width and shape of stellar absorption lines, reflects gas motions in stellar atmospheres. While in cool stars, it is thought to be caused by convection zones immediately beneath the stellar surface, the origin of macroturbulence in hot stars is still under discussion. Recent works established a correlation between the turbulent-to-total pressure ratio inside the envelope of stellar models and the macroturbulent velocities observed in corresponding Galactic stars. To probe this connection further, we evaluated the turbulent pressure that arises in the envelope convective zones of stellar models in the mass range 1-125 Msun based on the mixing-length theory and computed for metallicities of the Large and Small Magellanic Cloud. We find that the turbulent pressure contributions in models with these metallicities located in the hot high-luminosity part of the Hertzsprung-Russel (HR) diagram is lower than in similar models with solar metallicity, whereas the turbulent pressure in low-metallicity models populating the cool part of the HR-diagram is not reduced. Based on our models, we find that the currently available observations of hot massive stars in the Magellanic Clouds appear to support a connection between macroturbulence and the turbulent pressure in stellar envelopes. Multidimensional simulations of sub-surface convection zones and a larger number of high-quality observations are necessary to test this idea more rigorously.
△ Less
Submitted 2 August, 2016;
originally announced August 2016.
-
Models of low-mass helium white dwarfs including gravitational settling, thermal and chemical diffusion, and rotational mixing
Authors:
Alina Istrate,
Pablo Marchant,
Thomas M. Tauris,
Norbert Langer,
Richard J. Stancliffe,
Luca Grassitelli
Abstract:
A large number of extremely low-mass helium white dwarfs (ELM WDs) have been discovered in recent years. The majority of them are found in close binary systems suggesting they are formed either through a common-envelope phase or via stable mass transfer in a low-mass X-ray binary (LMXB) or a cataclysmic variable (CV) system. Here, we investigate the formation of these objects through the LMXB chan…
▽ More
A large number of extremely low-mass helium white dwarfs (ELM WDs) have been discovered in recent years. The majority of them are found in close binary systems suggesting they are formed either through a common-envelope phase or via stable mass transfer in a low-mass X-ray binary (LMXB) or a cataclysmic variable (CV) system. Here, we investigate the formation of these objects through the LMXB channel with emphasis on the proto-WD evolution in environments with different metallicities. We study, for the first time, the combined effects of rotational mixing and element diffusion (e.g. gravitational settling, thermal and chemical diffusion) on the evolution of proto-WDs and on the cooling properties of the resulting WDs. We present state-of-the-art binary stellar evolution models computed with MESA for metallicities between Z=0.0002 and Z=0.02, producing WDs with masses between 0.16-0.45 M$_{\odot}$. Our results confirm that element diffusion plays a significant role in the evolution of proto-WDs that experience hydrogen shell flashes. The occurrence of these flashes produces a clear dichotomy in the cooling timescales of ELM WDs, which has important consequences e.g. for the age determination of binary millisecond pulsars. Rotational mixing is found to counteract the effect of gravitational settling in the surface layers of young, bloated ELM proto-WDs and therefore plays a key role in determining their surface chemical abundances. We predict that these proto-WDs have helium-rich envelopes through a significant part of their lifetime, a crucial ingredient for understanding the newly observed ELM proto-WD pulsators. The hydrogen envelope at detachment, although small compared to the total mass of the WD, contains enough angular momentum such that the spin frequency of the resulting WD on the cooling track is well above the orbital frequency.
△ Less
Submitted 15 June, 2016;
originally announced June 2016.
-
Diagnostic of the unstable envelopes of Wolf-Rayet stars
Authors:
Luca Grassitelli,
Andre-Nicolas Chene',
Debashis Sanyal,
Norbert Langer,
Nicole St. Louis,
Joachim Bestenlehner,
Luca Fossati
Abstract:
The envelopes of stars near the Eddington limit are prone to various instabilities. A high Eddington factor in connection with the Fe 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. We determine the convective velocities…
▽ More
The envelopes of stars near the Eddington limit are prone to various instabilities. A high Eddington factor in connection with the Fe 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. 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 WR stars. We also investigate the effect of mass loss on the pulsations of our models. We approximated solar metallicity WR stars by models of mass-losing helium stars, and we characterized the properties of convection in the envelope adopting the standard MLT. Our results show the occurrence of sub-surface convective regions in all studied models. Small surface velocity amplitudes are predicted for models with masses below 10Msun. For models with M>10Msun, the surface velocity amplitudes are of the order of 10km/s. Moreover we find the occurrence of pulsations for stars in the mass range 9-14Msun, while mass loss appears to stabilize the more massive WR stars. We confront our results with observationally derived line variabilities of 17 WN stars. The data suggest variability to occur for stars above 10Msun, which is increasing linearly with mass above this value, in agreement with our results. We further find some of our models to be unstable to radial pulsations, and predict local magnetic fields of the order of hundreds of Gauss in WR stars more massive than 10Msun. 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 WR stars.
△ Less
Submitted 29 March, 2016;
originally announced March 2016.
-
Relating turbulent pressure and macroturbulence across the HR diagram with a possible link to γ-Dor stars
Authors:
Luca Grassitelli,
Luca Fossati,
Norbert Langer,
Andrea Miglio,
Alina Georgiana Istrate,
Debashis Sanyal
Abstract:
A significant fraction of the envelope of low- and intermediate-mass stars is unstable to convection, leading to sub-surface turbulent motion. Here, we consider and include the effects of turbulence pressure in our stellar evolution calculations. In search of an observational signature, we compare the fractional contribution of turbulent pressure to the observed macroturbulent velocities in stars…
▽ More
A significant fraction of the envelope of low- and intermediate-mass stars is unstable to convection, leading to sub-surface turbulent motion. Here, we consider and include the effects of turbulence pressure in our stellar evolution calculations. In search of an observational signature, we compare the fractional contribution of turbulent pressure to the observed macroturbulent velocities in stars at different evolutionary stages. We find a strong correlation between the two quantities, similar to what was previously found for massive OB stars. We therefore argue that turbulent pressure fluctuations of finite amplitude may excite high-order, high-angular degree stellar oscillations, which manifest themselves at the surface an additional broadening of the spectral lines, i.e., macroturbulence, across most of the HR diagram. When considering the locations in the HR diagram where we expect high-order oscillations to be excited by stochastic turbulent pressure fluctuations, we find a close match with the observational $γ$ Doradus instability strip, which indeed contains high-order, non-radial pulsators. We suggest that turbulent pressure fluctuations on a percentual level may contribute to the $γ$ Dor phenomenon, calling for more detailed theoretical modelling in this direction.
△ Less
Submitted 4 November, 2015;
originally announced November 2015.
-
Observational consequences of turbulent pressure in the envelopes of massive stars
Authors:
Luca Grassitelli,
Luca Fossati,
Sergio Simon-Diaz,
Norbert Langer,
Norberto Castro,
Debashis Sanyal
Abstract:
The major mass fraction of the envelope of hot luminous stars is radiatively stable. However, the partial ionisation of hydrogen, helium and iron gives rise to extended sub-surface convection zones in all of them. In this work, we investigate the effect of the pressure induced by the turbulent motion in these zones based on the mixing length theory, and search for observable consequences. We find…
▽ More
The major mass fraction of the envelope of hot luminous stars is radiatively stable. However, the partial ionisation of hydrogen, helium and iron gives rise to extended sub-surface convection zones in all of them. In this work, we investigate the effect of the pressure induced by the turbulent motion in these zones based on the mixing length theory, and search for observable consequences. We find that the turbulent pressure fraction can amount up to ~5% in OB supergiants, and to ~30% in cooler supergiants. The resulting structural changes are, however, not significantly affecting the evolutionary tracks compared to previous calculations. Instead, a comparison of macroturbulent velocities derived from high quality spectra of OB stars with the turbulent pressure fraction obtained in corresponding stellar models reveals a strong correlation of these two quantities. We discuss a possible physical connection, and conclude that turbulent pressure fluctuations may drive high-order oscillations, which - as conjectured earlier - manifest themselves as macroturbulence in the photospheres of hot luminous stars.
△ Less
Submitted 14 July, 2015;
originally announced July 2015.
-
Massive main sequence stars evolving at the Eddington limit
Authors:
Debashis Sanyal,
Luca Grassitelli,
Norbert Langer,
Joachim M. Bestenlehner
Abstract:
The evolution of massive stars even on the main sequence is not yet well understood. Due to the steep mass-luminosity relation, massive main sequence stars become very luminous. This brings their envelopes very close to the Eddington limit. We are analysing stellar evolutionary models in which the Eddington limit is reached and exceeded, and explore the rich diversity of physical phenomena which t…
▽ More
The evolution of massive stars even on the main sequence is not yet well understood. Due to the steep mass-luminosity relation, massive main sequence stars become very luminous. This brings their envelopes very close to the Eddington limit. We are analysing stellar evolutionary models in which the Eddington limit is reached and exceeded, and explore the rich diversity of physical phenomena which take place in their envelopes, and investigate their observational consequences. We use the grids of detailed stellar models by Brott et al. (2011) and Koehler et al. (2015), to investigate the envelope properties of core hydrogen burning massive stars. We find that at the stellar surface, the Eddington limit is almost never reached, even for stars up to 500 Msun. When an appropriate Eddington limit is defined locally in the stellar envelope, most stars more massive than 40 Msun actually exceed this limit, in particular in the partial ionization 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*10^5 Lsun i.e. stars above 40 Msun 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 to produce major LBV eruptions.
△ Less
Submitted 9 June, 2015;
originally announced June 2015.