Galactic Foregrounds Research Articles

Detecting a gravitational wave background from inflation with null energy condition violation: prospects for Taiji

The null energy condition (NEC) is a fundamental principle in general relativity, and its violation could leave discernible signatures in gravitational waves (GWs). A violation of the NEC during the primordial era would imprint a blue-tilted spectrum on the stochastic gravitational wave background (SGWB) at nanohertz frequencies, potentially accounting for the recently detected signal by pulsar timing arrays. Remarkably, models of NEC violation during inflation also predict a nearly scale-invariant GW spectrum in the millihertz frequency range, which could be detectable by upcoming space-based GW detectors such as Taiji. The observation of this distinctive spectrum would provide compelling evidence for new physics beyond the standard cosmological paradigm. In this study, we explore Taiji’s ability to detect an SGWB arising from NEC violation during inflation, considering various foregrounds and noise sources, including an extragalactic foreground from binary black hole mergers throughout the Universe, a galactic foreground from white dwarf binaries, and the intrinsic noise of the Taiji detector. Employing comprehensive Bayesian parameter estimation techniques to analyze simulated Taiji data, we demonstrate a remarkable improvement in precision of three orders of magnitude compared to the NANOGrav 15-year data set for measuring the tensor power spectrum amplitude, PT,2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$P_{T,2}$$\\end{document}, during the second inflationary stage. This substantial enhancement in measurement capabilities underscores Taiji’s potential as a powerful probe for investigating the NEC violation in the early Universe.

Gravitational waves from inflation in LISA: reconstruction pipeline and physics interpretation

Various scenarios of cosmic inflation enhance the amplitude of the stochastic gravitational wave background (SGWB) at frequencies detectable by the LISA detector. We develop tools for a template-based analysis of the SGWB and introduce a template databank to describe well-motivated signals from inflation, prototype their template-based searches, and forecast their reconstruction with LISA. Specifically, we classify seven templates based on their signal frequency shape, and we identifyrepresentative fundamental physics models leading to them. By running a template-based analysis, we forecast the accuracy with which LISA can reconstruct the template parameters of representative benchmark signals, with and without galactic and extragalactic foregrounds. We identify the parameter regions that can be probed by LISA within each template. Finally, we investigate how our signal reconstructions shed light on fundamental physics models of inflation: we discuss their impact for measurements of e.g., the couplings of inflationary axions to gauge fields; the graviton mass during inflation; the fluctuation seeds of primordialblack holes; the consequences of excited states during inflation, and the presence of small-scale spectral features.

Gravitational waves from first-order phase transitions in LISA: reconstruction pipeline and physics interpretation

We develop a tool for the analysis of stochastic gravitational wave backgrounds from cosmological first-order phase transitions with LISA: we initiate a template databank for these signals, prototype their searches, and forecast their reconstruction.The templates encompass the gravitational wave signals sourced by bubble collisions, sound waves and turbulence. Accounting for Galactic and extra-Galactic foregrounds, we forecast the region of the parameter space that LISA will reconstruct with better than ∼ 10% accuracy, if certain experimental and theoretical uncertainties are solved by the time LISA flies.We illustrate the accuracy with which LISA can reconstruct the parameters on a few benchmark signals, both in terms of the template parameters and the phase transition ones.To show the impact of the forecasts on physics beyond the Standard Model, we map the reconstructed benchmark measurements into the parameter spaces of the singlet extension of the Standard Model and of the classically conformal invariant U(1) B-L model.

A power spectrum approach to the search for axion-like particles from resolved galaxy clusters using CMB as a backlight

Axions or axion-like particles (ALPs) are hypothetical particles predicted by beyond standard model theories, which make one of the dark matter candidates. These particles can convert into photons and vice-versa in the presence of a magnetic field, with a probability decided by its coupling strength gaγ . One of the ways to detect these particles is by using the Cosmic Microwave Background (CMB) as a backlight. As the CMB photons pass through a galaxy cluster, they can get converted into ALPs in the mass range 10-15 eV to 10-11 eV through resonant conversion in the presence of cluster magnetic fields. This leads to a polarized spectral distortion (α-distortion) in the CMB as the photon polarization parallel to the magnetic field in the galaxy cluster is involved in the conversion. The fluctuations in the magnetic field and electron density in a galaxy cluster lead to spatially varying α-distortion around the cluster, with a power spectrum that is different from the lensed E-mode and B-mode CMB polarization power spectrum for the standard model of cosmology. By measuring the difference in the polarization power spectrum around a galaxy cluster from the all-sky signal, one can find new α-distortion in the sky. For the resolved galaxy clusters, if the redshift, electron density, and magnetic field profiles of the cluster can be constrained using optical, X-ray, and radio observations, one can measure the coupling strength gaγ from the ALP power spectrum. The contamination from CMB and galactic foregrounds such as synchrotron and dust can be mitigated by using multiple frequency bands by leveraging on the difference in the spectral shape of the signal from foregrounds. Using the ILC technique to clean the foregrounds, we show that the new power spectrum-based approach of the resolved galaxy clusters from upcoming CMB experiments such as Simons Observatory and CMB-S4 can detect (or put constraints) on the ALP-photon coupling strength of gaγ < 5.2 × 10-12 GeV-1 and gaγ < 3.6 × 10-12 GeV-1 at 95% C.I. respectively for ALPs of masses 10-13 eV or for smaller gaγ for lighter ALP masses.

