The characterization of the dust polarization foreground to the cosmic microwave background (CMB) is a necessary step toward the detection of the B-mode signal associated with primordial ...gravitational waves. We present a method to simulate maps of polarized dust emission on the sphere that is similar to the approach used for CMB anisotropies. This method builds on the understanding of Galactic polarization stemming from the analysis of Planck data. It relates the dust polarization sky to the structure of the Galactic magnetic field and its coupling with interstellar matter and turbulence. The Galactic magnetic field is modeled as a superposition of a mean uniform field and a Gaussian random (turbulent) component with a power-law power spectrum of exponent αM. The integration along the line of sight carried out to compute Stokes maps is approximated by a sum over a small number of emitting layers with different realizations of the random component of the magnetic field. The model parameters are constrained to fit the power spectra of dust polarization EE, BB, and TE measured using Planck data. We find that the slopes of the E and B power spectra of dust polarization are matched for αM = −2.5, an exponent close to that measured for total dust intensity but larger than the Kolmogorov exponent − 11/3. The model allows us to compute multiple realizations of the Stokes Q and U maps for different realizations of the random component of the magnetic field, and to quantify the variance of dust polarization spectra for any given sky area outside of the Galactic plane. The simulations reproduce the scaling relation between the dust polarization power and the mean total dust intensity including the observed dispersion around the mean relation. We also propose a method to carry out multifrequency simulations, including the decorrelation measured recently by Planck, using a given covariance matrix of the polarization maps. These simulations are well suited to optimize component separation methods and to quantify the confidence with which the dust and CMB B-modes can be separated in present and future experiments. We also provide an astrophysical perspective on our phenomenological modeling of the dust polarization spectra.
Context. The interstellar medium (ISM) is now widely acknowledged to display features ascribable to magnetized turbulence. With the public release of Planck data and the current balloon-borne and ...ground-based experiments, the growing amount of data tracing the polarized thermal emission from Galactic dust in the submillimetre provides choice diagnostics to constrain the properties of this magnetized turbulence. Aims. We aim to constrain these properties in a statistical way, focussing in particular on the power spectral index βB of the turbulent component of the interstellar magnetic field in a diffuse molecular cloud, the Polaris Flare. Methods. We present an analysis framework based on simulating polarized thermal dust emission maps using model dust density (proportional to gas density nH) and magnetic field cubes, integrated along the line of sight (LOS), and comparing these statistically to actual data. The model fields are derived from fractional Brownian motion (fBm) processes, which allows a precise control of their one- and two-point statistics. The parameters controlling the model are (1)–(2) the spectral indices of the density and magnetic field cubes, (3)–(4) the RMS-to-mean ratios for both fields, (5) the mean gas density, (6) the orientation of the mean magnetic field in the plane of the sky (POS), (7) the dust temperature, (8) the dust polarization fraction, and (9) the depth of the simulated cubes. We explore the nine-dimensional parameter space through a Markov chain Monte Carlo analysis, which yields best-fitting parameters and associated uncertainties. Results. We find that the power spectrum of the turbulent component of the magnetic field in the Polaris Flare molecular cloud scales with wavenumber as k−βB with a spectral index βB = 2.8 ± 0.2. It complements a uniform field whose norm in the POS is approximately twice the norm of the fluctuations of the turbulent component, and whose position angle with respect to the north-south direction is χ0 ≈−69°. The density field nH is well represented by a log-normally distributed field with a mean gas density 〈nH〉≈40 cm−3$\left\langle{n_{\mathrm{H}}}\right\rangle \approx 40\,\mathrm{cm}^{-3}$〈 nH 〉≈40 cm−3, a fluctuation ratio σnH/〈nH〉≈1.6$\sigma_{{n_{\mathrm{H}}}/\langle{n_{\mathrm{H}}}\rangle\approx 1.6$σnH/〈nH〉≈1.6 , and a power spectrum with an index βn=1.7−0.3+0.4 $\beta_n=1.7^{+0.4}_{-0.3}$βn=1.7−0.3+0.4 . We also constrain the depth of the cloud to be d ≈ 13 pc, and the polarization fraction p0 ≈ 0.12. The agreement between the Planck data and the simulated maps for these best-fitting parameters is quantified by a χ2 value that is only slightly larger than unity. Conclusions. We conclude that our fBm-based model is a reasonable description of the diffuse, turbulent, magnetized ISM in the Polaris Flare molecular cloud, and that our analysis framework is able to yield quantitative estimates of the statistical properties of the dust density and magnetic field in this cloud.
ABSTRACT
Anomalous Microwave Emission (AME) is a significant component of Galactic diffuse emission in the frequency range 10–$60\, \mathrm{GHz}$ and a new window into the properties of ...sub-nanometre-sized grains in the interstellar medium. We investigate the morphology of AME in the ≈10○ diameter λ Orionis ring by combining intensity data from the QUIJOTE experiment at 11, 13, 17, and $19\, \mathrm{GHz}$ and the C-Band All Sky Survey (C-BASS) at $4.76\, \mathrm{GHz}$, together with 19 ancillary data sets between 1.42 and $3000\, \mathrm{GHz}$. Maps of physical parameters at 1○ resolution are produced through Markov chain Monte Carlo (MCMC) fits of spectral energy distributions (SEDs), approximating the AME component with a lognormal distribution. AME is detected in excess of $20\, \sigma$ at degree-scales around the entirety of the ring along photodissociation regions (PDRs), with three primary bright regions containing dark clouds. A radial decrease is observed in the AME peak frequency from $\approx 35\, \mathrm{GHz}$ near the free–free region to $\approx 21\, \mathrm{GHz}$ in the outer regions of the ring, which is the first detection of AME spectral variations across a single region. A strong correlation between AME peak frequency, emission measure and dust temperature is an indication for the dependence of the AME peak frequency on the local radiation field. The AME amplitude normalized by the optical depth is also strongly correlated with the radiation field, giving an overall picture consistent with spinning dust where the local radiation field plays a key role.
ABSTRACT
Anomalous microwave emission (AME) is an important emission component between 10 and 60 GHz that is not yet fully understood. It seems to be ubiquitous in our Galaxy and is observed at a ...broad range of angular scales. Here we use the new QUIJOTE-MFI wide survey data at 11, 13, 17, and 19 GHz to constrain the AME in the Galactic plane (|b| < 10°) on degree scales. We built the spectral energy distribution between 0.408 and 3000 GHz for each of the 5309 0.9° pixels in the Galactic plane, and fitted a parametric model by considering five emission components: synchrotron, free–free, AME, thermal dust and CMB anisotropies. We show that not including QUIJOTE-MFI data points leads to the underestimation (up to 50 per cent) of the AME signal in favour of free–free emission. The parameters describing these components are then intercompared, looking for relations that help to understand AME physical processes. We find median values for the AME width, WAME, and for its peak frequency, νAME, respectively of $0.560^{+0.059}_{-0.050}$ and $20.7^{+2.0}_{-1.9}$ GHz, slightly in tension with current theoretical models. We find spatial variations throughout the Galactic plane for νAME, but only with reduced statistical significance. We report correlations of AME parameters with certain ISM properties, such as that between the AME emissivity (which shows variations with the Galactic longitude) and the interstellar radiation field, and that between the AME peak frequency and dust temperature. Finally, we discuss the implications of our results on the possible molecules responsible for AME.
ABSTRACT
We derive linearly polarized astrophysical component maps in the Northern Sky from the QUIJOTE-MFI data at 11 and 13 GHz in combination with the Wilkinson Microwave Anisotropy Probe K and Ka ...bands (23 and 33 GHz) and all Planck polarized channels (30–353 GHz), using the parametric component separation method B-SeCRET. The addition of QUIJOTE-MFI data significantly improves the parameter estimation of the low-frequency foregrounds, especially the estimation of the synchrotron spectral index, βs. We present the first detailed βs map of the Northern Celestial Hemisphere at a smoothing scale of 2°. We find statistically significant spatial variability across the sky. We obtain an average value of −3.08 and a dispersion of 0.13, considering only pixels with reliable goodness of fit. The power-law model of the synchrotron emission provides a good fit to the data outside the Galactic plane but fails to track the complexity within this region. Moreover, when we assume a synchrotron model with uniform curvature, cs, we find a value of cs = −0.0797 ± 0.0012. However, there is insufficient statistical significance to determine which model is favoured, either the power law or the power law with uniform curvature. Furthermore, we estimate the thermal dust spectral parameters in polarization. Our cosmic microwave background, synchrotron, and thermal dust maps are highly correlated with the corresponding products of the PR4 Planck release, although some large-scale differences are observed in the synchrotron emission. Finally, we find that the βs estimation in the high signal-to-noise synchrotron emission areas is prior-independent, while, outside these regions, the prior governs the βs estimation.
ABSTRACT
We present new intensity and polarization maps obtained with the QUIJOTE experiment towards the Galactic regions W49, W51 and IC443, covering the frequency range from 10 to 20 GHz at $\sim ...1\, \text{deg}$ angular resolution, with a sensitivity in the range 35–79 $\mu \text{K}\, \text{beam}^{-1}$ for total intensity and 13–23 $\mu \text{K}\, \text{beam}^{-1}$ for polarization. For each region, we combine QUIJOTE maps with ancillary data at frequencies ranging from 0.4 to 3000 GHz, reconstruct the spectral energy distribution and model it with a combination of known foregrounds. We detect anomalous microwave emission (AME) in total intensity towards W49 at 4.7σ and W51 at 4.0σ with peak frequencies $\nu _{\rm AME}=(20.0\pm 1.4)\, \text{GHz}$ and $\nu _{\rm AME}=(17.7\pm 3.6)\, \text{GHz}$, respectively; this is the first detection of AME towards W51. The contamination from ultracompact H ii regions to the residual AME flux density is estimated at 10 per cent in W49 and 5 per cent in W51, and does not rule out the AME detection. The polarized SEDs reveal a synchrotron contribution with spectral indices αs = −0.67 ± 0.10 in W49 and αs = −0.51 ± 0.07 in W51, ascribed to the diffuse Galactic emission and to the local supernova remnant, respectively. Towards IC443 in total intensity we measure a broken power-law synchrotron spectrum with cut-off frequency $\nu _{\rm 0,s}=(114\pm 73)\, \text{GHz}$, in agreement with previous studies; our analysis, however, rules out any AME contribution which had been previously claimed towards IC443. No evidence of polarized AME emission is detected in this study.
ABSTRACT
The QUIJOTE-MFI Northern Hemisphere Wide Survey has provided maps of the sky above declinations −30° at 11, 13, 17, and 19 GHz. These data are combined with ancillary data to produce ...Spectral Energy Distributions in intensity in the frequency range 0.4–3 000 GHz on a sample of 52 candidate compact sources harbouring anomalous microwave emission (AME). We apply a component separation analysis at 1° scale on the full sample from which we identify 44 sources with high AME significance. We explore correlations between different fitted parameters on this last sample. QUIJOTE-MFI data contribute to notably improve the characterization of the AME spectrum, and its separation from the other components. In particular, ignoring the 10–20 GHz data produces on average an underestimation of the AME amplitude, and an overestimation of the free–free component. We find an average AME peak frequency of 23.6 ± 3.6 GHz, about 4 GHz lower than the value reported in previous studies. The strongest correlation is found between the peak flux density of the thermal dust and of the AME component. A mild correlation is found between the AME emissivity (AAME/τ250) and the interstellar radiation field. On the other hand no correlation is found between the AME emissivity and the free–free radiation Emission Measure. Our statistical results suggest that the interstellar radiation field could still be the main driver of the intensity of the AME as regards spinning dust excitation mechanisms. On the other hand, it is not clear whether spinning dust would be most likely associated with cold phases of the interstellar medium rather than with hot phases dominated by free–free radiation.
ABSTRACT
The Andromeda Galaxy (M31) is the Local Group galaxy that is most similar to the Milky Way (MW). The similarities between the two galaxies make M31 useful for studying integrated properties ...common to spiral galaxies. We use the data from the recent QUIJOTE-MFI Wide Survey, together with new raster observations focused on M31, to study its integrated emission. The addition of raster data improves the sensitivity of QUIJOTE-MFI maps by almost a factor 3. Our main interest is to confirm if anomalous microwave emission (AME) is present in M31, as previous studies have suggested. To do so, we built the integrated spectral energy distribution of M31 between 0.408 and 3000 GHz. We then performed a component separation analysis taking into account synchrotron, free–free, AME, and thermal dust components. AME in M31 is modelled as a log-normal distribution with maximum amplitude, AAME, equal to 1.03 ± 0.32 Jy. It peaks at ${\nu _{\rm AME}}=17.2\pm 3.2{\rm \, GHz}{}$ with a width of WAME = 0.58 ± 0.16. Both the Akaike and Bayesian information criteria find the model without AME to be less than 1 per cent as probable as the one taking AME into consideration. We find that the AME emissivity per 100 $\mu$m intensity in M31 is ${\epsilon _{\rm AME}^{\rm 28.4\, GHz}}=9.6\pm 3.1\,\mu$K MJy−1 sr, similar to that of the MW. We also provide the first upper limits for the AME polarization fraction in an extragalactic object. M31 remains the only galaxy where an AME measurement has been made of its integrated spectrum.
ABSTRACT
The Haze is an excess of microwave intensity emission surrounding the Galactic Centre. It is spatially correlated with the γ-ray Fermi bubbles, and with the S-PASS radio polarization plumes, ...suggesting a possible common provenance. The models proposed to explain the origin of the Haze, including energetic events at the Galactic Centre and dark matter decay in the Galactic halo, do not yet provide a clear physical interpretation. In this paper, we present a reanalysis of the Haze including new observations from the Multi-Frequency Instrument (MFI) of the Q-U-I Joint TEnerife (QUIJOTE) experiment, at 11 and 13 GHz. We analyse the Haze in intensity and polarization, characterizing its spectrum. We detect an excess of diffuse intensity signal ascribed to the Haze. The spectrum at frequencies 11 GHz $\, \le \nu \le \,$ 70 GHz is a power law with spectral index βH = −2.79 ± 0.08, which is flatter than the Galactic synchrotron in the same region (βS = −2.98 ± 0.04), but steeper than that obtained from previous works (βH ∼ −2.5 at 23 GHz $\, \le \, \nu \le \,$ 70 GHz). We also observe an excess of polarized signal in the QUIJOTE-MFI maps in the Haze area. This is a first hint detection of polarized Haze, or a consequence of curvature of the synchrotron spectrum in that area. Finally, we show that the spectrum of polarized structures associated with Galactic Centre activity is steep at low frequencies (β ∼ −3.2 at 2.3 GHz ≤ ν ≤ 23 GHz), and becomes flatter above 11 GHz.
ABSTRACT
We use the new QUIJOTE-MFI wide survey (11, 13, 17, and 19 GHz) to produce spectral energy distributions (SEDs), on an angular scale of 1°, of the supernova remnants (SNRs) CTB 80, Cygnus ...Loop, HB 21, CTA 1, Tycho, and HB 9. We provide new measurements of the polarized synchrotron radiation in the microwave range. The intensity and polarization SEDs are obtained and modelled by combining QUIJOTE-MFI maps with ancillary data. In intensity, we confirm the curved spectra of CTB 80 and HB 21 with a break frequency νb at 2.0$^{+1.2}_{-0.5}$ and 5.0$^{+1.2}_{-1.0}$ GHz, respectively; and spectral indices above the break of $-0.6^{+0.04}_{-0.05}$ and $-0.86^{+0.04}_{-0.05}$. We provide constraints on the Anomalous Microwave Emission, suggesting that it is negligible towards these SNRs. From a simultaneous intensity and polarization fit, we recover synchrotron spectral indices as flat as −0.24, and the whole sample has a mean and scatter of −0.44 ± 0.12. The polarization fractions have a mean and scatter of 6.1 ± 1.9 per cent. When combining our results with the measurements from other QUIJOTE (Q-U-I JOint TEnerife CMB experiment) studies of SNRs, we find that radio spectral indices are flatter for mature SNRs, and particularly flatter for CTB 80 ($-0.24^{+0.07}_{-0.06}$) and HB 21 ($-0.34^{+0.04}_{-0.03}$). In addition, the evolution of the spectral indices against the SNRs age is modelled with a power-law function, providing an exponent −0.07 ± 0.03 and amplitude −0.49 ± 0.02 (at 10 kyr), which are conservative with respect to previous studies of our Galaxy and the Large Magellanic Cloud.