Serrin Overdetermined Problem

• The Serrin overdetermined problem is a rigidity question in elliptic PDEs where prescribing both Dirichlet and constant Neumann conditions forces the domain to be a ball.
• It employs methods like moving planes and Pohozaev-type identities, offering insights into shape optimization and quantitative stability.
• The theory extends to various settings—including nonlocal, quasilinear, and multi-phase problems—revealing deep connections between boundary data and domain geometry.

The Serrin overdetermined problem is the rigidity question asking when an elliptic boundary value problem with both Dirichlet and constant Neumann data can be solved on a domain. In its classical form, if a bounded $C^2$ domain $\Omega\subset\mathbb{R}^n$ admits $u\in C^2(\overline{\Omega})$ satisfying

$$\begin{cases} -\Delta u = 1 & \text{in } \Omega,\ u = 0 & \text{on } \partial\Omega,\ \partial_\nu u = c & \text{on } \partial\Omega, \end{cases}$$

then $\Omega$ is a ball and $u$ is the quadratic radial profile. The sign of $c$ depends on whether $\nu$ is taken as the inward or outward unit normal, so equivalent formulations differ by orientation. The problem is “overdetermined” because the Dirichlet problem already determines the solution, and the extra Neumann condition constrains the geometry of the domain. Over time, the subject has expanded from this prototype to semilinear, quasilinear, nonlocal, Hessian, multi-phase, manifold, rough-boundary, unbounded, and shape-optimization settings (Celentano et al., 2024, Cao et al., 3 Jun 2025).

1. Classical formulation and rigidity theorem

Serrin’s theorem is the benchmark result of the subject. In one standard normalization, if $\Omega\subset\mathbb{R}^n$ is connected, bounded, open, with $C^2$ boundary, and there exists $\Omega\subset\mathbb{R}^n$0 such that

$\Omega\subset\mathbb{R}^n$1

then $\Omega\subset\mathbb{R}^n$2 is a ball $\Omega\subset\mathbb{R}^n$3 and

$\Omega\subset\mathbb{R}^n$4

Equivalent versions appear with $\Omega\subset\mathbb{R}^n$5 or with the outward normal, changing only signs and normalization (Celentano et al., 2024, Cao et al., 3 Jun 2025).

The theorem extends beyond the torsion equation. Serrin’s original framework also covers semilinear problems of the form

$\Omega\subset\mathbb{R}^n$6

under structural assumptions on $\Omega\subset\mathbb{R}^n$7 and assuming that $\Omega\subset\mathbb{R}^n$8 has constant sign; the conclusion remains that $\Omega\subset\mathbb{R}^n$9 is a ball and $u\in C^2(\overline{\Omega})$0 is radial with monotone radial profile (Celentano et al., 2024, Fall et al., 2016). In this sense, the phrase “Serrin overdetermined problem” now designates a class of rigidity statements in which prescribing both $u\in C^2(\overline{\Omega})$1 and a flux-type boundary datum forces spherical symmetry.

The classical result also has a variational interpretation. In Euclidean space it characterizes stationary domains for torsional rigidity under volume-preserving deformations, and balls are precisely the stationary shapes in the one-phase setting (Fall et al., 2014). This variational viewpoint became one of the main routes by which overdetermined boundary conditions entered shape optimization.

2. Proof paradigms, weak formulations, and quantitative stability

The original proof uses the method of moving planes, inherited from Alexandrov’s reflection method for constant-mean-curvature hypersurfaces. In a fixed direction, one reflects caps of the domain across moving hyperplanes, compares the solution with its reflection, and uses the strong maximum principle together with Hopf’s boundary lemma and a refined corner-point lemma at the critical position. This remains the canonical qualitative proof strategy for many local Serrin-type problems (Celentano et al., 2024, Cao et al., 3 Jun 2025). A distinct route is Weinberger’s $u\in C^2(\overline{\Omega})$2-function method, based on an auxiliary quantity and Pohožaev-type identities, which is especially effective in the Laplacian torsion case (Cao et al., 3 Jun 2025).

A major recent development concerns rough boundaries. The classical theorem was long tied to $u\in C^2(\overline{\Omega})$3 domains, but it now extends to bounded indecomposable sets of finite perimeter satisfying a uniform upper density bound on the reduced boundary. In that formulation, the boundary condition is encoded distributionally as

$u\in C^2(\overline{\Omega})$4

with $u\in C^2(\overline{\Omega})$5 almost everywhere outside $u\in C^2(\overline{\Omega})$6. This framework covers bounded Lipschitz domains and even allows slit discontinuities, while preserving the conclusion that $u\in C^2(\overline{\Omega})$7 must be a ball (Figalli et al., 2024). The argument replaces classical boundary regularity with geometric-measure-theoretic analysis of reduced boundaries, Green functions, and a rough-domain version of Weinberger’s maximum-principle strategy.

