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:2108.12494 · Journal of Physics A Mathematical and Theoretical

Path integrals, complex probabilities and the discrete Weyl representation

Wayne Nicholas Polyzou

Published Jun 6, 2024• 4 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1088/1751-8121/ad5503

Related papers

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

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

Leonardo A. Pachon, Andres D. Gomez

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.

quant-pharXiv:2604.20776•Apr 23, 2026•0 citations

Scattering using real-time path integrals

Wayne Nicholas Polyzou, Ekaterina S. Nathanson

Background: Path integrals are a powerful tool for solving problems in quantum theory that are not amenable to a treatment by perturbation theory. Most path integral computations require an analytic continuation to imaginary time. While imaginary time treatments of scattering are possible, imaginary time is not a natural framework for treating scattering problems. More importantly, quantum algorithms for calculating path integrals require real-time evolution.Purpose: We test a recently introduced method for performing direct calculations of scattering observables using real-time path integrals in order to understand the challenges facing real-time path integral calculations of scattering observables.Method: The computations are based on a new interpretation of the path integral as the expectation value of a potential functional on cylinder sets of continuous paths with respect to a complex probability distribution. The method can in principle be applied to arbitrary short-range potentials.Results: The method is applied to compute matrix elements of M\o{}ller wave operators applied to narrow wave packets. These are used to calculate half-shell sharp-momentum transition matrix elements for one-dimensional potential scattering. The calculations for half-shell transition operator matrix elements converge to the numerical solution of the Lippmann-Schwinger equation.Conclusions: This work presents a proof in principle that scattering observables can be computed using real-time Feynman path integrals. While the computational method is not efficient, it can be improved. It provides a laboratory for studying quantum computational algorithms that are applicable to scattering problems.

arXiv:1712.00046•Physical Review C•Jun 4, 2020•2 citations

The Feynman path integral

Mark E. Tuckerman

Abstract The focus of Chapter 12 is the Feynman path integral formulations of quantum mechanics and quantum statistical mechanics. The connection between real-time quantum propagator and the density matrix is established via the Wick rotation, and a Trotter expansion is used to derive discrete path integral expressions for both. The continuous limit of the discrete path integral is derived to yield the functional integral formulation of the path integral. The construction of estimators and their averages for computing quantum observables is discussed. Next, the path integral is extended to incorporate Boltzmann, Fermi-Dirac, and Bose-Einstein statistics, and a useful recursive formula for the partition function is derived that accelerates convergence by building up N-body partition functions from those of smaller systems. Quantum free energies are discussed, and algorithms for the evaluation of path integrals using molecular dynamics or Monte Carlo are presented. The chapter concludes with a discussion of higher-order factorizations and ring-polymer contraction methods to reduce the computational overhead of path-integral calculations

•Aug 1, 2023•0 citations

Gaussian integrals. Algebraic preliminaries

Jean Zinn‐Justin

Abstract In this work, the perturbative aspects of quantum mechanics (QM) and quantum field theory (QFT), to a large extent, are studied with functional (path or field) integrals and functional techniques. This physics textbook thus begins with a discussion of algebraic properties of Gaussian measures, and Gaussian expectation values for a finite number of variables. The important role of Gaussian measures is not unrelated to the central limit theorem of probabilities, although the interesting physics is generally hidden in essential deviations from Gaussian distributions. A few algebraic identities about Gaussian expectation values, in particular Wick's theorem are recalled. Integrals over some type of formally complex conjugate variables, directly relevant for boson systems are defined. Fermion systems require the introduction of Grassmann or exterior algebras, and the corresponding generalization of the notions of differentiation and integration. Both for complex and Grassmann integrals, Gaussian integrals, and Gaussian expectation values are calculated, and generalized Wick's theorems proven. The concepts of generating functions and Legendre transformation are recalled.

•Apr 15, 2021•0 citations

Path integration with non-positive distributions and applications to the Schrödinger equation

Ekaterina S. Nathanson

