Polarization power spectra and dust cloud morphology

In the framework of studies of the CMB polarization and its Galactic foregrounds, the angular power spectra of thermal dust polarization maps have revealed an intriguing E/B asymmetry and a positive TE correlation. In interpretation studies of these observations, magnetized ISM dust clouds have been treated as filamentary structures only; however, sheet-like shapes are also supported by observational and theoretical evidence. In this work, we study the influence of cloud shape and its connection to the local magnetic field on angular power spectra of thermal dust polarization maps. We simulate realistic filament-like and sheet-like interstellar clouds, and generate synthetic maps of their thermal dust polarized emission using the software $Asterion$. We compute their polarization power spectra in multipole range $\ell \in [100,500]$ and quantify the E/B power asymmetry through the $R_{EB}$ ratio, and the correlation coefficient $r^{TE}$ between T and E modes. We quantify the dependence of $R_{EB}$ and $r^{TE}$ values on the offset angle (between longest cloud axis and magnetic field) and inclination angle (between line-of-sight and magnetic field) for both cloud shapes embedded either in a regular or a turbulent magnetic field. We find that both cloud shapes cover the same regions of the ($R_{EB}$, $r^{TE}$) parameter space. The dependence on inclination and offset angles are similar for both shapes although sheet-like structures generally show larger scatter. In addition to the known dependence on the offset angle, we find a strong dependence of $R_{EB}$ and $r^{TE}$ on the inclination angle. The fact that filament-like and sheet-like structures may lead to polarization power spectra with similar ($R_{EB}$, $r^{TE}$) values complicates their interpretation. In future analyses, this degeneracy should be accounted for as well as the connection to the magnetic field geometry.


Introduction
The study of polarization of the cosmic microwave background (CMB) is very important for the field of cosmology as it is the only foreseeable way to test observationally whether the Universe has undergone an early period of inflationary expansion.Indeed, gravitational waves generated during the inflation era would have left unambiguous imprints on CMB polarization.These imprints are the so-called B-modes in the polarization power spectra which, although difficult to detect, are in principle observable (Polnarev 1985).
Most of the cosmological information can be obtained from the CMB radiation through the analysis of its T , E, and B angular power spectra (intensity and linear polarization) and their correlation (e.g., Kamionkowski et al. 1997;Seljak & Zaldarriaga 1997;Hu & White 1997) in which the polarization is decomposed into a scalar, curl-free field, the E mode, and a pseudoscalar, curl field, the B mode. akonstantinou@physics.uoc.grpelgrims@physics.uoc.grGalactic dust, through which we observe the CMB, contaminates the cosmological signal as the grain thermal emission is also polarized and contains B-mode components as well.Nonspherical dust grains align their short axis with the magnetic field of the Galaxy which is present everywhere.As a result, they emit thermal radiation with their electric vector preferentially aligned with their long axis (e.g., Andersson et al. 2015).The polarized emission from Galactic dust grains, with polarization fraction about 20% (e.g., Benoît et al. 2004;Planck Collaboration XII 2020), higher than that of the CMB, which is about 5% (Bond & Efstathiou 1987;Seljak 1997;Hinshaw et al. 2003), always dominates the measured signal at frequency higher than ≈ 80 GHz (e.g., Planck Collaboration Int.XXX 2016).The B-mode component of Galactic thermal dust emission can be confused with the primordial signal of inflation (e.g., BICEP2 Collaboration et al. 2014;Planck Collaboration Int. XXX 2016;Planck Collaboration XI 2020).It is therefore essential to understand the emission and absorption properties of interstellar dust so as to remove its contribution to the map of the polarized primordial microwave sky.
In this context, the dust polarization is characterized by its T , E, and B representations rather than through its representation in Stokes parameters I, Q, and U. E and B are invariant under rotations making the E-B decomposition more convenient since it is unnecessary to define a reference frame.
The characterization of the dust polarization auto-and crossangular power spectra has led to the intriguing evidences, first reported in (Planck Collaboration Int. XXXVIII 2016), for an E/B asymmetry (R EB ∼ 2) and a positive T E correlation (r T E ∼ 0.36) that were not anticipated, at least for sky regions at intermediate and high Galactic latitudes in which a weak signal with random polarization orientations was expected from the dust.Dedicated analysis of the data has revealed that this -a priori unexpected-signal may likely originate from the observed correlation between the orientation of ridges seen in intensity maps and the sky-projected orientation of the integrated magnetic field (Planck Collaboration Int. XXXII 2016;Planck Collaboration Int. XXXV 2016;Planck Collaboration Int. XXXVIII 2016;Rotti & Huffenberger 2019;Soler 2019).Such a correlation between ISM clouds and magnetic field was known for some time for molecular clouds (e.g., Goodman et al. (1990)) and subsequently observed for clouds of the diffuse ISM seen in H I emission (McClure-Griffiths et al. 2006;Clark et al. 2014).This observed correlation is likely the result of the coupling of interstellar matter and the ambient magnetic field, which is expected in the framework of the magneto-hydrodynamic (MHD) description of ISM and it has been extensively discussed in the literature (e.g., Hennebelle 2013;Soler et al. 2013;Caldwell et al. 2017;Kandel et al. 2018;Kim et al. 2019).
To make progress in the characterization of Galactic foregrounds to the CMB polarization, it is necessary to consider the sky and the ISM structures as three-dimensional (3D) systems projected onto the sky (e.g., Tassis & Pavlidou 2015;Pelgrims et al. 2021a).In this context, dedicated studies have been carried out to estimate and characterize the result of projection effects of ISM structures to the measured polarization power spectra (e.g., Rotti & Huffenberger 2019;Huffenberger et al. 2020;Hervías-Caimapo et al. 2021).In these studies, the shape of ISM structures have been routinely, if not exclusively, assumed to be filamentary.The authors have indeed considered and modeled ISM clouds as being cylinders or prolate (elongated) spheroids.However, to date, there is no consensus that all ISM clouds have a filament-like morphology.On the contrary, numerical simulations show that turbulent flows tend to stretch and compress the ISM gas into sheet-like and filament-like structures, which appear elongated in column density maps because of projection effects (e.g., Hennebelle 2013).This shape ambiguity is not new and has been at the heart of a long lasting debate in the ISM literature (e.g., Kainulainen et al. 2016vs. Tritsis & Tassis 2018).Dense, molecular, self-gravitating clouds generally appear filamentary on the sky, giving rise to the picture of a cosmic web in which clouds fragment to give rise to star-forming regions at its nodes (e.g., Myers 2009).However, there are several observational pieces of evidence indicating the existence of clouds with sheet-like morphology that are necessary to explain light echoes of supernovae explosions and pulsar scattering data (e.g., Williamson 1972;Suntzeff et al. 1988;Spyromilio et al. 1995;Brisken et al. 2010;Yang et al. 2016).Such sheet-like clouds may result from large scale planar shocks when supernovae induced bubbles interact (Hartmann et al. 2001).In addition, the very fact that filamentary patterns are not all straight but may show large-scale curvature, as suggested in (Planck Collaboration Int.XXXVIII 2016), might also be indicative of sheet-like structures that would, for example, be part of expanding bub-bles, and would appear brighter where viewed tangentially (e.g., Tahani et al. 2022).Hence, while diverse observations of the diffuse ISM -from relatively low-resolution submillimeter observations by Planck to high-resolution infrared observations by Herschel, and of H I observation in the radio waves (e.g., Kalberla & Kerp 2009;Miville-Deschênes et al. 2010;André et al. 2014;Clark et al. 2014)-reveal its filamentary appearance, it might very well be that it is made of clouds having both filament-like and sheet-like shapes.
The preferred shape of ISM clouds in the diffuse ISM, along with their specific relation to the magnetic field orientation, is expected to depend on the exact balance of several factors: magnetic field strength, gravitational energy, and turbulent energy (e.g., Heiles & Crutcher 2005;Crutcher et al. 2010;Soler et al. 2013).However, it is very challenging to determine observationally the regime corresponding to any given ISM region.The characterization of relative orientations between apparent structures and projected magnetic field is an equally difficult task, with results that may also depend on the details of the specific method used in the analysis (e.g., Micelotta et al. 2021).
In this paper, motivated by the fact that the possibility for sheet-like clouds has been disregarded so far in CMBforeground characterization studies that want to account for 3D effects, we explore the possible effects from ISM-cloud morphology on the characterization of polarization power spectra.
Our main goal is to study how the shape of interstellar clouds affects the angular power spectra of emitted polarized radiation from dust.We set up a simulation-based experiment to investigate the degeneracy that cloud shapes might lead to in the properties of observed polarization power spectra.Our toy-models allow us to quantify the effect of the angle between the structuremajor axes and the local magnetic field orientation, as well as the effect of the viewing angle, i.e. the angle that the sight-lines make with the magnetic field lines.
We simulate interstellar clouds with both filamentary and sheet-like structure using the software Asterion that we developed and present in Sect.2, where we also produce synthetic polarization maps and compute their angular power spectra.In Sect. 3 we present the results of our exploration of the effect the two angles have on the power asymmetry between E and B modes, and on the correlation between T and E modes.We discuss our results in Sect. 4 and provide our conclusions in Sect. 5.

