Source author record

Max Tegmark

Max Tegmark appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

53works
20topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

53 published item(s)

preprint2022arXiv

Machine-learning hidden symmetries

We present an automated method for finding hidden symmetries, defined as symmetries that become manifest only in a new coordinate system that must be discovered. Its core idea is to quantify asymmetry as violation of certain partial differential equations, and to numerically minimize such violation over the space of all invertible transformations, parametrized as invertible neural networks. For example, our method rediscovers the famous Gullstrand-Painleve metric that manifests hidden translational symmetry in the Schwarzschild metric of non-rotating black holes, as well as Hamiltonicity, modularity and other simplifying traits not traditionally viewed as symmetries.

preprint2022arXiv

Pareto-optimal clustering with the primal deterministic information bottleneck

At the heart of both lossy compression and clustering is a trade-off between the fidelity and size of the learned representation. Our goal is to map out and study the Pareto frontier that quantifies this trade-off. We focus on the optimization of the Deterministic Information Bottleneck (DIB) objective over the space of hard clusterings. To this end, we introduce the primal DIB problem, which we show results in a much richer frontier than its previously studied Lagrangian relaxation when optimized over discrete search spaces. We present an algorithm for mapping out the Pareto frontier of the primal DIB trade-off that is also applicable to other two-objective clustering problems. We study general properties of the Pareto frontier, and we give both analytic and numerical evidence for logarithmic sparsity of the frontier in general. We provide evidence that our algorithm has polynomial scaling despite the super-exponential search space, and additionally, we propose a modification to the algorithm that can be used where sampling noise is expected to be significant. Finally, we use our algorithm to map the DIB frontier of three different tasks: compressing the English alphabet, extracting informative color classes from natural images, and compressing a group theory-inspired dataset, revealing interesting features of frontier, and demonstrating how the structure of the frontier can be used for model selection with a focus on points previously hidden by the cloak of the convex hull.

preprint2020arXiv

AI Feynman: a Physics-Inspired Method for Symbolic Regression

A core challenge for both physics and artificial intellicence (AI) is symbolic regression: finding a symbolic expression that matches data from an unknown function. Although this problem is likely to be NP-hard in principle, functions of practical interest often exhibit symmetries, separability, compositionality and other simplifying properties. In this spirit, we develop a recursive multidimensional symbolic regression algorithm that combines neural network fitting with a suite of physics-inspired techniques. We apply it to 100 equations from the Feynman Lectures on Physics, and it discovers all of them, while previous publicly available software cracks only 71; for a more difficult test set, we improve the state of the art success rate from 15% to 90%.

preprint2020arXiv

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

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

preprint2020arXiv

Pareto-optimal data compression for binary classification tasks

The goal of lossy data compression is to reduce the storage cost of a data set $X$ while retaining as much information as possible about something ($Y$) that you care about. For example, what aspects of an image $X$ contain the most information about whether it depicts a cat? Mathematically, this corresponds to finding a mapping $X\to Z\equiv f(X)$ that maximizes the mutual information $I(Z,Y)$ while the entropy $H(Z)$ is kept below some fixed threshold. We present a method for mapping out the Pareto frontier for classification tasks, reflecting the tradeoff between retained entropy and class information. We first show how a random variable $X$ (an image, say) drawn from a class $Y\in\{1,...,n\}$ can be distilled into a vector $W=f(X)\in \mathbb{R}^{n-1}$ losslessly, so that $I(W,Y)=I(X,Y)$; for example, for a binary classification task of cats and dogs, each image $X$ is mapped into a single real number $W$ retaining all information that helps distinguish cats from dogs. For the $n=2$ case of binary classification, we then show how $W$ can be further compressed into a discrete variable $Z=g_β(W)\in\{1,...,m_β\}$ by binning $W$ into $m_β$ bins, in such a way that varying the parameter $β$ sweeps out the full Pareto frontier, solving a generalization of the Discrete Information Bottleneck (DIB) problem. We argue that the most interesting points on this frontier are "corners" maximizing $I(Z,Y)$ for a fixed number of bins $m=2,3...$ which can be conveniently be found without multiobjective optimization. We apply this method to the CIFAR-10, MNIST and Fashion-MNIST datasets, illustrating how it can be interpreted as an information-theoretically optimal image clustering algorithm.

preprint2019arXiv

The role of artificial intelligence in achieving the Sustainable Development Goals

The emergence of artificial intelligence (AI) and its progressively wider impact on many sectors across the society requires an assessment of its effect on sustainable development. Here we analyze published evidence of positive or negative impacts of AI on the achievement of each of the 17 goals and 169 targets of the 2030 Agenda for Sustainable Development. We find that AI can support the achievement of 128 targets across all SDGs, but it may also inhibit 58 targets. Notably, AI enables new technologies that improve efficiency and productivity, but it may also lead to increased inequalities among and within countries, thus hindering the achievement of the 2030 Agenda. The fast development of AI needs to be supported by appropriate policy and regulation. Otherwise, it would lead to gaps in transparency, accountability, safety and ethical standards of AI-based technology, which could be detrimental towards the development and sustainable use of AI. Finally, there is a lack of research assessing the medium- and long-term impacts of AI. It is therefore essential to reinforce the global debate regarding the use of AI and to develop the necessary regulatory insight and oversight for AI-based technologies.

preprint2016arXiv

Research Priorities for Robust and Beneficial Artificial Intelligence

Success in the quest for artificial intelligence has the potential to bring unprecedented benefits to humanity, and it is therefore worthwhile to investigate how to maximize these benefits while avoiding potential pitfalls. This article gives numerous examples (which should by no means be construed as an exhaustive list) of such worthwhile research aimed at ensuring that AI remains robust and beneficial.

preprint2016arXiv

The Hydrogen Epoch of Reionization Array Dish II: Characterization of Spectral Structure with Electromagnetic Simulations and its science Implications

