Foreground Emission Research Articles

Determination of miscalibrated polarization angles from observed cosmic microwave background and foreground EB power spectra: Application to partial-sky observation

Abstract We study a strategy to determine miscalibrated polarization angles of cosmic microwave background (CMB) experiments using the observed $EB$ polarization power spectra of CMB and Galactic foreground emission. We apply the methodology of Y. Minami et al. (Prog. Theor. Exp. Phys. 2019, 083E02, 2019), developed for full-sky observations to ground-based experiments such as Simons Observatory. We take into account the $E$-to-$B$ leakage and $\ell$-to-$\ell$ covariance due to partial sky coverage using the public code NaMaster. We show that our method yields an unbiased estimate of miscalibrated angles. Our method also enables simultaneous determination of miscalibrated angles and the intrinsic $EB$ power spectrum of polarized dust emission when the latter is proportional to $\sqrt{C_\ell^{EE}C_\ell^{BB}}$ and $C_\ell^{BB}$ is proportional to $C_\ell^{EE}$.

Impact of polarized galactic foreground emission on CMB lensing reconstruction and delensing of B-modes

Next generation CMB experiments such as CMB-S4 aim at measuring the CMB lensing potential at sub-percent precision where most of the constraining power will come from CMB polarization. We investigate the prospects of achieving this goal in the presence of large-scale, diffuse galactic foreground emission by using non-Gaussian sky simulations and exploit multi-frequency information to clean those. We show that, while prior to foreground cleaning, cosmological parameter estimates from the contaminated lensing potential estimation with the quadratic estimator can be significantly biased, these can be successfully mitigated by applying a parametric foreground cleaning approach. We further observe no significant additional bias in the delensed B-mode power spectrum after applying foreground cleaning and are therefore able to obtain an unbiased measurement of the tensor-to-scalar ratio, r, after delensing with the quadratic estimate of the lensing potential.

Foreground modelling via Gaussian process regression: an application to HERA data

ABSTRACT The key challenge in the observation of the redshifted 21-cm signal from cosmic reionization is its separation from the much brighter foreground emission. Such separation relies on the different spectral properties of the two components, although, in real life, the foreground intrinsic spectrum is often corrupted by the instrumental response, inducing systematic effects that can further jeopardize the measurement of the 21-cm signal. In this paper, we use Gaussian Process Regression to model both foreground emission and instrumental systematics in ∼2 h of data from the Hydrogen Epoch of Reionization Array. We find that a simple co-variance model with three components matches the data well, giving a residual power spectrum with white noise properties. These consist of an ‘intrinsic’ and instrumentally corrupted component with a coherence scale of 20 and 2.4 MHz, respectively (dominating the line-of-sight power spectrum over scales k∥ ≤ 0.2 h cMpc−1) and a baseline-dependent periodic signal with a period of ∼1 MHz (dominating over k∥ ∼ 0.4–0.8 h cMpc−1), which should be distinguishable from the 21-cm Epoch of Reionization signal whose typical coherence scale is ∼0.8 MHz.

HIR4: cosmology from a simulated neutral hydrogen full sky using Horizon Run 4

ABSTRACT The distribution of cosmological neutral hydrogen will provide a new window into the large-scale structure of the Universe with the next generation of radio telescopes and surveys. The observation of this material, through 21 cm line emission, will be confused by foreground emission in the same frequencies. Even after these foregrounds are removed, the reconstructed map may not exactly match the original cosmological signal, which will introduce systematic errors and offset into the measured correlations. In this paper, we simulate future surveys of neutral hydrogen using the Horizon Run 4 (HR4) cosmological N-body simulation. We generate H i intensity maps from the HR4 halo catalogue, and combine with foreground radio emission maps from the Global Sky Model, to create accurate simulations over the entire sky. We simulate the H i sky for the frequency range 700–800 MHz, matching the sensitivity of the Tianlai pathfinder. We test the accuracy of the fastICA, PCA, and log-polynomial fitting foreground removal methods to recover the input cosmological angular power spectrum and measure the parameters. We show the effect of survey noise levels and beam sizes on the recovered the cosmological constraints. We find that while the reconstruction removes power from the cosmological 21 cm distribution on large scales, we can correct for this and recover the input parameters in the noise-free case. However, the effect of noise and beam size of the Tianlai pathfinder prevents accurate recovery of the cosmological parameters when using only intensity mapping information.

