Issue |
A&A
Volume 595, November 2016
|
|
---|---|---|
Article Number | A120 | |
Number of page(s) | 13 | |
Section | Stellar structure and evolution | |
DOI | https://doi.org/10.1051/0004-6361/201628770 | |
Published online | 10 November 2016 |
Dissecting a supernova impostor’s circumstellar medium: MUSEing about the SHAPE of η Carinae’s outer ejecta ⋆
1 ESO – European Organisation for Astronomical Research in the Southern Hemisphere, Alonso de Cordova 3107, Vitacura, Santiago de Chile Chile
e-mail: amehner@eso.org
2 Instituto de Astronomía, Universidad Nacional Autónoma de México, Apdo Postal 106, Ensenada 22800, Baja California, Mexico
3 School of Physics, Trinity College Dublin, Dublin 2, Ireland
4 ESO – European Organisation for Astronomical Research in the Southern Hemisphere, Karl-Schwarzschild-Straße 2, 85748 Garching, Germany
5 Department of Astronomy, University of Minnesota, Minneapolis, MN 55455, USA
6 School of Physics and Astronomy, The University of Leeds, Leeds, LS2 9JT, UK
Received: 22 April 2016
Accepted: 30 September 2016
Aims. The role of episodic mass loss is one of the outstanding questions in massive star evolution. The structural inhomogeneities and kinematics of their nebulae are tracers of their mass-loss history. We conduct a three-dimensional morpho-kinematic analysis of the ejecta of η Car outside its famous Homunculus nebula.
Methods. We carried out the first large-scale integral field unit observations of η Car in the optical, covering a field of view of 1′× 1′ centered on the star. Observations with the Multi Unit Spectroscopic Explorer (MUSE) at the Very Large Telescope (VLT) reveal the detailed three-dimensional structure of η Car’s outer ejecta. Morpho-kinematic modeling of these ejecta is conducted with the code SHAPE.
Results. The largest coherent structure in η Car’s outer ejecta can be described as a bent cylinder with roughly the same symmetry axis as the Homunculus nebula. This large outer shell is interacting with the surrounding medium, creating soft X-ray emission. Doppler velocities of up to 3000 km s-1 are observed. We establish the shape and extent of the ghost shell in front of the southern Homunculus lobe and confirm that the NN condensation can best be modeled as a bowshock in the orbital/equatorial plane.
Conclusions. The SHAPE modeling of the MUSE observations provides a significant gain in the study of the three-dimensional structure of η Car’s outer ejecta. Our SHAPE modeling indicates that the kinematics of the outer ejecta measured with MUSE can be described by a spatially coherent structure, and that this structure also correlates with the extended soft X-ray emission associated with the outer debris field. The ghost shell immediately outside the southern Homunculus lobe hints at a sequence of eruptions within the time frame of the Great Eruption from 1837–1858 or possibly a later shock/reverse shock velocity separation. Our 3D morpho-kinematic modeling and the MUSE observations constitute an invaluable dataset to be confronted with future radiation-hydrodynamics simulations. Such a comparison may shed light on the yet elusive physical mechanism responsible for η Car-like eruptions.
Key words: stars: individual: ηCarinae / stars: emission-line, Be / stars: evolution / stars: massive / stars: mass-loss / stars: winds, outflows
© ESO, 2016
1. Introduction
Very massive stars with initial masses above 100 M⊙ may lose much of their mass in violent pre-supernova (SN) eruptions, which profoundly affect the appearance of spectra and lightcurves in the subsequent SN (Smith 2014). In several cases, giant outbursts during the luminous blue variable (LBV) phase have been confused with low-luminosity SNe IIn (van Dyk & Matheson 2012). Eta Car’s 1843 “Great Eruption” is the most famous of such events and the star has become the proto-type of the so-called SN impostors (van Dyk 2005; Davidson & Humphreys 2012a). In a few cases, LBV-like outbursts have been followed by a genuine core-collapse SN (Smith et al. 2007; Pastorello et al. 2007; Kiewe et al. 2012; Ofek et al. 2014; Moriya et al. 2014).
Eta Car is a binary, and a very luminous and massive system (see Davidson & Humphreys 2012a, for a recent summary, and references therein). The star is surrounded by a complex circumstellar environment, which is among the angularly largest of all LBVs. Its famous bipolar nebula, called the Homunculus, was ejected in the 1840s during the Great Eruption. It has an extent of approximately 18′′ (0.2 pc) and contains a mass of 10−35 M⊙ (Smith & Ferland 2007). The Homunculus expands with velocities of up to 650 km s-1. Several models for its formation exist, which include stellar interior and atmospheric instabilities, a fast rotating star, and binary interaction (e.g., Langer et al. 1999; Dwarkadas & Balick 1998; Frank et al. 1995; González et al. 2010; Maeder & Desjacques 2001; Smith & Townsend 2007; Kashi & Soker 2010). A smaller bipolar nebula, called the Little Homunculus, was likely ejected in the 1890s and has approximately 1 M⊙ (Ishibashi et al. 2003; Smith et al. 2004). It is contained within the larger Homunculus.
Beyond the bipolar Homunculus there are additional nebulous features referred to as the outer ejecta with a total mass of >2−4 M⊙ (see Weis 2012, for a review and Fig. 1). The outer ejecta are distributed in a region of 60′′ (0.7 pc) in diameter and include a variety of structures of different sizes and morphologies (Walborn 1976; Weis et al. 1999). A number of puzzling structures exist, such as highly collimated filaments and a jet-like feature. A global bipolar expansion pattern of these outer ejecta was also recognized, but the morphology appears irregular. Kinematic analysis showed the bi-directional expansion of the outer ejecta with velocities between 400−900 km s-1 and the same symmetry axis as the Homunculus (Weis et al. 2001). An expansion shell in front of the southeast lobe of the Homunculus also matches the global expansion pattern (Currie et al. 2002).
The 19th century Great Eruption of η Car ejected more than 10 M⊙ at speeds of 650 km s-1, with a kinetic energy of approximately 1050 erg (Smith et al. 2003), compared to 1050 − 1051 erg in kinetic energy for a typical core-collapse SN. Some features in the outer ejecta are moving at up to 3500−6000 km s-1 (Smith 2008), which approximately doubles the total kinetic energy, and suggests that the eruption released a blast wave, which is led by a shock front of compressed gas. The high velocities cause strong interaction with the surrounding medium and give rise to soft X-ray emission. As a result, η Car’s outer shell mimics a low-energy SN remnant. The eruption may have been powered by a deep-seated explosion rivaling a SN, perhaps triggered by pulsational pair instability (Heger & Woosley 2002). Pulsational pair instability induces strong pulsations of the entire star and part of the outer envelope is ejected (Heger & Woosley 2002; Woosley et al. 2007; Waldman 2008). The number of pulsations, the time scales between them (several weeks to a few years), and the mass ejected, depend on the mass of the CO core (Yoshida et al. 2016). Eta Car’s brief brightening events during the Great Eruption may be associated with such pulsations (see Smith & Frew 2011, for the historical lightcurve). Most models for η Car’s Great Eruption, however, assume a radiation-driven super-Eddington wind (Humphreys & Davidson 1994; Davidson & Humphreys 1997; Shaviv 2000; Owocki et al. 2004; Davidson 2016).
While the Homunculus is mainly a reflection nebula, the outer ejecta are emission nebulae. Chemical composition analysis revealed an overabundance of nitrogen from CNO processed material in the outer ejecta (Davidson et al. 1982, 1986), supporting their formation during the evolved stellar phase. However, if the star undergoes tidal mixing in a binary or is a fast rotator, apparent overabundances may be due to enhanced mixing. Smith & Morse (2004) find that the more distant clumps have a lower nitrogen abundance. This could be explained by interactions of younger ejecta with unprocessed material from the previous stellar mass loss, but the outer ejecta might also be the relic of an earlier eruptive phase (González et al. 2004). A comprehensive model for the formation of the outer ejecta is still missing.
![]() |
Fig. 1 HST/ACS F658N filter image, reproduced from Weis et al. (2004). The field of view is 1.5′× 1.5′. Several distinct features are identified with the terminology used by Walborn (1976). The main filter transmission is between λ6549−6623 Å, ranging over 3000 km s-1 in velocity space. Eta Car’s outer ejecta have a large range of velocities and thus [N ii] λλ6548, 6583 and Hα emission contribute to the features observed here. |
We investigate the outer ejecta of η Car using data obtained with the Very Large Telescope Multi Unit Spectroscopic Explorer (VLT MUSE; Bacon et al. 2010). This work is complementary to recent work on the three-dimensional structure of the Homunculus (Steffen et al. 2014). Walborn (1976) was the first to review the detailed structure of η Car’s outer ejecta and his terminology is still in use today. Figure 1 provides a map to some of the features discussed in this paper. This narrow-band image does, however, not reveal the detailed morphology of the outer ejecta discussed in this paper, because of their large velocity range and line blending.
In Sect. 2 we describe the MUSE observations and the modeling with the three-dimensional morpho-kinematic code SHAPE (Steffen et al. 2011). In Sect. 3 we present our results on the structures of the outer ejecta, and explore the spatial correlation between them and the extended X-ray emission around η Car. In Sects. 4 and 5 we discuss our findings and present our conclusions.
2. Data and analysis
2.1. Observations
Eta Car was observed with VLT MUSE over three nights (2014 November 13, 2014 December 12, and 2015 January 8). MUSE is an Integral Field Spectrograph, composed of 24 integral field units (IFUs) that sample a continuous 1′× 1′ field of view. The instrument covers most of the optical and part of the near-infrared domain (4800−9300 Å) with a spectral resolving power of R ~ 1600 (blue) to R ~ 3600 (red), corresponding to a median velocity resolution of ~100 km s-1. The spatial sampling for MUSE’s currently offered wide field mode is 0.̋2. We obtained seeing-limited spatial resolutions between 0.̋8 and 1.̋3.
Eta Car’s circumstellar material with an extent of approximately 60′′, velocities of several hundred to a few thousand km s-1, and numerous optical emission lines is ideally suited for observations with MUSE. However, the bright central source causes strong artifacts (vertical and horizontal stripes, concentric rings) at all wavelengths, probably due to internal reflections in the spectrographs and the fore-optics module. These artifacts are located close to the central star and span mostly over the extent of the Homunculus. They are spatially stable (i.e., they occur in all velocity channels at the same position and, therefore, cannot be confused with intrinsic features of the nebulosities), but vary in intensity across the wavelength range of the reduced single MUSE data cube, where the intensity variations are driven by the spectrum of the bright offending source. We remove these artifacts from the different emission line channel maps presented here by subtracting median-averaged white light images constructed from nearby emission line-free regions. The rapid intensity variations (with wavelength) of the artifacts renders this procedure difficult and some residual artifacts remain. The data also show the usual diffraction spikes around bright sources. Our treatment of the data leaves residuals. However, because these are spatially stable, they can easily be identified and avoided in our analysis.
MUSE provides the first complete and velocity resolved set of optical line ratios and line centroids at each position of the large-scale nebula around η Car. Observations with exposure times from 0.02 s to 300 s were obtained with η Car in the center of the field of view. Additional exposures of 300 s and 600 s were obtained at offset positions 35′′ and 70′′ along the Homunculus axis and perpendicular to it. In this paper, we only use the observations centered on η Car, covering a field of view of 1′ × 1′. The channel maps shown in Fig. 2 are extracted from a 300 s exposure.
The data were reduced using version 1.1.0 of the MUSE ESO standard pipeline (Weilbacher et al. 2012). Bias, arc, and flatfield master calibration files were created using the default set of calibration exposures provided by ESO. These include an illumination correction flatfield taken within one hour of the science observations to correct for potential temperature variation in the illumination pattern of the slices. (The MUSE field of view is split into 24 slices, which are sent to separate IFUs. These 24 channels are further sliced into 48 15″ × 0.2″ slitlets.) All data cubes were resampled to .
QFitsView1 was used to extract channel maps, position-velocity diagrams, and circular aperture spectra. Figure 2 shows selected channel maps of the Hβ emission in η Car’s outer ejecta. Each panel displays a single wavelength channel of the MUSE data cube, corresponding to a narrow velocity range of approximately 80 km s-1. We analyzed the velocity space from − 6000 to 6000 km s-1 for a few selected lines. Hα and the nearby [N ii] lines at λλ6548, 6583 are very strong in emission in η Car’s outer ejecta and can have similar line peak strengths. Line blending due to the large radial velocities of the ejecta and multiple components per line of sight is a serious issue for the [N ii] and Hα lines. [N ii] λ6548 is − 675 km s-1 and [N ii] λ6583 is 945 km s-1 in velocity space from Hα, very common velocities in the outer ejecta. Clearly disentangling each of these lines is a challenging task. Here, we thus chose the Hβ emission line as the primary line to analyze the three-dimensional morpho-kinematic behavior of η Car’s outer ejecta. Another suitable, but fainter, emission line for this type of work is [S iii] λ9069. Both lines have little contamination from nearby lines of Fe ii, [Fe ii], Cr ii, and N i. The channel maps shown in Fig. 2 are presented with a square-root color stretch and were normalized in flux individually to highlight the location of the outer (and fainter) ejecta. Absolute brightness distribution is not required for the work presented here, which focuses on the structure and dynamics of the outer ejecta.
Unless mentioned otherwise, the velocities in the following sections are radial velocities and the distances are the values projected onto the plane of the sky. Velocities are not corrected for η Car’s systemic velocity of roughly −8 km s-1 (heliocentric; Davidson et al. 1997; Smith 2004). Throughout the paper we assume a distance to η Car of 2.35 kpc (Allen & Hillier 1993; Davidson et al. 2001; Smith 2006). Quoted wavelengths are air values.
2.2. Morpho-kinematic modeling with SHAPE
![]() |
Fig. 2 MUSE channel maps showing the observed Hβ emission in the MUSE field of view of 1′× 1′. Each panel shows a single channel of the MUSE data cube, corresponding to different velocities of the Hβ emission. The velocities are indicated in each panel in units of km s-1 (each channel map spans a range of approximately 80 km s-1, velocities are not corrected for η Car’s systemic heliocentric velocity of approximately −8 km s-1). The corresponding SHAPE model channel map is shown on the right side of each observational channel map (the grayscale only indicates the location of the emission and does not reproduce the brightness distribution), see text for details. The channel maps are irregularly spaced in velocity to show best the brightest emission features corresponding to the outer shell, the ghost shell, the bowshock, and the jet, identified by labels and arrows pointing toward the feature (see also Fig. 4 for their three-dimensional location). The jet is not modeled with SHAPE, since the origin of the emission is not clear, see Sect. 3.4. The central star and part of the Homunculus are saturated. The channel maps are presented with a square-root color stretch and normalized in flux individually to peak intensity. They only contain physical information on the structure and velocities, but do not permit channel-to-channel comparisons of the fluxes to be made. |
![]() |
Fig. 3 MUSE channel maps showing the observed Hβ emission in the MUSE field of view of 1′× 1′ as in Fig. 2, but overlayed with the brightness outlines of the model features. Only a few contours are shown to avoid obstructing the underlying observation. |
From the spatio-kinematic MUSE data, we construct a three-dimensional geometry of η Car’s outer ejecta with the morpho-kinematic modeling code SHAPE (Steffen et al. 2011). SHAPE has been used extensively to model the complex three-dimensional structure of planetary nebulae. Recently, Steffen et al. (2014) constructed the first three-dimensional model of the Homunculus based on kinematic molecular hydrogen data from VLT X-shooter data. The X-shooter data cover the entire Homunculus nebula, but none of the outer ejecta.
SHAPE enables the user to interactively optimize model parameters by direct comparison of model predictions with observations. The general modeling procedure is described in Steffen et al. (2011) and on the SHAPE support website2. The program consists of a three-dimensional modeling view in which the geometry and behavior of the model are defined, and a two-dimensional view where the simulated appearance of the model can be compared to observational data. SHAPE models are constructed interactively with three-dimensional structural mesh elements, such as spheres or cylinders, based on the user’s assumptions about the object’s geometry, and the emission and the velocity information extracted from observations. Small-scale structure can be added in a variety of ways by modifying the elementary meshes. The mesh models can also be combined with hydrodynamic simulations performed with the corresponding SHAPE module. We have used this technique to model the bowshock structure to the northeast of the Homunculus. The three-dimensional model is the input to a rendering module, which emulates how the object would be observed with an imager or spectrograph. The user can then interactively refine the model until it qualitatively fits the observational data. In the case of our MUSE data, we modeled several large-scale features, which were selected based on their apparent spatial association as a contiguous structure in the observational data. This includes the larger shells, the ghost shell, and the bow-shock-like feature, described in Sect. 3. Other, small-scale features were not included, because they do not provide additional constraints for our morpho-kinematic modeling. Quantitative fits of SHAPE models to our MUSE observations are not possible within the SHAPE tool. However, such an analysis is not required here, because we are limited by other factors described in the following paragraph. Given an assumption for the velocity field, the model solution is unique with an additional geometric constraint, such as an approximately circular cross-section, for example.
A cautionary note: reconstruction based on Doppler-velocity measurements of expanding nebulae often assumes homologous expansion, that is, radial expansion with the velocity increasing linearly with distance (). We adopt homologous expansion for the structures analyzed in this paper with the exception of the bowshock-like feature discussed in Sect. 3.4, where a hydrodynamical simulation was performed. This is because the events that produced the structures were short compared to the dynamic time scale of the expansion, thereby resembling more a ballistic blast wave than a continuous wind shock. A homologous flow may not be valid in detail or throughout the expanse. Ejecta may be slowed down when hitting older material and by interactions with the environment, or may have been accelerated by a blast wave. Interaction with preexisting density fluctuations alter the velocity field and distort the reconstructions to an extent that depends on the time-scale of the interaction (Steffen et al. 2009). Assuming an incorrect velocity field will lead to distortions in the reconstructed shape. Refinements of the three-dimensional structure model will require improvements of the constraints on the ejection time and velocity field. Our reconstruction of the outer ejecta thus only gives a suitable first order approximation. We urge the reader to keep this in mind.
3. Results
![]() |
Fig. 4 Three-dimensional model of η Car’s outer ejecta seen from Earth (left) and from the side (right) with the bipolar Homunculus nebula in the center (dark red and blue lobes). The model has been constructed with SHAPE and the Hβ emission and velocity (over a range of ± 3000 km s -1) information extracted from the MUSE data. The outer ejecta form a partial bent cylinder. The open cylinder parts are not observed in the MUSE data, i.e., we only trace the bottom of the shell moving toward us and the top of the shell moving away from us. The open cylinder parts are likely too faint, perhaps due to only weak interactions with less dense surrounding material. This large tube is probably bent or indented due to material being slowed down by interactions with the environment. The ghost shell (green) is in our line of sight to the southern Homunculus lobe and consists of two closely aligned shells. An additional half-shell (green mesh, complementing the bottom of the large cylindrical shell moving away from us) lies behind the northern Homunculus lobe. The bowshock-like feature is displayed in yellow. The colors are selected to distinguish the different shapes selected for the SHAPE modeling. |
![]() |
Fig. 5 Three-dimensional velocity structure of η Car’s outer ejecta seen from Earth (left) and from the side (right). The color coding is from −2500 km s-1 (blue/violet) to +2500 km s-1 (red). This SHAPE model shows only qualitatively where material is located. The colors of the velocity vectors and the vector lengths show roughly the direction and speed of the material. |
With MUSE IFU observations we uncover the apparently irregular structure of η Car’s outer ejecta and show for the first time their spatially contiguous large-scale morphology. In the past, narrow-band images excluded the high velocity ejecta and contained emission from more than one atomic transition. This made the detection of a contiguous morphology difficult. Broad-band images were dominated by the bright central source and the Homunculus nebula, and were not deep enough to show the structure of the outer ejecta. Unlike the bipolar Homunculus, η Car’s outer ejecta thus did not appear symmetric or connected (Fig. 1), but seemed to be composed of many irregularly shaped and randomly placed structures in an elliptically shaped region (described as filaments, jets, arcs, bullets or knots, and strings). By integrating the MUSE data over the HST WFPC2 filter transmission curves, we identify the same features observed in the HST images (filtered to match the angular resolution of MUSE) and can obtain their proper motions with a baseline of 17.5 yr by comparing 1997 HST images with our 2014 MUSE data. We use proper motion measurements to constrain the ejection time of the NN bow, one of the most prominent features in η Car’s outer ejecta and discussed in Sect. 3.4.
MUSE provides the first complete spatio-spectral sampling of η Car’s outer ejecta. The outer ejecta are indeed more clumpy than the material in the Homunculus, but most of the knots outline a coherent superstructure. Figures 2 and 3 show selected observed MUSE velocity channel maps and the corresponding SHAPE model channel maps. Rough morpho-kinematic modeling with SHAPE of the emission features visible in the MUSE position-velocity diagrams across the Hβ emission line reveal a continuous morphology of η Car’s outer ejecta that is not a simple bipolar shell. The three-dimensional model is shown in Figs. 4 and 5 and described in detail in the following subsections. The largest continuous structure surrounding η Car can be identified as the “cocoon” or “blast wave” described in Smith & Morse (2004) and Smith (2008). We also determine the approximate shape and extent of the so-called ghost shell in front of the southeast Homunculus lobe (Currie et al. 2002). In addition, we discuss the jet-like NN bow and condensation in the equatorial plane (Meaburn et al. 1993).
3.1. Prominent features in the outer ejecta
Before describing the main large-scale features in η Car’s outer ejecta (i.e., the outer shell, the ghost shell, and the NN bow) in the following subsections, we give a short description of the prominent features in Fig. 1 and how the MUSE data improves their interpretation. We follow the terminology used by Walborn (1976).
-
The E condensations are five discrete knots southeast of the Homunculus (Walborn 1976; Weis 2012). The MUSE data reveal the velocity structure of these condensations. They have a range of negative velocities. The condensation E5 has the fastest velocities, with speeds of v ~ − 600 km s-1 and consists of several knots and arc-like features. E2 and E4 can be observed up to velocities of v ~ − 450 km s-1. E2 and E3 can be better described as arcs than knots and are visible up to velocities of v ~ − 300 km s-1. They can be seen in several of the channel maps in Fig. 2. The E condensations are not part of the continuous shell discussed in Sect. 3.2, such that their nature remains particularly mysterious.
-
In narrow-band images, the S ridge appears to be the largest structure in η Car’s outer ejecta (Walborn 1976; Weis 2012). The S ridge is composed of many filaments and knots with velocities between v ~ − 300 km s-1 in the southwest and up to more than v ~ 1000 km s-1 in the northwest. MUSE data show that this ridge is part of the large outer shell discussed in Sect. 3.2. The S ridge and S condensation have not been modeled as separate features, because they appear to be part of the overall structure of the outer shell (Figs. 2 and 4).
-
The W arc (Walborn 1976; Weis 2012) is also part of the same outer shell structure and can be observed up to velocities of v ~ 2000 km s-1. The W condensation, on the other hand, is a distinct knot at smaller radial velocities of v ~ 350 km s-1.
-
The NN, NS condensations, and the NN bow are intriguing features discussed in detail in Sect. 3.4. They lie within the equatorial plane and can be observed over a wide range of negative and positive velocities (± 1000 km s-1).
-
Five long, highly collimated linear features, called strings (Weis et al. 1999), are perhaps some of the most interesting structures in the outer ejecta. They are found in the northwest and southeast, pointing radially away from the Homunculus (they are faint and thus not visible with the cuts used in Fig. 2). Their velocity structure is resolved in the MUSE data. The velocities increase with distance up to v ~ − 1000 km s-1. The MUSE data do not allow for greater constraint of their three-dimensional motion than previous work (Weis et al. 1999).
3.2. The outer shell
MUSE channel maps across Hβ reveal a large shell-like emission feature with velocities between − 2200 km s-1 southeast of the central source and + 2000 km s-1 northwest of the central source (Fig. 2). At a velocity of − 2200 km s-1, the shell has an extent of approximately 30′′ toward the southeast of the central star. For velocities between − 2200 km s-1 and 0 km s-1, the emission of this shell appears to gradually move toward the central source. The shell is oriented along a similar direction as the southern Homunculus lobe. At positive velocities, a complementary shell structure becomes evident northwest of the central star, which can be traced up to a distance of approximately 30′′ and velocities up to + 2000 km s-1. With deeper images, faster velocities may be detected and the structure may be found to extend beyond what is visible in our MUSE data set.
In the MUSE data we measure the radial velocities of the ejecta and the projected distances. We then have to adopt an ejection time in order to determine the distances of material from the star – unless we assume a priori a geometry. For the SHAPE modeling of this large shell-like structure seen in emission, we use a homologous expansion (i.e., material is neither decelerated nor accelerated) for the ejecta and set the ejection date at the time of the Great Eruption. The selected velocity field corresponds to an ejection of the material in 1836. Within the accuracy of the model and observations, an ejection date of 1843 could be used as well. We adopt the same age for the outer ejecta as for the Homunculus (i.e., ejected in the same physical event) for the following reason. If we assume a homologous flow with the geometrical constraint that the tube of the outer ejecta has roughly a circular cross-section, we find an age for the outer ejecta 1.4 times older than the Homunculus. But no major eruptions were reported in η Car in the two centuries before the Great Eruption (Smith & Frew 2011, though only few data points are recorded from the 17th century).
With SHAPE we derived that the emission of the outer shell can be matched with a bent cylinder geometry embedding the entire Homunculus nebula (Fig. 4). This structure is reminiscent of the shape of the planetary nebula M 2−9 (Clyne et al. 2015). The bending of the outer shell might have been caused by the interaction with a denser region of the ambient medium, slowing down the material. The reconstructed bend might be somewhat exaggerated, because the interaction with surrounding material might only have started recently compared to the age of the outer shell. The actual bending may therefore be less pronounced than depicted in the current reconstruction, but will increase with time. In Appendix A we show a SHAPE model, which assumes a scaled-up Homunculus geometry for the outer ejecta. It illustrates that a Homunculus geometry does not match the MUSE data.
The emission brightness is irregularly distributed over the extent of the cylinder. Brighter knots are observed and the cylinder has its strongest emission on opposite sides for its blue and red components (see the structure labeled “outer shell” in Fig. 2). The emission originates from the bottom half of the cylinder for the blue component (viewed from our line of sight), while most of the emission comes from the top half for the red component. These bright parts lie along the line where the cylinder wall intersects with the Homunculus axis. That we only observe a partial cylinder is probably because the other parts are too faint to be detected in our data, perhaps due to only weak interactions of the ejecta with less dense surrounding material in those directions. The shell depicted as a green mesh behind the northern Homunculus lobe may be completing the tube, but we cannot clearly establish this (Figs. 4 and 5).
The orientation of the cylinder spatially close to η Car is misaligned with respect to the Homunculus axis (Fig. 4), which has an inclination of 41° (Davidson et al. 2001; Smith 2006). Further away from η Car, the inclination of the cylinder is approximately aligned with the Homunculus axis. Projected on the plane of the sky, the cylinder is rotated by approximately −10° with an estimated error of approximately 3° with respect to the Homunculus. The outer ejecta are not a perfect cylinder, but locally the orientation varies. The difference in orientation with respect to the Homunculus could be a secondary effect of the interaction with an inhomogeneous environment.
The diameter of the cylinder section moving toward us varies between 1.0 × 1018 cm and 1.2 × 1018 cm (0.32–0.39 pc). The diameter of the cylinder section moving away from us is noticeably smaller at approximately 6.7 × 1017 cm (0.22 pc). The diameters vary by approximately 20% with position. The thickness of the shell is approximately 15% of its radius. The material behind the star probably runs into denser material and thus cannot expand as easily as in front of the star, resulting in the smaller cylinder diameter. Hβ emission tracing the component in front of the star is smooth and continuous, while the emission of the component behind the star is more irregular and clumpy. There is also more soft X-ray emission in the northwest than in the southeast of η Car. This is consistent with stronger interactions between the shell and the surrounding material on the red side. We will discuss this further in Sect. 3.5. Note that with SHAPE we can only model the brightness distribution of the ejecta, but do not obtain any information on the density distribution.
3.3. The ghost shell(s)
In the MUSE channel maps, two horn-like features can be seen southeast of η Car. The first horn-like feature becomes apparent at velocities of approximately − 900 km s-1 (labeled “ghost shell” in the −841 km s-1 observed and model channel maps in Fig. 2) and the second at velocities of approximately − 600 km s-1. These emission features can be modeled by two adjoining shells, which lie along our line of sight immediately in front of the southern Homunculus lobe and are moving toward us (green mesh in Fig. 4). Figure 2 shows a good fit between the MUSE observations and the SHAPE model. These two shells appear to align with the southern Homunculus lobe. Behind the northern Homunculus lobe, we identify a structure that extends approximately twice as far as the size of the Homunculus with projected velocities larger than + 400 km s-1 (green mesh in Fig. 4). It is not clear if this shell is the counterpart of the ghost shell in front of the southern lobe or part of the outer shell.
We identify the two shells in front of the southern Homunculus lobe with the ghost shell discussed in Currie et al. (2002). Currie et al. (2002) detected the ghost shell in emission in multiple Balmer lines and in forbidden lines such as [N ii], [S ii], and [Ar iii] at velocities between −675 km s-1 and −850 km s-1. It is probably also associated with the complex absorption structure of Ca H and Ca K lines (Davidson et al. 2001) and the −513 km s-1 absorption system in the near-ultraviolet (Gull et al. 2006). Currie et al. (2002) hypothesized that the form of the shell is an approximate sphere with radius ~11′′ that surrounds the Homunculus, but that the quasi-spherical approximation is not valid near the central star and in the northwest lobe region. They propose that the shell is the forward shock between the fast stellar wind of the Great Eruption and the older slow massive stellar wind, distorted by ejecta such as the equatorial disk.
MUSE enables us to show clearly the structure and extent of the ghost shell (Fig. 4). We identify two adjacent shells, which are located immediately outside the southeast lobe of the Homunculus. They have a similar shape as the southern lobe, but velocities several hundreds of km s-1 faster than the Homunculus material.
3.4. The NN bow and the “jet”
In the Homunculus equatorial plane, toward the northeast, is a remarkable structure, identified as NN bow and NN condensation (Fig. 1). Meaburn et al. (1993) interpreted this prominent knot as the interaction of a jet with the ambient gas. In the MUSE velocity channel maps, we observe a large feature at negative velocities that is reminiscent of a bowshock, but not the corresponding jet. The emission feature has the structure of a giant loop and has a radial velocity of approximately − 500 km s-1 at its furthest extent from η Car at the NN condensation. The entire structure has a large velocity range between − 850 km s-1 and − 200 km s-1. We measured the proper motion of the NN bow between July 1997 and December 2014 using a HST WFPC2 F656N image and our MUSE data integrated over the WFPC2 F656N filter transmission curve. The displacement of the tip of the NN bowshock (the NN condensation) in the 17.5 yr between these two epochs is 1.5−2′′ (corresponding to a transverse velocity of approximately 200 km s-1; compared to Kiminki et al. 2016, who find proper motions of 200 km s-1 to 1400 km s-1 for this region).
Instead of using simple structural elements, we model this bowshock-like feature using the numerical hydrodynamics extension of SHAPE (Steffen et al. 2013). The measured radial velocity and proper motion provide some constraints for the model parameters. We find that the observations can be matched by a dense jet with a spatial bowshock velocity of 1300 km s-1 running into an inhomogeneous medium (yellow mesh in Fig. 4). The jet originates 10° ± 10° from the plane of the equatorial skirt toward the blue lobe, roughly perpendicular to the Homunculus axis, and at approximately 25° ± 5° from the plane of the sky toward the observer. Its deprojected extent from the star is approximately 6.9 × 1017 cm (0.22 pc). Figure 2 shows the observed channel maps next to our model. We observe the bowshock, but not the corresponding jet in the MUSE data at negative velocities.
The region where this bowshock-like feature is observed is faint in X-ray emission, which suggests absorption of the soft X-rays by an extension of the equatorial disk or by the bowshock itself (Weis et al. 2004). The lack of soft X-rays may also be caused if the material was expanding into a low-density region. In this case, the bowshock-like structure could be explained by material expanding freely into the low-density ambient medium, in contrast to other outflow directions.
At the same spatial location, but at positive velocities, a very bright collimated emission resembling a jet is visible (Fig. 2). Velocities along this “jet” increase only slightly with distance to the central source. The emission brightness peaks at 500−600 km s-1. The emission line profile is very broad and extends over approximately 800 km s-1. A jet at positive velocities cannot cause a bowshock at negative velocities and the relation of this feature to the bowshock-like structure is thus questionable. The spectrum of this “jet” shows continuum emission and strong emission from hydrogen and helium lines present in the stellar spectrum, only red-shifted. Thus, this structure most likely results from processes including reflection and/or scattering of light from the central star by an extension of the equatorial skirt. The bright structures observed at negative and positive velocities may not be directly related, but the equatorial plane is a preferred direction for outflows and scattering processes. This structure has a possible counterpart toward the southwest with similar projected velocities.
3.5. Correlation of the outer ejecta with the soft X-ray emission
X-ray images obtained with Einstein, ROSAT, and Chandra show low-energy X-ray emission (0.2−1.5 keV) in a hook shape around η Car (Fig. 6; Seward et al. 1979; Weis et al. 2001). This soft X-ray emission is produced as ejected material runs into the surrounding gas and dust. Two brighter knots can be identified with the optical S ridge and the W arc, where also medium-energy X-ray emission (1.5−3.0 keV) is observed. There is very little soft X-ray emission in the southern part of the outer ejecta, where the density of optical clumps is also low. The central source and the Homunculus nebula are detected only in hard X-ray bands (3.0−8.0 keV) with the emission from the Homunculus being reflected X-ray emission (Corcoran et al. 2004). No hard X-ray emission is detected outside the Homunculus.
The X-ray brightness is related to the combination of the speed of the ejecta and the density of the ambient medium. In the past, it was thus surprising that optical images showed only few direct correlations with the X-ray images, though optical emission was observed in most of the regions where X-ray emission is found (Weis et al. 2001, 2004). On the other hand, the expansion velocities of the ejecta were found to be in good agreement with the X-ray brightness distribution and X-ray temperatures, that is, regions with higher X-ray emission have higher expansion velocities (Weis et al. 2001, 2004). The temperature of the X-ray gas is of the order of 0.65 keV, indicating velocities of the shocking gas on the order of 750 km s-1 (Weis et al. 2004). This is consistent with the velocity structure seen in the MUSE data.
![]() |
Fig. 6 Composite images displaying our SHAPE model together with the 0.5–1.2 keV Chandra X-ray image. The soft X-ray emission is produced as material runs into the nearby gas and dust. X-ray emission from the Homunculus is detected only in hard X-ray bands. The large-scale optical tube structure fits right inside the X-ray bubble. The SHAPE model contours only qualitatively show the outlines of the emission. The blue and red colors are to distinguish the different shapes. There is no physical information associated with them. For the image scale, see Fig. 7. |
The optical emission and velocity information from our MUSE data together with the SHAPE modeling let us predict where the X-ray emission (in three dimensions) should be and see where it actually is. With the complete spatio-spectral sampling at optical wavelengths in the MUSE data we are able to recognize the excellent correlation between the X-ray emission and the optical ejecta. The large outer shell surrounding η Car aligns tightly inside the X-ray bubble and some bright X-ray features can be associated with counterparts in the optical (Figs. 6 and 7).
Figure 7 shows the Hβ intensity contours in the different MUSE cube wavelength slices, color-coded as a function of the corresponding line of sight velocity of the emitting gas. These contours are overlaid on top of the 0.5−1.2 keV Chandra X-ray image. The analysis is restricted to the regions outside the Homunculus due to the artifacts close to the central star in the final datacube, which are described in Sect. 2. The figure shows the regions (purple shapes) that we mask out for this analysis. They contain the entire Homunculus, the concentric reflection rings around the central star, and several nearby stars. At negative velocities the intensity contours at different velocities clearly trace the expanding shell and the largest velocities align with the bright X-ray emission rim in the southeast. The bowshock is also distinctly resolved. The Hβ emission at positive velocities is fainter and even very weak artifacts in the data make it difficult to determine their correlation with the X-rays. We thus chose two different mask-out region sizes (Fig. 7, lower left and lower right) and different intensity levels to determine the Hβ intensity contours at different velocities. We find that for positive velocities also, the cylinder is clearly traced and the largest velocities align well with the bright X-ray hook. Note that the artifacts are primarily concentric reflection rings around the central source and horizontal and vertical stripes, which could be easily identified as artifacts in this analysis.
![]() |
Fig. 7 The 0.5–1.2 keV Chandra X-ray image (gray-scale image) overlaid with the Hβ velocity contours obtained from our MUSE data. Contours are fit in individual MUSE slices at an intensity level “Ic” in 10-20 erg/s/cm2/Å and color-coded as a function of velocity. The purple shapes indicate masked-out regions not used in this analysis that include the Homunculus, the artifacts from the optical ghosts described in Sect. 2, and several nearby stars. Because of the number of artifacts present and the faintness of the clumps, two different mask-out regions are used for the positive velocity contours. We also adjust the “Ic” value in each case. Upper left: extent of the soft X-ray emission surrounding η Car. The green box indicates the region shown in the other three panels. Upper right: contours of the blueshifted emission, highlighting the correlation of the optical emission with the soft X-rays in the southeast. The outer shell in the southeast coincides with the inner wall of the X-ray emission. Lower left and right: contours of the redshifted emission. The red part of the outer shell is fainter and it is difficult to determine the velocity contours of the ejecta due to artifacts, but the optical emission here also accurately traces the inner wall of the X-ray emission. |
The NN bow is not associated with strong soft or medium X-ray emission, which may point to a high foreground column density (Weis et al. 2004). The X-ray bridge between the northern and southern part of the X-ray nebula could be attributed to a large expanding disk, maybe an extension of the equatorial disk, which produces shocks at the outer edge (Duschl et al. 1995). The lack of soft X-ray emission from the NN bow may suggest that this feature lies in front of the disk and absorbs the X-rays (Weis et al. 2004). An alternative explanation for the absence of X-rays from the NN bow is that the material expands almost freely into a low density region (Weis et al. 2004).
4. Discussion
It is not yet resolved if η Car’s outer ejecta originate from eruptive events or high velocity winds prior to the Great Eruption, or if they were expelled during the Great Eruption. Eta Car’s history of sudden changes in its light curve and its spectrum, which have occurred quasi-periodically every ~50 yr since the Great Eruption (Humphreys et al. 2008; Davidson et al. 1999) may indicate that this eruption was not an isolated event. It is possible that eruptive events have occurred in the past, which could be held responsible for the formation of the outer ejecta. Several (theoretical) models aim to explain the origin of the Homunculus (for references see Sect. 1), but no working model exists yet for η Car’s outer ejecta. González et al. (2004) is the only example of a theoretical work which investigates if a high velocity pre-outburst wind just before the Great Eruption could have produced the ghost shell and the outer ejecta. Their work is inconclusive.
To better understand the nature of η Car’s Great Eruption, it is essential to determine the three-dimensional kinematics and geometry of the ejected material. The unique value of the MUSE data is that they enable us to sketch the three-dimensional structure of η Car’s outer ejecta. With the help of the observed velocities, kinematic bridges can be built between morphologically isolated islands.
4.1. The outer shell
Eta Car’s ejecta are chemically stratified, following a clear pattern of progressive nitrogen enrichment. Coincident with the soft X-ray shell are less nitrogen-rich ejecta than the material immediately outside the Homunculus nebula, and the gas beyond the X-ray shell has not been significantly processed through the CNO-burning cycle (Smith & Morse 2004; see also Tsuboi et al. 1997; Corcoran et al. 1998; Hamaguchi et al. 2007). The outer ejecta thus interact with the unprocessed material from previous stellar wind mass loss. This suggests that η Car has only been ejecting nitrogen-enriched material in the past few thousand years.
Fast material in η Car’s outer ejecta, in excess of 1000 km s-1, has been discussed since the late 1980s (e.g., Dufour 1989; Meaburn et al. 1996; Weis et al. 2001). Smith (2008) estimated velocities of up to 3500 − 6000 km s-1 and proposed that these high velocities indicate that the material has been accelerated by the pressure behind a “blast wave”, which originated from a deep-seated explosion. This blast wave involves only a small amount of high-velocity material originating from η Car’s SN impostor event and should not be confused with the blast waves seen in SN explosions. Rest et al. (2012) derived from light-echo spectra that η Car was very cool during the Great Eruption with an effective temperature of ~5000 K and argued that this supports the notion that a physical mechanism such as an energetic blast wave may be associated with the outer ejecta (Smith 2008), but see Davidson & Humphreys (2012b). The Great Eruption may thus have been powered by a deep-seated explosion. Candidates for an explosion mechanism could be a pulsational pair instability (Woosley et al. 2007) or another instability associated with nuclear burning in the last stages of evolution. These events are expected to occur only ~10−1000 yr before the final core collapse SN, as proposed for the progenitors of the extremely luminous SN 2006gy (Smith et al. 2007) and SN 2006jc (Pastorello et al. 2007).
Some models for η Car’s Great Eruption, however, assume a radiation-driven super-Eddington wind (Humphreys & Davidson 1994; Davidson & Humphreys 1997; Shaviv 2000; Owocki et al. 2004). A recent review by Davidson (2016) strongly advocates the idea that LBV giant eruptions are super-Eddington mass outflows. Owocki & Shaviv (2016) found that the low temperatures from light echoes of η Car’s Great Eruption are consistent with the very large mass-loss rates and luminosities estimated for this eruption epoch. However, Owocki & Shaviv (2016) also note that such a cool spectral temperature is not unique to a steady wind outflow model and cannot be used by itself to discriminate between the explosion and steady-wind scenarios.
As explained in Sect. 3.2, we can only roughly constrain the age of the outer shell surrounding η Car with the SHAPE modeling, because of the degeneracy between the velocity field and the linear distance of ejecta from the central star. Without any historical record of eruptions prior to the Great Eruption and assuming that material decelerates, it is most appropriate (for all intents and purposes) to presume that all outer shell material was ejected by the energy release that caused the Great Eruption. (Note, that observations in the 18th century were sparse and we cannot entirely rule out eruptions prior to the Great Eruption, see Sect. 3.2.) The shell has a thickness of approximately 15% of its radius, which implies an upper limit of 25 yr for its formation and is thus consistent with the formation during the Great Eruption (e.g., Smith & Frew 2011). The fact that the outer ejecta have a slightly different orientation than the Homunculus of ~10° does not indicate ejection in different directions and at different times. The cylindrical structure is also not perfectly straight and the orientation varies locally by a similar amount. Variations in the orientation could be due to interactions with the inhomogeneous environment. Without further study, the significance of these deviations cannot be assessed. Slower knots and arcs surrounding η Car, such as the E condensations, may be older, consistent with their less nitrogen-enriched chemical abundances.
Smith (2008) proposed a bipolar forward shock geometry similar to the Homunculus, but three to four times its size and speed for the shape of the outer ejecta. This would result in an up-scaled version of the Homunculus, yet our data are inconsistent with this idea. We find that the outer shell is not a scaled version of the Homunculus nebula. If we assume the same bipolar geometry as for the Homunculus, we cannot match the observed spatial distribution of the outer ejecta with a SHAPE model (Appendix A). The geometry we derive is instead more appropriately described by a bent tube (Fig. 4). Smith (2008) suggested that the interaction between the blast wave from η Car’s Great Eruption with 500−1000 yr old clumpy ejecta gives rise to the soft X-ray emission. We confirm a strong correlation between the outer shell and the soft X-ray emission with our MUSE data (Figs. 6 and 7, see also Weis et al. 2001). The outer shell tightly aligns inside the X-ray bubble and bright optical emission is associated with bright X-ray features.
We revise the ejecta velocities in Smith (2008), who estimated a range for the largest observed (deprojected) velocities in η Car’s outer ejecta between 3500 and 6000 km s-1, depending on the assumed geometry. From the derived three-dimensional geometry of the outer ejecta with the MUSE data and SHAPE modeling, we find that the largest deprojected velocities are approximately 4500 km s-1.
4.2. The ghost shell
Currie et al. (2002) proposed that the ghost shell in front of the southern Homunculus lobe is material associated with the forward shock between the fast stellar wind of the Great Eruption and the older, slower, massive wind. The reverse shock formed the leading edge of the Homunculus. This would be consistent with the lack of soft X-rays from the Homunculus, since the material would be expanding into a low density medium.
We find that the ghost shell consists of two aligned shells in front of the southern Homunculus lobe. The shells expand with velocities of only a few hundred km s-1 faster than the Homunculus material. A scenario that could accommodate this invokes several pulses of ejection during the Great Eruption (maybe a few years apart), during which material was expelled, forming the two ghost shells and the Homunculus nebula, but possibly also involving a later shock/reverse shock velocity separation.
4.3. The NN bow
The NN condensation and the associated NN bow consist of fast moving material (deprojected velocities are approximately 1300 km s-1, the deprojected extent from the star is approximately 0.22 pc) close to the plane of the equatorial skirt. Together with the measured proper motion between the 1997 HST WFPC2 images and the 2014 MUSE data, we estimate an ejection date approximately 170 yr ago, which (again) points to the Great Eruption. This suggests that the main structures in the ejecta surrounding η Car, that is, the Homunculus, the ghost shell, the NN bow, and the outer shell, are all products of the Great Eruption. It is not clear what mechanism could have ejected material with these velocities in a direction that is perpendicular to the Homunculus axis and located in the presumed orbital plane.
We have shown that a hydrodynamical model for a jet or “bullet” reproduces the morphology of the NN bow well. However, this does not mean that this feature is a jet, as suggested by Meaburn et al. (1993). There is no consensus on the mechanism for forming jets from eruptive stellar systems. For η Car, the interacting companion offers a plausible cause. The jet could have originated from an outflow following the accretion of material by the secondary star (Soker 2005; Akashi et al. 2013; Schreier & Soker 2016). In this case, however, a peculiar orientation of the accretion disk needs to be invoked. There might have been an explosive instability in the accretion disk around the secondary star. In this case, the lack of a counter jet is also more readily understood.
The NN bow resembles a slingshot prominence, that is, a large arc or loop. It is somewhat reminiscent of structures that can be found in the ejecta of the red supergiant VY CMa and the post-red supergiant IRC+10420 (e.g., Humphreys et al. 2007; Tiffany et al. 2010). In these objects the eruptive mass loss is probably due to convective activity and magnetic fields. However, while the NN bow has a similar morphology and spatial extent as the arcs seen in these objects, the observed velocities are 1−2 orders of magnitude larger. η Car is also at least an order of magnitude more luminous than VY CMa and IRC+10420. In contrast to η Car, the large arcs in VY CMa were ejected in different directions over several hundred years; in IRC+10420 the semi-circular arcs are equatorial. These objects are clearly very different from η Car, but perhaps the NN bow represents a counterpart of the arc-like, episodic eruptive events seen in these less luminous objects. The similar effective temperature reached by η Car during the Great Eruption (though not in hydrostatic equilibrium) may suggest that a similar physical mechanism is at play in both cases, and further theoretical work is warranted.
One could also imagine more exotic explanations for the NN bow. For example, Portegies Zwart & van den Heuvel (2016) discuss the Great Eruption in terms of a merger event in a triple system, with massive loops of ejected material. In their model, the Homunculus was produced prior to the merger by an extremely enhanced stellar wind, energized by tidal energy dissipation. The merger itself resulted in a massive asymmetric outflow (two loops) in the equatorial plane.
5. Conclusion
Three-dimensional information on the morphology and kinematics of η Car’s ejecta provides evidence for their ejection history and structure. In the past, η Car’s outer ejecta were seen as a conglomerate of individual small structures. This was regarded as evidence for very rapid and efficient fragmentation within the nebula or a mechanism which favors the ejection of clumps. Using MUSE IFU observations and the modeling tool SHAPE, we reconstructed the large-scale three-dimensional geometry of η Car’s outer ejecta and show their contiguous nature.
Our analysis suggests that the star is surrounded by a large bent partial cylinder (the outer shell), centered on the star and roughly aligned with the Homunculus. This structure fits tightly into the X-ray bubble. For the modeling with SHAPE, we have to make an assumption on the velocity field and on the ejection times of the ejecta. Based on several considerations, described in Sects. 3.2 and 4.3, the probable ejection time of the outer shell is during the Great Eruption and it is improbable that the expelled material is more than twice as old as the Great Eruption.
The existence of two ghost shells immediately outside the southern Homunculus lobe, and the NN outflow suggests a sequence of several outbursts during the Great Eruption, and/or later shock/reverse shock velocity separation of material. This confirms the notion of event-driven mass loss, breaks it down to smaller scales, and poses an invitation for theorists to search for an explanation.
In future work, it will be interesting to study in detail the correlation between the variations in X-ray temperature and the radial velocities of the ejecta. The MUSE data also provide for detailed velocity and abundance studies, which will shed further light on the chemical composition during different ejection episodes. However, the blending of lines due to the large velocity ranges of ejecta within the same line of sight makes this a difficult task. The MUSE data set and the SHAPE modeling show that three-dimensional geometry is needed to study η Car’s outer ejecta, against which future three-dimensional hydrodynamical and radiative transfer simulations should be compared to.
Note added in proof. During the second revision of this paper, Kiminki et al. (2016) was published, whose content is directly related to our paper. Kiminki et al. (2016) measured proper motions of approximately 800 features in η Car’s outer ejecta with a baseline of over 21 yr. They find that the velocities are consistent with features moving at constant velocities, as we have assumed in our modeling. The authors argue for two ejection dates for the outer ejecta prior to the Great Eruption:
-
1.
In the mid-1200s: The E and NNE condensation originated inthis event. These features are not part of our model and we agreethat they are from a previous mass loss episode.
-
2.
An intermediate eruption in the 16th century: The SE arc, the W condensation, and the NW condensation originated in this event. These three features are also not part of our model. They are certainly distinct to the features that comprise the coherent outer shell structure, which originated in the Great Eruption. However, we assume that the S ridge is associated to the Great Eruption.
Acknowledgments
This research has made use of NASA’s Astrophysics Data System Bibliographic Services, SAOImage DS9 (Joye & Mandel 2003), developed by Smithsonian Astrophysical Observatory, QFitsView (http://www.mpe.mpg.de/~ott/dpuser/qfitsview.html), Aladin sky atlas developed at CDS, Strasbourg Observatory, France (Bonnarel et al. 2000; Boch & Fernique 2014), Astropy, a community-developed core Python package for Astronomy (Astropy Collaboration et al. 2013), APLpy, an open-source plotting package for Python hosted at http://aplpy.github.com, and of Montage, funded by the National Science Foundation under Grant Number ACI-1440620, and was previously funded by the National Aeronautics and Space Administration’s Earth Science Technology Office, Computation Technologies Project, under Cooperative Agreement Number NCC5-626 between NASA and the California Institute of Technology. W.S. acknowledges support from grant IN101014 from UNAM-DGAPA-PAPIIT.
References
- Akashi, M. S., Kashi, A., & Soker, N. 2013, New A, 18, 23 [NASA ADS] [CrossRef] [Google Scholar]
- Allen, D. A., & Hillier, D. J. 1993, Proc. Astron. Soc. Aust., 10, 338 [NASA ADS] [CrossRef] [Google Scholar]
- Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., et al. 2013, A&A, 558, A33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bacon, R., Accardo, M., Adjali, L., et al. 2010, in SPIE Conf. Ser., 7735, 8 [Google Scholar]
- Boch, T., & Fernique, P. 2014, in Astronomical Data Analysis Software and Systems XXIII, eds. N. Manset, & P. Forshay, ASP Conf. Ser., 485, 277 [Google Scholar]
- Bonnarel, F., Fernique, P., Bienaymé, O., et al. 2000, A&AS, 143, 33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Clyne, N., Akras, S., Steffen, W., et al. 2015, A&A, 582, A60 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Corcoran, M. F., Petre, R., Swank, J. H., et al. 1998, ApJ, 494, 381 [NASA ADS] [CrossRef] [Google Scholar]
- Corcoran, M. F., Hamaguchi, K., Gull, T., et al. 2004, ApJ, 613, 381 [NASA ADS] [CrossRef] [Google Scholar]
- Currie, D. G., Dorland, B. N., & Kaufer, A. 2002, A&A, 389, L65 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Davidson, K. 2016, J. Phys. Conf. Ser., 728, 022008 [NASA ADS] [CrossRef] [Google Scholar]
- Davidson, K., & Humphreys, R. M. 1997, ARA&A, 35, 1 [NASA ADS] [CrossRef] [Google Scholar]
- Davidson, K., & Humphreys, R. M. 2012a, Eta Carinae and the Supernova Impostors, 384 [Google Scholar]
- Davidson, K., & Humphreys, R. M. 2012b, Nature, 486, E1 [NASA ADS] [CrossRef] [Google Scholar]
- Davidson, K., Walborn, N. R., & Gull, T. R. 1982, ApJ, 254, L47 [NASA ADS] [CrossRef] [Google Scholar]
- Davidson, K., Dufour, R. J., Walborn, N. R., & Gull, T. R. 1986, ApJ, 305, 867 [NASA ADS] [CrossRef] [Google Scholar]
- Davidson, K., Ebbets, D., Johansson, S., Morse, J. A., & Hamann, F. W. 1997, AJ, 113, 335 [NASA ADS] [CrossRef] [Google Scholar]
- Davidson, K., Gull, T. R., Humphreys, R. M., et al. 1999, AJ, 118, 1777 [NASA ADS] [CrossRef] [Google Scholar]
- Davidson, K., Smith, N., Gull, T. R., Ishibashi, K., & Hillier, D. J. 2001, AJ, 121, 1569 [NASA ADS] [CrossRef] [Google Scholar]
- Dufour, R. J. 1989, Rev. Mex. Astron. Astrofis., 18, 87 [Google Scholar]
- Duschl, W. J., Hofmann, K.-H., Rigaut, F., & Weigelt, G. 1995, in Rev. Mex. Astron. Conf. Ser. 2, eds. V. Niemela, N. Morrell, & A. Feinstein, 17 [Google Scholar]
- Dwarkadas, V. V., & Balick, B. 1998, AJ, 116, 829 [NASA ADS] [CrossRef] [Google Scholar]
- Frank, A., Balick, B., & Davidson, K. 1995, ApJ, 441, L77 [NASA ADS] [CrossRef] [Google Scholar]
- González, R. F., de Gouveia Dal Pino, E. M., Raga, A. C., & Velázquez, P. F. 2004, ApJ, 616, 976 [Google Scholar]
- González, R. F., Villa, A. M., Gómez, G. C., et al. 2010, MNRAS, 402, 1141 [NASA ADS] [CrossRef] [Google Scholar]
- Gull, T. R., Kober, G. V., & Nielsen, K. E. 2006, ApJS, 163, 173 [NASA ADS] [CrossRef] [Google Scholar]
- Hamaguchi, K., Petre, R., Matsumoto, H., et al. 2007, PASJ, 59, 151 [NASA ADS] [Google Scholar]
- Heger, A., & Woosley, S. E. 2002, ApJ, 567, 532 [NASA ADS] [CrossRef] [Google Scholar]
- Humphreys, R. M., & Davidson, K. 1994, PASP, 106, 1025 [NASA ADS] [CrossRef] [Google Scholar]
- Humphreys, R. M., Helton, L. A., & Jones, T. J. 2007, AJ, 133, 2716 [NASA ADS] [CrossRef] [Google Scholar]
- Humphreys, R. M., Davidson, K., & Koppelman, M. 2008, AJ, 135, 1249 [NASA ADS] [CrossRef] [Google Scholar]
- Ishibashi, K., Gull, T. R., Davidson, K., et al. 2003, AJ, 125, 3222 [NASA ADS] [CrossRef] [Google Scholar]
- Joye, W. A., & Mandel, E. 2003, in Astronomical Data Analysis Software and Systems XII, eds. H. E. Payne, R. I. Jedrzejewski, & R. N. Hook, ASP Conf. Ser., 295, 489 [Google Scholar]
- Kashi, A., & Soker, N. 2010, ApJ, 723, 602 [NASA ADS] [CrossRef] [Google Scholar]
- Kiewe, M., Gal-Yam, A., Arcavi, I., et al. 2012, ApJ, 744, 10 [NASA ADS] [CrossRef] [Google Scholar]
- Kiminki, M. M., Reiter, M., & Smith, N. 2016, MNRAS, 463, 845 [NASA ADS] [CrossRef] [Google Scholar]
- Langer, N., García-Segura, G., & Mac Low, M.-M. 1999, ApJ, 520, L49 [NASA ADS] [CrossRef] [Google Scholar]
- Maeder, A., & Desjacques, V. 2001, A&A, 372, L9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Meaburn, J., Gehring, G., Walsh, J. R., et al. 1993, A&A, 276, L21 [NASA ADS] [Google Scholar]
- Meaburn, J., Boumis, P., Walsh, J. R., et al. 1996, MNRAS, 282, 1313 [NASA ADS] [CrossRef] [Google Scholar]
- Moriya, T. J., Maeda, K., Taddia, F., et al. 2014, MNRAS, 439, 2917 [NASA ADS] [CrossRef] [Google Scholar]
- Ofek, E. O., Sullivan, M., Shaviv, N. J., et al. 2014, ApJ, 789, 104 [NASA ADS] [CrossRef] [Google Scholar]
- Owocki, S. P., & Shaviv, N. J. 2016, MNRAS, 462, 345 [NASA ADS] [CrossRef] [Google Scholar]
- Owocki, S. P., Gayley, K. G., & Shaviv, N. J. 2004, ApJ, 616, 525 [NASA ADS] [CrossRef] [Google Scholar]
- Pastorello, A., Smartt, S. J., Mattila, S., et al. 2007, Nature, 447, 829 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
- Portegies Zwart, S. F., & van den Heuvel, E. P. J. 2016, MNRAS, 456, 3401 [NASA ADS] [CrossRef] [Google Scholar]
- Rest, A., Prieto, J. L., Walborn, N. R., et al. 2012, Nature, 482, 375 [NASA ADS] [CrossRef] [Google Scholar]
- Schreier, R., & Soker, N. 2016, Res. Astron. Astrophys., 16, 001 [NASA ADS] [CrossRef] [Google Scholar]
- Seward, F. D., Forman, W. R., Giacconi, R., et al. 1979, ApJ, 234, L55 [NASA ADS] [CrossRef] [Google Scholar]
- Shaviv, N. J. 2000, ApJ, 532, L137 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
- Smith, N. 2004, MNRAS, 351, L15 [NASA ADS] [CrossRef] [Google Scholar]
- Smith, N. 2006, ApJ, 644, 1151 [NASA ADS] [CrossRef] [Google Scholar]
- Smith, N. 2008, Nature, 455, 201 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
- Smith, N. 2014, ARA&A, 52, 487 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
- Smith, N., & Ferland, G. J. 2007, ApJ, 655, 911 [NASA ADS] [CrossRef] [Google Scholar]
- Smith, N., & Frew, D. J. 2011, MNRAS, 415, 2009 [NASA ADS] [CrossRef] [Google Scholar]
- Smith, N., & Morse, J. A. 2004, ApJ, 605, 854 [NASA ADS] [CrossRef] [Google Scholar]
- Smith, N., & Townsend, R. H. D. 2007, ApJ, 666, 967 [NASA ADS] [CrossRef] [Google Scholar]
- Smith, N., Gehrz, R. D., Hinz, P. M., et al. 2003, AJ, 125, 1458 [NASA ADS] [CrossRef] [Google Scholar]
- Smith, N., Morse, J. A., Gull, T. R., et al. 2004, ApJ, 605, 405 [NASA ADS] [CrossRef] [Google Scholar]
- Smith, N., Li, W., Foley, R. J., et al. 2007, ApJ, 666, 1116 [NASA ADS] [CrossRef] [Google Scholar]
- Soker, N. 2005, ApJ, 635, 540 [NASA ADS] [CrossRef] [Google Scholar]
- Steffen, W., García-Segura, G., & Koning, N. 2009, ApJ, 691, 696 [NASA ADS] [CrossRef] [Google Scholar]
- Steffen, W., Koning, N., Wenger, S., Morisset, C., & Magnor, M. 2011, IEEE Transactions on Visualization and Computer Graphics, 17, 454 [Google Scholar]
- Steffen, W., Koning, N., Esquivel, A., et al. 2013, MNRAS, 436, 470 [NASA ADS] [CrossRef] [Google Scholar]
- Steffen, W., Teodoro, M., Madura, T. I., et al. 2014, MNRAS, 442, 3316 [NASA ADS] [CrossRef] [Google Scholar]
- Tiffany, C., Humphreys, R. M., Jones, T. J., & Davidson, K. 2010, AJ, 140, 339 [NASA ADS] [CrossRef] [Google Scholar]
- Tsuboi, Y., Koyama, K., Sakano, M., & Petre, R. 1997, PASJ, 49, 85 [NASA ADS] [Google Scholar]
- van Dyk, S. D. 2005, in The Fate of the Most Massive Stars, eds. R. Humphreys, & K. Stanek, ASP Conf. Ser., 332, 47 [Google Scholar]
- van Dyk, S. D., & Matheson, T. 2012, Eta Carinae and the Supernova Impostors, ASSL, 384, 249 [NASA ADS] [CrossRef] [Google Scholar]
- Walborn, N. R. 1976, ApJ, 204, L17 [NASA ADS] [CrossRef] [Google Scholar]
- Waldman, R. 2008, ApJ, 685, 1103 [NASA ADS] [CrossRef] [Google Scholar]
- Weilbacher, P. M., Streicher, O., Urrutia, T., et al. 2012, in SPIE Conf. Ser., 8451, 84510B [Google Scholar]
- Weis, K. 2012, in Ap&SSL 384, eds. K. Davidson, & R. M. Humphreys, 171 [Google Scholar]
- Weis, K., Duschl, W. J., & Chu, Y.-H. 1999, A&A, 349, 467 [NASA ADS] [Google Scholar]
- Weis, K., Duschl, W. J., & Bomans, D. J. 2001, A&A, 367, 566 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Weis, K., Corcoran, M. F., Bomans, D. J., & Davidson, K. 2004, A&A, 415, 595 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Woosley, S. E., Blinnikov, S., & Heger, A. 2007, Nature, 450, 390 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
- Yoshida, T., Umeda, H., Maeda, K., & Ishii, T. 2016, MNRAS, 457, 351 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: A resized Homunculus model for the outer shell
Figure A.1 shows a SHAPE model of a resized Homunculus geometry for the outer ejecta. This geometry assumption does not reproduce the observations nearly as well as the cylinder geometry shown in Fig. 2. By altering the scale and the speed (i.e., age) of the resized Homunculus model, the match with the observational data becomes only worse than the case shown here.
![]() |
Fig. A.1 Alternative SHAPE model of a resized Homunculus geometry compared to the observed MUSE channel maps of Hβ emission at different velocities in the MUSE field of view of 1′× 1′. This geometry assumption does not reproduce the observations nearly as well as the cylinder geometry shown in Fig. 2. With a given velocity field, there is not a great deal of freedom in the geometry. |
All Figures
![]() |
Fig. 1 HST/ACS F658N filter image, reproduced from Weis et al. (2004). The field of view is 1.5′× 1.5′. Several distinct features are identified with the terminology used by Walborn (1976). The main filter transmission is between λ6549−6623 Å, ranging over 3000 km s-1 in velocity space. Eta Car’s outer ejecta have a large range of velocities and thus [N ii] λλ6548, 6583 and Hα emission contribute to the features observed here. |
In the text |
![]() |
Fig. 2 MUSE channel maps showing the observed Hβ emission in the MUSE field of view of 1′× 1′. Each panel shows a single channel of the MUSE data cube, corresponding to different velocities of the Hβ emission. The velocities are indicated in each panel in units of km s-1 (each channel map spans a range of approximately 80 km s-1, velocities are not corrected for η Car’s systemic heliocentric velocity of approximately −8 km s-1). The corresponding SHAPE model channel map is shown on the right side of each observational channel map (the grayscale only indicates the location of the emission and does not reproduce the brightness distribution), see text for details. The channel maps are irregularly spaced in velocity to show best the brightest emission features corresponding to the outer shell, the ghost shell, the bowshock, and the jet, identified by labels and arrows pointing toward the feature (see also Fig. 4 for their three-dimensional location). The jet is not modeled with SHAPE, since the origin of the emission is not clear, see Sect. 3.4. The central star and part of the Homunculus are saturated. The channel maps are presented with a square-root color stretch and normalized in flux individually to peak intensity. They only contain physical information on the structure and velocities, but do not permit channel-to-channel comparisons of the fluxes to be made. |
In the text |
![]() |
Fig. 3 MUSE channel maps showing the observed Hβ emission in the MUSE field of view of 1′× 1′ as in Fig. 2, but overlayed with the brightness outlines of the model features. Only a few contours are shown to avoid obstructing the underlying observation. |
In the text |
![]() |
Fig. 4 Three-dimensional model of η Car’s outer ejecta seen from Earth (left) and from the side (right) with the bipolar Homunculus nebula in the center (dark red and blue lobes). The model has been constructed with SHAPE and the Hβ emission and velocity (over a range of ± 3000 km s -1) information extracted from the MUSE data. The outer ejecta form a partial bent cylinder. The open cylinder parts are not observed in the MUSE data, i.e., we only trace the bottom of the shell moving toward us and the top of the shell moving away from us. The open cylinder parts are likely too faint, perhaps due to only weak interactions with less dense surrounding material. This large tube is probably bent or indented due to material being slowed down by interactions with the environment. The ghost shell (green) is in our line of sight to the southern Homunculus lobe and consists of two closely aligned shells. An additional half-shell (green mesh, complementing the bottom of the large cylindrical shell moving away from us) lies behind the northern Homunculus lobe. The bowshock-like feature is displayed in yellow. The colors are selected to distinguish the different shapes selected for the SHAPE modeling. |
In the text |
![]() |
Fig. 5 Three-dimensional velocity structure of η Car’s outer ejecta seen from Earth (left) and from the side (right). The color coding is from −2500 km s-1 (blue/violet) to +2500 km s-1 (red). This SHAPE model shows only qualitatively where material is located. The colors of the velocity vectors and the vector lengths show roughly the direction and speed of the material. |
In the text |
![]() |
Fig. 6 Composite images displaying our SHAPE model together with the 0.5–1.2 keV Chandra X-ray image. The soft X-ray emission is produced as material runs into the nearby gas and dust. X-ray emission from the Homunculus is detected only in hard X-ray bands. The large-scale optical tube structure fits right inside the X-ray bubble. The SHAPE model contours only qualitatively show the outlines of the emission. The blue and red colors are to distinguish the different shapes. There is no physical information associated with them. For the image scale, see Fig. 7. |
In the text |
![]() |
Fig. 7 The 0.5–1.2 keV Chandra X-ray image (gray-scale image) overlaid with the Hβ velocity contours obtained from our MUSE data. Contours are fit in individual MUSE slices at an intensity level “Ic” in 10-20 erg/s/cm2/Å and color-coded as a function of velocity. The purple shapes indicate masked-out regions not used in this analysis that include the Homunculus, the artifacts from the optical ghosts described in Sect. 2, and several nearby stars. Because of the number of artifacts present and the faintness of the clumps, two different mask-out regions are used for the positive velocity contours. We also adjust the “Ic” value in each case. Upper left: extent of the soft X-ray emission surrounding η Car. The green box indicates the region shown in the other three panels. Upper right: contours of the blueshifted emission, highlighting the correlation of the optical emission with the soft X-rays in the southeast. The outer shell in the southeast coincides with the inner wall of the X-ray emission. Lower left and right: contours of the redshifted emission. The red part of the outer shell is fainter and it is difficult to determine the velocity contours of the ejecta due to artifacts, but the optical emission here also accurately traces the inner wall of the X-ray emission. |
In the text |
![]() |
Fig. A.1 Alternative SHAPE model of a resized Homunculus geometry compared to the observed MUSE channel maps of Hβ emission at different velocities in the MUSE field of view of 1′× 1′. This geometry assumption does not reproduce the observations nearly as well as the cylinder geometry shown in Fig. 2. With a given velocity field, there is not a great deal of freedom in the geometry. |
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.