Anderson impurity solver integrating tensor network methods with quantum computing

Francois Jamet,Lachlan P. Lindoy,Yannic Rath,Connor Lenihan,Abhishek Agarwal,Enrico Fontana,Fedor Simkovic IV,Baptiste Anselme Martin,Ivan Rungger

2024-10-21

Abstract:Solving the Anderson impurity model typically involves a two-step process, where one first calculates the ground state of the Hamiltonian, and then computes its dynamical properties to obtain the Green's function. Here we propose a hybrid classical/quantum algorithm where the first step is performed using a classical computer to obtain the tensor network ground state as well as its quantum circuit representation, and the second step is executed on the quantum computer to obtain the Green's function. Our algorithm exploits the efficiency of tensor networks for preparing ground states on classical computers, and takes advantage of quantum processors for the evaluation of the time evolution, which can become intractable on classical computers. We demonstrate the algorithm using 24 qubits on a quantum computing emulator for SrVO$_3$ with a multi-orbital Anderson impurity model within the dynamical mean field theory. The tensor network based ground state quantum circuit preparation algorithm can also be performed for up to 60 qubits with our available computing resources, while the state vector emulation of the quantum algorithm for time evolution is beyond what is accessible with such resources. We show that, provided the tensor network calculation is able to accurately obtain the ground state energy, this scheme does not require a perfect reproduction of the ground state wave function on the quantum circuit to give an accurate Green's function. This hybrid approach may lead to quantum advantage in materials simulations where the ground state can be computed classically, but where the dynamical properties cannot.

Quantum Physics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: in the multi - orbital Anderson impurity model (AIM), calculate the Green's function to obtain dynamic properties. Specifically, this problem involves two main steps:

1. **Calculate the ground state of the Hamiltonian**: This step requires obtaining the ground - state wave function and energy of the system.
2. **Calculate dynamic properties**: Based on the ground - state wave function, further calculate the dynamic properties of the system, such as the Green's function.

Traditional methods usually use the classical tensor network (TN) method to handle the first step, but for the time - evolution calculation in the second step, it may become intractable on a classical computer because as the system complexity increases, the computational resource requirements grow exponentially.

To solve this problem, this paper proposes a hybrid classical/quantum algorithm:

- The **classical computer** is used to obtain the ground - state wave function and its corresponding quantum - circuit representation through the tensor - network method.
- The **quantum computer** is used to perform time - evolution simulations, thereby calculating the Green's function.

This method combines the efficiency of tensor networks in preparing the ground state on a classical computer and the advantages of quantum processors in evaluating time - evolution, enabling more effective calculation of dynamic properties for complex systems that are difficult to handle on a classical computer.

### Formula Summary

#### Hamiltonian of the Anderson Impurity Model
The Hamiltonian \(H\) of the Anderson impurity model is expressed as follows:
\[
H=\sum_{ij}\sum_{\sigma}\epsilon_{ij}\hat{c}_{i\sigma}^{\dagger}\hat{c}_{j\sigma}+\sum_{i}U n_{i\uparrow}n_{i\downarrow}+\sum_{i > j}\sum_{\sigma}(U - 2J)n_{i\sigma}n_{j\bar{\sigma}}+(U - 3J)n_{i\sigma}n_{j\sigma}
\]
\[
-\sum_{i > j}\sum_{\sigma}J(\hat{c}_{i\sigma}^{\dagger}\hat{c}_{i\bar{\sigma}}\hat{c}_{j\bar{\sigma}}^{\dagger}\hat{c}_{j\sigma}-\hat{c}_{i\sigma}^{\dagger}\hat{c}_{i\bar{\sigma}}^{\dagger}\hat{c}_{j\sigma}\hat{c}_{j\bar{\sigma}})+\sum_{i}\sum_{j}\sum_{\sigma}V_{ij}(\hat{c}_{i\sigma}^{\dagger}\hat{d}_{j\sigma}+\hat{d}_{j\sigma}^{\dagger}\hat{c}_{i\sigma})+\sum_{ij}\sum_{\sigma}\epsilon_d^{ij}\hat{d}_{i\sigma}^{\dagger}\hat{d}_{j\sigma}
\]

#### Definition of the Green's Function
The zero - temperature large Green's function \(G^>_{\alpha\beta}(z)\) and the small Green's function \(G^<_{\alpha\beta}(z)\) are defined as follows respectively:
\[
G^>_{\alpha\beta}(z)=\langle GS|\hat{c}_\alpha(z - (\hat{H}-E_{GS}))^{-1}\hat{c}_\beta^{\dagger}|GS\rangle
\]
\[
G^<_{\alpha\beta}(z)=\langle GS|\hat{c}_\beta^{\dagger}(z + (