Pith. Sign in

REVIEW

The connection between massive galaxies and their dark matter haloes from $z=0.25$ to $z \sim 4$ from the combined COSMOS, XMM-LSS and ECDF-S fields

Astrophysics of Galaxies · Submitted 13 Oct 2026 · 155 references

Pith's verdict: Largest combined deep field analysis to date of the galaxy-halo connection at high mass and high redshift; central claim is plausible but N(z) validation is half-done.

T0 · AI review

How this review was made: An AI referee (deepseek/deepseek-v4-flash) read the paper on October 15, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 2 major objections · 5 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv ID2610.19030 (version 1)
Pith IDMBC67QGQ, Pith's permanent record of this paper and its review
CitationCatherine L. Hale, Joel Hamlett, Matt J. Jarvis, Natalia Stylianou, Peter W. Hatfield, Rohan G. Varadaraj, Tariq Yasin. "The connection between massive galaxies and their dark matter haloes from $z=0.25$ to $z \sim 4$ from the combined COSMOS, XMM-LSS and ECDF-S fields." arXiv:2610.19030 (2026). https://pith.science/paper/MBC67QGQ
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report
The paper

The paper claims

Massive galaxies stop growing in stars while their dark-matter haloes keep growing below redshift 1.

The reading

This paper combines optical and near-infrared data from three extragalactic fields covering 9.3 square degrees to study how massive galaxies populate dark matter haloes from redshift 0.25 to 4. By measuring the angular clustering of stellar-mass-selected galaxies and fitting a five-parameter halo occupation distribution model, the authors obtain the stellar-to-halo mass relation and its evolution. They find that for the most massive galaxies, the halo mass at which the stellar-to-halo mass ratio peaks increases from z=1 to z=0, while the corresponding peak stellar mass stays nearly constant. This indicates that in the last half of cosmic time, the stellar mass growth of the most massive central galaxies has become inefficient while their host haloes continue to grow hierarchically. The authors interpret this as evidence that feedback from active galactic nuclei – likely low-power radio jets – has shut down star formation in these galaxies.

What carries the argument

The five-parameter halo occupation distribution (HOD) model of Zheng et al. (2005), describing central and satellite galaxy occupation as a function of halo mass. The model is fitted to the angular two-point correlation function measured with the Landy–Szalay estimator, using a Markov chain Monte Carlo method that simultaneously fits the number density and clustering in each of three independent survey fields. The redshift distribution N(z) is constructed from photometric-redshift likelihoods with a magnitude-dependent prior and a Gaussian broadening of 0.04(1+z). The model projects the 3D galaxy power spectrum into the angular correlation function via the Limber approximation, allowing the extraction of the stellar-to-halo mass relation and derived quantities such as galaxy bias, satellite fraction, and peak halo/ stellar masses.

Extended reading notes

Core claim

Over the redshift range z=0.25 to z=4, the stellar-to-halo mass ratio for the most massive galaxy population peaks at z~1. Below z=1, the peak halo mass M_h^peak increases (from ~$10^{12}$.5 to ~$10^{13}$ M_sun in the highest stellar-mass bin), while the peak stellar mass M_*^peak in those haloes remains approximately constant at ~5×$10^{10}$ M_sun. This divergence implies that from z=1 to today, massive haloes grow through hierarchical accretion but their central galaxies do not build up additional stars. The satellite fraction also rises steadily since z~1, consistent with haloes growing via accretion of smaller systems whose central galaxies become satellites. The combined evidence points to a mechanism that suppresses star formation in the most massive central galaxies at late times, and the authors attribute this to radio-mode AGN feedback.

Load-bearing premise

The redshift distribution N(z) of each galaxy sample, constructed from photometric-redshift likelihoods with a magnitude-dependent prior and a Gaussian broadening of width 0.04(1+z), is assumed to be accurate; the authors acknowledge that misrepresentations in the shape of N(z), especially its wings, can cause up to 20–50% differences in the predicted angular correlation function, and that at z>1.5 the HOD parameter sigma_logM hits its upper prior boundary, indicating that the N(z) uncertainty is not fully resolved.

