🔢 mathematics

Spectral Theory of Fractional Cooperative Systems and Threshold Dynamics in Epidemic Models

This paper establishes a spectral theory for cooperative systems involving the spectral fractional Laplacian to derive a sharp principal eigenvalue criterion and variational characterizations, which are then applied to prove the existence, uniqueness, and threshold-type long-time dynamics of solutions in fractional epidemic models.

Original authors: Cong-Bang Trang, Hoang-Hung Vo

Published 2026-02-17

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Cong-Bang Trang, Hoang-Hung Vo

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 a bustling city where two types of residents are interacting: Bacteria (let's call them "The Invaders") and Infected Humans (let's call them "The Carriers"). The city is our domain, $\Omega$ .

In the old days, scientists modeled how these groups moved and spread using a simple rule: Diffusion. Think of this like a crowd of people walking randomly down a street. They bump into neighbors, move slowly, and spread out evenly. This is "local" movement.

But in the real world, things are messier. A person might take a bus across town, or a bird might fly hundreds of miles in one go. This is non-local movement. In mathematics, we model this with something called the Fractional Laplacian. Instead of just walking to the next house, an organism can "jump" to a distant house instantly.

This paper is a deep dive into understanding how these "jumping" movements affect the spread of an epidemic, specifically looking at the mathematical heartbeat of the system.

Here is the breakdown of their work, translated into everyday language:

1. The "Heartbeat" of the System: The Principal Eigenvalue

The authors are obsessed with finding a single number, which they call $\lambda_p$ (the principal eigenvalue).

The Analogy: Imagine the epidemic system is a giant, complex machine. The $\lambda_p$ is the machine's "tuning fork."
What it tells us:
- If $\lambda_p$ is negative, the machine is unstable in a "good" way for stopping the disease. It means the infection will die out.
- If $\lambda_p$ is positive, the machine is unstable in a "bad" way. It means the disease will take hold and become permanent (endemic).
The Challenge: Because the movement involves "jumps" (fractional diffusion) rather than just walking, the usual math tools don't work. The authors had to invent new ways to find this number and prove it exists, is unique, and behaves predictably.

2. The "Jump" vs. The "Walk" (Fractional Orders)

The paper studies how the "jumpiness" of the movement changes the outcome. They use a dial called $s$ (the fractional order).

$s$ close to 1: The movement is mostly like walking (local).
$s$ close to 0: The movement is wild, with massive jumps (non-local).

The Discovery: The authors found that if you turn up the "jumpiness" (change $s$ ), the "heartbeat" ( $\lambda_p$ ) changes smoothly. They proved that even with these wild jumps, the system still has a clear, predictable rhythm. They also looked at what happens if the "diffusion rate" (how fast they move) goes to zero or infinity, showing that the system settles into predictable patterns in both extremes.

3. The "Basic Reproduction Number" ( $R_0$ )

You've probably heard of $R_0$ in the news. It's the average number of people one sick person infects.

The Connection: The authors proved a direct link between their "heartbeat" ( $\lambda_p$ $λ_{p}$ ) and $R_0$ $R_{0}$ .
- If $\lambda_p < 0$ , then $R_0 > 1$ . The disease spreads.
- If $\lambda_p > 0$ , then $R_0 < 1$ . The disease dies out.
Why it matters: This gives scientists a powerful new tool. Instead of simulating the whole epidemic, they can just calculate this one number ( $\lambda_p$ ) to know if the disease will persist or vanish.

4. The "Sliding Method" (Finding the Solution)

One of the hardest parts of the paper is proving that if the disease does take hold, there is exactly one stable state where it lives forever (an "endemic equilibrium").

The Analogy: Imagine trying to fit a square peg in a round hole, but the hole is shifting. The authors used a technique called the "Sliding Method."
How it works: They imagine sliding a "ceiling" (a maximum possible infection level) down and a "floor" (a minimum infection level) up. They proved that no matter how you start the epidemic, the solution gets "squeezed" between the floor and the ceiling until it settles into a single, unique spot. It's like a ball rolling down a hill until it finds the only valley at the bottom.

