A publishing partnership

The following article is Open access

Near-infrared Spectroscopy of Dense Ejecta Knots in the Outer Eastern Area of the Cassiopeia A Supernova Remnant

, , , and

Published 2023 August 10 © 2023. The Author(s). Published by the American Astronomical Society.
, , Citation Bon-Chul Koo et al 2023 ApJ 953 131DOI 10.3847/1538-4357/acda2d

Download Article PDF
DownloadArticle ePub

You need an eReader or compatible software to experience the benefits of the ePub3 file format.

0004-637X/953/2/131

Abstract

The Cassiopeia A supernova remnant has a complex structure, manifesting the multidimensional nature of core-collapse supernova explosions. To further understand this, we carried out near-infrared multiobject spectroscopy on the ejecta knots located in the northeastern (NE) jet and Fe K plume regions, which are two distinct features in the outer eastern area of the remnant. Our study reveals that the knots exhibit varying ratios of [S ii] 1.03, [P ii] 1.189, and [Fe ii] 1.257 μm lines depending on their locations within the remnant, suggesting regional differences in elemental composition. Notably, the knots in the NE jet are mostly S-rich with weak or no [P ii] lines, implying that they originated below the explosive Ne-burning layer, consistent with the results of previous studies. We detected no ejecta knots exhibiting only [Fe ii] lines in the NE jet area that are expected in the jet-driven supernova explosion model. Instead, we discovered a dozen Fe-rich knots in the Fe K plume area. We propose that they are dense knots produced by a complete Si burning with α-rich freeze-out in the innermost region of the progenitor and ejected with the diffuse X-ray-emitting Fe ejecta but decoupled after crossing the reverse shock. In addition to these metal-rich ejecta knots, several knots emitting only He i 1.083 μm lines were detected, and their origin remains unclear. We also detected three extended H emission features of circumstellar or interstellar origin in this area and discuss their association with the supernova remnant.

Export citation and abstractBibTeXRIS

Original content from this work may be used under the terms of the Creative Commons Attribution 4.0 licence. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.

1. Introduction

The core-collapse supernova (SN) explosion is one of the outstanding problems in modern astrophysics. It is generally agreed that, for the majority of core-collapse SNe, the explosion is initiated by a core-bounce shock that becomes stalled but is revived by neutrino heating. Modern state-of-the-art 3D simulations show that this neutrino-driven SN explosion is basically multidimensional, facilitated by turbulent convection and hydrodynamic instabilities (e.g., see Janka 2017; Burrows & Vartanyan 2021). The nonuniform, inhomogeneous, and radially mixed/overturned distribution of metal-rich ejecta in SNe and young SN remnants (SNRs) provides strong observational evidence for it.

Cassiopeia A (Cas A) is a unique SNR where we can see the details of the spatial and velocity distribution of SN ejecta. It is young (∼350 yr; Thorstensen et al. 2001; Fesen et al. 2006) and nearby (3.4 kpc; Reed et al. 1995; Alarie et al. 2014). It is a remnant of a Type IIb SN with a probable progenitor mass of 15–25 M (Young et al. 2006; Krause et al. 2008; see also Koo & Park 2017 and references therein). The remnant is of shell type and has the basic appearance of a small, clumpy bright ring surrounded by a limb-brightened faint plateau. The bright ring has a radius of (1.7 pc) and is mainly composed of C- and O-burning ejecta material heated by reverse shock. The faint plateau extends to (2.5 pc), and its outer boundary defines the current position of the SN blast wave.

In detail, Cas A has a quite complex structure, manifesting the violent and asymmetric explosion. In the outer eastern area of the SNR, which is the area of primary interest in this study, there are two features distinguishable from the rest of the ejecta material. One is the "jet" structure in the northeastern (NE) area, and the other is the Fe K–bright "plume" (hereafter the NE jet and Fe K plume, respectively). The NE jet is a stream of compact optical knots confined into a narrow cone well outside the SN blast wave. It appears to have punched through the main ejecta shell, and the inferred ejection velocity is as large as 16,000 km s−1, which is three times faster than the bulk of the ejecta in the main ejecta shell (Kamper & van den Bergh 1976; Fesen & Gunderson 1996; Fesen 2001; Fesen & Milisavljevic 2016). The optical knots show strong O, S, and Ar emission lines (van den Bergh 1971; Fesen & Gunderson 1996; Fesen 2001; Fesen et al. 2006; Hammell & Fesen 2008; Fesen & Milisavljevic 2016). The NE jet structure was also observed by Chandra, which showed that the X-ray-emitting jet material is rich in O-burning and incomplete Si-burning materials (Si, S, Ar, and Ca) but poor in C-/Ne-burning products (O, Ne, and Mg; e.g., Hughes et al. 2000; Hwang et al. 2004; Vink 2004). The chemical composition of the jet material and the jet morphology suggest that the NE jet originated from a Si–S–Ar-rich layer deep inside the progenitor and penetrated through the outer N- and He-rich envelope (Fesen 2001; Fesen et al. 2006; Laming et al. 2006; Milisavljevic & Fesen 2013). The other feature is the Fe K plume detected by Chandra and XMM-Newton in the eastern area (Willingale et al. 2002; Hwang & Laming 2003; Hwang et al. 2004; Hwang & Laming 2012; Tsuchioka et al. 2022). This plume of hot gas, bright in the Fe K line, extends beyond the main ejecta shell, bright in the Si K and S K lines, and reaches the forward shock front. The velocity of the ejecta material in the outermost Fe K plume is >4500 km s−1, much higher than that of Si-/O-rich ejecta (Tsuchioka et al. 2022). The ejecta material in the Fe K plume has no S or Ar, and its Fe/Si abundance ratio is up to an order of magnitude higher than that expected in the incomplete Si-burning layer, indicating that it originated from a complete Si-burning process with α-rich freeze-out, where the burning products are almost exclusively 56Fe (Thielemann et al. 1996; Hwang & Laming 2003; see also Sato et al. 2021). DeLaney et al. (2010) proposed a model for the Fe K plume where the Fe ejecta pushed through the Si layer as a "piston," while Hwang & Laming (2003) and Milisavljevic & Fesen (2013) suggested that 56Ni bubbles could produce plumelike structures in Fe K emission, but the ionization age of the Fe K plume is inconsistent with a Ni bubble origin. We will explore the properties of the NE jet and Fe K plume in more detail in Section 5, where we discuss the results of our work in this paper. Here we only emphasize that the distribution of Fe ejecta has not been fully explored. The Fe K plume revealed by X-ray observations is hot and "diffuse" Fe-rich gas of density ≲10 cm−3 (Hwang & Laming 2003). Cooler (∼104 K) and presumably denser Fe-rich ejecta has not been explored. The presence or absence of Fe-rich ejecta material in the NE jet is also an important issue related to its origin, but it could not have been addressed in the previous optical studies (Fesen & Gunderson 1996; Fesen 2001; Fesen & Milisavljevic 2016).

Recently, Koo et al. (2018) obtained a long-exposure (∼10 hr) image of Cas A using the UKIRT 4 m telescope with a narrowband filter centered at [Fe ii] 1.644 μm emission. The band includes the [S i] 1.645 μm line, but it is faint and nonnegligible only for the unshocked ejecta in the inner area of the SNR (Koo et al. 2018; Raymond et al. 2018), so, in this paper, we refer to the image as the "deep [Fe ii] image." The [Fe ii] 1.644 μm line is one of the brightest near-infrared (NIR) lines in shocked gas, so it traces dense shocked material (e.g., Koo et al. 2016). The deep [Fe ii] image revealed numerous knots in the outer eastern area beyond the main ejecta shell, the proper motion of which indicates that they are mostly SN ejecta material. A comparison with a Hubble Space Telescope (HST) WFC3/IR F098M image, which is dominated by sulfur lines (Hammell & Fesen 2008; Fesen & Milisavljevic 2016), showed that the majority of the knots in the image correspond to S-rich ejecta knots in optical studies. But a considerable fraction (21%) of the knots does not have optical counterparts, and the knots around the Fe K plume in the eastern area appear to have high Fe abundance (Koo et al. 2018). Hence, the study of these knots could provide new insights on the origin of the NE jet and Fe K plume, as well as the explosion dynamics of the Cas A SN.

In this work, we present the results of NIR spectroscopic observations of the knots in the deep [Fe ii] image in the outer eastern area. In the NIR band (0.95–1.8 μm), in addition to the familiar He i 1.083 μm line and the forbidden lines of S ii, S iii, and Fe ii, there is the [P ii] 1.089 μm (1 D23 P2) line, which is the major NIR line of P ii. Phosphorus (31P) is an uncommon element with a cosmic abundance relative to H and Fe of 2.8 × 10−7 and 8.1 × 10−3 by number, respectively (Asplund et al. 2009). In young SNRs such as Cas A, where we can observe SN ejecta, however, it could be a major element. Phosphorus is produced in hydrostatic neon-burning shells in the pre-SN stage and also in explosive carbon- and neon-burning layers during SN explosion (Arnett 1996; Woosley et al. 2002). Therefore, the [P ii] 1.089 μm line could be strong in some ejecta material of young SNRs, and, together with the other strong lines of α elements (e.g., O, S, and Fe), it can be used for the study of the SN explosion dynamics (Koo et al. 2013).

The organization of the paper is as follows. In Section 2, we describe our NIR observations and explain the data reduction procedures. In Section 3, we explain how we obtain the spectral parameters of the knots for analysis. We identify individual knots in 2D spectra and extract their 1D spectra. We then derive the parameters of individual emission lines by performing single Gaussian fits. The derived fluxes are extinction-corrected using the Herschel SPIRE 250 μm data. In Section 4, we analyze the spectral properties of the ejecta knots. We find that there is significant variation in the ratios of bright lines among the knots in different areas. We classify the knots into different groups based on the bright line ratios and compare their spectral properties. We also search for optical counterparts of the knots and present the results. In Section 5, we discuss the origins of the knots and their connection to the explosion dynamics in Cas A. We focus on the knots in the NE jet and Fe K plume regions. We also discuss the properties of the extended H emission features detected in the outer eastern area and their association with the Cas A SNR. Finally, in Section 6, we summarize the paper.

2. Observations and Data Reduction

We carried out NIR spectroscopic observations of the SN ejecta knots in the outer eastern area of Cas A in 2017 September using the MMT/Magellan Infrared Spectrograph (MMIRS) attached to the MMT 6.5 m telescope. The MMIRS is equipped with a HAWAII-2RG 2048 × 2048 pixel HgCdTe detector with a pixel scale of 0farcs2012 pixel–1, which provides a 6farcm9 × 6farcm9 field of view for imaging observations (McLeod et al. 2012). In order to take J- and H-band spectra of the ejecta knots, we utilized the multiobject spectroscopy (MOS) mode with the J + zJ and H3000 + H (grism+filter) configurations that provide J-band (0.95–1.5 μm) and H-band (1.50–1.80 μm) spectra with a moderate spectral resolution.

In order to take the NIR spectra of the ejecta knots, we carefully selected the bright knots in the outer eastern area from the deep [Fe ii] image (Figure 1). Since the [Fe ii] image was taken in 2013, we estimated the positions of the knots in the 2017 epoch assuming that they have been expanding freely from the explosion center during the age of Cas A. We prepared three MOS masks, each of which includes 20, 17, and 15 MOS slits (Figure 1). Each slit has a width of 0farcs8 and a length of 10''. This "moderate" slit width provides a mean spectral resolving power of ∼8 Å for J-band spectra and ∼7 Å for H-band spectra, corresponding to R ∼ 1600 at 1.26 μm and R ∼ 2300 at 1.64 μm, respectively. At all slit positions in Figure 1, we obtained J-band spectra. For the slits in mask 1, we also took the H-band spectra. In order to subtract sky airglow emission lines, we utilized the ABBA observing sequence with dithering lengths of +2farcs6, −2farcs0, −2farcs0, and +2farcs6. The single exposure time per frame was 300 s, while the total effective exposure time per mask by the multiple frames was 40 or 60 minutes for the J-band spectra and 20 minutes for the H-band spectra. For the absolute flux calibration of the spectra, we also took the spectra of nearby A0V stars several times just before or after the target observations. The seeing during the observation was 0farcs7–1farcs2. The detailed observation log is given in Table 1.

Figure 1. Refer to the following caption and surrounding text.

Figure 1. Finding charts for the slit positions. The red bars represent the locations of the 10''-long MOS slits in three masks. The background is the deep [Fe ii] image taken in 2013 September (Koo et al. 2018). The green plus sign marks the SN expansion center, (23h23m27fs77 ± 0fs05, ), obtained by Thorstensen et al. (2001). The blue dashed line represents the P.A. of 78° measured from north to east at the explosion center, which is used to divide the knots into two groups: (a) NE jet area knots and (b) Fe K plume area knots. Note that the slit positions are shifted to the 2013 epoch to match the background image (see the text for more details). The 1D spectra presented in Figure 7 correspond to the four ejecta knots labeled in blue in the figure.

Standard image High-resolution image

Table 1. Log of MOS Spectroscopy

MaskNumber of SlitsGrismFilterP.A.Exposure Time
    (deg)(s)
