Path integral formulation of finite-dimensional quantum mechanics in discrete phase space

This paper develops an exact path integral formulation for finite-dimensional quantum systems in discrete phase space, deriving a sum-over-paths propagator that captures full entanglement dynamics through coherent contributions from all fluctuation sectors, thereby overcoming the limitations of single-sector approximations and providing a framework for semiclassical simulation and non-classicality characterization.

The Gist

Imagine you are trying to predict the future of a tiny, magical particle. In the world of big, everyday objects (like baseballs or planets), we use simple rules: if you know where a ball is and how fast it's moving, you can draw a single, smooth line showing exactly where it will be next. This is classical physics.

But in the quantum world (the world of atoms and electrons), things are weird. Particles don't just have one position; they exist in a "cloud" of possibilities. To track this cloud, scientists use a special map called the Wigner function. Think of this map not as a picture of where the particle is, but a weather map showing the probability of it raining (being found) in different places.

The Problem: The Map is Broken for Small Worlds

For a long time, scientists had a perfect way to draw this weather map for continuous things (like a ball rolling on a smooth table). They could use a "Path Integral" method. Imagine this as a way of calculating the future by summing up every single possible path the particle could take, not just the straight line. Some paths cancel each other out, and others add up to create the final result.

However, this method was missing for finite-dimensional systems—the tiny, discrete building blocks of quantum computers (like qubits and qutrits). These systems are like a chessboard with a fixed number of squares, rather than an infinite smooth floor. The old "smooth path" math didn't work on a grid.

The Solution: A New Way to Walk the Grid

The authors of this paper, Leonardo Pachón and Andrés Gómez, have built a new mathematical bridge. They created a Path Integral for the Quantum Grid.

Here is the analogy:

• The Old Way (Continuous): Imagine walking through a foggy forest. You can take any step, any angle. The math sums up all these infinite steps to predict where you end up.
• The New Way (Discrete): Imagine you are on a giant, glowing chessboard (the "discrete phase space"). You can only land on the squares. You can't step between them. The authors figured out how to sum up every possible route you could take across these specific squares to predict the future of the quantum system.

The Secret Ingredient: The "Ghost" Paths

The most fascinating part of their discovery is about fluctuations.

In their new formula, to get the correct answer, you have to sum up two types of paths:

1. The Main Path: The "average" route the particle seems to take.
2. The Ghost Paths (Fluctuations): These are the weird, jiggly, "what-if" paths that deviate from the average.

The authors discovered a shocking truth: If you ignore the Ghost Paths, the whole thing breaks.

• The Trap: If you only look at the Main Path (a method used in some current simulations called DTWA), you get a result that looks okay for a split second, but then it becomes nonsense. It predicts things that aren't real (like negative probabilities in a way that doesn't make sense) or it just says "the particle is equally likely to be anywhere," which is wrong.
• The Fix: The "Ghost Paths" are the secret sauce. They interfere with each other like waves in a pond. When you add them all up together, they cancel out the nonsense and reveal the true, complex behavior of the particle, including entanglement (where two particles become linked and share a fate).

Why This Matters: The Magic of Quantum Computers

The paper uses a specific example of a "Qutrit" (a quantum coin with three sides instead of two) to prove their point.

1. The Test: They watched two quantum coins interact.
2. The Result: They showed that to understand how these coins become "entangled" (linked together), you must include all the weird, fluctuating paths. If you try to simplify the math by ignoring the fluctuations, you lose the magic. You lose the ability to predict how the quantum computer will actually work.

The "Pseudo-Classical" Exception

There is one special case where the Ghost Paths disappear. If the quantum system is very simple (linear) and the timing is perfectly synchronized with the grid (like a clock ticking in perfect rhythm with the squares), the particle behaves like a classical object. It just hops from square to square in a straight line. This is the "boring" part of quantum mechanics where it acts like normal physics.

But the moment you introduce complexity or change the timing, the Ghost Paths wake up, and the system becomes truly quantum.

In a Nutshell

