Showing 1–19 of 19 results for author: Naess, S K

Cosmoglobe DR2. II. CIB monopole measurements from COBE-DIRBE through global Bayesian analysis

Authors: D. J. Watts, M. Galloway, E. Gjerløw, M. San, R. Aurlien, A. Basyrov, M. Brilenkov, H. K. Eriksen, U. Fuskeland, D. Herman, H. T. Ihle, J. G. S. Lunde, S. K. Næss, N. -O. Stutzer, H. Thommesen, I. K. Wehus

Abstract: We derive new constraints on the CIB monopole spectrum from reprocessed COBE-DIRBE sky maps with lower instrumental and astrophysical contamination than the legacy DIRBE maps. These maps have been generated through a global Bayesian analysis framework that simultaneously fits cosmological, astrophysical, and instrumental parameters, described in a series of papers referred to as Cosmoglobe DR2. We… ▽ More We derive new constraints on the CIB monopole spectrum from reprocessed COBE-DIRBE sky maps with lower instrumental and astrophysical contamination than the legacy DIRBE maps. These maps have been generated through a global Bayesian analysis framework that simultaneously fits cosmological, astrophysical, and instrumental parameters, described in a series of papers referred to as Cosmoglobe DR2. We have applied this method to the DIRBE Calibrated Individual Observations, complemented by selected HFI and FIRAS sky maps to break key astrophysical degeneracies, as well as WISE and Gaia compact object catalogs. In this paper, we focus on CIB constraints that result from this work. We report detections of an isotropic signal in five out of the ten DIRBE bands (1.25, 2.2, 3.5, 140, and 240 $μ$m). For the 2.2 $μ$m channel, we find an amplitude of $13\pm3\,\mathrm{nW\,m^{-2}\,sr^{-1}}$, 36% lower than reported from the official DIRBE release. For the 240 $μ$m channel, we find $10\pm3\,\mathrm{nW\,m^{-2}\,sr^{-1}}$, 26% lower than the official DIRBE release. We interpret these lower values as resulting from improved zodiacal light and Galactic foreground modeling. For the bands between 4.9 and 100 $μ$m, the presence of excess radiation in solar-centric coordinates reported in a companion paper precludes the definition of lower limits. However, the analysis still provides well-defined upper limits. For the 12 $μ$m channel, we find an upper 95% confidence limit of 103 $\mathrm{nW\,m^{-2}\,sr^{-1}}$, a factor of four lower than the corresponding legacy result of 468 $\mathrm{nW\,m^{-2}\,sr^{-1}}$. The results presented in this paper redefines the state-of-the-art CIB monopole constraints from COBE-DIRBE, and it provides a real-world illustration of the power of global end-to-end analysis of multiple complementary data sets which is the foundational idea of the Cosmoglobe project. △ Less

Submitted 3 June, 2024; originally announced June 2024.

Comments: 14 pages, 7 figures, submitted to A&A

CMB-HD: Astro2020 RFI Response

Authors: Neelima Sehgal, Simone Aiola, Yashar Akrami, Kaustuv moni Basu, Michael Boylan-Kolchin, Sean Bryan, Caitlin M Casey, Sébastien Clesse, Francis-Yan Cyr-Racine, Luca Di Mascolo, Simon Dicker, Thomas Essinger-Hileman, Simone Ferraro, George Fuller, Nicholas Galitzki, Dongwon Han, Matthew Hasselfield, Gil Holder, Bhuvnesh Jain, Bradley R. Johnson, Matthew Johnson, Pamela Klaassen, Amanda MacInnis, Mathew Madhavacheril, Philip Mauskopf , et al. (23 additional authors not shown)

Abstract: CMB-HD is a proposed ultra-deep (0.5 uk-arcmin), high-resolution (15 arcseconds) millimeter-wave survey over half the sky that would answer many outstanding questions in both fundamental physics of the Universe and astrophysics. This survey would be delivered in 7.5 years of observing 20,000 square degrees, using two new 30-meter-class off-axis cross-Dragone telescopes to be located at Cerro Toco… ▽ More CMB-HD is a proposed ultra-deep (0.5 uk-arcmin), high-resolution (15 arcseconds) millimeter-wave survey over half the sky that would answer many outstanding questions in both fundamental physics of the Universe and astrophysics. This survey would be delivered in 7.5 years of observing 20,000 square degrees, using two new 30-meter-class off-axis cross-Dragone telescopes to be located at Cerro Toco in the Atacama Desert. Each telescope would field 800,000 detectors (200,000 pixels), for a total of 1.6 million detectors. △ Less

Submitted 28 February, 2020; originally announced February 2020.

Comments: Response to request for information (RFI) by the Panel of Radio, Millimeter, and Submillimeter Observations from the Ground (RMS) of the Astro2020 Decadal Survey regarding the CMB-HD APC (arXiv:1906.10134). Note some text overlap with original APC. Note also detector count and cost have been reduced by 1/3, and observing time increased by 1/3 compared to original APC; science goals expanded

Constraints on the spectral index of polarized synchrotron emission from WMAP and Faraday-corrected S-PASS data

Authors: U. Fuskeland, K. J. Andersen, R. Aurlien, R. Banerji, M. Brilenkov, H. K. Eriksen, M. Galloway, E. Gjerløw, S. K. Næss, T. L. Svalheim, I. K. Wehus

