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

Hiranya V. Peiris

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

ResearcherUnclaimed source record
astro-ph.COhep-thastro-ph.IMhep-phgr-qcastro-phastro-ph.GAcond-mat.quant-gasInformation Theorymath.ITArtificial IntelligenceMachine Learning

Catalog footprint

What is connected

62works
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

62 published item(s)

preprint2024arXiv

Analog vacuum decay from vacuum initial conditions

Ultracold atomic gases can undergo phase transitions that mimic relativistic vacuum decay, allowing us to empirically test early-Universe physics in tabletop experiments. We investigate the physics of these analog systems, going beyond previous analyses of the classical equations of motion to study quantum fluctuations in the cold-atom false vacuum. We show that the fluctuation spectrum of this vacuum state agrees with the usual relativistic result in the regime where the classical analogy holds, providing further evidence for the suitability of these systems for studying vacuum decay. Using a suite of semiclassical lattice simulations, we simulate bubble nucleation from this analog vacuum state in a 1D homonuclear potassium-41 mixture, finding qualitative agreement with instanton predictions. We identify realistic parameters for this system that will allow us to study vacuum decay with current experimental capabilities, including a prescription for efficiently scanning over decay rates, and show that this setup will probe the quantum (rather than thermal) decay regime at temperatures $T\lesssim10\,\mathrm{nK}$. Our results help lay the groundwork for using upcoming cold-atom experiments as a new probe of nonperturbative early-Universe physics.

0 citations0 reviewsNo code linkedOpen access
preprint2024arXiv

Towards Accurate Field-Level Inference of Massive Cosmic Structures

