Spherical Hankel Functions

Many fundamental physical processes, from the decay of a subatomic particle to the transmission of a radio signal, involve waves radiating outwards in three-dimensional space. While the basic wave equation in spherical coordinates is well-known, its simplest solutions—spherical Bessel functions—describe standing waves, which don't carry energy away from a source. This creates a gap in our mathematical toolkit: how do we rigorously describe a purely outgoing or incoming traveling wave? The answer lies in a special class of functions known as spherical Hankel functions, which are ingeniously constructed to solve this very problem.

This article provides a comprehensive overview of these vital mathematical tools. In the "Principles and Mechanisms" chapter, we will delve into their origins, exploring how they are built from more basic functions and proving their physical meaning as representations of traveling waves. Following this, the "Applications and Interdisciplinary Connections" chapter will showcase their remarkable versatility, demonstrating how spherical Hankel functions form the backbone of theories in acoustics, quantum mechanics, scattering theory, and more.

Imagine you toss a pebble into a vast, still pond. Ripples spread out in perfect circles, carrying energy away from the point of impact. Now, picture this not on a flat surface, but in three-dimensional space. A tiny light bulb flashes, a speaker emits a "pop," or a subatomic particle decays. In each case, a wave radiates outwards in a spherical shell. How do we describe this fundamental process mathematically? The answer lies in a beautiful and surprisingly practical family of functions known as ​​spherical Hankel functions​​.

This chapter is a journey to understand these functions—not as abstract mathematical curiosities, but as the natural language for describing waves that travel. We will see how they are built, what they mean physically, and why they are indispensable in fields from quantum mechanics to antenna design.

Let's return to our pond analogy. The outgoing ripples are ​​traveling waves​​; they have a clear direction of motion. But you can also create waves in a bathtub by sloshing the water back and forth. The water level rises and falls in a fixed pattern, but the wave itself doesn't "go" anywhere. This is a ​​standing wave​​.

In physics, when we solve the fundamental wave equation (like the Schrödinger equation for a free particle or the Helmholtz equation for electromagnetic waves) in spherical coordinates, we arrive at a crucial equation for the radial part of the wave, known as the ​​spherical Bessel differential equation​​. Just as the equation for a simple harmonic oscillator gives us sines and cosines, this equation gives us a set of solutions.

