The Q & U Bolometric Interferometer for Cosmology, QUBIC, is an innovative experiment designed to measure the polarization of the Cosmic Microwave Background and in particular the signature left therein by the inflationary expansion of the Universe. The expected signal is extremely faint, thus extreme sensitivity and systematic control are necessary in order to attempt this measurement. QUBIC addresses these requirements using an innovative approach combining the sensitivity of Transition Edge Sensor cryogenic bolometers, with the deep control of systematics characteristic of interferometers. This makes QUBIC unique with respect to others classical imagers experiments devoted to the CMB polarization. In this contribution we report a description of the QUBIC instrument including recent achievements and the demonstration of the bolometric interferometry performed in lab. QUBIC will be deployed at the observation site in Alto Chorrillos, in Argentina at the end of 2019.

QUBIC (Q \& U Bolometric Interferometer for Cosmology) is an international ground-based experiment dedicated in the measurement of the polarized fluctuations of the Cosmic Microwave Background (CMB). It is based on bolometric interferometry, an original detection technique which combine the immunity to systematic effects of an interferometer with the sensitivity of low temperature incoherent detectors. QUBIC will be deployed in Argentina, at the Alto Chorrillos mountain site near San Antonio de los Cobres, in the Salta province. The QUBIC detection chain consists in 2048 NbSi Transition Edge Sensors (TESs) cooled to 350mK.The voltage-biased TESs are read out with Time Domain Multiplexing based on Superconducting QUantum Interference Devices (SQUIDs) at 1 K and a novel SiGe Application-Specific Integrated Circuit (ASIC) at 60 K allowing to reach an unprecedented multiplexing (MUX) factor equal to 128. The QUBIC experiment is currently being characterized in the lab with a reduced number of detectors before upgrading to the full instrument. I will present the last results of this characterization phase with a focus on the detectors and readout system.

In order to better analyse the polarization of the cosmic microwave background (CMB), which is dominated by emission from our Galaxy, we need tools that can detect residual foregrounds in cleaned CMB maps. Galactic foregrounds introduce statistical anisotropy and directionality to the polarization pseudo-vectors of the CMB, which can be investigated by using the D statistic of Bunn and Scott. This statistic is rapidly computable and capable of investigating a broad range of data products for directionality. We demonstrate the application of this statistic to detecting foregrounds in polarization maps by analysing the uncleaned Planck 2018 frequency maps. For the Planck 2018 CMB maps, we find no evidence for residual foreground contamination. In order to examine the sensitivity of the D statistic, we add a varying fraction of the polarized thermal dust and synchrotron foreground maps to the CMB maps and show the percent-level foreground contamination that would be detected with 95 percent confidence. We also demonstrate application of the D statistic to another data product by analysing the gradient of the minimum-variance CMB lensing potential map (i.e., the deflection angle) for directionality. We find no excess directionality in the lensing potential map when compared to the simulations provided by the Planck Collaboration.

The statistics of extremal points in the cosmic microwave background (CMB) temperature (hot and cold spots) have been well explored in the literature, and have been used to constrain models of the early Universe. Here, we extend the study of critical points in the CMB to the set that remains after removing extrema, namely the saddle points. We perform stacks of temperature and polarization about temperature saddle points in simulations of the CMB, as well as in data from the Planck satellite. We then compute the theoretical profile of saddle-point stacks, given the underlying power spectra of the CMB. As an example of the utility of such stacks, we constrain models of cosmic birefringence, and compare the constraining power of the saddle points with that of extremal points. We find that, in the specific example of birefringence, we can place tighter constraints using saddle points in our analysis than using extrema. In fact, we find our saddle-point analysis yields close to optimal constraints, as seen by comparing to a power spectrum analysis. We, therefore, suggest that stacking on saddle points may, in general, be a useful way of testing for non-standard physics effects that change the CMB power spectra.

In this paper we describe QUBIC, an experiment that will observe the polarized microwave sky with a novel approach, which combines the sensitivity of state-of-the art bolometric detectors with the systematic effects control typical of interferometers. QUBIC unique features are the so-called "self-calibration", a technique that allows us to clean the measured data from instrumental effects, and its spectral imaging power, i.e. the ability to separate the signal in various sub-bands within each frequency band. QUBIC will observe the sky in two main frequency bands: 150 GHz and 220 GHz. A technological demonstrator is currently under testing and will be deployed in Argentina during 2019, while the final instrument is expected to be installed during 2020.

QUBIC, the QU Bolometric Interferometer for Cosmology, is a novel forthcoming instrument to measure the B-mode polarization anisotropy of the Cosmic Microwave Background. The detection of the B-mode signal will be extremely challenging; QUBIC has been designed to address this with a novel approach, namely bolometric interferometry. The receiver cryostat is exceptionally large and cools complex optical and detector stages to 40 K, 4 K, 1 K and 350 mK using two pulse tube coolers, a novel 4He sorption cooler and a double-stage 3He/4He sorption cooler. We discuss the thermal and mechanical design of the cryostat, modelling and thermal analysis, and laboratory cryogenic testing.

We summarize the radio synchrotron background workshop that took place July 19-21, 2017 at the University of Richmond. This first scientific meeting dedicated to the topic was convened because current measurements of the diffuse radio monopole reveal a surface brightness that is several times higher than can be straightforwardly explained by known Galactic and extragalactic sources and processes, rendering it by far the least well understood photon background at present. It was the conclusion of a majority of the participants that the radio monopole level is at or near that reported by the ARCADE 2 experiment and inferred from several absolutely calibrated zero level lower frequency radio measurements, and unanimously agreed that the production of this level of surface brightness, if confirmed, represents a major outstanding question in astrophysics. The workshop reached a consensus on the next priorities for investigations of the radio synchrotron background.

