Asymptotic behavior of solutions to elliptic problems with Robin boundary conditions

1. Problem Statement

The paper investigates the asymptotic behavior of positive solutions $u_\beta$ to the semilinear elliptic equation with Robin boundary conditions as the parameter $\beta \to 0^+$. The problem is formulated as:

$$\begin{cases} -\Delta u = u^p & \text{in } \Omega, \\ u > 0 & \text{in } \Omega, \\ \frac{\partial u}{\partial \nu} + \beta u = 0 & \text{on } \partial\Omega, \end{cases}$$

where:

• $\Omega \subset \mathbb{R}^N$ ($N \ge 2$) is a bounded smooth domain.
• $\nu$ is the outward unit normal.
• $p \ge 0$ is the exponent of the nonlinearity.
• $\beta > 0$ is the Robin parameter.

Context:

• $\beta \to \infty$ corresponds to the Dirichlet problem ($u=0$ on $\partial\Omega$).
• $\beta \to 0$ corresponds to the Neumann problem ($\frac{\partial u}{\partial \nu} = 0$ on $\partial\Omega$).
• The case $p=1$ is treated as an eigenvalue problem: $-\Delta u = \lambda_{1,\beta} u$.

The primary goal is to characterize the precise asymptotic profile of $u_\beta$ as $\beta \to 0$ across different regimes of $p$ (sublinear, linear, superlinear, critical, and supercritical).

2. Methodology

The authors employ a combination of variational methods, maximum principles, elliptic regularity theory, and a specific auxiliary function transformation to analyze the solutions.

Key Technical Tools:

1. Auxiliary Function Transformation:
   The core of the analysis involves defining an auxiliary function $v_\beta$ to separate the "bulk" behavior from the boundary effects.

   • Let $c_\beta = \int_\Omega u_\beta^p \, dx$.
   • Define a constant $d_\beta = \frac{c_\beta}{\beta |\partial\Omega|}$.
   • Define $v_\beta = u_\beta - d_\beta$.
   • This transforms the original problem into a new equation for $v_\beta$ with a zero-mean boundary condition ($\int_{\partial\Omega} v_\beta \, dS = 0$).
   • The authors prove that $v_\beta$ remains bounded relative to $d_\beta$ (or converges to 0 after normalization), allowing them to approximate $u_\beta \approx d_\beta$.

2. Uniqueness Analysis:
   The paper establishes the uniqueness of positive solutions for specific ranges of $p$ and small $\beta$. This is done by assuming two distinct solutions $u$ and $v$, analyzing the difference $w = u-v$, and using the mean value theorem to derive a contradiction via the sign of the resulting linearized operator.

3. Existence in Supercritical Regime:
   For the supercritical case ($p > \frac{N+2}{N-2}$), standard variational methods fail due to the lack of compactness in Sobolev embeddings.

   • The authors restrict the domain to a ball ($\Omega = B_1(0)$).
   • They utilize known radial solutions to the problem in the whole space $\mathbb{R}^N$ (referencing Serrin and Zou).
   • A scaling argument is used to match the boundary condition of the ball problem, proving existence for small $\beta$.

4. Asymptotic Estimates:
   By analyzing the integral identities derived from the transformed equations and applying elliptic regularity estimates (e.g., $W^{2,p}$ bounds), the authors derive the leading-order terms of $d_\beta$ (and thus $u_\beta$) as $\beta \to 0$.

3. Key Contributions and Results

The paper provides a unified classification of the asymptotic behavior based on the exponent $p$:

A. Uniqueness (Theorem 1.1)

• Subcritical/Linear ($0 \le p \le 2^*-1, p \neq 1$): Unique positive solution for sufficiently small $\beta$.
• Linear ($p=1$): Unique positive eigenfunction for any $\beta$.
• Supercritical ($p > 2^*-1$): At most one positive solution for small $\beta$.

B. Asymptotic Behavior by Regime

1. Sublinear Case ($0 \le p < 1$): Blow-up

• Behavior: The solution $u_\beta$ blows up uniformly in $\Omega$ as $\beta \to 0$.
• Asymptotic Formula:
  $$u_\beta \sim \left( \frac{|\Omega|}{|\partial\Omega|} \right)^{\frac{1}{1-p}} \beta^{-\frac{1}{1-p}}$$
• Significance: The solution grows to infinity, and the rate depends explicitly on the ratio of the domain volume to the boundary surface area.

2. Linear Case ($p = 1$): Convergence to Constant

• Behavior: The solution (normalized eigenfunction) converges uniformly to the constant $1$.
• Asymptotic Formula:
  $$u_\beta \to 1 \quad \text{as } \beta \to 0$$
• Note: The associated eigenvalue $\lambda_{1,\beta}$ behaves like $\beta \frac{|\partial\Omega|}{|\Omega|}$.

3. Superlinear Case ($p > 1$): Convergence to Zero

• Behavior: The solution $u_\beta$ converges uniformly to zero as $\beta \to 0$.
• Asymptotic Formula:
  $$u_\beta \sim \left( \frac{|\partial\Omega|}{|\Omega|} \right)^{\frac{1}{p-1}} \beta^{\frac{1}{p-1}}$$
• Significance: Unlike the sublinear case, the solution vanishes. The decay rate is determined by the same geometric ratio but with a positive power of $\beta$.

C. Existence in Supercritical Regime (Theorem 1.4)

• For $p > \frac{N+2}{N-2}$, existence is generally unknown for arbitrary domains.
• Result: If $\Omega$ is a ball and $0 < \beta < \frac{2}{p-1}$, a radial positive solution exists.
• Method: Constructed by scaling a known radial solution from $\mathbb{R}^N$ to satisfy the Robin boundary condition on the ball.

4. Significance and Impact

1. Unified Framework: The paper successfully unifies the analysis of sublinear, linear, and superlinear regimes under the Robin boundary condition, revealing a sharp transition in behavior at $p=1$.
2. Geometric Dependence: It explicitly quantifies how the geometry of the domain (specifically the ratio $|\Omega|/|\partial\Omega|$) dictates the magnitude of the solution in the limit $\beta \to 0$.
3. Resolution of Open Problems:
   • It provides the exact leading-order asymptotics for the sublinear case, refining previous estimates.
   • It establishes existence results for the supercritical Robin problem on balls, a regime where standard variational techniques usually fail.
4. Methodological Innovation: The use of the auxiliary function $v_\beta = u_\beta - d_\beta$ effectively isolates the singular behavior of the solution, allowing for precise asymptotic expansions that would be difficult to obtain through direct energy estimates alone.

In summary, this work fills a gap in the literature regarding the Neumann limit ($\beta \to 0$) of Robin problems, providing rigorous proofs for existence, uniqueness, and precise asymptotic profiles across all non-negative exponents $p$.