We investigate the accuracy requirements for field-level inference of cluster and void masses using data from galaxy surveys. We introduce a two-step framework that takes advantage of the fact that cluster masses are determined by flows on larger scales than the clusters themselves. First, we determine the integration accuracy required to perform field-level inference of cosmic initial conditions on these large scales, by fitting to late-time galaxy counts using the Bayesian Origin Reconstruction from Galaxies (BORG) algorithm. A 20-step COLA integrator is able to accurately describe the density field surrounding the most massive clusters in the Local Super-Volume ($<135\,h^{-1}\,\mathrm{Mpc}$), but does not by itself lead to converged virial mass estimates. Therefore we carry out `posterior resimulations', using full $N$-body dynamics while sampling from the inferred initial conditions, and thereby obtain estimates of masses for nearby massive clusters. We show that these are in broad agreement with existing estimates, and find that mass functions in the Local Super-Volume are compatible with $Λ$CDM.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Discovering the building blocks of dark matter halo density profiles with neural networks

The density profiles of dark matter halos are typically modeled using empirical formulae fitted to the density profiles of relaxed halo populations. We present a neural network model that is trained to learn the mapping from the raw density field containing each halo to the dark matter density profile. We show that the model recovers the widely-used Navarro-Frenk-White (NFW) profile out to the virial radius, and can additionally describe the variability in the outer profile of the halos. The neural network architecture consists of a supervised encoder-decoder framework, which first compresses the density inputs into a low-dimensional latent representation, and then outputs $ρ(r)$ for any desired value of radius $r$. The latent representation contains all the information used by the model to predict the density profiles. This allows us to interpret the latent representation by quantifying the mutual information between the representation and the halos' ground-truth density profiles. A two-dimensional representation is sufficient to accurately model the density profiles up to the virial radius; however, a three-dimensional representation is required to describe the outer profiles beyond the virial radius. The additional dimension in the representation contains information about the infalling material in the outer profiles of dark matter halos, thus discovering the splashback boundary of halos without prior knowledge of the halos' dynamical history.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

New limits on light dark matter - proton cross section from the cosmic large-scale structure

We set the strongest limits to-date on the velocity-independent dark matter (DM) - proton cross section $σ$ for DM masses $m = 10\,\mathrm{keV}$ to $100\,\mathrm{GeV}$, using large-scale structure traced by the Lyman-alpha forest: e.g., a 95% lower limit $σ< 6 \times 10^{-30}\,\mathrm{cm}^2$, for $m = 100\,\mathrm{keV}$. Our results complement direct detection, which has limited sensitivity to sub-GeV DM. We use an emulator of cosmological simulations, combined with data from the smallest cosmological scales used to-date, to model and search for the imprint of primordial DM-proton collisions. Cosmological bounds are improved by up to a factor of 25.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Angular momentum evolution can be predicted from cosmological initial conditions

The angular momentum of dark matter haloes controls their spin magnitude and orientation, which in turn influences the galaxies therein. However, the process by which dark matter haloes acquire angular momentum is not fully understood; in particular, it is unclear whether angular momentum growth is stochastic. To address this question, we extend the genetic modification technique to allow control over the angular momentum of any region in the initial conditions. Using this technique to produce a sequence of modified simulations, we can then investigate whether changes to the angular momentum of a specified region in the evolved universe can be accurately predicted from changes in the initial conditions alone. We find that the angular momentum in regions with modified initial conditions can be predicted between 2 and 4 times more accurately than expected from applying tidal torque theory. This result is masked when analysing the angular momentum of haloes, because particles in the outskirts of haloes dominate the angular momentum budget. We conclude that the angular momentum of Lagrangian patches is highly predictable from the initial conditions, with apparent chaotic behaviour being driven by stochastic changes to the arbitrary boundary defining the halo.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Considerations for optimizing photometric classification of supernovae from the Rubin Observatory

The Vera C. Rubin Observatory will increase the number of observed supernovae (SNe) by an order of magnitude; however, it is impossible to spectroscopically confirm the class for all the SNe discovered. Thus, photometric classification is crucial but its accuracy depends on the not-yet-finalized observing strategy of Rubin Observatory's Legacy Survey of Space and Time (LSST). We quantitatively analyze the impact of the LSST observing strategy on SNe classification using simulated multi-band light curves from the Photometric LSST Astronomical Time-Series Classification Challenge (PLAsTiCC). First, we augment the simulated training set to be representative of the photometric redshift distribution per supernovae class, the cadence of observations, and the flux uncertainty distribution of the test set. Then we build a classifier using the photometric transient classification library snmachine, based on wavelet features obtained from Gaussian process fits, yielding similar performance to the winning PLAsTiCC entry. We study the classification performance for SNe with different properties within a single simulated observing strategy. We find that season length is important, with light curves of 150 days yielding the highest performance. Cadence also has an important impact on SNe classification; events with median inter-night gap <3.5 days yield higher classification performance. Interestingly, we find that large gaps (>10 days) in light curve observations do not impact performance if sufficient observations are available on either side, due to the effectiveness of the Gaussian process interpolation. This analysis is the first exploration of the impact of observing strategy on photometric supernova classification with LSST.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

General framework for cosmological dark matter bounds using $N$-body simulations

We present a general framework for obtaining robust bounds on the nature of dark matter using cosmological $N$-body simulations and Lyman-alpha forest data. We construct an emulator of hydrodynamical simulations, which is a flexible, accurate and computationally-efficient model for predicting the response of the Lyman-alpha forest flux power spectrum to different dark matter models, the state of the intergalactic medium (IGM) and the primordial power spectrum. The emulator combines a flexible parameterization for the small-scale suppression in the matter power spectrum arising in "non-cold" dark matter models, with an improved IGM model. We then demonstrate how to optimize the emulator for the case of ultra-light axion dark matter, presenting tests of convergence. We also carry out cross-validation tests of the accuracy of flux power spectrum prediction. This framework can be optimized for the analysis of many other dark matter candidates, e.g., warm or interacting dark matter. Our work demonstrates that a combination of an optimized emulator and cosmological "effective theories," where many models are described by a single set of equations, is a powerful approach for robust and computationally-efficient inference from the cosmic large-scale structure.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Strong bound on canonical ultra-light axion dark matter from the Lyman-alpha forest

We present a new bound on the ultra-light axion (ULA) dark matter mass $m_\text{a}$, using the Lyman-alpha forest to look for suppressed cosmic structure growth: a 95% lower limit $m_\text{a} > 2 \times 10^{-20}\,\text{eV}$. This strongly disfavors ($> 99.7\%$ credibility) the canonical ULA with $10^{-22}\,\text{eV} < m_\text{a} < 10^{-21}\,\text{eV}$, motivated by the string axiverse and solutions to possible tensions in the cold dark matter model. We strengthen previous equivalent bounds by about an order of magnitude. We demonstrate the robustness of our results using an optimized emulator of improved hydrodynamical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Target Neutrino Mass Precision for Determining the Neutrino Hierarchy

Recent works combining neutrino oscillation and cosmological data to determine the neutrino hierarchy found a range of odds in favour of the normal hierarchy. These results arise from differing approaches to incorporating prior knowledge about neutrinos. We develop a hierarchy-agnostic prior and show that the hierarchy cannot be conclusively determined with current data. The determination of the hierarchy is limited by the neutrino mass scale $Σ_ν$ measurement. We obtain a target precision of $σ(Σ_ν) = 0.014$ eV, necessary for conclusively establishing the normal hierarchy with future data.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Nonlinear Dynamics of the Cold Atom Analog False Vacuum

We investigate the nonlinear dynamics of cold atom systems that can in principle serve as quantum simulators of false vacuum decay. The analog false vacuum manifests as a metastable vacuum state for the relative phase in a two-species Bose-Einstein condensate (BEC), induced by a driven periodic coupling between the two species. In the appropriate low energy limit, the evolution of the relative phase is approximately governed by a relativistic wave equation exhibiting true and false vacuum configurations. In previous work, a linear stability analysis identified exponentially growing short-wavelength modes driven by the time-dependent coupling. These modes threaten to destabilize the analog false vacuum. Here, we employ numerical simulations of the coupled Gross-Pitaevski equations (GPEs) to determine the non-linear evolution of these linearly unstable modes. We find that unless a physical mechanism modifies the GPE on short length scales, the analog false vacuum is indeed destabilized. We briefly discuss various physically expected corrections to the GPEs that may act to remove the exponentially unstable modes. To investigate the resulting dynamics in cases where such a removal mechanism exists, we implement a hard UV cutoff that excludes the unstable modes as a simple model for these corrections. We use this to study the range of phenomena arising from such a system. In particular, we show that by modulating the strength of the time-dependent coupling, it is possible to observe the crossover between a second and first order phase transition out of the false vacuum.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A framework for testing isotropy with the cosmic microwave background

We present a new framework for testing the isotropy of the Universe using cosmic microwave background data, building on the nested-sampling ANICOSMO code. Uniquely, we are able to constrain the scalar, vector and tensor degrees of freedom alike; previous studies only considered the vector mode (linked to vorticity). We employ Bianchi type VII$_h$ cosmologies to model the anisotropic Universe, from which other types may be obtained by taking suitable limits. In a separate development, we improve the statistical analysis by including the effect of Bianchi power in the high-$\ell$, as well as the low-$\ell$, likelihood. To understand the effect of all these changes, we apply our new techniques to WMAP data. We find no evidence for anisotropy, constraining shear in the vector mode to $(σ_V/H)_0 < 1.7 \times 10^{-10}$ (95% CL). For the first time, we place limits on the tensor mode; unlike other modes, the tensor shear can grow from a near-isotropic early Universe. The limit on this type of shear is $(σ_{T,\rm reg}/H)_0 < 2.4 \times 10^{-7}$ (95% CL).

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Bursty star formation feedback and cooling outflows

We study how outflows of gas launched from a central galaxy undergoing repeated starbursts propagate through the circumgalactic medium (CGM), using the simulation code RAMSES. We assume that the outflow from the disk can be modelled as a rapidly moving bubble of hot gas at $\mathrm{\sim1\;kpc}$ above disk, then ask what happens as it moves out further into the halo around the galaxy on $\mathrm{\sim 100\;kpc}$ scales. To do this we run 60 two-dimensional simulations scanning over parameters of the outflow. Each of these is repeated with and without radiative cooling, assuming a primordial gas composition to give a lower bound on the importance of cooling. In a large fraction of radiative-cooling cases we are able to form rapidly outflowing cool gas from in situ cooling of the flow. We show that the amount of cool gas formed depends strongly on the 'burstiness' of energy injection; sharper, stronger bursts typically lead to a larger fraction of cool gas forming in the outflow. The abundance ratio of ions in the CGM may therefore change in response to the detailed historical pattern of star formation. For instance, outflows generated by star formation with short, intense bursts contain up to 60 per cent of their gas mass at temperatures $<5 \times 10^4\,\mathrm{K}$; for near-continuous star formation the figure is $\lesssim$ 5 per cent. Further study of cosmological simulations, and of idealised simulations with e.g., metal-cooling, magnetic fields and/or thermal conduction, will help to understand the precise signature of bursty outflows on observed ion abundances.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Designing and testing inflationary models with Bayesian networks

Even simple inflationary scenarios have many free parameters. Beyond the variables appearing in the inflationary action, these include dynamical initial conditions, the number of fields, and couplings to other sectors. These quantities are often ignored but cosmological observables can depend on the unknown parameters. We use Bayesian networks to account for a large set of inflationary parameters, deriving generative models for the primordial spectra that are conditioned on a hierarchical set of prior probabilities describing the initial conditions, reheating physics, and other free parameters. We use $N_f$--quadratic inflation as an illustrative example, finding that the number of $e$-folds $N_*$ between horizon exit for the pivot scale and the end of inflation is typically the most important parameter, even when the number of fields, their masses and initial conditions are unknown, along with possible conditional dependencies between these parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Gaussianisation for fast and accurate inference from cosmological data

We present a method to transform multivariate unimodal non-Gaussian posterior probability densities into approximately Gaussian ones via non-linear mappings, such as Box--Cox transformations and generalisations thereof. This permits an analytical reconstruction of the posterior from a point sample, like a Markov chain, and simplifies the subsequent joint analysis with other experiments. This way, a multivariate posterior density can be reported efficiently, by compressing the information contained in MCMC samples. Further, the model evidence integral (i.e. the marginal likelihood) can be computed analytically. This method is analogous to the search for normal parameters in the cosmic microwave background, but is more general. The search for the optimally Gaussianising transformation is performed computationally through a maximum-likelihood formalism; its quality can be judged by how well the credible regions of the posterior are reproduced. We demonstrate that our method outperforms kernel density estimates in this objective. Further, we select marginal posterior samples from Planck data with several distinct strongly non-Gaussian features, and verify the reproduction of the marginal contours. To demonstrate evidence computation, we Gaussianise the joint distribution of data from weak lensing and baryon acoustic oscillations (BAO), for different cosmological models, and find a preference for flat $Λ$CDM. Comparing to values computed with the Savage-Dickey density ratio, and Population Monte Carlo, we find good agreement of our method within the spread of the other two.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Hierarchical Bayesian inference of galaxy redshift distributions from photometric surveys

Accurately characterizing the redshift distributions of galaxies is essential for analysing deep photometric surveys and testing cosmological models. We present a technique to simultaneously infer redshift distributions and individual redshifts from photometric galaxy catalogues. Our model constructs a piecewise constant representation (effectively a histogram) of the distribution of galaxy types and redshifts, the parameters of which are efficiently inferred from noisy photometric flux measurements. This approach can be seen as a generalization of template-fitting photometric redshift methods and relies on a library of spectral templates to relate the photometric fluxes of individual galaxies to their redshifts. We illustrate this technique on simulated galaxy survey data, and demonstrate that it delivers correct posterior distributions on the underlying type and redshift distributions, as well as on the individual types and redshifts of galaxies. We show that even with uninformative priors, large photometric errors and parameter degeneracies, the redshift and type distributions can be recovered robustly thanks to the hierarchical nature of the model, which is not possible with common photometric redshift estimation techniques. As a result, redshift uncertainties can be fully propagated in cosmological analyses for the first time, fulfilling an essential requirement for the current and future generations of surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

How isotropic is the Universe?

A fundamental assumption in the standard model of cosmology is that the Universe is isotropic on large scales. Breaking this assumption leads to a set of solutions to Einstein's field equations, known as Bianchi cosmologies, only a subset of which have ever been tested against data. For the first time, we consider all degrees of freedom in these solutions to conduct a general test of isotropy using cosmic microwave background temperature and polarization data from Planck. For the vector mode (associated with vorticity), we obtain a limit on the anisotropic expansion of $(σ_V/H)_0 < 4.7 \times 10^{-11}$ (95% CI), which is an order of magnitude tighter than previous Planck results that used CMB temperature only. We also place upper limits on other modes of anisotropic expansion, with the weakest limit arising from the regular tensor mode, $(σ_{T,\rm reg}/H)_0<1.0 \times 10^{-6}$ (95% CI). Including all degrees of freedom simultaneously for the first time, anisotropic expansion of the Universe is strongly disfavoured, with odds of 121,000:1 against.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

How to quench a galaxy

We show how the interplay between active galactic nuclei (AGN) and merger history determines whether a galaxy quenches star formation at high redshift. We first simulate, in a full cosmological context, a galaxy of total dynamical mass $10^{12}\,M_{\odot}$ at $z=2$. Then we systematically alter the accretion history of the galaxy by minimally changing the linear overdensity in the initial conditions. This "genetic modification" approach allows the generation of three sets of $Λ$CDM initial conditions leading to maximum merger ratios of 1:10, 1:5 and 2:3 respectively. The changes leave the final halo mass, large scale structure and local environment unchanged, providing a controlled numerical experiment. Interaction between the AGN physics and mergers in the three cases lead respectively to a star-forming, temporarily-quenched and permanently-quenched galaxy. However the differences do not primarily lie in the black hole accretion rates, but in the kinetic effects of the merger: the galaxy is resilient against AGN feedback unless its gaseous disk is first disrupted. Typical accretion rates are comparable in the three cases, falling below $0.1\,M_{\odot}$ yr$^{-1}$, equivalent to around $2\%$ of the Eddington rate or $10^{-3}$ times the pre-quenching star formation rate, in agreement with observations. This low level of black hole accretion can be sustained even when there is insufficient dense cold gas for star formation. Conversely, supernova feedback is too distributed to generate outflows in high-mass systems, and cannot maintain quenching over periods longer than the halo gas cooling time.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Photometric Supernova Classification With Machine Learning

Automated photometric supernova classification has become an active area of research in recent years in light of current and upcoming imaging surveys such as the Dark Energy Survey (DES) and the Large Synoptic Survey Telescope, given that spectroscopic confirmation of type for all supernovae discovered will be impossible. Here, we develop a multi-faceted classification pipeline, combining existing and new approaches. Our pipeline consists of two stages: extracting descriptive features from the light curves and classification using a machine learning algorithm. Our feature extraction methods vary from model-dependent techniques, namely SALT2 fits, to more independent techniques fitting parametric models to curves, to a completely model-independent wavelet approach. We cover a range of representative machine learning algorithms, including naive Bayes, k-nearest neighbors, support vector machines, artificial neural networks and boosted decision trees (BDTs). We test the pipeline on simulated multi-band DES light curves from the Supernova Photometric Classification Challenge. Using the commonly used area under the curve (AUC) of the Receiver Operating Characteristic as a metric, we find that the SALT2 fits and the wavelet approach, with the BDTs algorithm, each achieves an AUC of 0.98, where 1 represents perfect classification. We find that a representative training set is essential for good classification, whatever the feature set or algorithm, with implications for spectroscopic follow-up. Importantly, we find that by using either the SALT2 or the wavelet feature sets with a BDT algorithm, accurate classification is possible purely from light curve data, without the need for any redshift information.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Robust forecasts on fundamental physics from the foreground-obscured, gravitationally-lensed CMB polarization

[Abridged] Recent results from the BICEP, Keck Array and Planck Collaborations demonstrate that Galactic foregrounds are an unavoidable obstacle in the search for evidence of inflationary gravitational waves in the cosmic microwave background (CMB) polarization. Beyond the foregrounds, the effect of lensing by intervening large-scale structure further obscures all but the strongest inflationary signals permitted by current data. With a plethora of ongoing and upcoming experiments aiming to measure these signatures, careful and self-consistent consideration of experiments' foreground- and lensing-removal capabilities is critical in obtaining credible forecasts of their performance. We investigate the capabilities of instruments such as Advanced ACTPol, BICEP3 and Keck Array, CLASS, EBEX10K, PIPER, Simons Array, SPT-3G and SPIDER, and projects as COrE+, LiteBIRD-ext, PIXIE and Stage IV, to clean contamination due to polarized synchrotron and dust from raw multi-frequency data, and remove lensing from the resulting co-added CMB maps (either using iterative CMB-only techniques or through cross-correlation with external data). Incorporating these effects, we present forecasts for the constraining power of these experiments in terms of inflationary physics, the neutrino sector, and dark energy parameters. Made publicly available through an online interface, this tool enables the next generation of CMB experiments to foreground-proof their designs, optimize their frequency coverage to maximize scientific output, and determine where cross-experimental collaboration would be most beneficial. We find that analyzing data from ground, balloon and space instruments in complementary combinations can significantly improve component separation performance, delensing, and cosmological constraints over individual datasets.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

SILC: a new Planck Internal Linear Combination CMB temperature map using directional wavelets

We present new clean maps of the CMB temperature anisotropies (as measured by Planck) constructed with a novel internal linear combination (ILC) algorithm using directional, scale-discretised wavelets --- Scale-discretised, directional wavelet ILC or SILC. Directional wavelets, when convolved with signals on the sphere, can separate the anisotropic filamentary structures which are characteristic of both the CMB and foregrounds. Extending previous component separation methods, which use the frequency, spatial and harmonic signatures of foregrounds to separate them from the cosmological background signal, SILC can additionally use morphological information in the foregrounds and CMB to better localise the cleaning algorithm. We test the method on Planck data and simulations, demonstrating consistency with existing component separation algorithms, and discuss how to optimise the use of morphological information by varying the number of directional wavelets as a function of spatial scale. We find that combining the use of directional and axisymmetric wavelets depending on scale could yield higher quality CMB temperature maps. Our results set the stage for the application of SILC to polarisation anisotropies through an extension to spin wavelets.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Spin-SILC: CMB polarisation component separation with spin wavelets

We present Spin-SILC, a new foreground component separation method that accurately extracts the cosmic microwave background (CMB) polarisation $E$ and $B$ modes from raw multifrequency Stokes $Q$ and $U$ measurements of the microwave sky. Spin-SILC is an internal linear combination method that uses spin wavelets to analyse the spin-2 polarisation signal $P = Q + iU$. The wavelets are additionally directional (non-axisymmetric). This allows different morphologies of signals to be separated and therefore the cleaning algorithm is localised using an additional domain of information. The advantage of spin wavelets over standard scalar wavelets is to simultaneously and self-consistently probe scales and directions in the polarisation signal $P = Q + iU$ and in the underlying $E$ and $B$ modes, therefore providing the ability to perform component separation and $E$-$B$ decomposition concurrently for the first time. We test Spin-SILC on full-mission Planck simulations and data and show the capacity to correctly recover the underlying cosmological $E$ and $B$ modes. We also demonstrate a strong consistency of our CMB maps with those derived from existing component separation methods. Spin-SILC can be combined with the pseudo- and pure $E$-$B$ spin wavelet estimators presented in a companion paper to reliably extract the cosmological signal in the presence of complicated sky cuts and noise. Therefore, it will provide a computationally-efficient method to accurately extract the CMB $E$ and $B$ modes for future polarisation experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The DESI Experiment Part I: Science,Targeting, and Survey Design

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$α$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The DESI Experiment Part II: Instrument Design

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= λ/Δλ$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Unbiased pseudo-Cl power spectrum estimation with mode projection

With the steadily improving sensitivity afforded by current and future galaxy surveys, a robust extraction of two-point correlation function measurements may become increasingly hampered by the presence of astrophysical foregrounds or observational systematics. The concept of mode projection has been introduced as a means to remove contaminants for which it is possible to construct a spatial map reflecting the expected signal contribution. Owing to its computational efficiency compared to minimum-variance methods, the sub-optimal pseudo-Cl (PCL) power spectrum estimator is a popular tool for the analysis of high-resolution data sets. Here, we integrate mode projection into the framework of PCL power spectrum estimation. In contrast to results obtained with optimal estimators, we show that the uncorrected projection of template maps leads to biased power spectra. Based on analytical calculations, we find exact closed-form expressions for the expectation value of the bias and demonstrate that they can be recast in a form that allows a numerically efficient evaluation, preserving the favorable O(l_max^3) time complexity of PCL estimator algorithms. Using simulated data sets, we assess the scaling of the bias with various analysis parameters and demonstrate that it can be reliably removed. We conclude that in combination with mode projection, PCL estimators allow for a fast and robust computation of power spectra in the presence of systematic effects -- properties in high demand for the analysis of ongoing and future large scale structure surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

3D weak lensing with spin wavelets on the ball

We construct the spin flaglet transform, a wavelet transform to analyze spin signals in three dimensions. Spin flaglets can probe signal content localized simultaneously in space and frequency and, moreover, are separable so that their angular and radial properties can be controlled independently. They are particularly suited to analyzing of cosmological observations such as the weak gravitational lensing of galaxies. Such observations have a unique 3D geometrical setting since they are natively made on the sky, have spin angular symmetries, and are extended in the radial direction by additional distance or redshift information. Flaglets are constructed in the harmonic space defined by the Fourier-Laguerre transform, previously defined for scalar functions and extended here to signals with spin symmetries. Thanks to various sampling theorems, both the Fourier-Laguerre and flaglet transforms are theoretically exact when applied to bandlimited signals. In other words, in numerical computations the only loss of information is due to the finite representation of floating point numbers. We develop a 3D framework relating the weak lensing power spectrum to covariances of flaglet coefficients. We suggest that the resulting novel flaglet weak lensing estimator offers a powerful alternative to common 2D and 3D approaches to accurately capture cosmological information. While standard weak lensing analyses focus on either real or harmonic space representations (i.e., correlation functions or Fourier-Bessel power spectra, respectively), a wavelet approach inherits the advantages of both techniques, where both complicated sky coverage and uncertainties associated with the physical modeling of small scales can be handled effectively. Our codes to compute the Fourier-Laguerre and flaglet transforms are made publicly available.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Constraints on primordial non-Gaussianity from 800,000 photometric quasars

We derive robust constraints on primordial non-Gaussianity (PNG) using the clustering of 800,000 photometric quasars from the Sloan Digital Sky Survey in the redshift range $0.5<z<3.5$. These measurements rely on the novel technique of {\it extended mode projection} to control the impact of spatially-varying systematics in a robust fashion, making use of blind analysis techniques. This allows the accurate measurement of quasar halo bias at the largest scales, while discarding as little as possible of the data. The standard local-type PNG parameters $f_\mathrm{NL}$ and $g_\mathrm{NL}$ both imprint a $k^{-2}$ scale-dependent effect in the bias. Constraining these individually, we obtain $-49<f_\mathrm{NL}<31$ and $-2.7\times10^5<g_\mathrm{NL}<1.9\times10^5$, while their joint constraints lead to $-105<f_\mathrm{NL}<72$ and $-4.0\times10^5<g_\mathrm{NL}<4.9\times10^5$ (all at 95% CL) . Introducing a running parameter $n_{f_\mathrm{NL}}$ to constrain $b(k) \propto k^{-2+n_{f_\mathrm{NL}}}$ and a generalised PNG amplitude $\tilde{f}_\mathrm{NL}$, we obtain $-45.5 \exp({3.7\, n_{f_\mathrm{NL}}}) < \tilde{f}_\mathrm{NL} < 34.4 \exp({3.3\, n_{f_\mathrm{NL}}})$ at 95% CL. These results incorporate uncertainties in the cosmological parameters, redshift distributions, shot noise, and the bias prescription used to relate the quasar clustering to the underlying dark matter. These are the strongest constraints obtained to date on PNG using a single population of large-scale structure tracers, and are already at the level of pre-{\it Planck} constraints from the cosmic microwave background. A conservative forecast for a {\it Large Synoptic Survey Telescope}-like survey incorporating mode projection yields $σ(f_\mathrm{NL}) \sim 5$ -- competitive with the {\it Planck} result -- highlighting the power of upcoming large scale structure surveys to probe the initial conditions of the universe.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Exploiting the full potential of photometric quasar surveys: Optimal power spectra through blind mitigation of systematics

We present optimal measurements of the angular power spectrum of the XDQSOz catalogue of photometric quasars from the Sloan Digital Sky Survey. These measurements rely on a quadratic maximum likelihood estimator that simultaneously measures the auto- and cross-power spectra of four redshift samples, and provides minimum-variance, unbiased estimates even at the largest angular scales. Since photometric quasars are known to be strongly affected by systematics such as spatially-varying depth and stellar contamination, we introduce a new framework of extended mode projection to robustly mitigate the impact of systematics on the power spectrum measurements. This technique involves constructing template maps of potential systematics, decorrelating them on the sky, and projecting out modes which are significantly correlated with the data. Our method is able to simultaneously process several thousands of nonlinearly-correlated systematics, and mode projection is performed in a blind fashion. Using our final power spectrum measurements, we find a good agreement with theoretical predictions, and no evidence for further contamination by systematics. Extended mode projection not only obviates the need for aggressive sky and quality cuts, but also provides control over the level of systematics in the measurements, enabling the search for small signals of new physics while avoiding confirmation bias.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Forecasting constraints from the cosmic microwave background on eternal inflation

We forecast the ability of cosmic microwave background (CMB) temperature and polarization datasets to constrain theories of eternal inflation using cosmic bubble collisions. Using the Fisher matrix formalism, we determine both the overall detectability of bubble collisions and the constraints achievable on the fundamental parameters describing the underlying theory. The CMB signatures considered are based on state-of-the-art numerical relativistic simulations of the bubble collision spacetime, evolved using the full temperature and polarization transfer functions. Comparing a theoretical cosmic-variance-limited experiment to the WMAP and Planck satellites, we find that there is no improvement to be gained from future temperature data, that adding polarization improves detectability by approximately 30%, and that cosmic-variance-limited polarization data offer only marginal improvements over Planck. The fundamental parameter constraints achievable depend on the precise values of the tensor-to-scalar ratio and energy density in (negative) spatial curvature. For a tensor-to-scalar ratio of $0.1$ and spatial curvature at the level of $10^{-4}$, using cosmic-variance-limited data it is possible to measure the width of the potential barrier separating the inflating false vacuum from the true vacuum down to $M_{\rm Pl}/500$, and the initial proper distance between colliding bubbles to a factor $π/2$ of the false vacuum horizon size (at three sigma). We conclude that very near-future data will have the final word on bubble collisions in the CMB.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Genetically modified halos: towards controlled experiments in $Λ$CDM galaxy formation

We propose a method to generate `genetically-modified' (GM) initial conditions for high-resolution simulations of galaxy formation in a cosmological context. Building on the Hoffman-Ribak algorithm, we start from a reference simulation with fully random initial conditions, then make controlled changes to specific properties of a single halo (such as its mass and merger history). The algorithm demonstrably makes minimal changes to other properties of the halo and its environment, allowing us to isolate the impact of a given modification. As a significant improvement over previous work, we are able to calculate the abundance of the resulting objects relative to the reference simulation. Our approach can be applied to a wide range of cosmic structures and epochs; here we study two problems as a proof-of-concept. First, we investigate the change in density profile and concentration as the collapse time of three individual halos are varied at fixed final mass, showing good agreement with previous statistical studies using large simulation suites. Second, we modify the $z=0$ mass of halos to show that our theoretical abundance calculations correctly recover the halo mass function. The results demonstrate that the technique is robust, opening the way to controlled experiments in galaxy formation using hydrodynamic zoom simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Inverted initial conditions: exploring the growth of cosmic structure and voids

