BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

John A. Peacock

John A. Peacock contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.COastro-ph.GAastro-ph.IMphysics.comp-phastro-ph.SRgr-qc

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
12works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

12 published item(s)

preprint2022arXiv

Accelerating astronomical and cosmological inference with Preconditioned Monte Carlo

We introduce Preconditioned Monte Carlo (PMC), a novel Monte Carlo method for Bayesian inference that facilitates efficient sampling of probability distributions with non-trivial geometry. PMC utilises a Normalising Flow (NF) in order to decorrelate the parameters of the distribution and then proceeds by sampling from the preconditioned target distribution using an adaptive Sequential Monte Carlo (SMC) scheme. The results produced by PMC include samples from the posterior distribution and an estimate of the model evidence that can be used for parameter inference and model comparison respectively. The aforementioned framework has been thoroughly tested in a variety of challenging target distributions achieving state-of-the-art sampling performance. In the cases of primordial feature analysis and gravitational wave inference, PMC is approximately 50 and 25 times faster respectively than Nested Sampling (NS). We found that in higher dimensional applications the acceleration is even greater. Finally, PMC is directly parallelisable, manifesting linear scaling up to thousands of CPUs.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

An empirical measurement of the Halo Mass Function from the combination of GAMA DR4, SDSS DR12, and REFLEX II data

