Serrin's Problem: Rigidity and Extensions

• Serrin’s problem is an overdetermined boundary value problem that investigates when a domain supports a Poisson-type equation with constant Dirichlet and Neumann conditions, resulting in a ball-shaped solution.
• It connects PDE symmetry, geometric analysis, and torsional rigidity using methods like moving planes, integral identities, and quantitative stability estimates to assess closeness to a ball.
• Recent advances extend the classical result to weak formulations, fully nonlinear variants, two-phase settings, and domains on manifolds or cones, deepening its impact in shape optimization and geometric analysis.

Serrin’s problem is the overdetermined boundary value problem that asks when a domain can support a solution of a Poisson-type equation together with both constant Dirichlet and constant Neumann data on the boundary. In its classical Euclidean form, one seeks a bounded domain $\Omega\subset\mathbb{R}^n$ and a function $u$ such that

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

or equivalently $u=0$ and $|\nabla u|=c$ on $\partial\Omega$, with $c>0$. Serrin’s theorem states that, for a bounded $C^2$ domain, the existence of such a solution forces $\Omega$ to be a ball and $u$ to be radial. This rigidity result became a prototype for overdetermined elliptic problems, linking PDE symmetry, geometric analysis, torsional rigidity, constant mean curvature phenomena, and shape optimization (Magnanini et al., 2017).

1. Classical formulation and core rigidity

The classical problem arises from the torsion equation for the Prandtl stress function. In the normalization used in several modern treatments, the torsion problem is

$u$0

and the overdetermined condition is $u$1 on $u$2. The condition is overdetermined because the Dirichlet problem already determines $u$3, so prescribing a Neumann datum on the whole boundary is not generically compatible. Serrin’s theorem identifies the unique compatible geometry: $u$4 must be a ball, and $u$5 is radial with respect to its center (Ciraolo et al., 2014).

This rigidity is classically proved by the method of moving planes, together with the strong maximum principle, Hopf’s lemma, and Serrin’s corner lemma. An alternative proof due to Weinberger uses a $u$6-function and integral identities. In the torsion normalization $u$7, one convenient reference radius is

$u$8

and the quadratic comparison function

$u$9

leads to the harmonic remainder $$\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, which plays a central role in integral-identity approaches to symmetry and stability (Magnanini et al., 2017).

The problem is tightly connected to torsional rigidity. For the torsion function $$\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, the torsional rigidity is represented by $$\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 or, depending on normalization, $$\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. In shape optimization, the condition that $$\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 be constant appears as the Euler–Lagrange condition for extremizers of torsional rigidity under a volume constraint (Cavallina, 2024).

2. Stability theory and quantitative symmetry

A major development after the qualitative theorem is quantitative stability: if the overdetermined condition is only approximately satisfied, then the domain is quantitatively close to a ball. A standard geometric measure of closeness uses inner and outer radii centered at a point $$\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,

$$\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

so that $$\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 measures the thickness of the smallest concentric annulus containing $$\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 (Magnanini et al., 2019).

For Serrin’s problem, a representative sharp-type estimate is

$$\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

with

$u=0$0

where $u=0$1 solves the torsion equation and $u=0$2 (Magnanini et al., 2019). In the $u=0$3 setting, the exponent is halved: $u=0$4 These estimates are described as nearly optimal; in dimension $u=0$5, the exponent $u=0$6 is optimal, as shown by explicit ellipse computations (Magnanini et al., 2019).

A related integral-identity framework yields improved stability simultaneously for Serrin’s problem and Alexandrov’s Soap Bubble theorem. In that approach, the interior defect

$u=0$7

is the common analytic quantity governing both the deviation of $u=0$8 from a constant and the deviation of the mean curvature $u=0$9 from its reference value $|\nabla u|=c$0. This leads, for example, to

$|\nabla u|=c$1

with $|\nabla u|=c$2, $|\nabla u|=c$3 arbitrarily close to $|\nabla u|=c$4, and $|\nabla u|=c$5 for $|\nabla u|=c$6 (Magnanini et al., 2017).

A distinct quantitative direction connects Serrin stability to dynamics. For the normalized overdetermined problem with $|\nabla u|=c$7 on $|\nabla u|=c$8, the defect

$|\nabla u|=c$9

controls the symmetric-difference distance to the optimal ball linearly under geometric assumptions such as a $\partial\Omega$0 boundary, an interior ball condition, and an $\partial\Omega$1-John condition (Feldman, 2017). That defect is exactly the energy dissipation rate for a quasi-static capillary drop model, which permits exponential convergence of regular solutions toward the spherical steady state (Feldman, 2017).

