Letters

PRESSURE EQUILIBRIUM BETWEEN THE LOCAL INTERSTELLAR CLOUDS AND THE LOCAL HOT BUBBLE

, , , , , , , , , , , , , , , , , and

Published 2014 July 28 © 2014. The American Astronomical Society. All rights reserved.
, , Citation S. L. Snowden et al 2014 ApJL 791 L14 DOI 10.1088/2041-8205/791/1/L14

2041-8205/791/1/L14

ABSTRACT

Three recent results related to the heliosphere and the local interstellar medium (ISM) have provided an improved insight into the distribution and conditions of material in the solar neighborhood. These are the measurement of the magnetic field outside of the heliosphere by Voyager 1, the improved mapping of the three-dimensional structure of neutral material surrounding the Local Cavity using extensive ISM absorption line and reddening data, and a sounding rocket flight which observed the heliospheric helium focusing cone in X-rays and provided a robust estimate of the contribution of solar wind charge exchange emission to the ROSAT All-Sky Survey 1/4 keV band data. Combining these disparate results, we show that the thermal pressure of the plasma in the Local Hot Bubble (LHB) is P/k = 10, 700 cm−3 K. If the LHB is relatively free of a global magnetic field, it can easily be in pressure (thermal plus magnetic field) equilibrium with the local interstellar clouds, eliminating a long-standing discrepancy in models of the local ISM.

Export citation and abstract BibTeX RIS

1. INTRODUCTION

The observed surface brightness of the soft X-ray diffuse background (SXRB) at 1/4 keV is closely linked to the local structure of the interstellar medium (ISM). There are shadowing features demonstrating the clear existence of foreground and background emission regions relative to many ISM clouds within a few hundred parsecs of the Sun (e.g., Snowden et al. 2000). The existence of a non-zero X-ray flux in the Galactic plane, where nearby H i column densities of a few times 1020 cm−2 are sufficient to be optically thick, further requires a local origin for that flux. Combining the knowledge of a local deficiency of neutral hydrogen (e.g., Knapp 1975), and the use of reddening and optical absorption line measurements to map the geometry of that deficiency (e.g., Frisch & York 1983; Paresce 1984), led to the identification of the Local Cavity (LC) in the H i of the Galactic disk near the Sun. The requirement for a local source of X-ray emission coupled with the existence of the LC led to the concept of the Local Hot Bubble (LHB), where the LC contains an X-ray emitting thermal plasma with a temperature of T  ∼  106 K. The LHB is confined within the LC but does not completely fill it, extends from ∼50–200 pc from the Sun, and is preferentially elongated toward high Galactic latitudes (see Sanders et al. 1977; Cox & Snowden 1986; Snowden et al. 1990). Reviews of the 1/4 keV SXRB can be found in McCammon & Sanders (1990) and Snowden (2002).

One well-known problem with the existence of the LHB was that the required thermal pressure of the plasma was much larger than the thermal pressure of the local interstellar clouds (LICs) surrounding the Sun (e.g., Jenkins 2009). While magnetic field pressure within the LICs could account for some of the difference, the required field strength for equilibrium was unreasonably large.

An additional complication to the LHB picture came with the identification of solar wind charge exchange (SWCX) emission as a possible source for a significant flux of diffuse X-rays at 1/4 keV. SWCX X-ray emission occurs when highly ionized metals in the solar wind interact with neutral material gaining an electron in a highly excited state, which then decays emitting a photon (or photons) with a characteristic energy of the ion. First identified in ROSAT observations of comets (Lisse et al. 1996; Cravens 1997), Cox (1998) and Freyberg (1998) suggested that SWCX could be the source for the non-cosmic background enhancements observed in the ROSAT All-Sky Survey (RASS; Snowden et al. 1995). Cravens et al. (2001) confirmed this by showing the strong correlation between solar wind flux and the background enhancement count rate.

The question was then: how much of the 1/4 keV SXRB originates in the LHB and how much originates as SWCX with interstellar neutrals flowing through the solar system (e.g., Lallement 2004; Welsh & Shelton 2009) or with exospheric material in Earth's magnetosheath (e.g., Robertson & Cravens 2003)? While the temporally varying SWCX component, with timescales from a few hours to a few days, was modeled and removed from the RASS, any longer-term variations (such as those arising as distributed emission through the heliosphere where the variation is smoothed out) could not, and so remain as an uncertain zero-level offset in the RASS maps. It was that unknown offset that might possibly have provided all of the observed local emission and thus obviated the need for the LHB.

