---
title: 'Bispherical Harmonics: Theory and Applications'
url: https://www.emergentmind.com/topics/bispherical-harmonics
type: topic
---

# Bispherical Harmonics: Theory and Applications

to=arxiv_search  天天中彩票和აjson
{"query":"bispherical harmonics arXiv", "max_results": 5}
to=search_arxiv 植物百科通.json
{"query":"bispherical harmonics", "max_results": 10}
Bispherical harmonics are harmonic constructions associated with two-point or two-center geometry, and the literature represented here uses the term in two technically distinct settings. In the angular-momentum formulation of isotropic $N$-point basis functions, the $N=2$ case gives the bipolar, or bispherical, spherical harmonics $\mathcal Y^{LM}_{l_1l_2}$ on $S^2\times S^2$, obtained by Clebsch–Gordan coupling of two ordinary spherical harmonics; in the isotropic sector, only the scalar $L=0$ component remains, and it reduces to a Legendre polynomial in $\hat r_1\!\cdot\!\hat r_2$. In the potential-theoretic setting, bispherical harmonics are the separated solutions of Laplace’s equation in bispherical coordinates, with angular dependence carried by associated Legendre functions and with internal and external branches adapted to two-sphere boundary-value problems. The two usages share the representation theory of $\mathrm{SO}(3)$ and the central role of addition theorems, but they arise from different analytic problems [2010.14418] [2303.02235].

## 1. Angular-momentum formulation and basic definition