Quantitative stability has developed in parallel with rigidity. In the fractional setting, for

$u\in C^2(\overline{\Omega})$8

with fractional normal derivative $u\in C^2(\overline{\Omega})$9, one has a stability estimate of the form

$$\begin{cases} -\Delta u = 1 & \text{in } \Omega,\ u = 0 & \text{on } \partial\Omega,\ \partial_\nu u = c & \text{on } \partial\Omega, \end{cases}$$0

where $$\begin{cases} -\Delta u = 1 & \text{in } \Omega,\ u = 0 & \text{on } \partial\Omega,\ \partial_\nu u = c & \text{on } \partial\Omega, \end{cases}$$1 measures the gap between inner and outer enclosing balls. For $$\begin{cases} -\Delta u = 1 & \text{in } \Omega,\ u = 0 & \text{on } \partial\Omega,\ \partial_\nu u = c & \text{on } \partial\Omega, \end{cases}$$2 domains this exponent improves to

$$\begin{cases} -\Delta u = 1 & \text{in } \Omega,\ u = 0 & \text{on } \partial\Omega,\ \partial_\nu u = c & \text{on } \partial\Omega, \end{cases}$$3

The proof develops a quantitative nonlocal moving-planes apparatus, including an antisymmetric barrier, a quantitative nonlocal maximum principle, and nonlocal versions of Hopf’s lemma and Serrin’s corner lemma (Dipierro et al., 2023).

The same stability theme appears in perturbative local theory. For small perturbations of the constant Neumann datum around the ball, one can solve the perturbed overdetermined problem by parameterizing the boundary as a normal graph over the sphere and applying a modified implicit function theorem for Banach triplets that tolerates loss of derivatives away from the base point. This yields existence, local uniqueness up to translation, and linear stability estimates for the shape of the domain (Gilsbach et al., 2021).

3. Generalized operators and ambient geometries

A central direction of generalization replaces the scalar Laplacian problem by more complex operators while retaining the overdetermined Dirichlet–Neumann structure. One extension treats a boundary-coupled elliptic system

$$\begin{cases} -\Delta u = 1 & \text{in } \Omega,\ u = 0 & \text{on } \partial\Omega,\ \partial_\nu u = c & \text{on } \partial\Omega, \end{cases}$$4

where each $$\begin{cases} -\Delta u = 1 & \text{in } \Omega,\ u = 0 & \text{on } \partial\Omega,\ \partial_\nu u = c & \text{on } \partial\Omega, \end{cases}$$5 with $$\begin{cases} -\Delta u = 1 & \text{in } \Omega,\ u = 0 & \text{on } \partial\Omega,\ \partial_\nu u = c & \text{on } \partial\Omega, \end{cases}$$6 locally Lipschitz and $$\begin{cases} -\Delta u = 1 & \text{in } \Omega,\ u = 0 & \text{on } \partial\Omega,\ \partial_\nu u = c & \text{on } \partial\Omega, \end{cases}$$7 nondecreasing, and the boundary nonlinearity $$\begin{cases} -\Delta u = 1 & \text{in } \Omega,\ u = 0 & \text{on } \partial\Omega,\ \partial_\nu u = c & \text{on } \partial\Omega, \end{cases}$$8 is nondecreasing in each variable and strictly increasing in at least one. Under these hypotheses, $$\begin{cases} -\Delta u = 1 & \text{in } \Omega,\ u = 0 & \text{on } \partial\Omega,\ \partial_\nu u = c & \text{on } \partial\Omega, \end{cases}$$9 is a ball and each $\Omega$0 is radially symmetric and strictly radially decreasing (Celentano et al., 2024). This is a genuine system analogue of Serrin’s theorem: the PDEs are decoupled in the interior, but the normal derivatives are coupled nonlinearly on the boundary.

Another broad extension concerns quasilinear, potentially degenerate elliptic equations,

$\Omega$1

with $\Omega$2 continuous, strictly increasing, and $\Omega$3. In arbitrary bounded domains with no regularity assumed on $\Omega$4, symmetry can still be recovered if $\Omega$5 and $\Omega$6 approach prescribed boundary values uniformly from the interior. The same continuous Steiner symmetrization method also gives a ring-domain analogue: in a connected domain $\Omega$7, the existence of a positive solution with appropriate two-sided overdetermined data forces both $\Omega$8 and $\Omega$9 to be balls, although not necessarily concentric (Cao et al., 3 Jun 2025). This sharply contrasts with the bounded simply connected case, where the entire domain is forced to be a single ball.