In 1948, Richard Feynman published the first paper on his new approach to non-relativistic quantum mechanics. Before Feynman's work there were two mathematical formulations of quantum mechanics. Schrödinger's formulation was based on PDE (the Schrödinger equation) and states representation by wave functions, so it was in the framework of analysis and differential equations. The other formulation was Heisenberg's matrix algebra. Initially, they were thought to be competing. The proponents of one claimed that the other was “ wrong. ” Within a couple of years, John von Neumann had proved that they are equivalent.\nAlthough Feynman's theory was not fundamentally new, it nonetheless offered an entirely fresh and different perspective: via a precise formulation of Bohr's correspondence principle, it made quantum mechanics similar to classical mechanics in a precise sense. In addition, Feynman's approach made it possible to explain physical experiments, and, via diagrams, link them directly to computations. What resulted was a very powerful device for computing energies and scattering amplitudes - the famous Feynman's diagrams.\nIn his formulation, Feynman aimed at representing the solution to the non-relativistic Schrödinger equation in the form of an “ average ” over histories or paths of a particle. This solution is commonly known as the Feynman path integral. It plays an important role in the theory but appears as a postulate based on intuition coming from physics rather than a justified mathematical object. This is why Feynman's vision has caught the attention of many mathematicians as well as physicists.\nThe papers of Gelfand, Cameron, and Nelson are among the first, and more substantial, attempts to supply Feynman's theory with a rigorous mathematical foundation. These attempts were followed by many others, but unfortunately all of them were not quite satisfactory. The difficulty comes from a need to define a measure on an infinite-dimensional space of continuous functions that represent all possible paths of a particle. This Feynman's measure has to produce an integral with the properties requested by Feynman. In particular, the expression for the Feynman measure has to involve the non-absolutely integrable Fresnel integrands. The non-absolute integrability of the Fresnel integrands makes the measure fail to be positive and to have the countably additive property. Thus, a well-defined measure in the case of the Feynman path integral does not exist.\nExtensive research has been done on the methods of relating the Feynman path integral to the integral with respect to the Wiener measure. The method of analytic continuation in mass defines the Feynman path integral as a certain limit of the Wiener integrals. Unfortunately, this method can be used as definition for only almost all values of the mass parameter in the Schrödinger equation. For physicists, this is not a satisfactory result and needs to be improved.\nIn this work we examine those questions which originally led to the Feynman path integral. By now we know that Feynman's “ dream ” cannot be realized as a positive and countably additive measure on the path-space. Here, we offer a new way out by modifying Feynman's question, and thereby achieving a solution to the Schrödinger equation via a different kind of averages in the path-space.\nWe give our version of the question that Feynman “ should have asked ” in order to realize the elusive path integral. In our formulation, we get a Feynman path integral as a limit of linear functionals, as opposed to the more familiar inductive limits of positive measures, traditionally used for constructing the Wiener measure, and related Gaussian families. We adapt here an approach pioneered by Patrick Muldowney. In it, Muldowney suggested a Henstock integration technique in order to deal with the non-absolute integrability of the kind of Fresnel integrals which we need in our solution to Feynman's question. By applying Henstock's theory to Fresnel integrals, we construct a complex-valued “ probability distribution functions ” on the path-space. Then we use this “ probability ” distribution function to define the Feynman path integral as an inductive limit. This establishes a mathematically rigorous Feynman limit, and at the same time, preserves Feynman's intuitive idea in resulting functional. In addition, our definition, and our solution, do not place any restrictions on any of the parameters in the Schrödinger equation, and have a potential to offer useful computational experiments, and other theoretical insights.

