arxiv: 2604.22200 · v3 · submitted 2026-04-24 · 🌌 astro-ph.GA · stat.AP

Recognition: unknown

Formalizing Galaxy Population Evolution: Drift and Mergers as Transport Processes on Manifolds

Tsutomu T. Takeuchi (Nagoya University , Institute of Statistical Mathematics)

Authors on Pith no claims yet

Pith reviewed 2026-05-08 11:01 UTC · model grok-4.3

classification 🌌 astro-ph.GA stat.AP

keywords galaxy evolutionprobability measuresmanifold learningluminosity functionstransport processesmergersdriftSchechter function

0 comments

The pith

Galaxy evolution is the transport of probability measures on a latent manifold, with observations as projections.

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

The paper proposes treating galaxy populations as probability densities evolving on a manifold of galaxy states. This unifies continuous changes like mass growth with sudden mergers as different types of transport on the same space. Observational quantities such as luminosity functions are seen as projections of this underlying dynamics rather than independent equations. Existing models for mass functions and coagulation emerge as special cases when the manifold is reduced. This view turns the problem of understanding galaxy evolution into recovering the hidden dynamics from data projections.

Core claim

Galaxy evolution is formulated as the time evolution of a probability measure on the galaxy manifold M. With galaxy states as latent variables θ in M and population density ρ(θ,t), the dynamics follow a general equation combining continuous transport and nonlocal jump processes. Luminosity and stellar mass functions arise as projections, and the framework recovers continuity equations and the Smoluchowski equation as limits.

What carries the argument

The general transport equation for the probability density ρ(θ,t) on the manifold, where drift represents continuous evolution and jumps represent mergers.

Load-bearing premise

Galaxy states can be represented by latent variables on a manifold such that all observational statistics are projections of the probability measure with no essential dynamical information lost.

What would settle it

Finding galaxy population statistics that cannot be reproduced as projections of any single evolving density on a manifold, such as luminosity functions whose time dependence requires independent equations unrelated to any underlying transport.

read the original abstract

Galaxy evolution is commonly described through the time evolution of observational statistics such as luminosity functions and stellar mass functions. However, these quantities are projections of an underlying multivariate galaxy state space rather than fundamental dynamical variables. We develop a unified framework in which galaxy evolution is formulated as the time evolution of a probability measure on the galaxy manifold. Representing galaxy states by latent variables $\theta\in\mathcal{M}$ and the population by a density $\rho(\theta,t)$, the evolution is governed by a general equation containing continuous transport and nonlocal jump processes. By reinterpreting manifold learning as the pushforward of measures, we distinguish observational, representation, and physical measures, and emphasize that manifold coordinates themselves need not carry direct physical meaning. In this picture, luminosity functions and stellar mass functions arise as projected observables of a single underlying dynamics, and generally do not form closed equations in observational space. The framework contains existing models as limiting cases: reduction to a single mass variable yields continuity-equation models, while additive post-merger states recover the Smoluchowski coagulation equation. We further show that luminosity-function evolution is naturally described within the Schechter family, whose apparent stability is interpreted as an effective consequence of projection. Since observables are projections of measures, inference of galaxy evolution becomes a statistical inverse problem of recovering manifold dynamics from data. This framework shifts the focus from fitting observed statistics directly to inferring the underlying state-space dynamics, thereby bridging manifold learning and physical theory.

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 recasts galaxy evolution as measure transport on a manifold with drift and jumps, recovering some standard models as limits, but stays conceptual without derivations or tests.

read the letter

The main thing here is a proposal to treat galaxy population changes as the evolution of a probability measure on a manifold, where continuous transport handles drift and nonlocal jumps handle mergers, with observed statistics like luminosity functions arising as projections of that measure. It distinguishes observational, representation, and physical measures and ties manifold learning to pushforwards of measures. Existing approaches appear as special cases: drop to one mass variable and you get continuity equations; allow additive post-merger states and you recover the Smoluchowski coagulation equation. The stability of Schechter forms is presented as a projection effect rather than something requiring separate tuning. The setup is internally consistent and avoids claiming to fit data directly, instead framing inference as an inverse problem of recovering manifold dynamics. That is the actual novelty in the combination, and it is a clean way to organize things that are usually handled piecemeal. The soft spot is the absence of the actual transport equation or any derivation showing how the general form reduces to those limits or keeps luminosity functions inside the Schechter family. Without explicit math, error estimates, or even a toy example, it is difficult to judge whether this changes modeling practice or mainly reorganizes existing ideas. The manifold assumption is stated as a choice, which is fair, but the payoff still depends on whether a useful latent space can be found and whether the dynamics are recoverable from projections. This is for theorists who work on formal descriptions of galaxy populations and are comfortable with measure theory. Someone running simulations or fitting catalogs directly will not get immediate tools or predictions. It deserves peer review because the structure holds together and the unification is new enough that referees can ask for the missing derivations and see whether it leads to anything testable.