We use time-domain electromagnetic simulations to determine the spectral characteristics of the Hydrogen Epoch of Reionization Arrays (HERA) antenna. These simulations are part of a multi-faceted campaign to determine the effectiveness of the dish's design for obtaining a detection of redshifted 21 cm emission from the epoch of reionization. Our simulations show the existence of reflections between HERA's suspended feed and its parabolic dish reflector that fall below -40 dB at 150 ns and, for reasonable impedance matches, have a negligible impact on HERA's ability to constrain EoR parameters. It follows that despite the reflections they introduce, dishes are effective for increasing the sensitivity of EoR experiments at relatively low cost. We find that electromagnetic resonances in the HERA feed's cylindrical skirt, which is intended to reduce cross coupling and beam ellipticity, introduces significant power at large delays ($-40$ dB at 200 ns) which can lead to some loss of measurable Fourier modes and a modest reduction in sensitivity. Even in the presence of this structure, we find that the spectral response of the antenna is sufficiently smooth for delay filtering to contain foreground emission at line-of-sight wave numbers below $k_\parallel \lesssim 0.2$ $h$Mpc$^{-1}$, in the region where the current PAPER experiment operates. Incorporating these results into a Fisher Matrix analysis, we find that the spectral structure observed in our simulations has only a small effect on the tight constraints HERA can achieve on parameters associated with the astrophysics of reionization.

preprint2015arXiv

Consciousness as a State of Matter

We examine the hypothesis that consciousness can be understood as a state of matter, "perceptronium", with distinctive information processing abilities. We explore five basic principles that may distinguish conscious matter from other physical systems such as solids, liquids and gases: the information, integration, independence, dynamics and utility principles. If such principles can identify conscious entities, then they can help solve the quantum factorization problem: why do conscious observers like us perceive the particular Hilbert space factorization corresponding to classical space (rather than Fourier space, say), and more generally, why do we perceive the world around us as a dynamic hierarchy of objects that are strongly integrated and relatively independent? Tensor factorization of matrices is found to play a central role, and our technical results include a theorem about Hamiltonian separability (defined using Hilbert-Schmidt superoperators) being maximized in the energy eigenbasis. Our approach generalizes Giulio Tononi's integrated information framework for neural-network-based consciousness to arbitrary quantum systems, and we find interesting links to error-correcting codes, condensed matter criticality, and the Quantum Darwinism program, as well as an interesting connection between the emergence of consciousness and the emergence of time.

preprint2015arXiv

Empirical Covariance Modeling for 21 cm Power Spectrum Estimation: A Method Demonstration and New Limits from Early Murchison Widefield Array 128-Tile Data

The separation of the faint cosmological background signal from bright astrophysical foregrounds remains one of the most daunting challenges of mapping the high-redshift intergalactic medium with the redshifted 21 cm line of neutral hydrogen. Advances in mapping and modeling of diffuse and point source foregrounds have improved subtraction accuracy, but no subtraction scheme is perfect. Precisely quantifying the errors and error correlations due to missubtracted foregrounds allows for both the rigorous analysis of the 21 cm power spectrum and for the maximal isolation of the "EoR window" from foreground contamination. We present a method to infer the covariance of foreground residuals from the data itself in contrast to previous attempts at a priori modeling. We demonstrate our method by setting limits on the power spectrum using a 3 h integration from the 128-tile Murchison Widefield Array. Observing between 167 and 198 MHz, we find at 95% confidence a best limit of Delta^2(k) < 3.7 x 10^4 mK^2 at comoving scale k = 0.18 hMpc^-1 and at z = 6.8, consistent with existing limits.

preprint2015arXiv

Mapmaking for Precision 21 cm Cosmology

In order to study the "Cosmic Dawn" and the Epoch of Reionization with 21 cm tomography, we need to statistically separate the cosmological signal from foregrounds known to be orders of magnitude brighter. Over the last few years, we have learned much about the role our telescopes play in creating a putatively foreground-free region called the "EoR window." In this work, we examine how an interferometer's effects can be taken into account in a way that allows for the rigorous estimation of 21 cm power spectra from interferometric maps while mitigating foreground contamination and thus increasing sensitivity. This requires a precise understanding of the statistical relationship between the maps we make and the underlying true sky. While some of these calculations would be computationally infeasible if performed exactly, we explore several well-controlled approximations that make mapmaking and the calculation of map statistics much faster, especially for compact and highly redundant interferometers designed specifically for 21 cm cosmology. We demonstrate the utility of these methods and the parametrized trade-offs between accuracy and speed using one such telescope, the upcoming Hydrogen Epoch of Reionization Array, as a case study.

preprint2014arXiv

Friendly Artificial Intelligence: the Physics Challenge

Relentless progress in artificial intelligence (AI) is increasingly raising concerns that machines will replace humans on the job market, and perhaps altogether. Eliezer Yudkowski and others have explored the possibility that a promising future for humankind could be guaranteed by a superintelligent "Friendly AI", designed to safeguard humanity and its values. I argue that, from a physics perspective where everything is simply an arrangement of elementary particles, this might be even harder than it appears. Indeed, it may require thinking rigorously about the meaning of life: What is "meaning" in a particle arrangement? What is "life"? What is the ultimate ethical imperative, i.e., how should we strive to rearrange the particles of our Universe and shape its future? If we fail to answer the last question rigorously, this future is unlikely to contain humans.

preprint2014arXiv

MITEoR: A Scalable Interferometer for Precision 21 cm Cosmology

We report on the MIT Epoch of Reionization (MITEoR) experiment, a pathfinder low-frequency radio interferometer whose goal is to test technologies that improve the calibration precision and reduce the cost of the high-sensitivity 3D mapping required for 21 cm cosmology. MITEoR accomplishes this by using massive baseline redundancy, which enables both automated precision calibration and correlator cost reduction. We demonstrate and quantify the power and robustness of redundancy for scalability and precision. We find that the calibration parameters precisely describe the effect of the instrument upon our measurements, allowing us to form a model that is consistent with $χ^2$ per degree of freedom < 1.2 for as much as 80% of the observations. We use these results to develop an optimal estimator of calibration parameters using Wiener filtering, and explore the question of how often and how finely in frequency visibilities must be reliably measured to solve for calibration coefficients. The success of MITEoR with its 64 dual-polarization elements bodes well for the more ambitious Hydrogen Epoch of Reionization Array (HERA) project and other next-generation instruments, which would incorporate many identical or similar technologies.

preprint2014arXiv