We introduce and explore "paired" cosmological simulations. A pair consists of an A and B simulation with initial conditions related by the inversion $δ_A(x, t_{initial})=-δ_B(x,t_{initial})$ (underdensities substituted for overdensities and vice versa). We argue that the technique is valuable for improving our understanding of cosmic structure formation. The A and B fields are by definition equally likely draws from ΛCDM initial conditions, and in the linear regime evolve identically up to the overall sign. As non-linear evolution takes hold, a region that collapses to form a halo in simulation A will tend to expand to create a void in simulation B. Applications include (i) contrasting the growth of A-halos and B-voids to test excursion-set theories of structure formation; (ii) cross-correlating the density field of the A and B universes as a novel test for perturbation theory; and (iii) canceling error terms by averaging power spectra between the two boxes. Generalizations of the method to more elaborate field transformations are suggested.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Unbiased methods for removing systematics from galaxy clustering measurements

Measuring the angular clustering of galaxies as a function of redshift is a powerful method for extracting information from the three-dimensional galaxy distribution. The precision of such measurements will dramatically increase with ongoing and future wide-field galaxy surveys. However, these are also increasingly sensitive to observational and astrophysical contaminants. Here, we study the statistical properties of three methods proposed for controlling such systematics -- template subtraction, basic mode projection, and extended mode projection -- all of which make use of externally supplied template maps, designed to characterise and capture the spatial variations of potential systematic effects. Based on a detailed mathematical analysis, and in agreement with simulations, we find that the template subtraction method in its original formulation returns biased estimates of the galaxy angular clustering. We derive closed-form expressions that should be used to correct results for this shortcoming. Turning to the basic mode projection algorithm, we prove it to be free of any bias, whereas we conclude that results computed with extended mode projection are biased. Within a simplified setup, we derive analytical expressions for the bias and discuss the options for correcting it in more realistic configurations. Common to all three methods is an increased estimator variance induced by the cleaning process, albeit at different levels. These results enable unbiased high-precision clustering measurements in the presence of spatially-varying systematics, an essential step towards realising the full potential of current and planned galaxy surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Considerations in the Interpretation of Cosmological Anomalies