We construct the halo mass function (HMF) from the GAMA galaxy group catalogue over the mass range 10^12.7M_sol to 10^15.5M_sol, and find good agreement with the expectation from LambdaCDM. In comparison to previous studies, this result extends the mass range over which the HMF has now been measured over by an order of magnitude. We combine the GAMA DR4 HMF with similar data from the SDSS DR12 and REFLEX II surveys, and fit a four-parameter Murray-Robotham-Power (MRP) function, valid at z~0.1, yielding: a density normalisation of: log10 (phi Mpc^3)=-3.96[+0.55,-0.82], a high mass turn-over of: log10(M/M_sol)=14.13[+0.43,-0.40], a low mass power law slope of: alpha=-1.68[+0.21,-0.24] , and a high mass softening parameter of: beta= 0.63[+0.25,-0.11]. If we fold in the constraint on Omega_M from Planck 2018 Cosmology, we are able to reduce these uncertainties further, but this relies on the assumption that the power-law trend can be extrapolated from 10^12.7M_sol to zero mass. Throughout, we highlight the effort needed to improve on our HMF measurement: improved halo mass estimates that do not rely on calibration to simulations; reduced halo mass uncertainties needed to mitigate the strong Eddington Bias that arises from the steepness of the HMF low mass slope; and deeper wider area spectroscopic surveys. To our halo mass limit of 10^12.7 M_sol, we are directly resolving (`seeing') 41+/-5 per cent of the total mass density, i.e. Omega_[M>12.7]=0.128+/-0.016, opening the door for the direct construction of 3D dark matter mass maps at Mpc resolution.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Extended Hernquist-Springel formalism for cosmic star formation

We present a revised and extended version of the analytic model for cosmic star formation originally given by Hernquist & Springel in 2003. The key assumption of this formalism is that star formation proceeds from cold gas, at a rate that is limited by an internal consumption timescale at early times, or by the rate of generation of gas via cooling at late times. These processes are analysed as a function of the mass of dark matter haloes and integrated over the halo population. We modify this approach in two main ways to make it more general: (1) halo collapse times are included explicitly, so that the behaviour is physically reasonable at late times; (2) allowance is made for a mass-dependent baryon fraction in haloes, which incorporates feedback effects. This model reproduces the main features of the observed baryonic Tully-Fisher relationship, and is consistent with observational estimates of the baryon mass fraction in the intergalactic medium. With minimal adjustment of parameters, our approach reproduces the observed history of cosmic star formation within a factor of two over the redshift range $0 < z < 10$. This level of agreement is comparable to that achieved by state-of-the-art cosmological simulations. Our simplified apparatus has pedagogical value in illuminating the results of such detailed calculations, and also serves as a means for rapid approximate exploration of non-standard cosmological models.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Galaxy And Mass Assembly (GAMA): Data Release 4 and the z < 0.1 total and z < 0.08 morphological galaxy stellar mass functions

In Galaxy And Mass Assembly Data Release 4 (GAMA DR4), we make available our full spectroscopic redshift sample. This includes 248682 galaxy spectra, and, in combination with earlier surveys, results in 330542 redshifts across five sky regions covering ~250deg^2. The redshift density, is the highest available over such a sustained area, has exceptionally high completeness (95 per cent to r_KIDS=19.65mag), and is well suited for the study of galaxy mergers, galaxy groups, and the low redshift (z<0.25) galaxy population. DR4 includes 32 value-added tables or Data Management Units (DMUs) that provide a number of measured and derived data products including GALEX, ESO KiDS, ESO VIKING, WISE and Herschel Space Observatory imaging. Within this release, we provide visual morphologies for 15330 galaxies to z<0.08, photometric redshift estimates for all 18million objects to r_KIDS~25mag, and stellar velocity dispersions for 111830 galaxies. We conclude by deriving the total galaxy stellar mass function (GSMF) and its sub-division by morphological class (elliptical, compact-bulge and disc, diffuse-bulge and disc, and disc only). This extends our previous measurement of the total GSMF down to 10^6.75 M_sol h^-2_70 and we find a total stellar mass density of rho_*=(2.97+/-0.04)x10^8 M_sol h_70 Mpc^-3 or Omega_*=(2.17+/-0.03)x10^-3 h^-1_70. We conclude that at z<0.1, the Universe has converted 4.9+/-0.1 per cent of the baryonic mass implied by Big Bang Nucleosynthesis into stars that are gravitationally bound within the galaxy population.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

pocoMC: A Python package for accelerated Bayesian inference in astronomy and cosmology

pocoMC is a Python package for accelerated Bayesian inference in astronomy and cosmology. The code is designed to sample efficiently from posterior distributions with non-trivial geometry, including strong multimodality and non-linearity. To this end, pocoMC relies on the Preconditioned Monte Carlo algorithm which utilises a Normalising Flow in order to decorrelate the parameters of the posterior. It facilitates both tasks of parameter estimation and model comparison, focusing especially on computationally expensive applications. It allows fitting arbitrary models defined as a log-likelihood function and a log-prior probability density function in Python. Compared to popular alternatives (e.g. nested sampling) pocoMC can speed up the sampling procedure by orders of magnitude, cutting down the computational cost substantially. Finally, parallelisation to computing clusters manifests linear scaling.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Understanding the relation between thermal Sunyaev-Zeldovich decrement and halo mass using the SIMBA and TNG simulations

The relation between the integrated thermal Sunyaev-Zeldovich (tSZ) $y$-decrement versus halo mass ($Y$--$M$) can potentially constrain galaxy formation models, if theoretical and observational systematics can be properly assessed. We investigate the $Y$--$M$ relation in the SIMBA and IllustrisTNG-100 cosmological hydrodynamic simulations, quantifying the effects of feedback, line-of-sight projection, and beam convolution. We find that SIMBA&#39;s AGN jet feedback generates strong deviations from self-similar expectations for the $Y$--$M$ relation, especially at $M_{500}<10^{13}M_{\odot}$. In SIMBA, this is driven by suppressed in-halo $y$ contributions owing to lowered halo baryon fractions. IllustrisTNG results more closely resemble SIMBA without jets. Projections of line-of-sight structures weaken these model differences slightly, but they remain significant -- particularly at group and lower halo masses. In contrast, beam smearing at $\textit{Planck}$ resolution makes the models indistinguishable, and both models appear to agree well with $\textit{Planck}$ data down to the lowest masses probed. We show that the arcminute resolution expected from forthcoming facilities would retain the differences between model predictions, and thereby provide strong constraints on AGN feedback.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Calibration of a star formation and feedback model for cosmological simulations with Enzo

We present results from seventy-one zoom simulations of a Milky Way-sized (MW) halo, exploring the parameter space for a widely-used star formation and feedback model in the {\tt Enzo} simulation code. We propose a novel way to match observations, using functional fits to the observed baryon makeup over a wide range of halo masses. The model MW galaxy is calibrated using three parameters: the star formation efficiency $\left(f_*\right)$, the efficiency of thermal energy from stellar feedback $\left(ε\right)$ and the region into which feedback is injected $\left(r\ {\rm and}\ s\right)$. We find that changing the amount of feedback energy affects the baryon content most significantly. We then identify two sets of feedback parameter values that are both able to reproduce the baryonic properties for haloes between $10^{10}\,\mathrm{M_\odot}$ and $10^{12}\,\mathrm{M_\odot}$. We can potentially improve the agreement by incorporating more parameters or physics. If we choose to focus on one property at a time, we can obtain a more realistic halo baryon makeup. We show that the employed feedback prescription is insensitive to dark matter mass resolution between $10^5\,{\rm M_\odot}$ and $10^7\,{\rm M_\odot}$. Contrasting both star formation criteria and the corresponding combination of optimal feedback parameters, we also highlight that feedback is self-consistent: to match the same baryonic properties, with a relatively higher gas to stars conversion efficiency, the feedback strength required is lower, and vice versa. Lastly, we demonstrate that chaotic variance in the code can cause deviations of approximately 10\% and 25\% in the stellar and baryon mass in simulations evolved from identical initial conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Multi-tracer extension of the halo model: probing quenching and conformity in eBOSS

We develop a new Multi-Tracer Halo Occupation Distribution (\texttt{MTHOD}) framework for the galaxy distribution and apply it to the extended Baryon Oscillation Spectroscopic Survey (eBOSS) final data between $z=0.7-1.1$. We obtain a best fit \mthod\, for each tracer and describe the host halo properties of these galaxies. The mean halo masses for LRGs, ELGs and QSOs are found to be $1.9 \times 10^{13} \msolaroh$, $1.1 \times 10^{12} \msolaroh$ and $5 \times 10^{12} \msolaroh$ respectively in the eBOSS data. We use the \texttt{MTHOD} framework to create mock galaxy catalogues and predict auto- and cross-correlation functions for all the tracers. Comparing these results with data, we investigate galactic conformity, the phenomenon whereby the properties of neighbouring galaxies are mutually correlated in a manner that is not captured by the basic halo model. We detect \textsl{1-halo} conformity at more than 3$σ$ statistical significance, while obtaining upper limit on \textsl{2-halo} conformity. We also look at the environmental dependence of the galaxy quenching efficiency and find that halo mass driven quenching successfully explains the behaviour in high density regions, but it fails to describe the quenching efficiency in low density regions. In particular, we show that the quenching efficiency in low density filaments is higher in the observed data, as compared to the prediction of the \mthod\ with halo mass driven quenching. The mock galaxy catalogue constructed in this paper is publicly available on https://www.roe.ac.uk/~salam/MTHOD/ .

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Primordial black hole merger rates: distributions for multiple LIGO observables

We have calculated the detectable merger rate of primordial black holes, as a function of the redshift, as well as the binary&#39;s mass ratio, total mass and chirp mass (observables that have not previously been explored in great detail for PBHs). We consider both the current and design sensitivity of LIGO and five different primordial black hole mass functions, as well as showing a comparison to a predicted astrophysical black hole merger rate. We show that the empirical preference for nearly equal-mass binaries in current LIGO/Virgo data can be consistent with a PBH hypothesis once observational selection effects are taken into account. However, current data do exclude some PBH mass distributions, and future data may be able to rule out the possibility that all observed BH mergers had a primordial origin.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Protoclusters at z=5.7: A view from the MultiDark galaxies

Protoclusters, which will yield galaxy clusters at lower redshift, can provide valuable information on the formation of galaxy clusters. However, identifying progenitors of galaxy clusters in observations is not an easy task, especially at high redshift. Different priors have been used to estimate the overdense regions that are thought to mark the locations of protoclusters. In this paper, we use mimicked Ly$α$-emitting galaxies at $z=5.7$ to identify protoclusters in the MultiDark galaxies, which are populated by applying three different semi-analytic models to the 1 $Gpc h^{-1}$ MultiDark Planck2 simulation. To compare with observational results, we extend the criterion 1 (a Ly$α$ luminosity limited sample), to criterion 2 (a match to the observed mean galaxy number density). To further statistically study the finding efficiency of this method, we enlarge the identified protocluster sample (criterion 3) to about 3500 at $z=5.7$ and study their final mass distribution. The number of overdense regions and their selection probability depends on the semi-analytic models and strongly on the three selection criteria (partly by design). The protoclusters identified with criterion 1 are associated with a typical final cluster mass of $2.82\pm0.92 \times 10^{15} M_\odot$ which is in agreement with the prediction (within $\pm 1 σ$) of an observed massive protocluster at $z=5.7$. Identifying more protoclusters allows us to investigate the efficiency of this method, which is more suitable for identifying the most massive clusters: completeness ($\mathbb{C}$) drops rapidly with decreasing halo mass. We further find that it is hard to have a high purity ($\mathbb{P}$) and completeness simultaneously.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Three Hundred Project: the stellar and gas profiles

Using the catalogues of galaxy clusters from The Three Hundred project, modelled with both hydrodynamic simulations, (Gadget-X and Gadget-MUSIC), and semi-analytic models (SAMs), we study the scatter and self-similarity of the profiles and distributions of the baryonic components of the clusters: the stellar and gas mass, metallicity, the stellar age, gas temperature, and the (specific) star formation rate. Through comparisons with observational results, we find that the shape and the scatter of the gas density profiles matches well the observed trends including the reduced scatter at large radii which is a signature of self-similarity suggested in previous studies. One of our simulated sets, Gadget-X, reproduces well the shape of the observed temperature profile, while Gadget-MUSIC has a higher and flatter profile in the cluster centre and a lower and steeper profile at large radii. The gas metallicity profiles from both simulation sets, despite following the observed trend, have a relatively lower normalisation. The cumulative stellar density profiles from SAMs are in better agreement with the observed result than both hydrodynamic simulations which show relatively higher profiles. The scatter in these physical profiles, especially in the cluster centre region, shows a dependence on the cluster dynamical state and on the cool-core/non-cool-core dichotomy. The stellar age, metallicity and (s)SFR show very large scatter, which are then presented in 2D maps. We also do not find any clear radial dependence of these properties. However, the brightest central galaxies have distinguishable features compared to the properties of the satellite galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Tomographic measurement of the intergalactic gas pressure through galaxy-tSZ cross-correlations

We cross-correlate maps of the thermal Sunyaev-Zeldovich (tSZ) Compton-$y$ parameter published by Planck with the projected distribution of galaxies in a set of low-redshift tomographic bins. We use the nearly full-sky 2MASS Photometric Redshift and WISE $\times$ SuperCOSMOS public catalogues, covering the redshift range $z\lesssim0.4$. Our measurements allow us to place constraints on the redshift dependence of the mass-observable relation for tSZ cluster count analyses in terms of the so-called &#39;hydrostatic mass bias&#39; parameter $1-b_{\rm H}$. These results can also be interpreted as measurements of the bias-weighted average gas pressure $\langle bP_e\rangle$ as a function of redshift, a quantity that can be related to the thermodynamics of gas inside haloes and used to constrain energy injection processes. We measure $1-b_{\rm H}$ with $\sim6\%$ precision in 6 equispaced redshift bins, and find no evidence for a redshift-dependent mass bias parameter, in agreement with previous analyses. Our mean value of $1-b_{\rm H} = 0.75\pm0.03$ is also in good agreement with the one estimated by the joint analysis of Planck cluster counts and CMB anisotropies calibrated with CMB lensing. Our measurements of $\langle bP_e\rangle$, at the level of $\sim10\%$ in each bin, are the most stringent constraints on the redshift dependence of this parameter to date, and agree well both with previous measurements and with theoretical expectations from shock-heating models.

0 citations0 reviewsNo code linkedOpen access