(1)(2)(3)(4)(5)(6)
Mask 120 J zJ 32300 × 8
  H3000 H 32300 × 4
Mask 217 J zJ 21300 × 12
Mask 315 J zJ 60300 × 12

Note. (1) Mask name, (2) number of MOS slits in each mask, (3) grism, (4) filter, (5) P.A. of slits measured counterclockwise from north to east on the plane of the sky, (6) exposure time per frame (s) × number of dithering.

Download table as:  ASCIITypeset image

For the data reduction, we utilized the MMIRS pipeline written in IDL (Chilingarian et al. 2015). The pipeline starts with nonlinearity correction of the detector and the cosmic-ray removal process, followed by dark subtraction and flat-fielding correction. The sky background including bright OH airglow emission lines was removed by subtracting the dithered frames. Wavelength calibration and distortion correction were done by using bright OH airglow emission lines falling in the J- and H-band spectra. The typical 1σ uncertainty of the wavelength solution is 0.5 Å for the J band and 0.3 Å for the H band, corresponding to a velocity uncertainty of ∼12 km s−1 at 1.26 μm and ∼5 km s−1 at 1.64 μm, respectively. The absolute photometric calibration was performed by applying the correction factor as a function of wavelength, derived by comparing the observed spectra of standard A0V stars with the Kurucz model spectrum; 5 thus, it also corrects telluric absorption. The uncertainty in the absolute photometric calibration has been estimated at ≲20%, while the relative line fluxes in each wave band should be quite robust.

3. Knot/Line Identification and Extinction Correction

We clearly detected a total of 67 knots emitting at least one emission line from the 52 MOS slits: 22 knots in mask 1, 26 knots in mask 2, and 19 knots in mask 3. No emission line was detected in four slits, while multiple velocity components were detected in 17 slits. Among the 67 knots, 35 knots are located in the NE jet area with a position angle (P.A.) of 55°–78° from the explosion center, whereas the remaining 32 knots are distributed in the Fe K plume area with P.A. > 78° (see Figure 1). Table 2 displays the positions of the 52 slits where the knots have been detected, and Table 3 lists the parameters of the identified knots, including their radial velocities (see below).

Table 2. Slit Positions

MaskSlit No. α (J2000) δ (J2000)
Mask 10523:24:00.342+58:51:06.45
 0623:23:59.797+58:50:53.89
 0723:23:57.769+58:50:49.28
 0823:23:54.621+58:50:50.54
 0923:23:55.148+58:50:32.56
 1023:23:52.626+58:50:22.26
 1123:23:48.962+58:50:26.84
 1223:23:52.200+58:49:57.42
 1323:23:50.194+58:49:40.02
 1423:23:49.642+58:49:24.08
 1523:23:49.462+58:49:04.85
 1623:23:48.859+58:48:54.93
 1723:23:48.487+58:48:43.22
 1823:23:45.436+58:48:32.51
 1923:23:43.979+58:48:19.23
 2023:23:42.877+58:48:05.38
 2123:23:38.578+58:48:12.57
 2223:23:43.680+58:47:33.33
 2423:23:38.731+58:47:29.37
 2823:23:37.742+58:46:34.72
Mask 20323:23:52.229+58:51:39.22
 0423:23:56.865+58:51:07.09
 0523:23:56.026+58:50:48.34
 0623:23:55.432+58:50:34.78
 0723:23:49.371+58:50:31.02
 0823:23:50.915+58:50:12.42
 0923:23:51.953+58:49:55.88
 1023:23:50.194+58:49:40.02
 1223:23:49.296+58:49:03.66
 1323:23:48.340+58:48:53.28
 1423:23:49.155+58:48:28.37
 1523:23:48.034+58:48:17.59
 1623:23:46.344+58:48:04.12
 1723:23:41.761+58:48:05.40
 1823:23:44.543+58:47:31.34
 1923:23:40.775+58:47:28.46
 2323:23:42.721+58:46:19.80
Mask 30723:23:59.854+58:50:50.70
 0823:23:58.676+58:50:42.03
 0923:23:57.739+58:50:31.51
 1023:23:55.808+58:50:35.53
 1123:23:56.335+58:50:05.89
 1223:23:54.448+58:50:07.00
 1323:23:51.720+58:50:21.20
 1423:23:52.053+58:49:53.03
 1523:23:48.541+58:50:00.63
 1623:23:50.035+58:49:09.96
 1723:23:48.548+58:48:59.68
 1823:23:46.346+58:48:54.93
 1923:23:40.964+58:49:08.44
 2223:23:38.661+58:48:22.95
 2323:23:37.549+58:48:13.42

Download table as:  ASCIITypeset image

Table 3. Physical Parameters of 67 Ejecta Knots

Knot No.ID α (J2000) δ (J2000)P.A. drad FWHM vrad NH KnotK2018 No. a
    (deg)(arcmin)(km s−1)(km s−1)(1022 cm−2)Type 
(1)(2)(3)(4)(5)(6)(7)(8)(8)(10)
12-07-223:23:49.49+58:50:33.458.33.30277 (8)+2336 (2)1.42S169B
22-04-123:23:56.87+58:51:07.158.64.41292 (34)+420 (15)1.85S182
32-07-123:23:49.37+58:50:31.058.73.27298 (2)+141 (1)1.42S169A
41-11-123:23:48.96+58:50:26.859.33.19318 (6)+272 (3)1.41S168
51-08-123:23:54.62+58:50:50.559.84.02249 (1)+124 (1)1.60S180
61-05-223:24:00.52+58:51:08.761.24.83247 (26)+172 (11)1.75S b
72-05-123:23:56.03+58:50:48.361.54.16214 (16)+428 (7)1.74S181
81-05-123:24:00.34+58:51:06.561.54.79264 (19)+748 (9)1.72S184
91-07-123:23:57.77+58:50:49.362.74.37297 (17)+832 (10)1.88S200
101-06-223:24:00.19+58:50:58.862.74.72174 (11)+452 (4)1.77S
111-06-123:23:59.80+58:50:53.963.34.63286 (3)+10 (2)1.73S207
123-13-123:23:51.72+58:50:21.263.73.46248 (5)+429 (3)1.55S188
132-06-123:23:55.43+58:50:34.863.83.99339 (7)+541 (4)1.58S195A
143-10-123:23:55.81+58:50:35.564.04.04369 (9)+519 (4)1.61S195B
153-07-123:23:59.85+58:50:50.764.04.62311 (7)+110 (4)1.73S206
161-09-123:23:55.15+58:50:32.664.13.94227 (8)+483 (3)1.58S194
171-10-123:23:52.62+58:50:22.364.33.57283 (3)+490 (2)1.66S189
183-08-123:23:58.68+58:50:42.064.84.42245 (28)+625 (10)1.73S201
192-08-123:23:50.92+58:50:12.465.23.30421 (61)+799 (26)1.49S187A
202-08-223:23:50.90+58:50:12.065.23.29271 (39)−230 (17)1.49S187B
213-15-223:23:49.10+58:50:03.166.03.02312 (28)+329 (14)1.46S b
223-15-123:23:48.54+58:50:00.666.12.94441 (16)+2485 (8)1.46S
233-09-123:23:57.74+58:50:31.566.34.24227 (10)−234 (4)1.66Fe197
243-12-123:23:54.45+58:50:07.069.43.69302 (6)−696 (3)1.64S220
251-12-123:23:52.20+58:49:57.470.23.36252 (3)−600 (2)1.65S214A
261-12-323:23:52.16+58:49:56.970.33.35241 (15)+3242 (8)1.65He
272-09-123:23:51.95+58:49:55.970.53.32317 (6)−515 (2)1.65S214B
283-11-123:23:56.34+58:50:05.970.93.91304 (11)−1199 (5)1.62S224
291-12-223:23:51.87+58:49:53.371.13.30256 (11)−2116 (6)1.65S b
303-11-223:23:55.84+58:50:03.771.13.84247 (9)+752 (4)1.62S222
313-14-123:23:52.05+58:49:53.071.33.32281 (3)−2125 (2)1.65S212
323-14-223:23:51.78+58:49:51.871.43.28217 (8)−3170 (4)1.62S b
332-09-223:23:51.72+58:49:51.271.63.27249 (9)−3152 (3)1.62He
341-13-123:23:50.19+58:49:40.073.83.02250 (8)−2017 (4)1.49S226
351-14-223:23:50.05+58:49:29.277.02.96187 (10)+2314 (6)1.48He
363-19-123:23:40.96+58:49:08.479.51.74358 (35)−1954 (14)1.57S229
373-19-223:23:40.90+58:49:08.179.61.73286 (22)−2411 (9)1.57S
383-16-223:23:49.86+58:49:09.283.42.88211 (14)−3312 (6)1.47He
391-15-123:23:49.46+58:49:04.984.72.82247 (13)−1286 (5)1.42Fe235
402-12-123:23:49.30+58:49:03.785.12.80268 (13)−1282 (5)1.42Fe234
412-12-223:23:49.23+58:49:02.485.52.79271 (23)−269 (11)1.42He
423-17-123:23:48.55+58:48:59.786.32.70210 (17)−1616 (8)1.39Fe232
432-13-323:23:48.58+58:48:58.186.92.70259 (17)−1987 (5)1.39Fe
443-18-123:23:46.35+58:48:54.987.82.41261 (14)−1319 (6)1.47S238
452-13-223:23:48.42+58:48:55.088.02.68235 (8)−1165 (3)1.37Fe
461-16-123:23:48.86+58:48:54.988.02.73242 (3)−1189 (2)1.39Fe241
472-13-123:23:48.34+58:48:53.388.62.66328 (4)−1159 (2)1.37Fe239
481-17-123:23:48.49+58:48:43.292.22.68161 (14)−903 (6)1.37Fe240
492-14-223:23:49.34+58:48:32.195.92.81227 (17)−1554 (9)1.34Fe251
501-18-123:23:45.44+58:48:32.597.02.30273 (3)+755 (2)1.44S244
512-14-123:23:49.15+58:48:28.497.22.79257 (6)−207 (2)1.32Fe250
522-14-323:23:49.14+58:48:28.297.22.79246 (2)+765 (1)1.32He
532-15-123:23:48.03+58:48:17.6101.42.68303 (10)−644 (4)1.36Fe259
542-15-223:23:47.96+58:48:16.1102.02.67257 (2)+1187 (1)1.36He
551-19-123:23:43.98+58:48:19.2103.42.16349 (9)+753 (5)1.41S270
562-16-123:23:46.34+58:48:04.1107.42.52230 (11)−1413 (6)1.41Fe272
572-17-223:23:41.93+58:48:08.8110.31.95280 (15)−2052 (7)1.43S
581-20-123:23:42.88+58:48:05.4110.52.09290 (7)−1651 (4)1.39S277
592-17-123:23:41.76+58:48:05.4112.01.95320 (35)−531 (14)1.42S b 275
602-17-323:23:41.70+58:48:04.3112.61.95405 (37)+560 (18)1.46S
613-23-123:23:37.55+58:48:13.4115.31.40375 (5)−674 (3)1.68S282
623-23-223:23:37.33+58:48:12.4116.51.38300 (8)−4004 (5)1.69S281
632-18-123:23:44.54+58:47:31.3120.92.53280 (22)−2023 (10)1.48Fe290
641-22-123:23:43.68+58:47:33.3121.62.42274 (40)−1725 (21)1.47Fe288
652-19-123:23:40.78+58:47:28.5128.72.16339 (13)−44 (7)1.42S293
661-24-123:23:38.73+58:47:29.4133.21.95436 (24)+46 (13)1.43S296
672-23-123:23:42.72+58:46:19.8142.13.16197 (3)+1957 (2)1.44Fe299

Notes. (1) Knot serial number in this work. (2) Mask-slit-number, where the number is the serial number of a knot detected in the slit. (3) and (4) Coordinates of the emission feature in the slit. (5) P.A. measured counterclockwise from north to east at the center of the SN explosion. (6) Radial distance from the center of the SN explosion. (7) and (8) FWHM and central LSR velocity of the knot obtained from a weighted average of the emission lines, where the numbers in parentheses are the 1σ statistical uncertainties from single Gaussian fittings (see Section 2). (9) Hydrogen column density derived from a Herschel SPIRE 250 μm image. (10) Knot type: S = S-rich, Fe = Fe-rich, He = He i. See Section 4.2 for an explanation of the classification of the knots. (11) Knot serial number in Koo et al. (2018).

a When there are multiple velocity components, they are labeled "A" and "B." b The [S ii] lines are detected, but [Fe ii] lines are not detected.

Download table as:  ASCIITypeset images: 1 2

