arxiv: v2 [astro-ph.ga] 23 Feb 2017

Transcription

1 Astronomy & Astrophysics manuscript no. GASS_ms2 c ESO 2018 Juy 13, 2018 odeing and simuation of arge-scae poarized dust emission over the southern Gaactic cap using the GASS Hi data T. Ghosh 1, 2, F. Bouanger 2, P. G. artin 3, A. Bracco 4, 2, F. Vansynge 2, J. Aumont 2, J. J. Bock 1, 5, O. Doré 1, 5, U. Haud 6, P.. W. Kabera 7, and P. Serra 1, 5 arxiv: v2 [astro-ph.ga] 23 Feb Caifornia Institute of Technoogy, Pasadena, Caifornia, U.S.A. 2 Institut d Astrophysique Spatiae, CNRS (UR8617) Université Paris-Sud 11, Bâtiment 121, Orsay, France 3 CITA, University of Toronto, 60 St. George St., Toronto, ON 5S 3H8, Canada 4 Laboratoire AI, IRFU/Service d Astrophysique - CEA/DS - CNRS - Université Paris Diderot, Bât. 709, CEA-Sacay, F Gif-sur-Yvette Cedex, France 5 Jet Propusion Laboratory, Caifornia Institute of Technoogy, 4800 Oak Grove Drive, Pasadena, Caifornia, U.S.A. 6 Tartu Observatory, Tõravere, Tartumaa, Estonia 7 Argeander-Institut für Astronomie, Universität Bonn, Auf dem Hüge 71, D Bonn, Germany Received...; accepted... ABSTRACT The Panck survey has quantified poarized Gaactic foregrounds and estabished that they are a main imiting factor in the quest for the cosmic microwave background (CB) B-mode signa induced by primordia gravitationa waves during cosmic infation. Accurate separation of the Gaactic foregrounds therefore binds this quest to our understanding of the magnetized interstear medium (IS). The two most reevant empirica resuts from anaysis of Panck data are ine of sight depoarization arising from fuctuations of the Gaactic magnetic fied orientation and aignment of fiamentary dust structures with the magnetic fied at high Gaactic atitude. Furthermore, Panck and Hi emission data in combination indicate that most of the fiamentary dust structures are in the cod neutra medium. The goa of this paper is to test whether these saient observationa resuts, taken together, can account fuy for the statistica properties of the dust poarization over a seected ow coumn density region comprising 34 % of the southern Gaactic cap (b 30 ). To do this, we construct a dust mode that incorporates Hi coumn density maps as tracers of the dust intensity structures and a phenomenoogica description of the Gaactic magnetic fied. By adjusting the parameters of the dust mode, we were abe to reproduce the Panck dust observations at 353 GHz in the seected region. Reaistic simuations of the poarized dust emission enabed by such a dust mode are usefu for testing the accuracy of component separation methods, studying non-gaussianity, and constraining the amount of decorreation with frequency. Key words. Interstear medium: dust, observations IS: structure IS: magnetic fieds poarization 1. Introduction An intense focus in current observationa cosmoogy is detection of the CB B-mode signa induced by primordia gravitationa waves during the infation era (Starobinsky 1979; Fabbri & Poock 1983; Abbott & Wise 1984). The BICEP2 experiment 1 puts an upper bound on the ampitude of the CB B-mode signa, parameterized by a tensorto-scaar ratio (r) at a eve of r < 0.09 (95 % confidence eve (BICEP2 and Keck Array Coaboration et a. 2016). Therefore, the major chaenge for a BICEP2-ike experiment with imited frequency coverage is to detect such an incrediby faint CB B-mode signa in the presence of foreground Gaactic contamination (Betoue et a. 2009). Corresponding author: tghosh@catech.edu 1 The dominant poarized foreground above 100 GHz come from therma emission by aigned aspherica dust grains (Panck Coaboration Int. XXII 2015; Panck Coaboration X 2016). Unike the CB, the poarized dust emission is distributed non-uniformy on the sky with varying coumn density over which are summed contributions from mutipe components with different dust composition, size, and shape and with different magnetic fied orientation (see reviews by Prunet & Lazarian 1999; Lazarian 2008). Acknowedging these compexities, the goa of this paper is to mode the principa effects, namey mutipe components with different magnetic fied orientations, to produce reaistic simuated maps of the poarized dust emission. Our knowedge and understanding of the dust poarization has improved significanty in the submiimetre and microwave range through expoitation of data from the Panck Artice number, page 1 of 19

2 A&A proofs: manuscript no. GASS_ms2 sateite 2 (Panck Coaboration I 2016; Panck Coaboration X 2016; Panck Coaboration Int. XIX 2015), which provides important empirica constraints on the modes. The Panck sateite has mapped the poarized sky at seven frequencies between 30 and 353 GHz (Panck Coaboration I 2016). The Panck maps are avaiabe in HEALPix 3 format (Górski et a. 2005) with resoution parameter abeed with the N side vaue. The genera statistica properties of the dust poarization in terms of anguar power spectrum at intermediate and high Gaactic atitudes from 100 to 353 GHz are quantified in Panck Coaboration Int. XXX (2016). Two unexpected resuts from the Panck observations were an asymmetry in the ampitudes of the anguar power spectra of the dust E- and B- modes (EE and BB spectra, respectivey) and a positive temperature E-mode (T E) correation at 353 GHz. Expanation of these important features was expored by Panck Coaboration Int. XXXVIII (2016) through a statistica study of the fiamentary structures in the Panck Stokes maps and Cadwe et a. (2016) discussed them in the context of magnetohydrodynamic interstear turbuence. They are a centra focus of the modeing in this paper. The Panck 353 GHz poarization maps (Panck Coaboration Int. XIX 2015) have the best signa-to-noise ratio and enabe study of the ink between the structure of the Gaactic magnetic fied (GF) and dust poarization properties. A arge scatter of the dust poarization fraction (p) for tota coumn density N H < cm 2 at ow and intermediate Gaactic atitudes is reported in Panck Coaboration Int. XIX (2015). Based on numerica simuations of anisotropic magnetohydrodynamic turbuence in the diffuse interstear medium (IS), it has been concuded that the arge scatter of p in the range cm 2 < N H < cm 2 comes mosty from fuctuations in the GF orientation aong the ine of sight (LOS) rather than from changes in grain shape or the efficiency of the grain aignment (Panck Coaboration Int. XX 2015). To account for the observed scatter of p in the high atitude southern Gaactic sky (b 60 ), a finite number of poarization ayers was introduced (Panck Coaboration Int. XLIV 2016). The orientation of arge-scae GF with respect to the pane of the sky (POS) aso pays a crucia roe in expaining the scatter of p (Panck Coaboration Int. XX 2015). The Panck poarization data revea a tight reationship between the orientation of the dust intensity structures and the magnetic fied projected on the pane of the sky (B POS, Panck Coaboration Int. XXXII 2016; Panck Coaboration Int. XXXV 2016; Panck Coaboration Int. XXXVIII 2016). In particuar, the dust intensity structures are preferentiay aigned with B POS in the diffuse IS (Panck Coaboration Int. XXXII 2016). Simiar aignment is aso reported between the Gaactic Arecibo L-Band Feed Array (GALFA, Peek et a. 2011) Hi fiaments and B POS derived either from staright poarization (Cark et a. 2014) 2 Panck ( is a project of the European Space Agency (ESA) with instruments provided by two scientific consortia funded by ESA member states and ed by Principa Investigators from France and Itay, teescope refectors provided through a coaboration between ESA and a scientific consortium ed and funded by Denmark, and additiona contributions from NASA (USA). 3 or Panck poarization data (Cark et a. 2015). Furthermore, a detaied a-sky study of the Hi fiaments, combining the Gaactic A Sky Survey (GASS, ccure-griffiths et a. 2009) and the Effesberg Bonn Hi Survey (EBHIS, Winke et a. 2016), shows that most of the fiaments occur in the cod neutra medium (CN, Kabera et a. 2016). The Hi emission in the CN phase has a ine profie with 1σ veocity dispersion 3 km s 1 (FWH 7 km s 1 ) in the soar neighbourhood (Heies & Troand 2003) and fiamentary structure on the sky, which might resut from projection effects of gas organized primariy in sheets (Heies & Troand 2005; Kabera et a. 2016). Aignment of CN structures with B POS from Panck has been reported towards the north eciptic poe, among the targeted fieds of the Green Bank Teescope Hi Intermediate Gaactic Latitude Survey (GHIGLS, artin et a. 2015). Such aignment can be induced by shear strain of gas turbuent veocities stretching matter and the magnetic fied in the same direction (Hennebee 2013; Inoue & Inutsuka 2016). Aignment between the CN structures and B POS was not incuded in the pre-panck dust modes, for exampe Panck Sky ode (PS, Deabrouie et a. 2013) and FG- Po (O Dea et a. 2012). The most recent version of post- Panck PS uses a fitered version of the Panck 353 GHz poarization data as the dust tempates (Panck Coaboration XII 2016). Therefore, in ow signa-to-noise regions the PS dust tempates are not reiabe representations of the true dust poarized sky. In Panck Coaboration Int. XXXVIII (2016) it was shown that preferentia aignment of the dust fiaments with B POS coud account for both the E-B asymmetry in the ampitudes of the dust D EE and D BB anguar power spectra 4 over the mutipoe range 40 < < 600 and the ampitude ratio D T E /D EE, as measured by Panck Coaboration Int. XXX (2016). Cark et a. (2015) showed that the aignment of the B POS with the Hi structures can expain the observed E-B asymmetry. Recenty, Vansynge et a. (2016) have extended the approach introduced in Panck Coaboration Int. XLIV (2016) using a parametric mode to account for the observed dust power spectra at intermediate and high Gaactic atitudes. Our work in this paper is compementary, incorporating additiona astrophysica data and insight and focusing on the ceanest sky area for CB B-mode studies from the southern hemisphere. In more detai, the goa of this paper is to test whether preferentia aignment of the CN structures with B POS together with fuctuations in the GF orientation aong the LOS can account fuy for a number of observed statistica properties of the dust poarization over the southern ow coumn density sky. These statistica properties incude the scatter of the poarization fraction p, the dispersion of the poarization ange ψ around the arge-scae GF, and the observed dust power spectra D EE, D BB, and D T E at 353 GHz. In doing so we extend the anaysis of Panck Coaboration Int. XXXVIII (2016) and Panck Coaboration Int. XLIV (2016) to sky areas in which the fiaments have very itte contrast with respect to the diffuse background emission. Within the southern Gaactic cap (defined here to be b 30 ), in which dust and gas are we correated (Bouanger et a. 1996; Panck Coaboration XI 2014) and the GASS Hi emission is a good tracer of the dust intensity structures (Panck Coaboration Int. XVII 2014), we 4 D EE,BB ( + 1)C EE,BB /(2π). Artice number, page 2 of 19

3 T. Ghosh et a.: odeing and simuation of poarized dust emission focus our study on the ow coumn density portion. We choose to mode the dust poarization sky at 1 resoution, which corresponds to retaining mutipoes 160 where the Gaactic dust dominates over the CB B-mode signa. Aso, most of the ground-based CB experiments are focusing at a degree scae for the detection of the recombination bump of the CB B-mode signa. Our dust poarization mode shoud be a usefu product for the CB community. It covers 34 % of the southern Gaactic cap, a 3500 deg 2 region containing the ceanest sky accessibe to ongoing ground and baoon-based CB experiments: ABS (Essinger-Hieman et a. 2010), Advanced ACTPo (Henderson et a. 2016), BICEP/Keck (Ogburn et a. 2012), CLASS (Essinger-Hieman et a. 2014), EBEX (Reichborn- Kjennerud et a. 2010), Poarbear (Kermish et a. 2012), Simons array (Suzuki et a. 2016), SPIDER (Cri et a. 2008), and SPT3G (Benson et a. 2014). This paper is organized as foows. Section 2 introduces the Panck observations and GASS Hi data used in our anaysis, aong with a seected sky region within the southern Gaactic cap, caed SGC34. In Sect. 3, we summarize the observationa dust poarization properties in SGC34 as seen by Panck. A mode to simuate dust poarization using the GASS Hi data and a phenomenoogica description of the GF is detaied in Sect. 4. The procedure to simuate the poarized dust sky towards the southern Gaactic cap and the choice of the dust mode parameters are discussed in Sect. 5. Section 5.4 presents the main resuts where we statisticay compare the dust mode with the observed Panck data. In Sect. 6, we make some predictions using our dust mode. We end with an astrophysica perspective in Sect. 7 and a short discussion and summary of our resuts in Sect Data sets used and region seection 2.1. Panck 353 GHz data In this paper we used the pubicy avaiabe 5 maps from the Panck 2015 poarization data reease (PR2), specificay at 353 GHz (Panck Coaboration I 2016). 6 The data processing and caibration of the Panck High Frequency Instrument (HFI) maps are described in Panck Coaboration VIII (2016). We used mutipe Panck 353 GHz poarization data sets: the two yeary surveys, Year 1 (Y1) and Year 2 (Y2), the two haf-mission maps, H1 and H2, and the two detector-set maps, DS1 and DS2, to expoit the statistica independence of the noise between them and to test for any systematic effects. 7 The 353 GHz poarization maps are corrected for zodiaca ight emission at the map-making eve (Panck Coaboration VIII 2016). These data sets are described in fu detai in Panck Coaboration VIII (2016). A of the 353 GHz data sets used in this study have an anguar resoution of FWH and are projected on a HEALPix grid with a resoution parameter N side = 2048 ( pixes). We corrected for the main systematic effect, that is eakage from intensity to poarization, using the goba tempate fit, as described in Panck Coaboration These maps are expressed in thermodynamic units (K CB). For this passband, mutipy by to convert to Jy sr 1. 7 The noise is partiay correated between the two haf-ring maps (Panck Coaboration VI 2014) and so those maps are not used in our anaysis. VIII (2016), which accounts for the monopoe, caibration, and bandpass mismatches. To increase the signa-to-noise ratio of the Panck poarization measurements over SGC34, we smoothed the Panck poarization maps to 1 resoution using a Gaussian approximation to the Panck beam and reprojected to a HEALPix grid with N side =128 (30 pixes). We did not attempt to correct for the CB poarization anisotropies at 1 resoution because they are negigiby sma compared to the dust poarization at 353 GHz. The dust intensity at 353 GHz, D 353, is cacuated from the atest pubicy avaiabe Panck therma dust emission mode (Panck Coaboration Int. XLVIII 2016), which is a modified backbody (BB) fit to maps at ν 353 GHz from which the cosmic infrared background anisotropies (CIBA) have been removed using the generaized inear interna combination (GNILC) method and zero-eve offsets (e.g. the CIB monopoe) have been removed by correation with Hi maps at high Gaactic atitude. By construction, this dust intensity map is corrected for zodiaca ight emission, CB anisotropies, and resoved point sources as we as the CIBA and the CIB monopoe. As a resut of the GNILC method and the BB fitting, the noise in the D 353 map is ower than in the Stokes I frequency map at 353 GHz. We aso smoothed the D 353 map to the common resoution of 1, taking into account the effective beam resoution of the map GASS Hi data We use the GASS Hi survey 8 (ccure-griffiths et a. 2009; Kabera et a. 2010; Kabera & Haud 2015) which mapped the southern sky (decination δ < 1 ) with the Parkes teescope at an anguar resoution of FWH The rms uncertainties of the brightness temperatures per channe are at a eve of 57 mk at δv = 1 km s 1. The GASS survey is the most sensitive and highest anguar resoution survey of Hi emission over the southern sky. The third data reease version has improved performance by the remova of residua instrumenta probems aong with stray radiation and radio-frequency interference, which can be major sources of error in Hi coumn density maps (Kabera & Haud 2015). The GASS Hi data are projected on a HEALPix grid with N side = 1024 (3. 75 pixes). For this opticay thin gas, the brightness temperature spectra T b of the GASS Hi data are decomposed into mutipe Gaussian components (Haud & Kabera 2007; Haud 2013) as T b = [ T0 i exp 1 ( vlsr v i ) 2 ] c 2 σ i, (1) i where the sum is over a Gaussian components for a LOS, v LSR is the veocity of the gas reative to the oca standard of rest (LSR), and T i 0, v i c, and σ i are the peak brightness temperature (in K), centra veocity, and 1σ veocity dispersion of the i th Gaussian component, respectivey. Athough the Gaussian decomposition does not provide a unique soution, it is correct empiricay. Haud (2013) used the GASS survey at N side =1024 and decomposed the observed Hi profies into Gaussians. In this study, we use an updated (but unpubished) Gaussian decomposition of the atest version of the GASS Hi survey 8 Artice number, page 3 of 19

4 A&A proofs: manuscript no. GASS_ms2 0 o 0 o 45 o 315 o 45 o 315 o 90 o -60 o 270 o 90 o -60 o 270 o -30 o -30 o 135 o 225 o 135 o 225 o 180 o 180 o 0 2x10 20 cm x10 20 cm -2 0 o 0 o 45 o 315 o 45 o 315 o 90 o -60 o 270 o 90 o -60 o 270 o -30 o -30 o 135 o 225 o 135 o 225 o 180 o 0 2x10 20 cm o Fig. 1. N Hi maps of the three Gaussian-based Hi phases: CN (top eft), UN (top right), and WN (bottom eft) over the portion of the southern Gaactic cap (b 30 ) covered by the GASS Hi survey, at 1 resoution. Cooured region corresponds to the seected region SGC34 (bottom right). Orthographic projection in Gaactic coordinates: atitudes and ongitudes marked by ines and circes respectivey. (third data reease). The intermediate and high veocity gas over the southern Gaactic cap has negigibe dust emission (see discussion in Panck Coaboration Int. XVII 2014) and in any case ines of sight with uncorreated dust and Hi emissions are removed in our masking process (Sect. 2.3). Broady speaking, the Gaussian Hi components represent the three conventiona phases of the IS. In order of increasing veocity dispersion these are: (1) a cod dense phase, the cod neutra medium (CN), with a temperature in the range K; (2) a thermay unstabe neutra medium (UN), with a temperature in the range K; and (3) a warm diffuse phase, the warm neutra medium (WN), with a temperature in the range K. Based on sensitive high-resoution Arecibo observations, Heies & Troand (2003) found that a significant fraction of the Hi emission is in the UN phase at Artice number, page 4 of 19

5 T. Ghosh et a.: odeing and simuation of poarized dust emission high Gaactic atitude. Numerica simuations of the IS suggest that the magnetic fied and dynamica processes ike turbuence coud drive Hi from the stabe WN and CN phases into the UN phase (Audit & Hennebee 2005; Hennebee & Iffrig 2014; Saury et a. 2014). The tota Hi coumn density aong the LOS, N Hi, in units of cm 2, is obtained by integrating the brightness temperature over veocity. We can separate the tota N Hi into the three different phases of the IS based on the veocity dispersion of the fitted Gaussian. Thus N Hi = 1.82 T b dv LSR [ = 1.82 T0 i exp 1 ( vlsr v i ) 2 ] c 2 σ i dv LSR i = π ( T i 0 σ i w c + T0 i σ i w u + T0 i σ i ) w w N c Hi + N u Hi + N w Hi, (2) where w c, w u, and w w are the weighting factors to seect the CN, UN, and WN components, respectivey, and NHi c, N Hi u, and N Hi w are the CN, UN, and WN coumn densities, respectivey. The weighting factors are defined as 1 [ ( )] σ < σ c 1 π σ σ w c = cos c g σ c σ c σ σ c (1 + g) w w = 0 σ > σ c (1 + g) 0 [ ( σ < σ u )] 1 cos u 1 2 π g σ σ σ u σ u σ < σ u (1 + g) 1 σ σ u (1 + g) w u = 1 w c w w. We use overapping weighting factors to avoid abrupt changes in the assignment of the emission from a given Hi phase to another, adopting g = 0.2. There is a partia correation between two neighbouring Hi components. To seect oca veocity gas, we incuded ony Gaussian components with v i c 50 km s 1. For our anaysis we aso smoothed these N Hi maps to the common resoution of 1, taking into account the effective beam resoution of the map. Figure 1 shows the separation of the tota N Hi into the three Gaussian-based Hi phases, CN, UN, and WN. The weights were cacuated with σ c = 7.5 km s 1 and σ u = 10 km s 1 as adopted in Sect The tota N Hi map obtained by summing these component phases is very cose to the LVC map in Fig. 1 of Panck Coaboration Int. XVII (2014) obtained by integrating over the veocity range defined by Gaactic rotation, with dispersion ony cm 2 in the difference or 0.8 % in the fractiona difference over the seected region SGC34 defined beow. To buid our CN map we use an upper imit to the veocity dispersion (σ c = 7.5 km s 1 ) arger than the conventiona ine width of cod Hi gas (Wofire et a. 2003; Heies & Troand 2005; Kabera et a. 2016). The ine width can be characterised by a Dopper temperature T D. Because the ine width incudes a contribution from turbuence, T D is an upper imit to the gas kinetic temperature. In Kabera et a. (2016), the median Dopper temperature of CN gas is about T D = 220 K, corresponding to a veocity dispersion of about σ c = 1.3 km s 1, whie the argest vaue of (3) T D = 1100 K corresponds to σ c = 3.0 km s 1. The vaue of σ c = 7.5 km s 1 in our mode foows from the mode fit to the Panck poarization data (the observed E B asymmetry and the D T E /D EE ratio over SGC34, see Sect. 5.2). Thus our CN component (NHi c ) incudes some UN gas. To investigate the impact of this on our modeing, we divided our CN component map into two parts corresponding to components with veocity dispersions above and beow σ c = 3.0 km s 1. We then repeated our data fitting with these two separate maps repacing our singe CN component map and found no significant difference in the mode resuts Seection of the region SGC34 For our anaysis we seect a ow dust coumn density region of the sky within the southern Gaactic cap. As a first cut, we retain ony the sky pixes that have N Hi beow a threshod vaue, N Hi cm 2. oecuar gas is known from UV observations to be negigibe at this threshod N Hi vaue (Gimon et a. 2006). We adopt a second mask, from Panck Coaboration Int. XVII (2014), to avoid sky pixes that fa off the main trend of gas-dust correation between the GASS Hi coumn density and Panck HFI intensity maps, in particuar appying a threshod of 0.21 Jy sr 1 (3σ cut, see Figure 4 of Panck Coaboration Int. XVII 2014) on the absoute vaue of the residua emission at 857 GHz after subtraction of the dust emission associated with Hi gas. We smooth the second mask to Gaussian 1 FWH resoution and downgrade to N side =128 HEALPix resoution. Within this sky area, we aso mask out a 10 radius patch centred around (, b) = ( , ) to avoid Hi emission from the ageanic stream, whose mean radia veocity for this specific sky patch is within the Gaactic range of veocities (Nidever et a. 2010). After masking, we are eft with a 3500 deg 2 region comprising 34 % of the southern Gaactic cap, referred to hereafter as SGC34. This is presented in Fig. 1. The mean N Hi over SGC34 is N Hi = cm 2. For the adopted vaues of σ c and σ u, the mean coumn densities of the three IS phases (Fig. 1) are NHi c = cm 2, NHi u = cm 2, and NHi w = cm 2. SGC34 has a mean dust intensity of D 353 = 54 kjy sr 1. Correating the tota N Hi map with the D 353 map, we find an emissivity (sope) 32 kjy sr 1 (10 20 cm 2 ) 1 at 353 GHz (Sect. 5.2) and an offset of just 3 kjy sr 1, which we use beow as a 1σ uncertainty on the zero-eve of the Gaactic dust emission. For computations of the anguar power spectrum, we apodized the seection mask by convoving it with a 2 FWH Gaussian. The effective sky coverage of SGC34 as defined by the mean sky coverage of the apodized mask is fsky eff = (8.5 %), meaning that the smoothing extends the sky region sighty but conserves the effective sky area Power spectrum estimator We use the pubicy avaiabe Xpo code (Tristram et a. 2005) to compute the binned anguar power spectra over a given mask. A top-hat binning in intervas of 20 over the range 40 < < 160 is appied throughout this anaysis. The upper cutoff max = 160 is set by 1 smoothed maps used throughout this paper. The ower cutoff min = 40 is cho- Artice number, page 5 of 19

6 A&A proofs: manuscript no. GASS_ms2 Pixe density Fig. 2. Normaized histograms of p 2 over SGC34 at 353 GHz. Panck data, p 2 d, bue inverted-trianges, with 1σ error bars as described in text. The pixe density on the y-axis is the number of pixes in each bin of p 2 divided by the tota number of pixes over SGC34. Dust mode, p 2 m, back circes connected by ine, with 1σ error bars computed from 100 onte-caro reaizations. The vertica dashed ine corresponds to the parameter p 0 = 18.5 % (p 2 0 = 0.034) of this dust mode. sen to avoid ow systematic effects present in the pubicy avaiabe Panck poarization data (Panck Coaboration Int. XLVI 2016). The Xpo code corrects for the incompete sky coverage, eakage from E- to B-mode poarization, beam smoothing, and pixe window function. Other codes such as PoSpice (Chon et a. 2004) and Xpure (Grain et a. 2009) aso compute anguar power spectra from incompete sky coverage. We prefer Xpo over other codes because it directy computes the anaytica error bars of the power spectrum without invoving any onte-caro simuations. The fina resuts of this paper do not depend on the Xpo code, as the same mean power spectra are returned by other power spectrum estimators. 3. Panck dust poarization observations over the seected region SGC34 In this section we derive the statistica properties of the Panck dust poarization over SGC34. These statistica properties incude normaized histograms from pixe space data, specificay of p 2 and of ψ with respect to the oca poarization ange of the arge-scae GF, and the poarized dust power spectra in the harmonic domain Poarization fraction The naive estimator of p, defined as p = Q 2 + U 2 /I, is a biased quantity. The bias in p becomes dominant in the ow signa-to-noise regime (Simmons & Stewart 1985). ontier et a. (2015) review the different agorithms to debias the p estimate. To avoid the debiasing probem, we choose to work with the unbiased quantity p 2 d derived from a combination of independent subsets of the Panck data at 353 GHz (see Equation (12) of Panck Coaboration Int. p 2 p 2 d p 2 m XLIV 2016): Q s 1 p 2 d,353 Q s2 d,353 d = + U s1 d,353 U s2 D d,353, (4) where the subscript d refers to the observed poarized dust emission from the Panck 353 GHz data (as distinct from m for the mode beow). Normaized histograms of p 2 d are computed over SGC34 from the cross-products of the three subsets of the Panck data where (s 1, s 2 ) = {(H1,H2), (Y1,Y2), (DS1,DS2)}. Figure 2 presents the mean of these three normaized histograms (bue inverted-trianges). The two main sources of uncertainty that bias the normaized histogram of p 2 d are data systematics and the uncertainty of the zero-eve of the Gaactic dust emission. The former is computed per histogram bin using the standard deviation of p 2 d from the three subsets of the Panck data. The atter is propagated by changing the dust intensity D 353 by ± 3 kjy sr 1. The tota 1σ error bar per histogram bin is the quadrature sum of these two uncertainties. The mean, median, and maximum vaues of the dust poarization fraction for the three subsets of the Panck data are presented in Tabe 1. The maximum vaue of p d at 1 resoution is consistent with the ones reported in Panck Coaboration Int. XLIV (2016) over the high atitude southern Gaactic cap and the Panck Coaboration Int. XIX (2015) vaue at ow and intermediate Gaactic atitudes. We define a sky fraction f(p d > 25 %) where the observed dust poarization fraction is arger than a given vaue of 25 % or p 2 d > In SGC34 this fraction is about 1 %. The sma scatter between the resuts from the three different subsets shows that the Panck measurements are robust against any systematic effects. The negative vaues of p 2 d computed from Eq. (4) arise from instrumenta noise present in the three subsets of the Panck data. The width of the p 2 d normaized histogram is narrower than the width reported by Panck Coaboration Int. XLIV (2016) over the high atitude southern Gaactic cap. We were abe to reproduce the broader distribution of p 2 d by adopting the same anaysis region as in Panck Coaboration Int. XLIV (2016), showing the sensitivity to the choice of sky region Poarization ange We use mutipe subsets of the Panck poarization data at 353 GHz to compute the normaized Stokes parameters q s d,353 = Qs d,353 /D 353 and u s d,353 = U s d,353 /D 353, where s = {DS1, DS2, H1, H2, Y1, Y2}. Using these data over SGC34, we fit the mean direction of the arge-scae GF, as defined by two coordinates 0 and b 0, using mode A of Panck Coaboration Int. XLIV (2016). The best-fit vaues of 0 and b 0 for each of the six independent subsets of the Panck data are quoted in Tabe 1. The mean vaues of 0 and b 0 and their dispersions from the independent subsets are ± 2. 0 and ± 1. 1, respectivey. Our error bars are smaer than those (5 ) reported by Panck Coaboration Int. XLIV (2016) because the atter incude the difference between fit vaues obtained with and without appying the bandpass mismatch (BP) correction described in Sect. A.3 of Panck Coaboration VIII (2016). The mean vaues of 0 and b 0 reported by Panck Coaboration Int. XLIV (2016) were derived using maps with no correction Artice number, page 6 of 19

7 T. Ghosh et a.: odeing and simuation of poarized dust emission Tabe 1. Dust poarization fraction p statistics at 353 GHz and the best-fit mean direction of the arge-scae GF over SGC34 using mutipe subsets of the Panck data. Dataset poarization fraction (p d ) f(p d > 25 %) subsets 0 b 0 median mean maximum Detsets 8.0 % 11.0% 33.4 % 1.2 % DS DS Hafissions 8.5 % 11.6 % 35.4 % 1.9 % H H Years 8.3 % 11.3 % 34.6 % 1.4 % YR YR average 8.3 % 11.3 % ± ± 1. 1 Pixe density ψ R [deg] Fig. 3. Simiar to Fig. 2, but for normaized histograms of ψ R, the poarization ange reative to the oca poarization ange of the arge-scae GF. The vertica dashed ine corresponds to a mode in which the turbuence is absent (B turb = 0 so that ψ R = 0 ). appied. The difference of 3. 5 reative to our mean vaue of 0 arises primariy from the BP correction used in our anaysis. The best-fit vaues of 0 and b 0 are used to derive a map of a mean-fied poarization ange ψ 0 (ˆn) (see Panck Coaboration Int. XLIV 2016 for the Q and U patterns associated with the mean fied). We then rotate each sky pixe with respect to the new reference ange ψ 0 (ˆn), using the reation [ ] Q R d,353 (ˆn) = Ud,353 R ( (ˆn) cos 2ψ0 (ˆn) sin 2ψ 0 (ˆn) sin 2ψ 0 (ˆn) cos 2ψ 0 (ˆn) ) [ ] Q s d,353 (ˆn) Ud,353 s (ˆn) ψ R d ψ R m, (5) where Q R d,353 and U d,353 R are the new rotated Stokes parameters. The butterfy patterns caused by the arge-scae GF towards the southern Gaactic cap (Panck Coaboration Int. XLIV 2016) are now removed from the dust poarization maps. Non-zero vaues of the poarization ange, ψd R, derived from the QR d,353 and U d,353 R maps via ψ R d (ˆn) = 1 2 tan 1 ( U R d,353(ˆn), Q R d,353(ˆn)), (6) Tabe 2. Fitted 353 GHz dust power spectra of the Panck data and of the dust mode over SGC34, for which the mean dust intensity D 353 is 54 ± 3 kjy sr 1. Parameter Panck 353 GHz data Dust mode α EE 2.37 ± ± 0.26 α BB 2.33 ± ± 0.38 α T E 2.07 ± ± 0.25 Assuming a common power aw exponent α XX = 2.3 χ 2 EE(N d.o.f. = 6) χ 2 BB(N d.o.f. = 6) χ 2 T E(N d.o.f. = 6) A EE [µk 2 CB] ± ± 1.59 A BB [µk 2 CB ] 8.85 ± ± 1.23 A T E [µk 2 CB] ± ± 2.76 A BB/A EE 0.51 ± ± 0.09 A T E/A EE 2.18 ± ± 0.23 resut from the dispersion of B POS around the mean direction of the arge-scae GF. The minus sign in Ud,353 R is necessary to produce anges in the IAU convention from HEALPix format Stokes maps in the COSO" convention as used for the Panck data. The mean normaized histogram of ψd R is computed from the mutipe subsets of the Panck data and is presented in Fig. 3 (bue inverted-trianges). The 1σ error bars on ψd R are derived from the standard deviation of the six independent measurements of the Panck data, as isted in Tabe 1. We fit a Gaussian to the normaized histogram of ψd R and find a 1σ dispersion of ± 0. 4 over SGC34. Our 1σ dispersion vaue of ψd R is sighty higher than the one found over the high atitude southern Gaactic cap (Panck Coaboration Int. XLIV 2016). This resut is not unexpected because SGC34 extends up to b 30 and adopting a mean direction of the arge-scae GF is ony an approximation Poarized dust power spectra We extend the anaysis by Panck Coaboration Int. XXX (2016) of the region LR24 (fsky eff = 0.24 or 24 % of the tota sky) to the smaer ower coumn density seected area SGC34 within it (fsky eff = 0.085). We compute the poarized dust power spectra D EE, D BB, and D T E by cross- Artice number, page 7 of 19

8 A&A proofs: manuscript no. GASS_ms2 D EE [µk 2 ] Panck dust D EE,353 GHz power-aw fit dust D EE,353 GHz power-aw fit dust D EE,150 GHz best-fit ΛCD CB D EE utipoe EE [µk 2 ] ode dust EE,353 GHz power-aw fit dust EE,353 GHz power-aw fit dust EE,150 GHz best-fit ΛCD CB D EE utipoe D BB [µk 2 ] Panck dust D BB,353 GHz power-aw fit dust D BB,353 GHz power-aw fit dust D,150 BB GHz CB D BB (r = 0.15) CB D BB (ensing) CB D BB (r = 0.15 pus ensing) utipoe BB [µk 2 ] ode dust BB,353 GHz power-aw fit dust BB,353 GHz power-aw fit dust BB,150 GHz CB D BB (r = 0.15) CB D BB (ensing) CB D BB (r = 0.15 pus ensing) utipoe D TE [µk 2 ] Panck dust D TE,353 GHz power-aw fit dust D TE,353 GHz power-aw fit dust D TE,150 GHz best-fit ΛCD CB D TE utipoe TE [µk 2 ] ode dust TE,353 GHz power-aw fit dust TE,353 GHz power-aw fit dust TE,150 GHz best-fit ΛCD CB D TE utipoe Fig. 4. Left coumn: Panck 353 GHz D EE (top row), D BB (midde row), and D T E (bottom row) power spectra in units of µk 2 CB computed over SGC34 using the subsets of the Panck data. See text re cacuation of 1σ error bars. The best-fit power aw assuming exponent 2.3 is potted for each power spectrum (dashed bue ine). For comparison, we aso show the Panck 2015 best-fit ΛCD expectation curves for the CB signa (Panck Coaboration XIII 2016), with the T E expectation shown as dashed where it is negative. Extrapoations of the dust power spectra to 150 GHz (see text) are potted as dashed red ines. Right coumn: Simiar to eft panes, but for the dust mode (Sect. 5.4). Artice number, page 8 of 19

9 T. Ghosh et a.: odeing and simuation of poarized dust emission correating the subsets of the Panck data at 353 GHz in the mutipoe range 40 < < 160 over SGC34. We use the same D 353 map for a subsets of the Panck data. The binning of the measured power spectra in six mutipoe bins matches that used in Panck Coaboration Int. XXX (2016). We take the mean power spectra computed from the cross-spectra of DetSets (DS1 DS2), Hafissions (H1 H2), and Years (Y1 Y2). The mean dust D EE is corrected for the CB contribution using the Panck bestfit ΛCD mode at the power spectrum eve, whie the mean D BB and D T E are kept unatered. Figure 4 shows that the contribution of the CB BB spectra (for r = 0.15 with ensing contribution) is negigiby sma as compared to the mean dust spectra D BB over SGC34. The mean D T E is unbiased because we use D 353 as a dust-ony Stokes I map. However, the unsubtracted CB E-modes bias the 1σ error of the dust D T E power spectrum by a factor [D T T (dust) D EE (CB)]/ν, where ν is associated with the effective number of degrees of freedom (Tristram et a. 2005). The mean dust power spectra D EE, D BB, and D T E of the Panck 353 GHz data are presented in the eft coumn of Fig. 4. The 1σ statistica uncertainties from noise (soid error bar in the eft coumn of Fig. 4) are computed from the anaytica approximation impemented in Xpo. The 1σ systematic uncertainties are computed using the standard deviation of the dust poarization power spectra from the three subsets of the Panck data (as isted in Tabe 1). The tota 1σ error bar per bin (dashed error bar in the eft coumn of Fig. 4) is the quadrature sum of the two. The measured Panck 353 GHz poarized dust power spectra D EE, D BB, and D T E are we described by power aws in mutipoe space, consistent with the one measured at intermediate and high Gaactic atitudes (Panck Coaboration Int. XXX 2016). Simiar to the anaysis in Panck Coaboration Int. XXX (2016), we fit the power spectra with the power-aw mode, D XX = A XX (/80) α XX +2, where A XX is the best-fit ampitude at = 80, α XX is the best-fit sope, and XX = {EE, BB, T E}. We restrict the fit to six band-powers in the range 40 < < 160 so that we can compare the Panck data directy with the ow-resoution dust mode derived using the GASS Hi data. The exponents α XX of the unconstrained fits are quoted in the top row of Tabe 2. Next, we perform the fit with a fixed exponent 2.3 for the appropriate α XX. The resuting power aws are shown in the eft coumn of Fig. 4 and the vaues of χ 2 (with number of degrees of freedom, N d.o.f = 6) and the best-fit vaues of A XX are presented in Tabe 2. The ratio of the dust B- and E-mode power ampitudes, roughy equa to a haf (Tabe 2), is maintained from LR24 to SGC34. The measured ampitude A BB = 8.85 µk 2 CB at = 80 is about 35% ower than the vaue expected from the empirica power-aw A BB D (Panck Coaboration Int. XXX 2016) appied to SGC34. For future work, it wi be important to understand the ampitude and variations of the dust B-mode signa in these ow coumn density regions. We find a significant positive T E correation, D T E, over SGC34; the ampitude A T E normaized at = 80 has more than 4 σ significance when the mutipoes between 40 and 160 are combined in this way. Simiary, a positive T E correation was reported over LR24 (see Figure B.1 of Panck Coaboration Int. XXX 2016). The non-zero positive T E correation is a direct consequence of the correation between the dust intensity structures and the orientation of the magnetic fied. Panck Coaboration Int. XXXVIII (2016) report that oriented stacking of the dust E-map over T peaks of interstear fiaments identified in the dust intensity map picks up a positive T E correation (see their Figure 10). For comparison in Fig. 4 we show the Panck 2015 bestfit ΛCD CB D EE, D BB (primordia tensor-to-scaar ratio r = 0.15 pus ensing contribution), and D T E expectation curves (Panck Coaboration XIII 2016). Then assuming that the dust spectra energy distribution (SED) foows a BB spectrum with a mean dust spectra index 1.59 ± 0.02 and a mean dust temperature of 19.6 K (Panck Coaboration Int. XXII 2015), we extrapoate the best-fit dust power spectra from 353 GHz to 150 GHz (scaing factor s 150/353 = ) in order to show the eve of the dust poarization signa (red dashed ines in Fig. 4) reative to the CB signa over SGC ode framework In this section we extend the formaism deveoped in Panck Coaboration Int. XLIV (2016) to simuate the dust poarization over SGC34 using a phenomenoogica description of the magnetic fied and the GASS Hi emission data. We introduce two main ingredients in the dust mode: (1) fuctuations in the GF orientation aong the LOS and (2) aignment of the CN structures with B POS. Based on the fiament study by Panck Coaboration Int. XXXVIII (2016), we anticipate that aignment of the CN structures with B POS wi account for the positive dust T E correation and the D T E 4.1. Formaism /D EE ratio over SGC34. We use the genera framework introduced by Panck Coaboration Int. XX (2015) to mode the dust Stokes parameters, I, Q, and U, for opticay thin regions at a frequency ν: ( I m,ν (ˆn) = [1 p 0 cos 2 γ i (ˆn) 2 )] S i 3 ν(ˆn) dτ ν Q m,ν (ˆn) = p 0 cos 2 γ i (ˆn) cos 2ψ i (ˆn)S ν(ˆn) i dτ ν (7) U m,ν (ˆn) = p 0 cos 2 γ i (ˆn) sin 2ψ i (ˆn)S ν(ˆn) i dτ ν, where (see Appendix B of Panck Coaboration Int. XX 2015 for detais) ˆn is the direction vector, S ν is the source function given by the reation S ν = n d B ν (T )C avg, τ ν is the optica depth, p 0 = p dust R is the product of the intrinsic dust poarization fraction" (p dust ) and the Rayeigh reduction factor (R, reated to the aignment of dust grains with respect to the GF), ψ is the poarization ange measured from Gaactic North, and γ is the ange between the oca magnetic fied and the POS. Again, ike in Sect. 3.2, the minus sign in U m,ν is necessary to work from the ange map ψ given in the IAU convention to produce HEALPix format Stokes maps in the COSO convention. We start with the same dust Stokes parameters equations as in the Vansynge et a. (2016) parametric dust mode. Artice number, page 9 of 19

10 A&A proofs: manuscript no. GASS_ms2 Foowing the approach in Panck Coaboration Int. XLIV (2016), we repace the integration aong the LOS over SGC34 by the sum over a finite number of ayers with different poarization properties. These ayers are a phenomenoogica means of accounting for the effects of fuctuations in the GF orientation aong the LOS, the first main ingredient of the dust mode. We interpret these ayers aong the LOS as three distinct phases of the IS. The Gaussian decomposition of tota N Hi into the CN, UN, and WN components based on their ine-widths (or veocity dispersion) is described in detai in Sect We can repace the integra of the source function over each ayer directy with the product of the coumn density maps of the three Hi components and their mean dust emissivities. Within these astrophysica approximations, we can modify Eq. (7) to be I m,ν (ˆn) = Q m,ν (ˆn) = U m,ν (ˆn) = 3 i=1 [1 p 0 ( cos 2 γ i (ˆn) 2 3 )] ɛ i ν N i Hi(ˆn) 3 p 0 cos 2 γ i (ˆn) cos 2ψ i (ˆn) ɛ i ν NHi(ˆn) i (8) i=1 3 p 0 cos 2 γ i (ˆn) sin 2ψ i (ˆn) ɛ i ν NHi(ˆn) i, i=1 where ɛ i ν is the mean dust emissivity for each of the three Hi components. In genera, the emissivity is a function not ony of frequency, but aso of sky position ˆn and different Hi phases. However, in SGC34 the overa emissivity D 353 /N Hi is quite uniform (see Sect. 5.2) and from a simutaneous fit of the D 353 map to a three component N Hi maps we find itte evidence, for our decomposition, that the emissivities of the individua components are different than the overa mean. Simiary, the effective degree of aignment p 0 = p dust R coud be a function of frequency and different for different Hi components. Equation (8) aso tes us that the product p 0 is degenerate with ɛ ν for the dust Stokes Q m,ν and U m,ν parameters. Such effects coud make the source functions in the ayers even more different from one another. In this paper we expore how modes based on simpifying assumptions, namey a constant p 0 and a frequency dependent ɛ ν common to the three Hi components, might despite the reduced fexibiity nevertheess provide a good description of the data. This use of the GASS N Hi data can be contrasted to the approach of Panck Coaboration Int. XLIV (2016) and Vansynge et a. (2016), where the source functions are assumed to be the same in each ayer. Furthermore, Panck Coaboration Int. XLIV (2016) choose to ignore the correction term invoving p 0 in the integrand for I m,ν (ˆn). This phenomenoogica dust mode based on the GASS Hi data is not unique. It is based on severa astrophysica assumptions as noted in this section. Within the same genera framework, we coud contempate repacing the Hi emission with three-dimensiona (3D) extinction maps, though no reiabe maps are currenty avaiabe over the SGC34 region (Green et a. 2015). Furthermore, the Hi emission traces the temperature or density structure of the diffuse IS, whereas the dust extinction traces ony the dust coumn density. Thus, using the Hi emission to mode the dust poarization is an interesting approach even if we had reiabe 3D extinction maps toward the southern Gaactic cap Structure of the Gaactic magnetic fied The Gaactic magnetic fied, B, is expressed as a vector sum of a mean arge-scae (ordered, B ord ) and a turbuent (random, B turb ) component (Jaffe et a. 2010) B(ˆn) = B ord (ˆn) + B turb (ˆn) = B ord ( ˆB ord (ˆn) + f ˆBturb (ˆn)), (9) where f is the standard deviation of the reative ampitude of the turbuent component B turb with respect to the mean arge-scae B ord. The butterfy patterns seen in the orthographic projections of the Panck poarization maps, Q d,353 and U d,353, centred at (, b) = (0, 90 ) are we fitted by a mean direction of the arge-scae GF over the southern Gaactic cap (Panck Coaboration Int. XLIV 2016). We assume that the mean direction of the arge-scae GF is sti a good approximation to fit the butterfy patterns seen in the Panck Stokes Q d,353 and U d,353 maps over SGC34 (Fig. 1). The unit vector ˆB ord is defined as ˆB ord = (cos 0 cos b 0, sin 0 cos b 0, sin b 0 ). We mode B turb with a Gaussian reaization on the sky with an underying power spectrum, C α, for 2 (Panck Coaboration Int. XLIV 2016). The correated patterns of B turb over the sky are needed to account for the correated patterns of dust p and ψ seen in Panck 353 GHz data (Panck Coaboration Int. XIX 2015). Foowing Fauvet et a. (2012), we ony consider isotropic turbuence in our anaysis. We produce a spatia distribution of B turb in the x, y, and z direction from the above power spectrum on the ceestia sphere of HEALPix resoution N side. We add B turb to B ord with a reative strength f as given in Eq Aignment of the CN structures with the magnetic fied We use the agebra given in Section 4.1 of Panck Coaboration Int. XLIV (2016) to compute cos 2 γ i for the three Hi components. For the UN and WN components, we assume that there is no T E correation and the anges ψ u and ψ w are computed directy from the tota B again using the agebra by Panck Coaboration Int. XLIV (2016). For the CN we foow a quite different approach to compute ψ c. We introduce a positive T E correation for the CN component through the poarization ange ψ c. The aignment between the CN structures and B POS is introduced to test whether it fuy accounts for the observed E-B power asymmetry or not. To do that, we assume that the masked N c Hi map is a pure E-mode map and define spin-2 maps Q T and U T as (Q T ± iu T )(ˆn) = =2 m= a c m ±2Y m (ˆn), (10) where a c m are the harmonic coefficients of the masked N c Hi map. In the fat-sky imit, Q T ( 2 x 2 y) 2 N c Hi and U T 2 x y 2 N c Hi, where 2 is the inverse Lapacian operator (Bowyer et a. 2011). We compute the ange ψ c using the reation, ψ c = (1/2) tan 1 ( U T, Q T ) (see Panck Coaboration XVI 2016 for detais). The poarization ange ψ c is the same as the ange between the major axis defined by a oca quadrature expansion of the N c Hi map and the horizonta axis. This idea of aigning perfecty the CN Artice number, page 10 of 19

11 T. Ghosh et a.: odeing and simuation of poarized dust emission structures with B POS is motivated from the Hi studies by Cark et a. (2015), artin et a. (2015), and Kabera et a. (2016), and Panck dust poarization studies by Panck Coaboration Int. XXXII (2016) and Panck Coaboration Int. XXXVIII (2016). This is the second main ingredient of the dust mode ode parameters The dust mode has the foowing parameters: ɛ ν, the mean dust emissivity per Hi component at frequency ν; p 0, the normaization factor to match the observed dust poarization fraction over the southern Gaactic cap; 0 and b 0, describing the mean direction of B ord in the region of the southern Gaactic cap; f, the reative strength of B turb with respect to B ord for the different components; α, the spectra index of B turb power spectrum for the different components; and σ c and σ u to separate the tota N Hi into the CN, UN, and WN components (phases). We treat the three Hi components with different density structures and magnetic fied orientations as three independent ayers. Because we incude aignment of the CN structures with B POS, we treat the CN poarization ayer separatey from the WN and UN poarization ayers. The vaue of f c for the CN component is different than the common f u/w for the UN and WN components. For the same reason, the spectra index of the turbuent power spectrum α c for the CN component coud be different than the common α u/w for the UN and WN components, but we take them to be the same, simpifying the mode and reducing the parameter space in our anaysis. We aso use the same vaues of ɛ ν and p 0 for each ayer. This gives a tota of nine parameters. In this anaysis, we do not optimize the nine parameters of the dust mode expicity. Instead, we determine five parameters based on various astrophysica constraints and optimize ony ɛ ν, p 0, f u/w, and σ c simutaneousy (Sect. 5.2). Our main goa is to test quantitativey whether the fuctuations in the GF orientation aong the LOS pus the preferentia aignment of the CN structures with B POS can account simutaneousy for both the observed dust poarization power spectra and the normaized histograms of p 2 and ψ R over SGC34. We review here the reationship between the nine dust mode parameters and the various dust observabes. The emissivity ɛ ν simpy converts the observed N Hi (10 20 cm 2 ) to the Panck intensity units (in µk CB ) (Panck Coaboration Int. XVII 2014). This parameter is optimized in Sect. 5.2, subject to a strong prior obtained by correating the map D 353 with the N Hi map. The parameter p 0 is constrained by the observed distribution of p 2 d at 353 GHz and the ampitude of the dust poarization power spectra D EE, D BB, and D T E. It is optimized using the atter in Sect The two mode parameters 0 and b 0 determine the orientation of the butterfy patterns seen in the orthographic projections of the Panck poarization maps, Q d,353 and U d,353 (Panck Coaboration Int. XLIV 2016). They are found by fitting these data (Sect. 3.2). The parameter f for each poarization ayer is reated to the dispersion of ψd R and p2 d (the higher the strength of the turbuent fied with respect to the ordered fied, the higher the dispersion of ψd R and p2 d ) and to the ampitude of the dust poarization power spectra. Our prescription for the aignment of CN structure and B POS subject to statistica constraints from poarization data determines f c u/w (Sect. 5.1); f is optimized using the dust poarization band powers in Sect The parameter α c is constrained by the aignment of the CN structures with B POS (Sect. 5.1). The parameter σ c contros the fraction of the Hi in the CN and hence sets the amount of E-B power asymmetry and the D T E /D EE ratio through the aignment of the CN structures with B POS. It is optimized in Sect Our dust mode is insensitive to the precise vaue of the parameter σ u, which is used to separate the UN and WN into distinct ayers. To break the degeneracy between different dust mode parameters and obtain a consistent and robust mode, we need to make use of dust poarization observabes in pixe space as we as in harmonic space. 5. Dust sky simuations We simuate the dust intensity and poarization maps at ony a singe Panck frequency 353 GHz. It woud be straightforward to extrapoate the dust mode to other microwave frequencies using a BB spectrum of the dust emission (Panck Coaboration Int. XVII 2014; Panck Coaboration Int. XXII 2015; Panck Coaboration Int. XXX 2016) or repacing the vaue of ɛ ν from Panck Coaboration Int. XVII (2014). However, such a simpe-minded approach woud not generate any decorreation of the dust poarization pattern between different frequencies, as reported in Panck Coaboration Int. L (2016) Constraining the turbuent GF for the CN component As discussed in Sect. 4.3, we adopted an aternative approach to compute the poarization ange ψ c for the CN component. At first sight, it might ook physicay inconsistent and therefore a main source of bias. However, we now show that the origina and aternative approaches ead to statisticay simiar resuts and so can be compared in order to constrain the parameters of B turb, that is α c and f c, for the CN component. To do this we first adopt the origina approach from Panck Coaboration Int. XLIV (2016) (and used for the UN and CN components). We generate Gaussian reaizations of B turb for a set of (α, f ) vaues. We assume a fixed mean direction of B ord as ( 0, b 0 ) = (73. 5, 23. 5) (Tabe 1) and add to it B turb to derive a tota B. We use the agebra given in Section 4.1 of Panck Coaboration Int. XLIV (2016) to compute the two quantities cos 2 γ and cos 2ψ for different vaues of α and f. Because the poarization anges γ and ψ are circuar quantities, we choose to work with quantities ike cos 2 γ and cos 2ψ as they appear in the Stokes Q and U parameters (Eq. 8). We then appy the Xpo routine to compute the anguar power spectra of these two quantities for each sky reaization over SGC34. The mean power spectra of cos 2 γ and cos 2ψ are computed by averaging 100 onte-caro reaizations and are presented in Fig. 5. The power spectra of the two poarization anges are we described by power-aws in mutipoe Artice number, page 11 of 19

12 A&A proofs: manuscript no. GASS_ms2 D(cos 2 γ) D(cos 2ψ) f = 0.3, α = 2.4 f = 0.4, α = 2.4 f = 0.4, α = 2.3 f = 0.4, α = 2.2 utipoe f = 0.3, α = 2.4 f = 0.4, α = 2.4 f = 0.4, α = 2.3 f = 0.4, α = 2.2 CN component utipoe Fig. 5. ean power spectra of cos 2 γ (top pane) and cos 2ψ (bottom pane) over SGC34 for a set of α and f parameters in the mutipoe range 40 < < 160. The 1σ error bars are computed from the standard deviations of the 100 onte-caro reaizations of the dust mode. The aignment of the CN structures with B POS (see Sect. 5.1 for detais) constrains the vaue of α c and f c for the CN component. space, D α+2, where α is the same as the spectra index of B turb power spectrum. The ampitudes and sopes of these two quantities are reated systematicay to the vaues of α and f. The higher the f vaue, the higher is the ampitude of the power spectra in Fig. 5. The bestfit sope of the power spectrum over the mutipoe range 40 < < 160 is the same as the spectra index of the B turb power spectrum for reevant f vaues. For each onte Caro reaization we aso cacuated ψ c for the CN component by the aternative approach and from that computed the power spectra of cos 2ψ for the CN component over SGC34. By comparing these resuts to the spectra from the first Panck Coaboration Int. XLIV (2016) approach, as in Fig. 5 (bottom pane), we can constrain the vaues of α c and f c. The typica vaues of αc and f c are 2.4 and 0.4, respectivey. One needs a onte- Caro approach to put reaistic error bars on the α c and f c parameters, which is beyond the scope of this paper. However, we do note that the vaue of α c is cose to the exponent of the anguar power spectrum of the NHi c map computed over SGC34, refecting the correation between the structure of the intensity map and the orientation of the fied. We stress that the vaues of α c and f c derived from Fig. 5 are vaid ony over SGC34. If we increase the sky fraction to b 30 covered by the GASS Hi survey, we get the same α c vaue, but a sighty higher f c vaue (f c = 0.5). This suggests that the reative strength of B turb with respect to B ord might vary across the sky Vaues of mode parameters Here we describe how the nine mode parameters were found. We found the mean direction of B ord to be 0 = and b 0 = from a simpe mode fit of normaized Stokes parameters for poarization over SGC34 (Sect. 3.2, Tabe 1). As described in Sect. 5.1, using NHi c as a tracer of the dust poarization ange gives the constraints α c = 2.4 and f c = 0.4. The spectra index αu/w of B turb for the UN and WN components is assumed to be same as the α c vaue of the CN component, for simpicity. Ony four mode parameters, ɛ 353, p 0, f u/w, and σ c, are adjusted by evauating the goodness of fit through a χ 2 -test, specificay χ 2 XX = and χ 2 = N pix i max min [ D XX [ D 353 XX ] ( ɛ 353, p 0, f u/w 2, σ c) (11) σ XX si m,353 ( ɛ 353, p 0, f u/w, σ c) o] 2, (12) where XX = {EE, BB, T E}, D XX, and σ XX are the mean observed dust poarization band powers and standard deviation at 353 GHz, respectivey, XX is the mode dust poarization power spectrum over SGC34, N pix is the tota number of pixes in SGC34, s and o are the sope and the offset of the T-T correation between the D 353 and I m,353 maps. We evauated χ 2 XX using the six band powers in the range 40 < < 160. For χ 2 minimization, we put a fat prior on the mean dust emissivity, ɛ 353 = (111 ± 22) µk CB (10 20 cm 2 ) 1 at 353 GHz, as obtained by correating N Hi with D 353 map over SGC34. Based on the observed distribution of the poarization fraction (Sect. 3.1), we put a fat prior on p 0 = (18 ± 4) %. To match the observed dust ampitude at 353 GHz, the vaue of s shoud be cose to 1 and is kept so by adjusting ɛ 353 during the iterative soution of the two equations. The typica mean vaues of ɛ 353, p 0, f u/w, and σ c that we obtained are µk CB (10 20 cm 2 ) 1, 18.5 %, 0.1, and 7.5 km s 1, respectivey, with χ 2 EE = 5.0, χ2 BB = 3.3, and χ 2 T E = 4.9 for six degrees of freedom (or band-powers). The precise optimization is not important and is beyond the scope of this paper. Artice number, page 12 of 19

13 T. Ghosh et a.: odeing and simuation of poarized dust emission D353 [µk] best-fit ine I m,353 [µk] Fig. 6. Correation of the Panck 353 GHz dust intensity map and the dust mode Stokes I map. The dashed ine is the best-fit reation. The mode is not particuary sensitive to the vaue of the fina parameter, σ u, but it is important to have significant coumn density in each of UN and WN so as to have contributions from three ayers; we chose σ u = 10 km s 1, for which the average coumn densities of the phases over SGC34 are rather simiar (Sect. 2.3). Our vaue of ɛ 353 differs sighty from the mean vaue ɛ 353 quoted in Tabe 2 of Panck Coaboration Int. XVII (2014), where the Panck intensity maps are correated with Hi LVC emission over many patches within Gaactic atitude b < 25. A difference in the mean dust emissivity is not unexpected because our study is focussed on the ow coumn density region (N Hi cm 2 ) within the Panck Coaboration Int. XVII (2014) sky region (N Hi cm 2 ). Fig. 6 shows the correation between the data D 353 and the mode I m,353 over SGC34. The parameters of the best fit ine are s = 1.0 (iterativey by construction) and o = 6.6 µk CB. We compute a residua map D 353 I m,353 over SGC34. The dispersion of the residua emission with respect to the tota D 353 emission is σ /σ D353 = This means that the residua emission accounts for ony 13 % (σ 2 /σ2 D 353 ) of the tota dust power at 353 GHz. We cacuate the Pearson correation coefficients over SGC34 between the residua map and the three Hi tempates, namey the CN, UN, and WN, finding 0.09, 0.17, and 0.13, respectivey. The residua might originate in pixedependent variations of ɛ 353 or possiby from diffuse ionized gas that is not spatiay correated with the tota N Hi map. Panck Coaboration Int. XVII (2014) expored these two possibiities (see their Appendix D) and were abe to reproduce the observed vaue of σ / D 353 at 353 GHz. We do not use the residua emission in our anaysis because our main goa is to reproduce the mean dust poarization properties over the SGC34 region. For future work, the pixe-dependent variation of ɛ 353, ike in the Panck Coaboration Int. XVII (2014) anaysis, coud be incuded so that one coud compare the refined dust mode with the observed dust poarization data for a specific sky patch (e.g. the BICEP2 fied). TT [µk 2 ] ode dust TT,353 GHz power-aw fit dust TT,353 GHz power-aw fit dust TT,150 GHz best-fit ΛCD CB D TT utipoe Fig. 7. Best-fit power aw dust mode T T power spectra at 353 GHz (back dashed ine) and extrapoated dust power spectra at 150 GHz (red dashed ine) as computed over SGC34. For comparison, we aso show the Panck 2015 best-fit ΛCD expectation curves of the CB signa onte-caro simuations Using the genera framework described in Sect. 4, we simuate a set of 100 poarized dust sky reaizations at 353 GHz using a set of mode parameters, as discussed in Sect Because the mode of the turbuent GF is statistica, we can simuate mutipe dust sky reaizations for a given set of parameters. The simuated dust intensity and poarization maps are smoothed to an anguar resoution of FWH 1 and projected on a HEALPix grid of N side =128. These dust simuations are vaid ony in the ow coumn density regions. To compare the dust mode with the Panck observations, we add reaistic noise simuations caed fu foca pane 8 or FFP8" (Panck Coaboration XII 2016). These noise simuations incude the reaistic noise correation between pair of detectors for a given boometer. For a given dust reaization, we produce two independent sampes of dust Stokes Q and U noisy mode maps as Q si m,353 (ˆn) = Q m,353(ˆn) + Q i n(ˆn), U si m,353 (ˆn) = U m,353(ˆn) + Q i n(ˆn), (13) where i = 1, 2 are two independent sampes, and Q i n and Un i are the statistica noise of the subsets of the Panck poarization maps (as isted in Tabe 1) Comparison of the dust mode with Panck observations In this section we compare the statistica properties of the dust poarization from the dust mode with the Panck 353 GHz observations over the seected region SGC34. We showed in Fig. 6 that the mode Stokes I353 m provides a good fit to the D 353 map in ampitude. Here we show that this is the case for the shape of the D T T power spectrum too. We compute the mode dust T T power spectrum, T T ( + 1)C,m T T /(2π), over SGC34 in the range 40 < < 600. The dust T T power spectrum is we represented with a power-aw mode, T T (/80) α T T +2, with Artice number, page 13 of 19

14 A&A proofs: manuscript no. GASS_ms2 a best-fit sope α T T = 2.59±0.02, as shown in Fig. 7. Our derived best-fit vaue of α T T is consistent with that measured in Panck Coaboration XVI (2014) and Panck Coaboration XIII (2016) at high Gaactic atitude. For comparison, we aso show the T T power spectrum of the CB signa for the Panck 2015 best-fit ΛCD mode (Panck Coaboration XIII 2016). The best-fit dust T T power spectrum extrapoated from 353 GHz to 150 GHz, scaing as in Fig. 4 and shown here as a dashed red ine, reveas that over SGC34 the ampitude of the dust emission at 150 GHz is ess than 1 % compared to the ampitude of the CB signa. Next we show that the dust mode is abe to reproduce the observed distribution on p 2 d and dispersion of ψr d over SGC34 shown in Figs. 2 and 3. We note that these statistics were not used to optimize the four mode parameters in Sect. 5.2, and so this serves as an important quaity check for the mode. We compute the square of the poarization fraction, p 2 m, for the dust mode Q s 1m,353 p 2 Q s2 m,353 m = + U s1 m,353 U s2 m,353 Im,353 2, (14) where the index m again refers to the dust mode. We use the two independent sampes to expoit the statistica independence of the noise between them. The mean normaized histogram of p 2 m and the associated 1σ error bars are computed over SGC34 using 100 onte-caro simuations. Returning to Fig. 2, we compare the normaized histogram of p 2 m (back circes) with the normaized histogram of p 2 d (bue inverted-trianges). The Panck instrumenta noise (FFP8 noise) added in the dust mode accounts nicey for the observed negative vaues of p 2 m and aso contributes the extension of the p 2 m distribution beyond the input vaue of p 2 0. The vaue of p 0 found in this paper is cose to the vaue of 19% deduced at 1 resoution over intermediate and ow Gaactic atitudes (Panck Coaboration Int. XIX 2015). We aso compute the dispersion of the poarization ange ψ R m for the dust mode. We foow the same procedure as discussed in Sect. 3.3 and compute the ange ψ R m = 1 2 tan 1 ( U R m,353, Q R m,353), (15) where Q R m,353 and Um,353 R are rotated Stokes parameters with respect to the oca direction of the arge-scae GF. We make the normaized histogram of ψm R over SGC34. The mean normaized histogram and associated 1σ error bars of ψm R for the dust mode are computed from 100 onte-caro simuations. Returning to Fig. 3, we compare the dispersion of ψm R (back points) with the dispersion of ψd R (bue points). The 1σ dispersion of ψm R derived from the mean of 100 dust mode reaizations is ± 0. 7, sighty arger than for the Panck data (15. 0 ± 0. 4, Sect. 3.2). This coud come from the simpe modeing assumptions used in our anaysis. In particuar, the histograms of p 2 d and ψr d depend on the dust modeing at ow mutipoes ( < 40), which are not constrained by our data fitting. Varying some of the mode parameters such as α u/w, different R for different IS phases, or introducing a ow- cutoff in the spectra index of the B turb (Cho & Lazarian 2002, 2010) might provide a better fit of these histograms to the Panck data. Due to the imited sky coverage over SGC34 and residua systematics in the pubicy avaiabe Panck PR2 data at ow mutipoes (Panck Coaboration Int. XLVI 2016), we are unabe to test some of these possibiities within our modeing framework. Finay, we compute power spectra of the dust mode over SGC34 using the cross-spectra of two independent sampes as described in Sects. 3.3 and 5.3. ean dust poarization power spectra of the noisy dust mode, EE, BB, and T E, over the mutipoe range 40 < < 160 are cacuated from 100 onte-caro simuations and presented in the right coumn of Fig. 4. The associated 1σ error bar per bin is aso derived from the onte-caro simuations. A of the dust poarization power spectra are we represented by a power aw in. We fit the mode spectra with a power-aw mode over the mutipoe range 40 < < 160. The best-fit vaues of the exponents and then of the ampitudes with exponent 2.3 are quoted in Tabe 2. Because of the χ 2 optimization of the four mode parameters (Sect. 5.2), these power spectra agree we with the observed power spectra over SGC Predictions from the dust mode The simuated dust emission maps with FFP8 noise approximate the Panck 353 GHz data we over SGC34. In this section we use noiseess dust simuations to make a few predictions concerning the foreground dust B-mode that are important for the search for the CB B-mode signa. These predictions incude the T B and EB correation for the dust emission and the statistica variance of the dust B-modes at the power spectrum eve. But first we compare the spatia distribution of Stokes I and the poarization intensity (P 2 ) of the noiseess dust mode with the Panck 353 GHz data over the southern Gaactic cap. The comparison between one reaization of the dust mode and the Panck data is presented in Fig. 8. The dust simuation is derived competey from the GASS Hi data and a phenomenoogica description of the argescae and turbuent components of the GF with parameters optimized over SGC34, a subregion of the SGC. We find a good match between the data and the dust mode over SGC34, outined by the back contour. Outside of SGC34, the brightest emission features appear stronger in the Panck 353 GHz data because the Hi-based dust mode does not account for dust emission associated with H 2 gas and with Hi gas that is too cod to resut in significant net emission Dust T B and EB correation The T B and EB cross-spectra vanish for the CB signa in the standard ΛCD mode (Zadarriaga & Sejak 1997). Here we test whether the corresponding dust mode T B and EB cross-spectra vanish or not in the ow coumn density region SGC34. To do this we simuate 100 onte- Caro reaizations of noiseess dust skies and compute the cross-spectra for each reaization. The mean mode dust EB spectrum over the mutipoe range 40 < < 160 is presented in Fig. 9. Simiar to Abitbo et a. (2016), we fit the EB power spectrum with a power-aw mode, EB = A EB (/80) αeb+2 with a fixed sope α EB = 2.3 and extrapoate it to 150 GHz using the scaing factor s 150/353. Our dust mode predicts A EB = 0.31 ± 0.13 µk 2 CB (2.4 σ eve) at 353 GHz. The Artice number, page 14 of 19

15 T. Ghosh et a.: odeing and simuation of poarized dust emission 0 o 0 o 45 o 315 o 45 o 315 o 90 o -60 o 270 o 90 o -60 o 270 o -30 o -30 o 135 o 225 o 135 o 225 o 180 o 180 o µkcb 0 o 0 o 45 o 315 o 45 o 315 o 90 o -60 o 270 o 90 o -60 o 270 o -30 o -30 o 135 o 225 o 135 o 225 o 180 o 180 o µk 2 CB Fig. 8. Left coumn: Orthographic projections of the Panck 353 GHz D 353 (Stokes I) map (top row) and the square of the poarization intensity P 2 (bottom row), at 1 resoution over the southern Gaactic cap covered by the GASS Hi survey, with the seected region SGC34 outined by a back contour. The same coordinate system is used as in Fig. 1. Right coumn: Simiar to eft panes, but for a reaization of the 353 GHz dust mode. Judged visuay, the dust mode is abe to reproduce roughy the observed dust poarization sky over SGC34, and indeed over the entire fied shown. best-fit vaue of A EB at = 80 transates into an EB ampitude of µk 2 CB at 150 GHz, consistent with the resuts presented in Tabe A1 of Abitbo et a. (2016) over the BICEP2 patch (1.3 % of the ceanest sky region). Because our dust simuations are noiseess, we can determine the ampitude of the EB spectrum with reativey sma error bars. This non-zero dust EB spectrum towards the SGC34 region can produce a spurious B-mode signa, if not taken into account in the sef-caibrated teescope poarization ange (Abitbo et a. 2016). Comparing our resut with Abitbo et a. (2016) indicates that the ampitude of spuri- Artice number, page 15 of 19