2. NEW TOOLS FOR THE LICs AND LHB

Three recent advances have shed considerable light on the relationship between the LICs and the LHB: (1) the Diffuse X-rays from the Local galaxy (DXL) sounding rocket experiment (Galeazzi et al. 2011) produced a measurement of the amount of SWCX emission from the heliosphere that contributes to the RASS 1/4 keV map (Galeazzi et al. 2014, hereafter G14); (2) the three-dimensional mapping of neutral material in the local ISM using absorption line and reddening data has greatly improved our knowledge of its distribution in the nearest ∼500 pc from the Sun (Lallement et al. 2014; Puspitarini et al. 2014); and (3) the penetration of the Voyager 1 spacecraft beyond the heliosphere has allowed a more reliable in situ measurement of the interstellar magnetic field within the LICs (Burlaga & Ness 2014).

2.1. X-Ray Data

G14 provides an estimate for the heliospheric and magnetospheric SWCX contribution to the RASS (Snowden et al. 1997) 1/4 keV band (R12 band) by fitting the observed count rate along the DXL scan path, the RASS count rate for the scan path, and the model SWCX emission from the heliosphere calculated for the sounding rocket and RASS observation geometries and solar wind conditions. The fractional contribution of SWCX emission to the RASS data along the DXL scan path is 0.40 ± 0.07 (combining systematic and statistical errors in quadrature), typically 140 RU (ROSAT units, 10−6 counts s−1 arcmin−2). Using the HEASARC X-ray Background Tool,10 a RASS spectrum was extracted for a 4fdg5 radius region centered on l  ∼  144°, b  ∼  0°, the direction of minimum surface brightness along the DXL scan path. The R12 count rate for the extraction is 339 ± 2 RU. (Note that 0.40 ± 0.07 of the extracted count rate is 136 ± 24 RU, consistent with the typical SWCX flux value of 140 RU stated in G14.) The spectrum was fit (using Xspec11) with a simple three-component model for the diffuse X-ray background consisting of (1) an unabsorbed thermal emission component representing the sum of emission from SWCX and the LHB, (2) an absorbed thermal component representing emission from more distant regions of the Milky Way, and (3) an absorbed power-law component representing the extragalactic background. The fitted parameters for the unabsorbed component were a temperature of kT = 0.102 ± 0.001 keV (T = (1.18 ± 0.01) × 106 K), and an Xspec normalization of $5.37^{+0.05}_{-0.10}\times 10^{-7}$ for an emission measure $n_e n_{\rm H} L = 2.59^{+0.02}_{-0.05}\times 10^{-3}$ cm−6 pc, where L is the path length through the emission in parsecs and ne and nH are the electron and proton densities in the plasma, respectively (assumed to be constant along the path length). Current APEC (Smith et al. 2001) thermal emission models with Anders & Grevesse (1989) abundances were used for the fit. The fitted absorption column density was ∼1.3 × 1021 cm−2 so the two absorbed spectral components have no contribution to the R12 band. While there is some dependence on the choice of thermal emission model, the APEC model was the most conservative as it led to the highest normalization and thermal pressure. There is also a dependence on assumed elemental abundances, however, using the Savage & Sembach (1996) warm cloud abundances led to a thermal pressure only ∼15% higher and a required LIC magnetic field pressure for equilibrium only ∼10% higher than that derived below.

Subtracting the SWCX contribution from the RASS yields 203 ± 24 RU for the LHB emission. Scaling the fitted Xspec emission measure by 203/339 to represent just the LHB emission, and assuming that the LHB and SWCX spectra are similar (demonstrated during the processing of the RASS data; Snowden et al. 1995), yields an inferred emission measure for the LHB component of nenpL = (1.55  ±  0.18) × 10−3 cm−6 pc.

2.2. Interstellar Thermal Pressure

From Lallement et al. (2014), the distance to the nearest wall of H i in the direction of l  ∼  144°, b  ∼  0° is ∼90 ± 15 pc (see Figure 1). However, part of that path length is filled by the LICs so we subtract an estimate of 5 pc (Frisch et al. 2011) leaving a plasma path length of 85 ± 15 pc (assuming that the error in the path length in the LICs is much smaller than the uncertainty in the distance to the edge of the LC). Again, assuming that ne and nH are constant along the line of sight and using Xspec's assumption that $n_e n_{\rm H} = n_e^2/1.2$, we have

Equation (1)

The thermal pressure of the plasma is then

Equation (2)

