The evolution of massive stars

The evolution of stars with masses between 15 M ₀ and 100M ₀ is considered. Stars in this mass range lose a considerable fraction of their matter during their evolution.The treatment of convection, semi-convection and the influence of mass loss by stellar winds at different evolutionary phases are analysed as well as the adopted opacities.Evolutionary sequences computed by various groups are examined and compared with observations, and the advanced evolution of a 15M ₀ and a 25M ₀ star from zero-age main sequence (ZAMS) through iron collapse is discussed.The effect of centrifugal forces on stellar wind mass loss and the influence of rotation on evolutionary models is examined. As a consequence of the outflow of matter deeper layers show up and when the mass loss rates are large enough layers with changed composition, due to interior nuclear reactions, appear on the surface.The evolution of massive close binaries as well during the phase of mass loss by stellar wind as during the mass exchange and mass loss phase due to Roche lobe overflow is treated in detail, and the value of the parameters governing mass and angular momentum losses are discussed.The problem of the Wolf-Rayet stars, their origin and the possibilities of their production either as single stars or as massive binaries is examined.Finally, the origin of X-ray binaries is discussed and the scenario for the formation of these objects (starting from massive ZAMS close binaries, through Wolf-Rayet binaries leading to OB-stars with a compact companion after a supernova explosion) is reviewed and completed, including stellar wind mass loss.     

Clusters of Galaxies: Setting the Stage

Clusters of galaxies are self-gravitating systems of mass ∼10¹⁴–10¹⁵ h ⁻¹ M_⊙ and size ∼1–3h ⁻¹ Mpc. Their mass budget consists of dark matter (∼80%, on average), hot diffuse intracluster plasma (≲20%) and a small fraction of stars, dust, and cold gas, mostly locked in galaxies. In most clusters, scaling relations between their properties, like mass, galaxy velocity dispersion, X-ray luminosity and temperature, testify that the cluster components are in approximate dynamical equilibrium within the cluster gravitational potential well. However, spatially inhomogeneous thermal and non-thermal emission of the intracluster medium (ICM), observed in some clusters in the X-ray and radio bands, and the kinematic and morphological segregation of galaxies are a signature of non-gravitational processes, ongoing cluster merging and interactions. Both the fraction of clusters with these features, and the correlation between the dynamical and morphological properties of irregular clusters and the surrounding large-scale structure increase with redshift. In the current bottom-up scenario for the formation of cosmic structure, where tiny fluctuations of the otherwise homogeneous primordial density field are amplified by gravity, clusters are the most massive nodes of the filamentary large-scale structure of the cosmic web and form by anisotropic and episodic accretion of mass, in agreement with most of the observational evidence. In this model of the universe dominated by cold dark matter, at the present time most baryons are expected to be in a diffuse component rather than in stars and galaxies; moreover, ∼50% of this diffuse component has temperature ∼0.01–1 keV and permeates the filamentary distribution of the dark matter. The temperature of this Warm-Hot Intergalactic Medium (WHIM) increases with the local density and its search in the outer regions of clusters and lower density regions has been the quest of much recent observational effort. Over the last thirty years, an impressive coherent picture of the formation and evolution of cosmic structures has emerged from the intense interplay between observations, theory and numerical experiments. Future efforts will continue to test whether this picture keeps being valid, needs corrections or suffers dramatic failures in its predictive power.     

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

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.     

Wolf-Rayet stars in starbursts

Let us suppose that it is possible observationally to determine the number ratio of WR to O stars in a starburst galaxy (cf. e. g. Vacca &; Conti 1992) and that one can also have some information on the way the different WR subtypes are distributed (number ratios as WN/WR, WNL/WR etc ...), the question is, what can we deduce from these values on the burst of star formation which gave birth to these WR stars? Is it possible for instance to constrain the age of the burst (i.e. the time elapsed since the beginning of the burst of star formation), its intensity (i.e. the ratio of the star formation rate during the burst to that before the burst) or the metallicity of the cloud from which the stars formed? We present here models of starbursts based on the most recent models for single stars computed by the Geneva group and show that the study of the WR population in a starburst provides very useful insights on the age of the burst and on the metallicity of the star forming zone.     

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     

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.     

A spectroscopic analysis of B stars in the SMC cluster NGC 330

