Abstract:Quantum phase estimation is one of the most powerful quantum primitives. This work proposes a new approach for the problem of multiple eigenvalue estimation: Quantum Multiple Eigenvalue Gaussian filtered Search (QMEGS). QMEGS leverages the Hadamard test circuit structure and only requires simple classical postprocessing. QMEGS is the first algorithm to simultaneously satisfy the following two properties: (1) It can achieve the Heisenberg-limited scaling without relying on any spectral gap assumption. (2) With a positive energy gap and additional assumptions on the initial state, QMEGS can estimate all dominant eigenvalues to $\epsilon$ accuracy utilizing a significantly reduced circuit depth compared to the standard quantum phase estimation algorithm. In the most favorable scenario, the maximal runtime can be reduced to as low as $\log(1/\epsilon)$. This implies that QMEGS serves as an efficient and versatile approach, achieving the best-known results for both gapped and gapless systems. Numerical results validate the efficiency of our proposed algorithm in various regimes.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to estimate multiple eigenvalues efficiently and accurately in quantum computing. Specifically, the paper proposes a new method - Quantum Multiple Eigenvalue Gaussian filtered Search (QMEGS) to solve the multi - eigenvalue estimation problem. This method aims to overcome the limitations of existing algorithms when dealing with gapless systems, and at the same time, it can also significantly reduce the circuit depth and improve the algorithm efficiency in the presence of an energy gap.
### Main contributions
1. **Heisenberg - limit scaling under the gapless assumption**: QMEGS is the first algorithm that can achieve Heisenberg - limit scaling without relying on any spectral gap assumptions.
2. **Efficient estimation when there is a positive energy gap**: When there is a positive energy gap and additional assumptions on the initial state, QMEGS can use a significantly reduced circuit depth to estimate all the main eigenvalues with an accuracy of \(\epsilon\).
3. **Short circuit depth**: In the most favorable case, the maximum running time can be reduced to as low as \(\log(1/\epsilon)\), which means that QMEGS can achieve the best - known results in both gapped and gapless systems.
4. **Robustness**: QMEGS can handle imperfect initial states, that is, the case where \(p_{\text{tail}}>0\).
### Method overview
The main steps of the QMEGS method include:
1. **Data generation**: Use the Hadamard test circuit to generate a data set \((t_n, Z_n)\), where \(Z_n\) is an unbiased estimate of \(\langle\psi|e^{-it_nH}|\psi\rangle\).
2. **Filtering and searching**: Define a Gaussian filtering function \(G(\theta)\), and find approximate values of the main eigenvalues through an iterative search process.
### Mathematical representation
- **Eigenvalues and eigenvectors**: Let \((\lambda_m,|\psi_m\rangle)\) represent the eigenvalue and eigenvector pairs of the Hamiltonian \(H\), and \(p_m = |\langle\psi_m|\psi\rangle|^2\) represents the overlap between the initial state and the \(m\)-th eigenvector.
- **Main eigenvalue condition**: Assume that there is a set of indices \(D\subset\{1, 2,\ldots,M\}\) such that \(p_{\text{min}}=\min_{i\in D}p_i > p_{\text{tail}}=\sum_{i\in D^c}p_i\), where \(D^c = \{1, 2,\ldots,M\}\setminus D\).
- **Gaussian filtering function**: Define the filtering function \(G(\theta)=\left|\frac{1}{N}\sum_{n = 1}^N Z_n\exp(i\theta t_n)\right|\).
- **Fourier transform**: Define the Fourier transform \(F(x)=\int_{-\infty}^{\infty}a(t)\exp(ixt)dt\) of the probability density function \(a(t)\).
### Algorithm complexity
- **Maximum running time**: \(T_{\text{max}}=\max_{1\leq n\leq N}|t_n|\).
- **Total running time**: \(T_{\text{total}}=\sum_{n = 1}^N|t_n|\).
### Experimental verification
The numerical results verify the efficiency of the QMEGS algorithm in various situations (gapped and gapless systems).
In conclusion, this paper solves the key problems in multi - eigenvalue estimation by proposing the QMEGS method, especially the efficient estimation in gapless systems and the short - circuit - depth optimization when there is an energy gap.