Numerics & synthetic data
2.1.Asterion 2.1.1.The software Asterion is a scientific tool to simulate the magnetized ISM of our Galaxy in 3D, including dust clouds and the magnetic field that permeates them, to directly visualize the data in 3D and interact with it into an immersed virtual environment.
Ultimately, we hope that Asterion will assist in the reconstruction of the 3D structure of the magnetized ISM of our Galaxy, a long-standing problem in the quest for the origin of the magnetic field in our Galaxy.Asterion relies on real time 3D visualization techniques with virtual reality capabilities that reach a high degree of performance.It allows us to render the details of the magnetized ISM and enables the user to interact with the simulated environment as it is done in video games.
Asterion is not meant to be another visualization and raytracing tool, with or without radiative transfer implemented, as already existing ones such as Polaris (Reissl et al. 2016).As-terion enables real-time reprocessing of portions of the visualized data at a higher resolution, allowing sub-grid details to be computed and added to the visualized data in real-time.Asterion is implemented in the Real-Time-Engine 'Unreal Engine 4' (UE4) 1 .
The core features of creating and visualizing 3D structures use the parallel-computation paradigm and utilize the benefits of modern Graphics Processing Units (GPU).A significant amount of hardware threads (Shaders) performs calculations independently and presents solutions at once, as 3D textures, available for further computation and visualization.
Within Asterion (currently at the prototype level), the user is immersed in a 3D environment that emulates the Galaxy.The user can fly through the Galactic space and visualize dust density distribution and magnetic field.For the purpose of this work the user can specify, through a set of parameters, the large-scale dust-density distribution model, the large-scale Galactic magnetic field model and, clouds of dust with non-trivial morphology following some orientation relation with the ambient magnetic field lines.The structures become accessible using interactive controls and volume-rendering techniques.The user can also specify a virtual-telescope (Observer) in the environment and output two-dimensional (2D) projections (polarization and column density maps) with given angular size.These polarization maps, intended to simulate observations of the polarized thermal dust emission in the sub-millimeter, encode information on the density and the magnetic field properties of the observed portion of the Galaxy.These outputs can then be analyzed externally as if they were actual observations.To allow for statistical analyses and to make Asterion produce massive amounts of selected outputs, input specifications in the form of setup files (CSV format) can be passed into Asterion, which runs over each and produces the corresponding output maps in an automated way.
In its current implementation, Asterion simulates both the Galaxy at the large scales and dust clouds with higher spatial resolution in a smaller volume.Both volumes are sampled by cubic grids made of 256 3 voxels.The low-resolution, large-scale, grid is 40 kpc on a side and the high-resolution grid has a size that can be adjusted by the user.Both grids are represented as Volume-Texture on the GPUs and, for each cell, local density and magnetic field as defined by model settings is rendered.
Within Asterion, both the location of the high-resolution observed box (OB) and the observer can be determined by the user who is free to move independently through space.The position of the observer is set in a heliocentric Cartesian coordinate system and the OB location is fixed through specification of the longitude, latitude and distance as seen from the observer.The outputs that are generated by the software and that are of interest in this study (see Sect. 2.1.2) are related to the polarized thermal dust emission and correspond to the view of the OB (containing simulated clouds at high-resolution) as seen by the observer.Only the volume spanned by the OB is mapped into the outputs.We make this choice to make it possible, in the future, to simulate large volumes with high-resolution capabilities by the production of mosaic in 3D.
1 UE4 is a complete suite of development tools that allows for cutting edge visualization and immersive virtual worlds, multi-platform deployment, asset and plugin marketplace, among other features (https: //www.unrealengine.com/en-US/unreal).

Polarized thermal dust emission model
For a given setup of the simulated magnetized ISM, synthetic observation of the thermal dust polarization and column density can be produced by the software.These come in the form of four maps: I, Q, U for the three Stokes parameters of the linear polarization of the thermal dust emission, and K for the total column density.These quantities are obtained through integration along the lines of sight through the OB according to the following equations: where r is the radial distance from the observer along a given line of sight with direction specified by n.The different terms in the equation are: p is a parameter dependent on dust polarization properties (grain cross sections and the degree of alignment with the magnetic field), taking a value of 0.25 in this work 2 n d (r, n) is the 3D dust grain density at position (r, n) -α(r, n) is the inclination angle between the magnetic field and the line of sight at (r, n) -ψ(r, n) is the local plane-of-sky polarization angle expressed as where B θ and B φ are the local transverse components of the magnetic field in the local spherical coordinate basis (e r , e θ , e φ ) with e θ pointing towards the South pole.ψ(r, n) is rotated 90 • from the position angle of the plane-of-sky component of the local magnetic field.
This formulation assumes that the emissivity of the dust grains is constant through the Galaxy despite known evidences for variations of the dust spectral energy distribution (e.g., Finkbeiner et al. 1999;Planck Collaboration XI 2014;Pelgrims et al. 2021a).
As we do not consider multi-frequency analysis in this study, we are not affected by such variations which, in any case, should be small since we only produce maps corresponding to small 3D volumes to achieve high spatial resolution.

Populating the Galactic space with dust
Within Asterion the Galaxy is populated by dust grains according to a smooth large-scale density distribution on top of which filaments sheets wi = 0.1 < l a t e x i t s h a 1 _ b a s e 6 4 = " I 9 0 n V o / N w 5 R t a v e O j r + p j M V q / Y s y 0 3 x r m 5 J A y 7 / H O c 0 a F S s 8 p 5 V u d g v V U + y U e e x h W 3 s 0 j w P U c U Z a q i T t 4 9 H P O H Z O D e k M T L G n 6 l G L t N s 4 t s y H j 4 A B O W R q A = = < / l a t e x i t > wi = 0.1 < l a t e x i t s h a 1 _ b a s e 6 4 = " I 9 0 n V o / N w 5 H E S 3 / u g W 9 n N i 7 x Y n W b q e K q d F f u b 9 1 h 7 q r u N 6 O 9 m X j 6 x E r f E / q W b Z v 5 X p 2 q R 6 O N Y 1 y C o p k g z q j q W u a S 6 K + r m 5 p e q J D l E x C n c o 3 h M m G n l t M + m 1 i S 6 d t V b R 8 f f d K Z i 1 Z 5 l u S n e 1 S 1 p w P b P c c 6 C R q V s H 5 Q r l 4 e l 6 m k 2 6 j x 2 s I t 9 m u c R q j h H D X X y 9 v G I J z w b F 4 Y 0 x s b k M 9 X I Z Z p t f F v G w w f 7 V p G k < / l a t e x i t > wi = 0.3 < l a t e x i t s h a 1 _ b a s e 6 4 = " 6 v 8 S q d R T 3 g 2 L g x p j I 3 J Z 6 o x l 2 l 2 8 G 0 Z D x 8 A J Z G m < / l a t e x i t > wc = 0.0005 < l a t e x i t s h a 1 _ b a s e 6 4 = " z 7 Y k K Y t T h k P s g L 5 9 2 j U 1 G 9 B e e i Z K z e g U n 9 6 Y l C b 2 S B N S X k x Y n m a q e K q c J f u b 9 0 R 5 y r u N 6 e 9 q r y G R t a v e O j r + p j M V q / Y s y 0 3 x r m 5 J A y 7 / H O c 0 a F S s 8 p 5 V u d g v V U + y U e e x h W 3 s 0 j w P U c U Z a q i T t 4 9 H P O H Z O D e k M T L G n 6 l G L t N s 4 t s y H j 4 A B O W R q A = = < / l a t e x i t > and Table 1 for parameter definitions and labels.)The circular regions displayed on those maps have an angular radius of 12.53 • .All maps share the same color scale, which is expressed in 10 21 cm −2 .
dust-cloud structures are added.The large-scale density distribution follows a parametric model that implements either the exponential disk model (ED) or the four logarithmic-spiral arms model (ARM4) discussed and adjusted to a full-sky Planck map in Pelgrims et al. (2021b).Their ARM4 model with best-fit parameter values is implemented as default in Asterion and we use that model in this work.The large-scale model fills the lowresolution grid which is then integrated to generate a dust column density map as seen from the Sun's location.A latitude profile is subsequently built and a normalization factor is computed such that it fits the well-known relation (Kulkarni & Heiles 1987;Heiles & Troland 2003): where N H I is the column density of neutral hydrogen atoms and b is Galactic latitude.The normalization factor fixes the proportionality in Eq. 4 leading the large-scale density distribution model to determine the number density of hydrogen atoms in any place in the Galaxy.The propagation of this normalization factor to Eqs. 1 to 3 implies that our polarization maps (I, Q, and U) are given in units of dust grain emissivity, as stated above.
The particle density assigned to the high-resolution box (i.e., the OB) is calculated from this large-scale model which is evaluated at the location of the center of the OB.Inside the OB, the total number of particles is determined by assuming that the matter in the OB is constant and therefore, scales with the OB volume.The total number of particles is then divided into two phases of the ISM which are the Warm Neutral Medium (WNM) and the Cold Neutral Medium (CNM) (e.g., Wolfire et al. 2003).In Asterion we assume that the WNM corresponds to a very diffuse component that follows the large-scale density distribution model, whereas the CNM forms the ISM structures (the clouds) which will appear brighter on the maps.The modeling of the clouds which we consider as being either filament-like or sheetlike structures is discussed in Sect.2.1.4.The relative contribution of CNM and WNM to build the column density is a free parameter in Asterion.We currently do not consider the possibility for a contribution from an unstable neutral medium phase (Ghosh et al. 2017;Adak et al. 2020).
Assuming a given number density of dust particles in clouds and given the volume of the dust clouds that we model, the number of clouds required to account for the mass of the CNM in the OB is automatically determined.The clouds are then placed randomly within the OB according to a uniform distribution.Both the particle number density and the morphology of the clouds are parameters that the user can tune.
The orientation of clouds in the 3D space are not random.It is now well established observationally that the main axes of the projected shapes of ISM structures show orientations that are not independent of the orientation of the ambient plane-of-sky component of the magnetic field.Both orientations appear preferentially parallel or perpendicular, depending on their column density (e.g., Planck Collaboration Int.XXXII 2016; Planck Collaboration Int.XXXV 2016; Clark 2017; Soler 2019) and this alignment and misalignment may lead to differences in characteristics of the polarization power spectra (Planck Collaboration Int.XXXVIII 2016; Huffenberger et al. 2020;Clark et al. 2021).To control this correlation, Asterion allows the user to vary the angle of each cloud in 3D with the local magnetic field.We refer to this angle as the offset angle which we denote ω.The implementation of magnetic field is discussed in Sect.2.1.5.The correlation between dust density structures and local magnetic field is further discussed in Sect.2.1.6.

Dust cloud structures
Asterion allows for the creation and visualization of dust clouds of the ISM that are either filament-like or sheet-like.The basic filaments sheets e 6 I 8 5 d 3 G 9 H e 1 1 4 h   regular geometric shapes, which are then perturbed, are cylinders and flat rectangular parallelepipeds for filaments and sheets, respectively.The user can control the shape parameters (size, axis ratio, etc) and the density of dust grains within the structure.Generally, input parameters control the ranges of possible values and, internally, the software seeks random numbers from a uniform distribution to create any given implementation.This feature ensures intrinsic scatter which account for the fact that all dust clouds in the ISM are not identical.These basic shapes are then perturbed in order to make their appearance more realistic, producing 'wiggles'.For this purpose, Asterion relies on the use of Perlin noise.Perlin noise is a type of gradient noise used to increase appearance of realism in computer graphics simulating textures (Perlin 1985).This technology is extensively used in video games and by the imaging industry.Two parameters control the properties of the noise, an amplitude (w i ) and a correlation length (w c ), that the user can change.
Finally, the density distribution within the OB is smoothed in 3D through a distance-weighted averaging of the densities from direct neighboring voxels.This smoothing is such that the total mass is conserved.It is applied to avoid sharp transitions between structures and their surroundings.
When the emission is integrated along the lines of sight (see Eq. 1), dust clouds can appear as illustrated in Fig. 1 where we show the intensity maps (Eq. 1) of a population of filaments (left panel) and sheets (right panel) obtained for different combinations of the wiggle parameters.In these sky maps, the major axes of all the structures are horizontal.In the case of sheet-like clouds, the orientation of the minor axes are not bounded, which makes it possible to observe a cloud either edge-on or face-on.Among maps of filaments (sheets), only the wiggle parameters change, the 3D positioning of the clouds and their main orientation remain unchanged.In Fig. 1, the wiggle amplitude is constant through columns and increases from left to right while the wiggle correlation length is constant through rows and increases from top to bottom.The larger w i and the smaller w c , the noisier the structures.

Magnetic field
A parametric model for the regular large-scale Galactic magnetic field is assumed in Asterion.Two models of the magnetic field are implemented: the Logarithmic Spiral Arm model of Page et al. (2007) with the best-fit parameters obtained by Pelgrims et al. (2021b) for a fit to dust polarization maps, and a four logarithmic-spiral arms model with plane-parallel field lines but with an azimuthal modulation of the strength that follow the large-scale dust density.The latter is the model that we use in this work despite the fact that it has never been tested against data.Our choice is only motivated by the fact that this model allows for an easier control of the relative orientation of the field lines with respect to the lines of sight of observation (see Sect. 2.2).

Mimicking MHD turbulence
Following the description above, our simplistic approach to model the magnetized ISM does not couple matter density distribution and magnetic field as MHD equations do.To implement such a coupling in our toy-model and allow for the magnetic field lines to depart from the large-scale regular model, we proceed as follows.
First, we add a stochastic component to the magnetic field in the OB.The stochastic component is a random realization of a Perlin-noise vector field.This choice has the advantage of creating anisotropic fluctuations in the magnetic field directions.Two parameters control the stochastic component of the magnetic field; the amplitude of the fluctuations (n i ) measured in µG and a correlation length (n c ).Then, we add the density structures (clouds).They are first placed and oriented as before and then modified (perturbed).To account for coupling we rely on the flux-freezing approximation (expected to be very accurate at the low density clouds we consider here for the high Galactic latitudes of interest), which dictates that matter can move freely along magnetic field lines.Therefore, a shift along the local (total) magnetic field line is imposed to each voxel of each cloud.The amplitude of the shifts and their spatial correlation is obtained from a scalar Perlin noise (the one that is used to generate the wiggles in the previous section) which is again controlled by the wiggle intensity (w i ) and wiggle correlation (w c ) parameters.Given our implementation, the structures can be squeezed or stretched along the magnetic field lines and compression and rarefaction can happen if the field lines converge or diverge.As a result, within individual structures, overdensities generally appear where magnetic field lines are more compressed as observed in the magnetized ISM (e.g., Heyvaerts & Priest 1983;Heiles & Troland 2003).
The effect of the stochastic component in the magnetic field is demonstrated in Fig. 2 where we show maps of the Stokes Q (Eq.2) of a population of filaments (left) and sheets (right) located in a circular region towards the North Galactic pole.The orientation of the regular magnetic field crosses the sky regions making an angle of 28 • with the horizontal which corresponds to the value of the pitch angle of the spiral arms in our model.In those panels, only the parameters of the stochastic magneticfield component vary; n i is constant in columns and increases from left to right with values of n i = 0.02, 0.05 and 0.08, and n c is constant in rows and increases from top to bottom with values of n c = 0.005, 0.015, 0.035.The parameters of the wiggles are fixed to (w i , w c ) = (0.2, 0.001).
The effects of varying the wiggle parameters to perturb the density structures in this coupling scheme are illustrated in Fig. 3. To produce those synthetic maps we proceed as for Fig. 1 (with (ω, α) = (0 • , 90 • )) and fix the realization of the stochastic component in the magnetic field (including random seed value) to be the same with parameters (n i , n c ) = (0.05, 0.0015).The amplitude of the perturbation (w i ) is constant in a column and increases from left to right with values of n i = 0.1, 0.2 and 0.3 and the correlation length (n c ) is constant in a row and increases from top to bottom with values of n c = 0.0005, 0.0010, 0.0015.Comparison of Figs. 1 and 3 reveals, as expected, that the coupling of structure perturbations to magnetic field lead to significant difference in structures' morphology as the wiggles are not isotropic anymore.As a side effect of our implementation of the coupling between magnetic field and matter density, structures may generally appear more aligned with the magnetic field than as imposed by the offset angle ω, and may appear more elongated than when turbulence is switched off.Both effects are in qualitative agreement with what MHD theory predicts (e.g., Brandenburg & Lazarian 2013).