The Simons Observatory: Combining cross-spectral foreground cleaning with multitracer B -mode delensing for improved constraints on inflation

The Simons Observatory (SO), due to start full science operations in early 2025, aims to set tight constraints on inflationary physics by inferring the tensor-to-scalar ratio r from measurements of cosmic microwave background (CMB) polarization B-modes. Its nominal design including three small-aperture telescopes (SATs) targets a precision σ(r=0)≤0.003 without delensing. Achieving this goal and further reducing uncertainties requires a thorough understanding and mitigation of other large-scale B-mode sources such as Galactic foregrounds and weak gravitational lensing. We present an analysis pipeline aiming to estimate r by including delensing within a cross-spectral likelihood, and demonstrate it for the first time on SO-like simulations accounting for various levels of foreground complexity, inhomogeneous noise and partial sky coverage. As introduced in an earlier SO delensing paper, lensing B-modes are synthesized using internal CMB lensing reconstructions as well as -like cosmic infrared background maps and LSST-like galaxy density maps. We then extend SO’s power-spectrum-based foreground-cleaning algorithm to include all auto- and cross-spectra between the lensing template and the SAT B-modes in the likelihood function. This allows us to constrain r and the parameters of our foreground model simultaneously. Within this framework, we demonstrate the equivalence of map-based and cross-spectral delensing and use it to motivate an optimized pixel-weighting scheme for power spectrum estimation. We start by validating our pipeline in the simplistic case of uniform foreground spectral energy distributions. In the absence of primordial B-modes, we find that the 1σ statistical uncertainty on r, σ(r), decreases by 37% as a result of delensing. Tensor modes at the level of r=0.01 are successfully detected by our pipeline. Even when using more realistic foreground models including spatial variations in the dust and synchrotron spectral properties, we obtain unbiased estimates of r both with and without delensing by employing the moment-expansion method. In this case, uncertainties are increased due to the higher number of model parameters, and delensing-related improvements range between 27% and 31%. These results constitute the first realistic assessment of the delensing performance at SO’s nominal sensitivity level. Published by the American Physical Society 2024

LiteBIRD science goals and forecasts: a full-sky measurement of gravitational lensing of the CMB

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 μ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 80% 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.

Impact of beam far side-lobe knowledge in the presence of foregrounds for LiteBIRD

We present a study of the impact of a beam far side-lobe lack of knowledge on the measurement of the Cosmic Microwave Background B-mode signal at large scale. Beam far side-lobes induce a mismatch in the transfer function of Galactic foregrounds between the dipole and higher multipoles which degrads the performances of component separation methods. This leads to foreground residuals in the CMB map. It is expected to be one of the main source of systematic effects in future CMB polarization observations. Thus, it becomes crucial for all-sky survey missions to take into account the interplays between beam systematic effects and all the data analysis steps. LiteBIRD is the ISAS/JAXA second strategic large-class satellite mission and is dedicated to target the measurement of CMB primordial B modes by reaching a sensitivity on the tensor-to-scalar ratio r of σ(r) ≤ 10-3 assuming r = 0. The primary goal of this paper is to provide the methodology and develop the framework to carry out the end-to-end study of beam far side-lobe effects for a space-borne CMB experiment. We introduce uncertainties in the beam model, and propagate the beam effects through all the steps of the analysis pipeline, most importantly including component separation, up to the cosmological results in the form of a bias δr. As a demonstration of our framework, we derive requirements on the calibration and modeling for the LiteBIRD's beams under given assumptions on design, simulation, component separation method and allocated error budget. In particular, we assume a parametric method of component separation with no mitigation of the far side-lobes effect at any stage of the analysis pipeline. We show that δ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 δ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 the level of ∼ 10-4, to achieve the required limit on the bias δr < 1.9 × 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, such as PICO or CMB-Bharat. We further discuss the limitations of this framework and potential extensions to circumvent them.

