BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Benjamin D. Wandelt

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

ResearcherUnclaimed source record
astro-ph.COastro-ph.IMgr-qcastro-phastro-ph.GAastro-ph.SRhep-phastro-ph.HEhep-thastro-ph.EPhep-exMachine Learning

Catalog footprint

What is connected

93works
12topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

93 published item(s)

preprint2022arXiv

Bayesian Control Variates for optimal covariance estimation with pairs of simulations and surrogates

Predictions of the mean and covariance matrix of summary statistics are critical for confronting cosmological theories with observations, not least for likelihood approximations and parameter inference. The price to pay for accurate estimates is the extreme cost of running $N$-body and hydrodynamics simulations. Approximate solvers, or surrogates, greatly reduce the computational cost but can introduce significant biases, for example in the non-linear regime of cosmic structure growth. We propose "CARPool Bayes", an approach to solve the inference problem for both the means and covariances using a combination of simulations and surrogates. Our framework allows incorporating prior information for the mean and covariance. We derive closed-form solutions for Maximum A Posteriori covariance estimates that are efficient Bayesian shrinkage estimators, guarantee positive semi-definiteness, and can optionally leverage analytical covariance approximations. We discuss choices of the prior and propose a simple procedure for obtaining optimal prior hyperparameter values with a small set of test simulations. We test our method by estimating the covariances of clustering statistics of GADGET-III $N$-body simulations at redshift $z=0.5$ using surrogates from a 100-1000$\times$ faster particle-mesh code. Taking the sample covariance from 15,000 simulations as the truth, and using an empirical Bayes prior with diagonal blocks, our estimator produces nearly identical Fisher matrix contours for $Λ$CDM parameters using only $15$ simulations of the non-linear dark matter power spectrum. In this case the number of simulations is so small that the sample covariance would be degenerate. We show cases where even with a naïve prior our method still improves the estimate. Our framework is applicable to a wide range of cosmological and astrophysical problems where fast surrogates are available.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Breaking baryon-cosmology degeneracy with the electron density power spectrum

Uncertain feedback processes in galaxies affect the distribution of matter, currently limiting the power of weak lensing surveys. If we can identify cosmological statistics that are robust against these uncertainties, or constrain these effects by other means, then we can enhance the power of current and upcoming observations from weak lensing surveys such as DES, Euclid, the Rubin Observatory, and the Roman Space Telescope. In this work, we investigate the potential of the electron density auto-power spectrum as a robust probe of cosmology and baryonic feedback. We use a suite of (magneto-)hydrodynamic simulations from the CAMELS project and perform an idealized analysis to forecast statistical uncertainties on a limited set of cosmological and physically-motivated astrophysical parameters. We find that the electron number density auto-correlation, measurable through either kinematic Sunyaev-Zel'dovich observations or through Fast Radio Burst dispersion measures, provides tight constraints on $Ω_{m}$ and the mean baryon fraction in intermediate-mass halos, $\bar{f}_{\mathrm{bar}}$. By obtaining an empirical measure for the associated systematic uncertainties, we find these constraints to be largely robust to differences in baryonic feedback models implemented in hydrodynamic simulations. We further discuss the main caveats associated with our analysis, and point out possible directions for future work.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

CARPool Covariance: Fast, unbiased covariance estimation for large-scale structure observables

The covariance matrix $\boldsymbolΣ$ of non-linear clustering statistics that are measured in current and upcoming surveys is of fundamental interest for comparing cosmological theory and data and a crucial ingredient for the likelihood approximations underlying widely used parameter inference and forecasting methods. The extreme number of simulations needed to estimate $\boldsymbolΣ$ to sufficient accuracy poses a severe challenge. Approximating $\boldsymbolΣ$ using inexpensive but biased surrogates introduces model error with respect to full simulations, especially in the non-linear regime of structure growth. To address this problem we develop a matrix generalization of Convergence Acceleration by Regression and Pooling (CARPool) to combine a small number of simulations with fast surrogates and obtain low-noise estimates of $\boldsymbolΣ$ that are unbiased by construction. Our numerical examples use CARPool to combine GADGET-III $N$-body simulations with fast surrogates computed using COmoving Lagrangian Acceleration (COLA). Even at the challenging redshift $z=0.5$, we find variance reductions of at least $\mathcal{O}(10^1)$ and up to $\mathcal{O}(10^4)$ for the elements of the matter power spectrum covariance matrix on scales $8.9\times 10^{-3}<k_\mathrm{max} <1.0$ $h {\rm Mpc^{-1}}$. We demonstrate comparable performance for the covariance of the matter bispectrum, the matter correlation function and probability density function of the matter density field. We compare eigenvalues, likelihoods, and Fisher matrices computed using the CARPool covariance estimate with the standard sample covariance estimators and generally find considerable improvement except in cases where $Σ$ is severely ill-conditioned.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Implicit Likelihood Inference of Reionization Parameters from the 21 cm Power Spectrum

The first measurements of the 21 cm brightness temperature power spectrum from the epoch of reionization will very likely be achieved in the near future by radio interferometric array experiments such as the Hydrogen Epoch of Reionization Array (HERA) and the Square Kilometre Array (SKA). Standard MCMC analyses use an explicit likelihood approximation to infer the reionization parameters from the 21 cm power spectrum. In this paper, we present a new Bayesian inference of the reionization parameters where the likelihood is implicitly defined through forward simulations using density estimation likelihood-free inference (DELFI). Realistic effects including thermal noise and foreground avoidance are also applied to the mock observations from the HERA and SKA. We demonstrate that this method recovers accurate posterior distributions for the reionization parameters, and outperforms the standard MCMC analysis in terms of the location and size of credible parameter regions. With the minutes-level processing time once the network is trained, this technique is a promising approach for the scientific interpretation of future 21 cm power spectrum observation data. Our code 21cmDELFI-PS is publicly available at this link.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Snowmass2021 CMB-HD White Paper

CMB-HD is a proposed millimeter-wave survey over half the sky that would be ultra-deep (0.5 uK-arcmin) and have unprecedented resolution (15 arcseconds at 150 GHz). Such a survey would answer many outstanding questions about the fundamental physics of the Universe. Major advances would be 1.) the use of gravitational lensing of the primordial microwave background to map the distribution of matter on small scales (k~10 h Mpc^(-1)), which probes dark matter particle properties. It will also allow 2.) measurements of the thermal and kinetic Sunyaev-Zel'dovich effects on small scales to map the gas density and velocity, another probe of cosmic structure. In addition, CMB-HD would allow us to cross critical thresholds: 3.) ruling out or detecting any new, light (< 0.1 eV) particles that were in thermal equilibrium with known particles in the early Universe, 4.) testing a wide class of multi-field models that could explain an epoch of inflation in the early Universe, and 5.) ruling out or detecting inflationary magnetic fields. CMB-HD would also provide world-leading constraints on 6.) axion-like particles, 7.) cosmic birefringence, 8.) the sum of the neutrino masses, and 9.) the dark energy equation of state. The CMB-HD survey would be delivered in 7.5 years of observing 20,000 square degrees of sky, using two new 30-meter-class off-axis crossed Dragone telescopes to be located at Cerro Toco in the Atacama Desert. Each telescope would field 800,000 detectors (200,000 pixels), for a total of 1.6 million detectors.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Cosmic Graph: Optimal Information Extraction from Large-Scale Structure using Catalogues

We present an implicit likelihood approach to quantifying cosmological information over discrete catalogue data, assembled as graphs. To do so, we explore cosmological parameter constraints using mock dark matter halo catalogues. We employ Information Maximising Neural Networks (IMNNs) to quantify Fisher information extraction as a function of graph representation. We a) demonstrate the high sensitivity of modular graph structure to the underlying cosmology in the noise-free limit, b) show that graph neural network summaries automatically combine mass and clustering information through comparisons to traditional statistics, c) demonstrate that networks can still extract information when catalogues are subject to noisy survey cuts, and d) illustrate how nonlinear IMNN summaries can be used as asymptotically optimal compressed statistics for Bayesian simulation-based inference. We reduce the area of joint $Ω_m, σ_8$ parameter constraints with small ($\sim$100 object) halo catalogues by a factor of 42 over the two-point correlation function, and demonstrate that the networks automatically combine mass and clustering information. This work utilises a new IMNN implementation over graph data in Jax, which can take advantage of either numerical or auto-differentiability. We also show that graph IMNNs successfully compress simulations away from the fiducial model at which the network is fitted, indicating a promising alternative to n-point statistics in catalogue simulation-based analyses.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Accurate precision Cosmology with redshift unknown gravitational wave sources

Gravitational waves can provide an accurate measurement of the luminosity distance to the source, but cannot provide the source redshift unless the degeneracy between mass and redshift can be broken. This makes it essential to infer the redshift of the source independently to measure the expansion history of the Universe. We show that by exploiting the clustering scale of the gravitational wave sources with galaxies of known redshift, we can infer the expansion history from redshift unknown gravitational wave sources. By using gravitational wave sources of unknown redshift that are detectable from the network of gravitational wave detectors with Advanced LIGO design sensitivity, we will be able to obtain accurate and precise measurements of the local Hubble constant, the expansion history of the universe, and the gravitational wave bias parameter, which captures the distribution of gravitational wave sources with respect to the redshift tracer distribution. While we showcase its application to low redshift gravitational waves, this technique will be applicable also to the high redshift gravitational wave sources detectable from Laser Interferometer Space Antenna (LISA), Cosmic Explorer (CE), and Einstein Telescope (ET). Moreover, this method will also be applicable to samples of supernovae and fast radio bursts with unknown or photometric redshifts.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Simulation-Based Inference of Reionization Parameters From 3D Tomographic 21 cm Lightcone Images

Tomographic three-dimensional 21 cm images from the epoch of reionization contain a wealth of information about the reionization of the intergalactic medium by astrophysical sources. Conventional power spectrum analysis cannot exploit the full information in the 21 cm data because the 21 cm signal is highly non-Gaussian due to reionization patchiness. We perform a Bayesian inference of the reionization parameters where the likelihood is implicitly defined through forward simulations using density estimation likelihood-free inference (DELFI). We adopt a trained 3D Convolutional Neural Network (CNN) to compress the 3D image data into informative summaries (DELFI-3D CNN). We show that this method recovers accurate posterior distributions for the reionization parameters. Our approach outperforms earlier analysis based on two-dimensional 21 cm images. In contrast, an MCMC analysis of the 3D lightcone-based 21 cm power spectrum alone and using a standard explicit likelihood approximation results in less accurate credible parameter regions than inferred by the DELFI-3D CNN, both in terms of the location and shape of the contours. Our proof-of-concept study implies that the DELFI-3D CNN can effectively exploit more information in the 3D 21 cm images than a 2D CNN or power spectrum analysis. This technique can be readily extended to include realistic effects and is therefore a promising approach for the scientific interpretation of future 21 cm observation data.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

SSSpaNG! Stellar Spectra as Sparse, data-driven, Non-Gaussian processes