We identified 29 emission lines in total in the J and H bands, including forbidden lines from C, Si, P, S, and Fe species and the He i 1.083 μm line. All 67 knots show at least one of the four strongest emission lines: He i 1.083 μm, the [S ii] 1.03 μm multiplet, [S iii] 0.983 μm, and [Fe ii] 1.257 μm. About one-third (27) of them also show a strong [P ii] 1.189 μm line, and they always accompany strong [S ii], [S iii], and [Fe ii] lines. Seven knots show only a He i 1.083 μm line with or without a very weak [C i] 0.985 μm line above a 3σ–4σ rms noise level, while five knots show only the [S iii] 0.953 and/or [S ii] 1.03 μm lines (see Table 3). Five out of 22 knots in mask 1 show the [Si i] 1.646 μm line in the H band together with bright [S ii], [S iii], [Fe ii], and [P ii] lines. Previous NIR studies reported the detection of the [P ii] 1.147 and [Si x] 1.430 μm lines in the main ejecta shell (Gerardy & Fesen 2001; Lee et al. 2017). In our observations, the [P ii] 1.147 μm line is hardly visible from individual spectra, partly because the line is located at the part of the spectra with higher noise. We do barely see the line, however, if we stack the spectra of the knots in the NE jet area. The [Si x] 1.430 μm line mostly falls outside the detector area.

We extracted 1D spectra of the 67 knots by averaging the spectra of a given region and performed single Gaussian fitting for all the detected emission lines to derive three physical parameters of the individual lines: central velocities, line widths, and line fluxes. Some emission lines, however, are very close or even blended with each other, so that the single Gaussian fitting could not be done without some additional constraints. Around 1.08 μm, for example, there are two strong emission lines from different species, [S i] 1.082 and He i 1.083 μm, and their wavelength difference is only 9 Å. The spectral resolving power in the J band is 8 Å, so in most cases, we could distinguish whether the line is [S i] or He i from a single Gaussian fitting. But when the lines are broad (FWHM > 16 Å), the two lines are blended, and we had to perform double Gaussian fitting with the central wavelengths and widths of the lines fixed. For the [S ii] 1.03 μm multiplets at 1.029, 1.032, 1.034, and 1.037 μm, we also fixed their wavelengths and line widths, as well as their flux ratios based on Tayal & Zatsarinny (2010). Table 4 lists the detected lines in the individual knots and their derived parameters. For the undetected emission lines, we estimated 1σ upper limits to their fluxes using the rms noise level at their expected wavelengths. All of the measured radial velocities have been corrected to the heliocentric reference frame.

Table 4. Observed NIR Line Parameters of 67 Ejecta Knots

KnotLine ID λrest Observed Flux a
No.Transition (lu)(μm)(10−17 erg s−1 cm−2)
(1)(2)(3)(4)
1[S iii] 3 P21 D2 0.9533721.9 (2.8)
1[C i] 3 P11 D2 0.98268·(2.4)
1[C i] 3 P21 D2 0.9853011.5 (1.7)
1[S ii] 2 D3/22 P3/2 1.0289615.0 (0.8)
1[S ii] 2 D5/22 P3/2 1.0323318.2 (0.0)
1[S ii] 2 D3/22 P1/2 1.0339211.3 (1.3)
1[S ii] 2 D5/22 P1/2 1.037334.9 (0.0)
1[S i] 3 P21 D2 1.0824114.8 (2.0)
1He i 3 S13 P0,1,2 1.08332257.0 (5.0)
1[S i] 3 P11 D2 1.13089·(0.7)
1[P ii] 3 P21 D2 1.18861·(1.0)
1[Fe ii] a6 D7/2a4 D5/2 1.24888·(21.3)
1[Fe ii] a6 D3/2a4 D1/2 1.25248·(2.4)
1[Fe ii] a6 D9/2a4 D7/2 1.257025.9 (1.8)
1[Fe ii] a6 D1/2a4 D1/2 1.27069·(3.3)
1[Fe ii] a6 D3/2a4 D3/2 1.27913·(2.1)
1[Fe ii] a6 D5/2a4 D5/2 1.29462·(2.4)
1[Fe ii] a6 D1/2a4 D3/2 1.29813·(4.5)
1[Fe ii] a6 D7/2a4 D7/2 1.32092·(1.2)
1[Fe ii] a6 D3/2a4 D5/2 1.32814·(0.9)
1[Fe ii] a4 F9/2a4 D5/2 1.53389·(·)
1[Fe ii] a4 F7/2a4 D3/2 1.59991·(·)
1[Si i] 3 P11 D2 1.60727·(·)
1[Fe ii] a4 F9/2a4 D7/2 1.64400·(·)
1[Si i] 3 P21 D2 1.64590·(·)
1[Fe ii] a4 F5/2a4 D1/2 1.66422·(·)
1[Fe ii] a4 F7/2a4 D5/2 1.67733·(·)
1[Fe ii] a4 F5/2a4 D3/2 1.71159·(·)
1[Fe ii] a4 F3/2a4 D1/2 1.74541·(·)

Note. (1) Knot serial number in this work, (2) line ID and its transition, (3) rest wavelength in vacuum, (4) observed line flux. The numbers in parentheses are the 1σ statistical uncertainties from single Gaussian fittings (see Section 2).

a If the line falls within the observed wavelength range but not detected, its flux is represented as ".(σ)" where σ is the rms uncertainty at the wavelength of the line. If the line is beyond the observed wavelength range, e.g., the lines in the H band for the knots in masks 2 or 3, the flux is represented as ".(.)".

Only a portion of this table is shown here to demonstrate its form and content. A machine-readable version of the full table is available.

Download table as:  DataTypeset image

For an analysis of line flux ratios, we need to apply an extinction correction to the observed fluxes. The extinction to Cas A is large and varies significantly over the field, e.g., AV = 5–15 (Eriksen et al. 2009; Lee et al. 2015; Koo et al. 2018, and references therein). A column density map of the foreground gas/dust across the Cas A SNR had been obtained from X-ray spectral analysis (Hwang & Laming 2012), but our target knots are mostly outside this column density map. For the knots with both [Fe ii] 1.26 and 1.64 μm lines detected, we can derive extinction by using the flux ratio of these two lines because they share the same upper level (see, e.g., Koo & Lee 2015; Lee et al. 2017). However, we only obtained H-band spectra for the knots in mask 1, and some knots do not show any detectable [Fe ii] lines. So, for the majority of the knots, we could not use this technique either. We instead used the Herschel SPIRE 250 μm data for the extinction correction. The extinction toward Cas A is mostly due to the interstellar medium (ISM) in the Perseus spiral arm, so that the dust emission at 250 μm is a good measure of the extinction to Cas A (e.g., De Looze et al. 2017; Zhou et al. 2018). Koo et al. (2018) showed that there is a good correlation between the Herschel SPIRE 250 μm surface brightness (S250 μm) and the X-ray-absorbing column density (NH,X-ray) toward dense circumstellar knots scattered over the remnant. Their plot is reproduced in Figure 2, where we see that the correlation is linear at NH,X-ray ≲ 2 × 1022 cm−2 and consistent with the 250 μm emission being from 20 K dust with the general ISM dust opacity, i.e., S250 μm/NH,X-ray ≈ 1.1 MJy/1020 cm−2. For higher column densities, S250 μm are considerably smaller than those expected from NH,X-ray, which could be either because some column densities are due to molecular clouds and the temperatures of dust associated with molecular clouds are lower, or because some column densities are due to H gas without dust (Hartmann et al. 1997; Koo et al. 2018). In the figure, we overplot our target knots where both [Fe ii] 1.26 and 1.64 μm lines have been detected with a high (>6) signal-to-noise ratio. There are 11 such knots, and their [Fe ii] 1.26/1.64 μm line ratios range from 0.57 to 0.75, corresponding to hydrogen column densities of NH = 1.2–1.8 × 1022 cm−2 adopting a theoretical line ratio of 1.36, which is the value suggested by Nussbaumer & Storey (1988) and Deb & Hibbert (2010), and the dust opacity of the general ISM (Draine 2003). We see that the column densities derived from the [Fe ii] line ratios are generally consistent with the 250 μm brightness, although there is a considerable uncertainty in the former; i.e., the Einstein A-coefficients of the two lines are uncertain, and the [Fe ii] 1.26/1.64 μm line ratios in the literature range from 0.98 to 1.49 (see Giannini et al. 2015; Koo & Lee 2015, and references therein). In this work, we have derived the absorbing column densities NH for the 67 knots using the Herschel SPIRE 250 μm image as described above and applied an extinction correction to the observed line fluxes by using the extinction curve of the general ISM (Draine 2003). The derived column densities for individual knots are listed in Table 3. Note that the uncertainty in the extinction correction will not significantly affect the results of our analysis using the J-band line ratios; e.g., an error of 5 × 1021 cm2 in NH would yield an error of ∼10% in the extinction-corrected ratios of the [S ii] 1.03 and [Fe ii] 1.257 μm line fluxes.

Figure 2. Refer to the following caption and surrounding text.

Figure 2. Hydrogen column density vs. 250 μm brightness toward the Cas A SNR. The NH is derived from either NIR [Fe ii] line ratios (black open circles; Section 3) or X-ray analysis (gray dots; Hwang & Laming 2012). The S250 μm is obtained from the Herschel SPIRE 250 μm image after subtracting the synchrotron continuum. The dashed lines show the expected S250 μm of the ISM of NH with dust at Td = 18, 20, and 22 K. The right y-axis shows the visual extinction (A V ) corresponding to NH for the general ISM, i.e., A V /NH = 1.87 × 1021 mag cm−2 (Draine 2003).

Standard image High-resolution image

In addition to the 67 knots, we also detected "extended H emission features" in three slits (2–14, 2–15, and 3–18) in the eastern outer region with a P.A. of ∼90°. They only show hydrogen lines, without any metallic emission lines. They will be discussed in Section 5.4, where we show that they are most likely arising from the ISM/circumstellar medium (CSM) rather than the SN ejecta.

4. Results

4.1. Line Ratios and Their Regional Characteristics

Figure 3 shows the locations of the observed knots and line flux ratios of three bright emission lines (i.e., the [S ii] 1.03, [P ii] 1.189, and He i 1.083 μm lines) to the [Fe ii] 1.257 μm line, hereafter R([S ii]/[Fe ii]), R([P ii]/[Fe ii]), and R(He i/[Fe ii]). In the left panels, the line flux ratios are plotted as a function of the knot's radial distance from the explosion center, while in the right panels, the locations of the knots are shown in the (R.A., decl.) plane. In all panels, the color of the symbols represents the flux ratios, i.e., red (blue) means low (high) ratios. The squares and circles indicate the knots in the NE jet and Fe K plume areas, respectively. The red and green dashed lines in the right panels indicate the nominal location of the forward and reverse shocks, respectively. We can see that the NE jet knots are located outside of the forward shock, while most of the knots in the Fe K plume area are located between the forward and reverse shocks, except for one at P.A. = 142°, which is apparently located outside the forward shock.

Figure 3. Refer to the following caption and surrounding text.

Figure 3. Identified NIR ejecta knots and their line ratios. The rows represent different line ratios: R([S ii]/[Fe ii]) at the top, R([P ii]/[Fe ii]) in the middle, and R(He i/[Fe ii]) at the bottom. The left column displays the line ratio against the knot's radial distance from the explosion center, while the right column shows the locations of the knots in the (R.A., decl.) plane. In all panels, the squares and circles indicate the NE jet and Fe K plume area knots, respectively, and their colors represent the flux ratios as in the color bars. The red and green dashed lines in the right panels represent the nominal location of the forward and reverse shocks at the 2017 epoch, and , respectively (Gotthelf et al. 2001; DeLaney & Rudnick 2003).

Standard image High-resolution image

The ratio R([S ii]/[Fe ii]) of the ejecta knots varies over almost 2 orders of magnitude, i.e., 0.2–70 (see top row of Figure 3). A noticeable feature is a group of knots with low R([S ii]/[Fe ii]) at . Most of them have a strong [Fe ii] 1.257 μm line without any detectable [S ii] 1.03 μm lines, so they are marked with upper limits in the figure. They are knots located around the forward shock front just outside the Fe K plume. Their R([S ii]/[Fe ii]) upper limits are more than 1 or 2 orders of magnitude lower than the R([S ii]/[Fe ii]) of the knots around the reverse shock at . There also appears to be a systematic variation in the R([S ii]/[Fe ii]) of the knots at r ≲ 2.3; i.e., the knots at larger radial distances have smaller ratios. The knots at have ratios of a few tens, while those at have ratios of ∼10. On the other hand, the NE jet knots, except for one, have R([S ii]/[Fe ii]) > 4, which is comparable with the knots around the reverse shock. There is no systematic variation in the R([S ii]/[Fe ii]) of the NE jet knots.

The ratio R([P ii]/[Fe ii]) shows a similar pattern to the R([S ii]/[Fe ii]); i.e., the knots at have a very low R([P ii]/[Fe ii]), and the knots at smaller radial distances have higher ratios (see middle row of Figure 3). Again, most of the knots around have no detectable [P ii] 1.189 μm line. Most of the NE jet knots also have low R([P ii]/[Fe ii]), substantially lower than the knots around the reverse shock. In particular, most of the knots at do not have a detectable [P ii] 1.188 μm line. The 2σ upper limits of these knots are 0.05–0.3, which are a factor of 3–15 smaller than the ratios of the knots around the reverse shock.

