Probing the cold magnetized Universe with SPICA-POL (B-BOP)

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.


Preface
The following set of articles describe in detail the science goals of the future Space Infrared telescope for Cosmology and Astrophysics (SPICA). The SPICA satellite will employ a 2.5-m telescope, actively cooled to below 8 K, and a suite of mid-to far-infrared spectrometers and photometric cameras, equipped with state-of-theart detectors. In particular, the SPICA Far Infrared Instrument (SAFARI) will be a grating spectrograph with low (R = 300) and medium (R = 3000-11000) resolution observing modes instantaneously covering the 35-230 µm wavelength range. The SPICA Mid-Infrared Instrument (SMI) will have three operating modes: a large field of view (12 ×10 ) low-resolution 17-36 µm spectroscopic (R = 50-120) and photometric camera at 34 µm, a medium resolution (R = 2000) grating spectrometer covering wavelengths of 18-36 µm and a high-resolution echelle module (R = 28000) for the 12-18 µm domain. A large-field-of-view (160 ×160 ) 1 , three-channel (100, 200, and 350 µm) polarimetric camera (B-BOP 2 ) will also be part of the instrument complement. These articles will focus on some of the major scientific questions that the SPICA mission aims to address; more details about the mission and instruments can be found in 1 Some other SPICA papers refer to this field of view as 80 ×80 , but it is 160 ×160 according to the latest design.

