The Local Bubble: a magnetic veil to our Galaxy

The magnetic field in the local interstellar medium does not follow the large-scale Galactic magnetic field. The local magnetic field has probably been distorted by the Local Bubble, a cavity of hot ionized gas extending all around the Sun and surrounded by a shell of cold neutral gas and dust. However, so far no conclusive association between the local magnetic field and the Local Bubble has been established. Here we develop an analytical model for the magnetic field in the shell of the Local Bubble, which we represent as an inclined spheroid, off-centred from the Sun. We fit the model to Planck dust polarized emission observations within 30 deg of the Galactic poles. We find a solution that is consistent with a highly deformed magnetic field, with significantly different directions towards the north and south Galactic poles. This work sets a methodological framework for modelling the three-dimensional (3D) structure of the magnetic field in the local interstellar medium, which is a most awaited input for large-scale Galactic magnetic field models.


Introduction
The interstellar medium (ISM) of our Galaxy is threaded by a pervasive magnetic field that plays an important role in many astrophysical processes (see e.g. Ferrière 2001, Heiles & Haverkorn 2012. Various models of the large-scale Galactic magnetic field have been developed throughout the years (e.g. Han et al. 2006, Sun & Reich 2010, Jansson & Farrar 2012, Jaffe et al. 2013, Terral & Ferrière 2017. These aim at describing the large-scale structure of the magnetic field, commonly divided into disk and halo components. The models are mostly constrained by observations of synchrotron emission and Faraday rotation. However, the observations include several local structures, which appear large on the sky because of their proximity, but are not part of the large-scale magnetic field. The most important such structure is the Local Bubble, an interstellar bubble located around the Sun and extending out to ∼ 60 pc towards the Galactic centre and 250 pc towards the north and south Galactic poles (Puspitarini et al. 2014). The Local Bubble is thought to be the result of supernova explosions, which swept out a cavity of hot ionized gas and pushed most of the evacuated matter, together with the frozen-in magnetic field, into a dense shell of cold neutral gas and dust (Cox & Reynolds 1987;Shelton 1998;Fuchs et al. 2006;Lallement et al. 2014).
Several observational studies have shown that the magnetic field in the local ISM does not follow the large-scale Galactic magnetic field (e.g. Heiles 1998, Leroy 1999, Santos et al. 2011, Frisch et al. 2012, Berdyugin et al. 2014, which points approximately towards ( , b) = (83 • , 0 • ) at the Sun's position (Heiles 1996). Here, and b denote Galactic longitude and latitude, respectively. These departures from the large-scale magnetic field were related to the Local Bubble shell by Leroy Email: m.alves@astro.ru.nl (1999), who, based on stellar polarization observations, found coherent magnetic field orientations over the distance range 60 -150 pc. Recently, Planck Collaboration et al. (2016c) analysed the Planck maps of dust polarized emission towards the southern Galactic cap (b < −60 • ), and found that these could be represented by a uniform magnetic field pointing in the direction ( , b) = (70 • , 24 • ), thus with a significant vertical component. At these high Galactic latitudes, the observed dust emission is most likely associated with the shell of the Local Bubble (Kos et al. 2014). Despite these different studies, no definitive relation between the deformed local magnetic field and the Local Bubble has so far been established.
In this letter, we present the first model of the magnetic field in the Local Bubble, based on a new analytical solution and on Planck observations of dust polarized emission. This work is a stepping stone towards modelling the 3D structure of the magnetic field in the Solar neighbourhood. Such a model will contribute to our understanding of the large-scale Galactic magnetic field and will also lead to a more accurate modelling of the Galactic polarized foreground emission, which contaminates observations of the B modes of the cosmic microwave background (BICEP2/Keck Collaboration et al. 2015) and of Hi emission from the epoch of reionization (Jelić et al. 2010).

The magnetic field in the Local Bubble
Observations of dust polarized emission trace the orientation of the magnetic field averaged along the line of sight, without providing any information on the magnetic field strength. Hence, we do not attempt to model variations in the magnetic field along the line of sight, nor are we concerned with the magnetic field strength. Instead, we assume that all the swept-up matter and field lines are squeezed into a thin shell (whose actual thickness Article number, page 1 of 6 arXiv:1803.05251v1 [astro-ph.GA] 14 Mar 2018 is irrelevant) that follows the surface of the Local Bubble. Accordingly, the swept-up magnetic field is required to be tangent throughout to the surface of the Local Bubble. To fully define the magnetic field orientation, we make the additional assumption that the expansion motions that created the Local Bubble are purely radial from a single explosion centre. However, these motions are not required to have spherical symmetry.
In view of the above assumptions, we now derive an analytical expression for the magnetic field, B, in a very thin shell surrounding a bubble of arbitrary shape (see Appendix A for more details). Consider an initially homogeneous medium with uniform magnetic field B 0 , and imagine that supernova explosions occurred at a point O of this medium. Let us define a cartesian coordinate system (x, y, z) centered on O, and let us denote the associated spherical coordinates by (r, θ, φ). Since the expansion motions driven by the explosion are assumed to be purely radial (i.e. along the unit vector e r ), the initial position of a particle presently at position r = r e r can be written as (1) To derive the general expression of the magnetic field deformed by the explosion, we use the vector potential, A, defined such that B = ∇ × A. In the frozen-in approximation, the present vector potential, A, is related to the initial vector potential, A 0 , through Since the initial magnetic field is uniform, we can adopt Introducing Eqs. (1) and (3) into Eq.
(2) and taking B = ∇ × A, we obtain a lengthy and unwieldy equation for B. However, under our assumption that all the field lines swept up by the explosion are confined into a very thin shell, the general expression of B can be greatly simplified to where R sh is the shell radius (a function of (θ, φ)), B 0t = B 0θ e θ + B 0φ e φ is the orthoradial (or transverse) component of the initial magnetic field, and ∇ t = (1/r) (∂/∂θ) e θ + (1/r sin θ) (∂/∂φ) e φ is the orthoradial component of the gradient operator. The factor (r 0 /r) in Eq. (4) arises from spherical divergence. For a very thin shell, the second factor is very large, but this is irrelevant here because we are only interested in the direction of the magnetic field, not in its strength.
For numerical convenience, we rewrite Eq. (4) in the form where n = (e r − ∇ t R sh )/||e r − ∇ t R sh || is the unit vector normal to the surface of the shell. The advantage of writing B in terms of the unit vector n is that, in contrast to ∇ t R sh , n is independent of the reference frame. Interestingly, Eq. (5) can also be written as which expresses the fact that the swept-up magnetic field is both tangent to the shell surface (as required) and contained in the plane of the initial magnetic field, B 0 , and the direction of the expansion motion, e r (as expected). Thus the direction of B depends only on the direction of B 0 and on the shape of the bubble.
As we aim to model the ordered magnetic field in the shell of the Local Bubble, we do not attempt to reproduce its complex structure (Capitanio et al. 2017;Liu et al. 2017). Rather, we consider that its large-scale shape can be captured by a generic geometrical form, for which we adopt a simple ellipsoid centred on the explosion centre. The latter is not required to coincide with the Sun's position. For a given initial magnetic field and a given bubble shape, the magnetic field in the shell has a welldetermined configuration, but our view of its projection onto the plane of the sky depends on the position of the explosion centre relative to the Sun. In practice, to fit our magnetic field model to the Planck observations, we need to perform a change of coordinates from the explosion centre frame, (x, y, z), where our model equations are written, to the Sun's frame, (x , y , z ), where the observations are made. Here, we orient these two frames such that the x-and x -axes point towards the Galactic centre, the yand y -axes point towards = 90 • , and the z-and z -axes point towards b = 90 • .
Altogether, our magnetic field model involves 12 free parameters: 9 describing the size, the orientation, and the position of the ellipsoid, 2 defining the direction of the initial magnetic field, and 1 relating the magnetic field orientation to the observed dust polarized emission. More specifically, the 3 parameters giving the size of the ellipsoid are its semi-axes, a ell , b ell , and c ell . The 3 parameters giving the orientation of the ellipsoid are the standard Euler angles, θ ell (nutation), ψ ell (precession), and φ ell (intrinsic), defined in the explosion centre frame. The 3 parameters giving the position of the ellipsoid are the coordinates (δ x , δ y , δ z ) of its centre in the Sun's frame. The 2 parameters defining the direction of the initial magnetic field B 0 are its Galactic longitude, 0 , and latitude, b 0 . The last parameter is the maximum dust polarization fraction, p 0 , i.e. the polarization fraction obtained when the magnetic field lies exactly in the plane of the sky.

Fit to Planck observations
We take the latest publicly available Planck maps of the Stokes parameters Q and U at 353 GHz (Planck Legacy Archive 1 , Planck Collaboration et al. 2016a). For the total intensity I we use the model map derived by Planck Collaboration et al. (2014), obtained by fitting observations that have been corrected for the cosmic microwave background anisotropies and for zodiacal light. These data are smoothed to 1 • angular resolution.
Our model concerns the magnetic field in the shell of the Local Bubble only, that is, it does not include any component from the large-scale Galactic magnetic field. Therefore, to limit the contribution from the Galactic disk to the observed dust polarized emission, we restrict our analysis to the Galactic polar caps, |b| > 60 • . Given this restriction, we do not attempt to fit all the parameters in the model. Instead, we present a first illustration of our method based on a limited choice for the input parameters. First, we neglect the vertical shift of the explosion centre relative to the Sun, that is, we set δ z = 0, based on the argument that the Sun is located near the Galactic plane (Joshi 2007) and that massive stars, the most likely supernova progenitors, are also distributed close to the Galactic plane (Bronfman et al. 2000). Second, we consider that the shape of the Local Bubble can be described by a spheroid, and accordingly we set a ell = b ell . The exact value of φ ell is thus irrelevant, and for simplicity we adopt [deg] 71 ± 11 −16 ± 7 12 ± 1 set=1 set=1 2.7 ± 0.6 216 ± 47 30 ± 50 set=0 0.17 ± 0.3 0.56 ± 0.2 set=0 φ ell = 0 • . Our restriction to a spheroid follows from the observation that the Local Bubble is elongated towards the Galactic poles (see Sect. 1), together with the limitations of our analysis, which relies exclusively on high-|b| observations and is, therefore, unable to provide separate constraints on the Local Bubble dimensions in the Galactic plane. Finally, as dust polarized emission observations do not give any information about distances, we work with normalized sizes, with our reference length set by a ell = 1.
Under the above assumptions, our magnetic field model is left with 8 free parameters out of the initial 12. We follow the methodology detailed in Planck Collaboration et al. (2016c) to produce all-sky mock Stokes parameter maps from our 3D magnetic field model. In particular, we compute normalized Stokes parameters, Q/I and U/I, which depend on the magnetic field orientation only, and then multiply them by the intensity map I. We use the python mpfit routine to derive the set of parameter values that best fit the observations. We take into account the noise in the Planck Q and U data, which we add in quadrature with a contribution from the turbulent magnetic field component, which is not accounted for in our model. The latter is estimated from the power spectra of the Planck dust polarization maps at 353 GHz using data simulations from Vansyngel et al. (2017). In the region of the sky under study, the mean contribution from the statistical noise is σ noise Q,U = 1.4 × 10 −3 MJy sr −1 and that from the turbulent magnetic field is σ turb Q,U /I = 0.055. Table 1 shows the results of the fit to the Planck observations, corresponding to a reduced χ 2 of 0.8. The fitted Q and U Stokes parameter maps are displayed in Fig. 1, along with the observational Planck maps. The uncertainties in the different parameters are derived from the standard deviations of the parameter values obtained from similar fits to 12 different Planck sub-sets of data, with and without spectral mismatch correction (Planck Collaboration ES 2015;Planck Collaboration et al. 2016a), and thus account for the uncertainties due to systematic effects in the observations. The latter are significantly larger than the uncertainties resulting from the fitting procedure. This is particularly the case for ψ ell and θ ell , whose uncertainties further reflect the degeneracy between these parameters and the position of the spheroid.
We find that the initial magnetic field points in the direction ( 0 , b 0 ) = (71 • , −16 • ), which is close to that of the large-scale Galactic magnetic field at the Sun's position, ( , b) = (83 • , 0 • ) (see Sect. 1). The shape of the Local Bubble is consistent with a prolate spheroid of aspect ratio 2.7, whose long axis points towards ( , b) = (216 • , 60 • ). The explosion centre is located along the direction = 73 • , similar to that of the initial magnetic field, 0 = 71 • , and at a normalized distance of 0.58 from the Sun, corresponding to a large proportion (58 %) of the semi-minor axis of the bubble. In other words, the Sun is significantly off-centre. Our model indicates that the magnetic field points on average towards ( , b) = (70 • ±11 • , 43 • ±8 • ) and (74 • ±8 • , −14 • ±18 • ) in the north and south Galactic polar caps, respectively. These val- Fig. 2. The shape of the Local Bubble in great-circle cuts through the Galactic poles along two longitudes: 0 (top) and 0 + 90 • (bottom). The blue curves correspond to our model and the green curves to the data from Liu et al. (2017). The latter are re-scaled by a factor of 1/(100 pc). The dashed blue curves represent the sky area where |b| < 60 • , which is not included in the analysis. ues are obtained by averaging the magnetic field vector within 30 • from the poles. We perform the same averaging for each of the magnetic field solutions obtained from the 12 different Planck data sub-sets, and take their standard deviations as the uncertainties in and b. The direction of the magnetic field in the southern polar cap turns out to be different from that derived by Planck Collaboration et al. (2016c), who found ( , b) = (70 • , 24 • ), although they obtained the same maximum dust polarization fraction, p 0 = 12 %. The discrepancy probably finds its roots in the fundamental difference between both approaches: while Planck Collaboration et al. (2016c) model the polarization data solely based on variations over the polar cap in the projection of a presumably uniform magnetic field, our model includes both projection effects and intrinsic variations in the magnetic field orientation. Since these variations rely on a physical bubble model, our results might be more reliable. In any case, we emphasise that the results are model dependent and that our new results should be taken with particular caution, given the narrow angular range of the observations retained in our analysis added to the stringent assumptions underlying our bubble model.

Discussion and conclusions
Fitting our model of the magnetic field in the shell of the Local Bubble to Planck observations of dust polarized emission, we obtain a solution with a good reduced χ 2 . Although this solution is only indicative, it strongly suggests that the local mag-netic field as probed by Planck observations is highly distorted. In particular, the directions of the magnetic field in the north and south Galactic polar caps are found to differ by 57 • . The discrepancy between our results and those of Planck Collaboration et al. (2016c), regarding the direction of the magnetic field in the south polar cap shows that it is essential to consider physically motivated models to reproduce the observations.
In order to effectively constrain our magnetic field model, it is necessary to extend the sky area over which the fit is performed. However, to do so we need to take into account the polarized signal arising from the Galactic disk and possibly from the halo. In principle, this can be done by combining our geometrical model of the Local Bubble with existing models of the large-scale Galactic magnetic field (e.g. Terral & Ferrière 2017). Further constraints come from X-ray observations and stellar astrometric data, with which the 3D structure of the Local Bubble can be estimated (Capitanio et al. 2017;Liu et al. 2017). These observations are not taken into account in the present analysis. As it turns out, we find that the shape of the Local Bubble derived from our fit differs from that inferred from X-ray data by Liu et al. (2017) (as illustrated in Fig. 2). There are several reasons for this discrepancy, including the assumption by Liu et al. (2017) that the measured X-ray intensity is simply proportional to the size of the bubble in the considered direction and our own assumption that the centre of the explosion coincides with the geometric centre of the bubble. When incorporating observations of the Local Bubble structure into our modelling, we will need to relax this latter assumption.
The local magnetic field is currently not taken into account in state-of-the-art Galactic magnetic field models. However, it is undoubtedly important to consider this local contribution in order to model polarization observations outside of the Galactic plane (as recently pointed out by Planck Collaboration et al. 2016b). Here we illustrate this point using dust polarized emission observations. Clearly the same idea could be applied to other tracers of the magnetic field, including Faraday rotation measures and perhaps even synchrotron emission data.