Upcoming million-star spectroscopic surveys have the potential to revolutionize our view of the formation and chemical evolution of the Milky Way. Realizing this potential requires automated approaches to optimize estimates of stellar properties, such as chemical element abundances, from the spectra. The volume and quality of the observations strongly motivate that these approaches should be data-driven. With this in mind, we introduce SSSpaNG: a data-driven non-Gaussian Process model of stellar spectra. We demonstrate the capabilities of SSSpaNG using a sample of APOGEE red clump stars, whose model parameters we infer via Gibbs sampling. Pooling information between stars to infer their covariance, we permit clear identification of the correlations between spectral pixels. Harnessing these correlations, we infer the true spectrum of each star, inpainting missing regions and denoising by a factor of at least 2 for stars with signal-to-noise of ~20. As we marginalize over the covariance matrix of the spectra, the effective prior on these true spectra is non-Gaussian and sparsifying, favouring typically small but occasionally large excursions from the mean. The high-fidelity inferred spectra produced will enable improved elemental abundance measurements for individual stars. Our model also allows us to quantify the information gained by observing portions of a star's spectrum, and thereby define the most mutually informative spectral regions. Using 25 windows centred on elemental absorption lines, we demonstrate that the iron-peak and alpha-process elements are particularly mutually informative for these spectra, and that the majority of information about a target window is contained in the 10-or-so most informative windows. Such mutual-information estimates have the potential to inform models of nucleosynthetic yields and the design of future observations.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Testing the general theory of relativity using gravitational wave propagation from dark standard sirens

Alternative theories of gravity predict modifications in the propagation of gravitational waves (GW) through space-time. One of the smoking-gun predictions of such theories is the change in the GW luminosity distance to GW sources as a function of redshift relative to the electromagnetic (EM) luminosity distance expected from EM probes. We propose a multi-messenger test of the theory of general relativity from the propagation of gravitational waves by combining EM and GW observations to resolve these issues from GW sources without EM counterparts (which are also referred to as dark standard sirens). By using the relation between the geometric distances accessible from baryon acoustic oscillation measurements, and luminosity distance measurements from the GW sources, we can measure any deviation from the general theory of relativity via the GW sources of unknown redshift that will be detectable by networks of GW detectors such as LIGO, Virgo, and KAGRA. Using this technique, the fiducial value of the frictional term can be measured to a precision $Ξ_0=0.98^{+0.04}_{-0.23}$ after marginalizing over redshift dependence, cosmological parameters, and GW bias parameters with $\sim 3500$ dark standard sirens of masses $30\,\rm M_\odot$ each distributed up to redshift $z=0.5$. For a fixed redshift dependence, a value of $Ξ_0=0.99^{+0.02}_{-0.02}$ can be measured with a similar number of dark sirens. Application of our methodology to the far more numerous dark standard sirens detectable with next generation GW detectors, such as LISA, Einstein Telescope and Cosmic Explorer, will allow achievement of higher accuracy than possible from use of bright standard sirens.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The Age of the Universe with Globular Clusters: Reducing Systematic Uncertainties

The dominant systematic uncertainty in the age determination of galactic globular clusters is the depth of the convection envelope of the stars. This parameter is partially degenerate with metallicity which is in turn degenerate with age. However, if the metal content, distance and extinction are known, the position and morphology of the red giant branch in a color-magnitude diagram are mostly sensitive to the value of the depth of the convective envelope. Therefore, using external, precise metallicity determinations this degeneracy and thus the systematic error in age, can be reduced. Alternatively, the morphology of the red giant branch of globular clusters color magnitude diagram can also be used to achieve the same. We demonstrate that globular cluster red giant branches are well fitted by values of the depth of the convection envelope consistent with those obtained for the Sun and this finding is robust to the adopted treatment of the stellar physics. With these findings, the uncertainty in the depth of the convection envelope is no longer the dominant contribution to the systematic error in the age determination of the oldest globular clusters, reducing it from $0.5$ to $0.23$ or $0.33$ Gyr, depending on the methodology adopted: i.e., whether resorting to external data (spectroscopic metallicity determinations) or relying solely on the morphology of the clusters's color-magnitude diagrams. This results in an age of the Universe $t_{\rm U}=13.5^{+0.16}_{-0.14} {\rm (stat.)} \pm 0.23(0.33) ({\rm sys.})$ at 68\% confidence level, accounting for the formation time of globular clusters and its uncertainty. An uncertainty of 0.27(0.36) Gyr if added in quadrature. This agrees well with $13.8 \pm 0.02$ Gyr, the cosmological model-dependent value inferred by the Planck mission assuming the $Λ$CDM model.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The GIGANTES dataset: precision cosmology from voids in the machine learning era

We present GIGANTES, the most extensive and realistic void catalog suite ever released -- containing over 1 billion cosmic voids covering a volume larger than the observable Universe, more than 20 TB of data, and created by running the void finder VIDE on QUIJOTE's halo simulations. The expansive and detailed GIGANTES suite, spanning thousands of cosmological models, opens up the study of voids, answering compelling questions: Do voids carry unique cosmological information? How is this information correlated with galaxy information? Leveraging the large number of voids in the GIGANTES suite, our Fisher constraints demonstrate voids contain additional information, critically tightening constraints on cosmological parameters. We use traditional void summary statistics (void size function, void density profile) and the void auto-correlation function, which independently yields an error of $0.13\,\mathrm{eV}$ on $\sum\,m_ν$ for a 1 $h^{-3}\mathrm{Gpc}^3$ simulation, without CMB priors. Combining halos and voids we forecast an error of $0.09\,\mathrm{eV}$ from the same volume. Extrapolating to next generation multi-Gpc$^3$ surveys such as DESI, Euclid, SPHEREx, and the Roman Space Telescope, we expect voids should yield an independent determination of neutrino mass. Crucially, GIGANTES is the first void catalog suite expressly built for intensive machine learning exploration. We illustrate this by training a neural network to perform likelihood-free inference on the void size function. Cosmology problems provide an impetus to develop novel deep learning techniques, leveraging the symmetries embedded throughout the universe from physical laws, interpreting models, and accurately predicting errors. With GIGANTES, machine learning gains an impressive dataset, offering unique problems that will stimulate new techniques.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A new probe of Axion-Like Particles: CMB polarization distortions due to cluster magnetic fields

We propose using the upcoming Cosmic Microwave Background (CMB) ground based experiments to detect the signal of ALPs (Axion like particles) interacting with magnetic fields in galaxy clusters. The conversion between CMB photons and ALPs in the presence of the cluster magnetic field can cause a polarized spectral distortion in the CMB around a galaxy cluster. The strength of the signal depends upon the redshift of the galaxy cluster and will exhibit a distinctive spatial profile around it depending upon the structure of electron density and magnetic field. This distortion produces a different shape from the other known spectral distortions like $y$-type and $μ$-type and hence are separable from the multi-frequency CMB observation. The spectrum is close to kinematic Sunyaev-Zeldovich (kSZ) signal but can be separated from it using the polarization information. For the future ground-based CMB experiments such as Simons Observatory and CMB-S4, we estimate the measurability of this signal in the presence of foreground contamination, instrument noise and CMB anisotropies. This new avenue can probe the photon-ALP coupling over the ALP mass range from $10^{-13}$ eV to $10^{-12}$ eV with two orders of magnitude better accuracy from CMB-S4 than the current existing bounds.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Bayesian delensing delight: sampling-based inference of the primordial CMB and gravitational lensing

The search for primordial gravitational waves in the Cosmic Microwave Background (CMB) will soon be limited by our ability to remove the lensing contamination to $B$-mode polarization. The often-used quadratic estimator for lensing is known to be suboptimal for surveys that are currently operating and will continue to become less and less efficient as instrumental noise decreases. While foregrounds can in principle be mitigated by observing in more frequency bands, progress in delensing hinges entirely on algorithmic advances. We demonstrate here a new inference method that solves this problem by sampling the exact Bayesian posterior of any desired cosmological parameters, of the gravitational lensing potential, and of the delensed CMB maps, given lensed temperature and polarization data. We validate the method using simulated CMB data with non-white noise and masking on up to 650\,deg$^2$ patches of sky. A unique strength of this approach is the ability to jointly estimate cosmological parameters which control both the primordial CMB and the lensing potential, which we demonstrate here for the first time by sampling both the tensor-to-scalar ratio, $r$, and the amplitude of the lensing potential, $A_ϕ$. The method allows us to perform the most precise check to-date of several important approximations underlying CMB-S4 $r$ forecasting, and we confirm these yield the correct expected uncertainty on $r$ to better than 10%.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

How to measure CMB spectral distortions with an imaging telescope

We propose the implementation of an imaging telescope in combination with an inter-frequency calibrator to measure the spectral shape of the microwave sky by exploiting the differences in the sky intensity between multiple pairs of frequency channels. By jointly sampling the cosmological and foreground parameters in a Bayesian framework for $600$ instrument configurations, we determine the minimum calibration accuracy required in order to obtain measurements of spectral distortions and simultaneously measure spectral and spatial fluctuations of the CMB. We demonstrate the feasibility of this technique for a CMB mission like PICO (Probe of Inflation and Cosmic Origins), and show that a $10$-$σ$ measurement of the $y$-distortion along with a two orders of magnitude improvement on the FIRAS (Far-Infrared Absolute Spectrophotometer) $μ$-distortion limit is feasible from this technique. We argue that longer term applications may be envisaged at even higher sensitivity, capable of attaining the $μ$-distortion that provide a robust prediction of the $Λ$CDM model and even primordial recombination lines of hydrogen and helium, in the context of the ESA (European Space Agency) Voyage 2035-2050 program.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Inferring the Age of the Universe with Globular Clusters

We present an estimate of the absolute age of 68 galactic globular clusters obtained by exploiting the distribution of stars in the full color-magnitude diagram. In particular, we jointly estimate the absolute age, distance, reddening, metallicity ([Fe/H]) and [$α$/Fe] of each cluster, imposing priors motivated by independent observations; we also estimate possible systematics from stellar modeling. Our derived distances for the globular cluster sample are in agreement with those obtained from GAIA using main-sequence dwarf stars (where available), and the inferred ages are in good agreement with those previously published. The novelty of our approach is that, with the adopted priors, we are able to estimate robustly these parameters from the globular cluster color-magnitude diagram. We find that the average age of the oldest globular clusters is $t_{\rm GC}=13.32 \pm 0.1 {\rm (stat.)} \pm 0.5 {\rm (sys.)}$, at 68\% confidence level, including systematic uncertainties from stellar modeling. These measurements can be used to infer the age of the Universe, largely independently of the cosmological parameters: we find an age of the Universe $t_{\rm U}=13.5^{+0.16}_{-0.14} {\rm (stat.)} \pm 0.5 ({\rm sys.})$ at 68\% confidence level, accounting for the formation time of globular clusters and its uncertainty. This value is compatible with $13.8 \pm 0.02$ Gyr, the cosmological model-dependent value inferred by the Planck mission assuming the $Λ$CDM model.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Multi-messenger tests of gravity with weakly lensed gravitational waves

General relativity (GR) predicts concordant trajectories for photons and gravitational waves (GW). We propose a new multi-messenger avenue (GW-CMB-CMB) to prove this aspect of fundamental physics by cross-correlating the GW signal of astrophysical origin with the lensing field derived from the cosmic microwave background (CMB). This new window will allow robust measurement of the prediction from GR with high signal-to-noise and will be able to unveil the true nature of gravity using the GW sources detected by missions such as the Laser Interferometer Space Antenna (LISA), Einstein Telescope and Cosmic Explorer.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Probing the theory of gravity with gravitational lensing of gravitational waves and galaxy surveys

