Active OB : Laboratories for Stellar and Circumstellar Physics ASP Conference Series, Vol. 361, 2007 S. Stefl,ˇ S. P. Owocki and A. T. Okazaki

Nebulae from Eruptions of Luminous Evolved Stars: , RY Scuti, and the LBVs

N. Smith CASA, University of Colorado, 389 UCB, Boulder, CO 80309

Abstract. The most prodigious mass loss for luminous hot stars occurs during the (LBV) phase in transition to a Wolf-Rayet . Most of the mass loss is a result of a few brief eruptions, rather than a steady wind. For the most luminous stars, these eruptions eject several M⊙ at once, accounting for a large fraction of their total post–main-sequence mass loss. The geometry of their nebulae in the young free expansion phase traces the roles of rotation and binary interactions. Our most observable example is the around η Car, while nebulae around the eclipsing binary RY Sct and other LBVs share similar but less extreme properties. Both η Car and RY Sct have nebulae less than 200 yrs old with pronounced axial symmetry.

1. Introduction: High Luminosities and Circumstellar Nebulae

In this contribution, I would like to explore the results of stellar activity in stars with temperatures similar to Be stars, but with much higher luminosities. Various types of stellar activity and non-spherical mass loss discussed at this meeting are present in Be stars largely as a consequence of lower geff induced by very rapid rotation rates. In high-luminosity B supergiants, especially the luminous blue variables (LBVs), extreme mass-loss and the effects of rotation follow from lower geff due to the proximity to the Eddington limit and near- critical rotation, even if the stars do not have high rotational velocities. The most spectacular manifestation of this is the short-duration outbursts of LBVs, which may dominate the post-MS mass loss in luminous stars. While we cannot study rotation rates directly (i.e. broadened lines) in these stars due to their dense winds, we can infer the importance of rotation from their circumstellar geometry. However, we must catch these nebulae soon after ejection, while still in the free-expansion phase. The most well-studied case is the 160-yr old around η Car, but I will also mention a few other relevant objects.

2. Eta Carinae and the Homunculus

The around η Car is perhaps the most recognizable bipo- lar nebula seen by HST (Morse et al. 1998), and its 19th century eruption is arguably the most extreme example of stellar activity, barring the final destruc- tion of a massive star. That outburst ejected more than 10 M⊙ of material and released more kinetic energy than SN1054 (Smith et al. 2003b). The mass loss rate was so high that it could not have been a conventional line-driven wind because the material was optically thick, and so the mass-loss was either a super- 200 Nebulae from Eruptions 201

Figure 1. Dust structure in the Homunculus. (a) 18 µm image. (b) Color temperature. (c) Warm dust column density at 8.8 µm. (d) Cool dust column density at 18 µm (see Smith et al. 2003b).

Figure 2. Spatio-kinematic structure of (a) [Fe ii] λ16435 and (b) H2 2.122 µm, with the long-slit aperture along the major axis of the Homunculus, but offset slightly from the star; from Smith (2002). 202 Smith

Figure 3. High resolution spectrum of [Fe ii] λ16435 in P Cygniprofiles’s nebula, with the long-slit aperture oriented north-south (Smith & Hartigan 2006).

Eddington continuum-driven wind (e.g., Owocki et al. 2004) or a hydrodynamic explosion (Arnett et al. 2005). The huge amount of mass ejected suggests that these outbursts may be important in . Recent study of the Homunculus has revealed a double-shell excitation structure (see Figs. 1 and 2), with a thin outer shell seen in cool dust and H2 emission, and a thicker inner shell seen in warmer dust and [Fe ii] emission (Smith et al. 2003b; Smith 2002, 2005). Most of the mass is contained in the thin and cool H2-emitting shell, whose bipolar shape and pinched waist point toward rapid rotation of the erupting star as the shaping mechanism (Smith 2002). The nebula is too dense and massive to have been shaped by a pre-existing torus, and proper motions of the expanding nebula are ballistic. Subsequent mass ejections have followed the same bipolar symmetry axis: for example, the 1890 eruption that formed the “Little Homunculus” (Ishibashi et al. 2003; Smith 2005) and the present-day stellar wind (Smith et al. 2003a). Eta Car is probably a binary system, but the specific role of the companion star is unclear and continues to be investigated. It is unlikely, however, that the putative companion star played a major role in shaping the ejecta of the 19th century eruption, since the kinetic energy of the Homunculus is greater than the gravitational binding energy of the orbit (Smith et al. 2003b). Of course, recurring close periastron passages in this eccentric system may have spun up the outer layers of the star (e.g., Smith et al. 2003a).