The ratio R(He i/[Fe ii]) between the knots differs by a factor of 103 (see bottom row of Figure 3). For most knots located at , the He i 1.083 μm line is not detected, and the upper limit of R(He i/[Fe ii]) is 0.1–1.0. However, knots of high ratio abruptly appear at , and the ratio reaches up to 100. Several knots in the NE jet area near also show exceptionally high R(He i/[Fe ii]). Seven knots were identified that exhibit only the He i 1.083 μm line with/without the very weak [C i] 0.985 μm line. These knots will be referred to as "He i knots" for the rest of the paper. These He i knots are detected mainly along the forward shock from the southern base of the NE jet to the Fe K plume area. It is worth mentioning that the radial velocities of most of the He i knots are very large, so they are not circumstellar material (see Section 5.2). In the NE jet, the knots near the tip of the jet have low (≲1.0) R(He i/[Fe ii]).

4.2. Classification of Knots and Their Physical Conditions

4.2.1. Knot Classification

In the previous section, we have found that the ejecta knots in different areas show distinct NIR spectroscopic properties. In this section, we classify the knots using three line ratios: R([S ii]/[Fe ii]), R([P ii]/[Fe ii]), and R(He i/[Fe ii]). We use the [S ii] 1.03 μm multiplet line instead of the [S iii] 0.953 μm line. The two lines have a good correlation, and the latter is usually brighter than the [S ii] 1.03 μm multiplet, i.e., F([S iii] 0.953)/F([S ii] 1.03) = 1.5 ± 0.6. But the [S iii] 0.953 μm line is very close to the blue limit of the spectroscope throughput, so the line becomes hardly detectable if it gets significantly blueshifted (e.g., if its radial velocity is less than −2000 km s−1). There are two knots emitting only weak [S iii] 0.953 μm lines, but the expected [S ii] 1.03 μm fluxes of those two knots from F([S iii] 0.953)/F([S ii] 1.03) = 1.5 are smaller than their 3σ rms noise levels, indicating that the nondetection of the [S ii] line is probably due to the low sensitivity of our observations. For these two knots, we use the expected [S ii] 1.03 μm flux for their classification for convenience.

Figure 4 is a diagram of R([P ii]/[Fe ii]) versus R([S ii]/[Fe ii]). The NE jet and Fe K plume area knots are marked by squares and circles, respectively. We first examine the R([S ii]/[Fe ii]) of the knots. The R([S ii]/[Fe ii]) of the knots ranges from 0.2 to 70. The Fe K plume area knots are scattered over the entire range, while the NE jet knots have relatively high ratios, i.e., R([S ii]/[Fe ii]) ≳ 4. Note that most of the Fe K plume area knots with R([S ii]/[Fe ii]) ≲ 4 do not have detectable [S ii] emission. Such separation of the knots in R([S ii]/[Fe ii]) is similar to the properties of the knots in the main ejecta shell. For the knots in the main ejecta shell, Lee et al. (2017) performed principal component analysis of their NIR spectral properties and showed that they could be classified into two groups: (1) S-rich knots showing strong [S ii] and [P ii] lines and (2) Fe-rich knots showing strong [Fe ii] lines. They found that the two groups were well divided by the ratio F([S ii] 1.03)/F([Fe ii] 1.644) of ∼5. Following Lee et al. (2017), we adopt the threshold of R([S ii]/[Fe ii]) ∼ 3.7, assuming an intrinsic ratio of 1.36 for F([Fe ii] 1.257)/F([Fe ii] 1.644) (Nussbaumer & Storey 1988; Deb & Hibbert 2010). We therefore divide the knots into two groups: S-rich knots with R([S ii]/[Fe ii]) ≥ 3.7 and Fe-rich knots with R([S ii]/[Fe ii]) < 3.7. In Figure 4, S- and Fe-rich knots are represented by blue and red symbols, respectively. Note that all Fe-rich knots except one are the knots in the Fe K plume area and that almost of all of them have no detectable [S ii] and [P ii] emission. The Fe-rich knots are mostly located around the outer boundary of the Fe K plume (see Figure 3), and, as we will discuss in Section 5.2, this has an interesting implication for the SN explosion dynamics. Another thing to note is that all NE jet knots except one are S-rich knots and that they have R([P ii]/[Fe ii]) ratios that are generally lower than those of the S-rich knots in the Fe K plume area (see also Figure 3). We will explore the spectral properties of the NE jet knots in detail in Section 5.1.

Figure 4. Refer to the following caption and surrounding text.

Figure 4. Classification of the NIR knots in the R([S ii]/[Fe ii]) vs. R([P ii]/[Fe ii]) plane. The squares and circles indicate the knots in the NE jet and Fe K plume areas, respectively. The arrows represent 2σ upper/lower limits. The black dashed line marks the boundary between the S-rich (blue) and Fe-rich (red) knots. The number of knots belonging to the two groups is given in parentheses. The number (44) of S-rich knots includes five knots that could not be shown in the figure because only [S ii] lines are detected. There are also seven He i knots where these metal lines are not detected (see Section 4.2.1).

Standard image High-resolution image

Figure 4 does not include all 67 ejecta knots because some knots do not exhibit [Fe ii] 1.257 μm (and [P ii] 1.189 μm) lines. This may appear strange because our target knots have been selected from the deep [Fe ii] image. But the observation was done with a 10'' slit in ABBA mode, so there could be non-[Fe ii] emission knots accidentally located within the slit. There are 12 knots with no detectable [Fe ii] 1.257 μm lines. Among these, four knots show only [S ii] 1.257 and/or [S iii] 0.953 μm lines, one knot shows only [S ii] 1.257 and He i 1.083 μm lines, and seven knots show only He i 1.083 μm lines. The former five knots belong to S-rich knots because the [Fe ii] 1.257 μm line has not been detected. (Their 2σ lower limits of R([S ii]/[Fe ii]) are ≥3.7.) But R([P ii]/[Fe ii]) is not available for them because both the [Fe ii] 1.257 and [P ii] 1.189 μm lines are not detected. Then there are seven He i knots that show only the He i 1.083 μm line with no [S ii] 1.03, [Fe ii] 1.257, or [P ii] 1.189 μm lines. Their He i 1.083 μm fluxes are comparable to those of the metal-rich ejecta knots, but their [S ii] 1.03, [Fe ii] 1.257, and [P ii] 1.189 μm fluxes are considerably lower than those of the other ejecta knots (Figure 5). We note that these He i knots are mostly an accidental detection and that the category of the He i knots is not very rigorously defined. Indeed, some spectra of S-/Fe-rich knots with high R(He i/[Fe ii]) if degraded could become indistinguishable from the He i knots. However, the number of such knots is less than a few and thus should not affect our discussion of He i knots. From Figure 5, we can also see that S- and Fe-rich knots have a comparable range of [Fe ii] 1.257 μm line strengths. This is not unexpected given that the knots that we observed are the bright ones in the deep [Fe ii] image. The two groups of knots, however, are clearly distinguished in their [S ii] line fluxes. In [P ii] emission, the Fe-rich knots generally have lower 1.189 μm fluxes than the S-rich knots, but not necessarily. A good fraction of the S-rich knots have [P ii] 1.189 μm fluxes as low as the Fe-rich knots. As we will see in Section 5.1, they are the knots near the tip of the NE jet.

Figure 5. Refer to the following caption and surrounding text.

Figure 5.  F([Fe ii]) vs. F([S ii]), F([P ii]), and F(He i). The knot symbols and their red and blue colors are the same as in Figure 4. The green symbols represent He i knots. The arrows represent 2σ upper/lower limits. The dashed line in the left panel marks the boundary between S- and Fe-rich knots. The dashed line in the right panel is the line at which R(He i/[Fe ii]) = 8. The He i knots have R(He i/[Fe ii]) higher than this.

Standard image High-resolution image

In summary, we have classified the 67 knots into three groups: S-rich (44), Fe-rich (16), and He i (7) knots. Figure 6 shows which group each individual observed knots belong to. We see that the majority (31/35) of the NE jet knots are S-rich knots. Only one knot is classified as Fe-rich, which has R([S II]/[Fe II])(=1.8) close to the threshold dividing S- and Fe-rich knots. For comparison, the majority (26/32) of the Fe K plume area knots are S- or Fe-rich with clear spatial separation between them; the S-rich knots are near the main ejecta ring within the SNR, while the Fe-rich knots are aligned in the north–south direction along the boundary of the SNR. Figure 6 also shows that those Fe-rich knots appear to be clustered in the position–velocity diagram, implying that they were ejected in a narrow cone. We will explore the kinematic properties of the Fe-rich knots in Section 5.2. The He i knots are detected along the forward shock front, from the southern base of the NE jet to the Fe K plume area. Most of the He i knots have large radial velocities, so they are different from the circumstellar knots that are bright in the He i 1.083 μm line (see Lee et al. 2017). We will explore their origin in Section 5.2. Table 5 summarizes the properties of the three groups, and Figure 7 shows the representative sample spectra of individual groups.

Figure 6. Refer to the following caption and surrounding text.

Figure 6. Spectral properties of the 67 ejecta knots shown in the (a) (R.A., decl.) and (b) (radial distance, radial velocity) planes. The spectral properties are color-coded: S-rich knots in blue and Fe-rich knots in red. The He i knots are shown in green. The background in panel (a) is a Chandra X-ray 4.2–6.4 keV image produced from the data taken in 2013–2017. The dashed lines and green plus sign are the same as in Figure 3. The gray dashed line in the right panel indicates the systematic velocity centroid of the remnant, vr = +800 km s−1 (Reed et al. 1995; DeLaney et al. 2010; Isensee et al. 2010; Milisavljevic & Fesen 2013).

Standard image High-resolution image
Figure 7. Refer to the following caption and surrounding text.

Figure 7. Sample 1D spectra of three ejecta knot groups: S-rich (2-06-1 and 1-18-1), Fe-rich (1-16-1), and He i (2-15-2). For the S-rich knots, separate sample spectra are shown for the NE jet and Fe K plume areas because the intensities of the [P ii] 1.189 μm line differ significantly between the two regions. The positions of the knots are marked in Figure 1. The spiky features in the spectra are the residuals of bright OH airglow emission lines.

Standard image High-resolution image

Table 5. Spectral Classification of the 67 NIR Ejecta Knots and Their Statistics

Knot GroupSpectral PropertiesNumber of Knots
  NE Jet AreaFe K Plume Area
(1)(2)(3)(4)
S-rich knotsStrong [S ii], [Fe ii] lines3113
 Strong [P ii] line (Fe K plume area); weak/no [P ii] line (NE jet area)  
  R([S ii]/[Fe ii]) > 3.7  
Fe-rich knotsStrong [Fe ii] line; weak or no [S ii], [P ii] lines115
  R([S ii]/[Fe ii]) < 3.7  
He i knotsStrong He i line; weak or no [C i] line34
  R(He i/[Fe ii]) > 8  
Total 3532

Note. (1) Knot group, (2) line ratios, (3) and (4) number of knots in the NE jet and Fe K plume areas.

Download table as:  ASCIITypeset image

4.2.2. Spectral and Physical Properties of Three Knot Groups

In this section, we examine the basic spectral properties of three knot groups and compare the densities of their line-emitting regions, looking for differences other than the bright line ratios.

Figure 8 compares the distributions of the radial velocities and line widths of three knot groups. The radial velocities of the S-rich knots strongly peak near 0 km s−1, which is mostly due to the knots in the NE jet. The radial velocities of the knots in the NE jet are mostly confined to a narrow range of vr = 0 to +800 km s−1, but there are also knots with high positive and negative velocities; e.g., the radial velocities of the knots in the southernmost stream are vr = −2100 to −520 km s−1. This is consistent with the previous result from optical studies that the bright knots in the jet streams lie close to the plane of the sky, although the knots in the NE jet region, including those in the jet base area in general, encompass a broad range of radial velocities (Fesen & Gunderson 1996; Fesen 2001; Milisavljevic & Fesen 2013). For comparison, the S- and Fe-rich knots in the Fe K plume area are mostly blueshifted. As already mentioned in Section 4.2.1, the Fe-rich knots are mostly located in the Fe K plume area, and their radial velocities are confined to a narrow range (−2000 to −200 km s−1). On the other hand, the radial velocities of the He i knots are spread over a broad velocity range, from −3500 to +3500 km s−1. In line width, the S-rich knots have a broad distribution centered at around +290 km s−1, while the Fe- and He-rich knots have a relatively narrow distribution centered at approximately +250 km s−1. The velocity widths may be considered as a characteristic shock speed for the knots. But as we will see in Section 4.3, the ejecta knots appear to have complex structures with subknots, so the shock speed could be substantially lower than this. Figure 8 also compares the distribution of [Fe ii] 1.257 μm line fluxes among the three knot groups. As we mentioned in Section 4.2.1, the S- and Fe-rich knots have a similar distribution of [Fe ii] 1.257 μm line fluxes, while the He i knots have the [Fe ii] 1.257 μm line undetected with an upper limit of mostly a few times 10−16 erg s−1 cm2 (see Figure 3). Note that the [Fe ii] 1.257 μm line flux in the table is the flux within the slit. The fluxes of the knots in the deep [Fe ii] image may be found in Koo et al. (2018).

Figure 8. Refer to the following caption and surrounding text.

Figure 8. Distribution of the spectral parameters of the knots in three knot groups. The columns represent different parameters: radial velocity, line width, and [Fe ii] 1.257 μm flux. The rows represent different knot groups. In each panel, the lighter histogram represents the knots in the NE jet area, and the darker histogram represents the knots in the Fe K plume area.

