Zeros of Bessel Functions

Bessel functions and their zeros represent a cornerstone of mathematical physics, yet they can often appear as abstract and esoteric concepts. These special numbers, which are the roots of the Bessel function equations, seem like mere mathematical curiosities at first glance. However, they are deeply woven into the fabric of the physical world, dictating the behavior of phenomena constrained by cylindrical symmetry. This article bridges the gap between their abstract definition and their tangible reality, revealing the beautiful and surprisingly simple rules that govern them and their ubiquitous role across science.

The following chapters will guide you on a journey of discovery. In "Principles and Mechanisms," we will demystify the mathematical properties of these zeros, uncovering their connections to familiar trigonometric functions, the elegant interplay between different function orders, and the universal rhythm they follow at large scales. Subsequently, in "Applications and Interdisciplinary Connections," we will see these abstract numbers come to life, exploring how they define everything from the sound of a drum and the resolving power of a telescope to the fundamental energy levels of quantum systems and the structure of the quantum vacuum itself.

Now that we have been introduced to the world of Bessel functions, let's pull back the curtain and have a look at the machinery inside. At first glance, these functions and their zeros can seem abstract and intimidating, a menagerie of squiggles defined by arcane formulas. But as we start to explore, we find that they are not so alien. In fact, they are deeply connected to concepts we already know and love, and they possess a hidden structure and rhythm that is both beautiful and profoundly useful.

Let's begin our journey by dispelling some of the mystery. You might think that finding the zeros of a Bessel function requires high-powered numerical methods. Sometimes it does, but in certain wonderful cases, the function reveals itself to be an old friend in disguise.

Consider the Bessel function of order $\nu = -1/2$, which we call $J_{-1/2}(x)$. Its definition involves an infinite series and the Gamma function, which sounds complicated enough. But through a bit of mathematical alchemy using standard identities, this function miraculously simplifies into something you've known for years. It turns out that for positive $x$, $J_{-1/2}(x)$ is just a cosine function with a decaying amplitude:

Suddenly, the problem of finding its zeros becomes trivial! Where does this function equal zero? Well, the $\sqrt{2/(\pi x)}$ part is never zero for a finite, positive $x$. So, the zeros must be exactly where the cosine function is zero. And we all know where that is: at $\pi/2$, $3\pi/2$, $5\pi/2$, and so on. We can write a simple formula for the $n$-th positive zero: it's just $(n - 1/2)\pi$.

This is a wonderful first step. It tells us that buried within the complex world of Bessel functions are connections to the simple, periodic world of trigonometry. It gives us a foothold, a sense of confidence that we can, in fact, understand these things intuitively.

Of course, not all Bessel functions are so simple. The famous $J_0(x)$, which describes the fundamental vibration of a circular drumhead, does not simplify to a basic trigonometric function. Its zeros—the circular lines on the drumhead that remain perfectly still while everything else vibrates—are at the seemingly random-looking positions $2.4048\dots$, $5.5201\dots$, $8.6537\dots$, and so on.

But even here, there is a beautiful, hidden order. The behavior of $J_0(x)$ is intimately tied to the behavior of its neighbor, $J_1(x)$. Think about a swinging pendulum. At the highest points of its arc, where it momentarily stops before changing direction, its velocity is zero. In the same way, the function $J_0(x)$ reaches its local maxima and minima—its peaks and troughs—at points where its "velocity," or derivative, is zero. And here is the magic: the derivative of $J_0(x)$ is nothing other than $-J_1(x)$.

This means that to find the peaks and valleys of the vibrating drumhead pattern described by $J_0(x)$, you just need to find the ​​zeros​​ of the $J_1(x)$ function. The zeros of one function dictate the extrema of another. It's a beautiful dance between the different orders. The nodal lines of the $J_1$ vibrational mode are precisely the anti-nodes (points of maximum displacement) for the $J_0$ mode. This is not a coincidence; it is a fundamental property that stems directly from the differential equation that all Bessel functions must obey.

Let's zoom out. What happens to the zeros far away from the origin? If we trace the wiggles of any Bessel function $J_\nu(x)$ for very large values of $x$, a remarkable simplification occurs. The frantic, complex oscillations near the center calm down and settle into a predictable, universal rhythm.

