---
title: Chu–Harrington Limit (CHL) Explained
url: https://www.emergentmind.com/topics/chu-harrington-limit-chl
type: topic
---

# Chu–Harrington Limit (CHL) Explained

Searching arXiv for recent and directly relevant papers on the Chu–Harrington limit and related updates.
The Chu–Harrington limit (CHL) denotes a family of size-constrained bounds on the radiation quality factor \(Q\) of passive radiators whose sources are confined within a smallest enclosing sphere of radius \(a\), and, by extension, on the fractional bandwidth of electrically small single-resonant devices. In the classical spherical formulation, the exterior field is expanded in TE/TM spherical modes with free-space wavenumber \(k=\omega/c\), and the bound is expressed in terms of the electrical size \(ka\). For the dominant dipolar modes and \(ka \ll 1\), the canonical forms are \(Q_{\min}^{\mathrm{TM1}} \gtrsim (ka)^{-3} + (ka)^{-1}\) and \(Q_{\min}^{\mathrm{TE1}} \gtrsim (ka)^{-3} + (ka)^{-1}\). Recent work shows that this classical interpretation remains accurate in the electrically small regime, while later analyses refine its meaning for two-dimensional radiators, moderate \(ka\), directivity-constrained radiation, singular plasmonic resonators, and naturally resonant emitters [2408.07085] [2107.10770] [2604.26980].

## 1. Classical definition and mathematical setting

Chu’s historical construction treats each spherical TE/TM mode of a source enclosed within \(r \le a\) by means of an exact external ladder circuit, and defines modal \(Q\) as the energy stored in reactive elements divided by radiated power. Harrington later extended the framework to directivity–bandwidth tradeoffs for bounded radiators. In the notation used in recent re-examinations, \(k=\omega/c\), \(a\) is the radius of the smallest sphere enclosing the radiator, \(n\) is the spherical harmonic order, and \(D\) is the directivity. Depending on convention, \(Q\) is written either as \(Q = 2\omega \max(W_e,W_m)/P\) in a fields-based form or as \(Q = \omega W_s/P\) in a single-resonance convention [2408.07085].

The standard CHL assumptions are specific. The radiator is linear, passive, time-invariant, and typically reciprocal; all sources and materials are contained within a sphere of radius \(a\); the exterior region is homogeneous, isotropic, and lossless free space; and the fields outside the sphere are represented by outgoing spherical waves. Under those assumptions, the lower bound applies to the radiation-limited quality factor \(Q_{\text{rad}}\), not directly to input impedance, matching-network complexity, or peak cross-section. Collin and Rothschild re-derived the classical values via field integrals that subtract radiative energy, so the circuit-based and field-based formulations coincide in the electrically small regime [2408.07085].

## 2. Canonical small-antenna formulas and Harrington’s directivity criterion

For the lowest-order spherical modes, the widely cited electrically small asymptotics are
\[
Q^{\mathrm{Chu}}_{\mathrm{TM10}} \approx \frac{1}{(ka)^3} + \frac{1}{ka},
\qquad
Q^{\mathrm{Chu}}_{\mathrm{TE10}} \approx \frac{1}{(ka)^3} + \frac{1}{ka}.
\]
For a combined self-resonant TE+TM dipole with equal powers and matched phases, the best order-wise scaling is commonly quoted as \(Q \approx 1/(ka)^3\). In single-resonant high-\(Q\) cases, the bandwidth proxy \(B_{3\text{dB}} \approx 2/Q\) is used, whereas for multi-resonant behavior or \(Q \lesssim 1\) that relation is not precise [2408.07085].

Harrington’s contribution addressed directivity under a bandwidth constraint. With the first \(N\) spherical harmonics, using both TE and TM families phased constructively, the maximum “normal directivity” is
\[
D_{\max}(N)=N^2+2N.
\]
With the identification \(N \approx ka\), this becomes
\[
D_{\mathrm{normal}}(ka) = (ka)^2 + 2ka.
\]
Within Harrington’s assumptions, \(D > (ka)^2 + 2ka\) implies “high \(Q\).” The classical picture therefore couples enclosure size, modal content, and attainable bandwidth: strong directivity growth requires access to higher spherical orders, but those higher orders are accompanied by larger reactive-energy burdens in the electrically small regime [2408.07085].