Referee Report

3 major / 2 minor

Summary. The paper proposes a unified framework for galaxy population evolution by modeling galaxy states as latent variables θ on a manifold M, with the population described by a probability density ρ(θ, t) whose time evolution is governed by a general equation incorporating continuous transport (drift) and nonlocal jump processes (mergers). It claims that standard models are recovered as limits (continuity equations for single-variable reduction; Smoluchowski coagulation for additive post-merger states), that luminosity and stellar-mass functions arise as projections of the underlying measure dynamics and remain within the Schechter family as an effective consequence, and that inference reduces to a statistical inverse problem of recovering manifold dynamics from observational projections.

Significance. If the explicit derivations and limiting-case recoveries hold, the framework would offer a substantive conceptual unification that distinguishes physical, representation, and observational measures, interprets the apparent stability of Schechter functions as a projection artifact, and bridges manifold-learning techniques with physical galaxy-evolution theory. The parameter-free axiomatic structure and explicit framing of the manifold assumption as a modeling choice rather than a derived result are positive features that could facilitate falsifiable tests once the transport equation is specified.

major comments (3)

Abstract: the central claim that 'luminosity-function evolution is naturally described within the Schechter family' and that 'existing models are recovered as limiting cases' is asserted without any explicit form of the general evolution equation for ρ(θ, t), without the projection operator, and without the derivation showing how the transport-plus-jump dynamics projects onto the Schechter functional form or recovers the continuity/Smoluchowski equations; this leaves the load-bearing unification claim without demonstrated support.
Abstract: the reduction 'to a single mass variable yields continuity-equation models' and 'additive post-merger states recover the Smoluchowski coagulation equation' is stated as a fact, yet no intermediate steps, coordinate choices, or limiting procedures are supplied, rendering verification of these recoveries impossible and undermining the assertion that the framework contains prior models as special cases.
Abstract: the weakest modeling assumption—that observational statistics arise strictly as projections of the measure on M 'without essential loss of dynamical information'—is flagged as a choice, but its quantitative consequences for the inverse-problem formulation (recoverability, uniqueness, or information loss under realistic projection kernels) are not examined, leaving the practical utility of the framework untested.

minor comments (2)

Abstract: the notation θ∈M and ρ(θ, t) is introduced without a brief statement of the dimension or topology of M, which would help readers assess the generality of the construction.
Abstract: the phrase 'manifold coordinates themselves need not carry direct physical meaning' is conceptually important but left without an illustrative example of a coordinate choice that is purely representational.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We are grateful to the referee for the thoughtful and constructive report, which has helped us identify areas where the presentation of our framework can be strengthened. We address each of the major comments in detail below and outline the revisions we will make to the manuscript.

read point-by-point responses

Referee: Abstract: the central claim that 'luminosity-function evolution is naturally described within the Schechter family' and that 'existing models are recovered as limiting cases' is asserted without any explicit form of the general evolution equation for ρ(θ, t), without the projection operator, and without the derivation showing how the transport-plus-jump dynamics projects onto the Schechter functional form or recovers the continuity/Smoluchowski equations; this leaves the load-bearing unification claim without demonstrated support.