Synthetic data: Input parameters
In this work, we want to study how the shape of ISM clouds and their alignment with the ambient magnetic field impact the measured polarization power spectra as well as to gauge the importance of the viewing angle.This goal, together with the requirement to observe sky patches covering at least one per cent of the full sky (see next section), leads us to favour certain observer-OB configurations within Asterion.
First, we want to study a portion of the ISM typical of the intermediate-to-high Galactic latitude sky.Therefore, we select a region placed at (l , b ) ≈ (180 • , 44 • ) so that the entire observed region, as seen from the Sun, is above latitude 30 • .We want to resolve density structures at the pc scale and therefore decide for an OB size of 200 pc (the voxels are 0.78 pc on a side).To achieve an observed sky area of 1% of the full sky, the circle defined by the intersection of the observation cone and the back side of the OB has to have an angular size of ≈ 25 • and, therefore, the center of the OB has to be placed at 350 pc from the observer.The distance between the Sun and the OB is d ≈ 370 pc and the observer will not be at Sun position.
Second, the fact that the observables returned by Asterion reflect only the portion of the ISM that is encompassed by the OB makes it possible to study the same ISM configuration from different viewpoints.We take advantage of this feature to control values of the inclination angle between the (regular) magnetic field orientation at the center of the OB and the line of sight toward the center of the OB (the angle α in Eqs. 1 to 3).We notice that in general α may vary across the outputs.The largescale regular magnetic field model implements magnetic field lines following a logarithmic-spiral pattern with a constant pitch angle of 28 • and with zero out-of-plane component.Therefore, we choose to place the observer at a height z = 257.1 pc (= d sin(b )) above the Galactic disk.We move the observer in this plane while keeping the same distance to the OB center in order to explore the effect of the inclination angle on the observations.In this case, all observations made by the observer are at latitude 0 • and moving the observer around the OB changes the longitude only.We sample the inclination angles from 90 • (magnetic field in the plane of the sky) to 0 • (magnetic field along the line of sight pointing to the center of the OB) with step of 15 • .The angle that the main axes of the clouds make with the orientation of the magnetic field is the other parameter that we want to explore.We vary this angle from 0 to 90 • in 15 • step.To allow us to focus only on the effects from the shape, the viewing angle and the misalignment between structures and magnetic field, we start our study focusing on simulating data without the addition of the stochastic component in the magnetic field and the correlations between density perturbations and magnetic field.Another set of synthetic maps is generated in Sect.3.4 where we study the effects from those variables.
To decide on a default setting for the large number of free parameters allowed by Asterion, we rely on observational constraints and on visual inspection of the resulting column density maps so that the latter appear as realistic as possible.
First, according to Heiles & Troland (2003), we consider that the CNM (the cloud structures) must account for 39% of the gas mass and that the remaining 61% is attributed to the WNM (the large-scale density distribution).Second, based on Spitzer (1978), we consider that the particle number density in clouds must be chosen randomly for each cloud in the range of 10 to 70 cm −3 .
Following Heiles (1976), the typical length of clouds is expected to be around 20 pc with a width in an approximate ratio of 10:1.We explored Asterion's parameter space around those values also varying the parameters describing the wiggles (the fluctuations in the density distribution within the clouds), and the values or range of values of all parameters we choose are listed in Table 1.The wiggle parameters are fixed to the same values for filaments and sheets.For filaments, the length is drawn from a uniform random distribution in the range 7 to 25 pc while the thickness varies uniformly from 1.5 to 4.5 pc.For sheets, the length is drawn from a uniform random distribution in the range from 5 to 20 pc, the width of the sheets takes random values in the range from 10 to 20, and the thickness (the height) in the range 1 to 3.
We adopt the observer-OB configurations presented above together with the chosen settings for the cloud morphology and filaments sheets wi = 0.1 < l a t e x i t s h a 1 _ b a s e 6 4 = " I 9 0 n V o / N w 5 H E S 3 / u g W 9 n N i 7 x Y n W b q e K q d F f u b 9 1 h 7 q r u N 6 O 9 m X j 6 x E r f E / q W b Z v 5 X p 2 q R 6 O N Y 1 y C o p k g z q j q W u a S 6 K + r m 5 p e q J D l E x C n c o 3 h M m G n l t M + m 1 i S 6 d t V b R 8 f f d K Z i 1 Z 5 l u S n e 1 S 1 p w P b P c c 6 C R q V s H 5 Q r l 4 e l 6 m k 2 6 j x 2 s I t 9 m u c R q j h H D X X y 9 v G I J z w b F 4 Y 0 x s b k M 9 X I Z Z p t f F v G w w f 7 V p G k < / l a t e x i t > wi = 0.2 < l a t e x i t s h a 1 _ b a s e 6 4 = " H / z q 4 + 7 o 7 y f w + I 4 F N q k j 0 H t r k H k 7 e M R T 3 g 2 L g x p j I 3 J Z 6 q R y z T b + L a M h w / 9 t p G l < / l a t e x i t > wi = 0.3 < l a t e x i t s h a 1 _ b a s e 6 4 = " 6 v 8 S q d B e e i Z K z e g U n 9 6 Y l C b 2 S B N S X k x Y n m a q e K q c J f u b 9 0 R 5 y r u N 6 e 9 q r y G  x Y n W b q e K q d F f u b 9 1 h 7 q r u N 6 O 9 m X j 6 x E r f E / q W b Z v 5 X p 2 q R 6 O N Y 1 y C o p k g z q j q W u a S 6 K + r m 5 p e q J D l E x C n c o 3 h M m G n l t M + m 1 i S 6 d t V b R 8 f f d K Z i 1 Z 5 l u S n e 1 S 1 p w P b P c c 6 C R q V s H 5 Q r l 4 e l 6 m k 2 6 j x 2 s I t 9 m u c R q j h H D X X y 9 v G I J z w b F 4 Y 0 x s b k M 9 X I Z Z p t f F v G w w f 7 V p G k < / l a t e x i t > wi = 0.2 < l a t e x i t s h a 1 _ b a s e 6 4 = " H / z q 4 + 7 o 7 y f w + I 4 F N q k j 0 H t r k H R T 3 g 2 L g x p j I 3 J Z 6 q R y z T b + L a M h w / 9 t p G l < / l a t e x i t > wi = 0.3 < l a t e x i t s h a 1 _ b a s e 6 4 = " 6 v 8 S q d K Y t T h k P s g L 5 9 2 j U 1 G 9 B e e i Z K z e g U n 9 6 Y l C b 2 S B N S X k x Y n m a q e K q c J f u b 9 0 R 5 y r u N 6 e 9 q r y G x A t f E / q W b Z v 5 X J 2 s R 6 O F U 1 e B R T Z F i Z H V M u 6 S q K / L m 5 p e q B D l E x E n c p X h M m C n l t M + m 0 i S q d t l b R 8 X f V K Z k 5  Table 1.Parameters and values kept constant in our toy-models for filaments and sheets.generate 10 realizations of the magnetized ISM for each view point by changing the random seed in Asterion.To summarize, for both filaments and sheets, we explore seven values for the inclination angle, seven values for the offset angle (between the structures' main axis and the magnetic field) and we generate 10 realizations for each couple of angles in order to infer the variance induced by the positioning of the clouds and the scatter from the specific choices of cloud size and shape and small-scale fluctuations.Therefore, we generate a set of 490 observations for each.An example of I, Q and U polarization maps is shown in Fig. 4 for clouds being filaments (top) or sheets (bottom) with major axis aligned with the ambient magnetic field lines; the latter making an angle of 90 • with the line of sight for the filament output and 60 • for the sheet output.Since the magnetic field runs horizontally, most of the polarization signal is observed in positive Stokes Q.However, given the large area of our output maps, some signal is also observed in Stokes U far from the center of the maps because of projection effects.Even for the case of a uniform magnetic field in the OB, the effective inclination and position angles of the magnetic field with respect to the line of sight and in the plane of the sky vary in the image.The smaller the inclination angle, the larger the variation in the position angle.

Formalism and main characteristics
The analyses of polarization power spectra are generally carried out using pseudo angular power spectra defined as: D XY = ( + 1) C XY /(2π) where C XY are the angular (auto-or cross-) power spectra and where X and Y refer either to T , E or B (e.g., Bracco et al. 2019b).
The E/B power asymmetry is measured through the R EB ratio, which is obtained by averaging the ratio of the auto-power spectra D EE and D BB over a specified multipole range where • stands for the mean over multipole bins.
To quantify the correlation between the T and E power spectra we use the normalized parameter r T E introduced by Caldwell et al. (2017).r T E takes values 1, −1 and 0 in case of perfect correlation, perfect anti-correlation and absence of correlation, respectively, and is defined as The correlation between the power spectra T and B discussed in Sect. 4 is also computed through the correlation coefficient r T B obtained by substituting E for B in Eq. 8. Additionally, visual inspection of the power spectra informs us that most of the T T , EE and BB auto-power spectra present an apparent power-law dependence in in the range ∈ [100, 500].Consequently, we decide to adjust a power-law power spectrum model of the form D XX = A XX 80 ( /80) a XX +2 , where X ∈ {E, B}, and to characterize the spectra through the values of A XX 80 and a XX .

Computation from synthetic maps
We use the Xpol code3 (Tristram et al. 2005) to compute the polarization power spectra of the dust polarized sky and account for incomplete sky coverage.The reliability of Xpol at estimating polarization power spectra has been tested for sky area larger or equal to one per cent of the full-sky4 .This limit fixes the minimum angular size of the output that we generate using Asterion.This choice, together with the spatial resolution that we want to achieve in the OB, further sets the size of the OB and its distance to the observer (see previous section).
In practice, we project Asterion's outputs onto a full-sky HEALPix map with resolution parameter N side = 2048.Most of the map is empty.We then define a circular mask of radius 11.5 • centered on the outputs' center.This choice allows us to retain only the sky area that is free of artifacts from ray-tracing the volume-limited cubic OB.The sightlines that are kept pass through both the front and back faces of the OB, as seen by the observer.To avoid power leakage we further apodize (smooth) the mask using a 0.5 degree FWHM beam.This mask is applied to Asterion's outputs when projected on HEALPix map.The effective sky fraction of the mask ( f sky ) is 1% of the full-sky.Using Monte Carlo simulations we checked that our mask does not introduce bias or leakage in the power spectrum estimation following the same validation procedure as in (Pelgrims et al. 2022).
For the computation of the power spectra there are some constraints imposed for the multipole moments.The maximum value of is defined by the effective angular resolution achieved in Asterion's outputs.Visual inspection of several power spectra informs us that above = 500 there is a sudden drop of power, reflecting the effects of the limited spatial resolution achieved in the OB and the smoothing of density structure.
The threshold at low multipole values comes from uncertainties in power spectrum estimation.
Since our synthetic maps do not contain observational noise, we consider only the sampling variance as a source of uncertainty in our power spectrum estimates.The latter is linked to the number of unmasked pixels and their spatial arrangement on the sky.The analytical estimate from the sampling variance is given by (Planck Collaboration Int.XXX 2016): where f sky and ∆ bin are the sky fraction and the width of the multipole bins used to estimate the power spectra, respectively.We choose to reject power spectrum estimates for multipole below 100 so as to ensure that C /σ C ≥ 3 always.
To summarize, we use Xpol to estimate the polarization power spectra from Asterion outputs in the multipole range ∈ [100, 500] adopting bins of width 20.We compute the T T , EE, BB and T E spectra for every output of our sample and characterize them through their R EB and r T E values.An example of such a set of polarization power spectra is shown in Fig. 4 along with the corresponding polarization maps.

