Normalized to: Javid, K.
[1]
oai:arXiv.org:1909.00029 [pdf] - 2047557
Physical modelling of galaxy clusters and Bayesian inference in
astrophysics
Submitted: 2019-08-30, last modified: 2020-02-11
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.
[2]
oai:arXiv.org:2002.04123 [pdf] - 2046555
Geometric nested sampling: sampling from distributions defined on
non-trivial geometries
Submitted: 2020-02-10
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.
[3]
oai:arXiv.org:1805.01991 [pdf] - 2045135
Physical modelling of galaxy clusters detected by Planck
Submitted: 2018-05-04, last modified: 2020-02-09
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 $\theta$ - integrated
Comptonisation parameter $Y$ posterior distributions, and the other uses a
`slicing function' to provide information on $\theta$ 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.
[4]
oai:arXiv.org:1901.09980 [pdf] - 2045153
Sunyaev--Zel'dovich profile fitting with joint AMI-Planck analysis
Submitted: 2019-01-28, last modified: 2020-02-09
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.
[5]
oai:arXiv.org:1905.09110 [pdf] - 2045172
Nested sampling on non-trivial geometries
Submitted: 2019-05-22, last modified: 2020-02-09
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}.
[6]
oai:arXiv.org:1809.03325 [pdf] - 2045142
Physical modelling of galaxy clusters using Einasto dark matter profiles
Submitted: 2018-09-10, last modified: 2020-02-09
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\sigma$ 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.
[7]
oai:arXiv.org:1805.01968 [pdf] - 2045134
Comparison of physical and observational galaxy cluster modelling
Submitted: 2018-05-04, last modified: 2020-02-09
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 $\theta$ and
integrated Comptonisation $Y$, derived from simulations, while OM II uses
separable priors on $\theta$ 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 $\theta$ - $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.
[8]
oai:arXiv.org:1707.04237 [pdf] - 1636634
A Digital Correlator Upgrade for the Arcminute MicroKelvin Imager
Hickish, Jack;
Razavi-Ghods, Nima;
Perrott, Yvette C.;
Titterington, David J.;
Carey, Steve H.;
Scott, Paul F.;
Grainge, Keith J. B.;
Scaife, Anna M. M.;
Alexander, Paul;
Saunders, Richard D. E.;
Crofts, Mike;
Javid, Kamran;
Rumsey, Clare;
Jin, Terry Z.;
Ely, John A.;
Shaw, Clive;
Northrop, Ian G.;
Pooley, Guy;
D'Alessandro, Robert;
Doherty, Peter;
Willatt, Greg P.
Submitted: 2017-07-13, last modified: 2018-02-16
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.
[9]
oai:arXiv.org:1709.05136 [pdf] - 1588429
AMI-CL J0300+2613: a Galactic anomalous-microwave-emission ring
masquerading as a galaxy cluster
Perrott, Yvette C.;
Cantwell, Therese M.;
Carey, Steve H.;
Elwood, Patrick J.;
Feroz, Farhan;
Grainge, Keith J. B.;
Green, David A.;
Hobson, Michael P.;
Javid, Kamran;
Jin, Terry Z.;
Pooley, Guy G.;
Razavi-Ghods, Nima;
Rumsey, Clare;
Saunders, Richard D. E.;
Scaife, Anna M. M.;
Schammel, Michel P.;
Scott, Paul F.;
Shimwell, Timothy W.;
Titterington, David J.;
Waldram, Elizabeth M.
Submitted: 2017-09-15
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.