Source author record

Philip Bull

Philip Bull 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

39works
6topics
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

39 published item(s)

preprint2026arXiv

Angular correlation functions of bright Lyman-break galaxies at $\mathbf{3 \lesssim z \lesssim 5}$

We investigate the clustering of Lyman-break galaxies at redshifts of 3 $\lesssim z \lesssim$ 5 within the COSMOS field by measuring the angular two-point correlation function. Our robust sample of $\sim$60,000 bright ($m_{\rm UV}\lesssim 27$) Lyman-break galaxies was selected based on spectral energy distribution fitting across 14 photometric bands spanning optical and near-infrared wavelengths. We constrained both the 1- and 2-halo terms at separations up to 300 arcsec, finding an excess in the correlation function at scales corresponding to $<20$ kpc, consistent with enhancement due to clumps in the same galaxy or interactions on this scale. We then performed Bayesian model fits on the correlation functions to infer the Halo Occupation Distribution parameters, star formation duty cycle, and galaxy bias in three redshift bins. We examined several cases where different combinations of parameters were varied, showing that our data can constrain the slope of the satellite occupation function, which previous studies have fixed. For an $M_{\rm{UV}}$-limited sub-sample, we found galaxy bias values of $b_g=3.18^{+0.14}_{-0.14}$ at $z\simeq3$, $b_g=3.58^{+0.27}_{-0.29}$ at $z\simeq4$, $b_g=4.27^{+0.25}_{-0.26}$ at $z\simeq5$. The duty cycle values are $0.62^{+0.25}_{-0.26}$, $0.40^{+0.34}_{-0.22}$, and $0.39^{+0.31}_{-0.20}$, respectively. These results suggest that, as the redshift increases, there is a slight decrease in the host halo masses and a shorter timescale for star formation in bright galaxies, at a fixed rest-frame UV luminosity threshold.

preprint2025arXiv

RHINO: A large horn antenna for detecting the 21cm global signal