In this paper we describe QUBIC, an experiment that takes up the challenge posed by the detection of primordial gravitational waves with a novel approach, that combines the sensitivity of state-of-the art bolometric detectors with the systematic effects control typical of interferometers. The so-called "self-calibration" is a technique deeply rooted in the interferometric nature of the instrument and allows us to clean the measured data from instrumental effects. The first module of QUBIC is a dual band instrument (150 GHz and 220 GHz) that will be deployed in Argentina during the Fall 2018.

The next generation of CMB experiments (CMB Stage-4) will produce a Sunyaev-Zel'dovich (SZ) cluster catalog containing $\sim10^{5}$ objects, two orders of magnitudes more than currently available. In this paper, we discuss the detectability of the polarized signal generated by scattering of the CMB quadrupole on the cluster electron gas using this catalog. We discuss the possibility of using this signal to measure the relationship between cluster optical depth and mass. We find that the area of observation of S4 maximizes the signal-to-noise (S/N) on the polarized signal but that this S/N is extremely small for an individual cluster, of order 0.5% for a typical cluster in our catalog, the main source of noise being the residual primordial E-mode signal. However, we find that the signal could be detected using the full cluster catalog and that the significance of the result will increase linearly with the size of the CMB S4 telescope mirror.

In order to draw scientific conclusions from observations of cosmic microwave background (CMB) polarization, it is necessary to separate the contributions of the E and B components of the data. For data with incomplete sky coverage, there are ambiguous modes, which can be sourced by either E or B signals. Techniques exist for producing "pure" E and B maps, which are guaranteed to be free of cross-contamination, although the standard method, which involves constructing an eigenbasis, has a high computational cost. We show that such pure maps can be thought of as resulting from the application of a Wiener filter to the data. This perspective leads to far more efficient methods of producing pure maps. Moreover, by expressing the idea of purification in the general framework of Wiener filtering (i.e., maximization of a posterior probability), it leads to a variety of generalizations of the notion of pure E and B maps, e.g., accounting for noise or other contaminants in the data as well as correlations with temperature anisotropy.

QUBIC is an instrument aiming at measuring the B mode polarisation anisotropies at medium scales angular scales (30-200 multipoles). The search for the primordial CMB B-mode polarization signal is challenging, because of many difficulties: smallness of the expected signal, instrumental systematics that could possibly induce polarization leakage from the large E signal into B, brighter than anticipated polarized foregrounds (dust) reducing to zero the initial hope of finding sky regions clean enough to have a direct primordial B-modes observation. The QUBIC instrument is designed to address all aspects of this challenge with a novel kind of instrument, a Bolometric Interferometer, combining the background-limited sensitivity of Transition-Edge-Sensors and the control of systematics allowed by the observation of interference fringe patterns, while operating at two frequencies to disentangle polarized foregrounds from primordial B mode polarization. Its characteristics are described in details in this Technological Design Report.

We examine the degree to which observations of large-scale cosmic microwave background (CMB) polarization can shed light on the puzzling large-scale power modulation in maps of CMB anisotropy. We consider a phenomenological model in which the observed anomaly is caused by modulation of large-scale primordial curvature perturbations, and calculate Fisher information and error forecasts for future polarization data, constrained by the existing CMB anisotropy data. Because a significant fraction of the available information is contained in correlations with the anomalous temperature data, it is essential to account for these constraints. We also present a systematic approach to finding a set of normal modes that maximize the available information, generalizing the well-known Karhunen-Loeve transformation to take account of the constraints from the temperature data. A polarization map covering at least $\sim 60\%$ of the sky should be able to provide a $3\sigma$ detection of modulation at the level favored by the temperature data. A significant fraction of the information in such a data set is contained in the single mode that optimally encapsulates the signal due to temperature-polarization correlation.

We describe a method for forecasting errors in interferometric measurements of polarization of the cosmic microwave background (CMB) radiation, based on the use of the Fisher matrix calculated from the visibility covariance and relation matrices. In addition to noise and sample variance, the method can account for many kinds of systematic error by calculating an augmented Fisher matrix, including parameters that characterize the instrument along with the cosmological parameters to be estimated. The method is illustrated with examples of gain errors and errors in polarizer orientation. The augmented Fisher matrix approach is applicable to a much wider range of problems beyond CMB interferometry.

We present in this paper a new Bayesian semi-blind approach for foreground removal in observations of the 21-cm signal with interferometers. The technique, which we call HIEMICA (HI Expectation-Maximization Independent Component Analysis), is an extension of the Independent Component Analysis (ICA) technique developed for two-dimensional (2D) CMB maps to three-dimensional (3D) 21-cm cosmological signals measured by interferometers. This technique provides a fully Bayesian inference of power spectra and maps and separates the foregrounds from signal based on the diversity of their power spectra. Only relying on the statistical independence of the components, this approach can jointly estimate the 3D power spectrum of the 21-cm signal and, the 2D angular power spectrum and the frequency dependence of each foreground component, without any prior assumptions about foregrounds. This approach has been tested extensively by applying it to mock data from interferometric 21-cm intensity mapping observations under idealized assumptions of instrumental effects. We also discuss the impact when the noise properties are not known completely. As a first step toward solving the 21 cm power spectrum analysis problem we compare the semi-blind HIEMICA technique with the commonly used Principal Component Analysis (PCA). Under the same idealized circumstances the proposed technique provides significantly improved recovery of the power spectrum. This technique can be applied straightforwardly to all 21-cm interferometric observations, including epoch of reionization measurements, and can be extended to single-dish observations as well.