where n is the total particle density (n  ∼  1.92ne). This is ∼60% the LHB thermal pressure estimated by Galeazzi et al. (2007), which well represents previous measurements that did not include the effect of SWCX emission. (Note that the evolution of thermal models can affect the pressure scaling, as well as the accuracy of assumed path length through the emitting plasma.) On the other hand, the value of PT/k = 10, 200 cm−3 K derived in Puspitarini et al. (2014) is in quite good agreement. From Frisch et al. (2011), the LIC has a thermal pressure of $P_T/k=3200^{+2000}_{-1100}$ cm−3 K, which implies a thermal pressure difference between the LICs and the LHB of Δ(PT/k) ${\sim }\,7400^{+1530}_{-2270}$ cm−3 K.

Figure 1.

Figure 1. Distribution of the neutral ISM in the direction of l  ∼  144°, b  ∼  0° from the analysis of Lallement et al. (2014) where purple indicates the denser material of the LC wall and the distances are in parsecs. The inversion mapping is based on reddening and absorption line data, and is currently limited by the number of sample stars to a resolution of ∼15 pc.

Standard image High-resolution image

2.3. Interstellar Magnetic Field Pressure

Assuming that the LHB plasma has been swept relatively free of magnetic fields, the magnetic pressure within the LICs may compensate for the difference in thermal pressures. (Also assuming that other contributors to the total interstellar pressure either similarly affect both the LICs and the LHB or are negligible.) The magnetic field pressure is PB/k = 288 B2 cm−3 K where B is in μG. The required LIC magnetic field for pressure equilibrium is then

Equation (3)

From Burlaga & Ness (2014), the interstellar magnetic field in the LIC ranges from 3.8 μG to 5.9 μG with an average of 4.9 ± 0.5μG.

3. DISCUSSION

Assuming that the magnetic field measured by Voyager 1 is representative of the LICs, the agreement between the thermal plus magnetic field pressure of the LICs and the thermal pressure of the LHB plasma alleviates a long-standing inconsistency in the LHB model leading to a self-consistent picture of the local ISM. The pressure equilibrium explains why there has been no evidence for either the expansion or contraction of the LICs (Frisch et al. 2011), and is consistent with the picture of a rather sedate and aged bubble of hot plasma. (Note, however, that Gry & Jenkins 2014 have a new interpretation that the LICs are actually one monolithic structure which is compressing in the ram and anti-ram directions and expanding perpendicularly to the flow.), The results of G14 indicate that ∼40% of the minimum observed 1/4 keV flux in the Galactic plane originates as SWCX emission within the heliosphere and magnetosheath but still requires X-ray emitting plasma within the LC, even at low Galactic latitudes (the required LC/LHB plasma at higher latitudes is relatively unchanged). In short, contrary to rumors of its demise (e.g., Welsh & Shelton 2009), the LHB still lives.

A recent paper by Smith et al. (2014) came to a similar conclusion about the continued existence of the LHB using a spectral analysis of Diffuse X-ray Spectrometer (DXS; Sanders et al. 2001) and Cosmic Hot Interstellar Plasma Spectrometer (Hurwitz et al. 2005) data. Using their own model for the SWCX spectrum plus a thermal emission component they fit the high spectral resolution DXS 1/4 keV band data deriving the result that 26% ± 4% of the observed emission was from the LHB (the DXS data were from near the Galactic plane). While the fraction of the observed emission from the LHB is less than half of the amount found by G14, the Smith et al. (2014) results suggest a thermal pressure considerably closer to that of the LICs, and thus less reliant on a differential in magnetic field strength for pressure equilibrium.

The local ISM is clearly more complicated than simply filling the LC with a hot plasma. For instance, there are directions where neither neutral material nor X-ray emitting plasma are observed, such as the Region of Bizarre Emptiness (RBE; Cox 1998) although a large filling factor of H ii region-like-gas, such as observed in the direction of β CMa (Gry et al. 1985) could be responsible.

While we report here the characterization of the LHB, SWCX, LIC, and thermal and magnetic field pressures based on the data in one direction, these results can reasonably be extended to the entire LHB. As mentioned above, the LHB thermal pressure that is derived here is in very good agreement with the thermal pressure derived in Puspitarini et al. (2014), which considered the LHB's extent in the LC. With their scaling, the plasma required for the LHB emission fits inside their three-dimensional mapping of the LC.

We are grateful to Ed Jenkins for interesting and helpful insights into these results. This work has been supported in part by NASA's award #NNX11AF04G.

Footnotes

Please wait… references are loading.
10.1088/2041-8205/791/1/L14