---
title: Ball Harmonics Expansion Overview
url: https://www.emergentmind.com/topics/ball-harmonics-expansion
type: topic
---

# Ball Harmonics Expansion Overview

Searching arXiv for recent and foundational papers on ball harmonics expansion.
Ball harmonics expansion denotes a family of separable representations for functions, kernels, or fields on a ball in which the angular dependence is expanded in spherical harmonics or related spherical modes, while the radial dependence is carried by basis functions determined by the underlying operator, boundary condition, or reproducing structure. In the Dirichlet-Laplacian setting on the real unit ball, the canonical basis consists of eigenfunctions $\psi_{l,m,n}(r,\theta,\phi)=C_{l,n}j_l(\alpha_{l,n}r)Y_l^m(\theta,\phi)$ with eigenvalues $\lambda_{l,n}=\alpha_{l,n}^2$ [2406.05922]. In broader usage, the same separation principle appears in shellwise spherical-harmonic reconstructions, in vector spherical harmonics, in Sobolev-orthogonal polynomial bases on the ball, and in reproducing-kernel expansions for $M$-harmonic and $H$-harmonic functions [1602.01301].

## 1. Core separation principle

At its most basic, ball harmonics expansion is the statement that a function on a ball can be decomposed into angular modes on the sphere and radial modes on the interval. In the single-particle imaging formulation, this is written as
\[
f(r,\theta,\phi)=\sum_{n,\ell,m} a_{n\ell m}R_{n\ell}(r)Y_{\ell m}(\theta,\phi),
\]
with admissible radial choices such as spherical Bessel $j_\ell$, Zernike radial polynomials, or Fourier–Bessel modes; the same source emphasizes that a shellwise discretization in radius yields a practical ball expansion sampled in the radial coordinate [1602.01301]. On the unit ball in $\mathbb{R}^3$, the Dirichlet Laplacian eigenbasis gives a distinguished version of the expansion that is orthonormal, rotation-steerable, and frequency-ordered [2406.05922]. On the unit disk, the analogous basis is the Fourier–Bessel basis, also referred to as the harmonics on the disk [2207.13674].

The literature does not impose a single universal radial factor. Instead, the radial part is adapted to the operator under study: Dirichlet eigenfunctions for the Laplacian, shellwise samples in diffraction reconstruction, Sobolev orthogonal polynomials for variational PDE discretizations, or hypergeometric radial factors in reproducing-kernel problems. This suggests that “ball harmonics expansion” is best understood as an operator-adapted angular–radial decomposition rather than a single fixed basis.

| Setting | Angular component | Radial or kernel component |
|---|---|---|
| Dirichlet Laplacian on $B\subset\mathbb{R}^3$ | $Y_l^m$ | $j_l(\alpha_{l,n}r)$ |
| Shellwise diffraction intensity | $Y_{\ell m}$ on each shell | discrete shells $k=k_s$ |
| Disk harmonics on $D\subset\mathbb{R}^2$ | $e^{im\theta}$ | $J_m(j_{m,k}r)$ |
| Complex-ball $M$-harmonic kernels | unitary spherical harmonics $H_{p,q}$ | $S_{p,q}(r)$ |

## 2. Spectral bases, orthogonality, and rotation

For the Dirichlet Laplacian on the unit ball, separation of variables yields the radial equation
\[
r^2R''(r)+2rR'(r)+[\lambda r^2-l(l+1)]R(r)=0,
\]
whose regular solutions are spherical Bessel functions $j_l(\sqrt{\lambda}r)$; the Dirichlet boundary condition forces $\sqrt{\lambda}$ to be a positive zero $\alpha_{l,n}$ of $j_l$ [2406.05922]. With
\[
\psi_{l,m,n}(r,\theta,\phi)=C_{l,n}j_l(\alpha_{l,n}r)Y_l^m(\theta,\phi),\qquad C_{l,n}=\frac{\sqrt{2}}{j_{l+1}(\alpha_{l,n})},
\]
the basis is orthonormal in $L^2(B)$, and the corresponding coefficients are obtained by projection against $\overline{\psi}_{l,m,n}$ [2406.05922]. In the more general radius-$R$ formulation, the orthonormal basis of $L^2(B_R)$ is
\[
B_{lmn}(r,\theta,\phi)=R_{ln}(r)Y_l^m(\theta,\phi),\qquad R_{ln}(r)=c_{ln}(R)j_l(\kappa_{ln}r/R),
\]
where $\kappa_{ln}$ is the $n$-th positive root of $j_l$ and $c_{ln}(R)=\sqrt{2}/(R^{3/2}|j_{l+1}(\kappa_{ln})|)$ [2509.01180].

