Search

Future cosmic microwave background (CMB) experiments are primarily targeting a detection of the primordial $B$-mode polarisation. The faintness of this signal requires exquisite control of systematic effects which may bias the measurements. In this work, we derive requirements on the relative calibration accuracy of the overall polarisation gain ($\Delta g_\nu$) for LiteBIRD experiment, through the application of the blind Needlet Internal Linear Combination (NILC) foreground-cleaning method. We find that minimum variance techniques, as NILC, are less affected by gain calibration uncertainties than a parametric approach, which requires a proper modelling of these instrumental effects. The tightest constraints are obtained for frequency channels where the CMB signal is relatively brighter (166 GHz channel, $\Delta {g}_\nu \approx 0.16 \%$), while, with a parametric approach, the strictest requirements were on foreground-dominated channels. We then propagate gain calibration uncertainties, corresponding to the derived requirements, into all frequency channels simultaneously. We find that the overall impact on the estimated $r$ is lower than the required budget for LiteBIRD by almost a factor $5$. The adopted procedure to derive requirements assumes a simple Galactic model. We therefore assess the robustness of obtained results against more realistic scenarios by injecting the gain calibration uncertainties, according to the requirements, into LiteBIRD simulated maps and assuming intermediate- and high-complexity sky models. In this case, we employ the so-called Multi-Clustering NILC (MC-NILC) foreground-cleaning pipeline and obtain that the impact of gain calibration uncertainties on $r$ is lower than the LiteBIRD gain systematics budget for the intermediate-complexity sky model. For the high-complexity case, instead, it would be necessary to tighten the requirements by a factor $1.8$., Comment: 29 pages, 11 figures

This work focuses on the study of the AME, an important emission mechanism between 10 and 60 GHz whose polarization properties are not yet fully understood, and is therefore a potential contaminant for future CMB polarization observations. We use new QUIJOTE-MFI maps 11, 13, 17 and 19 GHz, together with other public ancillary data including WMAP and Planck, to study the polarization properties of the AME in three Galactic regions: rho-Ophiuchi, Perseus and W43. We have obtained the SEDs for those three regions over the frequency range 0.4-3000 GHz, both in intensity and polarization. The intensity SEDs are well described by a combination of free-free emission, thermal dust, AME and CMB anisotropies. In polarization, we extracted the flux densities using all available data between 11 and 353 GHz. We implemented an improved intensity-to-polarization leakage correction that has allowed for the first time to derive reliable polarization constraints well below the 1% level from Planck-LFI data. A frequency stacking of maps in the range 10-60 GHz has allowed us to reduce the statistical noise and to push the upper limits on the AME polarization level. We have obtained upper limits on the AME polarization fraction of order <1% (95% confidence level) for the three regions. In particular we get Pi_AME < 1.1% (at 28.4 GHz), Pi_AME < 1.1% (at 22.8 GHz) and Pi_AME < 0.28% (at 33 GHz) in rho-Ophiuchi, Perseus and W43 respectively. At the QUIJOTE 17 GHz frequency band, we get Pi_AME< 5.1% for rho-Ophiuchi, Pi_AME< 3.5% for Perseus and Pi_AME< 0.85% for W43. Our final upper limits derived using the stacking procedure are Pi_AME < 0.58% for rho-Ophiuchi, Pi_AME < 1.64% for Perseus and Pi_AME < 0.31% for W43. Altogether, these are the most stringent constraints to date on the AME polarization fraction of these three star-forming regions., Comment: 26 pages. 15 figures. Submitted to A&A. Comments welcome

Large angular scale surveys in the absence of atmosphere are essential for measuring the primordial $B$-mode power spectrum of the Cosmic Microwave Background (CMB). Since this proposed measurement is about three to four orders of magnitude fainter than the temperature anisotropies of the CMB, in-flight calibration of the instruments and active suppression of systematic effects are crucial. We investigate the effect of changing the parameters of the scanning strategy on the in-flight calibration effectiveness, the suppression of the systematic effects themselves, and the ability to distinguish systematic effects by null-tests. Next-generation missions such as LiteBIRD, modulated by a Half-Wave Plate (HWP), will be able to observe polarisation using a single detector, eliminating the need to combine several detectors to measure polarisation, as done in many previous experiments and hence avoiding the consequent systematic effects. While the HWP is expected to suppress many systematic effects, some of them will remain. We use an analytical approach to comprehensively address the mitigation of these systematic effects and identify the characteristics of scanning strategies that are the most effective for implementing a variety of calibration strategies in the multi-dimensional space of common spacecraft scan parameters. We also present Falcons, a fast spacecraft scanning simulator that we developed to investigate this scanning parameter space.

We assess the capabilities of the LiteBIRD mission to map the hot gas distribution in the Universe through the thermal Sunyaev-Zeldovich (SZ) effect. Our analysis relies on comprehensive simulations incorporating various sources of Galactic and extragalactic foreground emission, while accounting for specific instrumental characteristics of LiteBIRD, such as detector sensitivities, frequency-dependent beam convolution, inhomogeneous sky scanning, and $1/f$ noise. We implement a tailored component-separation pipeline to map the thermal SZ Compton $y$-parameter over 98% of the sky. Despite lower angular resolution for galaxy cluster science, LiteBIRD provides full-sky coverage and, compared to the Planck satellite, enhanced sensitivity, as well as more frequency bands to enable the construction of an all-sky $y$-map, with reduced foreground contamination at large and intermediate angular scales. By combining LiteBIRD and Planck channels in the component-separation pipeline, we obtain an optimal $y$-map that leverages the advantages of both experiments, with the higher angular resolution of the Planck channels enabling the recovery of compact clusters beyond the LiteBIRD beam limitations, and the numerous sensitive LiteBIRD channels further mitigating foregrounds. The added value of LiteBIRD is highlighted through the examination of maps, power spectra, and one-point statistics of the various sky components. After component separation, the $1/f$ noise from LiteBIRD is effectively mitigated below the thermal SZ signal at all multipoles. Cosmological constraints on $S_8=\sigma_8\left(\Omega_{\rm m}/0.3\right)^{0.5}$ obtained from the LiteBIRD-Planck combined $y$-map power spectrum exhibits a 15% reduction in uncertainty compared to constraints from Planck alone. This improvement can be attributed to the increased portion of uncontaminated sky available in the LiteBIRD-Planck combined $y$-map., Comment: 38 pages, 13 figures, abstract shortened. Updated to match version accepted by JCAP

LiteBIRD, the next-generation cosmic microwave background (CMB) experiment, aims for a launch in Japan's fiscal year 2032, marking a major advancement in the exploration of primordial cosmology and fundamental physics. Orbiting the Sun-Earth Lagrangian point L2, this JAXA-led strategic L-class mission will conduct a comprehensive mapping of the CMB polarization across the entire sky. During its 3-year mission, LiteBIRD will employ three telescopes within 15 unique frequency bands (ranging from 34 through 448 GHz), targeting a sensitivity of 2.2\,$\mu$K-arcmin and a resolution of 0.5$^\circ$ at 100\,GHz. Its primary goal is to measure the tensor-to-scalar ratio $r$ with an uncertainty $\delta r = 0.001$, including systematic errors and margin. If $r \geq 0.01$, LiteBIRD expects to achieve a $>5\sigma$ detection in the $\ell=$2-10 and $\ell=$11-200 ranges separately, providing crucial insight into the early Universe. We describe LiteBIRD's scientific objectives, the application of systems engineering to mission requirements, the anticipated scientific impact, and the operations and scanning strategies vital to minimizing systematic effects. We will also highlight LiteBIRD's synergies with concurrent CMB projects., Comment: 23 pages, 9 figures, 1 table, SPIE Astronomical Telescopes + Instrumentation 2024