LiteBIRD science goals and forecasts: improving sensitivity to inflationary gravitational waves with multitracer delensing

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.

The Simons Observatory: Pipeline comparison and validation for large-scale B-modes

Context. The upcoming Simons Observatory Small Aperture Telescopes aim at achieving a constraint on the primordial tensor-to-scalar ratio r at the level of σ(r = 0)≲0.003, observing the polarized CMB in the presence of partial sky coverage, cosmic variance, inhomogeneous non-white noise, and Galactic foregrounds. Aims. We present three different analysis pipelines able to constrain r given the latest available instrument performance, and compare their predictions on a set of sky simulations that allow us to explore a number of Galactic foreground models and elements of instrumental noise, relevant for the Simons Observatory. Methods. The three pipelines employ different combinations of parametric and non-parametric component separation at the map and power spectrum levels, and use B-mode purification to estimate the CMB B-mode power spectrum. We applied them to a common set of simulated realistic frequency maps, and compared and validated them with focus on their ability to extract robust constraints on the tensor-to-scalar ratio r. We evaluated their performance in terms of bias and statistical uncertainty on this parameter. Results. In most of the scenarios the three methodologies achieve similar performance. Nevertheless, several simulations with complex foreground signals lead to a &gt; 2σ bias on r if analyzed with the default versions of these pipelines, highlighting the need for more sophisticated pipeline components that marginalize over foreground residuals. We show two such extensions, using power-spectrum-based and map-based methods, that are able to fully reduce the bias on r below the statistical uncertainties in all foreground models explored, at a moderate cost in terms of σ(r).

A stochastic gravitational wave background in LISA from unresolved white dwarf binaries in the Large Magellanic Cloud

ABSTRACT The Laser Interferometer Space Antenna (LISA) is expected to detect a wide variety of gravitational wave sources in the mHz band. Some of these signals will elude individual detection, instead contributing as confusion noise to one of several stochastic gravitational-wave backgrounds (SGWBs) – notably including the ‘Galactic foreground’, a loud signal resulting from the superposition of millions of unresolved double white dwarf binaries (DWDs) in the Milky Way. It is possible that similar, weaker SGWBs will be detectable from other DWD populations in the local Universe, including the Large Magellanic Cloud (LMC). We use the Bayesian LISA Inference Package (blip) to investigate the possibility of an anisotropic SGWB generated by unresolved DWDs in the LMC. To do so, we compute the LMC SGWB from a realistic DWD population generated via binary population synthesis, simulate 4 years of time-domain data with blip comprised of stochastic contributions from the LMC SGWB and the LISA detector noise, and analyse this data with blip’s spherical harmonic anisotropic SGWB search. We also consider the case of spectral separation from the Galactic foreground. We present the results of these analyses and show, for the first time, that the unresolved DWDs in the LMC will comprise a significant SGWB for LISA.

Foreground removal with ILC methods for AliCPT-1

One of the main goals of most future CMB experiments is the precise measurement of CMB B-mode polarization, whose major obstacle is the Galactic foregrounds. In this paper, we evaluate the foreground cleaning performance of the variants of the ILC method on partial sky B-modes and analyze the main sources of biases on the BB power spectrum. Specially, we compare the NILC, the cILC (in three domains) and the cMILC methods for AliCPT-1 simulations. We find that the cILC methods implemented in harmonic space and needlet space are both competent to clean different models of foregrounds, which bias the tensor-to-scalar ratio about 0.008 at maximum, and constrain the tensor-to-scalar ratio to r < 0.043(95%CL) for the AliCPT-1 configuration. We also note that the deviation of the estimated noise bias from the actual one for ILC, dubbed the noise bias error (NBE) in this paper, might make significant effects on the power spectrum for a small footprint and low signal-to-noise ratio CMB experiment. We finally obtain its relation with respect to the noise residual which fits well with the simulated results.

Mapping the intracluster medium in the era of high-resolution X-ray spectroscopy