Overcoming real-world obstacles in 21 cm power spectrum estimation: A method demonstration and results from early Murchison Widefield Array data

We present techniques for bridging the gap between idealized inverse covariance weighted quadratic estimation of 21 cm power spectra and the real-world challenges presented universally by interferometric observation. By carefully evaluating various estimators and adapting our techniques for large but incomplete data sets, we develop a robust power spectrum estimation framework that preserves the so-called "EoR window" and keeps track of estimator errors and covariances. We apply our method to observations from the 32-tile prototype of the Murchinson Widefield Array to demonstrate the importance of a judicious analysis technique. Lastly, we apply our method to investigate the dependence of the clean EoR window on frequency--especially the frequency dependence of the so-called "wedge" feature--and establish upper limits on the power spectrum from z = 6.2 to z = 11.7. Our lowest limit is Delta(k) < 0.3 Kelvin at 95% confidence at a comoving scale k = 0.046 Mpc^-1 and z = 9.5.

preprint2014arXiv

Sharpening the Second Law of Thermodynamics with the Quantum Bayes' Theorem

We prove a generalization of the classic Groenewold-Lindblad entropy inequality, combining decoherence and the quantum Bayes theorem into a simple unified picture where decoherence increases entropy while observation decreases it. This provides a rigorous quantum-mechanical version of the second law of thermodynamics, governing how the entropy of a system (the entropy of its density matrix, partial-traced over the environment and conditioned on what is known) evolves under general decoherence and observation. The powerful tool of spectral majorization enables both simple alternative proofs of the classic Lindblad and Holevo inequalities without using strong subadditivity, and also novel inequalities for decoherence and observation that hold not only for von Neumann entropy, but also for arbitrary concave entropies.

preprint2014arXiv

What Next-Generation 21 cm Power Spectrum Measurements Can Teach Us About the Epoch of Reionization

A number of experiments are currently working towards a measurement of the 21 cm signal from the Epoch of Reionization. Whether or not these experiments deliver a detection of cosmological emission, their limited sensitivity will prevent them from providing detailed information about the astrophysics of reionization. In this work, we consider what types of measurements will be enabled by a next-generation of larger 21 cm EoR telescopes. To calculate the type of constraints that will be possible with such arrays, we use simple models for the instrument, foreground emission, and the reionization history. We focus primarily on an instrument modeled after the $\sim 0.1~\rm{km}^2$ collecting area Hydrogen Epoch of Reionization Array (HERA) concept design, and parameterize the uncertainties with regard to foreground emission by considering different limits to the recently described "wedge" footprint in k-space. Uncertainties in the reionization history are accounted for using a series of simulations which vary the ionizing efficiency and minimum virial temperature of the galaxies responsible for reionization, as well as the mean free path of ionizing photons through the IGM. Given various combinations of models, we consider the significance of the possible power spectrum detections, the ability to trace the power spectrum evolution versus redshift, the detectability of salient power spectrum features, and the achievable level of quantitative constraints on astrophysical parameters. Ultimately, we find that $0.1~\rm{km}^2$ of collecting area is enough to ensure a very high significance ($\gtrsim30σ$) detection of the reionization power spectrum in even the most pessimistic scenarios. This sensitivity should allow for meaningful constraints on the reionization history and astrophysical parameters, especially if foreground subtraction techniques can be improved and successfully implemented.

preprint2013arXiv

A Fast Method for Power Spectrum and Foreground Analysis for 21 cm Cosmology

We develop and demonstrate an acceleration of the Liu & Tegmark quadratic estimator formalism for inverse variance foreground subtraction and power spectrum estimation in 21 cm tomography from O(N^3) to O(N log N), where N is the number of voxels of data. This technique makes feasible the megavoxel scale analysis necessary for current and upcoming radio interferometers by making only moderately restrictive assumptions about foreground models and survey geometry. We exploit iterative and Monte Carlo techniques and the symmetries of the foreground covariance matrices to quickly estimate the 21 cm brightness temperature power spectrum, P(k_parallel, k_perpendicular), the Fisher information matrix, the error bars, the window functions, and the bias. We also extend the Liu & Tegmark foreground model to include bright point sources with known positions in a way that scales as O[(N log N)(N point sources)] < O(N^5/3). As a first application of our method, we forecast error bars and window functions for the upcoming 128-tile deployment of the Murchinson Widefield Array, showing that 1000 hours of observation should prove sufficiently sensitive to detect the power spectrum signal from the Epoch of Reionization.

preprint2013arXiv

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.

preprint2013arXiv

Mapping our Universe in 3D with MITEoR

Mapping our universe in 3D by imaging the redshifted 21 cm line from neutral hydrogen has the potential to overtake the cosmic microwave background as our most powerful cosmological probe, because it can map a much larger volume of our Universe, shedding new light on the epoch of reionization, inflation, dark matter, dark energy, and neutrino masses. We report on MITEoR, a pathfinder low-frequency radio interferometer whose goal is to test technologies that greatly reduce the cost of such 3D mapping for a given sensitivity. MITEoR accomplishes this by using massive baseline redundancy both to enable automated precision calibration and to cut the correlator cost scaling from N^2 to NlogN, where N is the number of antennas. The success of MITEoR with its 64 dual-polarization elements bodes well for the more ambitious HERA project, which would incorporate many identical or similar technologies using an order of magnitude more antennas, each with dramatically larger collecting area.

preprint2012arXiv

Born in an Infinite Universe: a Cosmological Interpretation of Quantum Mechanics

We study the quantum measurement problem in the context of an infinite, statistically uniform space, as could be generated by eternal inflation. It has recently been argued that when identical copies of a quantum measurement system exist, the standard projection operators and Born rule method for calculating probabilities must be supplemented by estimates of relative frequencies of observers. We argue that an infinite space actually renders the Born rule redundant, by physically realizing all outcomes of a quantum measurement in different regions, with relative frequencies given by the square of the wave function amplitudes. Our formal argument hinges on properties of what we term the quantum confusion operator, which projects onto the Hilbert subspace where the Born rule fails, and we comment on its relation to the oft-discussed quantum frequency operator. This analysis unifies the classical and quantum levels of parallel universes that have been discussed in the literature, and has implications for several issues in quantum measurement theory. It also shows how, even for a single measurement, probabilities may be interpreted as relative frequencies in unitary (Everettian) quantum mechanics. We also argue that after discarding a zero-norm part of the wavefunction, the remainder consists of a superposition of indistinguishable terms, so that arguably "collapse" of the wavefunction is irrelevant, and the "many worlds" of Everett's interpretation are unified into one. Finally, the analysis suggests a "cosmological interpretation" of quantum theory in which the wave function describes the actual spatial collection of identical quantum systems, and quantum uncertainty is attributable to the observer's inability to self-locate in this collection.

