Towards a Unified View of Inhomogeneous Stellar Winds in Isolated Supergiant Stars and Supergiant High Mass X-Ray Binaries

Massive stars, at least \(\sim10\) times more massive than the Sun, have two key properties that make them the main drivers of evolution of star clusters, galaxies, and the Universe as a whole. On the one hand, the outer layers of massive stars are so hot that they produce most of the ionizing ultraviolet radiation of galaxies; in fact, the first massive stars helped to re-ionize the Universe after its Dark Ages. Another important property of massive stars are the strong stellar winds and outflows they produce. This mass loss, and finally the explosion of a massive star as a supernova or a gamma-ray burst, provide a significant input of mechanical and radiative energy into the interstellar space. These two properties together make massive stars one of the most important cosmic engines: they trigger the star formation and enrich the interstellar medium with heavy elements, that ultimately leads to formation of Earth-like rocky planets and the development of complex life. The study of massive star winds is thus a truly multidisciplinary field and has a wide impact on different areas of astronomy.In recent years observational and theoretical evidences have been growing that these winds are not smooth and homogeneous as previously assumed, but rather populated by dense “clumps”. The presence of these structures dramatically affects the mass loss rates derived from the study of stellar winds. Clump properties in isolated stars are nowadays inferred mostly through indirect methods (i.e., spectroscopic observations of line profiles in various wavelength regimes, and their analysis based on tailored, inhomogeneous wind models). The limited characterization of the clump physical properties (mass, size) obtained so far have led to large uncertainties in the mass loss rates from massive stars. Such uncertainties limit our understanding of the role of massive star winds in galactic and cosmic evolution.Supergiant high mass X-ray binaries (SgXBs) are among the brightest X-ray sources in the sky. A large number of them consist of a neutron star accreting from the wind of a massive companion and producing a powerful X-ray source. The characteristics of the stellar wind together with the complex interactions between the compact object and the donor star determine the observed X-ray output from all these systems. Consequently, the use of SgXBs for studies of massive stars is only possible when the physics of the stellar winds, the compact objects, and accretion mechanisms are combined together and confronted with observations.This detailed review summarises the current knowledge on the theory and observations of winds from massive stars, as well as on observations and accretion processes in wind-fed high mass X-ray binaries. The aim is to combine in the near future all available theoretical diagnostics and observational measurements to achieve a unified picture of massive star winds in isolated objects and in binary systems.     

Presupernova evolution of the most massive stars

Taking as example a 60M _⊙ star of solar metallicity, the state of the art of model calculations for very massive, from the main sequence to the supernova stage, is reviewed. It is argued that — due to the simple internal structure of Wolf-Rayet stars — the post main sequence evolutionary phases are currently those which are better understood. A brief discussion of the supernova outcome from very massive stars is given. Then, the more uncertain main sequence evolution is discussed. A first attempt to incorporate results about pulsational instabilities of very massive stars in stellar evolutionary calculations is performed. On its basis, a new type of evolutionary sequence for very massive stars is obtained, namely O-star → Of-star → H-rich WNL → LBV → H-poor WNL → WNE → WC → SN. This scenario is shown to correspond better to many observed properties of very massive stars than the standard one. It includes a model for the prototype LBV P Cygni.     

Massive stars: Setting the stage

The paper gives a summary of the situation mid-1993 of theory and observations regarding massive stars. I describe: stellar mass loss and its implications, pre-main-sequence evolution, the main sequence, problems of atmospheric instability, Luminous Blue Supergiants, Yellow Hypergiants, Wolf-Rayet stars and supernovae.     

The synthesis of helium and CNO isotopes in massive stars

We present helium and CNO isotopic yields for massive mass-losing stars in the initial mass range 15M M _i 50M . We investigate their dependence on assumptions about mass loss rates, internal mixing processes, and metallicity, and specify the contributions from stellar winds and from supernova ejecta.     

The population of massive stars in R136 from HST/FOC UV observations

