Improved quantum algorithm for calculating eigenvalues of differential operators and its application to estimating the decay rate of the perturbation distribution tail in stochastic inflation

Quantum algorithms for scientific computing and their applications have been actively studied. In this paper, we propose a quantum algorithm for estimating the first eigenvalue of a differential operator L on R^{d} and its application to cosmic inflation theory. A common approach for this eigenvalue...

Full description

Saved in:
Bibliographic Details
Main Authors: Koichi Miyamoto, Yuichiro Tada
Format: Article
Language:English
Published: American Physical Society 2025-06-01
Series:Physical Review Research
Online Access:http://doi.org/10.1103/wpnm-rlrl
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:Quantum algorithms for scientific computing and their applications have been actively studied. In this paper, we propose a quantum algorithm for estimating the first eigenvalue of a differential operator L on R^{d} and its application to cosmic inflation theory. A common approach for this eigenvalue problem involves applying the finite-difference discretization to L and computing the eigenvalues of the resulting matrix, but this method suffers from the curse of dimensionality, namely, the exponential complexity with respect to d. Our first contribution is the development of a quantum algorithm for this task, leveraging recent quantum singular value transformation-based methods. Given a trial function that overlaps well with the eigenfunction, our method runs with query complexity scaling as O[over ̃](d^{3}/ε^{2}) with d and estimation accuracy ε, which is polynomial in d and shows improvement over existing quantum algorithms. Then, we consider the application of our method to a problem in a theoretical framework for cosmic inflation known as stochastic inflation, specifically calculating the eigenvalue of the adjoint Fokker–Planck operator, which is related to the decay rate of the tail of the probability distribution for the primordial density perturbation. We numerically see that, in some cases, simple trial functions overlap well with the first eigenfunction, indicating our method is promising for this problem.
ISSN:2643-1564