Preprint

Large diffusion sorts a nonlinear heat equation into three outcomes

Preprint: A mathematical analysis gives a uniform rule for decay, threshold behavior and finite-time blow-up as diffusion grows.

A new mathematical preprint reports that sufficiently strong diffusion, the model’s spatial-smoothing effect, can sort every initial profile in a bounded H1 ball into one of three outcomes: finite-time blow-up, meaning the solution becomes unbounded in finite time; exponential convergence to zero; or global convergence to a positive constant called ξa. H1 is the function space used to control a profile and its first spatial derivatives.

From one value to a spatial field

The benchmark is the spatially homogeneous scalar equation, where all points carry the same value. For reaction parameter a greater than 1, zero is stable and the positive equilibrium ξa is unstable. Values below ξa converge to zero, while values above it blow up in finite time. The full analysis uses this simple phase portrait as its reference.

For each radius R, the theorem identifies a sufficiently large diffusivity level. Once diffusion is above that level, every profile with H1 norm at most R falls into exactly one of the three categories. The proof also supplies an explicit sufficient threshold for this dynamical classification when R is at least 1, with a constant that depends only on a and the spatial domain.

One part of the result turns the classification into a test based on the average. In the large-diffusion regime, finite-time blow-up occurs if and only if the solution’s spatial mean, its average over the domain, exceeds ξa at some time before the time up to which the solution exists.

The cost of controlling concentration

The paper’s second question is how the diffusion cost changes as the allowed H1 radius grows. Fix a positive gap δ below ξa, so the initial mean is kept that distance from the positive equilibrium. For every positive R, the uniform stabilization threshold is itself positive and finite. As R becomes large, the logarithm of that threshold is R squared divided by 8 pi, up to a remainder of logarithmic order in R. The worst-case diffusion requirement therefore has an exponential scale set by the square of R, with polynomial corrections in the two-sided bounds.

The paper attributes this coefficient to the sharp mean-zero Moser-Trudinger inequality, a concentration estimate that controls the exponential e to the power of u through the function’s spatial mean and gradient energy. In this analysis, that control fixes the leading 1 over 8 pi in the logarithmic growth law.

The lower bound is built in the opposite direction. Boundary-concentrating Moser profiles, combined with a localized Kaplan eigenfunction argument, create the obstruction that prevents a smaller uniform threshold from working across the whole class. The Neumann boundary permits this concentration mechanism, so the same leading exponential scale appears as a worst-case lower bound.

For subthreshold data, the theorem gives a stronger conclusion than eventual bounded behavior. If the initial H1 norm is within R and the initial spatial mean is at least δ below ξa, the theorem states that diffusivity above the sufficient stabilization threshold yields a global solution that converges exponentially to zero in both the L-infinity and H1 norms. The statement is uniform over the specified class, so the threshold is set by the radius and the mean gap rather than by one selected profile.

The gap itself is essential. When R exceeds ξa times the square root of the domain’s area, the uniform threshold tends to infinity as δ shrinks toward zero. In that radius range, no fixed finite diffusion level can cover profiles whose mean approaches ξa from below.

The middle outcome is also genuinely populated. A constructed family contains spatially nonconstant initial data that converge to ξa, and within that family a unique transition parameter separates blow-up from decay. The threshold is therefore not just a formal possibility reserved for perfectly uniform profiles.

A result with a defined mathematical boundary

The document is an arXiv preprint, arXiv:2608.28061 version 1, dated 28 August 2026. It treats deterministic initial profiles in H1 of the spatial domain with positive diffusivity and their unique maximal solutions, so the result is a theorem about a specified mathematical class rather than an empirical dataset.

Paper data and sources

Original title: Large-diffusion dynamics for a planar Neumann heat equation with exponential nonlinearity
Authors: Juneyoung Seo
Journal/Repository: arXiv
Status: Preprint, not yet peer-reviewed
First online: 2026-08-28
DOI: Not available
Original paper · Full text

Versions and corrections

  1. Published automatically after legal-source, freshness, evidence, and independent-verification gates passed.