The cross-correlation of gravitational wave strain with upcoming galaxy surveys probe theories of gravity in a new way. This method enables testing the theory of gravity by combining the effects from both gravitational lensing of gravitational waves and the propagation of gravitational waves in spacetime. We find that within 10 years, the combination of the Advanced-LIGO and VIRGO detector networks with planned galaxy surveys should detect weak gravitational lensing of gravitational waves in the low redshift Universe ($z<0.5$). With the next generation gravitational wave experiments such as Voyager, LISA, Cosmic-Explorer and Einstein Telescope, we can extend this test of the theory of gravity to larger redshifts by exploiting the synergies between electromagnetic wave and gravitational wave probes.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Super-resolution emulator of cosmological simulations using deep physical models

We present an extension of our recently developed Wasserstein optimized model to emulate accurate high-resolution features from computationally cheaper low-resolution cosmological simulations. Our deep physical modelling technique relies on restricted neural networks to perform a mapping of the distribution of the low-resolution cosmic density field to the space of the high-resolution small-scale structures. We constrain our network using a single triplet of high-resolution initial conditions and the corresponding low- and high-resolution evolved dark matter simulations from the Quijote suite of simulations. We exploit the information content of the high-resolution initial conditions as a well constructed prior distribution from which the network emulates the small-scale structures. Once fitted, our physical model yields emulated high-resolution simulations at low computational cost, while also providing some insights about how the large-scale modes affect the small-scale structure in real space.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Velocity correction for Hubble constant measurements from standard sirens

Gravitational wave (GW) sources are an excellent probe of the luminosity distance and offer a novel measure of the Hubble constant, $H_0$. This estimation of $H_0$ from standard sirens requires an accurate estimation of the cosmological redshift of the host galaxy of the GW source, after correcting for its peculiar velocity. Absence of an accurate peculiar velocity correction affects both the precision and accuracy of the measurement of $H_0$, particularly for nearby sources. We propose a framework to incorporate such a peculiar velocity correction for GW sources. A first implementation of our method to the event GW170817 combined with the Very Large Baseline Interferometry (VLBI) observation leads to a revised value of $H_0= 68.3^{+ 4.6}_{-4.5}$ km/s/Mpc. While this revision is minor, it demonstrates that our method makes it possible for obtaining an unbiased and accurate measurements of $H_0$ at the precision required for the standard siren cosmology.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Neural physical engines for inferring the halo mass distribution function

An ambitious goal in cosmology is to forward-model the observed distribution of galaxies in the nearby Universe today from the initial conditions of large-scale structures. For practical reasons, the spatial resolution at which this can be done is necessarily limited. Consequently, one needs a mapping between the density of dark matter averaged over ~Mpc scales, and the distribution of dark matter halos (used as a proxy for galaxies) in the same region. Here we demonstrate a method for determining the halo mass distribution function by learning the tracer bias between density fields and halo catalogues using a neural bias model. The method is based on the Bayesian analysis of simple, physically motivated, neural network-like architectures, which we denote as neural physical engines, and neural density estimation. As a result, we are able to sample the initial phases of the dark matter density field whilst inferring the parameters describing the halo mass distribution function, providing a fully Bayesian interpretation of both the initial dark matter density distribution and the neural bias model. We successfully run an upgraded BORG inference using our new likelihood and neural bias model with halo catalogues derived from full N-body simulations. We notice orders of magnitude improvement in modelling compared to classical biasing techniques.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Bayesian semi-blind component separation for foreground removal in interferometric 21-cm observations

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Constraints on Cosmology and Gravity from the Dynamics of Voids

The Universe is mostly composed of large and relatively empty domains known as cosmic voids, whereas its matter content is predominantly distributed along their boundaries. The remaining material inside them, either dark or luminous matter, is attracted to these boundaries and causes voids to expand faster and to grow emptier over time. Using the distribution of galaxies centered on voids identified in the Sloan Digital Sky Survey and adopting minimal assumptions on the statistical motion of these galaxies, we constrain the average matter content $Ω_\mathrm{m}=0.281\pm0.031$ in the Universe today, as well as the linear growth rate of structure $f/b=0.417\pm0.089$ at median redshift $\bar{z}=0.57$, where $b$ is the galaxy bias ($68\%$ C.L.). These values originate from a percent-level measurement of the anisotropic distortion in the void-galaxy cross-correlation function, $\varepsilon = 1.003\pm0.012$, and are robust to consistency tests with bootstraps of the data and simulated mock catalogs within an additional systematic uncertainty of half that size. They surpass (and are complementary to) existing constraints by unlocking cosmological information on smaller scales through an accurate model of nonlinear clustering and dynamics in void environments. As such, our analysis furnishes a powerful probe of deviations from Einstein's general relativity in the low-density regime which has largely remained untested so far. We find no evidence for such deviations in the data at hand.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Cosmological parameter constraints from CMB lensing with cosmic voids

We investigate the potential of using cosmic voids as a probe to constrain cosmological parameters through the gravitational lensing effect of the cosmic microwave background (CMB) and make predictions for the next generation surveys. By assuming the detection of a series of $\approx 5 - 10$ voids along a line of sight within a square-degree patch of the sky, we found that they can be used to break the degeneracy direction of some of the cosmological parameter constraints (for example $ω_b$ and $Ω_Λ$) in comparison with the constraints from random CMB skies with the same size area for a survey with extensive integration time. This analysis is based on our current knowledge of the average void profile and analytical estimates of the void number function. We also provide combined cosmological parameter constraints between a sky patch where series of voids are detected and a patch without voids (a randomly selected patch). The full potential of this technique relies on an accurate determination of the void profile to $\approx 10$% level. For a small-area CMB observation with extensive integration time and a high signal-to-noise ratio, CMB lensing with such series of voids will provide a complementary route to cosmological parameter constraints to the CMB observations. Example of parameter constraints with a series of five voids on a $1.0^{\circ} \times 1.0^{\circ}$ patch of the sky are $100ω_b = 2.20 \pm 0.27$, $ω_c = 0.120 \pm 0.022$, $Ω_Λ= 0.682 \pm 0.078$, $Δ_{\mathcal{R}}^2 = \left(2.22 \pm 7.79\right) \times 10^{-9}$, $n_s = 0.962 \pm 0.097$ and $τ= 0.925 \pm 1.747$ at 68% C.L.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Bayesian inference on the sphere beyond statistical isotropy

We present a general method for Bayesian inference of the underlying covariance structure of random fields on a sphere. We employ the Bipolar Spherical Harmonic (BipoSH) representation of general covariance structure on the sphere. We illustrate the efficacy of the method as a principled approach to assess violation of statistical isotropy (SI) in the sky maps of Cosmic Microwave Background (CMB) fluctuations. SI violation in observed CMB maps arise due to known physical effects such as Doppler boost and weak lensing; yet unknown theoretical possibilities like cosmic topology and subtle violations of the cosmological principle, as well as, expected observational artefacts of scanning the sky with a non-circular beam, masking, foreground residuals, anisotropic noise, etc. We explicitly demonstrate the recovery of the input SI violation signals with their full statistics in simulated CMB maps. Our formalism easily adapts to exploring parametric physical models with non-SI covariance, as we illustrate for the inference of the parameters of a Doppler boosted sky map. Our approach promises to provide a robust quantitative evaluation of the evidence for SI violation related anomalies in the CMB sky by estimating the BipoSH spectra along with their complete posterior.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Counting voids to probe dark energy

We show that the number of observed voids in galaxy redshift surveys is a sensitive function of the equation of state of dark energy. Using the Fisher matrix formalism we find the error ellipses in the $w_0-w_a$ plane when the equation of state of dark energy is assumed to be of the form $w_{CPL}(z)=w_0 +w_a z/(1+z)$. We forecast the number of voids to be observed with the ESA Euclid satellite and the NASA WFIRST mission, taking into account updated details of the surveys to reach accurate estimates of their power. The theoretical model for the forecast of the number of voids is based on matches between abundances in simulations and the analytical prediction. To take into account the uncertainties within the model, we marginalize over its free parameters when calculating the Fisher matrices. The addition of the void abundance constraints to the data from Planck, HST and supernova survey data noticeably tighten the $w_0-w_a$ parameter space. We thus quantify the improvement in the constraints due to the use of voids and demonstrate that the void abundance is a sensitive new probe for the dark energy equation of state.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Probing cosmology and gravity with redshift-space distortions around voids

Cosmic voids in the large-scale structure of the Universe affect the peculiar motions of objects in their vicinity. Although these motions are difficult to observe directly, the clustering pattern of their surrounding tracers in redshift space is influenced in a unique way. This allows to investigate the interplay between densities and velocities around voids, which is solely dictated by the laws of gravity. With the help of $N$-body simulations and derived mock-galaxy catalogs we calculate the average density fluctuations around voids identified with a watershed algorithm in redshift space and compare the results with the expectation from general relativity and the $Λ$CDM model. We find linear theory to work remarkably well in describing the dynamics of voids. Adopting a Bayesian inference framework, we explore the full posterior of our model parameters and forecast the achievable accuracy on measurements of the growth rate of structure and the geometric distortion through the Alcock-Paczynski effect. Systematic errors in the latter are reduced from $\sim15\%$ to $\sim5\%$ when peculiar velocities are taken into account. The relative parameter uncertainties in galaxy surveys with number densities comparable to the SDSS MAIN (CMASS) sample probing a volume of $1h^{-3}{\rm Gpc}^3$ yield $σ_{f/b}\left/(f/b)\right.\sim2\%$ ($20\%$) and $σ_{D_AH}/D_AH\sim0.2\%$ ($2\%$), respectively. At this level of precision the linear-theory model becomes systematics dominated, with parameter biases that fall beyond these values. Nevertheless, the presented method is highly model independent; its viability lies in the underlying assumption of statistical isotropy of the Universe.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Real-space density profile reconstruction of stacked voids

We present a non-parametric, model-independent method to reconstruct the spherical density profiles of void stacks in real space, without redshift-space distortions. Our method uses the expected spherical symmetry of stacked voids to build the shape of the spherical density profile of cosmic voids in real space without any assumption about the cosmological model. We test the reconstruction algorithm with both a toy model and a full dark matter simulation. We present the result for the simulation: the reconstruction of the spherical density profile for a simulated stacked void in real space. We also present a first application of the algorithm to reconstruct real cosmic void stacks density profiles in real space from the Sloan Digital Sky Survey (Sutter et al. 2012b). We discuss capabilities of the algorithm and possible future improvements. Reconstructed density profiles from real voids open the way to the study of the spherically averaged dynamical structure of voids.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

sCOLA: The N-body COLA Method Extended to the Spatial Domain

We present sCOLA -- an extension of the N-body COmoving Lagrangian Acceleration (COLA) method to the spatial domain. Similar to the original temporal-domain COLA, sCOLA is an N-body method for solving for large-scale structure in a frame that is comoving with observers following trajectories calculated in Lagrangian Perturbation Theory. Incorporating the sCOLA method in an N-body code allows one to gain computational speed by capturing the gravitational potential from the far field using perturbative techniques, while letting the N-body code solve only for the near field. The far and near fields are completely decoupled, effectively localizing gravity for the N-body side of the code. Thus, running an N-body code for a small simulation volume using sCOLA can reproduce the results of a standard N-body run for the same small volume embedded inside a much larger simulation. We demonstrate that sCOLA can be safely combined with the original temporal-domain COLA. sCOLA can be used as a method for performing zoom-in simulations. It also allows N-body codes to be made embarrassingly parallel, thus allowing for efficiently tiling a volume of interest using grid computing. Moreover, sCOLA can be useful for cheaply generating large ensembles of accurate mock halo catalogs required to study galaxy clustering. Surveys that will benefit the most are ones with large aspect ratios, such as pencil-beam surveys, where sCOLA can easily capture the effects of large-scale transverse modes without the need to substantially increase the simulated volume. As an illustration of the method, we present proof-of-concept zoom-in simulations using a freely available sCOLA-based N-body code.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Linear Perturbation Theory of Reionization in Position-Space: Cosmological Radiative Transfer Along the Light-Cone