•Jul 1, 2014•5 citations
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.
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.
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.
Background: Path integrals are a powerful tool for solving problems in quantum theory that are not amenable to a treatment by perturbation theory. Most path integral computations require an analytic continuation to imaginary time. While imaginary time treatments of scattering are possible, imaginary time is not a natural framework for treating scattering problems. More importantly, quantum algorithms for calculating path integrals require real-time evolution.Purpose: We test a recently introduced method for performing direct calculations of scattering observables using real-time path integrals in order to understand the challenges facing real-time path integral calculations of scattering observables.Method: The computations are based on a new interpretation of the path integral as the expectation value of a potential functional on cylinder sets of continuous paths with respect to a complex probability distribution. The method can in principle be applied to arbitrary short-range potentials.Results: The method is applied to compute matrix elements of M\o{}ller wave operators applied to narrow wave packets. These are used to calculate half-shell sharp-momentum transition matrix elements for one-dimensional potential scattering. The calculations for half-shell transition operator matrix elements converge to the numerical solution of the Lippmann-Schwinger equation.Conclusions: This work presents a proof in principle that scattering observables can be computed using real-time Feynman path integrals. While the computational method is not efficient, it can be improved. It provides a laboratory for studying quantum computational algorithms that are applicable to scattering problems.
Abstract The focus of Chapter 12 is the Feynman path integral formulations of quantum mechanics and quantum statistical mechanics. The connection between real-time quantum propagator and the density matrix is established via the Wick rotation, and a Trotter expansion is used to derive discrete path integral expressions for both. The continuous limit of the discrete path integral is derived to yield the functional integral formulation of the path integral. The construction of estimators and their averages for computing quantum observables is discussed. Next, the path integral is extended to incorporate Boltzmann, Fermi-Dirac, and Bose-Einstein statistics, and a useful recursive formula for the partition function is derived that accelerates convergence by building up N-body partition functions from those of smaller systems. Quantum free energies are discussed, and algorithms for the evaluation of path integrals using molecular dynamics or Monte Carlo are presented. The chapter concludes with a discussion of higher-order factorizations and ring-polymer contraction methods to reduce the computational overhead of path-integral calculations
Abstract In this work, the perturbative aspects of quantum mechanics (QM) and quantum field theory (QFT), to a large extent, are studied with functional (path or field) integrals and functional techniques. This physics textbook thus begins with a discussion of algebraic properties of Gaussian measures, and Gaussian expectation values for a finite number of variables. The important role of Gaussian measures is not unrelated to the central limit theorem of probabilities, although the interesting physics is generally hidden in essential deviations from Gaussian distributions. A few algebraic identities about Gaussian expectation values, in particular Wick's theorem are recalled. Integrals over some type of formally complex conjugate variables, directly relevant for boson systems are defined. Fermion systems require the introduction of Grassmann or exterior algebras, and the corresponding generalization of the notions of differentiation and integration. Both for complex and Grassmann integrals, Gaussian integrals, and Gaussian expectation values are calculated, and generalized Wick's theorems proven. The concepts of generating functions and Legendre transformation are recalled.
In 1948, Richard Feynman published the first paper on his new approach to non-relativistic quantum mechanics. Before Feynman's work there were two mathematical formulations of quantum mechanics. Schrödinger's formulation was based on PDE (the Schrödinger equation) and states representation by wave functions, so it was in the framework of analysis and differential equations. The other formulation was Heisenberg's matrix algebra. Initially, they were thought to be competing. The proponents of one claimed that the other was “ wrong. ” Within a couple of years, John von Neumann had proved that they are equivalent.\nAlthough Feynman's theory was not fundamentally new, it nonetheless offered an entirely fresh and different perspective: via a precise formulation of Bohr's correspondence principle, it made quantum mechanics similar to classical mechanics in a precise sense. In addition, Feynman's approach made it possible to explain physical experiments, and, via diagrams, link them directly to computations. What resulted was a very powerful device for computing energies and scattering amplitudes - the famous Feynman's diagrams.\nIn his formulation, Feynman aimed at representing the solution to the non-relativistic Schrödinger equation in the form of an “ average ” over histories or paths of a particle. This solution is commonly known as the Feynman path integral. It plays an important role in the theory but appears as a postulate based on intuition coming from physics rather than a justified mathematical object. This is why Feynman's vision has caught the attention of many mathematicians as well as physicists.\nThe papers of Gelfand, Cameron, and Nelson are among the first, and more substantial, attempts to supply Feynman's theory with a rigorous mathematical foundation. These attempts were followed by many others, but unfortunately all of them were not quite satisfactory. The difficulty comes from a need to define a measure on an infinite-dimensional space of continuous functions that represent all possible paths of a particle. This Feynman's measure has to produce an integral with the properties requested by Feynman. In particular, the expression for the Feynman measure has to involve the non-absolutely integrable Fresnel integrands. The non-absolute integrability of the Fresnel integrands makes the measure fail to be positive and to have the countably additive property. Thus, a well-defined measure in the case of the Feynman path integral does not exist.\nExtensive research has been done on the methods of relating the Feynman path integral to the integral with respect to the Wiener measure. The method of analytic continuation in mass defines the Feynman path integral as a certain limit of the Wiener integrals. Unfortunately, this method can be used as definition for only almost all values of the mass parameter in the Schrödinger equation. For physicists, this is not a satisfactory result and needs to be improved.\nIn this work we examine those questions which originally led to the Feynman path integral. By now we know that Feynman's “ dream ” cannot be realized as a positive and countably additive measure on the path-space. Here, we offer a new way out by modifying Feynman's question, and thereby achieving a solution to the Schrödinger equation via a different kind of averages in the path-space.\nWe give our version of the question that Feynman “ should have asked ” in order to realize the elusive path integral. In our formulation, we get a Feynman path integral as a limit of linear functionals, as opposed to the more familiar inductive limits of positive measures, traditionally used for constructing the Wiener measure, and related Gaussian families. We adapt here an approach pioneered by Patrick Muldowney. In it, Muldowney suggested a Henstock integration technique in order to deal with the non-absolute integrability of the kind of Fresnel integrals which we need in our solution to Feynman's question. By applying Henstock's theory to Fresnel integrals, we construct a complex-valued “ probability distribution functions ” on the path-space. Then we use this “ probability ” distribution function to define the Feynman path integral as an inductive limit. This establishes a mathematically rigorous Feynman limit, and at the same time, preserves Feynman's intuitive idea in resulting functional. In addition, our definition, and our solution, do not place any restrictions on any of the parameters in the Schrödinger equation, and have a potential to offer useful computational experiments, and other theoretical insights.