The most straightforward solutions are the ​​spherical Bessel functions of the first kind​​, denoted $j_l(kr)$, where $k$ is the wavenumber (related to the wave's energy or frequency) and $r$ is the distance from the origin. These functions are well-behaved everywhere, including at the center ($r=0$). But what kind of wave do they describe?

Let's look at the simplest case, the spherically symmetric wave with zero angular momentum ($l=0$). Here, $j_0(kr) = \frac{\sin(kr)}{kr}$. A sine wave is the quintessential standing wave, the result of two identical waves traveling in opposite directions. And indeed, this is precisely what a spherical Bessel function is. It represents a spherical standing wave, with probability or energy sloshing inwards and outwards in perfect balance, resulting in no net flow in either direction.

This is a beautiful and important solution, but it doesn't describe our outgoing ripple from the pebble. For that, we need to dig deeper.

How can we get a traveling wave from a standing wave? The answer is as simple as it is brilliant: a standing wave is just a superposition of two traveling waves. A sine function, for example, can be written using Euler's formula as $\sin(x) = \frac{e^{ix} - e^{-ix}}{2i}$. The term $e^{ix}$ represents a wave traveling in one direction, and $e^{-ix}$ represents a wave traveling in the opposite direction. Our goal is to isolate one of these components.

To do this, we need a second, independent solution to the spherical Bessel equation. This partner function is the ​​spherical Neumann function​​, $n_l(kr)$, which, unlike $j_l(kr)$, is singular at the origin (it blows up to infinity). By itself, it often seems unphysical. But its true purpose is revealed when we combine it with $j_l(kr)$ using the magic of complex numbers.

We define two new functions, the ​​spherical Hankel functions of the first and second kind​​:

$h_l^{(1)}(kr) = j_l(kr) + i n_l(kr)$ $h_l^{(2)}(kr) = j_l(kr) - i n_l(kr)$

This might look like a purely formal mathematical trick. But let's see what happens for the $l=0$ case. We have $j_0(x) = \frac{\sin(x)}{x}$ and $n_0(x) = -\frac{\cos(x)}{x}$. Let's build $h_0^{(1)}(x)$ where $x=kr$:

$h_0^{(1)}(x) = \frac{\sin(x)}{x} + i \left(-\frac{\cos(x)}{x}\right) = \frac{\sin(x) - i\cos(x)}{x}$

Using Euler's identity, we can see that $\sin(x) - i\cos(x) = -i(\cos(x) + i\sin(x)) = -i e^{ix}$. So, the result is astonishingly clean:

$h_0^{(1)}(x) = -i \frac{e^{ix}}{x}$

We have successfully isolated one of the traveling wave components! If we re-attach the standard time dependence $e^{-i\omega t}$, the full wave becomes proportional to $\frac{1}{r}e^{i(kr - \omega t)}$. This is the classic signature of a pure ​​outgoing​​ spherical wave. The amplitude decreases as $1/r$ (so the total energy passing through any spherical shell is constant), and the phase $kr-\omega t$ continuously increases for an observer moving outward with the wave.

By the same token, $h_0^{(2)}(x)$ simplifies to $i \frac{e^{-ix}}{x}$, which represents a pure ​​incoming​​ spherical wave converging on the origin.

This insight allows us to reinterpret the standing wave $j_0(x)$ in a new light. A little algebra shows that:

$j_0(x) = \frac{1}{2}h_0^{(1)}(x) + \frac{1}{2}h_0^{(2)}(x)$

This equation beautifully confirms our intuition: the standing wave $j_0(x)$ is an equal superposition of a purely outgoing wave and a purely incoming wave. The Hankel functions have successfully "unmixed" the two.

This story about "incoming" and "outgoing" waves is compelling, but can we prove it? In quantum mechanics, there is a rigorous way to determine the direction of flow of a particle's probability: the ​​probability current​​, $\vec{j}$. A positive radial current means there is a net flow of probability out of a region, while a negative current signifies a net flow in.

Let's consider a quantum particle whose state is described by a wavefunction proportional to a Hankel function, say $\psi(\vec{r}) \propto h_l^{(1)}(kr)$. If we calculate the total probability flux passing through a very large sphere surrounding the origin, the result is unambiguous. The net flux is strictly positive. This isn't just a mathematical convention; it is a direct physical consequence of the function's structure. The function $h_l^{(1)}(kr)$ describes a particle definitively moving away from the origin. Conversely, a state described by $h_l^{(2)}(kr)$ would yield a strictly negative flux, representing a particle moving toward the origin. This provides the ultimate physical justification for labeling these functions as representing outgoing and incoming waves.

So, why is this distinction so vital? Because in the real world, waves don't just exist in empty space; they interact with things. This is the phenomenon of ​​scattering​​. Imagine a plane wave of light (which can be described using the regular Bessel functions, $j_l$) hitting a tiny dust particle. The particle acts like a new, tiny source, radiating waves in all directions. This new wave is the ​​scattered wave​​.

What kind of wave must this be? The scattered wave is created by the particle, so it must carry energy away from the particle to observers far away. It must be a purely ​​outgoing​​ wave. Therefore, to describe the scattered field in any scattering problem—be it light off a sphere (Mie scattering), an electron off an atom, or radar waves off an airplane—we must use spherical Hankel functions of the first kind, $h_l^{(1)}(kr)$. Using a Bessel function $j_l(kr)$ would be physically incorrect, as it implies that the scattered wave is also bringing energy in from infinity, which makes no sense. The Hankel functions are not just an option; they are a necessity imposed by the physical reality of causality and energy flow.

While their physical meaning is paramount, it is worth appreciating the beautiful mathematical structure of these functions.

• ​​They are bona fide solutions:​​ Spherical Hankel functions are exact solutions to the same spherical Bessel differential equation that the Bessel and Neumann functions solve. We can verify this by direct, if tedious, calculation.
• ​​They form a complete basis:​​ The two types, $h_l^{(1)}$ and $h_l^{(2)}$, are linearly independent. This is confirmed by calculating their ​​Wronskian determinant​​, a mathematical tool that is non-zero for independent solutions. This means that any possible solution to the radial wave equation outside the origin can be written as a combination of incoming and outgoing Hankel functions.
• ​​They are not monstrously complex:​​ For low orders of $l$, they can be written explicitly using elementary functions like polynomials in $1/z$ and the exponential function $e^{iz}$. They are "special functions" not because they are esoteric, but because they are purpose-built for a special, and very common, job.
• ​​They have an orderly structure:​​ These functions obey a set of simple ​​recurrence relations​​ that connect functions of different orders and their derivatives. These relations provide a powerful computational toolkit, allowing for the elegant manipulation and simplification of seemingly complex expressions.

In the end, spherical Hankel functions are a perfect example of mathematics in service of physics. They arise from a fundamental physical question—how to describe a traveling wave in three dimensions—and through a touch of mathematical elegance involving complex numbers, they provide an answer that is not only correct but also deeply insightful, forming the bedrock of our understanding of scattering and radiation.

In our previous discussion, we met the spherical Hankel functions, $h_l^{(1)}(kr)$ and $h_l^{(2)}(kr)$. We saw them as nature's chosen solutions for the Helmholtz equation in spherical coordinates, describing waves that travel gracefully outwards toward infinity or converge inwards toward a central point. At first glance, this might seem like a neat mathematical trick, a specialized tool for a narrow class of problems. But nothing could be further from the truth. The universe is filled with things that are, to a good approximation, spherical, and it is replete with phenomena that involve waves radiating from or scattering off of these things.

The spherical Hankel functions, therefore, are not just abstract mathematical constructs; they are the very language we use to describe a startlingly broad array of physical realities. In this chapter, we will embark on a journey to see these functions in action. We will see how this single mathematical concept provides a unified framework for understanding everything from the sound of a drum to the color of the sky, from the quantum behavior of a particle to the tremors of an earthquake.

Let's begin with something we can all imagine: sound. What is the simplest possible source of sound? Perhaps a tiny, pulsating sphere, rhythmically expanding and contracting, pushing the air around it. This "breathing" sphere, or monopole source, sends out pressure waves in all directions. How does the pressure change as we move away from the source? The answer is described perfectly by the simplest spherical Hankel function, $h_0^{(1)}(kr)$. This function tells us that far from the source, the pressure amplitude of the outgoing sound wave decreases as $1/r$. This makes perfect intuitive sense: the energy of the wave spreads out over the surface of an ever-larger sphere, so its intensity must fall off as $1/r^2$, and its amplitude as $1/r$. The spherical Hankel function has this behavior built right in. This simple model is the foundation for understanding almost any compact sound source, from a small engine to a popping balloon, provided the source is much smaller than the wavelength of the sound it produces.

Of course, not all sources just "breathe." What if a sphere vibrates back and forth, like a tiny spherical bell? This is a dipole source, and it creates a more complex sound pattern. It pushes air forward in one direction while pulling it in from the opposite direction. This corresponds to the $l=1$ mode in our spherical expansion. The radial part of its wave is described not by $h_0^{(1)}(kr)$, but by $h_1^{(1)}(kr)$. By combining these fundamental modes—monopole ($l=0$), dipole ($l=1$), quadrupole ($l=2$), and so on, each with its corresponding spherical Hankel function—we can describe the sound field radiated by any vibrating object, no matter how complex its motion. This "partial wave expansion" allows us to calculate crucial physical quantities, like the total acoustic power radiated by a source, which tells us how efficient it is at turning vibration into sound.

The leap from the classical world of sound to the strange realm of quantum mechanics might seem vast, but the mathematical language remains the same. According to quantum theory, a particle like an electron is also a wave, described by a wavefunction $\psi$. The evolution of this wavefunction is governed by the Schrödinger equation which, for a free particle, is just the Helmholtz equation in disguise.

Imagine a particle scattering off an atom. After the interaction, the particle flies away. How do we describe this? As an outgoing matter wave. And what is the mathematical description of a spherically outgoing wave with a definite angular momentum $l$? It's precisely a spherical Hankel function, $h_l^{(1)}(kr)$. In this context, the square of the wavefunction, $|\psi|^2$, represents a probability density. An outgoing Hankel function therefore describes a net outward flow of probability—the particle is, quite literally, leaving the scene.

The power of this formalism becomes even more apparent when we consider phenomena like quantum tunneling. Consider a particle trying to get past a spherical potential energy barrier, like an electron interacting with a hollow, charged shell. The wavefunction outside the shell will be a superposition of an incoming wave (perhaps described by $h_l^{(2)}(kr)$) and an outgoing, scattered wave ($h_l^{(1)}(kr)$). By matching the wavefunction at the boundaries of the shell, we can connect the coefficients of the incoming and outgoing waves. A powerful tool called the transfer matrix method formalizes this process, allowing us to calculate exactly how much of the particle's wave is reflected and how much is transmitted through the barrier, all in the language of spherical Hankel functions and their relatives.

Perhaps the most celebrated application of spherical Hankel functions is in the theory of scattering. The fundamental question of scattering theory is: what happens when a wave hits an obstacle?

Let's picture a plane wave—be it an acoustic wave or a light wave—marching forward and encountering a sphere. The wave that emerges on the other side is a combination of the original, undisturbed plane wave and a new wave, the scattered wave, which emanates from the sphere. This scattered wave must, by its very nature, be composed of purely outgoing waves. Therefore, we can immediately write it down as a sum—a partial wave expansion—of terms containing spherical Hankel functions $h_l^{(1)}(kr)$ multiplied by angular functions $P_l(\cos\theta)$.

The "only" thing left to do is figure out the coefficients in this sum. And those coefficients are determined entirely by the physical properties of the sphere, which manifest as boundary conditions at its surface. For instance, if the sphere is "soft" (for an acoustic wave, this means the total pressure on its surface is zero), we can easily solve for the coefficients. In the limit of long wavelengths, where the sphere is tiny compared to the wave, this leads to a famous result: the total scattering cross-section is simply four times the geometrical cross-section of the sphere, $\sigma_{tot} = 4\pi a^2$. This framework is incredibly flexible; it can handle all sorts of boundary conditions that describe different physical interactions at the surface.

This exact approach, when applied to the full vector nature of electromagnetic waves, becomes the renowned Mie theory. It describes how a spherical particle—a water droplet in a cloud, a dust mote in the air, or a gold nanoparticle in a biological sample—scatters light. Each term in the expansion corresponds to a different multipole oscillation of charges and currents induced in the sphere by the incoming light wave. The coefficients, found by enforcing the continuity of the electric and magnetic fields at the sphere's surface, depend on the sphere's size, its refractive index, and its magnetic permeability. This theory is the reason the sky is blue and clouds are white, and it is a workhorse of modern optics, atmospheric science, and nanotechnology.

The unifying power of the spherical Hankel function doesn't stop at sound, particles, and light. Its reach extends into a diverse range of scientific and engineering disciplines.

Consider the field of solid mechanics or seismology. When an earthquake sends a compressional (P) wave through the Earth's crust and it encounters a pocket of magma or a different type of rock—an inclusion—it scatters. The scattered waves, which include both P-waves and shear (S) waves, radiate outwards from the inclusion. Once again, their radial behavior is described by spherical Hankel functions. By analyzing these scattered waves picked up by seismographs, we can deduce the properties of the Earth's interior. The same principles apply on a smaller scale in materials science for the non-destructive testing of materials, using ultrasound to find defects or characterize composites.

Finally, in theoretical physics, Hankel functions are indispensable for constructing Green's functions. A Green's function is a fantastically powerful tool; it represents the response of a system to a single, localized point source. Once you have the Green's function, you can find the response to any arbitrary source distribution. When working in a region with a spherical boundary—say, the space outside a sphere—the Green's function is built directly from spherical Hankel functions of both the first and second kind. They are pieced together in a clever way that respects both the boundary conditions on the sphere's surface and the radiation condition at infinity, giving us a complete solution to the wave equation for any source configuration.

From the simple to the complex, from the concrete to the abstract, we find the same mathematical protagonists playing a central role. The spherical Hankel functions, born from the simple requirement of describing waves traveling to or from a point, provide a deep and beautiful unity to our understanding of the physical world. They are a testament to how nature, in its boundless complexity, often relies on a few elegant, recurring mathematical ideas.