The linear perturbation theory of inhomogeneous reionization (LPTR) has been developed as an analytical tool for predicting the global ionized fraction and large-scale power spectrum of ionized density fluctuations during reionization. In the original formulation of the LPTR, the ionization balance and radiative transfer equations are linearized and solved in Fourier space. However, the LPTR's approximation to the full solution of the radiative transfer equation is not straightforward to interpret, since the latter is most intuitively conceptualized in position space. To bridge the gap between the LPTR and the language of numerical radiative transfer, we present a new, equivalent, position-space formulation of the LPTR that clarifies the approximations it makes and facilitates its interpretation. We offer a comparison between the LPTR and the excursion-set model of reionization (ESMR), and demonstrate the built-in capability of the LPTR to explore a wide range of reionization scenarios, and to go beyond the ESMR in exploring scenarios involving X-rays.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Using cosmic voids to distinguish f(R) gravity in future galaxy surveys

We use properties of void populations identified in $N$-body simulations to forecast the ability of upcoming galaxy surveys to differentiate models of f(R) gravity from \lcdm~cosmology. We analyze multiple simulation realizations, which were designed to mimic the expected number densities, volumes, and redshifts of the upcoming Euclid satellite and a lower-redshift ground-based counterpart survey, using the public {\tt VIDE} toolkit. We examine void abundances, ellipicities, radial density profiles, and radial velocity profiles at redshifts 1.0 and 0.43. We find that stronger f(R) coupling strengths eliminates small voids and produces voids up to $\sim 20\%$ larger in radius, leading to a significant tilt in the void number function. Additionally, under the influence of modified gravity, voids at all scales tend to be measurably emptier with correspondingly higher compensation walls. The velocity profiles reflect this, showing increased outflows inside voids and increased inflows outside voids. Using the void number function as an example, we forecast that future surveys can constrain the modified gravity coupling strength to $\sim 3 \times 10^{-5}$ using voids.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Using hybrid GPU/CPU kernel splitting to accelerate spherical convolutions

We present a general method for accelerating by more than an order of magnitude the convolution of pixelated functions on the sphere with a radially-symmetric kernel. Our method splits the kernel into a compact real-space component and a compact spherical harmonic space component. These components can then be convolved in parallel using an inexpensive commodity GPU and a CPU. We provide models for the computational cost of both real-space and Fourier space convolutions and an estimate for the approximation error. Using these models we can determine the optimum split that minimizes the wall clock time for the convolution while satisfying the desired error bounds. We apply this technique to the problem of simulating a cosmic microwave background (CMB) anisotropy sky map at the resolution typical of the high resolution maps produced by the Planck mission. For the main Planck CMB science channels we achieve a speedup of over a factor of ten, assuming an acceptable fractional rms error of order 1.e-5 in the power spectrum of the output map.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A measurement of the Alcock-Paczynski effect using cosmic voids in the SDSS

We perform an Alcock-Paczynski test using stacked cosmic voids identified in the SDSS Data Release 7 main sample and Data Release 10 LOWZ and CMASS samples. We find ~1,500 voids out to redshift $0.6$ using a heavily modified and extended version of the watershed algorithm ZOBOV, which we call VIDE (Void IDentification and Examination). To assess the impact of peculiar velocities we use the mock void catalogs presented in Sutter et al. (2013). We find a constant uniform flattening of 14% along the line of sight when peculiar velocities are included. This flattening appears universal for all void sizes at all redshifts and for all tracer densities. We also use these mocks to identify an optimal stacking strategy. After correcting for systematic effects we find that our Alcock-Paczynski measurement leads to a preference of our best-fit value of $Ω_{\rm M}\sim 0.15$ over $Ω_{\rm M} = 1.0$ by a likelihood ratio of 10. Likewise, we find a factor of $4.5$ preference of the likelihood ratio for a $Λ$CDM $Ω_{\rm M} = 0.3$ model and a null measurement. Taken together, we find substantial evidence for the Alcock-Paczynski signal in our sample of cosmic voids. Our assessment using realistic mocks suggests that measurements with future SDSS releases and other surveys will provide tighter cosmological parameter constraints. The void-finding algorithm and catalogs used in this work will be made publicly available at http://www.cosmicvoids.net.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

An optimal estimator for resonance bispectra in the CMB

We propose an (optimal) estimator for a CMB bispectrum containing logarithmically spaced oscillations. There is tremendous theoretical interest in such bispectra, and they are predicted by a plethora of models, including axion monodromy models of inflation and initial state modifications. The number of resolved logarithmical oscillations in the bispectrum is limited due to the discrete resolution of the multipole bispectrum. We derive a simple relation between the maximum number of resolved oscillations and the frequency. We investigate several ways to factorize the primordial bispectrum, and conclude that a one dimensional expansion in the sum of the momenta $\sum k_i = k_t$ is the most efficient and flexible approach. We compare the expansion to the exact result in multipole space and show for $ω_{\rm eff}=100$ that $\mathcal{O}(10^3)$ modes are sufficient for an accurate reconstruction. We compute the expected $σ_{f_{\rm NL}}$ and find that within an effective field theory (EFT) the overall signal to noise scales as $S/N\propto ω^{3/2}$. Using only the temperature data we find $S/N\sim\mathcal{O}(1-10^2)$ for the frequency domain set by the EFT.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Cosmology with Void-Galaxy Correlations

Galaxy bias, the unknown relationship between the clustering of galaxies and the underlying dark matter density field is a major hurdle for cosmological inference from large-scale structure. While traditional analyses focus on the absolute clustering amplitude of high-density regions mapped out by galaxy surveys, we propose a relative measurement that compares those to the underdense regions, cosmic voids. On the basis of realistic mock catalogs we demonstrate that cross correlating galaxies and voids opens up the possibility to calibrate galaxy bias and to define a static ruler thanks to the observable geometric nature of voids. We illustrate how the clustering of voids is related to mass compensation and show that volume-exclusion significantly reduces the degree of stochasticity in their spatial distribution. Extracting the spherically averaged distribution of galaxies inside voids from their cross correlations reveals a remarkable concordance with the mass-density profile of voids.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Cosmology: from theory to data, from data to theory

Cosmology has come a long way from being based on a small number of observations to being a data-driven precision science. We discuss the questions "What is observable?", "What in the Universe is knowable?" and "What are the fundamental limits to cosmological knowledge?". We then describe the methodology for investigation: theoretical hypotheses are used to model, predict and anticipate results; data is used to infer theory. We illustrate with concrete examples of principled analysis approaches from the study of cosmic microwave background anisotropies and surveys of large-scale structure, culminating in a summary of the highest precision probe to date of the physical origin of cosmic structures: the Planck 2013 constraints on primordial non-Gaussianity.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

First measurement of gravitational lensing by cosmic voids in SDSS

We report the first measurement of the diminutive lensing signal arising from matter underdensities associated with cosmic voids. While undetectable individually, by stacking the weak gravitational shear estimates around 901 voids detected in SDSS DR7 by Sutter et al. (2012a), we find substantial evidence for a depression of the lensing signal compared to the cosmic mean. This depression is most pronounced at the void radius, in agreement with analytical models of void matter profiles. Even with the largest void sample and imaging survey available today, we cannot put useful constraints on the radial dark-matter void profile. We invite independent investigations of our findings by releasing data and analysis code to the public at https://github.com/pmelchior/void-lensing

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Modeling cosmic void statistics

Understanding the internal structure and spatial distribution of cosmic voids is crucial when considering them as probes of cosmology. We present recent advances in modeling void density- and velocity-profiles in real space, as well as void two-point statistics in redshift space, by examining voids identified via the watershed transform in state-of-the-art $Λ$CDM n-body simulations and mock galaxy catalogs. The simple and universal characteristics that emerge from these statistics indicate the self-similarity of large-scale structure and suggest cosmic voids to be among the most pristine objects to consider for future studies on the nature of dark energy, dark matter and modified gravity.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the observability of coupled dark energy with cosmic voids

Taking N-body simulations with volumes and particle densities tuned to match the SDSS DR7 spectroscopic main sample, we assess the ability of current void catalogs (e.g., Sutter et al. 2012b) to distinguish a model of coupled dark matter-dark energy from ΛCDM cosmology using properties of cosmic voids. Identifying voids with the VIDE toolkit, we find no statistically significant differences in the ellipticities, but find that coupling produces a population of significantly larger voids, possibly explaining the recent result of Tavasoli et al. (2013). In addition, we use the universal density profile of Hamaus et al. (2014) to quantify the relationship between coupling and density profile shape, finding that the coupling produces broader, shallower, undercompensated profiles for large voids by thinning the walls between adjacent medium-scale voids. We find that these differences are potentially measurable with existing void catalogs once effects from survey geometries and peculiar velocities are taken into account.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Past and present cosmic structure in the SDSS DR7 main sample

We present a chrono-cosmography project, aiming at the inference of the four dimensional formation history of the observed large scale structure from its origin to the present epoch. To do so, we perform a full-scale Bayesian analysis of the northern galactic cap of the Sloan Digital Sky Survey (SDSS) Data Release 7 main galaxy sample, relying on a fully probabilistic, physical model of the non-linearly evolved density field. Besides inferring initial conditions from observations, our methodology naturally and accurately reconstructs non-linear features at the present epoch, such as walls and filaments, corresponding to high-order correlation functions generated by late-time structure formation. Our inference framework self-consistently accounts for typical observational systematic and statistical uncertainties such as noise, survey geometry and selection effects. We further account for luminosity dependent galaxy biases and automatic noise calibration within a fully Bayesian approach. As a result, this analysis provides highly-detailed and accurate reconstructions of the present density field on scales larger than $\sim~3$ Mpc$/h$, constrained by SDSS observations. This approach also leads to the first quantitative inference of plausible formation histories of the dynamic large scale structure underlying the observed galaxy distribution. The results described in this work constitute the first full Bayesian non-linear analysis of the cosmic large scale structure with the demonstrated capability of uncertainty quantification. Some of these results will be made publicly available along with this work. The level of detail of inferred results and the high degree of control on observational uncertainties pave the path towards high precision chrono-cosmography, the subject of simultaneously studying the dynamics and the morphology of the inhomogeneous Universe.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Searching for oscillations in the primordial power spectrum

A small deviation from scale invariance in the form of oscillations in the primordial correlation spectra has been predicted by various cosmological models. In this paper we review a recently developed method to search for these features in the data in a more effective way. By Taylor expanding the small features around the 'background' cosmology, we have shown we are able to improve the search for these features compared to previous analyses. In this short paper we will extend that work by combining this method with a multi nested sampler. We recover our previous findings and are able to do so in 192 CPU hours. We will also briefly discuss the possibility of a long wavelength feature in the data to alleviate some tension between CMB data and the LCDM+r concordance cosmology.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Searching for Oscillations in the Primordial Power Spectrum: Constraints from Planck (Paper II)

