What are spherical harmonics?

Spherical harmonics Yₗᵐ(θ,φ) = Pₗᵐ(cosθ)e^(jmφ) are the angular solutions to the Laplace equation in spherical coordinates. They form a complete orthogonal basis on the sphere, with ℓ = 0 representing monopoles, ℓ = 1 dipoles, ℓ = 2 quadrupoles, and so on. Any function on a sphere can be expanded in spherical harmonics.

What are Legendre polynomials used for?

Legendre polynomials Pₗ(cosθ) are the polar-angle solutions for axially symmetric problems (no φ-dependence). They appear in gravitational multipole expansions, electrostatic boundary problems, and quantum mechanical orbital angular momentum. Their orthogonality allows boundary conditions to be matched term-by-term.

How do you solve the Laplace equation in spherical coordinates?

Use separation of variables f = R(r)Θ(θ)Φ(φ), producing: trigonometric solutions in φ (integer m), associated Legendre functions in θ (integer ℓ), and power-law solutions rˡ and r^(−ℓ−1) in r. Apply boundary conditions to determine coefficients. For axial symmetry, only m = 0 terms survive.

Laplace Equation in Spherical Coordinates

Quick Answer

The Laplace equation in spherical coordinates is (1/r²)∂/∂r(r²∂f/∂r) + (1/(r²sinθ))∂/∂θ(sinθ ∂f/∂θ) + (1/(r²sin²θ))∂²f/∂φ² = 0. Separation of variables yields radial power functions rⁿ, Legendre polynomials Pₙ(cosθ) in the polar angle, and trigonometric functions in the azimuthal angle. Solutions describe gravitational potentials, electrostatic fields, and quantum orbitals. Explore Laplace equation tools at www.lapcalc.com.

Calculate this instantly on LAPLACE Calculator →

Laplace Equation in Spherical Coordinates: Full Expression

The Laplace equation in spherical coordinates (r,θ,φ) sets the spherical Laplacian to zero: (1/r²)∂/∂r(r²∂f/∂r) + (1/(r²sinθ))∂/∂θ(sinθ ∂f/∂θ) + (1/(r²sin²θ))∂²f/∂φ² = 0. This PDE governs steady-state phenomena with spherical symmetry: gravitational potential outside a planet, electrostatic potential around a charged sphere, and temperature distribution in a spherical shell. The equation's structure naturally decomposes into radial, polar, and azimuthal parts through separation of variables, producing three ODEs each solvable in closed form.

Key Formulas

Separation of Variables: Radial, Polar, and Azimuthal Solutions

Assuming f(r,θ,φ) = R(r)Θ(θ)Φ(φ) and substituting into the Laplacian equation in spherical coordinates produces three separated ODEs. The azimuthal equation Φ″ + m²Φ = 0 has solutions e^(±jmφ) or cos(mφ), sin(mφ) with integer m for single-valuedness. The polar equation becomes the associated Legendre equation with solutions Pₗᵐ(cosθ), requiring integer ℓ ≥ |m|. The radial equation r²R″ + 2rR′ − ℓ(ℓ+1)R = 0 is Euler-Cauchy with solutions rˡ (regular at origin) and r^(−ℓ−1) (singular at origin). The complete solution combines these as spherical harmonics Yₗᵐ(θ,φ).

Compute laplace equation in spherical coordinates Instantly

Get step-by-step solutions with AI-powered explanations. Free for basic computations.

Open Calculator

Spherical Harmonics and Legendre Polynomials

Solutions to the Laplace equation in spherical coordinates are built from spherical harmonics Yₗᵐ(θ,φ) = Pₗᵐ(cosθ)e^(jmφ), where Pₗᵐ are associated Legendre functions. For axially symmetric problems (no φ-dependence, m = 0), solutions simplify to f(r,θ) = Σ(Aₗrˡ + Bₗr^(−ℓ−1))Pₗ(cosθ), where Pₗ are ordinary Legendre polynomials: P₀ = 1, P₁ = cosθ, P₂ = (3cos²θ−1)/2. These polynomials form a complete orthogonal set on [−1,1], allowing any boundary condition to be expanded in a Legendre series with coefficients determined by orthogonality integrals.

Applications: Gravitational and Electrostatic Potentials

The Laplace equation in spherical coordinates governs gravitational and electrostatic potentials in charge-free or mass-free regions. The gravitational potential outside a spherically symmetric mass distribution is Φ = −GM/r, the ℓ = 0 term. Deviations from spherical symmetry add higher multipole terms: ℓ = 1 (dipole), ℓ = 2 (quadrupole), etc. Earth's gravitational field is expanded in spherical harmonics to model geoid variations. In electrostatics, the potential of a conducting sphere in a uniform field involves the ℓ = 1 Legendre polynomial, and the field of a dielectric sphere uses continuity conditions at the boundary to determine coefficients.

Connecting Spherical Laplace Equation to Laplace Transforms

While the Laplace equation ∇²f = 0 and the Laplace transform L{f(t)} are different mathematical objects, they combine powerfully for time-dependent problems in spherical geometry. The wave equation in spherical coordinates ∂²u/∂t² = c²∇²u can be Laplace-transformed in time, reducing to a modified Helmholtz equation (∇² − s²/c²)Ũ = −(s/c²)u(r,0) − (1/c²)∂u/∂t|₀ in the spatial variables. This is solved using modified spherical Bessel functions, and the time-domain solution is recovered by inverse Laplace transform. This combined approach handles transient wave propagation, seismic analysis, and acoustic radiation problems at www.lapcalc.com.

Frequently Asked Questions

It is ∇²f = 0 written in (r,θ,φ): (1/r²)∂/∂r(r²∂f/∂r) + (1/(r²sinθ))∂/∂θ(sinθ∂f/∂θ) + (1/(r²sin²θ))∂²f/∂φ² = 0. It governs potential fields with spherical symmetry, such as gravitational, electrostatic, and steady-state thermal problems.

Master Your Engineering Math

Join thousands of students and engineers using LAPLACE Calculator for instant, step-by-step solutions.

Start Calculating Free →