Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2604.20776

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

Leonardo A. Pachon, Andres D. Gomez

quant-phmath-phSubmitted Apr 22, 2026• Announced Apr 23, 2026• 0 citations• Revised Jun 28, 2026 (v2)
v1 · Apr 22, 2026v2 · Jun 28, 2026
arXivHTMLPDF

Abstract

We develop a path integral representation for the dynamics of quantum systems with a finite-dimensional Hilbert space, formulated entirely within a discrete phase space. Starting from the discrete Wigner function on ${\mathbb{Z}_d} \times {\mathbb{Z}_d}$ ($d$ an odd prime) and the associated Weyl transform built from generalized displacement operators, we derive an exact kernel that propagates the discrete Wigner function in time and, by iterating its composition law through a short-time approximation, obtain a sum-over-paths expression weighted by a discrete phase-space action -- the natural finite-dimensional counterpart of Marinov's functional. For Hamiltonians linear in the phase-space coordinates and at times strictly commensurate with the lattice, the fluctuation sum collapses to a deterministic shift, realizing the discrete analog of classical Hamiltonian flow. Applying the formalism to one and to two interacting qutrits ($d=3$), we show that the full entanglement dynamics -- captured by a closed-form linear entropy valid for all times -- requires the coherent contribution of all fluctuation sectors; the boundary-term (mean-field) sector alone fails to reproduce it. For a non-stabilizer Hamiltonian, where the short-time kernel is only approximate, the time-sliced path integral converges to the exact dynamics, including the dynamical generation of Wigner negativity. We discuss implications for the semiclassical simulation of many-body spin systems and for the characterization of non-classicality through Wigner negativity.

Comments
12 pages, numerical examples included

Related papers

Quantum circuit dynamics via path integrals: Is there a classical action for discrete-time paths?

Mark D Penney, Dax Enshan Koh, Robert W Spekkens

It is straightforward to compute the transition amplitudes of a quantum circuit using the sum-over-paths methodology when the gates in the circuit are balanced, where a balanced gate is one for which all non-zero transition amplitudes are of equal magnitude. Here we consider the question of whether, for such circuits, the relative phases of different discrete-time paths through the configuration space can be defined in terms of a classical action, as they are for continuous-time paths. We show how to do so for certain kinds of quantum circuits, namely, Clifford circuits where the elementary systems are continuous-variable systems or discrete systems of odd-prime dimension. These types of circuit are distinguished by having phase-space representations that serve to define their classical counterparts. For discrete systems, the phase-space coordinates are also discrete variables. We show that for each gate in the generating set, one can associate a symplectomorphism on the phase-space and to each of these one can associate a generating function, defined on two copies of the configuration space. For discrete systems, the latter association is achieved using tools from algebraic geometry. Finally, we show that if the action functional for a discrete-time path through a sequence of gates is defined using the sum of the corresponding generating functions, then it yields the correct relative phases for the path-sum expression. These results are likely to be relevant for quantizing physical theories where time is fundamentally discrete, characterizing the classical limit of discrete-time quantum dynamics, and proving complexity results for quantum circuits.

arXiv:1604.07452•New Journal of Physics•Feb 21, 2017•10 citations

Path integrals and the discrete Weyl representation

Wayne Nicholas Polyzou

This paper uses Schwinger's representation of quantum systems of a finite number of degrees of freedom to formulate discrete path integral representations of the dynamics. The representation starts with an observable with a finite number of distinct eigenvalues. Using that observable Schwinger constructs complementary pairs of finite dimensional unitary operators that are finite dimensional analogs of the irreducible Weyl algebra. These can be decomposed into products of irreducible sub-algebras that reduce to q-bit gates for the case of $2^N$ degrees of freedom. In the limit of a large number of discrete degrees of freedom these representations can be used to model quantum systems with continuous degrees of freedom. This limit recovers the continuum Weyl algebra. Path integrals can be formulated in the discrete case. They have the advantage that the number of paths for a finite number of time slices is finite. In this work the path integral is interpreted as the expectation value of a potential functional with respect to a complex probability distribution on the space of paths. An application to scattering from a short-range potential is given. Multi-resolution wavelet bases are used exactly represent a local $\phi^4(x)$ quantum field theory as a theory with an infinite number of discrete local modes. The discrete path integral is illustrated by computing time evolution when this theory is truncated to two coupled modes.