Standard image High-resolution image

One of the physical parameters of the knots that can be easily obtained is electron density. In the J and H bands, there are several [Fe ii] lines originating from levels with similar excitation energies whose ratios depend mainly on electron densities and can be used as a density tracer (e.g., see Koo et al.2016; Lee et al. 2017). In Figure 9, we plot the F([Fe ii] 1.295)/F([Fe ii] 1.257) of the knots, which is sensitive to electron density between 103 and 105 cm−3. The x-axes of the panels in the figure are the [S ii] 1.03, [P ii] 1.189, and He i 1.083 μm fluxes normalized by the [Fe ii] 1.257 μm flux. For about two-thirds of the knots, the [Fe ii] 1.295 μm line has not been detected, so their upper limits are plotted. For the knots with the [Fe ii] 1.295 μm line detected, the electron densities inferred from their [Fe ii] line ratios are mostly in the range of (1–4) × 104 cm−3, with a mean of 1.8 × 104 cm−3. The S-rich knots have a slightly lower density than the Fe-rich knots, i.e., 1.5 × 104 versus 2.5 × 104 cm−3. Note that this is the density of the shock-compressed, line-emitting region of the knots. The initial density of a knot before being shocked should be much lower (e.g., ∼102 cm−3; see Figure S3 of Koo et al. 2013). There is no clear association between the line ratios of the bright lines and F([Fe ii] 1.295)/F([Fe ii] 1.257), except for a moderate positive correlation (ρ = +0.33) observed for R(He i/[Fe ii]), which could possibly be because the He i 1.083 μm line emissivity increases with density due to the contribution from the collisional excitation from the lower level (23S; e.g., see Koo et al. 2023). For a given F([Fe ii] 1.295)/F([Fe ii] 1.257) ratio, the R([S ii]/[Fe ii]) of S- and Fe-rich knots differs by more than an order of magnitude. This suggests that the density of the shocked gas is not the primary factor causing the division between the two groups of ejecta knots.

Figure 9. Refer to the following caption and surrounding text.

Figure 9.  F([Fe ii] 1.295)/F([Fe ii] 1.257) vs. R([S ii]/[Fe ii]), R([P ii]/[Fe ii]), and R(He i/[Fe ii]). The knot symbols and colors are the same as in Figure 4. The horizontal lines mark the F([Fe ii] 1.295)/F([Fe ii] 1.257) ratios of ionized gas of ne = 103, 104, and 105 cm−3 at 104 K in collisional equilibrium. The parameter ρ in each panel is the correlation coefficient between two ratios.

Standard image High-resolution image

4.3. Optical Counterparts

Numerous dense ejecta knots have been identified around the Cas A SNR in previous optical studies, and it will be interesting to check if our NIR ejecta knots have optical counterparts. We use the optical ejecta knot catalog of Hammell & Fesen (2008), who identified 1825 compact optical knots that lie beyond a radial distance of from the explosion center in Hubble ACS/WFC images taken with three different broadband filters (i.e., F652W, F775W, and F850LP) in 2004. They were able to classify the optical knots into three groups based on their flux ratios in the three filters: (1) [N ii] knots dominated by [N ii] λ λ6548, 6583 emissions, (2) [O ii] knots dominated by [O ii] λ λ7319, 7330 emissions, and (3) fast-moving knot (FMK)–like knots displaying filter flux ratios suggestive of [S ii], [O ii], and [Ar iii] λ7135 emission line strengths similar to the FMKs found in the remnant's main ejecta shell (see also Fesen et al. 2006). Of the 1825 knots identified by them, 444 were [N ii] knots, 192 were [O ii] knots, and 1189 were FMK-like knots. Their spatial distributions are distinct; [N ii] knots are arranged in a broad shell around the remnant, [O ii] knots are clustered around the base of jets, and FMK-like knots are mainly confined to NE and SW jet areas (see also Fesen & Milisavljevic 2016). We have compared the spatial positions of our 67 NIR ejecta knots with those of the optical outer knots shifted to the 2017 epoch assuming free expansion. Considering the seeing and slit width of our observations, we have searched their optical counterparts within a 0farcs5 radius. Note that 12 out of the 67 ejecta knots are located inside , where Hammell & Fesen (2008) did not search optical knots, so we can check the counterparts for only 55 NIR knots.

We found that 38 out of the 55 NIR knots have optical counterparts. All of them have FMK-like knots as a counterpart, except for one that appears to have a [N ii] knot as a counterpart. No NIR knots have [O ii] knots as a counterpart. It is not unusual that the NIR knots that are identified from the ground-based, deep [Fe ii] image are resolved into multiple knots in high-resolution HST images, implying that the knots can have a complicated structure consisting of subknots. Figure 10 shows the optical counterparts of the sample of the observed NIR knots. Some NIR knots appear as single, very compact optical knots, but they often have a more complex structure that is resolved into smaller subcomponents. In some cases, an NIR knot is found to be a part of a larger, more extended structure in the optical. We further note that optical knots tend to cluster together with other knots of the same type. However, there are rare cases where different types of optical knots are found in close proximity to each other, especially in the main shell region. While these can complicate the identification of optical counterparts, we do not think the overall trend will be affected. Figure 11 shows the pie charts for the optical counterparts of the three groups of NIR knots. About three-fourths of the S- and Fe-rich knots have optical counterparts, while about half of the He i knots have optical counterparts. It is worth noticing that one-fourth of the S-rich knots have no optical counterpart, although, since the F850LP filter is sensitive to the [S ii] 1.03 and [S iii] 0.907, 0.953 μm lines, the S-rich knots emitting strong [S ii] and [S iii] lines should at least be detected in the deep F850LP image. The S-rich knots not detected in the deep optical image may have appeared after 2004, when the HST images were taken. It has been found that the knots in the NE jet and Fe K plume areas reveal substantial emission variations ("flickering") in an interval as short as 9 months (Fesen et al. 2011; see also Figure 10). Indeed, we were able to confirm that several knots (e.g., knots 2, 7, and 15) were not present in the 2004 HST image but appeared in a subsequent 2010 HST image.

Figure 10. Refer to the following caption and surrounding text.

Figure 10. Optical counterparts of the sample of NIR knots detected in slits 1–9, 1–16, 2–12, and 3–23. For each slit area, four images are shown: the deep [Fe ii] image and three HST images obtained with the F098M, F625W, and F850LP filters. North is up, and east is to the left. In each panel, the gray rectangle represents the location of the 10''-long MOS slit, and the circles mark the knot locations, the colors of which represent the knot groups: S-rich knots in blue, Fe-rich knots in red, and He i knots in green. Note that the positions of the slits and circles are shifted to the epochs of the background images, which are given at the top of the individual columns.

Standard image High-resolution image
Figure 11. Refer to the following caption and surrounding text.

Figure 11. Pie charts showing the optical counterparts of 55 NIR ejecta knots. The red and blue slices represent the percentages of the NIR knots with optical [N ii] and FMK-like knots as counterparts, respectively, while the green slice represents the percentage of the NIR knots without optical counterparts. Note that the statistics of S-rich knots includes only 32 knots out of 44. The remaining 12 S-rich knots are located in the inner region of the Cas A SNR, where the characteristics of the optical knots were not investigated by Hammell & Fesen (2008).

Standard image High-resolution image

We also found that 70% of the Fe-rich knots have FMK-like knots as their optical counterparts. This is inconsistent with our results that most Fe-rich knots do not have [S ii] 1.03 and [S iii] 0.953 μm lines. One possibility is that the flux detected in the F850LP image is not due to [S ii] or [S iii] lines but rather [Fe ii] lines. There are several bright [Fe ii] lines, e.g., 8619, 8894, and 9229 Å (Dennefeld 1986; Hurford & Fesen 1996; Koo et al. 2016), in the wavelength band (8000–10000 Å) of the F850LP filter, while there are few or no bright [Fe ii] lines in the wave bands of the F652W and F775W filters. For example, the expected [Fe ii] 8619/1.257 μm ratio ranges from 0.08 to 0.49 for electron densities of 103–105 cm−3 (Koo et al. 2016). Hence, the [Fe ii] lines can make a substantial contribution to the F850LP filter flux, and, together with the high interstellar extinction toward Cas A, the Fe-rich knots in our sample could have been classified as FMK-like knots in the optical band. Note that the majority of the optical knots in the catalog of Hammell & Fesen (2008) are very faint, with an F850LP flux of less than 5 × 10−17 erg cm−2 s−1, which is more than 2 orders of magnitude smaller than the [Fe ii] 1.257 μm fluxes of Fe-rich knots (Figure 5). Another possibility is that the knots have [S ii] and/or [S iii] lines but that the lines are "weak," i.e., R([S ii]/[Fe ii]) < 3.7, so they have been classified as Fe-rich knots in this work, but they were detected in the deep F850LP image and classified as FMK-like knots. Indeed, at least three of the Fe-rich knots with FMK-like knot counterparts do have [S ii] lines (see Table 6 in Section 5.2). Yet another possibility is a chance coincidence. The area where the Fe-rich knots have been detected is crowded with FMK-like knots (e.g., see Figure 1 of Hammell & Fesen 2008), so we cannot rule out the possibility that an FMK-like knot accidentally falls within the circle of 0farcs5 radius. But, since the distances between the matched knots are usually an order of magnitude smaller than the median distance between the optical knots around Fe-rich knots (i.e., a few arcseconds), most of the matches are probably genuine.

Table 6. Properties of Fe-rich Knots

KnotK2019 No. α (J2000) δ (J2000)Size μ F[Fe II]1.644 μm vr Optical
No.   (arcsec)(km s−1)(erg cm−2 s−1)(km s−1)   Counterpart
(1)(2)(3)(4)(5)(6)(7)(8)(9)(10)(11)(12)
2319723:23:57.39958:50:30.362.26 × 1.445.75E-16−2341.78≤0.13≤0.26FMK
3923523:23:49.20858:49:04.651.85 × 1.484.66E-16−1286≤1.72≤0.100.71FMK
4023423:23:49.04458:49:03.511.69 × 1.445.17E-16−1282≤0.40≤0.040.59FMK
4223223:23:48.30558:48:59.551.85 × 1.3172914.81E-16−1616≤0.48≤0.060.24FMK
4323:23:48.33658:48:57.600.77 × 0.779.30E-17−1987≤1.13≤0.20≤0.60
4523:23:51.23558:49:51.881.47 × 1.393.82E-16−1165≤0.62≤0.090.45
4624123:23:48.61758:48:54.883.29 × 2.3272212.98E-15−1189≤0.26≤0.020.34FMK
4723923:23:48.10658:48:53.213.41 × 2.6169892.89E-15−1159≤0.24≤0.020.38FMK
4824023:23:48.25358:48:43.331.26 × 0.841.90E-16−903≤1.21≤0.20≤0.31FMK
4925123:23:49.10958:48:32.441.21 × 0.821.63E-16−15543.130.350.86FMK
5125023:23:48.90958:48:28.592.26 × 1.769.13E-16−207≤0.39≤0.138.02
5325923:23:47.85058:48:17.755.01 × 2.8976403.36E-15−6441.070.110.30FMK
5627223:23:46.13658:48:04.651.46 × 1.223.66E-16−1413≤0.85≤0.09≤0.16
6329023:23:44.34258:47:32.271.22 × 0.891.87E-16−2023≤3.43≤0.37≤0.74
6428823:23:43.48358:47:33.935.86 × 1.4665801.10E-15−1725≤1.82≤0.24≤1.00FMK
6729923:23:42.55458:46:21.522.54 × 2.091.81E-151957≤0.18≤0.032.25FMK

Note. (1) Knot serial number in this work, (2) knot serial number in Koo et al. (2018), (3) and (4) central coordinates, (5) major and minor axes, (6) velocity of proper motion, (7) observed [Fe ii] 1.644 μm flux, (8) radial velocity, (9)–(11) line ratios, (12) optical counterpart (see Section 4.3). Columns (3)–(7) are from Koo et al. (2018).

Download table as:  ASCIITypeset image

Lastly, it is interesting that about half of the He i knots have FMK-like knots as their optical counterparts. This may indicate that He i knots are just S-rich knots with [S ii], [S iii], and [Fe ii] lines below our detection limit. But, since most of the He i knots are also located in the area crowded with FMK-like knots, a chance coincidence is possible. We will explore the nature of He i knots in Section 5.3.

5. Discussion

5.1. S-rich Ejecta Knots in the NE Jet