Medium resolution (2A/px) but high s/n spectra of approximately twenty of the brightest blue stars in the young open cluster NGC 330 in the SMC have been analyzed in order to determine their atmospheric parameters and the evolutionary status. Stellar parameters are determined by comparison with LTE and NLTE model atmosphere calculations and an HR diagram constructed. Luminosities of the sample stars lie in the range 4.0L _*/L )<5.0 and spectral types between O9 and late-B. The stars in our sample appear to define 4 groups: main-sequence B-stars (B2-B4), B-supergiants (B4) in a blue-loop phase of evolution, a small number of blue stragglers (O9-B0 near main-sequence stars) and a group of luminous giants (B1-B2) which reside in the so-called post main-sequence gap of the HR diagram. Furthermore, we have confirmed spectroscopically the very high incidence of Be stars in this cluster. Finally the almost complete absence of metal lines (at this resolution) is in keeping with the expected very low metallicity of the SMC.     

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.     

The Milky Way 3-Helium Abundance

We are making precise determinations of the abundance of the light isotope of helium, ³He. The ³He abundance in Milky Way sources impacts stellar evolution, chemical evolution, and cosmology. The abundance of ³He is derived from measurements of the hyperfine transition of ³He⁺ which has a rest wavelength of 3.46 cm (8.665 GHz). As with all the light elements, the present interstellar ³He abundance results from a combination of Big Bang Nucleosynthesis (BBNS) and stellar nucleosynthesis. We are measuring the ³He abundance in Milky Way H ii regions and planetary nebulae (PNe). The source sample is currently comprised of 60 H ii regions and 12 PNe. H ii regions are examples of zero-age objects that are young relative to the age of the Galaxy. Therefore their abundances chronicle the results of billions of years of Galactic chemical evolution. PNe probe material that has been ejected from low-mass (M≤ 2M _⊙) to intermediate-mass (M∼2–5M _⊙) stars to be further processed by future stellar generations. Because the Milky Way ISM is optically thin at centimeter wavelengths, our source sample probes a larger volume of the Galactic disk than does any other light element tracer of Galactic chemical evolution. The sources in our sample possess a wide range of physical properties (including object type, size, temperature, excitation, etc.). The ³He abundances we derive have led to what has been called “The ³He Problem”.     

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.     

Magnetic Fields in Massive Stars, Their Winds, and Their Nebulae

Massive stars are crucial building blocks of galaxies and the universe, as production sites of heavy elements and as stirring agents and energy providers through stellar winds and supernovae. The field of magnetic massive stars has seen tremendous progress in recent years. Different perspectives—ranging from direct field measurements over dynamo theory and stellar evolution to colliding winds and the stellar environment—fruitfully combine into a most interesting and still evolving overall picture, which we attempt to review here. Zeeman signatures leave no doubt that at least some O- and early B-type stars have a surface magnetic field. Indirect evidence, especially non-thermal radio emission from colliding winds, suggests many more. The emerging picture for massive stars shows similarities with results from intermediate mass stars, for which much more data are available. Observations are often compatible with a dipole or low order multi-pole field of about 1 kG (O-stars) or 300 G to 30?kG (Ap/Bp stars). Weak and unordered fields have been detected in the O-star ζ Ori A and in Vega, the first normal A-type star with a magnetic field. Theory offers essentially two explanations for the origin of the observed surface fields: fossil fields, particularly for strong and ordered fields, or different dynamo mechanisms, preferentially for less ordered fields. Numerical simulations yield the first concrete stable (fossil) field configuration, but give contradictory results as to whether dynamo action in the radiative envelope of massive main sequence stars is possible. Internal magnetic fields, which may not even show up at the stellar surface, affect stellar evolution as they lead to a more uniform rotation, with more slowly rotating cores and faster surface rotation. Surface metallicities may become enhanced, thus affecting the mass-loss rates.     

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.     

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.     

The stellar content of the Orion OB1 association

We present a photometric investigation, using the VBLUW system, of the stellar content of Orion OB1. Physical parameters (logg, logT _eff) for the stars are derived with the aid of model atmospheres. From these, visual extinctions, absolute magnitudes and distance moduli are derived. The distance moduli are used to determine membership for the stars in each of the subgroups and distances to the subgroups are calculated. The ages of the subgroups are derived through isochrone fitting and the IMF is derived for each subgroup. The energy deposited into the ISM through stellar winds and supernovae is calculated and compared to observed large scale features in the ISM around Orion OB1.     

Infrared irradiance calibration

