🔢 mathematics

Asymptotic large time behavior of singular entire solutions of the fast diffusion equation

This paper establishes the asymptotic large-time behavior of singular entire solutions to the fast diffusion equation in $\mathbb{R}^n \setminus \{0\}$ for $n \ge 3$ , demonstrating that initial data with singularities bounded by $|x|^{-\gamma}$ converge to a specific singular radially symmetric self-similar profile while also providing an asymptotic expansion of this profile near the origin.

Original authors: Kin Ming Hui, Jongmyeong Kim

Published 2026-08-14

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Kin Ming Hui, Jongmyeong Kim

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the universe as a giant, invisible ocean where invisible particles are constantly drifting, bumping into each other, and spreading out. Sometimes, they spread out slowly, like a drop of ink in cold water; other times, they rush out wildly, like a firecracker exploding in a room. Scientists call this movement "diffusion." When the particles are ordinary, like heat in a metal rod, they follow a predictable, gentle path. But when the particles interact in a special, tricky way—where the more crowded they get, the faster they scatter—the math gets wild. This is called "fast diffusion." It describes everything from how plasma behaves in stars to how impurities move through the silicon chips in your phone.

Now, picture a single point in this ocean where the density of particles is infinitely high, a "singularity." It's like a black hole for particles, a place where the rules seem to break down. For a long time, mathematicians have been trying to understand what happens to the flow of particles near these infinite points as time goes on. Do they settle into a calm, predictable pattern? Do they explode? Or do they dance in a chaotic, never-ending loop? Understanding this isn't just about abstract math; it helps us predict how things evolve over vast periods, from the cooling of the universe to the behavior of complex materials.

This paper dives deep into that chaotic dance. The authors, Kin Ming Hui and Jongmyeong Kim, are looking at a specific type of fast diffusion equation where the particles start with a very sharp, infinite spike at the center. They want to know: if you wait a very, very long time, does this messy, exploding mess of particles eventually settle down into a specific, recognizable shape?

Think of it like watching a chaotic crowd of people running away from a sudden noise in a stadium. At first, everyone is screaming and running in random directions. But if you zoom out and watch for hours, you might notice they start to form a specific, smooth wave pattern. The authors prove that for these "fast diffusion" equations, the chaotic particles do indeed settle into a very specific, elegant shape called a "self-similar solution." This shape is special because it looks the same no matter how much you zoom in or out, just like a fractal.

The paper does three main things. First, it proves that this specific, smooth shape actually exists and is the only one that fits the rules. Before this, other methods were used to find it, but the authors used a clever new trick called a "fixed point method." Imagine trying to find a specific spot on a map by folding the map over and over; eventually, the folds land on the exact same spot. That's the math they used to show that this unique shape is the only possible destination for the particles.

Second, they figured out exactly what happens right at the very center of the explosion, near the singularity. They calculated a detailed "asymptotic expansion," which is just a fancy way of saying they wrote down a precise recipe for how the particles behave as they get closer and closer to the infinite point. It's like having a high-definition map of the eye of a hurricane, showing exactly how the wind speed changes as you get right to the center.