Editorial extensions

If this is right

  • If true, the suppression of stellar-mass growth in massive haloes below z~1 implies that AGN feedback is the primary quenching mechanism for the most massive galaxies in the nearby Universe, consistent with semi-analytic models that couple radio-jet heating to the hot gas halo.
  • The increasing satellite fraction with decreasing redshift suggests that halo mass assembly at z<1 occurs primarily through the accretion of pre-formed galaxies rather than through in-situ star formation in the central galaxy.
  • The redshift z~1 marks the epoch at which the star-formation efficiency in massive haloes peaks; above z=1 the efficiency declines because haloes are still assembling rapidly, and below z=1 it declines because star formation shuts off.
  • The finding motivates targeted radio observations of massive galaxies at z<1 to quantify the duty cycle of low-power radio jets and their correlation with the star-formation history of the central galaxy.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The result provides an observational anchor for galaxy-formation simulations: any model that fails to reproduce the decoupling of stellar-mass and halo-mass growth at z<1 may be missing a key feedback process or overestimating the late-time gas accretion rate.
  • If the suppression is indeed due to AGN feedback, one expects that the most massive central galaxies at z<0.5 should host compact radio cores or show evidence of recent mechanical feedback in X-ray cavities; upcoming surveys with MeerKAT and the SKA can test this directly.
  • The combined-field approach reduces cosmic variance by a factor ~2 compared to single-field studies; future surveys such as Euclid and Rubin will be able to extend this analysis to even larger volumes and fainter stellar-mass thresholds, testing whether the effect holds for lower-mass galaxies as well.
Original abstract

Understanding how galaxies populate and grow within dark matter haloes is essential for linking the observed evolution of galaxies to the hierarchical growth of structure. We present a galaxy clustering analysis using optical and near-infrared data from the UltraVISTA and VIDEO surveys covering $9.3$ deg$^2$ across the COSMOS, XMM-LSS and ECDF-S fields, $\sim 5\times$ larger than previous comparable studies, providing sufficient volume to determine the halo occupation of the most massive galaxies to $z\sim 4$. We measure the angular two-point correlation function for stellar-mass-selected galaxies over $0.25<z<4$ and by combining the three fields, we reduce the cosmic variance by a factor of $\sim2$ compared to the individual fields. The stellar-to-halo mass ratio peaks at a halo mass of $M_{h}\approx10^{12} {\rm M}_{\odot}$ and is a maximum at $z\sim 1$. For the most massive galaxies, the peak halo mass ($M_{h}^{\rm peak}$) increases from $z=1\rightarrow 0$, while the peak stellar mass in such haloes remains approximately constant, indicating that at $z<1$ massive galaxy formation becomes less efficient. We also find a steady increase in the satellite fraction since $z\sim 1$, from which we infer that massive halos continue to grow hierarchically while their central galaxies experience suppressed stellar-mass growth. We suggest this is due to the central massive galaxies having their growth controlled by feedback from active galactic nuclei feedback, likely from low-power stochastic radio emission.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13.0 via TypeSafe API, question bank abstract-p1-20260923, 2026-10-15: 1 request, 5,987 input tokens.

What the abstract itself states, as probabilities: claims first 0.06 · states limitation 0.13 · human participants 0.02 · controlled comparison 0.07 · randomized assignment 0.01 · states quantitative result 0.98 · compares with prior results 0.71 · mentions code or data availability 0.02.

The abstract, sentence by sentence

arXiv listing abstract v1. Results and the main claim are underlined; .

