BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

W. J. Handley

W. J. Handley contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.COastro-ph.IMgr-qcphysics.comp-phastro-ph.GAastro-ph.SRmath.NANumerical Analysis

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

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

10 published item(s)

preprint2022arXiv

A Comprehensive Bayesian re-analysis of the SARAS2 data from the Epoch of Reionization

We present a Bayesian re-analysis of the sky-averaged 21-cm experimental data from SARAS2 using nested sampling implemented with polychord, spectrally smooth foreground modelling implemented with maxsmooth, detailed systematic modelling and rapid signal emulation with globalemu. Our analysis differs from previous analysis of the SARAS2 data through the use of a full Bayesian framework and separate modelling of the foreground and non-smooth systematics. We use the most up-to-date signal models including Lyman-$α$ and CMB heating parameterised by astrophysical parameters such as star formation efficiency, X-ray heating efficiency, minimal virial circular velocity of star forming galaxies, CMB optical depth and the low energy cutoff of the X-ray spectral energy distribution. We consider models with an excess radio background above the CMB produced via radio emission from early galaxies and parameterised by a radio production efficiency. A non-smooth systematic is identified and modelled as both a frequency damped sinusoid introduced by the electronics and separately from the sky. The latter is modulated by the total efficiency of the antenna and marginally favoured by the data. We consider three different models for the noise in the data. The SARAS2 constraints on individual astrophysical parameters are extremely weak however we identify classes of disfavoured signals. We weakly disfavour standard astrophysical models with high Lyman-$α$ fluxes and weak heating and more confidently disfavour exotic models with high Lyman-$α$ fluxes, low X-ray efficiencies and high radio production efficiencies in early galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Impact of the Primordial Stellar Initial Mass Function on the 21-cm Signal

Properties of the first generation of stars (Pop III), such as their initial mass function (IMF), are poorly constrained by observations and have yet to converge between simulations. The cosmological 21-cm signal of neutral hydrogen is predicted to be sensitive to Lyman-band photons produced by these stars, thus providing a unique way to probe the first stellar population. In this paper, we investigate the impacts of the Pop III IMF on the cosmic dawn 21-cm signal via the Wouthuysen-Field effect, Lyman-Werner feedback, Ly-alpha heating, and CMB heating. We calculate the emission spectra of star-forming halos for different IMFs by integrating over individual metal-free stellar spectra, computed from a set of stellar evolution histories and stellar atmospheres, and taking into account variability of the spectra with stellar age. Through this study, we therefore relax two common assumptions: that the zero-age main sequence emission rate of a Pop III star is representative of its lifetime mean emission rate, and that Pop III emission can be treated as instantaneous. Exploring a bottom-heavy, a top-heavy, and intermediate IMFs, we show that variations in the 21-cm signal are driven by stars lighter than 20 solar masses. For the explored models we find maximum relative differences of 59% in the cosmic dawn global 21-cm signal, and 131% between power spectra. Although this impact is modest, precise modelling of the first stars and their evolution is necessary for accurate prediction and interpretation of the 21-cm signal.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Improved cosmological fits with quantized primordial power spectra

We observationally examine cosmological models based on primordial power spectra with quantized wavevectors. Introducing a linearly quantized power spectrum with $k_0=3.225\times10^{-4}\mathrm{Mpc}^{-1}$ and spacing $Δk = 2.257 \times 10^{-4} \mathrm{Mpc}^{-1}$ provides a better fit to the Planck 2018 observations than the concordance baseline, with $Δχ^2 = -8.55$. Extending the results of Lasenby et al [1], we show that the requirement for perturbations to remain finite beyond the future conformal boundary in a universe containing dark matter and a cosmological constant results in a linearly quantized primordial power spectrum. It is found that the infrared cutoffs for this future conformal boundary quantized cosmology do not provide cosmic microwave background power spectra compatible with observations, but future theories may predict more observationally consistent quantized spectra.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Perturbations and the Future Conformal Boundary

