Plane Wave Expansion

In the world of physics, waves come in many forms, but two of the most fundamental are the directed plane wave and the expanding spherical wave. One describes linear, uniform motion, while the other describes influence radiating from a central point. While they seem entirely distinct, a remarkable and powerful concept known as the plane wave expansion reveals they are deeply connected. This expansion shows how any plane wave can be perfectly reconstructed from an infinite sum of spherical waves, providing a transformative tool for solving complex problems. It addresses the challenge of analyzing interactions that have a natural spherical symmetry, like a particle scattering off an atom, when the incoming particle is best described by a linear plane wave.

This article explores the plane wave expansion, a cornerstone of modern physics. In the first chapter, ​​Principles and Mechanisms​​, we will deconstruct the famous Rayleigh formula, examining the mathematical ingredients—spherical Bessel functions and Legendre polynomials—and the physical reasoning that makes it work. Following that, the chapter on ​​Applications and Interdisciplinary Connections​​ will demonstrate how this theoretical recipe becomes a practical toolkit, enabling breakthroughs in quantum scattering theory and providing the very foundation for our understanding of electrons in solids.

Imagine you are standing on the shore of a vast, calm lake. A long, straight ripple—a perfect plane wave—travels across the surface towards you. It seems to be the simplest kind of wave imaginable, a uniform front marching forward with a single, well-defined direction. Now, imagine you toss a small pebble into the water. It creates a completely different pattern: circular ripples that expand outwards from the point of impact. These are spherical waves, defined by their origin, spreading their influence equally in all directions. At first glance, these two phenomena—the directed plane wave and the isotropic spherical wave—seem like polar opposites. One is about direction; the other is about a center.

But what if I told you that the straight, directed plane wave is actually a magnificent illusion, a grand conspiracy? What if that simple, straight ripple could be perfectly reconstructed by adding up an infinite number of these circular, spherical waves, all centered on a single, arbitrary point in the water? This is not a trick; it is one of the most elegant and powerful ideas in wave physics, known as the ​​plane wave expansion​​. It allows us to take a problem about a particle moving in a straight line and re-imagine it from a perspective of angular motion, a trick that is the key to unlocking the secrets of quantum scattering, antenna design, and the behavior of electrons in crystals.

Let’s write down the magic formula. A plane wave traveling along the z-axis, described in quantum mechanics by the wavefunction $\psi(\mathbf{r}) = \exp(ikz)$, can be written as an infinite sum:

This is the famous ​​Rayleigh plane wave expansion​​. Let's not be intimidated by the symbols. Think of it as a recipe. On the left is the final dish: the simple plane wave. On the right are the ingredients we must mix together in precise amounts. Each term in the sum, indexed by the integer $l=0, 1, 2, \dots$, is a "partial wave", each with a definite ​​angular momentum​​, a concept central to quantum mechanics.

The ingredients for each partial wave are:

• ​​The Radial Ingredient, $j_l(kr)$:​​ These are the ​​spherical Bessel functions​​. They tell us how the amplitude of each partial wave changes as we move away from our chosen origin ($r=0$). They are the "spherical" part of our spherical waves.

• ​​The Angular Ingredient, $P_l(\cos\theta)$:​​ These are the ​​Legendre polynomials​​. They describe the shape of each partial wave, telling us how its strength varies with the angle $\theta$ from the z-axis. Unlike a simple spherical wave, these components are not the same in all directions.

• ​​The Seasoning, $i^l (2l+1)$:​​ These are the carefully chosen coefficients—the weights and phases—that ensure all the partial waves add up just right. The factor $i^l$ is particularly interesting; it's a phase shift that is crucial for the delicate interference that makes this all work.

So, the expansion tells us that a particle with definite linear momentum (moving along $z$) is actually a superposition of states with definite angular momentum $l$. Let’s taste some of these ingredients to see how they work.

What is the simplest, most fundamental component of a plane wave? Let's look at the first term in our series, the one with zero angular momentum ($l=0$). This is called the ​​s-wave​​. For $l=0$, the recipe simplifies dramatically. The Legendre polynomial is just $P_0(\cos\theta) = 1$, meaning this component has no angular dependence at all—it is perfectly isotropic. The coefficient is simply $i^0 (2 \cdot 0 + 1) = 1$. So the entire $l=0$ term is just the spherical Bessel function $j_0(kr)$.

And what is this function? It has a beautifully simple form:

