A mathematical analysis of a radially symmetric reaction–diffusion model finds that its modeled population can follow three sharply different long-term paths, depending on a boundary-density parameter δ: the range can expand without limit, settle at a finite size, or shrink toward zero. The result comes from a theorem-driven study of a deterministic model, not from observations of a real species.
The model tracks a species density u(t,r) across a radial range written as [0,h(t)]. Its free boundary is set by the density level δ, so changing δ changes the mathematical rule governing the edge of the modeled population. The analysis is carried out for spatial dimensions N≥2, with positive diffusion, initial radius and boundary parameter, and an initial function in the stated class.
Three outcomes at the boundary
When 0<δ<1, the model is in a spreading regime: h(t) tends to infinity, while the density approaches 1 locally and uniformly across every fixed radial region. In plain terms, the modeled range keeps advancing and the interior settles toward the model’s upper density level.
At the threshold δ=1, the model instead reaches a transition regime. The boundary approaches a finite positive limit, h∞, and the density approaches 1 uniformly throughout the changing interval [0,h(t)]. The range therefore remains bounded even as the density reaches the same limiting level used in the spreading case.
For δ>1, the outcome reverses: h(t) tends to 0, and the density approaches δ uniformly on [0,h(t)]. This is the model’s vanishing regime, in which the radial population range collapses toward the center.
A closer look at the advancing front
The analysis goes beyond classifying the three outcomes. In the spreading regime, the boundary speed h′(t) approaches a limiting speed c*, and the full density profile approaches a traveling shape, written q*(h(t)−r). A semi-wave is a profile that moves with the front while describing how density changes behind it; the study states that the relevant semi-wave pair is unique for 0<δ<1.
The front does not follow its leading linear motion exactly. Its position has a logarithmic correction: h(t)−c*t+c_N(δ)log t converges to a finite constant. That means the difference between the actual boundary and the main speed-times-time estimate changes systematically with log t, rather than disappearing as a fixed short-term error.
The coefficient c_N(δ) is positive and approaches d(N−1)/c0 as δ approaches 0. The authors identify that limit with the corresponding shift coefficient for the high-dimensional radial pushed-case Cauchy problem, linking the free-boundary calculation to that related model through a parameter limit.
What the result establishes—and what it does not
The first theorem establishes global existence and uniqueness: under the stated equations and initial-data assumptions, there is one solution pair (u,h) for every positive time. The local existence argument uses a contraction mapping to obtain a unique fixed point, while the long-time analysis uses upper and lower solutions together with zero-number theory to establish the boundary limit and the three-regime classification.
The study does not estimate a spreading speed from field or laboratory data, and it reports no statistical uncertainty or sample-based effect estimates. Its evidence supports conditional theorems for the specified monostable reaction term, radial geometry and initial-data class. Nonradial geometries and reaction terms outside that class are not addressed.
The δ→0 result is also a mathematical convergence statement. Under the stated C1 convergence of the initial data, the free-boundary solutions converge locally to the corresponding Cauchy-problem solution as δ approaches 0; the analysis does not directly study an original δ=0 version of the boundary condition.
A preprint with a narrowly defined scope
The document is identified in its front matter as an arXiv mathematics item. The authors report partial support from the Natural Science Foundation of China and the Natural Science Foundation of Shandong Province.
Paper data and sources
Original title: The high dimensional monostable reaction-diffusion equation with free boundary and radial symmetry
Authors: Hongkai Cao, Jingyi Cui, Chengzhe Tang, Xiaoyan Zhang
Journal/Repository: arXiv
Status: Preprint, not yet peer-reviewed
First online: 2026-08-26
DOI: Not available
Original paper · Full text