---
title: Serrin-Type Overdetermined Systems
url: https://www.emergentmind.com/topics/serrin-type-overdetermined-system
type: topic
---

# Serrin-Type Overdetermined Systems

A Serrin-type overdetermined system is an elliptic boundary value problem in which one prescribes, on the same boundary or on a distinguished part of it, both Dirichlet data and a Neumann-type datum. In its classical torsion form, the prototype is
\[
\Delta u=-1 \quad \text{in }\Omega,\qquad u=0 \quad \text{on }\partial\Omega,\qquad \partial_\nu u=c \quad \text{on }\partial\Omega,
\]
while semilinear, quasilinear, fully nonlinear, anisotropic, and geometric analogues replace the Laplacian, the source term, or the ambient space. The system is called overdetermined because a second-order elliptic equation ordinarily does not admit simultaneous global Dirichlet and Neumann prescriptions except in highly constrained geometries. The central theme of the subject is therefore rigidity: solvability often forces the domain to be a ball, a geodesic ball, a spherical sector, a Wulff shape, or another canonical symmetric object. At the same time, modern work shows that this rigidity can fail in unbounded, multiphase, weak, or high-dimensional settings, so the theory now sits at the intersection of elliptic PDE, geometric analysis, free-boundary theory, and shape optimization [1310.4528][2407.02293].

## 1. Canonical formulation and the classical rigidity paradigm

The classical bounded-domain model is Serrin’s theorem: if a bounded \(C^2\) domain admits a classical solution of the overdetermined Poisson problem, then the domain must be a Euclidean ball, and the solution is the corresponding radial quadratic profile. In the normalization used in the rough-domain formulation,
\[
\Delta u=-1 \quad \text{in }\Omega,\qquad u=0 \quad \text{on }\partial\Omega,\qquad \partial_\nu u=c>0 \quad \text{on }\partial\Omega,
\]
the profile is
\[
u(x)=\frac{R^2-|x-x_0|^2}{2n},\qquad R=nc,
\]
up to translation [2407.02293]. This is the PDE analogue of Alexandrov’s theorem characterizing compact embedded constant-mean-curvature hypersurfaces as spheres, and that analogy remains structurally important in later developments [1310.4528].

A standard semilinear generalization replaces the torsion equation by
\[
\Delta u+f(u)=0 \quad \text{in }\Omega,\qquad u=0 \quad \text{on }\partial\Omega,\qquad \partial_\nu u=c \quad \text{on }\partial\Omega,
\]
with \(f\) smooth and monostable:
\[
f(0)=f(1)=0,\qquad f(s)>0 \ \text{for } s\in(0,1),\qquad f'(1)<0.
\]
This includes the Fisher–Kolmogorov and Allen–Cahn nonlinearities. Under these hypotheses there is a unique increasing one-dimensional profile \(w\) solving
\[
w''+f(w)=0,\qquad w(0)=0,\qquad w(+\infty)=1,
\]
and the half-space solution \(u(x)=w(x_N)\) solves the overdetermined problem in \(\{x_N>0\}\) [1310.4528]. This one-dimensional model becomes the local normal profile in several gluing constructions.

The modern literature uses “Serrin-type” in a broader sense. The overdetermined datum may be prescribed only on a portion of the boundary, may be replaced by \(|\nabla u|=\text{const}\), may involve a nonlinear functional \(F(\partial_\nu u_1,\dots,\partial_\nu u_m)\), or may appear on an internal interface rather than the outer boundary [2402.09817][2311.18585]. The defining feature is not the precise PDE, but the coexistence of elliptic interior structure with boundary data that are too strong for generic domains.

## 2. Rigidity mechanisms and proof architectures

Two proof paradigms dominate the classical theory. The first is Serrin’s method of moving planes, which compares the solution with its reflections until a critical position is reached; internal tangency and corner-touching are then excluded by the strong maximum principle, Hopf’s lemma, and the corner lemma. This framework extends to semilinear systems
\[
-\Delta u_i=f_i(u_i)\quad\text{in }\Omega,\qquad
u_i=0\quad\text{on }\partial\Omega,\qquad
F\!\left(\frac{\partial u_1}{\partial \nu},\dots,\frac{\partial u_m}{\partial \nu}\right)=c\quad\text{on }\partial\Omega,
\]
provided \(F\) is nondecreasing in each variable and strictly increasing in at least one component; the conclusion is again that \(\Omega\) is a ball and each \(u_i\) is radial [2402.09817].