preprint2012arXiv

How unitary cosmology generalizes thermodynamics and solves the inflationary entropy problem

We analyze cosmology assuming unitary quantum mechanics, using a tripartite partition into system, observer and environment degrees of freedom. This generalizes the second law of thermodynamics to "The system's entropy can't decrease unless it interacts with the observer, and it can't increase unless it interacts with the environment." The former follows from the quantum Bayes Theorem we derive. We show that because of the long-range entanglement created by cosmological inflation, the cosmic entropy decreases exponentially rather than linearly with the number of bits of information observed, so that a given observer can reduce entropy by much more than the amount of information her brain can store. Indeed, we argue that as long as inflation has occurred in a non-negligible fraction of the volume, almost all sentient observers will find themselves in a post-inflationary low-entropy Hubble volume, and we humans have no reason to be surprised that we do so as well, which solves the so-called inflationary entropy problem. An arguably worse problem for unitary cosmology involves gamma-ray-burst constraints on the "Big Snap", a fourth cosmic doomsday scenario alongside the "Big Crunch", "Big Chill" and "Big Rip", where an increasingly granular nature of expanding space modifies our life-supporting laws of physics. Our tripartite framework also clarifies when it is valid to make the popular quantum gravity approximation that the Einstein tensor equals the quantum expectation value of the stress-energy tensor, and how problems with recent attempts to explain dark energy as gravitational backreaction from super-horizon scale fluctuations can be understood as a failure of this approximation.

preprint2012arXiv

Testing Multi-Field Inflation: A Geometric Approach

We develop an approach for linking the power spectra, bispectrum, and trispectrum to the geometric and kinematical features of multifield inflationary Lagrangians. Our geometric approach can also be useful in determining when a complicated multifield model can be well approximated by a model with one, two, or a handful of fields. To arrive at these results, we focus on the mode interactions in the kinematical basis, starting with the case of no sourcing and showing that there is a series of mode conservation laws analogous to the conservation law for the adiabatic mode in single-field inflation. We then treat the special case of a quadratic potential with canonical kinetic terms, showing that it produces a series of mode sourcing relations identical in form to that for the adiabatic mode. We build on this result to show that the mode sourcing relations for general multifield inflation are extension of this special case but contain higher-order covariant derivatives of the potential and corrections from the field metric. In parallel, we show how these interactions depend on the geometry of the inflationary Lagrangian and on the kinematics of the associated field trajectory. Finally, we consider how the mode interactions and effective number of fields active during inflation are reflected in the spectra and introduce a multifield consistency relation, as well as a multifield observable that can potentially distinguish two-field scenarios from scenarios involving three or more effective fields.

preprint2011arXiv

A Method for 21cm Power Spectrum Estimation in the Presence of Foregrounds

21cm tomography promises to be a powerful tool for estimating cosmological parameters, constraining the epoch of reionization, and probing the so-called dark ages. However, realizing this promise will require the extraction of a cosmological power spectrum from beneath overwhelmingly large sources of foreground contamination. In this paper, we develop a unified matrix-based framework for foreground subtraction and power spectrum estimation, which allows us to quantify the errors and biases that arise in the power spectrum as a result of foreground subtraction. We find that existing line-of-sight foreground subtraction proposals can lead to substantial mode-mixing as well as residual noise and foreground biases, whereas our proposed inverse variance foreground subtraction eliminates noise and foreground biases, gives smaller error bars, and produces less correlated measurements of the power spectrum. We also numerically confirm the intuitive belief in the literature that 21cm foreground subtraction is best done using frequency rather than angular information.

preprint2011arXiv

Galaxy Clustering in the Completed SDSS Redshift Survey: The Dependence on Color and Luminosity

We measure the luminosity and color dependence of galaxy clustering in the SDSS DR7 main galaxy sample, focusing on the projected correlation function w_p(r_p) of volume-limited samples. We interpret our measurements using halo occupation distribution (HOD) modeling assuming a Lambda-CDM cosmology. The amplitude of w_p(r_p) grows slowly with luminosity for L < L_* and increases sharply at higher luminosities, with bias factor b(>L)=1.06+0.23(L/L_*)^{1.12}. At fixed luminosity, redder galaxies have a stronger and steeper w_p(r_p), a trend that runs steadily from the bluest galaxies to the reddest galaxies. The individual luminosity trends for the red and blue galaxy populations are strikingly different. Blue galaxies show a slow but steady increase of w_p(r_p) with luminosity, at all scales. The large-scale clustering of red galaxies shows little luminosity dependence until a sharp increase at L > 4L_*, but the lowest luminosity red galaxies (0.04-0.25 L_*) show very strong clustering on scales r_p < 2 Mpc/h. Most of the observed trends can be naturally understood within the LCDM+HOD framework. The growth of w_p(r_p) with luminosity reflects an overall shift in the halo mass scale, in particular an increase in the minimum host halo mass Mmin. The mass at which a halo has, on average, one satellite galaxy brighter than L is M_1 ~ 17 Mmin(L) over most of the luminosity range. The growth and steepening of w_p(r_p) for redder galaxies reflects the increasing fraction of galaxies that are satellite systems in high mass halos instead of central systems in low mass halos, a trend that is especially marked at low luminosities. Our extensive measurements, provided in tabular form, will allow detailed tests of theoretical models of galaxy formation, a firm grounding of semi-empirical models of the galaxy population, and new cosmological tests.

preprint2011arXiv

How well can we measure and understand foregrounds with 21 cm experiments?