INTRODUCTION: SPICA AND THE NATURE OF COSMIC MAGNETISM
Alongside gravity, magnetic fields play a key role in the formation and evolution of a wide range of structures in the Universe, from galaxies to stars and planets. They simultaneously are an actor, an outcome, and a tracer of cosmic evolution. These three facets of cosmic magnetism are intertwined and must be thought of together. On one hand, the role magnetic fields play in the formation of stars and galaxies results from and traces their interplay with gas dynamics. On the other hand, turbulence is central to the dynamo processes that initially amplified cosmic magnetic fields and have since maintained their strength in galaxies across time (Brandenburg & Subramanian, 2005). A transfer from gas kinetic to magnetic energy inevitably takes place in turbulent cosmic flows, while magnetic fields act on gas dynamics through the Lorentz force. These physical couplings relate cosmic magnetism to structure formation in the Universe across time and scales, and make the observation of magnetic fields a tracer of cosmic evolution, which is today yet to be disclosed. Improving our observational understanding of cosmic magnetism on a broad range of physical scales is thus at the heart of the "Origins" big question and is an integral part of one of ESA's four Grand Science Themes ("Cosmic Radiation and Magnetism") as defined by the ESA High-level Science Policy Advisory Committee (HISPAC) in 2013. As often in Astrophysics, our understanding of the Universe is rooted in observations of the very local universe: the Milky Way and nearby galaxies. In the interstellar medium (ISM) of these galaxies, the magnetic energy is observed to be in rough equipartition with the kinetic (e.g. turbulent), radiative, and cosmic ray energies, all on the order of ∼ 1 eV cm −3 , suggesting that magnetic fields are a key player in the dynamics of the ISM (e.g. Draine, 2011). Their exact role in the formation of molecular clouds and stellar systems is not well understood, however, and remains highly debated (e.g. Crutcher, 2012). Interstellar magnetic fields also hold the key for making headway on other main issues in Astrophysics, including the dynamics and energetics of the multiphase ISM, the acceleration and propagation of cosmic rays, and the physics of stellar and back-hole feedback. Altogether, a broad range of science topics call for progress in our understanding of interstellar magnetic fields, which in turn motivates ambitious efforts to obtain relevant data (cf. Boulanger et al., 2018).
Observations of Galactic polarization are a highlight and a lasting legacy of the Planck space mission. Spectacular images combining the intensity of dust emission with the texture derived from polarization data have received world-wide attention and have become part of the general scientific culture (Planck 2015(Planck res. I, 2016. Beyond their popular impact, the Planck polarization maps represented an immense step forward for Galactic astrophysics (Planck 2018res. XII, 2019. Planck has paved the way for statistical studies of the structure of the Galactic magnetic field and its coupling with interstellar matter and turbulence, in the diffuse ISM and star-forming molecular clouds. SPICA, the Space Infrared Telescope for Cosmology and Astrophysics proposed to ESA as an M5 mission concept (Roelfsema et al., 2018), provides one of the best opportunities to take the next big leap forward and gain fundamental insight into the role of magnetic fields in structure formation in the cold Universe, thanks to the unprecedented sensitivity, angular resolution, and dynamic range of its far-infrared (far-IR) imaging polarimeter, B-BOP 2 (previously called SPICA-POL, for "SPICA polarimeter"). The baseline B-BOP instrument will allow simultaneous imaging observations in three bands, 100 µm, 200 µm, and 350 µm, with an individual pixel NEP < 3 × 10 −18 W Hz −1/2 , over an instantaneous field of view of ∼ 2.7 × 2.7 at resolutions of 9 , 18 , and 32 , respectively (Rodriguez et al., 2018). Benefiting from a 2.5-m space telescope cooled to < 8 K, B-BOP will be two to three orders of magnitude more sensitive than current or planned far-IR/submillimeter polarimeters (see § 2.4 below) and will produce far-IR dust polarization images at a factor 20-30 higher resolution than the Planck satellite. It will provide wide-field 100-350 µm polarimetric images in Stokes Q and U of comparable quality (in terms of resolution, signal-to-noise ratio, and both intensity and spatial dynamic ranges) to Herschel images in Stokes I.
The present paper gives an overview of the main science drivers for the B-BOP polarimeter and is complementary to the papers by, e.g., Spinoglio et al. (2017) and van der Tak et al. (2018) which discuss the science questions addressed by the other two instruments of SPICA, SMI (Kaneda et al., 2016) and SAFARI (Roelfsema et al., 2014), mainly through highly sensitive spectroscopy. The outline is as follows: Section 2 describes the prime science driver for B-BOP, namely high dynamic range polarimetric mapping of Galactic filamentary structures to unravel the role of magnetic fields in the star formation process. Section 3 introduces the contribution of B-BOP to the statistical characterization of magnetized interstellar turbulence. Section 4 and Section 5 emphasize the importance of B-BOP polarization observations for our understanding of the physics of protostellar dense cores and high-mass star protoclusters, respectively. Section 6 discusses dust polarization observations of galaxies, focusing mainly on nearby galaxies. Section 7 describes how multi-wavelength polarimetry with B-BOP can constrain dust models and the physics of dust grain alignment. Finally, Sections 8, 9, and 10 discuss three topics which, although not among the main ≥ 20 /sec ≥ 20 /sec ≥ 20 /sec † Surface brightness level in I to map Q, U at 5σ over 1 deg 2 in 10hr assuming 5% fractional polarization.
drivers of the B-BOP instrument, will significantly benefit from B-BOP observations, namely the study of the origin of cosmic rays and of their interaction with molecular clouds (Sect. 8), the detection of polarized far-IR dust emission from protoplanetary disks, thereby tracing magnetic fields in the inner layers of the disks (Sect. 9), and the (non-polarimetric) monitoring of protostars in the far-IR, i.e., close to the peak of their spectral energy distributions (SEDs), to provide direct constraints on the process of episodic protostellar accretion (Sect. 10). Section 11 concludes the paper.

MAGNETIC FIELDS AND STAR FORMATION IN FILAMENTARY CLOUDS
Understanding how stars form in the cold ISM of galaxies is central in Astrophysics. Star formation is both one of the main factors that drive the evolution of galaxies on global scales and the process that sets the physical conditions for planet formation on local scales. Star formation is also a complex, multi-scale process, involving a subtle interplay between gravity, turbulence, magnetic fields, feedback mechanisms. As a consequence, and despite recent progress, the basic questions of what regulates star formation in galaxies and what determines the mass distribution of forming stars (i,e. the stellar initial mass function or IMF) remain two of the most debated problems in Astronomy. Today, a popular school of thought for understanding star formation and these two big questions is the gravo-turbulent paradigm (e.g. Mac Low & Klessen, 2004;McKee & Ostriker, 2007;Padoan et al., 2014), whereby magnetized supersonic turbulence creates structure and seeds in interstellar clouds, which subsequently grow and collapse under the primary influence of gravity. A variation on this scenario is that of dominant magnetic fields in cloud envelopes, and a turbulence-enhanced ambipolar diffusion leading to gravity-dominated subregions (e.g., Li & Nakamura, 2004;Kudoh & Basu, 2008). Moreover, while the global rate of star formation in galaxies and the positions of galaxies in the Schmidt-Kennicutt diagram (e.g. Kennicutt & Evans, 2012) are likely controlled by macroscopic phenomena such as cosmic accretion, large-scale feedback, and large-scale turbulence (Sánchez Almeida et al., 2014), there is some evidence that the star formation efficiency in the dense molecular gas of galaxies is nearly universal 3 (e.g. Gao & Solomon, 2004;Lada et al., 2012) and primarily governed by the physics of filamentary cloud fragmentation on much smaller scales (e.g. André et al., 2014;Padoan et al., 2014). As argued in § 2.2 and § 2.3 below, magnetic fields are likely a key element of the physics behind the formation and fragmentation of filamentary structures in interstellar clouds.
Often ignored, strong, organized magnetic fields, in rough equipartition with the turbulent and cosmic ray energy densities, have been detected in the ISM of a large number of galaxies out to z = 2 (e.g. Bernet et al., 2008). Recent cosmological magneto-hydrodynamic (MHD) simulations of structure formation in the Universe suggest that magnetic-field strengths comparable to those measured in nearby galaxies ( < ∼ 10 µG) can be quickly built up in high-redshift galaxies (in << 1 Gyr), through the dynamo amplification of initially weak seed fields (e.g. Rieder & Teyssier, 2017;Marinacci et al., 2018). Magnetic fields are therefore expected to play a dynamically important role in the formation of giant molecular clouds (GMCs) on kpc scales within galaxies (e.g. Inoue & Inutsuka, 2012) and in the formation of filamentary structures leading to individual star formation on ∼1-10 pc scales within GMCs (e.g. Inutsuka et al., 2015;Inoue et al., 2018, see § 2.3 below). On dense core (≤ 0.1 pc) scales, the magnetic field and angular momentum of most protostellar systems are likely inherited from the processes of filament formation and fragmentation (cf. Misugi et al., in prep.). On even smaller (< 0.01 pc or < 2000 au) scales, magnetic fields are essential to solve the angular momentum problem of star formation, generate protostellar outflows, and control the formation of protoplanetary disks (e.g. Pudritz et al., 2007;Machida et al., 2008;Li et al., 2014).
In this context, B-BOP will be a unique tool for characterizing the morphology of magnetic fields on scales ranging from ∼ 0.01 pc to ∼ 1 kpc in Milky Way like galaxies. In particular, a key science driver for B-BOP is to clarify the role of magnetic fields in shaping the rich web of filamentary structures pervading the cold ISM, from the low-density striations seen in HI clouds and the outskirts of CO clouds (e.g. Clark et al., 2014;Kalberla et al., 2016;Goldsmith et al., 2008) to the denser molecular filaments within which most prestellar cores and protostars are forming according to Herschel results (see Fig. 1 and § 2.2 below).

Dust grain alignment
Polarization of background starlight from dichroic extinction produced by intervening interstellar dust has been known since the late 1940s (Hall, 1949;Hiltner, 1949). The analysis of the extinction data in polarization, in particular its variation with wavelength in the visible to near-UV, has allowed major discoveries regarding dust properties, in particular regarding the size distribution of dust. Like the first large-scale total intensity mapping in the far-IR that was provided by the Infrared Astronomical Satellite (IRAS) satellite data, extensive studies of polarized far-IR emission today bring the prospect of a new revolution in our understanding of dust physics. This endeavor includes pioneering observations with ground-based, balloon-borne, and space-borne facilities such as the very recent all-sky observations by the Planck satellite at 850 µm and beyond (e.g. Planck 2018 res. XII, 2019). However, polarimetric imaging of polarized dust continuum emission is still in its infancy and amazing improvements are expected in the next decades from instruments such as the Atacama Large Millimeter Array (ALMA) in the submillimeter and B-BOP in the far-IR. The initial discovery that starlight extinguished by intervening dust is polarized led to the conclusion that dust grains must be somewhat elongated and globally aligned in space in order to produce the observed polarized extinction. While the elongation of dust grains was not unexpected, coherent grain alignment over large spatial scales has been more difficult to explain. A very important constraint has come from recent measurements in emission with, e.g., the Archeops balloon-borne experiment (Benoît et al., 2004) and the Planck satellite (Planck int. res. XIX, 2015) which indicated that the polarization degree of dust emission can be as high as 20% in some regions of the diffuse ISM in the solar neighborhood. This requires more efficient dust alignment processes than previously anticipated (Planck 2018 res. XII, 2019).
The most widely accepted dust grain alignment theo-ries, already alluded to by Hiltner (1949), propose that alignment is with respect to the magnetic field that pervades the ISM. Rapidly spinning grains will naturally align their angular momentum with the magnetic field direction (Purcell, 1979;Lazarian & Draine, 1999), but the mechanism leading to such rapid spin remains a mystery. The formation of molecular hydrogen at the surface of dust grains could provide the required momentum (Purcell, 1979). Today's leading grain alignment theory is Radiative Alignment Torques (RATs) (Dolginov & Mitrofanov, 1976;Draine & Weingartner, 1996;Lazarian & Hoang, 2007;Hoang & Lazarian, 2016, and references therein), where supra-thermal spinup of irregularly-shaped dust grains results from their irradiation by an anisotropic radiation field (a process experimentally confirmed, see Abbas et al., 2004).

Probing magnetic fields with imaging polarimetry
In the conventional picture that the minor axis of elongated dust grains is aligned with the local direction of the magnetic field, mapping observations of linearlypolarized continuum emission at far-IR and submillimeter wavelengths are a powerful tool to measure the morphology and structure of magnetic field lines in starforming clouds and dense cores (cf. Matthews et al., 2009;Crutcher et al., 2004;Crutcher, 2012). A key advantage of this technique is that it images the structure of magnetic fields through an emission process that traces the mass of cold interstellar matter, i.e., the reservoir of gas directly involved in star formation. Indirect estimates of the plane-of-sky magnetic field strength B P OS can also be obtained using the Davis-Chandrasekhar-Fermi method (Davis, 1951;Chandrasekhar & Fermi, 1953): B P OS = α corr √ 4πρ δV /δΦ, where ρ is the gas density (which can be estimated to reasonable accuracy from Herschel column density maps, especially in the case of resolved filaments and cores -cf. Palmeirim et al., 2013;Roy et al., 2014), δV is the one-dimensional velocity dispersion (which can be estimated from line observations in an appropriate tracer such as N 2 H + for star-forming filaments and dense cores-e.g. André et al., 2007;Tafalla & Hacar, 2015), δΦ is the dispersion in polarization position angles directly measured in a dust polarization map, and α corr ≈ 0.5 is a correction factor obtained through numerical simulations (cf. Ostriker et al., 2001). Large-scale maps that resolve the above quantities over a large dynamic range of densities can be used to estimate the mass-to-flux ratio in different parts of a molecular cloud. This can test the idea that cloud envelopes may be magnetically supported and have a subcritical mass-to-flux ratio (Mouschovias & Ciolek, 1999;Shu et al., 1999). Recent applications of the Davis-Chandrasekhar-Fermi method using SCUBA2-POL 850 µm data taken as part of the BISTRO survey (Ward-Thompson et al., 2017) toward dusty molecular clumps in the Orion and Ophiuchus clouds are presented in Pattle et al. (2017), Kwon et al. (2018), and Soam et al. (2018). Refined estimates of both the mean and the turbulent component of B P OS can be derived from an analysis of the second-order angular structure function (or angular dispersion function) of observed polarization position angles < ∆Φ 2 (l) Houde et al., 2009). Alternatively, in localized regions where gravity dominates over MHD turbulence, the polarization-intensity gradient method can be used to obtain maps of the local magnetic field strength from maps of the misalignment angle δ between the local magnetic field (estimated from observed polarization position angles) and the local column density gradient (estimated from maps of total dust emission). Indeed, such δ maps provide information on the local ratio between the magnetic field tension force and the gravitational force (Koch et al., 2012(Koch et al., , 2014. Additionally, the paradigm of Alfvénic turbulence can be tested in dense regions where gravity dominates, in which the observed angular dispersion ∆Φ is expected to decrease in amplitude toward the center of dense cores where δV also decreases (Auddy et al., 2019).
Because the typical degree of polarized dust continuum emission is low (∼ 2%-5% -e.g. Matthews et al., 2009) and the range of relevant column densities spans three orders of magnitude from equivalent visual extinctions 4 A V ∼ 0.1 in the atomic medium to A V > 100 in the densest molecular filaments/cores, a systematic dust polarization study of the rich filamentary networks pervading nearby interstellar clouds and their connection to star formation requires a large improvement in sensitivity, mapping speed, and dynamic range over existing far-IR/submillimeter polarimeters. A big improvement in polarimetric mapping speed is also needed for statistical reasons. As only the plane-of-sky component of the magnetic field is directly accessible to dust continuum polarimetry, a large number of systems must be imaged in various Galactic environments before physically meaningful conclusions can be drawn statistically on the role of magnetic fields. As shown in § 2.4 below, the required step forward in performance can be uniquely provided by a large, cryogenically cooled space-borne telescope such as SPICA, which can do in far-IR polarimetric imaging what Herschel achieved in total-power continuum imaging.

Insights from Herschel and Planck: A filamentary paradigm for star formation?
The Herschel mission has led to spectacular advances in our knowledge of the texture of the cold ISM and its link with star formation. While interstellar clouds have been known to be filamentary for a long time (e.g. Schneider & Elmegreen, 1979;Bally et al., 1987;Myers, 2009, and references therein), Herschel imaging surveys have established the ubiquity of filaments on almost all length scales (∼ 0.5 pc to ∼ 100 pc) in the molecular clouds of the Galaxy and shown that this filamentary structure likely plays a key role in the star formation process (e.g. André et al., 2010;Henning et al., 2010;Molinari et al., 2010;Hill et al., 2011;Schisano et al., 2014;Wang et al., 2015). The interstellar filamentary structures detected with Herschel span broad ranges in length, central column density, and mass per unit length (e.g. Schisano et al., 2014;Arzoumanian et al., 2019). In contrast, detailed analysis of the radial column density profiles indicates that, at least in the nearby molecular clouds of the Gould Belt, Herschel filaments are characterized by a narrow distribution of inner widths with a typical value of ∼ 0.1 pc and a dispersion of less than a factor of 2, when the data are averaged over the filament crests (Arzoumanian et al., 2011. Independent studies of filament widths in nearby clouds have generally confirmed this result when using submillimeter continuum data (e.g. Koch & Rosolowsky, 2015;Salji et al., 2015;Rivera-Ingraham et al., 2016), even if factor of ∼ 2-4 variations around the mean inner width of ∼ 0.1 pc have been found along the main axis of a given filament (e.g. Juvela et al., 2012;Ysard et al., 2013). Measurements of filament widths obtained in molecular line tracers (e.g., Pineda et al., 2011;Fernández-López et al., 2014;Panopoulou et al., 2014;Hacar et al., 2018) have been less consistent with the Herschel dust continuum results of Arzoumanian et al. (2011Arzoumanian et al. ( , 2018, but this can be attributed to the lower dynamic range achieved by observations in any given molecular line tracer. Panopoulou et al. (2017) pointed out an apparent contradiction between the existence of a characteristic filament width and the essentially scale-free nature of the power spectrum of interstellar cloud images (well described by a single power law from ∼ 0.01 pc to ∼ 50 pc - Miville-Deschênes et al., 2010, but Roy et al. (2019) showed that there is no contradiction given the only modest area filling factors ( < ∼ 10%) and column density contrasts (≤ 100% in most cases) derived by Arzoumanian et al. (2019) for the filaments seen in Herschel images. While further high-resolution submillimeter continuum studies would be required to investigate whether the same result holds beyond the Gould Belt, the median inner width of ∼ 0.1 pc measured with Herschel appears to reflect the presence of a true common scale in the filamentary structure of nearby interstellar clouds. If confirmed, this result may have far-reaching consequences as it introduces a characteristic scale in a system generally thought to be chaotic and turbulent (i.e. largely scale-free -cf. Guszejnov et al., 2018). It may thus present a severe challenge in any attempt to interpret all ISM observations Multi-resolution column density map of the Taurus molecular cloud as derived from a combination of high-resolution (18 -36 HPBW) observations from the Herschel Gould Belt survey and low-resolution (5 HPBW) Planck data. The superimposed "drapery" pattern traces the magnetic-field orientation projected on the plane of the sky, as inferred from Planck polarization data at 850 µm (Planck int. res. XXXV, 2016). (b) Herschel/SPIRE 250 µm dust continuum image of the B211/B213 filament in the Taurus cloud (Palmeirim et al., 2013;Marsh et al., 2016). The superimposed blue dashed curves trace the magnetic field orientation projected on the plane of the sky, as inferred from P lanck dust polarization data at 850 µm (Planck int. res. XXXV, 2016). Note the presence of faint striations oriented roughly perpendicular to the main filament and parallel to the plane-of-sky magnetic field. (c) IRAM/NIKA1 1.2 mm dust continuum image of the central part of the Herschel field shown in (b) (effective HPBW resolution of 20 ), showing a chain of at least four equally-spaced dense cores along the B211/B213 filament (from Bracco et al., 2017). B-BOP can image the magnetic field lines at a factor 30 better resolution than Planck over the entire Taurus cloud (cf. panel a), probing scales from ∼ 0.01 pc to > 10 pc.
in terms of scale-free processes.
Another major result from Herschel studies of nearby clouds is that most (> 75%) prestellar cores and protostars are found to lie in dense, "supercritical" filaments above a critical threshold ∼ 16 M /pc in mass per unit length, equivalent to a critical threshold ∼ 160 M /pc 2 (A V ∼ 8) in column density or n H2 ∼ 2 × 10 4 cm −3 in volume density (André et al., 2010;Könyves et al., 2015;Marsh et al., 2016). A similar column density threshold for the formation of prestellar cores (at A V ∼ 5-10) had been suggested earlier based on ground-based millimeter and submillimeter studies (e.g. Onishi et al., 1998;Johnstone et al., 2004;Kirk et al., 2006), but without clear connection to filaments. Interestingly, a comparable threshold in extinction (at A V ∼ 8) has also been observed in the spatial distribution of young stellar objects (YSOs) with Spitzer (e.g. Heiderman et al., 2010;Lada et al., 2010;Evans et al., 2014).
Overall, the Herschel results support a filamentary paradigm for star formation in two main steps (e.g. André et al., 2014;Inutsuka et al., 2015): First, multiple large-scale compressions of interstellar material in supersonic turbulent MHD flows generate a cobweb of ∼ 0.1 pc-wide filaments in the cold ISM; second, the densest filaments fragment into prestellar cores (and subsequently protostars) by gravitational instability above the critical mass per unit length M line,crit = 2 c 2 s /G of nearly isothermal, cylinder-like filaments (see Fig. 1), where c s is the sound speed and G the gravitational constant. This paradigm differs from the classical gravo-turbulent picture in that it relies on the unique features of filamentary geometry, such as the existence of a critical line mass for nearly isothermal filaments (e.g. Inutsuka & Miyama, 1997, and references therein). The validity and details of the filamentary paradigm are strongly debated, however, and many issues remain open. For instance, according to some numerical simulations, the above two steps may not occur consecutively but simultaneously, in the sense that both filamentary structures and dense cores may grow in mass at the same time (e.g. Gómez & Vázquez-Semadeni, 2014;Chen & Ostriker, 2015). The physical origin of the typical ∼ 0.1 pc inner width of molecular filaments is also poorly understood and remains a challenge for numerical models (e.g. Padoan et al., 2001;Hennebelle, 2013;Smith et al., 2014;Ntormousi et al., 2016). Auddy et al. (2016) point out that magnetized filaments may actually be ribbon-like and quasi-equilibrium structures supported by the magnetic field, and therefore not have cylindrical symmetry. Regardless of any particular scenario, there is nevertheless little doubt after Herschel results that dense molecular filaments represent an integral part of the initial conditions of the bulk of star formation in our Galaxy.
As molecular filaments are known to be present in the Large Magellanic Cloud (LMC -Fukui et al., 2015), the proposed filamentary paradigm may have implications on galaxy-wide scales. Assuming that all filaments have similar inner widths, it has been argued that they may help to regulate the star formation efficiency in dense  Palmeirim et al. (2013) using the algorithm getfilaments (Men'shchikov, 2013). In this view, all transverse angular scales larger than 72 (or ∼ 0.05 pc) were filtered out to enhance the contrast of the small-scale structure. The color scale is in MJy/sr at 250 µm. The colored curves display the velocity-coherent fibers independently identified by Hacar et al. (2013) using N 2 H + /C 18 O observations. (b) MHD simulation of a collapsing/accreting filament performed by E. Ntormousi & P. Hennebelle with the adaptive mesh refinement (AMR) code RAMSES. Line-of-sight velocities (in km/s) after one free-fall time (∼ 0.9 Myr) are coded by colors. For clarity, only the dense gas with 10 4 cm −3 < n H 2 < 10 5 cm −3 is shown. Note the braid-like velocity structure and the morphological similarity with the fiber-like pattern seen in the B211/B213 observations on the left. Thanks to its high resolution and dynamic range, B-BOP can probe, for the first time, the geometry of the magnetic field within the dense system of fibers and the connection with the low-density striations in the ambient cloud. molecular gas (André et al., 2014), and that they may be responsible for a quasi-universal star formation law in the dense molecular ISM of galaxies (cf. Lada et al., 2012;Shimajiri et al., 2017), with possible variations in extreme environments such as the CMZ . In parallel, the Planck mission has led to major advances in our knowledge of the geometry of the magnetic field on large scales in the Galactic ISM. The first allsky maps of dust polarization provided by Planck at 850 µm have revealed a very organized magnetic field structure on > ∼ 1-10 pc scales in Galactic interstellar clouds (Planck int. res. XXXV, 2016, see Fig. 1a). The large-scale magnetic field tends to be aligned with low-density filamentary structures with subcritical line masses such as striations (see Fig. 1b) and perpendicular to dense star-forming filaments with supercritical line masses (Planck int. res. XXXII, 2016; Planck int. res. XXXV, 2016, see Figs. 1b & 1c), a trend also seen in optical and near-IR polarization observations (Chapman et al., 2011;Palmeirim et al., 2013;Panopoulou et al., 2016;Soler et al., 2016). There is also a hint from Planck polarization observations of the nearest clouds that the direction of the magnetic field may change within dense filaments from nearly perpendicular in the ambient cloud to more parallel in the filament interior (cf. Planck int. res. XXXIII, 2016). These findings suggest that magnetic fields are dynamically important and play a key role in the formation and evolution of filamentary structures in interstellar clouds, supporting the view that dense molecular filaments form by accumulation of interstellar matter along field lines.
The low resolution of Planck polarization data (10 at best or 0.4 pc in nearby clouds) is however insufficient to probe the organization of field lines in the ∼ 0.1 pc interior of filaments, corresponding both to the characteristic transverse scale of filaments (Arzoumanian et al., 2011 and to the scale at which fragmentation into prestellar cores occurs (cf. Tafalla & Hacar, 2015). Consequently, the geometry of the magnetic field within interstellar filaments and its effects on fragmentation and star formation are essentially unknown today.

Investigating the role of magnetic fields in the formation and evolution of molecular filaments with B-BOP
Improving our understanding of the physics and detailed properties of molecular filaments is of paramount importance as the latter are representative of the initial conditions of star formation in molecular clouds and GMCs 5 (see § 2.2 above). In particular, investigating how dense, "supercritical" molecular filaments can maintain a roughly constant ∼0.1 pc inner width and fragment into prestellar cores instead of collapsing radially to spindles is crucial to understanding star formation. The topology of magnetic field lines may be one of the key elements here. For instance, a longitudinal magnetic field can support a filament against radial collapse but not against fragmentation along its main axis, while a perpendicular magnetic field works against fragmentation and increases the critical mass per unit length but cannot prevent the radial collapse of a supercritical filament (e.g. Tomisaka, 2014;Hanawa et al., 2017). The actual topology of the field within molecular filaments is likely more complex and may be a combination of these two extreme configurations. One plausible evolutionary scenario, consistent with existing observations, is that star-forming filaments accrete ambient cloud material along field lines through a network of magnetically-dominated striations (e.g. Palmeirim et al., 2013;Cox et al., 2016;Shimajiri et al., 2019, see also Figs. 1b & 2a). Accretion-driven MHD waves may then generate a system of velocity-coherent fibers within dense filaments (Hacar et al., 2013(Hacar et al., , 2018Arzoumanian et al., 2013;Hennebelle & André, 2013, cf. Fig. 2) and the corresponding organization of magnetic field lines may play a central role in accounting for the roughly constant ∼ 0.1 pc inner width of starforming filaments as measured in Herschel observations (cf. § 2.2). Constraining this process further is key to understanding star formation itself, since filaments with supercritical masses per unit length would otherwise undergo rapid radial contraction with time, effectively preventing fragmentation into prestellar cores and the formation of protostars (e.g. Inutsuka & Miyama, 1997). Information on the geometry of magnetic field lines within star-forming filaments at A V > 8 is thus crucially needed, and can be obtained through 200-350 µm dust polarimetric imaging at high angular resolution with B-BOP. Large area coverage and both high angular resolution and high spatial dynamic range are needed to resolve the 0.1 pc scale by a factor ∼ 3-10 on one hand and to probe spatial scales from > 10 pc in the lowdensity striations of the ambient cloud (see Fig. 1), down to ∼ 0.01-0.03 pc for the fibers of dense filaments (see Fig. 2). In nearby Galactic regions (at d ∼ 150-500 pc), this corresponds to angular scales from > 5 deg or more down to ∼ 20 or less.
Low-density striations are remarkably ordered structures in an otherwise chaotic-looking turbulent medium. While the exact physical origin of both low-density striations (Heyer et al., 2016;Tritsis & Tassis, 2016Chen et al., 2017) and high-density fibers (e.g. Clarke et al., 2017;Zamora-Avilés et al., 2017) is not well understood and remains highly debated in the literature, there is little doubt that magnetic fields are involved. For instance, Tritsis & Tassis (2016) modeled striations as density fluctuations associated with magnetosonic waves in the linear regime (the column density contrast of observed striations does not exceed 25%). These waves are excited as a result of the passage of Alfvén waves, which couple to other MHD modes through phase mixing Right panel: Volume density image from the simulations. Left panel: Zoomed-in column density view of a single striation, showing the "sausage" instability setting in, with characteristic imprints in both the magnetic-field and the column-density distribution. In both panels, the drapery pattern traces the magnetic field lines and the mean direction of the magnetic field is indicated by a black arrow. The passage of Alfvén waves excites magnetosonic modes that create compressions and rarefactions (colorbar) along field lines, giving rise to striations. The simulated data in both panels have been convolved to an effective spatial resolution of 0.012 pc, corresponding to the 18 HPBW of B-BOP at 200 µm.
(see Fig. 3, right panel). In contrast, Chen et al. (2017) proposed that striations do not represent real density fluctuations, but are rather a line-of-sight column density effect in a corrugated layer forming in the dense postshock region of an oblique MHD shock. High-resolution polarimetric imaging data would be of great interest to set direct observational constraints and discriminate between these possible models. Specifically, the magnetosonic wave model predicts that a zoo of MHD wave effects should be observable in these regions. One of them, that linear waves in an isolated cloud should establish standing waves (normal modes) imprinted in the striations pattern, has recently been confirmed in the case of the Musca cloud . Other such effects include the "sausage" and "kink" modes (see Fig. 3, left panel), which are studied extensively in the context of heliophysics (e.g. Nakariakov et al., 2016), and which could open a new window to probe the local conditions in molecular clouds .
A first specific objective of B-BOP observations will be to test the hypothesis, tentatively suggested by Planck polarization results (cf. Planck int. res. XXXIII, 2016) that the magnetic field may become nearly parallel to the long axis of star-forming filaments in their dense interiors at scales < 0.1 pc, due to, e.g., gravitational or turbulent compression (see Fig. 4) and/or reorientation of oblique shocks in magnetized colliding flows (Fogerty et al., 2017). A change of field orientation inside dense star-forming filaments is also predicted by numerical MHD simulations in which gravity dominates and the magnetic field is dragged by gas flowing along the filament axis ( observed in the velocity field of some massive infrared dark filaments (Peretto et al., 2014). An alternative topology for the field lines within dense molecular filaments often advocated in the literature is that of helical magnetic fields wrapping around the filament axis (e.g. Fiege & Pudritz, 2000;Stutz & Gould, 2016;Tahani et al., 2018). As significant degeneracies exist between different models because only the plane-of-sky magnetic field is directly accessible to dust polarimetry (cf. Reissl et al., 2018;Tomisaka, 2015), discriminating between these various magnetic topologies will require sensitive imaging observations of large samples of molecular filaments for which the distribution of viewing angles may be assumed to be essentially random. One advantage of the model of oblique MHD shocks (e.g. Chen & Ostriker, 2014;Inoue et al., 2018;Lehmann & Wardle, 2016) is that it could potentially explain both how dense filaments maintain a roughly constant ∼ 0.1 pc width while evolving (cf. Seifried & Walch, 2015) and why the observed spacing of prestellar cores along the filaments is significantly shorter than the characteristic fragmentation scale of 4 × the filament diameter expected in the case of non-magnetized nearly isothermal gas cylinders (e.g. Inutsuka & Miyama, 1992;Nakamura et al., 1993;Kainulainen et al., 2017).
A second specific objective of B-BOP observations will be to better characterize the transition column density at which a switch occurs between filamentary structures primarily parallel to the magnetic field (at low N H2 ) and filamentary structures preferentially perpendicular to the magnetic field (at high N H2 ) (see Planck int. res. XXXV, 2016, and § 2.2 above). Based on a detailed anal-ysis of numerical MHD simulations, Soler & Hennebelle (2017) postulated that this transition column density depends primarily on the strength of the magnetic field in the parent molecular cloud and therefore constitutes a key observable piece of information. Moreover, Chen et al. (2016) showed that, in their colliding flow MHD simulations, the transition occurs where the ambient gas is accelerated gravitationally from sub-Alfvénic to super-Alfvénic speeds. They also concluded that the nature of the transition and the 3D magnetic field morphology in the super-Alfvénic region can be constrained from the observed polarization fraction and dispersion of polarization angles in the plane of the sky, which provides information on the tangledness of the field.
As a practical illustration of what could be achieved with B-BOP, a reference polarimetric imaging survey would map, in Stokes I, Q, U at 100 µm, 200 µm, 350 µm, the same ∼ 500 deg 2 area in nearby interstellar clouds imaged by Herschel in Stokes I at 70-500 µm as part of the Gould Belt, HOBYS, and Hi-GAL surveys (André et al., 2010;Motte et al., 2010;Molinari et al., 2010). To first order, the gain in sensitivity of B-BOP over SPIRE & PACS on Herschel would compensate for the low degree of polarization (only a few %) and make it possible to obtain Q and U maps of polarized dust emission with a signal-to-noise ratio similar to the Herschel images in Stokes I. Assuming the B-BOP performance parameters given in Table 1 (see also Table 4 of Roelfsema et al., 2018, and Table 1 of Rodriguez et al., 2018) and an integration time of ∼ 2 hr per square degree, such a survey would reach a signal-to-noise ratio of 7 in Q, U intensity at both 200 µm and 350 µm in low column density areas with A V ∼ 0.2 (corresponding to the diffuse, cold ISM), for a typical polarization fraction of 5% and a typical dust temperature of T d ∼ 15 K. The same survey would reach a signal-to-noise ratio of 5 in Q, U at 100 µm down to A V ∼ 1. The entire survey of ∼ 500 deg 2 would require ∼ 1500 hr of telescope time, including overheads. It would provide key information on the magnetic field geometry for thousands of filamentary structures spanning ∼ 3 orders of magnitude in column density from lowdensity subcritical filaments in the atomic (HI) medium at A V < 0.5 to star-forming supercritical filaments in the dense inner parts of molecular clouds at A V > 100.

Key advantages of B-BOP over other polarimetric facilities
Far-IR/submillimeter polarimetric imaging from space with B-BOP will have unique advantages, especially in terms of spatial dynamic range and surface brightness dynamic range. Studying the multi-scale physics of star formation within molecular filaments requires a spatial dynamic range of ∼ 1000 or more to simultaneously probe scales > 10 pc in the parent clouds down to ∼ 0.01 pc in the interior of star-forming filaments (corresponding to angular scales from ∼ 18 to > 5 deg in the nearest molecular clouds -see Fig. 1). Such a high spatial dynamic range was routinely achieved with Herschel in non-polarized imaging, but has never been obtained in ground-based submillimeter continuum observations. It will be achieved for the first time with B-BOP in polarized far-IR imaging. The angular resolution and surface brightness dynamic range of B-BOP will make it possible to resolve 0.1 pcwide filaments out to 350 pc and to image a few % polarized dust emission through the entire extent of nearby cloud complexes (cf. Fig. 1), from the low-density outer parts of molecular clouds (A V < ∼ 0.5) all the way to the densest filaments and cores (A V > 100). In comparison, Planck had far too low resolution (10 at best in polarization) to probe the magnetic field within dense 0.1 pc-wide filaments or detect faint 0.1 pc-wide striations. Near-infrared polarimetry cannot penetrate the dense inner parts of star-forming filaments, and ground-based or air-borne millimeter/submillimeter polarimetric instruments, such as SCUBA2-POL (the polarimeter for the Sub-millimetre Common-User Bolometer Array 2 -Friberg et al., 2016), NIKA2-POL (the polarization channel for the New IRAM KID Arrays 2 -cf. Adam et al., 2018), HAWC+ (the High-resolution Airborne Wideband Camera-Plus for SOFIA, the Stratospheric Observatory for Infrared Astronomy - Harper et al., 2018), or BLAST-TNG (the next generation Balloonborne Large Aperture Submillimeter Telescope for Polarimetry - Galitzki et al., 2014), will lack the required sensitivity and dynamic range in both spatial scales and intensity. Figure 5. Surface-brightness sensitivity of B-BOP for wide-field polarimetric imaging compared to other existing or planned polarimetric facilities. The total surface-brightness level required to detect polarization (i.e., Stokes parameters Q, U) with a signalto-noise ratio of 7 per resolution element (e.g. 9 pixel at 200µm for B-BOP) when mapping 1 deg 2 in 2 hr assuming 5% fractional polarization is plotted as a function of wavelength for each instrument (SOFIA-HAWC+, B-BOP, BLAST-TNG, CCAT-p, SCUBA2-POL, CSST, NIKA2-POL). For comparison, the typical surface-brightness level expected in total intensity from the diffuse outer parts of molecular clouds (A V = 1) is shown for two representative dust temperatures (T d = 10 K and T d = 14 K, assuming simple modified blackbody emission with a dust emissivity index β = 2), as well as the SED of the halo of the nearby galaxy M82 (cf. Galliano et al., 2008;Roussel et al., 2010).
More specifically, BLAST-Pol, the Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry operating at 250, 350, and 500 µm (Fissel et al., 2010(Fissel et al., , 2016, has only modest resolution (30 -1 ), sensitivity, and dynamic range. HAWC+, the far-infrared camera and polarimeter for SOFIA (Dowell et al., 2010); and BLAST-TNG (cf. Dober et al., 2014), both benefit from a larger 2.5-m primary mirror equivalent to that of SPICA and thus have comparable angular resolution, but are not cooled and therefore are two to three orders of magnitude less sensitive (Noise Equivalent Power NEP > 10 −16 W Hz −1/2 ) than B-BOP. Stated another way, the mapping speed 6 of B-BOP will be four to five orders of magnitude higher than that of HAWC+ or BLAST-TNG. Future ground-based submillimeter telescopes on high, dry sites such as CCAT-p (the Cerro Chajnantor Atacama Telescope, prime) and CSST (the Chajnantor Submillimeter Survey Telescope) will benefit from larger aperture sizes (6 m and 30 m, respectively) and will thus achieve higher angular resolution than SPICA at 350 µm, but will be limited in sensitivity by the atmospheric background load on the detectors and in spatial dynamic range by the need to remove atmospheric fluctuations. The performance and advantage of B-BOP over other instruments for wide-field dust polarimetric imaging are illustrated in Fig. 5.
Dust polarimetric imaging with ALMA at λ ∼ 0.8-3 mm will provide excellent sensitivity and resolution, but only on small angular scales (from ∼ 0.02 to ∼ 20 ). Indeed, even with additional observations with ACA (the ALMA Compact Array), the maximum angular scale recoverable by the ALMA interferometer remains smaller than ∼ 1 in total intensity and ∼ 20 in polarized emission (see ALMA technical handbook) 7 . This implies that ALMA polarimetry is intrinsically insensitive to all angular scales > 20 , corresponding to structures larger than 0.015-0.05 pc in nearby clouds. Using multiconfiguration imaging, ALMA can achieve a spatial dynamic range of ∼ 1000, comparable to that of Herschel or B-BOP, but only for relatively high surface brightness emission. Because ALMA can only image the sky at high resolution, it is indeed ∼ 2-3 orders of magnitude less sensitive to low surface brightness emission than a cooled single-dish space telescope such as SPICA. Expressed in terms of column density, this means that ALMA can only produce polarized dust continuum images of compact objects with N H2 > ∼ 10 23 cm −2 (such as substructure in distant, massive supercritical filamentssee Beuther et al., 2018) at significantly higher resolution (∼ 1 or better) than SPICA, while polarimetric imaging of extended, low column density structures down to N H2 > ∼ 5 × 10 20 cm −2 (such as subcritical filaments and striations) will be possible with B-BOP. Furthermore, the small size of the primary beam (∼ 0.3 -1 at 0.8-3 mm) makes mosaicing of wide (> 1 deg 2 ) fields impractical and prohibitive with ALMA. In practice, ALMA polarimetric studies of star-forming molecular clouds will provide invaluable insight into the role of magnetic fields within individual protostellar cores/disks and will be very complementary to, but will not compete with, the B-BOP observations discussed here which target the role of magnetic fields in the formation and evolution of filaments on larger scales.

THE TURBULENT MAGNETIZED INTERSTELLAR MEDIUM
Magnetic fields and turbulence are central to the dynamics and energetics of gas in galactic disks, but also in their halos and possibly in the cosmic web at large. These two intertwined actors of cosmic evolution coupled to gravity drive the formation of coherent structures from the warm and hot tenuous gas phases to the onset of star formation in molecular clouds.
Reaching a statistical description of turbulence in the magnetized ISM is an outstanding challenge, because its extreme characteristics may not be reproduced in laboratory experiments nor in numerical simulations. This challenge is of fundamental importance to Astrophysics, in particular to understand how galaxies and stars form, as well as the chemical evolution of matter in space. This section summarizes the contribution we expect B-BOP to bring to this ambitious endeavor.

Interstellar magnetic fields
Magnetic fields pervade the multi-phase ISM of galaxies. In the Milky Way, and more generally in local universe galaxies, the ordered (mean) and turbulent (random) components of interstellar magnetic fields are comparable and in near equipartition with turbulent kinetic energy (Heiles & Troland, 2005;. The galactic dynamo amplification is saturated, but exchanges between gas kinetic and magnetic energy still occur and are of major importance to gas dynamics. Magnetic fields are involved in the driving of turbulence and in the turbulent energy cascade (Subramanian, 2007). The two facets of interstellar turbulence: gas kinematics and magnetic fields, are dynamically so intertwined that they may not be studied independently of each other.
The multiphase magnetized ISM is far too complex to be described by an analytic theory. Our understanding in this research field follows from observations, MHD simulations and phenomenological models. MHD simulations allow us to quantify the non-linear ISM physics but within numerical constraints that limit their scope. They may guide the interpretation observations but alone they do not provide conclusive answers because they are very far from reproducing the high Reynolds (R e ) and magnetic Prandtl numbers (P m ) 8 of interstellar turbulence (Kritsuk et al., 2011). The fluctuation dynamo and shock waves contribute to produce highly intermittent magnetic fields where the field strength is enhanced in localized magnetic structures. The volume filling factor of these structures decreseases for increasing values of the magnetic Reynolds number R m = R e × P m (Schekochihin et al., 2002;Brandenburg & Subramanian, 2005). The inhomegeneity in the degree of magnetization of matter associated with intermittency is an essential facet of interstellar turbulence (Falgarone et al., 2015;Nixon & Pringle, 2019), which simulations miss because they are far from reproducing the interstellar values of R m . In this context, to make headway, we must follow an empirical approach where a statistical model of interstellar turbulence is inferred from observations.

The promise of B-BOP
B-BOP will image dust polarization with an unprecedented combination of sensitivity and angular resolution, providing a unique data set (Sect. 2.4) to characterize the magnetic facet of interstellar turbulence. This leap forward will open an immense discovery space, which will revolutionize our understanding of interstellar magnetic fields, and their correlation with matter and gas kinematics.
At 200 µm, for the SED of the diffuse ISM (Planck int. res. XVII, 2014), the surface brightness sensitivity of B-BOP is three orders of magnitude greater than that of the Planck 353 GHz all-sky map for deep imaging (10 hr per square degree) and a few hundred times better for faster mapping (2 hr per square degree). The analysis of dust polarization at high Galactic latitude with the Planck data is limited by sensitivity to an effective angular resolution of ∼ 1 • in the diffuse ISM and 10 in molecular clouds where the column density is larger than 10 22 H cm −2 (Planck 2018 res. XII, 2019), while B-BOP will map dust polarization with a factor ∼ 20-70 better resolution at 100-350 µm.
Dust polarization probes the magnetic field orientation in dust-containing regions, i.e., mostly in the cold and warm phases of the ISM, which account for the bulk of the gas mass, and hence of the dust mass. These ISM phases comprise the diffuse ISM and star-forming molecular clouds. They account for most of the gas turbulent kinetic energy in galaxies (Hennebelle & Falgarone, 2012). Thus, among the various means available to map the structure of interstellar magnetic fields, dust polarization is best suited to trace the dynamical coupling between magnetic fields, turbulence, and gravity in the ISM. This interplay is pivotal to ISM physics and star formation. It is also central to cosmic magnetism because it underlies dynamo processes (Subramanian, 2007).
Observations have so far taught us that magnetic fields are correlated with the structure of matter in both the diffuse ISM and in molecular clouds (Clark et al., 2014; Planck int. res. XXXII, 2016; Planck int. res. XXXV, 2016) but this correlation does not fully describe interstellar magnetism. Data must also be used to characterize the intermittent nature of interstellar magnetic fields.
The magnetic structures identified in MHD simulations ( Fig. 6) may be described as filaments, ribbons or sheets with at least one dimension commensurate with dissipation scales of turbulent and magnetic energy in shocks, current sheets or through ambipolar diffusion (Momferratos et al., 2014;Falgarone et al., 2015). While the viscous and Ohmic dissipation scales of turbulence are too small to be resolved by B-BOP, turbulence dissipation due to ion-neutral friction is expected to occur on typical scales between ∼ 0.03 pc and ∼ 0.3 pc (cf. Momferratos et al., 2014), which is well within the reach of B-BOP for matter in the local ISM.
Although dust polarization does not measure the field strength, the polarization angle may be used to map these magnetic structures, as illustrated in Fig. 6. The figure shows that the largest values of the increment of the polarization angle, ∆Φ, delineate structures that tend to follow those of intense dissipation of turbulent energy. B-BOP will allow us to identify magnetic structures such as those in Fig. 6 even if their transverse size is unresolved because i) they are highly elongated and ii) their spatial distribution in the ISM is fractal.
Regions of intermittency in interstellar turbulence correspond to rare events. Their finding requires obtaining large data sets combining brightness sensitivity and angular resolution, as illustrated by the CO observations with the Institut de Radioastronomie Millimétrique (IRAM) 30m telescope analyzed by Hily-Blant et al. (2008) and Hily-Blant & Falgarone (2009). B-BOP has the unique capability to extend these pioneering studies of the intermittency of gas kinematics to dust polarization observations, tracing the structure of magnetic fields (Sect. 2.1.2), with a comparable angular resolution. B-BOP will also greatly strengthen their statistical significance by covering a total sky area more than two orders of magnitude larger. B-BOP holds promises to reveal a rich array of magnetic structures, characterizing the intermittency of interstellar magnetic fields. Planck data, at a much coarser scale, gives a first insight at the expected outcome of the observations illustrated in Fig. 7. Magnetic structures will be identified in the data as locations where the probability distributions of the increments of the polarization angle, and the Stokes Q/I and U/I ratios, depart from Gaussian distributions. Compared to Planck, B-BOP will only map a small fraction of the sky (∼ 1% for nearby molecular clouds and diffuse ISM observed away from the Galactic plane) but it will probe the field structure on much smaller scales (by a factor 30 or more) where the surface density of magnetic structures is expected to be much larger.

Observing strategy
B-BOP will considerably expand our ability to map the structure of interstellar magnetic fields. These data will be complementary to a diverse array of polarization observations of the Galaxy.
Stellar polarization surveys will be combined with Gaia astrometry (e.g. Tassis et al., 2018) to build 3D maps of the magnetic fields in the Galaxy but with a rather coarse resolution, comparable to that of the density structure of the local ISM in Lallement et al. (2018).
Synchrotron observations at radio wavelengths with the Square Kilometer Array (SKA) and its precursors will probe the structure of magnetic fields (Dickinson et al., 2015;Haverkorn et al., 2015), in particular in ionized phases through Faraday rotation (Gaensler et al., 2011;Zaroubi et al., 2015). SKA will also provide Faraday rotation measurements toward ∼ 10 7 extragalactic sources (Johnston-Hollitt et al., 2015), which will be available for comparison with B-BOP dust polarization data as illustrated in the pioneering study of Tahani et al. (2018).
B-BOP will allow us to study interstellar turbulence over an impressive range of physical scales and astrophysical environments from the warm ISM phases to molecular clouds. Observations of nearby galaxies are best suited to probe the driving of turbulence in relation to galaxy dynamics (spiral structure, bars, galaxy interaction, outflows) and stellar feedback as discussed in Sect. 6. Galactic observations will probe the inertial range of turbulence over 4 to 5 orders of magnitude from the injection scales (∼ 100 pc -1 kpc) down to 0.01 pc. The smallest physical scales will be reached by observing interstellar matter nearest the Sun, away from the Galactic plane: the diffuse ISM at high Galactic latitudes and star-forming molecular clouds in the Gould Belt. These sky regions are best suited for the study of turbulence because the overlap of structures along the line of sight is minimized. The Gould Belt clouds are already part of the filament science case in Sect. 2. This survey will include star-forming clouds and diffuse clouds representative of the cold neutral medium. Deeper polarimetric imaging of high Galactic latitude fields (10 hr per deg 2 ), sampling regions of low gas column density (A V ∼ 0.1 to 0.3), will allow us to probe turbulence in the warm ISM phases. These deep imaging observations could potentially share the same fields as those used to carry out a SPICA-SMI cosmological survey. The size of the area that may be mapped to that depth (of order ∼ 100 deg 2 ) will be optimized with the needs of this survey. Altogether, we estimate that B-BOP will cover a total area of about 500 deg 2 away from the Galactic Plane including diffuse ISM fields at high Galactic latitudes, which will be available to study turbulence in diverse interstellar environments. At the angular resolution of B-BOP at 200 µm, these data correspond to a total of 2 × 10 7 polarization measurements. This number is 20 times larger than the statistics offered by the Planck polarization data. The gain in angular resolution and sensitivity is so large, that B-BOP will supersede Planck in terms of data statistics, even if the maps used cover only ∼ 1% of the sky.
A wealth of spectroscopic observations of HI and molecular gas species, tracing the gas density, column density and kinematics, will become available before the launch of B-BOP with SKA and its precursors (McClure-Griffiths et al., 2015), and the advent of powerful heterodyne arrays on millimeter ground-based telescopes, e.g. the Large Millimeter Telescope and the IRAM 30m telescope. Furthermore, we will be able to investigate the link between coherent magnetic structures and turbulent energy dissipation observing main ISM cooling lines from H 2 , CII, and OI with the SPICA mid and far-IR spectrometers SMI and SAFARI. These complementary data from SPICA and ground-based observatories will be combined to characterize the turbulent magnetized ISM statistically. The data analysis will rely on on-going progress in the development of statistical methods (e.g. Makarenko et al., 2018), which we will use to characterize the structure of interstellar magnetic fields and their correlation with gas density and velocity. This process will converge toward an empirical model of interstellar turbulence, which will be related to ISM physics comparing data and MHD simulations.

Current state of the art
In molecular clouds, protostellar dense cores are the "seeds" where the gravitational force proceeds to form stars. Class 0 objects are the youngest known accreting protostars: most of their mass is still in the form of a dense core/envelope (M env M ) and this phase is characterized by high accretion rates of gas from the dense core onto the central stellar object, accompanied by ejection of powerful highly collimated flows (André et al., 2000;Dunham et al., 2014).
How many stars can be formed out of a typical molecular cloud depends not only on the physical conditions in filamentary structures, but also on the detailed manner the gravitational collapse proceeds within individual protostellar cores (i.e. would it form a single/binary stellar system, one or several low-mass stars or high-mass stars?).
By the end of the protostellar phase, the star has gained most of its final mass: understanding the role of magnetic fields during the protostellar stage is therefore crucial to clarify how they affect some of the most remarkable features of the star formation process, such as the distribution of stellar masses, the stellar multiplicity, or the ability to host planet-forming disks (McKee & Ostriker, 2007;Li et al., 2014).
The development of numerical magnetohydrodynamical models describing the collapse of protostellar cores and the formation of low-mass stars, has opened new ways to explore in more details the physical processes responsible for the formation of solar-type stars. MHD models suggest that protostellar collapse proceeding with initially strong and well aligned magnetic field produces significantly different outcomes than hydrodynamical or weakly magnetized models (Fiedler & Mouschovias, 1993;Hennebelle & Fromang, 2008;Masson et al., 2016). For example, if the field is strong enough and well coupled to the core material, magnetic braking will regulate the formation of disks and multiple systems during the Class 0 phase. This has been the focus of recent studies ( because it could potentially explain the low-end of the size distribution of protostellar disks (e.g. Maury et al., 2010;Segura-Cox et al., 2018;Maury et al., 2019).
All protostellar cores are magnetized to some level and current observations suggest that at least in some cases the magnetic field at core scale is remarkably well organized, pointing toward scenarii with strong field even at the high column densities typical of protostellar cores (e.g., IRAS 4A, G31.41, G240.31, NGC 6334, L1157, B335: Girart et al., 2006Girart et al., , 2009Qiu et al., 2014;Li et al., 2015;Stephens et al., 2013;Galametz et al., 2018;Maury et al., 2018), while in other cases (e.g., Girart et al. 2013;Hull et al. 2017;Ching et al. 2017, and Fig. 8 for an example in the G14 massive star forming region) the core-scale magnetic field shows very complex morphology. In Fig. 8, for instance, it is noteworthy that the northern hub of the G14 region, with a more uniform magnetic field, has a lower level of fragmentation than the southern hub (that shows a more perturbed magnetic field). These observations suggest the field may remain organized at scales where collapse occurs in most solar-type progenitors, and also at least some of the massive cores Beuther et al., 2018). Current results may be biased, however, because present single-dish facilities selectively trace magnetic fields from the brightest regions within star-forming cores (dust polarization is only detected at the highest column densities, see Fig. 8).

Role of magnetic fields in controlling the typical outcome of protostellar collapse
B-BOP can perform statistical studies in unprecedentedly large samples of protostellar cores, testing for example whether the magnetic field in cores is directly inherited from their environment (if the magnetic field in low-density filamentary structures, see § 2.3, connects to the magnetic field in high-density cores, which is expected in the strong field case), or if the field in cores is disconnected from the local field in the progenitor cloud (weak field case). While the ALMA interferometer can only provide constraints on the magnetic field topology at the smallest scales (∼ 0.02 to ∼ 20 , i.e. < 5000 au in Gould Belt clouds), the polarization capabilities of other facilities probing larger spatial scales (SMA, NIKA2, and HAWC+) are severely sensitivity limited. Accordingly, studies linking cores and filaments can currently be carried out in bright, massive star-forming regions mostly (see Fig. 8), and only in a handful of nearby solar-type protostellar cores.
Some indications have been found, in small samples (< 20 objects), that the topology of the magnetic field at core scales may be linked to the distribution of angular momentum in solar-type cores, and hence that the magnetic field may be of paramount importance to set the initial conditions for the formation of protoplanetary disks and multiple systems (see Fig. 9 and Galametz et al. 2018;Segura-Cox et al. 2018). An example of the type of studies that B-BOP could extend to the full mass function of protostellar cores in a statistical fashion is shown in Fig. 8 and Fig. 9: these two figures illustrate the tentative link between the magnetic field topology at core scale and the disks and multiplicity fraction found within protostellar cores at smaller scale. The magnetic field properties found with B-BOP at dense core scales can be compared with the protostellar properties observed with interferometers at smaller scales to build correlation diagrams similar to the one shown in Fig. 9. In this way, B-BOP observations can test the hypothesis, tentatively suggested by current studies of the brightest protostars, that magnetic fields regulate the formation of disks and multiple systems during the main accretion phase. Observations of large samples of protostars could be carried out thanks to the sensitivity and spatial resolution of B-BOP, which is crucially needed not only to populate diagrams such as Fig. 9, but also because only statistics will allow us to solve the degeneracy induced by projection effects intrinsically linked to dust polarization (tracing only the magnetic field component in the plane of the sky). Moreover, B-BOP will provide information on the geometry of magnetic field lines across the full protostellar core mass distribution, probing different behaviors in different mass regimes, and potentially as a function of environment in different star-forming regions. Figure 9. Potential role of the magnetic field topology at core scales in the formation of disks and multiple systems. Top: Magnetic field (red/orange line segments, from dust polarization observations with the SMA at 850 µm) in two solar-type Class 0 protostellar cores (Galametz et al., 2018). The blue arrows indicate the jet/rotation axis of these cores, aligned with the core-scale magnetic field in L1157 (left), and mostly orthogonal to it in L1448N (right). Bottom: Level of core rotation (from kinematic observations at core scales, Yen et al. 2015 and Gaudel et al. in prep) as a function of the misalignment between the rotation axis and the magnetic field (observed at core scale with the SMA - Galametz et al. 2018) in a sample of nearby Class 0 sources. There is a hint that large misalignments of the magnetic field at core scales lead to sources with large rotational gradients and multiple systems at smaller scales (red symbols).
The angular resolution and surface brightness dynamic range of B-BOP will make it possible to resolve most ∼ 2000-20000 au protostellar cores in nearby starforming regions out to 250 pc. A wide-field B-BOP survey of all nearby clouds as envisaged in § 2.3 (∼ 2 hr per square degree) will map dust polarization (fraction > 1%) at core scales with signal-to-noise ratio > 7, in complete populations of > ∼ 1000 protostars (Class 0 and Class I) and their parent cloud/environment, from massive protostellar cores down to the low-mass progenitors of solar-type stars. In contrast, current millimeter/submillimeter polarimetric instruments, such as SCUBA2-POL, NIKA2-POL, SOFIA/HAWC+, or BLAST-TNG (cf. § 2.4) are limited to the subset of the ∼50-100 brightest cores, and without the important context provided by the magnetic field information in the parental clouds.

ROLE OF MAGNETIC FIELDS IN HIGH-MASS STAR AND CLUSTER FORMATION
As mentioned in § 2.2, supercritical molecular filaments are believed to be the preferred birthplaces of solartype stars. It is however unclear whether the filamentary paradigm -or an extension of it, based on unusually high line masses or levels of turbulence (e.g.  -also applies to high-mass star formation and can lead to the quasi-static formation and monolithic collapse of massive prestellar cores, then forming highmass protostars. The most recent observational results, partly obtained with Herschel, suggest that high-mass stars and stellar clusters form in denser, more dynamical filamentary structures called ridges 9 which exceed the critical line mass of an isothermal filament by up to ∼ two orders of magnitude (see, e.g., the review by Motte et al., 2018b). Such massive structures may originate from highly dynamical events at large scales like converging flows and cloud-cloud collisions, continuing on median scales through the global collapse and filament feeding of ridges. The role of magnetic fields in this scheme is poorly known and may be as crucial as for low-mass star formation.
In the hypothesis of large-scale cloud collapse, gravity overcomes the magnetic field support and the magnetic field follows the infall gas streams from cloud scales (∼100 pc) to accumulation points at scales between ∼1 pc to ∼0.1 pc, with a typical hourglass geometry toward these accumulation points (Girart et al., 2009;Cortes et al., 2016). Large-scale collapse leads to very dense, massive structures at pc scales, which are either spherical (hubs) or elongated (ridges) (Hartmann & Burkert, 2007;Schneider et al., 2010;Hill et al., 2011;Peretto et al., 2013). Pilot works with ground-based facilities (SMA) toward the DR21 ridge in Cygnus X Ching et al., 2017) show that the magnetic field is ordered at the scale of the ridge and mostly perpendicular to its main axis, as for low-mass supercritical filaments (cf. § 2.2 and Fig. 1), suggesting mass accumulation along field lines. However, while largescale collapse and strong ordered magnetic fields are probably a key ingredient, the detailed physical processes at the origin of ridges remain, for now, a mystery.
At some point, ridges fragment into hundreds of protostellar cores in local, short, but violent bursts of star formation, leading to exceptionally large instantaneous star formation rates (Nguyen Luong et al., 2011;Louvet et al., 2014). This clustered mode of fragmentation in ridges may differ in nature from the filamentary mode of fragmentation leading to low-mass star formation at significantly lower average densities. As a matter of fact, top-heavy core mass functions, overpopulated with highmass protostellar cores, begin to be found with ALMA in the massive, young ridges of the Galactic plane (Csengeri et al., 2017;Motte et al., 2018a;Cheng et al., 2018). The magnetic-field configuration (field topology and field strength) inside the hubs and ridges, at scales of a few 0.1 pc, may limit fragmentation (see Commerçon et al., 2011 for MHD simulations and Fig. 8 for recent observations) and favor the formation of massive protostellar cores against their low-mass counterparts. Dynamical processes associated with local accretion streams and global collapse may also favor the growth of high-mass protostellar cores due to competitive accretion (Smith et al., 2009). Elucidating the relative roles of -and coupling between -magnetic fields and dynamics, is therefore of crucial importance for understanding the origin of high-mass stars and their associated clusters. For example, if cluster-scale global collapse is required to form massive stars near the bottom of the gravitational potential well, the collapsing flow should drag the magnetic field on the cluster scale into a more or less radial configuration. If, on the other hand, only localized collapse of a pre-existing massive core is required to produce a massive star (McKee & Tan, 2003), the collapse-induced field distortion is expected to be limited to the smaller, core region.
Observationally, this requires probing the magnetic field configuration from cloud scales (∼100 pc) to protostellar scales (∼0.01 pc) within massive dense ridges/hubs (down to so-called "massive dense cores" or MDCs; 0.1-0.3pc; Motte et al., 2007;Bontemps et al., 2010). Ultimately, a spatial dynamic range 10 as high as 10 4 (from 100 pc to 0.01 pc) is thus needed. In nearby high-mass star-forming regions, located at ∼ 1 to 3 kpc, this translates to angular scales from a few (∼ 2-6) degrees down to ∼ 0.1 . While the 0.1 scale of individual pre-/protostellar cores in these regions is only reachable with SMA or ALMA, the inner scale of massive dense ridges/hubs or MDCs (∼ 0.2 pc, or ∼ 14 − 40 at ∼1-3 kpc distance) can be reached with the angular resolution of B-BOP. This MDC scale is of particular importance for high-mass star and cluster formation since high-mass protostellar cores appear to form in only a subset of MDCs, possibly those with a high level of magnetization (see, e.g., Motte et al., 2018b;Ching et al., 2018 and Fig. 8). With a typical spatial resolution of 10 , Planck polarization data already provide some indications on the magnetic-field geometry at scales between 100 pc and ∼ 3-9 pc, but Planck maps are strongly limited by the confusion arising from several layers of dust emission along the line of sight within the large beam. The spatial resolution of B-BOP is required to separate 10 The spatial dynamic range is defined as the ratio of the largest to the smallest spatial scale accessible to an instrument or an observation (see also Sect. 2.4).
the contributions of these layers and focus on the polarization signal from high-mass star-forming ridges, hubs, and MDCs. The high sensitivity of B-BOP will also be crucial to trace the magnetic field topology all the way to the outer environment of ridges and hubs, where the column density of dust reaches values below A V ∼ 1-2.

MAGNETIC FIELDS IN GALAXIES
Magnetic fields are an important agent that influences the structure and evolution of galaxies (e.g. Tabatabaei et al., 2016). The magnetic pressure in the ISM is comparable in magnitude to the thermal, turbulent, and cosmic ray pressures (e.g. Ferrière, 2001;Boulares & Cox, 1990;Beck, 2007), so the magnetic field contributes significantly to the total pressure which supports a galactic gas disk against gravity. The interplay between the magnetic field, gravity, and turbulence is central to the process of star formation (see, e.g., McKee & Ostriker, 2007;Hennebelle & Falgarone, 2012;Crutcher, 2012), both on the scale of individual stars and filaments (cf. § 2 and § 4), and for the formation of molecular clouds out of the magnetized diffuse interstellar gas (e.g. Körtgen et al., 2018). On even larger scales in galaxies, magnetic fields control the density and distribution of cosmic rays (e.g. Kotera & Olinto, 2011), mediate the spiral arm shock strength (e.g. Shetty & Ostriker, 2006), and may even modulate rotation of galaxy gas disks (e.g. Elstner et al., 2014) and quench high-mass star formation (e.g. Tabatabaei et al., 2018). Magnetic fields play an important role in launching galactic winds and outflows (Heesen et al., 2011), regulate gas kinematics at the diskhalo interface (Henriksen & Irwin, 2016), and ultimately connect galaxies to the intergalactic medium (Bernet et al., 2013).

Current observational status
Significant progress has been made in recent decades to characterize the interstellar magnetic fields of external galaxies, with measurements of the magnetic field strength and orientation obtained for about one hundred nearby galaxies (see the appendix of Beck & Wielebinski, 2013) 11 . This effort has established a broadly consistent picture of the large-scale (> 1 kpc) properties of galactic magnetic fields. As in the Milky Way, the interstellar magnetic field in external galaxies can be described as a combination of large-scale regular fields and small-scale turbulent fields. Observations of face-on spiral galaxies demonstrate that galaxies typically host spiral fields in the disk, with the observed large-scale magnetic field orientations appearing similar to the material spiral arms. The vertical structure of the field is more easily probed via observations of edge-on systems, which typically show an X-shaped structure such that the field tends to become more inclined (and eventually perpendicular) with increasing distance from the midplane. Field strengths vary, but are usually in the range of several to tens of µG, with roughly similar contributions from ordered and random field components (Beck & Wielebinski, 2013).
To date, magnetic field properties in external galaxies have mostly been investigated via observations of synchrotron emission at GHz frequencies (for a review, see Beck & Wielebinski, 2013;. At the same frequencies, Faraday rotation of background radio sources provides an alternative, direct determination of the direction and strength of magnetic field within galaxies. This technique has been used for detailed studies of the Milky Way's magnetic field (e.g. Terral & Ferrière, 2017;Mao et al., 2010), but its application to external galaxies has thus far been limited to the most nearby galaxies (which have a large angular size and thus a sufficient number of bright background sources, e.g. M31, LMC, Han et al., 1998;Gaensler et al., 2005). Faraday rotation measurements for a much larger sample of nearby galaxies is an important science driver for SKA (see e.g. . At other wavelengths, the magnetic field structure of a much smaller number of external galaxies has been surveyed using optical polarization e.g. the Magellanic Clouds (LMC and SMC; Mathewson & Ford 1970), NGC1068, and M51 (Scarrott et al., 1987). Wide-field imaging of polarized extinction at infrared (IR) wavelengths is a newer capability that has been used to probe the magnetic field structure in nearby edge-on galaxies, with results that are generally in good agreement with radio observations (e.g. Clemens et al., 2013;Montgomery & Clemens, 2014). The IR extinction technique is less suited to observations of face-on galaxies, due to the relatively short path length that can produce internal extinction (see e.g. Pavel & Clemens, 2012, for the case of M51).
B-BOP will probe the magnetic field structure via observations of dust polarization in emission. An important advantage of this technique compared to radio synchrotron observations is that it traces the magnetic field structure in the cold gas where star formation occurs, with minimal contamination from the warm ionized gas in the halo (e.g. Mao et al., 2015). Studies in nearby galaxies further suggest that the magnetic field is coupled to the interstellar gas independently of the star formation rate (e.g. Schinnerer et al., 2013;Tabatabaei et al., 2018), highlighting the importance of tracing the field in the neutral ISM. The SCUBA-POL camera on the James Clerk Maxwell Telescope (JCMT), operat-ing at 850 µm, obtained the first such dust polarization observations of external galaxies (e.g. Matthews et al., 2009), but was only able to access extremely bright extragalactic regions, such as the centre of the nearby starburst system M82 (e.g. Greaves et al., 2000). The Planck mission has recently provided all-sky measurements of the polarized submillimeter dust emission, but with an angular resolution (∼ 10 at 353 GHz) sufficient to resolve sub-kiloparsec scales only in the closest Local Group galaxies (< 1 Mpc). With ∼arcsecond resolution, a key opportunity for ALMA will be targeted imaging of the detailed magnetic field structure in extragalactic molecular clouds, but wide-field polarization surveys of nearby galaxies will remain impractical, due to prohibitive integration times for fields larger than a few square arcminutes.

Key opportunities for B-BOP on nearby galaxies
Mapping the structure of interstellar magnetic fields in the cold ISM of nearby galaxies is crucial to understand how magnetic fields influence gas dynamics in galaxies, and in particular the role of the field in regulating star formation, driving galactic outflows, and fuelling galactic nuclei. Observational studies of these processes in nearby galaxies complement Milky Way studies, which typically have superior spatial resolution, but may be limited by distance ambiguity and line-of-sight confusion. Among current and near-future facilities, only B-BOP will be able to make these measurements across a representative sample of external galaxies, probing a much wider range of ISM conditions than those encountered in the Milky Way, and to conduct spatially complete mapping of the field structure in Local Group targets.
To highlight B-BOP's unique capabilities for such an effort, Fig. 10 shows the estimated signal-to-noise ratio in polarized intensity for several iconic nearby galaxies at 100 µm (top row) and 200 µm (bottom row) after 2 hr on-source integration with B-BOP. These maps are constructed assuming the performance parameters given in Table 1 and a conservative polarization fraction of 1%. At the intrinsic resolution of the 100 µm band, the sensitivity is sufficient for tracing the detailed polarization structure of bright features such as spiral arms, bar dust lanes, and galaxy centers. For the galaxies in Fig. 10, several hundred independent 100 µm polarization vectors would be obtained. Measurements in the more sensitive 200 µm and 350 µm bands would essentially cover the entire galactic disk within 0.6 R 25 (where R 25 is the optical radius). The excellent signal-to-noise ratio that can be achieved with a modest integration time per galaxy means that B-BOP could conduct the first systematic survey of the polarized far-IR dust emission -and hence magnetic field structure -in ∼ 100 nearby galaxies. This is slightly larger than the combined sample of galaxies targeted by the VNGS and KINGFISH Herschel nearby galaxy projects, and would require only 200 hr of onsource observing time. In the remainder of this section, we highlight some of the potential science drivers for such a B-BOP nearby galaxy survey.

Testing and refining Galactic dynamo models
The currently favored paradigm for interstellar magnetic fields is that they are amplified by dynamo action. In this scenario, weak primordial fields in young galaxies are quickly amplified by a small-scale turbulent dynamo, which continuously supplies turbulent fields to the ISM after the formation of a galactic disk in 10 9 yr (Schleicher et al., 2010). The large-scale field is then amplified by the mean-field α − Ω dynamo effect (e.g. Ruzmaikin et al., 1988), whereby the combination of differential rotation of the galactic disk (Ω -effect) and helical turbulence (α -effect) presumably driven by supernova explosions (Ferrière & Schmitt, 2000), produce small-scale turbulent and organize some fraction of them into regular large-scale patterns.
The mean field dynamo is expected to generate a regular magnetic field with both poloidal and azimuthal components, and nearly all polarized synchrotron observations of face-on disk galaxies show a large-scale spiral pattern. To date, the magnetic field pitch angles p B and azimuthal structure that have been observed in nearby disk galaxies via observations of polarized radio synchrotron emission are broadly compatible with the predictions of mean field dynamo theory (Fletcher, 2010;Van Eck et al., 2015). Yet the precise nature of the magnetic field generated via the mean-field dynamo depends on properties of the host galaxy. For example, the rotation curve determines the shear strength in a differentially rotating galaxy disk, and hence how the azimuthal field component is generated from the poloidal field. The α-effect -by which a poloidal field component is generated from the azimuthal field -is thought to be powered by supernova explosions, which depend on a galaxy's star formation rate. As our knowledge of external galaxies grows, the logical next step is to refine dynamo models for specific galaxies to include all relevant observed galaxy properties -e.g. the ionized and molecular gas density distributions, rotation curve, star formation rate, gas inflow and outflow rates -and test the model predictions for individual galaxies against the observed properties of the magnetic field. The first attempt to do this systematically for a sample of galaxies (Van Eck et al., 2015) was hampered by inconsistencies in the available radio observations. A sample of galaxies observed with the same instrument at the same resolution and sensitivity is necessary to allow the details of dynamo theory, such as how the dynamo saturates, to Figure 10. Expected signal-to-noise in 100 µm (top row) and 200 µm (bottom row) polarized intensity after 2 hr on-source integration with B-BOP for four nearby (d < 10 Mpc) galaxies. The maps are constructed using Herschel data at 70µm and 250µm as input. We assume the B-BOP performance parameters given in Table 1, a dust spectral index of β = 1.9, and a conservative polarization fraction of 1%. The color scale, which runs between a signal-to-noise ratio of 1 and 100 (top row) or 1 and 1000 (bottom row), uses a logarithmic stretch, with green indicating a signal-to-noise of 10.
be tested against data. B-BOP's moderate resolution, full-disk sampling of the magnetic field structure across a statistically significant sample of nearby galaxies would provide precisely this test.

Magnetic fields and gas flows in barred and spiral galaxies
Mapping the structure of the magnetic field across a sample of nearby galaxies is needed to understand the typical dynamical importance of the field on galactic scales. Of particular interest is how gas flows in galaxies -e.g. gas streaming along spiral arms, inflow along bar dust lanes, and starburst-driven outflows -interact with the field. While independent estimates of the field strength will still be required, B-BOP observations at sub-kiloparsec resolution of the magnetic field structure and complementary spectral line data for tracing interstellar gas kinematics will be extremely valuable for studying the interplay between the field and motions within the cold gas reservoir across the local galaxy population.
In face-on disk galaxies, the large-scale field traced by radio polarization observations tends to follow a spiral pattern. This pattern is expected from mean field dynamo theory, and is not directly connected to a galaxy's baryonic (i.e. gas/stellar) spiral structure. Observations of spiral galaxies indeed show that the field pattern is not always spatially coincident with the spiral arms, and in several cases (most famously NGC6496 Beck, 2007) the ordered field pattern is most pronounced in the interarm region. Some of the large-scale field patterns in galaxies may be due to the combined action of shear and compression in the interstellar gas, which renders the turbulent field anisotropic (and hence ordered) over large scales (e.g. in M51, Fletcher et al., 2011;Mulcahy et al., 2014Mulcahy et al., , 2016. Current observations also suggest that the average pitch angle of the regular spiral field pattern is often similar to the pitch angle of the local spiral arm p M . This is not a direct prediction of mean-field dynamo theory, but would be expected if spiral shocks amplify the magnetic field component parallel to the shock. Significant discrepancies between p B and p M in the inter-arm region, as well as large azimuthal and radial variations in p B , are also observed, the origin of which are not yet well understood. The central regions of barred galaxies are the site of fast radial gas inflow, strong shocks, and intense star formation. Barred galaxies often show strong gas streaming along the shock fronts at the edge of bars, which develop because the gas is rotating faster than the bar pattern. Radio polarization observations of the prototypical barred galaxy NGC1097 (Beck et al., 2005) reveal strongly polarized emission along the bar with field orientations parallel to the gas streamlines. The observed polarization pattern suggests that the field is amplified and stretched by shear in the compression region, and that the field is frozen into the gas and aligned with the gas flow over a large part of the bar. If this result holds generally in barred galaxies, the polarization pattern in bars -especially using a tracer that preferentially  Figure 11. Expected signal-to-noise in polarized intensity for B-BOP observations of the LMC after 50 hr on-source integration. The maps are constructed using Herschel data at 100µm and 250µm as input (Meixner et al., 2013). We assume the B-BOP performance parameters given in Table 1, a dust spectral index of β = 1.9, and a conservative polarization fraction of 1%. The color scale, which runs between a signal-to-noise ratio of 1 and 100, uses a logarithmic stretch, with green indicating a signal-to-noise of 10. probes the dense interstellar gas -would provide important complementary information on the plane-of-sky gas flows to the line-of-sight kinematic information obtained from molecular emission lines. In combination with estimates for the magnetic field strength, information about the magnetic field structure in the central regions of barred galaxies would also provide useful constraints for models of AGN fuelling. One of the main problems in this area is to generate mass inflow rates that are compatible with the observed nuclear activity. Magnetic stress in circumnuclear rings (e.g. Beck et al., 1999) and fast MHD density waves (e.g. Lou et al., 2001) have been proposed as potential mechanisms, but current observational data for the field strength and structure in the inner regions of barred galaxies is not sufficient for a rigorous test of these models.

Magnetism in dwarf galaxies
Due to their slow rotation, the amplification of magnetic fields should be less efficient in dwarf galaxies. Yet observations of radio polarized intensity show that several nearby low-mass galaxies host large-scale ordered fields, e.g. the Magellanic Clouds, NGC4449 and IC10 (Chyży et al., 2003;Gaensler et al., 2005;Mao et al., 2012;Chyży et al., 2016;Heesen et al., 2018). Dwarf galaxies are also more likely to exhibit star formation powered outflows and galactic winds, due to their shallow gravitational potential. The magnetized nature of these outflows has been observed in some dwarf systems (Chyży et al., 2000;Kepley et al., 2010), consistent with some models of a cosmic ray driven dynamo (Siejkowski et al., 2014;Dubois & Teyssier, 2010). To date, all dwarf galaxies with detected ordered magnetic fields are starbursting, participating in a galaxy-galaxy interaction, and/or experiencing significant gas infall, suggesting the importance of enhanced turbulence for the magnetic field properties and evolution of these systems. B-BOP observations of a sample of local dwarf galaxies with a range of masses, interaction properties and star formation histories, would provide valuable input for theories for the amplification of magnetic fields in such systems, and their role in magnetizing the intergalactic medium (IGM).

The Magellanic Clouds
The Large and the Small Magellanic Cloud (LMC, SMC) are the closest gas-rich galaxies to the Milky Way. A B-BOP survey of the Magellanic Clouds would for the first time probe the magnetic field structure in the cold ISM across all spatial scales between the clumpy sub-structure within GMCs (∼ 2 pc) and the galactic disk (several kpc). Observations across such a large range of spatial scales are needed to decipher the dynamical importance of the magnetic field for the inherently hierarchical process of star formation, i.e. from the formation of GMCs out of the diffuse ISM, down to the formation of individual stars. Spatially complete surveys of dust emission in the Magellanic Clouds with ALMA are unfeasible due to their large angular size (∼ 50 and ∼ 10 deg 2 for the LMC and SMC respectively).
As an example of what could be achieved with B-BOP, Fig. 11 shows the estimated signal-to-noise ratio for a 50 hr polarimetric imaging survey of the LMC at 100 µm, 200 µm, and 350 µm. This hypothetical survey would achieve a signal-to-noise ratio of 3 for the polarized intensity at 100 µm for interstellar gas with column densities above ∼ 2.5 × 10 21 cm −2 (equivalent to A V ∼ 0.4 in the LMC- Weingartner & Draine, 2001). This sensitivity would yield ∼0.5 million independent measurements of the magnetic field orientation in the interstellar gas on spatial scales of 2 pc, including in the column density regime of the atomic-to-molecular phase transition. At 200 µm and 350 µm, a similar number of significant detections of the magnetic field orientation would be achieved in even more diffuse gas (∼ 1 × 10 21 cm −2 ≈ 0.15 mag). This represents ∼ two orders of magnitude increase in detail over measurements with Planck's 353 GHz channel in the Magellanic Clouds, and would provide the first spatially complete view of the parsec-scale magnetic field structure in the molecular gas reservoir of any galaxy.

Wavelength dependence of polarization in U/LIRGs and AGNs
The polarization of luminous external galaxies such as Luminous Infrared Galaxies (LIRGs), Ultraluminous Infrared Galaxies (ULIRGs) and active galactic nuclei (AGNs), in the far-IR and submillimeter can arise from synchrotron emission but also from emission or absorption by aligned dust grains in the optically thick clouds that surround young stars and AGN tori (e.g. Efstathiou et al., 1997;Aitken et al., 2002). Information in the far-IR and submillimeter can be combined with information at 10 µm and 18 µm as well as near-IR data from the ground to study the switch in position angle by about 90 deg that is predicted as polarization changes from dichroic absorption at shorter wavelengths to dichroic emission at longer ones. Several highly polarized galaxies in the mid-IR were found by Siebenmorgen & Efstathiou (2001) with Infrared Space Observatory (ISO) and more recently by Lopez- Rodriguez et al. (2018a) with CanariCam on the 10.4-m Gran Telescopio Canarias (GTC). This is a science area where significant progress can be achieved with B-BOP, which will provide sensitive polarization measurements at 100-350 µm for a large sample of luminous external galaxies. Such information is currently available for very few objects. In a recent study of the nearby radio galaxy Cygnus A using data from HAWC+ onboard SOFIA, Lopez- Rodriguez et al. (2018b) showed that this approach can be very useful for unravelling the polarization mechanisms in the infrared and submillimeter and providing an independent method of estimating the contributions of AGN tori and starbursts to the SEDs. Exploring the role of AGNs and star formation in galaxies is a scientific objective of wide interest. The opportunity to study multi-wavelength polarization with B-BOP will be complementary to other methods such as spectroscopy (e.g. González-Alfonso et al., 2017) and traditional SED fitting of the total emission (e.g. Gruppioni et al., 2017).

Distant galaxies and the potential detection of the Cosmic Infrared Background polarization
The build-up of coherent magnetic fields in galaxies and their persistence along cosmic evolution is being investigated with analytical models of the galactic dynamo (e.g. Rodrigues et al., 2019) and numerical simulations of galaxy formation (e.g. Martin-Alvarez et al., 2018). These studies suggest that the mean-field dynamo is effective early in the evolution of galaxies but, today, polarization data available to trace the redshift evolution of galactic magnetic fields are very scarce (Mao et al., 2017). While SKA holds exciting promises to extend observations of cosmic magnetism to the distant universe (e.g. Basu et al., 2018;Mao, 2018), we argue here that B-BOP can also uniquely contribute providing the first polarimetric extragalactic survey at far-IR wavelengths.
To quantify what could be achieved with B-BOP, we consider the point source sensitivity of a polarimetric extragalactic survey for an integration time of 10 hr per deg 2 (Table 1). At the detection limit of Herschel imaging surveys in total intensity, the signal-to-noise ratio of B-BOP in Stokes Q and U is ∼ 200 at 200 µm, and ∼ 100 at 100 µm and 350 µm. This sensitivity needs to be compared to the few existing values of the far-IR polarization fraction for galaxies as a whole.
The net polarization fraction resulting from the integrated emission of galaxies depends on the existence of a coherent mean magnetic field and on viewing angle. In disk galaxies, the polarization angle is aligned with the projection of the galaxy angular momentum vector on the plane-of-the-sky, and the polarization fraction increases from a face-on to an edge-on view. Integrating the Planck dust polarization maps at 353 GHz over a 20 • wide band centered on the Galactic plane, De Zotti et al. (2018) found a polarization fraction p = 2.7%. Within a simple model, p is expected to scale as sin i, where i is the inclination angle of the galaxy axis to the line of sight. For this scaling, the mean p fraction averaged over inclinations is 1.4%. This may be taken as a reference value for spiral galaxies like the Milky Way, but p is likely to be on average lower for distant infrared galaxies. Indeed, SOFIA polarization imaging of the two template starburst galaxies, M82 and NGC253, revealed regions with different polarization orientations, which tend to average out when computing the integrated polarized emission yielding an overall mean p ∼ 0.1% (Jones et al., 2019).
Even if B-BOP only detects polarized emission from only a small fraction of Herschel galaxies, the number of detections will be significant given the present dearth of such measurements. If the detections are numerous, the emission from galaxies could even limit the polarization sensitivity of B-BOP deep surveys. This is a possibility that needs to be assessed. Beyond the study of individual galaxies, we anticipate that the main outcome of a deep polarimetric extragalactic survey with B-BOP could follow from a statistical analysis of the data.
Statistical analysis is the reference in cosmology and much can be learned without detecting sources individually. In particular, the cross-correlation of surveys across the electromagnetic spectrum is a powerful means commonly used. In the far-IR, this is illustrated by the results obtained stacking Herschel data on positions of extragalactic sources in the near-and mid-IR. This approach has been successfully used to statistically identify sources accounting for the the bulk of Cosmic Infrared Background (CIB), although they were too faint to be detected individually (Béthermin et al., 2012;Viero et al., 2013). B-BOP, which will extend these studies to polarization, is uniquely suited to detect -or set tight constraints on -the CIB polarization. We note that the analysis of point sources circumvents the difficulty of separating the CIB from the foreground polarized emission of the diffuse Galactic ISM.
For polarization, data stacking needs to be oriented to align polarization vectors. This can be achieved by using, e.g., galaxy shapes measured from near-IR surveys at the appropriate angular resolution. Another interesting path will be the study of correlations between B-BOP extragalactic survey data and maps of the cosmic web inferred from weak-lensing surveys. The angular momenta of galaxies are not randomly oriented on the sky. The cosmic web environment has a strong influence on galaxy formation and evolution, and tidal gravitational fields tend to locally align the spins of dark-matter halos and galaxies. Such alignments, observed in dark-matter simulations, bear information on galaxy formation and evolution, as well as on the growth of structure in the Universe (Kirk et al., 2015). In this picture, low-mass haloes tend to acquire a spin parallel to cosmic web filaments, while the most massive haloes, which are typically the products of later mergers, have a spin perpendicular to filaments (e.g. Codis et al., 2012;Dubois et al., 2014). Quasar observations provide observational evidence of a correlation between the polarization orientation of galaxies and the large-scale structure of the Universe (Hutsemékers et al., 2014), but only for a small number of sources. B-BOP can uniquely contribute to characterizing this correlation for infrared-luminous galaxies.
The scientific goals outlined here are new and promising but still qualitative. Modelling is required to assess the scientific outcome of a deep polarimetric extragalactic survey with B-BOP and decide on the best observing strategy in terms of survey depth and sky coverage.

CONSTRAINING DUST PHYSICS
The polarization of thermal dust emission depends on the shape, size, composition, and alignment efficiency of dust grains, and also on the 3D structure of the Galactic magnetic field on the line of sight and within the instrument beam. Dust polarized emission can therefore bring specific constraints on the alignment mechanism of dust grains and possibly on the grain shape. Despite this complex nature, it can also be used to constrain the optical properties of aligned dust grains (emissivity, spectral index), which are large grains at thermal equilibrium (Planck int. res. XXI, 2015; Planck int. res. XXII, 2015).
As described below, B-BOP, with its high angular resolution and good wavelength coverage of the polarized dust SED, will be a key instrument to provide new constraints on grain alignment theories and inform the evolution of aligned grain properties from the diffuse ISM to the densest cloud cores.

Probing the grain alignment mechanism
Grain alignment is subordinate to various processes. First, grains must rotate supra-thermally 12 to be well aligned. Radiative torques, or chemical torques resulting from the formation of H 2 molecules on the surface of grains are good candidates for grain spin-up. Second, alignment torques of magnetic, radiative, or mechanical origin (Hoang & Lazarian, 2016) are needed to align the supra-thermally rotating dust grains. Compared to the time evolution of molecular clouds and cores, the alignment of grains along magnetic field lines by radiative torques is a fast process, with a timescale on the order of 10 3 -10 4 yr (Hoang & Lazarian, 2014) 13 , so that situations where grain alignment would be out-ofequilibrium can be safely ignored.
While the mechanisms of dust grain alignment are still debated (see, e.g., Sect. 2.1.1), polarization measurements at UV to optical wavelengths imply that the alignment efficiency of dust grains is sensitive to grain size. Such a behavior, well observed in the Mie regime where absorption and scattering are size-dependent (Bohren & Huffman, 1983), is however more difficult to extract from the polarized thermal SED because the dust temperature and dust spectral index depend more on the grain shape, internal structure, and composition (through its emissivity) than on the exact grain size.
As already mentioned in Sect. 2.1.1, the leading grain alignment theory is Radiative Alignment Torques (RATs) (e.g. Lazarian & Hoang, 2007, and references therein). The RATs alignment mechanism, if present, will lead to characteristic signatures in observations. For instance, in this theory, the alignment efficiency is directly dependent on the angle between the incident radiation field and the magnetic field direction (ψ). When dust is heated by a single nearby star or in starless dense cores where the field is strongly attenuated and anisotropic, the incident radiation field direction is well characterized. Since the magnetic field orientation projected on the plane of the sky can be determined from the polarized signal, mapping dust polarization in such regions can in principle be used to test alignment by RATs. This is J J Figure 12. Sketch of the geometry around a single star dominating the heating of the local ISM. The magnetic field direction is represented by the horizontal dashed lines. The aligned dust grains are sketched as prolate rotating parallelograms. If RATs dominate, dust alignment will be more efficient in regions with low values of ψ, the angle between the stellar radiation and magnetic field directions. Figure  illustrated in Fig. 12 which sketches the relative geometry of the radiation field and magnetic field in a region of the ISM where the radiation field is dominated by a single star. In such regions, the RATs alignment theory predicts a stronger alignment and therefore a higher polarization fraction where ψ is close to 0 deg.
Despite an expected clear signature, direct observational evidence of this angular effect has been scarce. The only positive detection reported in emission is by Vaillancourt & Andersson (2015) who detected a periodic modulation of the dust polarization fraction around the Becklin-Neugebauer Kleinmann-Low (BNKL) object in Orion OMC-1. Some authors also claim to have evidenced a correlation between dust temperature and polarization fraction, as expected for dust grains aligned through the RATs mechanism (Andersson et al., 2011;Matsumura et al., 2011). There are also a few reports that indicate a possible influence of H 2 formation (e.g. Andersson et al., 2013, in IC63). In any case, only a handful of cases have been investigated so far and there is certainly a bias in the literature for publishing detections rather than non-detections. Given the intrinsically tangled nature of magnetic field geometry, chance coincidences are very difficult or even impossible to exclude for these few isolated studies and a statistically representative study is clearly needed.
Such a study has not been carried out so far using Planck all-sky data, essentially because the number of interstellar regions where dust is directly and predominantly heated by a single star is very low at the Planck angular resolution. Attempts to unambiguously detect a statistical increase of the polarization fraction with dust temperature in the Planck data, which would also be attributable to radiation-enhanced spin up and alignment of dust grains, have not led to a strong conclusion. Planck 2018 res. XII (2019) showed that it is possible to disentangle, statistically, between what can be attributed to variations in grain alignment efficiency or grain properties, and what is due to line-of-sight and beam averaging of magnetic field structures. This study demonstrated that there is no strong variations in grain alignment efficiency in the diffuse ISM (up to a column density N H ∼ 2 × 10 22 cm −2 ), but, due to the low angular resolution of the Planck data, could not conclude in the case of the high-density ISM. Planck 2018 res. XII (2019) did not find any correlation either between dust temperature and polarization fraction in the diffuse ISM.
In conclusion, the analysis of the Planck all-sky data have so far not allowed to strongly confirm or rule out any specific grain alignment theory, but have provided an upper limit to the drop of alignement in the diffuse ISM.
Owing to its much higher angular resolution and sensitivity, B-BOP will allow us to systematically map the polarization of dust emission around thousands of individual stars heating the nearby ISM locally. The good coverage of the polarized SED will allow us to measure the temperature of aligned dust grains responsible for the polarized emission. Analysis of the polarization fraction as a function of the angle between the known radiation field and the magnetic field direction derived from polarization will allow us to test, for the first time, the RATs alignment theory with statistical significance, in a way that is not affected by local variations and the complexity of individual objects. At the same time, we will be able to detect, if present, the polarized emission resulting from the small temperature difference between grains heated face-on and edge-on immersed in an anisotropic radiation field (Onaka, 1995). This process, which is only efficient at short wavelength (λ ≤ 100 µm, Onaka, 2000), would appear in B-BOP data as a characteristic difference between the polarization fraction and angle measured at 100 µm and the ones measured by unaffected channels at longer wavelengths.
Altogether, the high resolution, sensitivity, and spectral coverage of B-BOP will set unprecedented, probably unexpected constraints on the physics of grain alignment in star-forming regions, a topic which Planck observations could hardly address.

Dust polarization as a proxy for dust evolution
The wavelength range covered by B-BOP will allow us to disentangle between various dust models. The Wien part of the polarized dust emission is currently not constrained. It is in this wavelength range that dust models present the strongest differences in spectral variations of P/I (Fig. 13), in particular between those where carbon grains are aligned and those where they are not (Draine & Hensley, 2013;Guillet et al., 2018). Moreover, as dust models now predict both emission and absorption properties of dust grain populations in polarization (Siebenmorgen et al., 2014;Draine & Hensley, 2017;Guillet et al., 2018), joint observations of common targets with B-BOP and survey experiments targeting extinction polarization of background stars, such as PASIPHAE , can be used to further test such models. High-resolution polarization observations with B-BOP will allow us to probe dust properties in dense environments and to further characterize dust evolution between diffuse and dense media (e.g. Köhler et al., 2015). Observationally, a number of studies based on emission and extinction data have provided evidence of grain growth within dense clouds. One of the main results is an increase of the far-IR/submillimeter emissivity by a typical factor of 2-3 compared to standard grains in the diffuse medium (cf. Stepnik et al., 2003;Planck early res. XXV, 2011;Ysard et al., 2013;Roy et al., 2013;Juvela et al., 2015), as predicted by calculations of aggregate optical properties (Ossenkopf & Henning, 1994;Köhler et al., 2012Köhler et al., , 2015. Another striking result is the so-called coreshine effect, i.e., enhanced mid-IR light scattering detected with Spitzer toward a number of dense cores Steinacker et al., 2010), implying the presence of larger (Steinacker et al., 2015) or, taking into account the change in dust optical properties, only moderately larger (Ysard et al., 2016), dust grains.
Stochastic emission by very small grains (a 10 nm) is known to contribute significantly to the 60 µm and 100 µm emission bands. This contribution is estimated to be on the order of 13% at 100 µm, and 45% at 60 µm (e.g. Jones et al., 2013). When the 100 µm band is used to derive the dust temperature 14 , this contamination affects the accuracy of mass determinations using dust continuum measurements. The formation of dust aggregates first removes the very small grains from the gas phase, as suggested by observations showing a significant decrease in the 60 µm emission (Laureijs et al., 1991;Bernard et al., 1999;Stepnik et al., 2003;Ysard et al., 2013) and as predicted by dust evolution models (Ossenkopf & Henning, 1994;Köhler et al., 2015). Because small grains are not aligned with the magnetic field, such contamination is absent from the polarized thermal emission SED 15 . As a consequence, the dust temperature derived 16 from the polarized dust SED that B-BOP can observe will only reflect the temperature of large aligned dust grains in the transition from the diffuse to the dense ISM, a constraint that will be used in addition to that inferred from the total intensity SED to study dust evolution processes.
Just like unpolarized emission, polarized dust emission will also probe variations in dust emissivity as expected from dust evolution in dense clouds. Dense environments could not be properly characterized in polarization at the low resolution of the Planck data. Unlike unpolarized emission, the anisotropic nature of polarized emission makes it sensitive to the grain shape. The formation of dust aggregates by grain-grain coagulation must have its counterpart in polarization, and B-BOP will detect signatures which will have to be analyzed through detailed modeling of the coagulation process. Here again, the combination of B-BOP data with the increasing amount of high-resolution polarization observations in the optical and the near-IR will provide strong constraints on dust evolution (Planck int. res. XXI, 2015; Planck 2018 res. XII, 2019).
Observations of total dust emission intensity at farinfrared and submillimeter wavelengths with Planck and Herschel have also brought surprises. One of them is evidence that the logarithmic slope of the dust emission SED at long wavelengths, often referred to as the dust emissivity index, β, exhibits significant variations at large scales across the Galaxy. The Planck all-sky data clearly show variations of β along the Galactic Plane, from very steep SEDs toward inner regions of the Milky Way to much flatter SEDs (β 1.5) toward the Milky Way anticenter (Planck 2013res. XI, 2014. This has also been confirmed in the far-infrared by the analysis of Hi-GAL data (Paradis et al., 2012). Even larger variations have been found in observations of external galaxies, with the SMC and LMC having β 1.3 and β 1.0, respectively (Planck early res. XVII, 2011). Such variations are observed in the Herschel data within individual nearby galaxies such as M31 (Smith et al., 2012) and M33 (Tabatabaei et al., 2014). The origin of these variations is currently unclear and three main classes of dust models have been proposed to explain them. The first type involves the mixing of different materials during the dust life-cycle Ysard et al., 2015). The second type of models invokes Two-Level-System (TLS) low-energy transitions in the amorphous material composing dust grains (Meny et al., 2007) as the cause for the flattening of the SED. The third type of models proposes that magnetic inclusions in dust grains (Draine & Hensley, 2013) could produce the observed variations (Draine & Hensley, 2012). Determining the origin of these variations is critical in many respects, not only to understand the dust cycle in the ISM, but also for accurate mass determinations from dust continuum measurements (which require good knowledge of the dust emissivity, its wavelength dependence, and its spatial variations).
In this domain again, extensive polarimetric imaging at far-IR wavelengths is likely to play a critical role in the future. Dust models based on dust evolution have not yet presented their predictions in polarization, but the other two classes of models mentioned above predict significantly different behaviors for the polarization fraction as a function of wavelength. TLS-based models essentially predict a flat spectrum for the polarization fraction, a prediction compatible with Planck observations. In contrast, metallic-inclusion models predict variations of the polarization fraction in the submillimeter (Draine & Hensley, 2013), which are not observed. It will be possible to evidence those distortions of the polarization SED by comparing far-IR B-BOP observations with existing submillimeter Planck data for the Magellanic clouds and polarization data obtained with new ground-based polarimetric facilities such as NIKA2-POL and SCUBA2-POL for Milky Way regions/sources and nearby galaxies. Correlating changes in the polarized SED with variations of β will allow us to constrain models of the submillimeter dust emissivity in a very unique way.

Toward a tentative detection of polarization by dust self-scattering in the densest cores
In the past decade, interesting constraints on grain sizes in dense clouds have come from the detection of the "coreshine" effect , which results from scattering of near-IR stellar photons by dust grains present in the cloud. This has been interpreted as evidence of grain growth (a = 1 µm, Steinacker et al., 2010) or grain compositional and structural evolution with a modest size increase (a < 0.5 µm, Ysard et al., 2016). More recently, it has been demonstrated that very large (a > 10 µm) dust grains are able to produce polarization by scattering thermal dust emission, a process that is called 'self-scattering'. This was first predicted to be observable in protoplanetary disks (Kataoka et al., 2015) and then confirmed by numerous ALMA observations (Kataoka et al., 2016;Yang et al., 2017;Girart et al., 2018). B-BOP, with its 100 µm polarized channel, would be able to detect and characterize the spectral dependence of polarization by scattering due to ∼ 15 µm dust grains (Kataoka et al., 2015), if present. For the effect to be observable, the thermal emission of dust must first present a quadrupolar anisotropy: scattering grains must receive more far-IR irradiation along one direction in the plane of the sky than along the perpendicular direction. Such a condition is naturally met in a dense protostellar core, or in the presence of density gradients. Second, high local densities (> 10 6 cm −3 ) must be present along the line of sight so that dust grains can have grown to the very large sizes (a ∼ λ/2π) needed for scattering to occur in the far-IR. Observing such a high-density medium should be feasible at 100 µm at the resolution of SPICA, but simulations of grain growth and polarization by scattering are needed to confirm this idea.
Polarization by self-scattering at 100 µm with B-BOP will most likely concern only a few lines of sight through the densest cores. Because polarization due to scattering sharply declines at wavelengths larger than the grain size, it will not alter the polarized emission from aligned dust grains at 200 µm and 350 µm used to trace the local magnetic field orientation in molecular clouds ( § 2 to §5).

MOLECULAR CLOUDS AND THE ORIGIN OF COSMIC RAYS
Polarimetric imaging of molecular clouds (MCs) 17 with B-BOP will also be very useful for the study of the origin of cosmic rays (CRs). Indeed, CRs pervade the whole galaxy, and their interaction with the dense gas of MCs has two important consequences. First, the interactions of high energy CRs (kinetic energy larger than a few hundred MeV) with the gas make MCs bright γ-ray sources. Second, CRs of low energy (≈ 1 − 100 MeV) are the only ionizing agents able to penetrate MCs and regulate the ionization fraction of MC dense cores. For these reasons, observations of enhanced levels of γ-ray emission or ionization rates from MCs reveal the presence of a CR accelerator in their vicinity (see review by Gabici & Montmerle, 2015). A first problem arises, that of the propagation of cosmic rays inside molecular clouds, for which the strength, and above all, the topology of the magnetic field around and inside them plays a central role. At high energies, CRs are unaffected by the magnetic field, so that they interact with all the gas (atomic as well as molecular): the γ-ray emissivity is simply proportional to the product (CR flux × cloud mass). In other words, for a given cloud mass, determining the γ-ray luminosity of a molecular cloud allows the local CR flux to be measured, irrespective of the magnetic field. Over galactic scales, it is well established that the CR flux deduced in this way is essentially uniform (Ackermann et al., 2011), which means that the CR diffusion away from their sources is efficient enough not to be sensitive to large-scale spatial features like the spiral arms. However, there may be γ-ray "hot spots" close to CR sources, and this is precisely what happens when a supernova remnant (SNR) collides or is located in the vicinity of a molecular cloud. Many examples of such SNR-MC associations are known (Gabici & Montmerle, 2015). The reason for the enhanced γray emission is that SNR shock waves accelerate CR in situ, via the so-called Diffuse Shock Acceleration, or DSA, mechanism (Drury, 1983), so this particular configuration can be considered as a "CR laboratory": If the shock-accelerated CRs are insensitive to the local magnetic fields when they reach high energies, the process by which they do, in other words the acceleration mechanism itself depends very much on it, and again on its topology close to the SNR shock. For instance, it is well known from theoretical DSA models that the acceleration efficiency depends strongly on the angle between the shock front and the local magnetic field lines (Caprioli & Spitkovsky, 2014). So, again, knowing the magnetic field topology on small spatial scales in MCs impacted by SNRs would allow a detailed study of the CR acceleration process in the vicinity of the shock front.
In particular, a new picture of the ISM in star-forming regions has now to be taken into account: in the standard picture summarized above, the gaseous medium in which the SNR shock propagates is assumed uniform. But as recent work has shown (see Sect. 2.2 and Fig. 1), the structure of MCs is not uniform, but filamentary, down to scales of parsecs (in length) and ∼ 0.1 pc (in width), i.e., precisely those that are accessible to B-BOP at distances ≈ 1 kpc. One of the key changes is then that the shock would cross the ambient magnetic field lines at all angles, and likely perturb them and enhance the turbulent component of the magnetic field, on spatial scales comparable to that of the filaments: for not-toodistant sources, B-BOP would then act as a "magnifyingglass" to study the shock-ambient gas interactions at unprecedented small spatial scales, and put entirely new, perhaps even unexpected, constraints on DSA models (e.g., see the CR escape issues raised by Malkov et al. 2013).
At low energies, the situation is markedly different, because the transport properties of CRs become very sensitive to a variety of processes governed by the magnetic field properties which may hamper the penetration of CRs into MCs, and reduce the rate at which the gas is ionized by these particles (Phan et al., 2018). Like γ-ray production at high CR energies, the ionization by low-energy CR can also be measured by way of infrared and millimeter-wave observations, which detect lines of various molecules and radicals (like H + 3 , HCO + , DCO + , etc. -cf. Indriolo et al., 2015). This has been done for many MCs in the Galaxy, but more recently also for SNR-MC collision regions (Gabici & Montmerle, 2015): here again an enhancement of MC ionization has been observed. The results tentatively suggest a proportionality between the SNR-accelerated high-energy and low-energy CR fluxes, constraining the acceleration mechanism, or a proportionality between the partially irradiated, ionized gas, and the fully irradiated, γ-ray emitting gas, or both.
More generally, both low-and high-energy CRs are affected in their propagation in the diffuse ISM by diffusion effects, which are still poorly known -and directly influenced by magnetic fields. The spatially average diffusion coefficient of CRs in the Galaxy is constrained by a number of observations, and is often assumed to be isotropic on large Galactic scales ( 100 pc). On the other hand, in order to explain a number of γ-ray observations of SNR-MC associations (characterized by spatial scales of ∼ 10 − 100 pc), a diffusion coefficient about two orders of magnitudes smaller (i.e., slower bulk propagation) than the average Galactic one needs to be assumed. However, such a discrepancy could be reconciled if CR diffusion is in fact anisotropic on such small scales (Nava & Gabici, 2013). An anisotropic diffusion is indeed expected for spatial scales smaller than the magnetic field coherence length (Malkov et al., 2013). Knowing the topology of the magnetic field in such regions is thus of paramount importance in order to interpret γ-ray observations correctly.
All of these issues can be illustrated by a recent study of the W28 SNR (cf. Fig. 14), a well-known example of an SNR-MC collision (Vaupré et al., 2014). This SNR is located in the Galactic plane, at d ≈ 2 kpc from the Figure 14. The region surrounding the W28 SNR (d ∼ 2 kpc; shock approximated by the dashed white circle), as seen in cold dust emission at 353 GHz by Planck (color image with background B-field "drapery" from polarization data), in TeV-GeV γ-rays (white contours), and CO (green areas, well correlated with the γ-ray sources - Aharonian et al., 2008,). The labels highlight the various CR processes discussed in the text, at high energies (HECR) and low energies (LECR). The Planck and B-BOP beams are indicated by a light green circle (label "P") and a dot (label "S"), respectively. Sun, with an estimated (very uncertain) age ≈ 10 4 yr. At this distance, the SNR apparent diameter (∼ 30 ) gives a linear diameter D ≈ 20 pc. An observation by the High Energy Stereoscopic System (HESS) Čerenkov telescope, in the TeV γ-ray range (Aharonian et al., 2008), covering a large field-of-view of ∼ 1.5 • × 1.5 • (with a resolution of ∼ 0.1 • ), has revealed a complex of several resolved γ-ray sources. One of the sources, which is spatially correlated with a part of the SNR shock outline, was also detected as a bright GeV source by the Fermi satellite, contrary to the other sources, which are either dimmer or undetected (Abdo et al., 2010). This multiple source was soon correlated with a complex of molecular clouds mapped in CO by the NANTEN 18 telescope, showing that the SNR was in fact colliding with the molecular cloud associated with the GeV-TeV source (a physical contact being confirmed by the existence of several OH masers), the other sources being away, far upstream of the SNR shock.
Calculations indicated a factor ≈ 100 enhancement of the local high-energy CR flux, qualitatively consistent with a local CR acceleration by the SNR shock. Using the IRAM 30-m telescope to observe various molecular and radical tracers (H 13 CO + , C 18 O, etc.) in the millimeter range, Vaupré et al. (2014) were also able to calculate the MC ionization rate ζ at several locations. They found ∼ 2 − 3 order of magnitude enhancements of ζ (or lower 18 NANTEN means "Southern Sky" in Japanese. limits) over its average Galactic value (ζ 0 ≈ 4 − 5 × 10 −17 erg s −1 -e.g. Indriolo et al., 2015), for the cloud correlated with GeV-TeV emission, i.e., indirect evidence for a similar enhancement of the low-energy CR flux, but no such enhancements for the clouds far upstream of the SNR shock. Within the "GeV-TeV bright" cloud, the measurements were separated by the IRAM telescope resolution, ∼ 12 , i.e., comparable to (only 1.5 times better than) the B-BOP resolution at 200 µm (or a linear scale ∼ 0.15 pc). The (projected) distance to the other clouds is ∼ 10 pc, and this is seen as the diffusion length for high-energy CRs (TeV CRs reaching the distant clouds before the GeV CRs).
Thus, the W28 SNR and its environment provide us with a case study with all the ingredients needed to improve our understanding of the origin of CRs, and their relation with magnetic fields down to scales ∼ 0.15 pc, i.e., much smaller than observable before: (i) CR acceleration by the SNR shock; (ii) diffusion of CRs between clouds as a function of energy; (iii) penetration of low-energy CRs in "average" clouds, irradiated only by ambient, galactic CRs (cf. Fig. 14).
For more distant sources, the "magnifying-glass" effect of B-BOP on small spatial scales would of course decrease, but an interesting link could then be established with the Čerenkov Telescope Array (CTA, presently under construction; Actis et al. 2011), which operates in the 20 GeV-300 TeV regime. Until CTA actually observes, it is difficult to make accurate predictions on how far SNR-MC systems such as W28 could be detected in γ-rays. For CR studies, the main point is not simply the detection, but the location of the emission with respect to the shock. For the moment, only two cases are known in which the γ-rays are clearly upstream of the shock: W28 (detected by CGRO, Fermi and HESS, so from GeV to TeV energies), and W44 (detected by Fermi only, so at GeV energies only). About 20% of HESS sources, and most of the SNRs detected by Fermi are SNR-MC systems (H.E.S.S. Collab. et al. 2018;Ackermann et al. 2011), but apart from W28 and W44 it is difficult to distinguish between upstream and downstream γ-ray emissions (or both). Taking into account that CTA is ∼ 10 times more sensitive than HESS, and taking W28 as a template, representative of a Galactic disk SNR-MC population, we estimate that ≈ ten W28-like sources could be possibly detected and be sufficiently resolved by CTA, hence be good targets for future B-BOP observations, which would in turn spur new theoretical work on CR acceleration on very small spatial scales not considered at present.

POLARIZED DUST EMISSION FROM PROTOPLANETARY DISKS
As already mentioned in Sect. 2, magnetic fields may regulate the gravitational collapse and fragmentation of prestellar dense cores, thereby influencing the overall star formation efficiency (Mouschovias & Ciolek, 1999;Dullemond et al., 2007;Crutcher, 2012). It is thus natural to expect that, during core collapse, magnetic fields can be dragged inward, leaving a remnant field in the protoplanetary disk formed subsequently. If protoplanetary disks are indeed (weakly) magnetized, then the MHD turbulence arising from magnetorotational instability (MRI) is thought to be the primary source of disk viscosity, a crucial driving force for disk evolution (e.g. disk accretion) and planet formation (Balbus & Hawley, 1998;Turner et al., 2014). Despite this general consensus, our knowledge about magnetic fields in disks is actually very limited and incomplete at this stage, largely due to the lack of observational constraints on magnetic field properties (geometry and strength) in protoplanetary disks. Polarimetric observations of thermal dust emission at centimeter or millimeter wavelengths with single-dish telescopes, such as the Caltech Submillimeter Observatory (CSO) and JCMT, or interferometric arrays, such as the Karl G. Jansky Very Large Array (JVLA), the SMA, the Berkeley-Illinois-Maryland Array (BIMA), and the Combined Array for Research in Millimeter-wave Astronomy (CARMA), have been used extensively to map magnetic field structure in YSOs at scales from ∼ 50 AU to thousands of AU (see Crutcher, 2012 for a review). However, due to the limited sensitivity and angular resolution offered by current facilities and the nature of centimeter/millimeter observations, most of these studies have been focused on magnetic fields in molecular clumps/cores (cf. § 4), or Class 0/Class I objects (e.g. Qiu et al., 2013;Rao et al., 2014;Liu et al., 2016), rather than classical protoplanetary disks around Class II objects. Using CARMA, Stephens et al. (2014) spatially resolved the HL Tau disk in polarized light at 1.3 mm, and their best-fit model suggested that the observation was consistent with a highly tilted (by ∼50 • from the disk plane), toroidal magnetic field threading the disk. This conclusion was challenged by follow-up studies, which showed that the 1.3 mm polarization of HL Tau could also arise solely from dust scattering as opposed to dichroic emission from elongated grains aligned with the magnetic field (Kataoka et al., 2015;Yang et al., 2016). However, more recent ALMA observations at 0.87, 1.3, and 3.1 mm indicated that dust scattering alone may not be able to explain all of the multi-wavelength polarization data (Kataoka et al., 2017;Stephens et al., 2017).
Recently, Li et al. (2016) have been able to highlight the signature of a magnetic field in the AB Aur protoplanetary disk at mid-IR wavelengths. Using observations of the AB Aur protoplanetary disk at 10.5 µm with the GTC/CanariCam imager and polarimeter, they detected a polarization pattern in the inner regions of the disk compatible with dichroic emission polarization produced by elongated grains aligned by a tilted poloidal magnetic field. The observed polarization level (2-3 %) was somewhat lower than that predicted by theory (Cho & Lazarian, 2007), although this is something naturally expected since the modeling assumes alignment efficiencies and intrinsic particles polarizability which are probably overestimated. At a wavelength of 10.5 µm where protoplanetary disks are optically thick, the observations are probing the disk properties down to depths corresponding to the τ =1 optical depth surface. This depth is relatively small (less than ∼10%) compared to the disk scale height. Longer wavelengths, up to about 200 µm where the disk can still be moderately optically thick up to large distances from the star, are emitted by cooler material located deeper within the disk. Thus, by measuring the level of polarization at wavelengths in the range 100-300 µm, we expect to be able to com-pare the levels of polarization as a function of optical thickness, thereby getting an indirect signature of the magnetic field at different depths within protoplanetary disks. This will provide constraints on the importance of MRI-induced turbulence. Together with 10 µm and ALMA similar types of observations, this will allow us to build a tomographic view of the magnetic field along the vertical profile of the disks. Such measurements would also have large impacts on our understanding of planet formation processes.
For the purpose of this paper, we used the same modeling approach as described in Li et al. (2016), with also the same disk parameters based on the example of AB Aur, in order to predict expected polarization levels levels at 100 µm, 200 µm, and 350 µm if observed by SPICA. Given the angular resolution of SPICA at these wavelengths, we do not expect, apart from exceptional cases, to angularly resolve the polarized emission, therefore we computed an integrated value, considering the object as unresolved.
In Fig. 15, we show the predictions of the model, integrated over the full spatial extent of the disk, for different disk inclinations with respect to the line of sight and for observing wavelengths within the B-BOP range. Two types of magnetic field configuration are considered, poloidal and toroidal, which are the simplest ones, and those also widely discussed in the literature. We can see from the simulations that the poloidal magnetic field configuration produces stronger integrated polarization signatures compared to the toroidal configuration.
As mentioned earlier, the origin of the dust continuum polarization on the disk scale is still uncertain, with potential contributions from scattering by large grains in addition to that from emission by magnetically (or radiatively) aligned grains. The B-BOP instrument will generally not (or barely be able to) resolve protoplanetary disks, which poses the problem of disentangling these various mechanisms. Fortunately, polarization by aligned dust grains and polarization by dust self-scattering have different dependences on wavelength and optical depth (Yang et al., 2017). SPICA will greatly extend the wavelength coverage of ALMA (from 870 µm to 100 µm), which will make it easier to disentangle the contributions from the competing mechanisms. Such an effort is a pre-requisite for using dust continuum polarization to probe both disk magnetic fields and grain growth, the crucial first step toward the formation of planetesimals and ultimately planets. Moreover, given the plan to image the whole extent of nearby star-forming regions with B-BOP (cf. end of § 2.3), several tens of protoplanetary disks will be detected in Stokes I, Q, U. It will therefore be possible to derive statistical trends about the presence of magnetic fields, and any bias can be controlled providing that the inclination and position angles of the disks are known.

VARIABILITY STUDIES OF PROTOSTARS IN THE FAR-INFRARED
At the core (≤ 0.1 pc) scale, the formation of a solar-type star is well understood as a continual mass assembly process whereby material in the protostellar envelope is accreted onto a circumstellar disk and then transported inward and onto the protostar via accretion columns (Hartmann et al., 2016). The observational evidence for the mass assembly rate of low-mass stars is provided by the lifetimes of the various stages and the bolometric luminosities, which are dominated by accretion energy at early times (e.g. Dunham et al., 2014, and references therein). These two quantifiable measures are in significant disagreement and circumstantial evidence exists for the episodic nature of mass assembly -bullets in outflows (e.g. Plunkett et al., 2015), FU and EX Ori phenomenon (Hartmann & Kenyon, 1996;Herbig, 2008), numerical calculations of disk transport (e.g. Armitage, 2015). Moreover, high spatial resolution images of young disks reveal macroscopic structure including rings (ALMA Partnership et al., 2015), spirals (Pérez et al., 2016), and fragmentation , suggesting that the transport of material through the disk is not a smooth and steady process.
Recently, significant variability has been detected in the submillimeter continuum emission of several nearby protostars (see, e.g., Mairs et al., 2017b;Yoo et al., 2017;Johnstone et al., 2018a), through an ongoing multi-year monitoring survey of eight nearby star-forming regions with JCMT at 850 µm (Herczeg et al., 2017). While uncertain due to small number statistics, it appears that roughly 10% of deeply embedded protostars vary over year timescales by around 10% at submillimeter wavelengths (Johnstone et al., 2018a). The dominant mode of variability uncovered by the survey is quasi-secular, with the protostellar brightness increasing or decreasing for extended -multi-year -periods (Mairs et al., 2017b;Johnstone et al., 2018a) and suggesting a link to non-steady accretion processes taking place within the circumstellar disk where the orbital timescales match those of the observed variability. These long timescales also allow for significant amplification of the overall change in submillimeter brightness after many years. One source, EC 53 in Serpens Main, has an eighteen-month quasi-periodic light curve (Yoo et al., 2017, and Fig. 16), previously identified through near-IR observations (Hodapp et al., 2012), which is likely due to periodic forcing by a longlived structure within the inner several AU region of the disk.
Stronger variability is expected at far-IR wavelengths where Class 0 and Class I YSOs have the peak of their SEDs and the envelope emission directly scales with the internal (accretion) luminosity of the underlying protostar (Dunham et al., 2008;Johnstone et al., 2013). At The typical uncertainty in a single measurement is ∼ 20 mJy (S/N ∼ 50) and the peak to peak brightness variation is almost 500 mJy. submillimeter wavelengths the emission typically scales with the envelope temperature and thus the far-IR signal is expected to be around four times larger such that a 10% variability in the submillimeter relates to a 40% variability in the far-IR. Contemporaneous monitoring of EC 53 in the near-IR and at 850/450 µm has confirmed that the longer submillimeter wavelength shows just such a diminished response to the underlying change in internal luminosity as proxied by near-IR observations (Yoo et al., in prep.). Thus, carefully calibrated monitoring of protostars with SPICA should uncover a significantly larger fraction of variables than the 10% obtained by the JCMT survey. Thanks to its high continuum sensitivity and mapping speed at wavelengths around the peak of protostellar SEDs, B-BOP, used as a total-power imager, will be ideal for monitoring hundreds of forming stars over multi-year epochs, allowing an unprecedented statistical determination of the variation in accretion on these timescales. Typical nearby deeply embedded protostars have far-IR brightnesses greater than ∼ 10 mJy and thus will be observed to a S/N > ∼ 100 by B-BOP in a fast scanning mode. As demonstrated for ground-based submillimeter observations (Mairs et al., 2017a), instrument stability will need to be carefully monitored in order to achieve precise relative flux calibration of a few percent between epochs. Additional critical requirements for B-BOP will be a large, 10 5 or higher, dynamic range and instrument robustness against extremely bright sources within the field.
Three interconnected monitoring surveys are envisioned. First, the bulk of the ∼ 1000 nearby, Gould Belt, deeply embedded protostars will be observed every six months while SPICA is in orbit, requiring coverage of ∼ 20 deg 2 (a modest twenty hours of observing per epoch). This will allow for a detailed statistical characterization of variability across multiple years. Additionally, a few carefully chosen nearby star formation fields, each roughly a square degree, will be observed weekly during their expected few month continuous observing window (for information on observing strategies for SPICA, see Roelfsema et al., 2018). For both of these nearby samples, an even larger number of Class II YSOs will be observable within each field. While these sources will be fainter at far-IR wavelengths as the emission probes the disks directly, the enhanced numbers will allow for a determination of the importance of variability throughout the evolution of a protostar. Finally, a sample of more distant high-mass star-forming regions should be observed yearly to search for rare, but extremely bright, bursts such as FU Ori events. While SPICA will not have the spatial resolution to separate individual protostars within these regions, evidence of a significant brightening can be easily followed-up with ground-based telescopes such as ALMA (see Hunter et al., 2017 for an example of a brightening in a high-mass star-forming region).

CONCLUDING REMARKS
Magnetic fields are a largely unexplored "dimension" of the cold Universe. While they are believed to be a key dark ingredient of the star formation process through most of Cosmic time, they remain very poorly constrained observationally, especially in the cold ISM of galaxies (e.g. Crutcher, 2012).
Benefiting from a cryogenic telescope, SPICA-POL or B-BOP will be two to three orders of magnitude more sensitive than existing or planned far-IR/submillimeter polarimeters (cf. Fig. 5) and will therefore lead to a quantum step forward in the area of far-IR dust polarimetric imaging, one of the prime observational techniques to probe the topology of magnetic fields in cold, mostly neutral environments. In particular, systematic polarimetric imaging surveys of Galactic molecular clouds and nearby galaxies with B-BOP have the potential to revolutionize our understanding of the origin and role of magnetic fields in the cold ISM of Milky-Way-like galaxies on scales from ∼ 0.01 pc to a few kpc. The three main science drivers for B-BOP are 1) probing how magnetic fields control the formation, evolution, and fragmentation of dusty molecular filaments (Sect. 2), thereby setting the initial conditions for individual protostellar collapse (Sects. 4 and 5); 2) characterizing the structure of both turbulent and regular magnetic fields in the cold ISM of nearby galaxies, including the Milky Way, and constraining galactic dynamo models (Sect. 3 and Sect. 6); and 3) testing models of dust grain alignment and informing dust physics (Sect. 7). Other science areas can be tackled with, or uniquely informed by, B-BOP observations, including the problem of the interaction of cosmic rays with molecular clouds (Sect. 8), the study of the magnetization of protoplanetary disks (Sect. 9), and the characterization of variable accretion in embedded protostars (Sect. 10). Last, but not least, the leap forward provided by B-BOP in far-IR imaging polarimetry will undoubtedly lead to unexpected discoveries, such as the potential detection of polarization from the CIB (Sect. 6.3).

ACKNOWLEDGEMENTS
This paper is dedicated to the memory of Bruce Swinyard, who initiated the SPICA project in Europe, but unfortunately died on 2015 May 22 at the age of 52. He was ISO-LWS calibration scientist, Herschel-SPIRE instrument scientist, first European PI of SPICA and first design lead of SAFARI. This work has received support from the European Research