“…For fermions in one spatial dimension, the integrability of the Fermi-Hubbard model facilitates the derivation of the exact evolution after a quench including effects such as "pre-thermalization" [13,50,51]. Again, in higher dimensions, appropriate approximations are necessary, such as a time-dependent Monte-Carlo method [52], time-dependent dynamical mean field theory [53][54][55][56][57], the Gutzwiller ansatz for fermions [58][59][60], the flow equation method [9,10,61,62], or effective quasi-particle methods [63].…”