This paper is like finding the instruction manual for navigating a quantum video game that only exists on a pixelated screen.

• Before: We tried to use rules for smooth, real-world physics on a pixelated screen, and it glitched.
• Now: We have a new rulebook that says, "To predict the future on this grid, you must sum up every possible zig-zag, jitter, and ghost path."
• The Lesson: You cannot simplify quantum mechanics by ignoring the "noise" (the fluctuations). That noise is actually the signal. It's what makes quantum computers powerful and what creates the "spooky" connections between particles.

This new tool will help scientists build better simulations for future quantum computers and understand exactly why quantum systems are so different from the world we see with our eyes.

Technical Summary

1. Problem Statement

While the phase-space formulation of quantum mechanics (specifically the Wigner-Weyl-Moyal framework) is well-established for continuous-variable systems, its extension to finite-dimensional Hilbert spaces (such as qubits and qudits) lacks a systematic dynamical treatment.

• The Gap: Existing discrete phase-space methods (e.g., the Discrete Truncated Wigner Approximation, DTWA) rely on sampling classical trajectories but lack an exact path integral representation with an explicit action functional.
• The Challenge: Constructing a path integral for discrete systems requires defining a "sum-over-paths" on a discrete toroidal lattice ($Z_d \times Z_d$) that correctly reproduces quantum dynamics, including entanglement and non-classical interference, without relying on a continuum limit that destroys the discrete algebraic structure.
• Motivation: A rigorous path integral is essential for developing semiclassical simulation methods for many-body spin systems and for understanding the resource theory of quantum computation (specifically the role of Wigner negativity and "magic" states).

2. Methodology

The authors construct the theory entirely within a discrete phase space defined over the finite field $Z_d$, where $d$ is an odd prime. The methodology proceeds in four logical stages:

A. Algebraic Foundation

• Discrete Fourier Transform (DFT): They define position-like ($\hat{x}$) and momentum-like ($\hat{p}$) operators related by a DFT.
• Generalized Displacement Operators: Using clock ($\hat{Z}$) and shift ($\hat{X}$) operators, they construct the displacement operators $\hat{D}(k, j) = \omega^{-kj/2} \hat{Z}^k \hat{X}^j$, which form an orthogonal basis for the operator space.
• Discrete Weyl Transform: They define the Weyl symbol of an operator $\hat{A}$ as $\tilde{A}_W(k, j) = \text{Tr}[\hat{A}\hat{D}(k, j)]$. The inverse transform maps these symbols back to operators.
• Discrete Wigner Function: Defined as the inverse 2D DFT of the Weyl symbol of the density matrix, ensuring real-valued marginals that reproduce quantum probabilities.

B. Derivation of the Evolution Kernel

• Twisted Convolution: The authors derive the composition rule for the product of operators in the Weyl representation (the discrete analog of the Moyal star product).
• Exact Propagator: By applying the composition rule to the unitary evolution $\hat{\rho}(t) = \hat{U}(t)\hat{\rho}(0)\hat{U}^\dagger(t)$, they derive an exact closed-form expression for the evolution kernel $G_W$ in the discrete phase space. This kernel relates the Wigner function at time $t$ to its initial state via a sum over an auxiliary fluctuation variable $\tilde{\mu}$.

C. Path Integral Construction

• Time Slicing: The time interval is partitioned into $N$ steps. The composition law of the kernel is iterated to create a sum over intermediate phase-space points.
• Short-Time Approximation: For small time steps $\tau$, the evolution operator's symbol is approximated by the exponential of the Hamiltonian symbol ($U_W \approx e^{-i H_W \tau}$).
• Sum-over-Paths: By iterating the short-time kernel, the authors derive a path integral representation where the propagator is a sum over piecewise-constant paths $\gamma$ (the trajectory) and fluctuation paths $\xi$ (analogous to the $\tilde{r}$ variable in continuous Marinov integrals).

D. The Discrete Action

