Line-intensity mapping is an emerging field of observational work, with strong potential to fit into a larger effort to probe large-scale structure and small-scale astrophysical phenomena using ...multiple complementary tracers. Taking full advantage of such complementarity means, in part, undertaking line-intensity surveys with galaxy surveys in mind. We consider the potential for detection of a cross-correlation signal between COMAP and blind surveys based on photometric redshifts (as in COSMOS) or based on spectroscopic data (as with the HETDEX survey of Ly emitters). We find that obtaining accuracy in redshifts and 10−4 sources per Mpc3 with spectroscopic redshift determination should enable a CO-galaxy cross spectrum detection significance at least twice that of the CO auto spectrum. Either a future targeted spectroscopic survey or a blind survey like HETDEX may be able to meet both of these requirements.
Abstract
We present the power spectrum methodology used for the first-season COMAP analysis, and assess the quality of the current data set. The main results are derived through the Feed–Feed ...Pseudo-Cross-Spectrum (FPXS) method, which is a robust estimator with respect to both noise modeling errors and experimental systematics. We use effective transfer functions to take into account the effects of instrumental beam smoothing and various filter operations applied during the low-level data processing. The power spectra estimated in this way have allowed us to identify a systematic error associated with one of our two scanning strategies, believed to be due to residual ground or atmospheric contamination. We omit these data from our analysis and no longer use this scanning technique for observations. We present the power spectra from our first season of observing, and demonstrate that the uncertainties are integrating as expected for uncorrelated noise, with any residual systematics suppressed to a level below the noise. Using the FPXS method, and combining data on scales
k
= 0.051–0.62 Mpc
−1
, we estimate
P
CO
(k) = −2. 7 ± 1.7 × 10
4
μ
K
2
Mpc
3
, the first direct 3D constraint on the clustering component of the CO(1–0) power spectrum in the literature.
Abstract
We describe the first-season CO Mapping Array Project (COMAP) analysis pipeline that converts raw detector readouts to calibrated sky maps. This pipeline implements four main steps: gain ...calibration, filtering, data selection, and mapmaking. Absolute gain calibration relies on a combination of instrumental and astrophysical sources, while relative gain calibration exploits real-time total-power variations. High-efficiency filtering is achieved through spectroscopic common-mode rejection within and across receivers, resulting in nearly uncorrelated white noise within single-frequency channels. Consequently, near-optimal but biased maps are produced by binning the filtered time stream into pixelized maps; the corresponding signal bias transfer function is estimated through simulations. Data selection is performed automatically through a series of goodness-of-fit statistics, including
χ
2
and multiscale correlation tests. Applying this pipeline to the first-season COMAP data, we produce a data set with very low levels of correlated noise. We find that one of our two scanning strategies (the Lissajous type) is sensitive to residual instrumental systematics. As a result, we no longer use this type of scan and exclude data taken this way from our Season 1 power spectrum estimates. We perform a careful analysis of our data processing and observing efficiencies and take account of planned improvements to estimate our future performance. Power spectrum results derived from the first-season COMAP maps are presented and discussed in companion papers.
COMAP Early Science. I. Overview Cleary, Kieran A.; Borowska, Jowita; Breysse, Patrick C. ...
The Astrophysical journal,
07/2022, Letnik:
933, Številka:
2
Journal Article
Recenzirano
Odprti dostop
Abstract
The CO Mapping Array Project (COMAP) aims to use line-intensity mapping of carbon monoxide (CO) to trace the distribution and global properties of galaxies over cosmic time, back to the ...Epoch of Reionization (EoR). To validate the technologies and techniques needed for this goal, a Pathfinder instrument has been constructed and fielded. Sensitive to CO(1–0) emission from
z
= 2.4–3.4 and a fainter contribution from CO(2–1) at
z
= 6–8, the Pathfinder is surveying 12 deg
2
in a 5 yr observing campaign to detect the CO signal from
z
∼ 3. Using data from the first 13 months of observing, we estimate
P
CO
(
k
) = −2.7 ± 1.7 × 10
4
μ
K
2
Mpc
3
on scales
k
= 0.051 −0.62 Mpc
−1
, the first direct three-dimensional constraint on the clustering component of the CO(1–0) power spectrum. Based on these observations alone, we obtain a constraint on the amplitude of the clustering component (the squared mean CO line temperature bias product) of
Tb
2
<
49
μ
K
2
, nearly an order-of-magnitude improvement on the previous best measurement. These constraints allow us to rule out two models from the literature. We forecast a detection of the power spectrum after 5 yr with signal-to-noise ratio (S/N) 9–17. Cross-correlation with an overlapping galaxy survey will yield a detection of the CO–galaxy power spectrum with S/N of 19. We are also conducting a 30 GHz survey of the Galactic plane and present a preliminary map. Looking to the future of COMAP, we examine the prospects for future phases of the experiment to detect and characterize the CO signal from the EoR.
BEYONDPLANCK Galloway, M.; Andersen, K. J.; Aurlien, R. ...
Astronomy and astrophysics (Berlin),
06/2023, Letnik:
675
Journal Article
Recenzirano
Odprti dostop
We describe the computational infrastructure for end-to-end Bayesian cosmic microwave background (CMB) analysis implemented by the BeyondPlanck Collaboration. The code is called
Commander3
. It ...provides a statistically consistent framework for global analysis of CMB and microwave observations and may be useful for a wide range of legacy, current, and future experiments. The paper has three main goals. Firstly, we provide a high-level overview of the existing code base, aiming to guide readers who wish to extend and adapt the code according to their own needs or re-implement it from scratch in a different programming language. Secondly, we discuss some critical computational challenges that arise within any global CMB analysis framework, for instance in-memory compression of time-ordered data, fast Fourier transform optimization, and parallelization and load-balancing. Thirdly, we quantify the CPU and RAM requirements for the current B
EYOND
P
LANCK
analysis, finding that a total of 1.5 TB of RAM is required for efficient analysis and that the total cost of a full Gibbs sample for LFI is 170 CPU-hrs, including both low-level processing and high-level component separation, which is well within the capabilities of current low-cost computing facilities. The existing code base is made publicly available under a GNU General Public Library (GPL) license.
Abstract
Line-intensity mapping observations will find fluctuations of integrated line emission are attenuated by varying degrees at small scales due to the width of the line emission profiles. This ...attenuation may significantly impact estimates of astrophysical or cosmological quantities derived from measurements. We consider a theoretical treatment of the effect of line broadening on both the clustering and shot-noise components of the power spectrum of a generic line-intensity power spectrum using a halo model. We then consider possible simplifications to allow easier application in analysis, particularly in the context of inferences that require numerous, repeated, fast computations of model line-intensity signals across a large parameter space. For the CO Mapping Array Project and the CO(1–0) line-intensity field at
z
∼ 3 serving as our primary case study, we expect a ∼10% attenuation of the spherically averaged power spectrum on average at relevant scales of
k
≈ 0.2–0.3 Mpc
−1
compared to ∼25% for the interferometric Millimetre-wave Intensity Mapping Experiment targeting shot noise from CO lines at
z
∼ 1–5 at scales of
k
≳ 1 Mpc
−1
. We also consider the nature and amplitude of errors introduced by simplified treatments of line broadening and find that while an approximation using a single effective velocity scale is sufficient for spherically averaged power spectra, a more careful treatment is necessary when considering other statistics such as higher multipoles of the anisotropic power spectrum or the voxel intensity distribution.
Abstract
We present the current state of models for the
z
∼ 3 carbon monoxide (CO) line intensity signal targeted by the CO Mapping Array Project (COMAP) Pathfinder in the context of its early ...science results. Our fiducial model, relating dark matter halo properties to CO luminosities, informs parameter priors with empirical models of the galaxy–halo connection and previous CO (1–0) observations. The Pathfinder early science data spanning wavenumbers
k
= 0.051–0.62 Mpc
−1
represent the first direct 3D constraint on the clustering component of the CO (1–0) power spectrum. Our 95% upper limit on the redshift-space clustering amplitude
A
clust
≲ 70
μ
K
2
greatly improves on the indirect upper limit of 420
μ
K
2
reported from the CO Power Spectrum Survey (COPSS) measurement at
k
∼ 1 Mpc
−1
. The COMAP limit excludes a subset of models from previous literature and constrains interpretation of the COPSS results, demonstrating the complementary nature of COMAP and interferometric CO surveys. Using line bias expectations from our priors, we also constrain the squared mean line intensity–bias product,
Tb
2
≲ 50
μ
K
2
, and the cosmic molecular gas density,
ρ
H2
< 2.5 × 10
8
M
⊙
Mpc
−3
(95% upper limits). Based on early instrument performance and our current CO signal estimates, we forecast that the 5 yr Pathfinder campaign will detect the CO power spectrum with overall signal-to-noise ratio of 9–17. Between then and now, we also expect to detect the CO–galaxy cross-spectrum using overlapping galaxy survey data, enabling enhanced inferences of cosmic star formation and galaxy evolution history.
BEYONDPLANCK Andersen, K. J.; Herman, D.; Aurlien, R. ...
Astronomy and astrophysics (Berlin),
06/2023, Letnik:
675
Journal Article
Recenzirano
Odprti dostop
We present the intensity foreground algorithms and model employed within the B
EYOND
P
LANCK
analysis framework. The B
EYOND
P
LANCK
analysis is aimed at integrating component separation and ...instrumental parameter sampling within a global framework, leading to complete end-to-end error propagation in the
Planck
Low Frequency Instrument (LFI) data analysis. Given the scope of the B
EYOND
P
LANCK
analysis, a limited set of data is included in the component separation process, leading to foreground parameter degeneracies. In order to properly constrain the Galactic foreground parameters, we improve upon the previous
Commander
component separation implementation by adding a suite of algorithmic techniques. These algorithms are designed to improve the stability and computational efficiency for weakly constrained posterior distributions. These are: (1) joint foreground spectral parameter and amplitude sampling, building on ideas from M
IRAMARE
; (2) component-based monopole determination; (3) joint spectral parameter and monopole sampling; and (4) application of informative spatial priors for component amplitude maps. We find that the only spectral parameter with a significant signal-to-noise ratio using the current B
EYOND
P
LANCK
data set is the peak frequency of the anomalous microwave emission component, for which we find
ν
p
= 25.3 ± 0.5 GHz; all others must be constrained through external priors. Future works will be aimed at integrating many more data sets into this analysis, both map and time-ordered based, thereby gradually eliminating the currently observed degeneracies in a controlled manner with respect to both instrumental systematic effects and astrophysical degeneracies. When this happens, the simple LFI-oriented data model employed in the current work will need to be generalized to account for both a richer astrophysical model and additional instrumental effects. This work will be organized within the Open Science-based C
OSMOGLOBE
community effort.
BEYONDPLANCK Svalheim, T. L.; Zonca, A.; Andersen, K. J. ...
Astronomy and astrophysics (Berlin),
07/2023, Letnik:
675
Journal Article
Recenzirano
Odprti dostop
We discuss the treatment of bandpass and beam leakage corrections in the Bayesian B
EYOND
P
LANCK
cosmic microwave background (CMB) analysis pipeline as applied to the
Planck
LFI measurements. As a ...preparatory step, we first applied three corrections to the nominal LFI bandpass profiles, including the removal of a known systematic effect in the ground measuring equipment at 61 GHz, along with a smoothing of standing wave ripples and edge regularization. The main net impact of these modifications is an overall shift in the 70 GHz bandpass of +0.6 GHz. We argue that any analysis of LFI data products, either from
Planck
or B
EYOND
P
LANCK
, should use these new bandpasses. In addition, we fit a single free bandpass parameter for each radiometer of the form Δ
i
= Δ
0
+
δ
i
, where Δ
0
represents an absolute frequency shift per frequency band and
δ
i
is a relative shift per detector. The absolute correction is only fitted at 30 GHz, with a full
χ
2
-based likelihood, resulting in a correction of Δ
30
= 0.24 ± 0.03 GHz. The relative corrections were fitted using a spurious map approach that is fundamentally similar to the method pioneered by the WMAP team, but excluding the introduction of many additional degrees of freedom. All the bandpass parameters were sampled using a standard Metropolis sampler within the main B
EYOND
P
LANCK
Gibbs chain and the bandpass uncertainties were thus propagated to all other data products in the analysis. In summary, we find that our bandpass model significantly reduces leakage effects. For beam leakage corrections, we adopted the official
Planck
LFI beam estimates without any additional degrees of freedom and we only marginalized over the underlying sky model. We note that this is the first time that leakage from beam mismatch has been included for
Planck
LFI maps.
BEYONDPLANCK Svalheim, T. L.; Andersen, K. J.; Aurlien, R. ...
Astronomy and astrophysics (Berlin),
06/2023, Letnik:
675
Journal Article
Recenzirano
Odprti dostop
Using the
Planck
Low Frequency Instrument (LFI) and WMAP data within the global Bayesian B
EYOND
P
LANCK
framework, we constrained the polarized foreground emission between 30 and 70 GHz. We ...combined, for the first time, full-resolution
Planck
LFI time-ordered data with low-resolution WMAP sky maps at 33, 40, and 61 GHz. The spectral parameters were fit with a likelihood defined at the native resolution of each frequency channel. This analysis represents the first implementation of true multi-resolution component separation applied to CMB observations for both amplitude and spectral energy distribution (SED) parameters. For the synchrotron emission, we approximated the SED as a power-law in frequency and we find that the low signal-to-noise ratio of the current data strongly limits the number of free parameters that can be robustly constrained. We partitioned the sky into four large disjoint regions (High Latitude; Galactic Spur; Galactic Plane; and Galactic Center), each associated with its own power-law index. We find that the High Latitude region is prior-dominated, while the Galactic Center region is contaminated by residual instrumental systematics. The two remaining regions appear to be signal-dominated, and for these we derive spectral indices of
β
s
Spur
= −3.17 ± 0.06 and β
s
Plane
= −3.03 ± 0.07, which is in good agreement with previous results. For the thermal dust emission, we assumed a modified blackbody model and we fit a single power-law index across the full sky. We find
β
d
= 1.64 ± 0.03, which is slightly steeper than the value reported in
Planck
HFI data, but still statistically consistent at the 2
σ
confidence level.