Luttinger liquid tensor network: sine versus tangent dispersion of massless Dirac fermions

V. A. Zakharov,S. Polla,A. Donís Vela,P. Emonts,M. J. Pacholski,J. Tworzydło,C. W. J. Beenakker

2024-09-19

Abstract:To apply the powerful many-body techniques of tensor networks to massless Dirac fermions one wants to discretize the $p\cdot\sigma$ Hamiltonian and construct a matrix-product-operator (MPO) representation. We compare two alternative discretization schemes, one with a sine dispersion, the other with a tangent dispersion, applied to a one-dimensional Luttinger liquid with Hubbard interaction. Both types of lattice fermions allow for an exact MPO representation of low bond dimension, so they are efficiently computable, but only the tangent dispersion gives a power law decay of the propagator in agreement with the continuum limit: The sine dispersion is gapped by the interactions, evidenced by an exponentially decaying propagator. Our construction of a tensor network with an unpaired Dirac cone works around the fermion-doubling obstruction by exploiting the fact that the \textit{nonlocal} Hamiltonian of tangent fermions permits a \textit{local} generalized eigenproblem.

Strongly Correlated Electrons,Mesoscale and Nanoscale Physics,Quantum Physics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: how to effectively discretize the Luttinger liquid model of a continuous system onto a one - dimensional lattice while maintaining the topologically protected properties of massless Dirac fermions, and construct a low - rank matrix product operator (MPO) representation for efficient numerical calculations using tensor network methods such as DMRG.

Specifically, the author compares two different discretization schemes: one is sine dispersion, and the other is tangent dispersion. Both of these schemes allow the construction of low - rank MPO representations, but only the tangent dispersion can maintain a power - law decaying propagator consistent with the continuous limit, while the sine dispersion will introduce an energy gap due to the presence of interactions, resulting in an exponentially decaying propagator.

### Main problem summary:
1. **Maintaining topological protection**: Avoid the fermion doubling problem during the discretization process to maintain the topologically protected properties of massless Dirac fermions.
2. **Efficient tensor network representation**: Construct a low - rank matrix product operator (MPO) representation so that efficient numerical calculations can be carried out using tensor network methods such as the density matrix renormalization group (DMRG).
3. **Verifying the effectiveness of discretization schemes**: Verify the performance of different discretization schemes in free fermion and systems containing Hubbard interactions by comparing with the analytical results of the continuous limit.

### Formula presentation:
- **Sine dispersion relation**:
\[
E(k)=\frac{\hbar v}{a}\sin(ak)
\]
- **Tangent dispersion relation**:
\[
E(k)=\frac{2\hbar v}{a}\tan\left(\frac{ak}{2}\right)
\]
- **Continuous limit of the propagator at zero temperature**:
\[
C_{\sigma}(x,x')=\frac{\sigma\hbar v}{2i\beta\sinh\left[\pi\left(\frac{\hbar v}{\beta}\right)(x - x')\right]}
\]
When \(\beta\rightarrow\infty\):
\[
C_{\sigma}(x,x')=\frac{\sigma}{2\pi i(x - x')}
\]

Through these studies, the author hopes to find a discretization method that can achieve efficient numerical simulation while maintaining physical properties.