The Laplace Equation and Harmonic Functions

⚡ Key Takeaways

One-line statement: Laplace's equation $\nabla^2 u = 0$ characterizes harmonic functions — the smoothest, most "balanced" scalar fields, equal at every point to the average of their surroundings. Key result in words: A twice-differentiable function whose Laplacian vanishes has no local maxima or minima in the interior, is infinitely differentiable, and is uniquely determined by its boundary values. Where it applies: Steady-state heat flow, electrostatics, gravitation, incompressible fluid flow, complex analysis, and probability (Brownian motion). Who/when: Pierre-Simon Laplace (1780s, celestial mechanics); foundations laid by Euler, Gauss, Green, Riemann, and Poincaré through the 19th century.

💡 The Core Insight

A harmonic function is one with zero net curvature in the averaged sense: at every interior point its value equals the average of its values on any surrounding sphere. Curvature in one direction is exactly cancelled by opposite curvature in the perpendicular directions, so $u$ "bows up" in some directions only by "bowing down" equally in others. This single averaging property — the mean value property — generates almost everything else: smoothness, the maximum principle, uniqueness, and rigidity.

DefinitionThe Laplacian and Laplace's Equation

For a scalar function $u:\Omega \subseteq \mathbb{R}^n \to \mathbb{R}$ that is twice continuously differentiable ( $u \in C^2$ ), the Laplacian is

1

\quad \nabla^2 u \;=\; \Delta u \;=\; \sum_{i=1}^{n}\frac{\partial^2 u}{\partial x_i^2} \;=\; \frac{\partial^2 u}{\partial x_1^2} + \cdots + \frac{\partial^2 u}{\partial x_n^2}.

Derivation hint: The Laplacian is the divergence of the gradient, $\Delta u = \nabla\cdot(\nabla u)$ ; expand $\nabla\cdot\big(\partial_{x_1}u,\dots,\partial_{x_n}u\big)$ component-by-component.

Laplace's equation is the homogeneous PDE

2

\quad \Delta u = 0 \quad \text{on } \Omega.

Derivation hint: Set the Laplacian to zero; this is the source-free special case of Poisson's equation $\Delta u = f$ with $f \equiv 0$ . A function satisfying [2] on an open set $\Omega$ is harmonic on $\Omega$ . Here $\Omega$ is the open domain, $\partial\Omega$ its boundary, and the symbols $x_i$ are Cartesian coordinates.

DefinitionLaplacian in Other Coordinate Systems

In polar coordinates $(r,\theta)$ on $\mathbb{R}^2$ :

3

\quad \Delta u = \frac{\partial^2 u}{\partial r^2} + \frac{1}{r}\frac{\partial u}{\partial r} + \frac{1}{r^2}\frac{\partial^2 u}{\partial \theta^2}.

Derivation hint: Apply the chain rule to $x = r\cos\theta$ , $y = r\sin\theta$ , compute $\partial_x, \partial_y$ in terms of $\partial_r, \partial_\theta$ , square and add; cross terms cancel. In spherical coordinates $(r,\theta,\varphi)$ on $\mathbb{R}^3$ (with $\theta$ the polar angle):

4

\quad \Delta u = \frac{1}{r^2}\frac{\partial}{\partial r}\!\left(r^2\frac{\partial u}{\partial r}\right) + \frac{1}{r^2\sin\theta}\frac{\partial}{\partial\theta}\!\left(\sin\theta\,\frac{\partial u}{\partial\theta}\right) + \frac{1}{r^2\sin^2\theta}\frac{\partial^2 u}{\partial\varphi^2}.

Derivation hint: Use $\Delta u = \frac{1}{\sqrt{g}}\partial_i(\sqrt{g}\,g^{ij}\partial_j u)$ with the spherical metric $g = \mathrm{diag}(1, r^2, r^2\sin^2\theta)$ .

TheoremMean Value Property

Let $u$ be harmonic on an open set $\Omega \subseteq \mathbb{R}^n$ . For any ball $\overline{B(x_0,r)} \subset \Omega$ , the value at the center equals the average over the sphere and over the ball:

5

\quad u(x_0) = \frac{1}{|\partial B(x_0,r)|}\int_{\partial B(x_0,r)} u\,dS = \frac{1}{|B(x_0,r)|}\int_{B(x_0,r)} u\,dV.

Derivation hint: Define $\phi(r)$ as the spherical average; differentiate under the integral, rescale to the unit sphere, and apply the divergence theorem to $\int_{B}\Delta u\,dV = 0$ to show $\phi'(r)=0$ . Conversely, any continuous function satisfying the mean value property on every ball in $\Omega$ is harmonic.

TheoremMaximum Principle

