Fully nonlinear logistic equations with sanctuary

1. Problem Statement

The paper investigates the existence, uniqueness, and asymptotic behavior of positive solutions to a stationary fully nonlinear logistic equation in a bounded domain $\Omega \subset \mathbb{R}^N$. The problem is formulated as:

$$\begin{cases} F(x, D^2u) + \mu u = k(x)u^p & \text{in } \Omega \\ u > 0 & \text{in } \Omega \\ u = 0 & \text{on } \partial\Omega \end{cases}$$

Key Components:

• Operator $F$: A uniformly elliptic, positively homogeneous (degree 1) operator acting on the Hessian $D^2u$. This generalizes the classical Laplacian $\Delta u$ to nonlinear diffusion scenarios.
• Parameter $\mu$: Represents the natural growth rate of the population.
• Nonlinearity $p$: $p > 1$, representing the logistic competition term.
• Reaction Term $k(x)$: A bounded continuous function representing environmental carrying capacity. The paper analyzes two distinct cases:
  1. Case (K1): $k(x) \geq k_0 > 0$ everywhere in $\Omega$ (standard logistic growth).
  2. Case (K2): $k(x) = 0$ in a subdomain $\Omega_0 \subset \subset \Omega$ (the "sanctuary"), while $k(x) > 0$ outside $\Omega_0$. This models a region where the population faces no competition or predation.
• Solution Concept: Solutions are understood in the viscosity sense, which is standard for fully nonlinear elliptic equations.

2. Methodology

The authors employ a combination of viscosity solution theory, spectral theory for nonlinear operators, and barrier methods.

• Principal Eigenvalue ($\lambda_1^+$): The central tool is the principal eigenvalue defined via the maximum principle (following Berestycki, Nirenberg, and Varadhan):
  $$\lambda_1^+(\Omega) := \sup \{ \mu \in \mathbb{R} : \exists \psi > 0 \text{ s.t. } F(x, D^2\psi) + \mu\psi \leq 0 \text{ in } \Omega \}$$
  This value acts as the critical threshold separating existence from non-existence.

• Comparison Principles: The authors establish and utilize strong comparison principles for viscosity sub- and supersolutions. A key technical achievement is proving that if $u_1 \geq u_2$ on the boundary and satisfy the respective inequalities, then $u_1 \geq u_2$ in the domain. They also prove a strong comparison principle: if $u_1 \geq u_2$ and they touch at an interior point, they are identical.

• Barrier Construction:

  • Subsolutions: Constructed using the principal eigenfunction $\phi_\Omega$ scaled by a small constant.
  • Supersolutions:
    • For Case (K1), a constant function suffices.
    • For Case (K2), the construction is non-trivial. The authors construct a composite supersolution that behaves like the eigenfunction near the sanctuary boundary and transitions to a specific barrier function near the outer boundary to handle the singularity at $\partial \Omega_0$.

• Monotone Iteration: Existence is proven by constructing monotone sequences of sub- and supersolutions that converge to the unique solution.

• Asymptotic Analysis: The behavior of solutions as $\mu$ approaches the boundaries of the existence interval is studied using compactness arguments, $C^{0,\alpha}$ estimates, and blow-up analysis.

3. Key Contributions and Results

A. Sharp Existence and Uniqueness Conditions (Theorem 1)

The paper provides precise conditions for the existence of a unique positive solution based on the relationship between $\mu$ and the principal eigenvalues of the domain and the sanctuary.

1. Case (K1) (No Sanctuary): A unique positive solution exists if and only if $\mu > \lambda_1^+(\Omega)$.
2. Case (K2) (With Sanctuary): A unique positive solution exists if and only if
   $$\lambda_1^+(\Omega) < \mu < \lambda_1^+(\Omega_0)$$
   • If $\mu \leq \lambda_1^+(\Omega)$, the growth rate is too low to sustain a population against the logistic term.
   • If $\mu \geq \lambda_1^+(\Omega_0)$, the growth rate is too high for the "sanctuary" (where $k=0$) to support a bounded solution; the population would blow up.

B. Asymptotic Behavior as $\mu \to \lambda_1^+(\Omega)^+$ (Theorem 2)

As the growth rate $\mu$ approaches the lower threshold from above:

• The solution $u_\mu$ converges uniformly to 0 in $\Omega$.
• The normalized solution $v_\mu = u_\mu / \|u_\mu\|_\infty$ converges uniformly to the principal eigenfunction $\phi_\Omega$ of the operator $F$ in $\Omega$.

C. Asymptotic Behavior as $\mu \to \lambda_1^+(\Omega_0)^-$ (Theorem 3)

As the growth rate $\mu$ approaches the upper threshold (the eigenvalue of the sanctuary) from below:

• Inside the Sanctuary ($\Omega_0$): The solution $u_\mu$ blows up (converges to $+\infty$) uniformly. The normalized solution converges to the principal eigenfunction $\phi_{\Omega_0}$.
• Outside the Sanctuary ($\Omega \setminus \Omega_0$): The solution converges locally uniformly to a function $U$. This $U$ is the minimal positive solution to the "blow-up" problem:
  $$\begin{cases} F(x, D^2u) + \lambda_1^+(\Omega_0)u = k(x)u^p & \text{in } \Omega \setminus \Omega_0 \\ u = 0 & \text{on } \partial\Omega \\ u = +\infty & \text{on } \partial\Omega_0 \end{cases}$$
  This result characterizes the "large solution" behavior where the population density becomes infinite at the boundary of the sanctuary.

4. Significance and Novelty

• Generalization to Fully Nonlinear Operators: While logistic equations with the Laplacian ($\Delta u$) are well-studied (e.g., by Ouyang, Du, Huang), this paper extends these results to fully nonlinear, non-variational operators. This is significant because many physical diffusion processes are nonlinear and do not admit a variational structure (Rayleigh quotient).
• New Proofs for Classical Cases: The authors emphasize that their approach, based on comparison principles and explicit barrier constructions, provides new proofs even for the classical linear case ($F=\Delta$). This offers a more unified and robust framework for understanding the role of the principal eigenvalue.
• Sanctuary Dynamics: The paper rigorously analyzes the "sanctuary" effect in a nonlinear setting. It demonstrates how the eigenvalue of the sanctuary ($\lambda_1^+(\Omega_0)$) acts as a hard upper bound for the growth parameter $\mu$ to maintain a bounded population distribution.
• Viscosity Solution Framework: The work contributes to the theory of viscosity solutions by establishing comparison principles and asymptotic behaviors for logistic-type equations with singular boundary conditions (infinite boundary values).

In summary, the paper establishes a complete existence theory for fully nonlinear logistic equations with a sanctuary, identifying the principal eigenvalues as the critical thresholds for existence and characterizing the precise blow-up behavior of solutions as these thresholds are approached.