Analysis
In this section, we explore the dependence of the values of R EB and r T E with respect to (i) the offset angle between the clouds' long axes and the magnetic field (ω) and (ii) the inclination angle of the magnetic field with the line of sight, for both filament-like and sheet-like clouds.We vary both angles from 0 • to 90 • with a step of 15 • .Since the exact values of R EB and r T E depend sensitively on the specific choice of parameters used to produce the synthetic maps, we are primarily interested in the trends that they exhibit as a function with ω and α.

R EB and r T E dependence on ω
Here we infer the dependence of R EB and r T E on the angle that the major axes of the structure makes in three dimensions with the local magnetic field orientation.
We first select the synthetic maps created with the magnetic field in the plane of the sky (α = 90 • ) and we sort them according to the offset angle that the major axes of the structures make with the local magnetic field.In the top row of Fig. 5, we show the dependence of R EB (left) and r T E (right) as a function of ω, for both shape families of clouds.
For filament-like clouds, we observe that R EB as a function of ω follows a parabolic trend starting from a maximum when clouds' major axes are perfectly aligned with the magnetic field (ω = 0 • ), then decreases to reach a minimum toward ω = 45 • and increases back for structures perpendicular to the magnetic field (ω = 90 • ).This trend reflects well the generally accepted picture according to which E modes are maximized for structures perfectly parallel or perpendicular to the magnetic field whereas B modes dominate when (projected) linear structures make an angle of 45 • with it (e.g.Zaldarriaga 2001;Huffenberger et al. 2020;Clark et al. 2021).Obtaining this result, for filaments, validates our modeling and analysis pipeline.Still for filament-like clouds, we observe that r T E as a function of ω follows a decreasing squared-cosine trend, even showing negative values for large offset angles (ω ≥ 60 • ).This trend is reminiscent of the result obtained for individual straight cylinders (Huffenberger et al. 2020) and is understood by the same argument as above.When E modes are reduced due to misalignment, the corresponding E map systematically looses its correlation with the T map.The two even get anti-correlated when the structures are perpendicular to the field since E modes become negative in such a configuration.
For sheet-like clouds, the same parabolic trend is recovered for R EB as a function of ω but with a smaller absolute difference between the extrema than in the case of filament-like structures.This is also true for the squared-cosine trend of r T E versus ω.We understand these observations by the fact that sheets have one additional degree of freedom in their orientation.The fact that they could be seen either face-on, edge-on or with any intermediate angles, even though the orientation of their principal axes is fixed as compared to the magnetic field orientation (here parallel to it), has the effect of reducing the imbalance between E and B modes.Even in the case the sheets' major axes are parallel to the magnetic field, when projected on the sky, wiggles of the sheet-like structures may show appreciable angles with the magnetic field lines, contributing to the production of B modes rather than E modes.This reduces both the R EB and r T E maximum values.
For both sheet-and filament-like structures, we observe that the absolute difference between extrema in the trend of R EB versus ω is reduced when the inclination angle decreases, for example when the field lines depart from the plane of the sky.This is again understood by projection effects and is further explored in Sect.3.2.
In the bottom row of Fig. 5, we show the scatter plots of r T E versus R EB obtained for filament-like clouds (left) and sheet-like clouds (right), for all offset and inclination angles.The scatter plots are color-coded according to the offset angles.These plots offer a different and generalized view of the trends observed in the top panels which were restricted to the case α = 90 • .
Particularly, for the case of filament-like clouds (left), we notice a tail in high R EB and r T E values which corresponds to small offset angles (ω 30 • ).Then, as ω increases up to ≈ 45 • , both the values of R EB and r T E decrease forming a roundish cluster centered on (R EB , r T E ) ≈ (1, 0.1).As structures' major axes become perpendicular to the field, another smaller tail appears with negative r T E and large R EB values.This general picture is also observed for sheet-like clouds (Fig. 5, bottom right panel) although the scatter of the data points is larger; as expected due to the increase of the degree of freedom.The comparison between the 2D distributions of (R EB , r T E ) pairs obtained from filamentand sheet-like clouds is further developed below.

R EB and r T E dependence on α
Here we infer the dependence of the measured R EB and r T E as a function of the angle between the magnetic field lines and the line of sight.
We first select the synthetic maps created with all the clouds having their major axis parallel to the magnetic field (ω = 0 • ) and we sort them according to the inclination angle of the field with respect to the line of sight.As in the previous subsection, we have for each α value 10 realizations with filament-like and sheet-like clouds of the same ISM volume that is observed from seven different observer positions (see Sect. 2.2).A different set of polarization power spectra corresponds to each realization.In the top row of Fig. 6, we show the dependence of the corresponding values of R EB (left) and r T E (right) as a function of α.
For filament-like clouds, the viewing-angle and projection effects on R EB and r T E are clearly shown.R EB decreases following a squared-sine trend from its maximum (at ≈ 1.8) obtained for lines of sight perpendicular to the field lines (α = 90 • ) to its minimum (at ≈ 1) obtained for lines of sight parallel to the field (α = 0 • ).As for r T E , it also decreases from its maximum when α ≈ 90 • to zero when α ≈ 0 • but following a sine trend.For lines of sight and field lines nearly parallel, the degree of polarization is low and any small variations generate polarization with random orientations.This generates equal amounts of E and B modes (thus R EB ≈ 1) with no particular correlation with intensity (thus r T E ≈ 0).Similar trends are observed for sheet-like clouds with, however, much reduced difference between extrema.As before, we interpret this reduction by the presence of wiggles and the extra freedom in orientation that sheet-like structures allow for.The maxima of both R EB and r T E also do not appear for α = 90 • but at lower values ( 60 • ).We speculate that this peculiarity may be related to projection effects of the 3D structures.
In the bottom row of Fig. 6, we show the same scatter plots as shown in the bottom row of Fig. 5 but color-coded by their corresponding inclination angle values.We observe that the data points corresponding to α = 90 • draw a croissant-like pattern (already outlined in the bottom panels of Fig. 5) in which R EB is maximized for ω = 0 • and 90 • and minimized for ω = 45 • , and that is tighter for filament-like clouds than for sheet-like clouds.Then, as α decreases, the distributions of (R EB , r T E ) pairs depart from the croissant shape and become rounder with a center at about (1, 0.1).The extreme cases are for α = 0 • for which the data points (all ω values included) show a very large scatter for both filament-like and sheet-like clouds.