We apply our recently developed code to search for resonance features in the Planck CMB temperature data. We search both for log spaced oscillations or linear spaced oscillations and compare our findings with results of our WMAP9 analysis and the Planck team analysis. While there are hints of log spaced resonant features present in the WMAP9 data, the significance of these features weaken with more data. With more accurate small scale measurements, we also find that the best fit frequency has shifted and the amplitude has been reduced. We confirm the presence of a several low frequency peaks, earlier identified by the Planck team, but with a better improvement of fit (delta chi^2 ~ 12). We further investigate this improvement by allowing the lensing potential to vary as well, showing mild correlation between the amplitude of the oscillations and the lensing amplitude. We find that the improvement of the fit increases even more (delta chi^2 ~ 14) for the low frequencies that modify the spectrum in a way that mimics the lensing effect. Since these features were not present in the WMAP data, they are primarily due to better measurements of Planck at small angular scales. For linear spaced oscillations we find a maximum delta chi^2 ~ 13 scanning two orders of magnitude in frequency space, and the biggest improvements are at extremely high frequencies. We recover a best fit frequency very close to the one found in WMAP9, which confirms that the fit improvement is driven by low l. Further comparisons with WMAP9 show Planck contains many more features, both for linear and log space oscillations, but with a smaller improvement of fit. We discuss the improvement as a function of the number of modes and study the effect of the 217 GHz map, which appears to drive most of the improvement for log spaced oscillations. We conclude that none of the detected features are statistically significant.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Semi-blind Bayesian inference of CMB map and power spectrum

We present a new blind formulation of the Cosmic Microwave Background (CMB) inference problem. The approach relies on a phenomenological model of the multi-frequency microwave sky without the need for physical models of the individual components. For all-sky and high resolution data, it unifies parts of the analysis that have previously been treated separately, such as component separation and power spectrum inference. We describe an efficient sampling scheme that fully explores the component separation uncertainties on the inferred CMB products such as maps and/or power spectra. External information about individual components can be incorporated as a prior giving a flexible way to progressively and continuously introduce physical component separation from a maximally blind approach. We connect our Bayesian formalism to existing approaches such as Commander, SMICA and ILC, and discuss possible future extensions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Simulation of the analysis of interferometric microwave background polarization data

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Sparse Inpainting and Isotropy

Sparse inpainting techniques are gaining in popularity as a tool for cosmological data analysis, in particular for handling data which present masked regions and missing observations. We investigate here the relationship between sparse inpainting techniques using the spherical harmonic basis as a dictionary and the isotropy properties of cosmological maps, as for instance those arising from cosmic microwave background (CMB) experiments. In particular, we investigate the possibility that inpainted maps may exhibit anisotropies in the behaviour of higher-order angular polyspectra. We provide analytic computations and simulations of inpainted maps for a Gaussian isotropic model of CMB data, suggesting that the resulting angular trispectrum may exhibit small but non-negligible deviations from isotropy.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Sparse sampling, galaxy bias, and voids

To study the impact of sparsity and galaxy bias on void statistics, we use a single large-volume, high-resolution N-body simulation to compare voids in multiple levels of subsampled dark matter, halo populations, and mock galaxies from a Halo Occupation Distribution model tuned to different galaxy survey densities. We focus our comparison on three key observational statistics: number functions, ellipticity distributions, and radial density profiles. We use the hierarchical tree structure of voids to interpret the impacts of sampling density and galaxy bias, and theoretical and empirical functions to describe the statistics in all our sample populations. We are able to make simple adjustments to theoretical expectations to offer prescriptions for translating from analytics to the void properties measured in realistic observations. We find that sampling density has a much larger effect on void sizes than galaxy bias. At lower tracer density, small voids disappear and the remaining voids are larger, more spherical, and have slightly steeper profiles. When a proper lower mass threshold is chosen, voids in halo distributions largely mimic those found in galaxy populations, except for ellipticities, where galaxy bias leads to higher values. We use the void density profile of Hamaus et al. (2014) to show that voids follow a self-similar and universal trend, allowing simple translations between voids studied in dark matter and voids identified in galaxy surveys. We have added the mock void catalogs used in this work to the Public Cosmic Void Catalog at http://www.cosmicvoids.net.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Testing cosmic geometry without dynamic distortions using voids

We propose a novel technique to probe the expansion history of the Universe based on the clustering statistics of cosmic voids. In particular, we compute their two-point statistics in redshift space on the basis of realistic mock galaxy catalogs and apply the Alcock-Paczynski test. In contrast to galaxies, we find void auto-correlations to be marginally affected by peculiar motions, providing a model-independent measure of cosmological parameters without systematics from redshift-space distortions. Because only galaxy-galaxy and void-galaxy correlations have been considered in these types of studies before, the presented method improves both statistical and systematic uncertainties on the product of angular diameter distance and Hubble rate, furnishing the potentially cleanest probe of cosmic geometry available to date.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The dark matter of galaxy voids

How do observed voids relate to the underlying dark matter distribution? To examine the spatial distribution of dark matter contained within voids identified in galaxy surveys, we apply Halo Occupation Distribution models representing sparsely and densely sampled galaxy surveys to a high-resolution N-body simulation. We compare these galaxy voids to voids found in the halo distribution, low-resolution dark matter, and high-resolution dark matter. We find that voids at all scales in densely sampled surveys - and medium- to large-scale voids in sparse surveys - trace the same underdensities as dark matter, but they are larger in radius by ~20%, they have somewhat shallower density profiles, and they have centers offset by ~0.4Rv rms. However, in void-to-void comparison we find that shape estimators are less robust to sampling, and the largest voids in sparsely sampled surveys suffer fragmentation at their edges. We find that voids in galaxy surveys always correspond to underdensities in the dark matter, though the centers may be offset. When this offset is taken into account, we recover almost identical radial density profiles between galaxies and dark matter. All mock catalogs used in this work are available at http://www.cosmicvoids.net.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Universal Density Profile for Cosmic Voids

We present a simple empirical function for the average density profile of cosmic voids, identified via the watershed technique in $Λ$CDM N-body simulations. This function is universal across void size and redshift, accurately describing a large radial range of scales around void centers with only two free parameters. In analogy to halo density profiles, these parameters describe the scale radius and the central density of voids. While we initially start with a more general four-parameter model, we find two of its parameters to be redundant, as they follow linear trends with the scale radius in two distinct regimes of the void sample, separated by its compensation scale. Assuming linear theory, we derive an analytic formula for the velocity profile of voids and find an excellent agreement with the numerical data as well. In our companion paper [Sutter et al., Mon. Not. R. Astron. Soc. 442, 462 (2014)] the presented density profile is shown to be universal even across tracer type, properly describing voids defined in halo and galaxy distributions of varying sparsity, allowing us to relate various void populations by simple rescalings. This provides a powerful framework to match theory and simulations with observational data, opening up promising perspectives to constrain competing models of cosmology and gravity.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

VIDE: The Void IDentification and Examination toolkit

We present VIDE, the Void IDentification and Examination toolkit, an open-source Python/C++ code for finding cosmic voids in galaxy redshift surveys and N-body simulations, characterizing their properties, and providing a platform for more detailed analysis. At its core, VIDE uses a substantially enhanced version of ZOBOV (Neyinck 2008) to calculate a Voronoi tessellation for estimating the density field and a performing a watershed transform to construct voids. Additionally, VIDE provides significant functionality for both pre- and post-processing: for example, vide can work with volume- or magnitude-limited galaxy samples with arbitrary survey geometries, or dark matter particles or halo catalogs in a variety of common formats. It can also randomly subsample inputs and includes a Halo Occupation Distribution model for constructing mock galaxy populations. VIDE uses the watershed levels to place voids in a hierarchical tree, outputs a summary of void properties in plain ASCII, and provides a Python API to perform many analysis tasks, such as loading and manipulating void catalogs and particle members, filtering, plotting, computing clustering statistics, stacking, comparing catalogs, and fitting density profiles. While centered around ZOBOV, the toolkit is designed to be as modular as possible and accommodate other void finders. VIDE has been in development for several years and has already been used to produce a wealth of results, which we summarize in this work to highlight the capabilities of the toolkit. VIDE is publicly available at http://bitbucket.org/cosmicvoids/vide public and http://www.cosmicvoids.net.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Voids in the SDSS DR9: observations, simulations, and the impact of the survey mask

We present and study cosmic voids identified using the watershed void finder VIDE in the Sloan Digital Sky Survey Data Release 9, compare these voids to ones identified in mock catalogs, and assess the impact of the survey mask on void statistics such as number functions, ellipticity distributions, and radial density profiles. The nearly 1,000 identified voids span three nearly volume-limited samples from redshift z = 0.43 to 0.7. For comparison we use 98 of the publicly available 2LPT-based mock galaxy catalogs of Manera et al., and also generate our own mock catalogs by applying a Halo Occupation Distribution model to an N-body simulation. We find that the mask reduces the number density of voids at all scales by a factor of three and slightly skews the relative size distributions. This engenders an increase in the mean ellipticity by roughly 30%. However, we find that radial density profiles are largely robust to the effects of the mask. We see excellent agreement between the data and both mock catalogs, and find no tension between the observed void properties and the properties derived from ΛCDM simulations. We have added the void catalogs from both data and mock galaxy populations discussed in this work to the Public Cosmic Void Catalog at http://www.cosmicvoids.net.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A response to arXiv:1310.2791: A self-consistent public catalogue of voids and superclusters in the SDSS Data Release 7 galaxy surveys

Recently, Nadathur & Hotchkiss (2013) submitted a paper discussing a new cosmic void catalog. This paper includes claims about the void catalog described in Sutter et al. (2012). In this note, we respond to those claims, clarify some discrepancies between the text of Sutter et al. (2012) and the most recent version of the catalog, and provide some comments on the differences between our catalog and that of Nadathur & Hotchkiss (2013). All updates and documentation for our catalog are available at http://www.cosmicvoids.net.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Bayesian Inference of Polarized CMB Power Spectra from Interferometric Data

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Fast Wiener filtering of CMB maps

We present the application of a new method to compute the Wiener filter solution of large and complex data sets. Contrary to the iterative solvers usually employed in signal processing, our algorithm does not require the use of preconditioners to be computationally efficient. The new scheme is conceptually very simple and therefore easy to implement, numerically absolutely stable, and guaranteed to converge. We introduce a messenger field to mediate between the different preferred bases in which signal and noise properties can be specified most conveniently, and rephrase the signal reconstruction problem in terms of this auxiliary variable. We demonstrate the capabilities of the algorithm by applying it to cosmic microwave background (CMB) radiation data obtained by the WMAP satellite.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Maximum likelihood analysis of systematic errors in interferometric observations of the cosmic microwave background

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-σ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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Methods for Bayesian power spectrum inference with galaxy surveys