Anomalies drive scientific discovery -- they are associated with the cutting edge of the research frontier, and thus typically exploit data in the low signal-to-noise regime. In astronomy, the prevalence of systematics --- both "known unknowns" and "unknown unknowns" --- combined with increasingly large datasets, the widespread use of ad hoc estimators for anomaly detection, and the "look-elsewhere" effect, can lead to spurious false detections. In this informal note, I argue that anomaly detection leading to discoveries of new physics requires a combination of physical understanding, careful experimental design to avoid confirmation bias, and self-consistent statistical methods. These points are illustrated with several concrete examples from cosmology.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Gravitational wave consistency relations for multifield inflation

We study the tensor spectral index $n_t$ and the tensor-to-scalar ratio $r$ in the simplest multifield extension to single-field, slow-roll inflation models. We show that multifield models with potentials $V \sim \sum_i λ_i |ϕ_i|^p$ have different predictions for $n_t/r$ than single-field models, even when all the couplings are equal $λ_i=λ_j$, due to the probabilistic nature of the fields' initial values. We analyze well-motivated prior probabilities for the $λ_i$ and initial conditions to make detailed predictions for the marginalized probability distribution of $n_t/r$. With $\mathcal O(100)$ fields and $p>3/4$, we find that $n_t/r$ differs from the single-field result of $n_t/r=-1/8$ at the 5$σ$ level. This gives a novel and testable prediction for the simplest multifield inflation models.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

