In this paper, we propose and analyze a nonlocal cooperative reaction--diffusion system with free boundaries and drift terms, motivated by directional epidemic spread. Lacking a variational structure but requiring sharper regularity of solutions, the model poses substantial analytical challenges compared with previous works~\cite{Du,Berestycki2016a,Berestycki2016b,Cao2019,NguyenVo2022,Tang2024a,Tang2024b}. We first establish the well-posedness of the local problem and the global existence and uniqueness of classical solutions in $C^1$ space.
We then study the associated nonlocal eigenvalue problem, proving the existence, simplicity, qualitative properties, and asymptotic behavior of the principal eigenvalue. The analysis employs Fredholm theory, the Crandall--Rabinowitz bifurcation theorem, and Hadamard-type derivative formulas to describe its parameter dependence and connection with the basic reproduction number~$R_0$.
Building on this spectral characterization, we show that the system admits a \emph{sharp vanishing--spreading dichotomy} in its long-term dynamics. When $R_0\le1$, all solutions vanish; for $R_0>1$, the outcome depends on the initial domain size~$h_0$ and the free-boundary expansion rate~$μ$. There exists a critical habitat length~$\mathcal L^\ast$ such that if $h_0<\mathcal L^\ast$, a threshold $\widehatμ>0$ separates vanishing ($μ\in(0,\widehatμ]$) from spreading ($μ>\widehatμ$). In the spreading regime, solutions converge to the unique positive steady state, while in the vanishing regime they decay uniformly to zero. These results provide a rigorous framework for the threshold dynamics of cooperative--advective nonlocal systems and offer mathematical insight for further studies in epidemic modeling, ecological invasion, and population dynamics.
academicThe eigentheory for nonlocal cooperative-advective system and its role in the study of free boundary system for directional epidemic models
- Paper ID: 2510.10024
- Title: The eigentheory for nonlocal cooperative-advective system and its role in the study of free boundary system for directional epidemic models
- Authors: Soufiane Bentout (University of Ain Temouchent), Hoang-Hung Vo (Saigon University)
- Classification: math.AP (Partial Differential Equations)
- Submission Date: October 14, 2025
- Paper Link: https://arxiv.org/abs/2510.10024
This paper investigates nonlocal cooperative reaction-diffusion systems driven by directional epidemic propagation, incorporating free boundaries and drift terms. The model presents substantial analytical challenges compared to prior work due to the absence of a variational structure combined with stringent regularity requirements for solutions. The authors first establish well-posedness of the local problem and global existence and uniqueness of classical solutions in C¹ spaces. Subsequently, they investigate the associated nonlocal eigenvalue problem, proving the existence, simplicity, qualitative properties, and asymptotic behavior of the principal eigenvalue. The analysis employs Fredholm theory, the Crandall-Rabinowitz bifurcation theorem, and Hadamard-type derivative formulas to characterize parameter dependence and its connection to the basic reproduction number R₀. Based on this spectral characterization, the authors establish a sharp vanishing-spreading dichotomy in the long-term dynamics: when R₀ ≤ 1, all solutions vanish; when R₀ > 1, the outcome depends on the initial domain size h₀ and the free boundary expansion rate μ.
- Threat of waterborne diseases: Waterborne diseases such as cholera, botulism, and giardiasis continue to pose significant threats to global public health, particularly in developing countries. It is estimated that waterborne diseases account for approximately 3.6% of the global disease burden (DALYs), causing approximately 1.5 million deaths annually.
- Limitations of existing models: Classical reaction-diffusion frameworks based on the Laplacian operator are inadequate for handling inherently nonlocal or directionally-biased movement. Traditional models often neglect long-distance transmission and directional drift driven by environmental flow or human activity.
- Mathematical challenges: The combination of nonlocal diffusion systems with drift terms results in the absence of a variational structure, making classical spectral theory methods difficult to apply directly, necessitating the development of new analytical techniques.
- Practical necessity: Accurate modeling of pathogen transmission in water bodies, accounting for environmental factors such as river flow and wind direction
- Theoretical challenges: Development of new mathematical theory for handling nonlocal cooperative systems with free boundary problems
- Applied value: Provision of mathematical foundations for epidemiology, biological invasion, and population dynamics research
- Well-posedness theory: In the absence of a variational structure, the authors establish a global well-posedness theory for nonlocal cooperative reaction-diffusion systems, proving existence, uniqueness, and global existence of classical solutions in C¹ spaces.
- Spectral theory development:
- Proof of existence and simplicity of the principal eigenvalue
- Establishment of continuity and monotonicity of the principal eigenvalue with respect to parameters
- Analysis of differentiability of the principal eigenvalue with respect to diffusion coefficients using Hadamard-type formulas
- Characterization of basic reproduction number: Establishment of equivalence between the sign of the principal eigenvalue and the basic reproduction number R₀: R₀ > 1 ⟺ λₚ < 0
- Vanishing-spreading dichotomy: Complete classification of long-term dynamics:
- When R₀ ≤ 1, all solutions vanish
- When R₀ > 1 and h₀ < L*, there exists a critical threshold μ̂ distinguishing vanishing and spreading behavior
- In the spreading case, solutions converge to a unique positive steady state
Consider the following nonlocal cooperative reaction-diffusion system:
ut=d1[∫g(t)h(t)J1(x−y)u(t,y)dy−u(t,x)]+pux−a(x)u(t,x)+H(v(t,x))vt=d2[∫g(t)h(t)J2(x−y)v(t,y)dy−v(t,x)]+qvx−b(x)v(t,x)+G(u(t,x))where:
- u(t,x), v(t,x) represent the pathogen density in the environment and infected individual density, respectively
- J₁, J₂ are symmetric diffusion kernels
- p, q are drift coefficients representing directional transmission
- H(v), G(u) are cooperative interaction terms satisfying concavity conditions
- g(t), h(t) is the time-varying spatial domain
The boundary motion is determined by the following equations:
h′(t)=μ[∫g(t)h(t)∫h(t)∞J1(x−y)u(t,x)dydx+ρ∫g(t)h(t)∫h(t)∞J2(x−y)v(t,x)dydx]g′(t)=−μ[∫g(t)h(t)∫−∞g(t)J1(x−y)u(t,x)dydx+ρ∫g(t)h(t)∫−∞g(t)J2(x−y)v(t,x)dydx]- Generalization of comparison principles: Establishment of new comparison principles for nonlocal diffusion-advection operators, overcoming technical difficulties posed by drift terms.
- Spectral analysis methods:
- Utilization of operator decomposition techniques: L = J - T, where J is a compact positive operator and T contains drift terms
- Analysis of spectral radius through properties of the resolvent R(ξ) = (T + ξI)⁻¹
- Application of the Krein-Rutman theorem to handle non-self-adjoint operators
- Application of Hadamard formulas: Establishment of differentiability of the principal eigenvalue with respect to diffusion coefficients:∂λ∗/∂dj=∫−ZZw∗j(x)[∫−ZZJj(x−y)φj(y)dy−φj(x)]dx
Under assumptions (J) and (GH), for any h₀ > 0, problem (1) admits a unique positive solution (u,v,g,h) defined for all t > 0.
The linearized eigenvalue problem admits a principal eigenvalue λ*, which is simple and corresponds to a strictly positive eigenfunction.
The principal eigenvalue λ*(Z) is strictly increasing and continuous with respect to domain size Z. There exists a critical domain size Z* such that:
- λ*(Z) < 0 when Z < Z*
- λ*(Z) > 0 when Z > Z*
- Z* < ∞ if and only if R₀ > 1
Assume R₀ > 1 and h₀ < L*. Then there exists a threshold μ̂ > 0 such that:
- Vanishing occurs when μ ∈ (0, μ̂]
- Spreading occurs and converges to a positive steady state when μ > μ̂
- Local existence: Application of the Banach fixed-point theorem through construction of a contraction mapping
- Global existence: Utilization of comparison principles and a priori estimates to prove solutions do not blow up in finite time
- Regularity: Establishment of C¹ regularity through the method of characteristics and integral equation techniques
- Operator decomposition: Transformation of the problem into spectral analysis of compact operators
- Monotonicity: Utilization of operator sequences and monotonicity of spectral radius
- Continuity: Establishment through continuous dependence of operator norms
- Vanishing case: Construction of supersolutions and application of comparison principles
- Spreading case: Combination of subsolution construction and asymptotic analysis
- Critical threshold: Determination through combination of continuity and monotonicity
- Provision of mathematical models for spatial transmission of waterborne diseases (e.g., cholera)
- Quantification of environmental factors' (river flow, wind direction) impact on epidemic propagation
- Theoretical basis for public health intervention strategies
- Simulation of invasive species dispersal in heterogeneous environments
- Analysis of habitat fragmentation effects on population dynamics
- Prediction of species distribution changes under climate change
- Strong theoretical innovation:
- First systematic study of free boundary problems for nonlocal cooperative-advective systems
- Development of new techniques for handling systems without variational structure
- Establishment of a complete spectral theory framework
- Advanced mathematical techniques:
- Skillful combination of Fredholm theory, bifurcation theory, and Hadamard formulas
- Innovative handling of technical difficulties posed by drift terms
- Rigorous mathematical proofs and comprehensive theoretical framework
- Completeness of results:
- Complete analysis from well-posedness to long-term behavior
- Sharp characterization of vanishing-spreading dichotomy
- Precise description of parameter dependence
- Significant practical value:
- Direct application to epidemiological and ecological problems
- Mathematical foundation for practical prediction and control
- Restrictions of model assumptions:
- Assumption of symmetric kernels, though asymmetry may exist in practice
- Constant drift coefficients, which may vary spatially in reality
- One-dimensional spatial setting; higher-dimensional generalizations require further study
- Absence of numerical verification:
- Pure theoretical analysis lacking numerical experiments
- No practical methods for parameter estimation provided
- Lack of comparison with real data
- Limited application scope:
- Primarily focused on cooperative systems; competitive systems require different treatment
- Assumption of environmental homogeneity; generalization to heterogeneous environments is complex
- Academic value: Significant contribution to nonlocal diffusion equation theory, expected to have broad influence in mathematical analysis
- Application prospects: Provision of new tools for epidemic modeling and ecological research with important practical value
- Future research directions: Opens multiple new research avenues including higher-dimensional generalizations, stochastic perturbations, and heterogeneous environments
- Epidemiological research: Spatial transmission modeling and prediction of waterborne diseases
- Ecological applications: Invasive species dispersal, population dynamics analysis
- Theoretical mathematics: Nonlocal diffusion equations, free boundary problems research
- Numerical analysis: Provision of theoretical foundation for development of related numerical methods
This paper achieves important breakthroughs in the study of free boundary problems for nonlocal cooperative-advective systems, establishing a comprehensive mathematical theoretical framework. Major achievements include well-posedness theory, spectral analysis, and complete characterization of long-term dynamics. This work is not only theoretically significant in mathematics but also provides powerful tools for practical applications in epidemiology and ecology.
Future research directions include:
- Generalization to higher-dimensional spaces
- Treatment of asymmetric kernels
- Effects of stochastic perturbations
- Development of numerical methods
- Validation with real data
This work lays a solid theoretical foundation for further development in related fields and is expected to have sustained and profound impact.