ABSTRACT High-resolution spectroscopy in soft X-rays will open a new window to map multiphase gas in galaxy clusters and probe physics of the intracluster medium (ICM), including chemical enrichment histories, circulation of matter and energy during large-scale structure evolution, stellar and black hole feedback, halo virialization, and gas mixing processes. An eV-level spectral resolution, large field of view, and effective area are essential to separate cluster emissions from the Galactic foreground and efficiently map the cluster outskirts. Several mission concepts that meet these criteria have been proposed recently, e.g. LEM, HUBS, and Super DIOS. This theoretical study explores what information on ICM physics could be recovered with such missions and the associated challenges. We emphasize the need for a comprehensive comparison between simulations and observations to interpret the high-resolution spectroscopic observations correctly. Using Line Emission Mapper (LEM) characteristics as an example, we demonstrate that it enables the use of soft X-ray emission lines (e.g. O vii/viii and Fe-L complex) from the cluster outskirts to measure the thermodynamic, chemical, and kinematic properties of the gas up to r200 and beyond. By generating mock observations with full backgrounds, analysing their images/spectra with observational approaches, and comparing the recovered characteristics with true ones from simulations, we develop six key science drivers for future missions, including the exploration of multiphase gas in galaxy clusters (e.g. temperature fluctuations, phase-space distributions), metallicity, ICM gas bulk motions and turbulence power spectra, ICM-cosmic filament interactions, and advances for cluster cosmology.

The European Low Frequency Survey

In this paper we present the European Low Frequency Survey (ELFS), a project that will enable the detection of primordial B-mode polarization by measuring the Galactic and extra-Galactic foregrounds in the 5– 120 GHz frequency window. Indeed, the main difficulty in measuring the Bmode 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.5m SA telescopes in Atacama, Chile (“ELFS on SA”). The receiver will be equipped with a fully digital backend based on the latest Xilinx RF System-on-Chip devices that will provide frequency resolution of 1MHz 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.

Measuring the CMB primordial B-modes with Bolometric Interferometry

The Q&amp;U Bolometric Interferometer for Cosmology (QL’BIC) is the first bolometric interferometer designed to measure the primordial B-mode polarization of the Cosmic Microwave Background (CMB). Bolometric interferometry is a novel technique that combines the sensitivity of bolometric detectors with the control of systematic effects that is typical of interferometry, both key features in the quest for the faint signal of the primordial B-modes. A unique feature is the so-called “spectral imaging”, i.e., the ability to recover the sky signal in several sub-bands within the physical band during data analysis. This feature provides an in-band spectral resolution of ∆v/v ~ 0.04 that is unattainable by a traditional imager. This is a key tool for controlling the Galactic foregrounds contamination. In this paper, we describe the principles of bolometric interferometry, the current status of the QU BIC experiment and future prospects.

Fitting and Comparing Galactic Foreground Models for Unbiased 21 cm Cosmology

Accurate detection of the cosmological 21 cm global signal requires galactic foreground models that can remove power over 106. Although foreground and global signal models unavoidably exhibit overlap in their vector spaces inducing bias error in the extracted signal, a second source of bias and error arises from inadequate foreground models, i.e., models that cannot fit spectra down to the noise level of the signal. We therefore test the level to which seven commonly employed foreground models—including nonlinear and linear forward models, polynomials, and maximally smooth polynomials—fit realistic simulated mock foreground spectra, as well as their dependence upon model inputs. The mock spectra are synthesized for an EDGES-like experiment and we compare all models’ goodness of fit and preference using a Kolmogorov–Smirnov (K-S) test of the noise-normalized residuals in order to compare models with differing, and sometimes indeterminable, degrees of freedom. For a single local sidereal time (LST) bin spectrum and p-value threshold of p = 0.05, the nonlinear forward model with four parameters is preferred (p = 0.99), while the linear forward model fits well with six to seven parameters (p = 0.94, 0.97, respectively). The polynomials and maximally smooth polynomials, like those employed by the EDGES and SARAS3 experiments, cannot produce good fits with five parameters for the experimental simulations in this work (p < 10−6). However, we find that polynomials with six parameters pass the K-S test (p = 0.4), although a nine-parameter fit produces the highest p-value (p ∼ 0.67). When fitting multiple LST bins simultaneously, we find that the linear forward model outperforms (a higher p-value) the nonlinear model for 2, 5, and 10 LST bins. Importantly, the K-S test consistently identifies best-fit and preferred models.

Detecting Cosmic 21 cm global signal using an improved polynomial fitting algorithm

ABSTRACT Detecting the cosmic 21 cm signal from epoch of reionization has always been a difficult task. Although the Galactic foreground can be regarded as a smooth power-law spectrum, due to the chromaticity of the antenna, additional structure will be introduced into the global spectrum, making the polynomial fitting algorithm perform poorly. In this paper, we introduce an improved polynomial fitting algorithm – the Vari-Zeroth-Order Polynomial (VZOP) fitting and use it to fit the simulation data. This algorithm is developed for the upcoming low-frequency anechoic chamber experiment, yet it is a general method suitable for application in any single antenna-based global 21 cm signal experiment. VZOP defines a 24-h averaged beam model that brings information about the antenna beam into the polynomial model. Assuming that the beam can be measured, VZOP can successfully recover the 21 cm absorption feature, even if the beam is extremely frequency-dependent. In real observations, due to various systematics, the corrected measured beam contains residual errors that are not completely random. Assuming the errors are frequency-dependent, VZOP is capable of recovering the 21 cm absorption feature even when the error reaches 10 per cent. Even in the most extreme scenario where the errors are completely random, VZOP can at least give a fitting result that is not worse than the common polynomial fitting. In conclusion, the fitting effect of VZOP depends on the structure of the error and the accuracy of the beam measurement.