Before it becomes a sensitive probe of the Epoch of Reionization, the Dark Ages, and fundamental physics, 21 cm tomography must successfully contend with the issue of foreground contamination. Broadband foreground sources are expected to be roughly four orders of magnitude larger than any cosmological signals, so precise foreground models will be necessary. Such foreground models often contain a large number of parameters, reflecting the complicated physics that governs foreground sources. In this paper, we concentrate on spectral modeling (neglecting, for instance, bright point source removal from spatial maps) and show that 21 cm tomography experiments will likely not be able to measure these parameters without large degeneracies, simply because the foreground spectra are so featureless and generic. However, we show that this is also an advantage, because it means that the foregrounds can be characterized to a high degree of accuracy once a small number of parameters (likely three or four, depending on one's instrumental specifications) are measured. This provides a simple understanding for why 21 cm foreground subtraction schemes are able to remove most of the contaminants by suppressing just a small handful of simple spectral forms. In addition, this suggests that the foreground modeling process should be relatively simple and will likely not be an impediment to the foreground subtraction schemes that are necessary for a successful 21 cm tomography experiment.

preprint2010arXiv

Axion Cosmology and the Energy Scale of Inflation

We survey observational constraints on the parameter space of inflation and axions and map out two allowed windows: the classic window and the inflationary anthropic window. The cosmology of the latter is particularly interesting; inflationary axion cosmology predicts the existence of isocurvature fluctuations in the CMB, with an amplitude that grows with both the energy scale of inflation and the fraction of dark matter in axions. Statistical arguments favor a substantial value for the latter, and so current bounds on isocurvature fluctuations imply tight constraints on inflation. For example, an axion Peccei-Quinn scale of 10^16 GeV excludes any inflation model with energy scale > 3.8*10^14 GeV (r > 2*10^(-9)) at 95% confidence, and so implies negligible gravitational waves from inflation, but suggests appreciable isocurvature fluctuations.

preprint2010arXiv

Likely values of the Higgs vev

We make an estimate of the likelihood function for the Higgs vacuum expectation value by imposing anthropic constraints on the existence of atoms while allowing the other parameters of the Standard Model to also be variable. We argue that the most important extra ingredients are the Yukawa couplings, and for the intrinsic distribution of Yukawa couplings we use the scale invariant distribution which is favored phenomenologically. The result is successful phenomenologically, favoring values close to the observed vev. We also discuss modifications that could change these conclusions. Our work supports the hypothesis that the anthropic constraints could be the origin of the small value of the Higgs vev.

preprint2010arXiv

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 δ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 τ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.

preprint2010arXiv

Omniscopes: Large Area Telescope Arrays with only N log N Computational Cost

We show that the class of antenna layouts for telescope arrays allowing cheap analysis hardware (with correlator cost scaling as N log N rather than N^2 with the number of antennas N) is encouragingly large, including not only previously discussed rectangular grids but also arbitrary hierarchies of such grids, with arbitrary rotations and shears at each level. We show that all correlations for such a 2D array with an n-level hierarchy can be efficiently computed via a Fast Fourier Transform in not 2 but 2n dimensions. This can allow major correlator cost reductions for science applications requiring exquisite sensitivity at widely separated angular scales, for example 21cm tomography (where short baselines are needed to probe the cosmological signal and long baselines are needed for point source removal), helping enable future 21cm experiments with thousands or millions of cheap dipole-like antennas. Such hierarchical grids combine the angular resolution advantage of traditional array layouts with the cost advantage of a rectangular Fast Fourier Transform Telescope. We also describe an algorithm for how a subclass of hierarchical arrays can efficiently use rotation synthesis to produce global sky maps with minimal noise and a well-characterized synthesized beam.

preprint2010arXiv

Precision Calibration of Radio Interferometers Using Redundant Baselines

Growing interest in 21 cm tomography has led to the design and construction of broadband radio interferometers with low noise, moderate angular resolution, high spectral resolution, and wide fields of view. With characteristics somewhat different from traditional radio instruments, these interferometers may require new calibration techniques in order to reach their design sensitivities. Self-calibration or redundant calibration techniques that allow an instrument to be calibrated off complicated sky emission structures are ideal. In particular, the large number of redundant baselines possessed by these new instruments makes redundant calibration an especially attractive option. In this paper, we explore the errors and biases in existing redundant calibration schemes through simulations, and show how statistical biases can be eliminated. We also develop a general calibration formalism that includes both redundant baseline methods and basic point source calibration methods as special cases, and show how slight deviations from perfect redundancy and coplanarity can be taken into account.

preprint2010arXiv

Solving the Corner-Turning Problem for Large Interferometers

The so-called corner turning problem is a major bottleneck for radio telescopes with large numbers of antennas. The problem is essentially that of rapidly transposing a matrix that is too large to store on one single device; in radio interferometry, it occurs because data from each antenna needs to be routed to an array of processors that will each handle a limited portion of the data (a frequency range, say) but requires input from each antenna. We present a low-cost solution allowing the correlator to transpose its data in real time, without contending for bandwidth, via a butterfly network requiring neither additional RAM memory nor expensive general-purpose switching hardware. We discuss possible implementations of this using FPGA, CMOS, analog logic and optical technology, and conclude that the corner turner cost can be small even for upcoming massive radio arrays.

preprint2010arXiv

Testing Two-Field Inflation

We derive semi-analytic formulae for the power spectra of two-field inflation assuming an arbitrary potential and non-canonical kinetic terms, and we use them both to build phenomenological intuition and to constrain classes of two-field models using WMAP data. Using covariant formalism, we first develop a framework for understanding the background field kinematics and introduce a "slow-turn" approximation. Next, we find covariant expressions for the evolution of the adiabatic/curvature and entropy/isocurvature modes, and we discuss how the mode evolution can be inferred directly from the background kinematics and the geometry of the field manifold. From these expressions, we derive semi-analytic formulae for the curvature, isocurvature, and cross spectra, and the spectral observables, all to second-order in the slow-roll and slow-turn approximations. In tandem, we show how our covariant formalism provides useful intuition into how the characteristics of the inflationary Lagrangian translate into distinct features in the power spectra. In particular, we find that key features of the power spectra can be directly read off of the nature of the roll path, the curve the field vector rolls along with respect to the field manifold. For example, models whose roll path makes a sharp turn 60 e-folds before inflation ends tend to be ruled out because they produce strong departures from scale invariance. Finally, we apply our formalism to confront four classes of two-field models with WMAP data, including doubly quadratic and quartic potentials and non-standard kinetic terms, showing how whether a model is ruled out depends not only on certain features of the inflationary Lagrangian, but also on the initial conditions. Ultimately, models must possess the right balance of kinematical and dynamical behaviors, which we capture in a set of functions that can be reconstructed from spectral observables.

preprint2009arXiv

Baryon Acoustic Oscillations in the Sloan Digital Sky Survey Data Release 7 Galaxy Sample

The spectroscopic Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) galaxy sample represents the final set of galaxies observed using the original SDSS target selection criteria. We analyse the clustering of galaxies within this sample, including both the Luminous Red Galaxy (LRG) and Main samples, and also include the 2-degree Field Galaxy Redshift Survey (2dFGRS) data. Baryon Acoustic Oscillations are observed in power spectra measured for different slices in redshift; this allows us to constrain the distance--redshift relation at multiple epochs. We achieve a distance measure at redshift z=0.275, of r_s(z_d)/D_V(0.275)=0.1390+/-0.0037 (2.7% accuracy), where r_s(z_d) is the comoving sound horizon at the baryon drag epoch, D_V(z)=[(1+z)^2D_A^2cz/H(z)]^(1/3), D_A(z) is the angular diameter distance and H(z) is the Hubble parameter. We find an almost independent constraint on the ratio of distances D_V(0.35)/D_V(0.2)=1.736+/-0.065, which is consistent at the 1.1sigma level with the best fit Lambda-CDM model obtained when combining our z=0.275 distance constraint with the WMAP 5-year data. The offset is similar to that found in previous analyses of the SDSS DR5 sample, but the discrepancy is now of lower significance, a change caused by a revised error analysis and a change in the methodology adopted, as well as the addition of more data. Using WMAP5 constraints on Omega_bh^2 and Omega_ch^2, and combining our BAO distance measurements with those from the Union Supernova sample, places a tight constraint on Omega_m=0.286+/-0.018 and H_0 = 68.2+/-2.2km/s/Mpc that is robust to allowing curvature and non-Lambda dark energy. This result is independent of the behaviour of dark energy at redshifts greater than those probed by the BAO and supernova measurements. (abridged)