The presence of astrophysical emissions in microwave observations forces us to perform component separation to extract the Cosmic Microwave Background (CMB) signal. However, even in the most optimistic cases, there are still strongly contaminated regions, such as the Galactic plane or those with emission from extragalactic point sources, which require the use of a mask. Since many CMB analyses, especially the ones working in harmonic space, need the whole sky map, it is crucial to develop a reliable inpainting algorithm that replaces the values of the excluded pixels by others statistically compatible with the rest of the sky. This is especially important when working with $Q$ and $U$ sky maps in order to obtain $E$- and $B$-mode maps which are free from $E$-to-$B$ leakage. In this work we study a method based on Gaussian Constrained Realizations (GCR), that can deal with both intensity and polarization. Several tests have been performed to asses the validation of the method, including the study of the one-dimensional probability distribution function (1-PDF), E- and B-mode map reconstruction, and power spectra estimation. We have considered two scenarios for the input simulation: one case with only CMB signal and a second one including also Planck PR4 semi-realistic noise. Even if we are limited to low resolution maps, $N_{\mathrm{side}} = $ 64 if $T$, $Q$ and $U$ are considered, we believe that this is a useful approach to be applied to future missions such as LiteBIRD, where the target are the largest scales., Comment: 26 pages, 24 figures. Minor changes in the text. New Appendix B. Updated to match the final version published in JCAP

We present detailed forecasts for the constraints on primordial magnetic fields (PMFs) that will be obtained with the LiteBIRD satellite. The constraints are driven by the effects of PMFs on the CMB anisotropies: the gravitational effects of magnetically-induced perturbations; the effects on the thermal and ionization history of the Universe; the Faraday rotation imprint on the CMB polarization; and the non-Gaussianities induced in polarization anisotropies. LiteBIRD represents a sensitive probe for PMFs and by exploiting all the physical effects, it will be able to improve the current limit coming from Planck. In particular, thanks to its accurate $B$-mode polarization measurement, LiteBIRD will improve the constraints on infrared configurations for the gravitational effect, giving $B_{\rm 1\,Mpc}^{n_{\rm B} =-2.9} < 0.8$ nG at 95% C.L., potentially opening the possibility to detect nanogauss fields with high significance. We also observe a significant improvement in the limits when marginalized over the spectral index, $B_{1\,{\rm Mpc}}^{\rm marg}< 2.2$ nG at 95% C.L. From the thermal history effect, which relies mainly on $E$-mode polarization data, we obtain a significant improvement for all PMF configurations, with the marginalized case, $\sqrt{\langle B^2\rangle}^{\rm marg}<0.50$ nG at 95% C.L. Faraday rotation constraints will take advantage of the wide frequency coverage of LiteBIRD and the high sensitivity in $B$ modes, improving the limits by orders of magnitude with respect to current results, $B_{1\,{\rm Mpc}}^{n_{\rm B} =-2.9} < 3.2$ nG at 95% C.L. Finally, non-Gaussianities of the $B$-mode polarization can probe PMFs at the level of 1 nG, again significantly improving the current bounds from Planck. Altogether our forecasts represent a broad collection of complementary probes, providing conservative limits on PMF characteristics that will be achieved with LiteBIRD., Comment: 51 pages, 24 figures, abstract shortened

The lack of power anomaly is an unexpected feature observed at large angular scales in the CMB maps produced by the COBE, WMAP and Planck satellites. This signature, which consists in a missing of power with respect to that predicted by the LCDM model, might hint at a new cosmological phase before the standard inflationary era. The main point of this paper is taking into account the latest Planck polarisation data to investigate how the CMB polarisation improves the understanding of this feature. With this aim, we apply to the latest Planck data, both PR3 (2018) and PR4 (2020) releases, a new class of estimators capable of evaluating this anomaly by considering temperature and polarisation data both separately and in a jointly way. This is the first time that the PR4 dataset has been used to study this anomaly. To critically evaluate this feature, taking into account the residuals of known systematic effects present in the Planck datasets, we analyse the cleaned CMB maps using different combinations of sky masks, harmonic range and binning on the CMB multipoles. Our analysis shows that the estimator based only on temperature data confirms the presence of a lack of power with a lower-tail-probability (LTP), depending on the component separation method, $\leq 0.33\%$ and $\leq 1.76\%$ for PR3 and PR4, respectively. To our knowledge, the LTP$\leq 0.33\%$ for the PR3 dataset is the lowest one present in the literature obtained from Planck 2018 data, considering the Planck confidence mask. We find significant differences between these two datasets when polarisation is taken into account. Moreover, we also show that for the PR3 dataset the inclusion of the subdominant polarisation information provides estimates that are less likely accepted in a LCDM cosmological model than the only-temperature analysis over the entire harmonic-range considered.

We present a study of the impact of an uncertainty in the beam far side-lobe knowledge on the measurement of the Cosmic Microwave Background $B$-mode signal at large scale. It is expected to be one of the main source of systematic effects in future CMB observations. Because it is crucial for all-sky survey missions to take into account the interplays between beam systematic effects and all the data analysis steps, the primary goal of this paper is to provide the methodology to carry out the end-to-end study of their effect for a space-borne CMB polarization experiment, up to the cosmological results in the form of a bias $\delta r$ on the tensor-to-scalar ratio $r$. LiteBIRD is dedicated to target the measurement of CMB primordial $B$ modes by reaching a sensitivity of $\sigma \left( r \right) \leq 10^{-3}$ assuming $r=0$. As a demonstration of our framework, we derive the relationship between the knowledge of the beam far side-lobes and the tentatively allocated error budget under given assumptions on design, simulation and component separation method. We assume no mitigation of the far side-lobes effect at any stage of the analysis pipeline. We show that $\delta r$ is mostly due to the integrated fractional power difference between the estimated beams and the true beams in the far side-lobes region, with little dependence on the actual shape of the beams, for low enough $\delta r$. Under our set of assumptions, in particular considering the specific foreground cleaning method we used, we find that the integrated fractional power in the far side-lobes should be known at a level as tight as $\sim 10^{-4}$, to achieve the required limit on the bias $\delta r < 1.9 \times 10^{-5}$. The framework and tools developed for this study can be easily adapted to provide requirements under different design, data analysis frameworks and for other future space-borne experiments beyond LiteBIRD.

We estimate the efficiency of mitigating the lensing $B$-mode polarization, the so-called delensing, for the $LiteBIRD$ experiment with multiple external data sets of lensing-mass tracers. The current best bound on the tensor-to-scalar ratio, $r$, is limited by lensing rather than Galactic foregrounds. Delensing will be a critical step to improve sensitivity to $r$ as measurements of $r$ become more and more limited by lensing. In this paper, we extend the analysis of the recent $LiteBIRD$ forecast paper to include multiple mass tracers, i.e., the CMB lensing maps from $LiteBIRD$ and CMB-S4-like experiment, cosmic infrared background, and galaxy number density from $Euclid$- and LSST-like survey. We find that multi-tracer delensing will further improve the constraint on $r$ by about $20\%$. In $LiteBIRD$, the residual Galactic foregrounds also significantly contribute to uncertainties of the $B$-modes, and delensing becomes more important if the residual foregrounds are further reduced by an improved component separation method., Comment: 21 pages, 7 figures

