Source author record

F. R. Pearce

F. R. Pearce appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

13works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

13 published item(s)

preprint2015arXiv

Galaxy Cluster Mass Reconstruction Project: II. Quantifying scatter and bias using contrasting mock catalogues

This article is the second in a series in which we perform an extensive comparison of various galaxy-based cluster mass estimation techniques that utilise the positions, velocities and colours of galaxies. Our aim is to quantify the scatter, systematic bias and completeness of cluster masses derived from a diverse set of 25 galaxy-based methods using two contrasting mock galaxy catalogues based on a sophisticated halo occupation model and a semi-analytic model. Analysing 968 clusters, we find a wide range in the RMS errors in log M200c delivered by the different methods (0.18 to 1.08 dex, i.e., a factor of ~1.5 to 12), with abundance matching and richness methods providing the best results, irrespective of the input model assumptions. In addition, certain methods produce a significant number of catastrophic cases where the mass is under- or over-estimated by a factor greater than 10. Given the steeply falling high-mass end of the cluster mass function, we recommend that richness or abundance matching-based methods are used in conjunction with these methods as a sanity check for studies selecting high mass clusters. We see a stronger correlation of the recovered to input number of galaxies for both catalogues in comparison with the group/cluster mass, however, this does not guarantee that the correct member galaxies are being selected. We do not observe significantly higher scatter for either mock galaxy catalogues. Our results have implications for cosmological analyses that utilise the masses, richnesses, or abundances of clusters, which have different uncertainties when different methods are used.

preprint2014arXiv

Galaxy Cluster Mass Reconstruction Project: I. Methods and first results on galaxy-based techniques

This paper is the first in a series in which we perform an extensive comparison of various galaxy-based cluster mass estimation techniques that utilise the positions, velocities and colours of galaxies. Our primary aim is to test the performance of these cluster mass estimation techniques on a diverse set of models that will increase in complexity. We begin by providing participating methods with data from a simple model that delivers idealised clusters, enabling us to quantify the underlying scatter intrinsic to these mass estimation techniques. The mock catalogue is based on a Halo Occupation Distribution (HOD) model that assumes spherical Navarro, Frenk and White (NFW) haloes truncated at R_200, with no substructure nor colour segregation, and with isotropic, isothermal Maxwellian velocities. We find that, above 10^14 M_solar, recovered cluster masses are correlated with the true underlying cluster mass with an intrinsic scatter of typically a factor of two. Below 10^14 M_solar, the scatter rises as the number of member galaxies drops and rapidly approaches an order of magnitude. We find that richness-based methods deliver the lowest scatter, but it is not clear whether such accuracy may simply be the result of using an over-simplistic model to populate the galaxies in their haloes. Even when given the true cluster membership, large scatter is observed for the majority non-richness-based approaches, suggesting that mass reconstruction with a low number of dynamical tracers is inherently problematic.

preprint2013arXiv

Brighter galaxy bias: underestimating the velocity dispersions of galaxy clusters

We study the systematic bias introduced when selecting the spectroscopic redshifts of brighter cluster galaxies to estimate the velocity dispersion of galaxy clusters from both simulated and observational galaxy catalogues. We select clusters with Ngal > 50 at five low redshift snapshots from a semi-analytic model galaxy catalogue, and from a catalogue of SDSS DR8 groups and clusters across the redshift range 0.021<z<0.098. We employ various selection techniques to explore whether the velocity dispersion bias is simply due to a lack of dynamical information or is the result of an underlying physical process occurring in the cluster, for example, dynamical friction. The velocity dispersions and stacked particle velocity distributions of the parent dark matter (DM) halos are compared to the corresponding cluster dispersions and galaxy velocity distribution. We find a clear bias between the halo and the semi-analytic galaxy cluster velocity dispersion on the order of sigma gal / sigma DM = 0.87-0.95 and a distinct difference in the stacked galaxy and DM particle velocity distribution. We identify a systematic underestimation of the velocity dispersions when imposing increasing absolute I-band magnitude limits. This underestimation is enhanced when using only the brighter cluster members for dynamical analysis on the order of 5-35%, indicating that dynamical friction is a serious source of bias when using galaxy velocities as tracers of the underlying gravitational potential. In contrast to the literature we find that the resulting bias is not only halo mass-dependent but that the nature of the dependence changes according to the galaxy selection strategy. We make a recommendation that, in the realistic case of limited availability of spectral observations, a strictly magnitude-limited sample should be avoided to ensure an unbiased estimate of the velocity dispersion.

preprint2010arXiv

The evolution of galaxy cluster X-ray scaling relations