Characterizing the Sensitivity of 40 GHz TES Bolometers for BICEP Array

The BICEP/Keck (BK) experiment aims to detect the imprint of primordial gravitational waves in the Cosmic Microwave Background polarization, which would be direct evidence of the inflation theory. While the tensor-to-scalar ratio has been constrained to be r_0.05 < 0.06 at 95% c.l., further improvements on this upper limit are hindered by polarized Galactic foreground emissions and removal of gravitational lensing polarization. The 30/40 GHz receiver of the BICEP Array (BA) will deploy at the end of 2019 and will constrain the synchrotron foreground with unprecedented accuracy within the BK sky patch. We will show the design of the 30/40 GHz detectors and test results summarizing its performance. The low optical and atmospheric loading at these frequencies requires our TES detectors to have low saturation power in order to be photon-noise dominated. To realize the low thermal conductivity required from a 250 mK base temperature, we developed new bolometer leg designs. We will present the relevant measured detector parameters: G, Tc, Rn, Psat , and spectral bands, and noise spectra. We achieved a per bolometer NEP including all noise components of 2.07E-17 W/sqrt(Hz), including an anticipated photon noise level 1.54E-17 W/sqrt(Hz).

Signals of axion like dark matter in time dependent polarization of light

We consider the search for axion-like particles (ALPs) by using time series data of the polarization angle of the light. If the condensation of an ALP plays the role of dark matter, the polarization plane of the light oscillates as a function of time and we may be able to detect the signal of the ALP by continuously observing the polarization. In particular, we discuss that the analysis of the Fourier-transformed data of the time-dependent polarization angle is powerful to find the signal of the ALP dark matter. We pay particular attention to the light coming from astrophysical sources such as protoplanetary disks, supernova remnants, the foreground emission of the cosmic microwave background, and so on. We show that, for the ALP mass of ∼10−22–10−19eV, ALP searches in the Fourier space may reach the parameter region which is unexplored by other searches yet.

Improved upper limits on the 21 cm signal power spectrum of neutral hydrogen at z ≈ 9.1 from LOFAR

ABSTRACTA new upper limit on the 21 cm signal power spectrum at a redshift of z ≈ 9.1 is presented, based on 141 h of data obtained with the Low-Frequency Array (LOFAR). The analysis includes significant improvements in spectrally smooth gain-calibration, Gaussian Process Regression (GPR) foreground mitigation and optimally weighted power spectrum inference. Previously seen ‘excess power’ due to spectral structure in the gain solutions has markedly reduced but some excess power still remains with a spectral correlation distinct from thermal noise. This excess has a spectral coherence scale of 0.25–0.45 MHz and is partially correlated between nights, especially in the foreground wedge region. The correlation is stronger between nights covering similar local sidereal times. A best 2-σ upper limit of $\Delta ^2_{21} \lt (73)^2\, \mathrm{mK^2}$ at $k = 0.075\, \mathrm{h\, cMpc^{-1}}$ is found, an improvement by a factor ≈8 in power compared to the previously reported upper limit. The remaining excess power could be due to residual foreground emission from sources or diffuse emission far away from the phase centre, polarization leakage, chromatic calibration errors, ionosphere, or low-level radiofrequency interference. We discuss future improvements to the signal processing chain that can further reduce or even eliminate these causes of excess power.

Updated Design of the CMB Polarization Experiment Satellite LiteBIRD

