Quantum Quench and Charge Oscillations in the SU(3) Hubbard Model: a Test of Time Evolving Block Decimation with general non-Abelian Symmetries

TL;DR

This paper develops a non-Abelian tensor network method to efficiently simulate quantum dynamics in the SU(3) Hubbard model, revealing how interactions affect charge relaxation and quantum oscillations.

Contribution

It introduces a general non-Abelian TEBD scheme leveraging symmetries to enhance computational efficiency and applies it to study post-quench dynamics in the SU(3) Hubbard model.

Findings

Interactions convert algebraic charge relaxation into exponential decay

Quantum oscillations are rapidly suppressed by interactions

The method enables large bond dimension simulations on standard hardware

Abstract

We introduce the notion of non-Abelian tensors, and use them to construct a general non-Abelian time evolving block decimation (NA-TEBD) scheme that uses an arbitrary number of Abelian and non-Abelian symmetries. Our approach increases the speed and memory storage efficiency of matrix product state based computations by several orders of magnitudes, and makes large bond dimensions accessible even on simple desktop architectures. We use it to study post-quench dynamics in the repulsive SU(3) Hubbard model, and to determine the time evolution of various local operators and correlation functions efficiently. Interactions turn algebraic charge relaxation into exponential, and suppress coherent quantum oscillations rapidly.