We derive and implement a full Bayesian large scale structure inference method aiming at precision recovery of the cosmological power spectrum from galaxy redshift surveys. Our approach improves over previous Bayesian methods by performing a joint inference of the three dimensional density field, the cosmological power spectrum, luminosity dependent galaxy biases and corresponding normalizations. We account for all joint and correlated uncertainties between all inferred quantities. Classes of galaxies with different biases are treated as separate sub samples. The method therefore also allows the combined analysis of more than one galaxy survey. In particular, it solves the problem of inferring the power spectrum from galaxy surveys with non-trivial survey geometries by exploring the joint posterior distribution with efficient implementations of multiple block Markov chain and Hybrid Monte Carlo methods. Our Markov sampler achieves high statistical efficiency in low signal to noise regimes by using a deterministic reversible jump algorithm. We test our method on an artificial mock galaxy survey, emulating characteristic features of the Sloan Digital Sky Survey data release 7, such as its survey geometry and luminosity dependent biases. These tests demonstrate the numerical feasibility of our large scale Bayesian inference frame work when the parameter space has millions of dimensions. The method reveals and correctly treats the anti-correlation between bias amplitudes and power spectrum, which are not taken into account in current approaches to power spectrum estimation, a 20 percent effect across large ranges in k-space. In addition, the method results in constrained realizations of density fields obtained without assuming the power spectrum or bias parameters in advance.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Probabilistic image reconstruction for radio interferometers

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Searching for Oscillations in the Primordial Power Spectrum: Perturbative Approach (Paper I)

In this first of two papers, we present a new method for searching for oscillatory features in the primordial power spectrum. A wide variety of models predict these features in one of two different flavors: logarithmically spaced oscillations and linearly spaced oscillations. The proposed method treats the oscillations as perturbations on top of the scale-invariant power spectrum, allowing us to vary all cosmological parameters. This perturbative approach reduces the computational requirements for the search as the transfer functions and their derivatives can be precomputed. We show that the most significant degeneracy in the analysis is between the distance to last scattering and the overall amplitude at low frequencies. For models with logarithmic oscillations, this degeneracy leads to an uncertainty in the phase. For linear spaced oscillations, it affects the frequency of the oscillations. In this first of two papers, we test our code on simulated Planck-like data, and show we are able to recover fiducial input oscillations with an amplitude of a few times order 10^{-2}. We apply the code to WMAP9-year data and confirm the existence of two intriguing resonant frequencies for log spaced oscillations. For linear spaced oscillations we find a single resonance peak. We use numerical simulations to assess the significance of these features and conclude that the data do not provide compelling evidence for the existence of oscillatory features in the primordial spectrum.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Systematic Effects in Interferometric Observations of the CMB Polarization

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, $|ε_{rms}| = 5 \times 10^{-4}$ for antenna coupling, $δ_{rms} \approx 0.7^\circ$ for pointing, $ζ_{rms} \approx 0.7^\circ$ for beam shape, and $μ_{rms} = 5 \times 10^{-4}$ for beam cross-polarization.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A first application of the Alcock-Paczynski test to stacked cosmic voids

We report on the first application of the Alcock-Paczynski test to stacked voids in spectroscopic galaxy redshift surveys.We use voids from the Sutter et al. (2012) void catalog, which was derived from the Sloan Digital Sky Survey Data Release 7 main sample and luminous red galaxy catalogs. The construction of that void catalog removes potential shape measurement bias by using a modified version of the ZOBOV algorithm and by removing voids near survey boundaries and masks. We apply the shape-fitting procedure presented in Lavaux & Wandelt (2012) to ten void stacks out to redshift z=0.36. Combining these measurements, we determine the mean cosmologically induced "stretch" of voids in three redshift bins, with 1-sigma errors of 5-15%. The mean stretch is consistent with unity, providing no indication of a distortion induced by peculiar velocities. While the statistical errors are too large to detect the Alcock-Paczynski effect over our limited redshift range, this proof-of-concept analysis defines procedures that can be applied to larger spectroscopic galaxy surveys at higher redshifts to constrain dark energy using the expected statistical isotropy of structures that are minimally affected by uncertainties in galaxy velocity bias.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A public void catalog from the SDSS DR7 Galaxy Redshift Surveys based on the watershed transform

We produce the most comprehensive public void catalog to date using the Sloan Digital Sky Survey Data Release 7 main sample out to redshift z=0.2 and the luminous red galaxy sample out to z=0.44. Using a modified version of the parameter-free void finder ZOBOV, we fully take into account the presence of the survey boundary and masks. Our strategy for finding voids is thus appropriate for any survey configuration. We produce two distinct catalogs: a complete catalog including voids near any masks, which would be appropriate for void galaxy surveys, and a bias-free catalog of voids away from any masks, which is necessary for analyses that require a fair sampling of void shapes and alignments. Our discovered voids have effective radii from 5 to 135 h^-1 Mpc. We discuss basic catalog statistics such as number counts and redshift distributions and describe some additional data products derived from our catalog, such as radial density profiles and projected density maps. We find that radial profiles of stacked voids show a qualitatively similar behavior across nearly two decades of void radii and throughout the full redshift range.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Accelerating convolutions on the sphere with hybrid GPU/CPU kernel splitting

We present a general method for accelerating by more than an order of magnitude the convolution of pixelated function on the sphere with a radially-symmetric kernel. Our method splits the kernel into a compact real-space, and a compact spherical harmonic space component that can then be convolved in parallel using an inexpensive commodity GPU and a CPU, respectively. We provide models for the computational cost of both real-space and Fourier space convolutions and an estimate for the approximation error. Using these models we can determine the optimum split that minimizes the wall clock time for the convolution while satisfying the desired error bounds. We apply this technique to the problem of simulating a cosmic microwave background sky map at the resolution typical of the high resolution maps of the cosmic microwave background anisotropies produced by the Planck space craft. For the main Planck CMB science channels we achieve a speedup of over a factor of ten, assuming an acceptable fractional rms error of order 10^-5 in the (power spectrum of the) output map.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Bayesian angular power spectrum analysis of interferometric data

We present a Bayesian angular power spectrum and signal map inference engine which can be adapted to interferometric observations of anisotropies inthe cosmic microwave background, 21 cm emission line mapping of galactic brightness fluctuations, or 21 cm absorption line mapping of neutral hydrogen in the dark ages. The method uses Gibbs sampling to generate a sampled representation of the angular power spectrum posterior and the posterior of signal maps given a set of measured visibilities in the uv-plane. We use a mock interferometric CMB observation to demonstrate the validity of this method in the flat-sky approximation when adapted to take into account arbitrary coverage of the uv-plane, mode-mode correlations due to observations on a finite patch, and heteroschedastic visibility errors. The computational requirements scale as O(n_p log n_p) where n_p measures the ratio of the size of the detector array to the inter-detector spacing, meaning that Gibbs sampling is a promising technique for meeting the data analysis requirements of future cosmology missions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Bayesian physical reconstruction of initial conditions from large scale structure surveys

We present a fully probabilistic, physical model of the non-linearly evolved density field, as probed by realistic galaxy surveys. Our model is valid in the linear and mildly non-linear regimes and uses second order Lagrangian perturbation theory to connect the initial conditions with the final density field. Our parameter space consists of the 3D initial density field and our method allows a fully Bayesian exploration of the sets of initial conditions that are consistent with the galaxy distribution sampling the final density field. A natural byproduct of this technique is an optimal non-linear reconstruction of the present density and velocity fields, including a full propagation of the observational uncertainties. A test of these methods on simulated data mimicking the survey mask, selection function and galaxy number of the SDSS DR7 main sample shows that this physical model gives accurate reconstructions of the underlying present-day density and velocity fields on scales larger than ~6 Mpc/h. Our method naturally and accurately reconstructs non-linear features corresponding to three-point and higher order correlation functions such as walls and filaments. Simple tests of the reconstructed initial conditions show statistical consistency with the Gaussian simulation inputs. Our test demonstrates that statistical approaches based on physical models of the large scale structure distribution are now becoming feasible for realistic current and future surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Efficient Wiener filtering without preconditioning

We present a new approach to calculate the Wiener filter solution of general data sets. It is trivial to implement, flexible, numerically absolutely stable, and guaranteed to converge. Most importantly, it does not require an ingenious choice of preconditioner to work well. The method is capable of taking into account inhomogeneous noise distributions and arbitrary mask geometries. It iteratively builds up the signal reconstruction by means of a messenger field, introduced to mediate between the different preferred bases in which signal and noise properties can be specified most conveniently. Using cosmic microwave background (CMB) radiation data as a showcase, we demonstrate the capabilities of our scheme by computing Wiener filtered WMAP7 temperature and polarization maps at full resolution for the first time. We show how the algorithm can be modified to synthesize fluctuation maps, which, combined with the Wiener filter solution, result in unbiased constrained signal realizations, consistent with the observations. The algorithm performs well even on simulated CMB maps with Planck resolution and dynamic range.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Fast calculation of the Fisher matrix for cosmic microwave background experiments

The Fisher information matrix of the cosmic microwave background (CMB) radiation power spectrum coefficients is a fundamental quantity that specifies the information content of a CMB experiment. In the most general case, its exact calculation scales with the third power of the number of data points N and is therefore computationally prohibitive for state-of-the-art surveys. Applicable to a very large class of CMB experiments without special symmetries, we show how to compute the Fisher matrix in only O(N^2 log N) operations as long as the inverse noise covariance matrix can be applied to a data vector in time O(l_max^3 log l_max). This assumption is true to a good approximation for all CMB data sets taken so far. The method takes into account common systematics such as arbitrary sky coverage and realistic noise correlations. As a consequence, optimal quadratic power spectrum estimation also becomes feasible in O(N^2 log N) operations for this large group of experiments. We discuss the relevance of our findings to other areas of cosmology where optimal power spectrum estimation plays a role.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Likelihood, Fisher information, and systematics of cosmic microwave background experiments

Every experiment is affected by systematic effects that hamper the data analysis and have the potential to ultimately degrade its performance. In the case of probes of the cosmic microwave background (CMB) radiation, a minimal set of issues to consider includes asymmetric beam functions, correlated noise, and incomplete sky coverage. Presuming a simplified scanning strategy that allows for an exact analytical treatment of the problem, we study the impact of systematic effects on the likelihood function of the CMB power spectrum. We use the Fisher matrix, a measure of the information content of a data set, for a quantitative comparison of different experimental configurations. In addition, for various power spectrum coefficients, we explore the functional form of the likelihood directly, and obtain the following results: The likelihood function can deviate systematically from a Gaussian distribution up to the highest multipole values considered in our analysis. Treated exactly, realistic levels of asymmetric beam functions and correlated noise do not by themselves decrease the information yield of CMB experiments nor do they induce noticeable coupling between multipoles. Masking large fractions of the sky, on the other hand, results in a considerably more complex correlation structure of the likelihood function. Combining adjacent power spectrum coefficients into bins can partially mitigate these problems.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Precision cosmography with stacked voids

We present a purely geometrical method for probing the expansion history of the Universe from the observation of the shape of stacked voids in spectroscopic redshift surveys. Our method is an Alcock-Paczynski (AP) test based on the average sphericity of voids posited on the local isotropy of the Universe. It works by comparing the temporal extent of cosmic voids along the line of sight with their angular, spatial extent. We describe the algorithm that we use to detect and stack voids in redshift shells on the light cone and test it on mock light cones produced from N-body simulations. We establish a robust statistical model for estimating the average stretching of voids in redshift space and quantify the contamination by peculiar velocities. Finally, assuming that the void statistics that we derive from N-body simulations is preserved when considering galaxy surveys, we assess the capability of this approach to constrain dark energy parameters. We report this assessment in terms of the figure of merit (FoM) of the dark energy task force and in particular of the proposed EUCLID mission which is particularly suited for this technique since it is a spectroscopic survey. The FoM due to stacked voids from the EUCLID wide survey may double that of all other dark energy probes derived from EUCLID data alone (combined with Planck priors). In particular, voids seem to outperform Baryon Acoustic Oscillations by an order of magnitude. This result is consistent with simple estimates based on mode-counting. The AP test based on stacked voids may be a significant addition to the portfolio of major dark energy probes and its potentialities must be studied in detail.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Robust, data-driven inference in non-linear cosmostatistics