The stability viewpoint also appears in a partial-data variant. If the equation is known only in $\partial\Omega$2, for some $\partial\Omega$3, and the overdetermined condition is imposed on $\partial\Omega$4, then $\partial\Omega$5 is quantitatively close to a ball when $\partial\Omega$6 or $\partial\Omega$7 is small. The estimates control an $\partial\Omega$8 pseudo-distance to a sphere, Fraenkel-type asymmetry, and the annular thickness $\partial\Omega$9, thereby recovering the classical Serrin theorem when $c>0$0 (Dipierro et al., 2020).

3. Regularity thresholds and weak formulations

The classical theorem assumes a $c>0$1 boundary. A major recent result shows that this regularity is far from optimal. In a weak formulation, one extends $c>0$2 by zero outside $c>0$3 and encodes the overdetermined condition distributionally as

$c>0$4

for a bounded indecomposable set of finite perimeter $c>0$5. Under the density bound

$c>0$6

the only possible solution set is a ball, and $c>0$7 is the explicit quadratic torsion function (Figalli et al., 2024). This applies in particular to bounded Lipschitz domains and settles affirmatively the open question of whether Serrin’s theorem remains valid in that class (Figalli et al., 2024).

The same framework also accommodates slit discontinuities. If $c>0$8 is an $c>0$9-rectifiable interior slit and

$C^2$0

together with the symmetric blow-up condition

$C^2$1

then the same spherical rigidity holds (Figalli et al., 2024).

In the plane, weak formulations can be recast as harmonic quadrature identities. For a Jordan domain $C^2$2 with rectifiable boundary, the identity

$C^2$3

is equivalent, on Jordan quasidisks, to the weak Serrin condition (Zhang, 3 May 2026). Within the class of rectifiable Jordan Smirnov domains, this identity forces $C^2$4 to be a disk (Zhang, 3 May 2026). The Smirnov assumption is sharp: there exist rectifiable, non-Smirnov Jordan domains satisfying the same quadrature identity, hence supporting a weak Serrin formulation without being disks (Zhang, 3 May 2026). This exhibits a precise boundary-regularity threshold in the planar weak theory.

4. Fully nonlinear, two-phase, and non-elliptic variants

A broad generalization replaces the Laplacian by a fully nonlinear Hessian operator. In the planar real-analytic setting, one considers

$C^2$5

where $C^2$6 is rotationally invariant and not locally zero. If $C^2$7 is smooth, bounded, and simply connected, and $C^2$8 is nontrivial, then $C^2$9 is a disk and $\Omega$0 is radial, even though $\Omega$1 is not assumed elliptic (Gálvez et al., 2019).

The proof mechanism is completely different from moving planes. Rotational invariance reduces the equation to a functional dependence between $\Omega$2 and $\Omega$3, yielding

$\Omega$4

and the analysis proceeds via the eigenline fields of the Hessian and a Poincaré–Hopf index argument on a simply connected planar domain (Gálvez et al., 2019). The result is sharp in two separate senses: it fails if $\Omega$5 is not simply connected, and it fails if $\Omega$6 is merely $\Omega$7 rather than real analytic (Gálvez et al., 2019).

A different extension is the two-phase Serrin problem, motivated by shape optimization for composite torsional rigidity. Here one has a core $\Omega$8 with piecewise-constant conductivity

$\Omega$9

and solves

$u$0

together with transmission conditions across $u$1 (Cavallina, 2024). In this setting, the outer shape $u$2 is itself unknown.

Several qualitative features distinguish the two-phase problem from the one-phase case. Critical shapes of the two-phase torsional rigidity under a volume constraint are exactly those satisfying the overdetermined condition on $u$3, but such critical shapes are never local minimizers (Cavallina, 2024). The same work proves that solutions have no tentacles, the outer boundary contains no flat parts, and if the outer boundary contains a spherical portion then the only possibility is concentric balls (Cavallina, 2024). A strong parameter-rigidity result also holds: if a given configuration solves the two-phase problem for two distinct conductivity values $u$4, then the configuration must be concentric balls (Cavallina, 2024).

5. Curved ambient spaces, cones, and weighted analogues

On Riemannian manifolds, Serrin-type problems take the form

$u$5

with $u$6. Under $u$7, a closed conformal vector field, and a curvature-weighted integral condition, one obtains rigidity: $u$8 is a metric ball and $u$9 is radial (Freitas et al., 2023). In the Einstein case $u$00, the curvature condition is automatic, so any positive solution of the overdetermined problem on such a manifold forces $u$01 to be a metric ball (Freitas et al., 2023).

The principal analytic tools in that setting are a new Pohozaev identity involving scalar curvature,

$u$02

and a generalized Weinberger $u$03-function,

$u$04

which is subharmonic under the Ricci lower bound (Freitas et al., 2023).