A most striking morphological feature in Cas A known from early optical studies is the NE jet structure (van den Bergh & Dodd 1970; Fesen & Gunderson 1996). It is composed of the streams of bright knots confined into a narrow cone well outside the northern SN shell that appears to have punched through the main ejecta shell. The jet structure extends to beyond the main shell, implying ejection velocities of up to 16,000 km s−1, which is three times faster than the bulk of the O- and S-rich ejecta in the main ejecta shell (Kamper & van den Bergh 1976; Fesen & Gunderson 1996; Fesen 2001; Fesen & Milisavljevic 2016). Later, the SW "counterjet" composed of numerous optical knots with comparable maximum expansion velocities was detected on the opposite side (Fesen 2001; Milisavljevic & Fesen 2013). The origin of the bipolar NE jet–SW counterjet structure is not fully understood, although it must have been generated very near the explosion center during the first seconds of the explosion (see below). Critical information about the origin of the jet comes from its chemical composition. Previous optical studies showed that the NE jet knots exhibit strong S, Ca, and Ar emission lines with comparatively weak oxygen lines (van den Bergh 1971; Fesen & Gunderson 1996; Fesen 2001; Hammell & Fesen 2008; Milisavljevic & Fesen 2013; Fesen & Milisavljevic 2016). The knots lying farther out and possessing the highest expansion velocities show no detectable emission lines other than those of [S ii] (Fesen & Milisavljevic 2016). The knots near the jet base, on the other hand, show O and N emission too (Fesen & Gunderson 1996). The jet structure has been observed in the X-ray and IR too. In the X-ray, the jet material shows strong Si, S, Ar, and Ca lines but not a particularly strong Fe K line (Hughes et al. 2000; Hwang et al. 2004; Vink 2004). The Fe K emission has been found to be strong in other areas of Cas A, particularly the eastern area (Hwang & Laming 2003, 2012; Picquenot et al. 2021; see Section 5.2). These X-ray studies showed that the elemental composition of the jet, enriched in the intermediate-mass elements with a limited amount of Fe, is similar to that of oxygen or incomplete Si burning (see also Ikeda et al. 2022). In the IR, the ringlike jet base lifted from the main ejecta shell is pronounced in Ar emission (DeLaney et al. 2010). These findings suggest that the jet originated in the deep interior of the progenitor star where Si, S, Ar, and Ca are nucleosynthesized and penetrated through the outer stellar envelope.

Figure 12 gives a detailed view of the NE jet structure in [Fe ii] 1.644 μm emission. In the image, the NE jet appears as streams of bright knots confined to a fanlike structure with an opening angle of 10° viewed at the explosion center. The central stream at P.A. ≃ 64° is most prominent, extending to (=4.6 pc) from the explosion center. In the northern area above the central stream, there are many bright knots scattered over the area that appear to be aligned along several stream lines, while, in the southern area below the central stream, there is a prominent but fuzzy stream of ejecta material possessing a finite width at P.A. ≃ 70°, and not many bright knots are in between. This jet morphology is very similar to what we have already seen in [S ii] λ λ6716, 6731 emission (Fesen & Gunderson 1996).

Figure 12. Refer to the following caption and surrounding text.

Figure 12. Detail of the NE jet area. The background is the deep [Fe ii] image. North is up, and east is to the left. The knots where NIR spectra have been obtained in this work are marked by squares, the colors of which represent the knot groups: S-rich knots in blue, Fe-rich knots in red, and He i knots in green.

Standard image High-resolution image

In Figure 12, we mark the ejecta knots observed in this work using different colors for S-rich, Fe-rich, and He i knots. We see that the NE jet knots are mostly S-rich (see also Figure 6). We have not detected an Fe-rich knot exhibiting only [Fe ii] lines. There is one knot classified as an Fe-rich knot located between the central and southern streams, but it exhibits [S ii] lines of moderate strength (R([S ii]/[Fe ii]) = 1.8). These results are consistent with the previous results from optical and X-ray studies that the jet material shows strong S lines but not particularly strong Fe lines. In the optical bands, the NE jet knots show strong O, S, and Ar lines like the other S-rich knots and are not clearly distinguished from the S-rich knots (FMKs) in the main ejecta shell, although, for example, a few outermost knots show stronger than usual [Ca ii] λ λ7291, 7324 lines that could be due to abundance differences (Fesen & Gunderson 1996). In the NIR, however, the NE jet knots are distinguished from the S-rich knots in the main ejecta shell in [P ii] emission. The majority of the dense knots in the NE jet show weak or no [P ii] lines, whereas the S-rich knots in the main ejecta shell show strong [P ii] lines. This is shown in Figure 13, where the line ratios of the NE jet knots are plotted as a function of distance from the explosion center (see also Figure 3). The figure shows that R([S ii]/[Fe ii]) remains roughly constant along the jet, whereas R([P ii]/[Fe ii]) and R([P ii]/[S ii]) appear to decrease systematically along the jet. Note that the [P ii] line has not been detected in the majority of the knots near the tip of the jet. For comparison, the three ratios for the S-rich knots in the main ejecta shell are (R([S ii]/[Fe ii]), R([P ii]/[Fe ii]), R([P ii]/[S ii])) = (14.2 ± 13.5, 1.46 ± 1.31, 0.117 ± 0.083) (see Figure 2 of Koo et al. 2013).

Figure 13. Refer to the following caption and surrounding text.

Figure 13. Line ratios of the NE jet knots as a function of distance from the explosion center. The S- and Fe-rich knots are denoted by blue and red squares, respectively, while open symbols indicate upper limits.

Standard image High-resolution image

The simplest interpretation of the low [P ii]/[Fe ii] and [P ii]/[S ii] ratios of the NE jet knots is that the amount of phosphorus present there is lower compared to the knots in the main shell. Although the knots in the NE jet and the main shell have different shock environments, the decrease in these ratios from the main shell to the tip of the NE jet while the [S ii]/[Fe ii] ratio remains constant suggests that the phosphorus abundance is probably the primary factor affecting the line ratios. Therefore, we consider the weakness of the [P ii] line in the NE jet knots as direct evidence that the NE jet originated from a layer with a high sulfur and low phosphorus content. And, since 15P is mainly produced in carbon- and neon-burning layers (Arnett 1996; Woosley et al. 2002), it implies that the NE jet was ejected from a region below the explosive Ne-burning layer, i.e., the O-burning or incomplete Si-burning layers. This is consistent with the result of previous X-ray studies, which showed that the elemental composition of diffuse ejecta material in the NE jet is similar to that of oxygen and incomplete Si burning (Vink 2004; Ikeda et al. 2022).

We now discuss the implication of our results on the physical origin of the NE jet. The origin of the NE jet and its role in the SN explosion have been subject to study since its discovery in the 1960s. Minkowski (1968) suggested that the jet structure (or the "flare" in his paper) could be a surviving part of a spherically symmetric shell due to nonspherical distribution of the ambient ISM (see also Fesen & Gunderson 1996). But it is now rather well established that the NE jet originated from a Si–S–Ar-rich layer deep inside the progenitor and penetrated through the outer N- and He-rich envelope (Fesen 2001; Fesen et al. 2006; Laming et al. 2006; Milisavljevic & Fesen 2013). The chemical composition of the jet material and the spatial correlation between the jet and the disrupted main ejecta shell are strong evidence for that. On its launching mechanism, there have been several theoretical propositions in relation to the core-collapse SN explosion models. In the classical "jet-induced" explosion model, where the magnetic field anchored to the collapsing core is amplified during the collapse and drives the SN explosion, a magnetocentrifugal jet can be ejected along the rotation axis (Khokhlov et al. 1999; Wheeler et al. 2002; Akiyama et al. 2003). Recent observational studies, however, suggest that the properties of the NE jet are not consistent with the jet-induced explosion model; there is no "pure" Fe ejecta in the jet (Ikeda et al. 2022), the opening angle is large, the kinetic energy is insufficient, and the kick direction of the neutron star kick is not along but perpendicular to the jet axis (Fesen 2001; Hwang et al. 2004; Laming et al. 2006; Milisavljevic & Fesen 2013; Fesen & Milisavljevic 2016; see Soker 2022). Our results are also in line with previous results and do not support the jet scenario. We also have not found dense Fe-rich knots around the tip of the NE jet in this study. The opening angle of the jet in Figure 12 is small (≲15°), but the total kinetic energy will be less than the previous estimate of 1 × 1050 erg of the optical knots in the NE jet (Fesen & Milisavljevic 2016). Hence, our results also do not support the jet-induced explosion model.

In the neutrino-driven explosion model, which is the modern paradigm of the core-collapse SN explosion, narrow high-velocity ejecta can be generated either during or after the explosion by the newly formed neutron star. During the explosion, high-speed jetlike fingers can be produced from the core by neutrino convection bubbles (e.g., Burrows et al. 1995; Janka & Mueller 1996). In particular, the interface between Si and O composition is unstable, seeded by the flow structures resulting from neutrino-driven convection, leading to a fragmentation of this shell into metal-rich clumps (Kifonidis et al. 2003). Orlando et al. (2016) showed that Si-rich, jetlike features such as the NE jet in Cas A could be reproduced by introducing overdense knots moving faster than the surrounding ejecta just outside the Fe core. But a recent 3D simulation for modeling the Cas A SNR as a remnant of a neutrino-driven SN, going from the core collapse to the fully developed remnant, could not reproduce a structure similar to the NE jet (Orlando et al. 2021). Another possibility is a hydrodynamic jet from a newly formed neutron star. In a standard neutrino-driven explosion, a fraction of ejecta falls back to the newly formed neutron star to form a disk where an MHD jet can be launched (Burrows et al. 2005, 2007; Burrows & Vartanyan 2021; Janka et al. 2022). But as far as we are aware, there are few theoretical predictions that can be compared with observations for such an MHD jet. So, for the neutrino-driven explosion model, it is still an open question whether the NE jet can be explained by this model.

5.2. Fe Ejecta Knots in the Fe K plume area

We discovered Fe-rich knots in the Fe K plume area (see Figure 14). This is the area where the Fe-rich plume with bright Fe K emission has been detected in X-rays (Willingale et al. 2002; Hwang & Laming 2003; Hwang et al. 2004; Hwang & Laming 2012; Tsuchioka et al. 2022). The ejecta material in this Fe K plume has no S or Ar, and its Fe/Si abundance ratio is up to an order of magnitude higher than expected in the incomplete Si-burning layer, indicating that it originated from a complete Si-burning process with α-rich freeze-out, where the burning products are almost exclusively 56Fe because free α particles are unable to reassemble to heavier elements on a hydrodynamic timescale (Arnett 1996; Woosley et al. 2002). Recently, Sato et al. (2021) detected Ti and Cr in this area, the ratio of which to Fe is consistent with explosive complete Si burning with α-rich freeze-out. Therefore, the material in the Fe K plume has likely been produced in the innermost, high-entropy region during the SN explosion, where the temperature is high and the density is low (e.g., Arnett 1996; Thielemann et al. 1996). In 2D projected maps, however, the Fe K plume is located much further outside the main shell, which is composed of ejecta synthesized from incomplete Si burning, e.g., Si, S, and Ar, which could be interpreted as an overturning of the ejecta layers (Hughes et al. 2000). But a 3D reconstruction of the ejecta distribution using infrared and Chandra X-ray Doppler velocity measurements showed that, in the main shell, the Fe K plume occupies a "hole" surrounded by a ring structure composed of Si and Ar (DeLaney et al. 2010; see also Milisavljevic & Fesen 2013). As an explanation of the observed morphology, DeLaney et al. (2010) proposed a model where an ejecta "piston" faster than the average ejecta has pushed through the Si layers, rather than interpreting it as an overturning of the layers. In this model, the Si–Ar ring structure in the main shell represents the circumference of the piston intersecting with the reverse shock. Recently, Tsuchioka et al. (2022) showed that the velocity of the ejecta in the outermost Fe K plume is >4500 km s−1, much higher than that of Si-/O-rich ejecta, suggesting that the ejecta piston was formed at the early stages of the remnant evolution or during the SN explosion, which is consistent with the results of numerical simulations. On the other hand, recent 3D numerical simulations demonstrated that Ni clumps, created during the initial stages of the SN explosion in the innermost region by the convective overturning in the neutrino-heating layer and the hydrodynamic instabilities, later decay to 56Fe and can push out the less dense ejecta to produce structures such as the Fe K plume (Orlando et al. 2016, 2021). But Hwang & Laming (2003) pointed out that the morphology and the "ionization age" (ne t, where ne is the electron density, and t is the time since the plasma was shock-heated) of the Fe ejecta in the Fe K plume are inconsistent with a Ni bubble origin.

Figure 14. Refer to the following caption and surrounding text.

Figure 14. Detail of the Fe K plume area in the deep [Fe ii] and Chandra X-ray Fe K emission images. North is up, and east is to the left. The knots where NIR spectra have been obtained in this work are marked by circles (P.A. ≥ 78°) or squares (P.A. < 78°), the colors of which represent the knot groups: S-rich knots in blue, Fe-rich knots in red, and He i knots in green. Note that the positions of the knot markers are shifted to the epochs of the background images, i.e., 2013 September for the deep [Fe ii] image and 2004 April for the Chandra image. The yellow contour marks the approximate boundary of the SNR in the radio in 2003 (DeLaney 2004), and it is left as a reference point.

Standard image High-resolution image