The theory also extends outside the uniformly elliptic regime. In dimension two, if $u$0 solves a rotationally invariant Hessian equation

$u$1

on a smooth bounded simply connected domain and satisfies

$u$2

then $u$3 is a disk and $u$4 is radial, even though $u$5 is not assumed elliptic (Gálvez et al., 2019). The analyticity and simple connectivity are essential: the statement fails on non-simply-connected domains and fails for merely $u$6 solutions.

The ambient geometry can also be changed. In convex cones, Serrin-type mixed problems with Dirichlet and constant Neumann data on the interior part of the boundary and homogeneous Neumann data on the conical part force the domain to be a spherical sector. This holds for a class of possibly degenerate operators of the form

$u$7

and, in space forms, for the operator $u$8, yielding geodesic spherical sectors (Ciraolo et al., 2018). On compact Riemannian manifolds, small Serrin-type domains exist as smooth perturbations of geodesic balls. Their centers are governed by scalar curvature: if $u$9 is a nondegenerate critical point of $c$0, the boundaries of the resulting domains form a smooth foliation near $c$1, and conversely any concentrating family of such domains must converge to a critical point of $c$2 (Fall et al., 2014).

4. Coupled, two-phase, and multi-phase problems

The two-phase Serrin problem replaces the Laplacian by a divergence-form operator with piecewise constant coefficient,

$c$3

in a pair of nested domains $c$4, together with

$c$5

In this setting one distinguishes an outer problem, where $c$6 is fixed and $c$7 is sought, and an inner problem, where $c$8 is fixed and the inclusion $c$9 is sought (Cavallina et al., 2021, Cavallina et al., 2018). The concentric configuration $\nu$0 always provides a radial solution, but the one-phase rigidity picture no longer persists globally.

Near a concentric two-phase configuration, local solvability can be proved by shape derivatives and an implicit function theorem. For fixed $\nu$1 and $\nu$2, if the contrast $\nu$3 avoids a finite resonance set $\nu$4, then every sufficiently small perturbation of the core induces a unique nearby perturbation of the outer boundary solving the two-phase overdetermined problem; numerically, this can be computed by minimizing a Kohn–Vogelius mismatch functional under a volume constraint (Cavallina et al., 2018). This already shows that the outer boundary need not remain spherical.

At resonant coefficient values, symmetry breaking occurs. Writing $\nu$5 near special values $\nu$6, the linearized operator loses invertibility in the spherical-harmonic mode of degree $\nu$7. Crandall–Rabinowitz bifurcation then yields nontrivial branches of solutions

$\nu$8

for $\nu$9, and

$\Omega\subset\mathbb{R}^n$0

for $\Omega\subset\mathbb{R}^n$1, where $\Omega\subset\mathbb{R}^n$2 parameterizes the outer boundary and $\Omega\subset\mathbb{R}^n$3 is a degree-$\Omega\subset\mathbb{R}^n$4 spherical harmonic (Cavallina et al., 2020). In this precise sense, the two-phase Serrin problem admits symmetry-breaking outer domains even with a spherical core.

Quantitative stability survives near the one-phase regime. If a two-phase configuration solves the overdetermined problem and either $\Omega\subset\mathbb{R}^n$5 is small or the inclusion volume $\Omega\subset\mathbb{R}^n$6 is small, then $\Omega\subset\mathbb{R}^n$7 is quantitatively close to a ball, with explicit power-law control of the radius gap $\Omega\subset\mathbb{R}^n$8 by $\Omega\subset\mathbb{R}^n$9 or $C^2$0, where

$C^2$1

These estimates imply nonexistence for the inner problem on a fixed non-ball when the contrast or inclusion size is too small (Cavallina et al., 2021).

The multi-phase problem exhibits a further structural change. For a layered conductivity

$C^2$2

one can ask whether constant normal derivatives of successive orders on the outer boundary force radiality. The answer depends sharply on the number of phases. For $C^2$3, a single overdetermined condition $C^2$4 on the outer boundary does not force concentric balls, but imposing both

$C^2$5

does, and this characterization is if and only if (Cavallina, 2023). For $C^2$6, the pattern collapses: there exist non-radial configurations for which

$C^2$7

so countably infinitely many overdetermined conditions on the outer boundary still do not force concentricity (Cavallina, 2023). This is one of the clearest demonstrations that Serrin-type rigidity is highly sensitive to how internal structure is encoded.

5. Unbounded, epigraph, periodic, and exterior configurations

On bounded domains, the classical theorem forces balls. On unbounded domains, that classification fails. For the original torsion equation

$C^2$8

there exist nontrivial unbounded domains that are periodic in some variables and radial in the others. These domains are of the form

$C^2$9