MultiModeCode: An efficient numerical solver for multifield inflation

We present MultiModeCode, a Fortran 95/2000 package for the numerical exploration of multifield inflation models. This program facilitates efficient Monte Carlo sampling of prior probabilities for inflationary model parameters and initial conditions and is the first publicly available code that can efficiently generate large sample-sets for inflation models with $\mathcal O(100)$ fields. The code numerically solves the equations of motion for the background and first-order perturbations of multi-field inflation models with canonical kinetic terms and arbitrary potentials, providing the adiabatic, isocurvature, and tensor power spectra at the end of inflation. For models with sum-separable potentials MultiModeCode also computes the slow-roll prediction via the $δN$ formalism for easy model exploration and validation. We pay particular attention to the isocurvature perturbations as the system approaches the adiabatic limit, showing how to avoid numerical instabilities that affect some other approaches to this problem. We demonstrate the use of MultiModeCode by exploring a few toy models. Finally, we give a concise review of multifield perturbation theory and a user's manual for the program.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

No new cosmological concordance with massive sterile neutrinos

It has been claimed recently that massive sterile neutrinos could bring about a new concordance between observations of the cosmic microwave background (CMB), the large-scale structure (LSS) of the Universe, and local measurements of the Hubble constant, $H_0$. We demonstrate that this apparent concordance results from combining datasets which are in significant tension, even within this extended model, possibly indicating remaining systematic biases in the measurements. We further show that this tension remains when the cosmological model is further extended to include significant tensor modes, as suggested by the recent BICEP2 results. Using the Bayesian evidence, we show that the minimal $Λ$CDM model is strongly favoured over its neutrino extensions by various combinations of datasets. Robust data combinations yield stringent limits of $\sum m_ν\lesssim0.3$ eV and $m_{ν,{\rm sterile}}^{\rm eff} \lesssim 0.3$ eV at $95\%$ CL for the sum of active and sterile neutrinos, respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On spin scale-discretised wavelets on the sphere for the analysis of CMB polarisation

A new spin wavelet transform on the sphere is proposed to analyse the polarisation of the cosmic microwave background (CMB), a spin $\pm 2$ signal observed on the celestial sphere. The scalar directional scale-discretised wavelet transform on the sphere is extended to analyse signals of arbitrary spin. The resulting spin scale-discretised wavelet transform probes the directional intensity of spin signals. A procedure is presented using this new spin wavelet transform to recover E- and B-mode signals from partial-sky observations of CMB polarisation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Simple predictions from multifield inflationary models

