Mono-Stable Reaction-Diffusion-Advection Equation

• Mono-stable reaction-diffusion-advection equations are parabolic PDEs featuring a unique globally attracting equilibrium, modeling processes like invasion and fixation.
• They combine local reaction kinetics, spatial diffusion, and directional advection to capture phenomena in ecology, combustion, pattern formation, and population genetics.
• Analytical approaches such as travelling wave analysis and geometric desingularisation reveal threshold speeds, front classifications, and regime behaviors.

A mono-stable reaction-diffusion-advection equation refers to a class of parabolic partial differential equations (PDEs) in which a density field (e.g., population, concentration, mass fraction) evolves under three fundamental processes: local reaction (mono-stable kinetics), spatial diffusion, and directional advection or convection. Mono-stability in the reaction term implies that the system possesses only one globally attracting nontrivial equilibrium, typically representing successful invasion or fixation by a species, chemical, or phase. These equations underpin a broad spectrum of phenomena in ecology, combustion, pattern formation, and population genetics.

1. Mathematical Formulation and Mono-Stability Criteria

The generic form for one-dimensional mono-stable reaction-diffusion-advection equations is:

$$g(u)\,u_\tau + f(u)\,u_x = [D(u)u_x]_x + \rho(u), \quad u\in[0,1]$$

Here, $u=u(\tau,x)$ denotes the density at time $\tau$ and position $x$; $g(u)$ is an accumulation term; $D(u)$ is a possible density-dependent diffusivity; $f(u)$ encodes advection/convection; and $\rho(u)$ is the reaction term.

Mono-stable reaction: The reaction function $\rho$ satisfies:

$$\rho(0) = \rho(1) = 0, \qquad \rho(u) > 0\;\;\forall\,u \in (0,1)$$

This criterion ensures that $u=1$ is globally stable, while $u=0$ is unstable (Cantarini et al., 14 Feb 2025). Classical examples include the Fisher-KPP kinetics $f(u) = u(1-u)$, widely employed to model population invasion and gene propagation.

2. Existence and Uniqueness of Travelling Wave Solutions

Travelling fronts, or monotone wave solutions, are central in analyzing invasion and saturation dynamics:

• Seek solutions of the form $u(\tau, x) = U(z)$, $z = x - c\tau$, leading to second-order ODEs for $U(z)$:

