pith. sign in

arxiv: 2605.00110 · v1 · submitted 2026-04-30 · 🧮 math.AP

Formation and Behavior of Dirac Singularities in the Parabolic-Elliptic Keller-Segel System in Dimensions ngeq 3

Gregor Fl\"uchter This is my paper

Pith reviewed 2026-05-09 20:37 UTC · model grok-4.3

classification 🧮 math.AP
keywords Keller-Segel systemDirac singularitymeasure-valued solutionsblow-upparabolic-ellipticradial symmetrymass concentration
0
0 comments X

The pith

Solutions of the parabolic-elliptic Keller-Segel system in dimensions n at least 3 develop a Dirac mass at the origin continuously after blow-up, with all mass eventually absorbed there.

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

The paper considers nonnegative radially symmetric solutions of the Keller-Segel system in a ball under Neumann boundary conditions. It extends the notion of solution past finite-time blow-up by constructing a minimal measure-valued solution that includes a Dirac delta at the origin whose mass θ(t) is increasing and right-continuous. For certain initial data θ(t) becomes positive at some finite time t0, grows strictly thereafter, and the full mass is asymptotically concentrated at the origin. This behavior differs from the two-dimensional case, where the Dirac mass appears instantaneously.

Core claim

Nonnegative radially symmetric solutions of the parabolic-elliptic Keller-Segel system in dimensions n≥3 admit a minimal extension as measure-valued solutions u(t)=θ(t)δ_0 + ρ(·,t) dx for t≥0, where ρ is integrable and θ is increasing and right-continuous. There exist initial data for which θ(t0)>0 at some t0>0, θ increases strictly on [t0,∞), and the entire mass is asymptotically absorbed at the origin.

What carries the argument

The minimal measure-valued solution u(t)=θ(t)δ_0 + ρ(·,t) dx, where the increasing function θ(t) tracks the accumulated singular mass at the origin.

If this is right

  • A class of initial data produces a positive Dirac mass at the origin after finite time.
  • The singular mass θ(t) increases continuously and strictly after its first appearance.
  • The full initial mass concentrates at the Dirac singularity as time tends to infinity.
  • Formation of the singularity occurs gradually rather than by an instantaneous jump.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The gradual accumulation may reflect different scaling of the aggregation mechanism compared with two dimensions.
  • Numerical methods for the system in higher dimensions may need to track measure concentrations explicitly beyond the classical blow-up time.
  • The construction might extend to other boundary conditions or to systems with additional nonlinearities that preserve radial symmetry.

Load-bearing premise

Solutions are nonnegative and radially symmetric under homogeneous Neumann boundary conditions in a ball, allowing concentration specifically at the origin.

What would settle it

A radially symmetric nonnegative solution that blows up in finite time but for which the associated measure-valued extension has θ(t)=0 for all t, or for which θ(t) does not absorb the full mass asymptotically.

read the original abstract

We consider nonnegative radially symmetric solutions of the parabolic-elliptic Keller-Segel system \begin{align*} \left\lbrace \begin{array}{r@{}l@{\quad}l} &u_t=\Delta u-\nabla \cdot \big(u\nabla v\big),\\ &0=\Delta v -\mu + u , \\ \end{array}\right. \end{align*} where $\mu$ is the spatial average of $u$, under homogeneous Neumann boundary conditions in a ball in $\mathbb R^n$ for $n\geq 3$. In two dimensions, it is well established that solutions blowing up in finite time converge to a Dirac profile in the vague topology. In contrast, for $n\geq 3$, blow-up solutions with finite existence time do not appear to exhibit such concentration behavior. By generalizing to measure-valued solutions corresponding to accumulated densities of $u$, we extend the analysis beyond the blow-up time. Within this framework, we establish the existence of a minimal solution \[ u(t)=\theta(t)\delta_0 + \rho(\cdot,t) dx, \qquad t \geq 0, \] where $\rho$ is integrable and $\theta$ is increasing and right-continuous. We further construct a class of initial data for which $\theta(t_0)>0$ for some $t_0>0$, thereby establishing the formation of a Dirac mass at the origin. Unlike in the case $n=2$, the singular mass does not jump to a positive level instantaneously; instead, $\theta$ becomes positive continuously. Moreover, $\theta$ is strictly increasing on $[t_0,\infty)$, and the entire mass is asymptotically absorbed at the origin.

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

2 major / 2 minor

Summary. The manuscript considers nonnegative radially symmetric solutions of the parabolic-elliptic Keller-Segel system in a ball in R^n (n≥3) under homogeneous Neumann boundary conditions. It extends the analysis past finite-time blow-up by constructing minimal measure-valued solutions of the form u(t)=θ(t)δ_0 + ρ(·,t) dx, where ρ is integrable and θ is increasing and right-continuous. For a class of initial data, the paper proves that θ(t0)>0 for some finite t0>0 (with continuous onset), that θ is strictly increasing on [t0,∞), and that the entire mass is asymptotically absorbed at the origin. This contrasts with the instantaneous Dirac formation known in 2D.