Ambiguity from the shapes of dust clouds
It is obvious that sheet-like structures may appear filament-like if they are seen nearly edge-on, that is, if their minor axes is close to the plane of the sky.Such an example is shown in Fig. 4.However, if dust clouds of the ISM are sheet-like rather than filament- e F e 3 p A H b P 8 c 5 C e q l o r 1 T L J 3 u F s o H 6 a i z 2 M A m t m m e e y j j G B X U y P s K j 3 j C s 3 F i 3 B h 3 x v 1 n q p F J N e v 4 t o y H D 7 j 5 k Z w = < / l a t e x i t > R EB < l a t e x i t s h a 1 _ b a s e 6 4 = " p c I j W K m M T w s P w a W a I s 2 o 6 n j m k u q u q J t b X 6 q S 5 B A R p 3 C P 4 j F h r p X T P l t a k + j a V W + Z j r / p T M W q P c 9 y U 7 y r W 9 K A y z / H O Q u a l V L 5 o F Q 5 P y x W a 9 m o 8 9 j B L v Z p n k e o 4 g x 1 N M h 7 i E c 8 4 d l o G R P j z r j / T D V y m W Y b 3 5 b x 8 A F p f J Y u < / l a t e x i t > R EB < l a t e x i t s h a 1 _ b a s e 6 4 = " p c I j W K m M T w s P w a W a I s 2 o 6 n j m k u q u q J t b X 6 q S 5 B A R p 3 C P 4 j F h r p X T P l t a k + j a V W + Z j r / p T M W q P c 9 y U 7 y r W 9 K A y z / H O Q u a l V L 5 o F Q 5 P y x W a 9 m o 8 9 j B L v Z p n k e o 4 g x 1 N M h 7 i E c 8 4 d l o G R P j z r j / T D V y m W Y b 3 5 b x 8 A F p f J Y u < / l a t e x i t > r T E < l a t e x i t s h a 1 _ b a s e 6 4 = " 1 U g 3 A R h 6 A w 6 J 4 X 9 3 H 2 b S 0 I x X 5 v s = " > A A A C y X i c j V H L S s N A F D 2 N r 1 p f V Z d u g k V w V Z I q 6 L I o g u C m Q l 9 Q q y T T a Y 2 m S Z x M x F q 6 8 g f c 6 o + J f 6 B / 4 Z 0 x B b W I T k h y 5 t x 7 z s y 9 1 4 1 8 L 5 a W 9 Z o x p q Z n Z u e y 8 7 m F x a X l l f z q W j 0 O E 8 F 4 j Y V + K J q u E 3 P f C 3 h N e t L n z U h w p + / 6 v O F e H 6 p 4 4 5 a L 2 A u D q h x E v N 1 3 e o H X 9 Z g j i a q L 8 2 H 1 a H S R L 1 h F S y 9 z E t g p K C B d l T D / g j N 0 e F e 3 p A H b P 8 c 5 C e q l o r 1 T L J 3 u F s o H 6 a i z 2 M A m t m m e e y j j G B X U y P s K j 3 j C s 3 F i 3 B h 3 x v 1 n q p F J N e v 4 t o y H D 7 j 5 k Z w = < / l a t e x i t > r T E < l a t e x i t s h a 1 _ b a s e 6 4 = " 1 U g e F e 3 p A H b P 8 c 5 C e q l o r 1 T L J 3 u F s o H 6 a i z 2 M A m t m m e e y j j G B X U y P s K j 3 j C s 3 F i 3 B h 3 x v 1 n q p F J N e v 4 t o y H D 7 j 5 k Z w = < / l a t e x i t > F ilaments < l a t e x i t s h a 1 _ b a s e 6 4 = " Z z J d 0 w a W a I s 2 o 6 n j m k u q u q J t b X 6 q S 5 B A R p 3 C P 4 j F h r p X T P l t a k + j a V W + Z j r / p T M W q P c 9 y U 7 y r W 9 u a S 6 K + r m 5 p e q J D l E x C n c o 3 h M 2 N P K S Z 9 N r U l 0 7 a q 3 j o 6 / 6 U z F q r 2 X 5 a Z 4 V 7 e k A d s / x z k N G p W y v V e u n O 2 X q k f Z q P P Y w j Z 2 a Z 4 H q O I E N d T J + w q P e M K z c W p c G 7 f G 3 W e q k c s 0 m / i 2 j I c P w m m R o A = = < / l a t e x i t > ↵ = 90 < l a t e x i t s h a 1 _ b a s e 6 4 = " p f k T W L q T T a l E 6 X 1 J + 8 x J Y I 0 w B e M = " > e y z a z W J r d 3 0 l t n 4 q 8 0 0 r N n z L D f F m 7 k l D b j y f Z w / Q W O t X F k v r x 1 s l K q 7 2 a j z W M I y V m m e m 6 h i H z X U y f s S D 3 j E k 9 N y b p x b 5 + 4 j 1 c l l m k V 8 W c 7 9 O 8 / L l X 0 = < / l a t e x i t > ↵ = 90 < l a t e x i t s h a 1 _ b a s e 6 4 = " p f k T W L q T T a l E 6 X 1 J + 8 x J Y I 0 w B e M = " > α dependence of R EB and r T E r T E < l a t e x i t s h a 1 _ b a s e 6 4 = " 1 U g 3 A R h 6 A w 6 J 4 X 9 3 H 2 b S 0 I x X 5 v s = " > A A A C y X i c j V H L S s N A F D 2 N r 1 p f V Z d u g k V w V Z I q 6 L I o g u C m Q l 9 Q q y T T a Y 2 m S Z x M x F q 6 8 g f c 6 o + J f 6 B / 4 Z 0 x B b W I T k h y 5 t x 7 z s y 9 1 4 1 8 L 5 a W 9 Z o x p q Z n Z u e y 8 7 m F x a X l l f z q W j 0 O E 8 F 4 j Y V + K J q u E 3 P f C 3 h N e t L n z U h w p + / 6 v O F e H 6 p 4 4 5 a L 2 A u D q h x E v N 1 3 e o H X 9 Z g j i a q L 8 2 H 1 a H S R L 1 h F S y 9 z E t g p K C B d l T D / g j N 0   n j m k u q u q J t b X 6 q S 5 B A R p 3 C P 4 j F h r p X T P l t a k + j a V W + Z j r / p T M W q P c 9 y U 7 y r W 9 K A y z / H O Q u a l V L 5 o F Q 5 P y x W a 9 m o 8 9 j B L v Z p n k e o 4 g x 1 N M h 7 i E c 8 4 d l o G R P j z r j / T D V y m W Y b 3 5 b x 8 A F p f J Y u < / l a t e x i t > r T E < l a t e x i t s h a 1 _ b a s e 6 4 = " 1 U g 3 A R h 6 A w 6 J 4 X 9 3 H 2 b S 0 I x X 5 v s = " > A A A C y X i c j V H L S s N A F D 2 N r 1 p f V Z d u g k V w V Z I q 6 L I o g u C m Q l 9 Q q y T T a Y 2 m S Z x M x F q 6 8 g f c 6 o + J f 6 B / 4 Z 0 x B b W I T k h y 5 t x 7 z s y 9 1 4 1 8 L 5 a W 9 Z o x p q Z n Z u e y 8 7 m F x a X l l f z q W j 0 O E 8 F 4 j Y V + K J q u E 3 P f C 3 h N e t L n z U h w p + / 6 v O F e H 6 p 4 4 5 a L 2 A u D q h x E v N 1 3 e o H X 9 Z g j i a q L 8 2 H 1 a H S R L 1 h F S y 9 z E t r T E < l a t e x i t s h a 1 _ b a s e 6 4 = " 1 U g 3 A R h 6 A w 6 J 4 X 9 3 H 2 b S 0 I x X 5 v s = " > A A A C y X i c j V H L S s N A F D 2 N r 1 p f V Z d u g k V w V Z I q 6 L I o g u C m Q l 9 Q q y T T a Y 2 m S Z x M x F q 6 8 g f c 6 o + J f 6 B / 4 Z 0 x B b W I T k h y 5 t x 7 z s y 9 1 4 1 8 L 5 a W 9 Z o x p q Z n Z u e y 8 7 m F x a X l l f z q W j 0 O E 8 F 4 j Y V + K J q u E 3 P f C 3 h N e t L n z U h w p + / 6 v O F e H 6 p 4 4 5 a L 2 A u D q h x E v N 1 3 e o H X 9 Z g j i a q L 8 2 H 1 a H S R L 1 h F S y 9 z E t g p K C B d l T D / g j N 0 e 1 S 1 p w P b P c U 6 D R q V s 7 5 U r Z / u l 6 l E 2 6 j y 2 s I 1 d m u c B q j h B D X X y v s I j n v B s n B r X x q 1 x 9 5 l q 5 D L N J r 4 t 4 + E D u y W R n Q = = < / l a t e x i t > != 0 < l a t e x i t s h a 1 _ b a s e 6 4 = " m 7 b k Z s S 6 u t l B o B 3 s S q 0 6 k w r n p i r 6 5 8 6 U q R Q 4 p c R p 3 K S 4 J M 6 M c 9 d k x m s z U r n v r m / i b y d S s 3 r M i N 8 e 7 v i U N 2 P s 5 z n F Q 2 6 l 4 u 5 W d 8 7 3 y 0 X E x 6 j l s Y B P b N M 9 9 H O E U Z 6 i S d x + P e M K z V b d u r T v r / j P V m i g 0 6 / i 2 r I c P J W K V P Q = = < / l a t e x i t > Fig. 6.Dependence of R EB and r T values on α, the angle that makes the magnetic field lines with the line of sight.Top row shows R EB (left) and r T E (right) versus α for ω = 0 • for both filaments (blue) and sheets (orange).Plotting conventions are the same as in Fig. 5 except that the color-coded in the bottom row is according to α.
like structures, and given that generally the 3D orientation of structures does not depend on the observer position (us), there is no reason why all the minor axes of ISM clouds would lie in the plane of the sky 5 .As a consequence, due to projection effects, ripples in the sheets may produce filament-like patterns (ridges) uncorrelated with the projected orientation of the principal axes of the sheets, or of the magnetic field.Therefore, this could induce different patterns in the polarization maps and therefore, different characteristics of the (polarization) power spectra.
We studied the dependence of the summarizing characteristics of the polarization power spectra on the offset angle and on the inclination (viewing) angle for outputs generated out of both, filament-like and sheet-like clouds.Beside the overall similar trend observed for both shape types, we showed that sheetlike structures lead to a generally looser correlation of the R EB and r T E values than filaments.This reflects the increase of degree of freedom in the projected orientation of substructures in sheets than in filaments.In addition, reaching simultaneously large R EB and large positive r T E values for sheets appears less common.This increase of the possible ways to project a sheet on the sky results in the reduction of the E to B power asymmetry and the correlation between T and E.
However, despite these small differences, we find that the distributions of (R EB , r T E ) pairs from maps of filament-like and sheet-like clouds largely overlap.This is best seen in Fig. 7 where we show contour plots of the distributions of points on the (R EB , r T E ) plane for both shape categories on a same figure.The three embedded contours enclose about 10, 70 and 90 per cent of the data points.The overlap is striking and the larger scatter of data points from outputs made out of sheet-like clouds is evident.
Figure 7 makes it clear that the shape of the cloud structures cannot be distinguished by the unique consideration of the R EB and r T E values.The observation of an (R EB , r T E ) pair appears almost as likely coming from either family of shapes.This is even more true if one relaxes the constraints of the morphology (axis ratios) and size of the clouds, which we have kept fixed in this study (see Sect. 2.2).
Aside from the R EB and r T E characteristics, power spectra of outputs made from filaments and sheet-like structures are similar.This can already be inferred from the right column of Fig. 4 where, for those examples, it is seen that amplitudes and slopes of spectra have similar values.This is better demonstrated in Fig. 8 where we present histograms of the spectral indices of power-law fits to the auto-power spectrum of E and B modes, for filament (left) and sheet-like (right) clouds.No substantial difference can be spotted between cloud shapes.
Based on these considerations, we argue that interpreting polarization power spectra solely in terms of filament characteristics and their relationship to the ambient magnetic field, is risky because of this degeneracy on the intrinsic shape of the clouds.

Turbulent ISM and magnetic field
The toy-model simulations of the magnetized ISM corresponding to the polarization maps used in the analysis so far are somewhat unrealistic in the sense that they lack known physics.In our model, the orientations of density caustics and magnetic field 5 We note that even in the scenario where the Sun is embedded in the Local Bubble and the magnetic field lines follow the surface of the shell of the Local Bubble (as discussed in Sect.4), it is difficult to find a physical argument that would make the minor axes of all the sheets to lie in the plane of the sky.are not correlated at small scales although we expect them to be, at least to some extent, due to flux-freezing (Heiles & Crutcher 2005).These simple, but non-trivial simulations have allowed us to study in detail the dependence of the summary characteristics of the polarization power spectra on the inclination angle, the offset angle and the type of dust cloud morphology.Here, we want to study the effects of coupling density structures to a nonregular magnetic field on our results.To do so, we exploit the possibility in Asterion to add a stochastic component to the regular magnetic field and to generate density perturbations along the local (perturbed) magnetic field lines (see Sect. 2.1.6).Example of outputs obtained with such coupling were presented in Figs. 2 and 3 for Stokes Q maps and K maps, respectively.For the purpose of our analysis we fix the parameters for the stochastic component in magnetic field and for the coupled wiggles such that (i) the outputs are visually similar to the dust polarization sky at high-Galactic latitudes and (ii) the characteristics of polarization power spectra of maps made out of filaments are close to observed values when ω = 0 • and α = 90 • .This choice, leads to the set of parameters (n i , n c ) = (0.05, 0.0015) and (w i , w c ) = (0.2, 0.001) which we do not change while exploring the effect from varying the ω and α angles as in Sect.3.With those values we notice that generally the simulated density structures coupled to the magnetic field show their effective major axes more aligned with the local mean magnetic field (ω eff 45 • ) than set by the input ω values.This is a desirable feature motivated by observational evidences that clouds of the diffuse ISM, both H I fibers and dust clouds, tend to be preferentially aligned with the ambient magnetic field (McClure-Griffiths et al. 2006;Clark et al. 2014;Clark et al. 2015;Martin et al. 2015;Planck Collaboration Int. XXXII 2016;Planck Collaboration Int. XXXV 2016;Planck Collaboration Int. XXXVIII 2016;Soler 2019).As a result, and as illustrated in Fig. A.1, we find that R EB remains globally larger than unity and r T E remains positive while varying ω.This is the case for both cloud shapes.Comparing Figs.A.1 and A.2 with Figs. 5 and 6, we find that the different scalings present for the range ω 30 • are severely attenuated, if not suppressed, when density structures are coupled to the magnetic field and that R EB and r T E values of outputs from sheet-like structures are in general closer to outputs from filament-like structures.The main effects of coupling the density perturbations and the magnetic field is that 2D distributions of R EB and r T E from the two types of shapes do overlap even more than in the absence of such coupling and that negative r T E values are much less frequent.This is clearly seen in Fig. 9.

Discussion
In this paper we do not attempt to model the actual dust polarization sky observed at submillimeter wavelengths, for example by Planck at 353 GHz.Rather, our focus is to explore whether the assumed shape of the ISM clouds may lead to erroneous conclusions in the interpretation of dust polarization power spectra.
In order to model the dust polarization sky and in particular the measured polarization power spectra, it would be necessary to sample a large volume of the Galactic space with our OB, generating a 3D mosaic, possibly in a way similar to what has been proposed by Hervías-Caimapo et al. (2021).However, our analysis shows that the viewing angle -with which the portion of the ISM is observed-is significant in determining the ratio and correlation between power spectrum modes.We find an extra dependence on α than that proposed by Huffenberger et al. (2020) for r T E only.Our results make it clear that the overall geometrical arrangement of the magnetized ISM surrounding the observer (us) leaves its marks on the power spectra and in particular on the R EB and r T E parameters.We see it unlikely that the geometry dependence averages out when considering large swaths of the sky.This is reinforced by the fact that all viewing angles are not equally likely, especially at high Galactic latitudes.The reason is that the Sun sits in the Local Bubble, a cavity of hot plasma that was presumably created by supernova explosions, which is surrounded by a shell of cold neutral gas and dust with a radius between 100 to 300 pc, depending on observation direction.The polarized signal observed in the polar caps (|b| > 60 • ) is dominated by the emission from the shell of the Local Bubble (Skalidis & Pelgrims 2019).The formation of this structure has shaped the local, magnetized ISM surrounding the Sun (Alves et al. 2018;Pelgrims et al. 2020) and has left the magnetic field to lie mostly in the plane of the sky at high |b| but may show strong departures in some places (Pelgrims et al. 2020).In light of the highlighted α-dependence that we find, we caution that the exact morphology of the magnetic field in regions dominating the polarization signal needs to be accounted for to model large areas of the sky.This would involve the production of 3D mosaics, and the computation of polarization power spectra to be compared to actual observations.This is a highly complex task, which is beyond the scope of this paper.
Recently, Pelgrims et al. (2022) showed from maps synthesized from Galaxy-sized MHD simulations that the specific observer location has a strong effect on the characterization of the polarization power spectra.In light of our results, it would be interesting to figure out how much of the reported variance can be attributed to the peculiar 3D geometry of the part of the magnetic field that is imprinted in maps and thus, in power spectra.
On the other hand, Bracco et al. (2019b) have studied polarization power spectra from Planck data at intermediate-to-low Galactic latitudes in sky patches of the same size as ours.They found that the R EB and r T E parameters span large ranges of values.With our toy-model simulation-based study we show that a significant source of scatter in those quantities may actually come from the geometrical arrangement of the magnetized ISM, in addition to the scatter that may arise from the fact that the offset angle in different clouds is probably not constant but rather follows some distribution (see e.g., Huffenberger et al. 2020).Our toy-models, however, cannot be simply compared to the results of their observational study since our synthetic maps do not account for the full cone of observation, which becomes important at low latitudes, and also because they pre-processed the polarization sky maps in order to get rid of the large-scale Galactic gradient.We do not know how this processing affects their results.
Furthermore, in an attempt to model the actual sky, it would probably be relevant to consider that dust clouds can be described by a mixed population of sheet-like and filament-like structures with varying axis ratios.In this study, we did not explore such a possibility because we want to quantify the possible cloud-shape degeneracy of the R EB and r T E values.
Currently, our results are limited by the achieved spatial resolution that is fixed by the minimum angular size of the sky patch that Xpol can handle, and the size of Asterion grid (256 3 ).We achieve a spatial resolution of ∼ 0.78 pc which is more than twice the thickness obtained for actual sheet clouds of the neutral ISM ( 0.3 pc according to Kalberla et al. 2016).
Another source of limitation of our study comes from the limited portion along distance that is mapped.This reduces the number of structures that appear on our maps as compared to real observations.This does not affect our analysis since our main results relate to the (relative) comparison of power spectrum characteristics between (i) different cloud shapes, (ii) different offset angles, and (iii) different inclination angles.Besides the special attention paid to the E/B power asymmetry and the correlation between T and E modes which were detected with high significance in dust polarization maps from Planck, a more marginal signal has been reported between the T and B modes (Planck Collaboration Int. XXX 2016;Planck Collaboration XI 2020).Other recent studies (Weiland et al. 2020;Clark et al. 2021) relied on a more accurate version of the Planck data and on external data sets to prove that the weak positive T B signal (a mark of parity violation in the dust polarization sky) is a real property of the Galactic emission and infer its origin.
A detailed study of the effect of cloud morphology on the T B signal (e.g., Bracco et al. 2019a) and a careful analysis of the dependence on ω and α along with the scale dependence of the signal is not the focus of this paper and we postpone this to future work.Meanwhile, we show in Fig. 10, that the correlation coefficient r T B between T and B modes measured from filamentlike or sheet-like structures does not strongly differ, at least when the scatter from varying ω and or α is taken into account.All distributions are centered on zero and, as with R EB and r T E , we find that sheet-like structures lead to broader distributions than filament-like structures do.