We discuss two projects in non-linear cosmostatistics applicable to very large surveys of galaxies. The first is a Bayesian reconstruction of galaxy redshifts and their number density distribution from approximate, photometric redshift data. The second focuses on cosmic voids and uses them to construct cosmic spheres that allow reconstructing the expansion history of the Universe using the Alcock-Paczynski test. In both cases we find that non-linearities enable the methods or enhance the results: non-linear gravitational evolution creates voids and our photo-z reconstruction works best in the highest density (and hence most non-linear) portions of our simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Bayesian inference from photometric redshift surveys

We show how to enhance the redshift accuracy of surveys consisting of tracers with highly uncertain positions along the line of sight. Photometric surveys with redshift uncertainty delta_z ~ 0.03 can yield final redshift uncertainties of delta_z_f ~ 0.003 in high density regions. This increased redshift precision is achieved by imposing an isotropy and 2-point correlation prior in a Bayesian analysis and is completely independent of the process that estimates the photometric redshift. As a byproduct, the method also infers the three dimensional density field, essentially super-resolving high density regions in redshift space. Our method fully takes into account the survey mask and selection function. It uses a simplified Poissonian picture of galaxy formation, relating preferred locations of galaxies to regions of higher density in the matter field. The method quantifies the remaining uncertainties in the three dimensional density field and the true radial locations of galaxies by generating samples that are constrained by the survey data. The exploration of this high dimensional, non-Gaussian joint posterior is made feasible using multiple-block Metropolis-Hastings sampling. We demonstrate the performance of our implementation on a simulation containing 2.0 x 10^7 galaxies. These results bear out the promise of Bayesian analysis for upcoming photometric large scale structure surveys with tens of millions of galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Probability Distribution for Non-Gaussianity Estimators

One of the principle efforts in cosmic microwave background (CMB) research is measurement of the parameter fnl that quantifies the departure from Gaussianity in a large class of non-minimal inflationary (and other) models. Estimators for fnl are composed of a sum of products of the temperatures in three different pixels in the CMB map. Since the number ~Npix^2 of terms in this sum exceeds the number Npix of measurements, these ~Npix^2 terms cannot be statistically independent. Therefore, the central-limit theorem does not necessarily apply, and the probability distribution function (PDF) for the fnl estimator does not necessarily approach a Gaussian distribution for N_pix >> 1. Although the variance of the estimators is known, the significance of a measurement of fnl depends on knowledge of the full shape of its PDF. Here we use Monte Carlo realizations of CMB maps to determine the PDF for two minimum-variance estimators: the standard estimator, constructed under the null hypothesis (fnl=0), and an improved estimator with a smaller variance for |fnl| > 0. While the PDF for the null-hypothesis estimator is very nearly Gaussian when the true value of fnl is zero, the PDF becomes significantly non-Gaussian when |fnl| > 0. In this case we find that the PDF for the null-hypothesis estimator fnl_hat is skewed, with a long non-Gaussian tail at fnl_hat > |fnl| and less probability at fnl_hat < |fnl| than in the Gaussian case. We provide an analytic fit to these PDFs. On the other hand, we find that the PDF for the improved estimator is nearly Gaussian for observationally allowed values of fnl. We discuss briefly the implications for trispectrum (and other higher-order correlation) estimators.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Fast CMB lensing using statistical interpolation on the sphere

We describe a accurate and fast pixel-based statistical method to interpolate fields of arbitrary spin on the sphere. We call this method Fast and Lean Interpolation on the Sphere (FLINTS). The method predicts the optimal interpolated values based on the theory of isotropic Gaussian random fields and provides an accurate error estimate at no additional cost. We use this method to compute lensed Cosmic Microwave Background (CMB) maps precisely and quickly, achieving a relative precision of 0.02% at a HEALPix resolution of Nside=4096, for a bandlimit of l_max=4096 in the same time it takes to simulate the original, unlensed CMB map. The method is suitable for efficient, distributed memory parallelization. The power spectra of our lensed maps are accurate to better than 0.5% at l=3000 for the temperature, the E and B mode of the polarization. As expected theoretically, we demonstrate that, on realistic cases, this method is between two to three orders of magnitude more precise than other known interpolation methods for the same computational cost.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Local non-Gaussianity in the Cosmic Microwave Background the Bayesian way

We introduce an exact Bayesian approach to search for non-Gaussianity of local type in Cosmic Microwave Background (CMB) radiation data. Using simulated CMB temperature maps, the newly developed technique is compared against the conventional frequentist bispectrum estimator. Starting from the joint probability distribution, we obtain analytic expressions for the conditional probabilities of the primordial perturbations given the data, and for the level of non-Gaussianity, f_nl, given the data and the perturbations. We propose Hamiltonian Monte Carlo sampling as a means to derive realizations of the primordial fluctuations from which we in turn sample f_nl. Although being computationally expensive, this approach allows us to exactly construct the full target posterior probability distribution. When compared to the frequentist estimator, applying the Bayesian method to Gaussian CMB maps provides consistent results. For the analysis of non-Gaussian maps, however, the error bars on f_nl do not show excess variance within the Bayesian framework. This finding is of particular relevance in the light of upcoming high precision CMB measurements obtained by the Planck satellite mission.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

More on crinkles in the last scattering surface

Inhomogeneous recombination can give rise to perturbations in the electron number density which can be a factor of five larger than the perturbations in baryon density. We do a thorough analysis of the second order anisotropies generated in the cosmic microwave background (CMB) due to perturbations in the electron number density. We show that solving the second order Boltzmann equation for photons is equivalent to solving the first + second order Boltzmann equations and then taking the second order part of the solution. We find the approximate solution to the photon Boltzmann hierarchy in l modes and show that the contributions from inhomogeneous recombination to the second order monopole, dipole and quadrupole are numerically small. We also point out that perturbing the electron number density in the first order tight coupling and damping solutions for the monopole, dipole and quadrupole is not equivalent to solving the second order Boltzmann equations for inhomogeneous recombination. Finally we confirm our result in a previous paper that inhomogeneous recombination gives rise to a local type non-Gaussianity parameter f_{NL}~ -1. The signal to noise for the detection of the temperature bispectrum generated by inhomogeneous recombination is ~ 1 for an ideal full sky experiment measuring modes up to l_{max}=2500.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

O-V-S-Z and friends: Non-Gaussianity from inhomogeneous reionization

We calculate the cosmic microwave background (CMB) bispectrum due to inhomogeneous reionization. We calculate all the terms that can contribute to the bispectrum that are products of first order terms on all scales in conformal Newtonian gauge. We also correctly account for the de-correlation between the matter density and initial conditions using perturbation theory up to third order. We find that the bispectrum is of local type as expected. For a reasonable model of reionization, in which the Universe is completely ionized by redshift z_{ri} ~ 8 with optical depth to the last scattering surface τ_0=0.087 the signal to noise for detection of the CMB temperature bispectrum is S/N ~ 0.1 and confusion in the estimation of primordial non-Gaussianity is f_{NL} ~ -0.1. For an extreme model with z_{ri} ~ 12.5, τ_0=0.14 we get S/N ~ 0.5 and f_{NL} ~ -0.2.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Precision cosmology with voids: definition, methods, dynamics

We propose a new definition of cosmic voids based on methods of Lagrangian orbit reconstruction as well as an algorithm to find them in actual data called DIVA. Our technique is intended to yield results which can be modeled sufficiently accurately to create a new probe of precision cosmology. We then develop an analytical model of the ellipticity of voids found by our method based on Zel'dovich approximation. We measure in N-body simulation that this model is precise at the 0.1% level for the mean ellipticity of voids of size greater than ~4 Mpc/h. We estimate that at this scale, we are able to predict the ellipticity with an accuracy of 0.02. Finally, we compare the distribution of void shapes in N-body simulation for two different equations of state w of the dark energy. We conclude that our method is far more accurate than Eulerian methods and is therefore promising as a precision probe of dark energy phenomenology.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Primordial Non-Gaussianity in the Cosmic Microwave Background

In the last few decades, advances in observational cosmology have given us a standard model of cosmology. We know the content of the universe to within a few percent. With more ambitious experiments on the way, we hope to move beyond the knowledge of what the universe is made of, to why the universe is the way it is. In this review paper we focus on primordial non-Gaussianity as a probe of the physics of the dynamics of the universe at the very earliest moments. We discuss 1) theoretical predictions from inflationary models and their observational consequences in the cosmic microwave background (CMB) anisotropies; 2) CMB--based estimators for constraining primordial non-Gaussianity with an emphasis on bispectrum templates; 3) current constraints on non-Gaussianity and what we can hope to achieve in the near future; and 4) non-primordial sources of non-Gaussianities in the CMB such as bispectrum due to second order effects, three way cross-correlation between primary-lensing-secondary CMB, and possible instrumental effects.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Probing local non-Gaussianities within a Bayesian framework

Aims: We outline the Bayesian approach to inferring f_NL, the level of non-Gaussianity of local type. Phrasing f_NL inference in a Bayesian framework takes advantage of existing techniques to account for instrumental effects and foreground contamination in CMB data and takes into account uncertainties in the cosmological parameters in an unambiguous way. Methods: We derive closed form expressions for the joint posterior of f_NL and the reconstructed underlying curvature perturbation, Phi, and deduce the conditional probability densities for f_NL and Phi. Completing the inference problem amounts to finding the marginal density for f_NL. For realistic data sets the necessary integrations are intractable. We propose an exact Hamiltonian sampling algorithm to generate correlated samples from the f_NL posterior. For sufficiently high signal-to-noise ratios, we can exploit the assumption of weak non-Gaussianity to find a direct Monte Carlo technique to generate independent samples from the posterior distribution for f_NL. We illustrate our approach using a simplified toy model of CMB data for the simple case of a 1-D sky. Results: When applied to our toy problem, we find that, in the limit of high signal-to-noise, the sampling efficiency of the approximate algorithm outperforms that of Hamiltonian sampling by two orders of magnitude. When f_NL is not significantly constrained by the data, the more efficient, approximate algorithm biases the posterior density towards f_NL = 0.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The primordial non-Gaussianity of local type (f_NL) in the WMAP 5-year data: the length distribution of CMB skeleton

We present skeleton studies of non-Gaussianity in the CMB temperature anisotropy observed in the WMAP5 data. The local skeleton is traced on the 2D sphere by cubic spline interpolation which leads to more accurate estimation of the intersection positions between the skeleton and the secondary pixels than conventional linear interpolation. We demonstrate that the skeleton-based estimator of non-Gaussianity of the local type (f_NL) - the departure of the length distribution from the corresponding Gaussian expectation - yields an unbiased and sufficiently converged f_NL-likelihood. We analyse the skeleton statistics in the WMAP5 combined V- and W-band data outside the Galactic base-mask determined from the KQ75 sky-coverage. The results are consistent with Gaussian simulations of the the best-fitting cosmological model, but deviate from the previous results determined using the WMAP1 data. We show that it is unlikely that the improved skeleton tracing method, the omission of Q-band data, the modification of the foreground-template fitting method or the absence of 6 extended regions in the new mask contribute to such a deviation. However, the application of the Kp0 base-mask in data processing does improve the consistency with the WMAP1 results. The f_NL-likelihoods of the data are estimated at 9 different smoothing levels. It is unexpected that the best-fit values show positive correlation with the smoothing scales. Further investigation argues against a point-source or goodness-of-fit explanation but finds that about 30% of either Gaussian or f_NL samples having better goodness-of-fit than the WMAP5 show a similar correlation. We present the estimate f_NL=47.3+/-34.9 (1sigma error) determined from the first four smoothing angles and f_NL=76.8+/-43.1 for the combination of all nine. The former result may be overestimated at the 0.21sigma-level because of point sources.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Voids as a Precision Probe of Dark Energy