This is a brief introduction to general relativity, designed for both students and teachers of the subject. While there are many excellent expositions of general relativity, few adequately explain the geometrical meaning of the basic equation of the theory: Einstein's equation. Here we give a simple formulation of this equation in terms of the motion of freely falling test particles. We also sketch some of its consequences, and explain how the formulation given here is equivalent to the usual one in terms of tensors. Finally, we include an annotated bibliography of books, articles and websites suitable for the student of relativity.

We present results from an end-to-end simulation pipeline interferometric observations of cosmic microwave background polarization. We use both maximum-likelihood and Gibbs sampling techniques to estimate the power spectrum. In addition, we use Gibbs sampling for image reconstruction from interferometric visibilities. The results indicate the level to which various systematic errors (e.g., pointing errors, gain errors, beam shape errors, cross- polarization) must be controlled in order to successfully detect and characterize primordial B modes as well as other scientific goals. In addition, we show that Gibbs sampling is an effective method of image reconstruction for interferometric data in other astrophysical contexts.

We present a novel, general-purpose method for deconvolving and denoising images from gridded radio interferometric visibilities using Bayesian inference based on a Gaussian process model. The method automatically takes into account incomplete coverage of the uv-plane, signal mode coupling due to the primary beam, and noise mode coupling due to uv sampling. Our method uses Gibbs sampling to efficiently explore the full posterior distribution of the underlying signal image given the data. We use a set of widely diverse mock images with a realistic interferometer setup and level of noise to assess the method. Compared to results from a proxy for point source- based CLEAN method we find that in terms of RMS error and signal-to-noise ratio our approach performs better than traditional deconvolution techniques, regardless of the structure of the source image in our test suite. Our implementation scales as O(np log np), provides full statistical and uncertainty information of the reconstructed image, requires no supervision, and provides a robust, consistent framework for incorporating noise and parameter marginalizations and foreground removal.

QUBIC is a unique instrument that crosses the barriers between classical imaging architectures and interferometry taking advantage from both for high sensitivity and systematics mitigation. The scientific target is the detection of the primordial gravitational waves imprint on the Cosmic Microwave Background which are the proof of inflation, holy grail of modern cosmology. In this paper, we show the latest advances in the development of the architecture and the sub-systems of the first module of this instrument to be deployed in Dome Charlie Concordia base - Antarctica in 2015.

We investigate the impact of instrumental systematic errors in interferometric measurements of the cosmic microwave background (CMB) temperature and polarization power spectra. We simulate interferometric CMB observations to generate mock visibilities and estimate power spectra using the statistically optimal maximum likelihood technique. We define a quadratic error measure to determine allowable levels of systematic error that do not induce power spectrum errors beyond a given tolerance. As an example, in this study we focus on differential pointing errors. The effects of other systematics can be simulated by this pipeline in a straightforward manner. We find that, in order to accurately recover the underlying B-modes for r=0.01 at 28<l<384, Gaussian-distributed pointing errors must be controlled to 0.7^\circ rms for an interferometer with an antenna configuration similar to QUBIC, in agreement with analytical estimates. Only the statistical uncertainty for 28<l<88 would be changed at ~10% level. With the same instrumental configuration, we find the pointing errors would slightly bias the 2-\sigma upper limit of the tensor-to-scalar ratio r by ~10%. We also show that the impact of pointing errors on the TB and EB measurements is negligibly small.

The detection of the primordial $B$-mode spectrum of the polarized cosmic microwave background (CMB) signal may provide a probe of inflation. However, observation of such a faint signal requires excellent control of systematic errors. Interferometry proves to be a promising approach for overcoming such a challenge. In this paper we present a complete simulation pipeline of interferometric observations of CMB polarization, including systematic errors. We employ two different methods for obtaining the power spectra from mock data produced by simulated observations: the maximum likelihood method and the method of Gibbs sampling. We show that the results from both methods are consistent with each other, as well as, within a factor of 6, with analytical estimates. Several categories of systematic errors are considered: instrumental errors, consisting of antenna gain and antenna coupling errors, and beam errors, consisting of antenna pointing errors, beam cross-polarization and beam shape (and size) errors. In order to recover the tensor-to-scalar ratio, $r$, within a 10% tolerance level, which ensures the experiment is sensitive enough to detect the $B$-signal at $r=0.01$ in the multipole range $28 < \ell < 384$, we find that, for a QUBIC-like experiment, Gaussian-distributed systematic errors must be controlled with precisions of $|g_{rms}| = 0.1$ for antenna gain, $|\epsilon_{rms}| = 5 \times 10^{-4}$ for antenna coupling, $\delta_{rms} \approx 0.7^\circ$ for pointing, $\zeta_{rms} \approx 0.7^\circ$ for beam shape, and $\mu_{rms} = 5 \times 10^{-4}$ for beam cross-polarization.

Detection of B-mode polarization of the cosmic microwave background (CMB) radiation is one of the frontiers of observational cosmology. Because they are an order of magnitude fainter than E-modes, it is quite a challenge to detect B-modes. Having more manageable systematics, interferometers prove to have a substantial advantage over imagers in detecting such faint signals. Here, we present a method for Bayesian inference of power spectra and signal reconstruction from interferometric data of the CMB polarization signal by using the technique of Gibbs sampling. We demonstrate the validity of the method in the flat-sky approximation for a simulation of an interferometric observation on a finite patch with incomplete uv-plane coverage, a finite beam size and a realistic noise model. With a computational complexity of O(n^{3/2}), n being the data size, Gibbs sampling provides an efficient method for analyzing upcoming cosmology observations.

Separation of the B component of a cosmic microwave background (CMB) polarization map from the much larger E component is an essential step in CMB polarimetry. For a map with incomplete sky coverage, this separation is necessarily hampered by the presence of "ambiguous" modes which could be either E or B modes. I present an efficient pixel-space algorithm for removing the ambiguous modes and separating the map into "pure" E and B components. The method, which works for arbitrary geometries, does not involve generating a complete basis of such modes and scales the cube of the number of pixels on the boundary of the map.