A central structural property is steerability. Under a rotation $R\in SO(3)$, spherical harmonics transform by Wigner $D$-matrices,
\[
Y_{\ell m}(R\Omega)=\sum_{m'=-\ell}^{\ell}D^\ell_{mm'}(R)Y_{\ell m'}(\Omega),
\]
so the ball-harmonic coefficients transform block-diagonally in $\ell$ and the radial index [1602.01301]. In the Dirichlet-Laplacian basis used for subtomogram alignment, this becomes
\[
c'_{lmn}=\sum_{m'=-l}^{l}D^{(l)}_{mm'}(R)c_{lm'n},
\]
which is the basis of exact rotational matching in coefficient space [2509.01180]. The disk analogue is even simpler: rotating an image by angle $\phi$ multiplies each Fourier–Bessel coefficient by $e^{-im\phi}$ [2207.13674].

Vector-valued variants preserve the same separation principle. For divergence-free fields on the unit ball in $\mathbb{R}^3$, the expansion is organized by vector spherical harmonics, and the toroidal mode is singled out by the identity
\[
\nabla\cdot\big(f(r)\mathbf{Y}^{n,3}_\ell\big)=0.
\]
The resulting divergence-free vectorial ball PSWFs are
\[
\Psi^{\alpha,n}_{k,\ell}(x;c)=(x\times\nabla)\psi^{\alpha,n}_{k,\ell}(x;c),
\]
so the scalar basis and the divergence-free vector basis are linked by the solenoidal operator $x\times\nabla$ [2004.03888].

## 3. Kernel-theoretic expansions on real and complex balls

Ball harmonics expansion also appears in reproducing-kernel theory, where the object being expanded is a kernel rather than a generic $L^2$ function. On the unit ball
\[
B^n=\{z\in\mathbb{C}^n:\|z\|<1\},
\]
the invariant Laplacian
\[
\Delta=4(1-\|z\|^2)\sum_{j,k=1}^n(\delta_{jk}-\overline z_j z_k)\frac{\partial^2}{\partial z_j\partial\overline z_k}
\]
defines the class of $M$-harmonic functions, and the restriction representation of $U(n)$ on $L^2(S^{2n-1},d\sigma)$ decomposes multiplicity-free as
\[
L^2(S^{2n-1},d\sigma)=\bigoplus_{p,q\ge0}H_{p,q},
\]
where $H_{p,q}$ consists of restrictions of harmonic polynomials homogeneous of degree $p$ in $z$ and degree $q$ in $\overline z$ [2208.07358]. The associated Poisson–Szegő kernel admits the expansion
\[
P(r\zeta,\eta)=\sum_{p,q\ge0}S_{p,q}(r)H_{p,q}(\langle\zeta,\eta\rangle),
\]
and the $M$-harmonic Szegő kernel satisfies
\[
K_{Sz}(z,w)=\sum_{p,q\ge0}S_{p,q}(|z|)S_{p,q}(|w|)H_{p,q}(\langle\zeta_z,\zeta_w\rangle),
\]
with an explicit closed form in terms of the four-variable hypergeometric function $FD_1$ [2208.07358]. The same work gives a general series for weighted $M$-harmonic Bergman kernels and argues against a simple closed form in higher dimensions by exhibiting coefficients involving $\zeta(3)$ already in the case $(n,s)=(2,0)$ [2208.07358].

A related but distinct real-ball theory concerns $H$-harmonic functions, i.e. functions annihilated by the hyperbolic Laplacian on $B^n\subset\mathbb{R}^n$. There the Szegő kernel is
\[
K_h(x,y)=\int_{\partial B^n}P_h(\eta,x)P_h(\eta,y)\,d\sigma(\eta),
\]
and it admits an explicit series expansion via Exton’s triple hypergeometric function $X_9$, as well as a finite sum representation in Gauss ${}_2F_1$ functions [2510.11821]. In the weighted Bergman setting,
\[
K_h^s(x,y)=\sum_{m=0}^\infty \frac{S_m(|x|^2)S_m(|y|^2)}{I_m(s)}\,Z_m(x,y),
\]
which realizes the kernel as a spherical/radial decomposition into zonal harmonics and radial factors [2510.11821].

