Comparison methods for semilinear elliptic problems on Riemannian manifolds with a Ricci lower bound

1. Problem Statement

The paper addresses the semilinear Dirichlet problem on a bounded domain $\Omega$ within an $n$-dimensional Riemannian manifold $(M^n, g)$:
$$\begin{cases} \Delta u + f(u) = 0 & \text{in } \Omega, \\ u > 0 & \text{in } \Omega, \\ u = 0 & \text{on } \partial\Omega, \end{cases}$$
where $\Delta$ is the Laplace–Beltrami operator and $f$ is a locally Lipschitz nonlinearity.

The study is conducted under a Ricci curvature lower bound:
$$\text{Ric}_g \geq (n-1)k g$$
for some constant $k \in \mathbb{R}$. The primary goal is to establish sharp pointwise gradient comparisons between a solution $u$ on a general manifold and a "model" solution defined on a specific warped product space. The authors aim to derive geometric consequences, including isoperimetric inequalities and estimates for the location of the maximum set ("hot spots"), and to construct non-radial solutions on spheres.

2. Methodology

The authors extend the comparison framework previously developed for Euclidean and spherical settings to general Riemannian manifolds with Ricci bounds. The methodology relies on three main pillars:

A. Model Manifolds and Radial Profiles

The authors introduce Model Manifolds $(M^n_k(F), g_k)$, which are warped products of the form:
$$(I_k \times F, dr^2 + s_k(r)^2 g_F)$$
where $F$ is a compact $(n-1)$-dimensional manifold, and $s_k(r)$ is a warping function determined by the curvature bound $k$ (involving $\sinh$, $r$, or $\sin$).

• Reduction to ODE: By assuming radial symmetry (invariance under the fiber $F$), the PDE reduces to an Ordinary Differential Equation (ODE) for the profile $U(r)$:
  $$U''(r) + (n-1)\cot_k(r)U'(r) + f(U(r)) = 0$$
  where $\cot_k(r) = s_k'(r)/s_k(r)$.

B. Admissibility and $\tau$-Functions

To ensure the existence of suitable comparison profiles, the authors define $k$-admissible nonlinearities. A function $f$ is $k$-admissible if the associated ODE solutions with Cauchy data $(R, M)$ (where $R$ is the core radius and $M$ the maximum) are positive between their zeros and strictly monotone away from the maximum.

• $\tau$-functions: The authors define normalized boundary response functions $\tau^\pm_{M,k,f}(R)$, which relate the squared gradient of the model solution at the boundary to that of a centered solution. These functions are used to match the boundary gradient of the actual solution $u$ to a specific model solution.

C. The Comparison Scheme

1. Selection: For a connected component $U$ of $\Omega \setminus \text{Max}(u)$, the authors identify a boundary component $\Gamma$ with the maximal gradient. They use the normalized gradient value $\tau(\Gamma)$ to select a corresponding model pair $(\bar{U}, \bar{u})$ from the family of radial solutions.
2. Pseudo-radial Map: A map $\Psi: U \to \mathbb{R}$ is constructed such that $u = \bar{U} \circ \Psi$. This allows the comparison of geometric quantities on level sets of $u$ with those of the radial model.
3. Maximum Principle: The core of the proof involves constructing a specific auxiliary function $F_\beta$ involving the difference of squared gradients ($W = |\nabla u|^2$ and $\bar{W} = |\nabla \bar{u}|^2 \circ \Psi$). By deriving a specific elliptic inequality for $F_\beta$ (using the Bochner formula and the Ricci lower bound), they apply the maximum principle to prove $W \leq \bar{W}$.

3. Key Contributions and Results

Theorem A: Sharp Pointwise Gradient Comparison

The central result establishes that for any point $p$ in a component $U$:
$$|\nabla u(p)|^2 \leq |\nabla \bar{u}(\Psi(p))|^2$$
Rigidity: If equality holds at a single point, then:

1. The Ricci curvature satisfies $\text{Ric}(\nabla u, \nabla u) = (n-1)k |\nabla u|^2$ on $U$.
2. The metric on $U$ splits as a warped product isometric to the model.
3. The solution $(U, u)$ is congruent to the model pair $(\bar{U}, \bar{u})$.

Theorem B: Isoperimetric-Type Inequality

By combining the gradient bound with the coarea formula, the authors derive an explicit lower bound for the ratio of the volume of the domain to the area of the maximum set (or top level set):
$$\frac{\mathcal{H}^n(U)}{\mathcal{H}^{n-1}(\Gamma_M)} \geq \frac{\mathcal{H}^n_k(\bar{U})}{\mathcal{H}^{n-1}_k(\bar{\Gamma})}$$
Equality holds if and only if the solution is isometric to the model.

Theorem C: Quantitative Hot-Spot Localization

The paper provides a quantitative estimate for the distance between the maximum set of $u$ and the boundary $\partial\Omega$. Under convexity assumptions on $\Omega$ and specific conditions on $f$ (including $f(x) \leq nkx + 1$), the distance is bounded below by a function of the inradius and model parameters.
Rigidity Case: If the bound is sharp, the manifold must be a space form of curvature $k$, $f$ must be affine ($f(x) = nkx+1$), and the domain must be a geodesic ball (recovering Serrin's classical result for $k=0$).

Construction of Exotic Solutions (Section 4)

The authors demonstrate that on the sphere $S^n$, isoparametric foliations (level sets of Cartan-Münzner polynomials) can be used to construct non-rotational (exotic) $f$-extremal domains.

• These solutions are not radially symmetric but possess constant Neumann data on their boundaries.
• The construction descends to smooth quotients $S^n/G$ (e.g., Lens spaces, Projective spaces) whenever the isometry group $G$ preserves the foliation and acts freely.

4. Significance and Impact

• Generalization of Serrin's Theorem: The work extends the classical symmetry results of Serrin (which state that overdetermined problems on $\mathbb{R}^n$ imply spherical symmetry) to manifolds with Ricci curvature bounds, providing a unified framework for rigidity.
• Quantitative Estimates: Unlike previous qualitative results, this paper provides explicit, quantitative bounds for gradients, isoperimetric ratios, and the location of hot spots, dependent on the curvature bound $k$ and the nonlinearity $f$.
• New Geometric Tools: The introduction of the $\tau$-functions and the pseudo-radial map offers a robust mechanism for comparing solutions on general manifolds with model spaces, bypassing the need for explicit global symmetries in the domain.
• Exotic Solutions: The construction of non-rotational solutions on spheres and their quotients challenges the intuition that overdetermined problems always imply radial symmetry, highlighting the role of the ambient geometry (specifically isoparametric foliations) in determining solution structure.
• Applications: The results have implications for the study of static Einstein equations with positive cosmological constants, torsion problems, and the critical catenoid conjecture.

In summary, the paper establishes a powerful comparison principle for semilinear elliptic equations on manifolds with Ricci lower bounds, yielding sharp geometric inequalities and rigidity theorems while simultaneously constructing a rich class of non-symmetric solutions on spheres.