Global solution curves in harmonic parameters, and multiplicity of solutions

Technical Summary: Global Solution Curves in Harmonic Parameters, and Multiplicity of Solutions

Problem Statement
The paper investigates the semilinear Dirichlet problem on a smooth, bounded domain $\Omega \subset \mathbb{R}^m$:
$$\Delta u + g(u) = f(x), \quad x \in \Omega, \quad u = 0 \text{ on } \partial\Omega.$$
While the problem appears parameter-free, the author introduces "hidden" parameters based on the harmonic decomposition of the source term $f(x)$ and the solution $u(x)$. Specifically, $f(x)$ is decomposed as $f(x) = \sum_{i=1}^n \mu_i \phi_i + e(x)$, where $\phi_k$ are eigenfunctions of the Laplacian and $e(x)$ is orthogonal to the first $n$ eigenfunctions. Similarly, the solution is decomposed as $u(x) = \sum_{i=1}^n \xi_i \phi_i + U(x)$, with $U \perp \phi_i$.

The central objective is to analyze the solution set by studying the mapping from the "signature" of the solution $\xi = (\xi_1, \dots, \xi_n)$ to the corresponding parameters $\mu = (\mu_1, \dots, \mu_n)$. The paper focuses particularly on the principal solution curve $\mu_1 = \mu_1(\xi_1)$, where only the first harmonic is varied, to determine the multiplicity of solutions.

Methodology
The primary methodology involves a Lyapunov-Schmidt reduction combined with continuation in harmonic parameters.

1. Decomposition: The problem is split into a finite-dimensional part (governing the coefficients $\xi_i$ and $\mu_i$) and an infinite-dimensional part (governing the orthogonal component $U(x)$).
2. Continuation: Using the Implicit Function Theorem, the author constructs continuous solution curves $(u(x), \mu)(\xi)$ by varying the signature $\xi$. This is achieved by embedding the problem into a family parameterized by $k \in [0,1]$ or directly by $\xi$.
3. A Priori Estimates: A critical component of the analysis is establishing uniform bounds on the orthogonal component $U(x)$ in $H^2(\Omega)$ (or $H^1_0(\Omega)$). These bounds ensure that solution curves do not blow up at finite parameter values, allowing for global continuation.
4. Asymptotic Analysis: For the one-dimensional case, the stationary phase method is employed to derive precise asymptotic formulas for the solution curves as the harmonic amplitude $\xi_1 \to \pm \infty$.

Key Contributions and Results

• Global Solution Curves and Multiplicity:
  Under the condition $g'(u) < \lambda_2$ (where $\lambda_2$ is the second eigenvalue of the Laplacian), the paper proves that all solutions of the problem lie on a single, continuous solution curve parameterized by $\xi_1$. The multiplicity of solutions for a given $f(x)$ is determined by the number of times the horizontal line $\mu_1 = \text{const}$ intersects this curve.

  • If $\mu_1(\xi_1) \to \infty$ as $\xi_1 \to \pm \infty$, the curve has a global minimum $\mu_0$. The problem has zero solutions if $\mu_1 < \mu_0$, at least one if $\mu_1 = \mu_0$, and at least two if $\mu_1 > \mu_0$.
  • This result extends the work of Amann and Hess [1] by not requiring the existence of limits for $g(u)/u$ at infinity, and recovers the Ambrosetti-Prodi theorem [3] (via Berger and Podolak [6]) when $g'(u)$ is strictly monotone.

• Specific Equation Classes:

  • Cubic Nonlinearity: For $\Delta u + \lambda u - u^3 = \mu_1 \phi_1 + e(x)$ with $\lambda \in (0, \lambda_2)$, the solution curve is shown to be either monotone or S-shaped. If $\lambda \le \lambda_1$, the curve is monotone (unique solution). If $\lambda \in (\lambda_1, \lambda_2)$ and $e(x)$ is small, the curve makes at least two turns.
  • Resonant Problems: For $\Delta u + \lambda_1 u + g(u) = \mu_1 \phi_1 + e(x)$, the paper establishes the existence of a continuous solution curve exhausting the solution set, generalizing previous results by de Figueiredo, Ni, Iannacci, Nkashama, and Ward to cases where $\mu_1 \neq 0$.

• One-Dimensional Oscillatory Equations:
  The paper provides detailed results for the 1D resonant problem $u'' + \lambda_1 u + h(u)\sin u = \mu_1 \phi_1 + e(x)$.

  • Case $p \in (1/2, 1)$: If $h(u) \sim u^p$, the paper proves the existence of infinitely many solutions for any $\mu_1 \in \mathbb{R}$. An asymptotic formula is derived: $\mu_1(\xi_1) \sim \frac{2\sqrt{2}}{\sqrt{\pi}} \xi_1 \sin(\xi_1 - \pi/4) h(\xi_1)$. The amplitude of oscillation grows unbounded, ensuring infinite intersections with any horizontal line.
  • Case $p \in (0, 1/2)$: The problem has infinitely many solutions only when $\mu_1 = 0$. For $\mu_1 \neq 0$, there are at most finitely many solutions. The solution exhibits a universal asymptotic behavior $u(x) \sim \xi_1 \phi_1 + E(x)$, where $E(x)$ is the unique solution to the linearized problem with zero first harmonic.
  • Higher Eigenvalues: Similar results are obtained for resonance at higher eigenvalues $k$, showing that for $\mu_k = 0$, there are infinitely many solutions, while for large $|\mu_k|$, no solutions exist.

• Numerical Implementation:
  The author emphasizes that the approach of following solution curves in harmonic parameters is well-suited for numerical computation. The paper presents numerical computations (using a program based on Newton's method) that illustrate the theoretical curves and confirm the accuracy of the derived asymptotic formulas.

Significance
The paper claims that its primary significance lies in providing a unified framework for understanding solution multiplicity through the geometry of solution curves in harmonic parameter space. By demonstrating that all solutions lie on a single continuous curve (under specific spectral conditions), the complex problem of counting solutions is reduced to analyzing the shape of a planar curve $\mu_1(\xi_1)$.

The approach offers two distinct advantages:

1. Theoretical: It generalizes classical results (Ambrosetti-Prodi, Amann-Hess) by relaxing conditions on the nonlinearity $g(u)$, specifically removing the requirement for limits at infinity.
2. Computational: It facilitates the numerical computation of all solutions via curve-following algorithms, which are described as efficient and accurate. The paper explicitly notes that this method avoids the potential pitfalls of traditional bifurcation continuation (which may involve multiple turns) by using the harmonic parameter $\xi_1$ as a global parameter.

The work does not propose new experimental applications but rather refines the analytical and numerical tools available for studying semilinear elliptic problems, particularly in the context of resonance and multiplicity.