We explore the capability of measuring lensing signals in $LiteBIRD$ full-sky polarization maps. With a $30$ arcmin beam width and an impressively low polarization noise of $2.16\,\mu$K-arcmin, $LiteBIRD$ will be able to measure the full-sky polarization of the cosmic microwave background (CMB) very precisely. This unique sensitivity also enables the reconstruction of a nearly full-sky lensing map using only polarization data, even considering its limited capability to capture small-scale CMB anisotropies. In this paper, we investigate the ability to construct a full-sky lensing measurement in the presence of Galactic foregrounds, finding that several possible biases from Galactic foregrounds should be negligible after component separation by harmonic-space internal linear combination. We find that the signal-to-noise ratio of the lensing is approximately $40$ using only polarization data measured over $90\%$ of the sky. This achievement is comparable to $Planck$'s recent lensing measurement with both temperature and polarization and represents a four-fold improvement over $Planck$'s polarization-only lensing measurement. The $LiteBIRD$ lensing map will complement the $Planck$ lensing map and provide several opportunities for cross-correlation science, especially in the northern hemisphere.

We study the possibility of using the $LiteBIRD$ satellite $B$-mode survey to constrain models of inflation producing specific features in CMB angular power spectra. We explore a particular model example, i.e. spectator axion-SU(2) gauge field inflation. This model can source parity-violating gravitational waves from the amplification of gauge field fluctuations driven by a pseudoscalar "axionlike" field, rolling for a few e-folds during inflation. The sourced gravitational waves can exceed the vacuum contribution at reionization bump scales by about an order of magnitude and can be comparable to the vacuum contribution at recombination bump scales. We argue that a satellite mission with full sky coverage and access to the reionization bump scales is necessary to understand the origin of the primordial gravitational wave signal and distinguish among two production mechanisms: quantum vacuum fluctuations of spacetime and matter sources during inflation. We present the expected constraints on model parameters from $LiteBIRD$ satellite simulations, which complement and expand previous studies in the literature. We find that $LiteBIRD$ will be able to exclude with high significance standard single-field slow-roll models, such as the Starobinsky model, if the true model is the axion-SU(2) model with a feature at CMB scales. We further investigate the possibility of using the parity-violating signature of the model, such as the $TB$ and $EB$ angular power spectra, to disentangle it from the standard single-field slow-roll scenario. We find that most of the discriminating power of $LiteBIRD$ will reside in $BB$ angular power spectra rather than in $TB$ and $EB$ correlations., Comment: 22 pages, 13 figures. Submitted to JCAP

In this paper we present the European Low Frequency Survey (ELFS), a project that will enable foregrounds-free measurements of primordial $B$-mode polarization to a level 10$^{-3}$ by measuring the Galactic and extra-Galactic emissions in the 5--120\,GHz frequency window. Indeed, the main difficulty in measuring the B-mode polarization comes not just from its sheer faintness, but from the fact that many other objects in the Universe also emit polarized microwaves, which mask the faint CMB signal. The first stage of this project will be carried out in synergy with the Simons Array (SA) collaboration, installing a 5.5--11 GHz coherent receiver at the focus of one of the three 3.5\,m SA telescopes in Atacama, Chile ("ELFS on SA"). The receiver will be equipped with a fully digital back-end based on the latest Xilinx RF System-on-Chip devices that will provide frequency resolution of 1\,MHz across the whole observing band, allowing us to clean the scientific signal from unwanted radio frequency interference, particularly from low-Earth orbit satellite mega-constellations. This paper reviews the scientific motivation for ELFS and its instrumental characteristics, and provides an update on the development of ELFS on SA., Comment: to appear in Proc. of the mm Universe 2023 conference, Grenoble (France), June 2023, published by F. Mayet et al. (Eds), EPJ Web of conferences, EDP Sciences

In 2015, the Planck Collaboration released an all-sky map of the thermal Sunyaev-Zeldovich (SZ) effect, obtained by implementing the Needlet Internal Linear Combination (NILC) method on the Planck PR2 data. The quality of the Planck data has significantly improved since then. The Planck PR4 data release offers upgraded full-sky maps in the LFI and HFI frequency bands with improved systematics and sensitivity. We present a new all-sky thermal SZ Compton $y$-parameter map derived from the Planck PR4 data using NILC and highlight improvements, particularly in noise reduction and handling residual foreground contamination. The PR4 NILC Compton $y$-parameter map has been made publicly available to support further analyses., Comment: 6 pages, 6 figures; Minor changes to match the version published in Proc. of the mm Universe 2023 conference, Grenoble (France), June 2023, published by F. Mayet et al. (Eds), EPJ Web of conferences, EDP Sciences

In this article we examine recent developments in the research area concerning the creation of end-to-end models for the complete optimization of measuring instruments. The models we consider rely on differentiable programming methods and on the specification of a software pipeline including all factors impacting performance -- from the data-generating processes to their reconstruction and the extraction of inference on the parameters of interest of a measuring instrument -- along with the careful specification of a utility function well aligned with the end goals of the experiment. Building on previous studies originated within the MODE Collaboration, we focus specifically on applications involving instruments for particle physics experimentation, as well as industrial and medical applications that share the detection of radiation as their data-generating mechanism., Comment: 70 pages, 17 figures. To be submitted to journal

One of the foundations of the Standard Model of Cosmology is statistical isotropy, which can be tested, among other probes, through the study of the Cosmic Microwave Background (CMB). However, a hemispherical power asymmetry on large scales has been reported for WMAP and Planck data by different works. The statistical significance is above 3${\sigma}$ for temperature, suggesting a directional dependence of the local power spectrum, and thus a feature beyond the ${\Lambda}$CDM model. With the third release of the Planck data (PR3), a new analysis was performed including the E-mode polarization maps, finding an asymmetry at a modest level of significance. In this work, we perform an asymmetry analysis in intensity and polarization maps for the latest Planck processing pipeline (PR4). We obtain similar results to those obtained with PR3, with a slightly lower significance (2.8% for the Sevem method) for the amplitude of the E-mode local variance dipole as well as a significant variability with the considered mask. In addition, a hint of a possible T-E alignment between the asymmetry axes is found at the level of $\sim$ 5%. For the analysis, we have implemented an alternative inpainting approach in order to get an accurate reconstruction of the E-modes. More sensitive all-sky CMB polarization data, such as those expected from the future LiteBIRD experiment, are needed to reach a more robust conclusion on the possible existence of deviations from statistical isotropy in the form of a hemispherical power asymmetry., Comment: 26 pages, 17 figure. Minor changes in the discussion, conclusions and results unchanged. Updated to match the final version published in JCAP