We use numerical simulations to investigate, for the first time, the joint effect of feedback from supernovae (SNe) and active galactic nuclei (AGN) on the evolution of galaxy cluster X-ray scaling relations. Our simulations are drawn from the Millennium Gas Project and are some of the largest hydrodynamical N-body simulations ever carried out. Feedback is implemented using a hybrid scheme, where the energy input into intracluster gas by SNe and AGN is taken from a semi-analytic model of galaxy formation. This ensures that the source of feedback is a population of galaxies that closely resembles that found in the real universe. We show that our feedback model is capable of reproducing observed local X-ray scaling laws, at least for non-cool core clusters, but that almost identical results can be obtained with a simplistic preheating model. However, we demonstrate that the two models predict opposing evolutionary behaviour. We have examined whether the evolution predicted by our feedback model is compatible with observations of high-redshift clusters. Broadly speaking, we find that the data seems to favour the feedback model for z<0.5, and the preheating model at higher redshift. However, a statistically meaningful comparison with observations is impossible, because the large samples of high-redshift clusters currently available are prone to strong selection biases. As the observational picture becomes clearer in the near future, it should be possible to place tight constraints on the evolution of the scaling laws, providing us with an invaluable probe of the physical processes operating in galaxy clusters.

preprint2006arXiv

A Heating Model for the Millennium Gas Run

The comparison between observations of galaxy clusters thermo-dynamical properties and theoretical predictions suggests that non-gravitational heating needs to be added into the models. We implement an internally self-consistent heating scheme into GADGET-2 for the third (and fourth) run of the Millennium gas project (Pearce et al. in preparation), a set of four hydrodynamical cosmological simulations with N=2(5x10^8) particles and with the same volume (L=500 h-1 Mpc) and structures as the the N-body Millennium Simulation (Springel et al. 2005). Our aim is to reproduce the observed thermo-dynamical properties of galaxy clusters.

preprint2003arXiv

Stable clustering, the halo model and nonlinear cosmological power spectra

We present the results of a large library of cosmological N-body simulations, using power-law initial spectra. The nonlinear evolution of the matter power spectra is compared with the predictions of existing analytic scaling formulae based on the work of Hamilton et al. The scaling approach has assumed that highly nonlinear structures obey `stable clustering' and are frozen in proper coordinates. Our results show that, when transformed under the self-similarity scaling, the scale-free spectra define a nonlinear locus that is clearly shallower than would be required under stable clustering. Furthermore, the small-scale nonlinear power increases as both the power-spectrum index n and the density parameter Omega decrease, and this evolution is not well accounted for by the previous scaling formulae. This breakdown of stable clustering can be understood as resulting from the modification of dark-matter haloes by continuing mergers. These effects are naturally included in the analytic `halo model' for nonlinear structure; using this approach we are able to fit both our scale-free results and also our previous CDM data. This approach is more accurate than the commonly-used Peacock--Dodds formula and should be applicable to more general power spectra. Code to evaluate nonlinear power spectra using this method is available from http://as1.chem.nottingham.ac.uk/~res/software.html Following publication, we will make the power-law simulation data available through the Virgo website http://www.mpa-garching.mpg.de/Virgo

preprint2000arXiv

Collisional versus Collisionless Dark Matter

We compare the structure and substructure of dark matter halos in model universes dominated by collisional, strongly self interacting dark matter (SIDM) and collisionless, weakly interacting dark matter (CDM). While SIDM virialised halos are more nearly spherical than CDM halos, they can be rotationally flattened by as much as 20% in their inner regions. Substructure halos suffer ram-pressure truncation and drag which are more rapid and severe than their gravitational counterparts tidal stripping and dynamical friction. Lensing constraints on the size of galactic halos in clusters are a factor of two smaller than predicted by gravitational stripping, and the recent detection of tidal streams of stars escaping from the satellite galaxy Carina suggests that its tidal radius is close to its optical radius of a few hundred parsecs --- an order of magnitude smaller than predicted by CDM models but consistent with SIDM. The orbits of SIDM satellites suffer significant velocity bias $σ_{SIDM}/σ_{CDM}=0.85$ and are more circular than CDM, $β_{SIDM}} \approx 0.5$, in agreement with the inferred orbits of the Galaxy's satellites. In the limit of a short mean free path, SIDM halos have singular isothermal density profiles, thus in its simplest incarnation SIDM is inconsistent with galactic rotation curves.

preprint1998arXiv

Evolution of structure in cold dark matter universes