The sky-averaged brightness temperature of the 21cm line from neutral hydrogen provides a sensitive probe of the thermal state of the intergalactic medium, particularly before and during Cosmic Dawn and the Epoch of Reionisation. This `global signal' is faint, on the order of tens to hundreds of millikelvin, and spectrally relatively smooth, making it exceedingly difficult to disentangle from foreground radio emission and instrumental artefacts. In this paper, we introduce RHINO, an experiment based around a large horn antenna operating from 60-85 MHz. Horn antennas are highly characterisable and provide excellent shielding from their immediate environment, which are potentially decisive advantages when it comes to the beam measurement and modelling problems that are particularly challenging for this kind of experiment. The system also includes a novel continuous wave calibration source to control correlated gain fluctuations, allowing continuous monitoring of the overall gain level without needing to rapidly switch between the sky and a calibration source. Here, we describe the basic RHINO concept, including the antenna design, EM simulations, and receiver electronics. We use a basic simulation and analysis pipeline to study the impact of the limited bandwidth on recovery of physical 21cm global signal model parameters, and discuss a basic calibration scheme that incorporates the continuous wave signal. Finally, we report on the current state of a scaled-down prototype system under construction at Jodrell Bank Observatory.

preprint2022arXiv

Automated Detection of Antenna Malfunctions in Large-N Interferometers: A Case Study with the Hydrogen Epoch of Reionization Array

We present a framework for identifying and flagging malfunctioning antennas in large radio interferometers. We outline two distinct categories of metrics designed to detect outliers along known failure modes of large arrays: cross-correlation metrics, based on all antenna pairs, and auto-correlation metrics, based solely on individual antennas. We define and motivate the statistical framework for all metrics used, and present tailored visualizations that aid us in clearly identifying new and existing systematics. We implement these techniques using data from 105 antennas in the Hydrogen Epoch of Reionization Array (HERA) as a case study. Finally, we provide a detailed algorithm for implementing these metrics as flagging tools on real data sets.

preprint2022arXiv

Characterizing line-of-sight variability of polarized dust emission with future CMB experiments

Galactic dust emission is often accounted for in cosmic microwave background (CMB) analyses by fitting a simple two-parameter modified blackbody (MBB) model in each pixel, which nominally accounts for the temperature and opacity of the dust. While this may be a good approximation for individual dust clouds, typically a number of such clouds are found along each line of sight and within each angular pixel, resulting in a superposition of their spectra. In this paper, we study the effects of this superposition on pixel-based foreground fitting strategies by modelling the spectral energy distribution (SED) in each pixel as the integral of individual MBB spectra over various physically-motivated statistical distributions of dust cloud properties. We show that fitting these SEDs with the simple two-parameter MBB model generally results in unbiased estimates of the CMB Stokes Q and U amplitudes in each pixel, unless there are significant changes in both the dust SED and polarization angle along the line of sight, in which case significant ($ > 10σ$) biases are observed in an illustrative model. We also find that the best-fit values of the dust temperature, $T_d$, and spectral index, $β_d$, are significantly biased away from the mean/median of the corresponding statistical distributions when the distributions are broad, suggesting that MBB model fits can give an unrepresentative picture of the physical properties of the dust at microwave wavelengths if not interpreted carefully. Using a Fisher matrix analysis, we also determine the experimental sensitivity required to recover the parameters of the $T_d$ and $β_d$ distributions themselves by fitting a probabilistic MBB model, finding that only the parameters of broad distributions can be measured by SED fitting on a single line of sight.

preprint2022arXiv

Measurements of the diffuse Galactic synchrotron spectral index and curvature from MeerKLASS pilot data

21cm intensity mapping experiments are bringing an influx of high spectral resolution observational data in the $\sim100$ MHz $- 1$ GHz regime. We use pilot $971-1075$ MHz data from MeerKAT in single-dish mode, recently used to test the calibration and data reduction scheme of the upcoming MeerKLASS survey, to probe the spectral index of diffuse synchrotron emission below 1 GHz within $145^{\circ} < α< 180^{\circ}$, $-1^{\circ} < δ< 8^{\circ}$. Through comparisons with data from the OVRO Long Wavelength Array and the Maipu and MU surveys, we find an average spectral index of $-2.75 < β< -2.71$ between 45 and 1055 MHz. By fitting for spectral curvature with a spectral index of the form $β+ c \, {\rm{ln}}(ν/ 73~{\rm MHz})$, we measure $β= -2.55 \pm 0.13$ and $c = -0.12 \pm 0.05$ within our target field. Our results are in good agreement (within $1σ$) with existing measurements from experiments such as ARCADE2 and EDGES. These results show the calibration accuracy of current data and demonstrate that MeerKLASS will also be capable of achieving a secondary science goal of probing the interstellar medium.

preprint2021arXiv

First Results from HERA Phase I: Upper Limits on the Epoch of Reionization 21 cm Power Spectrum

We report upper-limits on the Epoch of Reionization (EoR) 21 cm power spectrum at redshifts 7.9 and 10.4 with 18 nights of data ($\sim36$ hours of integration) from Phase I of the Hydrogen Epoch of Reionization Array (HERA). The Phase I data show evidence for systematics that can be largely suppressed with systematic models down to a dynamic range of $\sim10^9$ with respect to the peak foreground power. This yields a 95% confidence upper limit on the 21 cm power spectrum of $Δ^2_{21} \le (30.76)^2\ {\rm mK}^2$ at $k=0.192\ h\ {\rm Mpc}^{-1}$ at $z=7.9$, and also $Δ^2_{21} \le (95.74)^2\ {\rm mK}^2$ at $k=0.256\ h\ {\rm Mpc}^{-1}$ at $z=10.4$. At $z=7.9$, these limits are the most sensitive to-date by over an order of magnitude. While we find evidence for residual systematics at low line-of-sight Fourier $k_\parallel$ modes, at high $k_\parallel$ modes we find our data to be largely consistent with thermal noise, an indicator that the system could benefit from deeper integrations. The observed systematics could be due to radio frequency interference, cable sub-reflections, or residual instrumental cross-coupling, and warrant further study. This analysis emphasizes algorithms that have minimal inherent signal loss, although we do perform a careful accounting in a companion paper of the small forms of loss or bias associated with the pipeline. Overall, these results are a promising first step in the development of a tuned, instrument-specific analysis pipeline for HERA, particularly as Phase II construction is completed en route to reaching the full sensitivity of the experiment.

preprint2021arXiv

Observing relativistic features in large-scale structure surveys -- I: Multipoles of the power spectrum

Planned efforts to probe the largest observable distance scales in future cosmological surveys are motivated by a desire to detect relic correlations left over from inflation, and the possibility of constraining novel gravitational phenomena beyond General Relativity (GR). On such large scales, the usual Newtonian approaches to modelling summary statistics like the power spectrum and bispectrum are insufficient, and we must consider a fully relativistic and gauge-independent treatment of observables such as galaxy number counts in order to avoid subtle biases, e.g. in the determination of the $f_{\rm NL}$ parameter. In this work, we present an initial application of an analysis pipeline capable of accurately modelling and recovering relativistic spectra and correlation functions. As a proof of concept, we focus on the non-zero dipole of the redshift-space power spectrum that arises in the cross-correlation of different mass bins of dark matter halos, using strictly gauge-independent observable quantities evaluated on the past light cone of a fully relativistic N-body simulation in a redshift bin $1.7 \le z \le 2.9$. We pay particular attention to the correct estimation of power spectrum multipoles, comparing different methods of accounting for complications such as the survey geometry (window function) and evolution/bias effects on the past light cone, and discuss how our results compare with previous attempts at extracting novel GR signatures from relativistic simulations.

preprint2021arXiv

Validation of the HERA Phase I Epoch of Reionization 21 cm Power Spectrum Software Pipeline

We describe the validation of the HERA Phase I software pipeline by a series of modular tests, building up to an end-to-end simulation. The philosophy of this approach is to validate the software and algorithms used in the Phase I upper limit analysis on wholly synthetic data satisfying the assumptions of that analysis, not addressing whether the actual data meet these assumptions. We discuss the organization of this validation approach, the specific modular tests performed, and the construction of the end-to-end simulations. We explicitly discuss the limitations in scope of the current simulation effort. With mock visibility data generated from a known analytic power spectrum and a wide range of realistic instrumental effects and foregrounds, we demonstrate that the current pipeline produces power spectrum estimates that are consistent with known analytic inputs to within thermal noise levels (at the 2 sigma level) for k > 0.2 h/Mpc for both bands and fields considered. Our input spectrum is intentionally amplified to enable a strong `detection' at k ~0.2 h/Mpc -- at the level of ~25 sigma -- with foregrounds dominating on larger scales, and thermal noise dominating at smaller scales. Our pipeline is able to detect this amplified input signal after suppressing foregrounds with a dynamic range (foreground to noise ratio) of > 10^7. Our validation test suite uncovered several sources of scale-independent signal loss throughout the pipeline, whose amplitude is well-characterized and accounted for in the final estimates. We conclude with a discussion of the steps required for the next round of data analysis.

preprint2020arXiv

Absolute Calibration Strategies for the Hydrogen Epoch of Reionization Array and Their Impact on the 21 cm Power Spectrum

