-
A 15.5 GHz detection of the galaxy cluster minihalo in RXJ1720.1+2638
Authors:
Yvette C. Perrott,
Pedro Carvalho,
Patrick J. Elwood,
Keith J. B. Grainge,
David A. Green,
Kamran Javid,
Terry Z. Jin,
Clare Rumsey,
Richard D. E. Saunders
Abstract:
RXJ1720.1+2638 is a cool-core, 'relaxed-appearing' cluster with a minihalo previously detected up to 8.4 GHz, confined by X-ray-detected cold fronts. We present observations of the minihalo at 13 - 18 GHz with the Arcminute Microkelvin Imager telescope, simultaneously modelling the Sunyaev-Zel'dovich signal of the cluster in conjunction with Planck and Chandra data in order to disentangle the non-…
▽ More
RXJ1720.1+2638 is a cool-core, 'relaxed-appearing' cluster with a minihalo previously detected up to 8.4 GHz, confined by X-ray-detected cold fronts. We present observations of the minihalo at 13 - 18 GHz with the Arcminute Microkelvin Imager telescope, simultaneously modelling the Sunyaev-Zel'dovich signal of the cluster in conjunction with Planck and Chandra data in order to disentangle the non-thermal emission of the minihalo. We show that the previously-reported steepening of the minihalo emission at 8.4 GHz is not supported by the AMI data and that the spectrum is consistent with a single power-law up to 18 GHz. We also show the presence of a larger-scale component of the minihalo extending beyond the cold fronts. Both of these observations could be explained by the 'hadronic' or 'secondary' mechanism for the production of relativistic electrons, rather than the currently-favoured 're-acceleration' mechanism and/or multiple episodes of jet activity from the active galactic nucleus in the brightest cluster galaxy.
△ Less
Submitted 16 September, 2021;
originally announced September 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.
-
Compromise-free Bayesian neural networks
Authors:
Kamran Javid,
Will Handley,
Mike Hobson,
Anthony Lasenby
Abstract:
We conduct a thorough analysis of the relationship between the out-of-sample performance and the Bayesian evidence (marginal likelihood) of Bayesian neural networks (BNNs), as well as looking at the performance of ensembles of BNNs, both using the Boston housing dataset. Using the state-of-the-art in nested sampling, we numerically sample the full (non-Gaussian and multimodal) network posterior an…
▽ More
We conduct a thorough analysis of the relationship between the out-of-sample performance and the Bayesian evidence (marginal likelihood) of Bayesian neural networks (BNNs), as well as looking at the performance of ensembles of BNNs, both using the Boston housing dataset. Using the state-of-the-art in nested sampling, we numerically sample the full (non-Gaussian and multimodal) network posterior and obtain numerical estimates of the Bayesian evidence, considering network models with up to 156 trainable parameters. The networks have between zero and four hidden layers, either $\tanh$ or $ReLU$ activation functions, and with and without hierarchical priors. The ensembles of BNNs are obtained by determining the posterior distribution over networks, from the posterior samples of individual BNNs re-weighted by the associated Bayesian evidence values. There is good correlation between out-of-sample performance and evidence, as well as a remarkable symmetry between the evidence versus model size and out-of-sample performance versus model size planes. Networks with $ReLU$ activation functions have consistently higher evidences than those with $\tanh$ functions, and this is reflected in their out-of-sample performance. Ensembling over architectures acts to further improve performance relative to the individual BNNs.
△ Less
Submitted 13 June, 2020; v1 submitted 25 April, 2020;
originally announced April 2020.
-
Geometric nested sampling: sampling from distributions defined on non-trivial geometries
Authors:
Kamran Javid
Abstract:
Metropolis Hastings nested sampling evolves a Markov chain, accepting new points along the chain according to a version of the Metropolis Hastings acceptance ratio, which has been modified to satisfy the nested sampling likelihood constraint. The geometric nested sampling algorithm I present here is based on the Metropolis Hastings method, but treats parameters as though they represent points on c…
▽ More
Metropolis Hastings nested sampling evolves a Markov chain, accepting new points along the chain according to a version of the Metropolis Hastings acceptance ratio, which has been modified to satisfy the nested sampling likelihood constraint. The geometric nested sampling algorithm I present here is based on the Metropolis Hastings method, but treats parameters as though they represent points on certain geometric objects, namely circles, tori and spheres. For parameters which represent points on a circle or torus, the trial distribution is "wrapped" around the domain of the posterior distribution such that samples cannot be rejected automatically when evaluating the Metropolis ratio due to being outside the sampling domain. Furthermore, this enhances the mobility of the sampler. For parameters which represent coordinates on the surface of a sphere, the algorithm transforms the parameters into a Cartesian coordinate system before sampling which again makes sure no samples are automatically rejected, and provides a physically intuitive way of the sampling the parameter space.
△ Less
Submitted 10 February, 2020;
originally announced February 2020.
-
Physical modelling of galaxy clusters and Bayesian inference in astrophysics
Authors:
Kamran Javid
Abstract:
I compare the mass values obtained with data taken from the Arcminute Microkelvin Imager (AMI) radio interferometer system and from the Planck satellite. The former of these uses a Bayesian analysis pipeline that parameterises a cluster in terms of its physical quantities, and models the dark matter \& baryonic components of a cluster using Navarro-Frenk-White (NFW) and generalised-NFW profiles re…
▽ More
I compare the mass values obtained with data taken from the Arcminute Microkelvin Imager (AMI) radio interferometer system and from the Planck satellite. The former of these uses a Bayesian analysis pipeline that parameterises a cluster in terms of its physical quantities, and models the dark matter \& baryonic components of a cluster using Navarro-Frenk-White (NFW) and generalised-NFW profiles respectively. I also analyse simulated AMI data with input values based on PwS mass estimates. I then compare three cluster models using AMI data for the 54 cluster sample. The two observational models considered only model the gas content of the cluster. To compare the physical and observational models I consider their posterior parameter estimates, including the calculation of a metric defined between two probability distributions. The models' fit to the cluster data is evaluated by looking at the Bayesian evidence values. Improvements to the physical modelling of galaxy clusters are then considered, either by relaxing some of the assumptions underlying the physical model, or by introducing a new profile for the dark matter component of clusters. The final part of the cluster analysis work focuses on Bayesian analysis using a joint likelihood function of data from both AMI and the Planck satellite simultaneously. Finally, a new Bayesian inference algorithm based on nested sampling is presented. The algorithm, named the "geometric nested sampler", is an adaption of the Metropolis-Hastings nested sampler and makes use of the geometrical interpretation of sets of parameters to sample from their domains efficiently. The geometric nested sampler is tested on several toy models as well as a model representing the emission of gravitational waves from binary black hole mergers.
△ Less
Submitted 11 February, 2020; v1 submitted 30 August, 2019;
originally announced September 2019.
-
Nested sampling on non-trivial geometries
Authors:
Kamran Javid
Abstract:
Metropolis nested sampling evolves a Markov chain from a current livepoint and accepts new points along the chain according to a version of the Metropolis acceptance ratio modified to satisfy the likelihood constraint, characteristic of nested sampling algorithms. The geometric nested sampling algorithm we present here is a based on the Metropolis method, but treats parameters as though they repre…
▽ More
Metropolis nested sampling evolves a Markov chain from a current livepoint and accepts new points along the chain according to a version of the Metropolis acceptance ratio modified to satisfy the likelihood constraint, characteristic of nested sampling algorithms. The geometric nested sampling algorithm we present here is a based on the Metropolis method, but treats parameters as though they represent points on certain geometric objects, namely circles, tori and spheres. For parameters which represent points on a circle or torus, the trial distribution is `wrapped' around the domain of the posterior distribution such that samples cannot be rejected automatically when evaluating the Metropolis ratio due to being outside the sampling domain. Furthermore, this enhances the mobility of the sampler. For parameters which represent coordinates on the surface of a sphere, the algorithm transforms the parameters into a Cartesian coordinate system before sampling which again makes sure no samples are automatically rejected, and provides a physically intutive way of the sampling the parameter space.
We apply the geometric nested sampler to two types of toy model which include circular, toroidal and spherical parameters. We find that the geometric nested sampler generally outperforms \textsc{MultiNest} in both cases. \\ %We also apply the algorithm to a gravitational wave detection model which includes circular and spherical parameters, and find that the geometric nested sampler and \textsc{MultiNest} appear to perform equally well as one another. Our implementation of the algorithm can be found at \url{https://github.com/SuperKam91/nested_sampling}.
△ Less
Submitted 9 February, 2020; v1 submitted 22 May, 2019;
originally announced May 2019.
-
Sunyaev--Zel'dovich profile fitting with joint AMI-Planck analysis
Authors:
Yvette C. Perrott,
Kamran Javid,
Pedro Carvalho,
Patrick J. Elwood,
Michael P. Hobson,
Anthony N. Lasenby,
Malak Olamaie,
Richard D. E. Saunders
Abstract:
We develop a Bayesian method of analysing Sunyaev-Zel'dovich measurements of galaxy clusters obtained from the Arcminute Microkelvin Imager (AMI) radio interferometer system and from the Planck satellite, using a joint likelihood function for the data from both instruments. Our method is applicable to any combination of Planck data with interferometric data from one or more arrays. We apply the an…
▽ More
We develop a Bayesian method of analysing Sunyaev-Zel'dovich measurements of galaxy clusters obtained from the Arcminute Microkelvin Imager (AMI) radio interferometer system and from the Planck satellite, using a joint likelihood function for the data from both instruments. Our method is applicable to any combination of Planck data with interferometric data from one or more arrays. We apply the analysis to simulated clusters and find that when the cluster pressure profile is known a-priori, the joint dataset provides precise and accurate constraints on the cluster parameters, removing the need for external information to reduce the parameter degeneracy. When the pressure profile deviates from that assumed for the fit, the constraints become biased. Allowing the pressure profile shape parameters to vary in the analysis allows an unbiased recovery of the integrated cluster signal and produces constraints on some shape parameters, depending on the angular size of the cluster. When applied to real data from Planck-detected cluster PSZ2 G063.80+11.42, our method resolves the discrepancy between the AMI and Planck $Y$-estimates and usefully constrains the gas pressure profile shape parameters at intermediate and large radii.
△ Less
Submitted 9 February, 2020; v1 submitted 28 January, 2019;
originally announced January 2019.
-
Physical modelling of galaxy clusters using Einasto dark matter profiles
Authors:
Kamran Javid,
Yvette C. Perrott,
Clare Rumsey,
Richard D. E. Saunders
Abstract:
We derive a model for Sunyaev--Zel'dovich data from a galaxy cluster which uses an Einasto profile to model the cluster's dark matter component. This model is similar to the physical models for clusters previously used by the Arcminute Microkelvin Imager (AMI) consortium, which model the dark matter using a Navarro-Frenk-White (NFW) profile, but the Einasto profile provides an extra degree of free…
▽ More
We derive a model for Sunyaev--Zel'dovich data from a galaxy cluster which uses an Einasto profile to model the cluster's dark matter component. This model is similar to the physical models for clusters previously used by the Arcminute Microkelvin Imager (AMI) consortium, which model the dark matter using a Navarro-Frenk-White (NFW) profile, but the Einasto profile provides an extra degree of freedom. We thus present a comparison between two physical models which differ only in the way they model dark matter: one which uses an NFW profile (PM I) and one that uses an Einasto profile (PM II). We illustrate the differences between the models by plotting physical properties of clusters as a function of cluster radius. We generate AMI simulations of clusters which are \textit{created} and \textit{analysed} with both models. From this we find that for 14 of the 16 simulations, the Bayesian evidence gives no preference to either of the models according to the Jeffreys scale, and for the other two simulations, weak preference in favour of the correct model. However, for the mass estimates obtained from the analyses, the values were within $1σ$ of the input values for 14 out of 16 of the clusters when using the correct model, but only in 6 out of 16 cases when the incorrect model was used to analyse the data. Finally we apply the models to real data from cluster A611 obtained with AMI, and find the mass estimates to be consistent with one another except in the case of when PM II is applied using an extreme value for the Einasto shape parameter.
△ Less
Submitted 9 February, 2020; v1 submitted 10 September, 2018;
originally announced September 2018.
-
Physical modelling of galaxy clusters detected by Planck
Authors:
Kamran Javid,
Malak Olamaie,
Yvette C. Perrott,
Pedro Carvalho,
Keith J. B. Grainge,
Michael P. Hobson,
Clare Rumsey,
Richard D. E. Saunders
Abstract:
We present a comparison of mass estimates for $54$ galaxy cluster candidates from the second Planck catalogue (PSZ2) of Sunyaev-Zel'dovich sources. We compare the mass values obtained with data taken from the Arcminute Microkelvin Imager (AMI) radio interferometer system and from the Planck satellite. The former of these uses a Bayesian analysis pipeline that parameterises a cluster in terms of it…
▽ More
We present a comparison of mass estimates for $54$ galaxy cluster candidates from the second Planck catalogue (PSZ2) of Sunyaev-Zel'dovich sources. We compare the mass values obtained with data taken from the Arcminute Microkelvin Imager (AMI) radio interferometer system and from the Planck satellite. The former of these uses a Bayesian analysis pipeline that parameterises a cluster in terms of its physical quantities, and models the dark matter & baryonic components of a cluster using NFW and GNFW profiles respectively. Our mass estimates derived from Planck data are obtained from the results of the Bayesian detection algorithm PowellSnakes (PwS), are based on the methodology detailed in the PSZ2 paper, and produce two sets of mass estimates; one estimate is calculated directly from the angular radius $θ$ - integrated Comptonisation parameter $Y$ posterior distributions, and the other uses a `slicing function' to provide information on $θ$ based on X-ray measurements and previous Planck mission samples. We find that for $37$ of the clusters, the AMI mass estimates are lower than both values obtained from Planck data. However the AMI and slicing function estimates are within one combined standard deviation of each other for $31$ clusters. We also generate cluster simulations based on the slicing-function mass estimates, and analyse them in the same way as we did the real AMI data. We find that inclusion in the simulations of radio-source confusion & CMB noise and measurable radio-sources causes AMI mass estimates to be systematically low.
△ Less
Submitted 9 February, 2020; v1 submitted 4 May, 2018;
originally announced May 2018.
-
Comparison of physical and observational galaxy cluster modelling
Authors:
Kamran Javid,
Yvette C. Perrott,
Michael P. Hobson,
Malak Olamaie,
Clare Rumsey,
Richard D. E. Saunders
Abstract:
We present a comparison between three cluster models applied to data obtained by the Arcminute Microkelvin Imager radio interferometer system. The physical model (PM) parameterises a cluster in terms of its physical quantities to model the dark matter and baryonic components of the cluster using NFW and GNFW profiles respectively. The observational models (OM I and OM II) model only the gas conten…
▽ More
We present a comparison between three cluster models applied to data obtained by the Arcminute Microkelvin Imager radio interferometer system. The physical model (PM) parameterises a cluster in terms of its physical quantities to model the dark matter and baryonic components of the cluster using NFW and GNFW profiles respectively. The observational models (OM I and OM II) model only the gas content of the cluster. The two OMs vary only in the priors they use in Bayesian inference: OM I has a joint prior on angular radius $θ$ and integrated Comptonisation $Y$, derived from simulations, while OM II uses separable priors on $θ$ and $Y$ which are based on calculations of the physical model. For the comparison we consider a sample of $54$ clusters which are a subsample of the second Planck catalogue of Sunyaev-Zel'dovich sources. We first compare the $Y$ estimates of the three models, and find that the PM generally yields lower estimates relative to the OMs. We then compute the Earth Mover's Distance between the $θ$ - $Y$ posterior distributions obtained from each model for each cluster, and find that the two models which are most discrepant are PM and OM I. Finally, we compare the Bayesian evidence values obtained from each model for each cluster. OM I generally provides the best fit to the data but not at a statistically significant level, according to the Jeffreys scale. The highest evidence ratio obtained is actually in favour of the PM over OM I.
△ Less
Submitted 9 February, 2020; v1 submitted 4 May, 2018;
originally announced May 2018.
-
AMI-CL J0300+2613: a Galactic anomalous-microwave-emission ring masquerading as a galaxy cluster
Authors:
Yvette C. Perrott,
Therese M. Cantwell,
Steve H. Carey,
Patrick J. Elwood,
Farhan Feroz,
Keith J. B. Grainge,
David A. Green,
Michael P. Hobson,
Kamran Javid,
Terry Z. Jin,
Guy G. Pooley,
Nima Razavi-Ghods,
Clare Rumsey,
Richard D. E. Saunders,
Anna M. M. Scaife,
Michel P. Schammel,
Paul F. Scott,
Timothy W. Shimwell,
David J. Titterington,
Elizabeth M. Waldram
Abstract:
The Arcminute Microkelvin Imager (AMI) carried out a blind survey for galaxy clusters via their Sunyaev-Zel'dovich effect decrements between 2008 and 2011. The first detection, known as AMI-CL J0300+2613, has been reobserved with AMI equipped with a new digital correlator with high dynamic range. The combination of the new AMI data and more recent high-resolution sub-mm and infra-red maps now show…
▽ More
The Arcminute Microkelvin Imager (AMI) carried out a blind survey for galaxy clusters via their Sunyaev-Zel'dovich effect decrements between 2008 and 2011. The first detection, known as AMI-CL J0300+2613, has been reobserved with AMI equipped with a new digital correlator with high dynamic range. The combination of the new AMI data and more recent high-resolution sub-mm and infra-red maps now shows the feature in fact to be a ring of positive dust-correlated Galactic emission, which is likely to be anomalous microwave emission (AME). If so, this is the first completely blind detection of AME at arcminute scales.
△ Less
Submitted 15 September, 2017;
originally announced September 2017.
-
A Digital Correlator Upgrade for the Arcminute MicroKelvin Imager
Authors:
Jack Hickish,
Nima Razavi-Ghods,
Yvette C. Perrott,
David J. Titterington,
Steve H. Carey,
Paul F. Scott,
Keith J. B. Grainge,
Anna M. M. Scaife,
Paul Alexander,
Richard D. E. Saunders,
Mike Crofts,
Kamran Javid,
Clare Rumsey,
Terry Z. Jin,
John A. Ely,
Clive Shaw,
Ian G. Northrop,
Guy Pooley,
Robert D'Alessandro,
Peter Doherty,
Greg P. Willatt
Abstract:
The Arcminute Microkelvin Imager (AMI) telescopes located at the Mullard Radio Astronomy Observatory near Cambridge have been significantly enhanced by the implementation of a new digital correlator with 1.2 MHz spectral resolution. This system has replaced a 750-MHz resolution analogue lag-based correlator, and was designed to mitigate the effects of radio frequency interference, particularly fro…
▽ More
The Arcminute Microkelvin Imager (AMI) telescopes located at the Mullard Radio Astronomy Observatory near Cambridge have been significantly enhanced by the implementation of a new digital correlator with 1.2 MHz spectral resolution. This system has replaced a 750-MHz resolution analogue lag-based correlator, and was designed to mitigate the effects of radio frequency interference, particularly from geostationary satellites that contaminate observations at low declinations. The upgraded instrument consists of 18 ROACH2 Field Programmable Gate Array platforms used to implement a pair of real-time FX correlators -- one for each of AMI's two arrays. The new system separates the down-converted RF baseband signal from each AMI receiver into two 2.3 GHz-wide sub-bands which are each digitized at 5-Gsps with 8 bits of precision. These digital data streams are filtered into 2048 frequency channels and cross-correlated using FPGA hardware, with a commercial 10 Gb Ethernet switch providing high-speed data interconnect. Images formed using data from the new digital correlator show over an order of magnitude improvement in dynamic range over the previous system. The ability to observe at low declinations has also been significantly improved.
△ Less
Submitted 16 February, 2018; v1 submitted 13 July, 2017;
originally announced July 2017.