We present an analysis of the clustering evolution of dark matter in four cold dark matter (CDM) cosmologies. We use a suite of high resolution, 17-million particle, N-body simulations which sample volumes large enough to give clustering statistics with unprecedented accuracy. We investigate both a flat and an open model with Omega_0=0.3, and two models with Omega=1, one with the standard CDM power spectrum and the other with the same power spectrum as the Omega_0=0.3 models. The amplitude of primordial fluctuations is set so that the models reproduce the observed abundance of rich galaxy clusters by the present day. The mass 2-point correlation function and power spectrum of all the simulations differ significantly from those of the observed galaxy distribution, in both shape and amplitude. Thus, for any of these models to provide an acceptable representation of reality, the distribution of galaxies must be biased relative to the mass in a non-trivial, scale-dependent, fashion. In the Omega=1 models the required bias is always greater than unity, but in the Omega_0=0.3 models an "antibias" is required on scales smaller than \sim 5\hmpc. The mass correlation functions in the simulations are well fit by recently published analytic models. The velocity fields are remarkably similar in all the models, whether they be characterised as bulk flows, 1-particle or pairwise velocity dispersions. This similarity is a direct consequence of our adopted normalisation. The small-scale pairwise velocity dispersion of the dark matter is somewhat larger than recent determinations from galaxy redshift surveys, but the bulk-flows predicted by our models are broadly in agreement with most available data.

preprint1997arXiv

The structure of galaxy clusters in different cosmologies

We investigate the internal structure of clusters of galaxies in high-resolution N-body simulations of 4 different cosmologies. There is a higher proportion of disordered clusters in critical-density than in low-density universes, although the structure of relaxed clusters is very similar in each. Crude measures of substructure, such as the shift in the position of the centre-of-mass as the density threshold is varied, can distinguish the two in a sample of just 20 or so clusters; it is harder to differentiate between clusters in open and flat models with the same density parameter. Most clusters are in a quasi-steady state within the virial radius and are well-described by the density profile of Navarro, Frenk & White (1995).

preprint1996arXiv

The Virgo consortium: simulations of dark matter and galaxy clustering

We report on work in progress by the Virgo consortium, a collaboration set up to carry out large simulations of the formation of galaxies and large-scale structure exploiting the latest generation of parallel supercomputers. We show results of $256^3$ particle N-body simulations of the clustering evolution of dark matter in four cold dark matter models with different cosmological parameters. The high resolution and large volume of these simulations allows us to determine reliably the mass autocorrelation function for pair separations in the range $40\hkpc$ to $20\hmpc$. Comparison of these with the observed galaxy correlation function shows that for any of these models to be viable, the distribution of galaxies must be biased relative to the distribution of mass in a non-trivial, scale-dependent fashion. In particular, low $Ω_0$ models require the galaxies to be more ``weakly'' clustered than the mass at small and intermediate pair separations. Simulations which include the evolution of gas show that cold gas knots form with approximately the abundance expected on theoretical grounds, although a few excessively massive objects grow near the centres of rich clusters. The locations where these cold gas knots form are, in general, biased relative to the distribution of mass in a scale-dependent way. Some of these biases have the required sign but they are, for the most part, weaker than is necessary for agreement with observations. The antibias present in our low $Ω_0$ N-body/SPH simulation appears to be related to the merging and disruption of galaxies in rich clusters.

preprint1993arXiv

Mergers of Systems Containing Gas

Several simple mergers between model galaxy clusters containing a mixture of gas and dark matter are examined, testing the coupling of the gas to the underlying collisionless material. The gas is shocked, irreversibly dissipating the energy fed into it by the collisionless component and forms a resolved constant-density core. For the dark matter, however, admixture of phase space vacuum is not very efficient and a constant-density core is not produced. In the final state the central gas has little residual kinetic energy, indicating that streaming motions do not help to support the gas.

preprint1993arXiv

On the probability of major-axis precession in triaxial ellipsoidal potentials

Orbits in triaxial ellipsoidal potentials precess about either the major or minor axis of the ellipsoid. In standard perturbation theory it can be shown that a circular orbit will precess about the minor axis if its angular momentum vector lies in a region bounded by two great circles which pass through the intermediate axis and which are inclined with minimum separation $i_T$ from the minor axis. We test the accuracy of the standard formula for $i_T$ by performing orbit integrations to determine $i_S$, the simulated turnover angle corresponding to $i_T$. We reach two principal conclusions: (i) $i_S$ is usually greater than $i_T$, by as much as 12 degrees even for moderate triaxialities, $A/1.2<B<C/0.8$. This reduces the expected frequency of polar rings. (ii) $i_S$ is not a single, well-defined number but can vary by a few degrees depending upon the initial phase of the orbit. This means that there is a reasonable probability for capture of gas onto orbits which precess about both axes. Interactions can then lead to substantial loss of angular momentum and subsequent infall to the galactic centre.