preprint2009arXiv

Cosmological Constraints from the Clustering of the Sloan Digital Sky Survey DR7 Luminous Red Galaxies

We present the power spectrum of the reconstructed halo density field derived from a sample of Luminous Red Galaxies (LRGs) from the Sloan Digital Sky Survey Seventh Data Release (DR7). The halo power spectrum has a direct connection to the underlying dark matter power for k <= 0.2 h/Mpc, well into the quasi-linear regime. This enables us to use a factor of ~8 more modes in the cosmological analysis than an analysis with kmax = 0.1 h/Mpc, as was adopted in the SDSS team analysis of the DR4 LRG sample (Tegmark et al. 2006). The observed halo power spectrum for 0.02 < k < 0.2 h/Mpc is well-fit by our model: chi^2 = 39.6 for 40 degrees of freedom for the best fit LCDM model. We find Ω_m h^2 * (n_s/0.96)^0.13 = 0.141^{+0.009}_{-0.012} for a power law primordial power spectrum with spectral index n_s and Ω_b h^2 = 0.02265 fixed, consistent with CMB measurements. The halo power spectrum also constrains the ratio of the comoving sound horizon at the baryon-drag epoch to an effective distance to z=0.35: r_s/D_V(0.35) = 0.1097^{+0.0039}_{-0.0042}. Combining the halo power spectrum measurement with the WMAP 5 year results, for the flat LCDM model we find Ω_m = 0.289 +/- 0.019 and H_0 = 69.4 +/- 1.6 km/s/Mpc. Allowing for massive neutrinos in LCDM, we find \sum m_ν < 0.62 eV at the 95% confidence level. If we instead consider the effective number of relativistic species Neff as a free parameter, we find Neff = 4.8^{+1.8}_{-1.7}. Combining also with the Kowalski et al. (2008) supernova sample, we find Ω_{tot} = 1.011 +/- 0.009 and w = -0.99 +/- 0.11 for an open cosmology with constant dark energy equation of state w.

preprint2009arXiv

The Second Law and Cosmology

I use cosmology examples to illustrate that the second law of thermodynamics is not old and tired, but alive and kicking, continuing to stimulate interesting research on really big puzzles. The question "Why is the entropy so low?" (despite the second law) suggests that our observable universe is merely a small and rather uniform patch in a vastly larger space stretched out by cosmological inflation. The question "Why is the entropy so high" (compared to the complexity required to describe many candidate "theories of everything") independently suggests that physical reality is much larger than the part we can observe.

preprint2005arXiv

Joint Efficient Dark-energy Investigation (JEDI): a Candidate Implementation of the NASA-DOE Joint Dark Energy Mission (JDEM)

JEDI will probe dark energy in three independent ways by measuring the expansion history of the universe: (1) using type Ia supernovae as cosmological standard candles over a range of distances, (2) using baryon oscillations as a cosmological standard ruler over a range of cosmic epochs, (3) mapping the weak gravitational lensing distortion by foreground galaxies of the images of background galaxies at different distances. JEDI will unravel the nature of dark energy with accuracy and precision. JEDI is a 2m-class space telescope with the unique ability of simultaneous wide-field imaging (0.8-4.2 micron in five bands) and multi-slit spectroscopy (0.8-3.2 micron) with a field of view of 1 square degree. What makes JEDI efficient is its ability to simultaneously obtain high signal-to-noise ratio, moderate resolution slit spectra for all supernovae and ~ 5000 galaxies in its wide field of view, and to combine imaging and spectroscopy so that the appropriate balance of time is devoted to each. JEDI will measure the cosmic expansion history H(z) as a free function to < 2% accuracy in redshift bins of 0.2-0.3. Assuming a flat universe and σ(Ω_m)=0.01 (0.03), JEDI could measure a departure from a vanilla ΛCDM model (w_0=-1, w'=0) with σ(w_0)=0.013 (0.031) and σ(w')=0.046 (0.063). JEDI will obtain the well-sampled lightcurves in Z, J, H, K, L bands and spectra of ~ 14,000 type Ia supernovae with redshifts ranging from 0 to 1.7; the redshifts of ~ 10-100 million galaxies to H ~ 23 and z ~ 4 over 1000-10,000 square degrees; and measurements of the shapes of galaxies over 1000-10,000 square degrees in Z,J,H,K,L for > 10^9 galaxies to H \~ 25.