We discuss absolute calibration strategies for Phase I of the Hydrogen Epoch of Reionization Array (HERA), which aims to measure the cosmological 21 cm signal from the Epoch of Reionization (EoR). HERA is a drift-scan array with a 10 degree wide field of view, meaning bright, well-characterized point source transits are scarce. This, combined with HERA's redundant sampling of the uv plane and the modest angular resolution of the Phase I instrument, make traditional sky-based and self-calibration techniques difficult to implement with high dynamic range. Nonetheless, in this work we demonstrate calibration for HERA using point source catalogues and electromagnetic simulations of its primary beam. We show that unmodeled diffuse flux and instrumental contaminants can corrupt the gain solutions, and present a gain smoothing approach for mitigating their impact on the 21 cm power spectrum. We also demonstrate a hybrid sky and redundant calibration scheme and compare it to pure sky-based calibration, showing only a marginal improvement to the gain solutions at intermediate delay scales. Our work suggests that the HERA Phase I system can be well-calibrated for a foreground-avoidance power spectrum estimator by applying direction-independent gains with a small set of degrees of freedom across the frequency and time axes.

preprint2019arXiv

Core Cosmology Library: Precision Cosmological Predictions for LSST

The Core Cosmology Library (CCL) provides routines to compute basic cosmological observables to a high degree of accuracy, which have been verified with an extensive suite of validation tests. Predictions are provided for many cosmological quantities, including distances, angular power spectra, correlation functions, halo bias and the halo mass function through state-of-the-art modeling prescriptions available in the literature. Fiducial specifications for the expected galaxy distributions for the Large Synoptic Survey Telescope (LSST) are also included, together with the capability of computing redshift distributions for a user-defined photometric redshift model. A rigorous validation procedure, based on comparisons between CCL and independent software packages, allows us to establish a well-defined numerical accuracy for each predicted quantity. As a result, predictions for correlation functions of galaxy clustering, galaxy-galaxy lensing and cosmic shear are demonstrated to be within a fraction of the expected statistical uncertainty of the observables for the models and in the range of scales of interest to LSST. CCL is an open source software package written in C, with a python interface and publicly available at https://github.com/LSSTDESC/CCL.

preprint2019arXiv

Mitigating Internal Instrument Coupling II: A Method Demonstration with the Hydrogen Epoch of Reionization Array

We present a study of internal reflection and cross coupling systematics in Phase I of the Hydrogen Epoch of Reionization Array (HERA). In a companion paper, we outlined the mathematical formalism for such systematics and presented algorithms for modeling and removing them from the data. In this work, we apply these techniques to data from HERA's first observing season as a method demonstration. The data show evidence for systematics that, without removal, would hinder a detection of the 21 cm power spectrum for the targeted EoR line-of-sight modes in the range 0.2 < k_parallel < 0.5\ h^-1 Mpc. After systematic removal, we find we can recover these modes in the power spectrum down to the integrated noise-floor of a nightly observation, achieving a dynamic range in the EoR window of 10^-6 in power (mK^2 units) with respect to the bright galactic foreground signal. In the absence of other systematics and assuming the systematic suppression demonstrated here continues to lower noise levels, our results suggest that fully-integrated HERA Phase I may have the capacity to set competitive upper limits on the 21 cm power spectrum. For future observing seasons, HERA will have upgraded analog and digital hardware to better control these systematics in the field.

preprint2018arXiv

Cosmology with Phase 1 of the Square Kilometre Array; Red Book 2018: Technical specifications and performance forecasts

We present a detailed overview of the cosmological surveys that will be carried out with Phase 1 of the Square Kilometre Array (SKA1), and the science that they will enable. We highlight three main surveys: a medium-deep continuum weak lensing and low-redshift spectroscopic HI galaxy survey over 5,000 sqdeg; a wide and deep continuum galaxy and HI intensity mapping survey over 20,000 sqdeg from z = 0.35 - 3; and a deep, high-redshift HI intensity mapping survey over 100 sqdeg from z = 3 - 6. Taken together, these surveys will achieve an array of important scientific goals: measuring the equation of state of dark energy out to z ~ 3 with percent-level precision measurements of the cosmic expansion rate; constraining possible deviations from General Relativity on cosmological scales by measuring the growth rate of structure through multiple independent methods; mapping the structure of the Universe on the largest accessible scales, thus constraining fundamental properties such as isotropy, homogeneity, and non-Gaussianity; and measuring the HI density and bias out to z = 6. These surveys will also provide highly complementary clustering and weak lensing measurements that have independent systematic uncertainties to those of optical surveys like LSST and Euclid, leading to a multitude of synergies that can improve constraints significantly beyond what optical or radio surveys can achieve on their own. This document, the 2018 Red Book, provides reference technical specifications, cosmological parameter forecasts, and an overview of relevant systematic effects for the three key surveys, and will be regularly updated by the Cosmology Science Working Group in the run up to start of operations and the Key Science Programme of SKA1.

preprint2016arXiv

Beyond $Λ$CDM: Problems, solutions, and the road ahead

Despite its continued observational successes, there is a persistent (and growing) interest in extending cosmology beyond the standard model, $Λ$CDM. This is motivated by a range of apparently serious theoretical issues, involving such questions as the cosmological constant problem, the particle nature of dark matter, the validity of general relativity on large scales, the existence of anomalies in the CMB and on small scales, and the predictivity and testability of the inflationary paradigm. In this paper, we summarize the current status of $Λ$CDM as a physical theory, and review investigations into possible alternatives along a number of different lines, with a particular focus on highlighting the most promising directions. While the fundamental problems are proving reluctant to yield, the study of alternative cosmologies has led to considerable progress, with much more to come if hopes about forthcoming high-precision observations and new theoretical ideas are fulfilled.

