How to make maps from CMB data without losing information

The next generation of CMB experiments can measure cosmological parameters with unprecedented accuracy - in principle. To achieve this in practice when faced with such gigantic data sets, elaborate data analysis methods are needed to make it computationally feasible. An important step in the data pipeline is to make a map, which typically reduces the size of the data set my orders of magnitude. We compare ten map-making methods, and find that for the Gaussian case, both the method used by the COBE DMR team and various variants of Wiener filtering are optimal in the sense that the map retains all cosmological information that was present in the time-ordered data (TOD). Specifically, one obtains just as small error bars on cosmological parameters when estimating them from the map as one could have obtained by estimating them directly from the TOD. The method of simply averaging the observations of each pixel (for total-power detectors), on the contrary, is found to generally destroy information, as does the maximum entropy method and most other non-linear map-making techniques. Since it is also numerically feasible, the COBE method is the natural choice for large data sets. Other lossless (e.g. Wiener-filtered) maps can then be computed directly from the COBE method map.Comment: Minor revisions to match published version. 12 pages, with 1 figure included. Color figure and links at http://www.sns.ias.edu/~max/mapmaking.html (faster from the US), from http://www.mpa-garching.mpg.de/~max/mapmaking.html (faster from Europe) or from [email protected]

Removing point sources from CMB maps

For high-precision cosmic microwave background (CMB) experiments, contamination from extragalactic point sources is a major concern. It is therefore useful to be able to detect and discard point source contaminated pixels using the map itself. We show that the sensitivity with which this can be done can often be greatly improved (by factors between 2.5 and 18 for the upcoming Planck mission) by a customized hi-pass filtering that suppresses fluctuations due to CMB and diffuse galactic foregrounds. This means that point source contamination will not severely degrade the cleanest Planck channels unless current source count estimates are off by more than an order of magnitude. A catalog of around 40,000 far infra-red sources at 857 GHz may be a useful by-product of Planck.Comment: 4 pages, with 2 figures included. Minor revisions to match accepted version. Color figure and links at http://www.sns.ias.edu/~max/cleaning.html (faster from the US), from http://www.mpa-garching.mpg.de/~max/cleaning.html (faster from Europe) or from [email protected], and Angelica's foreground links at http://www.sns.ias.edu/~angelica/foreground.htm

Cosmic Microwave Background Anisotropy Observing Strategy Assessment

I develop a method for assessing the ability of an instrument, coupled with an observing strategy, to measure the angular power spectrum of the cosmic microwave background (CMB). It allows for efficient calculation of expected parameter uncertainties. Related to this method is a means of graphically presenting, via the ``eigenmode window function'', the sensitivity of an observation to different regions of the spectrum, which is a generalization of the traditional practice of presenting the trace of the window function. I apply these techniques to a balloon-borne bolometric instrument to be flown this summer (MSAM2). I find that a smoothly scanning secondary is better than a chopping one and that, in this case, a very simple analytic formula provides a good (40\% or better) approximation to expected power spectrum uncertainties.Comment: Substantial revisions, LaTeX 15 pages including 3 figure

Non-Gaussianity in Two-Field Inflation

We derive semi-analytic formulae for the local bispectrum and trispectrum in general two-field inflation and provide a simple geometric recipe for building observationally allowed models with observable non-Gaussianity. We use the \delta N formalism and the transfer function formalism to express the bispectrum almost entirely in terms of model-independent physical quantities. Similarly, we calculate the trispectrum and show that the trispectrum parameter \tau NL can be expressed entirely in terms of spectral observables, which provides a new consistency relation unique to two-field inflation. We show that in order to generate observably large non-Gaussianity during inflation, the sourcing of curvature modes by isocurvature modes must be extremely sensitive to the initial conditions, and that the amount of sourcing must be moderate in order to avoid excessive fine-tuning. Under some minimal assumptions, we argue that the first condition is satisfied only when neighboring trajectories through the two-dimensional field space diverge during inflation. Geometrically, this means that the inflaton must roll along a ridge in the potential V for some time during inflation and that its trajectory must turn slightly (but not too sharply) in field space. Therefore, it follows that two-field scenarios with attractor solutions necessarily produce small non-Gaussianity. This explains why it has been so difficult to achieve large non-Gaussianity in two-field inflation, and why it has only been achieved in a narrow class of models like hybrid inflation and certain product potentials where the potential and/or the initial conditions are fine-tuned. Some of our conclusions generalize qualitatively to general multi-field inflation.Comment: Discussion improved, gNL formula and extra figure included, typos corrected, references added. 18 pages, 2 figure

Experimental neutrino physics

The current experimental status of neutrino physics is reviewed. It contains the evidences for a non-vanishing neutrino rest mass from neutrino oscillation searches. In addition an outlook is given on determining the various mixing matrix elements and mass differences more precisely with new experiments. Of special interest is the value of the mixing angle \theta_{13} determining the possibility of detecting leptonic CP violation in the future. The prospect for absolute mass measurements using beta and double beta decay as well as cosmological observations is presented.Comment: 11 pages, 11 figures, Inv. talk presented at the DPF 2004 meeting of the APS, Riverside, Aug. 200