Conclusions
The search for the primordial B modes in the polarization of CMB radiation calls for the understanding and characterization of the Galactic foregrounds in order to disentangle their contribution to the maps of the polarized sky.
In this project, we investigated whether the morphology of the dust clouds has an impact on the statistical characterization of the polarized emission.We simulated interstellar clouds with filament-like and sheet-like shapes and produced corresponding synthetic polarization maps, using the software Asterion.Then, we computed their polarization power spectra in multipole range ∈ [100, 500] and we focused our study on the power asymmetry between E and B modes, and on the cross correlation between T and E modes, through the parameters R EB and r T E , respectively.We explored the dependence of R EB and r T E on the inclination angle (α, between magnetic field and the line of sight), and the offset angle (ω, between the longest cloud axis and the magnetic field) for both types of shape.Our study provides important insights on the interpretation of the polarization power spectra of dust polarization in terms of ISM properties and could help in the characterization of Galactic foregrounds to the CMB polarization.
Our results for filaments are consistent with theoretical predictions and observations.For filamentary structures aligned with the magnetic field (ω ≈ 0 • ), something that seems to be true in nature for the diffuse ISM, the power ratio R EB can be as large as ∼ 2, in agreement with Planck measurements (e.g., Planck Collaboration Int.XXX 2016).Moreover, for regular magnetic field and large inclination angles, where polarization is maximized, E-mode power of structures is maximized for offset angles ω = 0 • and 90 • and B-mode power is maximized for ω = 45 • .This behaviour is in agreement with theoretical expectations Zaldarriaga (2001) and previous studies (Rotti & Huffenberger 2019;Huffenberger et al. 2020).As a result, R EB as a function of ω draws a parabola with minimum at ω ≈ 45 • , where B modes slightly dominate E modes.In parallel, the correlation coefficient r T E goes from its maximum at ω = 0 • to its minimum at ω = 90 • .These trends (with little variance from specific ISM realizations) are similarly observed for other inclination angle values, though with lower difference between extrema as α decreases.
To the best of our knowledge, this is the first time that an analysis of polarization power spectra obtained from sheet-like interstellar clouds has been performed.We found that the parameters R EB and r T E measured from maps made from sheet-like structures show similar trends as a function of ω and α angles to those observed for filaments.However, the amplitude between extrema is generally smaller and the variance with respect to the specific ISM realizations is larger.We understood these characteristics in terms of the extra degree of freedom that sheet-like structures possess related to the orientation of their corrugations with the long axis.
Varying only the inclination and offset angles, we found that sheet-like structures and filament-like structures lead to polarization power spectra that cannot be distinguished in the plane of (R EB , r T E ) parameters.
Consequently, we argued that measured R EB and r T E values cannot be used alone to discriminate among characteristics of ISM structures and that this degeneracy should be accounted for in the interpretation and modeling of observational data.We checked that our main conclusions remain valid when a stochastic component is added to the large-scale magnetic field and, most importantly, when density structure perturbations are added such that the orientations of density caustics and magnetic field are correlated at small scales as we expect them to be due to MHD physics.
Our analysis also highlights and quantifies the dependence of the E/B asymmetry and T E correlation on the inclination angle, a dependence that is changed by the inclusion of fluctuations in the magnetic field lines and their coupling to density perturbations but that preserves its significance.This strongly suggests that the geometrical arrangement of the large scale magnetic field as viewed by the observer is a critical factor shaping polarization power spectra.
This was already suggested by Bracco et al. (2019a) and could explain the large scatter on R EB and r T E values observed from different sky patches (Bracco et al. 2019b) and the strong cosmic variance reported from MHD simulation-based studies (Kim et al. 2019;Pelgrims et al. 2022).
Finally, and for completeness, we also examined the correlation coefficient between the T and B modes (r T B ) and observe that this feature of the power spectrum cannot be used to discriminate cloud shape types either.Instead, we find that variation in ω and α generates significant dispersion in the values of r T B .Further work will be needed to investigate the possible scale dependence of the T B signal as reported by real sky observation of dust polarization.
7 1 1 d P x N Z y p W 7 V m a m + B d 3 Z I G b P 8 c 5 z y o l 4 r 2 U b F U P S 6 U z 9 N R Z 7 G H f R z S P E 9 Q x i U q q J F 3 h E c 8 4 d m o G m P j z r j / T D U y q W Y X 3 5 b x 8 A G 7 4 5 J M < / l a t e x i t > wc = 0.0020 < l a t e x i t s h a 1 _ b a s e 6 4 = " i j 7 k p J t l 9 l L p P 6 O d D a G s a B m T G q w r n r r 6 P i b z l S s 2 r M 0 N 8 G 7 u i U N 2 P 4 5 z n l Q L x X t o 2 K p e l w o n 6 e j z m I P + z i k e Z 6 g jE t U U C P v C I 9 4 w r N R N c b G n X H / m W p k U s 0 u v i 3 j 4 Q O + R J J N < / l a t e x i t > wi = 0.5 < l a t e x i t s h a 1 _ b a s e 6 4 = " 4 c J w K / 1 d J M j h 1 u Y p k 8 3 F S T K L E d w = " > A A A C z H i c j V H L S s N A F D 2 N r 1 p f V Z d u g k V w F d K q 6 E Y o u n E l Fe x D a i n J d F o H 8 y K Z K K V 0 6 w + 4 1 e 8 S / 0 D / w j t j C m o R n Z D k z L n 3 n J l 7 r x t 5 I p G 2 / Z o z Z m b n 5 h f y i 4 W l 5 Z X V t e L 6 R i M J 0 5 j x O g u 9 M G 6 5 T s I 9 E f C 6 F N L j r S j m j u 9 6 v O n e n q p 4 8 4 7 H i Q i D S z m M e M d 3 B o H o C + Z I o q 7 u u 8 I 8 N m 3 r o F s s 2 Z a t l z k N y h k o I V u 1 s d 3 B o H o C + Z I o q 7 v u 8 I 8 M a 2 y 3 S 2 W r L K l l z k L 7 A y U k K 1 a W H z B D X o I w Z D C B 0 c A S d i D g 4 S e N m x Y i I j r Y E x c T E j o O M c E B d K m l M U p w y F 2 S N 8 B 7 d o Z G 9 B e e S Z a z e g U j 9 6 Y l C b 2 S B N S X k X 1 9 Y 3 N g t b 2 4 0 k T G P G 6 y z 0 w r j l O g n 3 R M D r U k i P t 6 K Y O 7 7 r 8 a Y 7 P F P x 5 h2 P E x E G V 3 I U 8 Y 7 v D A L R F 8 y R R F 3 f d 4 V 5 Y t p W p V s o 2 p a t l z k L S h k o I l u 1 s P C C G / Q Q g i G F D 4 4 A k r A H B w k 9 b Z R g I y K u g z F x M S G h 4 x w T 5 E m b U h a n D I f Y I X 0 H t G tn b E B 7 5 Z l o N a N T P H p j U p r Y J 0 1 I e T F h d Z q p 4 6 l 2 V u x v 3 m P t q e 4 2 o r + b e f n E S t w S + 5 d u m v l f n a p F o o 9 j X Y O g r a y u r W / k N 7 d q S Z j G j F d Z 6 I d x w 3 U S 7 n s B r w p P + L w R x d w Z u j 6 v u 4 M L G a / f 8 D j x w u B K j C P e H j r 9 w O t 5 z a t e x i t > wc = 0.0010 < l a t e x i t s h a 1 _ b a s e 6 4 = " + z p K c w x j l U A l h y M Z h 3 TE / u w v J m c = " > A A A C z 3 i c j V H L S s N A F D 2 N r 1 p f V Z d u g k V w V Z I q 6 E Y o u n H Z g n 1 A W 0 o y n d b Q N A m T i a U U x a 0 / 4 F b / S v w D / Q v v j C m o R X R C k j P n 3 n N m 7 r 1 u 5 H u x t K z X j L G w u L S 8 k l 3 N r a 1 v b G 7 l t 3 f q c Z g I x m s s 9 E P R d J 2 Y + 1 7 A a 9 K T P m 9 G g j s j 1 + c N d 3 i h 4 o 0 b L m I v D K 7 k J O K d k T M I v L 7H H E l U e 9 x l 5 p l p F S 3 L t r r 5 g g J q m f P A T k E B 6 a q E + R e 0 0 U M I h g Q j c A S Q h H 0 4 i O l p w Y a F i L g O p s Q J Q p 6 O c 9 w i R 9 q E s j h l O M Q O 6 T u g X S t l A 9 o r z 1 i r G Z 3 i 0 y t I a e K A N C H l C c L q N F P H E + 2 s 2 N + 8 p 9 p T 3 W 1 C f z f 1 G h E r c U 3 s X 7 p Z 5 n 9 1 q h a J P k 5 1 D R 7 V F G l G V c d S l 0 R 3 R d 3 c / F K V J I e I O I V 7 F B e E m V b O + m x q T a x r V 7 1 1 d P x N Z y p W 7 V m a m + B d 3 Z I G b P 8 c 5 z y o l 4 r 2 U b F U P S 6 U z 9 N R Z 7 G H f R z S P E 9 Q x i U q q J F 3 h E c 8 4 d m o G m P j z r j / T D U y q W Y X 3 5 b x 8 A G 7 4 5 J M < / l a t e x i t > wc = 0.0020 < l a t e x i t s h a 1 _ b a s e 6 4 = " i j 7 k p J t l 9 l L p P 6 O d D a G s a B m T G q w 6 2 c 9 d n U m l j X r n r r 6 P i b z l S s 2 r M 0 N 8 G 7 u i U N 2 P 4 5 z n l Q L x X t o 2 K p e l w o n 6 e j z m I P + z i k e Z 6 g j E t U UC P v C I 9 4 w r N R N c b G n X H / m W p k U s 0 u v i 3 j 4 Q O + R J J N < / l a t e x i t > wi = 0.5 < l a t e x i t s h a 1 _ b a s e 6 4 = " 4 c J w K / 1 d J M j h 1 u Y p k 8 3 F S T K L E d w = " > A A A C z H i c j V H L S s N A F D 2 N r 1 p f V Z d u g k V w F d K q 6 E Y o u n E l Fe x D a i n J d F o H 8 y K Z K K V 0 6 w + 4 1 e 8 S / 0 D / w j t j C m o R n Z D k z L n 3 n J l 7 r x t 5 I p G 2 / Z o z Z m b n 5 h f y i 4 W l 5 Z X V t e L 6 R i M J 0 5 j x O g u 9 M G 6 5 T s I 9 E f C 6 F N L j r S j m j u 9 6 v O n e n q p 4 8 4 7 H i Q i D S z m M e M d 3 B o H o C + Z I o q 7 u u 8 I 8 N m 3 r o F s s 2 Z a t l z k N y h k o I V u 1 s