Impact of Galactic dust non-Gaussianity on searches for B-modes from inflation

ABSTRACT A key challenge in the search for primordial B-modes is the presence of polarized Galactic foregrounds, especially thermal dust emission. Power-spectrum-based analysis methods generally assume the foregrounds to be Gaussian random fields when constructing a likelihood and computing the covariance matrix. In this paper, we investigate how non-Gaussianity in the dust field instead affects CMB and foreground parameter inference in the context of inflationary B-mode searches, capturing this effect via modifications to the dust power-spectrum covariance matrix. For upcoming experiments such as the Simons Observatory, we find no dependence of the tensor-to-scalar ratio uncertainty $\sigma (r)$ on the degree of dust non-Gaussianity or the nature of the dust covariance matrix. We provide an explanation of this result, noting that when frequency decorrelation is negligible, dust in mid-frequency channels is cleaned using high-frequency data in a way that is independent of the spatial statistics of dust. We show that our results hold also for non-zero levels of frequency decorrelation that are compatible with existing data. We find, however, that neglecting the impact of dust non-Gaussianity in the covariance matrix can lead to inaccuracies in goodness-of-fit metrics. Care must thus be taken when using such metrics to test B-mode spectra and models, although we show that any such problems can be mitigated by using only cleaned spectrum combinations when computing goodness-of-fit statistics.

Mitigating the effect of 1/f noise on the detection of the H i intensity mapping power spectrum from single-dish measurements

ABSTRACT We present and compare several methods to mitigate time-correlated (1/f) noise within the H i intensity mapping component of the MeerKAT Large Area Synoptic Survey (MeerKLASS). By simulating scan strategies, the H i signal, foreground emissions, white and correlated noise, we assess the ability of various data-processing pipelines to recover the power spectrum of H i brightness temperature fluctuations. We use MeerKAT pilot data to assess the level of 1/f noise expected for the MeerKLASS survey and use these measurements to create realistic levels of time-correlated noise for our simulations. We find the time-correlated noise component within the pilot data to be between 10 and 20 times higher than the white noise level at the scale of $k = 0.04 \, {\rm {Mpc}}^{-1}$. Having determined that the MeerKAT 1/f noise is partially correlated across all the frequency channels, we employ Singular Value Decomposition (SVD) as a technique to remove both the 1/f noise and Galactic foregrounds but find that over-cleaning results in the removal of H i power at large (angular and radial) scales; a power loss of 40 per cent is seen for a 3-mode SVD clean at the scale of $k = 0.04 \, {\rm {Mpc}}^{-1}$. We compare the impact of map-making using weighting by the full noise covariance (i.e. including a 1/f component), as opposed to just a simple unweighted binning, finding that including the time-correlated noise information reduces the excess power added by 1/f noise by up to 30 per cent.

Characterizing cosmic birefringence in the presence of Galactic foregrounds and instrumental systematic effects

Characterizing cosmic birefringence in the presence of Galactic foregrounds and instrumental systematic effects

Dark matter search with CMB: a study of foregrounds

The energy injected from dark matter annihilation and decay processes potentially raises the ionisation of the intergalactic medium and leaves visible footprints on the anisotropy maps of the cosmic microwave background (CMB). Galactic foregrounds emission in the microwave bands contaminate the CMB measurement and may affect the search for dark matter's signature. In this paper, we construct a full CMB data and foreground simulation based on the design of the next-generation ground-based CMB experiments. The foreground residual after the components separation on maps is fully considered in our data analysis, accounting for various contamination from the emission of synchrotron, thermal dust, free-free and spinning dust. We analyse the corresponding sensitivity on dark matter parameters from the temperature and polarization maps, and we find that the CMB foregrounds leave a non-zero yet controllable impact on the sensitivity. Comparing with statistics-only analysis, the CMB foreground residual leads to a factor of at most 19% weakening on energy-injection constraints, depending on the specific dark matter process and experimental configuration. Strong limits on dark matter annihilation rate and decay lifetime can be expected after foreground subtraction.