We explore whether multifield inflationary models make unambiguous predictions for fundamental cosmological observables. Focusing on $N$-quadratic inflation, we numerically evaluate the full perturbation equations for models with 2, 3, and $\mathcal{O}(100)$ fields, using several distinct methods for specifying the initial values of the background fields. All scenarios are highly predictive, with the probability distribution functions of the cosmological observables becoming more sharply peaked as $N$ increases. For $N=100$ fields, 95% of our Monte Carlo samples fall in the ranges $n_s \in (0.9455,0.9534)$; $α\in (-9.741,-7.047)\times 10^{-4}$; $r\in(0.1445,0.1449)$; and $r_\mathrm{iso} \in (0.02137,3.510)\times 10^{-3}$ for the spectral index, running, tensor-to-scalar ratio, and isocurvature-to-adiabatic ratio, respectively. The expected amplitude of isocurvature perturbations grows with $N$, raising the possibility that many-field models may be sensitive to post-inflationary physics and suggesting new avenues for testing these scenarios.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Simulating the universe(s) II: phenomenology of cosmic bubble collisions in full General Relativity

Observing the relics of collisions between bubble universes would provide direct evidence for the existence of an eternally inflating Multiverse; the non-observation of such events can also provide important constraints on inflationary physics. Realizing these prospects requires quantitative predictions for observables from the properties of the possible scalar field Lagrangians underlying eternal inflation. Building on previous work, we establish this connection in detail. We perform a fully relativistic numerical study of the phenomenology of bubble collisions in models with a single scalar field, computing the comoving curvature perturbation produced in a wide variety of models. We also construct a set of analytic predictions, allowing us to identify the phenomenologically relevant properties of the scalar field Lagrangian. The agreement between the analytic predictions and numerics in the relevant regions is excellent, and allows us to generalize our results beyond the models we adopt for the numerical studies. Specifically, the signature is completely determined by the spatial profile of the colliding bubble just before the collision, and the de Sitter invariant distance between the bubble centers. The analytic and numerical results support a power-law fit with an index $1< κ\lesssim 2$. For collisions between identical bubbles, we establish a lower-bound on the observed amplitude of collisions that is set by the present energy density in curvature.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Simulating the universe(s): from cosmic bubble collisions to cosmological observables with numerical relativity

The theory of eternal inflation in an inflaton potential with multiple vacua predicts that our universe is one of many bubble universes nucleating and growing inside an ever-expanding false vacuum. The collision of our bubble with another could provide an important observational signature to test this scenario. We develop and implement an algorithm for accurately computing the cosmological observables arising from bubble collisions directly from the Lagrangian of a single scalar field. We first simulate the collision spacetime by solving Einstein's equations, starting from nucleation and ending at reheating. Taking advantage of the collision's hyperbolic symmetry, simulations are performed with a 1+1-dimensional fully relativistic code that uses adaptive mesh refinement. We then calculate the comoving curvature perturbation in an open Friedmann-Robertson-Walker universe, which is used to determine the temperature anisotropies of the cosmic microwave background radiation. For a fiducial Lagrangian, the anisotropies are well described by a power law in the cosine of the angular distance from the center of the collision signature. For a given form of the Lagrangian, the resulting observational predictions are inherently statistical due to stochastic elements of the bubble nucleation process. Further uncertainties arise due to our imperfect knowledge about inflationary and pre-recombination physics. We characterize observational predictions by computing the probability distributions over four phenomenological parameters which capture these intrinsic and model uncertainties. This represents the first fully-relativistic set of predictions from an ensemble of scalar field models giving rise to eternal inflation, yielding significant differences from previous non-relativistic approximations. Thus, our results provide a basis for a rigorous confrontation of these theories with cosmological data.

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
preprint2013arXiv

(Lack of) Cosmological evidence for dark radiation after Planck

We use Bayesian model comparison to determine whether extensions to Standard-Model neutrino physics -- primarily additional effective numbers of neutrinos and/or massive neutrinos -- are merited by the latest cosmological data. Given the significant advances in cosmic microwave background (CMB) observations represented by the Planck data, we examine whether Planck temperature and CMB lensing data, in combination with lower redshift data, have strengthened (or weakened) the previous findings. We conclude that the state-of-the-art cosmological data do not show evidence for deviations from the standard cosmological model (which has three massless neutrino families). This does not mean that the model is necessarily correct -- in fact we know it is incomplete as neutrinos are not massless -- but it does imply that deviations from the standard model (e.g., non-zero neutrino mass) are too small compared to the current experimental uncertainties to be inferred from cosmological data alone.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Estimating the large-scale angular power spectrum in the presence of systematics: a case study of Sloan Digital Sky Survey quasars

The angular power spectrum is a powerful statistic for analysing cosmological signals imprinted in the clustering of matter. However, current galaxy and quasar surveys cover limited portions of the sky, and are contaminated by systematics that can mimic cosmological signatures and jeopardise the interpretation of the measured power spectra. We provide a framework for obtaining unbiased estimates of the angular power spectra of large-scale structure surveys at the largest scales using quadratic estimators. The method is tested by analysing the 600 CMASS mock catalogues constructed by Manera et al. (2013) for the Baryon Oscillation Spectroscopic Survey (BOSS). We then consider the Richards et al. (2009) catalogue of photometric quasars from the Sixth Data Release (DR6) of the Sloan Digital Sky Survey (SDSS), which is known to include significant stellar contamination and systematic uncertainties. Focusing on the sample of ultraviolet-excess (UVX) sources, we show that the excess clustering power present on the largest-scales can be largely mitigated by making use of improved sky masks and projecting out the modes corresponding to the principal systematics. In particular, we find that the sample of objects with photometric redshift $1.3 < z_p < 2.2$ exhibits no evidence of contamination when using our most conservative mask and mode projection. This indicates that any residual systematics are well within the statistical uncertainties. We conclude that, using our approach, this sample can be used for cosmological studies.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Hierarchical Bayesian Detection Algorithm for Early-Universe Relics in the Cosmic Microwave Background

A number of theoretically well-motivated additions to the standard cosmological model predict weak signatures in the form of spatially localized sources embedded in the cosmic microwave background (CMB) fluctuations. We present a hierarchical Bayesian statistical formalism and a complete data analysis pipeline for testing such scenarios. We derive an accurate approximation to the full posterior probability distribution over the parameters defining any theory that predicts sources embedded in the CMB, and perform an extensive set of tests in order to establish its validity. The approximation is implemented using a modular algorithm, designed to avoid a posteriori selection effects, which combines a candidate-detection stage with a full Bayesian model-selection and parameter-estimation analysis. We apply this pipeline to theories that predict cosmic textures and bubble collisions, extending previous analyses by using: (1) adaptive-resolution techniques, allowing us to probe features of arbitrary size, and (2) optimal filters, which provide the best possible sensitivity for detecting candidate signatures. We conclude that the WMAP 7-year data do not favor the addition of either cosmic textures or bubble collisions to the standard cosmological model, and place robust constraints on the predicted number of such sources. The expected numbers of bubble collisions and cosmic textures on the CMB sky within our detection thresholds are constrained to be fewer than 4.0 and 5.2 at 95% confidence, respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Is there evidence for additional neutrino species from cosmology?

It has been suggested that recent cosmological and flavor-oscillation data favor the existence of additional neutrino species beyond the three predicted by the Standard Model of particle physics. We apply Bayesian model selection to determine whether there is indeed any evidence from current cosmological datasets for the standard cosmological model to be extended to include additional neutrino flavors. The datasets employed include cosmic microwave background temperature, polarization and lensing power spectra, and measurements of the baryon acoustic oscillation scale and the Hubble constant. We also consider other extensions to the standard neutrino model, such as massive neutrinos, and possible degeneracies with other cosmological parameters. The Bayesian evidence indicates that current cosmological data do not require any non-standard neutrino properties.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Bayesian Analysis of Inflation III: Slow Roll Reconstruction Using Model Selection

