-
Resonant or asymmetric: The status of sub-GeV dark matter
Authors:
Sowmiya Balan,
Csaba Balázs,
Torsten Bringmann,
Christopher Cappiello,
Riccardo Catena,
Timon Emken,
Tomás E. Gonzalo,
Taylor R. Gray,
Will Handley,
Quan Huynh,
Felix Kahlhoefer,
Aaron C. Vincent
Abstract:
Sub-GeV dark matter (DM) particles produced via thermal freeze-out evade many of the strong constraints on heavier DM candidates but at the same time face a multitude of new constraints from laboratory experiments, astrophysical observations and cosmological data. In this work we combine all of these constraints in order to perform frequentist and Bayesian global analyses of fermionic and scalar s…
▽ More
Sub-GeV dark matter (DM) particles produced via thermal freeze-out evade many of the strong constraints on heavier DM candidates but at the same time face a multitude of new constraints from laboratory experiments, astrophysical observations and cosmological data. In this work we combine all of these constraints in order to perform frequentist and Bayesian global analyses of fermionic and scalar sub-GeV DM coupled to a dark photon with kinetic mixing. For fermionic DM, we find viable parameter regions close to the dark photon resonance, which expand significantly when including a particle-antiparticle asymmetry. For scalar DM, the velocity-dependent annihilation cross section evades the strongest constraints even in the symmetric case. Using Bayesian model comparison, we show that both asymmetric fermionic DM and symmetric scalar DM are preferred over symmetric fermionic DM due to the reduced fine-tuning penalty. Finally, we explore the discovery prospects of near-future experiments both in the full parameter space and for specific benchmark points. We find that the most commonly used benchmark scenarios are already in tension with existing constraints and propose a new benchmark point that can be targeted with future searches.
△ Less
Submitted 27 May, 2024;
originally announced May 2024.
-
Costless correction of chain based nested sampling parameter estimation in gravitational wave data and beyond
Authors:
Metha Prathaban,
Will Handley
Abstract:
Nested sampling parameter estimation differs from evidence estimation, in that it incurs an additional source of error. This error affects estimates of parameter means and credible intervals in gravitational wave analyses and beyond, and yet, it is typically not accounted for in standard error estimation methods. In this paper, we present two novel methods to quantify this error more accurately fo…
▽ More
Nested sampling parameter estimation differs from evidence estimation, in that it incurs an additional source of error. This error affects estimates of parameter means and credible intervals in gravitational wave analyses and beyond, and yet, it is typically not accounted for in standard error estimation methods. In this paper, we present two novel methods to quantify this error more accurately for any chain based nested sampler, using the additional likelihood calls made at runtime in producing independent samples. Using injected signals of black hole binary coalescences as an example, we first show concretely that the usual error estimation method is insufficient to capture the true error bar on parameter estimates. We then demonstrate how the extra points in the chains of chain based samplers may be carefully utilised to estimate this error correctly, and provide a way to check the accuracy of the resulting error bars. Finally, we discuss how this error affects $p$-$p$ plots and coverage assessments.
△ Less
Submitted 25 April, 2024;
originally announced April 2024.
-
Simulation-Based Inference of the sky-averaged 21-cm signal from CD-EoR with REACH
Authors:
Anchal Saxena,
P. Daniel Meerburg,
Christoph Weniger,
Eloy de Lera Acedo,
Will Handley
Abstract:
The redshifted 21-cm signal from the Cosmic Dawn and Epoch of Reionization carries invaluable information about the cosmology and astrophysics of the early Universe. Analyzing the data from a sky-averaged 21-cm signal experiment typically involves navigating through an intricate parameter space to accurately address various factors such as foregrounds, beam uncertainties, ionospheric distortions,…
▽ More
The redshifted 21-cm signal from the Cosmic Dawn and Epoch of Reionization carries invaluable information about the cosmology and astrophysics of the early Universe. Analyzing the data from a sky-averaged 21-cm signal experiment typically involves navigating through an intricate parameter space to accurately address various factors such as foregrounds, beam uncertainties, ionospheric distortions, and receiver noise for the search of the cosmological 21-cm signal. The traditional likelihood-based sampling methods for modeling these effects could become computationally demanding for such highly complex models, which makes it infeasible to include physically motivated 21-cm signal models in the analysis. Moreover, the inference with these traditional methods is driven by the assumed functional form of the likelihood function. This work demonstrates how Simulation-Based Inference through Truncated Marginal Neural Ratio Estimation (TMNRE) can naturally handle these issues at a significantly reduced computational cost than the likelihood-based methods. We estimate the posterior distribution on our model parameters with TMNRE for simulated mock observations, composed of beam-weighted foregrounds, physically motivated 21-cm signal, and radiometric noise. We find that maximizing the information content by simultaneously analyzing the data from multiple time slices and antennas significantly improves the parameter constraints and leads to a better exploration of the cosmological signal. We discuss the application of TMNRE for the current configuration of the REACH experiment and demonstrate how it can be utilized to explore potential avenues for REACH. The method presented here can be easily generalized for any sky-averaged 21-cm signal experiment.
△ Less
Submitted 21 March, 2024;
originally announced March 2024.
-
On the Constraints on Superconducting Cosmic Strings from 21-cm Cosmology
Authors:
T. Gessey-Jones,
S. Pochinda,
H. T. J. Bevins,
A. Fialkov,
W. J. Handley,
E. de Lera Acedo,
S. Singh,
R. Barkana
Abstract:
Constraints on the potential properties of superconducting cosmic strings provide an indirect probe of physics beyond the standard model at energies inaccessible to terrestrial particle colliders. In this study, we perform the first joint Bayesian analysis to extract constraints on superconducting cosmic strings from current 21-cm signal measurements while accounting rigorously for the uncertainti…
▽ More
Constraints on the potential properties of superconducting cosmic strings provide an indirect probe of physics beyond the standard model at energies inaccessible to terrestrial particle colliders. In this study, we perform the first joint Bayesian analysis to extract constraints on superconducting cosmic strings from current 21-cm signal measurements while accounting rigorously for the uncertainties in foregrounds and high redshift astrophysics. We include the latest publicly available 21-cm power spectrum upper limits from HERA, 21-cm global signal data from SARAS 3, and the synergistic probe of the unresolved X-ray background in our final analysis. This paper thus constitutes the first attempt to use 21-cm power spectrum data to probe cosmic strings. In contrast to previous works, we find no strong constraints can be placed on superconducting cosmic strings from current 21-cm measurements. This is because of uncertainties in the X-ray emission efficiency of the first galaxies, with X-ray emissivities greater than $3 \times 10^{40}$erg s$^{-1}$ M$_{\odot}^{-1}$ yr able to mask the presence of cosmic strings in the 21-cm signal. We conclude by discussing the prospects for future constraints from definitive 21-cm signal measurements and argue that the recently proposed soft photon heating should be cause for optimism due to its potential to break degeneracies that would have otherwise made the signatures of cosmic strings difficult to distinguish from those of astrophysical origin.
△ Less
Submitted 15 February, 2024; v1 submitted 14 December, 2023;
originally announced December 2023.
-
Non-Gaussian Likelihoods for Type Ia Supernovae Cosmology: Implications for Dark Energy and $H_0$
Authors:
Toby Lovick,
Suhail Dhawan,
Will Handley
Abstract:
The latest improvements in the scale and calibration of Type Ia supernovae catalogues allow us to constrain the specific nature and evolution of dark energy through its effect on the expansion history of the universe. We present the results of Bayesian cosmological model comparison on the SNe~Ia catalogue Pantheon+, where Flat $Λ$CDM is preferred by the data over all other models and we find moder…
▽ More
The latest improvements in the scale and calibration of Type Ia supernovae catalogues allow us to constrain the specific nature and evolution of dark energy through its effect on the expansion history of the universe. We present the results of Bayesian cosmological model comparison on the SNe~Ia catalogue Pantheon+, where Flat $Λ$CDM is preferred by the data over all other models and we find moderate evidence ($Δ\log \mathcal{Z} \sim 2.5$) to reject a number of the alternate dark energy models. The effect of peculiar velocity corrections on model comparison is analysed, where we show that removing the peculiar velocity corrections results in a varying fit on non-$Λ$CDM parameters. As well as comparing cosmological models, the Bayesian methodology is extended to comparing the scatter model of the data, testing for non-gaussianity in the Pantheon+ Hubble residuals. We find that adding a scale parameter to the Pantheon+ covariances, or alternately using a multivariate Student's t-distribution fits the data better than the fiducial analysis, producing a cosmology independent evidence increase of $Δ\log \mathcal{Z} = 2.29 $ and $2.46$ respectively. This improved treatment of the scatter decreases the uncertainty in the constraint on the Hubble constant, finding $H_0 = 73.67 \pm 0.99 $ km s$^{-1}$ Mpc$^{-1}$, in $ 5.7 σ$ tension with Planck. We also explore $M_B$ transition models as a potential solution for the Hubble tension, finding no evidence to support these models among the SNe data.
△ Less
Submitted 4 December, 2023;
originally announced December 2023.
-
aeons: approximating the end of nested sampling
Authors:
Zixiao Hu,
Artem Baryshnikov,
Will Handley
Abstract:
This paper presents analytic results on the anatomy of nested sampling, from which a technique is developed to estimate the run-time of the algorithm that works for any nested sampling implementation. We test these methods on both toy models and true cosmological nested sampling runs. The method gives an order-of-magnitude prediction of the end point at all times, forecasting the true endpoint wit…
▽ More
This paper presents analytic results on the anatomy of nested sampling, from which a technique is developed to estimate the run-time of the algorithm that works for any nested sampling implementation. We test these methods on both toy models and true cosmological nested sampling runs. The method gives an order-of-magnitude prediction of the end point at all times, forecasting the true endpoint within standard error around the halfway point.
△ Less
Submitted 30 November, 2023;
originally announced December 2023.
-
Kernel-, mean- and noise-marginalised Gaussian processes for exoplanet transits and $H_0$ inference
Authors:
Namu Kroupa,
David Yallup,
Will Handley,
Michael Hobson
Abstract:
Using a fully Bayesian approach, Gaussian Process regression is extended to include marginalisation over the kernel choice and kernel hyperparameters. In addition, Bayesian model comparison via the evidence enables direct kernel comparison. The calculation of the joint posterior was implemented with a transdimensional sampler which simultaneously samples over the discrete kernel choice and their h…
▽ More
Using a fully Bayesian approach, Gaussian Process regression is extended to include marginalisation over the kernel choice and kernel hyperparameters. In addition, Bayesian model comparison via the evidence enables direct kernel comparison. The calculation of the joint posterior was implemented with a transdimensional sampler which simultaneously samples over the discrete kernel choice and their hyperparameters by embedding these in a higher-dimensional space, from which samples are taken using nested sampling. Kernel recovery and mean function inference were explored on synthetic data from exoplanet transit light curve simulations. Subsequently, the method was extended to marginalisation over mean functions and noise models and applied to the inference of the present-day Hubble parameter, $H_0$, from real measurements of the Hubble parameter as a function of redshift, derived from the cosmologically model-independent cosmic chronometer and $Λ$CDM-dependent baryon acoustic oscillation observations. The inferred $H_0$ values from the cosmic chronometers, baryon acoustic oscillations and combined datasets are $H_0= 66 \pm 6\, \mathrm{km}\,\mathrm{s}^{-1}\,\mathrm{Mpc}^{-1}$, $H_0= 67 \pm 10\, \mathrm{km}\,\mathrm{s}^{-1}\,\mathrm{Mpc}^{-1}$ and $H_0= 69 \pm 6\, \mathrm{km}\,\mathrm{s}^{-1}\,\mathrm{Mpc}^{-1}$, respectively. The kernel posterior of the cosmic chronometers dataset prefers a non-stationary linear kernel. Finally, the datasets are shown to be not in tension with $\ln R=12.17\pm 0.02$.
△ Less
Submitted 12 February, 2024; v1 submitted 7 November, 2023;
originally announced November 2023.
-
Balancing ACT: weighing prior dependency and global tensions of DR6 lensing with other datasets
Authors:
A. N. Ormondroyd,
W. J. Handley,
M. P. Hobson,
A. N. Lasenby
Abstract:
We provide a complementary nested sampling analysis for the Atacama Cosmology Telescope lensing data release 6. This allows the quantification of global consistency statistics between ACT lensing and alternative datasets. In the context of flat $Λ$CDM, we find no inconsistency between ACT, Baryonic Acoustic Oscillations, Planck anisotropies, weak lensing datasets, or NPIPE lensing. As part of our…
▽ More
We provide a complementary nested sampling analysis for the Atacama Cosmology Telescope lensing data release 6. This allows the quantification of global consistency statistics between ACT lensing and alternative datasets. In the context of flat $Λ$CDM, we find no inconsistency between ACT, Baryonic Acoustic Oscillations, Planck anisotropies, weak lensing datasets, or NPIPE lensing. As part of our analysis, we also investigate the effect of the prior widths used in the ACT analysis and find that the headline results are quantitatively but not qualitatively affected by the chosen priors. We use both Bayes factors and the suspiciousness statistic to quantify the possibility of tension, and find suspiciousness unsuitable in the case of strong agreement between ACT DR6 and NPIPE. Nested sampling provides a competitive alternative to Metropolis Hastings and we recommend it be used alongside existing analyses. We release the chains and plotting source for the analysis using anesthetic.
△ Less
Submitted 3 December, 2023; v1 submitted 12 October, 2023;
originally announced October 2023.
-
Analytic Approximations for the Primordial Power Spectrum with Israel Junction Conditions
Authors:
D. D. Dineen,
W. J. Handley
Abstract:
This work compares cosmological matching conditions used in approximating generic pre-inflationary phases of the universe. We show that the joining conditions for primordial scalar perturbations assumed by Contaldi et al. are inconsistent with the physically motivated Israel junction conditions, however; performing general relativistic matching with the aforementioned constraints results in unreal…
▽ More
This work compares cosmological matching conditions used in approximating generic pre-inflationary phases of the universe. We show that the joining conditions for primordial scalar perturbations assumed by Contaldi et al. are inconsistent with the physically motivated Israel junction conditions, however; performing general relativistic matching with the aforementioned constraints results in unrealistic primordial power spectra. Eliminating the need for ambiguous matching, we look at an alternative semi-analytic model for producing the primordial power spectrum allowing for finite duration cosmological phase transitions.
△ Less
Submitted 13 May, 2024; v1 submitted 27 September, 2023;
originally announced September 2023.
-
Fully Bayesian Forecasts with Evidence Networks
Authors:
T. Gessey-Jones,
W. J. Handley
Abstract:
Sensitivity forecasts inform the design of experiments and the direction of theoretical efforts. To arrive at representative results, Bayesian forecasts should marginalize their conclusions over uncertain parameters and noise realizations rather than picking fiducial values. However, this is typically computationally infeasible with current methods for forecasts of an experiment's ability to disti…
▽ More
Sensitivity forecasts inform the design of experiments and the direction of theoretical efforts. To arrive at representative results, Bayesian forecasts should marginalize their conclusions over uncertain parameters and noise realizations rather than picking fiducial values. However, this is typically computationally infeasible with current methods for forecasts of an experiment's ability to distinguish between competing models. We thus propose a novel simulation-based methodology capable of providing expedient and rigorous Bayesian model comparison forecasts without relying on restrictive assumptions.
△ Less
Submitted 23 May, 2024; v1 submitted 13 September, 2023;
originally announced September 2023.
-
Receiver design for the REACH global 21-cm signal experiment
Authors:
Nima Razavi-Ghods,
Ian L. V. Roque,
Steven H. Carey,
John A. Ely,
Will Handley,
Alessio Magro,
Riccardo Chiello,
Tian Huang,
P. Alexander,
D. Anstey,
G. Bernardi,
H. T. J. Bevins,
J. Cavillot,
W. Croukamp,
J. Cumner,
E. de Lera Acedo,
D. I. L. de Villiers,
A. Fialkov,
T. Gessey-Jones,
Q. Gueuning,
A. T. Josaitis,
G. Kulkarni,
S. A. K. Leeney,
R. Maiolino,
P. D. Meerburg
, et al. (13 additional authors not shown)
Abstract:
We detail the the REACH radiometric system designed to enable measurements of the 21-cm neutral hydrogen line. Included is the radiometer architecture and end-to-end system simulations as well as a discussion of the challenges intrinsic to highly-calibratable system development. Following this, we share laboratory results based on the calculation of noise wave parameters utilising an over-constrai…
▽ More
We detail the the REACH radiometric system designed to enable measurements of the 21-cm neutral hydrogen line. Included is the radiometer architecture and end-to-end system simulations as well as a discussion of the challenges intrinsic to highly-calibratable system development. Following this, we share laboratory results based on the calculation of noise wave parameters utilising an over-constrained least squares approach demonstrating a calibration RMSE of 80 mK for five hours of integration on a custom-made source with comparable impedance to that of the antenna used in the field. This paper therefore documents the state of the calibrator and data analysis in December 2022 in Cambridge before shipping to South Africa.
△ Less
Submitted 14 July, 2023; v1 submitted 30 June, 2023;
originally announced July 2023.
-
Signatures of Cosmic Ray Heating in 21-cm Observables
Authors:
T. Gessey-Jones,
A. Fialkov,
E. de Lera Acedo,
W. J. Handley,
R. Barkana
Abstract:
Cosmic rays generated by supernovae carry away a significant portion of the lifetime energy emission of their parent star, making them a plausible mechanism for heating the early universe intergalactic medium (IGM). Following a review of the existing literature on cosmic ray heating, we develop a flexible model of this heating mechanism for use in 3D semi-numerical 21-cm signal simulations and con…
▽ More
Cosmic rays generated by supernovae carry away a significant portion of the lifetime energy emission of their parent star, making them a plausible mechanism for heating the early universe intergalactic medium (IGM). Following a review of the existing literature on cosmic ray heating, we develop a flexible model of this heating mechanism for use in 3D semi-numerical 21-cm signal simulations and conduct the first investigations of the signatures it imprints on the 21-cm power spectrum and tomographic maps. We find that cosmic ray heating of the IGM is short-ranged, leading to heating clustered around star-forming sites, and a sharp contrast between heated regions of 21-cm emission and unheated regions of absorption. This contrast results in greater small-scale power for cosmic ray heated scenarios compared to what is found for X-ray heating, thus suggesting a way to test the nature of IGM heating with future 21-cm observations. Finally, we find an unexpectedly rich thermal history in models where cosmic rays can only escape efficiently from low-mass halos, such as in scenarios where these energetic particles originate from population III star supernovae remnants. The interplay of heating and the Lyman-Werner feedback in these models can produce a local peak in the IGM kinetic temperature and, for a limited parameter range, a flattened absorption trough in the global 21-cm signal.
△ Less
Submitted 3 October, 2023; v1 submitted 14 April, 2023;
originally announced April 2023.
-
Joint analysis constraints on the physics of the first galaxies with low frequency radio astronomy data
Authors:
Harry T. J. Bevins,
Stefan Heimersheim,
Irene Abril-Cabezas,
Anastasia Fialkov,
Eloy de Lera Acedo,
William Handley,
Saurabh Singh,
Rennan Barkana
Abstract:
The first billion years of cosmic history remains largely unobserved. We demonstrate, using a novel machine learning technique, how combining upper limits on the spatial fluctuations in the 21-cm signal with observations of the sky-averaged 21-cm signal from neutral hydrogen can improve our understanding of this epoch. By jointly analysing data from SARAS3 (redshift $z\approx15-25$) and limits fro…
▽ More
The first billion years of cosmic history remains largely unobserved. We demonstrate, using a novel machine learning technique, how combining upper limits on the spatial fluctuations in the 21-cm signal with observations of the sky-averaged 21-cm signal from neutral hydrogen can improve our understanding of this epoch. By jointly analysing data from SARAS3 (redshift $z\approx15-25$) and limits from HERA ($z\approx8$ and $10$), we show that such a synergetic analysis provides tighter constraints on the astrophysics of galaxies 200 million years after the Big Bang than can be achieved with the individual data sets. Although our constraints are weak, this is the first time data from a sky-averaged 21-cm experiment and power spectrum experiment have been analysed together. In synergy, the two experiments leave only $64.9^{+0.3}_{-0.1}$% of the explored broad theoretical parameter space to be consistent with the joint data set, in comparison to $92.3^{+0.3}_{-0.1}$% for SARAS3 and $79.0^{+0.5}_{-0.2}$% for HERA alone. We use the joint analysis to constrain star formation efficiency, minimum halo mass for star formation, X-ray luminosity of early emitters and the radio luminosity of early galaxies. The joint analysis disfavours at 68% confidence a combination of galaxies with X-ray emission that is $\lesssim 33$ and radio emission that is $\gtrsim 32$ times as efficient as present day galaxies. We disfavour at $95$% confidence scenarios in which power spectra are $\geq126$ mK$^{2}$ at $z=25$ and the sky-averaged signals are $\leq-277$ mK.
△ Less
Submitted 18 December, 2023; v1 submitted 9 January, 2023;
originally announced January 2023.
-
Sky-averaged 21-cm signal extraction using multiple antennas with an SVD framework: the REACH case
Authors:
Anchal Saxena,
P. Daniel Meerburg,
Eloy de Lera Acedo,
Will Handley,
Léon V. E. Koopmans
Abstract:
In a sky-averaged 21-cm signal experiment, the uncertainty on the extracted signal depends mainly on the covariance between the foreground and 21-cm signal models. In this paper, we construct these models using the modes of variation obtained from the Singular Value Decomposition of a set of simulated foreground and 21-cm signals. We present a strategy to reduce this overlap between the 21-cm and…
▽ More
In a sky-averaged 21-cm signal experiment, the uncertainty on the extracted signal depends mainly on the covariance between the foreground and 21-cm signal models. In this paper, we construct these models using the modes of variation obtained from the Singular Value Decomposition of a set of simulated foreground and 21-cm signals. We present a strategy to reduce this overlap between the 21-cm and foreground modes by simultaneously fitting the spectra from multiple different antennas, which can be used in combination with the method of utilizing the time dependence of foregrounds while fitting multiple drift scan spectra. To demonstrate this idea, we consider two different foreground models (i) a simple foreground model, where we assume a constant spectral index over the sky, and (ii) a more realistic foreground model, with a spatial variation of the spectral index. For the simple foreground model, with just a single antenna design, we are able to extract the signal with good accuracy if we simultaneously fit the data from multiple time slices. The 21-cm signal extraction is further improved when we simultaneously fit the data from different antennas as well. This improvement becomes more pronounced while using the more realistic mock observations generated from the detailed foreground model. We find that even if we fit multiple time slices, the recovered signal is biased and inaccurate for a single antenna. However, simultaneously fitting the data from different antennas reduces the bias and the uncertainty by a factor of 2-3 on the extracted 21-cm signal.
△ Less
Submitted 6 April, 2023; v1 submitted 14 December, 2022;
originally announced December 2022.
-
SuperNest: accelerated nested sampling applied to astrophysics and cosmology
Authors:
Aleksandr Petrosyan,
William James Handley
Abstract:
We present a method for improving the performance of nested
sampling as well as its accuracy. Building on previous work by
Chen et al., we show that posterior repartitioning
may be used to reduce the amount of time nested sampling spends in
compressing from prior to posterior if a suitable ``proposal''
distribution is supplied. We showcase this on a cosmological example
with a Gaussian…
▽ More
We present a method for improving the performance of nested
sampling as well as its accuracy. Building on previous work by
Chen et al., we show that posterior repartitioning
may be used to reduce the amount of time nested sampling spends in
compressing from prior to posterior if a suitable ``proposal''
distribution is supplied. We showcase this on a cosmological example
with a Gaussian posterior, and release the code as an LGPL licensed,
extensible Python package
https://gitlab.com/a-p-petrosyan/sspr.
△ Less
Submitted 4 December, 2022;
originally announced December 2022.
-
Astrophysical Constraints from the SARAS3 non-detection of the Cosmic Dawn Sky-Averaged 21-cm Signal
Authors:
H. T. J. Bevins,
A. Fialkov,
E. de Lera Acedo,
W. J. Handley,
S. Singh,
R. Subrahmanyan,
R. Barkana
Abstract:
Observations of the redshifted 21-cm line of atomic hydrogen have resulted in several upper limits on the 21-cm power spectrum and a tentative detection of the sky-averaged signal at $z\sim17$. Made with the EDGES Low-Band antenna, this claim was recently disputed by the SARAS3 experiment, which reported a non-detection and is the only available upper limit strong enough to constrain cosmic dawn a…
▽ More
Observations of the redshifted 21-cm line of atomic hydrogen have resulted in several upper limits on the 21-cm power spectrum and a tentative detection of the sky-averaged signal at $z\sim17$. Made with the EDGES Low-Band antenna, this claim was recently disputed by the SARAS3 experiment, which reported a non-detection and is the only available upper limit strong enough to constrain cosmic dawn astrophysics. We use these data to constrain a population of radio-luminous galaxies $\sim 200$ million years after the Big Bang ($z\approx 20$). We find, using Bayesian data analysis, that the data disfavours (at 68% confidence) radio-luminous galaxies in dark matter halos with masses of $4.4\times10^{5}$ M$_\odot \lesssim M \lesssim 1.1\times10^{7}$M$_\odot$ (where $M_\odot$ is the mass of the Sun) at $z = 20$ and galaxies in which $>5$% of the gas is converted into stars. The data disfavour galaxies with radio luminosity per star formation rate of $L_\mathrm{r}/\mathrm{SFR} \gtrsim 1.549 \times 10^{25}$ W Hz$^{-1}$M$_\odot^{-1}$ yr at 150 MHz, a thousand times brighter than today, and, separately, a synchrotron radio background in excess of the CMB by $\gtrsim 6%$ at 1.42 GHz.
△ Less
Submitted 1 December, 2022;
originally announced December 2022.
-
Quantum initial conditions for curved inflating universes
Authors:
Mary I. Letey,
Zakhar Shumaylov,
Fruzsina J. Agocs,
Will J. Handley,
Michael P. Hobson,
Anthony N. Lasenby
Abstract:
We discuss the challenges of motivating, constructing, and quantizing a canonically normalized inflationary perturbation in spatially curved universes. We show that this has historically proved challenging due to the interaction of nonadiabaticity with spatial curvature. We construct a novel curvature perturbation that is canonically normalized in the sense of its equation of motion and is unique…
▽ More
We discuss the challenges of motivating, constructing, and quantizing a canonically normalized inflationary perturbation in spatially curved universes. We show that this has historically proved challenging due to the interaction of nonadiabaticity with spatial curvature. We construct a novel curvature perturbation that is canonically normalized in the sense of its equation of motion and is unique up to a single scalar parameter. With this construction it becomes possible to set initial conditions invariant under canonical transformations, overcoming known ambiguities in the literature. This corrected quantization has potentially observational consequences via modifications to the primordial power spectrum at large angular scales, as well as theoretical implications for quantization procedures in curved cosmologies filled with a scalar field.
△ Less
Submitted 8 July, 2024; v1 submitted 30 November, 2022;
originally announced November 2022.
-
Bayesian approach to radio frequency interference mitigation
Authors:
S. A. K. Leeney,
W. J. Handley,
E. de Lera Acedo
Abstract:
Interfering signals such as Radio Frequency Interference from ubiquitous satellite constellations are becoming an endemic problem in fields involving physical observations of the electromagnetic spectrum. To address this we propose a novel data cleaning methodology. Contamination is simultaneously flagged and managed at the likelihood level. It is modeled in a Bayesian fashion through a piecewise…
▽ More
Interfering signals such as Radio Frequency Interference from ubiquitous satellite constellations are becoming an endemic problem in fields involving physical observations of the electromagnetic spectrum. To address this we propose a novel data cleaning methodology. Contamination is simultaneously flagged and managed at the likelihood level. It is modeled in a Bayesian fashion through a piecewise likelihood that is constrained by a Bernoulli prior distribution. The techniques described in this paper can be implemented with just a few lines of code.
△ Less
Submitted 19 February, 2024; v1 submitted 28 November, 2022;
originally announced November 2022.
-
A General Bayesian Framework to Account for Foreground Map Errors in Global 21-cm Experiments
Authors:
Michael Pagano,
Peter Sims,
Adrian Liu,
Dominic Anstey,
Will Handley,
Eloy De Lera Acedo
Abstract:
Measurement of the global 21-cm signal during Cosmic Dawn (CD) and the Epoch of Reionization (EoR) is made difficult by bright foreground emission which is 2-5 orders of magnitude larger than the expected signal. Fitting for a physics-motivated parametric forward model of the data within a Bayesian framework provides a robust means to separate the signal from the foregrounds, given sufficient info…
▽ More
Measurement of the global 21-cm signal during Cosmic Dawn (CD) and the Epoch of Reionization (EoR) is made difficult by bright foreground emission which is 2-5 orders of magnitude larger than the expected signal. Fitting for a physics-motivated parametric forward model of the data within a Bayesian framework provides a robust means to separate the signal from the foregrounds, given sufficient information about the instrument and sky. It has previously been demonstrated that, within such a modelling framework, a foreground model of sufficient fidelity can be generated by dividing the sky into $N$ regions and scaling a base map assuming a distinct uniform spectral index in each region. Using the Radio Experiment for the Analysis of Cosmic Hydrogen (REACH) as our fiducial instrument, we show that, if unaccounted-for, amplitude errors in low-frequency radio maps used for our base map model will prevent recovery of the 21-cm signal within this framework, and that the level of bias in the recovered 21-cm signal is proportional to the amplitude and the correlation length of the base-map errors in the region. We introduce an updated foreground model that is capable of accounting for these measurement errors by fitting for a monopole offset and a set of spatially-dependent scale factors describing the ratio of the true and model sky temperatures, with the size of the set determined by Bayesian evidence-based model comparison. We show that our model is flexible enough to account for multiple foreground error scenarios allowing the 21-cm sky-averaged signal to be detected without bias from simulated observations with a smooth conical log spiral antenna.
△ Less
Submitted 1 November, 2023; v1 submitted 18 November, 2022;
originally announced November 2022.
-
The REACH radiometer for detecting the 21-cm hydrogen signal from redshift 7.5 to 28
Authors:
E. de Lera Acedo,
D. I. L. de Villiers,
N. Razavi-Ghods,
W. Handley,
A. Fialkov,
A. Magro,
D. Anstey,
H. T. J. Bevins,
R. Chiello,
J. Cumner,
A. T. Josaitis,
I. L. V. Roque,
P. H. Sims,
K. H. Scheutwinkel,
P. Alexander,
G. Bernardi,
S. Carey,
J. Cavillot,
W. Croukamp,
J. A. Ely,
T. Gessey-Jones,
Q. Gueuning,
R. Hills,
G. Kulkarni,
R. Maiolino
, et al. (9 additional authors not shown)
Abstract:
Observations of the 21-cm line from primordial hydrogen promise to be one of the best tools to study the early epochs of the Universe: the Dark Ages, the Cosmic Dawn, and the subsequent Epoch of Reionization. In 2018, the EDGES experiment caught the attention of the cosmology community with a potential detection of an absorption feature in the sky-averaged radio spectrum centred at 78 MHz. The fea…
▽ More
Observations of the 21-cm line from primordial hydrogen promise to be one of the best tools to study the early epochs of the Universe: the Dark Ages, the Cosmic Dawn, and the subsequent Epoch of Reionization. In 2018, the EDGES experiment caught the attention of the cosmology community with a potential detection of an absorption feature in the sky-averaged radio spectrum centred at 78 MHz. The feature is deeper than expected, and, if confirmed, would call for new physics. However, different groups have re-analyzed the EDGES data and questioned the reliability of the signal. The Radio Experiment for the Analysis of Cosmic Hydrogen (REACH) is a sky-averaged 21-cm experiment aiming at improving the current observations by tackling the issues faced by current instruments related to residual systematic signals in the data. The novel experimental approach focuses on detecting and jointly explaining these systematics together with the foregrounds and the cosmological signal using Bayesian statistics. To achieve this, REACH features simultaneous observations with two different antennas, an ultra wideband system (redshift range 7.5 to 28), and a receiver calibrator based on in-field measurements. Simulated observations forecast percent-level constraints on astrophysical parameters, potentially opening up a new window to the infant Universe.
△ Less
Submitted 13 October, 2022;
originally announced October 2022.
-
Use of Time Dependent Data in Bayesian Global 21cm Foreground and Signal Modelling
Authors:
Dominic Anstey,
Eloy de Lera Acedo,
Will Handley
Abstract:
Global 21cm cosmology aims to investigate the cosmic dawn and epoch of reionisation by measuring the sky averaged HI absorption signal, which requires, accurate modelling of, or correction for, the bright radio foregrounds and distortions arising from chromaticity of the antenna beam. We investigate the effect of improving foreground modelling by fitting data sets from many observation times simul…
▽ More
Global 21cm cosmology aims to investigate the cosmic dawn and epoch of reionisation by measuring the sky averaged HI absorption signal, which requires, accurate modelling of, or correction for, the bright radio foregrounds and distortions arising from chromaticity of the antenna beam. We investigate the effect of improving foreground modelling by fitting data sets from many observation times simultaneously in a single Bayesian analysis, fitting for the same parameter set by performing these fits on simulated data. We find that for a hexagonal dipole antenna, this simultaneous fitting produces a significant improvement in the accuracy of the recovered 21cm signal, relative to fitting a time average of the data. Furthermore, the recovered models of the foreground are also seen to become more accurate by up to a factor of $\sim$2-3 relative to time averaged fitting. For a less chromatic log spiral antenna, no significant improvement in signal recovery was found by this process. However, the modelling of the foregrounds was still significantly improved. We also investigate extending this technique to fit multiple data sets from different antennae simultaneously for the same parameters. This is also found to improve both 21cm signal and foreground modelling, to a higher degree than fitting data set from multiple times from the same antenna.
△ Less
Submitted 16 January, 2023; v1 submitted 10 October, 2022;
originally announced October 2022.
-
Marginal Bayesian Statistics Using Masked Autoregressive Flows and Kernel Density Estimators with Examples in Cosmology
Authors:
Harry Bevins,
Will Handley,
Pablo Lemos,
Peter Sims,
Eloy de Lera Acedo,
Anastasia Fialkov
Abstract:
Cosmological experiments often employ Bayesian workflows to derive constraints on cosmological and astrophysical parameters from their data. It has been shown that these constraints can be combined across different probes such as Planck and the Dark Energy Survey and that this can be a valuable exercise to improve our understanding of the universe and quantify tension between multiple experiments.…
▽ More
Cosmological experiments often employ Bayesian workflows to derive constraints on cosmological and astrophysical parameters from their data. It has been shown that these constraints can be combined across different probes such as Planck and the Dark Energy Survey and that this can be a valuable exercise to improve our understanding of the universe and quantify tension between multiple experiments. However, these experiments are typically plagued by differing systematics, instrumental effects and contaminating signals, which we collectively refer to as `nuisance' components, that have to be modelled alongside target signals of interest. This leads to high dimensional parameter spaces, especially when combining data sets, with > 20 dimensions of which only around 5 correspond to key physical quantities. We present a means by which to combine constraints from different data sets in a computationally efficient manner by generating rapid, reusable and reliable marginal probability density estimators, giving us access to nuisance-free likelihoods. This is possible through the unique combination of nested sampling, which gives us access to Bayesian evidences, and the marginal Bayesian statistics code MARGARINE. Our method is lossless in the signal parameters, resulting in the same posterior distributions as would be found from a full nested sampling run over all nuisance parameters, and typically quicker than evaluating full likelihoods. We demonstrate our approach by applying it to the combination of posteriors from the Dark Energy Survey and Planck.
△ Less
Submitted 24 November, 2022; v1 submitted 23 July, 2022;
originally announced July 2022.
-
Nested sampling for physical scientists
Authors:
Greg Ashton,
Noam Bernstein,
Johannes Buchner,
Xi Chen,
Gábor Csányi,
Andrew Fowlie,
Farhan Feroz,
Matthew Griffiths,
Will Handley,
Michael Habeck,
Edward Higson,
Michael Hobson,
Anthony Lasenby,
David Parkinson,
Livia B. Pártay,
Matthew Pitkin,
Doris Schneider,
Joshua S. Speagle,
Leah South,
John Veitch,
Philipp Wacker,
David J. Wales,
David Yallup
Abstract:
We review Skilling's nested sampling (NS) algorithm for Bayesian inference and more broadly multi-dimensional integration. After recapitulating the principles of NS, we survey developments in implementing efficient NS algorithms in practice in high-dimensions, including methods for sampling from the so-called constrained prior. We outline the ways in which NS may be applied and describe the applic…
▽ More
We review Skilling's nested sampling (NS) algorithm for Bayesian inference and more broadly multi-dimensional integration. After recapitulating the principles of NS, we survey developments in implementing efficient NS algorithms in practice in high-dimensions, including methods for sampling from the so-called constrained prior. We outline the ways in which NS may be applied and describe the application of NS in three scientific fields in which the algorithm has proved to be useful: cosmology, gravitational-wave astronomy, and materials science. We close by making recommendations for best practice when using NS and by summarizing potential limitations and optimizations of NS.
△ Less
Submitted 31 May, 2022;
originally announced May 2022.
-
Cosmological constraints on decaying axion-like particles: a global analysis
Authors:
Csaba Balázs,
Sanjay Bloor,
Tomás E. Gonzalo,
Will Handley,
Sebastian Hoof,
Felix Kahlhoefer,
Marie Lecroq,
David J. E. Marsh,
Janina J. Renk,
Pat Scott,
Patrick Stöcker
Abstract:
Axion-like particles (ALPs) decaying into photons are known to affect a wide range of astrophysical and cosmological observables. In this study we focus on ALPs with masses in the keV-MeV range and lifetimes between $10^4$ and $10^{13}$ seconds, corresponding to decays between the end of Big Bang Nucleosynthesis and the formation of the Cosmic Microwave Background (CMB). Using the CosmoBit module…
▽ More
Axion-like particles (ALPs) decaying into photons are known to affect a wide range of astrophysical and cosmological observables. In this study we focus on ALPs with masses in the keV-MeV range and lifetimes between $10^4$ and $10^{13}$ seconds, corresponding to decays between the end of Big Bang Nucleosynthesis and the formation of the Cosmic Microwave Background (CMB). Using the CosmoBit module of the global fitting framework GAMBIT, we combine state-of-the-art calculations of the irreducible ALP freeze-in abundance, primordial element abundances (including photodisintegration through ALP decays), CMB spectral distortions and anisotropies, and constraints from supernovae and stellar cooling. This approach makes it possible for the first time to perform a global analysis of the ALP parameter space while varying the parameters of $Λ$CDM as well as several nuisance parameters. We find a lower bound on the ALP mass of around $m_a > 300\,\text{keV}$, which can only be evaded if ALPs are stable on cosmological timescales. Future observations of CMB spectral distortions with a PIXIE-like mission are expected to improve this bound by two orders of magnitude.
△ Less
Submitted 19 December, 2022; v1 submitted 26 May, 2022;
originally announced May 2022.
-
Marginal Post Processing of Bayesian Inference Products with Normalizing Flows and Kernel Density Estimators
Authors:
Harry T. J. Bevins,
William J. Handley,
Pablo Lemos,
Peter H. Sims,
Eloy de Lera Acedo,
Anastasia Fialkov,
Justin Alsing
Abstract:
Bayesian analysis has become an indispensable tool across many different cosmological fields including the study of gravitational waves, the Cosmic Microwave Background and the 21-cm signal from the Cosmic Dawn among other phenomena. The method provides a way to fit complex models to data describing key cosmological and astrophysical signals and a whole host of contaminating signals and instrument…
▽ More
Bayesian analysis has become an indispensable tool across many different cosmological fields including the study of gravitational waves, the Cosmic Microwave Background and the 21-cm signal from the Cosmic Dawn among other phenomena. The method provides a way to fit complex models to data describing key cosmological and astrophysical signals and a whole host of contaminating signals and instrumental effects modelled with `nuisance parameters'. In this paper, we summarise a method that uses Masked Autoregressive Flows and Kernel Density Estimators to learn marginal posterior densities corresponding to core science parameters. We find that the marginal or 'nuisance-free' posteriors and the associated likelihoods have an abundance of applications including; the calculation of previously intractable marginal Kullback-Leibler divergences and marginal Bayesian Model Dimensionalities, likelihood emulation and prior emulation. We demonstrate each application using toy examples, examples from the field of 21-cm cosmology and samples from the Dark Energy Survey. We discuss how marginal summary statistics like the Kullback-Leibler divergences and Bayesian Model Dimensionalities can be used to examine the constraining power of different experiments and how we can perform efficient joint analysis by taking advantage of marginal prior and likelihood emulators. We package our multipurpose code up in the pip-installable code margarine for use in the wider scientific community.
△ Less
Submitted 18 December, 2023; v1 submitted 25 May, 2022;
originally announced May 2022.
-
Finite inflation in curved space
Authors:
Lukas T. Hergt,
Fruzsina J. Agocs,
Will J. Handley,
Michael P. Hobson,
Anthony N. Lasenby
Abstract:
We investigate the effects of non-zero spatial curvature on cosmic inflation in the light of cosmic microwave background (CMB) anisotropy measurements from the Planck 2018 legacy release and from the 2015 observing season of BICEP2 and the Keck Array. Even a small percentage of non-zero curvature today would significantly limit the total number of e-folds of the scale factor during inflation, rend…
▽ More
We investigate the effects of non-zero spatial curvature on cosmic inflation in the light of cosmic microwave background (CMB) anisotropy measurements from the Planck 2018 legacy release and from the 2015 observing season of BICEP2 and the Keck Array. Even a small percentage of non-zero curvature today would significantly limit the total number of e-folds of the scale factor during inflation, rendering just-enough inflation scenarios with a kinetically dominated or fast-roll stage prior to slow-roll inflation more likely. Finite inflation leads to oscillations and a cutoff towards large scales in the primordial power spectrum and curvature pushes them into the CMB observable window. Using nested sampling, we carry out Bayesian parameter estimations and model comparisons taking into account constraints from reheating and horizon considerations. We confirm the preference of CMB data for closed universes with Bayesian odds of over $100:1$ and with a posterior on the curvature density parameter of $Ω_{K,0}=-0.051\pm0.017$ for a curvature extension of LCDM and $Ω_{K,0}=-0.031\pm0.014$ for Starobinsky inflation. Model comparisons of various inflation models give similar results as for flat universes with the Starobinsky model outperforming most other models.
△ Less
Submitted 15 May, 2022;
originally announced May 2022.
-
Bayesian evidence-driven likelihood selection for sky-averaged 21-cm signal extraction
Authors:
K. H. Scheutwinkel,
W. Handley,
E. de Lera Acedo
Abstract:
We demonstrate that the Bayesian evidence can be used to find a good approximation of the ground truth likelihood function of a dataset, a goal of the likelihood-free inference (LFI) paradigm. As a concrete example, we use forward modelled sky-averaged 21-cm signal antenna temperature datasets where we artificially inject noise structures of various physically motivated forms. We find that the Gau…
▽ More
We demonstrate that the Bayesian evidence can be used to find a good approximation of the ground truth likelihood function of a dataset, a goal of the likelihood-free inference (LFI) paradigm. As a concrete example, we use forward modelled sky-averaged 21-cm signal antenna temperature datasets where we artificially inject noise structures of various physically motivated forms. We find that the Gaussian likelihood performs poorly when the noise distribution deviates from the Gaussian case e.g. heteroscedastic radiometric or heavy-tailed noise. For these non-Gaussian noise structures, we show that the generalised normal likelihood is on a similar Bayesian evidence scale with comparable sky-averaged 21-cm signal recovery as the ground truth likelihood function of our injected noise. We therefore propose the generalised normal likelihood function as a good approximation of the true likelihood function if the noise structure is a priori unknown.
△ Less
Submitted 20 March, 2023; v1 submitted 9 April, 2022;
originally announced April 2022.
-
Bayesian evidence-driven diagnosis of instrumental systematics for sky-averaged 21-cm cosmology experiments
Authors:
K. H. Scheutwinkel,
E. de Lera Acedo,
W. Handley
Abstract:
We demonstrate the effectiveness of a Bayesian evidence-based analysis for diagnosing and disentangling the sky-averaged 21-cm signal from instrumental systematic effects. As a case study, we consider a simulated REACH pipeline with an injected systematic. We demonstrate that very poor performance or erroneous signal recovery is achieved if the systematic remains unmodelled. These effects include…
▽ More
We demonstrate the effectiveness of a Bayesian evidence-based analysis for diagnosing and disentangling the sky-averaged 21-cm signal from instrumental systematic effects. As a case study, we consider a simulated REACH pipeline with an injected systematic. We demonstrate that very poor performance or erroneous signal recovery is achieved if the systematic remains unmodelled. These effects include sky-averaged 21-cm posterior estimates resembling a very deep or wide signal. However, when including parameterised models of the systematic, the signal recovery is dramatically improved in performance. Most importantly, a Bayesian evidence-based model comparison is capable of determining whether or not such a systematic model is needed as the true underlying generative model of an experimental dataset is in principle unknown. We, therefore, advocate a pipeline capable of testing a variety of potential systematic errors with the Bayesian evidence acting as the mechanism for detecting their presence.
△ Less
Submitted 23 September, 2022; v1 submitted 9 April, 2022;
originally announced April 2022.
-
Cosmology Intertwined: A Review of the Particle Physics, Astrophysics, and Cosmology Associated with the Cosmological Tensions and Anomalies
Authors:
Elcio Abdalla,
Guillermo Franco Abellán,
Amin Aboubrahim,
Adriano Agnello,
Ozgur Akarsu,
Yashar Akrami,
George Alestas,
Daniel Aloni,
Luca Amendola,
Luis A. Anchordoqui,
Richard I. Anderson,
Nikki Arendse,
Marika Asgari,
Mario Ballardini,
Vernon Barger,
Spyros Basilakos,
Ronaldo C. Batista,
Elia S. Battistelli,
Richard Battye,
Micol Benetti,
David Benisty,
Asher Berlin,
Paolo de Bernardis,
Emanuele Berti,
Bohdan Bidenko
, et al. (178 additional authors not shown)
Abstract:
In this paper we will list a few important goals that need to be addressed in the next decade, also taking into account the current discordances between the different cosmological probes, such as the disagreement in the value of the Hubble constant $H_0$, the $σ_8$--$S_8$ tension, and other less statistically significant anomalies. While these discordances can still be in part the result of system…
▽ More
In this paper we will list a few important goals that need to be addressed in the next decade, also taking into account the current discordances between the different cosmological probes, such as the disagreement in the value of the Hubble constant $H_0$, the $σ_8$--$S_8$ tension, and other less statistically significant anomalies. While these discordances can still be in part the result of systematic errors, their persistence after several years of accurate analysis strongly hints at cracks in the standard cosmological scenario and the necessity for new physics or generalisations beyond the standard model. In this paper, we focus on the $5.0\,σ$ tension between the {\it Planck} CMB estimate of the Hubble constant $H_0$ and the SH0ES collaboration measurements. After showing the $H_0$ evaluations made from different teams using different methods and geometric calibrations, we list a few interesting new physics models that could alleviate this tension and discuss how the next decade's experiments will be crucial. Moreover, we focus on the tension of the {\it Planck} CMB data with weak lensing measurements and redshift surveys, about the value of the matter energy density $Ω_m$, and the amplitude or rate of the growth of structure ($σ_8,fσ_8$). We list a few interesting models proposed for alleviating this tension, and we discuss the importance of trying to fit a full array of data with a single model and not just one parameter at a time. Additionally, we present a wide range of other less discussed anomalies at a statistical significance level lower than the $H_0$--$S_8$ tensions which may also constitute hints towards new physics, and we discuss possible generic theoretical approaches that can collectively explain the non-standard nature of these signals.[Abridged]
△ Less
Submitted 24 April, 2022; v1 submitted 11 March, 2022;
originally announced March 2022.
-
Impact of the Primordial Stellar Initial Mass Function on the 21-cm Signal
Authors:
T. Gessey-Jones,
N. S. Sartorio,
A. Fialkov,
G. M. Mirouh,
M. Magg,
R. G. Izzard,
E. de Lera Acedo,
W. J. Handley,
R. Barkana
Abstract:
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 invest…
▽ More
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.
△ Less
Submitted 21 July, 2022; v1 submitted 4 February, 2022;
originally announced February 2022.
-
A Comprehensive Bayesian re-analysis of the SARAS2 data from the Epoch of Reionization
Authors:
H. T. J. Bevins,
E. de Lera Acedo,
A. Fialkov,
W. J. Handley,
S. Singh,
R. Subrahmanyan,
R. Barkana
Abstract:
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…
▽ More
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.
△ Less
Submitted 20 May, 2022; v1 submitted 27 January, 2022;
originally announced January 2022.
-
Primordial power spectra from $k$-inflation with curvature
Authors:
Zakhar Shumaylov,
Will Handley
Abstract:
We investigate the primordial power spectra for general kinetic inflation models that support a period of kinetic dominance in the case of curved universes. We present derivations of the Mukhanov-Sasaki equations with a nonstandard scalar kinetic Lagrangian which manifests itself through the inflationary sound speed $c_s^2$. We extend the analytical approximations exploited in Contaldi et al [1] a…
▽ More
We investigate the primordial power spectra for general kinetic inflation models that support a period of kinetic dominance in the case of curved universes. We present derivations of the Mukhanov-Sasaki equations with a nonstandard scalar kinetic Lagrangian which manifests itself through the inflationary sound speed $c_s^2$. We extend the analytical approximations exploited in Contaldi et al [1] and Thavanesan et al [2] to general kinetic Lagrangians and show the effect of $k$-inflation on the primordial power spectra for models with curvature. In particular, the interplay between sound speed and curvature results in a natural low wavenumber cutoff for the power spectra in the case of closed universes. Using the analytical approximation, we further show that a change in the inflationary sound speed between different epochs of inflation results in non-decaying oscillations in the resultant power spectra for the comoving curvature perturbation.
△ Less
Submitted 24 June, 2022; v1 submitted 14 December, 2021;
originally announced December 2021.
-
Rescuing Palindromic Universes with Improved Recombination Modelling
Authors:
Metha Prathaban,
Will Handley
Abstract:
We explore the linearly quantised primordial power spectra associated with palindromic universes. Extending the results of Lasenby et al. [1] and Bartlett et al. [2], we improve the modelling of recombination and include higher orders in the photonic Boltzmann hierarchy. In so doing, we find that the predicted power spectra become largely consistent with observational data. The improved recombinat…
▽ More
We explore the linearly quantised primordial power spectra associated with palindromic universes. Extending the results of Lasenby et al. [1] and Bartlett et al. [2], we improve the modelling of recombination and include higher orders in the photonic Boltzmann hierarchy. In so doing, we find that the predicted power spectra become largely consistent with observational data. The improved recombination modelling involves developing further techniques for dealing with the future conformal boundary, by integrating the associated perturbation equations both forwards and backwards in conformal time. The resulting wavevector quantisation gives a lowest allowed wavenumber ${k_0 = 9.93 \times 10^{-5} \textrm{Mpc}^{-1}}$ and linear spacing ${Δk = 1.63 \times 10^{-4} \textrm{Mpc}^{-1}}$, providing fits consistent with observational data equivalent in quality to the $Λ$CDM model.
△ Less
Submitted 10 June, 2022; v1 submitted 29 November, 2021;
originally announced November 2021.
-
Radio antenna design for sky-averaged 21 cm cosmology experiments: the REACH case
Authors:
J. Cumner,
E. De Lera Acedo,
D. I. L. de Villiers,
D. Anstey,
C. I. Kolitsidas,
B. Gurdon,
N. Fagnoni,
P. Alexander,
G. Bernardi,
H. T. J. Bevins,
S. Carey,
J. Cavillot,
R. Chiello,
C. Craeye,
W. Croukamp,
J. A. Ely,
A. Fialkov,
T. Gessey-Jones,
Q. Gueuning,
W. Handley,
R. Hills,
A. T. Josaitis,
G. Kulkarni,
A. Magro,
R. Maiolino
, et al. (13 additional authors not shown)
Abstract:
Following the reported detection of an absorption profile associated with the 21~cm sky-averaged signal from the Cosmic Dawn by the EDGES experiment in 2018, a number of experiments have been set up to verify this result. This paper discusses the design process used for global 21~cm experiments, focusing specifically on the Radio Experiment for the Analysis of Cosmic Hydrogen (REACH). This experim…
▽ More
Following the reported detection of an absorption profile associated with the 21~cm sky-averaged signal from the Cosmic Dawn by the EDGES experiment in 2018, a number of experiments have been set up to verify this result. This paper discusses the design process used for global 21~cm experiments, focusing specifically on the Radio Experiment for the Analysis of Cosmic Hydrogen (REACH). This experiment will seek to understand and compensate for systematic errors present using detailed modelling and characterization of the instrumentation. There is detailed the quantitative figures of merit and numerical modelling used to assist the design process of the REACH dipole antenna (one of the 2 antenna designs for REACH Phase I). This design process produced a 2.5:1 frequency bandwidth dipole. The aim of this design was to balance spectral smoothness and low impedance reflections with the ability to describe and understand the antenna response to the sky signal to inform the critically important calibration during observation and data analysis.
△ Less
Submitted 12 January, 2023; v1 submitted 21 September, 2021;
originally announced September 2021.
-
Informing antenna design for sky-averaged 21-cm experiments using a simulated Bayesian data analysis pipeline
Authors:
Dominic Anstey,
John Cumner,
Eloy de Lera Acedo,
Will Handley
Abstract:
Global 21cm experiments aim to measure the sky averaged HI absorption signal from cosmic dawn and the epoch of reionisation. However, antenna chromaticity coupling to bright foregrounds can introduce distortions into the observational data of such experiments. We demonstrate a method for guiding the antenna design of a global experiment through data analysis simulations. This is done by performing…
▽ More
Global 21cm experiments aim to measure the sky averaged HI absorption signal from cosmic dawn and the epoch of reionisation. However, antenna chromaticity coupling to bright foregrounds can introduce distortions into the observational data of such experiments. We demonstrate a method for guiding the antenna design of a global experiment through data analysis simulations. This is done by performing simulated observations for a range of inserted 21cm signals, then attempting to identify the signals with a data analysis pipeline. We demonstrate this method on five antennae that were considered as potential designs for the Radio Experiment for the Analysis of Cosmic Hydrogen (REACH); a conical log spiral antenna, an inverted conical sinuous antenna and polygonal-, square- and elliptical-bladed dipoles. We find that the log spiral performs significantly better than the other antennae tested, able to correctly and confidently identify every inserted 21cm signal. In second place is the polygonal dipole antenna, which was only unable to detect signals with both very low amplitudes of 0.05K and low centre frequency of 80MHz. The conical sinuous antenna was found to perform least accurately, only able to detect the highest amplitude 21cm signals, and even then with biases. We also demonstrate that, due to the non-trivial nature of chromatic distortion and the processes of correcting for it, these are not the results that could have been expected superficially from the extent of chromatic variation in each antenna.
△ Less
Submitted 11 November, 2021; v1 submitted 18 June, 2021;
originally announced June 2021.
-
Thermal WIMPs and the Scale of New Physics: Global Fits of Dirac Dark Matter Effective Field Theories
Authors:
The GAMBIT Collaboration,
Peter Athron,
Neal Avis Kozar,
Csaba Balázs,
Ankit Beniwal,
Sanjay Bloor,
Torsten Bringmann,
Joachim Brod,
Christopher Chang,
Jonathan M. Cornell,
Ben Farmer,
Andrew Fowlie,
Tomás E. Gonzalo,
Will Handley,
Felix Kahlhoefer,
Anders Kvellestad,
Farvah Mahmoudi,
Markus T. Prim,
Are Raklev,
Janina J. Renk,
Andre Scaffidi,
Pat Scott,
Patrick Stöcker,
Aaron C. Vincent,
Martin White
, et al. (2 additional authors not shown)
Abstract:
We assess the status of a wide class of WIMP dark matter (DM) models in light of the latest experimental results using the global fitting framework $\textsf{GAMBIT}$. We perform a global analysis of effective field theory (EFT) operators describing the interactions between a gauge-singlet Dirac fermion and the Standard Model quarks, the gluons and the photon. In this bottom-up approach, we simulta…
▽ More
We assess the status of a wide class of WIMP dark matter (DM) models in light of the latest experimental results using the global fitting framework $\textsf{GAMBIT}$. We perform a global analysis of effective field theory (EFT) operators describing the interactions between a gauge-singlet Dirac fermion and the Standard Model quarks, the gluons and the photon. In this bottom-up approach, we simultaneously vary the coefficients of 14 such operators up to dimension 7, along with the DM mass, the scale of new physics and several nuisance parameters. Our likelihood functions include the latest data from $\mathit{Planck}$, direct and indirect detection experiments, and the LHC. For DM masses below 100 GeV, we find that it is impossible to satisfy all constraints simultaneously while maintaining EFT validity at LHC energies. For new physics scales around 1 TeV, our results are influenced by several small excesses in the LHC data and depend on the prescription that we adopt to ensure EFT validity. Furthermore, we find large regions of viable parameter space where the EFT is valid and the relic density can be reproduced, implying that WIMPs can still account for the DM of the universe while being consistent with the latest data.
△ Less
Submitted 13 November, 2021; v1 submitted 3 June, 2021;
originally announced June 2021.
-
Nested sampling for frequentist computation: fast estimation of small $p$-values
Authors:
Andrew Fowlie,
Sebastian Hoof,
Will Handley
Abstract:
We propose a novel method for computing $p$-values based on nested sampling (NS) applied to the sampling space rather than the parameter space of the problem, in contrast to its usage in Bayesian computation. The computational cost of NS scales as $\log^2{1/p}$, which compares favorably to the $1/p$ scaling for Monte Carlo (MC) simulations. For significances greater than about $4σ$ in both a toy p…
▽ More
We propose a novel method for computing $p$-values based on nested sampling (NS) applied to the sampling space rather than the parameter space of the problem, in contrast to its usage in Bayesian computation. The computational cost of NS scales as $\log^2{1/p}$, which compares favorably to the $1/p$ scaling for Monte Carlo (MC) simulations. For significances greater than about $4σ$ in both a toy problem and a simplified resonance search, we show that NS requires orders of magnitude fewer simulations than ordinary MC estimates. This is particularly relevant for high-energy physics, which adopts a $5σ$ gold standard for discovery. We conclude with remarks on new connections between Bayesian and frequentist computation and possibilities for tuning NS implementations for still better performance in this setting.
△ Less
Submitted 13 January, 2022; v1 submitted 27 May, 2021;
originally announced May 2021.
-
GLOBALEMU: A novel and robust approach for emulating the sky-averaged 21-cm signal from the cosmic dawn and epoch of reionisation
Authors:
H. T. J. Bevins,
W. J. Handley,
A. Fialkov,
E. de Lera Acedo,
K. Javid
Abstract:
Emulation of the Global (sky-averaged) 21-cm signal with neural networks has been shown to be an essential tool for physical signal modelling. In this paper we present globalemu, a Global 21-cm signal emulator that uses redshift as a character defining variable alongside a set of astrophysical parameters to estimate the signal brightness temperature. Combined with physically-motivated data pre-pro…
▽ More
Emulation of the Global (sky-averaged) 21-cm signal with neural networks has been shown to be an essential tool for physical signal modelling. In this paper we present globalemu, a Global 21-cm signal emulator that uses redshift as a character defining variable alongside a set of astrophysical parameters to estimate the signal brightness temperature. Combined with physically-motivated data pre-processing this makes for a reliable and fast emulator that is relatively insensitive to the network design. globalemu can emulate a high resolution signal in 1.3 ms in comparison to 133 ms, a factor of 102 improvement, when using the existing public state of the art 21cmGEM. We illustrate, with the standard astrophysical models used to train 21cmGEM, that globalemu is almost twice as accurate and for a test set of $\approx1,700$ signals we achieve a mean RMSE of 2.52 mK across the band z=7-28 ($\approx$10 per cent the expected noise of the Radio Experiment for the Analysis of Cosmic Hydrogen (REACH)). The models are parameterised by the star formation efficiency, $f_*$, minimum virial circular velocity, $V_c$, X-ray efficiency, $f_X$, CMB optical depth, $τ$, the slope and low energy cut off of the X-ray spectral energy density, $α$ and $ν_\mathrm{min}$, and the mean free path of ionizing photons, $R_\mathrm{mfp}$. globalemu provides a flexible framework for easily emulating updated simulations of the Global signal and in addition the neutral fraction history. The emulator is pip installable and available at: \url{https://github.com/htjb/globalemu}. globalemu will be used extensively by the REACH collaboration.
△ Less
Submitted 15 October, 2021; v1 submitted 9 April, 2021;
originally announced April 2021.
-
Constraining Quantum Initial Conditions before Inflation
Authors:
T. Gessey-Jones,
W. J. Handley
Abstract:
We theoretically and observationally investigate different choices of initial conditions for the primordial mode function that are imposed during an epoch preceding inflation. By deriving predictions for the observables resulting from several alternate quantum vacuum prescriptions we show some choices of vacua are theoretically observationally distinguishable from others. Comparing these predictio…
▽ More
We theoretically and observationally investigate different choices of initial conditions for the primordial mode function that are imposed during an epoch preceding inflation. By deriving predictions for the observables resulting from several alternate quantum vacuum prescriptions we show some choices of vacua are theoretically observationally distinguishable from others. Comparing these predictions to the Planck 2018 observations via a Bayesian analysis shows no significant evidence to favour any of the quantum vacuum prescriptions over the others. In addition we consider frozen initial conditions, representing a white-noise initial state at the big-bang singularity. Under certain assumptions the cosmological concordance model and frozen initial conditions are found to produce identical predictions for the cosmic microwave background anisotropies. Frozen initial conditions may thus provide an alternative theoretic paradigm to explain observations that were previously understood in terms of the inflation of a quantum vacuum.
△ Less
Submitted 7 September, 2021; v1 submitted 7 April, 2021;
originally announced April 2021.
-
Perturbations and the Future Conformal Boundary
Authors:
A. N. Lasenby,
W. J. Handley,
D. J. Bartlett,
C. S. Negreanu
Abstract:
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…
▽ More
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.
△ Less
Submitted 20 April, 2022; v1 submitted 5 April, 2021;
originally announced April 2021.
-
Improved cosmological fits with quantized primordial power spectra
Authors:
D. J. Bartlett,
W. J. Handley,
A. N. Lasenby
Abstract:
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…
▽ More
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.
△ Less
Submitted 20 April, 2022; v1 submitted 5 April, 2021;
originally announced April 2021.
-
Nested sampling with any prior you like
Authors:
Justin Alsing,
Will Handley
Abstract:
Nested sampling is an important tool for conducting Bayesian analysis in Astronomy and other fields, both for sampling complicated posterior distributions for parameter inference, and for computing marginal likelihoods for model comparison. One technical obstacle to using nested sampling in practice is the requirement (for most common implementations) that prior distributions be provided in the fo…
▽ More
Nested sampling is an important tool for conducting Bayesian analysis in Astronomy and other fields, both for sampling complicated posterior distributions for parameter inference, and for computing marginal likelihoods for model comparison. One technical obstacle to using nested sampling in practice is the requirement (for most common implementations) that prior distributions be provided in the form of transformations from the unit hyper-cube to the target prior density. For many applications - particularly when using the posterior from one experiment as the prior for another - such a transformation is not readily available. In this letter we show that parametric bijectors trained on samples from a desired prior density provide a general-purpose method for constructing transformations from the uniform base density to a target prior, enabling the practical use of nested sampling under arbitrary priors. We demonstrate the use of trained bijectors in conjunction with nested sampling on a number of examples from cosmology.
△ Less
Submitted 28 June, 2021; v1 submitted 24 February, 2021;
originally announced February 2021.
-
Bayesian evidence for the tensor-to-scalar ratio $r$ and neutrino masses $m_ν$: Effects of uniform vs logarithmic priors
Authors:
Lukas T. Hergt,
Will J. Handley,
Michael P. Hobson,
Anthony N. Lasenby
Abstract:
We review the effect that the choice of a uniform or logarithmic prior has on the Bayesian evidence and hence on Bayesian model comparisons when data provide only a one-sided bound on a parameter. We investigate two particular examples: the tensor-to-scalar ratio $r$ of primordial perturbations and the mass of individual neutrinos $m_ν$, using the cosmic microwave background temperature and polari…
▽ More
We review the effect that the choice of a uniform or logarithmic prior has on the Bayesian evidence and hence on Bayesian model comparisons when data provide only a one-sided bound on a parameter. We investigate two particular examples: the tensor-to-scalar ratio $r$ of primordial perturbations and the mass of individual neutrinos $m_ν$, using the cosmic microwave background temperature and polarisation data from Planck 2018 and the NuFIT 5.0 data from neutrino oscillation experiments. We argue that the Kullback-Leibler divergence, also called the relative entropy, mathematically quantifies the Occam penalty. We further show how the Bayesian evidence stays invariant upon changing the lower prior bound of an upper constrained parameter. While a uniform prior on the tensor-to-scalar ratio disfavours the $r$-extension compared to the base LCDM model with odds of about 1:20, switching to a logarithmic prior renders both models essentially equally likely. LCDM with a single massive neutrino is favoured over an extension with variable neutrino masses with odds of 20:1 in case of a uniform prior on the lightest neutrino mass, which decreases to roughly 2:1 for a logarithmic prior. For both prior options we get only a very slight preference for the normal over the inverted neutrino hierarchy with Bayesian odds of about 3:2 at most.
△ Less
Submitted 23 February, 2021;
originally announced February 2021.
-
When tension is just a fluctuation: How noisy data affect model comparison
Authors:
B. Joachimi,
F. Köhlinger,
W. Handley,
P. Lemos
Abstract:
Summary statistics of the likelihood, such as the Bayesian evidence, offer a principled way of comparing models and assessing tension between, or within, the results of physical experiments. Noisy realisations of the data induce scatter in these model comparison statistics. For a realistic case of cosmological inference from large-scale structure we show that the logarithm of the Bayes factor atta…
▽ More
Summary statistics of the likelihood, such as the Bayesian evidence, offer a principled way of comparing models and assessing tension between, or within, the results of physical experiments. Noisy realisations of the data induce scatter in these model comparison statistics. For a realistic case of cosmological inference from large-scale structure we show that the logarithm of the Bayes factor attains scatter of order unity, increasing significantly with stronger tension between the models under comparison. We develop an approximate procedure that quantifies the sampling distribution of the evidence at small additional computational cost and apply it to real data to demonstrate the impact of the scatter, which acts to reduce the significance of any model discrepancies. Data compression is highlighted as a potential avenue to suppressing noise in the evidence to negligible levels, with a proof of concept on Planck cosmic microwave background data.
△ Less
Submitted 29 April, 2021; v1 submitted 18 February, 2021;
originally announced February 2021.
-
The HARPS search for southern extra-solar planets XLV. Two Neptune mass planets orbiting HD 13808: a study of stellar activity modelling's impact on planet detection
Authors:
E. Ahrer,
D. Queloz,
V. M. Rajpaul,
D. Ségransan,
F. Bouchy,
R. Hall,
W. Handley,
C. Lovis,
M. Mayor,
A. Mortier,
F. Pepe,
S. Thompson,
S. Udry,
N. Unger
Abstract:
We present a comprehensive analysis of 10 years of HARPS radial velocities of the K2V dwarf star HD 13808, which has previously been reported to host two unconfirmed planet candidates. We use the state-of-the-art nested sampling algorithm PolyChord to compare a wide variety of stellar activity models, including simple models exploiting linear correlations between RVs and stellar activity indicator…
▽ More
We present a comprehensive analysis of 10 years of HARPS radial velocities of the K2V dwarf star HD 13808, which has previously been reported to host two unconfirmed planet candidates. We use the state-of-the-art nested sampling algorithm PolyChord to compare a wide variety of stellar activity models, including simple models exploiting linear correlations between RVs and stellar activity indicators, harmonic models for the activity signals, and a more sophisticated Gaussian process regression model. We show that the use of overly-simplistic stellar activity models that are not well-motivated physically can lead to spurious `detections' of planetary signals that are almost certainly not real. We also reveal some difficulties inherent in parameter and model inference in cases where multiple planetary signals may be present. Our study thus underlines the importance both of exploring a variety of competing models and of understanding the limitations and precision settings of one's sampling algorithm. We also show that at least in the case of HD 13808, we always arrive at consistent conclusions about two particular signals present in the RV, regardless of the stellar activity model we adopt; these two signals correspond to the previously-reported though unconfirmed planet candidate signals. Given the robustness and precision with which we can characterize these two signals, we deem them secure planet detections. In particular, we find two planets orbiting HD 13808 at distances of 0.11, 0.26 AU with periods of 14.2, 53.8 d, and minimum masses of 11, 10 Earth masses.
△ Less
Submitted 5 February, 2021;
originally announced February 2021.
-
Nonlinear Hamiltonian analysis of new quadratic torsion theories Part I. Cases with curvature-free constraints
Authors:
W. E. V. Barker,
A. N. Lasenby,
M. P. Hobson,
W. J. Handley
Abstract:
It was recently found that, when linearised in the absence of matter, 58 cases of the general gravitational theory with quadratic curvature and torsion are (i) free from ghosts and tachyons and (ii) power-counting renormalisable. We inspect the nonlinear Hamiltonian structure of the eight cases whose primary constraints do not depend on the curvature tensor. We confirm the particle spectra and uni…
▽ More
It was recently found that, when linearised in the absence of matter, 58 cases of the general gravitational theory with quadratic curvature and torsion are (i) free from ghosts and tachyons and (ii) power-counting renormalisable. We inspect the nonlinear Hamiltonian structure of the eight cases whose primary constraints do not depend on the curvature tensor. We confirm the particle spectra and unitarity of all these theories in the linear regime. We uncover qualitative dynamical changes in the nonlinear regimes of all eight cases, suggesting at least a broken gauge symmetry, and possibly the activation of negative kinetic energy spin-parity sectors and acausal behaviour. Two of the cases propagate a pair of massless modes at the linear level, and were interesting as candidate theories of gravity. However, we identify these modes with vector excitations, rather than the tensor polarisations of the graviton. Moreover, we show that these theories do not support a viable cosmological background.
△ Less
Submitted 7 January, 2021;
originally announced January 2021.
-
Simple and statistically sound recommendations for analysing physical theories
Authors:
Shehu S. AbdusSalam,
Fruzsina J. Agocs,
Benjamin C. Allanach,
Peter Athron,
Csaba Balázs,
Emanuele Bagnaschi,
Philip Bechtle,
Oliver Buchmueller,
Ankit Beniwal,
Jihyun Bhom,
Sanjay Bloor,
Torsten Bringmann,
Andy Buckley,
Anja Butter,
José Eliel Camargo-Molina,
Marcin Chrzaszcz,
Jan Conrad,
Jonathan M. Cornell,
Matthias Danninger,
Jorge de Blas,
Albert De Roeck,
Klaus Desch,
Matthew Dolan,
Herbert Dreiner,
Otto Eberhardt
, et al. (50 additional authors not shown)
Abstract:
Physical theories that depend on many parameters or are tested against data from many different experiments pose unique challenges to statistical inference. Many models in particle physics, astrophysics and cosmology fall into one or both of these categories. These issues are often sidestepped with statistically unsound ad hoc methods, involving intersection of parameter intervals estimated by mul…
▽ More
Physical theories that depend on many parameters or are tested against data from many different experiments pose unique challenges to statistical inference. Many models in particle physics, astrophysics and cosmology fall into one or both of these categories. These issues are often sidestepped with statistically unsound ad hoc methods, involving intersection of parameter intervals estimated by multiple experiments, and random or grid sampling of model parameters. Whilst these methods are easy to apply, they exhibit pathologies even in low-dimensional parameter spaces, and quickly become problematic to use and interpret in higher dimensions. In this article we give clear guidance for going beyond these procedures, suggesting where possible simple methods for performing statistically sound inference, and recommendations of readily-available software tools and standards that can assist in doing so. Our aim is to provide any physicists lacking comprehensive statistical training with recommendations for reaching correct scientific conclusions, with only a modest increase in analysis burden. Our examples can be reproduced with the code publicly available at https://doi.org/10.5281/zenodo.4322283.
△ Less
Submitted 11 April, 2022; v1 submitted 17 December, 2020;
originally announced December 2020.
-
Bayesian noise wave calibration for 21-cm global experiments
Authors:
I. L. V. Roque,
W. J. Handley,
N. Razavi-Ghods
Abstract:
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 mac…
▽ More
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.
△ Less
Submitted 1 June, 2021; v1 submitted 27 November, 2020;
originally announced November 2020.
-
Quantifying Ionospheric Effects on Global 21-cm Observations
Authors:
Emma Shen,
Dominic Anstey,
Eloy de Lera Acedo,
Anastasia Fialkov,
Will Handley
Abstract:
We modelled the two major layer of Earth's ionosphere, the F-layer and the D-layer, by a simplified spatial model with temporal variance to study the chromatic ionospheric effects on global 21-cm observations. From the analyses, we found that the magnitude of the ionospheric disruptions due to ionospheric refraction and absorption can be greater than the expected global 21-cm signal, and the varia…
▽ More
We modelled the two major layer of Earth's ionosphere, the F-layer and the D-layer, by a simplified spatial model with temporal variance to study the chromatic ionospheric effects on global 21-cm observations. From the analyses, we found that the magnitude of the ionospheric disruptions due to ionospheric refraction and absorption can be greater than the expected global 21-cm signal, and the variation of its magnitude can differ, depending on the ionospheric conditions. Within the parameter space adopted in the model, the shape of the global 21-cm signal is distorted after propagating through the ionosphere, while its amplitude is weakened. It is observed that the ionospheric effects do not cancel out over time, and thus should be accounted for in the foreground calibration at each timestep to account for the chromaticity introduced by the ionosphere.
△ Less
Submitted 20 November, 2020;
originally announced November 2020.
-
Nested sampling with plateaus
Authors:
Andrew Fowlie,
Will Handley,
Liangliang Su
Abstract:
It was recently emphasised by Riley (2019); Schittenhelm & Wacker (2020) that that in the presence of plateaus in the likelihood function nested sampling (NS) produces faulty estimates of the evidence and posterior densities. After informally explaining the cause of the problem, we present a modified version of NS that handles plateaus and can be applied retrospectively to NS runs from popular NS…
▽ More
It was recently emphasised by Riley (2019); Schittenhelm & Wacker (2020) that that in the presence of plateaus in the likelihood function nested sampling (NS) produces faulty estimates of the evidence and posterior densities. After informally explaining the cause of the problem, we present a modified version of NS that handles plateaus and can be applied retrospectively to NS runs from popular NS software using anesthetic. In the modified NS, live points in a plateau are evicted one by one without replacement, with ordinary NS compression of the prior volume after each eviction but taking into account the dynamic number of live points. The live points are replenished once all points in the plateau are removed. We demonstrate it on a number of examples. Since the modification is simple, we propose that it becomes the canonical version of Skilling's NS algorithm.
△ Less
Submitted 24 February, 2021; v1 submitted 26 October, 2020;
originally announced October 2020.