Finally, they showed that if you start with a cloud of particles that roughly matches this infinite spike (but isn't perfect), and you let time run forward, the cloud will eventually morph and settle down to match that perfect, smooth shape they found earlier. They proved this by showing that the difference between the messy starting cloud and the perfect final shape shrinks and shrinks until it disappears, like a rough stone being smoothed by a river until it becomes a perfect pebble.

The authors are very sure about these findings. They didn't just guess or run computer simulations; they provided a rigorous mathematical proof. They ruled out the idea that the particles might behave in some other, weirder way or that there might be multiple different shapes they could settle into. Instead, they showed that the universe, in this specific mathematical scenario, has a very strict preference: it always wants to settle into this one unique, self-similar pattern. So, even if the beginning is a chaotic mess of infinite density, the long-term future is surprisingly orderly and predictable.

Technical Summary: Asymptotic Large Time Behavior of Singular Entire Solutions of the Fast Diffusion Equation

Problem Statement
This paper investigates the asymptotic large-time behavior of singular solutions to the fast diffusion equation (FDE):
$u_t = \Delta(u^m/m)$
defined in the domain $(\mathbb{R}^n \setminus \{0\}) \times (0, \infty)$ , where $n \ge 3$ and $0 < m < \frac{n-2}{n}$ . The study focuses on solutions that blow up at the origin for all $t > 0$ . Specifically, the authors analyze initial data $u_0(x)$ satisfying the growth condition $A_1|x|^{-\gamma} \le u_0(x) \le A_2|x|^{-\gamma}$ for constants $A_2 > A_1 > 0$ .

The paper addresses two primary regimes for the exponent $\gamma$ :

The existence and uniqueness of singular forward self-similar solutions of the form $V(x,t) = t^{-\alpha}f(t^{-\beta}x)$ .
The asymptotic convergence of general singular solutions to these self-similar profiles as $t \to \infty$ .

A significant portion of the work extends previous results (specifically those of Hui and Kim [23]) from the range $\frac{2}{1-m} < \gamma < n$ to the extended range $n \le \gamma < \frac{n-2}{m}$ .

Methodology
The authors employ a combination of fixed-point arguments, asymptotic analysis, and weighted $L^1$ -contraction techniques.

Fixed Point Method for Existence and Uniqueness:
To prove the existence of singular radially symmetric forward self-similar solutions, the authors reduce the partial differential equation to an ordinary differential equation (ODE) for the profile function $f$ . By introducing logarithmic variables ( $s = \log r$ ) and transforming the dependent variable, the problem is converted into finding a solution to a specific ODE system (equation 2.7). The authors construct a Banach space and define a mapping $\Phi$ on a closed subset $D_{b_1}$ . They demonstrate that this mapping is a contraction, thereby establishing the existence and uniqueness of the solution via the Contraction Mapping Theorem. This provides a new, direct proof of results previously obtained by Hui and Kim [23] using an inversion method.
Asymptotic Expansion near the Origin:
To characterize the behavior of the singular profile $f$ near the origin, the authors utilize a modification of techniques from Hui and Kim [24]. They define a rescaled function $w(\rho)$ and derive a second-order asymptotic expansion. By analyzing the derivatives of $w$ at $\rho=0$ , they establish that the solution can be extended to a $C^2$ function at the origin, providing explicit formulas for the first and second derivatives in terms of the parameters $m, n, \alpha, \beta$ .
Weighted $L^1$ -Contraction and Convergence:
To study the large-time behavior, the authors introduce a specific non-negative radially symmetric weight function $\phi_\mu$ (constructed in Section 1) which satisfies $\Delta \phi_\mu \le 0$ . They prove a weighted $L^1$ -contraction principle (Theorem 1.4), showing that the weighted $L^1$ distance between two solutions is non-increasing in time.
Using this contraction property, combined with the strong contraction principle for rescaled functions (Lemma 4.2) and compactness arguments (Lemmas 4.3 and 4.4), they demonstrate that the rescaled solution $\tilde{u}(y, \tau) = t^\alpha u(t^\beta y, t)$ converges uniformly on compact subsets of $\mathbb{R}^n \setminus \{0\}$ and in the weighted $L^1$ space to a specific singular self-similar solution $f_{\lambda_0}$ as $\tau = \log t \to \infty$ .

Key Contributions and Results

Theorem 1.1 (Existence and Uniqueness): The paper establishes the existence and uniqueness of a singular radially symmetric solution $f$ to the self-similar profile equation satisfying specific boundary conditions at the origin ( $\lim_{|x|\to 0} |x|^{\gamma}f(x) = \eta$ ) and at infinity ( $\lim_{|x|\to \infty} |x|^{\frac{n-2}{m}}f(x) = D(\eta)$ ). This result is proven via a new fixed-point approach.
Theorem 1.2 (Asymptotic Expansion): The authors derive a precise second-order asymptotic expansion of the singular self-similar solution $f$ near the origin. They show that the rescaled profile $w(\rho)$ belongs to $C^2([0, \infty))$ and provide explicit expressions for $w(0)$ , $w_\rho(0)$ , and $w_{\rho\rho}(0)$ .
Theorem 1.4 (Weighted $L^1$ -Contraction): A contraction inequality is proven for the difference of two solutions in a weighted $L^1$ space, utilizing the constructed weight function $\phi_\mu$ . This is a crucial tool for stability analysis.
Theorem 1.5 (Asymptotic Large Time Behavior): The main result extends the asymptotic convergence theorem to the case where the initial data decay rate $\gamma$ satisfies $n \le \gamma < \frac{n-2}{m}$ . The paper proves that under mild additional assumptions on the initial data (specifically $u_0 - A_0|x|^{-\gamma} \in L^1(\phi_\mu)$ ), the rescaled solution converges to a forward singular self-similar solution determined by the constant $A_0$ .

Significance and Claims
The paper claims to provide a "new direct proof" of the existence and uniqueness of singular self-similar solutions, offering an alternative to the inversion method used in prior literature [23]. Furthermore, it extends the known asymptotic behavior results to a previously unaddressed range of the singularity exponent $\gamma$ (specifically $n \le \gamma < \frac{n-2}{m}$ ). The authors emphasize that understanding these singular self-similar solutions is essential for characterizing the behavior of general singular solutions that blow up at the origin. The work relies on and cites previous foundational studies on the fast diffusion equation, porous medium equations, and Yamabe flow (e.g., works by Daskalopoulos, del Pino, King, Sesum, Vazquez, and Hui), positioning its results as a rigorous extension of the existing theory regarding singular solutions in the subcritical fast diffusion regime.

More like this