PAC in DESI. I. Galaxy Stellar Mass Function into the $10^{6}{\rm M}_{\odot}$ Frontier

A. de la Macorra; Arjun Dey; B. A. Weaver; Carlos S. Frenk; C. Poppett; D. Bianchi; D. Brooks; D. Sprayberry; Durham; E. Gazta\~naga

arxiv: 2503.01948 · v4 · submitted 2025-03-03 · 🌌 astro-ph.GA · astro-ph.CO

PAC in DESI. I. Galaxy Stellar Mass Function into the 10⁶{rm M}_(odot) Frontier

Kun Xu (U Penn , Durham , SJTU) , Y. P. Jing , Shaun Cole , Carlos S. Frenk , Sownak Bose , Willem Elbers

show 31 more authors

Wenting Wang Yirong Wang Samuel Moore J. Aguilar S. Ahlen D. Bianchi D. Brooks T. Claybaugh A. de la Macorra Arjun Dey J. E. Forero-Romero E. Gazta\~naga S. Gontcho A Gontcho G. Gutierrez K. Honscheid M. Ishak T. Kisner S. E. Koposov M. Landriau L. Le Guillou R. Miquel J. Moustakas C. Poppett F. Prada I. P\'erez-R\`afols G. Rossi E. Sanchez D. Sprayberry G. Tarl\'e B. A. Weaver H. Zou

This is my paper

Pith reviewed 2026-05-23 01:18 UTC · model grok-4.3

classification 🌌 astro-ph.GA astro-ph.CO

keywords galaxy stellar mass functionPAC methodDESIDECaLSSchechter functionblue galaxiesred galaxieslow-mass galaxies

0 comments

The pith

Improved PAC method derives galaxy stellar mass functions down to 10^{5.3} solar masses for blue galaxies and 10^{6.3} for red ones.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper refines the PAC method to measure the excess surface density of photometric galaxies around spectroscopic tracers without redshift bins. It combines these measurements with wp from the DESI BGS sample to obtain the galaxy stellar mass function over the full range from 10^{5.3} to 10^{11.5} solar masses for blue galaxies and 10^{6.3} to 10^{11.9} for red galaxies. Under the assumption that bias depends mainly on mass and colour, the blue function follows a single Schechter form while the red follows a double Schechter form, with the red population becoming dominant below 10^{7.6} solar masses. This reaches the dwarf regime and opens the way for stage-IV surveys to map the entire local galaxy population.

Core claim

By improving the PAC method and combining n2 wp with wp measurements from the BGS sample, assuming galaxy bias is primarily determined by stellar mass and colour, the galaxy stellar mass functions are derived down to 10^{5.3} M_⊙ for blue galaxies and 10^{6.3} M_⊙ for red galaxies, while also setting lower limits for smaller masses. The blue and red GSMFs are well described by single and double Schechter functions, respectively, with low-mass end slopes of α_blue=-1.54^{+0.02}_{-0.02} and α_red=-2.50^{+0.08}_{-0.08}, resulting in the dominance of red galaxies below 10^{7.6} M_⊙.

What carries the argument

The PAC method, which measures the excess surface density n2 wp of photometric objects with specific properties around spectroscopic tracers.

Load-bearing premise

Galaxy bias is primarily determined by stellar mass and colour.

What would settle it

A complete spectroscopic census of the number density of galaxies with stellar masses between 10^{5.5} and 10^{6.5} solar masses in a local volume would falsify the derived GSMFs if the measured densities deviate from the Schechter extrapolations.

read the original abstract