Authors: We agree that the abstract, constrained by length, presents the claims at a high level without the supporting equations or derivations. The full manuscript provides these in detail: the general evolution equation for ρ(θ, t) is introduced in Section 2 as ∂_t ρ = -∇·(v ρ) + ∫ [K(θ,θ') ρ(θ') - K(θ',θ) ρ(θ)] dθ', where v is the drift velocity and K the jump kernel. The projection operator is defined in Section 3 as the marginalization or pushforward onto observable coordinates. The projection onto luminosity functions and their evolution within the Schechter family is derived in Section 4 by assuming specific forms for the transport that preserve the functional form under projection. Similarly, the limiting cases are shown in Section 5. To make this more accessible, we will revise the abstract to include a concise statement of the general equation and explicitly reference the sections containing the derivations. This revision will be made in the updated manuscript. revision: yes
Referee: Abstract: the reduction 'to a single mass variable yields continuity-equation models' and 'additive post-merger states recover the Smoluchowski coagulation equation' is stated as a fact, yet no intermediate steps, coordinate choices, or limiting procedures are supplied, rendering verification of these recoveries impossible and undermining the assertion that the framework contains prior models as special cases.

Authors: The referee correctly notes that the abstract does not supply the intermediate steps. These are provided in the body of the paper. In Section 5.1, we choose coordinates where the manifold is reduced to mass m by marginalizing over other latent variables, leading to the continuity equation ∂_t n(m,t) + ∂_m (v_m n) = 0 after appropriate limits on the jump term. For the Smoluchowski case in Section 5.2, we assume the post-merger state is additive in mass, θ'' = θ + θ', and the kernel K becomes the coagulation kernel, recovering the standard integral equation. To facilitate verification, we will add an appendix with the explicit coordinate transformations and limiting procedures. We believe this addresses the concern while maintaining the manuscript's focus. revision: yes
Referee: Abstract: the weakest modeling assumption—that observational statistics arise strictly as projections of the measure on M 'without essential loss of dynamical information'—is flagged as a choice, but its quantitative consequences for the inverse-problem formulation (recoverability, uniqueness, or information loss under realistic projection kernels) are not examined, leaving the practical utility of the framework untested.

Authors: We acknowledge that while the manuscript identifies the inverse problem in Section 6, it does not provide a quantitative analysis of information loss or uniqueness under specific projection kernels. This is a valid point, and a full treatment would require specifying observational kernels and deriving bounds, which is beyond the scope of the current work but important for practical applications. We will add a paragraph in the discussion section outlining the conditions under which the dynamics are recoverable (e.g., when the projection is invertible or when additional priors are used) and noting potential information loss as a direction for future research. This will be a partial revision, as a complete quantitative study is left for subsequent papers. revision: partial

Circularity Check

0 steps flagged

No significant circularity: general mathematical framework with independent content

full rationale

The paper formulates galaxy evolution as the time evolution of a probability measure on a manifold, governed by a general transport equation with continuous drift and nonlocal jumps. Existing models are recovered as explicit limiting cases (single-variable continuity equation; Smoluchowski coagulation), and luminosity functions are shown to arise as projections rather than closed dynamical variables. No step reduces a claimed prediction or uniqueness result to a fitted parameter, self-citation, or definitional tautology; the central construction is presented as a unifying ansatz whose consequences are derived from the stated measure-transport structure. The inverse-problem interpretation follows directly from the projection distinction without circular closure.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The framework rests on standard concepts from measure theory and manifold learning with no new free parameters or invented physical entities listed in the abstract.

axioms (2)

domain assumption Galaxy states can be represented by latent variables θ belonging to a manifold M
Explicitly stated when introducing ρ(θ,t) on M
domain assumption Observational statistics such as luminosity functions are projections of the underlying probability measure
Stated when explaining that luminosity functions arise as projected observables

pith-pipeline@v0.9.0 · 5567 in / 1366 out tokens · 63411 ms · 2026-05-08T11:01:23.212353+00:00 · methodology

discussion (0)

Forward citations

Cited by 1 Pith paper

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

A Measure-Theoretic Transport Formulation of Galaxy Evolution on the Galaxy Manifold: Geometric Constraints
physics.gen-ph 2026-04 unverdicted novelty 6.0

Galaxy evolution is cast as a geometrically constrained reaction-transport process on probability measures, using Wasserstein distance and CD(K,∞) conditions to enforce energy dissipation and interaction closure.

Reference graph

Works this paper leans on

20 extracted references · 7 canonical work pages · cited by 1 Pith paper · 1 internal anchor

[1]

Gott, III and M.J

J.R. Gott, III and M.J. Rees,A theory of galaxy formation and clustering., A&A45(1975) 365

1975
[2]

Press and P

W.H. Press and P. Schechter,Formation of Galaxies and Clusters of Galaxies by Self-Similar Gravitational Condensation, ApJ187(1974) 425

1974
[3]

J.R. Bond, S. Cole, G. Efstathiou and N. Kaiser,Excursion Set Mass Functions for Hierarchical Gaussian Fluctuations, ApJ379(1991) 440

1991
[4]

White and C.S

S.D.M. White and C.S. Frenk,Galaxy Formation through Hierarchical Clustering, ApJ379 (1991) 52

1991
[5]

Silk,Protogalaxy interactions in newly formed clusters: galaxy luminosities, colors, and intergalactic gas., ApJ220(1978) 390

J. Silk,Protogalaxy interactions in newly formed clusters: galaxy luminosities, colors, and intergalactic gas., ApJ220(1978) 390

1978
[6]

Silk and S.D

J. Silk and S.D. White,The development of structure in the expanding universe., ApJ223 (1978) L59

1978
[7]

Smoluchowski,Versuch einer mathematischen theorie der koagulationskinetik kolloider lösungen,Zeitschrift für Physikalische Chemie92(1917) 129

M. Smoluchowski,Versuch einer mathematischen theorie der koagulationskinetik kolloider lösungen,Zeitschrift für Physikalische Chemie92(1917) 129

1917
[8]

Schechter,An analytic expression for the luminosity function for galaxies., ApJ203(1976) 297

P. Schechter,An analytic expression for the luminosity function for galaxies., ApJ203(1976) 297

1976
[9]

Neistein, F.C

E. Neistein, F.C. van den Bosch and A. Dekel,Natural downsizing in hierarchical galaxy formation, MNRAS372(2006) 933 [astro-ph/0605045]. – 29 –

work page arXiv 2006
[10]

Neistein and A

E. Neistein and A. Dekel,Merger rates of dark matter haloes, MNRAS388(2008) 1792 [0802.0198]

work page arXiv 2008
[11]

Drory and M

N. Drory and M. Alvarez,The Contribution of Star Formation and Merging to Stellar Mass Buildup in Galaxies, ApJ680(2008) 41 [0803.1489]

work page arXiv 2008
[12]

Peng, S.J

Y.-j. Peng, S.J. Lilly, K. Kovač, M. Bolzonella, L. Pozzetti, A. Renzini et al.,Mass and Environment as Drivers of Galaxy Evolution in SDSS and zCOSMOS and the Origin of the Schechter Function, ApJ721(2010) 193

2010
[13]

Brosche,The Manifold of Galaxies

P. Brosche,The Manifold of Galaxies. Galaxies with known Dynamical Parameters,Astronomy and Astrophysics23(1973) 259

1973
[14]

Djorgovski,Galaxy manifolds and galaxy formation, inMorphological and Physical Classification of Galaxies, G

S. Djorgovski,Galaxy manifolds and galaxy formation, inMorphological and Physical Classification of Galaxies, G. Longo, M. Capaccioli and G. Busarello, eds., (Dordrecht), pp. 337–356, Springer Netherlands, 1992

1992
[15]

Takeuchi,Applications of Big Data and Machine Learning in Galaxy Formation and Evolution, Series in Astronomy and Astrophysics, CRC Press, Boca Raton (2025)

T.T. Takeuchi,Applications of Big Data and Machine Learning in Galaxy Formation and Evolution, Series in Astronomy and Astrophysics, CRC Press, Boca Raton (2025)

2025
[16]

Takeuchi, S

T.T. Takeuchi, S. Cooray and R. Kano,Galaxy Evolution with Manifold Learning,Entropy28 (2026) 288(28pp)

2026
[17]

A Measure-Theoretic Transport Formulation of Galaxy Evolution on the Galaxy Manifold: Geometric Constraints

T.T. Takeuchi,A measure-theoretic transport formulation of galaxy evolution on the galaxy manifold: Geometric constraints,arXiv(2026) [2604.24823]

work page internal anchor Pith review Pith/arXiv arXiv 2026
[18]

Cowie, A

L.L. Cowie, A. Songaila, E.M. Hu and J.G. Cohen,New Insight on Galaxy Formation and Evolution From Keck Spectroscopy of the Hawaii Deep Fields, AJ112(1996) 839 [astro-ph/9606079]

work page arXiv 1996
[19]

Takeuchi, K

T.T. Takeuchi, K. Yoshikawa and T.T. Ishii,Tests of Statistical Methods for Estimating Galaxy Luminosity Function and Applications to the Hubble Deep Field, ApJS129(2000) 1 [astro-ph/0003127]

work page arXiv 2000
[20]

Johnston,Shedding light on the galaxy luminosity function, A&A Rev.19(2011) 41 [1106.2039]

R. Johnston,Shedding light on the galaxy luminosity function, A&A Rev.19(2011) 41 [1106.2039]. – 30 –

work page arXiv 2011