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.

The atomic interstellar medium (ISM) is observed to be full of linear structures that are referred to as "fibers". Fibers exhibit similar properties to linear structures found in molecular clouds known as striations. Suggestive of a similar formation mechanism, both striations and fibers appear to be ordered, quasi-periodic, and well-aligned with the magnetic field. The prevailing formation mechanism for striations involves the excitation of fast magnetosonic waves. Based on this theoretical model, and through a combination of velocity centroids and column density maps, Tritsis et al. (2018) developed a method for estimating the plane-of-sky (POS) magnetic field from molecular cloud striations. We apply this method in two H\textsc{I} clouds with fibers along the same line-of-sight (LOS) towards the ultra-high-energy cosmic-ray (UHECR) hotspot, at the boundaries of Ursa Major. For the cloud located closer to Earth, where Zeeman observations from the literature were also available, we find general agreement in the distributions of the LOS and POS components of the magnetic field. We find relatively large values for the total magnetic field (ranging from $\sim$$\rm{10}$ to $\sim$$\rm{20} ~\rm{\mu G}$) and an average projection angle with respect to the LOS of $\sim$ 50$^\circ$. For the cloud located further away, we also find a large value for the POS component of the magnetic field of $15^{+8}_{-3}~\rm{\mu G}$. We discuss the potential of our new magnetic-field tomography method for large-scale application. We consider the implications of our findings for the accuracy of current reconstructions of the Galactic magnetic field and on the propagation of UHECR through the ISM.

Recent studies of the diffuse parts of molecular clouds have revealed the presence of parallel, ordered low-density filaments termed striations. Flows along magnetic field lines, Kelvin-Helmholtz instabilities and hydromagnetic waves are amongst the various formation mechanisms proposed. Through a synergy of observational, numerical and theoretical analysis, previous studies singled out the hydromagnetic waves model as the only one that can account for the observed properties of striations. Based on the predictions of that model, we develop here a method for measuring the temporal evolution of striations through a combination of molecular and dust continuum observations. Our method allows us to not only probe temporal variations in molecular clouds but also estimate the strength of both the ordered and fluctuating components of the magnetic field projected on the plane-of-the-sky. We benchmark our new method against chemical and radiative transfer effects through two-dimensional magnetohydrodynamic simulations coupled with non-equilibrium chemical modelling and non-local thermodynamic equilibrium line radiative transfer. We find good agreement between theoretical predictions, simulations and observations of striations in the Taurus molecular cloud. We find a value of $\rm{27 \pm 7} ~\rm{\mu G}$ for the plane-of-sky magnetic field, in agreement with previous estimates via the Davis-Chandrasekhar-Fermi method, and a ratio of fluctuating to ordered component of the magnetic field of $\sim$ 10\%.

We describe PyRaTE, a new, non-local thermodynamic equilibrium (non-LTE) line radiative transfer code developed specifically for post-processing astrochemical simulations. Population densities are estimated using the escape probability method. When computing the escape probability, the optical depth is calculated towards all directions with density, molecular abundance, temperature and velocity variations all taken into account. A very easy-to-use interface, capable of importing data from simulations outputs performed with all major astrophysical codes, is also developed. The code is written in Python using an `embarrassingly parallel' strategy and can handle all geometries and projection angles. We benchmark the code by comparing our results with those from RADEX (van der Tak et al. 2007) and against analytical solutions and present case studies using hydrochemical simulations. The code is available on GitHub (https://github.com/ArisTr/PyRaTE).

Stars and planets are formed inside dense interstellar molecular clouds, by processes imprinted on the 3-dimensional (3D) morphology of the clouds. Determining the 3D structure of interstellar clouds remains challenging, due to projection effects and difficulties measuring their extent along the line of sight. We report the detection of normal vibrational modes in the isolated interstellar cloud Musca, allowing determination of the 3D physical dimensions of the cloud. Musca is found to be vibrating globally, with the characteristic modes of a sheet viewed edge-on, not a filament as previously supposed. We reconstruct the physical properties of Musca through 3D magnetohydrodynamic simulations, reproducing the observed normal modes and confirming a sheet-like morphology.

Dust continuum and molecular observations of the low column density parts of molecular clouds have revealed the presence of elongated structures which appear to be well aligned with the magnetic field. These so-called striations are usually assumed to be streams that flow towards or away from denser regions. We perform ideal magnetohydrodynamic (MHD) simulations adopting four models that could account for the formation of such structures. In the first two models striations are created by velocity gradients between ambient, parallel streamlines along magnetic field lines. In the third model striations are formed as a result of a Kelvin-Helmholtz instability perpendicular to field lines. Finally, in the fourth model striations are formed from the nonlinear coupling of MHD waves due to density inhomogeneities. We assess the validity of each scenario by comparing the results from our simulations with previous observational studies and results obtained from the analysis of CO (J = 1 - 0) observations from the Taurus molecular cloud. We find that the first three models cannot reproduce the density contrast and the properties of the spatial power spectrum of a perpendicular cut to the long axes of striations. We conclude that the nonlinear coupling of MHD waves is the most probable formation mechanism of striations.

We present a new method for assessing the intrinsic 3D shape of prestellar cores from molecular column densities. We have employed hydrodynamic simulations of contracting, isothermal cores considering three intrinsic geometries: spherical, cylindrical/filamentary and disk-like. We have coupled our hydrodynamic simulations with non-equilibrium chemistry. We find that a) when cores are observed very elongated (i.e. for aspect ratios $\le$ 0.15) the intrinsic 3D geometry can be probed by their 2D molecular emission maps, since these exhibit significant qualitative morphological differences between cylindrical and disk-like cores. Specifically, if a disk-like core is observed as a filamentary object in dust emission, then it will be observed as two parallel filaments in $\rm{N_2H^{+}}$; b) for cores with higher aspect ratios (i.e. 0.15 $\sim$ 0.9) we define a metric $\Delta$ that quantifies whether a molecular column density profile is centrally peaked, depressed or flat. We have identified one molecule ($\rm{CN}$) for which $\Delta$ as a function of the aspect ratio probes the 3D geometry of the core; and c) for cores with almost circular projections (i.e. for aspect ratios $\sim$ 1), we have identified three molecules ($\rm{OH}$, $\rm{CO}$ and $\rm{H_2CO}$) that can be used to probe the intrinsic 3D shape by close inspection of their molecular column density radial profiles. We alter the temperature and the cosmic-ray ionization rate and demonstrate that our method is robust against the choice of parameters.

We revisit the relation between magnetic-field strength ($B$) and gas density ($\rho$) for contracting interstellar clouds and fragments (or, cores), which is central in observationally determining the dynamical importance of magnetic fields in cloud evolution and star formation. Recently, it has been claimed that a relation $B \propto \rho^{2/3} $ is statistically preferred over $B \propto \rho^{1/2}$ in molecular clouds, when magnetic field detections and nondetections from Zeeman observations are combined. This finding has unique observational implications on cloud and core geometry: The relation $B \propto \rho^{2/3} $ can only be realized under spherical contraction. However, no indication of spherical geometry can be found for the objects used in the original statistical analysis of the $B-\rho$ relation. We trace the origin of the inconsistency to simplifying assumptions in the statistical model used to arrive at the $B\propto \rho^{2/3}$ conclusion and to an underestimate of observational uncertainties in the determination of cloud and core densities. We show that, when these restrictive assumptions are relaxed, $B \propto \rho^{1/2}$ is the preferred relation for the (self-gravitating) molecular-cloud data, as theoretically predicted four decades ago.