If $u$ is harmonic on a bounded open set $\Omega$ and continuous on $\overline{\Omega}$ , then $u$ attains both its maximum and minimum on the boundary:

6

\quad \max_{\overline{\Omega}} u = \max_{\partial\Omega} u, \qquad \min_{\overline{\Omega}} u = \min_{\partial\Omega} u.

Derivation hint: If $u$ attained an interior maximum $M$ at $x_0$ , the mean value property [5] forces $u \equiv M$ on a surrounding sphere; a connectedness argument propagates $u\equiv M$ throughout $\Omega$ . The strong maximum principle states that if a non-constant harmonic function on a connected domain attains an interior maximum, that is impossible — interior extrema occur only for constants.

TheoremUniqueness for the Dirichlet Problem

The boundary value problem $\Delta u = 0$ in $\Omega$ , $u = g$ on $\partial\Omega$ (with $g$ continuous, $\Omega$ bounded) has at most one solution. Derivation hint: If $u_1, u_2$ both solve it, then $w = u_1 - u_2$ is harmonic with $w=0$ on $\partial\Omega$ ; the maximum principle [6] forces $\max w = \min w = 0$ , so $w\equiv 0$ .

ProofMean Value Property from Laplace's Equation

Step 1. Fix $x_0 \in \Omega$ and define the spherical average

\phi(r) = \frac{1}{|\partial B(x_0,r)|}\int_{\partial B(x_0,r)} u(y)\,dS(y).

Step 2. Change variables to the unit sphere via $y = x_0 + r z$ , $|z| = 1$ , $dS(y) = r^{n-1}dS(z)$ :

\phi(r) = \frac{1}{|\partial B(0,1)|}\int_{|z|=1} u(x_0 + rz)\,dS(z).

Step 3. Differentiate under the integral sign:

\phi'(r) = \frac{1}{|\partial B(0,1)|}\int_{|z|=1} \nabla u(x_0 + rz)\cdot z\,dS(z).

Here $z$ is the outward unit normal of the unit sphere, so $\nabla u \cdot z = \partial u/\partial \nu$ is the normal derivative.

Step 4. Convert back to the sphere of radius $r$ and apply the divergence theorem:

\phi'(r) = \frac{1}{|\partial B(x_0,r)|}\int_{\partial B(x_0,r)} \frac{\partial u}{\partial \nu}\,dS = \frac{1}{|\partial B(x_0,r)|}\int_{B(x_0,r)} \Delta u\,dV.

Step 5. Since $u$ is harmonic, $\Delta u = 0$ , hence $\phi'(r) = 0$ . Thus $\phi$ is constant in $r$ .

Step 6. As $r \to 0^+$ , continuity gives $\phi(r) \to u(x_0)$ . Therefore $\phi(r) = u(x_0)$ for all admissible $r$ , proving the spherical mean value property [5]. Integrating the spherical version over $r$ yields the solid (ball) version. $\blacksquare$

∎

ProofSmoothness of Harmonic Functions

Step 1. Let $u$ satisfy the mean value property. Choose a smooth radial mollifier $\eta_\varepsilon$ supported in $B(0,\varepsilon)$ with $\int \eta_\varepsilon = 1$ and depending only on $|x|$ .

Step 2. Form the convolution $u_\varepsilon = u * \eta_\varepsilon$ , which is automatically $C^\infty$ wherever defined.

Step 3. Because $\eta_\varepsilon$ is radial and $u$ satisfies the mean value property over every sphere, the weighted average $u_\varepsilon(x)$ collapses to $u(x)$ itself: $u_\varepsilon = u$ on the interior at distance $> \varepsilon$ from $\partial\Omega$ .

Step 4. Hence $u$ agrees locally with a $C^\infty$ function, so $u \in C^\infty(\Omega)$ . In fact $u$ is real-analytic. $\blacksquare$

∎

PL

Pierre-Simon Laplace

1749–1827

Laplace introduced the operator now bearing his name in his monumental Traité de mécanique céleste, where he studied the gravitational potential of celestial bodies. He showed that outside any mass distribution the potential $V$ satisfies $\Delta V = 0$ , reducing the calculation of gravitational fields to solving a single elegant equation. Famously, when Napoleon asked why his treatise made no mention of God, Laplace is said to have replied, "I had no need of that hypothesis." His work on potential theory, probability, and celestial mechanics shaped the analytic foundations of mathematical physics for the next century.

🧠 Intuition

