A novel estimator of the polarization amplitude from normally distributed Stokes parameters

We propose a novel estimator of the polarization amplitude from a single measurement of its normally distributed $(Q,U)$ Stokes components. Based on the properties of the Rice distribution and dubbed 'MAS' (Modified ASymptotic), it meets several desirable criteria:(i) its values lie in the whole positive region; (ii) its distribution is continuous; (iii) it transforms smoothly with the signal-to-noise ratio (SNR) from a Rayleigh-like shape to a Gaussian one; (iv) it is unbiased and reaches its components' variance as soon as the SNR exceeds 2; (v) it is analytic and can therefore be used on large data-sets. We also revisit the construction of its associated confidence intervals and show how the Feldman-Cousins prescription efficiently solves the issue of classical intervals lying entirely in the unphysical negative domain. Such intervals can be used to identify statistically significant polarized regions and conversely build masks for polarization data. We then consider the case of a general $[Q,U]$ covariance matrix and perform a generalization of the estimator that preserves its asymptotic properties. We show that its bias does not depend on the true polarization angle, and provide an analytic estimate of its variance. The estimator value, together with its variance, provide a powerful point-estimate of the true polarization amplitude that follows an unbiased Gaussian distribution for a SNR as low as 2. These results can be applied to the much more general case of transforming any normally distributed random variable from Cartesian to polar coordinates.Comment: Accepted by MNRA

Simulated CII observations for SPICA/SAFARI

