Multi-Wavelength Observations and Modeling of Loop I †

The article aims to overview the origin and current dynamical state of a giant structure on the northern galactic sky—the radio Loop I extending from ≈−45◦ to ≈+45◦ in longitudes and up to ≈80◦ in latitudes over the Galactic center (GC). The main issue addressed here is a description of possible sources of mass and energy able to build up the Loop I and associated structures seen in X-ray, 21 cm, far infrared (FIR), and maintain them on long timescales. This region of the sky is highly crowded, such that contaminations from many projected structures can be tangled, and not always current direct observations look sufficient to disentangle them. At such conditions indirect arguments based on analysis of underlying star formation (SF) rate, morphological features in radio, X-ray and FIR may be important for understanding the origin of Loop I. Simple estimates show that the observed rather weak SF rate is able to create and maintain Loop I, and under certain circumstances can provide the observed east-west asymmetry. However, an explanation of an apparent coexistence of morphologically similar HI and FIR filaments close to Loop I is challenging, indicating that most likely they may belong to the foreground. Recently discovered absorptions in diffuse interstellar bands seem to confirm this picture.


Introduction
Brief history.The existence of a bright radio-spur originating near the Galactic plane (b 10 • ) and extending up to latitudes b ∼ +85 • -referred to as the North Polar Spur (NPS)-has became known since the pioneering observations at Jodrell Bank in late 1950s, and was clearly articulated first apparently by Hanbury Brown et al. [1].Haslam et al. [2] concluded that the NPS is the brightest north-east part of the circle-like structure with a diameter ∼100 • defined later as Loop I.They also resolved a narrow neck connecting the NPS to the plane at (l ≈ 30 • , b ≈ 13 • ) and several filamentary structures, what might hint on dynamical interaction of Loop I and the Galactic disk.Later, an intensive study of radio sky with higher resolution has revealed several loop-like structures defined as Loop II, III and IV [3][4][5].
Concepts.In spite of a 60 year long history of theoretical understanding of these radio loops and all related structures discovered in the area later, is still illusive.Historically two different points of view are established by the end of 1980s, and to a certain extent they still remain in competition: one suggests predominantly local origin of an expanding superbubble driven by stellar activity in nearby OB associations and local spiral arms within the solid angle ∼π pointing towards l ∼ 0 • , the other point of view is that the NPS and related structures originate near the Galactic center (GC).
Sco-Cen-centric model.The first point of view is based on a simple and naturally attractive assumption that connects the four round-shaped loops with large faint supernova remnants in a relatively nearby interstellar medium (ISM) [2,4], and/or with a joint envelope from stellar wind of Sco-Cen OB association.This point of view lies close to the concept of giant HI supershells ubiquitously produced in the ISM disk by OB-associations with total energy release up to 10 SNe [6,7], see also theoretical arguments in [8], and later development in [9].Sco-Cen is the nearest such region of a recent active star formation [10], and given its current position falls into the direction of the center of Loop I, it seems likely to be a very natural source of energy for maintaining this structure.
Line of sight spread model.A complementary model was drawn by Sofue [12] who associated the spur-like structures with synchrotron 'banks' elevated up to ∼1 kpc above the local (within ∼3 kpc) spiral arms.In this scenario, a slight (20 to 30 • ) inclination of the spurs in the anticenter direction is explained by the ram pressure of an outward blowing gas or magnetohydrodymic waves from the Galactic center [12,13].Later on such kind of structures have been described numerically in 2D MHD as hydraulic jumps combined with shocks from spiral density waves-the resulted plum-like vertical structures extending up to ∼0.5-1 kpc [14,15].
It is clear that the whole variety of structural features observed toward the Galactic center stems from crowded events along the line of sight.The problem is therefore to disentangle superposing contributions and identify the driving sources.
With rather uncertain distances inferred from the galactic rotation curve within |l| ≤ 30 • , a reliable complementary tool for distance measurements is 3D-reddening measurements [16], (see, also [17]).In such circumstances, a reasonable compromise for understanding a qualitative picture can be found in estimates of timescales, gas and dust mass and energy budget, synchrotron emission (magnetic field, spectral index), morphological features of the structures etc.The following issues are therefore critical for understanding the origin of Loop I:

1.
Energetics: whether the underlying star formation and supernova (SN) rates during the history of Loop I formation are sufficient to feed it?2.
Mass budget: whether the gas mass confined into Loop I (accounting in particular, X-ray gas in the NPS) is consistent with the overall evolutionary scenario?(see also [18]).

3.
Dust mass budget: whether the dust mass as inferred from infrared observations fits this evolutionary path, particularly when the morphology of dust distribution and its temperature with possible dust sputtering in high temperature environment of the NPS are accounted for?4.
HI and its kinematics: what would be the fraction of atomic hydrogen in the expanding shell associated with Loop I and whether it is consistent with observational data?Whether the velocity field in 21 cm toward the Galactic center adequately displays the expanding supershell associated with Loop I? (see also [19]).

5.
Synchrotron: are characteristics of synchrotron emission in Loop I (spectral index, magnetic field) are compatible with its evolutionary scenario?(see also [20]).