Recent developments of transition-edge sensors (TESs), based on extensive experience in ground-based experiments, have been making the sensor techniques mature enough for their application on future satellite cosmic microwave background (CMB) polarization experiments. LiteBIRD is in the most advanced phase among such future satellites, targeting its launch in Japanese Fiscal Year 2027 (2027FY) with JAXA’s H3 rocket. It will accommodate more than 4000 TESs in focal planes of reflective low-frequency and refractive medium-and-high-frequency telescopes in order to detect a signature imprinted on the CMB by the primordial gravitational waves predicted in cosmic inflation. The total wide frequency coverage between 34 and 448 GHz enables us to extract such weak spiral polarization patterns through the precise subtraction of our Galaxy’s foreground emission by using spectral differences among CMB and foreground signals. Telescopes are cooled down to 5 K for suppressing thermal noise and contain polarization modulators with transmissive half-wave plates at individual apertures for separating sky polarization signals from artificial polarization and for mitigating from instrumental 1/f noise. Passive cooling by using V-grooves supports active cooling with mechanical coolers as well as adiabatic demagnetization refrigerators. Sky observations from the second Sun–Earth Lagrangian point, L2, are planned for 3 years. An international collaboration between Japan, the USA, Canada, and Europe is sharing various roles. In May 2019, the Institute of Space and Astronautical Science, JAXA, selected LiteBIRD as the strategic large mission No. 2.

Comparing simulated 26Al maps to gamma-ray measurements

Context. The diffuse gamma-ray emission of 26Al at 1.8 MeV reflects ongoing nucleosynthesis in the Milky Way and traces massive-star feedback in the interstellar medium due to its 1 Myr radioactive lifetime. The morphology and dynamics of the interstellar medium are investigated in astrophysics through 3D hydrodynamic simulations in fine detail as there are few suitable astronomical probes available. Aims. We aim to compare a galactic-scale hydrodynamic simulation of the Galaxy’s interstellar medium, including feedback and nucleosynthesis, with gamma-ray data on 26Al emission in the Milky Way, extracting constraints that are only weakly dependent on the particular realisation of the simulation or Galaxy structure. Methods. Due to constraints and biases in both the simulations and the gamma-ray observations, such comparisons are not straightforward. For a direct comparison, we performed maximum likelihood fits of both simulated sky maps and observation-based maximum entropy maps to measurements using INTEGRAL/SPI. In order to study general morphological properties, we compare the scale heights of 26Al emission produced by the simulation to INTEGRAL/SPI measurements. Results. The direct comparison shows that the simulation describes the observed inner Galaxy well, however it differs significantly from the observed full-sky emission morphology. Comparing the scale height distribution, we see similarities for small-scale height features and a mismatch at larger-scale heights. We attribute this to prominent foreground emission sites which are not captured by the simulation.

An Empirical Model of Energetic Neutral Atom Imaging of the Heliosphere and Its Implications for Future Heliospheric Missions at Great Heliocentric Distances

Abstract Several concepts for heliospheric missions operating at heliocentric distances far beyond Earth orbit are currently investigated by the scientific community. The mission concept of the Interstellar Probe, e.g., aims at reaching a distance of 1000 au away from the Sun within this century. This would allow the coming generation to obtain a global view of our heliosphere from an outside vantage point by measuring the energetic neutral atoms (ENAs) originating from the various plasma regions. It would also allow for direct sampling of the unperturbed interstellar medium, as well as for many observation opportunities beyond heliospheric science, such as visits to Kuiper Belt objects, a comprehensive view on the interplanetary dust populations, and infrared astronomy free from the foreground emission of the zodiacal cloud. In this study, we present a simple empirical model of ENAs from the heliosphere and derive basic requirements for ENA instrumentation on board a spacecraft at great heliocentric distances. We consider the full energy range of heliospheric ENAs from 10 eV to 100 keV because each part of the energy spectrum has its own merits for heliospheric science. To cover the full ENA energy range, two or three different ENA instruments are needed. Thanks to parallax observations, some insights about the nature of the IBEX ribbon and the dimensions of the heliosphere can already be gained by ENA imaging from a few au heliocentric distance. To directly reveal the global shape of the heliosphere, measurements from outside the heliosphere are, of course, the best option.

Mitigating Internal Instrument Coupling for 21 cm Cosmology. I. Temporal and Spectral Modeling in Simulations

Abstract We study the behavior of internal signal chain reflections and antenna cross coupling as systematics for 21 cm cosmological surveys. We outline the mathematics for how these systematics appear in interferometric visibilities and describe their phenomenology. We then describe techniques for modeling and removing these systematics without attenuating the 21 cm signal in the data. This has critical implications for low-frequency radio surveys aiming to characterize the 21 cm signal from the Epoch of Reionization (EoR) and Cosmic Dawn, as systematics can cause bright foreground emission to contaminate the EoR window and prohibit a robust detection. We also quantify the signal loss properties of the systematic modeling algorithms, and show that our techniques demonstrate resistance against EoR signal loss. In a companion paper, we demonstrate these methods on data from the Hydrogen Epoch of Reionization Array as a proof-of-concept.