For large $x$, every Bessel function $J_\nu(x)$ starts to look like a simple, decaying cosine wave:

Look at what this tells us! The zeros of the function, for large $x$, will be approximately where the cosine term is zero. The argument of the cosine is $x$ plus some constant phase shift. Since the zeros of $\cos(\theta)$ are separated by $\pi$, this means that the spacing between consecutive large zeros of any Bessel function $J_\nu(x)$ approaches ​​$\pi$​​.

This is a spectacular result. The specific order $\nu$ and the complex behavior near the origin fade into irrelevance at large distances. Out on the horizon, all Bessel functions march to the beat of the same drum, a rhythm governed by the fundamental constant $\pi$. This asymptotic simplicity is not just a mathematical curiosity; it's what allows engineers and physicists to understand and predict high-frequency wave phenomena in cylindrical structures, from optical fibers to acoustic ducts. Of course, this spacing is only approximately $\pi$. More advanced analysis reveals that the spacing itself approaches $\pi$ with correction terms that get smaller and smaller as the zeros get larger, in a very precise and predictable way.

So the zeros are regularly spaced far out. But what about their collective properties? Let's ask a strange-sounding question. If we take all the positive zeros of a Bessel function, $j_{\nu,1}, j_{\nu,2}, j_{\nu,3}, \dots$, and add up the reciprocal of their squares, what do we get?

This seems like a task for a supercomputer, adding up an infinite list of strange numbers. Yet, remarkably, there is an exact, elegant answer that we can find with a piece of paper and some clever thinking. The trick lies in a powerful idea from complex analysis. Just as a polynomial like $x^2 - 5x + 6$ can be defined by its coefficients or, equally well, by its roots $(x-2)(x-3)$, many important functions can be described in two ways: one is their ​​power series​​ (like a Taylor series), and the other is an ​​infinite product​​ built from their zeros. This is the essence of the Weierstrass factorization theorem.

By writing down the first few terms of the power series for a function related to $J_\nu(z)$ and comparing them to the expansion of its infinite product, we can perform a sort of mathematical miracle. The coefficient of the $z^2$ term in both expansions must be the same. One expression for this coefficient contains the sum we want, and the other contains the order $\nu$. Equating them yields a stunningly simple formula that holds for any $\nu > -1$:

For the drumhead function $J_0(x)$, this sum is exactly $1/4$. For $J_1(x)$, it's $1/8$. There is an incredible, hidden rigidity in the placement of these zeros. Their values are not arbitrary; they are constrained in such a way that the sum of their inverse squares is this simple, rational number. It's a deep statement about the function's structure, revealed by looking at it from two different points of view simultaneously. This tells us something profound: the rate at which the zeros spread out is precisely controlled, ensuring that the series $\sum 1/|z_n|^2$ converges, while the simpler series $\sum 1/|z_n|$ does not.

Is this just a trick of the complex analysis trade? Or does it represent a deeper physical truth? Let's call a second witness to the stand: the theory of waves and Fourier analysis.

Imagine you tap a drumhead, setting it to an initial shape—say, you displace the whole surface by a uniform height of 1. This shape can be represented as a sum of the drum's natural vibration modes, which are precisely the Bessel functions $J_0(\alpha_n x/L)$, where $\alpha_n$ are the zeros of $J_0$. This is a ​​Fourier-Bessel series​​.

A fundamental principle in physics, related to the conservation of energy, is ​​Parseval's theorem​​. It states that the total energy of the initial shape must equal the sum of the energies contained in all of its vibrational modes. We can calculate the "energy" of our initial flat displacement on the left side of an equation. On the right side, we can write down the sum of the "energies" of all the Bessel modes needed to build it.

When we carry out this calculation for a constant displacement $f(x)=1$, the coefficients of the expansion depend on the zeros $\alpha_n$. When we plug everything into Parseval's identity and simplify, the dust settles to reveal an equation that relates the energy to our mysterious sum. The result?