A related line of work treats convex cones in warped product manifolds. For a sector-like domain $u$05 in a convex cone $u$06, the overdetermined problem becomes

$u$07

where $u$08 and $u$09 (Araújo et al., 9 Jan 2025). Under a Ricci lower bound and a compatibility condition involving the scalar curvature and the closed conformal field $u$10, the only possible domains are intersections of geodesic balls with the cone (Araújo et al., 9 Jan 2025). In Einstein warped products the compatibility condition is automatic (Araújo et al., 9 Jan 2025).

The same setting supports cone analogues of Alexandrov’s Soap Bubble theorem and of the Heintze–Karcher inequality, characterizing $u$11 among sector-like domains. In the Euclidean cone case with drift Laplacian

$u$12

for a homogeneous weight $u$13, the weighted overdetermined problem likewise forces the domain to be an intersection of a ball with the cone (Araújo et al., 9 Jan 2025).

Warped products without cone structure were also treated by a direct Weinberger $u$14-function method. For

$u$15

in a warped product manifold with $u$16, one defines

$u$17

Under suitable assumptions on the warping function $u$18, or under a compatibility condition

$u$19

the domain is a metric ball and $u$20 is radial (Farina et al., 2019). In model manifolds, the argument becomes strong enough to show that the metric inside the relevant ball must actually be the space-form metric of curvature $u$21 (Farina et al., 2019).

6. Unbounded, periodic, and planar classification phenomena

Bounded domains are rigid in the classical theorem, but unbounded domains exhibit a richer geometry. One construction produces nontrivial unbounded periodic domains $u$22 of the form

$u$23

where $u$24, $u$25 is even in each $u$26, $u$27-periodic in each $u$28, and invariant under permutations of the $u$29. These domains bifurcate from straight cylinders $u$30, and for small $u$31 one has

$u$32

with $u$33 orthogonal to the first harmonics (Fall et al., 2016). On each such $u$34, the overdetermined problem

$u$35

admits a solution (Fall et al., 2016). These domains are also uniquely self-Cheeger relative to a period cell, with relative Cheeger constant $u$36 (Fall et al., 2016).

A different unbounded phenomenon occurs for semilinear problems. For all $u$37, there exist smooth entire epigraphs

$u$38

which are not half-spaces and for which

$u$39

has a positive bounded solution (Pino et al., 2013). This gives a negative answer, in high dimension, to the Berestycki–Caffarelli–Nirenberg question for epigraphs (Pino et al., 2013). The construction uses large dilations of nontrivial minimal or constant-mean-curvature hypersurfaces and gluing methods, thereby reinforcing the conceptual analogy between Serrin’s problem and Alexandrov’s theorem on constant mean curvature surfaces (Pino et al., 2013).

The sphere provides another setting where non-classical domains appear. On $u$40, the spherical overdetermined problem

$u$41

has, besides geodesic balls and straight tubular neighborhoods of the equator, families of nontrivial Serrin domains bifurcating from symmetric tubular neighborhoods (Fall et al., 2016). These domains are parameterized by perturbations

$u$42

where $u$43 is an axially symmetric spherical harmonic of degree $u$44 and $u$45 is orthogonal to $u$46 (Fall et al., 2016). They are the first examples of Serrin domains in $u$47 not bounded by geodesic spheres (Fall et al., 2016).

In the planar doubly connected and periodic two-boundary setting, a recent classification identifies all smooth ring domains and periodic bands that solve Serrin’s classical problem with locally constant boundary data as algebro-geometric potentials of the mKdV hierarchy (Cerezo et al., 14 Jan 2026). If $u$48 is the developing map of the domain, the quantity

$u$49

satisfies a finite-gap relation

$u$50

where $u$51 are the mKdV hierarchy operators (Cerezo et al., 14 Jan 2026). This organizes planar Serrin ring domains and periodic bands into finite-dimensional complexity levels indexed by the spectral genus $u$52. At the first nontrivial, elliptic level, the theory produces a global one-parameter family of periodic Serrin bands interpolating between a flat band and a chain of tangent disks, together with, for each $u$53, a two-dimensional moduli space $u$54 of non-radial Serrin ring domains with dihedral symmetry of order $u$55 (Cerezo et al., 14 Jan 2026).

Taken together, these developments show that Serrin’s problem is not a single rigidity statement but a broad research program. In the bounded Euclidean one-phase case it characterizes balls. In rough domains it remains rigid down to finite-perimeter settings with density control. In real-analytic fully nonlinear planar problems it survives beyond ellipticity. On manifolds and cones it interacts with curvature, closed conformal vector fields, and Reilly-type formulas. In two-phase media, periodic settings, spheres, and planar doubly connected domains, it reveals a substantial flexibility structured by spectral, geometric, and integrable mechanisms.