One of the major challenges of modern cosmology is the detection of B-mode polarization anisotropies in the CMB. These originate from tensor fluctuations of the metric produced during the inflationary phase. Their detection would therefore constitute a major step towards understanding the primordial Universe. The expected level of these anisotropies is however so small that it requires a new generation of instruments with high sensitivity and extremely good control of systematic effects. We propose the QUBIC instrument based on the novel concept of bolometric interferometry, bringing together the sensitivity advantages of bolometric detectors with the systematics effects advantages of interferometry. Methods: The instrument will directly observe the sky through an array of entry horns whose signals will be combined together using an optical combiner. The whole set-up is located inside a cryostat. Polarization modulation will be achieved using a rotating half-wave plate and interference fringes will be imaged on two focal planes (separated by a polarizing grid) tiled with bolometers. We show that QUBIC can be considered as a synthetic imager, exactly similar to a usual imager but with a synthesized beam formed by the array of entry horns. Scanning the sky provides an additional modulation of the signal and improve the sky coverage shape. The usual techniques of map-making and power spectrum estimation can then be applied. We show that the sensitivity of such an instrument is comparable with that of an imager with the same number of horns. We anticipate a low level of beam-related systematics thanks to the fact that the synthesized beam is determined by the location of the primary horns. Other systematics should be under good control thanks to an autocalibration technique, specific to our concept, that will permit the accurate determination of most of the systematics parameters.

Several anomalies appear to be present in the large-angle cosmic microwave background (CMB) anisotropy maps of WMAP, including the alignment of large-scale multipoles. Models in which isotropy is spontaneously broken (e.g., by a scalar field) have been proposed as explanations for these anomalies, as have models in which a preferred direction is imposed during inflation. We examine models inspired by these, in which isotropy is broken by a multiplicative factor with dipole and/or quadrupole terms. We evaluate the evidence provided by the multipole alignment using a Bayesian framework, finding that the evidence in favor of the model is generally weak. We also compute approximate changes in estimated cosmological parameters in the broken-isotropy models. Only the overall normalization of the power spectrum is modified significantly.

Several claims have been made of anomalies in the large-angle properties of the cosmic microwave background anisotropy as measured by WMAP. In most cases, the statistical significance of these anomalies is hard or even impossible to assess, due to the fact that the statistics used to quantify the anomalies were chosen a posteriori. On the other hand, the possibility of detecting new physics on the largest observable scales is so exciting that, in my opinion, it is worthwhile to examine the claims carefully. I will focus on three particular claims: the lack of large-angle power, the north-south power asymmetry, and multipole alignments. In all cases, the problem of a posteriori statistics can best be solved by finding a new data set that probes similar physical scales to the large-angle CMB. This is a difficult task, but there are some possible routes to achieving it.

Bolometric Interferometry is a technology currently under development that will be first dedicated to the detection of B-mode polarization fluctuations in the Cosmic Microwave Background. A bolometric interferometer will have to take advantage of the wide spectral detection band of its bolometers in order to be competitive with imaging experiments. A crucial concern is that interferometers are presumed to be importantly affected by a spoiling effect known as bandwidth smearing. In this paper, we investigate how the bandwidth modifies the work principle of a bolometric interferometer and how it affects its sensitivity to the CMB angular power spectra. We obtain analytical expressions for the broadband visibilities measured by broadband heterodyne and bolometric interferometers. We investigate how the visibilities must be reconstructed in a broadband bolometric interferometer and show that this critically depends on hardware properties of the modulation phase shifters. Using an angular power spectrum estimator accounting for the bandwidth, we finally calculate the sensitivity of a broadband bolometric interferometer. A numerical simulation has been performed and confirms the analytical results. We conclude (i) that broadband bolometric interferometers allow broadband visibilities to be reconstructed whatever the kind of phase shifters used and (ii) that for dedicated B-mode bolometric interferometers, the sensitivity loss due to bandwidth smearing is quite acceptable, even for wideband instruments (a factor 2 loss for a typical 20% bandwidth experiment).

A common belief about big-bang cosmology is that the cosmological redshift cannot be properly viewed as a Doppler shift (that is, as evidence for a recession velocity), but must be viewed in terms of the stretching of space. We argue that, contrary to this view, the most natural interpretation of the redshift is as a Doppler shift, or rather as the accumulation of many infinitesimal Doppler shifts. The stretching-of-space interpretation obscures a central idea of relativity, namely that it is always valid to choose a coordinate system that is locally Minkowskian. We show that an observed frequency shift in any spacetime can be interpreted either as a kinematic (Doppler) shift or a gravitational shift by imagining a suitable family of observers along the photon's path. In the context of the expanding universe the kinematic interpretation corresponds to a family of comoving observers and hence is more natural.

How did the universe evolve? The fine angular scale (l>1000) temperature and polarization anisotropies in the CMB are a Rosetta stone for understanding the evolution of the universe. Through detailed measurements one may address everything from the physics of the birth of the universe to the history of star formation and the process by which galaxies formed. One may in addition track the evolution of the dark energy and discover the net neutrino mass. We are at the dawn of a new era in which hundreds of square degrees of sky can be mapped with arcminute resolution and sensitivities measured in microKelvin. Acquiring these data requires the use of special purpose telescopes such as the Atacama Cosmology Telescope (ACT), located in Chile, and the South Pole Telescope (SPT). These new telescopes are outfitted with a new generation of custom mm-wave kilo-pixel arrays. Additional instruments are in the planning stages.

