- Split View
-
Views
-
Cite
Cite
Nathan Smith, W. David Arnett, John Bally, Adam Ginsburg, Alexei V. Filippenko, The ring nebula around the blue supergiant SBW1: pre-explosion snapshot of an SN 1987A twin, Monthly Notices of the Royal Astronomical Society, Volume 429, Issue 2, 21 February 2013, Pages 1324–1341, https://doi.org/10.1093/mnras/sts418
- Share Icon Share
Abstract
SBW1 is a B-type supergiant surrounded by a ring nebula that is a nearby twin of SN 1987A's progenitor and its circumstellar ring. We present images and spectra of SBW1 obtained with the Hubble Space Telescope (HST), the Spitzer Space Telescope and Gemini South. HST images of SBW1 do not exhibit long Rayleigh–Taylor (RT) fingers, which are presumed to cause the ‘hotspots’ in the SN 1987A ring when impacted by the blast wave, but instead show a geometrically thin (ΔR/R ≲ 0.05) clumpy ring. The radial mass distribution and size scales of inhomogeneities in SBW1’s ring closely resemble those in the SN 1987A ring, but the more complete disc expected to reside at the base of the RT fingers is absent in SBW1. This structure may explain why portions of the SN 1987A ring between the hotspots have not yet brightened, more than 15 years after the first hotspots appeared. The model we suggest does not require a fast wind colliding with a previous red supergiant wind, because a slowly expanding equatorial ring may be ejected by a rotating blue supergiant star or in a close binary system. More surprisingly, high-resolution images of SBW1 also reveal diffuse emission filling the interior of the ring seen in Hα and in thermal-infrared (IR) emission; ∼190 K dust dominates the 8–20 μm luminosity (but contains only 10−5 M⊙ of dust). Cooler (∼85 K) dust resides in the equatorial ring itself (and has a dust mass of at least 5 × 10−3 M⊙). Diffuse emission extends inward to ∼1 arcsec from the central star, where a paucity of Hα and IR emission suggests an inner hole excavated by the B-supergiant wind. We propose that diffuse emission inside the ring arises from an ionized flow of material photoevaporated from the dense ring, and its pressure prevents the B-supergiant wind from advancing in the equatorial plane. This inner emission could correspond to a structure hypothesized to reside around Sk−69°202 that was never directly detected. If this interpretation is correct, it would suggest that photoionization can play an important dynamical role in shaping the ring nebula, and we speculate that this might help explain the origin of the polar rings around SN 1987A. In effect, the photoevaporative flow shields the outer bipolar nebula at low latitudes, whereas the blue supergiant wind expands freely out the poles and clears away the polar caps of the nebula; the polar rings reside at the intersection of these two zones.
INTRODUCTION
SN 1987A in the Large Magellanic Cloud (LMC) is one of those vexing examples of Murphy's law where the nearest and best-studied supernova (SN) in modern times also appears unique when compared to the population of known extragalactic core-collapse supernovae (SNe). It is an oddball, consisting of a peculiar SN type having a slow rise to peak, attributed to the relatively compact nature of its (also unusual and unexpected) blue supergiant (BSG) progenitor (Arnett 1987, 1989; Arnett et al. 1989). SN 1987A is the only SN near enough to obtain a clear picture of its spatially resolved circumstellar nebula, but its bizarre triple-ring nebula still defies adequate explanation (Luo & McCray 1991; Blondin & Lundqvist 1993; Martin & Arnett 1995; Collins et al. 1999; Morris & Podsiadlowski 2009), and does not conform to the hourglass structures commonly seen in bipolar planetary nebulae and bipolar nebulae around massive stars.
There are, however, a few objects known that appear similar to SN 1987A in that they have prominent equatorial ring nebulae with bipolar lobes or rings surrounding BSG central stars. In this paper, we present a detailed study of the recently discovered ring nebula SBW1 (Smith, Bally & Walawender 2007). Two other known bipolar ring nebulae around massive stars are HD 168625 and Sher 25, which are discussed in detail elsewhere (Brandner et al. 1997; Smith 2007). Of these three, SBW1 appears to most closely resemble the nebula around SN 1987A and the physical parameters of its progenitor star.
Before massive stars explode as SNe, they can shed considerable mass as they attempt to remove their outer H-rich envelope. If they fail to shed their H envelopes, they will typically remain as red supergiants (RSGs) and will die as SNe of Type II. One expects single stars in the lower initial mass ranges of core-collapse SNe (8–20 M⊙) to retain their H envelopes and explode as RSGs, and this expectation is largely confirmed by pre-explosion detections of progenitor stars of SNe II-P (e.g. Smartt 2009; Leonard 2011, and references therein). SN 1987A challenged our view of stellar evolution because it was an SN II from an explosion of a BSG, not a RSG, with an initial mass of 18–20 M⊙ (see Arnett et al. 1989, and references therein). The reason that the progenitor was a BSG is still unclear, but various scenarios involving close binary evolution, binary mergers, rapid rotation, and enhanced mass loss have been suggested.
Pre-explosion data for SN1987A's progenitor star Sk−69°202 establish that it appeared to be a fairly normal B3 I supergiant (Rousseau et al. 1978; Walborn et al. 1989). Therefore, our conjectures about the progenitor and its pre-SN mass loss depend heavily on studies of the remarkable ring nebula surrounding the SN, made famous in early Hubble Space Telescope (HST) images (e.g. Burrows et al. 1995; Plait et al. 1995). Kinematic studies of the nebula's expansion indicate that it was ejected by the progenitor star roughly 104 yr before exploding (Meaburn et al. 1995; Crotts & Heathcote 2000). The ejection and shaping mechanisms of this nebula are intimately linked to the star's peculiar evolution just before explosion, but our understanding of that process is still tenuous. The poor understanding of how the ring nebula formed has become a more pressing problem in recent years. We are now lucky enough to witness a spectacular collision as the blast wave of the SN overtakes the ring nebula ejected by the progenitor, predicted shortly after the discovery of the nebula (Luo & McCray 1991).
The SN 1987A blast wave first began to collide with the dense circumstellar ring in 1997, heralded by the appearance of new ‘hot spots’ in the ring (Michael et al. 1998, 1999; Sonneborn et al. 1998). These hotspots had broader linewidths than the rest of the ring, confirming that they were bright because a shock was being driven into the dense clumps. These clumps were thought to be the ends of long ‘fingers’ created by Rayleigh–Taylor (RT) instabilities at the contact discontinuity between the slow RSG wind and the fast BSG wind. Since the hotspots first appeared in 1997, many more spots have brightened all around the ring (Sugerman et al. 2002), although the hotspots have not yet merged into a contiguous bright ring as one might expect when the blast wave catches up to gas in between the RT fingers.
In the decade between explosion and the start of this collision with the ring, the blast wave was expanding through the relatively low density region interior to the ring. Radio emission and hydrodynamic models suggest that the BSG wind had a surprisingly low mass-loss rate of ∼10−7 M⊙ yr−1 (Blondin & Lundqvist 1993; Stavely-Smith et al. 1993; Chevalier & Dwarkadas 1995; Martin & Arnett 1995). The expansion rate of the blast wave through this low-density progenitor wind was fast at first, but then slowed (Gaensler et al. 2000), attributed to the shock running into a higher-density H ii region caused by photoionized material from the dense RSG wind (Chevalier & Dwarkadas 1995; Meyer 1997). While the existence of this H ii region can account for some observed characteristics of the blast wave expansion, emission from this feature itself was not directly observed before it was hit by the blast wave. One of the key results from our analysis below is that such a feature is seen directly in the similar nebula around SBW1.
A different approach to shed light on the pre-SN evolution of SN 1987A is to study nearby analogues of SN 1987A's progenitor that have not yet exploded. As noted above, three possible cousins of SN 1987A's progenitor are known: Sher 25, HD 168625, and SBW1. Of these, SBW1 appears the most similar in terms of the nebular structure and the luminosity of the central star, but it is not as extensively studied as the other two. Smith et al. (2007) first discovered SBW1 and performed the initial study of its nebula. While it is seen projected in the Carina Nebula star-forming region, its positive radial velocity suggests that it is actually located at a much larger distance behind the Carina Nebula; a distance of ∼7 kpc rather than the well-established 2.3 kpc distance to Carina (Smith 2006) also provides a better match between the expected luminosity from the B1.5 Iab spectral type and the observed magnitude and relatively low reddening. This distance makes the luminosity of SBW1 comparable to that of Sk−69°202, suitable for an 18–25 M⊙ progenitor star. From our more detailed analysis below, we find that SBW1 is indeed a virtual twin of the progenitor of SN 1987A and its circumstellar environment.
OBSERVATIONS
HST/WFC3 images
Following discovery in our ground-based observations (Smith et al. 2007), we imaged SBW1 using the newly installed WFC3 camera onboard HST on 2009 Dec. 8. (ut dates are used throughout this paper; see Table 1.) We used three filters: F502N sampling [O iii] λ5007, F656N sampling Hα, and F658N sampling [N ii] λ6584. We employed standard image-reduction techniques, and the resulting monochromatic F502N, F656N and F658N images are shown in Figs 1(a)–(c), respectively, whereas a colour composite of the three is shown in Fig. 1(d).
The HST/WFC3 images detect no [O iii] emission from the ring nebula itself. This suggests that there is no nearby source of hard-ultraviolet (UV) photons (>35 eV) that can ionize O+ to O++ in the circumstellar gas. There are apparently plenty of hard ionizing photons from the much more massive early O-type stars in proximate regions of the Carina Nebula (also evidenced by the bright diffuse [O iii] emission seen around SBW1), so this provides yet another argument that SBW1 is not actually located within the Carina Nebula, but is instead far behind it and seen in projection (Smith et al. 2007). Similarly, we detect no features in absorption associated with the circumstellar dust in SBW1, contrary to expectations for dust features in silhouette against a background screen of an H ii region.
The Hα and [N ii] images reveal emission structures that are essentially identical. In a F656N − F658N difference image (not shown), the nebula vanishes almost completely, except for a small wisp of emission in the diffuse inner part of the ring at the 1 per cent level, which could easily be attributable to some foreground absorption or emission in the Carina Nebula. As we will see below, our analysis of the continuum-subtracted emission in HST Space Telescope Imaging Spectrograph (STIS) spectra shows no variation in the Hα/[N ii] flux ratio across the ring, consistent with flux ratios measured in images (STIS spectra are preferred for this comparison, since scattered starlight can be subtracted). Thus, the gas in the dense and thin equatorial ring, the gas in the outer bipolar regions and the more diffuse gas filling the interior of the ring are likely to all have the same relative N/H abundance. This is important for some models of the formation of the nebula (see below).
The most interesting results of the HST imaging are (1) the detailed structure in the equatorial ring itself, seen at roughly seven times the effective spatial resolution (not angular resolution) provided by HST images of the more distant SN 1987A ring and (2) diffuse emission structures inside the ring, which are associated with hot dust features (see below). Both of these provide important clues about the structure and origin of the nebula into which the SN 1987A blast wave has been expanding. The specific structures and implications are discussed in Sections 3 and 4.
HST/STIS spectra
The SBW1 ring was also observed with HST/STIS on 2010 May 18 (ut; see Table 1). The 52 × 0.2 arcsec2 slit aperture was oriented along position angle +37°, at two offset positions on either side of the central star as shown in Fig. 2.
Fig. 3(a) indicates that there is a slight velocity gradient across the ring, such that the NE side of the equatorial ring is redshifted by 10–20 km s−1, as was seen in the higher-resolution echelle spectra of Smith et al. (2007). Also, the entire ring centroid is redshifted by ∼20 km s−1 (the dashed line shows v = 0 for Hα), confirming the systemic-velocity measurement by Smith et al. (2007). As noted by those authors, this favours a large Galactic distance on the far side of the Carina spiral arm at 6–7 kpc (similar to well-known supergiants like AG Car and HR Car, as well as the cluster Westerlund 2, etc.), rather than a closer distance that would put SBW1 inside the Carina Nebula at 2.3 kpc.
Fig. 3(a) reveals no continuum emission in the interior parts of the ring. This means that even though there is dust located there (as indicated by mid-IR thermal emission; see below), the dust does not contribute enough scattered light to affect the Hα HST/WFC3 image. We can therefore assume that gas in the ring's interior is ionized, and that the corresponding Hα emission measure provides information about the electron density in that region (the [S ii] lines are detected at too low signal-to-noise ratio in our STIS spectra to use their flux ratio as an electron-density diagnostic).
The continuum-subtracted intensity tracings in Fig. 3(b) reveal no significant difference in the Hα/[N ii] flux ratio across the ring, consistent with WFC3 imaging as noted above. This suggests that nitrogen abundances and ionization/excitation conditions are the same in the dense equatorial ring and in the more diffuse gas that fills its interior. It is therefore likely that the roughly solar values for the nebular N abundance derived by Smith et al. (2007) from ground-based spectra with lower spatial resolution apply across the entire nebula.
Spitzer IRAC and MIPS images
The position of SBW1 was observed as part of the Spitzer Space Telescope (Spitzer) survey of star formation in the Carina Nebula (P.I.: Smith; see Smith et al. 2010b) using both the Infrared Array Camera (IRAC; Fazio et al. 2004) and the Mid-Infrared Photometer for Spitzer (MIPS; Rieke et al. 2004). Fig. 4 shows individual images in the four IRAC bands, the 24 μm MIPS image, and a colour composite of the IRAC images (note that the longer-wavelength 70 and 160 μm MIPS images were saturated). More details about the observations and data reduction are provided by Smith et al. (2010b) and Povich et al. (2011).
SBW1 is clearly detected in all five filters observed by Spitzer, and photometry is listed in Table 2. In IRAC Bands 1–3, the point-like central star dominated the total flux from SBW1, although low-level extended emission is seen in all three filters consistent with a few per cent of the total flux, and the relative contribution of the extended emission increases towards longer wavelengths. The spatial extent of this extended emission in IRAC Bands 1–3 is consistent with emission from the dense equatorial ring seen in HST images, and the emission mechanism may be a combination of polycyclic aromatic hydrocarbon (PAH) emission, free–free and scattered starlight. No spectrum of the 3–6 μm IR emission is available, however.
ut date . | Tel./instr. . | Filter . | Exp. . |
---|---|---|---|
2009 Dec. 08 | HST/WFC3 | F502N | 2 × 210 s |
2009 Dec. 08 | HST/WFC3 | F656N | 4 × 210 s |
2009 Dec. 08 | HST/WFC3 | F658N | 4 × 240 s |
2010 May 18 | HST/STIS | G750M/6581 | 2 × 2070 s |
2010 Mar. 29 | Gemini South/T-ReCS | 8.8 μm | 4 × 1570 s |
2008 Mar. 30 | Gemini South/T-ReCS | 11.7 μm | 4 × 600 s |
2009 Jun. 09 | Gemini South/T-ReCS | 11.7 μm | 4 × 1577 s |
2009 Apr. 18 | Gemini South/T-ReCS | 18 μm | 2 × 2065 s |
ut date . | Tel./instr. . | Filter . | Exp. . |
---|---|---|---|
2009 Dec. 08 | HST/WFC3 | F502N | 2 × 210 s |
2009 Dec. 08 | HST/WFC3 | F656N | 4 × 210 s |
2009 Dec. 08 | HST/WFC3 | F658N | 4 × 240 s |
2010 May 18 | HST/STIS | G750M/6581 | 2 × 2070 s |
2010 Mar. 29 | Gemini South/T-ReCS | 8.8 μm | 4 × 1570 s |
2008 Mar. 30 | Gemini South/T-ReCS | 11.7 μm | 4 × 600 s |
2009 Jun. 09 | Gemini South/T-ReCS | 11.7 μm | 4 × 1577 s |
2009 Apr. 18 | Gemini South/T-ReCS | 18 μm | 2 × 2065 s |
ut date . | Tel./instr. . | Filter . | Exp. . |
---|---|---|---|
2009 Dec. 08 | HST/WFC3 | F502N | 2 × 210 s |
2009 Dec. 08 | HST/WFC3 | F656N | 4 × 210 s |
2009 Dec. 08 | HST/WFC3 | F658N | 4 × 240 s |
2010 May 18 | HST/STIS | G750M/6581 | 2 × 2070 s |
2010 Mar. 29 | Gemini South/T-ReCS | 8.8 μm | 4 × 1570 s |
2008 Mar. 30 | Gemini South/T-ReCS | 11.7 μm | 4 × 600 s |
2009 Jun. 09 | Gemini South/T-ReCS | 11.7 μm | 4 × 1577 s |
2009 Apr. 18 | Gemini South/T-ReCS | 18 μm | 2 × 2065 s |
ut date . | Tel./instr. . | Filter . | Exp. . |
---|---|---|---|
2009 Dec. 08 | HST/WFC3 | F502N | 2 × 210 s |
2009 Dec. 08 | HST/WFC3 | F656N | 4 × 210 s |
2009 Dec. 08 | HST/WFC3 | F658N | 4 × 240 s |
2010 May 18 | HST/STIS | G750M/6581 | 2 × 2070 s |
2010 Mar. 29 | Gemini South/T-ReCS | 8.8 μm | 4 × 1570 s |
2008 Mar. 30 | Gemini South/T-ReCS | 11.7 μm | 4 × 600 s |
2009 Jun. 09 | Gemini South/T-ReCS | 11.7 μm | 4 × 1577 s |
2009 Apr. 18 | Gemini South/T-ReCS | 18 μm | 2 × 2065 s |
Tel./instr. . | Filter/λ . | Fν . | ± . |
---|---|---|---|
… . | (name or μm) . | (Jy) . | (Jy) . |
2MASS | J/1.235 μm | 0.0959 | 0.0020 |
2MASS | H/1.662 μm | 0.0866 | 0.0019 |
2MASS | K/2.159 μm | 0.0697 | 0.0014 |
WISE | 3.4 μm | 0.0369 | 0.0014 |
WISE | 4.6 μm | 0.0256 | 0.0012 |
WISE | 12 μm | 0.815 | 0.038 |
WISE | 22 μm | 6.89 | 0.32 |
MSX | A/8.28 μm | 0.3034 | 0.0149 |
MSX | C/12.13 μm | 1.328 | 0.0969 |
MSX | D/14.65 μm | 1.710 | 0.115 |
MSX | E/21.34 μm | 6.415 | 0.398 |
AKARI | 9 μm | 0.3102 | 0.0086 |
AKARI | 18 μm | 7.849 | 0.484 |
Spitzer/IRAC | 3.6 μm | 0.039 | 0.0015 |
Spitzer/IRAC | 4.5 μm | 0.032 | 0.0016 |
Spitzer/IRAC | 5.8 μm | 0.027 | 0.0032 |
Spitzer/IRAC | 8.0 μm | 0.209 | 0.0229 |
Spitzer/MIPS | 24 μm | 7.688 | 1.23 |
Tel./instr. . | Filter/λ . | Fν . | ± . |
---|---|---|---|
… . | (name or μm) . | (Jy) . | (Jy) . |
2MASS | J/1.235 μm | 0.0959 | 0.0020 |
2MASS | H/1.662 μm | 0.0866 | 0.0019 |
2MASS | K/2.159 μm | 0.0697 | 0.0014 |
WISE | 3.4 μm | 0.0369 | 0.0014 |
WISE | 4.6 μm | 0.0256 | 0.0012 |
WISE | 12 μm | 0.815 | 0.038 |
WISE | 22 μm | 6.89 | 0.32 |
MSX | A/8.28 μm | 0.3034 | 0.0149 |
MSX | C/12.13 μm | 1.328 | 0.0969 |
MSX | D/14.65 μm | 1.710 | 0.115 |
MSX | E/21.34 μm | 6.415 | 0.398 |
AKARI | 9 μm | 0.3102 | 0.0086 |
AKARI | 18 μm | 7.849 | 0.484 |
Spitzer/IRAC | 3.6 μm | 0.039 | 0.0015 |
Spitzer/IRAC | 4.5 μm | 0.032 | 0.0016 |
Spitzer/IRAC | 5.8 μm | 0.027 | 0.0032 |
Spitzer/IRAC | 8.0 μm | 0.209 | 0.0229 |
Spitzer/MIPS | 24 μm | 7.688 | 1.23 |
Tel./instr. . | Filter/λ . | Fν . | ± . |
---|---|---|---|
… . | (name or μm) . | (Jy) . | (Jy) . |
2MASS | J/1.235 μm | 0.0959 | 0.0020 |
2MASS | H/1.662 μm | 0.0866 | 0.0019 |
2MASS | K/2.159 μm | 0.0697 | 0.0014 |
WISE | 3.4 μm | 0.0369 | 0.0014 |
WISE | 4.6 μm | 0.0256 | 0.0012 |
WISE | 12 μm | 0.815 | 0.038 |
WISE | 22 μm | 6.89 | 0.32 |
MSX | A/8.28 μm | 0.3034 | 0.0149 |
MSX | C/12.13 μm | 1.328 | 0.0969 |
MSX | D/14.65 μm | 1.710 | 0.115 |
MSX | E/21.34 μm | 6.415 | 0.398 |
AKARI | 9 μm | 0.3102 | 0.0086 |
AKARI | 18 μm | 7.849 | 0.484 |
Spitzer/IRAC | 3.6 μm | 0.039 | 0.0015 |
Spitzer/IRAC | 4.5 μm | 0.032 | 0.0016 |
Spitzer/IRAC | 5.8 μm | 0.027 | 0.0032 |
Spitzer/IRAC | 8.0 μm | 0.209 | 0.0229 |
Spitzer/MIPS | 24 μm | 7.688 | 1.23 |
Tel./instr. . | Filter/λ . | Fν . | ± . |
---|---|---|---|
… . | (name or μm) . | (Jy) . | (Jy) . |
2MASS | J/1.235 μm | 0.0959 | 0.0020 |
2MASS | H/1.662 μm | 0.0866 | 0.0019 |
2MASS | K/2.159 μm | 0.0697 | 0.0014 |
WISE | 3.4 μm | 0.0369 | 0.0014 |
WISE | 4.6 μm | 0.0256 | 0.0012 |
WISE | 12 μm | 0.815 | 0.038 |
WISE | 22 μm | 6.89 | 0.32 |
MSX | A/8.28 μm | 0.3034 | 0.0149 |
MSX | C/12.13 μm | 1.328 | 0.0969 |
MSX | D/14.65 μm | 1.710 | 0.115 |
MSX | E/21.34 μm | 6.415 | 0.398 |
AKARI | 9 μm | 0.3102 | 0.0086 |
AKARI | 18 μm | 7.849 | 0.484 |
Spitzer/IRAC | 3.6 μm | 0.039 | 0.0015 |
Spitzer/IRAC | 4.5 μm | 0.032 | 0.0016 |
Spitzer/IRAC | 5.8 μm | 0.027 | 0.0032 |
Spitzer/IRAC | 8.0 μm | 0.209 | 0.0229 |
Spitzer/MIPS | 24 μm | 7.688 | 1.23 |
Extending to longer wavelengths in IRAC Band 4 at ∼8 μm, the emission is qualitatively different. The emission is no longer dominated by an unresolved point source, but instead appears elongated by a few arcsec in the SE/NW direction, consistent with the major axis of the nebula. While clearly not a point source, the Band 4 emission is less spatially extended than the ring emission seen in Bands 1–3, although the poorer angular resolution of Spitzer at ∼8 μm does not clearly resolve the structure of the emitting region. This ∼8 μm emission likely arises from thermal emission from hot dust that fills the interior regions of the ring. This emission is seen more clearly in the Gemini/T-ReCS images presented below. It is also evident from an analysis of the spectral energy distribution (SED) in Section 2.5 that the ∼8 μm flux is dominated by the hot-dust component and not any stellar photospheric emission, consistent with this interpretation.
The MIPS 24 μm image of SBW1 is shown in Fig. 4(e). Examining the raw image, it is evident that the source is slightly elongated along the major axis of the nebula. We have enhanced the contrast of axisymmetric structure in this image by subtracting a smoothed version of the image from the original, and the elongated nature of the source is clear from the contours drawn in Fig. 4(e). The asymmetry is present on angular scales larger than the diffraction limit of ∼7 arcsec. This suggests that much of the 24 μm flux arises from a source larger than the two peaks of hot-dust emission that dominate the ∼8 μm flux (as seen in IRAC Band 4), since those features are separated by only ∼4 arcsec. Thus, the 24 μm flux arises largely from cooler dust in the outer equatorial ring, and not from the diffuse emission interior to the ring. This conclusion is supported by our analysis of higher-resolution ground-based mid-IR images in the next section.
Gemini South/T-ReCS images
We obtained images of SBW1 at 8.8, 11.7 and 18.0 μm on 2008 March 30, 2009 April 18 and June 9, and 2010 March 29 using T-ReCS mounted on the 8-m Gemini South telescope (see Table 1). T-ReCS was the facility mid-IR imager and spectrograph at Gemini South, with a 320 × 240 pixel Si:As IBC array, a pixel scale of 0.089 arcsec, and a resulting field of view of 28.5 × 21.4 arcsec2. The observations were taken with a 15.0 arcsec east–west chop throw. Individual sky-subtracted frames were then combined to make a co-added image in each filter. Fig. 5 shows the resulting co-added T-ReCS images at 8.8, 11.7 and 18.0 μm, compared to the HST/WFC3 Hα image on the same scale (Fig. 5 a). The left-hand column displays the images in false colour, whereas the right-hand column gives the same images with contours of the Hα emission superposed in order to show the relative positions of the ionized equatorial ring and the hot inner dust traced by mid-IR emission.
A point source at the location of the central star is clearly detected at 8.8 and 11.7 μm. Both the 8.8 and 11.7 μm images were taken under non-photometric weather conditions, and the 11.7 μm images were taken over two separate epochs, so we do not use them to derive absolute photometry. However, the resulting images can be used to provide a precise measurement of the relative contribution of the central point source to the total flux in each filter. We find that the central star contributes 2.5 ± 0.2 per cent of the total 8.8 μm flux measured in a 5-arcsec-radius circular aperture, and similarly, the central star contributes 0.6 ± 0.08 per cent of the total 11.7 μm flux measured in a 7-arcsec-radius circular aperture (a larger aperture was used to measure the 11.7 μm total flux because the 11.7 μm emission is more extended, with some contribution from the main equatorial ring). The central star is not detected in the 18 μm filter. The fractions of the total flux contributed by the central point source are useful for our analysis of the SED discussed below (Section 2.5).
The T-ReCS images also provide critical information about the spatial distribution of warm dust grains in the SBW1 nebula, which is unclear from the lower-resolution imaging with Spitzer. After separating out the central star, the spatially resolved extended structure seen in the mid-IR can be understood as two spatial components whose relative contribution to the total flux changes with wavelength: (1) emission from the thin equatorial ring with a semimajor axis of ∼6 arcsec, whose emission becomes relatively stronger with increasing wavelength, and (2) diffuse emission arising from dust distributed throughout the interior of the ring, but concentrated mainly in two arcs of emission located 2–3 arcsec to the SE and to the NW of the central star, whose contribution to the total flux decreases with increasing wavelength. Both the thin equatorial ring and the inner diffuse emission can be seen in the Hα image taken with HST. Comparing this Hα emission to the T-ReCS images (right-hand column of Fig. 5), it is evident that the equatorial ring contributes no detectable emission to the 8.8 μm image. The outer equatorial ring emission can be seen as a faint halo in the 11.7 μm image, as well as the outer boundary of the (noisy) 18 μm image. This suggests that the equatorial ring contains relatively cool dust. The inner double-peaked structure associated with the more diffuse emission inside the ring dominates the 8.8 and 11.7 μm images, but is less prominent at 18 μm. This suggests that the inner dust is relatively hot.
Tracings of the intensity along the major axis of SBW1 are shown in Fig. 6, comparing the HST Hα emission to the mid-IR emission in the three filters observed with T-ReCS at Gemini South (the intensities are scaled arbitrarily for display). These tracings confirm the impression from images that the 8.8 and 11.7 μm flux is dominated by two inner peaks of emission located within ± 3 arcsec from the star. The location of the peak of this emission and the profile shape in tracings at 8.8 and 11.7 μm match, indicating that there is no strong radial temperature gradient in the dust. It is also clear from the tracings in Fig. 6 that these two peaks of hot-dust emission coincide with a subtle enhancement of Hα emission. The 18 μm emission has a flatter distribution along the major axis of the nebula, with the flux decreasing only slightly from 2 to 6 arcsec away from the star. Like the equatorial ring seen in Hα by HST, the 18 μm emission drops off abruptly at ∼ 6.5 arcsec from the star (Fig. 6), providing strong evidence that cooler dust in the ring emits a considerable fraction of the ∼20 μm emission. Tracings at all three mid-IR wavelengths show a pronounced deficit of emission within 2 arcsec of the central star, revealing an inner region devoid of dust grains.
The SED
Fig. 7 shows the SED of SBW1 from optical through IR wavelengths using data from several sources. We obtained publically available IR photometry of SBW1 from the Infrared Science Archive (IRSA1), including measurements from the Two Micron All-Sky Survey (2MASS; Skrutskie et al. 2006), and from point-source catalogues of the Midcourse Space Experiment (MSX; Price 1995), Wide-field Infrared Survey Explorer (WISE; Wright et al. 2010) and AKARI satellites. Flux densities (in Jy) from these sources are listed in Table 2, and then converted to νFν in Fig. 7. We also show V- and R-band photometry from Smith et al. (2007) in Fig. 7.
Fig. 7 and Table 2 include photometry measured from our Spitzer IRAC and MIPS surveys of the Carina Nebula (Smith et al. 2010b). Filled blue circles in Fig. 7 correspond to total fluxes of SBW1 measured in a 12-arcsec-radius circular aperture in all five bands, while blue asterisks are the stellar flux in Bands 1, 2 and 3 measured in a 3-arcsec-radius aperture with corrections for point spread function (PSF) flux outside this radius.
Finally, the long vertical red/orange bars in Fig. 7 are derived from the fraction of the total flux contributed by the central star in our 8.8 and 11.7 μm T-ReCS images of SBW1 from Gemini South. The upper end of these bars is tied to the solid-curve fit to the SED (see below), and the bottom of each then represents the expected stellar flux at each wavelength.
As shown in Fig. 7, the observed optical and IR photometry of SBW1 can be approximated reasonably well with three simple emission components: (1) a 21 000 K blackbody corresponding to the B1.5 Iab central star (dashed curve), reddened by foreground extinction of E(B − V) = 1.15 mag; (2) a greybody (with an emissivity proportional to λ−1 at long wavelengths) representing warm dust at 190 K and (3) another greybody representing cooler dust at 85 K. At a distance of 7 kpc, these stellar-photosphere, hot-dust and cool-dust components have luminosities of 50 000 L⊙, 400 L⊙ and 2500 L⊙, respectively.
The sharp increase in total flux from 6 to 8 μm, as the SED transitions from a photosphere to warm-dust emission, is indicative of a cavity with little or no dust close to the star. The 190 K ‘hot’-dust component dominates the total flux from about 6 to 15 μm, and arises from the diffuse structures filling the interior of the ring, as seen in the double-peaked feature in the 8.8 and 11.7 μm images from Gemini South (Fig. 5). The cooler 85 K dust component dominates the total flux longward of about 18 μm, and this cooler dust appears to reside primarily in the equatorial ring. We cannot rule out the presence of some cooler dust that may produce excess far-IR luminosity (and much higher mass) from presently available data.
The typical uncertainty for this type of rough estimate of the dust mass is ± 30 per cent, dominated by assumptions about grain properties, as well as uncertainties in the ranges of temperatures that can fit the SED. We have assumed astronomical silicate as the nominal dust composition, which seems the most reasonable given that SBW1 is a massive evolved star with normal CNO chemical abundances (Smith et al. 2007). Moreover, IR spectra of SN 1987A obtained with Spitzer revealed strong silicate emission features (Bouchet et al. 2006). Unfortunately, we do not have direct observational constraints on the grain composition in SBW1. Had we instead adopted the assumption of small (a = 0.2 μm) graphite grains and calculated the dust mass following the method described by Smith & Gehrz (2005), we would have derived dust masses of Md(190) = 1.1 × 10−5 M⊙ and Md(85) = 8.7 × 10−3 M⊙ for the Td = 190 K and 85 K components, respectively. The mass for the warmer component agrees within the expected uncertainty with that derived above for silicates. The cooler component is a factor of ∼2 larger for carbon grains, but as noted above, carbon grain composition seems unlikely since SBW1 is not C-enriched. In any case, we regard the value of Md(85) = 4.7 × 10−3 M⊙ as a lower limit to the possible dust mass, since our observations do not constrain the SED longward of 24 μm, making it possible that a larger mass of cooler (Td < 85 K) dust might reside in the equatorial ring, or just outside it.
To extend this dust mass to a total nebular mass requires an assumption of the gas-to-dust mass ratio, which is poorly constrained but generally taken to be about 100:1 in massive-star nebulae (see e.g. Smith et al. 2003, and references therein). We thus infer a total gas mass of the order of 0.5–1 M⊙ for the ring nebula around SBW1, derived from the mass of the coolest dust we detect. Although the gas:dust ratio could be substantially larger for the hotter 190 K dust component that is mixed with ionized gas in the interior region of the ring, this component makes a negligible contribution to the total nebular mass. Interestingly, a mass of the order of 0.5–1 M⊙ for the SBW1 ring nebula is about equal to the mass of pre-SN ejecta surrounding SN 1987A. The directly observed mass of ionized gas in the SN 1987A ring is at least 0.04 M⊙, although this is a lower limit to the nebula mass because it is derived from emission lines and corresponds only to a thin skin at the inward-facing edge of the ring that was ionized by the SN; observations of light echoes reveal a much larger mass of material in the ring and outside it of ∼1.7 M⊙ (Sugerman et al. 2005).2
The cool-dust component that dominates the λ > 20 μm emission has a total luminosity that is ∼5 per cent of the original stellar luminosity. The vertical thickness of the dust in the equatorial ring must therefore be ∼10 per cent of its radius, or 6 × 1016 cm, based on the fraction of the total luminosity that it intercepts. The much lower luminosity of the hot-dust component, despite its closer separation from the central star, means that it is very optically thin, and so its relative luminosity cannot provide a meaningful constraint on the vertical thickness of this hot-dust component.
X-ray emission
The Carina Nebula was observed as part of a large programme with the Chandra X-ray Observatory, called the Chandra Carina Complex Project (Townsley et al. 2011). There is no X-ray source listed in the resulting catalogue of X-ray sources in this survey (Broos et al. 2011) within a radius of 15 arcsec from SBW1’s position. The quoted completeness limit of this survey would suggest an absorption-corrected upper limit to the X-ray luminosity of 1030.7 erg s−2 for sources within the Carina nebula region. It is, however, difficult to place a meaningful upper limit on the intrinsic soft X-ray luminosity of SBW1, since it is located at a distance of 7 kpc, far behind the Carina nebula (D = 2.3 kpc; Smith 2006). Thus, there could be a large absorption column, and its intrinsic soft X-ray luminosity could be substantially higher than this nominal upper limit.
OBSERVED MULTI-WAVELENGTH STRUCTURE
Fig. 8 shows the relative spatial distribution of ionized gas and warm dust in colour, with the Hα image tinted as blue/green and the mid-IR emission tinted as red/orange. This image encapsulates the basic multiwavelength structure of SBW1, and is useful in the discussion below. One of the curious things about Fig. 8 is that it runs counter to normal expectations that dusty regions will be farther from the source of radiation than the ionized gas. This and other issues are clarified below.
Detailed structure of the dense equatorial ring
The new HST/WFC3 images in Hα and [N ii] provide a detailed view of the morphology in the ring nebula around SBW1. To the extent that SBW1 is a suitable analogue, these images also provide our best view of the circumstellar environment around an object like SN 1987A. The most prominent feature in the images is a thin ring of emission, presumably in the equatorial plane. The new HST images also provide a better view of the outer rings or hourglass structure in SBW1, and they clarify the nature of the diffuse emission structures in the ring's interior, closer to the central star. Below we examine the structure of the ring in detail, and then conduct a comparison with the equatorial ring of SN 1987A. Structures inside the ring are discussed in the following section.
This estimate of the mass of ionized gas is much less than the expected total mass of H gas. The ring's dust mass of 5 × 10−3 M⊙ measured from the luminosity of ∼85 K dust in the SED (see Section 2.5; Fig. 7) would imply a total gas mass of roughly 0.5 M⊙ if the equatorial ring has a normal gas:dust mass ratio of 100:1. This implies that the ionization fraction of the ring is only a few per cent. The equatorial ring is therefore ionization bounded, with the location of the ionization front determined by the flux of ionizing photons from the central star as well as the initial density structure of the ring. The incomplete ionization of the ring has a strong impact on the shape of the rest of the nebula around SBW1, and we return to the issue in Sections 4.1 and 4.2.
A strong increase in the ionizing flux would be expected to ionize a larger fraction of the total mass in the ring, and would affect the apparent-brightness structure of the ring. This may account for some of the observed morphological differences between SBW1 and the ring of SN 1987A, because the ring around SN 1987A is viewed after a huge burst of ionizing radiation from the UV flash of the SN. For example, currently the diffuse gas that fills the interior of the SBW1 ring is brighter compared to the thin ring itself, making SBW1 appear more filled-in than SN 1987A (see Fig. 9). If the ionizing flux of SBW1 were to increase suddenly, we would expect the thin equatorial ring to get much brighter, because more of the high-density neutral gas in the ring would be ionized and the emission measure (proportional to n2e) would rise, whereas the interior of the ring would not brighten because it is already fully ionized. Given the difference in ionizing flux, however, the morphologies of the two rings are already strikingly similar.
Fig. 9 shows Hα images of the equatorial rings around both SBW1 and SN 1987A, as well as plots of the deprojected intensity around the equatorial ring for each source. These were produced by first fitting an ellipse to the HST image of each ring, and then ‘deprojecting’ the appearance of the source by stretching the image by the appropriate factor along the minor axis, in order to make that ellipse into a circle. The deprojected intensity plots in Fig. 9 are then a radial tracing of the intensity at all position angles around the ring. For SN 1987A, this is meant to duplicate the analysis by Sugerman et al. (2002), which provides a basis for comparison. We then performed the same analysis for SBW1. For SN 1987A we adopted a ring inclination of i = 43| ${.\!\!\!\!\!\!^{\circ}}$|8 (Sugerman et al. 2002), and for SBW1 we measured i = 50| ${.\!\!\!\!\!\!^{\circ}}$|6 ± 0| ${.\!\!\!\!\!\!^{\circ}}$|6, in agreement with an earlier estimate of i = 50| ${.\!\!\!\!\!\!^{\circ}}$|2 ± 1° from a ground-based image (Smith et al. 2007).
Figs 9(c) and (d) are surprisingly similar. If we characterize the ring structure as dominated by a number of beads along a string, we see that both the wiggles in the string (Δr/R) and rough angular scale (Δl/R) of the beads are similar in the two objects, where R is the average ring radius (0.21 pc for SN 1987A and 0.19 pc for SBW1; see Smith et al. 2007). The small variations in Δr are typically about ±4 per cent of R for SBW1 and ±6 per cent of R for SN 1987A, with a few features reaching inward to 85 per cent of R for both objects. The angular scale of clumps appears somewhat larger for SN 1987A in these figures, but this depends on the effective spatial resolution, which is slightly higher for SBW1 because it is closer (we have smoothed the SBW1 image by about a factor of 3, but it is a factor of 7 closer.). A typical separation between major clumps in the ring is 8° to 15° for both, with roughly 22 to 25 major clumps around each ring. Note that these figures are intended to exaggerate the differences in structure between the two rings, with only a portion of R being plotted.
The main interesting result of this analysis is that if we allow for some small differences in spatial resolution and the precise arrangement of blobs, the equatorial ring around SBW1 is for all practical purposes indistinguishable from that around SN 1987A, in terms of its basic physical attributes. Yet, SBW1 is a factor of 7 closer to us, and so in the original HST image we have seven times better effective spatial resolution and can see more details of the structure than for SN 1987A. The chief reason this is interesting is that SBW1 does not exhibit the long RT fingers envisioned for SN 1987A in order to explain the occurrence of hotspots when the tips of the RT fingers were hit by the blast wave. Instead, SBW1 appears to be a fragmented, clumpy ring with a small radial extent. Since RT fingers are expected to mark the contact discontinuity between the faster BSG wind and the slower and flattened RSG disc wind, the lack of such features carries implications for the formation mechanism of the ring: the ring is probably not formed by a fast BSG wind sweeping into an extended disc. Indeed, after analysing the structures interior to the ring in the next section, we conclude that the expanding BSG wind is not directly interacting with the dense equatorial ring, due to the increased pressure resulting from the photoionization of the ring. The clumpy dense ring must therefore have had a different origin, as we discuss later.
Finally, we note a possibly interesting feature of the ring morphology. The brightest part of the ring on the north-west edge, at position angle PA ≈ 310°, coincides with a pronounced kink in the ring, exhibiting the largest jump in deviations from the average ring radius (see Fig. 9c). The significance or cause of this structure is not obvious. However, a similar kink is seen to be associated with the brightest portions of the ring around SN 1987A as well, coincidentally also located at PA ≈ 310°. One can readily imagine how a single large departure from azimuthal symmetry could arise in a sudden ejection in a binary system, but this is harder to explain when the structure of the ring arises from colliding winds alone (this is discussed in more detail below).
Structures inside the ring
The most unexpected result of this study of SBW1 concerns the diffuse emission from structures in the interior of the equatorial ring. While ground-based Hα images revealed some diffuse emission apparently filling the ring interior (Smith et al. 2007), it was difficult to reliably separate contributions from the ring and the stellar PSF. The new HST/WFC3 F656N image shows that the distribution is not uniform across the ring; there are two regions of enhanced Hα emission located about 1 to 3 arcsec on either side of the central star, and there is a marked deficit of Hα emission within 1 to 2 arcsec of the star (see Figs 1 and 5).
More striking, however, are the high-resolution mid-IR images obtained with T-ReCS on Gemini South. In the mid-IR, we see two distinct peaks of warm-dust emission located 1 to 3 arcsec on either side of the central star, oriented along the minor axis of the ring (see Figs 5 and 8). These dust-emission peaks appear to be spatially coincident with the enhancements of Hα mentioned above, but have a much stronger contrast in the mid-IR. The two peaks are probably caused by limb brightening of a toroidal structure with a projected radius of ∼2 arcsec, and if they are located in the equatorial plane, these inner dust peaks occur at about 1/3 of the radius of the dense equatorial ring (which is ∼0.2 pc; Smith et al. 2007). The IR emission also shows a clear deficit within a few arcsec of the star, indicating an inner region that is relatively devoid of hot dust.
The double-peaked mid-IR emission has the strongest contrast at the shorter mid-IR wavelengths (8.8 and 11.3 μm), whereas the emission structure appears more uniform at the longest thermal-IR wavelength (18 μm) due to a rising contribution from the outer equatorial ring. This suggests that the double-peaked features have the hottest dust, consistent with their location closer to the star than the other structures. The analysis of the optical/IR SED of SBW1 (Fig. 7; Section 2.5) suggests the presence of two dust components emitting at ∼190 K and 85 K. The mid-IR images therefore indicate that the warmer 190 K component must be associated with the inner double-peaked emission, and that the cooler 85 K component is associated with cooler dust located in the dense equatorial ring.
In any case, the exact heating mechanism of the dust interior to the ring is less critical than the fact that it resides there. The existence of dust at this location is problematic for conventional interacting-wind models for the formation of a ring nebula like the one around SN 1987A (Blondin & Lundqvist 1993; Martin & Arnett 1995; Collins et al. 1999; Morris & Podsiadlowski 2009). These models predict that the RT instabilities in the ring itself mark the contact discontinuity between the BSG wind and a pre-existing disc wind. If so, then the volume interior to the ring should be filled with mass from the BSG wind. There will be a reverse shock somewhere between the star and the ring, but this is not expected to form dust. Steady BSG winds do not form dust on their own, and the SED and spectrum of the central source in SBW1 reveal no evidence for a dusty red giant star in a binary system. Thus, the dust in the interior of the ring around SBW1 must have some other origin, and suggests that a simple interacting-winds model incorporating only hydrodynamic effects must be rejected.
Implications for a possible origin of this structure are discussed in the next section, related to the proposed existence of an ionized portion of a progenitor RSG wind around SN 1987A (Chevalier & Dwarkadas 1995). In short, we propose that the dusty region interior to the SBW1 ring originates from photoionization and photoevaporation of neutral dusty gas in the dense equatorial ring, which then expands to partially fill the interior region of the ring. This dusty photoevaporative flow (DPF) meets the expanding BSG wind at a shock interface with a toroidal geometry inside the ring. A sketch of the proposed geometry for SBW1 is given in Fig. 10, and is discussed in more detail below.
DISCUSSION
Origin of the dusty H ii region inside the ring: the role of ionization and photoevaporation
In the previous sections, we have described how observed structural properties of the SBW1 ring seem to be in conflict with some aspects of hydrodynamical models for the formation of SN 1987A's nebula. Specifically, we highlighted the lack of long RT fingers and the existence of dust interior to the ring. Both challenge the standard interacting-winds picture where RT fingers form as a result of instabilities in the contact discontinuity between a fast BSG wind and a flattened disc-like RSG wind. The presence of dust far interior to the ring is extremely important, because it means that the BSG wind (presumably devoid of dust) has not yet reached the ring, and in fact has not even penetrated past about 1/3 the radius of the ring in the equatorial plane.
In this section we suggest a modification to this scenario, wherein ionization and evaporation of a dense neutral ring produce a photoevaporative flow which, in turn, profoundly influences the overall hydrodynamics and shaping of the nebula. These ionization effects have so far been neglected in numerical simulations, but we argue that they are essential, and our intention is that our discussion of the observations will inspire such a numerical study. Our proposed scenario follows the suggestion by Chevalier & Dwarkadas (1995), who invoked a dense H ii region arising from photoionized portions of a previous RSG wind in order to account for some aspects of the advancing blast wave of SN 1987A. The additions and modifications that we make to this scenario concern the dynamical influence of the ionized photoevaporative flow off the ring, the nature of the interaction between this flow and the BSG wind, and the origin of dust in the system (not necessarily from a RSG wind).
In Fig. 10 we show a sketch of the proposed geometry for the SBW1 ring and material inside it (this is a side view in the equatorial plane). In brief, the picture we propose is that the warm dust seen in Gemini South/T-ReCS images resides at or near the location of a shock front at the interface where the outflowing BSG wind meets the inward flow of dusty ionized gas that has been photoevaporated from the dense equatorial ring/torus. In a two-dimensional cross-section, this structure has a parabolic shape due to the divergent flow from the ring. This is analogous to the opening angle of two colliding winds in a binary system, except here the structure is toroidal rather than a cone, because the source of the flow is a dense ring rather than a companion star. In effect, this ionized photoevaporative flow off the dense equatorial ring is able to keep the BSG wind at bay in the equatorial plane, preventing it from reaching the equatorial ring so that no direct hydrodynamical interaction between the ring and the BSG wind is able to occur. Instead, ionizing photons that travel ahead of the shock have the most important dynamical effect because of the increase in pressure when the gas is ionized. The geometry sketched in Fig. 10 could arise from a sequence of events such as that described in the following section.
Formation of the observed structure in SBW1
In this section, we outline a possible sequence of events that may have led to the formation of the ring nebula around SBW1. We also discuss the most relevant physical mechanisms for each component of the nebula. The key difference between this and previously suggested scenarios is the hydrodynamic role of an ionized photoevaporative flow off the ring. A key point is that the ionized gas pressure is allowed to have a major impact on the hydrodynamics because of the slow (10–20 km s−1) initial expansion speed of the dense ring.
Episodic ring ejection as a BSG, not a RSG disc-wind
Several previous studies that sought to explain the unusual structure of the nebula around SN 1987A assumed some form of a fast BSG wind (300–400 km s−1) blowing into a slower RSG wind (5–10 km s−1). Different versions of this are seen in the more standard ‘hydrodynamics-only’ interacting-winds model (Luo & McCray 1991; Blondin & Lundqvist 1993; Martin & Arnett 1995), as well as in the H ii region model of Chevalier & Dwarkadas (1995). In the sections to follow, we argue that the structures observed around SBW1 seem to clearly support the model of Chevalier & Dwarkadas (1995), wherein the effects of photoionization are critical in producing a dense H ii region at low latitudes, but allowing the BSG wind to expand freely out the poles. However, one aspect in which our suggested scenario differs from that of Chevalier & Dwarkadas (1995) is in the assumption of a previous RSG wind, as discussed in this section.
The assumption of a previous RSG wind is motivated by the slow (∼10 km s−1) expansion speed in the equatorial ring of SN 1987A. However, this assumption also introduces significant difficulties, as the RSG wind must be concentrated in a thin disc with a very high equator-to-pole density contrast (Luo & McCray 1991; Blondin & Lundqvist 1993; Martin & Arnett 1995). Moreover, these hydrodynamic simulations predict an hourglass structure with complete polar caps, rather than empty polar rings. Since the extended envelope of any single RSG star will have a very low angular velocity, one must invoke some interaction with a companion star in a binary system (possibly even a merger) in any scenario involving a RSG (Collins et al. 1999; Morris & Podsiadlowski 2009).
Yet, if one must invoke a binary system, then a previous RSG phase is not really needed, because a RSG wind is not the only way to form a slowly expanding dusty equatorial ring, and, moreover, such a flattened disc-like outflow has never been observed around a RSG (binary or single). A different case that has well-established observational and theoretical precedent is the episodic ejection of a ring by a BSG. This could have been the result of a binary merger event (Morris & Podsiadlowski 2009), an episodic ejection event occurring during Roche lobe overflow (RLOF) as seen in the eclipsing binary RY Scuti (Smith et al. 2011), or an eruptive ejection from a rotating star (Smith & Townsend 2007; Shacham & Shaviv 2012). In either a binary merger event or accretion from a companion in RLOF, a large amount of angular momentum must be shed in order for material to be accreted on to the merger product or the mass-gaining star, respectively, and one thus expects the mechanical mass shedding to occur in the equator. If the mass loss is episodic, this brief equatorial ejection will form a ring rather than an extended disc because of the brief nature of the mass ejection. Shedding of mass at the equator in this way forms a nearly Keplerian ring with a relatively slow radial expansion speed, which can be much less than the escape speed from the surface of a BSG star.
A prime example of the scenario we suggest is seen in the Galactic eclipsing binary system RY Scuti, in which an OB supergiant binary system undergoing RLOF has recently ejected a ring/torus that is spatially resolved in HST images (Smith et al. 1999, 2002, 2011). The radial expansion speed of the equatorial ring/torus in RY Scuti is only 30–40 km s−1 (Smith et al. 2002), much less than the several hundred km s−1 escape velocity or wind speed from the O9/B0 supergiant stars. Moreover, this has occurred without a RSG phase, because the stars are too close together with an orbital period of only ∼11 d (Smith et al. 2002; Grundstrom et al. 2007, and references therein). RY Scuti is an important example, because measured proper motions of the very young (120 yr old) nebula show that the equatorial mass loss was the result of a pair of brief episodic ejections rather than a steady flow of mass from the equator (Smith et al. 2011).
Moreover, the nebula around RY Scuti reveals that the slow equatorial outflow led to the formation of significant amounts of dust (Gehrz et al. 2001), even though it is a binary composed of two hot OB supergiants. In fact, episodic or eruptive mass loss like that seen in luminous blue variables (LBVs) may be an essential ingredient for the formation of dust around BSG stars, as discussed in detail recently by Kochanek (2011). The class of B[e] supergiant stars provides many additional examples of BSG stars that have formed significant quantities of dust in a slowly expanding equatorial ring or torus (e.g. Zickgraf et al. 1996). Thus, the formation of a slow and dusty equatorial ring does not require a previous RSG phase. In fact, of the three known Galactic analogues of the nebula around SN 1987A (these are HD 168625, Sher 25 and SBW1), two (Sher 25 and SBW1) have chemical abundances that are inconsistent with a previous RSG phase (Smartt et al. 2002; Smith et al. 2007). The third analogue, HD 168526, is an LBV that has quite likely suffered previous eruptive mass loss (Smith 2007).
What about the observed clumpy structure of the SBW1 ring? As noted above, the dust residing far interior to the ring suggests that the BSG wind has not yet reached the radius of the dense equatorial ring, and as such, the clumps observed in the equatorial ring cannot be the result of RT instabilities at the contact discontinuity as required in interacting-winds models. We suggest instead that this regular clumpy structure may occur naturally due to thermal instability and fragmentation of an ejected ring, rather than from RT instabilities. In both SBW1 and SN 1987A, the measured ring expansion speed is 10–20 km s−1, so the initially ionized ring would be expected to fragment as it cools and recombines (and its sound speed drops well below 10 km s−1). Perhaps this fragmentation could lead to the fairly regular series of clumps around the ring in both objects (Fig. 9). This is, however, a minor point.
Subsequently, as the BSG star now seen at the centre of SBW1 recovers from the mass-loss event and presumably contracts to become a hotter BSG, its wind speed (escape speed) and its ionizing flux could both increase, and their ratio has an important influence on the resulting structure. The influence of photoionization of the ring ejecta is discussed next.
Photoionization of the ring, and a dusty photoevaporative flow
After ejecting a thin equatorial ring, the inner edge of the ring will be exposed to ionizing photons from the central BSG star. Chevalier & Dwarkadas (1995) predicted that photoionization had a strong impact on the structure interior to the SN 1987A ring, and we suggest a similar process in the case of the observed structure around SBW1.
The originally neutral gas in the dense equatorial ring is mixed with dust. This cool dust resides outside the ionization front in the ring, and is observed as the ∼85 K component in the IR SED of SBW1 (Fig. 7), which is spatially resolved to be coincident with the thin ring in our 18 μm Gemini/T-ReCS image (Fig. 5). After being struck by ionizing photons, the ionized photoevaporated material from the ring will likely entrain the dust with it as it expands into the interior of the ring. Thus, we refer to the resulting expansion as a DPF in Fig. 10. This dust in the DPF moves inward and eventually piles up at the shock front (discussed in the following section), giving rise to hot (∼190 K) dust emission peaks inside the ring that are spatially resolved in Gemini images (Figs 5 and 8). This hypothesis seems far more likely than any alternative explanation for the origin of dust seen interior to the ring in mid-IR images. Such alternatives would require that dust forms quickly in the fast BSG wind and survives passage through the strong (several 102 km s−1) reverse shock.
This physical situation we propose for SBW1 has an interesting analogue in the evaporating protoplanetary discs (the so-called ‘proplyds’) seen in HST images of the Orion nebula (see Bally et al. 1998; Johnstone, Hollenbach & Bally 1998). In these objects, UV radiation from the nearby O6 V star θ1C Ori causes a photoevaporative flow off the dense and dusty protoplanetary disc envelopes that are seen in silhouette. After passing through an ionization front, the ionized photoevaporative flow continues to expand in a divergent flow until it collides directly with the fast stellar wind of θ1C Ori. The proplyds located in the Trapezium very close to θ1C Ori show bright arcs of Hα emission marking the shock between the ionized photoevaporative flow and the stellar wind (Bally et al. 1998). Interestingly, these same Hα arcs are also very bright in the thermal-IR continuum in high-resolution 11.7 μm images (Smith et al. 2005), indicating the presence of hot dust. Since the wind of θ1C Ori is too fast, rarefied and hot to form dust on its own, the dust in these shocks must have been entrained in the ionized photoevaporative flow from the proplyds. This is akin to the situation we propose for SBW1, except that instead of an externally ionized disc envelope in the Orion proplyds, we have a thin ring illuminated from the inside.
Because of the thin-ring geometry of the neutral gas reservoir, the DPF will expand into the interior of the ring toward the star, but the pressure of the ionized gas will also cause it to expand over a range of latitudes above and below the equator. This will form a thicker torus-like geometry stretching from the equatorial plane to mid-latitudes, as depicted by the DPF in Fig. 10.
Collision of the DPF and the BSG wind
The DPF will expand away from the dense equatorial ring, primarily into its interior, until it collides with the outflow from the stellar wind of the central BSG. This collision between the BSG wind and the DPF will have three boundaries: the reverse shock that decelerates the BSG wind, the forward shock driven into the DPF and the contact discontinuity between them. These are labelled in Fig. 10. The zone between the forward shock and contact discontinuity is drawn as a thin region due to the high-density gas and slow velocity of the forward shock. The region between the contact discontinuity and reverse shock is thicker due to the higher speed and lower density of the BSG wind.
Since the H ii region is toroidal and caused by evaporation of a thin equatorial ring, we expect the density to drop at latitudes above and below the equator because the ionized photoevaporative flow diverges. If so, this would cause the location of the shock front to curve outward to larger radii, as depicted in Fig. 10. The variation with latitude is discussed further below.
In the sketch of the geometry shown in Fig. 10, mass in the outflowing BSG wind and gas from the ionized photoevaporative flow will both pile up at the shock front. After colliding, the post-shock gas should then flow away from the star in a non-radial trajectory that follows the geometry of the shock interface (small black arrows in Fig. 10). This is analogous to flow down the shock cone in a colliding-wind binary, but here the geometry is toroidal rather than a cone. The material in this post-shock flow is denser than gas in the H ii region or the BSG wind. It is this higher-density zone of ionized gas that probably dominates Hα emission because of the n2e dependence of recombination emission.
In fact, it seems likely that the dense gas in this post-shock flow has already been observed directly. There is evidence for the resulting non-radial flow in high-resolution, long-slit echelle spectra of Hα and [N ii] emission in SBW1 that we presented in a previous study (Smith et al. 2007), obtained with the EMMI spectrograph on the New Technology Telescope (NTT). In that paper, we were perplexed by the unusual position/velocity structures, but it now makes sense in the context of the geometry proposed above. In Fig. 11 we show the position–velocity diagram of the earlier EMMI echelle spectra of [N ii] λ6583, presented originally by Smith et al. (2007), as well as a sketch for how these velocity structures may arise from the non-radial flow at the shock interface. While a quantitative test of this idea will require numerical simulations of the hydrodynamics, the sketch of the non-radial flow trajectories in Fig. 11(c) provides a reasonable explanation for the observed position–velocity structure in Fig. 11(b).
Latitudinal structure
In both SBW1 and around SN 1987A, HST images do not reveal dense nebular material in the polar directions. Rather, the polar rings of SN 1987A are true hollow rings residing at mid-latitudes (Burrows et al. 1995; Plait et al. 1995). The same appears to be true based on our HST images of SBW1. The lack of any dense polar structures in SBW1 or in SN 1987A argues against a shaping mechanism like the one discussed by Morris & Podsiadlowski (2009), because that formation model predicts dense polar caps in the nebula. The same is true for most of the colliding-wind models that predict complete hourglass shapes (see above), not empty rings at mid-latitudes. If such polar caps existed, they would be impacted by the fast polar wind of the BSG, and would likely be swept into a thin dense shell, and this interface would therefore be easy to detect in deep images.
Instead, the empty polar regions favour a different scenario like that suggested by Chevalier & Dwarkadas (1995), which is similar to the picture we advocate here. Rather than a pre-existing disc wind, we have proposed that the central BSG star episodically ejects a dense ring, which is then photoionized as it expands. The ionized gas expands to form a thick torus around the ring, which provides a barrier for the BSG wind. The collision between the wind and the ionized torus (H ii region) creates a curved shock front. In this picture, all of the slow and dense material is confined to low latitudes, so that the BSG wind is able to expand with no obstruction out the poles.
Is the BSG wind spherical? Ejection of an equatorial ring requires a system with excess angular momentum, suggesting either a rapidly rotating star or a close binary system. Even in a binary system with RLOF, the mass-gainer star accretes both mass and angular momentum, and would be expected to be in very rapid rotation. The geometry of the stellar wind from a rapidly rotating hot star is complicated, because hot stars have line-driven winds. Rapidly rotating stars suffer gravity darkening (von Zeipel 1924), producing cooler equatorial zones and hotter poles. Such stars also have higher escape speeds out the poles than at the equator, so the polar winds are faster. The oblate structure of the star produces a larger cross-section in the polar directions, and also the hotter poles produce much more radiative flux (Frad ∝ T4). This is critical for a radiation-driven wind. This stronger radiative flux at the poles, combined with velocity-dependent forces, means that rapidly rotating stars have higher mass flux and higher wind speed out the poles (see Owocki, Cranmer & Blondin 1994; Owocki, Cranmer & Gayley 1996; Owocki & Gayley 1997; Owocki, Gayley & Cranmer 1998). This is for radiatively driven mass loss – different from the more intuitive case of mechanical mass loss from a spinning object, which should produce an equatorial flow (appropriate for RLOF).
Including a stronger polar wind for the present state of the central BSG would not change the qualitative picture in Fig. 10, because a stronger wind out the poles would still escape unobstructed. It may, however, make it easier to produce this shape and to reconcile the required values of |$\dot{M}$| with those of normal BSG stars. Note that for a B1.5 Iab star, the expected value of |$\dot{M}$| is ∼5 × 10−7 M⊙ yr−1, or perhaps a small factor less if we allow for clumping. This would be in rough agreement with our upper limit for SBW1 derived from IR observations, ∼3 × 10−7 M⊙ yr−1. However, for both SBW1 and SN 1987A, densities at the equator would seem to suggest a slightly lower value of |$\dot{M}$|, closer to 10−7 M⊙ yr−1. Moreover, our assumed wind speed of 300 km s−1 (the same value assumed by Chevalier & Dwarkadas 1995) is somewhat lower than the speeds we expect for a B1.5 supergiant, closer to 500–600 km s−1. Perhaps the factor of 2–3 discrepancy in |$\dot{M}$| and VW is not too concerning, but allowing the central BSG to have a bipolar wind – with higher |$\dot{M}$| and VW out the poles, and lower |$\dot{M}$| and VW at low latitudes near the equator – would make the agreement better.
A bipolar wind from the BSG might also help explain another observed peculiarity of SN 1987A. Studies of the evolution of light echoes from SN 1987A have revealed a large bipolar shell, much larger than the triple-ring system seen in HST images (Wampler et al. 1990; Crotts, Kunkel & Heathcote 1995; Sugerman et al. 2005). In the geometry advocated here, there are no polar caps in the inner nebula, so the strong BSG wind is free to escape out the poles to much larger radii. Doing so, a bipolar wind could sweep into the surrounding interstellar medium and create a much larger bipolar shell.
Origin of the triple rings around SN 1987A
We have noted multiple times that while interacting-winds models can reproduce structures resembling the equatorial rings around SN 1987A and SBW1 by forming an hourglass shape (Blondin & Lundqvist 1993; Martin & Arnett 1995; Collins et al. 1999), the polar rings of SN 1987A have actually been an enduring problem with no satisfactory explanation. The nebula around SBW1 does not exhibit polar rings that are as prominent as those around SN 1987A, but the detailed geometry we see in SBW1 may provide some intriguing clues about how the polar rings might have been formed.
The key difference between previously published interacting-winds models and the model advocated here is the existence of a dense DPF off the compact and dense equatorial ring (Fig. 10). We have argued that its existence in SBW1 is supported by the presence of dust interior to the thin Hα ring (Fig. 8), as well as the unusual kinematic structure of the emission lines in long-slit echelle spectra (Fig. 11). Chevalier & Dwarkadas (1995) argued for the existence of a similar structure around SN 1987A based on the properties of the advancing blast wave. This DPF creates a stand-off shock front (reverse shock) in the equator at RRS = 0.05 pc, but at latitudes above and below the equator the interface will take a parabolic shape in cross-section (Fig. 10).
How does the existence of this DPF and reverse shock modify the situation? In this scenario, one would still expect the formation of a thin hourglass-shaped nebula as the H ii region tries to expand into any neutral gas outside it (see Chevalier & Dwarkadas 1995). The fast BSG wind will escape in the polar direction and will sweep away any of this outer nebulosity, but it will be prevented from doing so at low latitudes. In effect, the DPF and reverse shock keep the BSG stellar wind at bay, and thereby shield the nebula at low latitudes. Therefore, we expect that the thin outer hourglass structure will remain at mid- and low-latitudes, as depicted in Fig. 10.
Now, recall that the collision between the DPF and the BSG stellar wind leads to a non-radial flow of dense gas down the shock front (Fig. 11). An interesting consequence of the proposed flow down the shock interface is that the densest gas will flow outward at a narrow range of latitudes above and below the equator. This flow may eventually intersect the hourglass structure, which in this picture is the outer boundary of the photoionized cavity – labelled as the ionization shock front in Fig. 10. Interestingly, the geometrical intersection of these two curved surfaces in three dimensions is a pair of plane-parallel rings above and below the equator. (This would still be true if the outer ionization shock front were an hourglass, a sphere or a cylinder.) If there is a density enhancement associated with this intersection, it may help explain the origin of the polar rings around SN 1987A and HD 168625.
Alternatively, a similar idea was suggested by Chiţâ et al. (2008) based on their hydrodynamic simulations of a bipolar BSG wind interacting with a thin spherical shell formed by the terminal shock of the previous RSG phase. They assumed that as the star evolves in a blue loop after the RSG phase, the star must spin up, and that the wind from the BSG will have a bipolar shape for the reasons noted in Section 4.2.4. These authors showed that the collision of a bipolar wind with a thin shell could lead to a pair of polar rings. In this scenario, the existence of the reverse shock could help to stabilize the latitude range where the bipolar BSG wind and thin shell intersect, so that the rings would not be as transient as in those simulations.
We hope that the observations and discussion presented here will motivate new simulations that include the effects of a DPF. This basic picture also allows for some differences among objects, which may help explain some of the diversity in the structures seen around SN 1987A, SBW1, Sher 25 and HD 168625. For SBW1 and SN 1987A, we expect that the initial structures of the two nebulae were very similar, but that in SN 1987A, the structure has been modified: the dense ring has become more fully ionized and the interior dust was probably vaporized by a UV flash from SN shock breakout. Even ignoring the SN UV flash, however, certain aspects of the structure can be adjusted. For example, other central stars may have different values of QH or |$\dot{M}$|, either of which may change the relative location of the shock between the BSG and the photoevaporative flow, and may affect the opening angle. This could have important implications for understanding when the SN 1987A blast wave encountered the H ii region (see Chevalier & Dwarkadas 1995; Gaensler et al. 2000). The stellar luminosity values for SBW1 and Sk−69°202 are very similar, so their mass-loss rates should be similar, but SBW1 is hotter (B1.5 instead of B3) and has a higher value of QH. This may strengthen the relative importance of the DPF somewhat in SBW1. Nevertheless, SBW1 and SN 1987A may trace one another more closely than in other BSG ring nebulae, such as Sher 25 and HD 168625, both of which have substantially higher stellar luminosities. Nevertheless, one could probably tune the basic picture presented here to be applicable to those two objects as well.
SUMMARY
In this paper, we have presented a series of multiwavelength images and spectra of the nebula around the BSG SBW1, which is seen in the Carina nebula. The observations included the first images and spectra of this target obtained with HST using WFC3 and STIS, as well as IR images from Spitzer and Gemini South/T-ReCS. The analysis of this multiwavelength data set has led to several conclusions about SBW1 and related systems, enumerated below.
HST/WFC3-UVIS F656N and F658N images of SBW1 reveal a thin, clumpy ring, but do not show the long RT fingers that have been invoked to explain the origin of the hot spots around SN 1987A when the tips of the RT fingers are hit by the blast wave. Instead, the structure is more akin to beads along a string with little dense matter in between. Instead of arising when a wind swept into an extended thin equatorial disc, as suggested by hydrodynamical modelling of SN 1987A, the observed structure may suggest an episodic ejection of a ring that fragmented into clumps as it cooled.
We noted that the episodic ejection of a slow, clumpy ring could result during a brief phase of RLOF, as seen in the RY Scuti eclipsing binary system. A previous RSG phase is not needed to explain the slow expansion speed of the ring.
The radius of the ring, size and angular scale of clumps in the ring, and deviations from the average radius (wiggles) seen in SBW1 are all an excellent match to the structures seen in SN 1987A's ring.
The radial velocity measured in our STIS spectra as well as the luminosity and extinction indicated by the SED both confirm the large distance of ∼7 kpc (Smith et al. 2007), meaning that SBW1 is far behind the Carina Nebula and seen in projection, rather than actually being inside it.
High-resolution images reveal complex diffuse emission filling the interior of the ring. Hα emission seen by HST uniformly fills much of the ring interior, with two enhancements in the Hα flux about 1–2 arcsec away from the star on both sides along the major axis. Much more importantly, mid-IR images obtained with Gemini South reveal strong dust continuum emission associated with these same positions. The warm dust emission residing inside the ring is critically important, because the dust cannot arise from the BSG wind.
The IR SED reveals two distinct dust–temperature components: T = 190 K dust with a mass of only 10−5 M⊙, and ∼85 K dust with a larger mass of 5 × 10−3 M⊙. This implies a total gas mass for the SBW1 ring of at least 0.5 M⊙. High-resolution mid-IR images from Gemini South confirm that the cooler 85 K component arises from dust in the dense equatorial ring, whereas the warmer 190 K dust comes from the double-peaked diffuse emission in the ring's interior.
Both Hα and mid-IR images reveal a clear deficit of diffuse emission within 1–2 arcsec of the central star, suggesting that a cavity has been cleared by the BSG wind.
We propose a model for the origin of the observed structures in SBW1 that departs markedly from the standard interacting-winds model for SN 1987A. Instead of the equatorial ring arising as a consequence of a fast wind sweeping into a slow extended disc wind, we propose that the ring was ejected in an episodic event and then photoionized by the BSG. When the dense neutral clumps in the ring get ionized, they produce an ionized photoevaporative flow that has an important hydrodynamic effect because of the pressure of the ionized gas. The double-peaked diffuse Hα and mid-IR emission inside the ring marks the location of a shock where the BSG wind collides with the DPF. This DPF structure is analogous to the ‘H ii region’ proposed by Chevalier & Dwarkadas (1995) to account for the expansion of the blast wave of SN 1987A.
In this model, the DPF keeps the BSG wind at bay in the equatorial plane, thereby shielding the ring and other nebular structures at low latitudes. The BSG can expand out the poles unobstructed.
When the DPF and BSG wind collide in the equator, they must flow downstream. The shock front probably takes on a curved shape, diverting the flow to a narrow range of mid-latitudes (see Figs 10 and 11).
We note some implications for the origin of structures around SN 1987A. If the BSG wind expands freely out the poles (i.e. there are no polar caps in the nebula), it may help explain the existence of a much larger bipolar shell seen with light echoes (Sugerman et al. 2005). Moreover, the non-radial post-shock flow that is confined to a narrow range of latitudes (point 10 above) may help explain the existence of the polar rings around SN 1987A when this flow intersects the dense shell at the outer boundary of the H ii region. We encourage further investigation of this idea using numerical simulations that account for the effects of a photoionized flow. As we note in the text, this photoionized flow is permitted to have an important dynamical effect because the initial expansion speed of the neutral dense ring is slow, comparable to the sound speed in the ionized gas.
Support was provided by the National Aeronautics and Space Administration (NASA) through HST grants GO-11637, GO-11977 and AR-12623 from the Space Telescope Science Institute, which is operated by AURA, Inc., under NASA contract NAS5-26555. Additional support for this work was provided by NASA through awards issued by JPL/Caltech as part of Spitzer programs GO-3420, GO-20452 and GO-30848. AVF is also grateful for financial assistance from NSF grant AST-1211916. This publication makes use of data products from 2MASS, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, funded by NASA and the National Science Foundation (NSF). We have also used data products from the WISE, which is a joint project of the University of California, Los Angeles and the Jet Propulsion Laboratory/California Institute of Technology, funded by NASA.
Based in part on observations made with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc. (AURA), under NASA contract NAS5-26555.
Based in part on observations obtained at the Gemini Observatory, which is operated by AURA, under a cooperative agreement with the National Science Foundation (NSF) on behalf of the Gemini partnership, which comprises the NSF (US), the Particle Physics and Astronomy Research Council (UK), the National Research Council (Canada), CONICYT (Chile), the Australian Research Council (Australia), CNPq (Brazil) and CONICET (Argentina).
Note that the lower dust mass and higher dust temperature derived from recent observations of the SN 1987A ring by Bouchet et al. (2006) correspond to a very different physical regime of dusty gas that is heated and partly destroyed by the SN shock, so these do not necessarily reflect the properties of dust originally surrounding the progenitor.