  The real time evolution of two pieces of quantum insulators, initially at
different temperatures, is studied when they are glued together. Specifically,
each subsystem is taken as a Bose-Hubbard model in a Mott insulator state. The
process of temperature equilibration via heat transfer is simulated in real
time using the Minimally Entangled Typical Thermal States algorithm. The
analytic theory based on quasiparticles transport is also given.
