Quantum Simulation of the Fokker-Planck Equation via Schrödingerization

Author(s)

,
&

Abstract

This paper studies a quantum simulation technique for solving the Fokker-Planck equation. Traditional semi-discretization methods often fail to preserve the underlying Hamiltonian dynamics and may even modify the Hamiltonian structure, particularly when incorporating boundary conditions. We address this challenge by employing the Schrödingerization method – it converts any linear partial and ordinary differential equation with non-Hermitian dynamics into systems of Schrödinger-type equations. It does so via the so-called warped phase transformation that maps the equation into one higher dimension. We explore the application in two distinct forms of the Fokker-Planck equation. For the conservation form, we show that the semi-discretization-based Schrödingerization is preferable, especially when dealing with non-periodic boundary conditions. Additionally, we analyze the Schrödingerization approach for unstable systems that possess positive eigenvalues in the real part of the coefficient matrix or differential operator. Our analysis reveals that the direct use of Schrödingerization has the same effect as a stabilization procedure. For the heat equation form, we propose a quantum simulation procedure based on the time-splitting technique, and give explicitly its corresponding quantum circuit. We discuss the relationship between operator splitting in the Schrödingerization method and its application directly to the original problem, illustrating how the Schrödingerization method accurately reproduces the time-splitting solutions at each step. Furthermore, we explore finite difference discretizations of the heat equation form using shift operators. Utilizing Fourier bases, we diagonalize the shift operators, enabling efficient simulation in the frequency space. Providing additional guidance on implementing the diagonal unitary operators, we conduct a comparative analysis between diagonalizations in the Bell and the Fourier bases, and show that the former generally exhibits greater efficiency than the latter.

Author Biographies

  • Shi Jin

    School of Mathematical Sciences, Institute of Natural Sciences, MOE-LSC, Shanghai Jiao Tong University, Shanghai, 200240, China

  • Nana Liu

    School of Mathematical Sciences, Institute of Natural Sciences, MOE-LSC, Shanghai Jiao Tong University, Shanghai, 200240, China

    University of Michigan-Shanghai Jiao Tong University Joint Institute, Shanghai 200240, China

  • Yue Yu

    Hunan Key Laboratory for Computation and Simulation in Science and Engineering, Key Laboratory of Intelligent Computing and Information Processing of Ministry of Education, National Center for Applied Mathematics in Hunan, School of Mathematics and Computational Science, Xiangtan University, Xiangtan, Hunan 411105, China

About this article

Abstract View

  • 372

Pdf View

  • 117

DOI

10.4208/cicp.OA-2024-0097

How to Cite

Quantum Simulation of the Fokker-Planck Equation via Schrödingerization. (2026). Communications in Computational Physics, 39(4), 969-1001. https://doi.org/10.4208/cicp.OA-2024-0097