Taking advantage of the reduced levels of noise and systematics in the data of the latest Planck release (PR4, also known as NPIPE), we construct a new all-sky Compton-$y$ parameter map (hereafter, $y$-map) of the thermal Sunyaev-Zeldovich (SZ) effect from the Planck PR4 data. A tailored Needlet Internal Linear Combination (NILC) pipeline, first validated on detailed sky simulations, is applied to the nine single-frequency Planck PR4 sky maps, ranging from $30$ to $857$ GHz, to produce the PR4 $y$-map over 98% of the sky. Using map comparisons, angular power spectra and one-point statistics we show that the PR4 NILC $y$-map is of improved quality compared to that of the previous PR2 release. The new $y$-map shows reduced levels of large-scale striations associated with $1/f$ noise in the scan direction. Regions near the Galactic plane also show lower residual contamination by Galactic thermal dust emission. At small angular scales, the residual contamination by thermal noise and cosmic infrared background (CIB) emission is found to be reduced by around 7% and 34%, respectively, in the PR4 $y$-map. The PR4 NILC $y$-map is made publicly available for astrophysical and cosmological analyses of the thermal SZ effect., Comment: 17 pages, 17 figures. Public data products available in External Maps section of Planck Legacy Archive (see https://pla.esac.esa.int/#maps ) and at Zenodo (see https://doi.org/10.5281/zenodo.7940376 ). Updated to match the version published by MNRAS

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, $A_{\rm AME}$, equal to $1.03\pm0.32$ Jy. It peaks at $\nu_{\rm AME}=17.2\pm3.2$ GHz with a width of $W_{\rm AME}=0.58\pm0.16$. 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. We find that the AME emissivity per 100 $\mu$m intensity in M31 is $\epsilon_{\rm AME}^{\rm 28.4\,GHz}=9.6\pm3.1$ $\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., Comment: 17 pages, 9 figures. Accepted for publication in MNRAS. QUIJOTE data maps available at https://research.iac.es/proyecto/quijote

Anomalous Microwave Emission (AME) is an important emission component between 10 and 60 GHz that is not yet fully understood. It seems to be ubiquituous 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^\circ$) on degree scales. We built the spectral energy distribution between 0.408 and 3000 GHz for each of the 5309 0.9$^\circ$ 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 %) 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, $W_{\rm AME}$, and for its peak frequency, $\nu_{\rm 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 $\nu_{\rm 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., Comment: 32 pages, 31 figures. Accepted for publication in MNRAS. Derived data products available at https://research.iac.es/proyecto/quijote

We propose a machine learning approach to the blind detection of extragalactic point sources on maps of the temperature anisotropies of the cosmic microwave background. Using realistic simulations of the microwave sky as seen by Planck, we train a convolutional neural network (CNN) that solves source detection as an image segmentation problem. We divide the sky into regions of progressively increasing Galactic foreground intensity and independently train specialized CNNs for each region. This strategy leads to promising levels of completeness and reliability, with our CNN substantially outperforming traditional detection methods like the matched filter in regions close to the Galactic plane., Comment: 3 pages, 1 figure. Proceedings of the "Machine Learning: a giant leap towards space discovery in the era of peta and exabyte scale surveys" Symposia of the 2022 Annual meeting of the European Astronomical Society. To be published in Memorie della SAIt

LiteBIRD is a planned JAXA-led CMB B-mode satellite experiment aiming for launch in the late 2020s, with a primary goal of detecting the imprint of primordial inflationary gravitational waves. Its current baseline focal-plane configuration includes 15 frequency bands between 40 and 402 GHz, fulfilling the mission requirements to detect the amplitude of gravitational waves with the total uncertainty on the tensor-to-scalar ratio, $\delta r$, down to $\delta r<0.001$. A key aspect of this performance is accurate astrophysical component separation, and the ability to remove polarized thermal dust emission is particularly important. In this paper we note that the CMB frequency spectrum falls off nearly exponentially above 300 GHz relative to the thermal dust SED, and a relatively minor high frequency extension can therefore result in even lower uncertainties and better model reconstructions. Specifically, we compare the baseline design with five extended configurations, while varying the underlying dust modeling, in each of which the HFT (High-Frequency Telescope) frequency range is shifted logarithmically towards higher frequencies, with an upper cutoff ranging between 400 and 600 GHz. In each case, we measure the tensor-to-scalar ratio $r$ uncertainty and bias using both parametric and minimum-variance component-separation algorithms. When the thermal dust sky model includes a spatially varying spectral index and temperature, we find that the statistical uncertainty on $r$ after foreground cleaning may be reduced by as much as 30--50 % by extending the upper limit of the frequency range from 400 to 600 GHz, with most of the improvement already gained at 500 GHz. We also note that a broader frequency range leads to better ability to discriminate between models through higher $\chi^2$ sensitivity. (abridged), Comment: 18 pages, 13 figures. Published in A&A

The bright polarized synchrotron emission, away from the Galactic plane, originates mostly from filamentary structures. We implement a filament finder algorithm which allows the detection of bright elongated structures in polarized intensity maps. We analyse the sky at 23 and 30 GHz as observed respectively by WMAP and Planck. We identify 19 filaments, 13 of which have been previously observed. For each filament, we study the polarization fraction, finding values typically larger than for the areas outside the filaments, excluding the Galactic plane, and a fraction of about 30% is reached in two filaments. We study the polarization spectral indices of the filaments, and find a spectral index consistent with the values found in previous analysis (about -3.1) for more diffuse regions. Decomposing the polarization signals into the $E$ and $B$ families, we find that most of the filaments are detected in $P_E$, but not in $P_B$. We then focus on understanding the statistical properties of the diffuse regions of the synchrotron emission at 23 GHz. Using Minkowski functionals and tensors, we analyse the non-Gaussianity and statistical isotropy of the polarized intensity maps. For a sky coverage corresponding to 80% of the fainter emission, and on scales smaller than 6 degrees ($\ell > 30$), the deviations from Gaussianity and isotropy are significantly higher than 3$\sigma$. The level of deviation decreases for smaller scales, however, it remains significantly high for the lowest analised scale ($\sim 1.5^\circ$). When 60% sky coverage is analysed, we find that the deviations never exceed 3$\sigma$. Finally, we present a simple data-driven model to generate non-Gaussian and anisotropic simulations of the synchrotron polarized emission. The simulations are fitted in order to match the spectral and statistical properties of the faintest 80% sky coverage of the data maps., Comment: 35 pages, 17 figures

We present the catalogue of Q-U-I JOint TEnerife (QUIJOTE) Wide Survey radio sources extracted from the maps of the Multi-Frequency Instrument compiled between 2012 and 2018. The catalogue contains 786 sources observed in intensity and polarization, and is divided into two separate sub-catalogues: one containing 47 bright sources previously studied by the \emph{Planck} collaboration and an extended catalogue of 739 sources either selected from the \emph{Planck} Second Catalogue of Compact Sources or found through a blind search carried out with a Mexican Hat 2 wavelet. A significant fraction of the sources in our catalogue (38.7 per cent) are within the $|b| \leq 20^\circ$ region of the Galactic plane. We determine statistical properties for those sources that are likely to be extragalactic. We find that these statistical properties are compatible with currently available models, with a $\sim$1.8 Jy completeness limit at 11 GHz. We provide the polarimetric properties of (38, 33, 31, 23) sources with P detected above the $99.99\%$ significance level at (11, 13, 17, 19) GHz, respectively. Median polarization fractions are in the $2.8$-$4.7$\% range in the 11-19 GHz frequency interval. We do not distinguish between Galactic and extragalactic sources here. The results presented here are consistent with those reported in the literature for flat- and steep-spectrum radio sources., Comment: 21 pages, 13 figures, 7 tables

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 WMAP 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 $\beta_s$ map of the Northern Celestial Hemisphere at a smoothing scale of $2^{\circ}$. 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, $c_s$, we find a value of $c_s = -0.0797 \pm 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 CMB, 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 $\beta_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., Comment: 25 pages, 29 figures, Accepted for publication in MNRAS

The QUIJOTE-MFI Northern Hemisphere Wide-Survey has provided maps of the sky above declinations $-30^\circ$ 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$^\circ$ 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 characterisation 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 $\pm$ 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 ($A_{\rm AME}/\tau_{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., Comment: 36 pages. Accepted for publication in MNRAS. Maps and derived data products available at https://research.iac.es/proyecto/quijote/

The Haze is an excess of microwave intensity emission surrounding the Galactic centre. It is spatially correlated with the $\gamma$-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 re-analysis 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 analyze 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$\,\leq\nu\leq\,$70 GHz is a power-law with spectral index $\beta^{\rm H}=-2.79\pm0.08$, which is flatter than the Galactic synchrotron in the same region ($\beta^{\rm S}=-2.98\pm0.04$), but steeper than that obtained from previous works ($\beta^{\rm H}\sim-2.5$ at 23$\,\leq\,\nu\leq\,$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 ($\beta \sim -3.2$ at 2.3 $\leq\nu\leq$ 23 GHz), and becomes flatter above 11 GHz., Comment: 31 pages. Accepted for publication in MNRAS

We present new intensity and polarisation 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 deg angular resolution, with a sensitivity in the range 35-79 ${\mu}$K/beam for total intensity and 13-23 ${\mu}$K/beam for polarisation. 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${\sigma}$ and W51 at 4.0${\sigma}$ with peak frequencies ${\nu}_{AME}$ = (20.0 $\pm$ 1.4) GHz and ${\nu}_{AME}$ = (17.7 $\pm$ 3.6) GHz respectively; this is the first detection of AME towards W51. The contamination from ultra-compact HII regions to the residual AME flux density is estimated at 10% in W49 and 5% in W51, and does not rule out the AME detection. The polarised SEDs reveal a synchrotron contribution with spectral indices ${\alpha}_s$ = -0.67 $\pm$ 0.10 in W49 and ${\alpha}_s$ = -0.51 $\pm$ 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}_{0,s}$ = (114 $\pm$ 73) GHz, in agreement with previous studies; our analysis, however, rules out any AME contribution which had been previously claimed towards IC443. No evidence of polarised AME emission is detected in this study., Comment: 29 pages, 11 figures. Accepted for publication in MNRAS

We present QUIJOTE intensity and polarization maps in four frequency bands centred around 11, 13, 17 and 19GHz, and covering approximately 29000 deg$^2$, including most of the Northern sky region. These maps result from 9000 h of observations taken between May 2013 and June 2018 with the first QUIJOTE instrument (MFI), and have angular resolutions of around $1^\circ$, and sensitivities in polarization within the range 35-40 $\mu$K per 1-degree beam, being a factor $\sim 2$-$4$ worse in intensity. We discuss the data processing pipeline employed, and the basic characteristics of the maps in terms of real space statistics and angular power spectra. A number of validation tests have been applied to characterise the accuracy of the calibration and the residual level of systematic effects, finding a conservative overall calibration uncertainty of 5%. We also discuss flux densities for four bright celestial sources (Tau A, Cas A, Cyg A and 3C274) which are often used as calibrators at microwave frequencies. The polarization signal in our maps is dominated by synchrotron emission. The distribution of spectral index values between the 11GHz and WMAP 23GHz map peaks at $\beta=-3.09$ with a standard deviation of 0.14. The measured BB/EE ratio at scales of $\ell=80$ is $0.26\pm 0.07$ for a Galactic cut $|b|>10^\circ$. We find a positive TE correlation for 11GHz at large angular scales ($\ell \lesssim 50$), while the EB and TB signals are consistent with zero in the multipole range $30 \lesssim \ell \lesssim 150$. The maps discussed in this paper are publicly available., Comment: 60 pages. Accepted for publication in MNRAS. Maps and derived data products available at https://research.iac.es/proyecto/quijote/

This paper presents a polarization calibration method applied to a microwave polarimeter demonstrator based on a near-infrared (NIR) frequency up-conversion stage that allows both optical correlation and signal detection at a wavelength of 1550 nm. The instrument was designed to measure the polarization of cosmic microwave background (CMB) radiation from the sky, obtaining the Stokes parameters of the incoming signal simultaneously, in a frequency range from 10 to 20 GHz. A linearly polarized input signal with a variable polarization angle is used as excitation in the polarimeter calibration setup mounted in the laboratory. The polarimeter systematic errors can be corrected with the proposed calibration procedure, achieving high levels of polarization efficiency (low polarization percentage errors) and low polarization angle errors. The calibration method is based on the fitting of polarization errors by means of sinusoidal functions composed of additive or multiplicative terms. The accuracy of the fitting increases with the number of terms in such a way that the typical error levels required in low-frequency CMB experiments can be achieved with only a few terms in the fitting functions. On the other hand, assuming that the calibration signal is known with the required accuracy, additional terms can be calculated to reach the error levels needed in ultrasensitive B-mode polarization CMB experiments., Comment: 17 pages, 7 figures, 6 tables

The polarization of the cosmic microwave background (CMB) can be used to search for parity-violating processes like that predicted by a Chern-Simons coupling to a light pseudoscalar field. Such an interaction rotates $E$ modes into $B$ modes in the observed CMB signal by an effect known as cosmic birefringence. Even though isotropic birefringence can be confused with the rotation produced by a miscalibration of the detectors' polarization angles the degeneracy between both effects is broken when Galactic foreground emission is used as a calibrator. In this work, we use realistic simulations of the High-Frequency Instrument of the Planck mission to test the impact that Galactic foreground emission and instrumental systematics have on the recent birefringence measurements obtained through this technique. Our results demonstrate the robustness of the methodology against the miscalibration of polarization angles and other systematic effects, like intensity-to-polarization leakage, beam leakage, or cross-polarization effects. However, our estimator is sensitive to the $EB$ correlation of polarized foreground emission. Here we propose to correct the bias induced by dust $EB$ by modeling the foreground signal with templates produced in Bayesian component-separation analyses that fit parametric models to CMB data. Acknowledging the limitations of currently available dust templates like that of the Commander sky model, high-precision CMB data and a characterization of dust beyond the modified blackbody paradigm are needed to obtain a definitive measurement of cosmic birefringence in the future., Comment: 36 pages, 19 figures, matches the version accepted for publication at JCAP

We propose a novel experiment, the Canfranc Axion Detection Experiment (CADEx), to probe dark matter axions with masses in the range 330-460 $\mu$eV, within the W-band (80-110 GHz), an unexplored parameter space in the well-motivated dark matter window of Quantum ChromoDynamics (QCD) axions. The experimental design consists of a microwave resonant cavity haloscope in a high static magnetic field coupled to a highly sensitive detecting system based on Kinetic Inductance Detectors via optimized quasi-optics (horns and mirrors). The experiment is in preparation and will be installed in the dilution refrigerator of the Canfranc Underground Laboratory. Sensitivity forecasts for axion detection with CADEx, together with the potential of the experiment to search for dark photons, are presented., Comment: 23 pages, 10 figures

The purpose of this work is to characterize the diffuse Galactic polarized synchrotron. We present EE, BB, and EB power spectra estimated cross-correlating Planck and WMAP polarization frequency maps at 23 and 30 GHz, for a set of six sky regions covering from 30% to 94% of the sky. The EE and BB angular power spectra show a steep decay of the spectral amplitude as a function of multipole, approximated by a power law with power indices around -2.9 for both components. The B/E ratio is about 0.22. The EB cross-component is compatible with zero at 1$\sigma$, with upper constraint on the EB/EE ratio of 1.2% at the 2$\sigma$ level. The recovered SED, in the frequency range 23-30 GHz, shows E and B power-law spectral indices compatible between themselves with a value of about -3., Comment: Contribution to the 2022 Cosmology session of the 56th Rencontres de Moriond

We search for the signature of parity-violating physics in the Cosmic Microwave Background using Planck polarization data from the Public Release 4 (PR4 or $\mathtt{NPIPE}$). For nearly full-sky data, we initially find a birefringence angle $\beta=0.30^\circ\pm0.11^\circ$ ($68\%$~C.L.). We also find that the values of $\beta$ decrease as we enlarge the Galactic mask, which can be interpreted as the effect of polarized foreground emission. We use two independent approaches to model this effect and mitigate its impact on $\beta$. Although results are promising, and the good agreement between both models is encouraging, we do not assign cosmological significance to the measured value of $\beta$ until we improve our knowledge of the foreground polarization. Acknowledging that the miscalibration of polarization angles is not the only instrumental systematic that can create spurious TB and EB correlations, we also perform a detailed study of $\mathtt{NPIPE}$ end-to-end simulations to prove that our measurements of $\beta$ are not significantly affected by any of the known systematics., Comment: 8 pages, 2 figures, contribution to the 2022 Cosmology session of the 56th Rencontres de Moriond

LiteBIRD, the Lite (Light) satellite for the study of B-mode polarization and Inflation from cosmic background Radiation Detection, is a space mission for primordial cosmology and fundamental physics. The Japan Aerospace Exploration Agency (JAXA) selected LiteBIRD in May 2019 as a strategic large-class (L-class) mission, with an expected launch in the late 2020s using JAXA's H3 rocket. LiteBIRD is planned to orbit the Sun-Earth Lagrangian point L2, where it will map the cosmic microwave background (CMB) polarization over the entire sky for three years, with three telescopes in 15 frequency bands between 34 and 448 GHz, to achieve an unprecedented total sensitivity of 2.2$\mu$K-arcmin, with a typical angular resolution of 0.5$^\circ$ at 100 GHz. The primary scientific objective of LiteBIRD is to search for the signal from cosmic inflation, either making a discovery or ruling out well-motivated inflationary models. The measurements of LiteBIRD will also provide us with insight into the quantum nature of gravity and other new physics beyond the standard models of particle physics and cosmology. We provide an overview of the LiteBIRD project, including scientific objectives, mission and system requirements, operation concept, spacecraft and payload module design, expected scientific outcomes, potential design extensions and synergies with other projects., Comment: 155 pages, accepted for publication in PTEP

A methodology to provide the polarization angle requirements for different sets of detectors, at a given frequency of a CMB polarization experiment, is presented. The uncertainties in the polarization angle of each detector set are related to a given bias on the tensor-to-scalar ratio $r$ parameter. The approach is grounded in using a linear combination of the detector sets to obtain the CMB polarization signal. In addition, assuming that the uncertainties on the polarization angle are in the small angle limit (lower than a few degrees), it is possible to derive analytic expressions to establish the requirements. The methodology also accounts for possible correlations among detectors, that may originate from the optics, wafers, etc. The approach is applied to the LiteBIRD space mission. We show that, for the most restrictive case (i.e., full correlation of the polarization angle systematics among detector sets), the requirements on the polarization angle uncertainties are of around 1 arcmin at the most sensitive frequency bands (i.e., $\approx 150$ GHz) and of few tens of arcmin at the lowest (i.e., $\approx 40$ GHz) and highest (i.e., $\approx 400$ GHz) observational bands. Conversely, for the least restrictive case (i.e., no correlation of the polarization angle systematics among detector sets), the requirements are $\approx 5$ times less restrictive than for the previous scenario. At the global and the telescope levels, polarization angle knowledge of a few arcmins is sufficient for correlated global systematic errors and can be relaxed by a factor of two for fully uncorrelated errors in detector polarization angle. The reported uncertainty levels are needed in order to have the bias on $r$ due to systematics below the limit established by the LiteBIRD collaboration., Comment: 26 pages, 6 figures. Minor changes to match the final version in JCAP

We search for the signature of parity-violating physics in the cosmic microwave background, called cosmic birefringence, using the Planck data release 4. We initially find a birefringence angle of $\beta=0.30\pm0.11$ (68% C.L.) for nearly full-sky data. The values of $\beta$ decrease as we enlarge the Galactic mask, which can be interpreted as the effect of polarized foreground emission. Two independent ways to model this effect are used to mitigate the systematic impact on $\beta$ for different sky fractions. We choose not to assign cosmological significance to the measured value of $\beta$ until we improve our knowledge of the foreground polarization., Comment: 7 pages, 1 figure, 1 table. (v2) References added. Accepted for publication in Physical Review Letters

We present constraints on the tensor-to-scalar ratio r using a combination of BICEP/Keck 2018 and Planck PR4 data allowing us to fit for r consistently with the six parameters of the $\Lambda$CDM model without fixing any of them. In particular, we are able to derive a constraint on the reionization optical depth $\tau$ and thus propagate its uncertainty onto the posterior distribution for r. While Planck sensitivity to r is no longer comparable with ground-based measurements, combining Planck with BK18 and BAO gives results consistent with r = 0 and tightens the constraint to r < 0.032., Comment: 7 pages, 7 figures, published in PRD

The new generation of CMB polarization experiments will reach limits of sensitivity never achieved before to detect the primordial B-mode signal. However, all these efforts will be futile if we lack a tight control of systematics. Here, we focus on the systematic that arises from the uncertainty on the calibration of polarization angles. Miscalibrated polarization angles induce a mixing of E- and B-modes that obscures the primordial B-mode signal. We introduce an iterative angular power spectra maximum likelihood-based method to calculate polarization angles from the multi-frequency signal. The basis behind this methodology grounds on nulling the EB power spectra. To simplify the likelihood, we assume that the rotation angles are small (<6 deg) and, the maximum likelihood solution for the rotation angles is obtained by applying an iterative process where the covariance matrix does not depend on the angles per iteration, i.e., the rotation angles are fixed to the estimated angles in the previous iteration. With these assumptions, we obtain an analytical linear system which leads to a very fast computational implementation. We show that with this methodology we are able to determine the rotation angle for each frequency with sufficiently good accuracy. To prove the latter point we perform component separation analyses using the parametric component separation method B-SeCRET with two different approaches. In the first approach we apply the B-SeCRET pipeline to the signal de-rotated with the estimation of the angles, while in the second, the rotation angles are treated as model parameters using the estimation of the angles as a prior information. We obtain that the rotation angles estimations improve after applying the second approach, and show that the systematic residuals due to the non-null calibration polarization angles are mitigated to the order of a 1% at the power spectrum level., Comment: 28 pages, 12 figures, submitted to JCAP

The purpose of this work is to characterize the diffuse Galactic polarized synchrotron, which is the dominant CMB foreground emission at low frequency. We present EE, BB, and EB power spectra estimated from polarization frequency maps at 23 and 30 GHz as observed respectively by the WMAP K-band and the Planck lowest frequency channel, for a set of six sky regions covering from 30% to 94% of the sky. We study the synchrotron polarization angular distribution and spectral energy distribution (SED) by means of the so-called pseudo-$C_\ell$ formalism, provided by the NaMaster package, in the multipole interval 30 $\leq$ $\ell$ $\leq$ 300. Best results are obtained cross-correlating Planck and WMAP data. The EE and BB angular power spectra show a steep decay of the spectral amplitude as a function of multipole, approximated by a power law $C^{EE,BB} \propto \ell^{\alpha_{EE,BB}}$, with $\alpha_{EE}= -2.95\pm0.04$ and $\alpha_{BB}=-2.85\pm0.14$. The B/E power asymmetry is proved with a B-to-E ratio, computed as the amplitude ratio at the pivot multipole $\ell = 80$, of 0.22$\pm$0.02. The EB cross-component is compatible with zero at 1$\sigma$, with an upper constraint on the EB/EE ratio of 1.2% at the 2$\sigma$ level. We show that the EE and BB power-law model with null EB cross-correlation describes reasonably well the diffuse synchrotron polarization emission for the full sky if the bright Galactic center and point sources are masked. The recovered SED shows power-law spectral indices $\beta_{EE}= -3.00 \pm 0.10$ and $\beta_{BB} = -3.05\pm0.36$ compatible between themselves, in the frequency range 23-30 GHz. Results also seem to indicate that the SED gets steeper from low to high Galactic latitude.

Map-making is an important step for the data analysis of Cosmic Microwave Background (CMB) experiments. It consists of converting the data, which are typically a long, complex and noisy collection of measurements, into a map, which is an image of the observed sky. We present in this paper a new map-making code named PICASSO (Polarization and Intensity CArtographer for Scanned Sky Observations), which was implemented to construct intensity and polarization maps from the Multi Frequency Instrument (MFI) of the QUIJOTE (Q-U-I Joint TEnerife) CMB polarization experiment. PICASSO is based on the destriping algorithm, and is suited to address specific issues of ground-based microwave observations, with a technique that allows the fit of a template function in the time domain, during the map-making step. This paper describes the PICASSO code, validating it with simulations and assessing its performance. For this purpose, we produced realistic simulations of the QUIJOTE-MFI survey of the northern sky (approximately $\sim 20,000\,$deg$^2$), and analysed the reconstructed maps with PICASSO, using real and harmonic space statistics. We show that, for this sky area, PICASSO is able to reconstruct, with high fidelity, the injected signal, recovering all the scales with $\ell>10$ in TT, EE and BB. The signal error is better than 0.001% at $20<\ell<200$. Finally, we validated some of the methods that will be applied to the real wide-survey data, like the detection of the CMB anisotropies via cross-correlation analyses. Despite that the implementation of PICASSO is specific for QUIJOTE-MFI data, it could be adapted to other experiments., Comment: 20 pages, 12 figures, accepted for publication in MNRAS

We present ECLIPSE (Efficient Cmb poLarization and Intensity Power Spectra Estimator), an optimized implementation of the Quadratic Maximum Likelihood (QML) method for the estimation of the power spectra of the Cosmic Microwave Background (CMB). This approach allows one to reduce significantly the computational costs associated to this technique, allowing to estimate the power spectra up to higher multipoles than previous implementations. In particular, for a resolution of $N_\mathrm{side}=64$, $\ell_{\mathrm{max}}=192$ and a typical Galactic mask, the number of operations can be reduced by approximately a factor of 1000 in a full analysis including intensity and polarization with respect to an efficient direct implementation of the method. In addition, if one is interested in studying only polarization, it is possible to obtain the power spectra of the E and B modes with a further reduction of computational resources without degrading the results. We also show that for experiments observing a small fraction of the sky, the Fisher matrix becomes singular and, in this case, the standard QML can not be applied. To solve this problem, we have developed a binned version of the method that is unbiased and of minimum variance. We also test the robustness of the QML estimator when the assumed fiducial model differs from that of the sky and show the performance of an iterative approach. Finally, we present a comparison of the results obtained by QML and a pseudo-$C_{\ell}$ estimator (NaMaster) for a next-generation satellite, showing that, as expected, QML produces significantly smaller errors at low multipoles. The ECLIPSE fast QML code developed in this work will be made publicly available., Comment: 51 pages, 17 figures, 6 tables

LiteBIRD is a JAXA-led Strategic Large-Class mission designed to search for the existence of the primordial gravitational waves produced during the inflationary phase of the Universe, through the measurements of their imprint onto the polarization of the cosmic microwave background (CMB). These measurements, requiring unprecedented sensitivity, will be performed over the full sky, at large angular scales, and over 15 frequency bands from 34GHz to 448GHz. The LiteBIRD instruments consist of three telescopes, namely the Low-, Medium- and High-Frequency Telescope (respectively LFT, MFT and HFT). We present in this paper an overview of the design of the Medium-Frequency Telescope (89-224GHz) and the High-Frequency Telescope (166-448GHz), the so-called MHFT, under European responsibility, which are two cryogenic refractive telescopes cooled down to 5K. They include a continuous rotating half-wave plate as the first optical element, two high-density polyethylene (HDPE) lenses and more than three thousand transition-edge sensor (TES) detectors cooled to 100mK. We provide an overview of the concept design and the remaining specific challenges that we have to face in order to achieve the scientific goals of LiteBIRD., Comment: SPIE Conference

LiteBIRD, the Lite (Light) satellite for the study of B-mode polarization and Inflation from cosmic background Radiation Detection, is a space mission for primordial cosmology and fundamental physics. JAXA selected LiteBIRD in May 2019 as a strategic large-class (L-class) mission, with its expected launch in the late 2020s using JAXA's H3 rocket. LiteBIRD plans to map the cosmic microwave background (CMB) polarization over the full sky with unprecedented precision. Its main scientific objective is to carry out a definitive search for the signal from cosmic inflation, either making a discovery or ruling out well-motivated inflationary models. The measurements of LiteBIRD will also provide us with an insight into the quantum nature of gravity and other new physics beyond the standard models of particle physics and cosmology. To this end, LiteBIRD will perform full-sky surveys for three years at the Sun-Earth Lagrangian point L2 for 15 frequency bands between 34 and 448 GHz with three telescopes, to achieve a total sensitivity of 2.16 micro K-arcmin with a typical angular resolution of 0.5 deg. at 100GHz. We provide an overview of the LiteBIRD project, including scientific objectives, mission requirements, top-level system requirements, operation concept, and expected scientific outcomes., Comment: 20 pages, 9 figures

LiteBIRD has been selected as JAXA's strategic large mission in the 2020s, to observe the cosmic microwave background (CMB) $B$-mode polarization over the full sky at large angular scales. The challenges of LiteBIRD are the wide field-of-view (FoV) and broadband capabilities of millimeter-wave polarization measurements, which are derived from the system requirements. The possible paths of stray light increase with a wider FoV and the far sidelobe knowledge of $-56$ dB is a challenging optical requirement. A crossed-Dragone configuration was chosen for the low frequency telescope (LFT : 34--161 GHz), one of LiteBIRD's onboard telescopes. It has a wide field-of-view ($18^\circ \times 9^\circ$) with an aperture of 400 mm in diameter, corresponding to an angular resolution of about 30 arcminutes around 100 GHz. The focal ratio f/3.0 and the crossing angle of the optical axes of 90$^\circ$ are chosen after an extensive study of the stray light. The primary and secondary reflectors have rectangular shapes with serrations to reduce the diffraction pattern from the edges of the mirrors. The reflectors and structure are made of aluminum to proportionally contract from warm down to the operating temperature at $5\,$K. A 1/4 scaled model of the LFT has been developed to validate the wide field-of-view design and to demonstrate the reduced far sidelobes. A polarization modulation unit (PMU), realized with a half-wave plate (HWP) is placed in front of the aperture stop, the entrance pupil of this system. A large focal plane with approximately 1000 AlMn TES detectors and frequency multiplexing SQUID amplifiers is cooled to 100 mK. The lens and sinuous antennas have broadband capability. Performance specifications of the LFT and an outline of the proposed verification plan are presented., Comment: 21 pages, 14 figures

We present constraints on the tensor-to-scalar ratio r using Planck data. We use the latest release of Planck maps (PR4), processed with the NPIPE code, which produces calibrated frequency maps in temperature and polarization for all Planck channels from 30 GHz to 857 GHz using the same pipeline. We computed constraints on r using the BB angular power spectrum, and we also discuss constraints coming from the TT spectrum. Given Planck's noise level, the TT spectrum gives constraints on r that are cosmic-variance limited (with $\sigma$(r)=0.093), but we show that the marginalized posterior peaks towards negative values of r at about the 1.2$\sigma$ level. We derived Planck constraints using the BB power spectrum at both large angular scales (the 'reionization bump') and intermediate angular scales (the 'recombination bump') from $\ell$=2 to 150, and find a stronger constraint than that from TT, with $\sigma$(r)=0.069. The Planck BB spectrum shows no systematic bias, and is compatible with zero, given both the statistical noise and the systematic uncertainties. The likelihood analysis using B modes yields the constraint r<0.158 at 95% confidence using more than 50% of the sky. This upper limit tightens to r<0.069 when Planck EE, BB, and EB power spectra are combined consistently, and it tightens further to r<0.056 when the Planck TT power spectrum is included in the combination. Finally, combining Planck with BICEP2/Keck 2015 data yields an upper limit of r<0.044., Comment: 11 pages, 10 figures, 8 appendix pages. Revised to match version accepted by Astronomy & Astrophysics

We describe an extension of the most recent version of the Planck Catalogue of Compact Sources (PCCS2), produced using a new multi-band Bayesian Extraction and Estimation Package (BeeP). BeeP assumes that the compact sources present in PCCS2 at 857 GHz have a dust-like spectral energy distribution, which leads to emission at both lower and higher frequencies, and adjusts the parameters of the source and its SED to fit the emission observed in Planck's three highest frequency channels at 353, 545, and 857 GHz, as well as the IRIS map at 3000 GHz. In order to reduce confusion regarding diffuse cirrus emission, BeeP's data model includes a description of the background emission surrounding each source, and it adjusts the confidence in the source parameter extraction based on the statistical properties of the spatial distribution of the background emission. BeeP produces the following three new sets of parameters for each source: (a) fits to a modified blackbody (MBB) thermal emission model of the source; (b) SED-independent source flux densities at each frequency considered; and (c) fits to an MBB model of the background in which the source is embedded. BeeP also calculates, for each source, a reliability parameter, which takes into account confusion due to the surrounding cirrus. We define a high-reliability subset (BeeP/base), containing 26 083 sources (54.1 per cent of the total PCCS2 catalogue), the majority of which have no information on reliability in the PCCS2. The results of the BeeP extension of PCCS2, which are made publicly available via the PLA, will enable the study of the thermal properties of well-defined samples of compact Galactic and extra-galactic dusty sources., Comment: 55 pages. Accepted for publication in A&A. The BeeP catalogue will be published in the Planck Legacy Archive (https://pla.esac.esa.int/pla)

We present the NPIPE processing pipeline, which produces calibrated frequency maps in temperature and polarization from data from the Planck Low Frequency Instrument (LFI) and High Frequency Instrument (HFI) using high-performance computers. NPIPE represents a natural evolution of previous Planck analysis efforts, and combines some of the most powerful features of the separate LFI and HFI analysis pipelines. The net effect of the improvements is lower levels of noise and systematics in both frequency and component maps at essentially all angular scales, as well as notably improved internal consistency between the various frequency channels. Based on the NPIPE maps, we present the first estimate of the Solar dipole determined through component separation across all nine Planck frequencies. The amplitude is ($3366.6 \pm 2.7$)$\mu$K, consistent with, albeit slightly higher than, earlier estimates. From the large-scale polarization data, we derive an updated estimate of the optical depth of reionization of $\tau = 0.051 \pm 0.006$, which appears robust with respect to data and sky cuts. There are 600 complete signal, noise and systematics simulations of the full-frequency and detector-set maps. As a Planck first, these simulations include full time-domain processing of the beam-convolved CMB anisotropies. The release of NPIPE maps and simulations is accompanied with a complete suite of raw and processed time-ordered data and the software, scripts, auxiliary data, and parameter files needed to improve further on the analysis and to run matching simulations., Comment: 97 pages, 93 figures and 16 tables, abstract abridged for arXiv submission, accepted for publication in A&A

Most of the CMB experiments proposed for the next generation aim to detect the Primordial Gravitational Wave Background (PGWB). The fulfillment of this objective depends on our capacity to separate Galactic foreground emissions and to \emph{delens} the secondary B-mode component induced by weak gravitational lensing. Focusing on the latter of these efforts, in this work we briefly review the basic aspects of lensing, and exhaustively compare the performance of current delensing methodologies and implementations within the Born approximation as a preparation for the analysis of the data to come in the following years. Two of the main conclusions that can be drawn from our study are that, for next-generation experiments, delensing efficiency will still be limited by the quality of the data itself rather than by the limitations of current delensing methodologies, and that template delensing within the antilensing approximation will be the optimal (balancing accuracy and computational cost) technique to employ. We then evaluate the delensing capabilities of future experiments (like the Simons Observatory, the CMB Stage-IV, or the LiteBIRD and PICO satellites) by applying that methodology onto numerical simulations of the typical CMB and lensing potential reconstructions that they are expected to produce, and quantify how internal and external delensing will help them to improve their sensitivity to detect the PGWB. We also consider the benefits that a joint analysis of their data would provide., Comment: 37 pages, 12 figures, submitted to JCAP

The largest temperature anisotropy in the cosmic microwave background (CMB) is the dipole, which has been measured with increasing accuracy for more than three decades, particularly with the Planck satellite. The simplest interpretation of the dipole is that it is due to our motion with respect to the rest frame of the CMB. Since current CMB experiments infer temperature anisotropies from angular intensity variations, the dipole modulates the temperature anisotropies with the same frequency dependence as the thermal Sunyaev-Zeldovich (tSZ) effect. We present the first, and significant, detection of this signal in the tSZ maps and find that it is consistent with direct measurements of the CMB dipole, as expected. The signal contributes power in the tSZ maps, which is modulated in a quadrupolar pattern, and we estimate its contribution to the tSZ bispectrum, noting that it contributes negligible noise to the bispectrum at relevant scales., Comment: 15 pages, 8 figures. Added references, small clarifying and language edits. All results remain the same