An analytical model for fully developed three-dimensional incompressible turbulence was recently proposed in the hydrodynamics community, based on the concept of multiplicative chaos. It consists of a random field represented by means of a stochastic integral, which, with only a few parameters, shares many properties with experimental and numerical turbulence, including in particular energy transfer through scales (the cascade) and intermittency (non-Gaussianity) which is most conveniently controlled with a single parameter. Here, we propose three models extending this approach to MHD turbulence. Our formulae provide physically motivated 3D models of a turbulent velocity field and magnetic field coupled together. Besides its theoretical value, this work is meant to provide a tool for observers: a dozen of physically meaningful free parameters enter the description, which is useful to characterize astrophysical data.

The Sun is embedded in the so-called Local Bubble (LB) -- a cavity of hot plasma created by supernova explosions and surrounded by a shell of cold, dusty gas. Knowing the local distortion of the Galactic magnetic field associated with the LB is critical for the modeling of interstellar polarization data at high Galactic latitudes. In this his paper, we relate the structure of the Galactic magnetic field on the LB scale to three-dimensional (3D) maps of the local interstellar medium (ISM). First, we extracted the geometry of the LB shell, its inner surface, in particular from 3D dust extinction maps of the local ISM. We expanded the shell inner surface in spherical harmonics, up to a variable maximum multipole degree, which enabled us to control the level of complexity for the modeled surface. Next, we applied an analytical model for the ordered magnetic field in the shell to the modeled shell surface. This magnetic field model was successfully fitted to the \textit{Planck} 353~GHz dust polarized emission maps over the Galactic polar caps. For each polar cap, the direction of the mean magnetic field derived from dust polarization (together with the prior that the field points toward longitude $90^\circ \pm 90^\circ$) is found to be consistent with the Faraday spectra of the nearby diffuse synchrotron emission. Our work presents a new approach to modeling the local structure of the Galactic magnetic field. We expect our methodology and our results to be useful both in modeling the local ISM as traced by its different components and in modeling the dust polarized emission, which is a long-awaited input for studies of the polarized foregrounds for cosmic microwave background.

Reliable identification of the origin of the high-energy non-thermal emission from the Galactic Center (GC) is not achievable without adequate consideration of the ambient conditions such as the magnetic field configuration or gas distribution. In a first step, we present a model that can explain the diffuse gamma-ray emission as measured by H.E.S.S. for small longitudes in the Galactic Center region but comes to grief with higher longitudes. The model is given via the solution of a transport equation that allows for a radial dependency of the mass distribution. In order to move from this semi-analytical approximation toward a full understanding of the PeVatron signature, we present a new 3D analytical model of the gas distribution in the Central Molecular Zone (CMZ). Furthermore, we derive for the first time a 3D model of the magnetic field configuration and strength in the CMZ, which is analytical and divergence-free. The model is built via a combination of a model for the diffuse inter-cloud medium, local molecular clouds and non-thermal filaments at which local information are based on investigations from previous works and the molecular gas density. It can be shown that without an efficient longitudinal CR entrapment, a single source at the center does not facilely suffice the diffuse gamma-ray detection. Further, we show that the new magnetic field model GBFD19 is compatible with recent polarization data and has a significant impact on the longitudinal profiles of CR propagation.

The magnetic field in the Galactic Center (GC) is challenging to measure in detail due to its large distance in combination with a difficult foreground estimate. More precise knowledge is, however, required in particular in the context of proper modeling of Cosmic Ray (CR) transport. A detailed model of the magnetic field in the GC region, therefore, is not only required in order to understand the origin of CRs but also for appropriate modeling of electron transport with the resulted synchrotron spectra, as well as background estimates for dark matter searches in this region. For this kind of modeling, the gas distribution is equally important, as the production of high-energy gamma-ray via hadronic interaction depends on it. In this work, we are putting together and interpreting available data from the Galactic Center region in order to model both the gas distribution and the magnetic field configuration. We derive the configuration and strength of the total magnetic field in the GC via a combination of an analytical and divergence-free model for the diffuse inter-cloud medium and information from local molecular clouds and non-thermal filaments which are based on investigations from previous works and the molecular gas density. We present a comprehensive model of the magnetic field in the central 200 pc of the Galaxy. It can be shown that the model is compatible with recent polarization data. In first test simulations of cosmic-ray transport, we can show that using this more detailed magnetic field configuration in comparison to large-scale magnetic field configurations that discard the direct information from the GC has a significant impact on the propagation profiles of cosmic-rays.