New ultraviolet (1300 A, 3400 A),HST FOC observations have been used to derive the UV color-magnitude diagram (CMD) of R136, with the main scientific goal of studying the upper end of the stellar mass function at ultraviolet wavelengths where the color degeneracy encountered in visual CMDs is less severe. The CMD has been compared to a set of theoretical isochrones, which have been computed using the latest generation of evolutionary models and model atmospheres for early type stars. Wolf-Rayet stars are included. Comparison of theTheoretical andobserved CMD suggests that there are no stars brighter than M₁₃₀–11. We use the observed main sequence turn-off and the known spectroscopic properties of the stellar population to derive constraints on the most probable age of R136. The presence of WNL stars and the lack of red supergiants suggests a most likely age of 3±1 Myr. A theoretical isochrone of 3±1 Myr is consistent with the observed stellar content of R136 if the most massive stars have initial masses around 50 M.Bases on Observations with the NASA/ESA Hubble Space Telescope, obtained at the STScI, which is operated by AURA, Inc., under NASA contract NAS5-26555.Astrophysics Division, Space Science Department, ESA     

Formation and evolution of low-mass Algols

We discuss the origin, evolution and fate of low-mass Algols (LMA) that have components with initial masses less than 2.5 M₀. The semi-major axes of orbits of pre-LMA do not exceed 20–25 R₀. The rate of formation of Algol-type stars is ～ 0.01/year. Magnetic stellar winds may be the factor that determines the evolution of LMA. Most LMA end their lives as double helium degenerate dwarfs with M₁/M₂ ～ 0.88 (like L870-2). Some of them even merge through angular momentum loss caused by gravitational waves.     

Massive star distribution in external galaxies and starburst regions

Counts of hot and luminous stars in a number of associations in the Galaxy and Magellanic Clouds enable one to directly investigate the numbers and types of massive stars. There seems to be little, if any, dependence of the slope of the Intial Mass Function, or theM _upper on the initial composition of the stars. Indirect estimates of numbers of massive stars in various more distant environments are reviewed and discussed within a framework of acalibration of the methods using the stellar census of 30 Doradus. Very young starbursts, containing large numbers of massive stars, seem to be composed of smaller sub-units similar or somewhat larger than that object. These units might be newly born globular clusters.     

Formation and evolution of low-mass Algols

We discuss the origin, evolution and fate of low-mass Algols (LMA) that have components with initial masses less than 2.5 M₀. The semi-major axes of orbits of pre-LMA do not exceed 20–25 R₀. The rate of formation of Algol-type stars is 0.01/year. Magnetic stellar winds may be the factor that determines the evolution of LMA. Most LMA end their lives as double helium degenerate dwarfs with M₁/M₂ 0.88 (like L870-2). Some of them even merge through angular momentum loss caused by gravitational waves.     

Massive Stars: Input Physics and Stellar Models

We present a general overview of the structure and evolution of massive stars of masses ≥12 M _⊙ during their pre-supernova stages. We think it is worth reviewing this topic owing to the crucial role of massive stars in astrophysics, especially in the evolution of galaxies and the universe. We have performed several test computations with the aim to analyze and discuss many physical uncertainties still encountered in massive-star evolution. In particular, we explore the effects of mass loss, convection, rotation, ¹²C(α,γ)¹⁶O reaction and initial metallicity. We also compare and analyze the similarities and differences among various works and ours. Finally, we present useful comments on the nucleosynthesis from massive stars concerning the s-process and the yields for ²⁶Al and ⁶⁰Fe.     

Observations of the atmospheres and winds of O-stars,LBVs and Wolf-Rayet stars

This review summarises recent studies of O-stars, Luminous Blue Variables (LBVs) and Wolf-Rayet (WR) stars, emphasising observations and analyses of their atmospheres and stellar winds yielding determinations of their physical and chemical properties. Studies of these stellar groups provide important tests of both stellar wind theory and stellar evolution models incorporating mass-loss effects. Quantitative analyses of O-star spectra reveal enhanced helium abundances in Of and many luminous O-supergiants, together with CNO anomalies in OBN and Ofpe/WN9 stars, indicative of evolved objects. Enhanced helium, and CNO-cycle products are observed in several LBVs, implying a highly evolved status, whilst for the WR stars there is strong evidence for the exposition of CNO-cycle products in WN stars, and helium-burning products in WC and WO stars. The observed wind properties and mass-loss rates derived for O-stars show, in general terms, good agreement with predictions from the latest radiation-driven wind models, although some discrepancies are apparent. Several LBVs show similar mass-loss rates at maximum and minimum states, contrary to previous expectations, with the mechanism responsible for the variability and outbursts remaining unclear. WR stars exhibit the most extreme levels of mass-loss and stellar wind momenta. Whilst alternative mass-loss mechanisms have been proposed, recent calculations indicate that radiation pressure alone may be sufficient, given the strong ionization stratification present in their winds.     