Modern cosmology has sharpened questions posed for millennia about the origin of our cosmic habitat. The age-old questions have been transformed into two pressing issues primed for attack in the coming decade: How did the Universe begin? and What physical laws govern the Universe at the highest energies? The clearest window onto these questions is the pattern of polarization in the Cosmic Microwave Background (CMB), which is uniquely sensitive to primordial gravity waves. A detection of the special pattern produced by gravity waves would be not only an unprecedented discovery, but also a direct probe of physics at the earliest observable instants of our Universe. Experiments which map CMB polarization over the coming decade will lead us on our first steps towards answering these age-old questions.

In future microwave background polarization experiments, particularly those that aim to characterize the B component, careful attention will have to be paid to the mixing of E and B components due to finite sky coverage and pixelization. Any polarization map can be decomposed into "pure E," "pure B," and "ambiguous" components. In practice, since the B component is expected to be much weaker than the E component, nearly all of the recoverable $B$ information is contained in the pure B component. The amount of B information lost to ambiguous modes can be estimated in simple ways from the survey geometry and pixelization. Separation of any given map into pure and ambiguous components can be done by finding a basis of pure and ambiguous modes, but it is often more efficient to "purify" the map directly in real space by solving a certain differential equation to find the ambiguous component. This method may be useful in conjunction with power spectrum estimation techniques such as the pseudo-Cl method.

Cosmic microwave background (CMB) polarimetry has the potential to provide revolutionary advances in cosmology. Future experiments to detect the very weak B mode signal in CMB polarization maps will require unprecedented sensitivity and control of systematic errors. Bolometric interferometry may provide a way to achieve these goals. In a bolometric interferometer (or other adding interferometer), phase shift sequences are applied to the inputs in order to recover the visibilities. Noise is minimized when the phase shift sequences corresponding to all visibilities are orthogonal. We present a systematic method for finding sequences that produce this orthogonality, approximately minimizing both the length of the time sequence and the number of discrete phase shift values required. When some baselines are geometrically equivalent, we can choose sequences that read out those baselines simultaneously, which has been shown to improve signal to noise ratio.

Several anomalies appear to be present in the large-angle cosmic microwave background (CMB) anisotropy maps of WMAP. One of these is a lack of large-scale power. Because the data otherwise match standard models extremely well, it is natural to consider perturbations of the standard model as possible explanations. We show that, as long as the source of the perturbation is statistically independent of the source of the primary CMB anisotropy, no such model can explain this large-scale power deficit. On the contrary, any such perturbation always reduces the probability of obtaining any given low value of large-scale power. We rigorously prove this result when the lack of large-scale power is quantified with a quadratic statistic, such as the quadrupole moment. When a statistic based on the integrated square of the correlation function is used instead, we present strong numerical evidence in support of the result. The result applies to models in which the geometry of spacetime is perturbed (e.g., an ellipsoidal Universe) as well as explanations involving local contaminants, undiagnosed foregrounds, or systematic errors. Because the large-scale power deficit is arguably the most significant of the observed anomalies, explanations that worsen this discrepancy should be regarded with great skepticism, even if they help in explaining other anomalies such as multipole alignments.

Polarization is the next frontier of CMB analysis, but its signal is dominated over much of the sky by foregrounds which must be carefully removed. To determine the efficacy of this cleaning it is necessary to have sensitive tests for residual foreground contamination in polarization sky maps. The dominant Galactic foregrounds introduce a large-scale anisotropy on to the sky, so it makes sense to use a statistic sensitive to overall directionality for this purpose. Here we adapt the rapidly computable D statistic of Bunn and Scott to polarization data, and demonstrate its utility as a foreground monitor by applying it to the low resolution WMAP 3-yr sky maps. With a thorough simulation of the maps' noise properties, we find no evidence for contamination in the foreground cleaned sky maps.

We search for viable f(R) theories of gravity, making use of the equivalence between such theories and scalar-tensor gravity. We find that models can be made consistent with solar system constraints either by giving the scalar a high mass or by exploiting the so-called chameleon effect. However, in both cases, it appears likely that any late-time cosmic acceleration will be observationally indistinguishable from acceleration caused by a cosmological constant. We also explore further observational constraints from, e.g., big bang nucleosynthesis and inflation.

Cosmic microwave background (CMB) polarization observations will require superb control of systematic errors in order to achieve their full scientific potential, particularly in the case of attempts to detect the B modes that may provide a window on inflation. Interferometry may be a promising way to achieve these goals. This paper presents a formalism for characterizing the effects of a variety of systematic errors on interferometric CMB polarization observations, with particular emphasis on estimates of the B-mode power spectrum. The most severe errors are those that couple the temperature anisotropy signal to polarization; such errors include cross-talk within detectors, misalignment of polarizers, and cross-polarization. In a B mode experiment, the next most serious category of errors are those that mix E and B modes, such as gain fluctuations, pointing errors, and beam shape errors. The paper also indicates which sources of error may cause circular polarization (e.g., from foregrounds) to contaminate the cosmologically interesting linear polarization channels, and conversely whether monitoring of the circular polarization channels may yield useful information about the errors themselves. For all the sources of error considered, estimates of the level of control that will be required for both E and B mode experiments are provided. Both experiments that interfere linear polarizations and those that interfere circular polarizations are considered. The fact that circular experiments simultaneously measure both linear polarization Stokes parameters in each baseline mitigates some sources of error.

Measurements of cosmic microwave background (CMB) anisotropies by interferometers offer several advantages over single-dish observations. The formalism for analyzing interferometer CMB data is well developed in the flat-sky approximation, valid for small fields of view. As the area of sky is increased to obtain finer spectral resolution, this approximation needs to be relaxed. We extend the formalism for CMB interferometry, including both temperature and polarization, to mosaics of observations covering arbitrarily large areas of the sky, with each individual pointing lying within the flat-sky approximation. We present a method for computing the correlation between visibilities with arbitrary pointing centers and baselines and illustrate the effects of sky curvature on the l-space resolution that can be obtained from a mosaic.