$(D(U)U')' + (c\,g(U) - f(U)) U' + \rho(U) = 0$

with boundary conditions $U(-\infty)=1$, $U(+\infty)=0$. Existence hinges on the sign structure and regularity of $g, D, f, \rho$ (Cantarini et al., 14 Feb 2025, Marcelli, 27 Dec 2025).

Wave existence is typically characterized by a threshold speed $c^*$: there exists a (unique, up to translation) monotone front if and only if $c \geq c^*$. In degenerate or singular cases, sharp fronts may arise, where the solution attains an equilibrium in finite time (Marcelli, 27 Dec 2025).

The presence of cut-offs in the reaction or advection terms (e.g., Heaviside function $H(u-\varepsilon)$) can modify uniqueness and produce regimes with logarithmic or algebraic corrections to front speed (Popovic et al., 2024).

3. Determination of Minimal Wave Speed and Regime Classification

Wave speed selection is governed by linear or nonlinear eigenvalue problems derived from the leading edge analysis. For classical FKPP-type systems (normalized diffusivity and linear accumulation):

$$\frac{\partial u}{\partial t} + k u \frac{\partial u}{\partial x} H(u-\varepsilon) = \frac{\partial^2 u}{\partial x^2} + u(1-u) H(u-\varepsilon)$$

The minimal speed $c_0$ is:

$$c_0 = \begin{cases} 2,& k \leq 2\quad(\text{pulled regime}) \ k/2 + 2/k,& k > 2\quad(\text{pushed regime}) \end{cases}$$

Cut-offs $\varepsilon$ modify this to $c(\varepsilon) = c_0 - \Delta c(\varepsilon)$ with $\Delta c(\varepsilon)$ algebraically small for pushed fronts and logarithmically small for pulled fronts (Popovic et al., 2024). In systems with heterogeneous reaction domains or forced advection, multiregime behavior arises:

• Spreading for $0 \leq c < c_*$
• Vanishing for $c \geq c^*$
• Intermediate thresholds permitting either outcome depending on initial data (Bouhours et al., 2020, Zhao et al., 2015)

4. Analytical Techniques: Singular Reductions and Geometric Desingularisation

The rigorous analysis of wave existence and profiles often employs:

• First-order singular BVPs: By exploiting monotonicity, the second-order ODE can be reformulated for $Z(u)=D(u)U'(z(u))$, yielding singular first-order differential equations with boundary conditions at equilibria. Existence and sharpness are classified via local analysis at $u=0,1$, often requiring integral tests on $\rho$ and $D$ (Cantarini et al., 14 Feb 2025, Marcelli, 27 Dec 2025).
• Geometric desingularisation (“blow-up”): For equations with non-smooth cut-offs (e.g., switching surfaces at $u=\varepsilon$), analysis proceeds by lifting to a higher-dimensional dynamical system (charts $K_1$, $K_2$) and matching invariant manifolds. Persistence and uniqueness of orbits follow from transversality arguments (Popovic et al., 2024).

These methods accommodate degenerate diffusion (e.g., $D(u)$ vanishing at endpoints), sign-changing coefficients, and general nonlinearities.

5. Front Classification: Classical vs. Sharp, Finite vs. Infinite Propagation

Wavefronts are categorized based on their approach to equilibria $u=0,1$:

• Classical fronts: $U'(z)\to 0$ as $z\to\pm\infty$; the profile decays/explodes exponentially.
• Sharp fronts: $U$ reaches an equilibrium at a finite $z$, possible for degenerate $D(u)$ or abrupt changes in coefficients (Cantarini et al., 14 Feb 2025, Marcelli, 27 Dec 2025).

Local criteria such as $\lim_{u\to 0^+}(d(u)/u)^{1/(p-1)}\rho(u)$ and corresponding integrals of $u/\rho(u)$ determine whether finite propagation occurs (front reaches zero at a finite position) (Marcelli, 27 Dec 2025).

Table: Front types and conditions

Front Type	Criterion at $u=0$ or $u=1$	Example Models
Classical	$D(u)>0$ near equilibrium	Standard FKPP, $D>0$
Sharp (finite time)	$D(u)<0$ or $D(0)=0$	Degenerate $D(u)$, $c^*>0$

6. Impact of Advection and Variants

Advection modifies both the minimal wave speed and the qualitative behavior of solutions. In models with explicit advection ($\beta u_x$), spreading-vanishing thresholds depend on the relation of $|\beta|$ to the intrinsic speed $c_0$: strong advection may induce virtual spreading or vanishing, where the advancing front decouples from long-term invasion (Zhao et al., 2015). Competing advective mechanisms may yield divergence in front width or pinning of speed regardless of diffusion or coupling (Kogan et al., 2015).

Environmental heterogeneity and moving reaction domains introduce forced speed constraints, producing regimes where initial conditions and shifting rates jointly determine long-term outcomes (Bouhours et al., 2020).

7. Extensions, Open Problems, and Applications

The general framework of mono-stable reaction-diffusion-advection equations supports extensive extensions:

• Nonlinear and density-dependent diffusion, including $p$-Laplacian operators (Marcelli, 27 Dec 2025)
• Mixed and free boundary problems modeling expanding habitats (Zhao et al., 2015)
• Systems with cut-offs and switching kinetics (Popovic et al., 2024)
• Multi-dimensional settings, including aggregative movement (accumulation terms with sign changes) (Cantarini et al., 14 Feb 2025)

Open directions include multi-dimensional fronts, time-varying environments, weak regularity for sharp fronts, and refined eigenvalue-based characterization of threshold speeds.

A plausible implication is that the interaction of mono-stable kinetics, heterogeneous advection, and diffusive degeneracy governs rich front dynamics that cannot be encapsulated by classical linear or homogeneous models. This suggests the necessity of singular, geometric, and regime-based techniques for analyzing long-term behavior in applied and theoretical settings.