A signature of the dark energy equation of state may be observed in the shape of voids. We estimate the constraints on cosmological parameters that would be determined from the ellipticity distribution of voids from future spectroscopic surveys already planned for the study of large scale structure. The constraints stem from the sensitivity of the distribution of ellipticity to the cosmological parameters through the variance of fluctuations of the density field smoothed at some length scale. This length scale can be chosen to be of the order of the comoving radii of voids at very early times when the fluctuations are Gaussian distributed. We use Fisher estimates to show that the constraints from void ellipticities are promising. Combining these constraints with other traditional methods results in the improvement of the Dark Energy Task Force Figure of Merit on the dark energy parameters by an order of hundred for future experiments. The estimates of these future constraints depend on a number of systematic issues which require further study using simulations. We outline these issues and study the impact of certain observational and theoretical systematics on the forecasted constraints on dark energy parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

21 cm radiation: A new probe of fundamental physics

New low frequency radio telescopes currently being built open up the possibility of observing the 21-cm radiation before the Epoch of Reionization in the future, in particular at redshifts 200 > z > 30, also known as the dark ages. At these high redshifts, Cosmic Microwave Background (CMB) radiation is absorbed by neutral hydrogen at its 21-cm hyperfine transition. This redshifted 21-cm signal thus carries information about the state of the early Universe and can be used to test fundamental physics. We study the constraints these observations can put on the variation of fundamental constants. We show that the 21-cm radiation is very sensitive to the variations in the fine structure constant and can in principle place constraints comparable to or better than the other astrophysical experiments (fractional change < 10^ {-5}). Making such observations will require radio telescopes of collecting area 10 - 10^6 sq. km compared to 1 sq. km of current telescopes. These observations will thus provide independent constraints on the fine structure constant at high redshifts, observations of quasars being the only alternative. More importantly the 21-cm absorption of CMB is the only way to probe the redshift range between recombination and reionization.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Cosmic Cartography of the Large-Scale Structure with Sloan Digital Sky Survey Data Release 6

We present the largest Wiener reconstruction of the cosmic density field made to date. The reconstruction is based on the Sloan Digital Sky Survey data release 6 covering the northern Galactic cap. We use a novel supersampling algorithm to suppress aliasing effects and a Krylov-space inversion method to enable high performance with high resolution. These techniques are implemented in the ARGO computer code. We reconstruct the field over a 500 Mpc cube with Mpc grid-resolution while accounting both for the angular and radial selection functions of the SDSS, and the shot noise giving an effective resolution of the order of ~10 Mpc. In addition, we correct for the redshift distortions in the linear and nonlinear regimes in an approximate way. We show that the commonly used method of inverse weighting the galaxies by the corresponding selection function heads to excess noise in regions where the density of the observed galaxies is small. It is more accurate and conservative to adopt a Bayesian framework in which we model the galaxy selection/detection process to be Poisson-binomial. This results in heavier smoothing in regions of reduced sampling density. Our results show a complex cosmic web structure with huge void regions indicating that the recovered matter distribution is highly non-Gaussian. Filamentary structures are clearly visible on scales up to ~20 Mpc. We also calculate the statistical distribution of density after smoothing the reconstruction with Gaussian kernels of different radii r_S and find good agreement with a log-normal distribution for ~10 Mpc < r_S < ~30 Mpc.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Crinkles in the last scattering surface: Non-Gaussianity from inhomogeneous recombination

The perturbations in the electron number density during recombination contributes to the Cosmic Microwave Background bispectrum through second order terms. Perturbations in the electron density can be a factor of ~5 larger than the baryon density fluctuations on large scales as shown in the calculations by Novosyadlyj. This raises the possibility that the contribution to bispectrum arising from perturbations in the optical depth may be non-negligible. We calculate this bispectrum and find it to peak for squeezed triangles and of peak amplitude of the order of primordial non-Gaussianity of local type with fNL of 0.05 ~ -1 depending on the l-modes being considered. This is because the shape of the bispectrum is different from the primordial one although it peaks for squeezed configurations, similar to the local type primordial non-Gaussianity.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

LSST Science Book, Version 2.0

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Cosmic (super)string constraints from 21 cm radiation

We calculate the contribution of cosmic strings arising from a phase transition in the early universe, or cosmic superstrings arising from brane inflation, to the cosmic 21 cm power spectrum at redshifts z > 30. Future experiments can exploit this effect to constrain the cosmic string tension Gu and probe virtually the entire brane inflation model space allowed by current observations. Although current experiments with a collecting area of ~ 1 km^2 will not provide any useful constraints, future experiments with a collecting area of 10^4-10^6 km^2 covering the cleanest 10% of the sky can in principle constrain cosmic strings with tension Gu > 10^(-10) to 10^(-12) (superstring/phase transition mass scale >10^13 GeV).

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Fast Estimator of Primordial Non-Gaussianity from Temperature and Polarization Anisotropies in the Cosmic Microwave Background II: Partial Sky Coverage and Inhomogeneous Noise

In our recent paper (Yadav et al. 2007) we described a fast cubic (bispectrum) estimator of the amplitude of primordial non-Gaussianity of local type, f_{NL}, from a combined analysis of the Cosmic Microwave Background (CMB) temperature and E-polarization observations. In this paper we generalize the estimator to deal with a partial sky coverage as well as inhomogeneous noise. Our generalized estimator is still computationally efficient, scaling as O(N^3/2) compared to the O(N^5/2) scaling of the brute force bispectrum calculation for sky maps with N pixels. Upcoming CMB experiments are expected to yield high-sensitivity temperature and E-polarization data. Our generalized estimator will allow us to optimally utilize the combined CMB temperature and E-polarization information from these realistic experiments, and to constrain primordial non-Gaussianity.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Interacting Quintessence, Cosmic Acceleration and the Coincidence Problem

Faced by recent evidence for a flat universe dominated by dark energy, cosmologists grapple with deep cosmic enigmas such as the cosmological constant problem, extreme fine-tuning and the cosmic coincidence problem. The extent to which we observe the dimming of distant supernovae suggests that the cosmic acceleration is as least as severe as in cosmological constant models. Extrapolating this to our cosmic future implies terrifying visions of either a cold and empty universe or an explosive demise in a ``Big Rip.'' We construct a class of dynamical scalar field models of dark energy and dark matter. Within this class we can explain why supernovae imply a cosmic equation of state $w\lesssim-1$, address fine tuning issues, protect the universe from premature acceleration and predict a constant fraction of dark energy to dark matter in the future (thus solving the coincidence problem), satisfy the dominant energy condition, and ensure that gravitationally bound objects remain so forever (avoid a Big Rip). This is achieved with a string theory inspired Lagrangian containing standard kinetic terms, exponential potentials and couplings, and parameters of order unity.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

MAGIC: Exact Bayesian Covariance Estimation and Signal Reconstruction for Gaussian Random Fields

In this talk I describe MAGIC, an efficient approach to covariance estimation and signal reconstruction for Gaussian random fields (MAGIC Allows Global Inference of Covariance). It solves a long-standing problem in the field of cosmic microwave background (CMB) data analysis but is in fact a general technique that can be applied to noisy, contaminated and incomplete or censored measurements of either spatial or temporal Gaussian random fields. In this talk I will phrase the method in a way that emphasizes its general structure and applicability but I comment on applications in the CMB context. The method allows the exploration of the full non-Gaussian joint posterior density of the signal and parameters in the covariance matrix (such as the power spectrum) given the data. It generalizes the familiar Wiener filter in that it automatically discovers signal correlations in the data as long as a noise model is specified and priors encode what is known about potential contaminants. The key methodological difference is that instead of attempting to evaluate the likelihood (or posterior density) or its derivatives, this method generates an asymptotically exact Monte Carlo sample from it. I present example applications to power spectrum estimation and signal reconstruction from measurements of the CMB. For these applications the method achieves speed-ups of many orders of magnitude compared to likelihood maximization techniques, while offering greater flexibility in modeling and a full characterization of the uncertainty in the estimates.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Statistical Challenges of Cosmic Microwave Background Analysis

The Cosmic Microwave Background (CMB) is an abundant source of cosmological information. However, this information is encoded in non-trivial ways in a signal that is difficult to observe. The resulting challenges in extracting this information from CMB data sets have created a new frontier. In this talk I will discuss the challenges of CMB data analysis. I review what cosmological information is contained in the CMB data and the problem of extracting it. CMB analyses can be divided into two types: ``canonical'' parameter extraction which seeks to obtain the best possible estimates of cosmological parameters within a pre-defined theory space and "hypothesis testing" which seeks to test the assumption on which the canonical tests rest. Both of these activities are fundamentally important. In addition to mining the CMB for cosmological information cosmologists would like to strengthen the analysis with data from other cosmologically interesting observations as well as physical constraints. This gives an opportunity 1) to test the results from these separate probes for concordance and 2) if concordance is established to sharpen the constraints on theory space by combining the information from these separate sources.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

Halo Properties in Cosmological Simulations of Self-Interacting Cold Dark Matter

We present a comparison of halo properties in cosmological simulations of collisionless cold dark matter (CDM) and self-interacting dark matter (SIDM) for a range of dark matter cross sections. We find, in agreement with various authors, that CDM yields cuspy halos that are too centrally concentrated as compared to observations. Conversely, SIDM simulations using a Monte Carlo N-body technique produce halos with significantly reduced central densities and flatter cores with increasing cross section. We introduce a concentration parameter based on enclosed mass that we expect will be straightforward to determine observationally, unlike that of Navarro, Frenk & White, and provide predictions for SIDM and CDM. SIDM also produces more spherical halos than CDM, providing possibly the strongest observational test of SIDM. We discuss our findings in relation to various relevant observations as well as SIDM simulations of other groups. Taking proper account of simulation limitations, we find that a dark matter cross section per unit mass of sigma_DM ~= 10^{-23}-10^{-24} cm^2/GeV is consistent with all current observational constraints.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

The Pseudo-$C_l$ method: Cosmic microwave background anisotropy power spectrum statistics for high precision cosmology

As the era of high precision cosmology approaches, the empirically determined power spectrum of the microwave background anisotropy $C_l$ will provide a crucial test for cosmological theories. We present an exact semi-analytic framework for the study of the ampling statistics of the $C_l$ resulting from observations with partial sky coverage and anisotropic noise distributions. This includes space-borne, air-borne and ground-based experiments. We apply this theory to demonstrate its power for constructing fast but unbiased approximate methods for the joint estimation of cosmological parameters. Further applications, such as a test for possible non-Gaussianity of the underlying theory and a ``poor man's power spectrum estimator'' are suggested. An appendix derives recursion relations for the efficient computation of the couplings between spherical harmonics on the cut sky.

0 citations0 reviewsNo code linkedOpen access