The second paradigm is Weinberger’s \(P\)-function method. On model manifolds for the torsion equation \(\Delta u=-1\), one introduces
\[
P:=m|\nabla u|^2+2u,
\]
uses the Bochner formula and Cauchy–Schwarz to show \(\Delta P\ge 0\), and combines subharmonicity with boundary information to force constancy of \(P\). In the Euclidean case this recovers Weinberger’s proof of Serrin’s theorem; on model manifolds an additional compatibility condition is required to obtain the same rigidity conclusion [1708.02032]. Variants of the same idea appear in rough domains, in Riemannian manifolds, in cones, in anisotropic problems, and in the Heisenberg group.

A third major mechanism is the Rellich–Pohozaev identity. For Hessian quotient equations and curvature quotient equations, the proofs proceed by pairing a \(P\)-function maximum principle with identities that force equality in Newton–Maclaurin inequalities. In the Euclidean Hessian quotient case, the model overdetermined problem is
\[
S_k(D^2u)=\frac{C_n^k}{C_n^l}S_l(D^2u)\quad\text{in }\Omega,\qquad
u=0,\quad u_\gamma=1\quad\text{on }\partial\Omega,
\]
and the combination of the Pohozaev identity and the \(P\)-function \(P=|Du|^2-2u\) yields the unit ball and the quadratic profile [2209.06268].

More recent work shows that rigidity can also be proved without moving planes. For obstacle problems with radial obstacle \(\psi\), the proof compares the solution in \(\Omega\) with radial obstacle-problem solutions in the largest and smallest concentric balls \(B_\rho(0)\subset\Omega\subset B_R(0)\), and uses superharmonicity to force \(\rho=R\) [2306.12124]. For nonsmooth bounded domains and ring-shaped domains with potentially degenerate quasilinear operators, continuous Steiner symmetrization is used instead. There the key functional is
\[
G(z):=\int_0^z g(s)\,ds,
\]
and local symmetry under continuous Steiner symmetrization ultimately forces the positivity set to be a ball, or both boundary components to be balls in the ring-shaped case [2506.02423].

## 3. Unbounded domains, high dimensions, and the failure of flatness

The sharpest departure from the classical bounded rigidity is the construction of nontrivial unbounded epigraphs. For every \(N\ge 9\) and every monostable \(f\) as above, there exist smooth entire epigraphs
\[
\Omega=\{x\in\mathbb R^N:\ x_N>F(x_1,\dots,x_{N-1})\},
\]
which are not half-spaces, yet support a positive bounded solution of
\[
\Delta u+f(u)=0 \quad \text{in }\Omega,\qquad
u=0 \quad \text{on }\partial\Omega,\qquad
\partial_\nu u=c \quad \text{on }\partial\Omega.
\]
The model boundary is a small perturbation of a large dilation of the Bombieri–De Giorgi–Giusti minimal graph in \(\mathbb R^9\), and the solution is asymptotically of the form \(u(x)\approx w(z)\), with \(z\) the signed distance to \(\partial\Omega\) [1310.4528]. This gives a negative answer, for epigraphs in dimensions \(N\ge 9\), to the Berestycki–Caffarelli–Nirenberg question asking whether such domains must be a half-space, a cylinder, or the complement of a cylinder.

The same work pushes the Alexandrov–Serrin analogy into lower-dimensional noncompact geometry. If a hypersurface \(T\) has suitable Jacobi nondegeneracy, a large dilation of \(T\) can be perturbed to a nearby boundary carrying a Serrin solution. This produces domains near large dilations of compact nondegenerate CMC hypersurfaces, periodic cylindrically bounded domains in \(\mathbb R^3\) with boundary close to Delaunay surfaces, and domains near minimal surfaces of finite total curvature in \(\mathbb R^3\), including the catenoid and Costa-type surfaces [1310.4528].