Understanding how galaxies populate and grow within dark matter haloes is essential for linking the observed evolution of galaxies to the hierarchical growth of structure. We present a galaxy clustering analysis using optical and near-infrared data from the UltraVISTA and VIDEO surveys covering $9.3$ deg$^2$ across the COSMOS, XMM-LSS and ECDF-S fields, $\sim 5\times$ larger than previous comparable studies, providing sufficient volume to determine the halo occupation of the most massive galaxies to $z\sim 4$. We measure the angular two-point correlation function for stellar-mass-selected galaxies over $0.25<z<4$ and by combining the three fields, we reduce the cosmic variance by a factor of $\sim2$ compared to the individual fields. The stellar-to-halo mass ratio peaks at a halo mass of $M_{h}\approx10^{12} {\rm M}_{\odot}$ and is a maximum at $z\sim 1$. For the most massive galaxies, the peak halo mass ($M_{h}^{\rm peak}$) increases from $z=1\rightarrow 0$, while the peak stellar mass in such haloes remains approximately constant, indicating that at $z<1$ massive galaxy formation becomes less efficient. We also find a steady increase in the satellite fraction since $z\sim 1$, from which we infer that massive halos continue to grow hierarchically while their central galaxies experience suppressed stellar-mass growth. We suggest this is due to the central massive galaxies having their growth controlled by feedback from active galactic nuclei feedback, likely from low-power stochastic radio emission.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 Understanding how galaxies populate and grow within dark matter haloes is essential for linking the observed evolution of galaxies to the hierarchical growth of structure. background 1.00 0.04 0.03 asserted 1.00
s1 We present a galaxy clustering analysis using optical and near-infrared data from the UltraVISTA and VIDEO surveys covering $9.3$ deg$^2$ across the COSMOS, XMM-LSS and ECDF-S fields, $\sim 5\times$ larger than previous comparable studies, providing sufficient volume to determine the halo occupation of the most massive galaxies to $z\sim 4$. contribution 0.54 0.96 0.04 asserted 1.00
s2 We measure the angular two-point correlation function for stellar-mass-selected galaxies over $0.25<z<4$ and by combining the three fields, we reduce the cosmic variance by a factor of $\sim2$ compared to the individual fields. method 0.89 0.98 0.04 asserted 1.00
s3 The stellar-to-halo mass ratio peaks at a halo mass of $M_{h}\approx10^{12} {\rm M}_{\odot}$ and is a maximum at $z\sim 1$. result 1.00 0.92 0.24 asserted 1.00
s4 For the most massive galaxies, the peak halo mass ($M_{h}^{\rm peak}$) increases from $z=1\rightarrow 0$, while the peak stellar mass in such haloes remains approximately constant, indicating that at $z<1$ massive galaxy formation becomes less efficient. result 0.90 0.92 0.31 asserted 0.58
s5 We also find a steady increase in the satellite fraction since $z\sim 1$, from which we infer that massive halos continue to grow hierarchically while their central galaxies experience suppressed stellar-mass growth. result 0.54 0.97 0.21 asserted 0.98
s6 We suggest this is due to the central massive galaxies having their growth controlled by feedback from active galactic nuclei feedback, likely from low-power stochastic radio emission. interpretation 1.00 0.56 0.18 hedged possibility 0.98

What kind of paper the classifier reports this is: empirical study 1.00 · evidence from analysis of existing data 0.91.

What was examined, and what was not

