On the role of continuum-driven eruptions in the evolution of very massive stars and Population III stars

We suggest that the mass lost during the evolution of very massive stars may be dominated by optically thick, continuum-driven outbursts or explosions, instead of by steady line-driven winds. In order for a massive star to become a WR star, it must shed its H envelope, but new estimates of the effects of clumping in winds indicate that line driving is vastly insufficient. We discuss massive stars above roughly 40-50 Msun, for which the best alternative is mass loss during brief eruptions of luminous blue variables (LBVs). Our clearest example of this phenomenon is the 19th century outburst of eta Car, when the star shed 12-20 Msun or more in less than a decade. Other examples are circumstellar nebulae of LBVs, extragalactic eta Car analogs (``supernova impostors''), and massive shells around SNe and GRBs. We do not yet fully understand what triggers LBV outbursts, but they occur nonetheless, and present a fundamental mystery in stellar astrophysics. Since line opacity from metals becomes too saturated, the extreme mass loss probably arises from a continuum-driven wind or a hydrodynamic explosion, both of which are insensitive to metallicity. As such, eruptive mass loss could have played a pivotal role in the evolution and fate of massive metal-poor stars in the early universe. If they occur in these Population III stars, such eruptions would profoundly affect the chemical yield and types of remnants from early SNe and hypernovae.Comment: 4 pages, 1 figure, accepted by ApJ Letter

Crossing the `Yellow Void' -- Spatially Resolved Spectroscopy of the Post- Red Supergiant IRC+10420 and Its Circumstellar Ejecta

IRC +10420 is one of the extreme hypergiant stars that define the empirical upper luminosity boundary in the HR diagram. During their post--RSG evolution, these massive stars enter a temperature range (6000-9000 K) of increased dynamical instability, high mass loss, and increasing opacity, a semi--forbidden region, that de Jager and his collaborators have called the `yellow void'. We report HST/STIS spatially resolved spectroscopy of IRC +10420 and its reflection nebula with some surprising results. Long slit spectroscopy of the reflected spectrum allows us to effectively view the star from different directions. Measurements of the double--peaked Halpha emission profile show a uniform outflow of gas in a nearly spherical distribution, contrary to previous models with an equatorial disk or bipolar outflow. Based on the temperature and mass loss rate estimates that are usually quoted for this object, the wind is optically thick to the continuum at some and possibly all wavelengths. Consequently the observed variations in apparent spectral type and inferred temperature are changes in the wind and do not necessarily mean that the underlying stellar radius and interior structure are evolving on such a short timescale. To explain the evidence for simultaneous outflow and infall of material near the star, we propose a `rain' model in which blobs of gas condense in regions of lowered opacity outside the dense wind. With the apparent warming of its wind, the recent appearance of strong emission, and a decline in the mass loss rate, IRC +10420 may be about to shed its opaque wind, cross the `yellow void', and emerge as a hotter star.Comment: To appear in the Astronomical Journal, August 200

The Missing Luminous Blue Variables and the Bistability Jump

We discuss an interesting feature of the distribution of luminous blue variables on the H-R diagram, and we propose a connection with the bistability jump in the winds of early-type supergiants. There appears to be a deficiency of quiescent LBVs on the S Dor instability strip at luminosities between log L/Lsun = 5.6 and 5.8. The upper boundary, is also where the temperature-dependent S Dor instability strip intersects the bistability jump at about 21,000 K. Due to increased opacity, winds of early-type supergiants are slower and denser on the cool side of the bistability jump, and we postulate that this may trigger optically-thick winds that inhibit quiescent LBVs from residing there. We conduct numerical simulations of radiation-driven winds for a range of temperatures, masses, and velocity laws at log L/Lsun=5.7 to see what effect the bistability jump should have. We find that for relatively low stellar masses the increase in wind density at the bistability jump leads to the formation of a modest to strong pseudo photosphere -- enough to make an early B-type star appear as a yellow hypergiant. Thus, the proposed mechanism will be most relevant for LBVs that are post-red supergiants. Yellow hypergiants like IRC+10420 and rho Cas occupy the same luminosity range as the ``missing'' LBVs, and show apparent temperature variations at constant luminosity. If these yellow hypergiants do eventually become Wolf-Rayet stars, we speculate that they may skip the normal LBV phase, at least as far as their apparent positions on the HR diagram are concerned.Comment: 20 pages, 4 figs, accepted by Ap

A New View of the Circumstellar Environment of SN 1987A

We summarize the analysis of a uniform set of both previously-known and newly-discovered scattered-light echoes, detected within 30" of SN 1987A in ten years of optical imaging, and with which we have constructed the most complete three-dimensional model of the progenitor's circumstellar environment. Surrounding the SN is a richly-structured bipolar nebula. An outer, double-lobed ``peanut,'' which we believe is the contact discontinuity between the red supergiant and main sequence winds, is a prolate shell extending 28 ly along the poles and 11 ly near the equator. Napoleon's Hat, previously believed to be an independent structure, is the waist of this peanut, which is pinched to a radius of 6 ly. Interior, the innermost circumstellar material lies along a cylindrical hourglass, 1 ly in radius and 4 ly long, which connects to the peanut by a thick equatorial disk. The nebulae are inclined 41o south and 8o east of the line of sight, slightly elliptical in cross section, and marginally offset west of the SN. The 3-D geometry of the three circumstellar rings is studied, suggesting the equatorial ring is elliptical (b/a<0.98), and spatially offset in the same direction as the hourglass. Dust-scattering models suggest that between the hourglass and bipolar lobes: the gas density drops from 1--3 cm^{-3} to >0.03 cm^{-3}; the maximum dust-grain size increases from ~0.2 micron to 2 micron; and the Si:C dust ratio decreases. The nebulae have a total mass of ~1.7 Msun, yielding a red-supergiant mass loss around 5*10^{-6} Msun yr^{-1}.Comment: Accepted for publication in ApJ 2/14/05. 16 pages in emualteapj forma

The Evolution of Supernovae in Circumstellar Wind Bubbles II: Case of a Wolf-Rayet star

(Abridged) Mass-loss from massive stars leads to the formation of circumstellar wind-blown bubbles surrounding the star, bordered by a dense shell. When the star ends its life in a supernova (SN) explosion, the resulting shock wave will interact with this modified medium. In a previous paper we discussed the basic parameters of this interaction. In this paper we go a step further and study the evolution of SNe in the wind blown bubble formed by a 35 \msun star that starts off as an O star, goes through a red supergiant phase, and ends its life as a Wolf-Rayet star. We model the evolution of the CSM and then the expansion of the SN shock wave within this medium. Our simulations clearly reveal fluctuations in density and pressure within the surrounding medium. The SN shock interacting with these fluctuations, and then with the dense shell surrounding the wind-blown cavity, gives rise to a variety of transmitted and reflected shocks in the wind bubble. The interactions between these various shocks and discontinuities is examined, and its effects on the X-ray emission is noted. Our simulations reveal the presence of several hydrodynamic instabilities. They show that the turbulent interior, coupled with the large fluctuations in density and pressure, gives rise to an extremely corrugated SN shock wave. The shock shows considerable wrinkles as it impacts the dense shell, and the impact occurs in a piecemeal fashion, with some parts of the shock wave interacting with the shell before the others. Therefore different parts of the shell will `light-up' at different times. The non-spherical nature of the interaction means that it will occur over a prolonged period of time, and the spherical symmetry of the initial shock wave is destroyed.Comment: 50 pages, 19 figures. Accepted to the Astrophysical Journal. For a version with the original high-resolution color figures please download from http://astro.uchicago.edu/~vikram/sncsm.htm