•Aug 27, 2021•0 citations

Quantum Phase Space Representations and Their Negativities

Ninnat Dangniam

A classical simulation scheme of quantum computation given a restricted set of states and measurements may be---occasionally, but only occasionally---interpreted naturally as a statistical simulation of positive quasi-probability distributions on a phase space. In this dissertation, we explore phase space representations for finite-dimensional quantum systems and their negativities beyond the usual analogues of the Wigner function. The first line of study focuses on a characterization tool for valid quasi-probability distributions of (possibly mixed) quantum states. A quantum generalization of Bochner's theorem from classical probability theory simultaneously characterizes both the set of valid Wigner functions and the subset of positive ones. We extend this theorem to discrete phase spaces based on projective representations of abelian groups, including the discrete Wigner functions of Gross (J. Math. Phys. 47, 122107 (2006)) and Gibbons et al. (Phys. Rev. A 70, 062101 (2004)). The rest of the dissertation is then dedicated to phase space representations covariant with respect to a general Lie group. This means, in particular, that operators that lie in the same orbit of the group have the same amount of negativity in their quasi-probabilities---a natural requirement if one would like to interpret negativity as non-classicality, and the group action as a classical dynamics. We show that the construction due to Brif and Mann (Phys. Rev. A 59, 971 (1999)) of group covariant phase space representations satisfying the so-called Stratonovich-Weyl correspondence is essentially unique, given the group and an input state in the relevant representation. As an application, we construct quasi-probability representations on a compact phase space of fermionic Gaussian states tailored to a classically simulatable problem of fermionic linear optics and find that, in the first non-trivial case of four fermionic modes, the quasi-probabilities of the Gaussian states exhibit a considerable amount of negativity.

•Jan 1, 2018•0 citations

Path integrals, complex probabilities and the discrete Weyl representation

Wayne Nicholas Polyzou

Abstract A discrete formulation of the real-time path integral as the expectation value of a functional of paths with respect to a complex probability on a sample space of discrete valued paths is explored. The formulation in terms of complex probabilities is motivated by a recent reinterpretation of the real-time path integral as the expectation value of a potential functional with respect to a complex probability distribution on cylinder sets of paths. The discrete formulation in this work is based on a discrete version of the Weyl algebra that can be applied to any observable with a finite number of outcomes. The origin of the complex probability in this work is the completeness relation. In the discrete formulation the complex probability exactly factors into products of conditional probabilities and exact unitarity is maintained at each level of approximation. The approximation of infinite dimensional quantum systems by discrete systems is discussed. The method is illustrated by applying it to scattering theory and quantum field theory. The implications of these applications for quantum computing is discussed.

arXiv:2108.12494•Journal of Physics A Mathematical and Theoretical•Jun 6, 2024•4 citations

Stationary Phase Method in Discrete Wigner Functions and Classical Simulation of Quantum Circuits

Lucas Kocia, Peter J. Love

One of the lowest-order corrections to Gaussian quantum mechanics in infinite-dimensional Hilbert spaces are Airy functions: a uniformization of the stationary phase method applied in the path integral perspective. We introduce a "periodized stationary phase method" to discrete Wigner functions of systems with odd prime dimension and show that the π8 gate is the discrete analog of the Airy function. We then establish a relationship between the stabilizer rank of states and the number of quadratic Gauss sums necessary in the periodized stationary phase method. This allows us to develop a classical strong simulation of a single qutrit marginal on t qutrit π8 gates that are followed by Clifford evolution, and show that this only requires 3t2+1 quadratic Gauss sums. This outperforms the best alternative qutrit algorithm (based on Wigner negativity and scaling as ∼30.8t for 10−2 precision) for any number of π8 gates to full precision.

arXiv:1810.03622•Quantum•Jul 5, 2021•8 citations

From the discrete Weyl–Wigner formalism for symmetric ordering to a number–phase Wigner function