Learning from observations of the microwave background at small angular scales

In this paper, we focus our attention on the following question: How well can we recover the power spectrum of the cosmic microwave background from the maps of a given experiment?. Each experiment is described by a a pixelization scale, a beam size, a noise level and a sky coverage. We use accurate numerical simulations of the microwave sky and a cold dark matter model for structure formation in the universe. Angular scales smaller than those of previous simulations are included. The spectrum obtained from the simulated maps is appropriately compared with the theoretical one. Relative deviations between these spectra are estimated. Various contributions to these deviations are analyzed. The method used for spectra comparisons is discussed.Comment: 15 pages (LATEX), 2 postcript figures, accepted in Ap

The power spectrum of galaxies in the 2dF 100k redshift survey

We compute the real-space power spectrum and the redshift-space distortions of galaxies in the 2dF 100k galaxy redshift survey using pseudo-Karhunen-Loeve eigenmodes and the stochastic bias formalism. Our results agree well with those published by the 2dFGRS team, and have the added advantage of producing easy-to-interpret uncorrelated minimum-variance measurements of the galaxy-galaxy, galaxy-velocity and velocity-velocity power spectra in 27 k-bands, with narrow and well-behaved window functions in the range 0.01h/Mpc < k < 0.8h/Mpc. We find no significant detection of baryonic wiggles, although our results are consistent with a standard flat Omega_Lambda=0.7 ``concordance'' model and previous tantalizing hints of baryonic oscillations. We measure the galaxy-matter correlation coefficient r > 0.4 and the redshift-distortion parameter beta=0.49+/-0.16 for r=1 (beta=0.47+/- 0.16 without finger-of-god compression). Since this is an apparent-magnitude limited sample, luminosity-dependent bias may cause a slight red-tilt in the power spectum. A battery of systematic error tests indicate that the survey is not only impressive in size, but also unusually clean, free of systematic errors at the level to which our tests are sensitive. Our measurements and window functions are available at http://www.hep.upenn.edu/~max/2df.html together with the survey mask, radial selection function and uniform subsample of the survey that we have constructed.Comment: Replaced to match accepted MNRAS version, with new radial/angular systematics plot and sigma8 typo corrected. High-res figures, power spectra, windows and our uniform galaxy subsample with mask at http://www.hep.upenn.edu/~max/2df.html or from [email protected]. 26 journal pages, 28 fig

An iterative destriping technique for diffuse background polarization data

We describe a simple but effective iterative procedure specifically designed to destripe Q and U Stokes parameter data as those collected by the SPOrt experiment onboard the International Space Station (ISS). The method is general enough to be useful for other experiments, both in polarization and total intensity. The only requirement for the algorithm to work properly is that the receiver knee frequency must be lower than the signal modulation frequency, corresponding in our case to the ISS orbit period. Detailed performances of the technique are presented in the context of the SPOrt experiment, both in terms of added rms noise and residual correlated noise.Comment: Accepted for publication by A&A (8 pages, 6 figures

Global 21cm signal experiments: a designer's guide

[Abridged] The spatially averaged global spectrum of the redshifted 21cm line has generated much experimental interest, for it is potentially a direct probe of the Epoch of Reionization and the Dark Ages. Since the cosmological signal here has a purely spectral signature, most proposed experiments have little angular sensitivity. This is worrisome because with only spectra, the global 21cm signal can be difficult to distinguish from foregrounds such as Galactic synchrotron radiation, as both are spectrally smooth and the latter is orders of magnitude brighter. We establish a mathematical framework for global signal data analysis in a way that removes foregrounds optimally, complementing spectra with angular information. We explore various experimental design trade-offs, and find that 1) with spectral-only methods, it is impossible to mitigate errors that arise from uncertainties in foreground modeling; 2) foreground contamination can be significantly reduced for experiments with fine angular resolution; 3) most of the statistical significance in a positive detection during the Dark Ages comes from a characteristic high-redshift trough in the 21cm brightness temperature; and 4) Measurement errors decrease more rapidly with integration time for instruments with fine angular resolution. We show that if observations and algorithms are optimized based on these findings, an instrument with a 5 degree beam can achieve highly significant detections (greater than 5-sigma) of even extended (high Delta-z) reionization scenarios after integrating for 500 hrs. This is in contrast to instruments without angular resolution, which cannot detect gradual reionization. Abrupt ionization histories can be detected at the level of 10-100's of sigma. The expected errors are also low during the Dark Ages, with a 25-sigma detection of the expected cosmological signal after only 100 hrs of integration.Comment: 34 pages, 30 figures. Replaced (v2) to match accepted PRD version (minor pedagogical additions to text; methods, results, and conclusions unchanged). Fixed two typos (v3); text, results, conclusions etc. completely unchange