How Massive Single Stars End their Life

How massive stars die -- what sort of explosion and remnant each produces -- depends chiefly on the masses of their helium cores and hydrogen envelopes at death. For single stars, stellar winds are the only means of mass loss, and these are chiefly a function of the metallicity of the star. We discuss how metallicity, and a simplified prescription for its effect on mass loss, affects the evolution and final fate of massive stars. We map, as a function of mass and metallicity, where black holes and neutron stars are likely to form and where different types of supernovae are produced. Integrating over an initial mass function, we derive the relative populations as a function of metallicity. Provided single stars rotate rapidly enough at death, we speculate upon stellar populations that might produce gamma-ray bursts and jet-driven supernovae.Comment: 24 pages, 9 figues, submitted to Ap

Presupernova Evolution of Rotating Massive Stars I: Numerical Method and Evolution of the Internal Stellar Structure

The evolution of rotating stars with zero-age main sequence (ZAMS) masses in the range 8 to 25 M_sun is followed through all stages of stable evolution. The initial angular momentum is chosen such that the star's equatorial rotational velocity on the ZAMS ranges from zero to ~ 70 % of break-up. Redistribution of angular momentum and chemical species are then followed as a consequence of rotationally induced circulation and instablities. The effects of the centrifugal force on the stellar structure are included. Uncertain mixing efficiencies are gauged by observations. We find, as noted in previous work, that rotation increases the helium core masses and enriches the stellar envelopes with products of hydrogen burning. We determine, for the first time, the angular momentum distribution in typical presupernova stars along with their detailed chemical structure. Angular momentum loss due to (non-magnetic) stellar winds and the redistribution of angular momentum during core hydrogen burning are of crucial importance for the specific angular momentum of the core. Neglecting magnetic fields, we find angular momentum transport from the core to the envelope to be unimportant after core helium burning. We obtain specific angular momenta for the iron core and overlaying material of 1E16...1E17 erg s. These values are insensitive to the initial angular momentum. They are small enough to avoid triaxial deformations of the iron core before it collapses, but could lead to neutron stars which rotate close to break-up. They are also in the range required for the collapsar model of gamma-ray bursts. The apparent discrepancy with the measured rotation rates of young pulsars is discussed.Comment: 62 pages, including 7 tables and 19 figures. Accepted by Ap