The gradual acceleration of the outflow of VX-Sagittarii: A stellar evolutionary effect?

VX-Sagittarii is a red supergiant with a superwind which is observed in several maser lines. They provide an evidence that the outflow velocity keeps growing considerably at large distance from the star. It is argued that this phenomenon can be explained by stellar evolutionary effects.As a rule, the outflow velocity for late type stars correlates with the mass loss rate and from that it is suggested that the mass loss rate was higher in the past and is decreasing now. The mass of VX Sagittarii can be estimated on this basis and is about 40–50M      

The evolution of massive stars to explosion

We review the possible evolutionary paths from massive stars to explosive endpoints as various types of supernovae associated with Population I and hence with massive stars: Type II-P, Type II-L, Type Ib, Type Ic, and the hybrid events SN 1987K and SN 1993J. We identify SN 1954A as another hybrid event from the evidence for both H and He in its spectrum with velocities nearly the same as SN 1983J. Evidence for ejected⁵⁶Ni mass of 0.07 M suggests that SN II-P underwent standard iron core collapse, not collapse of an O–Ne–Mg core nor thermonuclear explosion of a C–O core. Most SN II-P presumably arise in single stars or wide binaries of 10–20 M. There may be indirect evidence for duplicity in some cases in the form of strong Ba II lines, such as characterized SN 1987A. SN II-L are recognizably distinct from typical SN II-P and must undergo a significantly different evolution. Despite indications that SN II-L have small envelopes that may be helium enriched, they are also distinct from events like SN 1993J that must have yet again a different evolution. The SN II-L that share a common Luminosity seem to have ejected a small nickel mass and hence may come from stars with O–Ne–Mg cores. The amount of nickel ejected by the exceptionally bright events, SN 1980K and SN 1979C, remains controversial. SN Ib require the complete loss of the H envelope, either to a binary companion or to a wind. The few identified have relatively large ejecta masses. It is not clear what evolutionary processes distinguish SN Ib's evolving in binary systems from hybrid events that retain some H in the envelope. SN Ic events are both H and He deficient. Binary models that can account for transfer of an extended helium envelope from low mass helium cores, 2 to 4 M, imply C–O core masses that are roughly consistent with that deduced from the ejecta mass plus a neutron star, 2 to 3 M. It is possible that the hybrid events are the result of Roche lobe overflow and that the pure events, SN Ib or SN Ic, result from common envelope evolution.     

Nucleosynthesis in massive stars

We discuss three aspects of the nucleosynthesis in massive and intermediate-mass stars during their early evolutionary phases. These are related to the CNO abundances in giant or supergiant stars, to the²⁶Al yield from massive stars via stellar wind, and to the production of the s-process nuclei in massive stars.     

Observational Appearance of Black Holes in X-Ray Binaries and AGN

