🔬 optics

Chebyshev self-imaging from inverse-sampled angular spectra

This paper demonstrates that by sampling the angular spectrum in inverse powers of the mode index, conventional self-imaging can be programmably transformed into a "prime-power staircase" of revivals governed by the second Chebyshev function, thereby using free-space diffraction to physically read out the multiplicative structure of integers.

Original authors: Layton A. Hall, Murat Yessenov

Published 2026-08-25

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Layton A. Hall, Murat Yessenov

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

For nearly two hundred years, physicists have watched a peculiar trick played by light. When a beam of light passes through a regular grid, the waves that emerge do not simply spread out and blur; instead, they periodically snap back into their original shape, recreating the grid pattern at precise distances without any lenses or mirrors to help them. This phenomenon, known as self-imaging, has been a reliable rule of nature: the distance at which the image reappears depends only on the spacing of the grid, not on how many different colors or frequencies of light are mixed together. It is a predictable, rhythmic return to the start, governed by a simple arithmetic that treats all the waves in the beam as part of a single, synchronized family.

A team of researchers at Los Alamos National Laboratory and Harvard University has now shown that this rule is not fixed, but can be rewritten. By carefully arranging the light waves so that their frequencies follow a specific, unusual pattern, they turned the predictable return of the image into a complex calculation that reveals the hidden structure of numbers. In their experiments, the light no longer returns at a single, fixed distance. Instead, the distance it must travel to rebuild itself changes depending on exactly which waves are present in the beam. If the researchers add a new wave to the mix, the image might stay exactly where it was, or it might suddenly jump to a much farther point. The distance the light travels becomes a physical readout of the mathematical properties of the numbers that define the waves, specifically how those numbers are built from prime factors.

The researchers achieved this by changing how they sampled the light. In a standard setup, the light waves are spaced out in a simple, even sequence. In this new work, the team programmed a device to create a beam where the spacing between the waves followed a specific rule: the higher the number assigned to a wave, the more closely packed it became, but in a way that inverted the usual relationship. They then sent this specially crafted beam through free space and watched how it evolved. As the light traveled, the different waves within it drifted in and out of step with one another. In a normal beam, they all line up again at the same moment. In this engineered beam, each wave had its own unique schedule for lining up. The entire image could only reappear when every single wave in the group happened to be in the right place at the same time.

This requirement for total synchronization created a surprising result. The distance the light had to travel to reform the image was determined by the least common multiple of the numbers assigned to the waves. In mathematics, the least common multiple is the smallest number that can be divided evenly by a set of other numbers. For a beam containing waves labeled with the numbers 1 through 4, the light had to travel a distance corresponding to the number 12. If the researchers added a wave labeled 5, the distance jumped to 60. However, if they added a wave labeled 6, the distance did not change at all, because 6 is made of factors that were already present in the previous set. The image would only jump again when a new prime number or a new power of a prime number was introduced.

The team observed this behavior directly using a laser and a camera. They programmed the light to contain different sets of waves, starting with just a few and gradually adding more. When they looked at the beam as it traveled, they saw the image fade and then sharpen back into focus at specific points. For a set of waves labeled 2 through 5, the image reappeared at a distance of 60 units. When they added the wave labeled 6, the image reappeared at the exact same distance of 60 units, confirming that the new wave did not alter the schedule. When they added the wave labeled 7, the distance jumped to 420 units. The researchers measured these distances with high precision, finding that the points where the image returned matched the mathematical predictions perfectly. The distance between the start and the return of the image formed a "staircase" that stepped up only when the beam included a new prime number or a prime number raised to a power, such as 4 or 9.

This work demonstrates that the behavior of light can be made to reflect the multiplicative structure of integers, rather than just their additive properties. In the past, scientists have used light to study the additive nature of numbers, but this experiment shows that by changing the way the light is prepared, the laws of physics can be made to perform a different kind of arithmetic. The distance the light travels becomes a physical manifestation of the least common multiple, a fundamental concept in number theory. The researchers noted that this effect is not limited to light; similar principles could apply to other types of waves, including sound or even the waves associated with particles in quantum mechanics. They suggested that this approach could serve as a classical way to simulate the behavior of electrons in atoms, where the energy levels follow a similar inverse pattern.

The findings also highlight a new way to think about the limits of optical systems. Usually, adding more waves to a beam is seen as a way to increase the detail or resolution of an image. Here, adding more waves changes the fundamental distance at which the image can be reconstructed. The researchers found that as they added more waves, the distance required for the image to return grew exponentially. This means that while the effect is powerful, it is also constrained by the physical size of the setup and the precision of the equipment. To see the effect with a large number of waves, the light would need to travel a distance that quickly becomes impractical in a standard laboratory. However, the researchers pointed out that this same principle could be applied to light traveling through optical fibers, where the distance could be extended to kilometers, potentially allowing for the observation of these effects with much larger sets of numbers.

The experiment serves as a bridge between the physical world of waves and the abstract world of number theory. By simply rearranging the frequencies of light, the researchers turned a beam of light into a machine that calculates the least common multiple of a set of numbers. The image that appears at the end of the journey is not just a picture of the source; it is a confirmation of a mathematical truth. The light travels, drifts, and realigns, and in doing so, it reveals the hidden architecture of the integers. The researchers did not just observe a new optical effect; they showed that the rules governing wave propagation can be programmed to expose the deep, structural properties of mathematics, turning the act of watching light travel into a way of reading the prime factors of numbers.

Technical Summary: Chebyshev Self-Imaging from Inverse-Sampled Angular Spectra