preprint2016arXiv

Distinguishing screening mechanisms with environment-dependent velocity statistics

Alternative theories of gravity typically invoke an environment-dependent screening mechanism to allow phenomenologically interesting deviations from general relativity (GR) to manifest on larger scales, while reducing to GR on small scales. The observation of the transition from screened to unscreened behavior would be compelling evidence for beyond-GR physics. We show that pairwise peculiar velocity statistics, in particular the relative radial velocity dispersion, $σ_\parallel$ , can be used to observe this transition when they are binned by some measure of halo environment. We established this by measuring the radial velocity dispersion between pairs of halos in N-body simulations for three $f(R)$ gravity and four symmetron models. We developed an estimator involving only line-of-sight velocities to show that this quantity is observable, and binned the results in halo mass, ambient density, and the isolatedness of halos. Ambient density is found to be the most relevant measure of environment; it is distinct from isolatedness, and correlates well with theoretical expectations for the symmetron model. By binning $σ_\parallel$ in ambient density, we find a strong environment-dependent signature for the symmetron models, with the velocities showing a clear transition from GR to non-GR behavior. No such transition is observed for $f(R)$, as the relevant scales are deep in the unscreened regime. Observations of the relative radial velocity dispersion in forthcoming peculiar velocity surveys, if binned appropriately by environment, therefore offer a valuable way of detecting the screening signature of modified gravity.

preprint2016arXiv

Extending cosmological tests of General Relativity with the Square Kilometre Array

Tests of general relativity (GR) are still in their infancy on cosmological scales, but forthcoming experiments promise to greatly improve their precision over a wide range of distance scales and redshifts. One such experiment, the Square Kilometre Array (SKA), will carry out several wide and deep surveys of resolved and unresolved neutral hydrogen (HI) 21cm line-emitting galaxies, mapping a significant fraction of the sky from $0 \le z \lesssim 6$. I present forecasts for the ability of a suite of possible SKA HI surveys to detect deviations from GR by reconstructing the cosmic expansion and growth history. SKA Phase 1 intensity mapping surveys can achieve sub-1\% measurements of $fσ_8$ out to $z\approx 1$, with an SKA1-MID Band 2 survey out to $z \lesssim 0.6$ able to surpass contemporary spectroscopic galaxy surveys such as DESI and Euclid in terms of constraints on modified gravity parameters if challenges such as foreground contamination can be tackled effectively. A more futuristic Phase 2 HI survey of $\sim10^9$ spectroscopic galaxy redshifts would be capable of detecting a $\sim 2\%$ modification of the Poisson equation out to $z\approx 2$.

preprint2016arXiv

Reconstructing cosmic growth with kSZ observations in the era of Stage IV experiments

Future ground-based CMB experiments will generate competitive large-scale structure datasets by precisely characterizing CMB secondary anisotropies over a large fraction of the sky. We describe a method for constraining the growth rate of structure to sub-1% precision out to $z\approx 1$, using a combination of galaxy cluster peculiar velocities measured using the kinetic Sunyaev-Zel'dovich (kSZ) effect, and the velocity field reconstructed from galaxy redshift surveys. We consider only thermal SZ-selected cluster samples, which will consist of $\mathcal{O}(10^4-10^5)$ sources for Stage 3 and 4 CMB experiments respectively. Three different methods for separating the kSZ effect from the primary CMB are compared, including a novel blind "constrained realization" method that improves signal-to-noise by a factor of $\sim 2$ over a commonly-used aperture photometry technique. Measurements of the integrated tSZ $y$-parameter are used to break the kSZ velocity-optical depth degeneracy, and the effects of including CMB polarization and SZ profile uncertainties are also considered. A combination of future Stage 4 experiments should be able to measure the product of the growth and expansion rates, $α\equiv f H$, to better than 1% in bins of $Δz = 0.1$ out to $z \approx 1$ -- competitive with contemporary redshift-space distortion constraints from galaxy surveys.

preprint2016arXiv

Science Impacts of the SPHEREx All-Sky Optical to Near-Infrared Spectral Survey: Report of a Community Workshop Examining Extragalactic, Galactic, Stellar and Planetary Science

SPHEREx is a proposed SMEX mission selected for Phase A. SPHEREx will carry out the first all-sky spectral survey and provide for every 6.2" pixel a spectra between 0.75 and 4.18 $μ$m [with R$\sim$41.4] and 4.18 and 5.00 $μ$m [with R$\sim$135]. The SPHEREx team has proposed three specific science investigations to be carried out with this unique data set: cosmic inflation, interstellar and circumstellar ices, and the extra-galactic background light. It is readily apparent, however, that many other questions in astrophysics and planetary sciences could be addressed with the SPHEREx data. The SPHEREx team convened a community workshop in February 2016, with the intent of enlisting the aid of a larger group of scientists in defining these questions. This paper summarizes the rich and varied menu of investigations that was laid out. It includes studies of the composition of main belt and Trojan/Greek asteroids; mapping the zodiacal light with unprecedented spatial and spectral resolution; identifying and studying very low-metallicity stars; improving stellar parameters in order to better characterize transiting exoplanets; studying aliphatic and aromatic carbon-bearing molecules in the interstellar medium; mapping star formation rates in nearby galaxies; determining the redshift of clusters of galaxies; identifying high redshift quasars over the full sky; and providing a NIR spectrum for most eROSITA X-ray sources. All of these investigations, and others not listed here, can be carried out with the nominal all-sky spectra to be produced by SPHEREx. In addition, the workshop defined enhanced data products and user tools which would facilitate some of these scientific studies. Finally, the workshop noted the high degrees of synergy between SPHEREx and a number of other current or forthcoming programs, including JWST, WFIRST, Euclid, GAIA, K2/Kepler, TESS, eROSITA and LSST.

