Source author record

Hiranya Peiris

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

18works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

18 published item(s)

preprint2026arXiv

StAD: Stein Amortized Divergence for Fast Likelihoods with Diffusion and Flow

Diffusion and flow-based models are ubiquitously used for generative modelling and density estimation. They admit a deterministic probability flow ordinary differential equation (PF-ODE), analogous to continuous normalizing flows (CNFs), which describes the transport of the probability mass. Obtaining the likelihood from these models is of interest to many workflows, especially Bayesian analysis, and requires solving the trace of the Jacobian to compute the divergence of the learned PF-ODE, which is either $\mathcal{O}(D^2)$ to compute exactly or $\mathcal{O}(D)$ with a noisy estimate. We introduce StAD, a new distillation method to predict and learn the divergence of the PF-ODE using the Langevin-Stein operator without ever computing the Jacobian. We show that our method is competitive with the Hutchinson and Hutch++ on CIFAR-10, ImageNet and other density estimation tasks, consistently improving the variance and speed of the likelihood predictions compared to the Hutchinson. We additionally show our method will generalize to a varied class of generative models, and show that under some regularity conditions these learned vector fields can be made to satisfy the Stein class.

preprint2022arXiv

Hierarchical Bayesian inference of photometric redshifts with stellar population synthesis models

We present a Bayesian hierarchical framework to analyze photometric galaxy survey data with stellar population synthesis (SPS) models. Our method couples robust modeling of spectral energy distributions with a population model and a noise model to characterize the statistical properties of the galaxy populations and real observations, respectively. By self-consistently inferring all model parameters, from high-level hyper-parameters to SPS parameters of individual galaxies, one can separate sources of bias and uncertainty in the data.We demonstrate the strengths and flexibility of this approach by deriving accurate photometric redshifts for a sample of spectroscopically-confirmed galaxies in the COSMOS field, all with 26-band photometry and spectroscopic redshifts. We achieve a performance competitive with publicly-released photometric redshift catalogs based on the same data. Prior to this work, this approach was computationally intractable in practice due to the heavy computational load of SPS model calls; we overcome this challenge using with neural emulators. We find that the largest photometric residuals are associated with poor calibration for emission line luminosities and thus build a framework to mitigate these effects. This combination of physics-based modeling accelerated with machine learning paves the path towards meeting the stringent requirements on the accuracy of photometric redshift estimation imposed by upcoming cosmological surveys. The approach also has the potential to create new links between cosmology and galaxy evolution through the analysis of photometric datasets.

preprint2021arXiv

Optimization of the Observing Cadence for the Rubin Observatory Legacy Survey of Space and Time: a pioneering process of community-focused experimental design

Vera C. Rubin Observatory is a ground-based astronomical facility under construction, a joint project of the National Science Foundation and the U.S. Department of Energy, designed to conduct a multi-purpose 10-year optical survey of the southern hemisphere sky: the Legacy Survey of Space and Time. Significant flexibility in survey strategy remains within the constraints imposed by the core science goals of probing dark energy and dark matter, cataloging the Solar System, exploring the transient optical sky, and mapping the Milky Way. The survey's massive data throughput will be transformational for many other astrophysics domains and Rubin's data access policy sets the stage for a huge potential users' community. To ensure that the survey science potential is maximized while serving as broad a community as possible, Rubin Observatory has involved the scientific community at large in the process of setting and refining the details of the observing strategy. The motivation, history, and decision-making process of this strategy optimization are detailed in this paper, giving context to the science-driven proposals and recommendations for the survey strategy included in this Focus Issue.

preprint2021arXiv

The Impact of Observing Strategy on Cosmological Constraints with LSST

The generation-defining Vera C. Rubin Observatory will make state-of-the-art measurements of both the static and transient universe through its Legacy Survey for Space and Time (LSST). With such capabilities, it is immensely challenging to optimize the LSST observing strategy across the survey's wide range of science drivers. Many aspects of the LSST observing strategy relevant to the LSST Dark Energy Science Collaboration, such as survey footprint definition, single visit exposure time and the cadence of repeat visits in different filters, are yet to be finalized. Here, we present metrics used to assess the impact of observing strategy on the cosmological probes considered most sensitive to survey design; these are large-scale structure, weak lensing, type Ia supernovae, kilonovae and strong lens systems (as well as photometric redshifts, which enable many of these probes). We evaluate these metrics for over 100 different simulated potential survey designs. Our results show that multiple observing strategy decisions can profoundly impact cosmological constraints with LSST; these include adjusting the survey footprint, ensuring repeat nightly visits are taken in different filters and enforcing regular cadence. We provide public code for our metrics, which makes them readily available for evaluating further modifications to the survey design. We conclude with a set of recommendations and highlight observing strategy factors that require further research.