Accretion onto black holes powers most luminous compact sources in the Universe. Black holes are found with masses extending over an extraordinary broad dynamic range, from several to a few billion times the mass of the Sun. Depending on their position on the mass scale, they may manifest themselves as X-ray binaries or active galactic nuclei. X-ray binaries harbor stellar mass black holes—endpoints of the evolution of massive stars. They have been studied by X-ray astronomy since its inception in the early 60-ies, however, the enigma of the most luminous of them—ultra-luminous X-ray sources, still remains unsolved. Supermassive black holes, lurking at the centers of galaxies, are up to hundreds of millions times more massive and give rise to the wide variety of different phenomena collectively termed “Active Galactic Nuclei”. The most luminous of them reach the Eddington luminosity limit for a few billions solar masses object and are found at redshifts as high as z≥5–7. Accretion onto supermassive black holes in AGN and stellar- and (possibly) intermediate mass black holes in X-ray binaries and ultra-luminous X-ray sources in star-forming galaxies can explain most, if not all, of the observed brightness of the cosmic X-ray background radiation. Despite the vast difference in the mass scale, accretion in X-ray binaries and AGN is governed by the same physical laws, so a degree of quantitative analogy among them is expected. Indeed, all luminous black holes are successfully described by the standard Shakura-Sunyaev theory of accretion disks, while the output of low-luminosity accreting black holes in the form of mechanical and radiative power of the associated jets obeys to a unified scaling relation, termed as the “fundamental plane of black holes”. From that standpoint, in this review we discuss formation of radiation in X-ray binaries and AGN, emphasizing their main similarities and differences, and examine our current knowledge of the demographics of stellar mass and supermassive black holes.     

Tailored analyses of 24 Galactic WN stars

The fundamental properties of 24 Galactic WN stars are determined from analyses of their optical, UV and IR spectra using sophisticated model atmosphere codes (Hillier, 1987, 1990). Terminal velocities, stellar luminosities, temperatures, mass loss rates and abundances of hydrogen, helium, carbon, nitrogen and oxygen are determined. Stellar parameters are derived using diagnostic lines and interstellar reddenings found from fitting theoretical continua to observed energy distributions.Our results confirm that the parameters of WN stars span a large range in temperature (T_*=30–90,000 K), luminosity (log L_*/L=4.8–5.9), mass loss (M=0.9–12×10^–5 M yr^–1) and terminal velocity (v =630–3300 km s^–1). Hydrogen abundances are determined, and found to be low in WNEw and WNEs stars (<15% by mass) and considerable in most WNL stars (1–50%). Metal abundances are also determined with the nitrogen content found to lie in the range N/He=1–5×10^–3 (by number) for all subtypes, and C/N 0.02 in broad agreement with the predictions of Maeder (1991). Enhanced O/N and O/C is found for HD 104994 (WN3p) suggesting a peculiar evolutionary history. Our results suggest that single WNL+abs stars may represent an evolutionary stage immediately after the Of phase. Since some WNE stars exist with non-negligible hydrogen contents (e.g. WR136) evolution may proceed directly from WNL+abs to WNE in some cases, circumventing the luminous blue variable (LBV) or red supergiant (RSG) stage.     

Spectral analyses of Wolf-Rayet stars: Theory,results, conclusions

Stratified Non-LTE models for expanding atmospheres became available in the recent years. They are based on the idealizing assumptions of spherical symmetry, stationarity and radiative equilibrium. From a critical discussion we conclude that this standard model is basically adequate for describing real Wolf-Rayet atmospheres and hence can be applied for quantitative spectral analyses of their spectra.By means of these models, the fundamental parameters have been determined meanwhile for the majority of the known Galactic WR stars. Most of them populate a vertical strip in the Herzsprung-Russell diagram at effective temperatures of 35 kK, the luminosities ranging from 10^4.5 to 10^5.9 L . Only early-type WN stars with strong lines and WC stars are hotter. The chemical composition of WR atmospheres corresponds to nuclear-processed material (WN: hydrogen burning in the CNO cycle; WC: helium burning). Hydrogen is depleted but still detectable in the cooler part of the WN subclass.Different scenarios for the evolutionary formation of the Wolf-Rayet stars are discussed in the light of the empirical data provided from the spectral analyses. Post-red-supergiant evolution can principally explain the basic observational properties, except the rather low luminosities of a considerable fraction of WN stars. Among the alternative scenarios, close-binary evolution can theoretically produce the least-luminous WN stars. However, final conclusions about the evolutionary formation of the WR stars are not yet possible.     

SN 1987A and SN 1993J: Testing stellar evolution theory?