Abstract: We constrain the spectral index of polarized synchrotron emission, $β_s$, by correlating the recently released 2.3 GHz S-Band Polarization All Sky Survey (S-PASS) data with the 23 GHz 9-year Wilkinson Microwave Anisotropy Probe (WMAP) sky maps. We subdivide the S-PASS field, which covers the southern ecliptic hemisphere, into 95 $15^{\circ}\times15^{\circ}$ regions and estimate the spectral index… ▽ More We constrain the spectral index of polarized synchrotron emission, $β_s$, by correlating the recently released 2.3 GHz S-Band Polarization All Sky Survey (S-PASS) data with the 23 GHz 9-year Wilkinson Microwave Anisotropy Probe (WMAP) sky maps. We subdivide the S-PASS field, which covers the southern ecliptic hemisphere, into 95 $15^{\circ}\times15^{\circ}$ regions and estimate the spectral index of polarized synchrotron emission within each region using a simple but robust T-T plot technique. Three different versions of the S-PASS data are considered, corresponding to: no correction for Faraday rotation; Faraday correction based on the rotation measure model presented by the S-PASS team; or Faraday correction based on a rotation measure model presented by Hutschenreuter and Enßlin. We find that the correlation between S-PASS and WMAP is strongest when applying the S-PASS model. Adopting this correction model, we find that the mean spectral index of polarized synchrotron emission gradually steepens from $β_s\approx-2.8$ at low Galactic latitudes to $β_s\approx-3.3$ at high Galactic latitudes, in good agreement with previously published results. Finally, we consider two special cases defined by the BICEP2 and SPIDER fields and obtain mean estimates of $β_{BICEP2}=-3.22\pm0.06$ and $β_{SPIDER}=-3.21\pm0.03$, respectively. A comparison with a similar analysis performed in the 23-33 GHz range suggests a flattening of about $Δβ_s \sim 0.1 \pm 0.2$ from low to higher frequencies, but with no statistical significance due to high uncertainties. △ Less

Submitted 11 March, 2021; v1 submitted 12 September, 2019; originally announced September 2019.

Comments: 16 pages, 15 figures. Figure number 10 and 15 are new. Published in A&A

Journal ref: A&A 646, A69 (2021)

How to avoid X'es around point sources in maximum likelihood CMB maps

Authors: Sigurd K. Naess

Abstract: The maximum likelihood estimator for CMB map-making is optimal and unbiased as long as the data model is correct, but in practice it rarely is, with model errors including sub-pixel structure and instrumental problems like time-variable gain and pointing errors. In the presence of such errors, the solution is biased, with the local error in each pixel leaking outwards along the scanning pattern by… ▽ More The maximum likelihood estimator for CMB map-making is optimal and unbiased as long as the data model is correct, but in practice it rarely is, with model errors including sub-pixel structure and instrumental problems like time-variable gain and pointing errors. In the presence of such errors, the solution is biased, with the local error in each pixel leaking outwards along the scanning pattern by a noise correlation length. The most important sources of such leakage are strong point sources, and for common scanning patterns the leakage manifests as an X around each such source. I discuss why this happens, and present several old and new methods for mitigating and/or eliminating this leakage, along with a small stand-alone TOD simulator and map-maker in Python that implements them. △ Less

Submitted 24 November, 2019; v1 submitted 19 June, 2019; originally announced June 2019.

Comments: 16 pages, 4 figures, accepted for publication in JCAP

The Atacama Cosmology Telescope: Two-season ACTPol Extragalactic Point Sources and their Polarization properties

Authors: Rahul Datta, Simone Aiola, Steve K. Choi, Mark Devlin, Joanna Dunkley, Rolando Dunner, Patricio A. Gallardo, Megan Gralla, Mark Halpern, Matthew Hasselfield, Matt Hilton, Adam D. Hincks, Shuay-Pwu P. Ho, Johannes Hubmayr, Kevin M. Huffenberger, John P. Hughes, Arthur Kosowsky, Carlos H. Lopez-Caraballo, Thibaut Louis, Marius Lungu, Tobias Marriage, Loic Maurin, Jeff McMahon, Kavilan Moodley, Sigurd K. Naess , et al. (8 additional authors not shown)

Abstract: We report on measurements of the polarization of extragalactic sources at 148 GHz made during the first two seasons of the Atacama Cosmology Telescope Polarization (ACTPol) survey. The survey covered 680 deg$^{2}$ of the sky on the celestial equator. Polarization measurements of 169 intensity-selected sources brighter than 30 mJy, that are predominantly Active Galactic Nuclei, are presented. Above… ▽ More We report on measurements of the polarization of extragalactic sources at 148 GHz made during the first two seasons of the Atacama Cosmology Telescope Polarization (ACTPol) survey. The survey covered 680 deg$^{2}$ of the sky on the celestial equator. Polarization measurements of 169 intensity-selected sources brighter than 30 mJy, that are predominantly Active Galactic Nuclei, are presented. Above a total flux of 215 mJy where the noise bias removal in the polarization measurement is reliable, we detect 26 sources, 14 of which have a detection of linear polarization at greater than 3$σ_{p}$ significance. The distribution of the fractional polarization as a function of total source intensity is analyzed. Our result is consistent with the scenario that the fractional polarization of our measured radio source population is independent of total intensity down to the limits of our measurements and well described by a Gaussian distribution with a mean fractional polarization $p=0.028\pm$0.005 and standard deviation $σ_{\mathrm{p}}=0.054$, truncated at $p=0$. Extrapolating this model for the distribution of source polarization below the ACTPol detection threshold, we predict that one could get a clean measure of the E-mode polarization power spectrum of the microwave background out to $\ell=6000$ with 1 $μ$K-arcminute maps over 10$\%$ of the sky from a future survey. We also study the spectral energy distribution of the total and polarized source flux densities by cross-matching with low radio frequency catalogs. We do not find any correlation between the spectral indices for total flux and polarized flux. △ Less