preprint2020arXiv

GenetIC -- a new initial conditions generator to support genetically modified zoom simulations

We present genetIC, a new code for generating initial conditions for cosmological N-body simulations. The code allows precise, user-specified alterations to be made to arbitrary regions of the simulation (while maintaining consistency with the statistical ensemble). These "genetic modifications" allow, for example, the history, mass, or environment of a target halo to be altered in order to study the effect on their evolution. The code natively supports initial conditions with nested zoom regions at progressively increasing resolution. Modifications in the high-resolution region must propagate self-consistently onto the lower resolution grids; to enable this while maintaining a small memory footprint, we introduce a Fourier-space filtering approach to generating fields at variable resolution. Due to a close correspondence with modifications, constrained initial conditions can also be produced by genetIC (for example with the aim of matching structures in the local Universe). We test the accuracy of modifications performed within zoom initial conditions. The code achieves sub-percent precision, which is easily sufficient for current applications in galaxy formation.

preprint2020arXiv

SPECULATOR: Emulating stellar population synthesis for fast and accurate galaxy spectra and photometry

We present SPECULATOR - a fast, accurate, and flexible framework for emulating stellar population synthesis (SPS) models for predicting galaxy spectra and photometry. For emulating spectra, we use principal component analysis to construct a set of basis functions, and neural networks to learn the basis coefficients as a function of the SPS model parameters. For photometry, we parameterize the magnitudes (for the filters of interest) as a function of SPS parameters by a neural network. The resulting emulators are able to predict spectra and photometry under both simple and complicated SPS model parameterizations to percent-level accuracy, giving a factor of $10^3$-$10^4$ speed up over direct SPS computation. They have readily-computable derivatives, making them amenable to gradient-based inference and optimization methods. The emulators are also straightforward to call from a GPU, giving an additional order-of-magnitude speed-up. Rapid SPS computations delivered by emulation offers a massive reduction in the computational resources required to infer the physical properties of galaxies from observed spectra or photometry and simulate galaxy populations under SPS models, whilst maintaining the accuracy required for a range of applications.

preprint2014arXiv

Constraints on ionising photon production from the large-scale Lyman-alpha forest

Recent work has shown that the z~2.5 Lyman-alpha forest on large scales encodes information about the galaxy and quasar populations that keep the intergalactic medium photoionized. We present the first forecasts for constraining the populations with data from current and next-generation surveys. At a minimum the forest should tell us whether galaxies or, conversely, quasars dominate the photon production. The number density and clustering strength of the ionising sources might be estimated to sub-10% precision with a DESI-like survey if degeneracies (e.g., with the photon mean-free-path, small-scale clustering power normalization and potentially other astrophysical effects) can be broken by prior information. We demonstrate that, when inhomogeneous ionisation is correctly handled, constraints on dark energy do not degrade.

preprint2014arXiv

PRISM (Polarized Radiation Imaging and Spectroscopy Mission): An Extended White Paper

PRISM (Polarized Radiation Imaging and Spectroscopy Mission) was proposed to ESA in May 2013 as a large-class mission for investigating within the framework of the ESA Cosmic Vision program a set of important scientific questions that require high resolution, high sensitivity, full-sky observations of the sky emission at wavelengths ranging from millimeter-wave to the far-infrared. PRISM's main objective is to explore the distant universe, probing cosmic history from very early times until now as well as the structures, distribution of matter, and velocity flows throughout our Hubble volume. PRISM will survey the full sky in a large number of frequency bands in both intensity and polarization and will measure the absolute spectrum of sky emission more than three orders of magnitude better than COBE FIRAS. The aim of this Extended White Paper is to provide a more detailed overview of the highlights of the new science that will be made possible by PRISM

preprint2013arXiv

Constraining Monodromy Inflation

We use cosmic microwave background (CMB) data from the 9-year WMAP release to derive constraints on monodromy inflation, which is characterized by a linear inflaton potential with a periodic modulation. We identify two possible periodic modulations that significantly improve the fit, but it is unclear whether this improvement is associated with a "signal" or with the scatter in the measured angular power spectrum. The unmodulated potential is mildly favored by standard model selection criteria. A modulated inflationary potential can generate substantial primordial non-Gaussianity, of a specific and characteristic form. For the best-fit parameters to the WMAP angular power spectrum, the corresponding non-Gaussianity might be detectable in upcoming CMB data, allowing nontrivial consistency checks on the predictions of a modulated inflationary potential.