The resulting action functional $S_W[\gamma, \xi, t]$ is the discrete counterpart of Marinov's functional:
$$S_W = -\frac{4\pi}{d} \sum \Delta \gamma_i \wedge \xi_i + \frac{1}{\hbar} \sum [H_W(\gamma_i + \xi_i) - H_W(\gamma_{i-1} - \xi_i)] \tau$$
Here, $\wedge$ denotes the discrete symplectic product on $Z_d \times Z_d$.

3. Key Contributions

1. Exact Path Integral Formulation: The paper provides the first exact path integral representation for the dynamics of finite-dimensional quantum systems in discrete phase space. It expresses the propagator as a sum over paths weighted by a discrete action (Eq. 29-30).
2. Identification of Fluctuation Sectors: The authors rigorously demonstrate that the $\tilde{\mu} = 0$ sector (the "mean-field" or boundary term) is insufficient.
   • At finite time steps, this sector yields a non-real kernel.
   • In the continuum limit ($N \to \infty$), it collapses to a trivial uniform kernel that carries no dynamical information.
   • Conclusion: Genuine quantum dynamics, including entanglement, arise strictly from the coherent interference of all non-zero fluctuation sectors ($\tilde{\mu} \neq 0$).
3. Pseudo-Classical Regime: They identify a specific regime where the path integral simplifies to a deterministic shift. For Hamiltonians linear in phase-space coordinates and at times strictly commensurate with the lattice (Clifford-covariant regime), the fluctuation sum collapses, realizing a discrete analog of classical Hamiltonian flow.
4. Application to Entanglement: The formalism is applied to interacting qutrits, showing that the exact time evolution of the linear entropy (a measure of entanglement) requires the full path integral sum. Truncating to the mean-field sector fails to reproduce entanglement generation.

4. Results and Verification

• Single Qutrit ($d=3$): The authors verified the formalism for a single qutrit under a diagonal Hamiltonian. The path integral correctly reproduced the time-evolved Wigner function, including the emergence of Wigner negativity when the system evolved into a non-stabilizer state.
• Two Interacting Qutrits:
  • Entanglement Dynamics: For a system with Hamiltonian $\hat{H} = \hbar\chi (\hat{x}_1 \otimes \hat{x}_2)$, the authors derived a closed-form expression for the purity $\text{Tr}[\hat{\rho}_1^2]$ and linear entropy.
  • Numerical Confirmation: They numerically verified that the full path integral (summing over all $d^{2n}$ fluctuation vectors) reproduces the exact entanglement dynamics to machine precision.
  • Failure of Truncation: They demonstrated that retaining only the $\tilde{\mu}=0$ term fails to capture entanglement, confirming that quantum correlations in discrete phase space are inherently non-local in the fluctuation variables.
• Code Availability: All quantitative claims are supported by a public Python repository that reproduces the results for $d \in \{3, 5, 7, 11\}$.

5. Significance and Outlook

• Semiclassical Simulation: This work provides a rigorous foundation for improving the Discrete Truncated Wigner Approximation (DTWA). By incorporating higher-order quantum corrections via the fluctuation sectors, it offers a pathway to simulate many-body spin systems in regimes where exact diagonalization is impossible.
• Quantum Resource Theory: The formulation links the generation of Wigner negativity (a resource for quantum advantage) directly to the interference of fluctuation paths in the path integral. This offers a dynamical perspective on why stabilizer circuits are classically simulable (no negativity, deterministic shifts) while non-stabilizer circuits are not.
• Future Directions: The authors suggest extending the framework to:
  • Qubits ($d=2$): Requiring a modification of the displacement operators to handle the lack of a multiplicative inverse of 2.
  • Prime-Power Dimensions: Using Galois field constructions.
  • Open Systems: Describing decoherence in finite-dimensional systems.
  • Monte Carlo Methods: Developing efficient sampling techniques for the fluctuation paths to overcome the "sign problem" inherent in the oscillatory weights.

In summary, this paper successfully bridges the gap between continuous path integrals and discrete quantum mechanics, providing a powerful tool for analyzing and simulating the dynamics of finite-dimensional quantum systems, particularly regarding entanglement and non-classicality.