ABSTRACT

Recent observations have revealed interstellar features that apparently connect energetic activity in the central region of our Galaxy to its halo. The nature of these features, however, remains largely uncertain. We present a Chandra mapping of the central 2° × 4° field of the Galaxy, revealing a complex of X-ray-emitting threads plus plume-like structures emerging from the Galactic Centre (GC). This mapping shows that the northern plume or fountain is offset from a well-known radio lobe (or the GCL), which however may represent a foreground H ii region, and that the southern plume is well wrapped by a corresponding radio lobe recently discovered by MeerKAT. In particular, we find that a distinct X-ray thread, G0.17−0.41, is embedded well within a non-thermal radio filament, which is locally inflated. This thread with a width of ∼1.6 arcsec (FWHM) is ∼2.6 arcmin or 6 pc long at the distance of the GC and has a spectrum that can be characterized by a power law or an optically-thin thermal plasma with temperature ≳ 3 keV. The X-ray-emitting material is likely confined within a strand of magnetic field with its strength ≳ 1 mG, not unusual in such radio filaments. These morphological and spectral properties of the radio/X-ray association suggest that magnetic field re-connection is the energy source. Such re-connection events are probably common when flux tubes of antiparallel magnetic fields collide and/or become twisted in and around the diffuse X-ray plumes, representing blowout superbubbles driven by young massive stellar clusters in the GC. The understanding of the process, theoretically predicted in analog to solar flares, can have strong implications for the study of interstellar hot plasma heating, cosmic ray acceleration and turbulence.

1 INTRODUCTION

Because of its proximity (D ≈ 8 kpc; 1 arcsec = 0.039 pc), the Galactic Centre (GC) of our Galaxy provides a critical test bed for developing our understanding of astrophysical processes under similar conditions throughout the cosmic time (e.g. Morris & Serabyn 1996; Kruijssen & Longmore 2013; Bland-Hawthorn & Gerhard 2016; Armillotta, Krumholz & Di Teodoro 2020). The environment of the GC is characterized by the high temperature, density, turbulent velocity, and magnetic field of the interstellar medium (ISM), in addition to the strong gravitational tidal force around Sgr A* – the central massive black hole of our Galaxy. Although Sgr A* is quiescent at present (e.g. Wang et al. 2013), it was evidently active in recent past, which has greatly influenced the structure, kinematics, and/or ionization of the ISM of the Galaxy and beyond (e.g. Su, Slatyer & Finkbeiner 2010; Ponti et al. 2015; Bland-Hawthorn & Gerhard 2016; Di Teodoro et al. 2018; Koyama 2018; Ashley et al. 2020). Around Sgr A* are some of the most massive young stellar clusters in the Galaxy: Arches and Quintuplet clusters, as well as the GC cluster, in the age range of 2–7 Myr (e.g. Figer et al. 2002). Dense gas is found predominately in the so-called central molecular zone (CMZ), including a twist ring or streams of giant molecular clouds (GMCs) orbiting around Sgr A* (Kruijssen, Dale & Longmore 2015; Ginsburg et al. 2016; Barnes et al. 2017; Libralato et al. 2020). These clouds have a mean temperature of ∼60 − 120 K and is extremely turbulent (e.g. FWHM ∼20 − 50 km s−1 on scales of a few parsecs; Ginsburg et al. 2016). Bursts of star formation are ongoing, but only in the most massive clouds such as Sgr B2 (e.g. Barnes et al. 2017).

Recently, an apparent connection of such activity in the GC to large-scale energetic structures in the Galactic bulge or halo has been revealed in radio and X-ray. Suzaku and XMM–Newton observations have shown hot plasma plumes (or fountains), emanating from the center to the high Galactic latitudes (Nakashima et al. 2013, 2019; Ponti et al. 2019), while in the same general regions, the MeerKAT mapping has revealed exquisite details of radio-emitting features (Heywood et al. 2019). Most noticeable are large-scale diffuse radio lobes, as well as numerous narrow filaments or their bundles, predominately ‘radiating’ away from the most active portion of the CMZ. They are seen not only in the CMZ, where such features were mostly known (e.g. Yusef-Zadeh, Morris & Chance 1984; LaRosa, Lazio & Kassim 2001), but also abundant in central regions relatively far away from the Galactic plane. Some of the radio filaments, observed to be strongly polarized, are clearly synchrotron in origin and are sometimes called non-thermal radio filaments (NTFs). The formation mechanism of such NTFs remains is still greatly uncertain (e.g. Heywood et al. 2019; Yusef-Zadeh & Wardle 2019; Sofue 2020).

To understand the nature of these interstellar structures and their connection to Galaxy-wide energetic phenomena, as well as the activity in the GC, we have constructed a large-scale high-resolution X-ray map of the GC and its interface with the Galactic bulge, based on Chandra/ACIS-I observations. They include new observations to fill various gaps in the coverage of archival ones. With a total 5.6-Ms exposure of the included observations, this combined mapping, covering an ∼2° × 4° field (Fig. 1A), now allows us to detect discrete X-ray sources, study their population, distribution, and other properties; map diffuse hot plasma with minimal confusion; and cross-correlate multiwavelength objects/features. We here report initial results of the mapping, focused on diffuse X-ray-emitting features and their connection to radio structures in the GC/bulge interface.

Figure 1.

Chandra/ACIS-I survey of the GC/bulge field and a multiwavelength montage: panel (A) shows three-colour X-ray intensity mapping – 1–4 (red), 4–6 (green), and 6–9 keV (blue); panel (B) shows the diffuse 1–4 keV map, where detected discrete sources have been excised, although dust scattering haloes of bright X-ray binaries remain; and panel (C) shows the comparison of the diffuse 1–4 keV map (blue) with the MeerKAT 1.3 GHz (red, Heywood et al. 2019) and WISE 22 μm (green) images. All these plots are projected in the Galactic coordinates. Several major X-ray-bright supernova remnants (SNRs) are marked in panel (B), while additional key stellar and interstellar structures referred to in the text are labeled in panel (C).

The rest of this paper is organized as follows. In Section 2, we briefly describe our data reduction and analysis procedure. We present our results in Section 3 and discuss their implications for understanding the nature and origin of various interstellar structures in Section 4. In Section 5, we summarize our main conclusions, as well as results. In this paper, we adopt a distance of the GC to be 8 kpc (1 arcsec = 0.039 pc) and present error bars at the 90 per cent confidence level.

2 OBSERVATIONS AND DATA ANALYSIS

We use all available Chandra observations taken before 2020 and with the ACIS-I at the focus. They are reprocessed with the standard pipeline of the Chandra Interactive Analysis of Observations (ciao; version 4.11 and caldb 4.8.4.1). Additional data reduction steps include the subtraction of the non-X-ray contribution to the data, using the ACIS-I non-X-ray (stowed) background data base, the detection of discrete sources and their excision from the data (Wang 2004). Details of these steps are not important for this work and will be given in a later paper.