Scattering of cosmic microwave background (CMB) radiation in galaxy clusters induces a polarization signal proportional to the CMB quadrupole anisotropy at the cluster's location and look-back time. A survey of such remote quadrupole measurements provides information about large-scale cosmological perturbations. This paper presents a formalism for calculating the correlation function of remote quadrupole measurements in spherical harmonic space. The number of independent modes probed by both single-redshift and volume-limited surveys is presented, along with the length scales probed by these modes. In a remote quadrupole survey sparsely covering a large area of sky, the largest-scale modes probe the same length scales as the quadrupole but with much narrower Fourier-space window functions. The largest-scale modes are significantly correlated with the local CMB, but even when this correlation is projected out the largest remaining modes probe gigaparsec scales (comparable to the CMB at l=2-10) with narrow window functions. These modes may provide insight into the possible anomalies in the large-scale CMB anisotropy. At fixed redshift, the data from such a survey form an E-type spin-2 field on the sphere to a good approximation; the near-absence of B modes will provide a valuable check on systematic errors. A survey of only a few low-redshift clusters allows an independent reconstruction of the five coefficients of the local CMB quadrupole, providing a test for contamination in the WMAP quadrupole. The formalism presented here is also useful for analyzing smaller-scale surveys to probe the late integrated Sachs-Wolfe effect and hence the properties of dark energy

In a microwave background polarization map that covers only part of the sky, it is impossible to separate the E and B components perfectly. This difficulty in general makes it more difficult to detect the B component in a data set. Any polarization map can be separated in a unique way into "pure E," "pure B," and "ambiguous" components. Power that resides in the pure E (B) component is guaranteed to be produced by E (B) modes, but there is no way to tell whether the ambiguous component comes from E or B modes. A polarization map can be separated into the three components either by finding an orthonormal basis for each component, or directly in real space by using Green functions or other methods.

[NOTE: Previous versions of this paper (both on astro-ph and published in Phys. Rev. D) contain results that are in error. The power spectra C_l were normalized incorrectly by a factor of 2 pi. All observing times in detector-years in those versions are too large by a factor of 2 pi. The main place these numbers appear is on the vertical axes of Figures 4 and 5. Note that because all calculations were based on the same power spectra, all conclusions pertaining to comparisons of different techniques remain unchanged. This error has been corrected in the present version of the paper. An erratum is being sent to Phys. Rev. D. I apologize for the error.] Using a Fisher-matrix formalism, we calculate the required sensitivities and observing times for an experiment to measure the amplitudes of both E and B components as a function of sky coverage, taking full account of the fact that the two components cannot be perfectly separated in an incomplete sky map. We also present a simple approximation scheme that accounts for mixing of E and B components in computing predicted errors in the E-component power spectrum amplitude. In an experiment with small sky coverage, mixing of the two components increases the difficulty of detecting the subdominant B component by a factor of two or more in observing time; however, for larger survey sizes the effect of mixing is less pronounced. Surprisingly, mixing of E and B components can enhance the detectability of the E component by increasing the effective number of independent modes that probe this component

Separation of the E and B components of a microwave background polarization map or a weak lensing map is an essential step in extracting science from it, but when the map covers only part of the sky and/or is pixelized, this decomposition cannot be done perfectly. We present a method for decomposing an arbitrary sky map into a sum of three orthogonal components that we term ``pure E'', ``pure B'' and ``ambiguous''. This method is useful both for providing intuition for experimental design and for analyzing data sets in practice. We show how to find orthonormal bases for all three components in terms of bilaplacian eigenfunctions. The number of ambiguous modes is proportional to the length of the map boundary so fairly round maps are preferred. For real-world data analysis, we present a simple matrix eigenvalue method for calculating nearly pure E and B modes in pixelized maps. We find that the dominant source of leakage between E and B is aliasing of small-scale power caused by the pixelization. This problem can be eliminated by heavily oversampling the map, but is exacerbated by the fact that the E power spectrum is expected to be much larger than the B power spectrum and extremely blue. We found that a factor of 2 to 3 more pixels are needed in a polarization map to achieve the same level of contamination by aliased power than in a temperature map.

We use a frequentist statistical approach to set confidence intervals on the values of cosmological parameters using the MAXIMA-1 and COBE measurements of the angular power spectrum of the cosmic microwave background. We define a $\Delta \chi^{2}$ statistic, simulate the measurements of MAXIMA-1 and COBE, determine the probability distribution of the statistic, and use it and the data to set confidence intervals on several cosmological parameters. We compare the frequentist confidence intervals to Bayesian credible regions. The frequentist and Bayesian approaches give best estimates for the parameters that agree within 15%, and confidence interval-widths that agree within 30%. The results also suggest that a frequentist analysis gives slightly broader confidence intervals than a Bayesian analysis. The frequentist analysis gives values of \Omega=0.89{+0.26\atop -0.19}, \Omega_{\rm B}h^2=0.026{+0.020\atop -0.011} and n=1.02{+0.31\atop -0.10}, and the Bayesian analysis gives values of \Omega=0.98{+0.14\atop -0.19}, \Omega_{\rm B}h^2=0.0.029{+0.015\atop-0.010}, and $n=1.18{+0.10\atop -0.23}$, all at the 95% confidence level.

