Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userPhysical modelling of galaxy clusters and Bayesian inference in astrophysics
97
0
0.0
(
0
)
Authors
Kamran Javid
Ask ChatGPT about the research
No Arabic 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 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.
rate research
Read More
We present a large-scale Bayesian inference framework to constrain cosmological parameters using galaxy redshift surveys, via an application of the Alcock-Paczynski (AP) test. Our physical model of the non-linearly evolved density field, as probed by galaxy surveys, employs Lagrangian perturbation theory (LPT) to connect Gaussian initial conditions to the final density field, followed by a coordinate transformation to obtain the redshift space representation for comparison with data. We generate realizations of primordial and present-day matter fluctuations given a set of observations. This hierarchical approach encodes a novel AP test, extracting several orders of magnitude more information from the cosmological expansion compared to classical approaches, to infer cosmological parameters and jointly reconstruct the underlying 3D dark matter density field. The novelty of this AP test lies in constraining the comoving-redshift transformation to infer the appropriate cosmology which yields isotropic correlations of the galaxy density field, with the underlying assumption relying purely on the cosmological principle. Such an AP test does not rely explicitly on modelling the full statistics of the field. We verify in depth via simulations that this renders our test robust to model misspecification. This leads to another crucial advantage, namely that the cosmological parameters exhibit extremely weak dependence on the currently unresolved phenomenon of galaxy bias, thereby circumventing a potentially key limitation. This is consequently among the first methods to extract a large fraction of information from statistics other than that of direct density contrast correlations, without being sensitive to the amplitude of density fluctuations. We perform several statistical efficiency and consistency tests on a mock galaxy catalogue, using the SDSS-III survey as template.
We derive a model for Sunyaev--Zeldovich data from a galaxy cluster which uses an Einasto profile to model the clusters 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 $1sigma$ 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.
A new method is presented for modelling the physical properties of galaxy clusters. Our technique moves away from the traditional approach of assuming specific parameterised functional forms for the variation of physical quantities within the cluster, and instead allows for a free-form reconstruction, but one for which the level of complexity is determined automatically by the observational data and may depend on position within the cluster. This is achieved by representing each independent cluster property as some interpolating or approximating function that is specified by a set of control points, or nodes, for which the number of nodes, together with their positions and amplitudes, are allowed to vary and are inferred in a Bayesian manner from the data. We illustrate our nodal approach in the case of a spherical cluster by modelling the electron pressure profile Pe(r) in analyses both of simulated Sunyaev-Zeldovich (SZ) data from the Arcminute MicroKelvin Imager (AMI) and of real AMI observations of the cluster MACS J0744+3927 in the CLASH sample. We demonstrate that one may indeed determine the complexity supported by the data in the reconstructed Pe(r), and that one may constrain two very important quantities in such an analysis: the cluster total volume integrated Comptonisation parameter (Ytot) and the extent of the gas distribution in the cluster (rmax). The approach is also well-suited to detecting clusters in blind SZ surveys.
We present the first public release of our Bayesian inference tool, Bayes-X, for the analysis of X-ray observations of galaxy clusters. We illustrate the use of Bayes-X by analysing a set of four simulated clusters at z=0.2-0.9 as they would be observed by a Chandra-like X-ray observatory. In both the simulations and the analysis pipeline we assume that the dark matter density follows a spherically-symmetric Navarro, Frenk and White (NFW) profile and that the gas pressure is described by a generalised NFW (GNFW) profile. We then perform four sets of analyses. By numerically exploring the joint probability distribution of the cluster parameters given simulated Chandra-like data, we show that the model and analysis technique can robustly return the simulated cluster input quantities, constrain the cluster physical parameters and reveal the degeneracies among the model parameters and cluster physical parameters. We then analyse Chandra data on the nearby cluster, A262, and derive the cluster physical profiles. To illustrate the performance of the Bayesian model selection, we also carried out analyses assuming an Einasto profile for the matter density and calculated the Bayes factor. The results of the model selection analyses for the simulated data favour the NFW model as expected. However, we find that the Einasto profile is preferred in the analysis of A262. The Bayes-X software, which is implemented in Fortran 90, is available at http://www.mrao.cam.ac.uk/facilities/software/bayesx/.
We present cosmological constraints from measurements of the gas mass fraction, $f_{gas}$, for massive, dynamically relaxed galaxy clusters. Our data set consists of Chandra observations of 40 such clusters, identified in a comprehensive search of the Chandra archive, as well as high-quality weak gravitational lensing data for a subset of these clusters. Incorporating a robust gravitational lensing calibration of the X-ray mass estimates, and restricting our measurements to the most self-similar and accurately measured regions of clusters, significantly reduces systematic uncertainties compared to previous work. Our data for the first time constrain the intrinsic scatter in $f_{gas}$, $(7.4pm2.3)$% in a spherical shell at radii 0.8-1.2 $r_{2500}$, consistent with the expected variation in gas depletion and non-thermal pressure for relaxed clusters. From the lowest-redshift data in our sample we obtain a constraint on a combination of the Hubble parameter and cosmic baryon fraction, $h^{3/2}Omega_b/Omega_m=0.089pm0.012$, that is insensitive to the nature of dark energy. Combined with standard priors on $h$ and $Omega_b h^2$, this provides a tight constraint on the cosmic matter density, $Omega_m=0.27pm0.04$, which is similarly insensitive to dark energy. Using the entire cluster sample, extending to $z>1$, we obtain consistent results for $Omega_m$ and interesting constraints on dark energy: $Omega_Lambda=0.65^{+0.17}_{-0.22}$ for non-flat $Lambda$CDM models, and $w=-0.98pm0.26$ for flat constant-$w$ models. Our results are both competitive and consistent with those from recent CMB, SNIa and BAO data. We present constraints on models of evolving dark energy from the combination of $f_{gas}$ data with these external data sets, and comment on the possibilities for improved $f_{gas}$ constraints using current and next-generation X-ray observatories and lensing data. (Abridged)
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