We implement Slow Roll Reconstruction -- an optimal solution to the inverse problem for inflationary cosmology -- within ModeCode, a publicly available solver for the inflationary dynamics. We obtain up-to-date constraints on the reconstructed inflationary potential, derived from the WMAP 7-year dataset and South Pole Telescope observations, combined with large scale structure data derived from SDSS Data Release 7. Using ModeCode in conjunction with the MultiNest sampler, we compute Bayesian evidence for the reconstructed potential at each order in the truncated slow roll hierarchy. We find that the data are well-described by the first two slow roll parameters, εand η, and that there is no need to include a nontrivial ξparameter.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Determining the outcome of cosmic bubble collisions in full General Relativity

Cosmic bubble collisions provide an important possible observational window on the dynamics of eternal inflation. In eternal inflation, our observable universe is contained in one of many bubbles formed from an inflating metastable vacuum. The collision between bubbles can leave a detectable imprint on the cosmic microwave background radiation. Although phenomenological models of the observational signature have been proposed, to make the theory fully predictive one must determine the bubble collision spacetime, and thus the cosmological observables, from a scalar field theory giving rise to eternal inflation. Because of the intrinsically non-linear nature of the bubbles and their collision, this requires a numerical treatment incorporating General Relativity. In this paper, we present results from numerical simulations of bubble collisions in full General Relativity. These simulations allow us to accurately determine the outcome of bubble collisions, and examine their effect on the cosmology inside a bubble universe. We confirm the validity of a number of approximations used in previous analytic work, and identify qualitatively new features of bubble collision spacetimes. Both vacuum bubbles and bubbles containing a realistic inflationary cosmology are studied. We identify the constraints on the scalar field potential that must be satisfied in order to obtain collisions that are consistent with our observed cosmology, yet leave detectable signatures.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Fast Computation of Bispectrum Features with Generalized Slow Roll

We develop a fast technique based on the generalized slow roll (GSR) approach for computing the curvature bispectrum of inflationary models with features. We show that all triangle configurations can be expressed in terms of three simple integrals over the inflationary background with typical accuracy of better than ~20%. With a first order GSR approach the typical accuracy can be improved to better than the 5% level. We illustrate this technique with the step potential model that has been invoked to explain the WMAP temperature power spectrum glitches at ell ~ 20-40 and show that the maximum likelihood model falls short of observability by more than a factor of 100 in amplitude. We also explicitly demonstrate that the bispectrum consistency relation with the local slope of the power spectrum is satisfied for these models. In the GSR approach, the bispectrum arises from integrals of nearly the same function of the background slow-roll parameters as the power spectrum but with a stronger weight to the epoch before horizon crossing. Hence this technique enables reverse engineering of models with large bispectrum but small power spectrum features.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Avoiding bias in reconstructing the largest observable scales from partial-sky data

Obscuration due to Galactic emission complicates the extraction of information from cosmological surveys, and requires some combination of the (typically imperfect) modeling and subtraction of foregrounds, or the removal of part of the sky. This particularly affects the extraction of information from the largest observable scales. Maximum-likelihood estimators for reconstructing the full-sky spherical harmonic coefficients from partial-sky maps have recently been shown to be susceptible to contamination from within the sky cut, arising due to the necessity to band-limit the data by smoothing prior to reconstruction. Using the WMAP 7-year data, we investigate modified implementations of such estimators which are robust to the leakage of contaminants from within masked regions. We provide a measure, based on the expected amplitude of residual foregrounds, for selecting the most appropriate estimator for the task at hand. We explain why the related quadratic maximum-likelihood estimator of the angular power spectrum does not suffer from smoothing-induced bias.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Bayesian Analysis of Inflation: Parameter Estimation for Single Field Models

Future astrophysical datasets promise to strengthen constraints on models of inflation, and extracting these constraints requires methods and tools commensurate with the quality of the data. In this paper we describe ModeCode, a new, publicly available code that computes the primordial scalar and tensor power spectra for single field inflationary models. ModeCode solves the inflationary mode equations numerically, avoiding the slow roll approximation. It is interfaced with CAMB and CosmoMC to compute cosmic microwave background angular power spectra and perform likelihood analysis and parameter estimation. ModeCode is easily extendable to additional models of inflation, and future updates will include Bayesian model comparison. Errors from ModeCode contribute negligibly to the error budget for analyses of data from Planck or other next generation experiments. We constrain representative single field models (phi^n with n=2/3, 1, 2, and 4, natural inflation, and "hilltop" inflation) using current data, and provide forecasts for Planck. From current data, we obtain weak but nontrivial limits on the post-inflationary physics, which is a significant source of uncertainty in the predictions of inflationary models, while we find that Planck will dramatically improve these constraints. In particular, Planck will link the inflationary dynamics with the post-inflationary growth of the horizon, and thus begin to probe the "primordial dark ages" between TeV and GUT scale energies.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

First Observational Tests of Eternal Inflation

The eternal inflation scenario predicts that our observable universe resides inside a single bubble embedded in a vast inflating multiverse. We present the first observational tests of eternal inflation, performing a search for cosmological signatures of collisions with other bubble universes in cosmic microwave background data from the WMAP satellite. We conclude that the WMAP 7-year data do not warrant augmenting LCDM with bubble collisions, constraining the average number of detectable bubble collisions on the full sky to be less than 1.6 at 68% CL. Data from the Planck satellite can be used to more definitively test the bubble collision hypothesis.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

First Observational Tests of Eternal Inflation: Analysis Methods and WMAP 7-Year Results

In the picture of eternal inflation, our observable universe resides inside a single bubble nucleated from an inflating false vacuum. Many of the theories giving rise to eternal inflation predict that we have causal access to collisions with other bubble universes, providing an opportunity to confront these theories with observation. We present the results from the first observational search for the effects of bubble collisions, using cosmic microwave background data from the WMAP satellite. Our search targets a generic set of properties associated with a bubble collision spacetime, which we describe in detail. We use a modular algorithm that is designed to avoid a posteriori selection effects, automatically picking out the most promising signals, performing a search for causal boundaries, and conducting a full Bayesian parameter estimation and model selection analysis. We outline each component of this algorithm, describing its response to simulated CMB skies with and without bubble collisions. Comparing the results for simulated bubble collisions to the results from an analysis of the WMAP 7-year data, we rule out bubble collisions over a range of parameter space. Our model selection results based on WMAP 7-year data do not warrant augmenting LCDM with bubble collisions. Data from the Planck satellite can be used to more definitively test the bubble collision hypothesis.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Minimally Parametric Power Spectrum Reconstruction from the Lyman-alpha Forest

Current results from the Lyman alpha forest assume that the primordial power spectrum of density perturbations follows a simple power law form, with running. We present the first analysis of Lyman alpha data to study the effect of relaxing this strong assumption on primordial and astrophysical constraints. We perform a large suite of numerical simulations, using them to calibrate a minimally parametric framework for describing the power spectrum. Combined with cross-validation, a statistical technique which prevents over-fitting of the data, this framework allows us to reconstruct the power spectrum shape without strong prior assumptions. We find no evidence for deviation from scale-invariance; our analysis also shows that current Lyman alpha data do not have sufficient statistical power to robustly probe the shape of the power spectrum at these scales. In contrast, the ongoing Baryon Oscillation Sky Survey (BOSS) will be able to do so with high precision. Furthermore, this near-future data will be able to break degeneracies between the power spectrum shape and astrophysical parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

CMB Isotropy Anomalies and the Local Kinetic Sunyaev-Zel'dovich Effect