Large patterns could exist on the microwave sky as a result of various non-standard possibilities for the large-scale Universe -- rotation or shear, non-trivial topology, and single topological defects are specific examples. All-sky (or nearly all-sky) CMB data sets allow us, uniquely, to constrain such exotica, and it is therefore worthwhile to explore a wide range of statistical tests. We describe one such statistic here, which is based on determining gradients and is useful for assessing the level of `preferred directionality' or `stripiness' in the map. This method is more general than other techniques for picking out specific patterns on the sky, and it also has the advantage of being easily calculable for the mega-pixel maps which will soon be available. For the purposes of illustration, we apply this statistic to the four-year COBE DMR data. For future CMB maps we expect this to be a useful statistical test of the large scale structure of the Universe. In principle, the same statistic could also be applied to sky maps at other wavelengths, to CMB polarisation maps, and to catalogues of discrete objects. It may also be useful as a means of checking for residual directionality (e.g., from Galactic or ecliptic signals) in maps.

A recent paper (Nodland and Ralston, PRL 78, 3043, astro-ph/9704196) claims to have detected evidence for birefringence in the propagation of radio waves across cosmological distances. In order to assess the statistical significance of their results, the authors analyze simulated data sets, finding a stronger correlation in the real data than in the simulations. Unfortunately, the procedure used for generating the simulated data sets is based on an incorrect null hypothesis. Furthermore, the correct null hypothesis would lead to a stronger correlation in the simulated data sets, weakening the case for birefringence. We conclude that the paper's analysis does not provide statistically significant evidence of birefringence.

We investigate the behaviour of light rays in an open universe with a partly periodic horn topology. The geodesics can be solved exactly and the periodic topology creates characteristic new effects in temperature maps of the microwave background. Large flat regions appear in COBE sky maps even when multiple images of astronomical sources are unobservable.

We present an analysis of the four-year data from the COBE DMR experiment. We use a Karhunen-Loeve expansion of the pixel data to calculate the normalization and goodness-of-fit of a range of models of structure formation. This technique produces unbiased normalization estimates and is capable of determining the normalization of any particular model with a statistical uncertainty of 7%. We present a parameterization of the normalization and likelihood function which summarizes our results for a wide range of models. We use the COBE normalization to compute small-scale fluctuation amplitudes for a variety of models and discuss the implications of these results for theories of large-scale structure.

More than a dozen papers analyzing the COBE data have now appeared. We review the different techniques and compare them to a ``brute force" likelihood analysis where we invert the full 4038 x 4038 Galaxy-cut pixel covariance matrix. This method is optimal in the sense of producing minimal error bars, and is a useful reference point for comparing other analysis techniques. Our maximum-likelihood estimate of the spectral index and normalization are n=1.15 (0.95) and Q=18.2 (21.3) micro-Kelvin including (excluding) the quadrupole. Marginalizing over the normalization C_9, we obtain n=1.10 +/- 0.29 (n=0.90 +/- 0.32). When we compare these results with those of the various techniques that involve a linear ``compression'' of the data, we find that the latter are all consistent with the brute-force analysis and have error bars that are nearly as small as the minimal error bars. We therefore conclude that the data compressions involved in these techniques do indeed retain most of the useful cosmological information.

We supply fitting formulae enabling the normalization of slow-roll inflation models to the four-year COBE data. We fully include the effect of the gravitational wave modes, including the predicted relation of the amplitude of these modes to that of the density perturbations. We provide the normalization of the matter power spectrum, which can be directly used for large-scale structure studies. The normalization for tilted spectra is a special case. We also provide fitting functions for the inflationary energy scale of COBE-normalized models and discuss the validity of approximating the spectra by power-laws. In an Appendix, we extend our analysis to include models with a cosmological constant, both with and without gravitational waves.

We review the physical processes that are thought to produce anisotropy in the cosmic microwave background, focusing primarily (but not exclusively) on the effects of acoustic waves in the early Universe. We attempt throughout to supply an intuitive, physical picture of the key ideas and to elucidate the ways in which the predicted anisotropy depends on cosmological parameters such as Omega and h. The second half of these lectures is devoted to a discussion of microwave background data analysis techniques, with an emphasis on the analysis of the COBE DMR data. In particular, the Karhunen-Loeve method of data compression is described in detail.

Large-scale cosmic microwave background anisotropies in homogeneous, globally anisotropic cosmologies are investigated. We perform a statistical analysis in which the four-year data from the Cosmic Background Explorer satellite is searched for the specific anisotropy patterns predicted by these models and thereby set definitive upper limits on the amount of shear, $(\sigma/H)_0$ and vorticity, $(\omega/H)_0$, which are orders of magnitude stronger than previous constraints. We comment on how these results might impact our understanding of primordial global anisotropy.

With the advent of the COBE detection of fluctuations in the Cosmic Microwave Background radiation, the study of inhomogeneous cosmology has entered a new phase. It is now possible to accurately normalize fluctuations on the largest observable scales, in the linear regime. In this paper we present a model-independent method of normalizing theories to the full COBE data. This technique allows an extremely wide range of theories to be accurately normalized to COBE in a very simple and fast way. We give the best fitting normalization and relative peak likelihoods for a range of spectral shapes, and discuss the normalization for several popular theories. Additionally we present both Bayesian and frequentist measures of the goodness of fit of a representative range of theories to the COBE data.

We apply a Wiener filter to the two-year COBE DMR data. The resulting sky map has significantly reduced noise levels compared to the raw data: the most prominent hot and cold spots are significant at the 4-sigma level. Furthermore, the entire covariance matrix of the errors in the filtered sky map is known, and it is therefore possible to make constrained realizations of the microwave sky with the correct a posteriori probability distribution. The filtered DMR sky map is used to make predictions for the Tenerife experiment. Two prominent features are predicted in a region of the sky not yet analyzed by the Tenerife group. The presence of these features is a robust prediction of the standard cosmological paradigm; if these features are not observed, some of our fundamental assumptions must be incorrect.

