Laplacian

Laplacian operator in physical field equations, harmonic potentials, and diffusion processes.

What it is

The Laplacian is a second-order differential operator, denoted by $\nabla^2$ or $\Delta$, defined as the divergence of the gradient:

$\nabla^2 \phi = \vec{\nabla} \cdot (\vec{\nabla} \phi)$

In Cartesian coordinates it is the sum of pure unmixed second partial derivatives:

$\nabla^2 \phi = \frac{\partial^2 \phi}{\partial x^2} + \frac{\partial^2 \phi}{\partial y^2} + \frac{\partial^2 \phi}{\partial z^2}$

Physically, the Laplacian measures the local mean curvature or the difference between the central value of a function and its spherical average in an infinitesimal surrounding neighborhood:

$\nabla^2 \phi(P) = \lim_{R \to 0} \frac{6}{R^2} \left( \langle \phi \rangle_{S_R} - \phi(P) \right)$

• If $\nabla^2 \phi(P) > 0$: $\phi(P)$ is smaller than the average of its neighbors (local valley / concave upward).
• If $\nabla^2 \phi(P) < 0$: $\phi(P)$ is greater than the average of its neighbors (local peak / concave downward).
• If $\nabla^2 \phi(P) = 0$: $\phi(P)$ matches its neighborhood average identically. Solutions are called harmonic functions.

Why it appears in physics

The Laplacian is the mathematical engine of field equations:

• Poisson's equation: $\nabla^2 V = -\frac{\rho}{\varepsilon_0}$ (electrostatics) and $\nabla^2 \Phi = 4\pi G \rho$ (gravitation).
• Laplace's equation: $\nabla^2 \phi = 0$. Governs steady-state diffusion, vacuum potentials, and irrotational potential flows.
• Heat / diffusion equation: $\frac{\partial T}{\partial t} = \alpha \nabla^2 T$. Heat flows to eliminate local discrepancies from the mean.
• Wave equation: $\nabla^2 u = \frac{1}{c^2} \frac{\partial^2 u}{\partial t^2}$. Governs acoustics, light, and seismic waves.
• Schrödinger equation: $-\frac{\hbar^2}{2m} \nabla^2 \psi + V \psi = E \psi$. Quantum kinetic energy is proportional to the Laplacian of the wavefunction.

Visual or geometric intuition

Picture an elastic drum membrane:

• When undisturbed, its height profile satisfies Laplace's equation $\nabla^2 z = 0$: every point sits at the exact average of its surrounding ring.
• Under a transverse localized weight, curvature changes non-trivially ($\nabla^2 z \neq 0$), modeling Poisson's equation with a point load.

Minimal symbolic language

• Cartesian coordinates: $\nabla^2 = \partial_{xx} + \partial_{yy} + \partial_{zz}$.
• Cylindrical coordinates: $\nabla^2 \phi = \frac{1}{r} \frac{\partial}{\partial r}\left(r \frac{\partial \phi}{\partial r}\right) + \frac{1}{r^2} \frac{\partial^2 \phi}{\partial \theta^2} + \frac{\partial^2 \phi}{\partial z^2}$.
• Spherical coordinates: $\nabla^2 \phi = \frac{1}{r^2} \frac{\partial}{\partial r}\left(r^2 \frac{\partial \phi}{\partial r}\right) + \frac{1}{r^2 \sin\theta} \frac{\partial}{\partial \theta}\left(\sin\theta \frac{\partial \phi}{\partial \theta}\right) + \frac{1}{r^2 \sin^2\theta} \frac{\partial^2 \phi}{\partial \varphi^2}$.
• Vector Laplacian: $\nabla^2 \vec{A} = \vec{\nabla}(\vec{\nabla}\cdot\vec{A}) - \vec{\nabla}\times(\vec{\nabla}\times\vec{A})$.

What the result means

Non-zero $\nabla^2 \phi$ signals that $\phi$ deviates from local equilibrium. In diffusion, $\nabla^2 T > 0$ means the point will warm up; $\nabla^2 T < 0$ means it will cool down.

When not to use it

• Do not use Cartesian second derivatives in polar or spherical systems without metric scale factors.
• For vector fields in curvilinear coordinates, do not take the Laplacian component-by-component because unit vectors rotate with position. Use $\nabla^2 \vec{A} = \vec{\nabla}(\vec{\nabla}\cdot\vec{A}) - \vec{\nabla}\times(\vec{\nabla}\times\vec{A})$.

Operational summary