Several anomalies have been identified which may imply a breakdown of the statistical isotropy of the cosmic microwave background (CMB). In particular, an anomalous alignment of the quadrupole and octopole and a hemispherical power asymmetry have increased in significance as the data have improved. There have been several attempts to explain these observations which explore isotropy breaking mechanisms within the early universe, but little attention has been given to the possibility that these anomalies have their origin within the local universe. We explore such a mechanism by considering the kinetic Sunyaev-Zel'dovich effect due to a gaseous halo associated with the Milky Way. Considering several physical models of an anisotropic free electron optical depth contributed by such a halo, we find that the associated screening maps of the primordial anisotropies have the necessary orientations to affect the anomaly statistics very significantly, but only if the column density of free electrons in the halo is at least an order of magnitude higher than indicated by current observations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The cut-sky cosmic microwave background is not anomalous

The observed angular correlation function of the cosmic microwave background has previously been reported to be anomalous, particularly when measured in regions of the sky uncontaminated by Galactic emission. Recent work by Efstathiou et al. presents a Bayesian comparison of isotropic theories, casting doubt on the significance of the purported anomaly. We extend this analysis to all anisotropic Gaussian theories with vanishing mean (<delta T> = 0), using the much wider class of models to confirm that the anomaly is not likely to point to new physics. On the other hand if there is any new physics to be gleaned, it results from low-l alignments which will be better quantified by a full-sky statistic. We also consider quadratic maximum likelihood power spectrum estimators that are constructed assuming isotropy. The underlying assumptions are therefore false if the ensemble is anisotropic. Nonetheless we demonstrate that, for theories compatible with the observed sky, these estimators (while no longer optimal) remain statistically superior to pseudo-C_l power spectrum estimators.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Cosmological Inflation: Theory and Observations

In this article we review the theory of cosmological inflation with a particular focus on the beautiful connection it provides between the physics of the very small and observations of the very large. We explain how quantum mechanical fluctuations during the inflationary era become macroscopic density fluctuations which leave distinct imprints in the cosmic microwave background (CMB). We describe the physics of anisotropies in the CMB temperature and polarization and discuss how CMB observations can be used to probe the primordial universe.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Photometric constraints on white dwarfs and the identification of extreme objects

It is possible to reliably identify white dwarfs (WDs) without recourse to spectra, instead using photometric and astrometric measurements to distinguish them from Main Sequence stars and quasars. WDs' colours can also be used to infer their intrinsic properties (effective temperature, surface gravity, etc.), but the results obtained must be interpreted with care. The difficulties stem from the existence of a solid angle degeneracy, as revealed by a full exploration of the likelihood, although this can be masked if a simple best-fit approach is used. Conversely, this degeneracy can be broken if a Bayesian approach is adopted, as it is then possible to utilise the prior information on the surface gravities of WDs implied by spectroscopic fitting. The benefits of such an approach are particularly strong when applied to outliers, such as the candidate halo and ultra-cool WDs identified by Vidrih et al. (2007). A reanalysis of these samples confirms their results for the latter sample but suggests that that most of the halo candidates are thick disk WDs in the tails of the photometric noise distribution.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The Shape of the Primordial Power Spectrum: A Last Stand Before Planck

We present a minimally-parametric reconstruction of the primordial power spectrum using the most recent cosmic microwave background and large scale structure data sets. Our goal is to constrain the shape of the power spectrum while simultaneously avoiding strong theoretical priors and over-fitting of the data. We find no evidence for any departure from a power law spectral index. We also find that an exact scale-invariant power spectrum is disfavored by the data, but this conclusion is weaker than the corresponding result assuming a theoretically-motivated power law spectral index prior. The reconstruction shows that better data are crucial to justify the adoption of such a strong theoretical prior observationally. These results can be used to determine the robustness of our present knowledge when compared with forthcoming precision data from Planck.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Cosmological Constraints on Dissipative Models of Inflation

(Abridged) We study dissipative inflation in the regime where the dissipative term takes a specific form, Γ=Γ(ϕ), analyzing two models in the weak and strong dissipative regimes with a SUSY breaking potential. After developing intuition about the predictions from these models through analytic approximations, we compute the predicted cosmological observables through full numerical evolution of the equations of motion, relating the mass scale and scale of dissipation to the characteristic amplitude and shape of the primordial power spectrum. We then use Markov Chain Monte Carlo techniques to constrain a subset of the models with cosmological data from the cosmic microwave background (WMAP three-year data) and large scale structure (SDSS Luminous Red Galaxy power spectrum). We find that the posterior distributions of the dissipative parameters are highly non-Gaussian and their allowed ranges agree well with the expectations obtained using analytic approximations. In the weak regime, only the mass scale is tightly constrained; conversely, in the strong regime, only the dissipative coefficient is tightly constrained. A lower limit is seen on the inflation scale: a sub-Planckian inflaton is disfavoured by the data. In both weak and strong regimes, we reconstruct the limits on the primordial power spectrum and show that these models prefer a {\it red} spectrum, with no significant running of the index. We calculate the reheat temperature and show that the gravitino problem can be overcome with large dissipation, which in turn leads to large levels of non-Gaussianity: if dissipative inflation is to evade the gravitino problem, the predicted level of non-Gaussianity might be seen by the Planck satellite.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

On Minimally-Parametric Primordial Power Spectrum Reconstruction and the Evidence for a Red Tilt

The latest cosmological data seem to indicate a significant deviation from scale invariance of the primordial power spectrum when parameterized either by a power law or by a spectral index with non-zero "running". This deviation, by itself, serves as a powerful tool to discriminate among theories for the origin of cosmological structures such as inflationary models. Here, we use a minimally-parametric smoothing spline technique to reconstruct the shape of the primordial power spectrum. This technique is well-suited to search for smooth features in the primordial power spectrum such as deviations from scale invariance or a running spectral index, although it would recover sharp features of high statistical significance. We use the WMAP 3 year results in combination with data from a suite of higher resolution CMB experiments (including the latest ACBAR 2008 release), as well as large-scale structure data from SDSS and 2dFGRS. We employ cross-validation to assess, using the data themselves, the optimal amount of smoothness in the primordial power spectrum consistent with the data. This minimally-parametric reconstruction supports the evidence for a power law primordial power spectrum with a red tilt, but not for deviations from a power law power spectrum. Smooth variations in the primordial power spectrum are not significantly degenerate with the other cosmological parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Primordial Black Holes, Eternal Inflation, and the Inflationary Parameter Space after WMAP5

We consider constraints on inflation driven by a single, minimally coupled scalar field in the light of the WMAP5 dataset, as well as ACBAR and the SuperNova Legacy Survey. We use the Slow Roll Reconstruction algorithm to derive optimal constraints on the inflationary parameter space. The scale dependence in the slope of the scalar spectrum permitted by WMAP5 is large enough to lead to viable models where the small scale perturbations have a substantial amplitude when extrapolated to the end of inflation. We find that excluding parameter values which would cause the overproduction of primordial black holes or even the onset of eternal inflation leads to potentially significant constraints on the slow roll parameters. Finally, we present a more sophisticated approach to including priors based on the total duration of inflation, and discuss the resulting restrictions on the inflationary parameter space.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Deciphering Inflation with Gravitational Waves: Cosmic Microwave Background Polarization vs. Direct Detection with Laser Interferometers

A detection of the primordial gravitational wave background is considered to be the ``smoking-gun '' evidence for inflation. While super-horizon waves are probed with cosmic microwave background (CMB) polarization, the relic background will be studied with laser interferometers. The long lever arm spanned by the two techniques improves constraints on the inflationary potential and validation of consistency relations expected under inflation. If gravitational waves with a tensor-to-scalar amplitude ratio greater than 0.01 are detected by the CMB, then a direct detection experiment with a sensitivity consistent with current concept studies should be pursued vigorously. If no primordial tensors are detected by the CMB, a direct detection experiment to understand the simplest form of inflation must have a sensitivity improved by two to three orders of magnitude over current plans.

0 citations0 reviewsNo code linkedOpen access