Our detection of Fe-rich knots in the Fe K plume area has quite significant implications. Figure 14 shows the distribution of IR ejecta knots in the Fe K plume area. We see a group of Fe-rich knots aligned in the north–south direction along the forward shock front. They are spatially confined in the area immediately outside the boundary of the diffuse X-ray Fe ejecta. Their radial velocities are in a narrow range, from −2000 to −200 km s−1 (see Figure 8 and Table 6), so, in a position–velocity diagram, they appear to be clustered (see Figure 6(b)). The ejection angle of this Fe-rich ejecta knot cluster is 14° ± 5° from the plane of the sky toward us. (In Figure 6(b), we see two Fe-rich knots located far off the cluster. One of them is knot 23 in the NE jet area, and the other is knot 67 near the southern edge of the field. The latter has a radial velocity of +1960 km s−1, which is very different from that of the Fe-rich ejecta knot cluster. They are not considered to belong to the cluster.) The radial velocities of the X-ray-emitting diffuse ejecta in this area measured with the Chandra High Energy Transmission Grating are from −2600 to −500 km s−1, comparable with those of the Fe-rich knots (Lazendic et al. 2006; Rutherford et al. 2013). On the other hand, the proper motion of the Fe-rich knots has a velocity of 6600–7600 km s−1 (Table 6), which is considerably higher than that of the diffuse ejecta at the tip of the Fe K plume (i.e., 4500–6700 km s−1; Tsuchioka et al. 2022). Hence, the Fe-rich knots move a little faster than the diffuse X-ray-emitting Fe-rich ejecta. The spatial and kinematic relations strongly support the physical association between the dense Fe-rich knots and the X-ray-diffuse Fe ejecta. And the spectral properties of the Fe-rich knots suggest that they are also produced by explosive Si complete burning with α freeze-out as the diffuse Fe ejecta. The He i 1.083 μm lines detected in the majority of Fe-rich knots might be from the helium remains of the α-rich freeze-out process (Table 6). In α-rich freeze-out, the burning products are almost exclusively 56Fe, but a considerable amount of He could remain. The local He mass fraction can be as high as 0.2, depending on the model (Rauscher et al. 2002; Woosley et al. 2002; Thielemann et al. 2018). The detection of dense Fe ejecta knots implies that the Fe ejecta produced in the innermost region was very inhomogeneous.

Theoretically, dense metallic ejecta knots can be produced during the explosion. As mentioned above, numerical simulations have shown that when the innermost ejecta expanding aspherically encounters a composition shell interface, it is compressed to a dense shell, which fragments into small, dense metallic clumps by hydrodynamic instability (Kifonidis et al. 2000, 2003; Hammer et al. 2010; Wongwathanarat et al. 2017). It is uncertain whether the dense knots had initially higher velocities. However, it is possible that these dense knots were initially comoving with the surrounding diffuse ejecta, but once they encounter the reverse shock, the dense and diffuse ejecta are decoupled, i.e., the diffuse Fe ejecta are decelerated more, so the dense Fe clumps are located ahead of the diffuse Fe ejecta. This explanation has been proposed for dense Fe knots beyond the SNR boundaries detected in X-rays in young Type Ia SNRs (Wang & Chevalier 2001; Tsebrenko & Soker 2015). A caveat in this scenario is that, since the radioactive decay from 56Ni to 56Fe produces postexplosion energy, the Fe clumps are expected to expand unless they are optically thin to γ-ray emission and should be characterized by the diffuse morphology in the remnant phase (Blondin et al. 2001; Hwang & Laming 2003; Gabler et al. 2021). Indeed, in the interior of Cas A, the bubble-like structures in S-rich ejecta filled with 44Ti are detected, supporting this scenario (Milisavljevic & Fesen 2013; Grefenstette et al. 2017). Therefore, the dense compact Fe-rich knots discovered in Cas A need an explanation. A possible explanation is that the dense Fe-rich knots are not 56Fe decayed from 56Ni but 54Fe (Wang & Chevalier 2001). But 54Fe is not a dominant element (see Hwang & Laming 2003, and below). In core-collapse SNe particularly, 54Fe turns into 58Ni (stable nuclei of nickel) in the α-rich freeze-out region of the complete Si-burning zone, so it is only detected in the incomplete Si-burning layer, where the S and Si abundances are comparable with or much larger than the 54Fe and 56Ni/56Fe abundances (e.g., Thielemann et al. 1996, 2018). Therefore, the Fe in our Fe-rich knots is most likely 56Fe, not 54Fe. Another possibility, which was proposed by Hwang & Laming (2003) to explain the compact X-ray knots embedded in the Fe K plume, is that the optical depth of the knots is not large enough for the γ-rays from the radioactive decay from 56Ni to 56Fe to escape. Using a hydrodynamics model for ejecta evolution, they showed that, for a typical knot of diameter ∼3'' and electron density ∼10 cm−3, the optical depth of the γ-rays (∼1 MeV) becomes less than 1 about 5 days after explosion, which is much shorter than the 56Ni and 56Co lifetimes, i.e., 8.8 and 111.3 days, respectively. We consider that the Fe-rich knots discovered in our NIR study could be a denser and faster version of the X-ray knots. The NIR Fe-rich knots are dense (e.g., ≳102 cm−3; Section 4.2.2), and their apparent sizes measured with an ∼1'' resolution are comparable to the X-ray knots (Table 6). So, if we naively accept these numbers, the NIR Fe-rich knots are expected to be optically thick for γ-rays generated from the radioactive decay from 56Ni to 56Fe. However, the structure of the NIR Fe-rich knots is not resolved. It is possible that they consist of multiple compact subknots with relatively low column densities, enabling the γ-rays to escape primarily from the knot. The HST images in Figure 10 show that indeed, the intrinsic sizes of the NIR knots could be much smaller than their apparent sizes in the deep [Fe ii] image. Furthermore, the ejecta clumps are destroyed by the passage of the reverse shock (e.g., see Figure 4 of Slavin et al. 2020), so the initial sizes of the NIR and X-ray knots could have been considerably smaller. Hence, the Ni bubble effect may not be significant for either type of Fe-rich knot. We will conduct subsequent studies in a forthcoming paper.

Before leaving this section, it is worth noticing that our observation has been done toward bright knots in the deep [Fe ii] image obtained using a narrowband filter with a bandwidth of ±2800 km s−1 (see Section 2), so there is a possibility that there are additional Fe-rich knots outside this radial velocity range. However, previous studies have shown that the diffuse Fe ejecta typically have radial velocities in the range of vr = −3000 to −1000 km s−1 (DeLaney et al. 2010; Tsuchioka et al. 2022). Therefore, it is likely that the majority of the bright Fe-rich ejecta knots are included in the [Fe ii] image of Figure 14, and our result provides insight into the distribution of the Fe-rich ejecta knots.

5.3. Knots with Only He i 1.083 μm lines

We detected seven knots that show only He i 1.083 μm lines. The He i 1.083 μm line intensities of these He i knots are not particularly strong in comparison with other metal-rich ejecta knots (Figure 5). But, since the [Fe ii] 1.257 μm line is not detected, they are noticeable by their very high R(He i/[Fe ii]), which is 1–2 orders of magnitude higher than the typical S- and Fe-rich ejecta knots (see the bottom panel of Figure 4). Table 7 lists the knots with R(He i/[Fe ii]) greater than 8, where we can see that the majority are He i knots.

Table 7. Ejecta Knots with High R(He i/[Fe ii])

 Knot vr FHe I 1.083 μm KnotOptical
 No.(km s−1)(erg cm−2 s−1)    TypeCounterpart
 (1)(2)(3)(4)(5)(6)(7)(8)(9)
He I knots2632413.89E-15≥8.21≤0.55≤0.25HeFMK
 33−31526.43E-15≥27.40≤0.550.31HeFMK
 3523142.22E-15≥14.50≤0.39≤0.57He
 38−33121.34E-15≥14.70≤0.47≤0.23HeFMK
 41−2681.28E-15≥22.09≤0.35≤0.24He
 527658.87E-15≥83.32≤0.080.11He
 5411861.05E-14≥77.55≤0.070.07He
Other types of knots123363.09E-1417.2173.230.240.07SFMK
 51231.80E-1371.8233.102.170.08SFMK
 32−31697.67E-15≥46.20≥54.430.85≤0.09S
 51−2061.49E-14≤0.398.02≤0.050.13Fe

Note. (1) Knot serial number in this work, (2) radial velocity, (3) He i 1.083 μm flux, (4)–(7) line ratios where [C I] = [C I] 0.985 μm, (8) knot type, (9) optical counterpart (see Section 4.3).

Download table as:  ASCIITypeset image

He i knots are detected along the forward shock front, from the southern base of the NE jet to the Fe K plume area (see Figure 6). The radial velocities of the He i knots are spread over a very wide range, from −3300 to +3200 km s−1 (Figure 8 and Table 7). Since most of them have very high radial velocities, they are probably not circumstellar material. Their spectral properties are also different from the low-velocity circumstellar clumps in Cas A, which are called quasi-stationary flocculi (QSFs). The QSFs are He- and N-enriched circumstellar material ejected by the progenitor before the SN explosion and have recently been shocked by the SN blast wave (Peimbert & van den Bergh 1971; McKee & Cowie 1975; see also Koo et al.2023 and references therein). In the NIR, QSFs show strong He i lines, together with moderately strong [Fe ii] and [S ii] lines and weak Paγ and [C i] 0.985 μm lines (Chevalier & Kirshner 1978; Gerardy & Fesen 2001; Lee et al. 2017). The ratio of Paγ to He i 1.083 μm is typically ∼0.03 (see Koo et al. 2023). For comparison, all He i knots show strong (>10−15 erg cm−2 s−1) He i lines without [Fe ii] or [S ii] lines, and about half of them also show weak [C i] 0.985 μm lines.

What is the origin of the He i knots? The first possibility that we can consider is that they are the debris of the He-rich envelope material of the progenitor star expelled by the SN explosion, like the high-velocity N-rich knots (NKs) detected in optical studies (Fesen et al. 1987, 1988; Fesen & Becker 1991; Fesen 2001; Hammell & Fesen 2008; Fesen & Milisavljevic 2016). The spectra of NKs are dominated by strong [N ii] λ λ6548, 6583 lines in the wavelength region 4500–7200 Å. For most NKs, the Hα line is either much weaker than the [N ii] lines or undetected, so NKs are considered important evidence that the Cas A progenitor had a thin N-rich and H-poor envelope at the time of the SN explosion. A deep, ground-based [N ii] imaging observation detected 50 NKs, whereas HST observations revealed numerous NKs distributed around the remnant. In particular, in the eastern area where He i knots have been detected, many NKs had been detected (Figure 15; see also Figure 5 of Fesen 2001 or Figure 9(E) of Fesen & Milisavljevic 2016). The radial velocities of the bright NK knots in the eastern area are in the range of −2000 to +2500 km s−1 (Fesen 2001). So the association of He i knots with NKs seems plausible considering that He abundance might be enhanced in the CNO-processed, N-rich envelope of the progenitor star. Apparently, however, none of the He i knots has an optical NK counterpart (Table 7; see Section 4.3). Another difficulty in connecting He i knots to NKs is the nondetection of the He i λ5876 line in NKs (Fesen & Becker 1991; Fesen 2001). The recombination emissivity of the He i λ5876 line is a factor of 2–10 smaller than that of the He i 1.083 μm line, depending on the electron density (102–104 cm−3) in the Case B approximation (e.g., Draine 2003). The flux will be further reduced by interstellar extinction. Hence, the He i λ5876 line is expected to be much weaker than the He i 1.083 μm line, but it is not obvious if that can explain the nondetection of He i λ5876 lines. An interesting possibility is that the He i knots originated from the lower portion of the He-rich envelope of the progenitor star where N is depleted. If a large amount of carbon is dredged up by convection, nitrogen could be depleted because 14N mixed into the lower layer by convection participates the reaction. Such a 14N-depleted region formed by convective dredge-up is often found at the pre-SN stage in massive star models (Y. Jung et al 2023, in preparation). Observationally, our slits are positioned toward [Fe ii] line-emitting knots by design. So it is possible that we are picking up He-rich, N-depleted knots that are located closer to the explosion center than the optical NKs but coincidentally overlap with the same region of sky as the Fe-emitting knots (Figure 15). Future studies may further explore the possibility.

Figure 15. Refer to the following caption and surrounding text.

Figure 15. Distribution of optical NKs (yellow circles) in the catalog of Hammell & Fesen (2008) shown on the deep [Fe ii] and HST F625W images. North is up, and east is to the left. The circles/squares and their colors are the same as in Figure 14. The black bars in the deep [Fe ii] image represent the locations of the 10''-long MOS slits. In the F6250W image, only the bright NK knots with an F6250W flux larger than 6 × 10−17 erg cm2 s−1 are shown. Note that the positions of the slit and knot symbols are shifted to the epochs of the background images, i.e., 2013 September for the deep [Fe ii] image and 2004 March for the HST image.

Standard image High-resolution image

Another possibility is that they are metal-enriched SN ejecta with a relatively high He abundance. We note that some S- and Fe-rich knots have very high R(He i/[Fe ii]), as high as the He i knots, and that the majority of He i knots have [S ii] 1.03 and [C i] 0.985 μm flux upper limits that are not particularly strong (Table 7). This explanation also appears to be consistent with the idea that half of the He knots seem to be FMK-like knots in optical classifications (Figure 11). The overabundance of low atomic elements could be attributed to significant mixing among the nucleosynthetic layers, as pointed out in previous X-ray and optical studies (Fesen & Becker 1991; Fesen & Gunderson 1996). From optical spectroscopy, Fesen & Becker (1991) reported the detection of a hybrid ejecta knot showing O and S emission lines along with [N ii] and Hα emission lines in the NE jet region, and they called it a "mixed emission knot" (MEK). Additional optical studies toward the outer region of the remnant have further detected several MEKs, which suggest chemical mixing among the nucleosynthetic layers between the explosive O-burning layer and the N-rich photospheric envelope of the progenitor (Fesen & Gunderson 1996; Fesen 2001). Yet another possibility is that the He i knots are the remains of the complete Si burning with α-rich freeze-out in the innermost region of the progenitor as helium in Fe-rich knots (see also Lee et al. 2017). But it seems to be difficult for this scenario to explain the difference in radial velocities between the He i and Fe-rich knots.