3. P Cygniprofiles: A Counter-Example

P Cygniprofiles is the only other LBV in our that has been observed during a giant outburst, back in 1600 A.D. (Humphreys et al. 1999). Its nebula ejected during that eruption (Smith & Hartigan 2006) is older and less massive than η Car’s, and so is more prone to shaping by the ambient medium. Unlike the Homunculus of η Car, the nebula around P Cygniprofiles appears to be spherical. Nebulae from Eruptions 203

It is a challenge to image this nebula because of the bright central star, but the spherical shape can be seen in high-resolution spectra, like the [Fe ii] data in Figure 3 (Smith & Hartigan 2006) or in similar [Ni ii] or [N ii] spectra (Barlow et al. 1994). The spherical shape is in stark contrast to the bipolar Homunculus of η Car. This implies that the effects of rotation in shaping outflows will be more severe at higher luminosities closer to the Eddington limit, since P Cyg is far 5.8 less luminous than η Car (P Cygniprofiles is ∼10 L⊙ at D=1.7 kpc (Najarro et al. 1997). The mildly bipolar or elliptical shape of AG Car’s nebula would support this view, since its luminosity is intermediate between η Car and P Cyg. The ejection speed and radiative luminosity during P Cygni’s outburst were less than for η Car. Recently, (Smith & Hartigan 2006) estimated the mass of P Cygni’s nebula to be ∼0.1 M⊙, which shows that the total mass ejected, the mass loss rate, and the kinetic energy during this eruption were also far less extreme than for η Car (although interestingly similar to η Car’s 1890 outburst which created the Little Homunculus (Smith 2005; Ishibashi et al. 2003). This illustrates the likely range of mass ejection in LBV outbursts, since P Cyg is at the lower luminosity end of LBVs that are not post-RSGs, relevant to the distribution of LBVs on the HR diagram (e.g., Smith et al. 2004).

a) HST/WFPC2 b) VLA Hα 15 GHz (1997) (1992)

1" 1"

c) Keck/LWS d) Flux ratio 11.7 µm 11.7 µm (1999) Hα

2" 2"

Figure 4. RY Scuti’s nebula. (a) HST/WFPC2 image in Hα from (Smith et al. 2002). (b) 2 cm free-free continuum image taken with the VLA (Smith et al. 2001a; Gehrz et al. 1995)). (c) Dust emission at 11.7 µm (Gehrz et al. 2001); note the difference in size compared to Panel a. (d) Flux ratio of panels a and c, showing the spatial separation between ionized gas and dust. 204 Smith

Figure 5. Left: Long-slit HST/STIS spectrum of [N ii] λ6583 in RY Scuti’s nebula, with the slit along the equatorial (major) axis. This shows the uneven brightness distribution around the ring, indicating severe azimuthal asymme- try (Smith et al. 2002). Right: Proper motions of RY Scuti’s nebula, showing the increase in the diameter of the nebula (Smith et al. 2001a).

4. RY Scuti, Equatorial Rings, and SN1987a