Assessment of the Projection-induced Polarimetry Technique for Constraining the Foreground Spectrum in Global 21 cm Cosmology

Abstract Detecting the cosmological sky-averaged (global) 21 cm signal as a function of observed frequency will provide a powerful tool to study the ionization and thermal history of the intergalactic medium (IGM) in the early universe (∼400 million years after the big bang). The greatest challenge in conventional total-power global 21 cm experiments is the removal of the foreground synchrotron emission (∼103–104 K) to uncover the weak cosmological signal (tens to hundreds of millikelvin), especially because the intrinsic smoothness of the foreground spectrum is corrupted by instrumental effects. Although the EDGES (Experiment to Detect the Global EoR Signature) team has recently reported an absorption profile at 78 MHz in the sky-averaged spectrum, it is necessary to confirm this detection with an independent approach. The projection effect from observing anisotropic foreground source emission with a wide-view antenna pointing at the North Celestial Pole can induce a net polarization, referred to as the projection-induced polarization effect (PIPE). Due to Earth’s rotation, observations centered at the circumpolar region will impose a dynamic sky modulation on the net polarization’s waveforms that is unique to the foreground component. In this study, we review the implementation practicality and underlying instrumental effects of this new polarimetry-based technique with detailed numerical simulations and a test-bed instrument, the Cosmic Twilight Polarimeter. In addition, we explore a singular value decomposition–based analysis approach for separating the foreground and instrumental effects from the background global 21 cm signal using the sky-modulated PIPE.

Keck/OSIRIS IFU Detection of a z ∼ 3 Damped Lyα Host Galaxy∗

Abstract We present Keck/OSIRIS infrared IFU observations of the z = 3.153 sub-DLA DLA2233+131, previously detected in absorption to a background quasar and studied with single-slit spectroscopy and Potsdam Multi Aperture Spectrophotometer integral field spectroscopy (IFU). We used the Laser Guide Star Adaptive Optics and OSIRIS IFU to reduce the point-spread function of the background quasar to FWHM ∼ 0.″15 and marginally resolve extended, foreground DLA emission. We detect emission with a flux erg s−1 cm−2, as well as unresolved and Hβλ4861 emission. Using a composite spectrum over the emission region, we measure dynamical mass M ⊙. We made several estimates of star formation rate (SFR) using and Hβλ4861 emission, and measured a SFR of ∼7.1 − 13.6 M ⊙ yr−1. We map and Hβλ4861 emission and the corresponding velocity fields to search for signs of kinematic structure. These maps allow for a more detailed kinematic analysis than previously possible for this galaxy. While some regions show slightly red and blueshifted emission indicative of potential edge-on disk rotation, the data are insufficient to support this interpretation.

On the contamination of the global 21 cm signal from polarized foregrounds

Abstract Global (i.e. sky-averaged) 21 cm signal experiments can measure the evolution of the universe from the Cosmic Dawn to the Epoch of Reionization. These measurements are challenged by the presence of bright foreground emission that can be separated from the cosmological signal if its spectrum is smooth. This assumption fails in the case of single polarization antennas as they measure linearly polarized foreground emission - which is inevitably Faraday rotated through the interstellar medium. We investigate the impact of Galactic polarized foregrounds on the extraction of the global 21 cm signal through realistic sky and dipole simulations both in a low frequency band from 50 to 100 MHz, where a 21 cm absorption profile is expected, and in a higher frequency band (100 − 200 MHz). We find that the presence of a polarized contaminant with complex frequency structure can bias the amplitude and the shape of the reconstructed signal parameters in both bands. We investigate if polarized foregrounds can explain the unexpected 21 cm Cosmic Dawn signal recently reported by the EDGES collaboration. We find that unaccounted polarized foreground contamination can produce an enhanced and distorted 21 cm absorption trough similar to the anomalous profile reported by Bowman et al. (2018), and whose amplitude is in mild tension with the assumed input Gaussian profile (at ∼1.5σ level). Moreover, we note that, under the hypothesis of contamination from polarized foreground, the amplitude of the reconstructed EDGES signal can be overestimated by around 30%, mitigating the requirement for an explanation based on exotic physics.

