DMRG Formulation of the Orbital Hatsugai-Kohmoto model in 2D

The orbital extension of the Hatsugai-Kohmoto model (OHK) has recently been used to simulate Hubbard physics on an exactly solvable system. OHK cures many of the pathologies of the "band" HK model while providing a systematic scheme to connect it to the Hubbard model through the number of orbitals. We demonstrate how the OHK scheme can be pushed further by implementing with DMRG, utilizing the connection between OHK and Hubbard with twisted boundary conditions. Further, we make a powerful scaling argument to show that OHK with fixed orbital number converges to Hubbard more rapidly when increasing the spatial dimension. Thus, the OHK+DMRG scheme can be used as a powerful tool for simulating 2D and higher strongly correlated physics with controllable numerics.