pith. sign in

arxiv: 2605.30243 · v1 · pith:QWX45QWWnew · submitted 2026-05-28 · 🧮 math.DS · math.PR

Energetic characterisation of transient clustering dynamics in aggregation-diffusion systems

Nathalie Wehlitz , Richard Scherzer , Carsten Hartmann , Stefanie Winkelmann This is my paper

Pith reviewed 2026-06-29 00:13 UTC · model grok-4.3

classification 🧮 math.DS math.PR
keywords aggregation-diffusion systemstransient clusteringMcKean-Vlasov equationWasserstein gradient flownon-monotone clusteringenergetic mechanismsdensity peak heightinteracting particle systems
0
0 comments X

The pith

Numerical experiments identify alternating aggregation- and diffusion-dominated regimes during transient clustering in aggregation-diffusion systems.

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

The paper examines transient clustering in nonlocal aggregation-diffusion systems by exploiting the Wasserstein gradient-flow structure of the associated McKean-Vlasov equation on the torus. This structure encodes the thermodynamic competition between interaction-driven aggregation and entropy-driven diffusion. Numerical simulations for locally attractive kernels reveal that trajectories to equilibrium alternate between regimes where one mechanism or the other dominates, producing non-monotone clustering. The work further shows that observables such as density peak height are only partially coupled to the underlying energy changes and therefore do not fully characterise the macroscopic transport.

Core claim

Starting from a stochastic interacting particle system, the macroscopic McKean-Vlasov equation on the torus is analysed via its Wasserstein gradient-flow structure to study the competition between interaction energy and entropy. Numerical experiments with locally attractive kernels identify alternating aggregation- and diffusion-dominated transient regimes along paths to fixed equilibria; these are interpreted as non-monotone clustering. Clustering observables such as density peak height prove only partially coupled to the energetic mechanisms and therefore fail to uniquely characterise the relevant transport dynamics.

What carries the argument

The Wasserstein gradient-flow structure of the McKean-Vlasov equation, which variationally encodes the competition between interaction energy and entropy and thereby organises the transient dynamics.

Load-bearing premise

Numerical experiments on the McKean-Vlasov equation with locally attractive kernels are sufficient to identify alternating regimes and to establish that clustering observables are only partially coupled to energetic mechanisms.

What would settle it

A numerical trajectory in the same setting where the density peak height changes monotonically in lockstep with the rate of energy dissipation would falsify the claim of partial coupling.

Figures

Figures reproduced from arXiv: 2605.30243 by Carsten Hartmann, Nathalie Wehlitz, Richard Scherzer, Stefanie Winkelmann.