We get the exact same answer. This is no accident. It is a powerful confirmation of the correctness and consistency of our mathematical description of the world. One path, through the abstract world of infinite products in the complex plane, and another, through the physical world of wave energies and orthogonal functions, lead to the identical, beautiful result. The locations of the zeros of Bessel functions are not just mathematical artifacts; they are fundamental constants woven into the fabric of geometry and physics, revealed to us through the unified language of mathematics.

It is one of the most remarkable and beautiful aspects of physics that a single mathematical idea can reappear, as if by magic, in the most disparate corners of the universe. The set of numbers we call the zeros of Bessel functions provides a stunning example of this unity. At first glance, they are merely the points where a particular oscillating function crosses the zero line. But as we look closer, we find that these numbers are not just mathematical curiosities; they are nature's "allowed" values, the fingerprints left by physical laws when they are constrained by cylindrical symmetry. Let us take a journey through science and see where these special numbers emerge, from the familiar vibrations of our world to the deepest puzzles of the quantum realm.

Perhaps the most intuitive place to meet the zeros of Bessel functions is in the sound of a drum. Imagine striking a circular drumhead. It begins to vibrate, producing a sound. But unlike a simple guitar string, which can only vibrate up and down along a line, the drumhead is a two-dimensional surface. It can vibrate in wonderfully complex patterns of hills and valleys. The patterns that produce a clear, sustained tone are called "normal modes," and they are standing waves on the membrane's surface.

Because the edge of the drumhead is clamped down and cannot move, any possible standing wave pattern must have zero displacement at the circular boundary. The mathematical description of the wave's shape away from the center involves Bessel functions. The fixed-boundary rule forces a strict condition: the Bessel function describing the radial part of the wave must be equal to zero precisely at the radius of the drum. This means the argument of the function, which depends on the frequency of vibration, cannot be just anything. It must be one of the special values—a zero of the appropriate Bessel function. Each zero corresponds to an allowed vibrational mode, each with its own unique frequency and pattern. The fundamental tone corresponds to the first zero, and the overtones, which give the drum its rich character, correspond to the subsequent zeros. Furthermore, for some of these higher-frequency modes, there are not only nodal lines at the edge but also perfect circles of silence within the drumhead where the surface remains completely still. The radii of these silent nodal circles are determined by the intermediate zeros of the same Bessel function that governs that mode. In this way, the abstract list of zeros is transformed into the audible spectrum of a percussion instrument.

Let's now turn from waves of sound to waves of light. When light passes through any opening, it diffracts, spreading out and creating a pattern of interference. If the opening is a perfect circle, like the aperture of a camera or a telescope, the diffraction pattern is not a simple spot but a beautiful bullseye known as an Airy pattern. It consists of a bright central disk surrounded by a series of concentric, progressively fainter bright rings, separated by rings of pure darkness.

Where do these dark rings come from? They are the locations where light waves, having traveled from different parts of the aperture, arrive perfectly out of phase and cancel each other out completely. The mathematics of this interference shows that the angles at which these dark rings appear are governed by the zeros of the Bessel function of order one, $J_1(x)$. Each dark ring corresponds directly to one of the zeros. This is not just a laboratory curiosity; it has profound consequences. The famous Rayleigh criterion, which defines the resolving power of any optical instrument, is based on this very principle. It states that two distant stars are just barely distinguishable when the central bright peak of one star's Airy pattern falls directly on top of the first dark ring of the other. The position of this first dark ring is dictated by the first zero of $J_1(x)$. Thus, a fundamental limit on our ability to see the universe is written in the language of Bessel function zeros.

So far we have seen waves expanding from a source. But what if we want to channel them, to guide them from one place to another? This is the job of a waveguide, which is essentially a metal pipe used to transport high-frequency electromagnetic waves, like microwaves or signals in a particle accelerator. When we confine waves within a conducting cylinder, we once again impose strict boundary conditions. The electric and magnetic fields must behave in a prescribed way at the metallic walls.