These kernel-theoretic examples show that ball harmonics expansion is not restricted to eigenfunction expansions of the Laplacian. It also functions as a representation-theoretic device for decomposing reproducing kernels into spherical projectors and radial coefficients.

## 4. Fast transforms and PDE-adapted polynomial bases

A major recent development is the construction of fast transforms between Cartesian data and ball-harmonic coefficients. For voxelized data on an $N\times N\times N$ grid, fast expansion into the Dirichlet-Laplacian basis on the unit ball is achieved by combining a 3D NUFFT, spherical harmonic transforms on $S^2$, and radial interpolation from Chebyshev nodes to the Dirichlet spectrum $\{\alpha_{l,n}\}$ [2406.05922]. The resulting forward and inverse transforms satisfy relative $\ell^1$–$\ell^\infty$ accuracy $\varepsilon$ and run in
\[
O(N^3(\log N)^2+N^3|\log\varepsilon|^2),
\]
whereas naive dense application of the expansion operators costs $O(N^6)$ [2406.05922]. The 2D analogue, the Fast Disk Harmonics Transform, computes Fourier–Bessel coefficients from an $L\times L$ Cartesian image in $O(L^2\log L)$ operations and preserves the diagonal action of rotations and radial convolutions in coefficient space [2207.13674].

Spectral methods on the ball also motivate polynomial ball-harmonic systems that are not Dirichlet eigenfunctions. For the stationary Schrödinger equation
\[
-\Delta u+\lambda(1-\|x\|^2)^\kappa u=f(x)\quad\text{in }\mathbb{B}^d,\qquad u|_{\mathbb{S}^{d-1}}=0,
\]
a fully diagonalized method is obtained by working in the polynomial space
\[
\mathcal{R}_n^d=\{(1-\|x\|^2)p(x):p\in\Pi_n^d\}
\]
and constructing basis functions
\[
\Phi_{n,j,\nu}(x)=(1-\|x\|^2)\,q_j^{(\beta)}(2\|x\|^2-1)Y_\nu^{n-2j}(x),
\]
where the radial factors are Sobolev orthogonal polynomials associated with a univariate Sobolev inner product [2601.15911]. Because the basis is orthogonal in the variational inner product, the Galerkin matrix becomes diagonal mode-by-mode [2601.15911].

An earlier ball-polynomial construction uses the inner product
\[
\langle f,g\rangle_\mu^\lambda
=
\frac{1}{\omega_\mu}\int_{B^d}f(x)g(x)W_\mu(x)\,dx
+
\frac{\lambda}{\sigma_{d-1}}\int_{S^{d-1}}f(\xi)g(\xi)\,d\sigma(\xi),
\]
which adds a boundary term on the sphere [1511.07056]. The associated basis functions keep the spherical harmonics unchanged and modify only the radial Jacobi factor to a Jacobi-type polynomial
\[
q_k^{(\alpha,\beta,\lambda)}(t)=\Big[a_k-(1+t)\frac{d}{dt}\Big]P_k^{(\alpha,\beta)}(t),
\]
thereby producing orthogonal ball expansions adapted to the ball-plus-sphere geometry [1511.07056]. For $\mu=\tfrac12$, these polynomials satisfy a fourth-order PDE on the ball [1511.07056].

## 5. Imaging, alignment, and field simulation

In inverse problems, ball harmonics expansion is often used to decouple angular and radial degrees of freedom while preserving exact rotation laws. In single-particle imaging with X-ray lasers, the diffraction intensity on each shell is expanded as
\[
I(k,\Omega)=\sum_{\ell=0}^{L_{\max}(k)}\sum_{m=-\ell}^{\ell}c_{\ell m}(k)Y_{\ell m}(\Omega),
\]
with Friedel symmetry forcing odd $\ell$ to vanish, and the radial dependence handled by processing concentric shells $k=k_s$ independently [1602.01301]. The shell-by-shell EMC algorithm alternates expansion, maximization, and compression; for each shell, the dominant M-step cost scales as
\[
\mathcal{O}(L^3N_sK),
\]
and the paper reports, for example, a full reconstruction of the STNV virus with $40$ shells, $\Delta k=0.01\,\text{\AA}^{-1}$, $k_{\max}=0.42\,\text{\AA}^{-1}$, $1000$ patterns, and R-factors $0.152$, $0.143$, and $0.0979$ for truncated theory at $L=3,5,7$ respectively [1602.01301].