Significance. If the central claims hold, the work supplies a precise description of Dirac singularity formation and long-time mass absorption in higher-dimensional Keller-Segel models, where standard blow-up solutions do not concentrate to a Dirac measure in the usual sense. The minimal measure-valued solution framework is a technically useful extension that permits continuation beyond blow-up time and yields falsifiable predictions on the monotonicity and asymptotics of θ(t).

major comments (2)
  1. [Main existence theorem and construction of minimal solution] The existence proof for the minimal measure-valued solution u(t)=θ(t)δ_0 + ρ dx and the verification that it satisfies the system in the weak measure sense (particularly consistency of the elliptic equation for v with the accumulated density) are load-bearing for the main theorem; the provided abstract states the result but the detailed construction and consistency checks with the original PDE are not fully verifiable from the given material.
  2. [Section on properties of θ(t) and asymptotic behavior] The claim that θ becomes positive continuously at t0 (rather than jumping) and is strictly increasing on [t0,∞) with full asymptotic absorption relies on the radial symmetry and Neumann boundary conditions; these assumptions are central but require explicit confirmation that no mass escapes or that the comparison principle used for minimality holds uniformly past blow-up.
minor comments (2)
  1. [Abstract] The abstract equation block would be clearer if the definition of μ (spatial average of u) were recalled explicitly in the displayed system.
  2. [Introduction] Notation for the vague topology and the precise sense in which the measure-valued solution satisfies the parabolic equation should be stated once in the introduction for readers unfamiliar with the 2D literature.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and insightful comments on our manuscript. We address each major comment point by point below, providing clarifications and indicating revisions made to strengthen the presentation.

read point-by-point responses

  1. Referee: The existence proof for the minimal measure-valued solution u(t)=θ(t)δ_0 + ρ dx and the verification that it satisfies the system in the weak measure sense (particularly consistency of the elliptic equation for v with the accumulated density) are load-bearing for the main theorem; the provided abstract states the result but the detailed construction and consistency checks with the original PDE are not fully verifiable from the given material.

    Authors: The construction of the minimal measure-valued solution is carried out in detail in Section 3 via a limiting procedure applied to a family of regularized problems, with the weak formulation verified directly in the limit. The consistency of the elliptic equation is established by showing that the potential v satisfies the Poisson equation in the sense of distributions, where the Dirac contribution enters explicitly through integration by parts against test functions. To improve verifiability, we have expanded the relevant subsection with additional intermediate estimates and an explicit computation of the limit passage for the elliptic part, including a new remark confirming that no extraneous boundary or singular terms arise. revision: yes

  2. Referee: The claim that θ becomes positive continuously at t0 (rather than jumping) and is strictly increasing on [t0,∞) with full asymptotic absorption relies on the radial symmetry and Neumann boundary conditions; these assumptions are central but require explicit confirmation that no mass escapes or that the comparison principle used for minimality holds uniformly past blow-up.

    Authors: Radial symmetry together with the homogeneous Neumann boundary conditions ensures that the total mass is conserved in the measure-valued sense, as the radial flux vanishes at the boundary for all approximating solutions and passes to the limit. The comparison principle underlying minimality is proved in Lemma 4.2 for the regularized problems and extends uniformly past blow-up by the definition of the minimal solution as the pointwise infimum. We have revised Section 4 to include an explicit paragraph confirming mass conservation (no escape) and the uniform validity of the comparison, together with a short argument showing that the continuous onset of θ follows from the strict positivity of the regular density near the origin at the first blow-up time. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation self-contained

full rationale

The paper derives existence of minimal measure-valued solutions u(t) = θ(t)δ_0 + ρ dx and constructs initial data leading to Dirac mass formation directly from the parabolic-elliptic Keller-Segel PDE system under radial symmetry and Neumann conditions. The claims follow from analysis and extension beyond blow-up without fitting parameters, self-definitional reductions, or load-bearing self-citations; θ(t) properties emerge from the equations and construction rather than being presupposed.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 1 invented entities

The central claim rests on the parabolic-elliptic Keller-Segel equations with radial symmetry and the extension to measure-valued solutions; no free parameters are fitted to data.

axioms (2)
  • domain assumption Solutions are nonnegative and radially symmetric with homogeneous Neumann boundary conditions in a ball in R^n.
    Assumed throughout to allow concentration at the origin and simplify the analysis.
  • domain assumption The parabolic-elliptic system holds classically before blow-up and extends via measures after.
    Core modeling choice enabling continuation past finite existence time.