The HERA-19 Commissioning Array: Direction-dependent Effects

Abstract Foreground power dominates the measurements of interferometers that seek a statistical detection of highly-redshifted H i emission from the Epoch of Reionization (EoR). The chromaticity of the instrument creates a boundary in the Fourier transform of frequency (proportional to k ∥) between spectrally smooth emission, characteristic of the strong synchrotron foreground (the “wedge”), and the spectrally structured emission from H i in the EoR (the “EoR window”). Faraday rotation can inject spectral structure into otherwise smooth polarized foreground emission, which through instrument effects or miscalibration could possibly pollute the EoR window. For instruments pursuing a “foreground avoidance” strategy of simply measuring in the EoR window, and not attempting to model and remove foregrounds, as is the plan for the first stage of the Hydrogen Epoch of Reionization Array (HERA), characterizing the intrinsic instrument polarization response is particularly important. Using data from the HERA 19-element commissioning array, we investigate the polarization response of this new instrument in the power-spectrum domain. We perform a simple image-based calibration based on the unpolarized diffuse emission of the Global Sky Model, and show that it achieves qualitative redundancy between the nominally redundant baselines of the array and reasonable amplitude accuracy. We construct power spectra of all fully polarized coherencies in all pseudo-Stokes parameters, and discuss the achieved isolation of foreground power due to the intrinsic spectral smoothness of the foregrounds, the instrument chromaticity, and the calibration. We compare to simulations based on an unpolarized diffuse sky model and detailed electromagnetic simulations of the dish and feed, confirming that in Stokes I, the calibration does not add significant spectral structure beyond that expected from the interferometer array configuration and the modeled primary beam response. Furthermore, this calibration is stable over the 8 days of observations considered. Excess power is seen in the power spectra of the linear polarization Stokes parameters, which is not easily attributable to leakage via the primary beam, and results from some combination of residual calibration errors and actual polarized emission. Stokes V is found to be highly discrepant from the expectation of zero power, strongly pointing to the need for more accurate polarized calibration.

Simultaneous determination of the cosmic birefringence and miscalibrated polarization angles from CMB experiments

Abstract We show that the cosmic birefringence and miscalibrated polarization angles can be determined simultaneously by cosmic microwave background (CMB) experiments using the cross-correlation between $E$- and $B$-mode polarization data. This is possible because the polarization angles of the CMB are rotated by both the cosmic birefringence and miscalibration effects, whereas those of the Galactic foreground emission are rotated only by the latter. Our method does not require prior knowledge of the $E$- and $B$-mode power spectra of the foreground emission, but uses only the knowledge of the CMB polarization spectra. Specifically, we relate the observed $EB$ correlation to the difference between the observed$E$- and $B$-mode spectra in the sky, and use different multipole dependences of the CMB (given by theory) and foreground spectra (given by data) to derive the likelihood for the miscalibration angle $\alpha$ and the birefringence angle $\beta$. We show that a future satellite mission similar to LiteBIRD can determine $\beta$ with a precision of 10 arcmin.

Fermi-LAT Observations of γ-Ray Emission toward the Outer Halo of M31.

The Andromeda galaxy is the closest spiral galaxy to us and has been the subject of numerous studies. It harbors a massive dark matter halo, which may span up to ~600 kpc across and comprises ~90% of the galaxy's total mass. This halo size translates into a large diameter of 42° on the sky, for an M31-Milky Way (MW) distance of 785 kpc, but its presumably low surface brightness makes it challenging to detect with γ-ray telescopes. Using 7.6 yr of Fermi Large Area Telescope (Fermi-LAT) observations, we make a detailed study of the γ-ray emission between 1-100 GeV toward M31's outer halo, with a total field radius of 60° centered at M31, and perform an in-depth analysis of the systematic uncertainties related to the observations. We use the cosmic-ray propagation code GALPROP to construct specialized interstellar emission models to characterize the foreground γ-ray emission from the MW, including a self-consistent determination of the isotropic component. We find evidence for an extended excess that appears to be distinct from the conventional MW foreground, having a total radial extension upward of ~120-200 kpc from the center of M31. We discuss plausible interpretations of the excess emission, but emphasize that uncertainties in the MW foreground-and in particular, modeling of the H i-related components-have not been fully explored and may impact the results.