Let $Y_{lm}(\hat r)$ denote the standard spherical harmonics in the Condon–Shortley convention, normalized by
$$
\int d\hat{\bf r}\, Y_{lm}(\hat{\bf r})Y_{l'm'}^*(\hat{\bf r})
= \delta^K_{ll'}\,\delta^K_{mm'}.
$$
For a common rotation $R$, they transform through Wigner $D$-matrices as
$$
Y_{lm}(R\,\hat{\bf r})=\sum_{m'} D^l_{mm'}(R)\,Y_{lm'}(\hat{\bf r}).
$$
Successive Clebsch–Gordan coupling then defines coupled tensor products of spherical harmonics, and the $N=2$ specialization yields the bipolar spherical harmonics
$$
\mathcal Y^{LM}_{l_1l_2}(\hat r_1,\hat r_2)
=
\sum_{m_1,m_2}
C^{LM}_{l_1m_1\,l_2m_2}
Y_{l_1m_1}(\hat r_1)\,
Y_{l_2m_2}(\hat r_2).
$$
These objects transform as rank-$L$ irreducible tensors under simultaneous rotation of both arguments, so the rotationally invariant sector is singled out by $L=0$ [2010.14418].

This construction is the $N=2$ instance of the general isotropic basis
$$
\mathcal P_\Lambda(\hat{\bf r}_1,\ldots,\hat{\bf r}_N)
=
\sum_{m_1,\ldots,m_N}
C^\Lambda_M
\prod_{i=1}^N Y_{\Lambda_i m_i}(\hat{\bf r}_i),
$$
where the weight $C^\Lambda_M$ encodes a fixed coupling scheme through primary angular momenta $\Lambda_i$ and intermediate angular momenta such as $\Lambda_{12},\Lambda_{123},\ldots$. Under simultaneous rotation, the coupled product transforms as an irreducible tensor of total angular momentum $L$; isotropy therefore amounts to projection onto the scalar $L=0$ component.

## 2. Bipolar harmonics on $S^2\times S^2$

The bipolar basis $\{\mathcal Y^{LM}_{l_1l_2}\}$ is orthonormal on $S^2\times S^2$:
$$
\int d\hat{\bf r}_1\,d\hat{\bf r}_2\;
\mathcal Y^{LM}_{l_1l_2}(\hat r_1,\hat r_2)\,
\left[\mathcal Y^{L'M'}_{l_1'l_2'}(\hat r_1,\hat r_2)\right]^*
=
\delta^K_{l_1l_1'}\,
\delta^K_{l_2l_2'}\,
\delta^K_{LL'}\,
\delta^K_{MM'}.
$$
It is also complete for square-integrable functions on $S^2\times S^2$:
$$
f(\hat r_1,\hat r_2)
=
\sum_{l_1,l_2}\sum_{L,M}
f^{LM}_{l_1l_2}\,
\mathcal Y^{LM}_{l_1l_2}(\hat r_1,\hat r_2),
$$
with coefficients obtained by projection against $\left[\mathcal Y^{LM}_{l_1l_2}\right]^*$ and the measure $d\hat{\bf r}_1\,d\hat{\bf r}_2$ [2010.14418].

The basis is constrained by standard selection rules. The triangle rule requires
$$
|l_1-l_2|\le L\le l_1+l_2.
$$
Parity follows directly from the parity of ordinary spherical harmonics:
$$
\mathcal Y^{LM}_{l_1l_2}\to (-1)^{l_1+l_2}\mathcal Y^{LM}_{l_1l_2}
$$
under inversion $\hat r_i\to -\hat r_i$. In the conventions used in the isotropic $N$-point basis paper, parity-even combinations are real and parity-odd combinations are purely imaginary. The same parity bookkeeping extends to the general isotropic basis through the factor
$$
S(\Lambda)\equiv (-1)^{\sum_i \Lambda_i},
$$
where the sum runs over primary angular momenta only.

## 3. Scalar $L=0$ sector and reduction to Legendre form

For isotropic two-point functions, only the scalar bipolar harmonic survives:
$$
\mathcal Y^{00}_{l_1l_2}(\hat r_1,\hat r_2)
=
\sum_{m_1,m_2}
C^{00}_{l_1m_1\,l_2m_2}
Y_{l_1m_1}(\hat r_1)\,
Y_{l_2m_2}(\hat r_2).
$$
Coupling to $L=0$ forces $l_1=l_2\equiv l$. Using the special Clebsch–Gordan coefficient for coupling $l$ with $l$ to zero and the identity $Y_{l,-m}=(-1)^mY_{lm}^*$, one obtains
$$
\mathcal Y^{00}_{ll}(\hat r_1,\hat r_2)
=
(-1)^l\frac{\sqrt{2l+1}}{2l+1}
\sum_m Y_{lm}(\hat r_1)Y_{lm}^*(\hat r_2).
$$
The spherical harmonic addition theorem,
$$
P_l(\cos\gamma)
=
\frac{4\pi}{2l+1}
\sum_m Y_{lm}(\hat r_1)Y_{lm}^*(\hat r_2),
\qquad
\cos\gamma=\hat r_1\cdot \hat r_2,
$$
then gives the Legendre form
$$
\mathcal Y^{00}_{ll}(\hat r_1,\hat r_2)
=
\kappa_l\,P_l(\hat r_1\cdot \hat r_2),
\qquad
\kappa_l=(-1)^l\frac{\sqrt{2l+1}}{4\pi}.
$$
Equivalently, the isotropic basis elements can be written as
$$
\mathcal P_{ll}(\hat r_1,\hat r_2)
=
\frac{\sqrt{2l+1}}{4\pi}(-1)^l
P_l(\hat r_1\cdot \hat r_2).
$$
This is the precise sense in which isotropic bispherical harmonics collapse to Legendre polynomials [2010.14418].

The lowest modes make the reduction explicit:
$$
\mathcal Y^{00}_{00}(\hat r_1,\hat r_2)=\frac{1}{4\pi},
$$
$$
\mathcal Y^{00}_{11}(\hat r_1,\hat r_2)
=
-\frac{\sqrt{3}}{4\pi}\,
\hat r_1\!\cdot\!\hat r_2,
$$
$$
\mathcal Y^{00}_{22}(\hat r_1,\hat r_2)
=
\frac{3}{2}\frac{\sqrt{5}}{4\pi}
\left[
(\hat r_1\!\cdot\!\hat r_2)^2-\frac13
\right].
$$
These examples are consistent with $\mathcal Y^{00}_{ll}=\kappa_l P_l(\hat r_1\!\cdot\!\hat r_2)$.

## 4. General isotropic $N$-point basis, recoupling, and symmetry

The $N=2$ case sits inside a broader coupled-harmonic framework in which tripolar and quadrupolar harmonics appear for $N=3$ and $N=4$. For three directions, one coupling scheme is
$$
\bigl[[Y_{l_1}(\hat r_1)\otimes Y_{l_2}(\hat r_2)]_{l_{12}}\otimes Y_{l_3}(\hat r_3)\bigr]_{LM},
$$
and isotropy requires $L=0$, which in this scheme enforces $l_{12}=l_3$. For four directions, canonical schemes couple pairs such as $(l_1,l_2)$ and $(l_3,l_4)$ before final coupling to total $L$. Different coupling schemes are related by Wigner $6$-$j$ and $9$-$j$ symbols, with the standard Racah recoupling identity for three angular momenta reading
$$
([l_1\otimes l_2]_{l_{12}}\otimes l_3)_{LM}
=
\sum_{l_{23}}
(-1)^{l_1+l_2+l_3+L}
\sqrt{(2l_{12}+1)(2l_{23}+1)}
\left\{
\begin{array}{ccc}
l_1 & l_2 & l_{12}\\
l_3 & L   & l_{23}
\end{array}
\right\}
(l_1\otimes [l_2\otimes l_3]_{l_{23}})_{LM}.
$$
For four angular momenta, the analogous recouplings are controlled by $9$-$j$ symbols [2010.14418].

The full isotropic basis is orthonormal in both primary and intermediate labels:
$$
\int \Bigl[\prod_{i=1}^N d\hat{\bf r}_i\Bigr]\;
\mathcal P_\Lambda(\hat{\bf r}_1,\ldots,\hat{\bf r}_N)\,
\mathcal P_{\Lambda'}(\hat{\bf r}_1,\ldots,\hat{\bf r}_N)
=
\delta^K_{\Lambda_1\Lambda_1'}
\delta^K_{\Lambda_2\Lambda_2'}
\cdots
\delta^K_{\Lambda_{12}\Lambda_{12}'}
\delta^K_{\Lambda_{123}\Lambda_{123}'}
\cdots.
$$
Any isotropic function $F_{\rm iso}(\hat{\bf r}_1,\ldots,\hat{\bf r}_N)$ can therefore be expanded in the basis $\mathcal P_\Lambda$, and rotational averaging acts as a projector onto the total-$L=0$ sector. Reordering the arguments $\hat{\bf r}_i$ produces linear combinations of the canonical-order basis, again with coefficients expressible in terms of $6$-$j$ and $9$-$j$ symbols. The paper emphasizes that Yutsis diagrams streamline these manipulations by graphically encoding products of $3$-$j$ symbols and their phases. The chain weight $C^\Lambda_M$ itself can be written as a product of $3$-$j$ symbols with a phase
$$
\kappa=(\Lambda_{12}-m_{12})+(\Lambda_{123}-m_{123})+\cdots,
$$
which makes the algebraic structure of the basis explicit.

## 5. Bispherical coordinates and Laplace-separable harmonics

A second usage of bispherical harmonics arises in the separation of Laplace’s equation in bispherical coordinates. In the formulation summarized from Alexander–Cohl–Volkmer, bispherical coordinates $(\mu,\nu,\phi)$ with focal distance $a>0$ are defined by
$$
x=\frac{a\sinh\mu}{\cosh\mu-\cos\nu}\cos\phi,\qquad
y=\frac{a\sinh\mu}{\cosh\mu-\cos\nu}\sin\phi,\qquad
z=\frac{a\sin\nu}{\cosh\mu-\cos\nu},
$$
with $\mu\in(-\infty,\infty)$, $\nu\in(0,\pi)$, and $\phi\in(-\pi,\pi]$. The scale factors and volume element are
$$
h_\mu=h_\nu=\frac{a}{\cosh\mu-\cos\nu},\qquad
h_\phi=\frac{a\sinh\mu}{\cosh\mu-\cos\nu},
$$
$$
dV=\frac{a^3\sinh\mu}{(\cosh\mu-\cos\nu)^3}\,d\mu\,d\nu\,d\phi.
$$
With the prefactor
$$
U(\mu,\nu,\phi)
=
[\cosh\mu-\cos\nu]^{1/2}M(\mu)N(\nu)e^{im\phi},
$$
the Laplacian becomes separable, leading to the associated Legendre equations
$$
(1-u^2)N''(u)-2uN'(u)+
\left[
\ell(\ell+1)-\frac{m^2}{1-u^2}
\right]N(u)=0,
\qquad
u=\cos\nu,
$$
and
$$
(x^2-1)M''(x)+2xM'(x)-
\left[
\ell(\ell+1)-\frac{m^2}{x^2-1}
\right]M(x)=0,
\qquad
x=\cosh\mu.
$$
The separated solutions are
$$
N(\nu)=P_\ell^m(\cos\nu),\qquad
M(\mu)\in\{P_\ell^m(\cosh\mu),\,Q_\ell^m(\cosh\mu)\},
$$
with integer $\ell\ge |m|$ [2303.02235].

Internal and external bispherical harmonics are then represented as
$$
\Phi^{(\mathrm{int})}_{\ell m}(\mu,\nu,\phi)
=
N^{(\mathrm{int})}_{\ell m}
[\cosh\mu-\cos\nu]^{1/2}
P_\ell^m(\cos\nu)\,
P_\ell^m(\cosh\mu)\,
e^{im\phi},
$$
$$
\Phi^{(\mathrm{ext})}_{\ell m}(\mu,\nu,\phi)
=
N^{(\mathrm{ext})}_{\ell m}
[\cosh\mu-\cos\nu]^{1/2}
P_\ell^m(\cos\nu)\,
Q_\ell^m(\cosh\mu)\,
e^{im\phi}.
$$
Here $P_\ell^m(\cos\nu)$ is the Ferrers function on $(-1,1)$, $P_\ell^m(\cosh\mu)$ grows with $\mu$, and $Q_\ell^m(\cosh\mu)$ decays as $\mu\to +\infty$, so the latter is appropriate for external solutions. The angular factor obeys the standard orthogonality relation
$$
\int_0^\pi
P_\ell^m(\cos\nu)\,
P_{\ell'}^m(\cos\nu)\,
\sin\nu\,d\nu
=
\frac{2}{2\ell+1}
\frac{(\ell+m)!}{(\ell-m)!}
\delta_{\ell\ell'},
$$
and a convenient angular normalization is
$$
N^{(\mathrm{ang})}_{\ell m}
=
\left[
\frac{2\ell+1}{4\pi}
\frac{(\ell-m)!}{(\ell+m)!}
\right]^{1/2}.
$$

## 6. Green’s functions, limiting constructions, and applications

The potential-theoretic significance of bispherical harmonics is encoded in the Green’s expansion
$$
\frac{1}{\|x-x'\|}
=
[\cosh\mu-\cos\nu]^{1/2}
[\cosh\mu'-\cos\nu']^{1/2}
\sum_{\ell=0}^\infty
\sum_{m=-\ell}^{\ell}
\frac{(\ell-m)!}{(\ell+m)!}
e^{-(\ell+1/2)(\mu_>-\mu_<)}
P_\ell^m(\cos\nu)\,
P_\ell^m(\cos\nu')\,
e^{im(\phi-\phi')},
$$
where $\mu_<=\min\{\mu,\mu'\}$ and $\mu_>=\max\{\mu,\mu'\}$. Equivalently,
$$
G(x,x')
=
\frac{1}{4\pi\|x-x'\|}
=
\sum_{\ell=0}^\infty
\sum_{m=-\ell}^{\ell}
C_{\ell m}\,
\Phi^{(\mathrm{int})}_{\ell m}(\mu_<,\nu_<,\phi_<)\,
\Phi^{(\mathrm{ext})}_{\ell m}(\mu_>,\nu_>,\phi_>),
$$
with
$$
C_{\ell m}
=
\frac{1}{4\pi}\frac{(\ell-m)!}{(\ell+m)!}
$$
for the unnormalized choice $N^{(\mathrm{int})}_{\ell m}=N^{(\mathrm{ext})}_{\ell m}=1$. Only modes with the same $m$ couple when boundary data are matched; for real-valued potentials one may use $\cos(m(\phi-\phi'))$ forms or pair the $m$ and $-m$ coefficients [2303.02235].

Alexander–Cohl–Volkmer derive this bispherical expansion as the $k\to 0$ limit of the bi-cyclide harmonic expansion. In that limit, the Lamé–Wangerin functions reduce to associated Legendre functions, and the bi-cyclide Green’s series reduces term-by-term to the classical bispherical series with $\ell=m+n$. The paper also emphasizes the addition theorem obtained from the bi-cyclide expansion and its reduction to the bispherical addition theorem. In this sense, bispherical harmonics sit at the intersection of special-function theory, separation of variables, and two-center potential theory.

The two principal application domains reflected in the cited literature are correspondingly different. In cosmology and related isotropic statistics, the $N=2$ scalar bispherical harmonics reduce angular dependence to Legendre polynomials of $\hat r_1\!\cdot\!\hat r_2$, while the $N=3$ and higher isotropic bases provide orthonormal expansions for angular structure in higher-point correlation functions and allow products and permutations to be reduced with $6$-$j$ and $9$-$j$ symbols [2010.14418]. In mathematical physics, bispherical-coordinate harmonics are tailored to the electrostatics of two coaxial spheres, capacitance and multipole interactions, and hydrodynamic or Stokes-flow problems around rigid spheres; internal modes enforce regularity, external modes enforce decay, and the Kelvin transform interchanges the two types [2303.02235].

Source: https://www.emergentmind.com/topics/bispherical-harmonics