QUIJOTE scientific results - VIII. Diffuse polarized foregrounds from component separation with QUIJOTE-MFI

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, [beta]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 [beta]s estimation.We thank the staff of the Teide Observatory for invaluable assistance in the commissioning and operation of QUIJOTE. The QUIJOTE experiment is being developed by the Instituto de Astrofisica de Canarias (IAC), the Instituto de Fisica de Cantabria (IFCA), and the Universities of Cantabria, Manchester, and Cambridge. Partial financial support was provided by the Spanish Ministry of Science and Innovation under the projects AYA2007-68058-C03-01, AYA2007- 68058-C03-02, AYA2010-21766-C03-01, AYA2010-21766-C03-02, AYA2014-60438-P, ESP2015-70646-C2-1-R, AYA2017-84185-P, ESP2017-83921-C2-1-R, AYA2017-90675-REDC (co-funded with EU FEDER funds), PGC2018-101814-B-I00, PID2019-110610RBC21, PID2020-120514GB-I00, IACA13-3E-2336, IACA15-BE3707, EQC2018-004918-P, the Severo Ochoa Programs SEV-2015- 0548 and CEX2019-000920-S, the Maria de Maeztu Program MDM2017-0765, and by the Consolider-Ingenio project CSD2010-00064 (EPI: Exploring the Physics of Inflation). We acknowledge support from the ACIISI, Consejeria de Economia, Conocimiento y Empleo del Gobierno de Canarias, and the European Regional Development Fund (ERDF) under grant with reference ProID2020010108. This project has received funding from the European Union’s Horizon 2020 research and innovation programme under grant agreement number 687312 (RADIOFOREGROUNDS). EdlH acknowledges financial support from the Concepcion´ Arenal Programme of the Universidad de Cantabria. DT acknowledges the support from the Chinese Academy of Sciences (CAS) President’s International Fellowship Initiative (PIFI) with grant no. 2020PM0042. FP acknowledges support from the Spanish State Research Agency (AEI) under grant number PID2019-105552RB-C43. The authors acknowledge the computer resources, technical expertise, and assistance provided by the Spanish Supercomputing Network (RES) node at Universidad de Cantabria. Some of the presented results are based on observations obtained with Planck (http://www.esa.int/Planck), an ESA science mission with instruments and contributions directly funded by ESA Member States, NASA, and Canada. We acknowledge the use of the Legacy Archive for Microwave Background Data Analysis (LAMBDA) and the Planck Legacy Archive (PLA). Support for LAMBDA is provided by the NASA Office of Space Science. Some of the results in this paper have been derived using the HEALPIX package (Gorski ´ et al. 2005), and the HEALPY (Zonca et al. 2019), NUMPY (Harris et al. 2020), EMCEE (ForemanMackey et al. 2013), and MATPLOTLIB (Hunter 2007) PYTHON packages

Planck 2018 results. IV. Diffuse component separation

We present full-sky maps of the cosmic microwave background (CMB) and polarized synchrotron and thermal dust emission, derived from the third set of Planck frequency maps. These products have significantly lower contamination from instrumental systematic effects than previous versions. The methodologies used to derive these maps follow closely those described in earlier papers, adopting four methods (Commander, NILC, SEVEM, and SMICA) to extract the CMB component, as well as three methods (Commander, GNILC, and SMICA) to extract astrophysical components. Our revised CMB temperature maps agree with corresponding products in the Planck 2015 delivery, whereas the polarization maps exhibit significantly lower large-scale power, reflecting the improved data processing described in companion papers; however, the noise properties of the resulting data products are complicated, and the best available end-to-end simulations exhibit relative biases with respect to the data at the few percent level. Using these maps, we are for the first time able to fit the spectral index of thermal dust independently over 3 degree regions. We derive a conservative estimate of the mean spectral index of polarized thermal dust emission of beta_d = 1.55 +/- 0.05, where the uncertainty marginalizes both over all known systematic uncertainties and different estimation techniques. For polarized synchrotron emission, we find a mean spectral index of beta_s = -3.1 +/- 0.1, consistent with previously reported measurements. We note that the current data processing does not allow for construction of unbiased single-bolometer maps, and this limits our ability to extract CO emission and correlated components. The foreground results for intensity derived in this paper therefore do not supersede corresponding Planck 2015 products. For polarization the new results supersede the corresponding 2015 products in all respects

Planck 2018 results. XII. Galactic astrophysics using polarized dust emission

We present 353 GHz full-sky maps of the polarization fraction $p$, angle $\psi$, and dispersion of angles $S$ of Galactic dust thermal emission produced from the 2018 release of Planck data. We confirm that the mean and maximum of $p$ decrease with increasing $N_H$. The uncertainty on the maximum polarization fraction, $p_\mathrm{max}=22.0$% at 80 arcmin resolution, is dominated by the uncertainty on the zero level in total intensity. The observed inverse behaviour between $p$ and $S$ is interpreted with models of the polarized sky that include effects from only the topology of the turbulent Galactic magnetic field. Thus, the statistical properties of $p$, $\psi$, and $S$ mostly reflect the structure of the magnetic field. Nevertheless, we search for potential signatures of varying grain alignment and dust properties. First, we analyse the product map $S \times p$, looking for residual trends. While $p$ decreases by a factor of 3--4 between $N_H=10^{20}$ cm$^{-2}$ and $N_H=2\times 10^{22}$ cm$^{-2}$, $S \times p$ decreases by only about 25%, a systematic trend observed in both the diffuse ISM and molecular clouds. Second, we find no systematic trend of $S \times p$ with the dust temperature, even though in the diffuse ISM lines of sight with high $p$ and low $S$ tend to have colder dust. We also compare Planck data with starlight polarization in the visible at high latitudes. The agreement in polarization angles is remarkable. Two polarization emission-to-extinction ratios that characterize dust optical properties depend only weakly on $N_H$ and converge towards the values previously determined for translucent lines of sight. We determine an upper limit for the polarization fraction in extinction of 13%, compatible with the $p_\mathrm{max}$ observed in emission. These results provide strong constraints for models of Galactic dust in diffuse gas

Biases in the estimation of velocity dispersions and dynamical masses for galaxy clusters

Using a set of 73 numerically simulated galaxy clusters, we have characterised the statistical and physical biases for three velocity dispersion and mass estimators, namely biweight, gapper and standard deviation, in the small number of galaxies regime (Ngal ≤ 75), both for the determination of the velocity dispersion and the dynamical mass of the clusters via the σ–M relation. These results are used to define a new set of unbiased estimators, that are able to correct for those statistical biases. By applying these new estimators to a subset of simulated observations, we show that they can retrieve bias-corrected values for both the mean velocity dispersion and the mean mass

Biases in the estimation of velocity dispersions and dynamical masses for galaxy clusters

Using a set of 73 numerically simulated galaxy clusters, we have characterised the statistical and physical biases for three velocity dispersion and mass estimators, namely biweight, gapper and standard deviation, in the small number of galaxies regime (Ngal ≤ 75), both for the determination of the velocity dispersion and the dynamical mass of the clusters via the σ–M relation. These results are used to define a new set of unbiased estimators, that are able to correct for those statistical biases. By applying these new estimators to a subset of simulated observations, we show that they can retrieve bias-corrected values for both the mean velocity dispersion and the mean mass

Optical Identifications of High-Redshift Galaxy Clusters from the Planck Sunyaev–Zeldovich Survey

International audienceWe present the results of optical identifications and spectroscopic redshift measurements for galaxy clusters from the second Planck catalogue of Sunyaev–Zeldovich sources (PSZ2) located at high redshifts, z ≈ 0.7−0.9. We used the data of optical observations with the Russian–Turkish 1.5-mtelescope (RTT-150), the Sayan Observatory 1.6-m telescope, the Calar Alto 3.5-m telescope, and the 6-m SAO RAS telescope (BTA). The spectroscopic redshift measurements were obtained for seven galaxy clusters, including one cluster, PSZ2 G126.57+51.61, from the cosmological sample of the PSZ2 catalogue. In the central regions of two clusters, PSZ2 G069.39+68.05 and PSZ2 G087.39−34.58, we detected arcs of strong gravitational lensing of background galaxies, one of which is at redshift z = 4.262. The data presented below roughly double the number of known galaxy clusters in the second Planck catalogue of Sunyaev–Zeldovich sources at high redshifts, z ≈ 0.8

QUIJOTE scientific results -- XIII. Intensity and polarization study of supernova remnants in the QUIJOTE-MFI wide survey: CTB 80, Cygnus Loop, HB 21, CTA 1, Tycho and HB 9

International audienceWe 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 deg, 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. For each SNR, the intensity and polarization SEDs are obtained and modelled by combining QUIJOTE-MFI maps with ancillary data. In intensity, we confirm the curved power law spectra of CTB 80 and HB 21 with a break frequency $\nu_{\rm b}$ at 2.0$^{+1.2}_{-0.5}$ GHz and 5.0$^{+1.2}_{-1.0}$ GHz respectively; and spectral indices respectively below and above the spectral break of $-0.34\pm0.04$ and $-0.86\pm0.5$ for CTB 80, and $-0.24\pm0.07$ and $-0.60\pm0.05$ for HB 21. In addition, we provide upper limits on the Anomalous Microwave Emission (AME), suggesting that the AME contribution is negligible towards these remnants. 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\pm0.12$. The polarization fractions have a mean and scatter of $6.1\pm1.9$%. When combining our results with the measurements from other QUIJOTE 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\pm0.03$ and amplitude $-0.49\pm0.02$ (normalised at 10 kyr), which are conservative with respect to previous studies of our Galaxy and the Large Magellanic Cloud

QUIJOTE scientific results -- XIII. Intensity and polarization study of supernova remnants in the QUIJOTE-MFI wide survey: CTB 80, Cygnus Loop, HB 21, CTA 1, Tycho and HB 9

International audienceWe 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 deg, 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. For each SNR, the intensity and polarization SEDs are obtained and modelled by combining QUIJOTE-MFI maps with ancillary data. In intensity, we confirm the curved power law spectra of CTB 80 and HB 21 with a break frequency $\nu_{\rm b}$ at 2.0$^{+1.2}_{-0.5}$ GHz and 5.0$^{+1.2}_{-1.0}$ GHz respectively; and spectral indices respectively below and above the spectral break of $-0.34\pm0.04$ and $-0.86\pm0.5$ for CTB 80, and $-0.24\pm0.07$ and $-0.60\pm0.05$ for HB 21. In addition, we provide upper limits on the Anomalous Microwave Emission (AME), suggesting that the AME contribution is negligible towards these remnants. 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\pm0.12$. The polarization fractions have a mean and scatter of $6.1\pm1.9$%. When combining our results with the measurements from other QUIJOTE 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\pm0.03$ and amplitude $-0.49\pm0.02$ (normalised at 10 kyr), which are conservative with respect to previous studies of our Galaxy and the Large Magellanic Cloud

QUIJOTE Scientific Results -- XVII. Studying the Anomalous Microwave Emission in the Andromeda Galaxy with QUIJOTE-MFI

International audienceThe 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 a factor greater than 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, $A_{\rm AME}$, equal to $1.06\pm0.30$ Jy. It peaks at $\nu_{\rm AME}=17.28\pm3.08$ GHz with a width of $W_{\rm AME}=0.57\pm0.15$. Both the Akaike and Bayesian Information Criteria find the model without AME to be less than 1 % as probable as the one taking AME into consideration, thus strongly favouring the presence of AME in M31. We find that the AME emissivity in M31 is $\epsilon_{\rm AME}^{\rm 28.4\,GHz}=9.1\pm2.9$ $\mu$K/(MJy/sr), similar to that computed for 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

QUIJOTE Scientific Results -- XVII. Studying the Anomalous Microwave Emission in the Andromeda Galaxy with QUIJOTE-MFI

International audienceThe 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 a factor greater than 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, $A_{\rm AME}$, equal to $1.06\pm0.30$ Jy. It peaks at $\nu_{\rm AME}=17.28\pm3.08$ GHz with a width of $W_{\rm AME}=0.57\pm0.15$. Both the Akaike and Bayesian Information Criteria find the model without AME to be less than 1 % as probable as the one taking AME into consideration, thus strongly favouring the presence of AME in M31. We find that the AME emissivity in M31 is $\epsilon_{\rm AME}^{\rm 28.4\,GHz}=9.1\pm2.9$ $\mu$K/(MJy/sr), similar to that computed for 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