## 3. Two-dimensional CHL-type bounds for translationally invariant radiators

For translationally invariant objects, the 2021 analysis of singular plasmonic resonators derives two-dimensional Chu-type bounds directly rather than importing the three-dimensional spherical result. In the quasi-static limit for electric-kind radiators, with longitudinal currents and zero curl, the radiation \(Q\) can be written as the functional
\[
\mathcal{Q}
=
-\frac{8}{\pi}
\frac{\displaystyle \oint_{\partial} \sigma \oint_{\partial} \sigma \log \Delta r \, dl \, dl'}
{\displaystyle \oint_{\partial} \sigma \oint_{\partial} \sigma (\Delta r)^2 \, dl \, dl'}
\frac{1}{x^2},
\]
where \(\partial\) is the cross-section boundary, \(\sigma=\mathbf j \cdot \hat{\mathbf n}\) is the surface charge density per unit length, \(\Delta r = |r-r'|/\ell_c\), and \(x=(\omega/c_0)\ell_c\) with \(\ell_c\) the radius of the smallest enclosing circle [2107.10770].

Minimizing the quasi-static functional subject to fixed dipole moments yields geometry-dependent lower bounds in terms of polarizabilities. For the electric kind,
\[
(x^2Q)_{\min} = \frac{8}{\gamma_{\max,e}},
\]
where \(\gamma_{\max,e}\) is the largest eigenvalue of the scaled two-dimensional electric polarizability tensor. For the magnetic kind,
\[
(x^2Q)_{\min} = \frac{4}{\gamma_m},
\]
where \(\gamma_m\) is the scalar two-dimensional magnetic polarizability. For a circular cylinder, the electric and magnetic bounds coincide:
\[
Q_{\min} = \frac{4}{\pi}\frac{1}{x^2}.
\]
For a dimer of two circular cylinders in the touching limit, the minimum normalized value becomes
\[
(x^2Q)_{\min} \to \frac{48}{\pi^3},
\qquad
Q_{\min} = \frac{48}{\pi^3}\frac{1}{x^2}.
\]
These expressions make explicit that, in two dimensions, the minimum \(Q\) scales as \(x^{-2}\) rather than \((ka)^{-3}\). The paper attributes that difference to the logarithmic Green’s function and the altered stored-energy–radiation balance in two-dimensional electrodynamics [2107.10770].

## 4. Singular plasmonic dimers and the loss of the simple \(1/Q\)-bandwidth correspondence

The singular-plasmonics case study is a translationally invariant dimer of two identical circular cylinders of radius \(R\), separated by an edge-edge gap \(\delta\), and described by a Drude permittivity
\[
\varepsilon_R(\omega)=1-\frac{\omega_p^2}{\omega(\omega+i\nu)}.
\]
Under uniform plane-wave excitation polarized along the dimer axis, the electrostatic modal expansion is organized by eigenpermittivities
\[
\varepsilon_k^+ = -\coth(k\mu),
\qquad
\varepsilon_k^- = -\tanh(k\mu),
\qquad
\varepsilon_k^- = \frac{1}{\varepsilon_k^+},
\]
with
\[
\mu = \operatorname{arccosh}\!\left(1+\frac{\delta}{2R}\right)\approx \sqrt{\frac{\delta}{R}}
\quad \text{for } \delta/R \ll 1.
\]
All eigenvalues accumulate at \(\varepsilon_{\text{acc}}=-1\). Under the Drude map, the bright-mode resonance frequencies are
\[
\omega_k^+ = \frac{\omega_p}{\sqrt{1-\varepsilon_k^+}},
\qquad
\omega_{\text{acc}}=\frac{\omega_p}{\sqrt{2}}.
\]
As the gap shrinks, the bright resonances redshift and cluster ever more densely near \(\omega_{\text{acc}}\) [2107.10770].

The \(k\)-th bright mode has normalized dipole moment
\[
\mathbf P_k
=
\sqrt{
\frac{8\pi(\delta/R)}{\delta/R+4}\,k
}\,
e^{-k\mu}\,\hat{\mathbf x},
\]
and radiation quality factor
\[
Q_k^r
=
\frac{8}{|\varepsilon_k-1|}\frac{1}{|\mathbf P_k|^2}\frac{1}{x^2}
=
\frac{4+\delta/R}{2\pi k(\delta/R)}
\left(e^{2k\mu}-1\right)\frac{1}{x^2}.
\]
In the touching limit, the dipole moment shrinks and
\[
Q_k^r \approx \frac{4}{\pi}\left(\frac{\delta}{R}\right)^{-1/2}\frac{1}{x^2},
\]
so the radiation \(Q\) diverges rather than decreases. With Drude loss, the dissipation quality factor is
\[
Q_k^d=\frac{\omega_k^+}{\nu},
\qquad
\frac{1}{Q_k}=\frac{1}{Q_k^r}+\frac{1}{Q_k^d}.
\]
The associated critical-coupling locus is
\[
\mathcal R_k(\nu)
=
\sqrt{
\frac{2c_0^2}{\pi(\delta/R)(4+\delta/R)}
\frac{e^{2k\mu}-1}{k}
\frac{\nu}{(\omega_k^+)^3}
}.
\]
For fixed loss, increasing \(k\) or decreasing \(\delta/R\) drives modes from radiation-limited to dissipation-limited behavior [2107.10770].

The central result is that there is no contradiction with CHL. When peaks are isolated and under-damped, the fractional bandwidth satisfies \(\mathrm{FBW}\approx 1/Q_{\text{tot}}\), and for low-order well-separated peaks the inverse FWHM agrees with the analytic \(Q_k^r\). Near touching, however, the discrete spectrum turns into a dense continuum around the accumulation point. The broadband response over any finite interval then results from the superposition of many plasmon resonances, each individually obeying the two-dimensional Chu bound and in fact possessing large \(Q_{\text{rad}}\). A common misconception is that the broadband scattering of singular dimers implies an anomalously low single-mode \(Q\); the paper shows instead that the link between overall bandwidth and any one resonance’s \(Q\) fails once the peaks cease to be well separated [2107.10770].

## 5. Updated spherical theory beyond the classical Chu model

A 2024 reformulation reexamines spherical radiation through an equivalent one-dimensional transmission-line model that exactly models the fields. For each outward TM\(_{n,m}\) mode, the radial problem is represented by a lossless tapered waveguide with propagation constant and characteristic impedance
\[
\beta_n(r)=k\sqrt{1-\frac{n(n+1)}{(kr)^2}},
\qquad
Z^{\mathrm{TM}}(r)=\frac{\eta_0\beta_n(r)}{k}.
\]
This immediately yields a precise cutoff condition,
\[
kr=\sqrt{n(n+1)}.
\]
Below that radius, \(\beta_n\) is imaginary and the mode is evanescent; above it, \(\beta_n\) is real and the mode propagates. The same analysis introduces a position-dependent reflection coefficient \(\Gamma(r)\) from the outward wave impedance and the local characteristic impedance, and uses it to partition energy in each spherical shell into radiated and stored components [2408.07085].

The new energy partition is
\[
U_r(r)r^2d\Omega
=
\frac{P^+(r)-P^-(r)}{P^+(r)+P^-(r)}\,U(r)r^2d\Omega,
\qquad
U_s(r)r^2d\Omega
=
\frac{2|\Gamma(r)|^2}{1+|\Gamma(r)|^2}\,U(r)r^2d\Omega.
\]
Below cutoff, \(Z^{\mathrm{TM}}\) is purely imaginary, \(|\Gamma|=1\), \(U_r=0\), and all energy is stored. This leads to modal \(Q_n\) values that agree with Chu and Collin in the small-\(ka\) regime but deviate as \(ka\) approaches or exceeds the mode cutoff. For \(n=1\), the updated small-\(ka\) expansion matches Chu to the first two leading orders, yet at \(ka=1\) the new \(Q\) is approximately \(32\%\) smaller than Chu’s \(n=1\) value [2408.07085].

The same framework revises Harrington’s directivity analysis. For equal TE/TM pairs of order \(n\), with amplitudes \(a_n\) phased to add in the main beam,
\[
D(\{|a_n|\})
=
\frac{\left(\sum_n |a_n|\right)^2}
{\sum_n \frac{|a_n|^2}{2n+1}},
\qquad
Q(\{|a_n|\})
=
\frac{\sum_n \frac{|a_n|^2}{2n+1}Q_n}
{\sum_n \frac{|a_n|^2}{2n+1}}.
\]
Lagrange-multiplier optimization at fixed \(Q\) yields
\[
|a_n| \propto \frac{2n+1}{(n+1)(Q_n+\lambda)}.
\]
Unlike Harrington’s truncation at \(N\le ka\), this optimum excites all spherical harmonics, including \(n>ka\), with smoothly decaying weights. Quantitatively, for \(Q=1\) the maximum directivity is within a factor \(\approx 1.4\) of Harrington’s \(D_{\mathrm{normal}}=(ka)^2+2ka\) over \(ka>1\); for \(Q=3\) at \(ka=6.75\), \(D\approx 19.7\) dB, corresponding to \(>200\%\) aperture efficiency; and for large \(ka\) with \(Q=20\), the maximum aperture efficiency remains above unity under the updated \(Q_n\), whereas Chu-based modal \(Q_n\) would predict a decay toward zero [2408.07085].

## 6. Scattering-channel bounds, superscattering, and what CHL does not bound

CHL is frequently conflated with peak-scattering bounds, but the two are distinct. For three-dimensional free-space scattering by a passive object, the optical theorem and S-matrix unitarity impose a per-partial-wave limit
\[
\sigma_{\mathrm{sc},\max}^{(l)}
=
\frac{2l+1}{2\pi}\lambda^2
=
\frac{(2l+1)2\pi}{k^2}.
\]
For the dipolar channel \(l=1\),
\[
\sigma_{\mathrm{sc},\max}^{(l=1)}
=
\frac{3\lambda^2}{2\pi}
=
\frac{6\pi}{k^2}.
\]
This is a cap on peak scattering in one multipolar channel at one frequency. By contrast, CHL constrains minimum radiation \(Q\) and therefore bandwidth. A radiator can exceed the single-channel dipolar peak by combining multiple multipolar channels constructively and still fully satisfy CHL [2209.12627].

The flat split-ring superscatterer demonstrates that distinction explicitly. The reported device is a two-dimensional array of near-field coupled split-ring resonators optimized over a 19-dimensional search space. A genetically optimized six-element configuration, operated in the \(5\)–\(6\) GHz range, exhibits a measured peak total scattering cross-section of about \(40.7\ \mathrm{cm}^2\) at the optimized orientation, exceeding the single-channel dipolar bound at the measured frequency by approximately \(2.17\times\). The enhancement is produced by engineered overlap of multiple resonances, including MD, ED, EQ, MQ, and octupolar terms. The paper proposes an empirical “super-radiant criterion”
\[
\sigma_{\max}^{\mathrm{array}}
\lesssim
N\,\sigma_{\mathrm{sc},\max}^{(l=1)}
=
N\,\frac{3\lambda^2}{2\pi},
\]
together with a practical device-specific form,
\[
\frac{\sigma_{\max}^{\mathrm{array}}}{N\,\sigma_{\max}^{\mathrm{single\ CSRR}}}\lesssim 1,
\]
and reports Monte Carlo evidence that randomly distributed resonant dipoles do not exceed that normalized ceiling [2209.12627].

The result is consistent with CHL for two separate reasons. First, CHL does not constrain peak \(\sigma\) at a single frequency. Second, the aggregate array is not electrically small in Chu’s sense: with a practical circumscribing radius \(a \approx 28\)–\(30\) mm at \(f\approx 5\) GHz, \(ka \approx 2.7\)–\(3.1\), so the classical lower bound \(Q_{\min} \approx 1/(ka)^3 + 1/(ka)\) is weak. Even for a single constituent CSRR, with \(a \approx 3.5\) mm and \(ka \approx 0.38\) near \(5.2\) GHz, the simulated \(Q \approx 61\) remains comfortably above the Chu lower bound of approximately \(20\)–\(21\). A common misconception is therefore that surpassing a single-channel scattering limit constitutes a CHL violation; the data instead show that scattering-channel bounds and CHL regulate different observables [2209.12627].

## 7. Naturally resonant emitters, efficiency bounds, and broader scope

The CHL also extends to naturally resonant electrically small emitters under its standard assumptions. In that setting, \(Q\) is written as
\[
Q=\omega\frac{W_e+W_m}{P_{\mathrm{diss}}},
\qquad
P_{\mathrm{diss}}=P_{\mathrm{rad}}+P_{\mathrm{loss}},
\]
and near a single Lorentzian resonance the bandwidth quality factor satisfies
\[
Q_{\mathrm{bw}}=Q=\frac{f}{\Delta f}=\frac{\omega}{\Delta\omega}.
\]
For dipolar radiators,
\[
Q_{\mathrm{CHL}}(ka)=\frac{1}{(ka)^3}+\frac{1}{ka}+O(ka),
\]
with the deep-subwavelength approximation \(Q_{\mathrm{CHL}}\approx (ka)^{-3}\). Since \(Q_{\mathrm{rad}} \ge Q_{\mathrm{CHL}}\) and \(Q_{\mathrm{total}}=\eta Q_{\mathrm{rad}}\) for radiation efficiency \(\eta=P_{\mathrm{rad}}/(P_{\mathrm{rad}}+P_{\mathrm{loss}})\), one obtains the efficiency bound
\[
\eta \le \min\!\left(1,\frac{f/\Delta f}{Q_{\mathrm{CHL}}}\right).
\]
In the same limit, the paper rewrites this as a bound on radiated power density per unit input power and unit volume,
\[
\frac{(P_{\mathrm{rad}}/P_{\mathrm{in}})}{V}
\le
\frac{6\pi^2 f^4}{c^3\Delta f},
\]
with \(V=(4\pi/3)a^3\) [2604.26980].

This motivates a figure of merit
\[
\mathrm{FOM}
=
\frac{(P_{\mathrm{rad}}/P_{\mathrm{in}})\,c^3\Delta f}{6\pi^2 f^4 V},
\]
for which \(\mathrm{FOM}=1\) is the CHL-attainable ideal of a passive, single-resonance, lossless emitter. The same framework yields CHL-dictated constraints on atomic emitters. Using the transform-limited inequality
\[
Q_{\mathrm{bw}} \le 4\pi f\tau,
\]
with excited-state lifetime \(\tau\), CHL implies
\[
\tau \ge \frac{Q_{\mathrm{CHL}}(ka)}{4\pi f}.
\]
With the spontaneous-emission relation
\[
\Gamma_{\mathrm{rad}}=\frac{\omega^3|\mu|^2}{3\pi\varepsilon_0\hbar c^3},
\]
one obtains
\[
\Gamma_{\mathrm{rad}} \le \frac{4\pi f}{Q_{\mathrm{CHL}}(ka)},
\qquad
|\mu|
\le
\sqrt{\frac{6\pi^2\varepsilon_0\hbar c^3}{\omega^2 Q_{\mathrm{CHL}}(ka)}}.
\]
These results rest on the observation that natural resonance mechanisms change the internal realization of the source but not the exterior electromagnetic boundary-value problem on which CHL is based [2604.26980].

The scope of the limit is therefore broad but not unlimited. It remains applicable to passive, linear, time-invariant free-space emitters bounded by a sphere, regardless of whether resonance is achieved by conventional matching, plasmonic localization, mechanical strain, or quantum transitions. It ceases to apply in its standard form when the assumptions themselves are altered: increasing the effective bounding sphere through arrays or distributed apertures, introducing active or non-Foster circuitry, using parametric or nonlinear time variance, or coupling strongly to non-free-space propagating channels. Across these cases, the consistent theme is that CHL fixes a minimum radiation \(Q\) for the allowable external modal content; apparent exceptions arise when the relevant observable is not \(Q\), when many resonances overlap, or when the enclosure and channel assumptions are changed rather than violated [2604.26980].

Source: https://www.emergentmind.com/topics/chu-harrington-limit-chl