Matrix-product-state-based band-Lanczos solver for quantum cluster approaches

Sebastian Paeckel,Thomas Köhler,Salvatore R. Manmana,Benjamin Lenz

2023-10-17

Abstract:We present a matrix-product state (MPS) based band-Lanczos method as solver for quantum cluster methods such as the variational cluster approximation (VCA). While a naïve implementation of MPS as cluster solver would barely improve its range of applicability, we show that our approach makes it possible to treat cluster geometries well beyond the reach of exact diagonalization methods. The key modifications we introduce are a continuous energy truncation combined with a convergence criterion that is more robust against approximation errors introduced by the MPS representation and provides a bound to deviations in the resulting Green's function. The potential of the resulting cluster solver is demonstrated by computing the self-energy functional for the single-band Hubbard model at half filling in the strongly correlated regime, on different cluster geometries. Here, we find that only when treating large cluster sizes, observables can be extrapolated to the thermodynamic limit, which we demonstrate at the example of the staggered magnetization. Treating clusters sizes with up to $6\times 6$ sites we obtain excellent agreement with quantum Monte-Carlo results.

Strongly Correlated Electrons,Computational Physics,Quantum Physics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: when dealing with two - dimensional strongly correlated electron systems, traditional methods (such as the exact diagonalization method) are difficult to handle larger - sized quantum clusters due to limitations in computational complexity. Specifically, the article aims to develop a band - like Lanczos algorithm based on matrix product states (MPS) as a solver for quantum cluster methods (such as the variational cluster approximation VCA). By introducing continuous energy truncation and improved convergence criteria, this method can handle cluster geometries far beyond what the exact diagonalization method can achieve, and can calculate the self - energy functional of the single - band Hubbard model in the strongly correlated regime.

### Main problems
1. **Computational complexity**: Traditional methods (such as the exact diagonalization method) can only handle clusters of a limited size due to the exponential growth of the Hilbert space.
2. **Challenges of strongly correlated systems**: The study of quantum states in two - dimensional strongly correlated electron systems (such as high - temperature superconductors) is very difficult and requires more effective numerical methods to handle large - sized clusters.
3. **Limitations of existing methods**: Existing cluster embedding techniques (such as CPT, DMFT, etc.) are effective but have limitations when dealing with two - dimensional systems, especially when considering spectral functions at different wave vectors.

### Solutions
- **MPS - based band - like Lanczos algorithm**: Combine MPS representation and the band - like Lanczos method, and reduce the growth of entanglement entropy through continuous energy truncation, so that larger - sized clusters can be handled.
- **Improved convergence criteria**: Introduce the Hochbruck - Lubich criteria to provide a more robust convergence judgment standard and avoid numerical instability caused by MPS arithmetic and truncation errors.
- **Application examples**: By calculating the self - energy functional of the single - band Hubbard model at half - filling, the effectiveness of this method is demonstrated, especially when dealing with 6x6 - sized clusters, which is highly consistent with the quantum Monte Carlo results.

### Key formulas
1. **Hamiltonian of the Hubbard model**:
\[
\hat{H}=-t \sum_{\langle i,j \rangle, \sigma}\left(\hat{c}_{i,\sigma}^{\dagger} \hat{c}_{j,\sigma}+\text{h.c.}\right)+U \sum_{i} \hat{n}_{i,\uparrow} \hat{n}_{i,\downarrow}
\]
2. **Definition of the Green function**:
\[
G^{(e)}_{\mu\mu'}(\omega)=\langle\psi|\hat{c}_{\mu}(\hat{H}-E_{0}-\omega)^{-1} \hat{c}_{\mu'}^{\dagger}|\psi\rangle
\]
3. **Energy truncation formula**:
\[
|\psi\rangle=\sum_{n}\langle\tilde{E}_{n}|\psi\rangle|\tilde{E}_{n}\rangle\approx\sum_{\epsilon_{n}>\delta} \epsilon_{n}|\tilde{E}_{n}\rangle
\]
4. **Hochbruck - Lubich error bound**:
\[
\epsilon_{D} \leq \frac{40}{\tau} e^{-\tau / 4}\left(e^{\tau / 4D}\right)^{D}
\]

Through these improvements, this method can significantly expand the size of the clusters that can be processed while maintaining computational accuracy, providing a new tool for studying the properties of two - dimensional strongly correlated electron systems.