5. The Big Picture: Why This Matters

Before this paper, modeling epidemics with "jumping" movements (like long-distance travel or animal migration) was very difficult because the math was messy and unpredictable.

The Breakthrough: This paper provides the rulebook. It tells us:
1. The system always has a clear "heartbeat" ( $\lambda_p$ ).
2. We can predict exactly when the disease will die out or become permanent based on how fast things move and how "jumpy" they are.
3. Even with complex, long-range jumps, the epidemic eventually settles into a stable pattern, and that pattern is unique.

In Summary:
Think of this paper as building a weather forecast model for epidemics, but instead of wind and rain, it tracks "jumps" and "diffusion." The authors proved that even in a chaotic world where organisms jump long distances, the spread of disease follows a strict, predictable mathematical rhythm. If you know the rhythm (the eigenvalue), you know the future of the epidemic.

1. Problem Statement

The paper addresses two interconnected mathematical challenges:

Spectral Theory: The development of a rigorous spectral theory for a cooperative linear system involving spectral fractional Laplacians on bounded domains. The system is given by:
$\begin{cases} (-d_1\Delta_B)^{s_1}u + a_{11}(x)u + a_{12}(x)v = \lambda u, \\ (-d_2\Delta_B)^{s_2}v + a_{21}(x)u + a_{22}(x)v = \lambda v, \end{cases}$
where $B \in \{D, N\}$ denotes Dirichlet or Neumann boundary conditions, $d_i > 0$ are diffusion coefficients, $0 < s_i < 1$ are fractional orders, and the matrix $A=(a_{ij})$ satisfies cooperative conditions ( $a_{12}, a_{21} < 0$ for the operator $-A$ ).
Epidemic Modeling: The application of this spectral theory to analyze the long-time dynamics of a fractional reaction-diffusion epidemic model describing the interaction between a pathogen ( $u$ ) and an infected human population ( $v$ ):
$\begin{cases} u_t + (-d_1\Delta_N)^{s_1}u = -a(x)u + H(v), \\ v_t + (-d_2\Delta_N)^{s_2}v = -b(x)v + G(u), \end{cases}$
where $H$ and $G$ are nonlinear infection terms. The goal is to establish well-posedness, determine threshold conditions for disease persistence vs. extinction, and characterize the stability of steady states.

Key Difficulty: Unlike classical local operators, the spectral fractional Laplacian is non-local. This absence of a strong maximum principle in the classical sense and the specific regularity properties of fractional Sobolev spaces ( $H^s$ ) make standard variational and comparison techniques insufficient.

2. Methodology

The authors employ a combination of functional analysis, spectral theory, and dynamical systems techniques:

Spectral Fractional Laplacian Definition: They utilize the definition by Caffarelli and Stinga, where the operator is defined via the eigen-expansion of the underlying Laplacian ( $-\Delta_B$ ) on the domain $\Omega$ .
Krein-Rutman Theorem & Fredholm Theory: To prove the existence and simplicity of the principal eigenvalue ( $\lambda_p$ ), the authors construct a compact, positive operator on the space of continuous functions $C(\Omega)$ . They apply the Krein-Rutman theorem to establish the existence of a positive eigenfunction and the simplicity of the principal eigenvalue.
Variational Characterization: Under the symmetry condition $a_{12} = a_{21}$ , the problem is self-adjoint. The authors use the Lax-Milgram theorem and the Fredholm alternative to derive a variational characterization of $\lambda_p$ as the infimum of a Rayleigh quotient.
Asymptotic Analysis: They investigate the behavior of $\lambda_p$ $λ_{p}$ with respect to:
- Diffusion coefficients ( $d_i \to 0$ and $d_i \to \infty$ ).
- Fractional orders ( $s_i \to 0^+$ and $s_i \to 1^-$ ).
- Domain scaling ( $\Omega_l = l\Omega$ ).
- Techniques include holomorphic perturbation theory (Kato), Fatou's lemma, and comparison principles for sub/super-solutions.
Dynamical Systems Approach: For the epidemic model, they use:
- Semigroup Theory: To establish local existence and uniqueness.
- Comparison Principle: Proven via monotone dynamical systems theory to establish global existence and bounds.
- Sliding Method: Used to prove the uniqueness of the positive steady state.
- Basic Reproduction Number ( $R_0$ ): Defined via the spectral radius of the next-generation operator, linking $R_0$ to the sign of the principal eigenvalue $\lambda_p$ .