This is a profound result. A plane wave, whose very definition involves a single direction, contains within it a component that is a perfect spherical wave, oscillating like $\sin(kr)$ and fading with distance as $1/r$. It’s a standing wave in the radial direction, emanating from the origin. Finding a perfectly spherical wave hidden inside a perfectly linear one is the first hint of the deep unity this expansion reveals.

Before we go on, we must address a crucial point. When solving the underlying wave equation (the free-particle Schrödinger equation), mathematics actually provides two families of radial solutions: the well-behaved spherical Bessel functions, $j_l(kr)$, and their unruly cousins, the ​​spherical Neumann functions​​, $n_l(kr)$. Why do we only use the $j_l$ in our expansion? Here, physics imposes its authority on the mathematics. The Neumann functions have a fatal flaw: they all diverge, or "blow up," at the origin ($r=0$). Our plane wave, $\exp(ikz)$, is perfectly finite and smooth everywhere, including the origin. To build a smooth function, we can only use well-behaved building blocks. Including any Neumann functions would introduce an unphysical infinity at the origin, so they must be completely excluded. Physical reality is the ultimate arbiter.

What about the other terms in the series, the partial waves with non-zero angular momentum? For $l=1$ (the ​​p-wave​​), the angular part is $P_1(\cos\theta) = \cos\theta$. This shape is not a sphere; it looks like a dumbbell aligned along the z-axis, positive on one side and negative on the other. For $l=2$ (the ​​d-wave​​), the angular shape is $P_2(\cos\theta) = \frac{1}{2}(3\cos^2\theta - 1)$, which has lobes along the z-axis and a "belt" of opposite sign around its equator. Sometimes, it's more convenient to use a slightly different but equivalent set of angular functions called ​​spherical harmonics​​, $Y_{l,m}(\theta, \phi)$, which are directly related to the Legendre polynomials when the wave is symmetric around the z-axis ($m=0$).

Each partial wave is a product of its radial part $j_l(kr)$ and its unique angular shape $P_l(\cos\theta)$. So how can adding up all these strange, lobed shapes create a perfectly flat plane wave? It is a symphony of interference. In the forward direction ($\theta=0$), all the partial waves conspire to add up constructively, building the crest of the plane wave. In every other direction, their positive and negative lobes intricately cancel each other out, leading to zero amplitude. It is a breathtakingly delicate balance, an infinite orchestra of spherical waves playing in perfect harmony to produce the simple melody of a single, straight line. Far from the origin, each of these partial waves behaves like a simple sine wave with a phase shift that depends on $l$, written as $\frac{1}{kr}\sin(kr - l\pi/2)$. This behavior is essential for understanding how waves scatter off objects.

You might be wondering if those coefficients, $i^l (2l+1)$, were just guessed. They were not; they are a direct consequence of a powerful mathematical principle called ​​orthogonality​​.