The Photometric Objects Around Cosmic Webs (PAC) method integrates cosmological photometric and spectroscopic surveys, offering valuable insights into galaxy formation. PAC measures the excess surface density of photometric objects, $\bar{n}_2w_{\rm{p}}$, with specific physical properties around spectroscopic tracers. In this study, we improve the PAC method to make it more rigorous and eliminate the need for redshift bins. We apply the enhanced PAC method to the DESI Y1 BGS Bright spectroscopic sample and the deep DECaLS photometric sample, obtaining $\bar{n}_2w_{\rm{p}}$ measurements across the complete stellar mass range, from $10^{5.3}{\rm M}_{\odot}$ to $10^{11.5}{\rm M}_{\odot}$ for blue galaxies, and from $10^{6.3}{\rm M}_{\odot}$ to $10^{11.9}{\rm M}_{\odot}$ for red galaxies. We combine $\bar{n}_2w_{\rm{p}}$ with $w_{\rm{p}}$ measurements from the BGS sample, which is not necessarily complete in stellar mass. Assuming that galaxy bias is primarily determined by stellar mass and colour, we derive the galaxy stellar mass functions (GSMFs) down to $10^{5.3}{\rm M}_{\odot}$ for blue galaxies and $10^{6.3}{\rm M}_{\odot}$ for red galaxies, while also setting lower limits for smaller masses. The blue and red GSMFs are well described by single and double Schechter functions, respectively, with low-mass end slopes of $\alpha_{\rm{blue}}=-1.54^{+0.02}_{-0.02}$ and $\alpha_{\rm{red}}=-2.50^{+0.08}_{-0.08}$, resulting in the dominance of red galaxies below $10^{7.6}{\rm M}_{\odot}$. Stage-IV cosmological photometric surveys, capable of reaching 2-3 magnitudes deeper than DECaLS, present an opportunity to explore the entire galaxy population in the local universe with PAC. This advancement allows us to address critical questions regarding the nature of dark matter, the physics of reionization, and the formation of dwarf galaxies.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper gives new low-mass GSMF points from an updated PAC estimator on DESI Y1, but the faint-end results rest on the untested claim that bias depends only on mass and color.

read the letter

The main thing to know is that this work extends the local galaxy stellar mass function to 10^5.3 solar masses for blue galaxies and 10^6.3 for red ones by applying a redshift-bin-free version of the PAC method to DESI Y1 BGS data plus DECaLS photometry. They combine the measured excess surface density with wp from the spectroscopic sample to recover the GSMF, fit Schechter functions, and report that red galaxies dominate below 10^7.6 Msun. The slopes come out at -1.54 for blue and -2.50 for red, with lower limits set even fainter. This is new data in a mass range that is usually incomplete in direct surveys. The technical step of removing redshift bins from PAC is a practical improvement that makes the estimator cleaner for this kind of cross-survey work. The measurements themselves are the concrete output that people working on dwarf galaxies or reionization can actually use. The central assumption is that galaxy bias is set primarily by stellar mass and color. If that breaks down at low masses because of environment or assembly bias, the inferred densities shift and the red slope in particular becomes unreliable. The abstract states the assumption but gives no numbers on tests or the size of possible systematics, so the error budget looks incomplete. The quoted uncertainties on the slopes are probably statistical only. The paper cites the relevant DESI, DECaLS, and earlier PAC papers in the usual way. This is aimed at observers and modelers who need faint-end constraints or who want to adapt clustering methods to photometric samples. It has enough new measurements and a documented method that it deserves a serious referee, even though the bias modeling will need close checking. I would send it to peer review.

Referee Report

1 major / 2 minor

Summary. The manuscript presents an improved Photometric Objects Around Cosmic Webs (PAC) method that eliminates redshift bins. Applied to the DESI Y1 BGS Bright spectroscopic sample and DECaLS photometric sample, it measures the excess surface density n̄₂w_p across the full stellar mass range (10^{5.3}–10^{11.5} M_⊙ for blue galaxies; 10^{6.3}–10^{11.9} M_⊙ for red galaxies). Combining these with w_p measurements from the BGS sample and assuming galaxy bias depends primarily on stellar mass and colour, the authors derive the GSMFs down to 10^{5.3} M_⊙ (blue) and 10^{6.3} M_⊙ (red), fit single and double Schechter functions respectively, report low-mass slopes α_blue = −1.54^{+0.02}_{-0.02} and α_red = −2.50^{+0.08}_{-0.08}, and conclude that red galaxies dominate below 10^{7.6} M_⊙. The work also discusses prospects for Stage-IV surveys.

Significance. If the bias assumption is validated, the result extends GSMF measurements into the dwarf regime, providing new empirical constraints on galaxy formation, dark matter models, and reionization. The PAC improvement (removal of redshift bins) is a clear methodological advance that enables deeper probes with future photometric surveys. The combination of independent photometric and spectroscopic datasets is a strength when the modelling assumptions are robust.

major comments (1)