We conduct spatial and spectral analyses of individual X-ray features, as well as point-like sources detected in the survey. In this work, we focus on the analysis of an X-ray thread G0.17–0.41 (Fig. 2), which is probably the most outstanding vertical linear X-ray feature detected in the Chandra survey. The presence of this thread was first reported in an XMM–Newton survey of the GC (Ponti et al. 2015). But because of both the limited exposure and spatial resolution of this survey, little specifics can be learned of the thread. In the present Chandra survey, the thread is covered by three observations (Obs # 7157, 19448, and 20111) with a total exposure of 72 ks, all near the axis. We conduct both spatial and spectral analyses of G0.17–0.41. A fit to the count distribution in the thread give a position angle of 2|${_{.}^{\circ}}$|7 west from the Galactic north. The spectral analysis of G0.17–0.41, in particularly, is based on the on-thread data extracted from a 15 × 154 arcsec2 box and a local background estimated in the 100 × 210 arcsec2 box east of the thread. Spectral fitting results obtained are not sensitive to the exact choice of these boxes.

Figure 2.

High-resolution close-up of G0.17−0.41 via the Chandra/ACIS-I 2.5–6 keV event intensity distribution. Detected discrete sources are marked with squares.

3 RESULTS

We detect about 104 sources in the survey field with the local false detection probability P < 10−6. The detection limit depends on the source spectral shape, as well as the local background intensity, exposure, and point spread function, which vary greatly across the field. Nevertheless, we should typically detect sources with luminosities ≳ 1 × 1032 erg s−1, which include black hole/neutron star X-ray binaries, active or quiescent, young pulsars and their associated compact wind nebulae (e.g. Wang, Gotthelf & Lang (e.g. Wang, Gotthelf & Lang 2002a; Hong et al. 2009; Johnson, Dong & Wang 2009; Muno et al. 2009; Jonker et al. 2014). Fainter undetected sources are numerous stars, primarily cataclysmic variables and active binaries (CV + AB; Wang et al. 2002a; Revnivtsev et al. 2009). Their collective contribution should smoothly follow the old stellar distribution in the field and may dominate the ‘diffuse’ (detected source-excised) emission above ∼4 keV (more on this topic in Section 4.4.2). Here we focus on X-ray-emitting interstellar structures.

Figs 1(A) and (B) present an overview of the Chandra survey of the GC/bulge field. While detailed Chandra maps of the generally X-ray bright region along the Galactic plane have been presented previously (e.g. Wang et al. 2002a; Muno et al. 2009), our presentation here is optimized to show large-scale features. In this presentation, two outstanding plumes from the GC stand out, mostly in the 1–4 keV band. These plumes, detected previously in Suzaku and XMM–Newton observations (Nakashima et al. 2013, 2019; Ponti et al. 2019), are now covered fully at the arcsec resolution of Chandra and with little confusion with discrete sources. The appearance of the plumes is affected by foreground soft X-ray absorption, however, which is most severe at |b| ≲ 0|${_{.}^{\circ}}$|3. Here the presence of the plumes is also evidenced by enhanced diffuse 4–6 keV emission, which extends further to higher Galactic latitudes, even in the inner region of the north plume right above Sgr A*. Foreground absorption remains significant in this part of the plume. Accounting for the absorption effect, the X-ray emission associated with the plumes qualitatively decreases in its intensity and softens in its spectrum with the increasing distance away from the Galactic plane. Therefore, the plumes apparently link the high-energy activity within the CMZ to larger scale diffuse soft X-ray structures seen at |b| ≳ 1° (Fig. 1B), where the X-ray absorption is weak. Compared to the north plume, the south plume is substantially more collimated and centrally filled in soft X-ray and extends further away from the GC into the Galactic bulge.

Fig. 1(C) presents a multiwavelength comparison of various relevant structures seen in the GC/bulge field covered by the MeerKAT survey (Heywood et al. 2019). The southern radio lobe, newly discovered by the survey, encloses the south X-ray plume well. In contrast, the radio and X-ray structures on the Galactic north side are not correlated. The north X-ray plume is systematically offset to the Galactic east from the northern radio lobe, which is sometimes called Galactic Centre Lobe (GCL; e.g. Law et al. 2009). The GCL appears to be closely linked to a distinct warm dust ‘double helix’ nebula (Morris, Uchida & Do 2006), seen here in the 22- μm emission. The implications of these findings are discussed in Section 4.4.

Various small-scale X-ray-emitting linear features or threads are present in the Chandra data, especially in the inner GC region (e.g. Wang, Lu & Lang 2002b; Lu, Wang & Lang 2003; Johnson et al. 2009; Zhang et al. 2020). Some of these X-ray threads are associated with NTFs. While a systematic identification and analysis of such associations is ongoing, Figs 24 highlight an outstanding X-ray thread, G0.17−0.41, which is located quite off the CMZ and is nearly vertically oriented, with respect to the Galactic plane. Most interestingly, G0.17−0.41 is embedded well within an NTF, which appears to be locally swelled (Fig. 3B). Their physical association is thus evident.

Figure 3.

Adaptively smoothed Chandra/ACIS-I 2.5–6 keV intensity image of the G0.17−0.41 (panel A) and comparison with the MeerKAT 1.3-GHz image (panel B), where the X-ray intensity contours are at 3σ, σ5, 9σ, 17σ, and 33σ above the background.

The presence of G0.17−0.41 is most significant in the 2.5–6 keV band and does not seem to be related to any detected discrete sources, as shown in Fig. 2, which is nearly at the full resolution of Chandra. The thread is very narrow in its middle part and becomes dimer and wider on both sides. The image presented in Fig. 3 is adaptively smoothed to reach a non-X-ray background-removed signal-to-noise ratio of >4 across a larger field of view that is linearly twice larger than that covered by Fig. 2. The total length of the thread is about 2.6 arcmin or lx ≈ 6 pc at the distance of the GC. In both Figs 2 and  3, G0.17−0.41 shows substructures, although individually they are only marginally significant at ≲ 2σ levels. Fig. 4(A) shows a surface brightness profile from a cut perpendicular to the full length of G0.17−0.41 and is adaptively grouped to achieve S/N > 5 for each bin. We use a similar profile for a cut through the inner half of the length to estimate the effective width of the thread as ∼1.5 arcsec (FWHM).

Figure 4.

Panel (A): diffuse 2.5–6 keV intensity, averaged over the length of G0.17−0.41, as a function of the vertical distance from it; the positive direction is toward the Galactic (slightly south) west. The horizontal line marks the average off-G0.17−0.41 background intensity. Panel (B): ACIS-I spectrum of G0.17−0.41, together with the best-fitting power-law model (Table 1).

With a net count rate of (5.0 ± 0.4) × 10−3 counts s−1, the ACIS-I detection of G0.17−0.41 allows for simple characterizations of its spectral shape. Table 1 lists model parameters: a power law and an optically-thin thermal (APEC), assuming collisional ionization equilibrium. Solar metal abundances (Anders & Grevesse 1989) are assumed for both the hot plasma. Fig. 4(B) shows the power law fit to the X-ray spectrum of this thread, which is only slightly more favored statistically than the plasma characterization with the abundances fixed to the solar values (first column under APEC in Table 1). The APEC model fit can be marginally improved, if the abundance is allowed to vary, which suggests an abundance of <0.51 solar (at the 95 per cent confidence; second column under APEC in Table 1). This abundance is somewhat lower than what is expected for the GC, but is consistent with the value obtained for G359.55+0.16 (Yamauchi et al. 2014, see Section 4.1 for more discussion). The large X-ray-absorbing column density NH inferred from either model fits is consistent with the GC location of G0.17−0.41.

Table 1.

X-ray spectral fitting results of G0.17−0.41.

ParametersPower lawAPEC
NH(1022 cm−2)11 (5.4–17)5.2 (3.8–7.0)9.0(5.1–15)
Γ or kT (keV)3.0 (1.6–4.5)>7.93.2(1.5–14)
Z (solar abundance)1<0.51
K5.7 × 10−41.5 × 10−44.9 × 10−4
χ2/dof36.0/3239.4/3235.8/31
fx(10−13 erg s−1 cm−2)1.31.61.2
Lx(1033 erg s−1)2.61.82.1
ParametersPower lawAPEC
NH(1022 cm−2)11 (5.4–17)5.2 (3.8–7.0)9.0(5.1–15)
Γ or kT (keV)3.0 (1.6–4.5)>7.93.2(1.5–14)
Z (solar abundance)1<0.51
K5.7 × 10−41.5 × 10−44.9 × 10−4
χ2/dof36.0/3239.4/3235.8/31
fx(10−13 erg s−1 cm−2)1.31.61.2
Lx(1033 erg s−1)2.61.82.1

Notes. The parameter uncertain ranges are all presented at the 90 per cent confidence level. K is the model normalization. For the APEC modelling, two fits are conducted: one assumes that the metal abundance is solar, and the other lets it to be fitted. For each model fit, both flux fx and luminosity Lx are inferred from the best-fitting model in the 2–10 keV band. Lx further corrects for the absorption (i.e. NH = 0).

Table 1.

X-ray spectral fitting results of G0.17−0.41.

ParametersPower lawAPEC
NH(1022 cm−2)11 (5.4–17)5.2 (3.8–7.0)9.0(5.1–15)
Γ or kT (keV)3.0 (1.6–4.5)>7.93.2(1.5–14)
Z (solar abundance)1<0.51
K5.7 × 10−41.5 × 10−44.9 × 10−4
χ2/dof36.0/3239.4/3235.8/31
fx(10−13 erg s−1 cm−2)1.31.61.2
Lx(1033 erg s−1)2.61.82.1
ParametersPower lawAPEC
NH(1022 cm−2)11 (5.4–17)5.2 (3.8–7.0)9.0(5.1–15)
Γ or kT (keV)3.0 (1.6–4.5)>7.93.2(1.5–14)
Z (solar abundance)1<0.51
K5.7 × 10−41.5 × 10−44.9 × 10−4
χ2/dof36.0/3239.4/3235.8/31
fx(10−13 erg s−1 cm−2)1.31.61.2
Lx(1033 erg s−1)2.61.82.1

Notes. The parameter uncertain ranges are all presented at the 90 per cent confidence level. K is the model normalization. For the APEC modelling, two fits are conducted: one assumes that the metal abundance is solar, and the other lets it to be fitted. For each model fit, both flux fx and luminosity Lx are inferred from the best-fitting model in the 2–10 keV band. Lx further corrects for the absorption (i.e. NH = 0).

Adopting the optically thin thermal plasma characterization of the X-ray emission from G0.17−0.41, we can estimate the density as |$n \approx n_{\rm e} \approx 20 D_{\rm kpc}(K/V)^{0.5} f_{\rm h}^{-0.5} {\rm cm^{-3}} \approx 15 f_{\rm h}^{-0.5} {\rm cm^{-3}}$|⁠, where Dkpc ≈ 8 is the distance to G0.17−0.41, K is from Table 1, V ≈ 0.017 pc3 is the total volume, assuming a cylinder of diameter ≈0.06 pc (roughly the FWHM of the thread) and the above length estimate lx, and the effective volume filling factor fh of the hot plasma is probably least certain. The total thermal energy and the radiative cooling time are |$E_{\rm th} \sim 3 n_{\rm e} kT f_{\rm h}V\sim 3 \times 10^{47} f_{\rm h}^{0.5}$| erg and |$t_{\rm c} \sim 5 \times 10^{7} f_{\rm h}^{0.5}$| yr, where kT ∼ 10 keV is assumed. The dynamic age of the thread may be approximated as the sound-crossing time of the re-connection region: tdlx/cs ∼ 6 × 103(kT/10 keV)−0.5 yr, over which the radiative cooling is negligible. A similar age may also be obtained if the Alfvén speed, 2.2 × 103n−0.5BmG km s−1, is used.

4 DISCUSSION

Here we first compare G0.17−0.41 with other most relevant X-ray threads and NTFs observed in the GC, providing both a brief review of existing work and a context for later discussion, then explore potential origins of G0.17−0.41 and another similar X-ray thread/NTF association G359.55+0.16, and finally focus on the potential role of magnetic field re-connection events in regulating the ISM in the GC and beyond.

4.1 Comparison with other similar X-ray and/or radio features

Among the X-ray/NTR associations observed in the GC, G0.17−0.41 most resembles G359.55+0.16. Discovered by Lu et al. (2003) in an early Chandra GC survey, G359.55+0.16 also shows a linear X-ray structure of a length similar to G0.17−0.41’s, contains no point-like source, and is located quite far away from and oriented nearly perpendicular to the Galactic plane. The associated NTF represents a strand of flux tubes, consisting of two main parallel bundles. Radio polarization measurements show that the average magnetic field is orientated along these bundles (Yusef-Zadeh, Wardle & Parastaran 1997). With a deep exposure of Suzaku, Yamauchi et al. (2014) show that the X-ray spectrum of G359.55+0.16 exhibits a distinct Fe xxv-He α 6.7-keV emission line with an equivalent width of |$850^{+400}_{-390}$| eV. The spectrum can be well characterized by an optically-thin thermal plasma with the temperature and metal abundance of |$4.1^{+2.7}_{-1.8}$| keV and |$0.58^{+0.41}_{-0.32}$| solar, which are qualitatively similar to those obtained for G0.17−0.41 (Section 3). Because of these similarities between G0.17−0.41 and G359.55+0.16, we explore their origins together in the following. To do so, we also need to place them in a large context of the NTFs observed in the GC.

The NTFs, although most of which do not have any X-ray counterparts, share some common morphological characteristics, as well as radio spectral and polarization properties. One of the striking collective characteristics of the NTFs is their predominantly radial orientation and concentration in the GC. They appear mostly radiating from within ∼80 pc of Sgr A*, tracing the poloidal magnetic field component or the magnetosphere of the GC (Morris 2006). The component may represent protogalactic field that has been concentrated over the history of the Galaxy by mass inflow (Chandran, Cowley & Morris 2000) and/or outflows themselves under the influence of the Coriolis force and magnetic re-connection (e.g. Hanasz, Otmianowska-Mazur & Lesch 2002). Analyses of existing observations show that the NTFs likely represent the periphery of the ongoing outflows from the GC (Heywood et al. 2019; Yusef-Zadeh & Wardle 2019). Incidentally, G0.17−0.41 is at the east rim of the southern radio lobe, while G359.55+0.16 is at the western boundary of the north plume of X-ray emission. G359.55+0.16 also seems to connect between the CMZ and a warm dust emission region, AFGL 5376 (Fig. 1C; Uchida et al. 1994). Indeed, the NTFs all appear to be undergoing an interaction with a molecular cloud (e.g. Morris 2015).

The most tight constraint on the magnetic field strength inside NTFs relies on a simple dynamical argument, originally proposed by Yusef-Zadeh & Morris (1987a). They argue that the absence of any significant bending or distortion of an NTF against the strong ram pressure of its associated molecular cloud indicates a field strength of ≳ 1 mG. The fattening of the NTF associated with G0.17−0.41 suggests an overpressure produced by the X-ray thread. Adopting the optically thin thermal plasma characterization of this X-ray emission from G0.17−0.41 (Section 3) and assuming the balance between the outward thermal pressure and inward magnetic tension, we estimate the required total vertical field strength |$B_{\perp } \sim 3 f_{\rm h}^{-0.25}$| mG. This value is quite large and may not be provided by the fattening of the NTF alone. Other contributions to the confinement can also be important. We will come back to this point in Section 4.3.

The above discussion, together with our new results, now provide a context for us to explore the origins of such X-ray thread/NTF associations as G0.17−0.41 and G359.55+0.16.

4.2 Pulsar wind nebula interpretation

Based on the satisfactory power law characterization of the X-ray spectrum (Table 1), one may suggest that G0.17−0.41 could represent a pulsar wind nebula (PWN), confined by a strongly magnetized medium (Wang et al. 2002b; Bandiera 2008; Barkov et al. 2019; Barkov & Lyutikov 2019). In this scenario, the X-ray thread would be interpreted as the synchrotron emission of relativistic electrons (and positrons), while streaming mostly along magnetic lines. Indeed, the scenario has been proposed for multiple X-ray threads in the GC, with or without radio counterparts (Wang et al. 2002b; Johnson et al. 2009; Zhang et al. 2020), although none of these cases has been confirmed by the detection of a putative pulsar. The scenario is considered to be most plausible when an X-ray thread contains a point-like X-ray source (i.e. the putative pulsar) and/or exhibits evidence for synchrotron cooling (Wang et al. 2002b; Gaensler et al. 2004; Park et al. 2005; Wang, Lu & Gotthelf 2006; Zhang et al. 2020).

One of – if not the – best candidates for such a PWN interpretation in the GC is G0.13−0.11, which shows a distinct curved linear X-ray morphology and borders a bow-shaped radio feature (Wang et al. 2002b; Zhang et al. 2020). The properties of the X-ray and radio emissions and their morphological configuration strongly suggest that G0.13−0.11 results from a moving pulsar with the propagation of its wind material strongly regulated by surrounding organized magnetic field (Wang et al. 2002b; Zhang et al. 2020). Well embedded in the X-ray thread is a point-like source of L2 − 10 keV = 8.4 × 1032 erg s−1. The X-ray properties of both the source and the thread are consistent with the PWN interpretation (Wang et al. 2002b; Zhang et al. 2020). G0.13−0.11 is also apparently associated with a point-like TeV source, detected at a significance of 5.9σ in the HESS GC survey (Aharonian et al. 2006; H. E. S. S. Collaboration et al. 2018). The inverse-Compton origin of the TeV emission matches well with the synchrotron interpretation for the X-ray emission of G0.13−0.11. It may thus be considered as a reference source for us to evaluate the likelihood for G0.17−0.41 or G359.55+0.16 to be a PWN.

We find that an application of the PWN (or SNR) scenario to G0.17−0.41 or G359.55+0.16 meets several major difficulties:

  • There is no evidence for a point-like source as the putative pulsar in either X-ray thread. Assuming the X-ray versus spin-down luminosity (Lsd) correlation, logL2 − 10 keV = 1.34 logLsd − 15.34 (Possenti et al. 2002), one would expect a quite powerful pulsar with Lsd ∼ 3 × 1036 erg s−1. Such a source, if present in G0.17−0.41 for example, would be easily detected because its overall X-ray luminosity is about the same as that of G0.13−0.11.

  • The power law photon index of G0.17−0.41 (though with a large statistical uncertainty range; Table 1) is not typical for a PWN and is steeper than that of G0.13−0.11, for example, at the 95 per cent confidence level. For G359.55+0.16, the PWN scenario for its X-ray emission can be completely ruled out because the detection of the strong He-like Fe 6.7-keV emission line, which arises from an optically-thin thermal plasma.

  • The PWN scenario also has the difficulty to explain the X-ray thread length of G0.17−0.41. The X-ray synchrotron cooling time scale is |$t_{{\rm s},x} \sim 1 {\rm yr}B_{\rm mG}^{-1.5}E_{\rm keV}^{-0.5}$|⁠, where EkeV is the photon energy. Thus ts, x is too small to be compared with the time scale for the electrons to move cross the thread, which suggests that the particle acceleration cannot be in a point-like source, unless BmG is substantially less than 1, which is not favoured for such a bright NTF.

  • Because of the long radio-synchrotron cooling time scale, one may expect that a pulsar with a typical moving speed (a few × 102 km s−1) should leave an extended radio remnant or PWN of a fan-out morphology extending only on one side of the X-ray thread and at the opposite direction of the pulsar’s proper motion, as is the case for G0.13−0.11, but not for G0.17−0.41. The fattening of the NTF around G0.17−0.41 is rather symmetric and has a half-width of ≲ 1 pc, as seen in the 1.3-GHz image. If the radio synchrotron cooling dominates, then we may constrain the vertical speed of the pulsar crossing the NTF filament to be |$\sim 5 {\rm km\,s^{-1}} B_{\rm mG}^{3/2}$|⁠. This speed is much smaller than the velocity expected for a pulsar or even for a source with a relative orbital motion, which is common in the GC. Clearly, a small BmG value would make this problem worse. Thus, it is difficult to understand such a nearly perfect and symmetric alignment between the X-ray thread and the NTF, unless they are physically associated and are produced by the flux tube itself. A similar argument can also be made for G359.55+0.16.

  • There is no indication of any TeV emission enhancement at the location of G0.17−0.41 or G359.55+0.16 in the HESS GC survey (Aharonian et al. 2006). The diffuse TeV background at the locations of these threads, away from the CMZ, is much low, in comparison to the location of G0.13−0.11. The lack of the HESS detection indicates that the flux of a point-like source would be at least about 30 lower than that of G0.13−0.11. This limiting TeV flux needs to be explained if the X-ray flux of G0.17−0.41, for example, is to be interpreted as synchrotron emission.

Finally, could G0.17−0.41 and G359.55+0.16 be SNRs regulated by their strongly magnetized environments? Both, more so G0.17−0.41, are located quite off the CMZ, where recent star formation occurs. Correspondingly, the probability should be small to detect a young pulsar or an SNR at such a location, especially in coincidence with a probably distinct flux tube of strongly enhanced magnetic field. Moreover, our estimated thermal energy of G0.17−0.41 (Section 3; similar for G359.55+0.16) is far too small to be consistent with that expected for a typical SNR. Its age should not matter much here, because the thermal radiative cooling time is long. This energy estimate should also hardly depend on the thermal emission assumption. Fundamentally the pressure and volume of the X-ray thread are too limited to be consistent with an SNR interpretation. However, it may be possible that G0.17−0.41 or G359.55+0.16 represents just a small piece of an SNR, ejected from the CMZ along a magnetic flux tube, although the plausibility of this scenario is yet to be carefully examined, both theoretically and observationally.

4.3 Magnetic field re-connection scenario

This scenario provides a unified explanation for the X-ray thread/NTF associations of G0.17−0.41 and G359.55+0.16, in terms of their formation and emission. We consider that they represent interstellar magnetic field strands. Various models have been proposed for the formation of such structures, mostly related to the outflow or wind from the GC (e.g. LaRosa et al. 2004; Sofue, Kigure & Shibata 2005; Banda-Barragán et al. 2016; Yusef-Zadeh & Wardle 2019).

In the cometary model (e.g. Shore & LaRosa 1999), for example, NTFs are magnetized wakes generated by the interaction of molecular clouds with a global GC wind. Naturally, field lines of opposite directions are expected to meet at the rear of such clouds, forming re-connection-prone wakes – a situation similar to what happens in the heliospheric tail (e.g. Lazarian & Desiati 2010; Banda-Barragán et al. 2016). The lines could also be twisted, e.g. due to the spin of the molecular clouds and/or the different rotation of the CMZ (Yusef-Zadeh & Morris 1987b; Mezger, Duschl & Zylka 1996; LaRosa et al. 2004). The twisting leads to the formation of the tightened ropes of enhanced magnetic field flux tubes and eventually to the magnetic re-connection (e.g. Titov & Démoulin 1999; Tsap et al. 2020).

Re-connection (more precisely breaking and rejoining) events take place at interface between ionized or partially ionized gases with opposite magnetic field lines (e.g. Lazarian & Vishniac 1999; Bicknell & Li 2001; Zweibel & Yamada 2009; Lazarian et al. 2020). In such an event, the fractal tearing instability of field lines leads to the formation of plasmoids (magnetic islands) in the re-connection region (e.g. Furth, Killeen & Rosenbluth 1963). Both observations and simulations have shown that the local vertical field strength in re-connection regions become comparable to or even slightly greater than the large-scale mean field (e.g. Ieda et al. 1998; Lu et al. 2020). The total pressure inside a plasmoid is typically higher than that in the background because magnetic tension force acts to compress the plasmoid. The situation here is analogous to the sun, where the energy density in solar flare regions can be much stronger than the average value in the solar corona (e.g. Toriumi & Wang 2019). Magnetic re-connection is also a dynamic process. Electric field generated by the changing magnetic field and the collision among the plasmoids further accelerate particles, producing cosmic rays (CRs; e.g. Bicknell & Li 2001). Direct observational evidence for re-connection events is abundant in the solar corona and space weather (e.g. Toriumi & Wang 2019), but so far little in the ISM.

The production of diffuse hot plasma via re-connection in the interstellar space has been predicted theoretically, in analogy to the heating of the solar corona (e.g. Heyvaerts, Norman & Pudritz 1988; Raymond 1992; Tanuma et al. 2003; Florido-Llinas, Nieves-Chinchilla & Linton 2020), and has been used to explain diffuse X-ray emission observed in nearby galaxies (e.g. Weżgowiec et al. 2020), as well as in our Galaxy (e.g. Tanuma et al. 1999). MHD simulations in interstellar contexts (Tanuma et al. 1999, 2003; Hanasz et al. 2002) show that re-connection events can occur in expanding shells (e.g. due to Parker or magnetic buoyancy instability) and, under the influence of Coriolis force, produce ‘helical tubes’ of magnetic field lines, as well as in wakes of molecular clouds interacting with a global GC wind (e.g. Shore & LaRosa 1999; Banda-Barragán et al. 2016). This mechanism seems to naturally explain the close X-ray/NTF associations, as seen in G0.17−0.41 and G359.55+0.16, as well as their morphological and spectral properties. The 6.7-keV iron-line emission, in particular, is also most naturally explained by an interstellar re-connection event (Lu et al. 2003). G0.17−0.41 probably represents an even stronger case than G359.55+0.16 as a re-connection event in an interstellar magnetic field rope. Much of the X-ray emission may arise from plasmoids as indicated by the discrete substructures of G0.17−0.41 (Section 3). The small-scale radially inward tension force of magnetic loops around these plasmoids may largely responsible for the confinement of G0.17−0.41 and G359.55+0.16. This small-scale plasmoid confinement is in addition to the global inward tension of the curved field lines of the NTF around G0.17−0.41, as discussed at the end of Section 4.1.

Furthermore, a reconnection event is likely triggered by external overpressure of the ambient medium (e.g. due to the cloud–cloud collision and/or twisting of a magnetic rope), although no direct evidence is available for such a dynamical process around the G0.17−0.41/NTF association. Our argument for the reconnection scenario is that it provides a consistent interpretation for the X-ray properties of G0.17−0.41 and G359.55+0.16, as well as their physical association with the NTFs. While the X-ray threads likely represent the ongoing or recent reconnection sites, the associated NTFs are the illuminated parts of the underlying magnetic flux tubes, which could extend to much larger scales and even connect to features on the opposite side of the Galactic plane (e.g. Fig. 1C). Therefore, the X-ray thread/NTF associations such as G0.17−0.41 and G359.55+0.16 may provide potentially excellent laboratories to test our understanding of the re-connection astrophysics in the ISM. They allow for spatially resolved studies of the hot plasma production, as well as the particle acceleration and magnetic field structure, in the re-connection zone.

In the re-connection scenario, the X-ray emission directly traces the heated thermal plasma. According to Tanuma et al. (2003), the plasma is heated to a characteristic temperature of |$kT \sim 10 n^{-1}B_{\rm mG}^{2}$| keV. The density n of the gas accompanied with the magnetic field can vary greatly from one location to another in the GC. But most ubiqious is diffuse warm molecular gas of n ≲ 100 cm−3 (Oka et al. 2019) or more tenuous diffuse hot plasma. The estimated n and B values (Section 3 and S 4.1) are reasonably consistent with those required to produce kT ∼ 10-keV plasma via the re-connection. Using the results given in Section 3, we also infer the plasma heating rate of ∼Eth/td ∼ 2 × 1036 erg s−1. The total energy generation rate of the region should be larger than this value, accounting for the energy in the plasma motion and in the particle acceleration.

Different radio/X-ray associations may represent re-connection events at different evolutionary stages (Tanuma et al. 1999, 2003; Hanasz et al. 2002). The swelled radio morphology of G0.17−0.41 indicates that the X-ray-emitting material is dynamically important and is still shaping the surrounding magnetic field topologically, which may represent an early event stage (e.g. fig. 6 in Tanuma et al. 1999). In comparison, G359.55+0.16 may be a manifest of a later-stage re-connection event. The NTF morphology around this thread is topologically different from that associated with G0.17−0.41. G359.55+0.16 appears at the closest (possibly crossing) part of the two major bundles of the NTF. Off the X-ray thread, they branch out. The X-ray thread is associated with only one of them at present. There is no evidence for any significant local inflation of the associated NTF in G359.55+0.16, suggesting that magnetic field tension dominates the hot plasma pressure. Some of the plasma may have already escaped from the re-connection region, at least from parts of the NTF. Therefore, the re-connection scenario provides a consistent explanation of the X-ray properties of G0.17−0.41 and G359.55+0.16, as well as their association with the distinct NTFs.

Of course, more conclusive evidence will be needed to establish the re-connection scenario. Future more sensitive X-ray and radio observations will enable spatially resolved and even time-dependent spectroscopy of such re-connection phenomena. In particular, a deeper X-ray observation of G0.17−0.41 will better constrain the presence of of the 6.7-keV line, as expected for the re-connection scenario. Multiband high-resolution radio observations with polarization measurements will be extremely useful to explore the magnetic field structure, as well as the propagation properties of the CR electrons along the associated NTFs. Dedicated theoretical studies, including computational simulations, are needed to improve our understanding of the magnetic re-connection in various interstellar environments. For now, we can examine potential connections of re-connection events to other NTFs and to larger scale interstellar structures observed in the GC.

In addition to G0.17−0.41 and G359.55+0.16, some of the other NTFs in the GC may also arise from re-connection events. Their X-ray emission may just be too weak and/or diffuse to be detected individually. As shown in MHD simulations (Tanuma et al. 2003), the energy release profile of a re-connection event is strongly model-dependent and typically shows an overall fast-rising slow-decay pattern. The half-peak duration is on the order of 10 times the sound-crossing time td. Many re-connection events may not be able to heat plasma to high enough temperature for significant X-ray emission. Identified radio/X-ray filament associations likely represent only the tip of the ‘re-connection iceberg’. Its true magnitude may be better traced by NTFs. CR electrons diffuse along the field lines at the Alfvén speed over the characteristic synchrotron cooling time-scale ts, r (Morris 2006; Thomas, Pfrommer & Enßlin 2020). One then expects to detect radio emission off the re-connection sites on the spatial scale of 44(nBmGν GHz)−0.5 pc, which is comparable to the observed length of the longest filaments in the GC. The resultant NTFs should in general be more extended and live longer than the X-ray threads, again consistent with the observations. Therefore, we may reasonably assume that these re-connection events are intimately linked to energetic outflows as traced by the large-scale interstellar structures in the GC.

4.4 Re-connection and global ISM structures

To explore the potential links of magnetic field re-connection to interstellar structures observed in various wavelengths, we first need to determine whether or not and how they may be physically associated in the GC. Here we focus on the X-ray plumes and the radio lobes, as we have described in Section 3.

4.4.1 Origin of the X-ray plumes and their connection to the radio lobes

In the GC field, projection effect can be significant. Particularly relevant here is the GCL. It has recently been proposed to be a foreground H ii region (Nagoshi et al. 2019; Tsuboi et al. 2020). This proposition appears consistent with various observations: the detection of the associated H α emission (Law et al. 2009), the small systemic velocity gradient (between −4 and +10 km s−1) observed in radio recombination lines (Law et al. 2009; Sofue 2013; Alves et al. 2015; Nagoshi et al. 2019), and the lack of an X-ray counterpart (Fig. 1C). Even the brightest (western) part of the lobe shows no enhanced X-ray emission in either 1–4 or 4–6 keV band (Fig. 1A). While the X-ray plume seems to be rooted in the most active portion of the CMZ, including all the three massive-and-young star clusters, the GCL is offset to the Galactic west and is not associated with any evident energetic source in the GC (Fig. 1C). An ongoing multiwavelength study is placing tight constraints on the line-of-sight location of this H ii region and on a nearby young stellar cluster that may be the source for the ionizing radiation (Benjamin, in preparation and private communications). One may speculate that the H ii region is also responsible for much of the angular size and temporal pulse broadening observed in GC radio pulsars, especially the transient magnetar SGR J1745−2900, which is only 2.5 arcsec from Sgr A* (Bower et al. 2014; Sicheneder & Dexter 2017). Indeed, a single H ii region modeling of the broadening measurements places the radio scattering screen at 1.5–4.8 kpc away from Earth, probably in the Scutum spiral arm (Sicheneder & Dexter 2017). Therefore, we consider the GCL to be a foreground of the GC and is not physically related to the X-ray plumes.

The north X-ray plume (Fig. 1) does not seem to be enclosed by any coherent radio structure and may represent an aged outflow or fountain of hot plasma, which is apparently powered by some or all of the young massive stellar clusters. Unlike in the Galactic disc, the dynamics of the outflow in the GC may be strongly influenced by the motion of Arches and Quintuplet clusters, which orbit around Sgr A* with period ∼106 yr, shorter than their ages. In contrast, the physical association of the south X-ray plume with the newly discovered MeerKAT southern lobe (Heywood et al. 2019) appears secure. The radio lobe extends to the southern edge of the field in Fig. 1(C) and wraps around the X-ray plume well. The X-ray emission from the plume is clearly dominated by optically-thin thermal plasma in an overionized state (Nakashima et al. 2013). This state, together with the strongly centrally filled and collimated morphology, suggests that the plume represents a recent blowout from the high-pressure GC region, resulting in the overcooling of the enclosed plasma. The collimation of the hot plasma and its association with the radio lobe further indicate that the confinement of the plasma is largely due to magnetic field lines, which are likely anchored in the CMZ and are twisted because of its orbital motion.

4.4.2 Heating of diffuse hot plasma in the GC

Diffuse hot plasma is best traced by X-ray emission. However, there is an ongoing debate as to the nature of the X-ray emission from the GC. Relevant studies have largely been focused on the nature of the prevalent Fe xxv-He α (6.7-keV) line emission (e.g. Nishiyama et al. 2013; Koyama 2018). This line emission has been interpreted mostly as a superposition of unresolved point sources (CV + AB; Wang et al. 2002a), which is largely confirmed in a 1-Ms-deep Chandra exposure of an inner bulge field (l ∼ 0|${_{.}^{\circ}}$|08, b ∼ −1|${_{.}^{\circ}}$|42; Revnivtsev et al. 2009). But extending this interpretation to the GC has been challenged vigorously (e.g. Nishiyama et al. 2013). Suzaku observations have shown that the specific line emission per stellar mass in regions outside |b| ∼ 0.5 is consistent with the Galactic bulge/ridge value, but increases (by a factor of up to ∼5) with the decreasing |b| in the inner GC region. Metallicity increase toward the GC could be a factor, although the iron abundance there is on average only up to about twice solar (e.g. Nogueras-Lara et al. 2018). Alternatively, the enhanced line emission, as well as the rather flat X-ray spectral continuum (e.g. Wang et al. 2002a), may be largely due to the presence of diffuse plasma with a characteristic temperature of ∼7 keV. This scenario, if confirmed, would have profound implications for our understanding the heating and confining mechanisms in the GC (e.g. Nishiyama et al. 2013); a huge amount of magnetic energy (∼1055 erg) would be expected in the central ∼300 pc, comparable to the kinetic energy associated with the rotation of the gas in the CMZ (Ponti et al. 2015). There are also lines of evidence for a diffuse soft X-ray component that is most likely due to plasma with a characteristic temperature probably at ∼1 keV (e.g. Ponti et al. 2015; Yamauchi et al. 2018). Such a characterization, however, is sensitive to the spectral model assumed in the multiple-component decomposition of the observed X-ray emission.

The presence of a substantial diffuse hot plasma component in the GC now becomes clear, with the detection of the distinct interstellar structures such as the north and south X-ray plumes, as revealed in Section 3 (see also Nakashima et al. 2013, 2019; Ponti et al. 2019). These plumes apparently originate in the most active portion of the CMZ and provide links to global energetic structures in the bulge and halo of our Galaxy.

Much of the diffuse X-ray emission from the plumes and from the GC in general could be produced by re-connection events. The magnetic pinch confinement of super hot plasma in individual flux tubes helps to mitigate the energy supply problem (e.g. Nishiyama et al. 2013), as it is not freely escaping from the GC. As discussed in Section 4.3, most of such hot plasma-filled flux tubes may not be detected individually with existing data. Because of the long radiative cooling time, the plasma will eventually expand and cool down, converting some of the thermal energy to kinetic energy. Therefore, one may expect that this unresolved plasma contribution should show a softer X-ray spectral shape than that in such confined threads as G0.17−0.41 and probably G359.55+0.16, consistent with the observed spectral properties of the diffuse X-ray emission from our Chandra survey field. Therefore, re-connection events may collectively contribute significantly to the X-ray emission observed in the GC.

One may also expect that re-connection events should occur more frequently and vigorously in and near the CMZ, where the overall interstellar pressure is more extreme. This is not only because of strong energy feedback from massive stars and possibly Sgr A*, but also due to differential orbital and turbulent motion of magnetized ISM. Indeed, many NTFs have been observed in the CMZ. Some of them do have X-ray counterparts (e.g. Wang et al. 2002b; Lu et al. 2003; Johnson et al. 2009; Zhang et al. 2020), although detection or identification of individual re-connection events is more difficult because of the confusion with other X-ray-emitting structures and the overall enhanced background.

4.4.3 Re-connection-regulated global ISM

A key ingredient for the re-connection to be an important interstellar process is the presence of a sufficiently large collective interface where magnetic fields of opposite directions collide. Such a dynamical event can be realized in various ways in a strongly turbulent medium, driven by differential rotation, as well as recent energetic feedback from massive stars and perhaps Sgr A*. High-velocity collisions can occur not only in the CMZ, but also around it. Energetic outflows, as evidenced by the radio lobes and X-ray plumes, naturally stretch mostly toroidal magnetic field lines in the CMZ to high Galactic latitudes, leading to the formation neutral sheets between colliding media with different magnetic field orientations. Strongly enhanced flux tubes may also form out of swept-up material, e.g. due to collisions between poloidal magnetic fields and a jet-like outflow from Sgr A* (e.g. Sofue & Fujimoto 1987) or dense clouds (Serabyn & Morris 1994), stripping of gas off dense gas clouds (e.g. Shore & LaRosa 1999; Banda-Barragán et al. 2016), and the presence of various other obstacles (e.g. individual stellar wind bubbles Yusef-Zadeh & Wardle 2019). These processes and hence associated re-connection events tend to occur at the outer boundaries of global outflows. As pointed out in Section 4.1, the G0.17−0.41 and G359.55+0.16 and their associated NTFs indeed appear at the rims of the radio lobes or X-ray plumes. This should motive future observations of gas associated with the southern radio shell, for example, to study its kinematics and potential interplay with various NTFs. Moreover, flux tubes in outflowing plasma can be naturally produced by rotational motion of associated molecular clouds in the GC, as shown in simulations (e.g. Mezger et al. 1996). In fact, recent 3D simulations show that turbulence and re-connection in the ISM are inter-connected on all scales: re-connection is an intrinsic part of the turbulent cascade (Lazarian et al. 2020).

In a re-connection event, magnetic field lines break and reconnect, leading to their topological change. The resultant strong curved magnetic tension can slingshot plasma to the Alfvén speed. Such local kinetic energy injection could be an important driving force for high-turbulent motion, as seen in the ISM of the GC. The re-connection also helps to shape the global magnetic field topology of the ISM. As shown in MHD simulations (e.g. Hanasz et al. 2002), the formation of large-scale poloidal magnetic field flux tubes in a turbulent medium requires the re-connection. In analog to solar flares, we may call such re-connection phenomena, driven by underlying interstellar turbulent motion, as ‘galactic flares’ (Sturrock & Stern 1980).

In the extreme environment of the GC, magnetic field may play a key (if not dominant) role in regulating the thermal and dynamical properties, as well as the global structure of the ISM. CRs from re-connection events can transport released magnetic field energy into surrounding medium. Such ‘in-situ’ produced CRs could be responsible for the heating of the ubiqious diffuse warm molecular gas with T ∼ 200 K, as recently revealed by observations of H|$_{3}^+$| absorption lines at 3.5–4.0 μm (Wiener, Zweibel & Oh 2013; Oka et al. 2019). The magnetic field energy may also be stored in global outflows. In this case, later release of the energy via re-connection may greatly influence the emission and the overall energetics of outflows far away from the GC. Lines of evidence for various large-scale bipolar structures around the GC have been revealed across the spectrum: e.g. the so-called Fermi Bubbles in γ-ray (Su et al. 2010), the WMAP ‘haze’ in microwave, and large-scale diffuse emissions in (polarized) radio, mid-IR, and X-ray (e.g. Bland-Hawthorn & Cohen 2003; Carretti et al. 2013). Kinematical information has also been obtained from observations of UV absorption lines in the spectra of stars at various distances, as well as background AGNs, and from HI 21 cm line emission mapping, strongly suggesting that clouds with velocities up to ∼103 km s−1 have been accelerated out from the GC, probably episodically (e.g. Di Teodoro et al. 2018; Ashley et al. 2020). Indeed, the re-connection has been proposed as a mechanism for CR acceleration to explain non-thermal emission from such large-scale structures as the Fermi Bubbles (e.g. Su et al. 2010).

5 SUMMARY AND CONCLUSIONS

We have mapped a 2° × 4° field of the GC/bulge with Chandra/ACIS-I observations, based on a dedicated gap-filling program, as well as archival data. We present maps in various bands after cleaning artifacts (mostly readout streaks and some very extended dust scattering haloes of bright sources) and removing non-X-ray background contributions. Discrete sources detected in individual observations are merged together and are excised from the data to produce diffuse X-ray emission maps. We have conducted some preliminary analysis of these maps and a multiwavelength comparison, mainly with the recent MeerKAT survey of the GC, to provide insights into the nature of various X-ray and radio structures, as well as their relationship. Our main results and conclusions as follows:

  • A pair of diffuse X-ray-emitting plumes, apparently emanating from the GC, are revealed at arcsec resolution and with minimum confusion with discrete sources. The Galactic south plume is centrally filled and collimated and is well enclosed by a radio lobe discovered recently by MeerKAT. This association indicates that the diffuse hot plasma in the plume is confined by a strongly magnetized medium. In contrast, the Galactic north plume is spatially offset from the well-known GCL, consistent with the proposition that this radio lobe is a foreground H ii region. An accurate localization and mapping of this H ii region will be important to the determination of its role as a scatter screen of radio signals, especially those from distant pulsars (Sicheneder & Dexter 2017), and to a comprehensive study of interstellar structures associated with the GC.

  • We highlight the spatial and spectral properties of an X-ray thread G0.17−0.41, which is embedded in a locally swelled radio filament, which is oriented nearly vertically with regard to the Galactic plane and is quite far off the CMZ. The X-ray spectrum of this thread can be characterized by either a power law or an optically thin thermal plasma with a temperature ≳ 3 keV. The combination of the radio and X-ray characteristics suggests that the association represents an ongoing re-connection event of antiparallel magnetic field flux tubes. G0.17−0.41, together with another similar X-ray thread/NTF association G359.55+0.16, will potentially be excellent sites to seriously learn about the interstellar re-connection astrophysics.

  • Major re-connection events tend to occur at outer boundaries of expanding plumes or supershells of dense gas, where strong magnetic fields of different directions collide or get twisted. Detected X-ray thread/NTF associations such as G0.17−0.41 and G359.55+0.16 probably represent only the tip of the re-connection iceberg in the GC. Most re-connection events likely produce too faint and/or too diffuse X-ray emission to be detected individually with existing observations. Re-connection events may play a major role in producing high-temperature optically-thin thermal plasma, accelerating CRs, driving turbulence, and regulating global interstellar structures in and beyond the GC.

These results and conclusions demonstrate the importance of studying high-energy phenomena and processes in the GC and their impact on larger-scale environments over a broad range of spatial or time-scales. It is also clear that the GC is a truly complex system, involving not only interplay among various stellar and interstellar components, plus Sgr A*, but also inflows and outflows, multiple energy sources, as well as heating and cooling mechanisms. Confusion with foreground and potentially background and embedded features can also be serious. A comprehensive study of the GC with this complexity will require truly a multiwavelength approach, together with dedicated theoretical and computer simulations. Ultimately, what we learn from the GC ecosystem and its connection to larger scale structures will provide us with insights into the working of similar extreme regions in other galaxies.

ACKNOWLEDGEMENTS

We thank the referee for useful comments and appreciate the feedback offered by Mark Morris on the work, which was supported by the National Aeronautics and Space Administration through Chandra Award Number GO9-20023X issued by the Chandra X-ray Center, which is operated by the Smithsonian Astrophysical Observatory for and on behalf of the National Aeronautics Space Administration under contract NAS8-03060, and by the ADAP grant NNX17AL67G. This research has made use of data obtained from the Chandra Data Archive and software provided by the Chandra X-ray Center (CXC) in the application packages ciao. The data include all ACIS-I observations taken before 2020 and within the field of view of Fig. 1(A).

DATA AVAILABILITY

The Chandra data as described in Section 2 are available in the Chandra data archive (https://asc.harvard.edu/cda/). Processed data products underlying this paper will be shared on reasonable request to the author.

REFERENCES

Aharonian
F.
et al. ,
2006
,
Nature
,
439
,
695

Alves
M. I. R.
,
Calabretta
M.
,
Davies
R. D.
,
Dickinson
C.
,
Staveley-Smith
L.
,
Davis
R. J.
,
Chen
T.
,
Barr
A.
,
2015
,
MNRAS
,
450
,
2025

Anders
E.
,
Grevesse
N.
,
1989
,
Geochim. Cosmochim. Acta
,
53
,
197

Armillotta
L.
,
Krumholz
M. R.
,
Di Teodoro
E. M.
,
2020
,
MNRAS
,
493
,
5273

Ashley
T.
,
Fox
A. J.
,
Jenkins
E. B.
,
Wakker
B. P.
,
Bordoloi
R.
,
Lockman
F. J.
,
Savage
B. D.
,
Karim
T.
,
2020
,
ApJ
,
898
,
128

Banda-Barragán
W. E.
,
Parkin
E. R.
,
Federrath
C.
,
Crocker
R. M.
,
Bicknell
G. V.
,
2016
,
MNRAS
,
455
,
1309

Bandiera
R.
,
2008
,
A&A
,
490
,
L3

Barkov
M. V.
,
Lyutikov
M.
,
2019
,
MNRAS
,
489
,
L28

Barkov
M. V.
,
Lyutikov
M.
,
Klingler
N.
,
Bordas
P.
,
2019
,
MNRAS
,
485
,
2041

Barnes
A. T.
,
Longmore
S. N.
,
Battersby
C.
,
Bally
J.
,
Kruijssen
J. M. D.
,
Henshaw
J. D.
,
Walker
D. L.
,
2017
,
MNRAS
,
469
,
2263

Bicknell
G. V.
,
Li
J.
,
2001
,
ApJ
,
548
,
L69

Bland-Hawthorn
J.
,
Cohen
M.
,
2003
,
ApJ
,
582
,
246

Bland-Hawthorn
J.
,
Gerhard
O.
,
2016
,
ARA&A
,
54
,
529

Bower
G. C.
et al. ,
2014
,
ApJ
,
790
,
1

Carretti
E.
et al. ,
2013
,
Nature
,
493
,
66

Chandran
B. D. G.
,
Cowley
S. C.
,
Morris
M.
,
2000
,
ApJ
,
528
,
723

Di Teodoro
E. M.
,
McClure-Griffiths
N. M.
,
Lockman
F. J.
,
Denbo
S. R.
,
Endsley
R.
,
Ford
H. A.
,
Harrington
K.
,
2018
,
ApJ
,
855
,
33

Figer
D. F.
et al. ,
2002
,
ApJ
,
581
,
258

Florido-Llinas
M.
,
Nieves-Chinchilla
T.
,
Linton
M. G.
,
2020
,
Sol. Phys.
,
295
,
118

Furth
H. P.
,
Killeen
J.
,
Rosenbluth
M. N.
,
1963
,
Phys. Fluids
,
6
,
459

Gaensler
B. M.
,
van der Swaluw
E.
,
Camilo
F.
,
Kaspi
V. M.
,
Baganoff
F. K.
,
Yusef-Zadeh
F.
,
Manchester
R. N.
,
2004
,
ApJ
,
616
,
383

Ginsburg
A.
et al. ,
2016
,
A&A
,
586
,
A50

H. E. S. S. Collaboration
et al. .,
2018
,
A&A
,
612
,
A9

Hanasz
M.
,
Otmianowska-Mazur
K.
,
Lesch
H.
,
2002
,
A&A
,
386
,
347

Heyvaerts
J.
,
Norman
C.
,
Pudritz
R. E.
,
1988
,
ApJ
,
330
,
718

Heywood
I.
et al. ,
2019
,
Nature
,
573
,
235

Hong
J. S.
,
van den Berg
M.
,
Grindlay
J. E.
,
Laycock
S.
,
2009
,
ApJ
,
706
,
223

Ieda
A.
,
Machida
S.
,
Mukai
T.
,
Saito
Y.
,
Yamamoto
T.
,
Nishida
A.
,
Terasawa
T.
,
Kokubun
S.
,
1998
,
J. Geophys. Res.
,
103
,
4453

Johnson
S. P.
,
Dong
H.
,
Wang
Q. D.
,
2009
,
MNRAS
,
399
,
1429

Jonker
P. G.
et al. ,
2014
,
ApJS
,
210
,
18

Koyama
K.
,
2018
,
PASJ
,
70
,
R1

Kruijssen
J. M. D.
,
Longmore
S. N.
,
2013
,
MNRAS
,
435
,
2598
Kruijssen J. M. D., Longmore S. N., 2013, MNRAS, 435, 2598

Kruijssen
J. M. D.
,
Dale
J. E.
,
Longmore
S. N.
,
2015
,
MNRAS
,
447
,
1059

LaRosa
T. N.
,
Lazio
T. J. W.
,
Kassim
N. E.
,
2001
,
ApJ
,
563
,
163

LaRosa
T. N.
,
Nord
M. E.
,
Lazio
T. J. W.
,
Kassim
N. E.
,
2004
,
ApJ
,
607
,
302

Law
C. J.
,
Backer
D.
,
Yusef-Zadeh
F.
,
Maddalena
R.
,
2009
,
ApJ
,
695
,
1070

Lazarian
A.
,
Desiati
P.
,
2010
,
ApJ
,
722
,
188

Lazarian
A.
,
Vishniac
E. T.
,
1999
,
ApJ
,
517
,
700

Lazarian
A.
,
Eyink
G. L.
,
Jafari
A.
,
Kowal
G.
,
Li
H.
,
Xu
S.
,
Vishniac
E. T.
,
2020
,
Phys. Plasmas
,
27
,
012305

Libralato
M.
,
Fardal
M.
,
Lennon
D.
,
van der Marel
R. P.
,
Bellini
A.
,
2020
,
MNRAS
,
497
,
4733

Lu
F. J.
,
Wang
Q. D.
,
Lang
C. C.
,
2003
,
AJ
,
126
,
319

Lu
S.
,
Angelopoulos
V.
,
Artemyev
A. V.
,
Pritchett
P. L.
,
Sun
W. J.
,
Slavin
J. A.
,
2020
,
ApJ
,
900
,
145

Mezger
P. G.
,
Duschl
W. J.
,
Zylka
R.
,
1996
,
A&AR
,
7
,
289

Morris
M. R.
,
2015
,
Manifestations of the Galactic Center Magnetic Field
,
Springer International Publishing
,
Switzerland
, p.
391

Morris
M.
,
Serabyn
E.
,
1996
,
ARA&A
,
34
,
645

Morris
M.
,
Uchida
K.
,
Do
T.
,
2006
,
Nature
,
440
,
308

Muno
M. P.
et al. ,
2009
,
ApJS
,
181
,
110

Nagoshi
H.
et al. ,
2019
,
PASJ
,
71
,
80

Nakashima
S.
,
Nobukawa
M.
,
Uchida
H.
,
Tanaka
T.
,
Tsuru
T. G.
,
Koyama
K.
,
Murakami
H.
,
Uchiyama
H.
,
2013
,
ApJ
,
773
,
20

Nakashima
S.
,
Koyama
K.
,
Wang
Q. D.
,
Enokiya
R.
,
2019
,
ApJ
,
875
,
32

Nishiyama
S.
et al. ,
2013
,
ApJ
,
769
,
L28

Nogueras-Lara
F.
et al. ,
2018
,
A&A
,
620
,
A83

Oka
T.
,
Geballe
T. R.
,
Goto
M.
,
Usuda
T.
,
Benjamin
,
McCall
J.
,
Indriolo
N.
,
2019
,
ApJ
,
883
,
54

Park
S.
et al. ,
2005
,
ApJ
,
631
,
964

Ponti
G.
et al. ,
2015
,
MNRAS
,
453
,
172

Ponti
G.
et al. ,
2019
,
Nature
,
567
,
347

Possenti
A.
,
Cerutti
R.
,
Colpi
M.
,
Mereghetti
S.
,
2002
,
A&A
,
387
,
993

Raymond
J. C.
,
1992
,
ApJ
,
384
,
502

Revnivtsev
M.
,
Sazonov
S.
,
Churazov
E.
,
Forman
W.
,
Vikhlinin
A.
,
Sunyaev
R.
,
2009
,
Nature
,
458
,
1142

Serabyn
E.
,
Morris
M.
,
1994
,
ApJ
,
424
,
L91

Shore
S. N.
,
LaRosa
T. N.
,
1999
,
ApJ
,
521
,
587

Sicheneder
E.
,
Dexter
J.
,
2017
,
MNRAS
,
467
,
3642

Sofue
Y.
,
2013
,
PASJ
,
65
,
118

Sofue
Y.
,
2020
,
PASJ
,
72
,
L4

Sofue
Y.
,
Fujimoto
M.
,
1987
,
ApJ
,
319
,
L73

Sofue
Y.
,
Kigure
H.
,
Shibata
K.
,
2005
,
PASJ
,
57
,
L39

Sturrock
P. A.
,
Stern
R.
,
1980
,
ApJ
,
238
,
98

Su
M.
,
Slatyer
T. R.
,
Finkbeiner
D. P.
,
2010
,
ApJ
,
724
,
1044

Tanuma
S.
,
Yokoyama
T.
,
Kudoh
T.
,
Matsumoto
R.
,
Shibata
K.
,
Makishima
K.
,
1999
,
PASJ
,
51
,
161

Tanuma
S.
,
Yokoyama
T.
,
Kudoh
T.
,
Shibata
K.
,
2003
,
ApJ
,
582
,
215

Thomas
T.
,
Pfrommer
C.
,
Enßlin
T.
,
2020
,
ApJ
,
890
,
L18

Titov
V. S.
,
Démoulin
P.
,
1999
,
A&A
,
351
,
707

Toriumi
S.
,
Wang
H.
,
2019
,
Living Rev. Sol. Phys.
,
16
,
3

Tsap
Y.
,
Fedun
V.
,
Cheremnykh
O.
,
Stepanov
A.
,
Kryshtal
A.
,
Kopylova
Y.
,
2020
,
ApJ
,
901
,
99

Tsuboi
M.
,
Tsutsumi
T.
,
Kitamura
Y.
,
Miyawaki
R.
,
Miyazaki
A.
,
Miyoshi
M.
,
2020
,
PASJ
,
72
,
10

Uchida
K. I.
,
Morris
M. R.
,
Serabyn
E.
,
Bally
J.
,
1994
,
ApJ
,
421
,
505

Wang
Q. D.
,
2004
,
ApJ
,
612
,
159

Wang
Q. D.
,
Gotthelf
E. V.
,
Lang
C. C.
,
2002a
,
Nature
,
415
,
148

Wang
Q. D.
,
Lu
F.
,
Lang
C. C.
,
2002b
,
ApJ
,
581
,
1148

Wang
Q. D.
,
Lu
F. J.
,
Gotthelf
E. V.
,
2006
,
MNRAS
,
367
,
937

Wang
Q. D.
et al. ,
2013
,
Science
,
341
,
981

Weżgowiec
M.
,
Ehle
M.
,
Soida
M.
,
Dettmar
R. J.
,
Beck
R.
,
Urbanik
M.
,
2020
,
A&A
,
640
,
A109

Wiener
J.
,
Zweibel
E. G.
,
Oh
S. P.
,
2013
,
ApJ
,
767
,
87

Yamauchi
S.
,
Shimizu
M.
,
Nakashima
S.
,
Nobukawa
M.
,
Tsuru
T. G.
,
Koyama
K.
,
2014
,
PASJ
,
66
,
125

Yamauchi
S.
,
Shimizu
M.
,
Nobukawa
M.
,
Nobukawa
K. K.
,
Uchiyama
H.
,
Koyama
K.
,
2018
,
PASJ
,
70
,
82

Yusef-Zadeh
F.
,
Morris
M.
,
1987a
,
AJ
,
94
,
1178

Yusef-Zadeh
F.
,
Morris
M.
,
1987b
,
ApJ
,
322
,
721

Yusef-Zadeh
F.
,
Wardle
M.
,
2019
,
MNRAS
,
490
,
L1

Yusef-Zadeh
F.
,
Morris
M.
,
Chance
D.
,
1984
,
Nature
,
310
,
557

Yusef-Zadeh
F.
,
Wardle
M.
,
Parastaran
P.
,
1997
,
ApJ
,
475
,
L119

Zhang
S.
et al. ,
2020
,
ApJ
,
893
,
3

Zweibel
E. G.
,
Yamada
M.
,
2009
,
ARA&A
,
47
,
291

This article is published and distributed under the terms of the Oxford University Press, Standard Journals Publication Model (https://academic.oup.com/journals/pages/open_access/funder_policies/chorus/standard_publication_model)