For subtomogram alignment in cryo-ET, the Dirichlet-Laplacian basis on the ball is used directly. Rotational correlation between two volumes is expressed as
\[
S(R;s)=\sum_{l=0}^{L}\sum_{n=1}^{N_l}[c^A_{ln}(s)]^\dagger D^{(l)}(R)c^B_{ln},
\]
so rotation evaluation reduces to small matrix-vector contractions in $(l,n)$ blocks [2509.01180]. The reported implementation uses $L_{\max}=42$ and $L_{\text{set}}=\{7,12,33\}$, combines frequency marching with Newton refinement, and achieves at least $10\times$ speed-up compared to FRM-style exhaustive rotational matching while retaining identical accuracy to FRM and often sub-degree accuracy; a representative case achieves approximately $1.5^\circ$ orientation error [2509.01180].

In plasma simulation, the same angular–radial decomposition appears in vector form. The VSHPIC algorithm expands electric and magnetic fields as
\[
\mathbf{E}=\sum_{l=0}^{l_{\max}}\sum_{m=-l}^{l}\Big(E_{lm}^{(r)}\mathbf{\Upsilon}_{lm}+E_{lm}^{(\Psi)}\mathbf{\Psi}_{lm}+E_{lm}^{(\Phi)}\mathbf{\Phi}_{lm}\Big),
\]
and similarly for $\mathbf{B}$ [2402.05446]. Projection of Maxwell’s equations onto the vector spherical harmonics reduces the field solve to a set of 1D radial systems for each retained $(l,m)$ mode, while a staggered deposition scheme enforces the discrete continuity equation and achieves deviations smaller than $10^{-13}$ in double precision over long runs [2402.05446]. The code also reproduces the cold-plasma dispersion relation $\omega^2=\omega_p^2+c^2k^2$ from the transformed $E_{lm}^{(\Phi)}$ modes [2402.05446].

## 6. Special cases, misconceptions, and related expansions

Several special cases clarify what ball harmonics expansion is, and what it is not. On the unit disk, the Dirichlet eigenbasis reduces to the Fourier–Bessel system
\[
\psi_{m,k}(r,\theta)=c_{m,k}J_m(j_{m,k}r)e^{im\theta},
\]
which is orthogonal, frequency-ordered, and steerable; on the complex unit disc, the $M$-harmonic Szegő kernel reduces to the classical harmonic formula
\[
K_{Sz}(z,w)=\frac{1-|z|^2|w|^2}{2\pi|1-z\overline w|^2}
\]
when $n=1$ [2207.13674; 2208.07358]. These low-dimensional cases are much simpler than their higher-dimensional analogues, where multivariable hypergeometric structure and logarithmic-type boundary behavior appear [2208.07358].

A common misunderstanding is to identify ball harmonics expansion exclusively with regular solid harmonics $r^\ell Y_\ell^m$. The cited literature supports a broader view. Regular and irregular solid spherical harmonics are one important family, but Dirichlet eigenfunctions use $j_\ell(\alpha_{\ell n}r)$, vectorial PSWFs use toroidal vector spherical modes, Sobolev methods use univariate Sobolev orthogonal polynomials in $2r^2-1$, and reproducing-kernel theories use hypergeometric radial factors [2406.05922; 2004.03888; 2601.15911; 2208.07358]. A plausible implication is that the term is most coherent when read functionally: it denotes a decomposition on the ball into spherical angular modes and a problem-specific radial sector.

Related coordinate systems further delimit the scope of the method. Relationships between solid spherical and toroidal harmonics allow one to convert toroidal-harmonic representations into regular or irregular spherical expansions, with coefficients determined by forward-stable recurrences for $C_{nk}^m$ and $S_{nk}^m$ [1802.03484]. However, those expansions have restricted convergence domains, and the paper emphasizes that the relationships between toroidal and spherical harmonics are non-invertible in general because singularity geometry determines whether a convergent expansion exists [1802.03484]. This is an important caution: ball harmonics expansion is powerful when its radial basis matches the operator and singularity structure, but it is not basis-independent, nor does it erase domain-of-convergence issues.

Taken together, these developments place ball harmonics expansion at the intersection of harmonic analysis, spectral theory, numerical analysis, and computational imaging. Its unifying mechanism is separation into spherical angular modes and radial structure on the interval, while its diversity comes from the fact that the radial sector may encode Dirichlet spectra, Sobolev orthogonality, reproducing kernels, divergence constraints, or shellwise discretizations, depending on the problem class [2406.05922].

Source: https://www.emergentmind.com/topics/ball-harmonics-expansion