🔭 astrophysics

Wigner function shapelets: Symplectic representation of astronomical images

This paper introduces Wigner function shapelets (WFSs), a novel symplectic framework that extends traditional shapelet analysis to four-dimensional phase space using Laguerre-Gaussian modes and $\mathrm{Sp}(4,\mathbb{R})$ group theory to provide a symmetry-preserving, quantum-information-based representation of astronomical images with inherent resolution limits and sensitivity to coherent morphological structures.

Original authors: Shun Arai

Published 2026-07-10

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Shun Arai

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

Imagine you are looking at a galaxy through a telescope. Usually, astronomers take a picture and measure its shape by looking at the light in two separate ways: either by looking at the picture itself (configuration space) or by looking at the pattern of light waves (Fourier space). It's like trying to understand a song by only listening to the lyrics or only listening to the sheet music, but never hearing them together.

In this paper, Shun Arai suggests a new way to look at these cosmic pictures. He introduces a tool called Wigner function shapelets (WFSs). Think of this as a magical "phase-space" lens that lets you see the galaxy's image and its wave patterns simultaneously in a single, four-dimensional view.

The "Quantum" Camera

To understand this, imagine the galaxy isn't just a static picture, but a living, breathing object that exists in a special "phase space." In this space, every point has a location (where the light is) and a momentum (how the light is moving).

The paper proposes that we can describe the galaxy's shape using a mathematical object called the Wigner function. This function is like a "quasi-probability map." Unlike a normal map that only shows positive numbers (like how much light is there), this map can show negative numbers and wiggles. These wiggles are crucial because they hold secret information about how the light waves interfere with each other—information that gets lost if you just look at the brightness alone.

The Shape-Shifting Blocks

The authors build their new tool using "shapelets," which are like a set of Lego blocks made from Laguerre-Gaussian modes. You can think of these modes as the fundamental "notes" or "vibrations" that a galaxy image can be made of.

The Old Way: Previous methods used these blocks to build the image either in the "picture" world or the "wave" world, but not both at once.
The New Way (WFS): The paper constructs these blocks directly in the four-dimensional phase space. This allows the blocks to preserve a special symmetry called symplectic geometry.

Here is the cool part: The paper argues that the way gravity bends light (gravitational lensing) acts on these galaxy shapes exactly like quantum mechanics acts on particles. When a galaxy is stretched by the gravity of a massive cluster, it's like a quantum particle jumping between energy levels. The WFS tool is designed to catch these "jumps" perfectly because it respects the same mathematical rules (the SU(2) group) that govern both quantum particles and galaxy shapes.

The "Hopf Torus" Map

To make sense of this four-dimensional mess, the authors use a clever trick called Hopf fibration. Imagine the galaxy's data is wrapped around a donut (a torus). The paper shows that you can flatten this donut onto a simple two-dimensional map defined by two numbers: $Q_0$ (which represents the total energy or size of the shape) and $Q_2$ (which represents the spin or twist of the shape).

On this map, the galaxy's shape isn't just a blob; it has a "band structure" with specific patterns of positive and negative waves. The authors define a new observable, $W_{k\ell}(Q_0, Q_2)$ , which acts like a dial that lets you tune into specific "twists" or "spins" of the galaxy.

