[Paper Review] Convergence to equilibrium for positive solutions of some mutation-selection model

Jérôme Coville|arXiv (Cornell University)|Aug 29, 2013

Mathematical and Theoretical Epidemiology and Ecology Models24 references22 citations

TL;DR

This paper establishes global convergence to equilibrium for positive solutions of a nonlocal mutation-selection model with Neumann boundary conditions, using nonlinear relative entropy identities and orthogonal decomposition. It proves unique, globally stable positive steady states exist under blind competition (K independent of x), and for small perturbations of the kernel, the steady state remains positively asymptotically stable.

ABSTRACT

In this paper we are interested in the long time behaviour of the positive solutions of the mutation selection model with Neumann Boundary condition: $$ \\frac{\\partial u(x,t)}{dt}=u\\left[r(x)-\\int_{\\O}K(x,y)|u|^{p}(y)\\,dy\ ight]+\ abla\\cdot\\left(A(x)\ abla u(x)\ ight),\\qquad \ ext{in}\\quad \\R^+\ imes\\O$$ where $\\O\\subset \\R^N$ is a bounded smooth domain, $k(.,.) \\in C(\\bar \\O \ imes C(\\bar\\O), \\R), p\\ge 1$ and $A(x)$ is a smooth elliptic matrix. In a blind competition situation, i.e $K(x,y)=k(y)$, we show the existence of a unique positive steady state which is positively globally stable. That is, the positive steady state attracts all the possible trajectories initiated from any non negative initial datum. When $K$ is a general positive kernel, we also present a necessary and sufficient condition for the existence of a positive steady states. We prove also some stability result on the dynamic of the equation when the competition kernel $K$ is of the form $K(x,y)=k_0(y)+\\eps k_1(x,y)$. That is, we prove that for sufficiently small $\\eps$ there exists a unique steady state, which in addition is positively asymptotically stable. The proofs of the global stability of the steady state essentially rely on non-linear relative entropy identities and an orthogonal decomposition. These identities combined with the decomposition provide us some a priori estimates and differential inequalities essential to characterise the asymptotic behaviour of the solutions.

Motivation & Objective

• To analyze the long-time behavior of positive solutions in a nonlocal mutation-selection model without assuming small mutation rates.
• To establish the existence and global stability of a unique positive steady state under blind competition (K(x,y) = k(y)).
• To derive a necessary and sufficient condition for the existence of a positive steady state when the competition kernel K is general.
• To investigate the stability of the dynamic when the kernel is a small perturbation of a blind competition form.
• To develop and apply nonlinear relative entropy identities and orthogonal decomposition to characterize asymptotic behavior.

Proposed method

• Uses nonlinear relative entropy identities to derive a priori estimates and differential inequalities for solution decay.
• Applies an orthogonal decomposition of the solution in terms of the principal eigenfunction φ₁ to decouple dynamics.
• Employs Schauder parabolic estimates and bootstrap arguments to prove uniform boundedness and regularity of solutions.
• Utilizes the parabolic maximum principle and sub-solution techniques to ensure positivity and boundedness of solutions.
• Applies compactness arguments and diagonal extraction to pass to the limit in approximating sequences.
• Analyzes the system via energy-type estimates involving the L^p norm and the weight φ₁, leading to Lyapunov-type decay.

Experimental results

Research questions

• RQ1 Under what conditions does the mutation-selection model with Neumann boundary conditions converge globally to a unique positive steady state?
• RQ2 What is the necessary and sufficient condition for the existence of a positive steady state when the competition kernel K is general?
• RQ3 How does the stability of the steady state behave when the kernel is a small perturbation of a blind competition kernel?
• RQ4 Can global stability be established without assuming small mutation rates, contrary to typical adaptive dynamics frameworks?
• RQ5 What role does the principal eigenfunction φ₁ of the operator ∇·(A∇) + r play in characterizing the long-time behavior?

Key findings

• In the blind competition case (K(x,y) = k(y)), a unique positive steady state exists and attracts all non-negative initial data globally.
• The steady state is positively globally stable, meaning all positive solutions converge to it as t → ∞.
• For a general positive kernel K, a necessary and sufficient condition for the existence of a positive steady state is derived via spectral analysis.
• When K(x,y) = k₀(y) + εk₁(x,y) with small ε > 0, a unique positive steady state exists and is positively asymptotically stable.
• The relative entropy identity combined with orthogonal decomposition yields differential inequalities that imply exponential decay of the solution to the steady state.
• Uniform bounds on the L^p norm and C^{1,α} regularity of solutions are established independently of ε, enabling convergence via compactness.