SPICA, the cryogenic infrared space telescope recently pre-selected for a `Phase A' concept study as one of the three remaining candidates for ESA's fifth medium class (M5) mission, is foreseen to include a far-infrared polarimetric imager (SPICA-POL, now called B-BOP), which would offer a unique opportunity to resolve major issues in our understanding of the nearby, cold magnetized Universe. This paper presents an overview of the main science drivers for B-BOP, including high dynamic range polarimetric imaging of the cold interstellar medium (ISM) in both our Milky Way and nearby galaxies. Thanks to a cooled telescope, B-BOP will deliver wide-field 100-350 micron images of linearly polarized dust emission in Stokes Q and U with a resolution, signal-to-noise ratio, and both intensity and spatial dynamic ranges comparable to those achieved by Herschel images of the cold ISM in total intensity (Stokes I). The B-BOP 200 micron images will also have a factor ~30 higher resolution than Planck polarization data. This will make B-BOP a unique tool for characterizing the statistical properties of the magnetized interstellar medium and probing the role of magnetic fields in the formation and evolution of the interstellar web of dusty molecular filaments giving birth to most stars in our Galaxy. B-BOP will also be a powerful instrument for studying the magnetism of nearby galaxies and testing galactic dynamo models, constraining the physics of dust grain alignment, informing the problem of the interaction of cosmic rays with molecular clouds, tracing magnetic fields in the inner layers of protoplanetary disks, and monitoring accretion bursts in embedded protostars.

We present a statistical study of the relative orientation in the plane of the sky between interstellar magnetic fields and filaments hosting cold clumps. For the first time, we consider both the density of the environment and the density contrast between the filaments and their environment. Moreover, we geometrically distinguish between the clumps and the remaining portions of the filaments. We infer the magnetic field orientations in the filaments and in their environment from the Stokes parameters, assuming optically thin conditions. Thus, we analyze the relative orientations between filaments, embedded clumps, and internal and background magnetic fields, depending on the filament environment and evolutionary stages. We recover the previously observed trend for filaments in low column density environments to be aligned parallel to the background magnetic field; however, we find that this trend is significant only for low contrast filaments, whereas high contrast filaments tend to be randomly orientated with respect to the background magnetic field. Filaments in high column density environments do not globally show any preferential orientation, although low contrast filaments alone tend to have perpendicular relative orientation with respect to the background magnetic field. For a subsample of nearby filaments, for which volume densities can be derived, we find a clear transition in the relative orientation with increasing density, at $n_{\rm H} \sim 10^{3}~{\rm cm}^{-3}$, changing from mostly parallel to mostly perpendicular in the off-clump portions of filaments and from even to bimodal in the clumps. Our results confirm a strong interplay between interstellar magnetic fields and filaments during their formation and evolution.

We present 353 GHz full-sky maps of the polarization fraction $p$, angle $\psi$, and dispersion of angles $S$ of Galactic dust thermal emission produced from the 2018 release of Planck data. We confirm that the mean and maximum of $p$ decrease with increasing $N_H$. The uncertainty on the maximum polarization fraction, $p_\mathrm{max}=22.0$% at 80 arcmin resolution, is dominated by the uncertainty on the zero level in total intensity. The observed inverse behaviour between $p$ and $S$ is interpreted with models of the polarized sky that include effects from only the topology of the turbulent Galactic magnetic field. Thus, the statistical properties of $p$, $\psi$, and $S$ mostly reflect the structure of the magnetic field. Nevertheless, we search for potential signatures of varying grain alignment and dust properties. First, we analyse the product map $S \times p$, looking for residual trends. While $p$ decreases by a factor of 3--4 between $N_H=10^{20}$ cm$^{-2}$ and $N_H=2\times 10^{22}$ cm$^{-2}$, $S \times p$ decreases by only about 25%, a systematic trend observed in both the diffuse ISM and molecular clouds. Second, we find no systematic trend of $S \times p$ with the dust temperature, even though in the diffuse ISM lines of sight with high $p$ and low $S$ tend to have colder dust. We also compare Planck data with starlight polarization in the visible at high latitudes. The agreement in polarization angles is remarkable. Two polarization emission-to-extinction ratios that characterize dust optical properties depend only weakly on $N_H$ and converge towards the values previously determined for translucent lines of sight. We determine an upper limit for the polarization fraction in extinction of 13%, compatible with the $p_\mathrm{max}$ observed in emission. These results provide strong constraints for models of Galactic dust in diffuse gas.

Faraday tomography allows us to map diffuse polarized synchrotron emission from our Galaxy and use it to interpret the magnetic field in the interstellar medium (ISM). We have applied Faraday tomography to 60 observations from the LOFAR Two-meter Sky Survey (LOTSS) and produced a Faraday depth cube mosaic covering 568 square degrees at high Galactic latitudes, at 4.3' angular resolution and 1 rad m$^{-2}$ Faraday depth resolution, with a typical noise level of 50--100 $\mu$Jy per point spread function (PSF) per rotation measure spread function (RMSF) (40-80 mK RMSF$^{-1}$). While parts of the images are strongly affected by instrumental polarization, we observe diffuse polarized emission throughout most of the field, with typical brightness between 1 and 6 K RMSF$^{-1}$, and Faraday depths between $-7$ and +25 rad m$^{-2}$. We observed many new polarization features, some up to 15 degrees in length. These include two regions with very uniformly structured, linear gradients in the Faraday depth; we measured the steepness of these gradients as 2.6 and 13 rad m$^{-2}$ deg$^{-1}$. We also observed a relationship between one of the gradients and an HI filament in the local ISM. Other ISM tracers were also checked for correlations with our polarization data and none were found, but very little signal was seen in most tracers in this region. We conclude that the LOTSS data are very well suited for Faraday tomography, and that a full-scale survey with all the LOTSS data has the potential to reveal many new Galactic polarization features and map out diffuse Faraday depth structure across the entire northern hemisphere.

We present the first far infrared (FIR) dust emission polarization map covering the full extent Milky Way's Central molecular zone (CMZ). The data, obtained with the PILOT balloon-borne experiment, covers the Galactic Center region $-2\,^\circ<l<2\,^\circ$, $-4\,^\circ<b<3\,^\circ$ at a wavelength of 240 $\mu$m and an angular resolution $2.2\,'$. From our measured dust polarization angles, we infer a magnetic field orientation projected onto the plane of the sky that is remarkably ordered over the full extent of the CMZ, with an average tilt angle of $\simeq 22\,^\circ$ clockwise with respect to the Galactic plane. Our results confirm previous claims that the field traced by dust polarized emission is oriented nearly orthogonal to the field traced by GHz radio synchrotron emission in the Galactic Center region. The observed field structure is globally compatible with the latest Planck polarization data at 353 GHz and 217 GHz. Upon subtraction of the extended emission in our data, the mean field orientation that we obtain shows good agreement with the mean field orientation measured at higher angular resolution by the JCMT within the 20 km/s and 50 km/s molecular clouds. We find no evidence that the magnetic field orientation is related to the 100 pc twisted ring structure within the CMZ. We propose that the low polarization fraction in the Galactic Center region and the highly ordered projected field orientation can be reconciled if the field is strong, with a 3D geometry that is is mostly oriented $\simeq 15\,^\circ$ with respect to the line-of-sight towards the Galactic center. Assuming equipartition between the magnetic pressure and ram pressure, we obtain magnetic field strengths estimates as high as a few mG for several CMZ molecular clouds.

The study of polarized dust emission has become entwined with the analysis of the cosmic microwave background (CMB) polarization. We use new Planck maps to characterize Galactic dust emission as a foreground to the CMB polarization. We present Planck EE, BB, and TE power spectra of dust polarization at 353 GHz for six nested sky regions covering from 24 to 71 % of the sky. We present power-law fits to the angular power spectra, yielding evidence for statistically significant variations of the exponents over sky regions and a difference between the values for the EE and BB spectra. The TE correlation and E/B power asymmetry extend to low multipoles that were not included in earlier Planck polarization papers. We also report evidence for a positive TB dust signal. Combining data from Planck and WMAP, we determine the amplitudes and spectral energy distributions (SEDs) of polarized foregrounds, including the correlation between dust and synchrotron polarized emission, for the six sky regions as a function of multipole. This quantifies the challenge of the component separation procedure required for detecting the reionization and recombination peaks of primordial CMB B modes. The SED of polarized dust emission is fit well by a single-temperature modified blackbody emission law from 353 GHz to below 70 GHz. For a dust temperature of 19.6 K, the mean spectral index for dust polarization is $\beta_{\rm d}^{P} = 1.53\pm0.02 $. By fitting multi-frequency cross-spectra, we examine the correlation of the dust polarization maps across frequency. We find no evidence for decorrelation. If the Planck limit for the largest sky region applies to the smaller sky regions observed by sub-orbital experiments, then decorrelation might not be a problem for CMB experiments aiming at a primordial B-mode detection limit on the tensor-to-scalar ratio $r\simeq0.01$ at the recombination peak.

The "Square Kilometre Array" (SKA) is a large international radio telescope project characterised, as suggested by its name, by a total collecting area of approximately one square kilometre, and consisting of several interferometric arrays to observe at metric and centimetric wavelengths. The deployment of the SKA will take place in two sites, in South Africa and Australia, and in two successive phases. From its Phase 1, the SKA will be one of the most formidable scientific machines ever deployed by mankind, and by far the most impressive in terms of data throughput and required computing power. With the participation of almost 200 authors from forty research institutes and six private companies, the publication of this French SKA white paper illustrates the strong involvement in the SKA project of the French astronomical community and of a rapidly growing number of major scientific and technological players in the fields of Big Data, high performance computing, energy production and storage, as well as system integration.

As the next step toward an improved large scale Galactic magnetic field model, we present a simple comparison of polarised synchrotron and thermal dust emission on the Galactic plane. We find that the field configuration in our previous model that reproduces the polarised synchrotron is not compatible with the WMAP 94 GHz polarised emission data. In particular, the high degree of dust polarisation in the outer Galaxy (90deg < l < 270deg) implies that the fields in the dust-emitting regions are more ordered than the average of synchrotron-emitting regions. This new dust information allows us to constrain the spatial mixing of the coherent and random magnetic field components in the outer Galaxy. The inner Galaxy differs in polarisation degree and apparently requires a more complicated scenario than our current model. In the scenario that each interstellar component (including fields and now dust) follows a spiral arm modulation, as observed in external galaxies, the changing degree of ordering of the fields in dust-emitting regions may imply that the dust arms and the field component arms are shifted as a varying function of Galacto-centric radius. We discuss the implications for how the spiral arm compression affects the various components of the magnetised interstellar medium but conclude that improved data such as that expected from the Planck satellite will be required for a thorough analysis.

The polarization properties of radio sources at very low frequencies (<200 MHz) have not been widely measured, but the new generation of low-frequency radio telescopes, including the Low Frequency Array (LOFAR: a Square Kilometre Array Low pathfinder), now gives us the opportunity to investigate these properties. In this paper, we report on the preliminary development of a data reduction pipeline to carry out polarization processing and Faraday tomography for data from the LOFAR Two-meter Sky Survey (LOTSS) and present the results of this pipeline from the LOTSS preliminary data release region (10h45m - 15h30m right ascension, 45 - 57 degrees declination, 570 square degrees). We have produced a catalog of 92 polarized radio sources at 150 MHz at 4.3 arcminute resolution and 1 mJy rms sensitivity, which is the largest catalog of polarized sources at such low frequencies. We estimate a lower limit to the polarized source surface density at 150 MHz, with our resolution and sensitivity, of 1 source per 6.2 square degrees. We find that our Faraday depth measurements are in agreement with previous measurements and have significantly smaller errors. Most of our sources show significant depolarization compared to 1.4 GHz, but there is a small population of sources with low depolarization indicating that their polarized emission is highly localized in Faraday depth. We predict that an extension of this work to the full LOTSS data would detect at least 3400 polarized sources using the same methods, and probably considerably more with improved data processing.

The new generation of low-frequency radio telescopes, such as the Low Frequency Array (LOFAR: a Square Kilometre Array-low pathfinder), provides advancements in our capability of probing Galactic magnetism through low-frequency polarimetry. Maps of diffuse polarized radio emission and Faraday rotation can be used to infer properties of, and trace structure in, the magnetic fields in the ISM. However, to date very little of the sky has been probed at high angular and Faraday depth resolution. We observed a 5x5 degree region centred on the nearby galaxy IC342 using LOFAR in the frequency range 115-178 MHz at 4 arcmin resolution and performed Faraday tomography to detect foreground Galactic polarized synchrotron emission separated by Faraday depth (different amounts of Faraday rotation). Our Faraday depth cube shows rich polarized structure, with up to 30 K of polarized emission at 150 MHz. We detect two overlapping diffuse polarized features that are clearly separated in Faraday depth. Faraday-thick structures at such low frequencies would be too strongly depolarized to explain the observations and are therefore rejected. Only Faraday thin structures will not be strongly depolarized; producing such structures requires localized variations in the ratio of synchrotron emissivity to Faraday depth per unit distance, which can arise from several physical phenomena, such as a transition between regions of ionized and neutral gas. We conclude that the observed polarized emission is Faraday thin, and propose that the emission originates from two neutral clouds in the local ISM. We have modeled the Faraday rotation for this line of sight and estimated that the line of sight component of magnetic field of the local ISM for this direction varies between -0.86 and +0.12 uG. We propose that this may be a useful method for mapping magnetic fields within the local ISM.

Recent models for the large-scale Galactic magnetic fields in the literature have been largely constrained by synchrotron emission and Faraday rotation measures. We use three different but representative models to compare their predicted polarized synchrotron and dust emission with that measured by the Planck satellite. We first update these models to match the Planck synchrotron products using a common model for the cosmic-ray leptons. We discuss the impact on this analysis of the ongoing problems of component separation in the Planck microwave bands and of the uncertain cosmic-ray spectrum. In particular, the inferred degree of ordering in the magnetic fields is sensitive to these systematic uncertainties, and we further show the importance of considering the expected variations in the observables in addition to their mean morphology. We then compare the resulting simulated emission to the observed dust polarization and find that the dust predictions do not match the morphology in the Planck data but underpredict the dust polarization away from the plane. We modify one of the models to roughly match both observables at high latitudes by increasing the field ordering in the thin disc near the observer. Though this specific analysis is dependent on the component separation issues, we present the improved model as a proof of concept for how these studies can be advanced in future using complementary information from ongoing and planned observational projects.

To study the strength and structure of the magnetic field in the Galactic centre (GC) we measured Faraday rotation of the radio emission of pulsars which are seen towards the GC. Three of these pulsars have the largest rotation measures (RMs) observed in any Galactic object with the exception of Sgr A*. Their large dispersion measures, RMs and the large RM variation between these pulsars and other known objects in the GC implies that the pulsars lie in the GC and are not merely seen in projection towards the GC. The large RMs of these pulsars indicate large line-of-sight magnetic field components between ~ 16-33 microgauss; combined with recent model predictions for the strength of the magnetic field in the GC this implies that the large-scale magnetic field has a very small inclination angle with respect to the plane of the sky (~ 12 degrees). Foreground objects like the Radio Arc or possibly an ablated, ionized halo around the molecular cloud G0.11-0.11 could contribute to the large RMs of two of the pulsars. If these pulsars lie behind the Radio Arc or G0.11-0.11 then this proves that low-scattering corridors with lengths >~ 100 pc must exist in the GC. This also suggests that future, sensitive observations will be able to detect additional pulsars in the GC. Finally, we show that the GC component in our most accurate electron density model oversimplifies structure in the GC.

We study the statistical properties of interstellar dust polarization at high Galactic latitude, using the Stokes parameter Planck maps at 353 GHz. Our aim is to advance the understanding of the magnetized interstellar medium (ISM), and to provide a model of the polarized dust foreground for cosmic microwave background component-separation procedures. Focusing on the southern Galactic cap, we examine the statistical distributions of the polarization fraction ($p$) and angle ($\psi$) to characterize the ordered and turbulent components of the Galactic magnetic field (GMF) in the solar neighbourhood. We relate patterns at large angular scales in polarization to the orientation of the mean (ordered) GMF towards Galactic coordinates $(l_0,b_0)=(70^\circ \pm 5^\circ,24^\circ \pm 5^\circ)$. The histogram of $p$ shows a wide dispersion up to 25 %. The histogram of $\psi$ has a standard deviation of $12^\circ$ about the regular pattern expected from the ordered GMF. We use these histograms to build a phenomenological model of the turbulent component of the GMF, assuming a uniform effective polarization fraction ($p_0$) of dust emission. To model the Stokes parameters, we approximate the integration along the line of sight (LOS) as a sum over a set of $N$ independent polarization layers, in each of which the turbulent component of the GMF is obtained from Gaussian realizations of a power-law power spectrum. We are able to reproduce the observed $p$ and $\psi$ distributions using: a $p_0$ value of (26 $\pm$ 3)%; a ratio of 0.9 $\pm$ 0.1 between the strengths of the turbulent and mean components of the GMF; and a small value of $N$. We relate the polarization layers to the density structure and to the correlation length of the GMF along the LOS.

Galactic cosmic rays (CR) are particles presumably accelerated in supernova remnant shocks that propagate in the interstellar medium up to the densest parts of molecular clouds, losing energy and their ionisation efficiency because of the presence of magnetic fields and collisions with molecular hydrogen. Recent observations hint at high levels of ionisation and at the presence of synchrotron emission in protostellar systems, which leads to an apparent contradiction. We want to explain the origin of these CRs accelerated within young protostars as suggested by observations. Our modelling consists of a set of conditions that has to be satisfied in order to have an efficient CR acceleration through diffusive shock acceleration. We analyse three main acceleration sites, then we follow the propagation of these particles through the protostellar system up to the hot spot region. We find that jet shocks can be strong accelerators of CR protons, which can be boosted up to relativistic energies. Other promising acceleration sites are protostellar surfaces, where shocks caused by impacting material during the collapse phase are strong enough to accelerate CR protons. In contrast, accretion flow shocks are too weak to efficiently accelerate CRs. Though CR electrons are weakly accelerated, they can gain a strong boost to relativistic energies through re-acceleration in successive shocks. We suggest a mechanism able to accelerate both CR protons and electrons through the diffusive shock acceleration mechanism, which can be used to explain the high ionisation rate and the synchrotron emission observed towards protostellar sources. The existence of an internal source of energetic particles can have a strong and unforeseen impact on the ionisation of the protostellar disc, on the star and planet formation processes, and on the formation of pre-biotic molecules.

The quest for a B-mode imprint from primordial gravity waves on the polarization of the cosmic microwave background (CMB) requires the characterization of foreground polarization from Galactic dust. We present a statistical study of the filamentary structure of the 353 GHz Planck Stokes maps at high Galactic latitude, relevant to the study of dust emission as a polarized foreground to the CMB. We filter the intensity and polarization maps to isolate filaments in the range of angular scales where the power asymmetry between E-modes and B-modes is observed. Using the Smoothed Hessian Major Axis Filament Finder, we identify 259 filaments at high Galactic latitude, with lengths larger or equal to 2\deg\ (corresponding to 3.5 pc in length for a typical distance of 100 pc). These filaments show a preferred orientation parallel to the magnetic field projected onto the plane of the sky, derived from their polarization angles. We present mean maps of the filaments in Stokes I, Q, U, E, and B, computed by stacking individual images rotated to align the orientations of the filaments. Combining the stacked images and the histogram of relative orientations, we estimate the mean polarization fraction of the filaments to be 11 %. Furthermore, we show that the correlation between the filaments and the magnetic field orientations may account for the E and B asymmetry and the $C_{\ell}^{TE}/C_{\ell}^{EE}$ ratio, reported in the power spectra analysis of the Planck 353 GHz polarization maps. Future models of the dust foreground for CMB polarization studies will need to take into account the observed correlation between the dust polarization and the structure of interstellar matter.

The main signature of the interaction between cosmic rays and molecular clouds is the high ionisation degree. This decreases towards the densest parts of a cloud, where star formation is expected, because of energy losses and magnetic effects. However recent observations hint to high levels of ionisation in protostellar systems, therefore leading to an apparent contradiction that could be explained by the presence of energetic particles accelerated within young protostars. Our modelling consists of a set of conditions that has to be satisfied in order to have an efficient particle acceleration through the diffusive shock acceleration mechanism. We find that jet shocks can be strong accelerators of protons which can be boosted up to relativistic energies. Another possibly efficient acceleration site is located at protostellar surfaces, where shocks caused by impacting material during the collapse phase are strong enough to accelerate protons. Our results demonstrate the possibility of accelerating particles during the early phase of a proto-Solar-like system and can be used as an argument to support available observations. The existence of an internal source of energetic particles can have a strong and unforeseen impact on the star and planet formation process as well as on the formation of pre-biotic molecules.

Planck has mapped the polarized dust emission over the whole sky, making it possible to trace the Galactic magnetic field structure that pervades the interstellar medium (ISM). We combine polarization data from Planck with rotation measure (RM) observations towards a massive star-forming region, the Rosette Nebula in the Monoceros molecular cloud, to study its magnetic field structure and the impact of an expanding HII region on the morphology of the field. We derive an analytical solution for the magnetic field, assumed to evolve from an initially uniform configuration following the expansion of ionized gas and the formation of a shell of swept-up ISM. From the RM data we estimate a mean value of the line-of-sight component of the magnetic field of about 3microG (towards the observer) in the Rosette Nebula, for a uniform electron density of about 12cm-3. The dust shell that surrounds the Rosette HII region is clearly observed in the Planck intensity map at 353 GHz, with a polarization signal significantly different from that of the local background when considered as a whole. The Planck observations constrain the plane-of-the-sky orientation of the magnetic field in the Rosette's parent molecular cloud to be mostly aligned with the large-scale field along the Galactic plane. The Planck data are compared with the analytical model, which predicts the mean polarization properties of a spherical and uniform dust shell for a given orientation of the field. This comparison leads to an upper limit of about 45degr on the angle between the line of sight and the magnetic field in the Rosette complex, for an assumed intrinsic dust polarization fraction of 4%. This field direction can reproduce the RM values detected in the ionized region if the magnetic field strength in the Monoceros molecular cloud is in the range 6.5--9microG.

Within ten nearby (d < 450 pc) Gould Belt molecular clouds we evaluate statistically the relative orientation between the magnetic field projected on the plane of sky, inferred from the polarized thermal emission of Galactic dust observed by Planck at 353 GHz, and the gas column density structures, quantified by the gradient of the column density, $N_H$. The selected regions, covering several degrees in size, are analyzed at an effective angular resolution of 10' FWHM, thus sampling physical scales from 0.4 to 40 pc in the nearest cloud. The column densities in the selected regions range from $N_H \approx 10^{21}$ to $10^{23}$ cm$^{-2}$, and hence they correspond to the bulk of the molecular clouds. The relative orientation is evaluated pixel by pixel and analyzed in bins of column density using the novel statistical tool called "Histogram of Relative Orientations". Throughout this study, we assume that the polarized emission observed by Planck at 353 GHz is representative of the projected morphology of the magnetic field in each region, i.e., we assume a constant dust grain alignment efficiency, independent of the local environment. Within most clouds we find that the relative orientation changes progressively with increasing $N_H$, from preferentially parallel or having no preferred orientation to preferentially perpendicular. In simulations of magnetohydrodynamic turbulence in molecular clouds this trend in relative orientation is a signature of Alfv\'enic or sub-Alfv\'enic turbulence, implying that the magnetic field is significant for the gas dynamics at the scales probed by Planck. We compare the deduced magnetic field strength with estimates we obtain from other methods and discuss the implications of the Planck observations for the general picture of molecular cloud formation and evolution.

Recent Low Frequency Array (LOFAR) observations at 115-175 MHz of a field at medium Galactic latitudes (centered at the bright quasar 3C196) have shown striking filamentary structures in polarization that extend over more than 4 degrees across the sky. In addition, the Planck satellite has released full sky maps of the dust emission in polarization at 353GHz. The LOFAR data resolve Faraday structures along the line of sight, whereas the Planck dust polarization maps probe the orientation of the sky projected magnetic field component. Hence, no apparent correlation between the two is expected. Here we report a surprising, yet clear, correlation between the filamentary structures, detected with LOFAR, and the magnetic field orientation, probed by the Planck satellite. This finding points to a common, yet unclear, physical origin of the two measurements in this specific area in the sky. A number of follow-up multi- frequency studies are proposed to shed light on this unexpected finding.

Planck observations at 353GHz provide the first fully-sampled maps of the polarized dust emission towards interstellar filaments and their backgrounds. The polarization data provide insight on the structure of their magnetic field (B). We present the polarization maps of three nearby star forming filament of moderate column density (NH~10^22cm^-2): Musca, B211, and L1506. We use the spatial information to separate Stokes I, Q, and U of the filaments from those of their backgrounds, an essential step to measure the intrinsic polarization fraction (p) and angle (psi) of each emission component. We find that the polarization angles in the three filaments (psi_fil) are coherent along their lengths and not the same as in their backgrounds (psi_bg). The differences between psi_fil and psi_bg are 12deg, 6deg, and 54deg for Musca, B211, and L1506, respectively. These differences for Musca and L1506 are larger than the dispersions of psi, both along the filaments and in their backgrounds. The observed changes of psi are direct evidence for variations of the orientation of the plane of the sky (POS) projection of the B-field. As in previous studies, we find a decrease of several percent of p with NH. We show that the drop in p cannot be explained by random fluctuations of the orientation of B within the filaments because they are too small (sigma_psi<10deg). We recognize the degeneracy between dust alignment efficiency and the structure of B in causing variations in p, but we argue that the decrease of p from the backgrounds to the filaments results in part from depolarization associated with the 3D structure of B: both its orientation in the POS and with respect to the POS. We do not resolve the inner structure of the filaments, but at the smallest scales accessible with Planck (~0.2pc), the observed changes of psi and p hold information on the B-field structure within filaments.

The role of the magnetic field in the formation of the filamentary structures observed in the interstellar medium (ISM) is a debated topic. The Planck all-sky maps of linearly polarized emission from dust at 353GHz provide the required combination of imaging and statistics to study the correlation between the structures of the Galactic magnetic field and of interstellar matter, both in the diffuse ISM and in molecular clouds. The data reveal structures, or ridges, in the intensity map with counterparts in the Stokes Q and/or U maps. We focus on structures at intermediate and high Galactic latitudes with column density from $10^{20}$ to $10^{22}$ cm$^{-2}$. We measure the magnetic field orientation on the plane of the sky from the polarization data, and present an algorithm to estimate the orientation of the ridges from the dust intensity map. We use analytical models to account for projection effects. Comparing polarization angles on and off the structures, we estimate the mean ratio between the strengths of the turbulent and mean components of the magnetic field to be between 0.6 and 1.0, with a preferred value of 0.8. We find that the ridges are preferentially aligned with the magnetic field measured on the structures. This trend becomes more striking for increasing polarization fraction and decreasing column density. We interpret the increase of alignment with polarization fraction as a consequence of projections effects. The decrease of alignment for high column density is not due to a loss of correlation between the structures and the geometry of the magnetic field. In molecular complexes, we observe structures perpendicular to the magnetic field, which cannot be accounted for by projection effects. We discuss our results in the context of models and MHD simulations, which describe the formation of structures in the magnetized ISM.

Magnetic fields are an important ingredient of the interstellar medium (ISM). Besides their importance for star formation, they govern the transport of cosmic rays, relevant to the launch and regulation of galactic outflows and winds, which in turn are pivotal in shaping the structure of halo magnetic fields. Mapping the small-scale structure of interstellar magnetic fields in many nearby galaxies is crucial to understand the interaction between gas and magnetic fields, in particular how gas flows are affected. Elucidation of the magnetic role in, e.g., triggering star formation, forming and stabilising spiral arms, driving outflows, gas heating by reconnection and magnetising the intergalactic medium has the potential to revolutionise our physical picture of the ISM and galaxy evolution in general. Radio polarisation observations in the very nearest galaxies at high frequencies (>= 3 GHz) and with high spatial resolution (<= 5") hold the key here. The galaxy survey with SKA1 that we propose will also be a major step to understand the galactic dynamo, which is important for models of galaxy evolution and for astrophysical magnetohydrodynamics in general. Field amplification by turbulent gas motions, which is crucial for efficient dynamo action, has been investigated so far only in simulations, while compelling evidence of turbulent fields from observations is still lacking.

Polarized emission observed by Planck HFI at 353 GHz towards a sample of nearby fields is presented, focusing on the statistics of polarization fractions $p$ and angles $\psi$. The polarization fractions and column densities in these nearby fields are representative of the range of values obtained over the whole sky. We find that: (i) the largest polarization fractions are reached in the most diffuse fields; (ii) the maximum polarization fraction $p_\mathrm{max}$ decreases with column density $N_\mathrm{H}$ in the more opaque fields with $N_\mathrm{H} > 10^{21}\,\mathrm{cm}^{-2}$; and (iii) the polarization fraction along a given line of sight is correlated with the local spatial coherence of the polarization angle. These observations are compared to polarized emission maps computed in simulations of anisotropic magnetohydrodynamical (MHD) turbulence in which we assume a uniform intrinsic polarization fraction of the dust grains. We find that an estimate of this parameter may be recovered from the maximum polarization fraction $p_\mathrm{max}$ in diffuse regions where the magnetic field is ordered on large scales and perpendicular to the line of sight. This emphasizes the impact of anisotropies of the magnetic field on the emerging polarization signal. The decrease of the polarization fraction with column density in nearby molecular clouds is well reproduced in the simulations, indicating that it is essentially due to the turbulent structure of the magnetic field: an accumulation of variously polarized structures along the line of sight leads to such an anti-correlation. In the simulations, polarization fractions are also found to anti-correlate with the angle dispersion function $\mathcal{S}$. [abridged]

This paper presents the large-scale polarized sky as seen by Planck HFI at 353 GHz, which is the most sensitive Planck channel for dust polarization. We construct and analyse large-scale maps of dust polarization fraction and polarization direction, while taking account of noise bias and possible systematic effects. We find that the maximum observed dust polarization fraction is high (pmax > 18%), in particular in some of the intermediate dust column density (AV < 1mag) regions. There is a systematic decrease in the dust polarization fraction with increasing dust column density, and we interpret the features of this correlation in light of both radiative grain alignment predictions and fluctuations in the magnetic field orientation. We also characterize the spatial structure of the polarization angle using the angle dispersion function and find that, in nearby fields at intermediate latitudes, the polarization angle is ordered over extended areas that are separated by filamentary structures, which appear as interfaces where the magnetic field sky projection rotates abruptly without apparent variations in the dust column density. The polarization fraction is found to be anti-correlated with the dispersion of the polarization angle, implying that the variations are likely due to fluctuations in the 3D magnetic field orientation along the line of sight sampling the diffuse interstellar medium.We also compare the dust emission with the polarized synchrotron emission measured with the Planck LFI, with low-frequency radio data, and with Faraday rotation measurements of extragalactic sources. The two polarized components are globally similar in structure along the plane and notably in the Fan and North Polar Spur regions. A detailed comparison of these three tracers shows, however, that dust and cosmic rays generally sample different parts of the line of sight and confirms that much of the variation observed in the Planck data is due to the 3D structure of the magnetic field.

PRISM (Polarized Radiation Imaging and Spectroscopy Mission) was proposed to ESA in May 2013 as a large-class mission for investigating within the framework of the ESA Cosmic Vision program a set of important scientific questions that require high resolution, high sensitivity, full-sky observations of the sky emission at wavelengths ranging from millimeter-wave to the far-infrared. PRISM's main objective is to explore the distant universe, probing cosmic history from very early times until now as well as the structures, distribution of matter, and velocity flows throughout our Hubble volume. PRISM will survey the full sky in a large number of frequency bands in both intensity and polarization and will measure the absolute spectrum of sky emission more than three orders of magnitude better than COBE FIRAS. The aim of this Extended White Paper is to provide a more detailed overview of the highlights of the new science that will be made possible by PRISM

We seek to understand the propagation mechanisms of positrons in the interstellar medium (ISM). This understanding is a key to determine whether the spatial distribution of the annihilation emission observed in our Galaxy reflects the spatial distribution of positron sources and, therefore, makes it possible to place constraints on the origin of positrons. We review the different processes that are likely to affect the transport of positrons in the ISM. These processes fall into three broad categories: scattering off magnetohydrodynamic waves, collisions with particles of the interstellar gas and advection with large-scale fluid motions. We assess the efficiency of each process and describe its impact on the propagation of positrons. We also develop a model of positron propagation, based on Monte-Carlo simulations, which enable us to estimate the distances traveled by positrons in the different phases of the ISM. We find that low-energy (< 10 MeV) positrons generally have negligible interactions with magnetohydrodynamic waves, insofar as these waves are heavily damped. Positron propagation is mainly controlled by collisions with gas particles. Under these circumstances, positrons can travel huge distances (up to 30 kpc/n_H for 1 MeV positrons) along magnetic field lines before annihilating.