It is also worth noticing that Figure 6 does not represent the overall distribution of He i knots. Since our observation was done toward ejecta knots bright in [Fe ii] emission, the detection of He i knots was unexpected. The He i knots were detected because they fell on the slit by chance. In order to see the overall distribution of the He i knots, which is essential for understanding their origin, we need He i 1.083 μm narrowband observations of the Cas A remnant.

5.4. Extended H Emission Features Associated with East Cloud

In addition to the 67 SN ejecta knots, we also detected three extended H emission features in the outer eastern region (Section 3), and their physical properties are listed in Table 8. They show only unresolved (FWHM ≲9 Å) H i lines, e.g., Paγ 1.094 and Paβ 1.282 μm lines, without any emission lines from metallic elements, and have radial velocities between −20 and −30 km s−1, indicating that they are of CSM/ISM origin. We can estimate the line-of-sight extinction to the H-emitting gas using the flux ratio of the H i emission lines. All three features show the Paβ line, and the brightest one (H2) also shows the Paγ line. The observed H i line ratio, F(Paγ)/F(Paβ), is 0.41 ± 0.09. By comparing with the unreddened line ratio of 0.56 at Te = 104 K and ne = 104 cm−3 (Hummer & Storey 1987), we obtained an NH of cm−2, which corresponds to an AV of mag. The extinction-corrected F([Fe ii] 1.257)/F(Paβ) is <0.2 (2σ upper limit), which is at least an order of magnitude smaller than that in the typical shock-excited CSM/ISM detected in the several Galactic SNRs (2–20; Koo & Lee 2015, and references therein). This low F([Fe ii] 1.257)/F(Paβ) indicates that the H i lines are associated with a photoionized gas.

Table 8. Physical Parameters of Three H Emission Features

Feature α (J2000) δ (J2000)P.A. drad FWHM vLSR F(Paγ)F(Paβ)
Name  (deg)(arcmin)(Å)(km s−1)(10−17 erg s−1 cm−2)
(1)(2)(3)(4)(5)(6)(7)(8)(9)
H123:23:45.81+58:48:52.588.72.348.8 (1.8)−29 (17)⋯(1.1)5.1 (1.3)
H223:23:49.14+58:48:28.097.32.798.7 (0.6)−19 (6)6.1 (1.2)14.6 (1.4)
H323:23:47.94+58:48:15.7102.12.677.0 (1.0)−28 (10)⋯(0.8)6.6 (1.2)

Note. (1) Feature name. (2) and (3) Central coordinates of the feature. (4) and (5) P.A. and radial distance of the feature from the center of the SN explosion. (6) and (7) FWHM and LSR velocity of the Paβ line. The numbers in parentheses are their 1σ statistical uncertainties from single Gaussian fittings. They do not include the systematic uncertainties from the absolute wavelength calibration (∼12 km s−1; see Section 2). (8) and (9) Observed fluxes of the Paγ and Paβ lines. The numbers in parentheses are their 1σ statistical uncertainties from single Gaussian fittings. They do not include the systematic uncertainties from the absolute photometric calibration (∼20%; see Section 2).

Download table as:  ASCIITypeset image

We noticed that the positions of the extended H emission features are coincident with the region where a bright, diffuse optical cloud is detected (Figure 16). The diffuse cloud in the eastern region (hereafter "East Cloud"), located 2farcm7 east of the center of the remnant, was first identified in an early optical photograph taken by R. Minkowski (Minkowski 1968; van den Bergh 1971), and it has a triangular shape with a size of (see also the right panel of Figure 16). In a deep Hα narrowband image, however, it turned out to be more extended and diffuse (see Figure 10 in Fesen 2001 or Figure 2 in Weil et al. 2020). The follow-up optical spectroscopy showed that the East Cloud emits a strong Hα emission line accompanied by relatively weak [O iii] λ λ4959, 5007, [N i] λ λ6548, 6583, and [S ii]λ λ 6716, 6725 lines (Fesen et al. 1987; Weil et al. 2020).

Figure 16. Refer to the following caption and surrounding text.

Figure 16. (a) Continuum-subtracted Hα image of the eastern area of Cas A. The symbols and lines are the same as in Figure 1. (b) Enlarged view of the boxed area in panel (a). The blue circles represent the central positions of the H emission features: H1, H2, and H3 in Table 8. The image was produced from the Isaac Newton Telescope Wide Field Survey data obtained in 2005 October.

Standard image High-resolution image

The origin of the East Cloud, as well as its physical association with Cas A, has been controversial. The early optical studies suggested that it is a small, diffuse H ii region adjacent to the remnant (Minkowski 1968; van den Bergh 1971). However, the lack of OB stars around the eastern cloud suggested that it is a diffuse CSM/ISM near Cas A excited by the UV/X-ray emission from the SN outburst (Peimbert 1971; van den Bergh 1971; Weil et al. 2020). In this scenario, the East Cloud could be either the diffuse ISM surrounding the progenitor star (Minkowski 1968; Peimbert & van den Bergh 1971; Peimbert 1971; van den Bergh 1971) or the diffuse CSM blown out from the progenitor star during its red supergiant phase (Fesen et al. 1987; Chevalier & Oishi 2003; Weil et al. 2020). The peculiar structure of the East Cloud, its projected proximity to the NE jet, and its distinct spectroscopic properties relative to the surrounding diffuse emissions are indirect evidence suggesting the physical association of the East Cloud with the remnant (Fesen et al. 1987; Fesen & Gunderson 1996; Fesen et al. 2006; Weil et al. 2020), while the presence of brightened ejecta knots matching the East Cloud's emission structure provides strong evidence that the East Cloud is physically adjacent to Cas A (Weil et al. 2020).

According to our spectroscopy, the mean vLSR of the Paβ lines weighted by 1/ (where σv is the 1σ uncertainty of the central velocity) is −22 ± 5 km s−1 (see Table 8). Considering the systematic uncertainty from the wavelength calibration (∼12 km s−1; see Section 2), the vLSR of the East Cloud is −22 ± 13 km s−1. This is consistent with the previous result: vLSR = −20 ± 75 km s−1 from a few bright optical lines (Fesen et al. 1987) and vLSR = −13 ± 10 km s−1 from the Hα lines (Alarie et al. 2014). Previous radio observations, however, showed that Cas A is located at the far side of the Perseus spiral arm, implying vLSR ∼ −50 km s−1 for Cas A (Zhou et al. 2018, and references therein). Furthermore, recent NIR observations for the southern optical QSF (knot 24) reported the detection of narrow (FWHM ∼ 8 km s−1) [Fe ii] lines from unshocked pristine CSM with vLSR ≃ −50 km s−1 (Koo et al. 2020), which is consistent with the systematic velocity center vLSR ≃ −50 km s−1 of Cas A. So the velocity of the East Cloud (vLSR ∼ −22 km s−1) is considerably different from the systemic velocity of the Cas A SNR. This might indicate that the East Cloud is not physically associated with the remnant but rather is located ∼2 kpc from us (more than 1 kpc closer than the remnant), assuming the flat Galactic rotation model with the IAU standard (R0 = 8.5 kpc and v0 = 220 km s−1). Relatively low line-of-sight extinction of the East Cloud (A V ∼ 4 mag) derived from the H i lines compared to that toward the remnant (A V = 6–10 mag; see Figure 2) seems to support this possibility. On the other hand, if the East Cloud is the CSM material ejected from the progenitor star, as suggested by Chevalier & Oishi (2003) and Weil et al. (2020), the velocity difference (∼30 km s−1) should represent the line-of-sight component of the ejection velocity of the CSM. And, considering that the East Cloud is located at the eastern outer edge of the remnant, the ejection velocity should have been much higher than 30 km s−1 or the typical wind velocity of red supergiant stars, e.g., 10–20 km s−1 (Smith 2014).

6. Summary

The Cas A SNR has quite a complex structure, manifesting the violent and asymmetric explosion. Two representative features are the NE jet and the Fe K plume in the outer eastern area of the SNR. The NE jet is a stream of ejecta material dominated by intermediate-mass elements beyond the SN blast wave, and the Fe K plume is a plume of X-ray-emitting hot gas dominated by Fe ejecta outside the main ejecta shell. These two features are expanding much faster than the main ejecta shell of the SNR, suggesting turbulent convection and hydrodynamic instabilities in the early stages of the SN explosion. We carried out NIR (0.95–1.75 μm) MOS spectroscopy of the NE jet and Fe K plume regions of the Cas A SNR using MMIRS. In the 2D spectra of 52 MOS slits, which were positioned on the bright knots in the deep [Fe ii] 1.64 μm image of Koo et al. (2018), a total of 67 knots have been identified. All knots show at least one of the following strong lines: [S iii] 0.983, He i 1.083, [S ii] 1.03, and [Fe ii] 1.257 μm. And about one-third of the knots also show a strong [P ii] 1.189 μm line. We find that the knots in different areas show distinctively different ratios of these lines, suggesting their different elemental composition. The NIR lines are emitted from shocked gas, so their intensities depend on shock environments, e.g., the density and velocity of the knots, density contrast between the knot and the ambient gas, ambient pressure, as well as the elemental composition of the knots. In metal-rich ejecta knots, the elemental composition also profoundly affects the cooling function and therefore the physical structure of the shocked gas (e.g., Raymond 2019). Hence, it would be a formidable task to derive the elemental composition from just the NIR spectra. In this work, we simply classify the knots into three groups based on the relative strengths of the [S ii], [Fe ii], and He i lines, i.e., S-rich, Fe-rich, and He i knots, and explore the origins of these knots and their connection to the explosion dynamics of the Cas A SN. We summarize our main results as follows.

  • 1.  
    The NE jet is dominated by S-rich knots. There are no Fe-rich knots without [S ii] lines. The knots show weak or no [P ii] lines, which clearly differentiates them from the S-rich ejecta in the main ejecta shell that have strong [P ii] lines. The low P abundance inferred from their low [P ii]/[Fe ii] line ratios indicates that these S-rich knots were produced below the explosive Ne-burning layer, which is consistent with the results of previous studies. Our results do not support the jet-induced explosion model, but it is also not clear if the NE jet can be explained by the neutrino-driven explosion model.
  • 2.  
    In the Fe K plume area, along the forward shock front just outside the boundary of the diffuse X-ray-emitting Fe ejecta, we detected Fe-rich knots showing only [Fe ii] lines with or without the He i line. These NIR Fe-rich ejecta knots are expanding with velocities considerably higher than the X-ray Fe ejecta. The spatial and kinematic relations support the physical association of these dense NIR Fe-rich knots with the X-ray-diffuse Fe ejecta produced by explosive complete Si burning with α freeze-out. We suggest that the initial density distribution of the Fe ejecta produced in the innermost region was very inhomogeneous and that the dense knots were ejected with the diffuse Fe ejecta but decoupled after crossing the reverse shock.
  • 3.  
    We also detected several He i knots emitting only He i 1.083 μm lines with or without very weak [C i] 0.985 μm lines. They are detected along the forward shock front from the southern base of the NE jet to the Fe K plume area. The origin of these He-rich knots is unclear. They are likely the debris of the He-rich layer above the carbon–oxygen core of the progenitor expelled during the SN explosion, but they could also be metal-enriched SN ejecta with a relatively high He abundance or the remains of the explosive complete Si burning in α-rich freeze-out. The detection of He i knots was unexpected because our MMT observations were directed toward ejecta knots bright in [Fe ii] emission. So an imaging observation is needed to reveal the distribution of the He i knots.
  • 4.  
    In addition to the 67 SN ejecta knots, we also detected three extended H emission features associated with the diffuse cloud, the East Cloud, in the eastern area well beyond the SNR boundary. They show only narrow H recombination lines with low line-of-sight velocities (∼−20 km s−1), indicating that the lines arise from photoionized CSM/ISM. Their velocities are substantially different from the systematic velocity of Cas A (vLSR ∼ −50 km s−1), and this implies that, if the East Cloud is CSM ejected from the progenitor star of the Cas A SN, the ejection velocity should have been much higher than 30 km s−1 considering its location.

Acknowledgments

This work was supported by the K-GMT Science Program (PID: MMT-2017B-4) of the Korea Astronomy and Space Science Institute (KASI). B.C.K. was supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Science, ICT and Future Planning (2020R1A2B5B01001994). The observations reported here were obtained at the MMT Observatory, a joint facility of the Smithsonian Institution and the University of Arizona.

Facilities: MMT - MMT at Fred Lawrence Whipple Observatory, UKIRT - United Kingdom Infrared Telescope.

Software: IDL, Python.

Footnotes

10.3847/1538-4357/acda2d
undefined