We investigate the case of CII 158 micron observations for SPICA/SAFARI using a three-dimensional magnetohydrodynamical (MHD) simulation of the diffuse interstellar medium (ISM) and the Meudon PDR code. The MHD simulation consists of two converging flows of warm gas (10,000 K) within a cubic box 50 pc in length. The interplay of thermal instability, magnetic field and self-gravity leads to the formation of cold, dense clumps within a warm, turbulent interclump medium. We sample several clumps along a line of sight through the simulated cube and use them as input density profiles in the Meudon PDR code. This allows us to derive intensity predictions for the CII 158 micron line and provide time estimates for the mapping of a given sky area.Comment: 4 pages, 5 figures, to appear in the proceedings of the workshop "The Space Infrared Telescope for Cosmology & Astrophysics: Revealing the Origins of Planets and Galaxies" (July 2009, Oxford, United Kingdom

Polarization measurements analysis II. Best estimators of polarization fraction and angle

With the forthcoming release of high precision polarization measurements, such as from the Planck satellite, it becomes critical to evaluate the performance of estimators for the polarization fraction and angle. These two physical quantities suffer from a well-known bias in the presence of measurement noise, as has been described in part I of this series. In this paper, part II of the series, we explore the extent to which various estimators may correct the bias. Traditional frequentist estimators of the polarization fraction are compared with two recent estimators: one inspired by a Bayesian analysis and a second following an asymptotic method. We investigate the sensitivity of these estimators to the asymmetry of the covariance matrix which may vary over large datasets. We present for the first time a comparison among polarization angle estimators, and evaluate the statistical bias on the angle that appears when the covariance matrix exhibits effective ellipticity. We also address the question of the accuracy of the polarization fraction and angle uncertainty estimators. The methods linked to the credible intervals and to the variance estimates are tested against the robust confidence interval method. From this pool of estimators, we build recipes adapted to different use-cases: build a mask, compute large maps, and deal with low S/N data. More generally, we show that the traditional estimators suffer from discontinuous distributions at low S/N, while the asymptotic and Bayesian methods do not. Attention is given to the shape of the output distribution of the estimators, and is compared with a Gaussian. In this regard, the new asymptotic method presents the best performance, while the Bayesian output distribution is shown to be strongly asymmetric with a sharp cut at low S/N.Finally, we present an optimization of the estimator derived from the Bayesian analysis using adapted priors

First detection of [N II] 205 micrometer absorption in interstellar gas

We present high resolution [NII] 205 micrometer ^3P_1-^3P_0 spectra obtained with Herschel-HIFI towards a small sample of far-infrared bright star forming regions in the Galactic plane: W31C (G10.6-0.4), W49N (G43.2-0.1), W51 (G49.5-0.4), and G34.3+0.1. All sources display an emission line profile associated directly with the HII regions themselves. For the first time we also detect absorption of the [NII] 205 micrometer line by extended low-density foreground material towards W31C and W49N over a wide range of velocities. We attribute this absorption to the warm ionised medium (WIM) and find N(N^+)\approx 1.5x10^17 cm^-2 towards both sources. This is in agreement with recent Herschel-HIFI observations of [CII] 158 micrometer, also observed in absorption in the same sight-lines, if \approx7-10 % of all C^+ ions exist in the WIM on average. Using an abundance ratio of [N]/[H] = 6.76x10^-5 in the gas phase we find that the mean electron and proton volume densities are ~0.1-0.3 cm^-3 assuming a WIM volume filling fraction of 0.1-0.4 with a corresponding line-of-sight filling fraction of 0.46-0.74. A low density and a high WIM filling fraction are also supported by RADEX modelling of the [NII] 205 micrometer absorption and emission together with visible emission lines attributed mainly to the WIM. The detection of the 205 micrometer line in absorption emphasises the importance of a high spectral resolution, and also offers a new tool for investigation of the WIM.Comment: 7 pages, 4 figures, accepted for publication in Astronomy & Astrophysics, 11 June 201

Magnetic field morphology in nearby molecular clouds as revealed by starlight and submillimetre polarization

Within four nearby (d < 160 pc) molecular clouds, we statistically evaluate the structure of the interstellar magnetic field, projected on the plane of the sky and integrated along the line of sight, as inferred from the polarized thermal emission of Galactic dust observed by Planck at 353 GHz and from the optical and NIR polarization of background starlight. We compare the dispersion of the field orientation directly in vicinities with an area equivalent to that subtended by the Planck effective beam at 353 GHz (10') and using the second-order structure functions of the field orientation angles. We find that the average dispersion of the starlight-inferred field orientations within 10'-diameter vicinities is less than 20 deg, and that at these scales the mean field orientation is on average within 5 deg of that inferred from the submillimetre polarization observations in the considered regions. We also find that the dispersion of starlight polarization orientations and the polarization fractions within these vicinities are well reproduced by a Gaussian model of the turbulent structure of the magnetic field, in agreement with the findings reported by the Planck collaboration at scales greater than 10' and for comparable column densities. At scales greater than 10', we find differences of up to 14.7 deg between the second-order structure functions obtained from starlight and submillimetre polarization observations in the same positions in the plane of the sky, but comparison with a Gaussian model of the turbulent structure of the magnetic field indicates that these differences are small and are consistent with the difference in angular resolution between both techniques.Comment: 15 pages, 10 figures, submitted to A&

On the use of fractional Brownian motion simulations to determine the 3D statistical properties of interstellar gas

Based on fractional Brownian motion (fBm) simulations of 3D gas density and velocity fields, we present a study of the statistical properties of spectro-imagery observations (channel maps, integrated emission, and line centroid velocity) in the case of an optically thin medium at various temperatures. The power spectral index gamma_W of the integrated emission is identified with that of the 3D density field (gamma_n) provided the medium's depth is at least of the order of the largest transverse scale in the image, and the power spectrum of the centroid velocity map is found to have the same index gamma_C as that of the velocity field (gamma_v). Further tests with non-fBm density and velocity fields show that this last result holds, and is not modified either by the effects of density-velocity correlations. A comparison is made with the theoretical predictions of Lazarian & Pogosyan (2000).Comment: 28 pages, 14 figures, accepted for publication in ApJ. For preprint with higher-resolution figures, see http://www.cita.utoronto.ca/~mamd/miville_fbm2003.pd

MESMER: MeerKAT Search for Molecules in the Epoch of Reionization

[Abridged] Observations of molecular gas at all redshifts are critical for measuring the cosmic evolution in molecular gas density and understanding the star-formation history of the Universe. The 12CO molecule (J=1-0 transition = 115.27 GHz) is the best proxy for extragalactic H2, which is the gas reservoir from which star formation occurs, and has been detected out to z~6. Typically, redshifted high-J lines are observed at mm-wavelengths, the most commonly targeted systems exhibiting high SFRs (e.g. submm galaxies), and far-IR-bright QSOs. While the most luminous objects are the most readily observed, detections of more typical galaxies with modest SFRs are essential for completing the picture. ALMA will be revolutionary in terms of increasing the detection rate and pushing the sensitivity limit down to include such galaxies, however the limited FoV when observing at such high frequencies makes it difficult to use ALMA for studies of the large-scale structure traced out by molecular gas in galaxies. This article introduces a strategy for a systematic search for molecular gas during the EoR (z~7 and above), capitalizing on the fact that the J=1-0 transition of 12CO enters the upper bands of cm-wave instruments at high-z. The FoV advantage gained by observing at such frequencies, coupled with modern broadband correlators allows significant cosmological volumes to be probed on reasonable timescales. In this article we present an overview of our future observing programme which has been awarded 6,500 hours as one of the Large Survey Projects for MeerKAT, the forthcoming South African SKA pathfinder instrument. Its large FoV and correlator bandwidth, and high-sensitivity provide unprecedented survey speed for such work. An existing astrophysical simulation is coupled with instrumental considerations to demonstrate the feasibility of such observations and predict detection rates.Comment: 7 pages, 4 figures, to appear in the proceedings of "Astronomy with megastructures: Joint science with the E-ELT and SKA", 10-14 May 2010, Crete, Greece (Eds: Isobel Hook, Dimitra Rigopoulou, Steve Rawlings and Aris Karastergiou

Galactic interstellar filaments as probed by LOFAR and Planck

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

UV-driven chemistry in simulations of the interstellar medium. I. Post-processed chemistry with the Meudon PDR code

Our main purpose is to estimate the effect of assuming uniform density on the line-of-sight in PDR chemistry models, compared to a more realistic distribution for which total gas densities may well vary by several orders of magnitude. A secondary goal of this paper is to estimate the amount of molecular hydrogen which is not properly traced by the CO (J = 1 -> 0) line, the so-called "dark molecular gas". We use results from a magnetohydrodynamical (MHD) simulation as a model for the density structures found in a turbulent diffuse ISM with no star-formation activity. The Meudon PDR code is then applied to a number of lines of sight through this model, to derive their chemical structures. It is found that, compared to the uniform density assumption, maximal chemical abundances for H2, CO, CH and CN are increased by a factor 2 to 4 when taking into account density fluctuations on the line of sight. The correlations between column densities of CO, CH and CN with respect to those of H2 are also found to be in better overall agreement with observations. For instance, at N(H2) > 2.10^{20} cm-2, while observations suggest that d[log N(CO)]=d[log N(H2)] = 3.07 +/- 0.73, we find d[log N(CO)]=d[log N(H2)] =14 when assuming uniform density, and d[log N(CO)]=d[log N(H2)] = 5.2 when including density fluctuations.Comment: 14 pages, 16 figures, accepted for publication in Astronomy & Astrophysic

Planck intermediate results. XLI. A map of lensing-induced B-modes

The secondary cosmic microwave background (CMB) $B$-modes stem from the post-decoupling distortion of the polarization $E$-modes due to the gravitational lensing effect of large-scale structures. These lensing-induced $B$-modes constitute both a valuable probe of the dark matter distribution and an important contaminant for the extraction of the primary CMB $B$-modes from inflation. Planck provides accurate nearly all-sky measurements of both the polarization $E$-modes and the integrated mass distribution via the reconstruction of the CMB lensing potential. By combining these two data products, we have produced an all-sky template map of the lensing-induced $B$-modes using a real-space algorithm that minimizes the impact of sky masks. The cross-correlation of this template with an observed (primordial and secondary) $B$-mode map can be used to measure the lensing $B$-mode power spectrum at multipoles up to $2000$. In particular, when cross-correlating with the $B$-mode contribution directly derived from the Planck polarization maps, we obtain lensing-induced $B$-mode power spectrum measurement at a significance level of $12\,\sigma$, which agrees with the theoretical expectation derived from the Planck best-fit $\Lambda$CDM model. This unique nearly all-sky secondary $B$-mode template, which includes the lensing-induced information from intermediate to small ($10\lesssim \ell\lesssim 1000$) angular scales, is delivered as part of the Planck 2015 public data release. It will be particularly useful for experiments searching for primordial $B$-modes, such as BICEP2/Keck Array or LiteBIRD, since it will enable an estimate to be made of the lensing-induced contribution to the measured total CMB $B$-modes.Comment: 20 pages, 12 figures; Accepted for publication in A&A; The B-mode map is part of the PR2-2015 Cosmology Products; available as Lensing Products in the Planck Legacy Archive http://pla.esac.esa.int/pla/#cosmology; and described in the 'Explanatory Supplement' https://wiki.cosmos.esa.int/planckpla2015/index.php/Specially_processed_maps#2015_Lensing-induced_B-mode_ma