Maciej Przanowski, Jaromir Tosiek

The general Weyl–Wigner formalism in finite dimensional phase spaces is investigated. Then this formalism is specified to the case of symmetric ordering of operators in an odd-dimensional Hilbert space. A respective Wigner function on the discrete phase space is found and the limit, when the dimension of Hilbert space tends to infinity, is considered. It is shown that this limit gives the number–phase Wigner function in quantum optics. Analogous results for the “almost” symmetric ordering in an even-dimensional Hilbert space are obtained. Relations between the discrete Wigner functions introduced in our paper and some other discrete Wigner functions appearing in literature are studied.

arXiv:1609.03751•Journal of Mathematical Physics•Oct 1, 2017•12 citations
We develop a path integral representation for the dynamics of quantum systems with a finite-dimensional Hilbert space, formulated entirely within a discrete phase space. Starting from the discrete Wigner function on Zd×Zd{\mathbb{Z}_d} \times {\mathbb{Z}_d}Zd​×Zd​ (ddd an odd prime) and the associated Weyl transform built from generalized displacement operators, we derive an exact kernel that propagates the discrete Wigner function in time and, by iterating its composition law through a short-time approximation, obtain a sum-over-paths expression weighted by a discrete phase-space action -- the natural finite-dimensional counterpart of Marinov's functional. For Hamiltonians linear in the phase-space coordinates and at times strictly commensurate with the lattice, the fluctuation sum collapses to a deterministic shift, realizing the discrete analog of classical Hamiltonian flow. Applying the formalism to one and to two interacting qutrits (d=3d=3d=3), we show that the full entanglement dynamics -- captured by a closed-form linear entropy valid for all times -- requires the coherent contribution of all fluctuation sectors; the boundary-term (mean-field) sector alone fails to reproduce it. For a non-stabilizer Hamiltonian, where the short-time kernel is only approximate, the time-sliced path integral converges to the exact dynamics, including the dynamical generation of Wigner negativity. We discuss implications for the semiclassical simulation of many-body spin systems and for the characterization of non-classicality through Wigner negativity.
It is straightforward to compute the transition amplitudes of a quantum circuit using the sum-over-paths methodology when the gates in the circuit are balanced, where a balanced gate is one for which all non-zero transition amplitudes are of equal magnitude. Here we consider the question of whether, for such circuits, the relative phases of different discrete-time paths through the configuration space can be defined in terms of a classical action, as they are for continuous-time paths. We show how to do so for certain kinds of quantum circuits, namely, Clifford circuits where the elementary systems are continuous-variable systems or discrete systems of odd-prime dimension. These types of circuit are distinguished by having phase-space representations that serve to define their classical counterparts. For discrete systems, the phase-space coordinates are also discrete variables. We show that for each gate in the generating set, one can associate a symplectomorphism on the phase-space and to each of these one can associate a generating function, defined on two copies of the configuration space. For discrete systems, the latter association is achieved using tools from algebraic geometry. Finally, we show that if the action functional for a discrete-time path through a sequence of gates is defined using the sum of the corresponding generating functions, then it yields the correct relative phases for the path-sum expression. These results are likely to be relevant for quantizing physical theories where time is fundamentally discrete, characterizing the classical limit of discrete-time quantum dynamics, and proving complexity results for quantum circuits.
This paper uses Schwinger's representation of quantum systems of a finite number of degrees of freedom to formulate discrete path integral representations of the dynamics. The representation starts with an observable with a finite number of distinct eigenvalues. Using that observable Schwinger constructs complementary pairs of finite dimensional unitary operators that are finite dimensional analogs of the irreducible Weyl algebra. These can be decomposed into products of irreducible sub-algebras that reduce to q-bit gates for the case of 2N2^N2N degrees of freedom. In the limit of a large number of discrete degrees of freedom these representations can be used to model quantum systems with continuous degrees of freedom. This limit recovers the continuum Weyl algebra. Path integrals can be formulated in the discrete case. They have the advantage that the number of paths for a finite number of time slices is finite. In this work the path integral is interpreted as the expectation value of a potential functional with respect to a complex probability distribution on the space of paths. An application to scattering from a short-range potential is given. Multi-resolution wavelet bases are used exactly represent a local ϕ4(x)\phi^4(x)ϕ4(x) quantum field theory as a theory with an infinite number of discrete local modes. The discrete path integral is illustrated by computing time evolution when this theory is truncated to two coupled modes.
A classical simulation scheme of quantum computation given a restricted set of states and measurements may be---occasionally, but only occasionally---interpreted naturally as a statistical simulation of positive quasi-probability distributions on a phase space. In this dissertation, we explore phase space representations for finite-dimensional quantum systems and their negativities beyond the usual analogues of the Wigner function. The first line of study focuses on a characterization tool for valid quasi-probability distributions of (possibly mixed) quantum states. A quantum generalization of Bochner's theorem from classical probability theory simultaneously characterizes both the set of valid Wigner functions and the subset of positive ones. We extend this theorem to discrete phase spaces based on projective representations of abelian groups, including the discrete Wigner functions of Gross (J. Math. Phys. 47, 122107 (2006)) and Gibbons et al. (Phys. Rev. A 70, 062101 (2004)). The rest of the dissertation is then dedicated to phase space representations covariant with respect to a general Lie group. This means, in particular, that operators that lie in the same orbit of the group have the same amount of negativity in their quasi-probabilities---a natural requirement if one would like to interpret negativity as non-classicality, and the group action as a classical dynamics. We show that the construction due to Brif and Mann (Phys. Rev. A 59, 971 (1999)) of group covariant phase space representations satisfying the so-called Stratonovich-Weyl correspondence is essentially unique, given the group and an input state in the relevant representation. As an application, we construct quasi-probability representations on a compact phase space of fermionic Gaussian states tailored to a classically simulatable problem of fermionic linear optics and find that, in the first non-trivial case of four fermionic modes, the quasi-probabilities of the Gaussian states exhibit a considerable amount of negativity.
Abstract A discrete formulation of the real-time path integral as the expectation value of a functional of paths with respect to a complex probability on a sample space of discrete valued paths is explored. The formulation in terms of complex probabilities is motivated by a recent reinterpretation of the real-time path integral as the expectation value of a potential functional with respect to a complex probability distribution on cylinder sets of paths. The discrete formulation in this work is based on a discrete version of the Weyl algebra that can be applied to any observable with a finite number of outcomes. The origin of the complex probability in this work is the completeness relation. In the discrete formulation the complex probability exactly factors into products of conditional probabilities and exact unitarity is maintained at each level of approximation. The approximation of infinite dimensional quantum systems by discrete systems is discussed. The method is illustrated by applying it to scattering theory and quantum field theory. The implications of these applications for quantum computing is discussed.
One of the lowest-order corrections to Gaussian quantum mechanics in infinite-dimensional Hilbert spaces are Airy functions: a uniformization of the stationary phase method applied in the path integral perspective. We introduce a "periodized stationary phase method" to discrete Wigner functions of systems with odd prime dimension and show that the π8 gate is the discrete analog of the Airy function. We then establish a relationship between the stabilizer rank of states and the number of quadratic Gauss sums necessary in the periodized stationary phase method. This allows us to develop a classical strong simulation of a single qutrit marginal on t qutrit π8 gates that are followed by Clifford evolution, and show that this only requires 3t2+1 quadratic Gauss sums. This outperforms the best alternative qutrit algorithm (based on Wigner negativity and scaling as ∼30.8t for 10−2 precision) for any number of π8 gates to full precision.
The general Weyl–Wigner formalism in finite dimensional phase spaces is investigated. Then this formalism is specified to the case of symmetric ordering of operators in an odd-dimensional Hilbert space. A respective Wigner function on the discrete phase space is found and the limit, when the dimension of Hilbert space tends to infinity, is considered. It is shown that this limit gives the number–phase Wigner function in quantum optics. Analogous results for the “almost” symmetric ordering in an even-dimensional Hilbert space are obtained. Relations between the discrete Wigner functions introduced in our paper and some other discrete Wigner functions appearing in literature are studied.