Fig. 1 .
Fig. 1.Column density (K) maps obtained for filament-like (left) and sheet-like (right) density structures with different wiggle parameters.The wiggle intensity (w i ) is constant in columns and increases from left to right and the wiggle correlation length (w c ) is constant in rows and increases from top to bottom.In our settings, w i and w c have units of [pc/(100 × 200)].The shape parameters of filaments and sheets are the same in each panel and read as follows (in pc).Filaments: L ∈ [7, 25] and R ∈ [1.5, 4.5]; Sheets: L ∈ [5, 20], R 12 ∈ [10, 20] and R 13 ∈ [1, 3] (see Sect. 2.1.4andTable1for parameter definitions and labels.)The circular regions displayed on those maps have an angular radius of 12.53 • .All maps share the same color scale, which is expressed in 10 21 cm −2 .
y c s 9 0 b o p 3 e U s a s P 1 z n L O g X i r a B 8 V S 9 b B Q P t O j z m I H u 9 i n e R 6 j j A t U U C P v C I 9 4 w r N R N W 6 M O + P + M 9 X I a M 0 2 v i 3 j 4 Q O y A J J I < / l a t e x i t > nc = 0.0035 < l a t e x i t s h a 1 _ b a s e 6 4 = " Q a Fig. 2. Stokes Q polarization maps obtained for filaments (left) and sheets (right) with different parameter values for the stochastic component of the magnetic field B sto .The component amplitude (n i ) is constant in column and increases from left to right and the correlation length (n c ) is constant in rows and increases from top to bottom.The circular regions, here centered on the North Galactic pole, have an angular radius of 12.53 • .Shape parameters, including those controlling the wiggles, are the same as in Fig. 3 and are kept fixed between panels.All maps share the same color scale.
1 e N M d n q l 4 8 4 7 H i Q i D K z m K e M d 3 B o H o C + Z I o q 7 v u 8 I 8 M a 2 y 3 S 2 W r L K l l z k L 7 A y U k K 1 a W H z B D X o I w Z D C B 0 c A S d i D g 4 S e N m x Y i I j r Y E x c T E j o O M c E B d K m l M U p w y F 2 S N 8 B 7 d o Z G 9 B e e S Z a z e g U j 9 6 Y l C b 2 S B N S X k X 1 9 Y 3 N g t b 2 4 0 k T G P G 6 y z 0 w r j l O g n 3 R M D r U k i P t 6 K Y O 7 7 r 8 a Y 7 P F P x 5 h2 P E x E G V 3 I U 8 Y 7 v D A L R F 8 y R R F 3 f d 4 V 5 Y t p W p V s o 2 p a t l z k L S h k o I l u 1 s P C C G / Q Q g i G F D 4 4 A k r A H B w k 9 b Z R g I y K u g z F x M S G h 4 x w T 5 E m b U h a n D I f Y I X 0 H t G tn b E B 7 5 Z l o N a N T P H p j U p r Y J 0 1 I e T F h d Z q p 4 6 l 2 V u x v 3 m P t q e 4 2 o r + b e f n E S t w S + 5 d u m v l f n a p F o o 9 j X Y O g r k 7 e M R T 3 g 2 L g x p j I 3 J Z 6 o x l 2 l 2 8 G 0 Z D x 8 A J Z G m < / l a t e x i t > wc = 0.0005 < l a t e x i t s h a 1 _ b a s e 6 4 = " z 7 Y k C k j P n 3 n N m 7 r 1 u 5 H u J s K z X j D E 3 v 7 C 4 l F 3 O r a y u r W / k N 7 d q S Z j G j F d Z 6 I d x w 3 U S 7 n s B r w p P + L w R x d w Z u j 6 v u 4 M L G a / f 8 D j x w u B K j C P e H j r 9 w O t 5 z 5 Z 7 p 3 B T v 8 p Y 0 Y P v n O G d B 7 a B o H x Y P K k e F 0 r k e d R Y 7 2 M U + z f M E J V y i j C p 5 R 3 j E E 5 6 N i j E y 7 o z 7 z 1 Q j o z X b + L a M h w / F Y p J Q < / l a t e x i t > wc = 0.0010 < l a t e x i t s h a 1 _ b a s e 6 4 = " + z p K c w x j l U A l h yM Z h 3 T E / u w v J m c = " > A A A C z 3 i c j V H L S s N A F D 2 N r 1 p f V Z d u g k V w V Z I q 6 E Y o u n H Z g n 1 A W 0 o y n d b Q N A m T i a U U x a 0 / 4 F b / S v w D / Q v v j C m o R X R C k j P n 3 n N m 7 r 1 u 5 H u x t K z X j L G w u L S 8 k l 3 N r a 1 v b G 7 l t 3 f q c Z g I x m s s 9 E P R d J 2 Y + 1 7 A a 9 K T P m 9 G g j s j 1 + c N d 3 i h 4 o 0 b L m I v D K 7 k J O K d k T M I v L 7H H E l U e 9 x l 5 p l p F S 3 L t r r 5 g g J q m f P A T k E B 6 a q E + R e 0 0 U M I h g Q j c A S Q h H 0 4 i O l p w Y a F i L g O p s Q J Q p 6 O c 9 w i R 9 q E s j h l O M Q O 6 T u g X S t l A 9 o r z 1 i r G Z 3 i 0 y t I a e K A N C H l C c L q N F P H E + 2 s 2 N + 8 p 9 p T 3 W 1 C f z f 1 G h E r c U 3 s X 7 p Z 5 n 9 1 q h a J P k 5 1 D R 7 V F G l G V c d S l 0 R 3 R d 3 c / F K V J I e I O I V 7 F B e E m V b O + m x q T a x r V 7 1 1 d P x N Z y p W 7 V m a m + B d 3 Z I G b P 8 c 5 z y o l 4 r 2 U b F U P S 6 U z 9 N R Z 7 G H f R z S P E 9 Q x i U q q J F 3 h E c 8 4 d m o G m P j z r j / T D U y q W Y X 3 5 b x 8 A G 7 4 5 J M < / l a t e x i t > wc = 0.0015 < l a t e x i t s h a 1 _ b a s e 6 4 = " 6 1 1 e N M d n q l 4 8 4 7 H i Q i D K z m K e M d 3 B o H o C + Z I o q 7 v u 8 I 8 M a 2 y 3 S 2 W r L K l l z k L 7 A y U k K 1 a W H z B D X o I w Z D C B 0 c A S d i D g 4 S e N m x Y i I j r Y E x c T E j o O M c E B d K m l M U p w y F 2 S N 8 B 7 d o Z G 9 B e e S Z a z e g U j 9 6 Y l C b 2 S B N S X k X 1 9 Y 3 N g t b 2 4 0 k T G P G 6 y z 0 w r j l O g n 3 R M D r U k i P t 6 K Y O 7 7 r 8 a Y 7 P F P x 5 h2 P E x E G V 3 I U 8 Y 7 v D A L R F 8 y R R F 3 f d 4 V 5 Y t p W p V s o 2 p a t l z k L S h k o I l u 1 s P C C G / Q Q g i G F D 4 4 A k r A H B w k 9 b Z R g I y K u g z F x M S G h 4 x w T 5 E m b U h a n D I f Y I X 0 H t G t n b E B 7 5 Z lo N a N T P H p j U p r Y J 0 1 I e T F h d Z q p 4 6 l 2 V u x v 3 m P t q e 4 2 o r + b e f n E S t w S + 5 d u m v l f n a p F o o 9 j X Y O g m i L N q O p Y 5 p L q r q i b m1 + q k u Q Q E a d w j + I x Y a a V 0 z 6 b W p P o 2 l V v H R 1 / 0 5 m K V X u W 5 a Z 4 V 7 e k A Z d + j n M W N M p W q W K V L w + L 1 d N s 1 D n s Y g 8 H N M 8 j V H G O G u r k 7 e MR T 3 g 2 L g x p j I 3 J Z 6 o x l 2 l 2 8 G 0 Z D x 8 A J Z G m < / l a t e x i t > wc = 0.0005 < l a t e x i t s h a 1 _ b a s e 6 4 = " z 7 Y k C k j P n 3 n N m 7 r 1 u 5 H u J s K z X j D E 3 v 7 C 4 l F 3 O r a y u r W / k N 7 d q S Z j G j F d Z 6 I d x w 3 U S 7 n s B r w p P + L w R x d w Z u j 6 v u 4 M L G a / f 8 D j x w u B K j C P e H j r 9 w O t 5 z Z 7 p 3 B T v 8 p Y 0 Y P v n O G d B 7 a B o H x Y P K k e F 0 r k e d R Y 7 2 M U + z f M E J V y i j C p 5 R 3 j E E 5 6 N i j E y 7 o z 7 z 1 Q j o z X b + L a M h w / F Y p J Q < / l a t e x i t > wc = 0.0010 < l a t e x i t s h a 1 _ b a s e 6 4 = " + z p K c w x j l U A l h yM Z h 3 T E / u w v J m c = " > A A A C z 3 i c j V H L S s N A F D 2 N r 1 p f V Z d u g k V w V Z I q 6 E Y o u n H Z g n 1 A W 0 o y n d b Q N A m T i a U U x a 0 / 4 F b / S v w D / Q v v j C m o R X R C k j P n 3 n N m 7 r 1 u 5 H u x t K z X j L G w u L S 8 k l 3 N r a 1 v b G 7 l t 3 f q c Z g I x m s s 9 E P R d J 2 Y + 1 7 A a 9 K T P m 9 G g j s j 1 + c N d 3 i h 4 o 0 b L m I v D K 7 k J O K d k T M I v L 7H H E l U e 9 x l 5 p l p F S 3 L t r r 5 g g J q m f P A T k E B 6 a q E + R e 0 0 U M I h g Q j c A S Q h H 0 4 i O l p w Y a F i L g O p s Q J Q p 6 O c 9 w i R 9 q E s j h l O M Q O 6 T u g X S t l A 9 o r z 1 i r G Z 3 i 0 y t I a e K A N C H l C c L q N F P H E + 2 s 2 N + 8 p 9 p T 3 W 1 C f z f 1 G h E r c U 3 s X 7 p Z 5 n 9 1 q h a J P k 5 1 D R 7 V F G l G V c d S l 0 R 3 R d 3 c / F K V J I e I O I V 7 F B e E m V b O + m x q T a x r V 7 1 1 d P x N Z y p W 7 V m a m + B d 3 Z I G b P 8 c 5 z y o l 4 r 2 U b F U P S 6 U z 9 N R Z 7 G H f R z S P E 9 Q x i U q q J F 3 h E c 8 4 d m o G m P j z r j / T D U y q W Y X 3 5 b x 8 A G 7 4 5 J M < / l a t e x i t > wc = 0.0015 < l a t e x i t s h a 1 _ b a s e 6 4 = " 6 1 Fig. 3. Same as for Fig. 1 but for the case of non-regular magnetic field and density perturbation coupled to the local magnetic field 3D orientation.The realization of the magnetic-field stochastic component is the same for each panel and corresponds to (n i , n c ) = (0.05, 0.0015).