preprint2012arXiv

Bayesian Analysis of Inflation II: Model Selection and Constraints on Reheating

We discuss the model selection problem for inflationary cosmology. We couple ModeCode, a publicly-available numerical solver for the primordial perturbation spectra, to the nested sampler MultiNest, in order to efficiently compute Bayesian evidence. Particular attention is paid to the specification of physically realistic priors, including the parametrization of the post-inflationary expansion and associated thermalization scale. It is confirmed that while present-day data tightly constrains the properties of the power spectrum, it cannot usefully distinguish between the members of a large class of simple inflationary models. We also compute evidence using a simulated Planck likelihood, showing that while Planck will have more power than WMAP to discriminate between inflationary models, it will not definitively address the inflationary model selection problem on its own. However, Planck will place very tight constraints on any model with more than one observationally-distinct inflationary regime -- e.g. the large- and small-field limits of the hilltop inflation model -- and put useful limits on different reheating scenarios for a given model.

preprint2010arXiv

The Herschel-SPIRE Legacy Survey (HSLS): the scientific goals of a shallow and wide submillimeter imaging survey with SPIRE

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

preprint2009arXiv

Lecture notes on the physics of cosmic microwave background anisotropies

We review the theory of the temperature anisotropy and polarization of the cosmic microwave background (CMB) radiation, and describe what we have learned from current CMB observations. In particular, we discuss how the CMB is being used to provide precise measurements of the composition and geometry of the observable universe, and to constrain the physics of the early universe. We also briefly review the physics of the small-scale CMB fluctuations generated during and after the epoch of reionization, and which are the target of a new breed of arcminute-resolution instruments.

preprint2006arXiv

Considerations in optimizing CMB polarization experiments to constrain inflationary physics

We quantify the limiting factors in optimizing current-technology cosmic microwave background (CMB) polarization experiments in order to learn about inflationary physics. We consider space-based, balloon-borne and ground-based experiments. We find that foreground contamination and residuals from foreground subtraction are ultimately the limiting factors in detecting a primordial gravity wave signal. For full-sky space-based experiments, these factors hinder the detection of tensor-to-scalar ratios of r < 10^{-3} on large scales, while for ground-based experiments these factors impede the ability to apply delensing techniques. We consider ground-based/balloon-borne experiments of currently planned or proposed designs and find that it is possible for a value of r=0.01 to be measured at ~ 3-sigma level. A small space-based CMB polarization experiment, with current detector technology and full sky coverage, can detect r ~ 10^{-3} at the ~ 3-sigma level, but a markedly improved knowledge of polarized foregrounds is needed. We advocate using as wide a frequency coverage as possible in order to carry out foreground subtraction at the percent level, which is necessary to measure such a small primordial tensor amplitude. To produce a clearly detectable (>3-sigma) tensor component in a realistic CMB experiment, inflation must either involve large-field variations, Δϕ>~ 1 or multi-field/hybrid models. Hybrid models can be easily distinguished from large-field models due to their blue scalar spectral index. Therefore, an observation of a tensor/scalar ratio and n < 1 in future experiments with the characteristics considered here may be an indication that inflation is being driven by some physics in which the inflaton cannot be described as a fundamental field.

preprint2006arXiv

Implications of a Running Spectral Index for Slow Roll Inflation

We analyze the weak (2 sigma) evidence for a running spectral index seen in the three-year WMAP dataset and its implications for single field, slow roll inflation. We assume that the running is comparable to the central value found from the WMAP data analysis, and use the Hubble Slow Roll formalism to follow the evolution of the slow roll parameters. For all parameter choices consistent with a large, negative running, single field, slow roll inflation lasts less than 30 efolds after CMB scales leave the horizon. Thus, a definitive observation of a large negative running would imply that any inflationary phase requires multiple fields or the breakdown of slow roll. Alternatively, if single field, slow roll inflation is sources the primordial fluctuations, we can expect the observed running to move much closer to zero as the CMB is measured more accurately at small angular scales.

preprint2006arXiv

Recovering the Inflationary Potential and Primordial Power Spectrum With a Slow Roll Prior: Methodology and Application to WMAP 3 Year Data