3. Key Contributions and Results

A. Spectral Theory of the Cooperative System

Existence and Simplicity: Theorem 1 establishes that the system admits a simple principal eigenvalue $\lambda_p$ with a strictly positive eigenfunction $\phi_1$ (under conditions A1-A3).
Variational Characterization: If $a_{12} = a_{21}$ , $\lambda_p$ is characterized variationally as the minimum of the energy functional $J(u)$ .
Maximum Principles: The authors establish weak and strong maximum principles. Specifically, $\lambda_p > 0$ is equivalent to the validity of the weak maximum principle for the operator.
Parameter Dependence:
- Diffusion ( $d$ ): $\lambda_p(d)$ $λ_{p} (d)$ is increasing and concave with respect to diffusion rates.
  - As $d \to 0$ , $\lambda_p \to \min_{x \in \Omega} \lambda(A(x))$ .
  - As $d \to \infty$ (Neumann), $\lambda_p \to \lambda(\bar{A})$ (the principal eigenvalue of the averaged matrix).
- Fractional Order ( $s$ ): $\lambda_p(s)$ $λ_{p} (s)$ is analytic.
  - As $s \to 1^-$ , $\lambda_p$ converges to the principal eigenvalue of the classical local system.
  - As $s \to 0^+$ , $\lambda_p$ converges to the principal eigenvalue of a matrix involving the projection operator $P_0$ .
- Domain Scaling: For Dirichlet conditions, as the domain expands ( $l \to \infty$ ), $\lambda_p \to \lambda(A)$ ; as it shrinks ( $l \to 0$ ), $\lambda_p \to \infty$ .

B. Epidemic Model Dynamics

Well-Posedness: Theorem 5 proves the global existence, uniqueness, and boundedness of solutions for the fractional epidemic system for non-negative initial data.
Threshold Dynamics: Theorem 17 establishes a sharp threshold based on the basic reproduction number $R_0$ $R_{0}$ :
- $R_0 > 1 \iff \lambda_p < 0$ .
- $R_0 < 1 \iff \lambda_p > 0$ .
Steady State Existence (Theorem 4):
- If $R_0 > 1$ , there exists a unique positive endemic equilibrium.
- If $R_0 \le 1$ , no non-trivial non-negative steady state exists (only the disease-free equilibrium).
Long-Time Behavior (Theorem 5):
- Persistence: If $R_0 > 1$ , solutions converge globally to the unique endemic equilibrium.
- Extinction: If $R_0 < 1$ , solutions converge globally to the disease-free equilibrium (0,0).

4. Significance

Theoretical Advancement: This work fills a gap in the spectral theory of fractional cooperative systems. While spectral theory for single fractional equations is well-developed, the extension to coupled systems with general boundary conditions and the rigorous derivation of the principal eigenvalue's properties (simplicity, variational forms, parameter dependence) is novel.
Bridging Local and Non-Local: The paper provides a unified framework connecting classical local epidemic models (where $s=1$ ) with non-local fractional models, showing how fractional diffusion alters persistence thresholds and asymptotic behaviors.
Application to Real-World Phenomena: By modeling anomalous diffusion (Lévy flights) in epidemic spread, the results offer a more accurate mathematical description of disease transmission in heterogeneous environments compared to classical Brownian motion models.
Methodological Rigor: The paper successfully adapts classical tools (Krein-Rutman, sliding method, comparison principles) to the non-local setting, providing a blueprint for future studies on fractional reaction-diffusion systems in biology and physics.

In summary, the paper establishes a robust spectral foundation for fractional cooperative systems and demonstrates its critical utility in determining the fate of infectious diseases, proving that the principal eigenvalue serves as the definitive threshold for disease persistence in non-local diffusion environments.