Submitted 9 November, 2018; v1 submitted 5 November, 2018; originally announced November 2018.

Cosmic Microwave Background Mapmaking with a Messenger Field

Authors: Kevin M. Huffenberger, Sigurd K. Næss

Abstract: We apply a messenger field method to solve the linear minimum-variance mapmaking equation in the context of Cosmic Microwave Background (CMB) observations. In simulations, the method produces sky maps that converge significantly faster than those from a conjugate gradient descent algorithm with a diagonal preconditioner, even though the computational cost per iteration is similar. The messenger me… ▽ More We apply a messenger field method to solve the linear minimum-variance mapmaking equation in the context of Cosmic Microwave Background (CMB) observations. In simulations, the method produces sky maps that converge significantly faster than those from a conjugate gradient descent algorithm with a diagonal preconditioner, even though the computational cost per iteration is similar. The messenger method recovers large scales in the map better than conjugate gradient descent, and yields a lower overall $χ^2$. In the single, pencil beam approximation, each iteration of the messenger mapmaking procedure produces an unbiased map, and the iterations become more optimal as they proceed. A variant of the method can handle differential data or perform deconvolution mapmaking. The messenger method requires no preconditioner, but a high-quality solution needs a cooling parameter to control the convergence. We study the convergence properties of this new method, and discuss how the algorithm is feasible for the large data sets of current and future CMB experiments. △ Less

Submitted 26 June, 2017; v1 submitted 4 May, 2017; originally announced May 2017.

Comments: 12 pages, 6 figures, revision includes new sections that treat differential data. Submitted to ApJ

The Q/U Imaging ExperimenT: Polarization Measurements of the Galactic Plane at 43 and 95 GHz

Authors: QUIET Collaboration, T. M. Ruud, U. Fuskeland, I. K. Wehus, M. Vidal, D. Araujo, C. Bischoff, I. Buder, Y. Chinone, K. Cleary, R. N. Dumoulin, A. Kusaka, R. Monsalve, S. K. Naess, L. B. Newburgh, R. A. Reeves, J. T. L. Zwart, L. Bronfman, R. D. Davies, R. Davis, C. Dickinson, H. K. Eriksen, T. Gaier, J. O. Gundersen, M. Hasegawa , et al. (17 additional authors not shown)

Abstract: We present polarization observations of two Galactic plane fields centered on Galactic coordinates (l,b)=(0 deg,0 deg) and (329 deg, 0 deg) at Q- (43 GHz) and W-band (95 GHz), covering between 301 and 539 square degrees depending on frequency and field. These measurements were made with the QUIET instrument between 2008 October and 2010 December, and include a total of 1263 hours of observations.… ▽ More We present polarization observations of two Galactic plane fields centered on Galactic coordinates (l,b)=(0 deg,0 deg) and (329 deg, 0 deg) at Q- (43 GHz) and W-band (95 GHz), covering between 301 and 539 square degrees depending on frequency and field. These measurements were made with the QUIET instrument between 2008 October and 2010 December, and include a total of 1263 hours of observations. The resulting maps represent the deepest large-area Galactic polarization observations published to date at the relevant frequencies with instrumental rms noise varying between 1.8 and 2.8 uK deg, 2.3-6 times deeper than corresponding WMAP and Planck maps. The angular resolution is 27.3' and 12.8' FWHM at Q- and W-band, respectively. We find excellent agreement between the QUIET and WMAP maps over the entire fields, and no compelling evidence for significant residual instrumental systematic errors in either experiment, whereas the Planck 44 GHz map deviates from these in a manner consistent with reported systematic uncertainties for this channel. We combine QUIET and WMAP data to compute inverse-variance-weighted average maps, effectively retaining small angular scales from QUIET and large angular scales from WMAP. From these combined maps, we derive constraints on several important astrophysical quantities, including a robust detection of polarized synchrotron spectral index steepening of ~0.2 off the plane, as well as the Faraday rotation measure toward the Galactic center (RM=-4000 +/- 200 rad m^-2), all of which are consistent with previously published results. Both the raw QUIET and the co-added QUIET+WMAP maps are made publicly available together with all necessary ancillary information. △ Less

Submitted 11 August, 2015; originally announced August 2015.