We introduce a new method for applying an inflationary prior to a cosmological dataset that includes relations between observables at arbitrary order in the slow roll expansion. The process is based on the inflationary flow equations, and the slow roll parameters appear explicitly in the cosmological parameter set. We contrast our method to other ways of imposing an inflationary prior on a cosmological dataset, and argue that this method is ideal for use with heterogeneous datasets. In particular, it would be well suited to exploiting any direct detection of fundamental tensor modes by a BBO-style mission. To demonstrate the practical use of this method we apply it to the WMAPI+All dataset, and the newly released WMAPII dataset on its own and together with the SDSS data. We find that all basic classes of single field inflationary models are still allowed at the 1-2sigma level, but the overall parameter space is sharply constrained. In particular, we find evidence that the combination of WMAPII+SDSS is sensitive to effects arising from terms that are quadratic in the two leading-order slow roll parameters.

preprint2006arXiv

Slow Roll Reconstruction: Constraints on Inflation from the 3 Year WMAP Dataset

We study the constraints on the inflationary parameter space derived from the 3 year WMAP dataset using ``slow roll reconstruction'', using the SDSS galaxy power spectrum to gain further leverage where appropriate. This approach inserts the inflationary slow roll parameters directly into a Monte Carlo Markov chain estimate of the cosmological parameters, and uses the inflationary flow hierarchy to compute the parameters' scale-dependence. We work with the first three parameters (epsilon, eta and xi) and pay close attention to the possibility that the 3 year WMAP dataset contains evidence for a ``running'' spectral index, which is dominated by the xi term. Mirroring the WMAP team's analysis we find that the permitted distribution of xi is broad, and centered away from zero. However, when we require that inflationary parameters yield at least 30 additional e-folds of inflation after the largest observable scales leave the horizon, the bounds on xi tighten dramatically. We make use of the absence of an explicit pivot scale in the slow roll reconstruction formalism to determine the dependence of the computed parameter distributions on the pivot. We show that the choice of pivot has a significant effect on the inferred constraints on the inflationary variables, and the spectral index and running derived from them. Finally, we argue that the next round of cosmological data can be expected to place very stringent constraints on the region of parameter space open to single field models of slow roll inflation.

preprint2005arXiv

Boundary Effective Field Theory and Trans-Planckian Perturbations: Astrophysical Implications

We contrast two approaches to calculating trans-Planckian corrections to the inflationary perturbation spectrum: the New Physics Hypersurface [NPH] model, in which modes are normalized when their physical wavelength first exceeds a critical value, and the Boundary Effective Field Theory [BEFT] approach, where the initial conditions for all modes are set at the same time, and modified by higher dimensional operators enumerated via an effective field theory calculation. We show that these two approaches -- as currently implemented -- lead to radically different expectations for the trans-Planckian corrections to the CMB and emphasize that in the BEFT formalism we expect the perturbation spectrum to be dominated by quantum gravity corrections for all scales shorter than some critical value. Conversely, in the NPH case the quantum effects only dominate the longest modes that are typically much larger than the present horizon size. Furthermore, the onset of the breakdown in the standard inflationary perturbation calculation predicted by the BEFT formalism is likely to be associated with a feature in the perturbation spectrum, and we discuss the observational signatures of this feature in both CMB and large scale structure observations. Finally, we discuss possible modifications to both calculational frameworks that would resolve the contradictions identified here.

preprint2005arXiv

Observing Trans-Planckian Signatures in the Cosmic Microwave Background

We examine the constraints cosmological observations can place on any trans-Planckian corrections to the primordial spectrum of perturbations underlying the anisotropies in the Cosmic Microwave Background. We focus on models of trans-Planckian physics which lead to a modulated primordial spectrum. Rather than looking at a generic modulated spectrum, our calculations are based on a specific model, and are intended as a case study for the sort of constraints one could hope to apply on a well-motivated model of trans-Planckian physics. We present results for three different approaches -- a grid search in a subset of the overall parameter space, a Fisher matrix estimate of the likely error ellipses, and a Monte Carlo Markov Chain fit to a simulated CMB sky. As was seen in previous analyses, the likelihood space has multiple peaks, and we show that their distribution can be reproduced via a simple semi-analytic argument. All three methods lead to broadly similar results. We vary 10 cosmological parameters (including two related to the trans-Planckian terms) and show that the amplitude of the tensor perturbations is directly correlated with the detectability of any trans-Planckian modulation. We argue that this is likely to be true for any trans-Planckian modulation in the paradigm of slow-roll inflation. For the specific case we consider, we conclude that if the tensor to scalar ratio, $r \sim 0.15$, the ratio between the inflationary Hubble scale $H$, and the scale of new physics $M$ has to be on the order of 0.004 if the modulation is detectable at the 2$σ$ level. For a lower value of $r$, the bound on $H/M$ becomes looser.