Figure 1
Figure 1. Figure 1: Pure diffusion-dominated relaxation to the homogeneous equilibrium. Numerical solution of the PDE (2) for σ = 1.1 starting in ρ0 ∼ Nper(0, 0.5 2 ). (a) Trajectory snapshots of the particle density ρ(x, t), (b) evolution of peak height over time, (c) evolution of free energy over time. The dynamics are diffusion-dominated for all times. allows for loss of stability of the homogeneous equilibrium and the eme… view at source ↗
Figure 2
Figure 2. Figure 2: Pure aggregation-dominated relaxation to the clustered equilibrium. Numerical solution of the PDE (2) for σ = 0.5 starting in ρ0 ∼ Nper(0, 0.5 2 ). (a) Trajectory snapshots of the particle density ρ(x, t), (b) evolution of peak height over time, (c) evolution of free energy over time. The dynamics are aggregation-dominated for all times. mean 0 and variance 0.5 2 , ρ0 ∼ Nper(0, 0.5 2 ). The noise strength … view at source ↗
Figure 3
Figure 3. Figure 3: Alternating transient regimes during relaxation to the ho￾mogeneous equilibrium. Numerical solution of the PDE (2) for σ = 0.838 starting in ρ0 ∼ Nper(0, 0.5 2 ), see Example 1. (a) Trajectory snapshots of the particle density ρ(x, t), (b) evolution of peak height over time, (c) evolution of free energy over time. The dynamics exhibit both aggregation-dominated regimes (0.4 ≲ t ≲ 3.2) as well as diffusion-… view at source ↗
Figure 4
Figure 4. Figure 4: Chemical potential and flux. (a) Chemical potential µ = δF δρ and (b) flux J = −ρ∇µ at time t = 0 for Example 1. Although the entropic contri￾bution induces outward transport in low-density regions, the density-weighted flux is dominated by the highly concentrated central region, leading to an initial accumulation of mass near the center. As a consequence, mass initially accumulates at the center despite t… view at source ↗
Figure 5
Figure 5. Figure 5: Alternating transient regimes during relaxation to the clus￾tered equilibrium. Numerical solution of the PDE (2) for σ = 0.65 starting in ρ0 ∼ 1 2 (Nper(−0.5, 0.2 2 ) + Nper(0.5, 0.2 2 )). (a) Trajectory snapshots of the particle density ρ(x, t), (b) evolution of peak height over time, (c) evolution of free energy over time. The dynamics shows both diffusion-dominated regimes (1.0 ≲ t ≲ 4.7) as well as agg… view at source ↗
Figure 6
Figure 6. Figure 6: Aggregation-dominated relaxation to the clustered equi￾librium. Numerical solution of the PDE (2) for σ = 0.838 starting in ρ0 ∼ Nper(0, 0.4 2 ). (a) Trajectory snapshots of the particle density ρ(x, t), (b) evolution of peak height over time, (c) evolution of free energy over time. With d dtFint < 0 and d dtFent > 0 for almost all t (except a short cooperative regime in the beginning for t ≲ 0.08), the dy… view at source ↗
Figure 7
Figure 7. Figure 7: Diffusion-dominated relaxation to the homogeneous equilib￾rium (but non-monotonicity in peak height). Numerical solution of the PDE (2) for σ = 0.838 starting in ρ0 ∼ Nper(0, 0.6 2 ). (a) Trajectory snap￾shots of the particle density ρ(x, t), (b) evolution of peak height over time, (c) evolution of free energy over time. With d dtFint > 0 and d dtFent < 0 for all t, the dynamics are purely diffusion-domina… view at source ↗
Figure 8
Figure 8. Figure 8: Time evolution of the second moment (left panel) and transition path [PITH_FULL_IMAGE:figures/full_fig_p016_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Hegselmann–Krause: Alternating transient regimes during relaxation to the homogeneous equilibrium. Numerical solution of the PDE (2) for U(x) = 1 2 (|x| 2 −1)1[0,R0](|x|) with R0 = 0.5 and σ = 0.485 starting in ρ0 ∼ Nper(0, 0.5 2 ). (a) Trajectory snapshots of the particle density ρ(x, t), (b) evolution of peak height over time, (c) evolution of free energy over time. The dynamics shows both aggregation-do… view at source ↗
read the original abstract

We investigate transient clustering dynamics in nonlocal aggregation-diffusion systems from an energetic perspective. Starting from a stochastic interacting particle system, we study the associated macroscopic McKean-Vlasov equation on the torus and exploit its Wasserstein gradient-flow structure to analyse the thermodynamic competition between interaction-driven aggregation and entropy-driven diffusion. Through numerical experiments for locally attractive interaction kernels, we identify alternating aggregation- and diffusion-dominated transient regimes along trajectories converging to fixed equilibria. These dynamics can be interpreted as a form of non-monotone clustering behaviour. Moreover, we demonstrate that clustering observables, such as the density peak height, are only partially coupled to the underlying energetic mechanisms and therefore do not uniquely characterise the relevant macroscopic transport dynamics. Our results highlight the role of the variational structure not only for equilibrium analysis, but also as a framework for understanding transient clustering phenomena in interacting particle systems.

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.

Referee Report

1 major / 0 minor

Summary. The manuscript investigates transient clustering dynamics in nonlocal aggregation-diffusion systems modeled by the McKean-Vlasov equation on the torus. It exploits the Wasserstein gradient-flow structure to analyze the competition between interaction-driven aggregation and entropy-driven diffusion. Numerical experiments with locally attractive interaction kernels identify alternating aggregation- and diffusion-dominated transient regimes along trajectories to equilibria (interpreted as non-monotone clustering), and demonstrate that observables such as density peak height are only partially coupled to the underlying energetic mechanisms.

Significance. If the numerical observations are robust, the work extends the application of variational structures from equilibrium analysis to transient dynamics, offering a framework for understanding non-monotone clustering in interacting particle systems. The energetic perspective on the partial decoupling of clustering observables from energy dissipation is a potentially useful contribution to the study of aggregation-diffusion models.

major comments (1)
  1. [Abstract] Abstract (numerical experiments paragraph): the central observational claims rest on numerical identification of alternating regimes and partial coupling of observables, yet no details are provided on discretization scheme, time-stepping, error controls, kernel definitions, or quantitative thresholds for regime identification. This renders the claims unverifiable from the given information and is load-bearing for the main results.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their detailed review and valuable feedback on our manuscript. We appreciate the recognition of the potential significance of our work on transient clustering dynamics. Below, we provide a point-by-point response to the major comment.

read point-by-point responses

  1. Referee: [Abstract] Abstract (numerical experiments paragraph): the central observational claims rest on numerical identification of alternating regimes and partial coupling of observables, yet no details are provided on discretization scheme, time-stepping, error controls, kernel definitions, or quantitative thresholds for regime identification. This renders the claims unverifiable from the given information and is load-bearing for the main results.

    Authors: We agree with the referee that the abstract's description of the numerical experiments lacks the necessary details on the discretization, time-stepping, and regime identification criteria, which are essential for verifying the central claims. While these details are elaborated in the main text (specifically in the numerical methods section), we acknowledge that the abstract should be self-contained to a greater extent for such observational results. We will revise the abstract to incorporate a brief mention of the numerical scheme (e.g., finite volume discretization with IMEX time stepping), the kernel specifications, and the quantitative criteria used for identifying aggregation- vs. diffusion-dominated regimes (based on the relative contributions to the energy dissipation). This revision will be made in the next version of the manuscript. revision: yes

Circularity Check

0 steps flagged

No significant circularity

full rationale

The paper's central claims rest on numerical experiments performed on the standard McKean-Vlasov equation and its established Wasserstein gradient-flow structure, which is drawn from prior literature rather than derived or fitted within this work. The identification of alternating aggregation/diffusion regimes and the partial coupling of clustering observables to energetic mechanisms are presented as direct observational results from those simulations, without any reduction of the reported dynamics to quantities defined by self-citations, fitted parameters renamed as predictions, or ansatzes smuggled via the authors' own prior work. The energetic interpretation simply invokes the known competition between interaction and entropy terms. No load-bearing step in the derivation chain reduces to its own inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Abstract-only review; the paper invokes the McKean-Vlasov equation and its Wasserstein gradient-flow structure as the foundation for energetic analysis. No free parameters, invented entities, or additional ad-hoc axioms are mentioned.

axioms (1)
  • domain assumption The macroscopic dynamics are governed by the McKean-Vlasov equation possessing a Wasserstein gradient-flow structure that encodes the competition between interaction energy and entropy.
    Stated in the abstract as the starting point for the energetic perspective and numerical study.

pith-pipeline@v0.9.1-grok · 5676 in / 1400 out tokens · 31973 ms · 2026-06-29T00:13:54.071350+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

34 extracted references · 21 canonical work pages · 4 internal anchors

  1. [1]

    Z. P. Adams, M. Engel, and R. S. Gvalani. Separation of time scales in weakly interacting diffusions. Arch. Rational Mech. Anal., 250 0 (3): 0 33, 2026. doi:10.1007/s00205-026-02180-w

    work page doi:10.1007/s00205-026-02180-w 2026

  2. [2]

    Ambrosio, N

    L. Ambrosio, N. Gigli, and G. Savaré. Gradient Flows. In Metric Spaces and in the Space of Probability Measures. Lectures in Mathematics. ETH Zürich. Birkhäuser Basel, 2nd edition, 2008

    2008

  3. [3]

    Arnold and J

    A. Arnold and J. Erb. Sharp entropy decay for hypocoercive and non-symmetric F okker- P lanck equations with linear drift. ASC Report (TU Wien), 29: 0 1--45, 2014. URL http://hdl.handle.net/20.500.12708/28336

    2014

  4. [4]

    Bailo, J

    R. Bailo, J. A. Carrillo, and J. Hu. Fully discrete positivity-preserving and energy-dissipating schemes for aggregation-diffusion equations with a gradient-flow structure. Commun. Math. Sci., 18 0 (5): 0 1259--1303, 2020. doi:10.4310/CMS.2020.v18.n5.a5

    work page doi:10.4310/cms.2020.v18.n5.a5 2020

  5. [5]

    Bellomo, D

    N. Bellomo, D. Clarke, L. Gibelli, P. Townsend, and B. Vreugdenhil. Human behaviours in evacuation crowd dynamics: F rom modelling to ``big data''' toward crisis management. Phys. Life Rev., 18: 0 1--21, 2016

    2016

  6. [6]

    J. A. Carrillo, M. Fornasier, G. Toscani, and F. Vecil. Particle, kinetic, and hydrodynamic models of swarming, pages 297--336. Birkh \"a user Boston, Boston, 2010. doi:10.1007/978-0-8176-4946-3_12

    work page doi:10.1007/978-0-8176-4946-3_12 2010

  7. [7]

    J. A. Carrillo, A. Chertock, and Y. Huang. A finite-volume method for nonlinear nonlocal equations with a gradient flow structure. Commun. Comput. Phys., 17 0 (1): 0 233--258, 2015. doi:10.4208/cicp.160214.010814a

    work page doi:10.4208/cicp.160214.010814a 2015

  8. [8]

    J. A. Carrillo, K. Craig, and Y. Yao. Aggregation--Diffusion Equations: D ynamics, Asymptotics, and Singular Limits , pages 65--108. Springer International Publishing, Cham, 2019. doi:10.1007/978-3-030-20297-2_3

    work page doi:10.1007/978-3-030-20297-2_3 2019

  9. [9]

    J. A. Carrillo, R. Gvalani, G. Pavliotis, and A. Schlichting. Long-time behaviour and phase transitions for the McKean -- Vlasov equation on the torus. Arch. Ration. Mech. Anal., 235 0 (1): 0 635--690, 2020. doi:10.1007/s00205-019-01430-4

    work page doi:10.1007/s00205-019-01430-4 2020

  10. [10]

    Chayes and V

    L. Chayes and V. Panferov. The McKean--Vlasov equation in finite volume. J. Stat. Phys., 138 0 (1): 0 351--380, 2010

    2010

  11. [11]

    D. A. Dawson. Critical dynamics and fluctuations for a mean-field model of cooperative behavior. J. Stat. Phys., 31 0 (1): 0 29--85, 1983

    1983

  12. [12]

    D'Orsogna, Y.-L

    M. D'Orsogna, Y.-L. Chuang, A. Bertozzi, and L. Chayes. Self-propelled particles with soft-core interactions: P atterns, stability, and collapse. Phys. Rev. Lett., 96: 0 104302, 2006. doi:10.1103/PhysRevLett.96.104302

    work page doi:10.1103/physrevlett.96.104302 2006

  13. [13]

    Garnier, G

    J. Garnier, G. Papanicolaou, and T.-W. Yang. Consensus convergence with stochastic effects. Vietnam J. Math., 45: 0 51--75, 2017. doi:10.1007/s10013-016-0190-2

    work page doi:10.1007/s10013-016-0190-2 2017

  14. [14]

    G \"a rtner

    J. G \"a rtner. On the McKean‐Vlasov limit for interacting diffusions. Math. Nachr., 137 0 (1): 0 197--248, 1988. doi:10.1002/mana.19881370116

    work page doi:10.1002/mana.19881370116 1988

  15. [15]

    Formation of clusters and coarsening in weakly interacting diffusions

    N. Gerber, R. Gvalani, M. Hairer, G. Pavliotis, and A. Schlichting. Formation of clusters and coarsening in weakly interacting diffusions. Working paper or preprint, 2025. URL https://arxiv.org/abs/2510.17629

    work page internal anchor Pith review Pith/arXiv arXiv 2025

  16. [16]

    Guillin, P

    A. Guillin, P. L. Bris, and P. Monmarché. Uniform in time propagation of chaos for the 2D vortex model and other singular stochastic systems. J. Eur. Math. Soc., 27 0 (6): 0 2359--2386, 2025

    2025

  17. [17]

    Hegselmann and U

    R. Hegselmann and U. Krause. Opinion dynamics and bounded confidence: Models , analysis and simulation. J. Artif. Soc. Soc. Simul., 5 0 (3), 2002

    2002

  18. [18]

    Helfmann, N

    L. Helfmann, N. Djurdjevac Conrad, A. Djurdjevac, S. Winkelmann, and C. Sch \"u tte. From interacting agents to density-based modeling with stochastic PDEs . Commun. Appl. Math. Comput. Sci., 16 0 (1): 0 1--32, 2021. doi:10.2140/camcos.2021.16.1

    work page doi:10.2140/camcos.2021.16.1 2021

  19. [19]

    Jabin and Z

    P.-E. Jabin and Z. Wang. Quantitative estimates of propagation of chaos for stochastic systems with W ^ -1, kernels. Invent. Math., 214 0 (1): 0 523--591, 2018. doi:10.1007/s00222-018-0808-y

    work page doi:10.1007/s00222-018-0808-y 2018

  20. [20]

    Jordan, D

    R. Jordan, D. Kinderlehrer, and F. Otto. The variational formulation of the Fokker--Planck equation. SIAM J. Math. Anal., 29 0 (1): 0 1--17, 1998. doi:10.1137/S0036141096303359

    work page doi:10.1137/s0036141096303359 1998

  21. [21]

    Leimkuhler, R

    B. Leimkuhler, R. Lohmann, G. A. Pavliotis, and P. A. Whalley. Cluster formation in diffusive systems. Working paper or preprint, 2025. URL https://arxiv.org/abs/2510.25034

    work page internal anchor Pith review arXiv 2025

  22. [22]

    A. J. Leverentz, C. M. Topaz, and A. J. Bernoff. Asymptotic dynamics of attractive--repulsive swarms. SIAM J. Appl. Dyn. Syst., 8 0 (3): 0 880--908, 2009. doi:10.1137/090749037

    work page doi:10.1137/090749037 2009

  23. [23]

    P. Li, Z. Li, H. Zhang, and J. Bian. On the generalization properties of diffusion models. Adv. Neural Inf. Process. Syst., 36: 0 2097--2127, 2023

    2097

  24. [24]

    Neureither and C

    L. Neureither and C. Hartmann. Time scales and exponential trend to equilibrium: G aussian model problems. In International workshop on Stochastic Dynamics out of Equilibrium, pages 391--410. Springer, 2017

    2017

  25. [25]

    Rouwhorst, C

    J. Rouwhorst, C. Ness, S. Stoyanov, A. Zaccone, and P. Schall. Nonequilibrium continuous phase transition in colloidal gelation with short-range attraction. Nat. Commun., 11 0 (1): 0 3558, 2020

    2020

  26. [26]

    Santambrogio

    F. Santambrogio. Optimal Transport for Applied Mathematicians. Calculus of Variations, PDEs and Modeling. Progress in Nonlinear Differential Equations and Their Applications. Birkhäuser Cham, 1st edition, 2015

    2015

  27. [27]

    Sznitman

    A.-S. Sznitman. Topics in propagation of chaos. In Ecole d'Et \'e de Probabilit \'e s de Saint-Flour XIX --- 1989 , pages 165--251. Springer Berlin Heidelberg, 1991

    1989

  28. [28]

    Taboada-Serrano, C.-J

    P. Taboada-Serrano, C.-J. Chin, S. Yiacoumi, and C. Tsouris. Modeling aggregation of colloidal particles. Curr. Opin. Colloid Interface Sci., 10 0 (3-4): 0 123--132, 2005

    2005

  29. [29]

    C. M. Topaz, A. J. Bernoff, S. Logan, and W. Toolson. A model for rolling swarms of locusts. Eur. Phys. J. Spec. Top., 157: 0 93--109, 2008. doi:10.1140/epjst/e2008-00633-y

    work page doi:10.1140/epjst/e2008-00633-y 2008

  30. [30]

    Wehlitz, M

    N. Wehlitz, M. Sadeghi, A. Montefusco, C. Sch\" u tte, G. A. Pavliotis, and S. Winkelmann. Approximating particle-based clustering dynamics by stochastic PDE s. SIAM J. Appl. Dyn. Syst., 24 0 (2): 0 1231--1250, 2025. doi:10.1137/24M1676661

    work page doi:10.1137/24m1676661 2025

  31. [31]

    Data-driven Reduction of Transfer Operators for Particle Clustering Dynamics

    N. Wehlitz, G. A. Pavliotis, C. Schütte, and S. Winkelmann. Data-driven reduction of transfer operators for particle clustering dynamics. Working paper or preprint, 2026. URL https://arxiv.org/abs/2601.02932

    work page internal anchor Pith review Pith/arXiv arXiv 2026

  32. [32]

    S. Yang, W. Chu, and P. G. Kevrekidis. Noisy nonlocal aggregation model with gradient flow structures. Working paper or preprint, 2026. URL https://arxiv.org/abs/2602.03654

    work page arXiv 2026

  33. [33]

    Zimper, A

    S. Zimper, A. Djurdjevac, C. Hartmann, C. Schütte, and N. Djurdjevac Conrad . Mean-field optimal control with stochastic leaders. Working paper or preprint, 2025. URL https://arxiv.org/abs/2512.19201

    work page arXiv 2025

  34. [34]

    Clustering in co-evolving opinion dynamics: reduced SPDE models

    S. Zimper, N. Djurdjevac Conrad , F. Cornalba, and A. Djurdjevac. Clustering in co-evolving opinion dynamics: reduced SPDE models. Working paper or preprint, 2026. URL https://arxiv.org/abs/2604.27961

    work page internal anchor Pith review Pith/arXiv arXiv 2026