Motivated by recent studies of the one-bubble inflationary scenario, simple open cold dark matter models are tested for consistency with cosmological observations. The initial perturbation spectrum is derived by solving for the evolution of fluctuations in an open inflationary stage. A likelihood analysis is performed for the Cosmic Microwave Background anisotropies using the two-year COBE DMR data and considering models based on both the Bunch-Davies and conformal vacua. Having normalized the perturbation spectrum to fit the COBE data, we reconsider the validity of the open model from the view point of cosmic structure formation. Open models may be severely constrained by the COBE likelihood analysis. In particular, small values of $\Omega_0$ are ruled out in the Bunch-Davies case: we find that $\Omega_0\ge 0.34$ at 95\% confidence for this model.

We consider Bayesian constraints on standard isocurvature baryon models from the slope and normalization of the anisotropy power spectrum detected by the COBE DMR experiment in their two year maps. In conjunction with either the amplitude of matter fluctuations $\sigma_8$ {\it or} its slope, all open models are ruled out at greater than 95\% confidence, whereas cosmological constant $\Lambda$ dominated models are constrained to be highly ionized. By including the COBE FIRAS 95\% confidence upper limit on spectral distortion, we further reduce the available parameter space for $\Lambda$ models by excluding these highly ionized models. These constraints define a single class of barely viable standard models which makes definite and testable predictions for degree scale anisotropies and large scale structure.

We use a Maximum Entropy technique to reconstruct a map of the microwave sky near the star Gamma Ursae Minoris, based on data from flights 2, 3 and 4 of the Millimeter-wave Anisotropy eXperiment (MAX).

We compare the two-year COBE DMR sky maps with the predictions of cosmological-constant cold dark matter models. Using a Bayesian analysis, we find that the most likely value of the cosmological constant in such a model is Lambda = 0. The data set an upper limit on Lambda of 0.78 (0.85) at 90% confidence, and 0.86 (0.92) at 95% confidence with (without) the quadrupole anisotropy.

The COBE detection of CMB anisotropies provides the best way of fixing the amplitude of fluctuations on the largest scales. This normalization is usually given for an n=1 spectrum, including only the anisotropy caused by the Sachs- Wolfe effect. This is certainly not a good approximation for a model containing any reasonable amount of baryonic matter. In fact, even tilted S-W spectra are not a good fit to models like CDM. Here we normalize standard CDM (sCDM) to the 2-year COBE data, and quote the best amplitude in terms of the conventionally used measures of power. We also give normalizations for some specific variants of this standard model, and we indicate how the normalization depends on the assumed values of n, Omega_B and H_0. For sCDM we find <Q>=19.9\pm1.5uK, corresponding to sigma_8=1.34\pm0.10, with the normalization at large scales being B=(8.16\pm1.04)\times10^5 (Mpc/h)^4, and other numbers given in the Table. The measured rms temperature fluctuation smoothed on 10deg is a little low relative to this normalization. This is mainly due to the low quadrupole in the data: when the quadrupole is removed, the measured value of sigma(10) is quite consistent with the best-fitting <Q>. The use of <Q> should be preferred over sigma(10), when its value can be determined for a particular theory, since it makes full use of the data.

We derive an optimal linear filter to suppress the noise from the COBE DMR sky maps for a given power spectrum. We then apply the filter to the first-year DMR data, after removing pixels within $20^\circ$ of the Galactic plane from the data. The filtered data have uncertainties 12 times smaller than the noise level of the raw data. We use the formalism of constrained realizations of Gaussian random fields to assess the uncertainty in the filtered sky maps. In addition to improving the signal-to-noise ratio of the map as a whole, these techniques allow us to recover some information about the CMB anisotropy in the missing Galactic plane region. From these maps we are able to determine which hot and cold spots in the data are statistically significant, and which may have been produced by noise. In addition, the filtered maps can be used for comparison with other experiments on similar angular scales.

A formalism is presented that allows cosmological experiments to be tested for consistency, and allows a simple frequentist interpretation of the resulting significance levels. As an example of an application, this formalism is used to place constraints on bulk flows of galaxies using the results of the microwave background anisotropy experiments COBE and SP91, and a few simplifying approximations about the experimental window functions. It is found that if taken at face value, with the quoted errors, the recent detection by Lauer and Postman of a bulk flow of 689 km/s on scales of 150$h^{-1}$Mpc is inconsistent with SP91 at a 95\% confidence level within the framework of a Cold Dark Matter (CDM) model. The same consistency test is also used to place constraints that are completely model-independent, in the sense that they hold for any power spectrum whatsoever --- the only assumption being that the random fields are Gaussian. It is shown that the resulting infinite-dimensional optimization problem reduces to a set of coupled non-linear equations that can readily be solved numerically. Applying this technique to the above-mentioned example, we find that the Lauer and Postman result is inconsistent with SP91 even if no assumptions whatsoever are made about the power spectrum.

We analyze the Gamma Ursae Minoris (GUM) and Mu Pegasii (MuP) scans of the Millimeter-wave Anisotropy eXperiment in the context of cold dark matter (CDM) models of structure formation, paying particular attention to the two-dimensional nature of the GUM scan. If all of the structure in the (foreground subtracted) data is attributed to cosmic microwave background anisotropy, then there is a detection in each scan. For a standard CDM model, the amplitudes of the signals are individually compatible with the COsmic Background Explorer measurement, but marginally inconsistent with each other.

We present results on the consistency of standard cold dark matter (CDM) models with both the COBE normalization and the data from the 4th channel of the 9-point South Pole scan. We find that CDM models are consistent with both experiments, a conclusion which is at odds with some other analyses. This is partly due to our careful treatment of the temperature autocorrelation function, but also derives from a realization that the statistical conclusions depend strongly on assumptions about the prior distribution for the normalization $Q$.