6.
Time scales: whether typical time scales: dynamical, radiative cooling, and dust destruction, do fit the evolutionary scenario?7.
Morphology and asymmetry: is Loop I (and the X-ray NPS) physically and genetically connected to the Fermi Bubbles?
Understanding these questions will provide generic qualitative constraints on evolutionary scenarios and physical mechanisms.In this review we will discuss energetics and mass budget of Loop I and their correspondence to the observed star formation rate in the underlying disk as the principal energy source.We focus on the morphological consistency of Loop I as seen in X-ray, synchrotron, HI and FIR emissions.Recently discovered diffuse interstellar bands within 2-3 kpc from Sun and their relation to Loop I are also briefly discussed.A systematic analysis of the origin and the distance to the NPS and its relation to the Fermi bubbles based on X-ray and gamma-ray data is given in the review by Kataoka et al. [18] in the same special issue of the Galaxies journal.

Local Hot Bubble and the Supershell from Sco-Cen
Within the first scenario with a bubble produced by explosions of a few SNe and stellar wind in the Sco-Cen OB association in the last 3-5 Myr the required energetics is of the order of Ė ≈ 7 × 10 37 erg s −1 , and the total energy released Ėdt ≈ 6 × 10 51 erg = 6 SNe [21,22].One of the key arguments in [22] scenario is the morphology of a dark ring-like structure in the ROSAT low energy band E = 1/4 keV (0.1-0.5 keV) which was attributed to absorption by thick HI ring formed as an overlap of the supershell from Sco-Cen OB association and the Local Hot Bubble (LHB) in which the Solar system is embedded.
The HI column density in the ring is N(HI) ≈ 7 × 10 20 cm −2 as measured by Copernicus, IUE and EUVE ultraviolet (UV) satellites (see in [22]).This number is at least half an order of magnitude higher than would be expected: with the assumed radius of the supershell from the Sco-Cen of R ss = 150 pc and the mean ambient density n 0 ∼ 1 cm −3 (see [23]), the swept up column density is N(HI) ∼ n 0 R ss /3 ∼ 2 × 10 20 cm −2 .Similar numbers for n 0 and the swept-up column density can be inferred from more recent measurements of HI and H 2 distribution in the inner galaxy (galactocentric radii R ≤ 8.5 kpc) by Marasco et al. [24], see their Figure 11 and Table 1.
Moreover, it has to be taken into account that the center of Loop I lies at around b ≈ 24 • , and thus the median vertical position of the driving SN explosions in the Sco-Cen OB should be z SN ≈ 70 pc, i.e., nearly 1.5 of the H 2 scale height (which is ≈50 pc, see top panel of Figure 10 in [24]).One would expect that at such circumstances the shock wave from the SNe will propagate predominantly upwards and hardly can break through the midplane of the gaseous disk to push sufficient amount of H 2 gas towards the LHB.It is also worth noting in this connection that the effective (Gaussian) scale height of neutral ISM gas-HI along with H 2 , is nearly flat h x ≈ 75 pc within R = 2-10 kpc as derived from observations of ≈100 Galactic point X-ray sources [25], thus making this inconsistency even more severe.Note also that if the X-ray absorbing shell toward the NPS was produced by collective SN explosions in ambient gas with n 0 ∼ 1 cm −3 , the net radiative energy losses would increase the required SN energy by an order of magnitude Ėdt ∼ 60 SNe.
In this context a very important conclusion has been reached in numerical simulations of the timing of 60 Fe isotop sediments in ocean by [26]: the mean ambient gas density (being assumed homogeneous) within the whole volume including LHB and Loop I may not exceed 0.3 cm −3 , otherwise the observed overabundant 60 Fe would not appear in the Earth's crust and ocean.If this the case, the discrepancy between the column density in the LHB-Loop I interface layer increases further: the estimated N(HI) ≈ 7 × 10 19 cm −2 is an order of magnitude lower than the observed towards the center of Loop I [22].In principle, the expansion law of SN remnants in a non-homogeneous (clumpy) medium with the mean density n = f v n 0 ( f v 1 is the volume filling factor of clumps, n 0 n the gas density in clumps) depends on the mean density approximately as n −1/5 Korolev et al. [27], such that the overall dynamical behavior remains close to what takes place in a homogeneous medium with density n = n .Moreover, measurements of polarization within the ring in [28] show a strong asymmetry between east and west parts and do not confirm it to be a unique ring-like structure.
An additional difficulty is related to the ages of stellar populations in the Sco-Cen association-the median age is ≈15 Myr with a large scatter between the subgroups: τ US = 13 ± 4 Myr in Upper Scorpius, τ UCL = 16 ± 5 Myr in Upper Centaurus Lupus, and τ LSS = 17 ± 5 Myr in Lower Centaurus Crux subgroups [29].Moreover, very recently Donaldson et al. [30] found a hierarcical decrease of ages in the US association on lower scales with τ sg 1 = 4.9 ± 0.8 Myr and τ sg 2 = 8.9 ± 0.9 Myr.At such circumstances a very rough estimate of the number of SNe what could have exploded 3-5 Myr ago is 3-10, which is apparently less than needed. 1  The HI ring-the interface between the Sco-Cen superbubble and the LHB, shadows the NPS emission in the ROSAT low energy E = 1/4 keV X-ray band.There is also shadowing in the higher energy ROSAT band (R5, 0.89 keV) seen as a sharp drop in X-ray counts at the latitudes b > 10 • , which can be attributed to absorption in the Aquila rift located much farther (at = 0.642 ± 0.174 kpc) [32].By assuming the NPS to be a part of a shell structure Sofue [32] put a lower limit of the distance to its center at 1.1 kpc, (see more detailed discussion in [18]).
Even stronger constraints have been inferred in [16] from an unprecedented analysis of absorption measurements of ROSAT R4, R5 and R6 (mean energies 0.725, 0.885, and 1.145 keV, correspondingly) and XMM-Newton data towards the terminus of NPS, accompanied with measurements of interstellar optical reddening within 3D Pan-STARRS survey and a fine spaced radio CO survey of the terminus region.One of the important consequences from scanning of the XMM-Newton surface brightness and hardness ratio, is that the NPS terminus extends lower toward the Galactic plane than might be expected if the absorption is attributed only to the local ring between the Sco-Cen superbubble and the LHB.Indeed, the regions contributing to X-ray hardening on lowest parts of the NPS terminus scanned in [16], likely lie beyond the Aquila rift, as suggested first by Sofue [32].The absolute lower limit on distance to the NPS boundary is set to 260 pc from comparing the HI absorption toward 3 stars in the NPS field and XMM absorbing column density.A similar lower boundary (300 pc) is derived from comparison of reddening to 19 dedicated stars in the XMM field and of dust maps for a larger sample of Pan-STARRS stars from [17].An additional interesting feature of the reddening measured by Pan-STARRS is its deficiency from 0.6 to 4 kpc, as clearly seen from Figure 11 in [16].Indeed, with the assumed scale height of 130 pc the inferred column density N(H) at low latitudes (b = 6 • -7 • ) should have increased between 0.6 and 4 kpc from 40 to 60 (×10 20 cm −2 ), while it remains nearly constant.This observation shifts the upper estimate on distance to NPS to the Planck limit ≈4 kpc derived from comparison of Pan-STARRS and Planck reddenings [33].
Radio polarimetry of the inner Galaxy (10 • < l < 34 • ) at 2.3 and 4.8 GHz placed the NPS terminus at a distance 2-3 kpc [34].Moreover, in [34] they found also the rotation of the position angle of magnetic lines at 2.3 and 4.8 GHz ∆φ ≈ +100 • towards a particular sight-line (l, b) = (22 • .7, 4 • .0)with the highest polarization.Sofue [32] derived from here the lower estimate of the distance towards the NPS terminus ∼5 kpc.
This line of arguments converges to the second scenario.The very first detection of a structure expected from wind-like outflow (the Galactic Center Lobe-GCL) of the very galactic center, has been reported by Sofue & Handa [11] on the base of 10 GHz observations.The east-side and west-side spurs restricting this "outflow" were detected earlier at lower frequencies (see references in [11]).The spectrum is derived to be thermal-with the slope α = −0.1, the structure therefore represents HII gas with electron density n e ∼ 5 cm −3 and temperature T ≈ 5000 K.The gas mass is M HII ≈ 4 × 10 5 M , and the required energy is E ≤ 10 51 erg [11].More precise estimate would increase the energy by at least 10 if radiation losses are taken into account (see Figure 1) provided the energy is injected by explosions. 1 Similar estimate is obtained recently from a detailed analysis of Sco-Cen star formation history [31].The volume filling factor of gas with given temperatures (see legend) for the same SN rates.From [35].
Later on Sofue [36] identifies in the 20-cm VLA image of the Galactic center arc a set of concentric shell structures 20 west-side and 5 souther of the GCL, and attributes them to a recurrent star bust events during the last few Myr.It seems quite natural under such circumstances to assume explicitly that star formation activity in the Galactic center not only drives the small-scale structures observed in the very narrow central part, but might have provided energy for the larger scale structure like Loop I, as mentioned first by Oort [37] and by Sofue & Handa [11], Sofue [38].Bland-Hawthorn & Cohen [39] combined 10 GHz data with mid-infrared (4.3-21 µm) image of a bipolar structure by the Midcourse Space Experiment (MSX), and extended them to much larger scales from (l to demonstrate their physical connection with the NPS (see their Figure 6).

Stellar Activity in the Central Molecular Zone and Loop I
The principal energy source for Loop I (and its hot X-ray NPS part) can apparently be connected only with stellar activity, and more explicitly with SN explosions; the luminosity of the supermassive black hole is too weak: ∼10 −9 of the Eddington luminosity, which is far from being sufficient. 2herefore the key parameter determining energetics of Loop I is the SF rate in the central region within |l| < 1 • -the central molecular zone (CMZ) [40].When averaged over the CMZ the current SF rate per unit area is Σ sf ≈ (0.09-0.11)M kpc −2 yr −1 (see Figure 16 in [41]), which gives overall SFR ≈ (0.009-0.017)M yr −1 .
It is worthwhile to note though that the optimistic conclusions by Sofue [36] and Bland-Hawthorn & Cohen [39] are drawn from a simplified hydrodynamical model of the central outflow, with the total energy injection from a few 10 55 ergs [39] to 3 × 10 56 ergs [38]. 3The SF rate in the central zone has a large uncertainty with estimates4 ranging from 0.04 to ∼ 1.5 M yr −1 (see [39]).A more detailed state-of-the art simulation by Sarkar et al. [42] has confirmed such a conclusion: in order for the central activity to produce a large scale Ω-shaped wind-driven supershell with radial and vertical sizes comparable to the observed in the NPS, the central SF rate has to be Ṁ ≈ 0.5 M yr −1 , with the estimated age of 24 Myr.This is at least factor 3 in excess of the observed SF rate: Spitzer and MSX data lead Yusef-Zadeh et al. [41] to conclude that the peak SF rate in the central 200 pc × 50 pc region traced by 10 5 years young stellar objects (YSO) is 0.14 M yr −1 .This estimate is likely an upper limit because not in all cases the sources from Yusef-Zadeh et al. [41] are indeed YSO.For instance, ref. [43] suggests that main sequence stars in the CMZ might "masquerade" as YSOs, which can cause a threefold overestimate of the SFR in this area.On a longer timescale an estimate 0.007 M yr −1 is inferred from the nonthermal radio emission.Overall, Yusef-Zadeh et al. [41] assume for a recent SFR 0.04-0.08M yr −1 , while over the Hubble time the following strict upper limit is put 0.15 M yr −1 .
More recent estimates of the SFR ≈ 0.09 M yr −1 within the last 0.1-3 Myr were derived from Spitzer and Herschel data by fitting fluxes at 24 µm and 70 µm, and the bolometric spectrum in the range 5.8 to 500 µm [44].
Therefore, following these estimates and numerical models the observed SF rate in the Galactic center seems to be weaker than needed to produce large scale Loop I structure within the Galacto-centric scenario, unless the initial mass function of stellar population in the CMZ is a heavy-top type (see discussion in [45]).

Continuous Wind vs. Instant Implosions
A possible solution of this inconsistency can be sought in changing the energy injection regime in the GC.An important detail of the numerical model by Sarkar et al. [42] is a continuous injection of energy and mass in a luminosity-driven wind regime.In a uniform (non-stratified) medium the wind regime is energetically the most efficient.The wind regime suggests frequent ("coherent") SN explosions [46,47], in which case a smaller fraction of energy becomes lost, as obviously seen in Figure 1 [35]: at higher SN rates-ν sn ≥ 10 −9 yr −1 pc −3 , the fraction of total explosion energy retained at high temperatures approaches asymptotically E/E t 0.3.In a spherically symmetric model by Sharma et al. [48] the transition to the wind regime occurs when energy is injected by individual SN events with a rate higher than a certain critical value, such that thermolization occurs before a considerable fraction of energy is lost radiatively, see Figure 2. The differences in numerical values of SN rates and the fractions of lost energy between [35] and [48] are because in [35] SNe explode randomly in the volume 200 3 pc 3 , while in [48] all energy in injected into the central numerical sphere with radius of 1 to 2 pc.Density profiles vs normalized radius of the remnant formed by sequential explosions of SNe at r = 0 in a small injection sphere with r ej , see details in [48].Energy is injected by N OB with a rate t −1 sn = N OB /t OB , t OB = 30 Myr as assumed in [48]; energy is injected in the form of kinetic energy (KE), thermal energy (TE), or the bubble is fed by the mechanical luminosity It is seen that for rare explosion events each SN ejecta expands into a very low-density cavity-it is nearly an expansion into vacuum.For higher N OB the profile approach the Chevalier-Clegg wind solution [49], where every portion of energy injection has to share it with the whole bubble through thermolization.From [48].
Dynamics change in a stratified medium.In this case, after a few SNe explosions, the collective bubble breaks out and continues growing predominantly upward.As a result, all subsequent explosions expand nearly into vacuum with minimal energy losses.A brief episode of a considerable radiative energy loss may occur in the initial stages, when the shock front propagating downward partly enters denser layers of the ISM disk where the radiative loss rate increases.It is seen in Figure 3: left panel depicts thermal energy, in the right panel the fraction of thermal energy is shown for a growing superbubble from clustered SNe explosions in the Galactic center.SNe explode one by one spread randomly in a volume of stellar cluster with radius r c = 30 pc; the rates are 10 −3 yr −1 (upper curves in the left panel), and 10 −4 yr −1 (lower curves).The star formation rates are ≈0.08 M yr −1 and 0.008 M yr −1 correspondingly, with a Kroupa IMF.On Figure 4 below the bubbles are shown for these two models: left panel is for 10 −4 SNe per yr at t = 4 Myr, while right panel is for 10 −3 SNe per yr at t = 2 Myr [50].
With a constant SF (and SN) rate the top parts of superbubbles grow approximately as z ∝ t 0.5 and will reach the height of 10 kpc within 50 Myr for the lower rate, and within 25 Myr for the higher one (see in [50]).One may conclude therefore that even such an anemic SF rate 0.04-0.08M yr −1 as measured in [41] can support large scale structures like Loop I and X-ray NPS when a regime of discrete energy injection from SN explosions is accounted for [50].In an alternative model with two sequential star burts separated by around 10-15 Myr are responsible for the observed Loop I and NPS [18].

Mass Budget
Transporting mass from the ISM disk into halo is an issue deserving a separate discussion.Generally it is thought that fountaining SN remnants efficiently elevate interstellar gas.However, as shown by Mac Low et al. [51] only at most 3% of the mass from an exponential or a gaussian interstellar disk is indeed elevated, while the rest remains in wall-like shells with a vertical extent of a few scale heights.This kind of a shell is observed in the terminus of the NPS in 21 cm-a "crater" with a height 450 pc [52].On smaller scales there is a 200-pc radius molecular cylinder reported by Sofue [53], with a vertical FHM length ∼60 pc-nearly twice of the CMZ scale height, might be due to a more recent explosion event in the Galactic center.
Therefore, when the mass budget of extraplanar structures, like NPS, is concerned the two sources are relevant: mass-loading by SNe explosions in the Galactic center or its disk, and partly the mass of a swept-up halo gas with a shallow density profile.It is exactly what happens in 2D numerical models by Sarkar et al. [42] and Sofue et al. [54]: the total mass injected by SNe is M inj ≈ 4 × 10 6 M and the swept-up halo gas in 2π shell is M sw ∼ 5 × 10 6 M in [42].while M inj ≈ 10 7 M and nearly the same for the swept-up mass in [54].It fully corresponds to the assumed SFR: 0.5M yr −1 in [42], and ∼3M yr −1 in [54], with the vertical density profiles with n d = 3 cm −3 and n h = 2.2 × 10 −3 cm −3 for the disk and halo gas in [42], while n d = 2 cm −3 (for gaussian z-profile) and n h = 2 × 10 −2 cm −3 in [54].It is also valid in the 3D model of a growing superbubble with low SF rate shown in Figure 4: gas mass in the shells (≈ 3 × 10 5 M on the left, and ≈ 4 × 10 5 M on the right) consists mostly on the swept-up halo gas, while SNe ejecta contributes little.This is consistent with the conclusion by Kataoka et al. [18] that the X-ray emitting gas in the NPS is a halo gas compressed by shocks coming from the GC.

M X -The Hot Gas Mass in the NPS
Mass of X-ray gas is with Ω nps 0.3 being the solid angle of the NPS (of ∆l ∼ 20 • in width, and ∆b ∼ 80 • in length), D 10 , the distance to the base normalized by 10 kpc, H 1 is the thickness of the emitting layer in units of 1 pc.
H and n e are connected through the emission measure EM ∼ 0.1 pc cm −6 [55].Depending on the NPS morphology H can vary from around the projected thickness of the NPS to the chord along the line of sight perpendicular to the projected NPS.It can be estimated as αD, where α varies from ≈0.25 in the first case at b 30 • , to around 0.6 for the second case.With a given emission measure the mass of X-ray gas in the NPS is then For α ≈ 0.25-0.6 and D 10 = 0.8 this estimate is in accord with the above numbers shown for numerical models by Sarkar et al. [42] and Sofue et al. [54] when only the swept-up mass is taken into account.The corresponding thermal energy E X ∼ 10 55 erg would be injected by SN explosions corresponding to star formation with the rate ∼0.04 M yr −1 in ∼30η −1 X Myr, where η X > ∼ 0.3 is the heating efficiency of collective SN explosions.

NPS Morphology
It is obvious that the observed asymmetry of the X-ray NPS and 408 MHz Loop I, along with all features seen in HI maps, dust emission, Planck dust optical depth-their concentration mostly towards the eastern part, is an essential and important characteristic from the point of view of their origin.In the galacto-centric scenario the structure which manifests through these features, requires huge energy-E ∼ 10 54 -10 56 erg and mass-M ∼ 10 7 M , which suggests averaging over large scales in the central Galactic region, and therefore a rather symmetric configuration. 5For mean surface density in the Galactic center Σ ∼ 10 3 M pc −2 [58], one can get the gas mass around 3 × 10 7 M within the radius 100 pc.To explain the north-east asymmetry of the NPS within the Galacto-centric model, 5 All numerical models of galacto-centric winds deal with axially symmetric density distributions, besides small-scale deviations from the symmetry, like in models of galactic wind filamentation [56,57].
one has to assume a highly energetic flow which would elevate selectively a third of gas mass of the Galactic central region up to heights z ∼ 10 kpc north-east keeping its integrity.It looks very unlikely provided the central SN activity is rather symmetric, because the flow prefers to extend towards low density regions with local velocity v ∝ ρ −1/2 .In addition, Rayleigh-Taylor and Kelvin-Helmholtz instabilities act on timescales shorter than that needed for the flow to carry this mass up to ∼10 kpc.
One of possible solutions here might be connected with highly nonhomogeneous distribution of star formation and SN rate in the central molecular zone (see review in [45]).In this case an occasionally east-biased local starburst in the Centrale molecular zone can initiate an outflow mostly towards the north-east side of the GC (see examples in [50]).

HI Around NPS
In spite of a general morphological similarity of all structures around Galactic symmetry axis: NPS, synchrotron (e.g., 408 MHz) emission, HI, Planck dust optical depth at 353 GHz τ 353 , several features deserve a closer look.
Soft X-ray as seen from ROSAT, Suzaku and XMM-Newton reveals high-temperature plasma: T 0.3 keV (see in [59], for more recent discussion see [18]), which suggests that the gas is highly ionized by a high-velocity shock wave(s) from SNe explosions in the center.In this scenario it looks unlikely that the observed HI loop-like gas is directly connected to the NPS and Loop I. Figure 5 shows a 3D numerical model of HI column density map of the shell being swept-up by central SNe explosions with the rate 10 −4 yr −1 [50].HI shell can be recognized only up to heights z < ∼ 0.6 kpc, while at higher z strong shock waves completely converge atomic hydrogen into fully ionized-only a few tiny clumps and filaments survive up to z ∼ 1 kpc.Please note that the postshock gas does not cool radiatively at z > 0.5 kpc because the cooling time is longer than the dynamical one (see below).This circumstance excludes possible formation of clouds via thermal instability.
It is important though to stress here that this model implies a homogeneous vertically stratified density.However, the existence of clumpy and molecular outflows (e.g., in M82) suggests a nonhomogeneous ambient gas distribution before the starburst (see, e.g., filamentation model by Tanner et al. [57]).Though at higher heights of z > 0.5 kpc, clumps and filaments are usually getting erased (see in [56,57,60]).From this point of view the presence of HI high-velocity clumps in a conical region towards the Galactic center at z > ∼ 1.5 kpc (b 10 • ) demonstrated by McCluer-Griffits [19] looks exceptionally interesting.For conditions in the NPS, gas remains highly ionized deeper behind the shock front because its radiative cooling time t c ∼ k B T/Λ(T)n e > ∼ 100 Myr for n e < ∼ 0.01 cm −3 (see in [50]).This might mean that the HI low latitude tongue around the polar axis is most likely a front-side feature not connected with the shock wave presumably responsible for creating the NPS.Another feature is the low column density tongue (log 10 N HI 20.5, and the corresponding mass of M HI ∼ 5 × 10 7 M ) offset north-east of the NPS, and extending upwards closer to the polar axis.It continues farther north-west beyond the polar axis along Loop I (see Figure 6).The overall morphological picture of NPS, radio-continuum Loop I and HI loop looks like if they are physically connected.However, if HI north-east filament is connected with the X-ray NPS, it precedes the shock front and thus seems to stay not in the right place.One issue related to this HI feature is though worth to speculate: isolated HI structures in the ISM being exposed to the diffuse Lyman continuum, are expected to be bounded by HII envelopes, which might be seen in Hα emission.A rough estimate of the Hα surface brightness in the rims is [50], assuming the structure of a cylindrical filament and exposure by diffuse Lyman continuum flux F Ly−c ∼ 2 × 10 6 phot cm −2 s −1 [64,65].Therefore, when such a filamentary structure lies at D ∼ 10 kpc the rim brightness falls below the detection limit.However, if they have local origin and are at D ∼ 100 pc the ionized rims around them might be recognizable.Such rims are not seen though at high latitudes in Hα all-sky surveys and are not obviuosly present around the HI loop close to the NPS [66][67][68], indicating probably that this tongue is located at high heights closer to the polar Galactic axis.Thus the origin of this HI loop is not quite clear.
Abundance Pattern in the NPS Miller et al. [59] reported that heavy elements (C, O, Ne, Mg and Fe) are significantly depleted in the NPS down to less than 0.5 solar, while N is highly enhanced: N/O 4 as compared to the solar N/O ratio.This observation may indicate a contribution from an elder population of asymptotic giant branch stars.Miller et al. [59] conclude that such an enhancement can be understood in the Sco-Cen scenario as due to a swept-up ISM gas enriched on longer time by Pop I stellar population.In the Galacto-centric scenario such an explanation of overabundant N can be problematic: with the total X-ray mass in the NPS of M X ∼ 0.4 × 10 7 M (see, Section 2.2.1) the mass of nitrogen in it would be of M X,N ∼ 6 × 10 3 M while the expected mass of oxygen from the GC is M X,O > ∼ 10 4 M , assuming M O > ∼ 1M oxygen from a single core-collapse SN [69].Another important issue concerns the presence of dust in an immediate neighbourhood of the X-ray NPS gas (see lower panels in Figure 6).The optical depth τ 353 at the lower left panel of Figure 6 shows a clearly seen filament following the HI loop from north-east to north-west.Its eastern part lies ≈10 • easter of the X-ray NPS (before the front, if the NPS traces hot gas behind the shock moving eastward).The total gas mass in this tongue estimated from τ 353 ∼ (0.6-1) × 10 −5 , after integration over the solid angle Ω ∼ 0.1 for this dust structure 6 is: M HI ∼ 6 × 10 7 D 2 10 M , which is consistent with the above estimated HI mass in the slightly north-east offset 21 cm filament coincident with the τ 353 one.
Given that the dust-to-gas mass ratio is similar to the one in the Galactic disc ζ d = 140 −1 , the total amount of dust in this feature can be as M d ∼ 4 × 10 5 D 2 10 M .One can conclude from this that the dominant source of dust mass in τ 353 filament can be only the dust in swept-up gas, which therefore is supposed to have the dust mass fraction not less than in the ISM disk.
The two fllaments-HI and τ 353 represent apparently a separate physically connected gas structure standing aside the X-ray gas of NPS.Moreover, as mentioned above, their position with respect to X-ray gas does not look consistent if the latter is driven by a shock wave from the Galactic center, and is hard to be interpreted in the Galacto-centric scenario.

Dust Temperature
The dust temperature distribution around the X-ray component of NPS (see the lower right panel in Figure 6) displays a hot (T d ≈ 23 K) dust tongue nearly coincident with the NPS.If the tongue is physically connected to the NPS, then dust temperature T d is determined by collisional heating from hot plasma-the equilibrium value is around T d ∼ 18n 1/6  0.01 T 1/4 0.3 a −1/6 0.1 K (n = 0.01n 0.01 cm −3 -plasma density in NPS, T = 0.3T 0.3 keV, its temperature, a = 0.1a 0.1 µm-dust grain radius): particles with a = 0.01-0.1µmfit well the whole range of dust temperature in the NPS.
However, when dust grains are maintained by thermal collisions at T d ∼ 20 K, their characteristic survival time against thermal sputtering is only few to ten Myr for silicates and carbonaceous grains, respectively, which is obviously shorter than the dynamical time for the NPS.From this point of view this dust may be a residual fraction of dust being destroyed in the hot hostile environment of X-ray plasma.However, with the mean flux density F ν ∼ 1 MJy at 353 GHz from the region occupied by the T d ∼ 20 K dust (solid angle Ω ∼ 0.05) within the NPS X-ray hot gas, one can estimate the corresponding dust mass M d ∼ 3 × 10 4 D 2 10 M .The dust mass fraction within X-ray NPS gas is for the assumed numbers.This ratio seems to be at least factor of 3 higher than would be expected for a 30 Myr old NPS X-ray region, suggesting that two additional parts of hot (T d ≈ 23 K) dust lie in the same projected area as the NPS. 6 For this estimate we used dust extinction model σ ext (λ > 30 µm) ≈ 10 −26 cm 2 per H atom at λ ∼ 0.1 mm, as accepted for the diffuse ISM at high galactic latitudes in [70].

Synchrotron
Supernova remnants and their ensembles growing within the ISM disk (i.e., within 1 to 2 scale heights) have similar physical characteristics at their shells: gas density, temperature, cosmic ray density.For such remnants the Σ − D relation falls on to the general trend (see, discussion in [71][72][73]), and the spectral index of their synchrotron emission is expected to be generally invariant along the shell (see [72]).Remnants of larger scales produced by a collection of SNe and consequently extending upward higher than 2 to 3 scale heights are expected to become individualized by their large scale gradient of the spectral index α.Such a gradient emerges because of several reasons: (i) relativistic electrons lose their energy ∝E 2 , such that the characteristic cooling time for them t cool ∝ E −1 ; and (ii) this effect is enhanced for relativistic electons being elevated in a vertically growing superbubble, in which an additional adiabatic cooling due to diverging velocity field in the superbubble shift their energy spectrum towards lower energies (see in [74]).One can expect therefore the Loop I spectral index to grow upwards, i.e., the spectrum to get steeper along the Loop I from its terminus at low latitudes to its polar parts.
Observations of polarized emission in the east northern part of Loop I [73] do not confirm this conclusion: going from the northest (field 1, Figure 6 in [73]) through the middle (field 2) and to the lowest (field 3) part α changes from 0.26 to 0.12 and to 1.3, respectively (see upper panel on Figure 7 in [73]).At the same time, global maps of synchrotron intensity from 10 MHz to 2.3 GHz obviously show softening in the spectrum from the central midplane area towards the north-north-east direction [75].Similar behavior is found in the wind flow of M82 with strong variations of α along the outflow: from 0.6 in the central region to 1.3 on scales of 0.5 kpc outward [76].Sychrotron spectral index in halos of edge-on galaxies do also reveal a considerable increase in the vertical direction consistent with a strong (velocities 100 to 250 km s −1 ) advective transport [77].

DIBs in the NPS
Lan et al. [78] scanned around 5000 deg 2 of the northern Galactic hemisphere mostly at high latitudes (l > ∼ 30 • ) including the NPS domain, in order to find diffuse interstellar bands (DIB) in the Milky Way halo.The most interesting result of their study connected to the question of interest-the origin of Loop I and related structures, is that the strongest DIB absorptions (W ∼ 0.5-0.8Å) are located exactly in the HI filaments close to the NPS X-ray part from east and west sides, and demonstrate a very similar morphology seen in Figure 7, left panel.This observation is reflected in a very tight correlation between their DIBs absorptions and the total-integrated over the integrated velocity range from −450 km s −1 to +400 km s −1 of atomic hydrogen in the LAN survey [63]: W DIB ∝ N α HI N µ H 2 with W DIB being the equivalent width of a DIB feature, α ∼ 1, |µ| 1 (see Figures 11 and 12, and Equation (4) in [78]).
A striking feature is that the typical distance to SDSS stars in this coverage is 2 to 3 kpc as seen from Figure 7, right panel.This finding definitely demonstrates that high-latitude HI tongues, and possibly connected to them gas traced by τ 353 , cannot be directly associated with the NPS X-ray gas provided it locates closer to the Galactic center.Otherwise they all: HI, τ 353 and X-ray, are placed nearby at D < ∼ 2-3 kpc, and represent a part of a huge (quasi)conical shell structure in the foreground.

Conclusions
Current observations and models of Loop I and related structures can be summarized as follows: • Energetic events in the Galactic center and the very existence of the Fermi-bubbles make it very attractive to think, that all this is produced by a collective SN activity in the center.The conclusion of whether a weak central SF rate-SRF < ∼ 0.04 M yr −1 is sufficient for maintaining these structures, crucially depends on the numerical model: luminosity-driven wind scenarios require an order of magnitude higher SFR, while randomly spread isolated SNe can apparently do the work with a lower SFR.

•
The mass of HI filament inferred from 21 cm column density, and the gas mass of a coincidental τ 353 filament north-east of the X-ray component of NPS are nearly equal, which would mean that they are physically connected.However, their position with respect to the X-ray NPS indicates most likely that they represent warm gas isolated of the X-ray gas.

•
Mass budget and energetics of the X-ray gas of NPS can be maintained by the central star formation of SFR ∼ 0.04 M yr −1 during ∼30η −1 X Myr, within the model of randomly positioned isolated SN explosions as the energy source.However, the dust mass in the X-ray NPS as inferred from far infrared emission at 353 GHz is factor of 3 higher than would be expected from a more than 30 Myr old NPS.

•
An enhanced N/O ratio in the NPS reported by Miller et al. [59] does not look to match conditions in the Galacto-centric model.It is easer understood within the Sco-Cen scenario.

•
Loop I synchrotron emission and, in particular, softening of its spectrum toward lower frequencies (10 to 820 MHz) in comparison to higher ones (1.4 to 2.3 GHz) is consistent most likely with the galacto-centric scenario.

•
On the other hand, discovery of DIB SDSS absorptions in spectra of stars within 2 to 3 kpc sheds a new light on the problem, because the most strong DIB absorptions seem to outline the same regions on sky containing the Loop I and related HI, τ 353 and the X-ray structures.This may indicate that either these structures are all placed relatively locally, within 2-3 kpc, or otherwise one has to conclude that not all of them are physically connected.For instance, synchrotron Loop I and X-ray NPS gas might arise due to SN activity in the GC, while HI and τ 353 filaments and plumes along with DIB absorptions represent frontside structures within 2-3 kpc.
The overall conclusion is therefore that even though some components of Loop I seem to be connected with activity of the Galactic center, including relativistic electrons shining in synchrotron and perhaps the hot NPS X-ray gas, there are contaminations from structures lying closer to us, and projecting onto structures over the GC.New techniques, such as 3D measurements of interstellar dust redenning or 21 cm velocity field towards Loop I, are to be applied in order to disentangle projection effects.

Figure 1 .
Figure 1.(Left) Evolution of the heating efficiency-a fraction of total energy injected by explosions retained in form of thermal energy in the phase with a given temperature (as indicated in the legend), top to bottom: SN rate 10 −3 , 10 −4 , 5 × 10 −5 , ∼3 × 10 −5 , 2.5 × 10 −5 and 10 −5 yr −1 in the computational volume 200 3 pc 3 -the volume nearly equal to the CMZ.(Right)The volume filling factor of gas with given temperatures (see legend) for the same SN rates.From[35].

Figure 2 .
Figure2.Density profiles vs normalized radius of the remnant formed by sequential explosions of SNe at r = 0 in a small injection sphere with r ej , see details in[48].Energy is injected by N OB with a rate t −1 sn = N OB /t OB , t OB = 30 Myr as assumed in[48]; energy is injected in the form of kinetic energy (KE), thermal energy (TE), or the bubble is fed by the mechanical luminosity L ej = N OB E 0 /t sn , E 0 = 1 B. It is seen that for rare explosion events each SN ejecta expands into a very low-density cavity-it is nearly an expansion into vacuum.For higher N OB the profile approach the Chevalier-Clegg wind solution[49], where every portion of energy injection has to share it with the whole bubble through thermolization.From[48].

Figure 3 .
Figure 3. (Left)Thermal energy in a superbubble from sequential SNe explosions in a stratified ISM in the Galactic center; black lines depict the total thermal energy, red lines-thermal energy of gas with T > 10 6 K, upper curves are for the SN rate 10 −3 yr −1 , lower are for 10 −4 yr −1 , SNe explode in a cluster of radius r c = 30 pc, gas scale height is 100 pc, the halo gas component has constant density n h = 10 −4 cm −3 .(Right) The ratio of thermal to total injected energy for all three model in the (Left) panel.The flattening of the energy growth is clearly seen in the (Left) panel between 0.3 and 0.8 Myr for lower SN rate, and between 0.2 and 0.6 Myr for higher SN rate.After the flattening radiative losses in the disk become inefficient compared to continuously injected energy, and thermal energy grows almost linearly.Later on, at around 1 to 4 Myr radiation increases in the shell when the downward parts of the shock wave cover a considerable fraction of the ISM disk.From[50].

Figure 5 .
Figure 5. HI distribution in an expanding superbubble in the model shown in Figure 4 at the age 4 Myr (the central region is zoomed in to compare to Figure 4), color bar shows HI column density from log 10 N HI = 15 to 18 cm −2 .(From [50]).

Figure 7 .
Figure 7. (Left) The full sky map of cumulative DIB absorptions (combined in total of 20 DIB lines) smoothed with a 5 • -gaussian kernel; the map includes regions with E(B − V) > 0.02; (Right) the map showing median distance to the SDSS stars used in the coverage[78].Clearly seen is that in the area of HI tongues neighbouring the X-ray NPS, the median distance is fairly flat of 2.5 kpc with a weak patchy deviations.Both maps are from[78].