The 21 cm Power Spectrum from the Cosmic Dawn: First Results from the OVRO-LWA

The 21\,cm transition of neutral hydrogen is opening an observational window into the cosmic dawn of the universe---the epoch of first star formation. We use 28\,hr of data from the Owens Valley Radio Observatory Long Wavelength Array (OVRO-LWA) to place upper limits on the spatial power spectrum of 21\,cm emission at $z \approx 18.4$ ($\Delta_{21} \lesssim 10^4\,\text{mK}$), and within the absorption feature reported by the EDGES experiment (Bowman et al. 2018). In the process we demonstrate the first application of the double Karhunen-Lo\`{e}ve transform for foreground filtering, and diagnose the systematic errors that are currently limiting the measurement. We also provide an updated model for the angular power spectrum of low-frequency foreground emission measured from the northern hemisphere, which can be used to refine sensitivity forecasts for next-generation experiments.

Joint estimation of the Epoch of Reionization power spectrum and foregrounds

ABSTRACTThe power spectrum of redshifted 21 cm emission brightness temperature fluctuations is a powerful probe of the Epoch of Reionization (EoR). However, bright foreground emission presents a significant impediment to its unbiased recovery from interferometric data. We estimate the power spectrum within a Bayesian framework and demonstrate that incorporating a priori knowledge of the spectral structure of foregrounds in the large spectral scale component of the data model enables significantly improved modelling of the foregrounds without increasing the model complexity. We explore two astrophysically motivated parametrizations of the large spectral scale model: (i) a constant plus power-law model of the form $q_{0}+q_{1}(\nu /\nu _{0})^{b_{1}}$ for two values of b1: b1 = 〈β〉GDSE and b1 = 〈β〉EGS, the mean spectral indices of the Galactic diffuse synchrotron emission and extragalactic source foreground emission, respectively; and (ii) a constant plus double power-law model of the form $q_{0}+q_{1}(\nu /\nu _{0})^{b_{1}}+q_{2}(\nu /\nu _{0})^{b_{2}}$ with b1 = 〈β〉GDSE and b2 = 〈β〉EGS. We estimate the EoR power spectrum from simulated interferometric data consisting of an EoR signal, Galactic diffuse synchrotron emission, extragalactic sources, and diffuse free–free emission from the Galaxy. We show that, by jointly estimating a model of the EoR signal with the constant plus double power-law parametrization of the large spectral scale model, unbiased estimates of the EoR power spectrum are recoverable on all spatial scales accessible in the data set, including on the large spatial scales that were found to be contaminated in earlier work.

Impact of foregrounds on H i intensity mapping cross-correlations with optical surveys

ABSTRACT The future of precision cosmology could benefit from cross-correlations between intensity maps of unresolved neutral hydrogen (H i) and more conventional optical galaxy surveys. A major challenge that needs to be overcome is removing the 21cm foreground emission that contaminates the cosmological H i signal. Using N-body simulations, we simulate H i intensity maps and optical catalogues that share the same underlying cosmology. Adding simulated foreground contamination and using state-of-the-art reconstruction techniques, we investigate the impacts that 21cm foregrounds and other systematics have on these cross-correlations. We find that the impact a Fast Independent Component Analysis 21cm foreground clean has on the cross-correlations with spectroscopic optical surveys with well-constrained redshifts is minimal. However, problems arise when photometric surveys are considered: We find that a redshift uncertainty σz ≥ 0.04 causes significant degradation in the cross-power spectrum signal. We diagnose the main root of these problems, which relates to arbitrary amplitude changes along the line of sight in the intensity maps caused by the foreground clean and suggest solutions that should be applicable to real data. These solutions involve a reconstruction of the line-of-sight temperature means using the available overlapping optical data along with an artificial extension to the H i data through redshift to address edge effects. We then put these solutions through a further test in a mock experiment that uses a clustering-based redshift estimation technique to constrain the photometric redshifts of the optical sample. We find that with our suggested reconstruction, cross-correlations can be utilized to make an accurate prediction of the optical redshift distribution.