Imagine a thin metal plate whose edge is held at fixed temperatures. After a long time the interior settles into a steady state: heat no longer flows on net, and the temperature distribution is harmonic. At any interior point the temperature equals the average of its neighbors — if it were higher, heat would flow out and cool it; if lower, heat would flow in and warm it. Equilibrium is the mean value property. This also explains why harmonic functions have no interior bumps. A local hot spot would be hotter than its surrounding average, violating [5]. So the hottest and coldest points always sit on the boundary. A taut soap film stretched across a wire loop is the same idea: with no pressure difference, the film's height is harmonic and the surface has no interior peaks or pits — it is a saddle everywhere or flat.

Figure

A harmonic saddle: curvature in

x

cancels curvature in

y

ExampleSteady Temperature on a Disk

Given: A circular plate of radius $R = 2$ m. The boundary temperature is held at $g(\theta) = 100\cos\theta$ °C. Find the steady-state temperature $u(r,\theta)$ inside, and its value at the center and at $(r,\theta) = (1, 0)$ .

Step 1 — General harmonic solution (Poisson/Fourier on a disk). Bounded harmonic functions on a disk have the form

u(r,\theta) = a_0 + \sum_{n=1}^{\infty} r^n\big(a_n\cos n\theta + b_n\sin n\theta\big),

because each $r^n\cos n\theta$ and $r^n\sin n\theta$ satisfies [3].

Step 2 — Match boundary data. At $r = R = 2$ we need $u(2,\theta) = 100\cos\theta$ . Only the $n=1$ cosine term survives:

a_1 R^1 \cos\theta = 100\cos\theta \;\Rightarrow\; a_1 \cdot 2 = 100 \;\Rightarrow\; a_1 = 50.

All other coefficients vanish.

Step 3 — Write the solution.

u(r,\theta) = 50\,r\cos\theta = 50\,x.

Step 4 — Evaluate. Center $(r=0)$ : $u = 0$ °C — consistent with the mean value property, since $\frac{1}{2\pi}\int_0^{2\pi}100\cos\theta\,d\theta = 0$ . At $(r,\theta)=(1,0)$ : $u = 50(1)\cos 0 = 50$ .

\boxed{u(r,\theta) = 50\,r\cos\theta \;\text{°C}, \quad u(\text{center}) = 0\ \text{°C}, \quad u(1,0) = 50\ \text{°C}}

ExampleVerifying and Building Harmonic Functions

Given: Test whether $u(x,y) = e^x\cos y$ is harmonic, and find its harmonic conjugate $v$ so that $f = u + iv$ is holomorphic.

Step 1 — Compute second derivatives.

u_{xx} = e^x\cos y, \qquad u_{yy} = -e^x\cos y.

Step 2 — Sum them.

\Delta u = u_{xx} + u_{yy} = e^x\cos y - e^x\cos y = 0. \;\checkmark \text{ Harmonic.}

Step 3 — Find conjugate via Cauchy–Riemann. Require $v_y = u_x = e^x\cos y$ and $v_x = -u_y = e^x\sin y$ . Integrate the first: $v = e^x\sin y + h(x)$ . Differentiate: $v_x = e^x\sin y + h'(x)$ . Compare with $v_x = e^x\sin y$ , so $h'(x)=0$ , $h = C$ .

Step 4 — Identify the holomorphic function.

f(z) = e^x\cos y + i\,e^x\sin y = e^x(\cos y + i\sin y) = e^{x+iy} = e^z.

\boxed{u = e^x\cos y \text{ is harmonic};\quad v = e^x\sin y + C;\quad f(z) = e^z}

⚠️ Common mistake

" $\Delta u = 0$ just means $u$ is a constant or a plane." False. Harmonic functions form an infinite-dimensional space: $x^2 - y^2$ , $xy$ , $e^x\cos y$ , $\ln r$ , $\arctan(y/x)$ , and infinitely many spherical harmonics are all non-trivial harmonic functions. What is true is that a harmonic function cannot have interior strict extrema.

"Harmonic implies bounded or implies decaying." Not on unbounded domains. On $\mathbb{R}^2$ , $u = x$ is harmonic and unbounded. Liouville's theorem only forbids functions that are both harmonic and bounded on all of $\mathbb{R}^n$ — those must be constant.

"The Laplacian is just the second derivative." In one dimension yes, but in higher dimensions it is a sum of second derivatives. The cancellation between directions (positive in one, negative in another) is precisely what makes the harmonic constraint non-trivial.

CounterexampleHarmonic in 2D but Not in 3D

The function $u(x,y) = \ln\sqrt{x^2+y^2} = \ln r$ is harmonic on $\mathbb{R}^2 \setminus \{0\}$ :

\Delta(\ln r) = \frac{d^2}{dr^2}\ln r + \frac{1}{r}\frac{d}{dr}\ln r = -\frac{1}{r^2} + \frac{1}{r}\cdot\frac{1}{r} = 0.