Think of the Legendre polynomials, $P_l(\cos\theta)$, as a set of "pure shapes" on a sphere, much like the pure frequencies of a musical instrument. Just as any complex sound can be decomposed into a sum of pure frequencies (a Fourier series), any reasonable function of the angle $\theta$ can be decomposed into a sum of these pure Legendre polynomials. The "orthogonality" property is the tool that lets us do this. It states that if you multiply two different Legendre polynomials ($P_l$ and $P_{l'}$) and integrate over all angles, the result is exactly zero. They are independent, in a geometric sense.

To find how much of the "shape" $P_l$ is contained within the plane wave $\exp(ikz)$, we "project" the plane wave onto that shape. In mathematical terms, we multiply $\exp(ikr\cos\theta)$ by $P_l(\cos\theta)$ and integrate over all angles. Thanks to orthogonality, all the other partial waves vanish in this process, leaving us with only the coefficient for the $l$-th wave. This procedure confirms the coefficients in our recipe and, in a beautiful reversal of logic, provides a powerful integral representation for the spherical Bessel functions themselves:

The physics of the plane wave expansion gives us a profound mathematical identity for free!

This expansion is far more than a mathematical curiosity. It is a workhorse of modern physics and engineering. In ​​quantum scattering​​, a particle (an incoming plane wave) strikes a target. The target's force is typically short-ranged and spherically symmetric. The expansion allows us to analyze the interaction not as a whole, but one partial wave at a time. The target might only affect the s-wave ($l=0$) and p-wave ($l=1$), leaving the higher-$l$ waves untouched. By measuring how each partial wave is deflected (phase-shifted), we can deduce the nature of the force, a technique that was essential in unraveling the structure of the atomic nucleus.

Furthermore, the expansion of the simple plane wave $\exp(ikz)$ becomes a building block for more complex waves. For example, a standing wave like $\cos(kz)$ can be seen as the sum of two plane waves traveling in opposite directions, $\cos(kz) = (\exp(ikz) + \exp(-ikz))/2$. By using our recipe on each piece, we can find the partial wave expansion for the standing wave. We can even handle more exotic functions by manipulating the series mathematically, for instance, finding the expansion for a wave like $\cos\theta \cdot \exp(ikz)$ or for the derivative of a wave.

The idea of expanding functions into plane waves is a cornerstone of ​​solid-state physics​​, where it forms the basis of some of the most powerful methods for calculating the properties of materials. The journey from a simple, straight ripple to an infinite orchestra of spherical waves is a testament to the hidden unity and surprising beauty that mathematics brings to our understanding of the physical world.

At first glance, the plane wave expansion seems like a curious exercise in making things more complicated. We begin with the simplest of all waves—a perfectly flat, infinitely broad plane wave like $e^{i\mathbf{k} \cdot \mathbf{r}}$, the very picture of uniform, directed motion—and we exchange it for an infinite sum of spherical waves, each described by a tangle of Bessel functions and spherical harmonics. Why would we trade such elegant simplicity for this apparent complexity?

The answer, as is so often the case in physics, lies in a profound shift in perspective. The plane wave expansion is not just a mathematical identity; it is the key that unlocks the relationship between two fundamental kinds of motion: moving in a straight line and spinning in a circle. It is a bridge between the worlds of linear momentum and angular momentum. Once we cross this bridge, we find we have a powerful tool that can be used to probe the secrets of the atom, to design new materials from the ground up, and to connect seemingly disparate fields of science. The machinery is so robust that it can even be turned to solve problems in pure mathematics, such as taming ferocious-looking infinite series of special functions with surprising ease.

In quantum mechanics, a plane wave $e^{i\mathbf{k} \cdot \mathbf{r}}$ describes a particle with a perfectly defined linear momentum, $\mathbf{p} = \hbar\mathbf{k}$. If you measure its linear momentum, you will get this value, every single time. But what if you ask a different question: what is the particle's angular momentum with respect to some origin? Naively, an object moving in a straight line that doesn't pass through the origin should have a well-defined angular momentum. But a plane wave extends through all of space; it isn't a localized particle. The quantum answer is subtle and beautiful, and it is given to us directly by the plane wave expansion.

The expansion tells us that a particle in a plane-wave state is simultaneously in a superposition of all possible (non-negative integer) orbital angular momentum states. The spherical waves in the expansion, characterized by the spherical harmonics $Y_{lm}(\hat{\mathbf{r}})$, are states of definite angular momentum, with quantum number $l$ for the total magnitude and $m$ for its projection on an axis. When we measure the angular momentum of our plane-wave particle, we could find it to have $l=0$, or $l=1$, or $l=2$, and so on, each with a specific, calculable probability. The expansion provides the exact recipe for these probabilities.

This becomes even more concrete when we consider a state formed by the superposition of two plane waves, for example, one traveling along the z-axis and one along the x-axis. Such a state has a more complex spatial structure, and consequently, a different distribution of angular momentum content. By applying the expansion to each plane wave component and adding them together, we can precisely calculate the relative probability of measuring different angular momentum outcomes, such as finding the particle in a state with $(l=2, m=0)$ versus a state with $(l=2, m=2)$. This ability to translate between the language of linear motion and rotational motion is not merely a mathematical convenience; it is a cornerstone of how we analyze and interpret quantum phenomena.

Perhaps the most classic and powerful application of the plane wave expansion is in the theory of scattering. Imagine you want to know the shape and structure of an object hidden in a dark room. A good way to do this is to shout and listen carefully to the echoes. The incoming sound wave from your shout is uniform and planar, but the echo that returns is a complex wave, shaped and sculpted by its interaction with the object.

In physics, this is precisely what we do in a scattering experiment. We fire a beam of particles (like electrons, neutrons, or protons), which can be described as a plane wave, at a target (like an atomic nucleus or a crystal). The particles scatter off the target, and we place detectors at various angles to measure how many particles come out in each direction. The plane wave expansion gives us the language to understand this process.

The incident plane wave is decomposed into an infinite series of incoming spherical waves. Each of these spherical waves interacts with the target potential, which modifies it, typically by shifting its phase. The scattered wave is the sum of all these modified, outgoing spherical waves. The amount by which each spherical wave's phase is shifted, known as the phase shift $\delta_l$, is the fingerprint of the interaction. A strong, repulsive potential will produce a different set of phase shifts than a weak, attractive one.

The end result is that the outgoing pattern of scattered particles—the "echo"—is an interference pattern created by the superposition of all the phase-shifted spherical waves. If we are given the set of phase shifts, which might be calculated from a theoretical model of the target potential, we can use the expansion to reconstruct the entire angular distribution of scattered particles, a quantity known as the differential cross section, $d\sigma/d\Omega$. The interference between the different partial waves, each corresponding to a different angular momentum $l$, often leads to beautiful diffraction-like oscillations in the scattering pattern, which can be compared directly with experimental data. The same framework allows us to analyze simpler wave patterns, like a standing wave $\cos(kz)$, and determine its composition in terms of spherical waves, which is essential for understanding resonant scattering and particles confined in cavities.

The true power of a fundamental idea is often revealed when it is used to understand complex, collective phenomena. So it is with the plane wave expansion, which provides the very foundation for our modern understanding of electrons in crystalline solids. A crystal is a vast, regular array of atoms, containing an astronomical number of electrons, all interacting with each other and with the periodic lattice of atomic nuclei. Describing this system from first principles seems like an utterly hopeless task.

Yet, here a variation of the plane wave idea comes to the rescue in the form of Bloch's Theorem. This theorem, a direct consequence of the crystal's periodic symmetry, states that the wavefunction of an electron moving in the crystal's periodic potential is not just any complicated function. It must take the form of a plane wave, $e^{i\mathbf{k} \cdot \mathbf{r}}$, modulated by a function, $u_{n\mathbf{k}}(\mathbf{r})$, that has the exact same periodicity as the crystal lattice itself.

This is a monumental simplification. It means we don't have to solve for the electron's wavefunction throughout the entire infinite crystal. Instead, we only need to solve for the periodic part, $u_{n\mathbf{k}}(\mathbf{r})$, within a single, tiny repeating unit of the crystal—the primitive unit cell. The overall behavior is then determined by the "crystal momentum" vector $\mathbf{k}$. This masterstroke transforms the impossible problem of an infinite system into a manageable problem solved within a single unit cell, which is then repeated for a representative sample of $\mathbf{k}$-vectors in the Brillouin zone. This block-diagonalization of the problem is precisely what makes modern computational materials science feasible.

And how do we describe the wavefunction within that unit cell? With another expansion of plane waves! Because the function $u_{n\mathbf{k}}(\mathbf{r})$ is periodic, it can be represented by a Fourier series, which is an expansion in plane waves whose wavevectors are the reciprocal lattice vectors of the crystal. This "plane-wave basis set" is the engine behind Density Functional Theory (DFT), the workhorse method for calculating the electronic structure and properties of almost any material imaginable.

This powerful framework can be further specialized to tackle problems in surface science. In techniques like Low-Energy Electron Diffraction (LEED), scientists shoot a beam of low-energy electrons at a surface to determine the precise arrangement of the atoms there. The theoretical modeling of this process is a beautiful application of the plane wave expansion. An incoming electron plane wave is modeled as it scatters from a two-dimensional periodic layer of atoms. The expansion is used twice: once to express the incoming plane wave as a sum of spherical waves centered on an atom, and again to express the spherical wave scattered from that atom as a sum of new plane waves traveling in all the allowed diffracted directions. This allows physicists to construct a "scattering matrix" for a single atomic layer, a crucial building block for simulating the entire experiment and deciphering the surface structure.

Of course, in a real crystal or on a surface, the electron waves scatter from atoms located at many different positions, not just at the origin of our coordinate system. The formalism is flexible enough to handle this. The plane wave expansion can be generalized to expand a wave around any arbitrary point in space or to transform the expansion coefficients when the physical system is more naturally described in a rotated frame of reference.

From the quantum dance of a single particle to the collective electronic symphony of a solid, the plane wave expansion is far more than a mathematical formula. It is a fundamental statement about the dual nature of motion, a versatile tool for calculation, and a conceptual bridge that unifies vast and seemingly disconnected fields of physics. It teaches us that sometimes, the most powerful way to understand a simple thing is to see it as the sum of many other simple things, viewed from just the right perspective.