Solutions of Laplces's equation are also called harmonic. One special case are the spherical harmonics. In 3D space they are simply known as spherical harmonics and are closely related to the Legendre Polynomials, which are the solution with "magnetic quantum number", [itex]m=0[/itex]. The other functions are called the associated Legendre functions. With the factor [itex]\exp(\mathrm{i} \varphi)[/itex] and normalized this gives the spherical harmonics, [itex]\mathrm{Y}_{lm}(\vartheta,\varphi)[/itex] which are a complete set on the Hilbert space, [itex]\mathrm{L}^2(\mathrm{S}_1)[/itex], i.e., all functions that are defined on the unit sphere in [itex]\mathbb{R}^3[/itex] that are square integrable. For the spherical harmonics one has
[tex]\int_{\mathrm{S}_1} \mathrm{d} \Omega \mathrm{Y}_{lm}^*(\vartheta,\varphi)\mathrm{Y}_{l'm'}(\vartheta,\varphi)=\int_0^{\pi} \mathrm{d} \vartheta \int_0^{2 \pi} \mathrm{d} \varphi \sin \vartheta \mathrm{Y}_{lm}^*(\vartheta,\varphi)\mathrm{Y}_{l'm'}(\vartheta,\varphi)=\delta_{ll'} \delta_{mm'}.[/tex]
Any harmonic function on [itex]\mathbb{R}^3[/itex] then can be expanded in the sense of convergence with respect to the corresponding Hilbert-space norm by
[tex]f(\vec{x})=\sum_{l=0}^{\infty} \sum_{m=-l}^{l} \left [f_{1lm}(r) r^l + f_{2lm} \frac{1}{r^{l+1}} \right ]\mathrm{Y}_{lm}.[/tex]
In arbitrary dimensions the analogues of the spherical harmonics are known as Gegenbauer functions.