The concordance model of cosmology predicts a universe which finishes in a finite amount of conformal time at a future conformal boundary. We show that for particular cases we study, the background variables and perturbations may be analytically continued beyond this boundary and that the "end of the universe" is not necessarily the end of their physical development. Remarkably, these theoretical considerations of the end of the universe might have observable consequences today: perturbation modes consistent with these boundary conditions have a quantised power spectrum which may be relevant to features seen in the large scale cosmic microwave background. Mathematically these cosmological models may either be interpreted as a palindromic universe mirrored in time, a reflecting boundary condition, or a double cover, but are identical with respect to their observational predictions and stand in contrast to the predictions of conformal cyclic cosmologies.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Bayesian noise wave calibration for 21-cm global experiments

Detection of millikelvin-level signals from the 'Cosmic Dawn' requires an unprecedented level of sensitivity and systematic calibration. We report the theory behind a novel calibration algorithm developed from the formalism introduced by the EDGES collaboration for use in 21-cm experiments. Improvements over previous approaches are provided through the incorporation of a Bayesian framework and machine learning techniques such as the use of Bayesian evidence to determine the level of frequency variation of calibration parameters that is supported by the data, the consideration of correlation between calibration parameters when determining their values and the use of a conjugate-prior based approach that results in a fast algorithm for application in the field. In self-consistency tests using empirical data models of varying complexity, our methodology is used to calibrate a 50 $Ω$ ambient-temperature load. The RMS error between the calibration solution and the measured temperature of the load is 8 mK, well within the 1$σ$ noise level. Whilst the methods described here are more applicable to global 21-cm experiments, they can easily be adapted and applied to other applications, including telescopes such as HERA and the SKA.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

maxsmooth: Rapid Maximally Smooth Function Fitting With Applications in Global 21-cm Cosmology

Maximally Smooth Functions (MSFs) are a form of constrained functions in which there are no inflection points or zero crossings in high order derivatives. Consequently, they have applications to signal recovery in experiments where signals of interest are expected to be non-smooth features masked by larger smooth signals or foregrounds. They can also act as a powerful tool for diagnosing the presence of systematics. The constrained nature of MSFs makes fitting these functions a non-trivial task. We introduce maxsmooth, an open source package that uses quadratic programming to rapidly fit MSFs. We demonstrate the efficiency and reliability of maxsmooth by comparison to commonly used fitting routines and show that we can reduce the fitting time by approximately two orders of magnitude. We introduce and implement with maxsmooth Partially Smooth Functions, which are useful for describing elements of non-smooth structure in foregrounds. This work has been motivated by the problem of foreground modelling in 21-cm cosmology. We discuss applications of maxsmooth to 21-cm cosmology and highlight this with examples using data from the Experiment to Detect the Global Epoch of Reionization Signature (EDGES) and the Large-aperture Experiment to Detect the Dark Ages (LEDA) experiments. We demonstrate the presence of a sinusoidal systematic in the EDGES data with a log-evidence difference of $86.19\pm0.12$ when compared to a pure foreground fit. MSFs are applied to data from LEDA for the first time in this paper and we identify the presence of sinusoidal systematics. maxsmooth is pip installable and available for download at: https://github.com/htjb/maxsmooth

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Addressing $H_0$ tension with emergent dark radiation in unitary gravity