But the same logarithm of the radial coordinate is not harmonic in 3D — there the fundamental solution is $1/r$ instead:

\Delta\!\left(\frac{1}{r}\right) = 0 \text{ on } \mathbb{R}^3\setminus\{0\}, \quad \text{whereas} \quad \Delta(\ln r) = \frac{1}{r^2} \neq 0 \text{ in 3D}.

Lesson: The form of the fundamental solution depends critically on dimension — $\ln r$ in 2D, $r^{2-n}$ in dimension $n \ge 3$ . Naively reusing a 2D harmonic function in 3D fails.

connection

Laplace's equation is the hub of an enormous web of mathematics:

Complex Analysis: The real and imaginary parts of every holomorphic function are harmonic, linked by the Cauchy-Riemann Equations. Conformal maps preserve harmonicity, giving a powerful method to solve boundary value problems.
Poisson's Equation: The inhomogeneous version $\Delta u = f$ models fields with sources; Laplace's equation is the source-free case $f=0$ .
<a class="wiki-link" href="/wiki/the-heat-equation">The Heat Equation</a>: $u_t = \Delta u$ describes diffusion; its time-independent (steady-state) solutions are exactly the harmonic functions.
<a class="wiki-link" href="/wiki/ndli-iit-kgp-ph11001-physics-the-wave-equation">The Wave Equation</a>: $u_{tt} = c^2\Delta u$ shares the Laplacian but is hyperbolic, behaving very differently (no maximum principle, finite propagation speed).
<a class="wiki-link" href="/wiki/fourier-series">Fourier Series</a> and <a class="wiki-link" href="/wiki/separation-of-variables">Separation of Variables</a>: The principal tools for constructing explicit harmonic solutions in rectangles, disks, and spheres.
Brownian Motion: A harmonic function's value at a point equals the expected boundary value reached by a random walker started there — a probabilistic incarnation of the mean value property.

🔬 Where this appears

⏳ History

- 1752 — Euler writes down what is essentially the Laplacian condition for the velocity potential of incompressible fluid flow.

- 1782–1785 — Laplace systematically uses

\Delta V = 0

for gravitational potentials in his celestial mechanics; the equation acquires his name.

- 1813 — Poisson generalizes to

\Delta V = -4\pi\rho

(Poisson's equation), correctly handling regions containing mass or charge.

- 1828 — George Green introduces the Green's function and potential-theoretic identities (Green's theorems), enabling boundary integral solutions.

- 1839–1850s — Gauss develops the theory of the potential and the mean value / maximum principle; Dirichlet poses the boundary value problem now bearing his name.

- 1851 — Riemann uses the Dirichlet principle (minimizing the energy integral) to assert existence of harmonic functions, fueling his theory of complex functions.

- 1870 — Weierstrass exposes a gap in the Dirichlet principle (the minimizer need not exist), triggering decades of rigor in the calculus of variations.

- 1887–1900 — Schwarz, Neumann, and Poincaré develop rigorous existence methods (alternating method, sweeping-out / balayage).

- 1900–1923 — Hilbert rehabilitates the Dirichlet principle; Perron (1923) gives the elegant subharmonic-function method that resolves existence for general domains.

open-problem

While the theory of harmonic functions is classically complete, several deep and active threads remain:

Regularity of boundaries. The Dirichlet problem is solvable for nice domains, but boundary behavior can fail at irregular points (e.g., the tip of a sharp inward cusp — the classical Lebesgue spine). Characterizing solvability via the Wiener criterion connects to potential theory and capacity, and analogues in non-smooth and fractal domains are still studied.
Higher-order and nonlinear analogues. The biharmonic equation $\Delta^2 u = 0$ (plate bending) lacks a maximum principle, and the p-Laplacian $\nabla\cdot(|\nabla u|^{p-2}\nabla u) = 0$ generalizes harmonicity nonlinearly with subtle regularity theory.
Geometric harmonic maps. Harmonic maps between Riemannian manifolds generalize harmonic functions; existence, uniqueness, and singularity structure (bubbling) remain research frontiers.
Discrete and probabilistic harmonicity. On graphs and lattices, discrete harmonic functions connect to random walks, electrical networks, and the convergence of numerical schemes; sharp results on infinite graphs are ongoing.
Unique continuation. How strongly a harmonic function's local behavior controls it globally (and quantitative estimates thereof) ties into open questions in spectral geometry and nodal-set bounds (Yau's conjecture territory).

The Laplace Equation and Harmonic Functions

Formal Statement

Why It Holds / Derivation

Building Intuition

Worked Examples

Common Misconceptions

Connections

Real-World Applications

Historical Development

Subtleties & Limits

See Also