with $\Omega\subset\mathbb{R}^n$00 even, $\Omega\subset\mathbb{R}^n$01-periodic, and permutation-invariant in the periodic variables, and they bifurcate from generalized cylinders or slabs. In each period cell they are uniquely self-Cheeger relative to that cell (Fall et al., 2016). This gives a direct counterpoint to the bounded-domain ball theorem.

Epigraphs form a different unbounded class. For

$\Omega\subset\mathbb{R}^n$02

if $\Omega\subset\mathbb{R}^n$03 is bounded from below and $\Omega\subset\mathbb{R}^n$04 is a classical solution of

$\Omega\subset\mathbb{R}^n$05

then rigidity can return. For specified Allen–Cahn-type nonlinearities, one obtains that $\Omega\subset\mathbb{R}^n$06 must be an affine half-space and $\Omega\subset\mathbb{R}^n$07 must be one-dimensional in all dimensions $\Omega\subset\mathbb{R}^n$08; for locally Lipschitz $\Omega\subset\mathbb{R}^n$09, analogous results hold in dimensions $\Omega\subset\mathbb{R}^n$10; and when $\Omega\subset\mathbb{R}^n$11, a new monotonicity theorem $\Omega\subset\mathbb{R}^n$12 is proved in any dimension $\Omega\subset\mathbb{R}^n$13 under $\Omega\subset\mathbb{R}^n$14 hypotheses on the epigraph (Beuvin et al., 7 Feb 2025). These results partially answer the Berestycki–Caffarelli–Nirenberg question in the epigraph setting.

Exterior domains introduce yet another geometry. For Hessian equations in $\Omega\subset\mathbb{R}^n$15,

$\Omega\subset\mathbb{R}^n$16

with

$\Omega\subset\mathbb{R}^n$17

and prescribed quadratic asymptotics at infinity, two regimes appear. If the asymptotic quadratic form is isotropic, $\Omega\subset\mathbb{R}^n$18, there is a threshold $\Omega\subset\mathbb{R}^n$19 such that for $\Omega\subset\mathbb{R}^n$20 the unique bounded domain is a ball and the solution is radial; if $\Omega\subset\mathbb{R}^n$21, no strictly convex solution exists (Wang et al., 2024). By contrast, for general anisotropic $\Omega\subset\mathbb{R}^n$22, there exists a unique bounded domain and a unique strictly convex solution, but the domain need not be spherical. The construction proceeds through a Legendre transform to a dual Hessian-quotient problem on the complement of a ball (Wang et al., 2024). This anisotropic exterior theory differs sharply from the bounded Serrin problem, where anisotropy of this type does not produce non-spherical domains.

6. Variational interpretations and current frontiers

A recurring principle is that overdetermined boundary conditions are Euler–Lagrange conditions for constrained shape problems. For the boundary-coupled system setting, critical domains under a volume constraint for

$\Omega\subset\mathbb{R}^n$23

satisfy the boundary relation

$\Omega\subset\mathbb{R}^n$24

which is exactly a generalized Neumann condition for a two-component system. The resulting Serrin-type theorem implies that every such critical shape is a ball (Celentano et al., 2024). On compact manifolds, the same variational viewpoint links small-volume torsional rigidity to scalar curvature: the small-volume profile is governed at first correction order by $\Omega\subset\mathbb{R}^n$25, and stationary small domains concentrate near critical points of $\Omega\subset\mathbb{R}^n$26 (Fall et al., 2014).

The perturbative theory around the ball has also become precise. For the classical torsion problem with perturbed Neumann datum

$\Omega\subset\mathbb{R}^n$27

one can parameterize the boundary by a normal graph $\Omega\subset\mathbb{R}^n$28 over the sphere and solve the nonlinear boundary-shape equation by a modified implicit function theorem that handles derivative loss except at the reference ball. This yields local existence of a unique domain up to translation and an optimal linear estimate

$\Omega\subset\mathbb{R}^n$29

together with a detailed spectral description of the translation kernel through spherical harmonics and the Dirichlet-to-Neumann operator (Gilsbach et al., 2021). This is a local, constructive complement to global rigidity theorems.

Several frontiers remain open or are explicitly suggested by the current literature. Natural problems include extending system results to PDE-coupled rather than only boundary-coupled systems, treating quasilinear or anisotropic systems with generalized Neumann data, weakening positivity and regularity assumptions in moving-planes arguments, classifying all unbounded or exterior domains that support Serrin-type solutions, and determining sharp thresholds in anisotropic exterior Hessian problems (Celentano et al., 2024, Beuvin et al., 7 Feb 2025, Wang et al., 2024). The modern theory therefore presents two simultaneous themes: a persistent rigidity principle—extra boundary data often forces canonical geometry—and a growing catalogue of precisely identified regimes in which that rigidity either weakens or breaks.