Problem Statement
Conventional self-imaging, exemplified by the Talbot effect, relies on the quadratic phase accumulation of spatial frequency components during paraxial propagation. In this regime, a transversely periodic field reconstructs itself at a fixed distance (the Talbot length, $z_T$ ) and its integer multiples. This revival distance is independent of the field's bandwidth because the harmonic phases ( $\propto m^2$ ) are commensurate; all modes rephase simultaneously at integer multiples of a single period. While previous work has utilized this quadratic phase kernel to evaluate Gauss sums and factor integers via interferometry, the underlying propagation arithmetic remains fixed and quadratic. The authors ask whether the arithmetic of self-imaging itself can be engineered to expose or compute number-theoretic structures directly, rather than merely reading them against a fixed kernel.

Methodology
The authors propose a method to alter the fundamental arithmetic of self-imaging by engineering the angular spectrum of the wavefield. Instead of the standard linear spacing of spatial frequencies, they sample the angular spectrum in inverse powers of the mode index $m$ .

The field is modeled as a superposition of $M$ transverse plane waves:
$E(x,z) = \sum_{m=1}^{M} a_m e^{ik_x(m)x} e^{-ik_x^2(m)z/2k_o}$
where the transverse spatial frequency is defined as $k_x(m) = k_L m^{s/2}$ .

Conventional Case ( $s > 0$ ): For $s=2$ (Talbot) or $s=1$ (Montgomery effect), the rephasing periods nest within the unit interval, leading to a universal revival distance independent of $M$ .
Inverse-Sampled Case ( $s = -r$ ): The authors set $s = -r$ (where $r$ is a positive integer), resulting in $k_x(m) = k_L m^{-r/2}$ . In this configuration, the phase of mode $m$ evolves as $\exp(-2\pi i \zeta / m^r)$ , where $\zeta = z/z_T$ .

Under this inverse sampling, the rephasing period for mode $m$ becomes proportional to $m^r$ . Exact self-imaging ( $E(x,z) = E(x,0)$ ) occurs only when the propagation distance $\zeta$ is a common multiple of all individual periods $\{m^r\}_{m=1}^M$ . Consequently, the first exact revival distance is determined by the least common multiple (LCM) of the mode indices raised to the power $r$ :
$z_{rev}(M) = z_T \cdot \text{lcm}(1, 2, \dots, M)^r$

Key Contributions and Results
The paper demonstrates that this engineering transforms the revival distance from a fixed constant into a programmable quantity governed by the multiplicative structure of integers.

Chebyshev Revivals and the Prime-Power Staircase:
The revival coordinate is expressed as $z_{rev}/z_T = e^{r\psi(M)}$ , where $\psi(M)$ is the second Chebyshev function ( $\psi(M) = \ln(\text{lcm}(1, \dots, M))$ ). The authors observe that the revival distance remains constant as modes are added, except when the mode index $M$ reaches a new prime power $p^a$ . At these specific indices, the distance jumps by a factor of $p^r$ . This creates a "prime-power staircase" in propagation distance, turning free-space diffraction into a physical readout of the multiplicative structure of integers.
Divisor Revivals:
Between exact revivals, the field exhibits "divisor revivals." At intermediate integer planes $\zeta$ , subsets of modes whose indices divide $\zeta$ rephase, forming partial self-images. The on-axis intensity $F(\zeta)$ displays a hierarchy of peaks corresponding to the divisor lattice of the integers up to $M$ . The first exact revival is distinct from these partial peaks as it is the first point where all modes rephase simultaneously.
Experimental Demonstration:
Using a continuous-wave laser ( $\lambda = 635$ nm) and a phase-only Spatial Light Modulator (SLM), the authors programmed inverse-sampled spectra ( $r=1$ ) and measured the axial intensity.
- For consecutive spectra $S = \{2, \dots, M\}$ , they observed the first revivals at $\zeta = 6, 12, 60, 60$ for $M=3, 4, 5, 6$ . The plateau between $M=5$ and $M=6$ confirms that adding a composite number (6) does not change the LCM, whereas adding a prime (5) caused a fivefold increase.
- Sparse spectra (e.g., $S=\{4,5,6,7\}$ ) were shown to achieve the same revival distance as the full set $\{1, \dots, 7\}$ , demonstrating that the recurrence is programmable based on the prime-power content of the occupied bandwidth, not just the bandwidth size.
- Measurements for $r=2$ (inverse-square sampling) were also presented, showing jumps of $p^2$ and confirming the squared divisor lattice structure.

Significance
The paper establishes that the revival distance of a wave is not an intrinsic, fixed property of propagation but a programmable arithmetic quantity. By shaping modal phase rates, the authors shift the governing principle from the additive commensurability of the Talbot effect to the multiplicative structure of the integers.

The significance lies in three main areas:

Fundamental Physics: It inverts the arithmetic of self-imaging, deferring exact revivals to the LCM of mode periods rather than a universal constant.
Number-Theoretic Readout: The propagation distance itself performs a primality and prime-power test on the mode index. The "carpet" of intensity between launch and revival physically renders the divisor lattice of the integers.
Quantum Analogy: The $r=2$ case ( $k_x \propto m^{-1}$ ) corresponds to the energy spectrum of the Coulomb problem ( $E_n \propto -1/n^2$ ). The authors suggest these optical revivals serve as a classical emulator for Rydberg wave-packet dynamics, where the exact arithmetic limit replaces the approximate Taylor expansions typically used in quantum wave-packet studies.

The authors note that while the revival distance grows exponentially with bandwidth ( $z_{rev} \sim e^{rM}$ ), limiting the accessible $M$ in free space, the construction is platform-independent. They suggest dispersive fiber optics could extend the accessible range, allowing the LCM of integers up to $\sim 10^4$ to be encoded in propagation length. The work reframes the finite spectral cutoff not as a truncation error, but as a programmable arithmetic register.

Technical Summary: Chebyshev Self-Imaging from Inverse-Sampled Angular Spectra

More like this