preprint2016arXiv

Spatial curvature endgame: Reaching the limit of curvature determination

Current constraints on spatial curvature show that it is dynamically negligible: $|Ω_{\rm K}| \lesssim 5 \times 10^{-3}$ (95% CL). Neglecting it as a cosmological parameter would be premature however, as more stringent constraints on $Ω_{\rm K}$ at around the $10^{-4}$ level would offer valuable tests of eternal inflation models and probe novel large-scale structure phenomena. This precision also represents the "curvature floor", beyond which constraints cannot be meaningfully improved due to the cosmic variance of horizon-scale perturbations. In this paper, we discuss what future experiments will need to do in order to measure spatial curvature to this maximum accuracy. Our conservative forecasts show that the curvature floor is unreachable - by an order of magnitude - even with Stage IV experiments, unless strong assumptions are made about dark energy evolution and the $Λ$CDM parameter values. We also discuss some of the novel problems that arise when attempting to constrain a global cosmological parameter like $Ω_{\rm K}$ with such high precision. Measuring curvature down to this level would be an important validation of systematics characterisation in high-precision cosmological analyses.

preprint2015arXiv

A CMB Gibbs sampler for localized secondary anisotropies

As well as primary fluctuations, CMB temperature maps contain a wealth of additional information in the form of secondary anisotropies. Secondary effects that can be identified with individual objects, such as the thermal and kinetic Sunyaev-Zel'dovich (SZ) effects due to galaxy clusters, are difficult to unambiguously disentangle from foreground contamination and the primary CMB however. We develop a Bayesian formalism for rigorously characterising anisotropies that are localised on the sky, taking the TSZ and KSZ effects as an example. Using a Gibbs sampling scheme, we are able to efficiently sample from the joint posterior distribution for a multi-component model of the sky with many thousands of correlated physical parameters. The posterior can then be exactly marginalised to estimate properties of the secondary anisotropies, fully taking into account degeneracies with the other signals in the CMB map. We show that this method is computationally tractable using a simple implementation based on the existing Commander component separation code, and also discuss how other types of secondary anisotropy can be accommodated within our framework.

preprint2015arXiv

A systematic study of Lyman-Alpha transfer through outflowing shells: Model parameter estimation