RY Scuti is a luminous WR+OB eclipsing binary with an 11 day period and a separation of about 0.4 AU (Smith et al. 2002; Gehrz et al. 1995). What makes it interesting for trying to understand how binary stars shape outflows is its unusual circumstellar nebula (see Fig. 4). RY Scuti’s nebula is an equatorial dust torus, with the inner edge of the torus ionized by the central star (1995 Gehrz et al. 2001; Smith et al. 2002). The ionized nebula is especially weird, as it seems to show a pair of plane-parallel ionized rings with the same radius above and below the equator, while kinematics in spectra (Fig. 5) reveal severe azimuthal asymmetry (Smith et al. 2002). This geometry is likely to be a direct imprint of the mass-ejection geometry, since the expanding nebula was ejected sometime in the 19th century (Smith et al. 2001a, Fig. 5). The total mass of ionized gas in RY Scuti’s nebula is about 0.003 M⊙ (Smith et al. 2002). Thus, −1 the mass-loss rate to create the ionized nebula is 0.003/(∆tyr) M⊙ yr . For −1 ∆tyr∼<11 days, that’s M˙ ∼>0.1 M⊙ yr , which is extremely high. The mass in RY Scuti’s nebula is concentrated near the equatorial plane, implying that mass loss through the outer Lagrangian points in this over-contact system is the dominant mechanism shaping the outflow. We can then attempt to understand the unique double-ring structure in the nebula as a consequence of close binary interaction, although no satisfactory explanation exists. A similar conundrum exists for the formation of the famous triple-ring neb- ula around SN1987a. The current paradigm to explain SN1987a’s nebula is that a fast blue supergiant wind expanded into a previously-deposited RSG wind with an equatorial density enhancement, causing the formation of a swept-up Nebulae from Eruptions 205 equatorial ring and bipolar nebula (Blondin & Lundvist 1993; Martin & Arnett 1995). The strong density contrast in the RSG wind would have required a rapidly rotating RSG, which in turn, would probably require a binary merger. The problem is that simulations justifying this scenario explain the triple rings as the limb-brightened edges of an hourglass-shaped bipolar nebula, whereas the nebula really is three distinct rings and not a limb-brightened hourglass (Sug- erman et al. 2005). Furthermore, the kinematics of the nebula suggest that the three rings are coeval, with an age of ∼104 (Crotts & Heathcote 2000). While I don’t presently offer a specific alternative, I simply note that there’s room for additional work to understand the formation of the nebula ejected by SN1987a’s blue supergiant progenitor, and that the geometries seen in LBVs may point in a fruitful direction. This paradigm of rapid rotation arising from a red supergiant swallowing a companion will also have trouble for the ring and bipolar nebula around the blue supergiant , sometimes described as a “twin” of SN1987a, because Sher 25 is too luminous to have been a RSG if it is at the same distance as NGC 3603 (Smartt et al. 2002).

5. The Role of Extreme Mass-Loss Events in Stellar Evolution

Although its physical properties are extreme, η Car is not alone in having giant LBV outbursts. In addition to P Cyg’s eruption in 1600 A.D., a number of these events have been observed in other : SN1954j (or V12 in NGC2403, Tammann & Sandage 1968; Smith et al. 2001b), SN1961v (van Dyk et al. 2000), and a number of more recent “ impostors” that have been classified as Type IIn SNe (e.g., Filippenko 2005). Many Galactic LBVs and candidates show prominent ring nebulae that probably resulted from similar outbursts. Although it may be difficult to draw conclusions about the roles of rotation and binary interactions from these older nebulae and unresolved extragalactic systems, the distribution of nebular mass and corresponding stellar luminosity will be important for understanding the role of this extreme mass-loss activity in the broader context of stellar evolution. Acknowledgments. I thank the SOC and LOC for an enjoyable and mem- orable meeting, and for partial financial support. My research is supported by NASA through grant HF-01166 from the Space Telescope Science Institute.

References