invented entities (1)
  • Minimal measure-valued solution with time-dependent Dirac mass θ(t) at the origin no independent evidence
    purpose: To describe the accumulated density and extend the solution beyond the blow-up time while satisfying the system in a weak sense.
    Constructed in the paper as the key object whose existence and properties are proved.

pith-pipeline@v0.9.0 · 5618 in / 1657 out tokens · 73413 ms · 2026-05-09T20:37:36.690016+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

30 extracted references · 30 canonical work pages

  1. [1]

    Bai and M

    X. Bai and M. Zhou. On the blow-up profile of keller–segel–patlak system.Mathematische Annalen, 392:313 – 337,

    work page

  2. [2]

    URL:https://api.semanticscholar.org/CorpusID:276091231

    work page

  3. [3]

    Bellomo, A

    N. Bellomo, A. Bellouquid, Y. Tao, and M. Winkler. Toward a mathematical theory of keller–segel models of pattern formation in biological tissues.Mathematical Models and Methods in Applied Sciences, 25(09):1663–1763, 2015. arXiv:https://doi.org/10.1142/S021820251550044X,doi:10.1142/S021820251550044X

    work page doi:10.1142/s021820251550044x 2015

  4. [4]

    Bellomo, N

    N. Bellomo, N. Outada, J. Soler, Y. Tao, and M. Winkler. Chemotaxis and cross-diffusion models in complex environments: models and analytic problems toward a multiscale vision.Math. Models Methods Appl. Sci., 32(4):713– 792, 2022.doi:10.1142/S0218202522500166

    work page doi:10.1142/s0218202522500166 2022

  5. [5]

    Finite-timeblow-upinadegeneratechemotaxissystemwithfluxlimitation.Transactions of the American Mathematical Society, Series B, 4:31–67, 06 2017.doi:10.1090/btran/17

    N.BellomoandM.Winkler. Finite-timeblow-upinadegeneratechemotaxissystemwithfluxlimitation.Transactions of the American Mathematical Society, Series B, 4:31–67, 06 2017.doi:10.1090/btran/17

    work page doi:10.1090/btran/17 2017

  6. [6]

    P. Biler. Existence and nonexistence of solutions for a model of gravitational interaction of particles iii.Colloquium Mathematicum, 68:229–239, 01 1995.doi:10.4064/cm-68-2-229-239

    work page doi:10.4064/cm-68-2-229-239 1995

  7. [7]

    P. Biler. Local and global solvability of some systems modelling chemotaxis.Advances in Mathematical Sciences and Applications, 8, 01 1998

    work page 1998

  8. [8]

    P. Biler. Radially symmetric solutions of a chemotaxis model in the plane - the supercritical case.Banach Center Publications, 81:31–42, 01 2008.doi:10.4064/bc81-0-2

    work page doi:10.4064/bc81-0-2 2008

  9. [9]

    Blanchet, J

    A. Blanchet, J. Dolbeault, and B. Perthame. Two-dimensional keller-segel model: Optimal critical mass and quali- tative properties of the solutions, in "electron.Electronic Journal of Differential Equations, 44, 04 2006

    work page 2006

  10. [10]

    Cieślak and M

    T. Cieślak and M. Winkler. Finite-time blow-up in a quasilinear system of chemotaxis.Nonlinearity, 21:1057 – 1076,

    work page

  11. [11]

    URL:https://api.semanticscholar.org/CorpusID:120503706

    work page

  12. [12]

    Daners and P

    D. Daners and P. K. Medina. Abstract evolution equations, periodic problems and applications. 1992. URL: https://api.semanticscholar.org/CorpusID:116694258. 41

    work page 1992

  13. [13]

    Methods Appl

    K. Fujie, M. Winkler, and T. Yokota. Boundedness of solutions to parabolic–elliptic keller–segel systems with signal-dependent sensitivity.Mathematical Methods in the Applied Sciences, 38(6):1212–1224, 2015. URL:https://onlinelibrary.wiley.com/doi/abs/10.1002/mma.3149,arXiv:https://onlinelibrary.wiley.com/ doi/pdf/10.1002/mma.3149,doi:10.1002/mma.3149

    work page doi:10.1002/mma.3149 2015

  14. [14]

    Herrero and Juan J

    M.HerreroandJ.Velázquez. Singularitypatternsinachemotaxismodelmath.Mathematische Annalen, 306:583–623, 09 1996.doi:10.1007/BF01445268

    work page doi:10.1007/bf01445268 1996

  15. [15]

    M. A. Herrero, E. Medina, and J. J. L. Velázquez. Finite-time aggregation into a single point in a reaction-diffusion system.Nonlinearity, 10(6):1739–1754, 1997.doi:10.1088/0951-7715/10/6/016

    work page doi:10.1088/0951-7715/10/6/016 1997

  16. [16]

    Horstmann

    D. Horstmann. From 1970 until present: the Keller-Segel model in chemotaxis and its consequences I.Jahresbericht der Deutschen Mathematiker-Vereinigung, 105(3):103–165, 2003

    work page 1970

  17. [17]

    Jäger and S

    W. Jäger and S. Luckhaus. On explosions of solutions to a system of partial differential equations modelling chemo- taxis.Transactions of the American Mathematical Society, 329:819, 1992.doi:10.2307/2153966

    work page doi:10.2307/2153966 1992

  18. [18]

    N. I. Kavallaris and P. Souplet. Grow-up rate and refined asymptotics for a two-dimensional patlak–keller–segel model in a disk.SIAM Journal on Mathematical Analysis, 40(5):1852–1881, 2009.arXiv:https://doi.org/10. 1137/080722229,doi:10.1137/080722229

    work page doi:10.1137/080722229 2009

  19. [19]

    E. F. Keller and L. A. Segel. Initiation of slime mold aggregation viewed as an instability.Journal of Theoretical Biology, 26:399–415, 1970.doi:10.1016/0022-5193(70)90092-5

    work page doi:10.1016/0022-5193(70)90092-5 1970

  20. [20]

    O. A. Ladyzhenskaja, V. A. Solonnikov, and N. N. Ural’ceva.Linear and Quasilinear Equations of Parabolic Type, volume 23 ofTranslations of Mathematical Monographs. American Mathematical Society, Providence, RI, 1968

    work page 1968

  21. [21]

    T. Nagai. Blow-up of radially symmetric solutions to a chemotaxis system.Adv. Math. Sci. Appl., 5(2):581–601, 1995

    work page 1995

  22. [22]

    Nanjundiah

    V. Nanjundiah. Chemotaxis, signal relaying and aggregation morphology.Journal of Theoretical Biology, 42(1):63– 105, 1973. Conjecture that Dirac singularities can form in KS on page 90. URL:https://www.sciencedirect.com/ science/article/pii/0022519373901495,doi:10.1016/0022-5193(73)90149-5

    work page doi:10.1016/0022-5193(73)90149-5 1973

  23. [23]

    Osaki and A

    K. Osaki and A. Yagi. Finite dimensional attractor for one-dimensional Keller-Segel equations.Funkcialaj Ekvacioj, Volume 44, No.3, pp. 441-470, 2001

    work page 2001

  24. [24]

    T. Senba. Chemotactic collapse of radial solutions to jäger-luckhaus system.Advances in Mathematical Sciences and Applications, 14, 01 2004

    work page 2004

  25. [25]

    T. Senba. Blowup behavior of radial solutions to jäger-luckhaus system in high dimensional domains.Funkcialaj Ekvacioj-serio Internacia - FUNKC EKVACIOJ-SER INT, 48:247–271, 08 2005.doi:10.1619/fesi.48.247

    work page doi:10.1619/fesi.48.247 2005

  26. [26]

    Senba and T

    T. Senba and T. Suzuki. Chemotactic collapse in a parabolic-elliptic system of mathematical biology.Advances in Differential Equations, 6, 01 2001.doi:10.57262/ade/1357141500

    work page doi:10.57262/ade/1357141500 2001

  27. [27]

    Souplet and M

    P. Souplet and M. Winkler. Blow-up profiles for the parabolic-elliptic keller-segel system in dimensionsn≥3. Communications in Mathematical Physics, 367, 02 2018.doi:10.1007/s00220-018-3238-1

    work page doi:10.1007/s00220-018-3238-1 2018

  28. [28]

    Finite-timeblow-upinthehigher-dimensionalparabolic–parabolickeller–segelsystem.Journal de Mathé- matiques Pures et Appliquées, 100(5):748–767, 2013

    M.Winkler. Finite-timeblow-upinthehigher-dimensionalparabolic–parabolickeller–segelsystem.Journal de Mathé- matiques Pures et Appliquées, 100(5):748–767, 2013. URL:https://www.sciencedirect.com/science/article/pii/ S0021782413000287,doi:10.1016/j.matpur.2013.01.020

    work page doi:10.1016/j.matpur.2013.01.020 2013

  29. [29]

    M. Winkler. How unstable is spatial homogeneity in Keller-Segel systems? A new critical mass phenomenon in two- and higher-dimensional parabolic-elliptic cases.Math. Ann., 373(3-4):1237–1282, 2019.doi:10.1007/ s00208-018-1722-8

    work page 2019

  30. [30]

    Winkler and K

    M. Winkler and K. C. Djie. Boundedness and finite-time collapse in a chemotaxis system with volume-filling effect. Nonlinear Analysis: Theory, Methods & Applications, 72(2):1044–1064, 2010. URL:https://www.sciencedirect. com/science/article/pii/S0362546X09009389,doi:10.1016/j.na.2009.07.045. 42

    work page doi:10.1016/j.na.2009.07.045 2010