What This Tool Can Do (and Can't Do)

The paper suggests that this method is incredibly powerful for several reasons, but it is careful to note what is proven and what is still being worked out:

It Keeps All the Secrets: Because the Wigner function keeps both the picture and the wave info together, it preserves information that other methods throw away.
It Handles Noise Like a Quantum Computer: The paper treats noise and blurring (like the atmosphere making stars twinkle) as "quantum channel operations." This is a fancy way of saying they can use math to describe how the image gets messed up and potentially fix it, similar to how quantum computers correct errors.
It Spies on the Universe's Secrets:
- Weak Lensing: The authors show how to calculate exactly how much a galaxy's shape changes due to cosmic shear (the stretching of space). They derive formulas that link the galaxy's "spin" modes to the shear.
- Parity Violation: They suggest a way to look for "parity violation" (a breaking of left-right symmetry in the universe) by checking if the galaxy shapes have a specific kind of "handedness" that shouldn't exist in a normal universe. This is done by looking at the correlation between two different galaxy images.

The "What Ifs" and Limits

The paper is very clear about what it is not yet.

It is not a finished product for daily use: The authors state that while the math is solid and the formulas are derived, they have not yet tested this on real telescope data or complex simulations to see if it works better than current methods in practice. They explicitly say they will "validate our methodology in practical image analysis" in the future.
It doesn't solve the blurring problem instantly: While the math describes how to handle blurring (like the Point Spread Function or PSF), the paper notes that you still need to choose the right "resolution" parameter, called $\lambda$ . If you pick $\lambda$ too small, you just see noise; if it's too big, you miss the details. The paper suggests a way to pick $\lambda$ based on the telescope's limits (like the size of the telescope mirror or the "seeing" from the atmosphere), estimating that for a diffraction-limited image, $\lambda$ is roughly 1.37 (a dimensionless number).
It's not a magic eraser: The paper admits that information loss due to noise and blurring is fundamental. You can't get back information that was erased, but this method helps you see exactly what is lost and what can be recovered.

The Bottom Line

Shun Arai's paper is a theoretical blueprint. It builds a bridge between the world of quantum physics and the world of astronomy, suggesting that if we look at galaxy images through the lens of "phase space" and "quantum channels," we might find new, sharper ways to measure the universe's shape and test the laws of physics.

The authors are confident in the math—they have derived the formulas for how these shapes respond to gravity and how they handle noise. However, they are cautious about the results: this is a new framework waiting to be tested on real data. It's like inventing a new type of microscope; the theory says it should see things no one else can, but we haven't looked through it yet to see the actual cells. The paper aims to provide the "blueprint" for this new tool, hoping that future astronomers will use it to unlock secrets about dark matter, cosmic expansion, and the fundamental nature of reality.

Technical Summary: Wigner Function Shapelets: Symplectic Representation of Astronomical Images

Problem Statement
Conventional shapelet formalisms for analyzing astronomical images operate separately in either configuration space (using Hermite–Gaussian modes) or Fourier space (using Laguerre–Gaussian modes). While effective, these approaches treat the image as a function of position or frequency alone, potentially losing the joint information encoded in the phase space. Furthermore, standard analyses often struggle to simultaneously preserve symplectic geometry, handle systematics (like Point Spread Function blurring) as algebraic operations, and capture topological or non-local correlations inherent in galaxy morphology. The paper addresses the need for a unified framework that represents astronomical images directly in a four-dimensional phase space, preserving the symplectic structure and enabling a more robust decomposition of image information.

Methodology
The authors introduce Wigner Function Shapelets (WFSs), a basis set designed to expand the Wigner function of an astronomical image. The methodology proceeds through the following steps:

Phase Space Representation: The image field $\psi(\theta)$ is mapped to a Wigner function $W(\theta, p)$ in a four-dimensional phase space defined by the transverse coordinate $\theta$ and conjugate momentum $p$ . This space is equipped with a symplectic structure $Sp(4, \mathbb{R})$ . The Wigner function is treated as a quasiprobability distribution that encodes full image information, including interference patterns and coherence, which are marginalized out in standard intensity or power spectrum analyses.
Basis Construction: The WFSs are constructed as the cross-Wigner functions of Laguerre–Gaussian (LG) modes. These modes form an orthogonal and complete basis for the Hilbert-Schmidt space of the image. The expansion coefficients are derived using the Moyal/Plancherel identity.
Symplectic and Group-Theoretic Decomposition: The authors utilize the $SU(2)$ irreducible representations of the two-mode oscillator to decompose the phase space. By employing the Hopf fibration, the four-dimensional phase space is mapped to a two-dimensional space of invariants $(Q_0, Q_2)$ (representing harmonic energy and axial angular momentum) and a two-torus of angles $(\phi_u, \phi_v)$ .
Observable Definition: A scalar function $W_{k\ell}(Q_0, Q_2)$ is introduced, obtained by averaging the Wigner function over the Hopf torus with specific winding numbers $(k, \ell)$ . This reduces the four-dimensional data to a "band structure" on the $(Q_0, Q_2)$ plane, where the indices $k$ and $\ell$ correspond to morphology and rotation modes, respectively.
Systematics and Noise Modeling: Systematics (such as PSF blurring) and noise are modeled as Completely Positive Trace-Preserving (CPTP) maps (quantum channels). The paper distinguishes between deterministic PSFs (coherent blurring) and stochastic PSFs (random displacements), analyzing their distinct effects on the morphology and rotation charges of the WFS.

Key Contributions and Results

Symplectic Covariance: The WFS formalism provides a representation that is covariant under linear canonical (symplectic) transformations, including translations, rotations, shears, and squeezes. This allows for a natural description of image deformations induced by gravitational lensing.
Analytic Closed Forms: The paper derives closed analytic expressions for the Wigner function of LG modes and their cross-terms. These expressions are factorized into products of one-mode kernels, explicitly showing dependence on the Hopf torus variables.
Selection Rules and Symmetries: The authors establish selection rules based on statistical isotropy and parity symmetry.
- Statistical Isotropy: Enforces diagonal selection rules ( $k=\ell$ ) for the Wigner-torus coefficients in an ensemble average.
- Parity Symmetry: Requires that modes with odd $k+\ell$ vanish. Consequently, non-zero off-diagonal or odd-sum modes in an ensemble can serve as sensitive probes for parity violation (e.g., cosmic birefringence).
- Chiral Symmetry: Links the symmetry of the coefficients under the exchange of indices and the sign flip of $Q_2$ .
Response to Weak Lensing: The paper derives the linear response of the WFS to weak gravitational lensing (shear and flexion). The response is expressed as a coupling between the WFS coefficients and the lensing potential, governed by specific selection rules ( $\Delta s = \pm 2$ for shear, $\Delta s = \pm 1, \pm 3$ for flexion).
Connection to Existing Methods: The work establishes analytic correspondences between WFS, BiPoSHs (Bipolar Spherical Harmonics), and the Fourier Power Function Shapelets (FPFS), demonstrating that WFS generalizes these existing methods within a phase-space framework.
Information Recovery: Using quantum information theory concepts, the paper defines the entanglement fidelity of image recovery. It shows that the WFS basis diagonalizes the fidelity kernel, allowing for the identification of which morphology and rotation modes can be optimally recovered from noisy, blurred data.

Significance and Claims
The paper claims to provide a formal baseline for applying quantum information theory to astronomical image analysis. Its primary significance lies in:

Unified Framework: It unifies configuration and Fourier space analyses into a single symplectic phase-space representation, capturing correlations inaccessible to either domain alone.
Topological Sensitivity: By utilizing Hopf fibrations, the method exposes topological and non-local information (phase coherence) in galaxy shapes, which is crucial for detecting subtle cosmological imprints like parity violation.
Algebraic Treatment of Systematics: It reframes image degradation (PSF, noise) as quantum channel operations, offering a rigorous algebraic approach to deconvolution and error propagation that respects the underlying symplectic geometry.
Foundation for Future Analysis: The authors position WFS not as a completed solution for immediate data processing, but as a theoretical foundation. They acknowledge that practical implementation (e.g., optimal recovery algorithms) requires further development, drawing parallels to the rapid advancement of quantum image processing in other fields. The paper aims to facilitate the validation and benchmarking of advanced analysis methods for upcoming high-resolution surveys (e.g., Euclid, Roman, Rubin Observatory) by providing a symmetry-preserving reference model.