These boundary conditions act as a filter, allowing only certain wave patterns, or "modes," to propagate. For one class of modes, the Transverse Magnetic (TM) modes, the physics is similar to the drumhead: the longitudinal electric field must be zero at the walls, a condition satisfied only if a parameter related to the wave's frequency corresponds to a zero of a Bessel function $J_m(x)$. For another class, the Transverse Electric (TE) modes, the boundary condition is different: it requires the rate of change (the derivative) of the longitudinal magnetic field to be zero at the walls. This leads to a new set of allowed modes, determined not by the zeros of $J_m(x)$ itself, but by the zeros of its derivative, $J'_m(x)$. Engineers use these principles to design communication systems, selecting waveguide dimensions and materials to ensure that only desired frequencies and modes can travel through them, a task where knowledge of both types of Bessel zeros is indispensable.

The leap from classical waves to quantum mechanics is a giant one, yet our Bessel functions come right along with us. According to quantum theory, a particle like an electron is also a wave—a wave of probability. If we confine a particle to a cylindrical region, say a tiny carbon nanotube, we have created a "quantum drum." The particle's wavefunction must vanish at the impenetrable walls, just like the displacement of the drumhead at its rim.

Unsurprisingly, the Schrödinger equation for this situation gives solutions involving Bessel functions. The boundary conditions once again quantize the system. The particle cannot have just any energy; it is restricted to a discrete set of allowed energy levels. And what determines these levels? The squares of the zeros of Bessel functions. A fascinating twist arises if we thread a line of magnetic flux down the center of the cylinder. A particle moving in this space experiences the Aharonov-Bohm effect: its energy levels shift even though it never touches the region with the magnetic field. This mysterious, non-local quantum interaction manifests itself mathematically by changing the order of the Bessel function whose zeros determine the energy, providing a deep connection between geometry, electromagnetism, and quantum mechanics.

The influence of our Bessel zeros extends to the most modern and exotic frontiers of physics. In quantum field theory, the vacuum is not an empty void but a seething cauldron of "virtual" particle-antiparticle pairs that flicker in and out of existence. If we confine this vacuum inside a cylinder, the boundary conditions restrict the modes of these virtual fluctuations. This restriction creates a real, measurable energy density in the vacuum known as the Casimir effect. In a simplified but insightful model of this phenomenon, the vacuum energy density on the cylinder's axis is proportional to an infinite sum involving the inverse squares of all the zeros of the $J_0$ Bessel function. Miraculously, this infinite sum converges to the simple, elegant value of $\frac{1}{4}$, allowing for a concrete prediction of a force from "nothing".

Equally fascinating is the behavior of a Bose-Einstein Condensate (BEC), a state of matter where millions of atoms cool down to such a low temperature that they behave as a single, macroscopic quantum entity. If you rotate a bucket of this quantum fluid, it doesn't form a smooth whirlpool like water. Instead, it develops quantized vortices—tiny, stable whirlpools where the angular momentum is a multiple of a fundamental unit. It turns out that there is a critical angular velocity at which it becomes energetically favorable for the first vortex to appear at the center of the condensate. This critical velocity is determined by the energy difference between the non-rotating ground state and the first rotating state. These energies, in turn, are governed by the zeros of the Bessel functions $J_0$ and $J_1$, respectively. The birth of a quantum whirlpool is thus heralded by a competition between two of nature's "allowed" values.

Given how ubiquitous these zeros are, one might wonder if they can be used in other mathematical contexts, for instance, in numerical methods. The technique of Gaussian quadrature is a powerful way to approximate definite integrals by sampling a function at a small number of cleverly chosen points. These points are the roots of a special class of functions called orthogonal polynomials. Could the zeros of a Bessel function serve as these optimal sampling points? The answer, surprisingly, is no. There is a deep mathematical reason: the roots of consecutive orthogonal polynomials in a sequence must strictly interlace one another. The zeros of a single Bessel function, taken in increasing order, do not satisfy this interlacing property with respect to each other. This shows that Bessel functions, while powerful, belong to a different mathematical family than the classical orthogonal polynomials. It highlights that the applicability of a mathematical tool is as much about the structures it doesn't fit as the ones it does.

From the tangible sound of a drum to the intangible energy of the quantum vacuum, the zeros of Bessel functions thread a unifying needle through the fabric of physics. They are a powerful reminder that the universe, in all its complexity, often speaks in a surprisingly simple and unified mathematical language.