sources read
  • arXiv listing abstract v1
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 1. Spatial distribution of galaxies in each field after masking. It is important to note that the three fields are not shown to scale with each other, since the two VIDEO fields (XMM-LSS and ECDF-S) are about three times larger than the COSMOS field. S field and any unmatched COSMOS and XMM-LSS sources, we used LePHARE to get photometric redshift measurements, which uses SED fitting techniques and 𝜒 2 minimisation to best match the observed photometry with a library of redshifted galaxy templates. We utilise the COSMOS SED input template provided by LePHARE, which consists of 32…
Figure 2. Colour-colour diagrams on all three fields with 𝐺 − 𝐼 from HSC versus 𝐽 − 𝐾 𝑠 from VISTA. The plots depict sources with 𝐾 𝑠 < 25 in COSMOS and 𝐾 𝑠 < 23.8 in XMM-LSS and ECDF-S. The dashed line indicates the stellar locus fit in the optical/NIR data and the solid line denotes the stellar cut separating the stars at the bottom from the galaxies at the top. are then removed from our sample (see
Figure 3. Stellar mass against the redshift of galaxies in COSMOS, XMM-LSS and ECDF-S after applying masks, a 𝐾 𝑠 magnitude cut, star and point source elimination. The plots are coloured by the number of galaxies in each field, indicated in the colour-bar. The log10 (𝑀lim) lines represent the stellar mass completeness limits for each field as described in Section 2.4.4. The green boxes depict the stellar mass and redshift sample selections used in subsequent analysis. These plots also show how galaxies from each field contribute to different parts of the stellar mass-redshift parameter space.…
Figure 4. Dark matter cosmic variance, 𝜎dm, as a function of mean redshift for the three fields used. The different markers represent the different redshift bin sizes Δ𝑧 that match the 𝑧-bins used in subsequent analyses. 4.2 Cosmic Variance Accurately determining the number density of galaxies is fundamen￾tal for understanding the galaxy-halo relationship, and it directly in￾forms the HOD modelling. However, number density measurements are subject to uncertainties arising from the statistical uncertainty that occurs when observations are limited to a finite volume of the universe. This…
Figure 5. The normalised 𝑁 (𝑧) redshift distributions (solid lines) of the COSMOS, XMM-LSS and ECDF-S fields and COMBINED are shown for each redshift bin used in this study. The 𝑁 (𝑧) distributions broadened by a 0.04(1 + 𝑧) Gaussian that are used during the modelling are also shown in dashed lines for each field. photo-𝑧 estimation code for each galaxy. In this framework, each galaxy is assigned a likelihood function, 𝐿(𝒐|𝑧), which represents the probability of observing its set of fluxes, 𝒐, at a given redshift. By stacking these likelihoods for all galaxies in a selected sample, we get:…
Figure 6. The observed angular two-point correlation functions for COSMOS (in red), XMM-LSS (in blue), and ECDF-S (in purple) fields are shown as circular points for galaxies with stellar mass threshold log10 (𝑀★/M⊙ ) > 9.75 (left block of 4 panels) in four redshift bins: 0.25 < 𝑧 < 0.5 (top left), 0.5 < 𝑧 < 0.75 (top right), 0.75 < 𝑧 < 1 (bottom left) and 1 < 𝑧 < 1.5 (bottom right) and for log10 (𝑀★/M⊙ ) > 10.5 (right block of four panels), for 1.5 < 𝑧 < 2 (top left), 2 < 𝑧 < 2.5 (top right), 2.5 < 𝑧 < 3 (bottom left) and 3 < 𝑧 < 3.5 (bottom right). The combined modelled correlation function…
Figure 7. The observed angular two-point correlation functions for COSMOS (in red), XMM-LSS (in blue), and ECDF-S (in purple) fields are shown as circular points for galaxies at 0.75 < 𝑧 < 1 in four stellar mass thresholds (left block of four panels): log10 (𝑀★/𝑀⊙ ) > 9.25 (top left panel), log10 (𝑀★/M⊙ ) > 10.0 (top right panel), log10 (𝑀★/M⊙ ) > 10.5 (bottom left panel) and log10 (𝑀★/M⊙ ) > 11.0 (bottom right panel), and for galaxies at 1.5 < 𝑧 < 2 in four stellar mass thresholds (right block of four panels): log10 (𝑀★/M⊙ ) > 10.0 (top left panel), log10 (𝑀★/M⊙ ) > 10.25 (top right panel),…
Figure 8. The evolution of the HOD parameters 𝑀1 (typical halo mass for a halo to host one satellite galaxy) and 𝑀min (minimum halo mass for a halo to host a central galaxy) as a function of stellar mass threshold for each redshift bin. 109 1010 1011 Stellar Mass [M ] 101 102 M1/Mmin 0.25 < z < 0.5 0.5 < z < 0.75 0.75 < z < 1 1 < z < 1.5 1.5 < z < 2 2 < z < 2.5 2.5 < z < 3 3 < z < 3.5 3.5 < z < 4
Figure 9. The evolution of the HOD parameters’ ratio 𝑀1/𝑀min as a function of stellar mass threshold for each redshift bin. tively
Figure 10. The evolution of the satellite slope 𝛼 as a function of weighted mean redshift, for different stellar mass thresholds. with halo mass. The total log–normal scatter 𝜎log 𝑀, comprises of an intrinsic component, 𝜎 𝑖 log 𝑀∗ , and a stellar mass measure￾ment error component, 𝜎 𝑚 log 𝑀 , arising from uncertainties in redshift, photometry, and stellar mass modelling, and it can be expressed as(𝜎log 𝑀) 2 = (𝜎 𝑖 log 𝑀 ) 2 + (𝜎 𝑚 log 𝑀 ) 2 (Leauthaud et al. 2011). While the measurement error can differ from study to study due to dif￾ferences in data and stellar mass estimation methods, the…
Figure 11. The evolution of galaxy bias as a function of stellar mass threshold and redshift bin, denoted in the top right legend. The top panel shows our galaxy bias results for each redshift bin. The bottom panel shows a comparison of our results (solid lines) with other papers, presented in the bottom legend, with their equivalent shapes. 6 DERIVED QUANTITIES 6.1 Galaxy Bias The bias (see Equation 12) relates the galaxy overdensity to the dark matter overdensity and can equivalently be defined as the ratio of their correlation functions as: 𝑏𝑔 = 𝛿𝑔 𝛿DM = √︄ 𝜉𝑔 𝜉DM (22) where, 𝛿𝑔 and 𝛿DM…
Figure 12. The fractional cosmic variance (solid lines) and Poisson uncertainties (dashed lines) as a function of redshift are shown for each stellar mass range across the three survey fields. The final plot on the right displays the combined uncertainties for the three fields, where the Poisson errors were combined in quadrature and the cosmic variances were combined using the volume-weighted equation 23 (Moster et al. 2011). Due to its smaller size, we see that COSMOS is more affected by cosmic variance. By including two other larger independent fields, we reduce cosmic variance by a factor…
Figure 13. The evolution of the satellite fraction as a function of stellar mass threshold for each redshift bin, denoted in the top right legend. The top panel shows our satellite fraction results for each redshift bin. The bottom panel shows a comparison of our results (solid lines) with other studies, presented in the bottom legend, with their equivalent shapes. fractions reported by Zaidi et al. 2024 (0.2 < 𝑧 < 0.5, 0.5 < 𝑧 < 1, 1 < 𝑧 < 1.25) rather than the lower values provided by Hatfield et al. 2016 (0.5 < 𝑧 < 0.75, 1 < 𝑧 < 1.25) and McCracken et al. 2015 (1.5 < 𝑧 < 2). Our results at…
Figure 14. The stellar-to-halo mass ratio (SHMR = 𝑀★/𝑀ℎ), for different redshifts (𝑧-bins noted in each legend) and for galaxies above a certain stellar mass (i.e. each panel shows the SHMR of galaxies above a different stellar mass threshold, log10 (𝑀★/𝑀⊙ ) > 9.25, 9.5, 9.75, 10, 10.25, 10.5). The SMHR is separated into contributions from central galaxies (dashed lines), satellite galaxies (dotted lines) and total galaxies (solid lines). On the right y-axis, we present the corresponding integrated star formation efficiency SFE = (𝑀★/𝑀ℎ )/ 𝑓𝑏. The grey shaded area in each plot highlights…
Figure 15. The peak SHMR of the total galaxies, SHMRpeak total (top panel) and the corresponding peaks of halo mass, 𝑀 peak ℎ (central panel) and stellar mass, 𝑀 peak ★ (bottom panel), for different stellar mass thresholds, as a function of the weighted mean redshift. apparent out to 𝑧 ∼ 2.5 for the log10 (𝑀★/M⊙) ≳ 10.25 threshold sample. In the bottom panel of

Figures reproduced from arXiv:2610.19030 by the authors.

References

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.