Fig. 4 .
Fig. 4. From left to right: Intensity (I), Q and U polarization maps and polarization power spectra obtained for filament-like (top) and sheet-like (bottom) clouds.The values of quantities extracted from the power spectra are R EB = 1.69, r T E = 0.65, a T T = −2.47, a EE = −2.45, a BB = −2.47 for filaments and R EB = 1.36, r T E = 0.61, a T T = −2.19, a EE = −2.40, a BB = −1.91 for sheets.
Fig. 5. Dependence of R EB and r T E values on ω, the misalignment angle between magnetic field's orientation and structures' main axes for both filament and sheet structures' shapes.Top row shows R EB (left) and r T E (right) versus ω for α = 90 • for both filaments (blue) and sheets (orange).The shaded areas include 68% of the data around the median (thick lines) obtained through 10 random realizations of the same setup.Bottom row shows scatter plots of (R EB , r T E ) pairs for filaments (left) and sheets (right) color-coded according to the ω value and obtained for the complete set of synthetic map.For visual reference, the black star in the bottom row indicate the fiducial Planck values (R EB = 1.89, r T E = 0.36).
C X J I S J O 4 Q 7 F B W G m l e M + m 1 o T 6 9 p V b x 0 d f 9 O Z i l V 7 l u Y m e F e 3 p A H b P 8 c 5 C e q l o r 1 T L J 3 u F s o H 6 a i z 2 M A m t m m e e y j j G B X U y P s K j 3 j C s 3 F i 3 B h 3 x v 1 n q p F J N e v 4 t o y H D 7 j 5 k Z w = < / l a t e x i t > R EB < l a t e x i t s h a 1 _ b a s e 6 4 = " p c I j W K m M T w s P P r M o f c w B T A s 4 F c 8 D n 3 n j N z 7 3 U i z 0 2 k b b / m j L n 5 h c W l / H J h Z X V t f c P c 3 G o m Y R p z 0 e C h F 8 Z t h y X C c w P R k K 7 0 R D u K B f M d T 7 S c w Y m K t 4 Y i T t w w u J S j S F z 7 7 C Z w + y 5 n k q i u a V 7 5 T N 5 y 5 o 0 v J t 3 x a W 3 S N Y t 2 y d b L m g X l D B S R r X p o v u A K P Y T g S O F D I I A k 7 I E h o a e D M m x E x F 1 j T F x M y N V x g Q k K p E 0 p S 1 A G I 3 Z A 3 x v a d T I 2 o L 3 y T L S a 0 y k e v T E p L e y R J q S 8 m L A 6 z d L x V D s r 9 j f v s f Z U d x v R 3 8 m 8 f G I l b o n 9 S z f N / K 9 O 1 S L R x 7 G uw a W a I s 2 o 6 n j m k u q u q J t b X 6 q S 5 B A R p 3 C P 4 j F h r p X T P l t a k + j a V W + Z j r / p T M W q P c 9 y U 7 y r W 9 K A y z / H O Q u a l V L 5 o F Q 5 P y x W a 9 m o 8 9 j B L v Z p n k e o 4 g x 1 N M h 7 i E c 8 4 d l o G R P j z r j / T D V y m W Y b 3 5 b x 8 A F p f J Y u < / l a t e x i t > R EB < l a t e x i t s h a 1 _ b a s e 6 4 = " p c I j W K m M T w s P P r M o f c w B T A s 4 F c 8 D n 3 n j N z 7 3 U i z 0 2 k b b / m j L n 5 h c W l / H J h Z X V t f c P c 3 G o m Y R p z 0 e C h F 8 Z t h y X C c w P R k K 7 0 R D u K B f M d T 7 S c w Y m K t 4 Y i T t w w u J S j S F z 7 7 C Z w + y 5 n k q i u a V 7 5 T N 5 y 5 o 0 v J t 3 x a W 3 S N Y t 2 y d b L m g X l D B S R r X p o v u A K P Y T g S O F D I I A k 7 I E h o a e D M m x E x F 1 j T F x M y N V x g Q k K p E 0 p S 1 A G I 3 Z A 3 x v a d T I 2 o L 3 y T L S a 0 y k e v T E p L e y R J q S 8 m L A 6 z d L x V D s r 9 j f v s f Z U d x v R 3 8 m 8 f G I l b o n 9 S z f N / K 9 O 1 S L R x 7 G u w a W a I s 2 o 6 g p K C B d l T D / g j N 0 E I I h Q R 8 c A S R h H w 5 i e l q w Y S E i r o 0 h c Y K Q p + M c I + R I m 1 A W p w y H 2 G v 6 9 m j X S t m A 9 s o z 1 m pG p / j 0 C l K a 2 C J N S H m C s D r N 1 P F E O y v 2 N + + h 9 l R 3 G 9 D f T b 3 6 x E p c E v u X b p z 5 X 5 2 q R a K L f V 2 D R z V F m l H V s d Q l 0 V 1 R N z e / VC X J I S J O 4 Q 7 F B W G m l e M + m 1 o T 6 9 p V b x 0 d f 9 O Z i l V 7 l u Y m e F e 3 p A H b P 8 c 5 C e q l o r 1 T L J 3 u F s o H 6 a i z 2 M A m t m m e e y j j G B X U y P s K j 3 j C s 3 F i 3 B h 3 x v 1 n q p F J N e v 4 t o y H D 7 j 5 k Z w = < / l a t e x i t > C X J I S J O 4 Q 7 F B W G m l e M + m 1 o T 6 9 p V b x 0 d f 9 O Z i l V 7 l u Y m e F e 3 p A H b P 8 c 5 C e q l o r 1 T L J 3 u F s o H 6 a i z 2 M A m t m m e e y j j G B X U y P s K j 3 j C s 3 F i 3 B h 3 x v 1 n q p F J N e v 4 t o y H D 7 j 5 k Z w = < / l a t e x i t > F ilaments < l a t e x i t s h a 1 _ b a s e 6 4 = " Z z J d 0 I b 5 Z d g K 6 a k 2 + 8 X V m 5 r p s B 4 8 0 M e U J w u o 0 W 8 c z 7 a z Y 3 7 z H 2 l P d bU R / z 3 i F x E r c E f u X b p r 5 X 5 2 q R a K P Y 1 0 D p 5 o S z a j q f O O S 6 a 6 o m 9 t f q p L k k B C n c I / i g r C v l d M + 2 1 q T 6 t p V b 1 0 d f 9 O Z i l V 7 3 + R m e F e 3 p A F X f o 5 z F j S r 5 c p B u X p 5 W K q d m F H n s Y N d 7 N M 8 j 1 D D O e p o k H e I R z z h 2 b q w p D W 2 J p + p V s 5 o t v F t W Q 8 f L v i S / A = = < / l a t e x i t >R EB < l a t e x i t s h a 1 _ b a s e 6 4 = " p c I j W K m M T w s P P r M o f c w B T A s 4 F c 8 D n 3 n j N z 7 3 U i z 0 2 k b b / m j L n 5 h c W l / H J h Z X V t f c P c 3 G o m Y R p z 0 e C h F 8 Z t h y X C c w P R k K 7 0 R D u K B f M d T 7 S c w Y m K t 4 Y i T t w w u J S j S F z 7 7 C Z w + y 5 n k q i u a V 7 5 T N 5 y 5 o 0v J t 3 x a W 3 S N Y t 2 y d b L m g X l D B S R r X p o v u A K P Y T g S O F D I I A k 7 I E h o a e D M m x E x F 1 j T F x M y N V x g Q k K p E 0 p S 1 A G I 3 Z A 3 x v a d T I 2 o L 3 y T L S a 0 y k e v T E p L e y R J q S 8 m L A 6 z d L x V D s r 9 j f v s f Z U d x v R 3 8 m 8 f G I l b o n 9 S z f N / K 9 O 1 S L R x 7 G uw a W a I s 2 o 6 n j m k u q u q J t b X 6 q S 5 B A R p 3 C P 4 j F h r p X T P l t a k + j a V W + Z j r / p T M W q P c 9 y U 7 y r W9 K A y z / H O Q u a l V L 5 o F Q 5 P y x W a 9 m o 8 9 j B L v Z p n k e o 4 g x 1 N M h 7 i E c 8 4 d l o G R P j z r j / T D V y m W Y b 3 5 b x 8 A F p f J Y u < / l a t e x i t > Sheets < l a t e x i t s h a 1 _ b a s e 6 4 = " e R j s x U T z l a c I S r U S d M 8 H a Y 5 K X Q 0 = " > A A A C y X i c j V H L S s N A F D 2 N r 1 p f V Z d u g k V w V Z I q 6 L L o R n B T 0 T 6 g F k n S a T s 2 L 5 O J W I s r f 8 C t / p j 4 B / o X 3 h l H U I v o h C R n z r 3 n z N x 7 3 d j n q b C s l 5 w x N T 0 z O 5 e f L y w s L i 2 v F F f X G m m U J R 6 r e 5 E f J S 3 X S Z n P Q 1 Y X X P i s F S f M C V y f N d 3 h o Y w 3 r 1 m S 8 i g 8 E 6 O Y d Q K n H / I e 9 x x B V O N 0 w J h I L 4 o l q 2 y p Z U 4 C W 4 M S 9 K p F x W e c o 4 s I H j I E Y A g h C P t w k N L T h g 0 L M X E d j I l L C H E V Z 7 h D g b Q Z Z T H K c I g d 0 r d P u 7 Z m Q 9 p L z 1 S p P T r F p z c h p Y k t 0 k S U l x C W p 5 k q n i l n y f 7 m P V a e 8 m 4 j + r v a K y B W Y E D s X 7 r P z P / q Z C 0 C P e y r G j j V F C t G V u d p l 0 x 1 R d 7 c / F K V I I e Y O I m 7 F E 8 I e 0 r 5 2 W d T a V J V u + y t o + K v K l O y c u / p 3 A x v 8 p Y 0 Y P v n O C d B o 1 K 2 d 8 q V k 9 1 S 9 U C P O o 8 N b G K b 5 r m H K o 5 Q Q 5 2 8 L / G A R z w Z x 8 a V c W P c f q Q a O a 1 Z x 7 d l 3 L 8 D 0 r W R p w = = < / l a t e x i t > ↵ < l a t e x i t s h a 1 _ b a s e 6 4 = " s u 9 O u v c t / 8 / 1 x j w 1 F x p 9 X 1 R W 9 Z c = " > A A A C y X i c j V H L T s J A F D 3 U F + I L d e m m k Z i 4 I i 2 a 6 J L o x s Q N J g I m Q M x0 G K D S l + 3 U i M S V P + B W f 8 z 4 B / o X 3 h l L o h K j 0 7 Q 9 c + 4 9 Z + b e 6 0 S e m 0 j L e s 0 Z M 7 N z 8 w v 5 x c L S 8 s r q W n F 9 o 5 G E a c x F n Y d e G F 8 4 L B G e G 4 i 6 d r n p i r 6 5 8 6 U q R Q 4 p c R p 3 K S 4 J M 6 M c 9 d k x m s z U r n v r m / i b y d S s 3 r M i N 8 e 7 v i U N 2 P s 5 z n F Q 2 6 l 4 u 5 W d 8 7 3 y 0 X E x 6 j l s Y B P b N M 9 9 H O E U Z 6 i S d x + P e M K z V b d u r T v r / j P V m i g 0 6 / i 2 r I c P J W K V P Q = = < / l a t e x i t > != 0 < l a t e x i t s h a 1 _ b a s e 6 4 = " m 7 b k Z s S 6 u t l B o B 3 s S q 0 6 k w b o W s

Fig. 7 .
Fig. 7. Contour plots of the 2D distributions of (R EB , r T E ) from filament-like (blue) and sheet-like (orange) clouds.The embedded contours enclose about 10, 70 and 90 per cent of the data points.To guide the eye, the black star indicates the fiducial Planck values (R EB = 1.89, r T E = 0.36).

Fig. 8 .
Fig. 8. Histograms of spectral indices of the power-law power spectra fitted to the E (orange) and B (green) auto-power spectra obtained for filament-like clouds (left) and sheet-like clouds (right); the full sample is included.

Fig. 9 .
Fig.9.Same as Fig.7but with noisy magnetic field and correlated perturbations in density structures.

Fig. 10 .
Fig.10.Histograms of r T B obtained for filament-like clouds (blue) and sheet-like clouds (orange) for regular magnetic field (left) and with "turbulence" switched on (right); the full sample is included.