Curator's Take
This article presents a clever quantum algorithm for solving the Fokker-Planck equation, which governs how probability distributions evolve over time in systems with both drift and random fluctuations - a fundamental problem in fields from finance to robotics. The researchers' key insight is encoding probability densities directly into quantum amplitudes, allowing them to represent exponentially large state spaces with a linear number of qubits, potentially offering dramatic speedups for tracking complex, high-dimensional probability distributions. While they had to introduce a creative workaround for the diffusion term using a "Wick rotation" technique, the approach opens new possibilities for quantum-enhanced Bayesian filtering and state estimation in real-world applications. This work bridges quantum computing with classical control theory and could eventually enable quantum sensors and autonomous systems to process uncertainty more efficiently than classical methods allow.
— Mark Eatherly
Summary
We propose a gate-based quantum algorithm for the prediction step of Bayesian state estimation based on the Fokker-Planck equation on a discretized position-velocity state space. The probability density is encoded in the amplitudes of a quantum state, enabling a compact representation of high-dimensional distributions. Exploiting the circulant structure of finite-difference operators, the evolution is realized in the spectral domain using quantum Fourier transforms and phase rotations. A key result is that the drift component can be implemented exactly in amplitude space, leading to an accurate reproduction of the classical transport dynamics. In contrast, the diffusion term does not admit a linear representation in amplitude space due to the nonlinear relation between probability density and wave function. To enable a quantum implementation, we introduce a unitary surrogate based on a Wick rotation, transforming diffusion into a dispersive phase evolution. This yields a fully unitary propagation that can be implemented efficiently on a gate-based quantum computer. The proposed method is evaluated numerically for different scenarios and shows strong agreement with the exact solution of the Fokker-Planck equation. The approach demonstrates the potential of quantum computing for Bayesian state estimation, as the representable state space grows exponentially with the number of qubits. This allows the efficient representation and propagation of probability densities that would otherwise require complex tensor decompositions on classical hardware, making the method a promising candidate for high-dimensional filtering problems.