Comments: 22 pages, 18 figures; accepted for publication in ApJ; data are available on Lambda (http://lambda.gsfc.nasa.gov)

Journal ref: ApJ 811, 89 (2015)

The Q/U Imaging Experiment: Polarization Measurements of Radio Sources at 43 and 95 GHz

Authors: QUIET Collaboration, K. M. Huffenberger, D. Araujo, C. Bischoff, I. Buder, Y. Chinone, K. Cleary, A. Kusaka, R. Monsalve, S. K. Næss, L. B. Newburgh, R. Reeves, T. M. Ruud, I. K. Wehus, J. T. L. Zwart, C. Dickinson, H. K. Eriksen, T. Gaier, J. O. Gundersen, M. Hasegawa, M. Hazumi, A. D. Miller, S. J. E. Radford, A. C. S. Readhead, S. T. Staggs , et al. (2 additional authors not shown)

Abstract: We present polarization measurements of extragalactic radio sources observed during the Cosmic Microwave Background polarization survey of the Q/U Imaging Experiment (QUIET), operating at 43 GHz (Q-band) and 95 GHz (W-band). We examine sources selected at 20 GHz from the public, $>$40 mJy catalog of the Australia Telescope (AT20G) survey. There are $\sim$480 such sources within QUIET's four low-fo… ▽ More We present polarization measurements of extragalactic radio sources observed during the Cosmic Microwave Background polarization survey of the Q/U Imaging Experiment (QUIET), operating at 43 GHz (Q-band) and 95 GHz (W-band). We examine sources selected at 20 GHz from the public, $>$40 mJy catalog of the Australia Telescope (AT20G) survey. There are $\sim$480 such sources within QUIET's four low-foreground survey patches, including the nearby radio galaxies Centaurus A and Pictor A. The median error on our polarized flux density measurements is 30--40 mJy per Stokes parameter. At S/N $> 3$ significance, we detect linear polarization for seven sources in Q-band and six in W-band; only $1.3 \pm 1.1$ detections per frequency band are expected by chance. For sources without a detection of polarized emission, we find that half of the sources have polarization amplitudes below 90 mJy (Q-band) and 106 mJy (W-band), at 95% confidence. Finally, we compare our polarization measurements to intensity and polarization measurements of the same sources from the literature. For the four sources with WMAP and Planck intensity measurements $>1$ Jy, the polarization fraction are above 1% in both QUIET bands. At high significance, we compute polarization fractions as much as 10--20% for some sources, but the effects of source variability may cut that level in half for contemporaneous comparisons. Our results indicate that simple models---ones that scale a fixed polarization fraction with frequency---are inadequate to model the behavior of these sources and their contributions to polarization maps. △ Less

Submitted 24 November, 2014; originally announced December 2014.

Comments: 16 pages, 10 figures. Submitted to ApJ

Spatial variations in the spectral index of polarized synchrotron emission in the 9 yr WMAP sky maps

Authors: U. Fuskeland, I. K. Wehus, H. K. Eriksen, S. K. Næss

Abstract: We estimate the spectral index, beta, of polarized synchrotron emission as observed in the 9 yr WMAP sky maps using two methods, linear regression ("T-T plot") and maximum likelihood. We partition the sky into 24 disjoint sky regions, and evaluate the spectral index for all polarization angles between 0 deg and 85 deg in steps of 5. Averaging over polarization angles, we derive a mean spectral ind… ▽ More We estimate the spectral index, beta, of polarized synchrotron emission as observed in the 9 yr WMAP sky maps using two methods, linear regression ("T-T plot") and maximum likelihood. We partition the sky into 24 disjoint sky regions, and evaluate the spectral index for all polarization angles between 0 deg and 85 deg in steps of 5. Averaging over polarization angles, we derive a mean spectral index of beta_all-sky=-2.99+-0.01 in the frequency range of 23-33 GHz. We find that the synchrotron spectral index steepens by 0.14 from low to high Galactic latitudes, in agreement with previous studies, with mean spectral indices of beta_plane=-2.98+-0.01 and beta_high-lat=-3.12+-0.04. In addition, we find a significant longitudinal variation along the Galactic plane with a steeper spectral index toward the Galactic center and anticenter than toward the Galactic spiral arms. This can be well modeled by an offset sinusoidal, beta(l)=-2.85+0.17sin(2l-90). Finally, we study synchrotron emission in the BICEP2 field, in an attempt to understand whether the claimed detection of large-scale B-mode polarization could be explained in terms of synchrotron contamination. Adopting a spectral index of beta=-3.12, typical for high Galactic latitudes, we find that the most likely bias corresponds to about 2% of the reported signal (r=0.003). The flattest index allowed by the data in this region is beta=-2.5, and under the assumption of a straight power-law frequency spectrum, we find that synchrotron emission can account for at most 20% of the reported BICEP2 signal. △ Less

Submitted 1 August, 2014; v1 submitted 21 April, 2014; originally announced April 2014.

Comments: 11 pages, 9 figures, updated to match version published in ApJ

Journal ref: ApJ, 790, 104 (2014)

Lensing simulations by Taylor expansion -- not so inefficient after all

Authors: Sigurd K. Næss, Thibaut Louis

Abstract: Cosmic Microwave Background lensing simulation by Taylor expansion has long been considered impractical due to slow convergence, but a recent flat-sky implementation shows that a simple trick eliminates this problem, making Taylor lensing a fast and simple lensing algorithm for the flat sky. Here we generalize the method to the full sky, and study its convergence and performance relative to a comm… ▽ More Cosmic Microwave Background lensing simulation by Taylor expansion has long been considered impractical due to slow convergence, but a recent flat-sky implementation shows that a simple trick eliminates this problem, making Taylor lensing a fast and simple lensing algorithm for the flat sky. Here we generalize the method to the full sky, and study its convergence and performance relative to a commonly used numerical code, Lenspix, with extensive benchmarks of both. Compared to the flat sky case, the method takes a speed hit due to the slow speed of spherical harmonic transforms compared to fast Fourier transforms, resulting in speeds of 1/3 to 2/3 of Lenspix for similar accuracy. △ Less

Submitted 3 August, 2013; v1 submitted 2 July, 2013; originally announced July 2013.

Comments: 7 pages, 6 figures, submitted to JCAP

Journal ref: JCAP 1309 (2013) 001

CMB likelihood approximation for banded probability distributions

Authors: E. Gjerløw, K. Mikkelsen, H. K. Eriksen, K. M. Górski, G. Huey, J. B. Jewell, S. K. Næss, G. Rocha, D. S. Seljebotn, I. K. Wehus

Abstract: We investigate sets of random variables that can be arranged sequentially such that a given variable only depends conditionally on its immediate predecessor. For such sets, we show that the full joint probability distribution may be expressed exclusively in terms of uni- and bivariate marginals. Under the assumption that the CMB power spectrum likelihood only exhibits correlations within a banded… ▽ More We investigate sets of random variables that can be arranged sequentially such that a given variable only depends conditionally on its immediate predecessor. For such sets, we show that the full joint probability distribution may be expressed exclusively in terms of uni- and bivariate marginals. Under the assumption that the CMB power spectrum likelihood only exhibits correlations within a banded multipole range, Δl, we apply this expression to two outstanding problems in CMB likelihood analysis. First, we derive a statistically well-defined hybrid likelihood estimator, merging two independent (e.g., low- and high-l) likelihoods into a single expression that properly accounts for correlations between the two. Applying this expression to the WMAP likelihood, we verify that the effect of correlations on cosmological parameters in the transition region is negligible in terms of cosmological parameters for WMAP; the largest relative shift seen for any parameter is 0.06σ. However, because this may not hold for other experimental setups (e.g., for different instrumental noise properties or analysis masks), but must rather be verified on a case-by-case basis, we recommend our new hybridization scheme for future experiments for statistical self-consistency reasons. Second, we use the same expression to improve the convergence rate of the Blackwell-Rao likelihood estimator, reducing the required number of Monte Carlo samples by several orders of magnitude, and thereby extend it to high-l applications. △ Less

Submitted 1 April, 2013; originally announced April 2013.

Comments: 9 pages, 7 figures, 1 table, submitted to ApJ

Grid-based exploration of cosmological parameter space with Snake

Authors: K. Mikkelsen, S. K. Næss, H. K. Eriksen

Abstract: We present a fully parallelized grid-based parameter estimation algorithm for investigating multidimensional likelihoods called Snake, and apply it to cosmological parameter estimation. The basic idea is to map out the likelihood grid-cell by grid-cell according to decreasing likelihood, and stop when a certain threshold has been reached. This approach improves vastly on the "curse of dimensionali… ▽ More We present a fully parallelized grid-based parameter estimation algorithm for investigating multidimensional likelihoods called Snake, and apply it to cosmological parameter estimation. The basic idea is to map out the likelihood grid-cell by grid-cell according to decreasing likelihood, and stop when a certain threshold has been reached. This approach improves vastly on the "curse of dimensionality" problem plaguing standard grid-based parameter estimation simply by disregarding grid-cells with negligible likelihood. The main advantages of this method compared to standard Metropolis-Hastings MCMC methods include 1) trivial extraction of arbitrary conditional distributions; 2) direct access to Bayesian evidences; 3) better sampling of the tails of the distribution; and 4) nearly perfect parallelization scaling. The main disadvantage is, as in the case of brute-force grid-based evaluation, a dependency on the number of parameters, N_par. One of the main goals of the present paper is to determine how large N_par can be, while still maintaining reasonable computational efficiency; we find that N_par = 12$ is well within the capabilities of the method. The performance of the code is tested by comparing cosmological parameters estimated using Snake and the WMAP-7 data with those obtained using CosmoMC, the current standard code in the field. We find fully consistent results, with similar computational expenses, but shorter wall time due to the perfect parallelization scheme. △ Less

Submitted 13 November, 2012; originally announced November 2012.

Comments: 5 pages, 4 figures, submitted to ApJ

The QUIET Instrument

Authors: QUIET Collaboration, C. Bischoff, A. Brizius, I. Buder, Y. Chinone, K. Cleary, R. N. Dumoulin, A. Kusaka, R. Monsalve, S. K. Naess, L. B. Newburgh, G. Nixon, R. Reeves, K. M. Smith, K. Vanderlinde, I. K. Wehus, M. Bogdan, R. Bustos, S. E. Church, R. Davis, C. Dickinson, H. K. Eriksen, T. Gaier, J. O. Gundersen, M. Hasegawa , et al. (30 additional authors not shown)

Abstract: The Q/U Imaging ExperimenT (QUIET) is designed to measure polarization in the Cosmic Microwave Background, targeting the imprint of inflationary gravitational waves at large angular scales (~ 1 degree). Between 2008 October and 2010 December, two independent receiver arrays were deployed sequentially on a 1.4 m side-fed Dragonian telescope. The polarimeters which form the focal planes use a highly… ▽ More The Q/U Imaging ExperimenT (QUIET) is designed to measure polarization in the Cosmic Microwave Background, targeting the imprint of inflationary gravitational waves at large angular scales (~ 1 degree). Between 2008 October and 2010 December, two independent receiver arrays were deployed sequentially on a 1.4 m side-fed Dragonian telescope. The polarimeters which form the focal planes use a highly compact design based on High Electron Mobility Transistors (HEMTs) that provides simultaneous measurements of the Stokes parameters Q, U, and I in a single module. The 17-element Q-band polarimeter array, with a central frequency of 43.1 GHz, has the best sensitivity (69 uK sqrt(s)) and the lowest instrumental systematic errors ever achieved in this band, contributing to the tensor-to-scalar ratio at r < 0.1. The 84-element W-band polarimeter array has a sensitivity of 87 uK sqrt(s) at a central frequency of 94.5 GHz. It has the lowest systematic errors to date, contributing at r < 0.01. The two arrays together cover multipoles in the range l= 25-975. These are the largest HEMT-based arrays deployed to date. This article describes the design, calibration, performance of, and sources of systematic error for the instrument. △ Less

Submitted 31 July, 2012; v1 submitted 23 July, 2012; originally announced July 2012.

Journal ref: The Astrophysical Journal, Volume 768, Number 1 (2013), 28 pages

Second Season QUIET Observations: Measurements of the CMB Polarization Power Spectrum at 95 GHz

Authors: QUIET Collaboration, D. Araujo, C. Bischoff, A. Brizius, I. Buder, Y. Chinone, K. Cleary, R. N. Dumoulin, A. Kusaka, R. Monsalve, S. K. Næss, L. B. Newburgh, R. Reeves, I. K. Wehus, J. T. L. Zwart, L. Bronfman, R. Bustos, S. E. Church, C. Dickinson, H. K. Eriksen, T. Gaier, J. O. Gundersen, M. Hasegawa, M. Hazumi, K. M. Huffenberger , et al. (26 additional authors not shown)

Abstract: The Q/U Imaging ExperimenT (QUIET) has observed the cosmic microwave background (CMB) at 43 and 95GHz. The 43-GHz results have been published in QUIET Collaboration et al. (2011), and here we report the measurement of CMB polarization power spectra using the 95-GHz data. This data set comprises 5337 hours of observations recorded by an array of 84 polarized coherent receivers with a total array se… ▽ More The Q/U Imaging ExperimenT (QUIET) has observed the cosmic microwave background (CMB) at 43 and 95GHz. The 43-GHz results have been published in QUIET Collaboration et al. (2011), and here we report the measurement of CMB polarization power spectra using the 95-GHz data. This data set comprises 5337 hours of observations recorded by an array of 84 polarized coherent receivers with a total array sensitivity of 87 uK sqrt(s). Four low-foreground fields were observed, covering a total of ~1000 square degrees with an effective angular resolution of 12.8', allowing for constraints on primordial gravitational waves and high-signal-to-noise measurements of the E-modes across three acoustic peaks. The data reduction was performed using two independent analysis pipelines, one based on a pseudo-Cl (PCL) cross-correlation approach, and the other on a maximum-likelihood (ML) approach. All data selection criteria and filters were modified until a predefined set of null tests had been satisfied before inspecting any non-null power spectrum. The results derived by the two pipelines are in good agreement. We characterize the EE, EB and BB power spectra between l=25 and 975 and find that the EE spectrum is consistent with LCDM, while the BB power spectrum is consistent with zero. Based on these measurements, we constrain the tensor-to-scalar ratio to r=1.1+0.9-0.8 (r<2.8 at 95% C.L.) as derived by the ML pipeline, and r=1.2+0.9-0.8 (r<2.7 at 95% C.L.) as derived by the PCL pipeline. In one of the fields, we find a correlation with the dust component of the Planck Sky Model, though the corresponding excess power is small compared to statistical errors. Finally, we derive limits on all known systematic errors, and demonstrate that these correspond to a tensor-to-scalar ratio smaller than r=0.01, the lowest level yet reported in the literature. △ Less

Submitted 30 January, 2013; v1 submitted 20 July, 2012; originally announced July 2012.

Comments: 10 pages, 6 figures, 3 tables, submitted to ApJ, This paper should be cited as "QUIET Collaboration (2012)." v2: updated to reflect published version

Journal ref: ApJ 760, 145 (2012)

Bayesian noise estimation for non-ideal CMB experiments

Authors: I. K. Wehus, S. K. Næss, H. K. Eriksen

Abstract: We describe a Bayesian framework for estimating the time-domain noise covariance of CMB observations, typically parametrized in terms of a 1/f frequency profile. This framework is based on the Gibbs sampling algorithm, which allows for exact marginalization over nuisance parameters through conditional probability distributions. In this paper we implement support for gaps in the data streams and ma… ▽ More We describe a Bayesian framework for estimating the time-domain noise covariance of CMB observations, typically parametrized in terms of a 1/f frequency profile. This framework is based on the Gibbs sampling algorithm, which allows for exact marginalization over nuisance parameters through conditional probability distributions. In this paper we implement support for gaps in the data streams and marginalization over fixed time-domain templates, and also outline how to marginalize over confusion from CMB fluctuations, which may be important for high signal-to-noise experiments. As a by-product of the method, we obtain proper constrained realizations, which themselves can be useful for map making. To validate the algorithm, we demonstrate that the reconstructed noise parameters and corresponding uncertainties are unbiased using simulated data. The CPU time required to process a single data stream of 100 000 samples with 1000 samples removed by gaps is 3 seconds if only the maximum posterior parameters are required, and 21 seconds if one also want to obtain the corresponding uncertainties by Gibbs sampling. △ Less

Submitted 6 October, 2011; originally announced October 2011.

Comments: 8 pages, 4 figures, submitted to ApJS

Report number: Imperial/TP/2011/ikw/02

Application of the Kolmogorov-Smirnov test to CMB data: Is the universe really weakly random?

Authors: Sigurd Kirkevold Næss

Abstract: A recent application of the Kolmogorov-Smirnov test to the WMAP 7 year W-band maps claims evidence that the CMB is "weakly random", and that only 20% of the signal can be explained as a random Gaussian field. I here repeat this analysis, and in contrast to the original result find no evidence for deviation from the standard ΛCDM model. Instead, the results of the original analysis are consistent w… ▽ More A recent application of the Kolmogorov-Smirnov test to the WMAP 7 year W-band maps claims evidence that the CMB is "weakly random", and that only 20% of the signal can be explained as a random Gaussian field. I here repeat this analysis, and in contrast to the original result find no evidence for deviation from the standard ΛCDM model. Instead, the results of the original analysis are consistent with not properly taking into account the correlations of the ΛCDM power spectrum. △ Less

Submitted 25 May, 2011; originally announced May 2011.

Comments: 4 pages, 5 figures, submitted to A&A

First Season QUIET Observations: Measurements of CMB Polarization Power Spectra at 43 GHz in the Multipole Range 25 <= ell <= 475

Authors: QUIET Collaboration, C. Bischoff, A. Brizius, I. Buder, Y. Chinone, K. Cleary, R. N. Dumoulin, A. Kusaka, R. Monsalve, S. K. Næss, L. B. Newburgh, R. Reeves, K. M. Smith, I. K. Wehus, J. A. Zuntz, J. T. L. Zwart, L. Bronfman, R. Bustos, S. E. Church, C. Dickinson, H. K. Eriksen, P. G. Ferreira, T. Gaier, J. O. Gundersen, M. Hasegawa , et al. (26 additional authors not shown)

Abstract: The Q/U Imaging ExperimenT (QUIET) employs coherent receivers at 43GHz and 95GHz, operating on the Chajnantor plateau in the Atacama Desert in Chile, to measure the anisotropy in the polarization of the CMB. QUIET primarily targets the B modes from primordial gravitational waves. The combination of these frequencies gives sensitivity to foreground contributions from diffuse Galactic synchrotron ra… ▽ More The Q/U Imaging ExperimenT (QUIET) employs coherent receivers at 43GHz and 95GHz, operating on the Chajnantor plateau in the Atacama Desert in Chile, to measure the anisotropy in the polarization of the CMB. QUIET primarily targets the B modes from primordial gravitational waves. The combination of these frequencies gives sensitivity to foreground contributions from diffuse Galactic synchrotron radiation. Between 2008 October and 2010 December, >10,000hours of data were collected, first with the 19-element 43GHz array (3458hours) and then with the 90-element 95GHz array. Each array observes the same four fields, selected for low foregrounds, together covering ~1000deg^2. This paper reports initial results from the 43GHz receiver which has an array sensitivity to CMB fluctuations of 69uK sqrt(s). The data were extensively studied with a large suite of null tests before the power spectra, determined with two independent pipelines, were examined. Analysis choices, including data selection, were modified until the null tests passed. Cross correlating maps with different telescope pointings is used to eliminate a bias. This paper reports the EE, BB and EB power spectra in the multipole range ell=25-475. With the exception of the lowest multipole bin for one of the fields, where a polarized foreground, consistent with Galactic synchrotron radiation, is detected with 3sigma significance, the E-mode spectrum is consistent with the LCDM model, confirming the only previous detection of the first acoustic peak. The B-mode spectrum is consistent with zero, leading to a measurement of the tensor-to-scalar ratio of r=0.35+1.06-0.87. The combination of a new time-stream double-demodulation technique, Mizuguchi-Dragone optics, natural sky rotation, and frequent boresight rotation leads to the lowest level of systematic contamination in the B-mode power so far reported, below the level of r=0.1 △ Less

Submitted 20 July, 2012; v1 submitted 14 December, 2010; originally announced December 2010.

Comments: 19 pages, 14 figures, higher quality figures are available at http://quiet.uchicago.edu/results/index.html; Fixed a typo and corrected statistical error values used as a reference in Figure 14, showing our systematic uncertainties (unchanged) vs. multipole; Revision to ApJ accepted version, this paper should be cited as "QUIET Collaboration et al. (2011)"

Journal ref: Astrophys.J.741:111,2011

Fast and precise map-making for massively multi-detector CMB experiments

Authors: D. Sutton, J. A. Zuntz, P. G. Ferreira, M. L. Brown, H. K. Eriksen, B. R. Johnson, A. Kusaka, S. K. Naess, I. K. Wehus

Abstract: Future cosmic microwave background (CMB) polarisation experiments aim to measure an unprecedentedly small signal - the primordial gravity wave component of the polarisation field B-mode. To achieve this, they will analyse huge datasets, involving years worth of time-ordered data (TOD) from massively multi-detector focal planes. This creates the need for fast and precise methods to complement the M… ▽ More Future cosmic microwave background (CMB) polarisation experiments aim to measure an unprecedentedly small signal - the primordial gravity wave component of the polarisation field B-mode. To achieve this, they will analyse huge datasets, involving years worth of time-ordered data (TOD) from massively multi-detector focal planes. This creates the need for fast and precise methods to complement the M-L approach in analysis pipelines. In this paper, we investigate fast map-making methods as applied to long duration, massively multi-detector, ground-based experiments, in the context of the search for B-modes. We focus on two alternative map-making approaches: destriping and TOD filtering, comparing their performance on simulated multi-detector polarisation data. We have written an optimised, parallel destriping code, the DEStriping CARTographer DESCART, that is generalised for massive focal planes, including the potential effect of cross-correlated TOD 1/f noise. We also determine the scaling of computing time for destriping as applied to a simulated full-season data-set for a realistic experiment. We find that destriping can out-perform filtering in estimating both the large-scale E and B-mode angular power spectra. In particular, filtering can produce significant spurious B-mode power via EB mixing. Whilst this can be removed, it contributes to the variance of B-mode bandpower estimates at scales near the primordial B-mode peak. For the experimental configuration we simulate, this has an effect on the possible detection significance for primordial B-modes. Destriping is a viable alternative fast method to the full M-L approach that does not cause the problems associated with filtering, and is flexible enough to fit into both M-L and Monte-Carlo pseudo-Cl pipelines. △ Less

Submitted 20 May, 2010; v1 submitted 14 December, 2009; originally announced December 2009.

Comments: 16 pages, 14 figures. MNRAS accepted. Typos corrected and computing time/memory requirement orders-of-magnitude numbers in section 4 replaced by precise numbers

CMBPol Mission Concept Study: Prospects for polarized foreground removal

Authors: J. Dunkley, A. Amblard, C. Baccigalupi, M. Betoule, D. Chuss, A. Cooray, J. Delabrouille, C. Dickinson, G. Dobler, J. Dotson, H. K. Eriksen, D. Finkbeiner, D. Fixsen, P. Fosalba, A. Fraisse, C. Hirata, A. Kogut, J. Kristiansen, C. Lawrence, A. M. Magalhaes, M. A. Miville-Deschenes, S. Meyer, A. Miller, S. K. Naess, L. Page , et al. (6 additional authors not shown)

Abstract: In this report we discuss the impact of polarized foregrounds on a future CMBPol satellite mission. We review our current knowledge of Galactic polarized emission at microwave frequencies, including synchrotron and thermal dust emission. We use existing data and our understanding of the physical behavior of the sources of foreground emission to generate sky templates, and start to assess how wel… ▽ More In this report we discuss the impact of polarized foregrounds on a future CMBPol satellite mission. We review our current knowledge of Galactic polarized emission at microwave frequencies, including synchrotron and thermal dust emission. We use existing data and our understanding of the physical behavior of the sources of foreground emission to generate sky templates, and start to assess how well primordial gravitational wave signals can be separated from foreground contaminants for a CMBPol mission. At the estimated foreground minimum of ~100 GHz, the polarized foregrounds are expected to be lower than a primordial polarization signal with tensor-to-scalar ratio r=0.01, in a small patch (~1%) of the sky known to have low Galactic emission. Over 75% of the sky we expect the foreground amplitude to exceed the primordial signal by about a factor of eight at the foreground minimum and on scales of two degrees. Only on the largest scales does the polarized foreground amplitude exceed the primordial signal by a larger factor of about 20. The prospects for detecting an r=0.01 signal including degree-scale measurements appear promising, with 5 sigma_r ~0.003 forecast from multiple methods. A mission that observes a range of scales offers better prospects from the foregrounds perspective than one targeting only the lowest few multipoles. We begin to explore how optimizing the composition of frequency channels in the focal plane can maximize our ability to perform component separation, with a range of typically 40 < nu < 300 GHz preferred for ten channels. Foreground cleaning methods are already in place to tackle a CMBPol mission data set, and further investigation of the optimization and detectability of the primordial signal will be useful for mission design. △ Less

Submitted 2 November, 2009; v1 submitted 24 November, 2008; originally announced November 2008.

Comments: 42 pages, 14 figures, Foreground Removal Working Group contribution to the CMBPol Mission Concept Study, v2, matches AIP version

Journal ref: AIP Conf.Proc.1141:222-264,2009