Researcher profile

Steven Murray

Steven Murray contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 19 - Baseline
5works
0followers
2topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

5 published item(s)

preprint2022arXiv

Machine learning astrophysics from 21 cm lightcones: impact of network architectures and signal contamination

Imaging the cosmic 21 cm signal will map out the first billion years of our Universe. The resulting 3D lightcone (LC) will encode the properties of the unseen first galaxies and physical cosmology. Here, we build on previous work using neural networks (NNs) to infer astrophysical parameters directly from 21 cm LC images. We introduce recurrent neural networks (RNNs), capable of efficiently characterizing the evolution along the redshift axis of 21 cm LC images. Using a large database of simulated cosmic 21 cm LCs, we compare the relative performance in parameter estimation of different network architectures. These including two types of RNNs, which differ in their complexity, as well as a more traditional convolutional neural network (CNN). For the ideal case of no instrumental effects, our simplest and easiest to train RNN performs the best, with a mean squared parameter estimation error (MSE) that is lower by a factor of $\ge 2$ compared with the other architectures studied here, and a factor of $\ge 8$ lower than the previously-studied CNN. We also corrupt the cosmic signal by adding noise expected from a 1000 h integration with the Square Kilometre Array, as well as excising a foreground-contaminated 'horizon wedge'. Parameter prediction errors increase when the NNs are trained on these contaminated LC images, though recovery is still good even in the most pessimistic case (with $R^2 \ge 0.5-0.95$). However, we find no notable differences in performance between network architectures on the contaminated images. We argue this is due to the size of our data set, highlighting the need for larger data sets and/or better data augmentation in order to maximize the potential of NNs in 21 cm parameter estimation.

preprint2020arXiv

The Impact of Realistic Foreground and Instrument Models on 21cm Epoch of Reionization Experiments

Predictions for the ability of 21-cm interferometric experiments to discriminate Epoch of Reionization (EoR) signal models are typically limited by the simplicity of data models, whereby foreground signals and characteristics of the instrument are often simplified or neglected.To move towards more realistic scenarios, we explore the effects of applying more realistic foreground and instrument models to the 21cm signal, and the ability to estimate astrophysical parameters with these additional complexities. We use a highly-optimized version of \textsc{21cmFAST}, integrated into \textsc{21cmMC}, to generate lightcones of the brightness temperature fluctuation for Bayesian parameter estimation. We include a statistical point-source foreground model and an instrument model based on the Murchison Widefield Array (MWA) scaled in observation time to have an effective sensitivity similar to the future Square Kilometre Array (SKA). We also extend the current likelihood prescription to account for the presence of beam convolution and foregrounds, the 2-Dimensional Power Spectrum (PS), and the correlation of PS modes. We use frequency bands between 150 and 180 MHz to constrain the ionizing efficiency ($ζ$), the minimum virial temperature of halos ($T_{\mathrm{vir}}$), the soft X-ray emissivity per unit Star Formation Rate (SFR) ($L_X/SFR$ ), and the X-ray energy threshold ($E_0$). We find that the inclusion of realistic foregrounds and instrumental components biases the parameter constraints due to unaccounted for cross-power between the EoR signal, foregrounds and thermal noise. This causes estimates of $ζ$ to be biased by up to $5σ$ but the estimates of $T_{vir}$, L$_X$/SFR and E$_0$ remain unaffected and are all within $1σ$.

preprint2015arXiv

nIFTy Cosmology: Galaxy/halo mock catalogue comparison project on clustering statistics

We present a comparison of major methodologies of fast generating mock halo or galaxy catalogues. The comparison is done for two-point and the three-point clustering statistics. The reference catalogues are drawn from the BigMultiDark N-body simulation. Both friend-of-friends (including distinct halos only) and spherical overdensity (including distinct halos and subhalos) catalogs have been used with the typical number density of a large-volume galaxy surveys. We demonstrate that a proper biasing model is essential for reproducing the power spectrum at quasilinear and even smaller scales. With respect to various clustering statistics a methodology based on perturbation theory and a realistic biasing model leads to very good agreement with N-body simulations. However, for the quadrupole of the correlation function or the power spectrum, only the method based on semi-N-body simulation could reach high accuracy (1% level) at small scales, i.e., r<25 Mpc/h or k>0.15 h/Mpc. Full N-body solutions will remain indispensable to produce reference catalogues. Nevertheless, we have demonstrated that the far more efficient approximate solvers can reach a few percent accuracy in terms of clustering statistics at the scales interesting for the large-scale structure analysis after calibration with a few reference N-body calculations. This makes them useful for massive production aimed at covariance studies, to scan large parameter spaces, and to estimate uncertainties in data analysis techniques, such as baryon acoustic oscillation reconstruction, redshift distortion measurements, etc.

preprint2013arXiv

HMFcalc: An Online Tool for Calculating Dark Matter Halo Mass Functions

The dark matter halo mass function (HMF) is a characteristic property of cosmological structure formation models, quantifying the number density of dark matter haloes per unit mass in the Universe. A key goal of current and planned large galaxy surveys is to measure the HMF and to use it to test theories of dark matter and dark energy. We present a new web application for calculating the HMF -- the frontend HMFcalc and the engine hmf. HMFcalc has been designed to be flexible, efficient and easy to use, providing observational and theoretical astronomers alike with the means to explore standard functional forms of the HMF or to tailor their own. We outline the theoretical background needed to compute the HMF, we show how it has been implemented in hmf, and finally we provide worked examples that illustrate HMFcalc&#39;s versatility as an analysis tool.

preprint2013arXiv

How well do we know the Halo Mass Function?

The parameters governing the standard ΛCold Dark Matter cosmological model have been constrained with unprecedented accuracy by precise measurements of the cosmic microwave background by the Wilkinson Microwave Anisotropy Probe (WMAP) and Planck satellites. Each new data release has refined further our knowledge of quantities -- such as the matter density parameter Ω_M -- that are imprinted on the dark matter halo mass function (HMF), a powerful probe of dark matter and dark energy models. In this letter we trace how changes in the cosmological parameters over the last decade have influenced uncertainty in our knowledge of the HMF. We show that this uncertainty has reduced significantly since the 3rd WMAP data release, but the rate of this reduction is slowing. This is limited by uncertainty in the normalisation σ_8, whose influence is most pronounced at the high mass end of the mass function. Interestingly, we find that the accuracy with which we can constrain the HMF in terms of the cosmological parameters has now reached the point at which it is comparable to the scatter in HMF fitting functions. This suggests that the power of the HMF as a precision probe of dark matter and dark energy hinges on more accurate determination of the theoretical HMF. Finally, we assess prospects of using the HMF to differentiate between Cold and Warm Dark Matter models based on ongoing improvements in measurements of Ω_M, and we comment briefly on optimal survey strategies for constraining dark matter and dark energy models using the HMF.