Nearby supernovae like SN 1987A and SN 1993J provide valuable constraints on the late evolution of massive stars. For this purpose, we review evolutionary models for the progenitor of SN 1987A and confront them with five observational/theoretical tests we devised. We show that single-star models (with the possible exception of rapid-rotation models) fail at least two of these tests, while two binary models (accretion and merger models) are consistent with all available constraints. We conclude that it is most likely that the progenitor of SN 1987A had a binary companion, either at the time of the explosion or at least in the not-too-distant past, and that SN 1987A should therefore not be used to calibrate single stellar evolution theory. For SN 1993J, we infer from the presupernova photometry and the early light curve that its progenitor was a 15M star that lost almost all of its hydrogen-rich envelope prior to the supernova. This seems to require that the progenitor underwent stable case C mass transfer. We discuss future observational tests of binary models for both supernovae.     

Mass Measurements of Stellar and Intermediate-Mass Black Holes

We discuss the method, and potential systematic effects therein, used for measuring the mass of stellar-mass black holes in X-ray binaries. We restrict our discussion to the method that relies on the validity of Kepler’s laws; we refer to this method as the dynamical method. We briefly discuss the implications of the mass distribution of stellar-mass black holes and provide an outlook for future measurements. Further, we investigate the evidence for the existence of intermediate-mass black holes i.e. black holes with masses above 100 M_⊙, the limit to the black hole mass that can be produced by stellar evolution in the current Universe.     

P Cygni: An Extraordinary Luminous Blue Variable

P Cygni is a prototype for understanding mass loss from massive stars. This textbook star is known first of all because of two great eruptions in the 17th century. In the first half of this century it has given its name to a class of stars which are characterized by spectral lines consisting of nearly undisplaced emissions accompanied by a blue-displaced absorption component. This characteristic P Cygni-type profile betrays the presence of a stellar wind, but P Cygni's wind is quite unlike that of other hot supergiants. P Cygni was the first star that showed the effects of stellar evoluton from a study of its photometric history. It shares some common properties with the so-called Luminous Blue Variables. However, P Cygni is a unique object. This review deals with P Cygni's photometric properties, its circumstellar environment - including infrared and radio observations - and its optical and ultraviolet spectrum. Smaller sections deal with P Cygni's wind structure and evolution. This revised version was published online in August 2006 with corrections to the Cover Date.     

A gentle process for the formation of Algols

This work is concerned with binary systems that we call moderately close. These are systems in which the primary (by which we mean the initially more massive star) fills its Roche lobe when it is on the giant branch with a deep convective envelope but before helium ignition (late case B). We find that if the mass ratio q(= M ₁/M ₂) < q _crit = 0.7 when the primary fills its Roche lobe positive feedback will lead to a rapid hydrodynamic phase of mass transfer which will probably lead to common envelope evolution and thence to either coalescence or possibly to a close binary in a planetary nebula. Although most Algols have probably filled their Roche lobes before evolving off the main-sequence we find that some could not have and are therefore moderately close. Since rapid overflow is unlikely to lead to an Algol-like system there must be some way of avoiding it. The most likely possibility is that the primary can lose sufficient mass to reduce q below q _crit before overflow begins. Ordinary mass loss rates are insufficient but evidence that enhanced mass loss does take place is provided by RS CVn systems that have inverted mass ratios but have not yet begun mass transfer. We postulate that the cause of enhanced mass loss lies in the heating of the corona by by magnetic fields maintained by an — dynamo which is enhanced by tidal effects associated with corotation. In order to model the the effects of enhanced mass loss we ignore the details and adopt an empirical approach calibrating a simple formula with the RS CVn system Z Her. Using further empirical relations (deduced from detailed stellar models) that describe the evolution of red giants we have investigated the effect on a large number of systems of various initial mass ratios and periods. These are notable in that some systems can now enter a much gentler Algol-like overflow phase and others are prevented from transferring mass altogether. We have also investigated the effects of enhanced angular momentum loss induced by corotation of the wind in the strong magnetic fields and consider this in relation to observed period changes. We find that a typical moderately close Algol-like system evolves through an RS CVn like system and then possibly a symbiotic state before becoming an Algol and then goes on through a red giant-white dwarf state which may become symbiotic before ending up as a double white dwarf system in either a close or wide orbit depending on how much mass is lost before the secondary fills its Roche lobe.