Exterior-domain versions introduce a different kind of flexibility. For the \(k\)-Hessian equation
\[
\sigma_k(\lambda(D^2u))=\frac{d!}{(d-k)!k!}\quad\text{in }\mathbb R^d\setminus\Omega,
\]
with interior-boundary conditions
\[
u=0,\qquad \frac{\partial u}{\partial \nu}=1 \quad\text{on }\partial\Omega,
\]
and prescribed quadratic asymptotics at infinity, one obtains a free-boundary problem in which the bounded obstacle \(\Omega\) is determined by the far-field data. In the isotropic case \(A=I\), the unique domain is a ball centered at \(-b\), but for anisotropic \(A\neq I\) the exterior domain need not be a ball [2412.11841]. This suggests that overdetermined rigidity in exterior geometries is governed jointly by boundary data and asymptotic structure, rather than by boundary data alone.

## 4. Ambient geometry: manifolds, cones, and sub-Riemannian models

Once the Euclidean ambient space is replaced, the rigidity object changes with the geometry. On rotationally symmetric model manifolds
\[
(M^m,g)=([0,R)\times\mathbb S^{m-1}/\sim,\ dr^2+\sigma^2(r)g_{\mathbb S^{m-1}}),
\]
the overdetermined torsion problem
\[
\Delta u=-1 \quad \text{in }\Omega,\qquad
u=0 \quad \text{on }\partial\Omega,\qquad
\partial_\nu u=c \quad \text{on }\partial\Omega
\]
forces \(\Omega\) to be an Euclidean ball centered at the pole, with
\[
u(r)=\frac{1}{2m}(\rho^2-r^2),
\]
provided \(\Omega\) contains the pole, the Ricci curvature is nonnegative in the model sense, \(\sigma'>0\), and the solution satisfies the compatibility condition
\[
\int_\Omega \frac{(\sigma''\,\sigma^{m-1})'}{\sigma^{m-1}}\,u^2 \le 0.
\]
In Euclidean space this condition is automatic because \(\sigma''=0\); in curved settings it is a genuine restriction [1708.02032].

In convex cones, the natural rigidity object is a spherical sector. For sector-like domains \(\Omega\subset E\subset\mathbb R^N\), the partially overdetermined problem
\[
L_f u=-1 \quad \text{in }\Omega,\qquad
u=0,\ \partial_\nu u=-c \quad \text{on }\Gamma,\qquad
\partial_\nu u=0 \quad \text{on }\Gamma_1\setminus\{0\},
\]
with
\[
L_f u=\operatorname{div}\!\big(f'(|\nabla u|)\nabla u\big),
\]
implies
\[
\Omega=E\cap B_R(x_0),
\]
so the domain is the intersection of the cone with a ball. The same conclusion holds in space forms for the equation \(\Delta u+NK\,u=-1\), yielding geodesic-ball sectors [1806.08553]. Related partially overdetermined results in convex cones and outside cones, especially in dimension \(2\), identify ball sectors and exterior ball sectors by a \(P\)-function and integral identities [2009.00280].

In more general Riemannian manifolds, one can recover metric balls under structural assumptions on the ambient geometry. If \((M^n,g)\) admits a conformal vector field and
\[
\Delta_g u + nk\,u = -1 \quad \text{in }\Omega,\qquad
u = 0 \quad \text{on }\partial\Omega,\qquad
|\nabla u| = c \quad \text{on }\partial\Omega,
\]
then a Pohozaev-type identity combined with Weinberger’s \(P\)-function
\[
P(u)=|\nabla u|^2+2u+ku^2
\]
forces \(\Omega\) to be a metric ball and \(u\) to be radial, provided \(\operatorname{Ric}_g\ge (n-1)k\,g\) and a compatibility condition involving the scalar curvature and the conformal factor holds [2405.17312].

In the Heisenberg group \(\mathbb H^n\), Euclidean balls are replaced by gauge balls. The family
\[
\Delta_H u=(Q+a-2)F_a \quad \text{in }\Omega,\qquad
u=0 \quad \text{on }\partial\Omega,\qquad
|\nabla_H u|=c\,F_a^{1/2}\quad \text{on }\partial\Omega,
\]
with
\[
F_a(\xi)=|x|^2p(\xi)^{a-4}=|\nabla_H p(\xi)|^2p(\xi)^{a-2},
\]
characterizes gauge balls within classes of cylindrically symmetric or toric symmetric domains; the model solution on \(B_R(0)\) is
\[
u_a(\xi)=\frac{p(\xi)^a-R^a}{a},\qquad c=R^{a-1}
\]
[2310.10389].

| Setting | Overdetermined problem | Rigidity object |
|---|---|---|
| Model manifolds | \(\Delta u=-1,\ u=0,\ \partial_\nu u=c\) | Euclidean ball centered at the pole |
| Convex cones / space forms | Partial Dirichlet–Neumann overdetermination | Spherical sector or geodesic-ball sector |
| Riemannian manifolds with conformal field | \(\Delta_g u+nku=-1,\ u=0,\ |\nabla u|=c\) | Metric ball |
| Heisenberg group | \(\Delta_H u=(Q+a-2)F_a,\ u=0,\ |\nabla_H u|=cF_a^{1/2}\) | Gauge ball |

## 5. Transmission media, free boundaries, and coupled variants

The two-phase Serrin problem replaces the Laplacian by a divergence-form operator with piecewise constant conductivity,
\[
-\operatorname{div}(\sigma_c\nabla u)=1,\qquad
\sigma_c(x)=1+(c-1)\chi_D(x),
\]
on a bounded domain \(\Omega\) containing an inclusion \(D\). The overdetermined condition is
\[
u=0 \quad \text{on }\partial\Omega,\qquad \partial_n u=c_0 \quad \text{on }\partial\Omega,\qquad
c_0=-\frac{|\Omega|}{|\partial\Omega|}.
\]
Near concentric balls \(D=B_R\), \(\Omega=B_1\), the outer problem is locally solvable and unique for \(\sigma_c\notin\Sigma\), where \(\Sigma=\{s(k):k=1,2,\dots\}\) is a finite resonance set defined explicitly in terms of \(R\), \(N\), and the spherical-harmonic degree \(k\). The proof uses shape derivatives and the implicit function theorem, and the numerical treatment minimizes a Kohn–Vogelius functional by a steepest descent algorithm [1811.07156].

At the exceptional conductivities \(\sigma_c=s(m)\), the linearized operator loses invertibility and symmetry breaking occurs. There are infinitely many bifurcating branches of nontrivial solutions \((B_R,\Omega_{g(\varepsilon)})\), one for each critical mode \(m\), with
\[
g(\varepsilon)=\varepsilon \cos(m\theta)+o(\varepsilon)\quad\text{in }N=2,
\qquad
g(\varepsilon)=\varepsilon Y_m(\theta)+o(\varepsilon)\quad\text{in }N\ge 3.
\]
Thus the two-phase problem admits nonradial outer boundaries even near a radial configuration [2001.10212].

Quantitative stability results show that this flexibility is limited near the one-phase regime. If \(v\) solves the one-phase torsion problem in \(\Omega\), then
\[
\rho_e-\rho_i \le C_1\|\partial_n v-c_0\|_{L^2(\partial\Omega)}^{\tau_N},
\]
with
\[
\tau_2=1,\qquad \tau_3=1-\theta \ \text{for any }\theta>0,\qquad \tau_N=\frac{2}{N-1}\quad (N\ge 4).
\]
For the two-phase problem this yields
\[
\rho_e-\rho_i \le C_2 |c-1|^{\tau_N},\qquad
\rho_e-\rho_i \le C_3 |D|^{\tau_N/2},
\]
so small contrast or small inclusion forces the outer boundary to be quantitatively close to a ball [2107.05889].

Other variants shift the overdetermination to incomplete or free-boundary settings. If the torsion equation is known only on \(\Omega\setminus\overline\omega\),
\[
\Delta u=1 \quad \text{in }\Omega\setminus\overline\omega,\qquad
u=0,\ u_\nu=c \quad \text{on }\partial\Omega,
\]
then a small unknown subdomain \(\omega\) still forces \(\Omega\) to be close to the sphere \(B_{Nc}(z)\), in boundary pseudodistance, Fraenkel-type asymmetry, and shell thickness \(\rho_e-\rho_i\) [2005.04859]. For obstacle problems,
\[
\min\{-\Delta u,\ u-\psi\}=0 \quad \text{in }\Omega,\qquad
u=0 \quad \text{on }\partial\Omega,\qquad
\partial_\nu u=c \quad \text{on }\partial\Omega,
\]
radial obstacle data imply that \(\Omega\) is a ball centered at the origin and \(u\) is radial; an approximate Neumann condition yields \(R-\rho\le K\varepsilon\) [2306.12124].

Partially overdetermined mixed problems in the half-space characterize capillary spherical caps. For
\[
\bar\Delta f=1 \quad\text{in }\Omega,\qquad
f=0 \quad\text{on }\Sigma,\qquad
\bar\nabla_{\bar N}f=c \quad\text{on }T,\qquad
f_\nu\equiv c_0 \quad\text{on }\Sigma,
\]
one obtains
\[
\Omega=B_{(n+1)c_0}(z)\cap\mathbb R^{n+1}_+,
\qquad
f(x)=\frac{|x-z|^2-(n+1)^2c_0^2}{2(n+1)},
\]
and the contact angle satisfies \(\cos\theta=-c/c_0\) [2311.18585].

## 6. Weak formulations, rough domains, and sharp thresholds

A major recent development is the extension of Serrin rigidity to rough sets. For a bounded indecomposable set of finite perimeter \(\Omega\subset\mathbb R^n\) satisfying a uniform upper density bound on the reduced boundary, if
\[
u\in W^{1,2}(\mathbb R^n),\qquad
u=0 \ \text{a.e. in }\mathbb R^n\setminus\Omega,
\]
and
\[
\Delta u = c\,\mathcal H^{n-1}\!\llcorner \partial^*\Omega - \mathbf 1_\Omega\,dx
\]
in the distributional sense, then \(\Omega\) is a ball and \(u\) is exactly the Serrin profile. The result applies in particular to Lipschitz domains, answers the previously unresolved Lipschitz case affirmatively, and even extends to certain slit-domain formulations [2407.02293].

An anisotropic Lipschitz-domain version replaces the Euclidean Laplacian by
\[
\Delta_H u=\operatorname{div}(H(Du)DH(Du)),
\]
where \(H\) is the Wulff potential associated with a bounded ellipsoid \(K\). If
\[
u \in W^{1,2}(\mathbb R^n),\qquad
u=0\ \text{a.e. in }\mathbb R^n\setminus\Omega,\qquad
\Delta_H u = \mathbf c\,\mathscr H^{n-1}|_{\partial^*\Omega} - \mathbf 1_\Omega\,dx,
\]
and \(u\) satisfies the gradient nondegeneracy and \(L^2\)-Dini-VMO assumptions near \(\partial^*\Omega\), then
\[
u(x)=\frac{r^2-H_*^2(x)}{2n}
\]
and \(\Omega\) is homothetic to \(K\) [2509.05155].

In the plane, the weak formulation reveals a sharp regularity threshold. For a bounded Jordan domain with rectifiable boundary, the weak Serrin condition is equivalent to the harmonic quadrature identity
\[
\int_\Omega h\,dA = c\int_{\partial\Omega} h\,ds
\qquad \forall\, h\in C(\overline\Omega),\ \Delta h=0 \text{ in }\Omega.
\]
If the domain is Smirnov, this identity forces \(\Omega\) to be a disk. But there also exist rectifiable Jordan non-Smirnov domains satisfying the same identity, showing that the Smirnov hypothesis is sharp [2605.02034]. This resolves a common misconception: weak Serrin-type data do not by themselves force disks among all rectifiable Jordan domains; an additional factorization condition on the conformal map is essential.

The subject has also acquired a constructive perturbative branch. For the perturbed planar problem
\[
-\Delta u+\mu w u=1 \quad \text{in }\Omega_\mu,\qquad
u=0,\quad \partial_\nu u=\frac12 \quad \text{on }\partial\Omega_\mu,
\]
with
\[
\Omega_\mu=\{(r,\theta):0\le r\le 1+\mu g(\theta)\},
\]
analytic expansions
\[
u=\sum_{k\ge 0}\mu^k u_k,\qquad
w=\sum_{k\ge 1}\mu^{k-1}w_k
\]
produce infinitely many approximants whose PDE and boundary defects are super-polynomially small in \(\mu\), via an optimal truncation inspired by Nekhoroshev theory [2603.19446]. This suggests a complementary viewpoint: beyond exact rigidity the overdetermined structure also supports highly accurate constructive near-solutions.

Taken together, these results show that Serrin-type overdetermined systems no longer denote a single theorem but a large rigidity–flexibility program. In smooth bounded Euclidean domains they still single out balls; in curved, anisotropic, and sub-Riemannian geometries they characterize the corresponding canonical shapes; in weak and rough settings they require refined measure-theoretic or conformal hypotheses; and in high-dimensional, multiphase, or asymptotic regimes they can support genuinely nonclassical geometries [1310.4528][2605.02034].

Source: https://www.emergentmind.com/topics/serrin-type-overdetermined-system