preprint2003arXiv

CMBfit: Rapid WMAP likelihood calculations with normal parameters

We present a method for ultra-fast confrontation of the WMAP cosmic microwave background observations with theoretical models, implemented as a publicly available software package called CMBfit, useful for anyone wishing to measure cosmological parameters by combining WMAP with other observations. The method takes advantage of the underlying physics by transforming into a set of parameters where the WMAP likelihood surface is accurately fit by the exponential of a quartic or sextic polynomial. Building on previous physics based approximations by Hu et.al., Kosowsky et.al. and Chu et.al., it combines their speed with precision cosmology grade accuracy. A Fortran code for computing the WMAP likelihood for a given set of parameters is provided, pre-calibrated against CMBfast, accurate to Delta lnL ~ 0.05 over the entire 2sigma region of the parameter space for 6 parameter ``vanilla'' Lambda CDM models. We also provide 7-parameter fits including spatial curvature, gravitational waves and a running spectral index.

preprint2003arXiv

The Large-Scale Polarization of the Microwave Foreground

Most of the useful information about inflationary gravitational waves and reionization is on large angular scales where Galactic foreground contamination is the worst, so a key challenge is to model, quantify and remove polarized foregrounds. We use the Leiden radio surveys to quantify the polarized synchrotron radiation at large angular scales, which is likely to be the most challenging polarized contaminant for the WMAP satellite. We find that the synchrotron E- and B-contributions are equal to within 10% from 408-820MHz with a hint of E-domination at higher frequencies. We quantify Faraday Rotation & Depolarization effects and show that they cause the synchrotron polarization percentage to drop both towards lower frequencies and towards lower multipoles.

preprint1998arXiv

Cosmic Complementarity: H_0 and Omega_m from Combining CMB Experiments and Redshift Surveys

We show that upcoming CMB satellite experiments and large redshift surveys can be used together to yield 5% determinations of H_0 and Omega_m, an order of magnitude improvement over CMB data alone. CMB anisotropies provide the sound horizon at recombination as a standard ruler. For reasonable baryon fractions, this scale is imprinted on the galaxy power spectrum as a series of spectral features. Measuring these features in redshift space determines the Hubble constant, which in turn yields Omega_m once combined with CMB data. Since the oscillations in both power spectra are frozen in at recombination, this test is insensitive to low-redshift cosmology.

preprint1998arXiv

Is "the theory of everything'' merely the ultimate ensemble theory?

We discuss some physical consequences of what might be called ``the ultimate ensemble theory'', where not only worlds corresponding to say different sets of initial data or different physical constants are considered equally real, but also worlds ruled by altogether different equations. The only postulate in this theory is that all structures that exist mathematically exist also physically, by which we mean that in those complex enough to contain self-aware substructures (SASs), these SASs will subjectively perceive themselves as existing in a physically ``real'' world. We find that it is far from clear that this simple theory, which has no free parameters whatsoever, is observationally ruled out. The predictions of the theory take the form of probability distributions for the outcome of experiments, which makes it testable. In addition, it may be possible to rule it out by comparing its a priori predictions for the observable attributes of nature (the particle masses, the dimensionality of spacetime, etc) with what is observed.

preprint1997arXiv

Analyzing redshift surveys to measure the power spectrum on large scales

Upcoming large redshift surveys potentially allow precision measurements of the galaxy power spectrum. To accurately measure P(k) on the largest scales, comparable to the depth of the survey, it is crucial that finite volume effects are accurately corrected for in the data analysis. Here we derive analytic expressions for the one such effect that has not previously been worked out exactly: that of the so-called integral constraint. We also show that for data analysis methods based on counts in cells, multiple constraints can be included via simple matrix operations, thereby rendering the results less sensitive to galactic extinction and misestimates of the shape of the radial selection function.

preprint1997arXiv

Forecasting Cosmic Parameter Errors from Microwave Background Anisotropy Experiments

Accurate measurements of the cosmic microwave background (CMB) anisotropies with an angular resolution of a few arcminutes can be used to determine fundamental cosmological parameters such as the densities of baryons, cold and hot dark matter, and certain combinations of the cosmological constant and the curvature of the Universe to percent-level precision. Assuming the true theory is a variant of inflationary cold dark matter cosmologies, we calculate the accuracy with which cosmological parameters can be determined by the next generation of CMB satellites, MAP and Planck. We pay special attention to: (a) the accuracy of the computed derivatives of the CMB power spectrum C_L; (b) the number and choices of parameters; (c) the inclusion of prior knowledge of the values of various parameters.

preprint1997arXiv

Is lensing of point sources a problem for future CMB experiments?

Weak gravitational lensing from large-scale structure enhances and reduces the fluxes from extragalactic point sources with an rms amplitude of order 15%. In cosmic microwave background (CMB) experiments, sources exceeding some flux threshold are removed, which means that lensing will modulate the brightness map of the remaining unresolved sources. Since this mean brightness is of order 100 micro-Kelvin at 30 GHz for a reasonable flux cut, one might be concerned that this modulation could cause substantial problems for future CMB experiments. We present a detailed calculation of this effect and, fortunately, find that its power spectrum is always smaller than the normal point source power spectrum. Thus although this effect should be taken into account when analysing future high-precision CMB measurements, it will not substantially reduce the accuracy with which cosmological parameters can be measured.

preprint1997arXiv

On the dimensionality of spacetime