Arnett, D., Meakin, C., & Young, P. 2005, ASP Conf. Ser. 332, 75 Barlow, M.J., Drew, J.E., Meaburn, J., & Massey, R.M. 1994, MNRAS, 268, L29 Blondin, J.M., & Lundqvist 1993, ApJ, 405, 337 Crotts, A.P.S., & Heathcote, S.R. 2000, ApJ, 528, 426 Filippenko, A. V. 2005, ASP Conf. Ser. 332, 33 Gehrz, R. D., Smith, N., Jones, B., Puetter, R., & Yahil, A. 2001, ApJ, 559, 395 Gehrz, R. D., Hayward, T. L., Houck, J. B., et al. 1995, ApJ, 439, 417 Humphreys, R. M., Davidson, K., & Smith, N. 1999, PASP, 111, 1124 Ishibashi, K., Gull, T. R., Davidson, K., et al. 2003, AJ, 125, 3222 Martin, C. L., & Arnett, D. 1995, ApJ, 447, 378 Morse, J. A., Davidson, K., B., John, et al. 1998, AJ, 116, 24443 Najarro, F., Hillier, D. J., & Stahl, O. 1997, A&A, 326, 1117 206 Smith

Owocki, S. P., Gayley, K., & Shaviv, N. 2004, ApJ, 616, 525 Smartt, S., Lennon, D. J., Kudritzki, R. P., et al. 2002, A&A, 391, 979 Smith, N. 2002, MNRAS, 337, 1252 Smith, N. 2005, MNRAS, 357, 1330 Smith, N., Gehrz, R. D.; Stahl, O.. Balick, B., Kaufer, A. 2002, ApJ, 578, 464 Smith, N., Davidson, K., Gull, T. R., Ishibashi, K., Hillier, D. J. 2003a, ApJ, 586, 432 Smith, N., Gehrz, R. D., Hinz, P. M., et al. 2003b, AJ, 125, 1458 Smith, N., Gehrz, R. D., & Goss, W. M. 2001a, AJ, 122, 2700 Smith, N., & Hartigan, P. 2006, ApJ, in press (astro-ph/0510836) Smith, N., Morse, J. A., Collins, N., & Gull, T. R. 2004, ApJ, 610, L105 Smith, N., Humphreys, R.M., & Gehrz, R.D. 2001b, PASP, 113, 692 Smith, N., Vink, J., & de Koter, A. 2004, ApJ, 615, 475 Sugerman, B.E.K., et al. 2005, ApJ, 627, 888 Tammann, G.A., & Sandage, A. 1968, ApJ, 151, 825 van Dyk, S.D., Peng, C. Y., King, J. Y., et al. 2000, PASP, 112, 1532

Discussion

D. Baade: What do we know about the duration of the outburst of RY Scuti? Is it highly unlikely that it lasted for less than an orbital period? N. Smith: The orbital period is about 11 days, so the dominant mass ejection episode(s) would need to have been very short indeed. If the ionized nebula we see today had been ejected in a few days time, the instantaneous mass-loss rate −1 would need to have been about 0.1 M⊙ yr . This seems rather high, but still less than η Car’s. For such a high mass-loss rate, though, the star should have been very bright in the 19th century when the nebula was ejected. O. Chesneau: Concerning RY Scuti, in PNe or Herbig AeBe disk simulations, the disk flaring shadow can easily mimic rings in an edge-on disk. N. Smith: Yes, I think an illumination effect might be an alternative explanation for the appearance of the “rings” around RY Sct; here it would need to be an unresolved circumbinary disk casting a Lyman continuum shadow, giving the illusion of a gap between the rings. That type of effect won’t explain the rings around SN1987a, though, so we still need a mechanism to eject parallel rings. R. Townsend: How do the rings of SN1987a move? Do they expand radially, or perpendicular to the rotation axis? N. Smith: When nebulae reach these distances from the star, all the motion is essentially radial. Observed Doppler shifts suggest that the 3 rings are coeval. N. Miller: You briefly alluded to recent data regarding the binarity of η Car. I was wondering if you would elaborate on that. N. Smith: I was trying to avoid this. It’s been an ongoing debate, to say the least...to me, the most convincing evidence for binarity is the X-ray variations studied by Corcoran et al., which look a lot like a colliding-wind binary. The pictures I showed in my talk were UV images taken with HST, which show suspicious moving shadows that change their direction before and after perias- tron passages. Those sharp directional changes are more easily explained with a binary system (see Smith et al. 2004).