[Abstract] The central claims—the reported GSMF slopes, the numerical values of the number densities at low mass, and the red-galaxy dominance below 10^{7.6} M_⊙—rest on the assumption that galaxy bias is determined primarily by stellar mass and colour (explicitly stated in the abstract). No quantitative robustness tests, sensitivity analyses to environment or assembly bias, or completeness-correction validation at M_* ≲ 10^7 M_⊙ are described; if additional bias dependencies are present, the inferred n(M,colour) and therefore α_red would shift systematically. This is load-bearing for the extrapolation below spectroscopic completeness.

minor comments (2)

The abstract states that lower limits are set for masses below the quoted ranges but does not indicate the method or uncertainty treatment used for these limits.
Notation such as n̄₂w_p and the selection factor should be defined explicitly on first use in the main text.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their constructive review and positive assessment of the methodological advance. We address the single major comment below and will revise the manuscript to strengthen the presentation of our key assumption.

read point-by-point responses

Referee: [Abstract] The central claims—the reported GSMF slopes, the numerical values of the number densities at low mass, and the red-galaxy dominance below 10^{7.6} M_⊙—rest on the assumption that galaxy bias is determined primarily by stellar mass and colour (explicitly stated in the abstract). No quantitative robustness tests, sensitivity analyses to environment or assembly bias, or completeness-correction validation at M_* ≲ 10^7 M_⊙ are described; if additional bias dependencies are present, the inferred n(M,colour) and therefore α_red would shift systematically. This is load-bearing for the extrapolation below spectroscopic completeness.

Authors: We agree that the assumption of bias depending primarily on stellar mass and colour is central and load-bearing for the low-mass extrapolation. This assumption is standard in the literature on colour-dependent clustering and is explicitly stated, but we acknowledge the absence of dedicated quantitative robustness tests in the current manuscript. In revision we will add a new subsection (likely in Section 5 or 6) that (i) reviews supporting evidence from prior clustering studies, (ii) performs sensitivity analyses by varying the bias model (e.g., allowing a modest assembly-bias term or environment dependence) and propagating the effect on α_red and the 10^{7.6} M_⊙ crossover, and (iii) validates the PAC completeness corrections at M_* ≲ 10^7 M_⊙ using mock catalogues that include realistic bias variations. These additions will quantify the systematic uncertainty and make the extrapolation more transparent without altering the core results. revision: yes

Circularity Check

0 steps flagged

No circularity: GSMF derived from independent measurements plus explicit external assumption

full rationale

The paper states its central step explicitly: it measures n̄₂w_p via the improved PAC method on DESI+DECaLS data, combines with w_p from the BGS sample, and obtains n(M,colour) via the relation that divides by b(M,colour) under the stated assumption that bias depends primarily on stellar mass and colour. This assumption is introduced as an input rather than derived inside the analysis, and no equation or section shows the reported α_blue, α_red, or low-mass dominance reducing by construction to quantities already fitted from the same data. No self-citations, fitted-input-as-prediction, or ansatz smuggling appear in the provided text. The result is therefore an extrapolation under a modeling choice, not a tautology.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central claim rests on one domain assumption about galaxy bias and on the validity of the improved PAC estimator; the Schechter parameters are fitted results rather than free parameters introduced by the authors.

free parameters (1)

Schechter function parameters (alpha, phi*, M*)
These are obtained by fitting the derived GSMF data points; they are outputs, not inputs chosen to make the derivation work.

axioms (1)

domain assumption galaxy bias is primarily determined by stellar mass and colour
This assumption is invoked to convert the measured PAC and wp signals into a GSMF.

pith-pipeline@v0.9.0 · 6150 in / 1423 out tokens · 75690 ms · 2026-05-23T01:18:27.890191+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

Assuming that galaxy bias is primarily determined by stellar mass and colour, we derive the galaxy stellar mass functions (GSMFs) down to 10^{5.3} M_⊙ … with low-mass end slopes of α_blue = −1.54 … and α_red = −2.50 …
IndisputableMonolith/Foundation/ArithmeticFromLogic.lean LogicNat ≃ Nat recovery unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

The blue and red GSMFs are well described by single and double Schechter functions, respectively

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

J-PLUS: The stellar mass function of quiescent and star-forming galaxies at 0.05 <= z <= 0.2
astro-ph.GA 2025-09 accept novelty 4.0

J-PLUS DR3 yields Schechter-function stellar mass functions for quiescent and star-forming galaxies at 0.05 <= z <= 0.2 that match prior work, with quiescent galaxies showing higher characteristic mass and flatter fai...