Some superstring theories have more than one effective low-energy limit, corresponding to classical spacetimes with different dimensionalities. We argue that all but the 3+1-dimensional one might correspond to ``dead worlds'', devoid of observers, in which case all such ensemble theories would actually predict that we should find ourselves inhabiting a 3+1-dimensional spacetime. With more or less than one time-dimension, the partial differential equations of nature would lack the hyperbolicity property that enables observers to make predictions. In a space with more than three dimensions, there can be no traditional atoms and perhaps no stable structures. A space with less than three dimensions allows no gravitational force and may be too simple and barren to contain observers.

preprint1997arXiv

The Interpretation of Quantum Mechanics: Many Worlds or Many Words?

As cutting-edge experiments display ever more extreme forms of non-classical behavior, the prevailing view on the interpretation of quantum mechanics appears to be gradually changing. A (highly unscientific) poll taken at the 1997 UMBC quantum mechanics workshop gave the once all-dominant Copenhagen interpretation less than half of the votes. The Many Worlds interpretation (MWI) scored second, comfortably ahead of the Consistent Histories and Bohm interpretations. It is argued that since all the above-mentioned approaches to nonrelativistic quantum mechanics give identical cookbook prescriptions for how to calculate things in practice, practical-minded experimentalists, who have traditionally adopted the ``shut-up-and-calculate interpretation'', typically show little interest in whether cozy classical concepts are in fact real in some untestable metaphysical sense or merely the way we subjectively perceive a mathematically simpler world where the Schrodinger equation describes everything - and that they are therefore becoming less bothered by a profusion of worlds than by a profusion of words. Common objections to the MWI are discussed. It is argued that when environment-induced decoherence is taken into account, the experimental predictions of the MWI are identical to those of the Copenhagen interpretation except for an experiment involving a Byzantine form of ``quantum suicide''. This makes the choice between them purely a matter of taste, roughly equivalent to whether one believes mathematical language or human language to be more fundamental.

preprint1996arXiv

A method for extracting maximum resolution power spectra from microwave sky maps

A method for extracting maximal resolution power spectra from microwave sky maps is presented and applied to the 2 year COBE data, yielding a power spectrum that is consistent with a standard n=1, Q=20 micro-Kelvin model. By using weight functions that fall off smoothly near the galactic cut, it is found that the spectral resolution Δl can be more than doubled at l=15 and more than tripled at l=20 compared to simply using galaxy-cut spherical harmonics. For a future high-resolution experiment with reasonable sky coverage, the resolution around the CDM Doppler peaks would be enhanced by a factor of about 100, down to Δl\approx 1, thus allowing spectral features such as the locations of the peaks to be determined with great accuracy. The reason that the improvement is so large is basically that functions with a sharp edge at the galaxy cut exhibit considerable "ringing" in the Fourier domain, whereas smooth functions do not. The method presented here is applicable to any survey geometry, chopping strategy and exposure pattern whatsoever. The so called signal-to-noise eigenfunction technique is found to be a special case, corresponding to ignoring the width of the window functions.

preprint1996arXiv

A method for subtracting foregrounds from multi-frequency CMB sky maps

An improved method for subtracting contaminants from Cosmic Microwave Background (CMB) sky maps is presented, and used to estimate how well future experiments will be able to recover the primordial CMB fluctuations. We find that the naive method of subtracting foregrounds (such as dust emission, synchrotron radiation, free-free-emission, unresolved point sources, etc) on a pixel by pixel basis can be improved by more than an order of magnitude by taking advantage of the correlation of the emission in neighboring pixels. The optimal multi-frequency subtraction method improves on simple pixel-by-pixel subtraction both by taking noise-levels into account, and by exploiting the fact that most contaminants have angular power spectra that differ substantially from that of the CMB. The results are natural to visualize in the two-dimensional plane with axes defined by multipole l and frequency v. We present a brief overview of the geography of this plane, showing the regions probed by various experiments and where we expect contaminants to dominate. We illustrate the method by estimating how well the proposed ESA COBRAS/SAMBA mission will be able to recover the CMB fluctuations against contaminating foregrounds.

preprint1996arXiv

An icosahedron-based method for pixelizing the celestial sphere

For power spectrum estimation it's important that the pixelization of a CMB sky map be smooth and regular to high degree. With this criterion in mind the ``COBE sky cube" was defined. This paper has as central theme to further improve on this elegant scheme which uses a cube as projective base -- here an icosahedron is used in its place. Although the sky cube is excellent, a further reduction of 20 percent of the number of pixels can be obtained while the pixel distance is maintained, and without any degradation of accuracy for integration. The pixels are rounder in this scheme, where they are hexagonal rather than square, and the faces are small in this implementation which simplifies area-equalization. The reason distortion is lessened is that the faces are smaller and therefore more flat. To use the method, you can get a FORTRAN code from the Internet.

preprint1995arXiv

ESTIMATING MICROWAVE POWER SPECTRA

A new method for estimating the power spectrum $C_l$ from cosmic microwave background (CMB) maps was recently presented by the author and applied to the 2 year COBE data. It was found that the spectral resolution $Δl$ for COBE could be more than doubled at $l=15$, thereby revealing previously unresolved features in the power spectrum. Whereas that paper was rather technical, the present paper shows that all qualitative features of this method can be understood from a simple analogy with quantum mechanics.

preprint1995arXiv

Using the kinematic Sunyaev-Zeldovich effect to determine the peculiar velocities of clusters of galaxies

We have investigated the possibility of inferring peculiar velocities for clusters of galaxies from the Doppler shift of scattered cosmic microwave background (CMB) photons. We find that if the core radius of the gas distribution or the beam size of the instrument is larger than 3-7 arcminutes, then the maximum attainable signal-to-noise ratio is determined by confusion with primary fluctuations. For smaller angular scales, ``cosmic confusion'' is less important and instrumental noise and/or foreground emission will be the limiting factor. For a cluster with the optical depth of the Coma cluster and for an optimal filtering technique, typical one-sigma errors span the wide range from 400 to 1600 km/s, depending on the cosmological model, the resolution of the instrument and the core radius of the cluster. The results have important implications for the design of future high-resolution surveys of the CMB. Individual peculiar velocities will be measurable only for a few fast moving clusters at intermediate redshift unless cosmic fluctuations are smaller than most standard cosmological scenarios predict. However, a reliable measurement of bulk velocities of ensembles of X-ray bright clusters will be possible on very large scales (100-500 Mpc/h).