Abstract:We propose quantum algorithms for projective ground-state preparation and calculations of the many-body Green's functions directly in frequency domain. The algorithms are based on the linear combination of unitary (LCU) operations and essentially only use quantum resources. To prepare the ground state, we construct the operator ${\exp}(-\tau \hat{H}^2)$ using Hubbard-Stratonovich transformation by LCU and apply it on an easy-to-prepare initial state. Our projective state preparation procedure saturates the near-optimal scaling, $\mathcal{O}(\frac{1}{\gamma\Delta} \log \frac{1}{\gamma\eta})$, of other currently known algorithms, in terms of the spectral-gap lower bound $\Delta$, the additive error $\eta$ in the state vector, and the overlap lower bound $\gamma$ between the initial state and the exact ground state. It is straightforward to combine our algorithm with spectral-gap amplification technique to achieve quadratically improved scaling $\mathcal{O}(1/\sqrt{\Delta})$ for ground-state preparation of frustration-free Hamiltonians, which we demonstrate with numerical results of the $q$-deformed XXZ chain. To compute the Green's functions, including the single-particle and other response functions, we act on the prepared ground state with the retarded resolvent operator $R(\omega + i\Gamma; \hat{H})$ in the LCU form derived from the Fourier-Laplace integral transform (FIT). Our resolvent algorithm has $\mathcal{O}(\frac{1}{\Gamma^2} \log\frac{1}{\Gamma\epsilon})$ complexity scaling for the frequency resolution $\Gamma$ of the response functions and the targeted error $\epsilon$, while classical algorithms for FIT usually have polynomial scaling over the error $\epsilon$. To illustrate the complexity scaling of our algorithms, we provide numerical results for their application to the paradigmatic Fermi-Hubbard model on a one-dimensional lattice with different numbers of sites.

Global fermionic mode optimization via swap gates

Two dimensional quantum lattice models via mode optimized hybrid CPU-GPU density matrix renormalization group method

Optimizing Strongly Interacting Fermionic Hamiltonians

Efficient Representation of Gaussian Fermionic Pure States in Non-Computational Bases

Efficient conversion from fermionic Gaussian states to matrix product states

Optimizing sparse fermionic Hamiltonians

Expanding the reach of quantum optimization with fermionic embeddings

Quantum Algorithms for Ground-State Preparation and Green's Function Calculation

A recipe for local simulation of strongly-correlated fermionic matter on quantum computers: the 2D Fermi-Hubbard model

Ever more optimized simulations of fermionic systems on a quantum computer

Quantum Computed Green's Functions using a Cumulant Expansion of the Lanczos Method

Optimal Particle-Conserved Linear Encoding for Practical Fermionic Simulation

Discovering optimal fermion-qubit mappings through algorithmic enumeration

Low-depth simulations of fermionic systems on square-grid quantum hardware

High-efficiency quantum Monte Carlo algorithm for extracting entanglement entropy in interacting fermion systems

From bosons and fermions to spins: A multi-mode extension of the Jordan-Schwinger map

Low-depth circuit ansatz for preparing correlated fermionic states on a quantum computer

Many-Body Quantum States with Exact Conservation of Non-Abelian and Lattice Symmetries through Variational Monte Carlo

Gate Efficient Composition of Hamiltonian Simulation and Block-Encoding with its Application on HUBO, Fermion Second-Quantization Operators and Finite Difference Method

Quantum Computation and Simulation using Fermion-Pair Registers