Infrared astronomical measurements are calibrated against reference sources, usually primary standard stars that are, in turn, calibrated either by direct or indirect means. A direct calibration compares the star with a certified source, typically a blackbody. Indirect methods extrapolate a direct measurement of the flux at one wavelength to the flux at another. Historically, α Lyr (Vega) has been used as the primary standard as it is bright, easily accessible from the northern hemisphere, and is well calibrated in the visual. Until recently, the direct absolute infrared calibrations of α Lyr and those derived from the absolute solar flux scaled to the observed spectral energy distributions of solar type stars increasingly diverged with wavelength from those obtained using a model atmosphere to extrapolate the absolute visual flux of Vega into the infrared. The exception is the direct calibration by the 1996/97 Midcourse Space Experiment of the absolute fluxes for a number of the commonly used infrared standard stars, including Vega.In the mid-1980s, the Air Force Geophysics Laboratory began a program that led to the establishment of a network of stars with which to calibrate infrared space-based sensors. α Lyr and a CMa were adopted as the fundamental references and the absolute 1.2 to 35 µm infrared spectral energy distributions for the 616 secondary standard stars in the network were derived through spectral and photometric comparisons with the primary standards. The stars are also used for calibration at ground-based infrared observatories. For applications in which the network stars may not be bright enough, particularly at the longer infrared wavelengths, planets and the larger asteroids are used. Planets and asteroids move and rather sophisticated thermal modeling of the bodies is required to predict the disk-integrated brightness at a specific time with reasonable accuracy. The Infrared Space Observatory applied such a sophisticated ‘thermo-physical’ model to the largest asteroids to support calibration of the sensors to a claimed accuracy of within 5%. The AFRL program also created a spectral atlas of the brightest stars in the sky that, although they are variable, may be used for calibration if the large(r) attendant uncertainties are acceptable.This revised version was published online in July 2005 with a corrected cover date.     

Infrared Irradiance Calibration

Infrared astronomical measurements are calibrated against reference sources, usually primary standard stars that are, in turn, calibrated either by direct or indirect means. A direct calibration compares the star with a certified source, typically a blackbody. Indirect methods extrapolate a direct measurement of the flux at one wavelength to the flux at another. Historically, α Lyr (Vega) has been used as the primary standard as it is bright, easily accessible from the northern hemisphere, and is well calibrated in the visual. Until recently, the direct absolute infrared calibrations of α Lyr and those derived from the absolute solar flux scaled to the observed spectral energy distributions of solar type stars increasingly diverged with wavelength from those obtained using a model atmosphere to extrapolate the absolute visual flux of Vega into the infrared. The exception is the direct calibration by the 1996/97 Midcourse Space Experiment of the absolute fluxes for a number of the commonly used infrared standard stars, including Vega.In the mid-1980s, the Air Force Geophysics Laboratory began a program that led to the establishment of a network of stars with which to calibrate infrared space-based sensors. α Lyr and a CMa were adopted as the fundamental references and the absolute 1.2 to 35 µm infrared spectral energy distributions for the 616 secondary standard stars in the network were derived through spectral and photometric comparisons with the primary standards. The stars are also used for calibration at ground-based infrared observatories. For applications in which the network stars may not be bright enough, particularly at the longer infrared wavelengths, planets and the larger asteroids are used. Planets and asteroids move and rather sophisticated thermal modeling of the bodies is required to predict the disk-integrated brightness at a specific time with reasonable accuracy. The Infrared Space Observatory applied such a sophisticated ‘thermo-physical’ model to the largest asteroids to support calibration of the sensors to a claimed accuracy of within 5%. The AFRL program also created a spectral atlas of the brightest stars in the sky that, although they are variable, may be used for calibration if the large(r) attendant uncertainties are acceptable.     

Fundamental parameters of the W Serpentis stars

Application of digital cross-correlation spectroscopy to the spectra of the W Serpentis binaries SX Cas and RX Cas has allowed an accurate determination of the orbits and rotations of the (mass-losing) K-subgiant secondary components. The distortion of the primary radial-velocity curves due to the influence of the prominent accretion disks in these systems has been modelled to first order. This enables us to estimate k ₁, and thereby the mass ratio q 0.30, to within ± 20%. The absolute radii of the secondaries are derived independently from the observed rotations and periods, assuming synchronous rotation. They show that the stars fill their Roche lobes, or at least very nearly so. Rough fits to the available photometry shows both primaries to be unevolved mid-B stars; that in RX Cas appears completely obscured by the disk. Preliminary spectroscopic data for W Ser and W Cru show some promise for similar analyses of these systems.