Mott physics in the half-filled Hubbard model on a family of vortex-full square lattices

We study the half-filled Hubbard model on a one-parameter family of vortex-full square lattices ranging from the isotropic case to weakly coupled Hubbard dimers. The ground-state phase diagram consists of four phases: A semimetal and a band insulator which are connected to the weak-coupling limit, and a magnetically ordered Neel phase and a valence bond solid (VBS) which are linked to the strong-coupling Mott limit. The phase diagram is obtained by quantum Monte Carlo (QMC) and continuous unitary transformations (CUTs). The CUT is performed in a two-step process: Nonperturbative graph-based CUTs are used in the Mott insulating phase to integrate out charge fluctuations. The resulting effective spin model is tackled by perturbative CUTs about the isolated dimer limit yielding the breakdown of the VBS by triplon condensation. We find three scenarios when varying the interaction for a fixed anisotropy of hopping amplitudes: (i) one direct phase transition from Neel to semimetal, (ii) two phase transitions VBS to Neel and Neel to semimetal, or (iii) a smooth crossover from VBS to the band insulator. Our results are consistent with the absence of spin-liquid phases in the whole phase diagram.