To compute $\nabla^2 \phi$: evaluate second partial derivatives and sum them. If $\nabla^2 \phi = 0$, the function obeys the mean-value property and has no local extrema in the interior (Earnshaw's theorem).

Procedure

Before calculating

Confirm whether the input is a scalar field $\phi$ or a vector field $\vec{A}$ and select the optimal coordinate system.

Step 1: Identify the mathematical object

Formulate the target equation: Laplace ($\nabla^2 \phi = 0$), Poisson, diffusion, or Helmholtz.

Step 2: Choose the operation

Select the appropriate Laplacian formula:

• Cartesian: $\partial_{xx} \phi + \partial_{yy} \phi + \partial_{zz} \phi$.
• Cylindrical: $\frac{1}{r}\partial_r(r\partial_r \phi) + \frac{1}{r^2}\partial_\theta^2 \phi + \partial_z^2 \phi$.
• Spherical: $\frac{1}{r^2}\partial_r(r^2\partial_r \phi) + \frac{1}{r^2\sin\theta}\partial_\theta(\sin\theta\partial_\theta \phi) + \frac{1}{r^2\sin^2\theta}\partial_\varphi^2 \phi$.

Step 3: Conserve units and signs

Perform inner first derivatives, multiply by coordinate scale weights, then evaluate outer derivatives.

Step 4: Calculate

Combine all terms and simplify. Substitute spatial coordinates for numerical evaluations.

Step 5: Interpret the result

Interpret physical significance: $\nabla^2 \phi = 0$ indicates harmonic equilibrium; non-zero values reveal local source density or diffusion rates.

Self-check checklist

• Is the output for a scalar field a scalar?
• Were curvilinear scale weights differentiated correctly?
• Are dimensions $[\phi]/[\text{m}^2]$?

Guided example

Physical situation

An infinitely long coaxial cable consists of an inner core of radius $a = 1.0 \text{ mm}$ held at voltage $V_a = 100 \text{ V}$ and an outer shield of radius $b = 5.0 \text{ mm}$ grounded at $V_b = 0 \text{ V}$. Space between ($a < r < b$) is vacuum.

1. Solve Laplace's equation $\nabla^2 V = 0$ for $V(r)$.
2. Calculate electric field $\vec{E}(r) = -\vec{\nabla} V$.

Mathematical translation

By cylindrical symmetry: $V = V(r)$. $\nabla^2 V = \frac{1}{r} \frac{d}{dr}\left( r \frac{dV}{dr} \right) = 0$

Step-by-step calculation

1. Integrate ODE: $r \frac{dV}{dr} = C_1 \implies \frac{dV}{dr} = \frac{C_1}{r} \implies V(r) = C_1 \ln(r) + C_2$

2. Impose boundary conditions: $V(a) = V_a$ and $V(b) = 0$: $C_1 = -\frac{V_a}{\ln(b/a)}, \quad C_2 = \frac{V_a \ln(b)}{\ln(b/a)}$ $V(r) = V_a \frac{\ln(b/r)}{\ln(b/a)}$

3. Numerical evaluation ($V_a = 100\text{ V}, b/a = 5$): $V(r) = 100 \frac{\ln(0.005/r)}{\ln(5)} \approx 62.133 \ln(0.005/r) \quad [\text{V}]$

4. Electric field: $\vec{E}(r) = -\frac{dV}{dr} \hat{e}_r = \frac{V_a}{r \ln(b/a)} \hat{e}_r = \frac{62.133}{r} \hat{e}_r \quad [\text{V/m}]$ At $r = a = 10^{-3}\text{ m}$: $E(a) \approx 62.13 \text{ kV/m}$.

Reading the result

Potential falls off logarithmically and field scales as $1/r$. Zero Laplacian confirms charge-free dielectric space.

Common error avoided

Using Cartesian $\frac{d^2 V}{dr^2} = 0$, which yields an invalid linear profile.

Inverse example

Given result

A spherical cloud yields a constant Laplacian $\nabla^2 V = -K$ inside $r < R$.

What physical question it can answer

What is the internal charge density? By Poisson's equation $\nabla^2 V = -\rho / \varepsilon_0 \implies \rho = K \varepsilon_0$. The constant Laplacian proves the sphere contains a perfectly uniform charge distribution.

Common errors

Error 1: Conflating the Laplacian with squared gradient magnitude

Confusing $\nabla^2 \phi$ with $\|\vec{\nabla}\phi\|^2$. The Laplacian is a linear sum of second derivatives, whereas the squared gradient norm is a nonlinear product of squared first derivatives.

Error 2: Omitting curvilinear metric curvature terms

Evaluating the spherical radial Laplacian as merely $\partial_{rr} \phi$. The correct radial term is $\partial_{rr}\phi + \frac{2}{r}\partial_r\phi$. Neglecting $\frac{2}{r}\partial_r$ invalidates the physical solution.

Error 3: Taking component-wise Laplacians in curvilinear vector fields

Writing $\nabla^2 \vec{A} = \sum_i (\nabla^2 A_i) \hat{e}_i$ in spherical or cylindrical frames. Curvilinear unit vectors rotate with coordinates; one must use $\nabla^2 \vec{A} = \vec{\nabla}(\vec{\nabla}\cdot\vec{A}) - \vec{\nabla}\times(\vec{\nabla}\times\vec{A})$.

Error 4: Expecting local extrema inside harmonic domains

Searching for electrostatic traps for free charged particles in vacuum. Earnshaw's theorem states that because $\nabla^2 V = 0$, stable isolated potential minima cannot exist in source-free space.

Error 5: Surprised that linear functions have zero Laplacian

Overlooking that any linear function $a x + b y + c z + d$ trivially satisfies Laplace's equation.

Quick detection rule

If a Laplacian formula contains cross derivatives $\partial_{xy}$ in Cartesian frames, or fails to carry physical units of $[\phi]/[\text{m}^2]$, it is mathematically incorrect.

Physical applications

1. Electrostatics and gravity: Poisson and Laplace formulations

• $\nabla^2 V = -\rho/\varepsilon_0$ (Coulomb potential).
• $\nabla^2 \Phi = 4\pi G \rho$ (Newtonian gravity).
• $\nabla^2 V = 0$ in vacuum (Earnshaw's theorem forbidding static charge traps).

2. Quantum mechanics: Schrödinger wave mechanics

The kinetic energy operator is $\hat{T} = -\frac{\hbar^2}{2m} \nabla^2$. The Laplacian dictates hydrogen atomic orbitals, solid-state electronic band structures, and tunneling.

3. Heat conduction and molecular diffusion

Fourier's heat equation $\partial_t T = \alpha \nabla^2 T$ and Fick's second law $\partial_t c = D \nabla^2 c$ describe relaxation toward maximum entropy.

4. Acoustics, optics, and electrodynamic waves

The d'Alembert wave equation $\nabla^2 u - \frac{1}{c^2} \partial_{tt} u = 0$ governs sound propagation, seismic waves, and electromagnetic light waves.

5. Elasticity and plate mechanics

Transverse deflections of elastic thin plates under load satisfy the biharmonic equation $\nabla^4 w = \nabla^2(\nabla^2 w) = q/D$.

Practice

Problem 1: Mean value theorem and harmonic functions

Test if $V(x,y) = e^x \cos(y) + 3xy$ is harmonic by evaluating $\nabla^2 V$.

Solution:

$\partial_{xx} V = e^x \cos y, \quad \partial_{yy} V = -e^x \cos y$ $\nabla^2 V = e^x \cos y - e^x \cos y = 0$ The function is harmonic on $\mathbb{R}^2$. By Earnshaw's theorem, it admits no isolated local maxima or minima.

Problem 2: Poisson's equation and electron cloud density

Electrostatic potential is $V(r) = \frac{e}{4\pi\varepsilon_0} \frac{e^{-2r/a_0}}{r}$. Compute $\nabla^2 V$ for $r > 0$ and find $\rho(r)$.

Solution:

Using $\nabla^2 V = \frac{1}{r}\frac{d^2}{dr^2}(rV)$: $r V = \frac{e}{4\pi\varepsilon_0} e^{-2r/a_0} \implies \frac{d^2(rV)}{dr^2} = \frac{4}{a_0^2} \frac{e}{4\pi\varepsilon_0} e^{-2r/a_0}$ $\nabla^2 V = \frac{e}{\pi\varepsilon_0 a_0^2} \frac{e^{-2r/a_0}}{r} \implies \rho(r) = -\varepsilon_0 \nabla^2 V = -\frac{e}{\pi a_0^2 r} e^{-2r/a_0}$

Problem 3: Steady-state conduction in a spherical shell

Shell of radii $R_1 = 0.10 \text{ m}$ and $R_2 = 0.20 \text{ m}$ has boundary temperatures $T(R_1) = 100^\circ\text{C}$ and $T(R_2) = 20^\circ\text{C}$. Find $T(r)$ and evaluate at midpoint $r = 0.15 \text{ m}$.

Solution:

$\nabla^2 T = 0 \implies T(r) = -A/r + B$ $-10A + B = 100, \quad -5A + B = 20 \implies A = -16, \quad B = -60$ $T(r) = \frac{16}{r} - 60$ $T(0.15) = \frac{16}{0.15} - 60 \approx 46.67^\circ\text{C}$