We propose a one-parameter extension to $Λ$CDM, expected to strongly affect cosmological tensions. An effective dark radiation component in the early universe redshifts away as hot dark matter, then quintessence, tracking the dominant equation-of-state parameter and leaving a falsifiable torsion field in the current epoch. This picture results from a new Poincaré gauge theory (PGT), one of the most promising among the latest batch of 58 PGTs found to be both power-counting renormalisable and free from ghosts and tachyons. We systematically categorise the cosmologies of 33 of these PGTs, as special cases of the most general parity-preserving, Ostrogradsky-stable PGT with a purely Yang-Mills action. The theory we consider contains two propagating massless gravitons, which may be $J^P=2^+$ (long-range gravitation and gravitational waves). A conspiracy among the coupling constants eliminates the spatial curvature $k\in\{\pm 1,0\}$ from the field equations. We show that this `$k$-screening' is not restricted to conformal gravity theories. The flat Friedmann equations are then emergent, with potentially tension-resolving freedom at the early scale-invariant epoch that reliably gives way to an attractor-like state of modern $Λ$CDM evolution. We compare with related theories and promising special cases, such as $k$-screened theories with negative-definite effective $k$, and more traditional theories with effective $Λ$ and a $J^P=0^-$ massive graviton (dark matter candidate). As a bonus, we analyse similarly constrained actions in the new extended Weyl gauge theory (eWGT). We show that in cosmology, PGT and eWGT span exactly the same classical phenomenology up to a linear map between their coupling constants, hinting at a deeper relationship between the two.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Dense output for highly oscillatory numerical solutions

We present a method to construct a continuous extension (otherwise known as dense output) for a numerical routine in the special case of the numerical solution being a scalar-valued function exhibiting rapid oscillations. Such cases call for numerical routines that make use of the known global behaviour of the solution, one example being methods using asymptotic expansions to forecast the solution at each step of the independent variable. An example is oscode, numerical routine which uses the Wentzel-Kramers-Brillouin (WKB) approximation when the solution oscillates rapidly and otherwise behaves as a Runge-Kutta (RK) solver. Polynomial interpolation is not suitable for producing the solution at an arbitrary point mid-step, since efficient numerical methods based on the WKB approximation will step through multiple oscillations in a single step. Instead we construct the continuous solution by extending the numerical quadrature used in computing a WKB approximation of the solution with no additional evaluations of the differential equation or terms within, and provide an error estimate on this dense output. Finally, we draw attention to previous work on the continuous extension of Runge-Kutta formulae, and construct an extension to a RK method based on Gauss--Lobatto quadrature nodes, thus describing how to generate dense output from each of the methods underlying oscode.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Quantum initial conditions for inflation and canonical invariance

We investigate the transformation of initial conditions for primordial curvature perturbations under two types of transformations of the associated action: simultaneous redefinition of time and the field to be quantised, and the addition of surface terms. The latter encompasses all canonical transformations, whilst the time- and field-redefinition is a distinct, non-canonical transformation since the initial and destination systems use different times. Actions related to each other via such transformations yield identical equations of motion and preserve the commutator structure. They further preserve the time-evolution of expectation values of quantum operators unless the vacuum state also changes under the transformation. These properties suggest that it is of interest to investigate vacuum prescriptions that also remain unchanged under canonical transformations. We find that initial conditions derived via minimising the vacuum expectation value of the Hamiltonian and those obtained using the Danielsson vacuum prescription are not invariant under these transformations, whereas those obtained by minimising the local energy density are. We derive the range of physically distinct initial conditions obtainable by Hamiltonian diagonalisation, and illustrate their effect on the scalar primordial power spectrum and the Cosmic Microwave Background under the just enough inflation model. We also generalise the analogy between the dynamics of a quantum scalar field on a curved, time-dependent spacetime and the gauge-invariant curvature perturbation. We argue that the invariance of the vacuum prescription obtained by minimising the renormalised stress--energy tensor should make it the preferred procedure for setting initial conditions for primordial perturbations. All other procedures reviewed in this work yield ambiguous initial conditions, which is problematic both in theory and practice.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

An efficient method for solving highly oscillatory ordinary differential equations with applications to physical systems

We present a novel numerical routine (oscode) with a C++ and Python interface for the efficient solution of one-dimensional, second-order, ordinary differential equations with rapidly oscillating solutions. The method is based on a Runge-Kutta-like stepping procedure that makes use of the Wentzel-Kramers-Brillouin (WKB) approximation to skip regions of integration where the characteristic frequency varies slowly. In regions where this is not the case, the method is able to switch to a made-to-measure Runge-Kutta integrator that minimises the total number of function evaluations. We demonstrate the effectiveness of the method with example solutions of the Airy equation and an equation exhibiting a burst of oscillations, discussing the error properties of the method in detail. We then show the method applied to physical systems. First, the one-dimensional, time-independent Schrödinger equation is solved as part of a shooting method to search for the energy eigenvalues for a potential with quartic anharmonicity. Then, the method is used to solve the Mukhanov-Sasaki equation describing the evolution of cosmological perturbations, and the primordial power spectrum of the perturbations is computed in different cosmological scenarios. We compare the performance of our solver in calculating a primordial power spectrum of scalar perturbations to that of BINGO, an efficient code specifically designed for such applications.

0 citations0 reviewsNo code linkedOpen access