Issue 
A&A
Volume 611, March 2018



Article Number  L5  
Number of page(s)  6  
Section  Letters to the Editor  
DOI  https://doi.org/10.1051/00046361/201832637  
Published online  05 April 2018 
Letter to the Editor
The Local Bubble: a magnetic veil to our Galaxy
^{1}
Department of Astrophysics/IMAPP, Radboud University,
PO Box 9010,
6500 GL
Nijmegen, The Netherlands
email: m.alves@astro.ru.nl
^{2}
Institut de Recherche en Astrophysique et Planétologie (IRAP), CNRS, Université de Toulouse, CNRS,
9 Avenue du Colonel Roche,
BP 44346,
31028
Toulouse Cedex 4, France
^{3}
École normale supérieure/LERMA, Observatoire de Paris, Sorbonne Université, Université PSL, CNRS,
Paris, France
Received:
13
January
2018
Accepted:
6
March
2018
The magnetic field in the local interstellar medium does not follow the largescale 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, offcentred from the Sun. We fit the model to Planck dust polarized emission observations within 30^{°} 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 threedimensional (3D) structure of the magnetic field in the local interstellar medium, which is a most awaited input for largescale Galactic magnetic field models.
Key words: ISM: bubbles / ISM: magnetic fields / ISM: individual objects: Local Bubble
© ESO 2018
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0;), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 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 largescale 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 largescale 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 largescale 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 frozenin 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 largescale 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 largescale magnetic field were related to the Local Bubble shell by Leroy (1999), who, based on stellar polarization observations, found coherent magnetic field orientations over the distance range 60–150 pc. Recently, Planck Collaboration Int. XLIV (2016) 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 largescale 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).
2 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 sweptup matter and field lines are squeezed into a thin shell (whose actual thickness is irrelevant) that follows the surface of the Local Bubble. Accordingly, the sweptup 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 imaginethat 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 frozenin approximation, the present vector potential, A, is related to the initial vector potential, A_{0} , through (2)
Since the initial magnetic field is uniform, we can adopt (3)
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 (4)
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∕rsinθ) (∂∕∂ϕ) 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 (5)
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 (6)
which expresses the fact that the sweptup magnetic field is bothtangent 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 directionof B_{0} and on the shape of the bubble.
As we aim to model the ordered magnetic field in the shell of the LocalBubble, we do not attempt to reproduce its complex structure (Capitanio et al. 2017; Liu et al. 2017). Rather, we consider that its largescale 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 y and 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 semiaxes, a_{ell} , b_{ell} , andc_{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.
3 Fit to Planck observations
We takethe latest publicly available Planck maps of the Stokes parameters Q and U at 353 GHz (Planck Legacy Archive^{1}, Planck Collaboration VIII 2016). For the total intensity I we use the model map derived by Planck Collaboration XI (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 LocalBubble only, that is, it does not include any component from the largescale 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 ϕ_{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 highb 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 Int. XLIV (2016) to produce allsky 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 pythonmpfit 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 MJy sr^{−1} and that from the turbulent magnetic field is .
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 subsets of data, with and without spectral mismatch correction (Planck Collaboration ES 2015; Planck Collaboration VIII 2016), 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 largescale 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 semiminor axis of the bubble. In other words, the Sun is significantly offcentre.
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 values 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 subsets, 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 Int. XLIV (2016), 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 Int. XLIV (2016) 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.
Results of fits of our magnetic field model to Planck polarization observations.
Fig. 1 Observed (top) and fitted (bottom) Q (left) and U (right) Stokes parameter maps in orthographic projection, in units of MJy sr^{−1} . Each pair of maps shows the sky divided into its two Galactic hemispheres: north (left) and south (right). The vertical full line in each map indicates ℓ = 0^{°}. The grey regions correspond to the areas that are masked. The texture seen in the model maps, overlaid on the regular patterns, arises from the intensity I. Planck observations show additional structure from the turbulent component of the magnetic field, not included in our model. 
4 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 magnetic 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 Int. XLIV (2016), 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 largescale Galactic magnetic field (e.g. Terral & Ferrière 2017). Further constraints come from Xray 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 Xray 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 Xray 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 stateoftheart Galactic magnetic field models. However, it is undoubtedly important to consider this local contribution in order to model polarization observations outside of the Galacticplane (as recently pointed out by Planck Collaboration Int. XLII 2016). 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.
Fig. 2 The shape of the Local Bubble in greatcircle 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 rescaled by a factor of 1/(100 pc). The dashed blue curves represent the sky area where b < 60^{°}, which is notincluded in the analysis. 
Acknowledgements
We thank the referee for the useful comments. We also thank W. Liu for kindly providing the allsky map of the structure of the Local Bubble. The results in this paper were derived using the HEALPix package (Górski et al. 2005) and the python mpfit routine written by M. Cappellari. This research is supported by the Agence Nationale de la Recherche (project BxB: ANR17CE310022).
Appendix A The magnetic field model
In this appendix, we describe in more detail some of the steps involved in the derivation of the magnetic field model presented in Sect. 2, in particular, the justification of Eq. (2) and how to go from Eq. (2) to Eq. (4). Except for a few new physical quantities defined below, all other quantities are as defined in Sect. 2.
Our magnetic field model relies on the frozenin approximation, which implies that the magnetic field is coupled to the matter, thatis, magnetic field lines follow the matter in its motions. Generally, the frozenin approximation is given by the induction equation in the limit of zero magnetic diffusion: (A.1)
where B = ∇×A is the magnetic field, A the vector potential, and V the velocity field. “Decurling” Eq. (A.1) leads to the evolution equation for A: (A.2)
with f an arbitrary scalar function. Choosing f = −V ⋅A, we can rewrite Eq. (A.2) in the form (A.3)
Finally, integrating Eq. (A.3) in time leads to Eq. (2) (also see Parker 1970).
The present vector potential is obtained by introducing Eq. (1) (valid for purely radial motions) and Eq. (3) (valid for a uniform initial magnetic field) into Eq. (2): (A.4)
The final magnetic field then follows from B = ∇×A: (A.6)
Now, assume that all the matter that was initially interior to the current shell has been swept up into an infinitely thin shell presently located at r = R_{sh}(θ, ϕ). Then, for r ≤ R_{sh}, the initial radius can be written as (A.7)
and its spatial derivatives as (A.8)
where h is the step function and δ the delta function. Noting that ∂r_{0}∕∂r ≫ r_{0}∕r in the shell, we introduce Eqs. (A.5), (A.7), and (A.8) into Eq. (A.6) and retain only the leading terms in ∂r_{0}∕∂r, whereupon the first term of Eq. (A.6) becomes (A.9)
the second term becomes (A.10)
and their sum gives Eq. (4).
References
 Berdyugin, A., Piirola, V., & Teerikorpi, P. 2014, A&A, 561, A24 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 BICEP2/KeckCollaboration, Planck Collaboration, Ade, P. A. R., et al. 2015, Phys. Rev. Lett., 114, 101301 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Bronfman, L., Casassus, S., May, J., & Nyman, L.Å. 2000, A&A, 358, 521 [Google Scholar]
 Capitanio, L., Lallement, R., Vergely, J. L., Elyajouri, M., & MonrealIbero, A. 2017, A&A, 606, A65 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Cox, D. P., & Reynolds, R. J. 1987, ARA&A, 25, 303 [NASA ADS] [CrossRef] [Google Scholar]
 Ferrière, K. M. 2001, Rev. Mod. Phys., 73, 1031 [NASA ADS] [CrossRef] [Google Scholar]
 Frisch, P. C., Andersson, B.G., Berdyugin, A., et al. 2012, ApJ, 760, 106 [NASA ADS] [CrossRef] [Google Scholar]
 Fuchs, B., Breitschwerdt, D., de Avillez, M. A., Dettbarn, C., & Flynn, C. 2006, MNRAS, 373, 993 [NASA ADS] [CrossRef] [Google Scholar]
 Górski, K. M., Hivon, E., Banday, A. J., et al. 2005, ApJ, 622, 759 [NASA ADS] [CrossRef] [Google Scholar]
 Han, J. L., Manchester, R. N., Lyne, A. G., Qiao, G. J., & van Straten W. 2006, ApJ, 642, 868 [NASA ADS] [CrossRef] [Google Scholar]
 Heiles, C. 1996, ApJ, 462, 316 [NASA ADS] [CrossRef] [Google Scholar]
 Heiles, C. 1998, in IAU Colloq. 166: The Local Bubble and Beyond, eds. D. Breitschwerdt, M. J. Freyberg, & J. Truemper (Berlin: Springer Verlag), Lect. Notes Phys., 506, 229 [NASA ADS] [Google Scholar]
 Heiles, C., & Haverkorn, M. 2012, Space Sci. Rev., 166, 293 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Jaffe, T. R., Ferrière, K. M., Banday, A. J., et al. 2013, MNRAS, 431, 683 [NASA ADS] [CrossRef] [Google Scholar]
 Jansson, R., & Farrar, G. R. 2012, ApJ, 757, 14 [NASA ADS] [CrossRef] [Google Scholar]
 Jelić, V., Zaroubi, S., Labropoulos, P., et al. 2010, MNRAS, 409, 1647 [NASA ADS] [CrossRef] [Google Scholar]
 Joshi, Y. C. 2007, MNRAS, 378, 768 [NASA ADS] [CrossRef] [Google Scholar]
 Kos, J., Zwitter, T., Wyse, R., et al. 2014, Science, 345, 791 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Lallement, R., Vergely, J.L., Valette, B., et al. 2014, A&A, 561, A91 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Leroy, J. L. 1999, A&A, 346, 955 [NASA ADS] [Google Scholar]
 Liu, W., Chiao, M., Collier, M. R., et al. 2017, ApJ, 834, 33 [NASA ADS] [CrossRef] [Google Scholar]
 Parker, E. N. 1970, ApJ, 162, 665 [NASA ADS] [CrossRef] [Google Scholar]
 Planck Collaboration XI. 2014, A&A, 571, A11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration ES. 2015, http://wiki.cosmos.esa.int/planckpla/index.php/Main_Page (ESA) [Google Scholar]
 Planck Collaboration VIII. 2016, A&A, 594, A8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration Int. XLII. 2016, A&A, 596, A103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration Int. XLIV. 2016, A&A, 596, A105 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Puspitarini, L., Lallement, R., Vergely, J.L., & Snowden, S. L. 2014, A&A, 566, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Santos, F. P., Corradi, W., & Reis, W. 2011, ApJ, 728, 104 [NASA ADS] [CrossRef] [Google Scholar]
 Shelton, R. L. 1998, ApJ, 504, 785 [NASA ADS] [CrossRef] [Google Scholar]
 Sun, X.H., & Reich, W. 2010, Res. Astron. Astrophys., 10, 1287 [NASA ADS] [CrossRef] [Google Scholar]
 Terral, P., & Ferrière, K. 2017, A&A, 600, A29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Vansyngel, F., Boulanger, F., Ghosh, T., et al. 2017, A&A, 603, A62 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
All Tables
All Figures
Fig. 1 Observed (top) and fitted (bottom) Q (left) and U (right) Stokes parameter maps in orthographic projection, in units of MJy sr^{−1} . Each pair of maps shows the sky divided into its two Galactic hemispheres: north (left) and south (right). The vertical full line in each map indicates ℓ = 0^{°}. The grey regions correspond to the areas that are masked. The texture seen in the model maps, overlaid on the regular patterns, arises from the intensity I. Planck observations show additional structure from the turbulent component of the magnetic field, not included in our model. 

In the text 
Fig. 2 The shape of the Local Bubble in greatcircle 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 rescaled by a factor of 1/(100 pc). The dashed blue curves represent the sky area where b < 60^{°}, which is notincluded in the analysis. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.