Comparative study of the density matrix embedding theory for the Hubbard models

TL;DR

This paper evaluates the density matrix embedding theory (DMET) for Hubbard models, emphasizing the entanglement spectrum as a key criterion for assessing the quality of local approximations and phase transitions.

Contribution

It introduces the use of the entanglement spectrum as a sensitive measure to evaluate DMET's effectiveness in capturing local physics and phase transitions in Hubbard models.

Findings

Entanglement spectrum is more sensitive than energy and double occupancy.

Symmetric potentials reproduce the ES of phases connected to noninteracting limits.

ES reveals Mott and symmetry-breaking transitions through singularities.

Abstract

We examine the performance of the density matrix embedding theory (DMET) recently proposed in [G. Knizia and G. K.-L. Chan, Phys. Rev. Lett. 109, 186404 (2012)]. The core of this method is to find a proper one-body potential that generates a good trial wave function for projecting a large scale original Hamiltonian to a local subsystem with a small number of basis. The resultant ground state of the projected Hamiltonian can locally approximate the true ground state. However, the lack of the variational principle makes it difficult to judge the quality of the choice of the potential. Here we focus on the entanglement spectrum (ES) as a judging criterion; accurate evaluation of the ES guarantees that the corresponding reduced density matrix well reproduces all physical quantities on the local subsystem. We apply the DMET to the Hubbard model on the one-dimensional chain, zigzag chain, and triangular lattice and test several variants of potentials and cost functions. It turns out that ES serves as a more sensitive quantity than the energy and double occupancy to probe the quality of the DMET outcomes. A symmetric potential reproduces the ES of the phase that continues from a noninteracting limit. The Mott transition as well as symmetry-breaking transitions can be detected by the singularities in the ES. However, the details of the ES in the strongly interacting parameter region depends much on these variants, meaning that the present DMET algorithm allowing for numerous variant is insufficient to fully characterize the particular phases that require characterization by the ES.