Outflows promote the escape of Lyman-$α$ (Ly$α$) photons from dusty interstellar media. The process of radiative transfer through interstellar outflows is often modelled by a spherically symmetric, geometrically thin shell of gas that scatters photons emitted by a central Ly$α$ source. Despite its simplified geometry, this `shell model' has been surprisingly successful at reproducing observed Ly$α$ line shapes. In this paper we perform automated line fitting on a set of noisy simulated shell model spectra, in order to determine whether degeneracies exist between the different shell model parameters. While there are some significant degeneracies, we find that most parameters are accurately recovered, especially the HI column density ($N_{\rm HI}$) and outflow velocity ($v_{\rm exp}$). This work represents an important first step in determining how the shell model parameters relate to the actual physical properties of Ly$α$ sources. To aid further exploration of the parameter space, we have made our simulated model spectra available through an interactive online tool.

preprint2015arXiv

Cosmology with a SKA HI intensity mapping survey

HI intensity mapping (IM) is a novel technique capable of mapping the large-scale structure of the Universe in three dimensions and delivering exquisite constraints on cosmology, by using HI as a biased tracer of the dark matter density field. This is achieved by measuring the intensity of the redshifted 21cm line over the sky in a range of redshifts without the requirement to resolve individual galaxies. In this chapter, we investigate the potential of SKA1 to deliver HI intensity maps over a broad range of frequencies and a substantial fraction of the sky. By pinning down the baryon acoustic oscillation and redshift space distortion features in the matter power spectrum -- thus determining the expansion and growth history of the Universe -- these surveys can provide powerful tests of dark energy models and modifications to General Relativity. They can also be used to probe physics on extremely large scales, where precise measurements of spatial curvature and primordial non-Gaussianity can be used to test inflation; on small scales, by measuring the sum of neutrino masses; and at high redshifts where non-standard evolution models can be probed. We discuss the impact of foregrounds as well as various instrumental and survey design parameters on the achievable constraints. In particular we analyse the feasibility of using the SKA1 autocorrelations to probe the large-scale signal.

preprint2015arXiv

Cross correlation surveys with the Square Kilometre Array

By the time that the first phase of the Square Kilometre Array is deployed it will be able to perform state of the art Large Scale Structure (LSS) as well as Weak Gravitational Lensing (WGL) measurements of the distribution of matter in the Universe. In this chapter we concentrate on the synergies that result from cross-correlating these different SKA data products as well as external correlation with the weak lensing measurements available from CMB missions. We show that the Dark Energy figures of merit obtained individually from WGL/LSS measurements and their independent combination is significantly increased when their full cross-correlations are taken into account. This is due to the increased knowledge of galaxy bias as a function of redshift as well as the extra information from the different cosmological dependences of the cross-correlations. We show that the cross-correlation between a spectroscopic LSS sample and a weak lensing sample with photometric redshifts can calibrate these same photometric redshifts, and their scatter, to high accuracy by modelling them as nuisance parameters and fitting them simultaneously cosmology. Finally we show that Modified Gravity parameters are greatly constrained by this cross-correlations because weak lensing and redshift space distortions (from the LSS survey) break strong degeneracies in common parameterisations of modified gravity.

preprint2015arXiv

Cross-correlating 21cm intensity maps with Lyman Break Galaxies in the post-reionization era

We investigate the cross-correlation between the spatial distribution of Lyman Break Galaxies (LBGs) and the 21cm intensity mapping signal at $z\sim[3-5]$. At these redshifts, galactic feedback is supposed to only marginally affect the matter power spectrum, and the neutral hydrogen distribution is independently constrained by quasar spectra. Using a high resolution N-body simulation, populated with neutral hydrogen a posteriori, we forecast for the expected LBG-21cm cross-spectrum and its error for a 21cm field observed by the Square Kilometre Array (SKA1-LOW and SKA1-MID), combined with a spectroscopic LBG survey with the same volume. The cross power can be detected with a signal-to-noise ratio (SNR) up to ~10 times higher (and down to ~4 times smaller scales) than the 21cm auto-spectrum for this set-up, with the SNR depending only very weakly on redshift and the LBG population. We also show that while both the 21cm auto- and LBG-21cm cross-spectra can be reliably recovered after the cleaning of smooth-spectrum foreground contamination, only the cross-power is robust to problematic non-smooth foregrounds like polarized synchrotron emission.

preprint2015arXiv

Euclid & SKA Synergies

Over the past few years two of the largest and highest fidelity experiments conceived have been approved for construction: Euclid is an ESA M-Class mission that will map three-quarters of the extra galactic sky with Hubble Space Telescope resolution optical and NIR imaging, and NIR spectroscopy, its scientific aims (amongst others) are to create a map of the dark Universe and to determine the nature of dark energy. The Square Kilometre Array (SKA) has similar scientific aims (and others) using radio wavelength observations. The two experiments are synergistic in several respects, both through the scientific objectives and through the control of systematic effects. SKA Phase-1 and Euclid will be commissioned on similar timescales offering an exciting opportunity to exploit synergies between these facilities.

preprint2015arXiv

Foreground Subtraction in Intensity Mapping with the SKA

21cm intensity mapping experiments aim to observe the diffuse neutral hydrogen (HI) distribution on large scales which traces the Cosmic structure. The Square Kilometre Array (SKA) will have the capacity to measure the 21cm signal over a large fraction of the sky. However, the redshifted 21cm signal in the respective frequencies is faint compared to the Galactic foregrounds produced by synchrotron and free-free electron emission. In this article, we review selected foreground subtraction methods suggested to effectively separate the 21cm signal from the foregrounds with intensity mapping simulations or data. We simulate an intensity mapping experiment feasible with SKA phase 1 including extragalactic and Galactic foregrounds. We give an example of the residuals of the foreground subtraction with a independent component analysis and show that the angular power spectrum is recovered within the statistical errors on most scales. Additionally, the scale of the Baryon Acoustic Oscillations is shown to be unaffected by foreground subtraction.

preprint2015arXiv

HI galaxy simulations for the SKA: number counts and bias

This chapter describes the assumed specifications and sensitivities for HI galaxy surveys with SKA1 and SKA2. It addresses the expected galaxy number densities based on available simulations as well as the clustering bias over the underlying dark matter. It is shown that a SKA1 HI galaxy survey should be able to find around $5\times 10^6$ galaxies over 5,000 deg$^2$ (up to $z\sim 0.8$), while SKA2 should find $\sim 10^9$ galaxies over 30,000 deg$^2$ (up to $z\sim 2.5$). The numbers presented here have been used throughout the cosmology chapters for forecasting.

preprint2015arXiv

Late-time cosmology with 21cm intensity mapping experiments

We present a framework for forecasting cosmological constraints from future neutral hydrogen intensity mapping experiments at low to intermediate redshifts. In the process, we establish a simple way of comparing such surveys with optical galaxy redshift surveys. We explore a wide range of experimental configurations and assess how well a number of cosmological observables (the expansion rate, growth rate, and angular diameter distance) and parameters (the densities of dark energy and dark matter, spatial curvature, the dark energy equation of state, etc.) will be measured by an extensive roster of upcoming experiments. A number of potential contaminants and systematic effects are also studied in detail. The overall picture is encouraging -- if autocorrelation calibration can be controlled to a sufficient level, Phase I of the SKA should be able to constrain the dark energy equation of state about as well as a DETF Stage IV galaxy redshift survey like Euclid, in roughly the same timeframe.

preprint2015arXiv

Measuring baryon acoustic oscillations with future SKA surveys

The imprint of baryon acoustic oscillations (BAO) in large-scale structure can be used as a standard ruler for mapping out the cosmic expansion history, and hence for testing cosmological models. In this article we briefly describe the scientific background to the BAO technique, and forecast the potential of the Phase 1 and 2 SKA telescopes to perform BAO surveys using both galaxy catalogues and intensity mapping, assessing their competitiveness with current and future optical galaxy surveys. We find that a 25,000 sq. deg. intensity mapping survey on a Phase 1 array will preferentially constrain the radial BAO, providing a highly competitive 2% constraint on the expansion rate at z ~ 2. A 30,000 sq. deg. galaxy redshift survey on SKA2 will outperform all other planned experiments for z < 1.4.

preprint2015arXiv

Measuring redshift-space distortions with future SKA surveys

The peculiar motion of galaxies can be a particularly sensitive probe of gravitational collapse. As such, it can be used to measure the dynamics of dark matter and dark energy as well the nature of the gravitational laws at play on cosmological scales. Peculiar motions manifest themselves as an overall anisotropy in the measured clustering signal as a function of the angle to the line-of-sight, known as redshift-space distortion (RSD). Limiting factors in this measurement include our ability to model non-linear galaxy motions on small scales and the complexities of galaxy bias. The anisotropy in the measured clustering pattern in redshift-space is also driven by the unknown distance factors at the redshift in question, the Alcock-Paczynski distortion. This weakens growth rate measurements, but permits an extra geometric probe of the Hubble expansion rate. In this chapter we will briefly describe the scientific background to the RSD technique, and forecast the potential of the SKA phase 1 and the SKA2 to measure the growth rate using both galaxy catalogues and intensity mapping, assessing their competitiveness with current and future optical galaxy surveys.

preprint2015arXiv

Observational signatures of modified gravity on ultra-large scales

Extremely large surveys with future experiments like Euclid and the SKA will soon allow us to access perturbation modes close to the Hubble scale, with wavenumbers $k \sim \mathcal{H}$. If a modified gravity theory is responsible for cosmic acceleration, the Hubble scale is a natural regime for deviations from General Relativity (GR) to become manifest. The majority of studies to date have concentrated on the consequences of alternative gravity theories for the subhorizon, quasi-static regime, however. In this paper we investigate how modifications to the gravitational field equations affect perturbations around the Hubble scale. We choose functional forms to represent the generic scale-dependent behaviour of gravity theories that modify GR at long wavelengths, and study the resulting deviations of ultra large-scale relativistic observables from their GR behaviour. We find that these are small unless modifications to the field equations are drastic. The angular dependence and redshift evolution of the deviations is highly parameterisation- and survey-dependent, however, and so they are possibly a rich source of modified gravity phenomenology if they can be measured.

preprint2015arXiv

Synergy between the Large Synoptic Survey Telescope and the Square Kilometre Array

We provide an overview of the science benefits of combining information from the Square Kilometre Array (SKA) and the Large Synoptic Survey Telescope (LSST). We first summarise the capabilities and timeline of the LSST and overview its science goals. We then discuss the science questions in common between the two projects, and how they can be best addressed by combining the data from both telescopes. We describe how weak gravitational lensing and galaxy clustering studies with LSST and SKA can provide improved constraints on the causes of the cosmological acceleration. We summarise the benefits to galaxy evolution studies of combining deep optical multi-band imaging with radio observations. Finally, we discuss the excellent match between one of the most unique features of the LSST, its temporal cadence in the optical waveband, and the time resolution of the SKA.

preprint2015arXiv

Ultra large-scale cosmology in next-generation experiments with single tracers

Future surveys of large-scale structure will be able to measure perturbations on the scale of the cosmological horizon, and so could potentially probe a number of novel relativistic effects that are negligibly small on sub-horizon scales. These effects leave distinctive signatures in the power spectra of clustering observables and, if measurable, would open a new window on relativistic cosmology. We quantify the size and detectability of the effects for the most relevant future large-scale structure experiments: spectroscopic and photometric galaxy redshift surveys, intensity mapping surveys of neutral hydrogen, and radio continuum surveys. Our forecasts show that next-generation experiments, reaching out to redshifts $z\simeq 4$, will not be able to detect previously-undetected general-relativistic effects by using individual tracers of the density field, although the contribution of weak lensing magnification on large scales should be clearly detectable. We also perform a rigorous joint forecast for the detection of primordial non-Gaussianity through the excess power it produces in the clustering of biased tracers on large scales, finding that uncertainties of $σ(f_{\rm NL})\sim 1-2$ should be achievable. We study the level of degeneracy of these large-scale effects with several tracer-dependent nuisance parameters, quantifying the minimal priors on the latter that are needed for an optimal measurement of the former. Finally, we discuss the systematic effects that must be mitigated to achieve this level of sensitivity, and some alternative approaches that should help to improve the constraints. The computational tools developed to carry out this study, which requires the full-sky computation of the theoretical angular power spectra for $\mathcal{O}(100)$ redshift bins, as well as realistic models of the luminosity function, are publicly available.

preprint2015arXiv

Weighing neutrinos with cosmic neutral hydrogen

We investigate the signatures left by massive neutrinos on the spatial distribution of neutral hydrogen (HI) in the post-reionization era by running hydrodynamic simulations that include massive neutrinos as additional collisionless particles. We find that halos in massive/massless neutrino cosmologies host a similar amount of neutral hydrogen, although for a fixed halo mass, on average, the HI mass increases with the sum of the neutrino masses. Our results show that HI is more strongly clustered in cosmologies with massive neutrinos, while its abundance, $Ω_{\rm HI}(z)$, is lower. These effects arise mainly from the impact of massive neutrinos on cosmology: they suppress both the amplitude of the matter power spectrum on small scales and the abundance of dark matter halos. Modelling the HI distribution with hydrodynamic simulations at $z > 3$, and a simple analytic model at $z<3$, we use the Fisher matrix formalism to conservatively forecast the constraints that Phase 1 of the Square Kilometre Array (SKA) will place on the sum of neutrino masses, $M_ν\equiv Σm_ν$. We find that with 10,000 hours of interferometric observations at $3 \lesssim z \lesssim 6$ from a deep and narrow survey with SKA1-LOW, the sum of the neutrino masses can be measured with an error $σ(M_ν)\lesssim0.3$ eV (95% CL). Similar constraints can be obtained with a wide and deep SKA1-MID survey at $z \lesssim 3$, using the single-dish mode. By combining data from MID, LOW, and Planck, plus priors on cosmological parameters from a Stage IV spectroscopic galaxy survey, the sum of the neutrino masses can be determined with an error $σ(M_ν)\simeq0.06$ eV (95% CL).

preprint2014arXiv

Blind foreground subtraction for intensity mapping experiments

We make use of a large set of fast simulations of an intensity mapping experiment with characteristics similar to those expected of the Square Kilometre Array (SKA) in order to study the viability and limits of blind foreground subtraction techniques. In particular, we consider different approaches: polynomial fitting, principal component analysis (PCA) and independent component analysis (ICA). We review the motivations and algorithms for the three methods, and show that they can all be described, using the same mathematical framework, as different approaches to the blind source separation problem. We study the efficiency of foreground subtraction both in the angular and radial (frequency) directions, as well as the dependence of this efficiency on different instrumental and modelling parameters. For well-behaved foregrounds and instrumental effects we find that foreground subtraction can be successful to a reasonable level on most scales of interest. We also quantify the effect that the cleaning has on the recovered signal and power spectra. Interestingly, we find that the three methods yield quantitatively similar results, with PCA and ICA being almost equivalent.

preprint2014arXiv

Quintessence in a quandary: prior dependence in dark energy models

The archetypal theory of dark energy is quintessence: a minimally coupled scalar field with a canonical kinetic energy and potential. By studying random potentials we show that quintessence imposes a restricted set of priors on the equation of state of dark energy. Focusing on the commonly-used parametrisation, $w(a)\approx w_0+w_a(1-a)$, we show that there is a natural scale and direction in the $(w_0, w_a)$ plane that distinguishes quintessence as a general framework. We calculate the expected information gain for a given survey and show that, because of the non-trivial prior information, it is a function of more than just the figure of merit. This allows us to make a quantitative case for novel survey strategies. We show that the scale of the prior sets target observational requirements for gaining significant information. This corresponds to a figure of merit FOM$\gtrsim 200$, a requirement that future galaxy redshift surveys will meet.

preprint2013arXiv

What if Planck's Universe isn't flat?

Inflationary theory predicts that the observable Universe should be very close to flat, with a spatial-curvature parameter |Omega_K| < 10^-4. The WMAP satellite currently constrains |Omega_K| < 0.01, and the Planck satellite will be sensitive to values near 10^-3. Suppose that Planck were to find Omega_K to be non-zero at this level. Would this necessarily be a serious problem for inflation? We argue that an apparent departure from flatness could be due either to a local (wavelength comparable to the observable horizon) inhomogeneity, or a truly superhorizon departure from flatness. If there is a local inhomogeneity, then secondary CMB anisotropies distort the CMB frequency spectrum at a level potentially detectable by a next-generation experiment. We discuss how these spectral distortions would complement constraints on the Grishchuk-Zel'dovich effect from the CMB power spectrum at large angular scales in discovering the source of the departure from flatness.

preprint2012arXiv

Local and non-local measures of acceleration in cosmology

Current cosmological observations, when interpreted within the framework of a homogeneous and isotropic Friedmann-Lemaitre-Robertson-Walker (FLRW) model, strongly suggest that the Universe is entering a period of accelerating expansion. This is often taken to mean that the expansion of space itself is accelerating. In a general spacetime, however, this is not necessarily true. We attempt to clarify this point by considering a handful of local and non-local measures of acceleration in a variety of inhomogeneous cosmological models. Each of the chosen measures corresponds to a theoretical or observational procedure that has previously been used to study acceleration in cosmology, and all measures reduce to the same quantity in the limit of exact spatial homogeneity and isotropy. In statistically homogeneous and isotropic spacetimes, we find that the acceleration inferred from observations of the distance-redshift relation is closely related to the acceleration of the spatially averaged universe, but does not necessarily bear any resemblance to the average of the local acceleration of spacetime itself. For inhomogeneous spacetimes that do not display statistical homogeneity and isotropy, however, we find little correlation between acceleration inferred from observations and the acceleration of the averaged spacetime. This shows that observations made in an inhomogeneous universe can imply acceleration without the existence of dark energy.

preprint2012arXiv

The isotropic blackbody CMB as evidence for a homogeneous universe

The question of whether the Universe is spatially homogeneous and isotropic on the largest scales is of fundamental importance to cosmology, but has not yet been answered decisively. Surprisingly, neither an isotropic primary CMB nor combined observations of luminosity distances and galaxy number counts are sufficient to establish such a result. The inclusion of the Sunyaev-Zel'dovich effect in CMB observations, however, dramatically improves this situation. We show that even a solitary observer who sees an isotropic blackbody CMB can conclude that the universe is homogeneous and isotropic in their causal past when the Sunyaev-Zel'dovich effect is present. Critically, however, the CMB must either be viewed for an extended period of time, or CMB photons that have scattered more than once must be detected. This result provides a theoretical underpinning for testing the Cosmological Principle with observations of the CMB alone.

preprint2012arXiv

The kSZ effect as a test of general radial inhomogeneity in LTB cosmology

The apparent accelerating expansion of the Universe, determined from observations of distant supernovae, and often taken to imply the existence of dark energy, may alternatively be explained by the effects of a giant underdense void if we relax the assumption of homogeneity on large scales. Recent studies have made use of the spherically-symmetric, radially-inhomogeneous Lemaitre-Tolman-Bondi (LTB) models to derive strong constraints on this scenario, particularly from observations of the kinematic Sunyaev-Zel'dovich (kSZ) effect which is sensitive to large scale inhomogeneity. However, most of these previous studies explicitly set the LTB 'bang time' function to be constant, neglecting an important freedom of the general solutions. Here we examine these models in full generality by relaxing this assumption. We find that although the extra freedom allowed by varying the bang time is sufficient to account for some observables individually, it is not enough to simultaneously explain the supernovae observations, the small-angle CMB, the local Hubble rate, and the kSZ effect. This set of observables is strongly constraining, and effectively rules out simple LTB models as an explanation of dark energy.