“…The Bose-Hubbard model (BHM), known as one of the simplest models that captures the essence of the superfluid-Mott insulator transition [1], has given rise to a plethora of studies with cold bosonic atoms in optical lattices [2][3][4][5] in which the ratio between hopping and onsite interaction is widely tunable. Extensions to multicomponent BHMs enable studies of polaron physics [6,7] and quantum magnetism [8][9][10][11], and by adding a dynamical coupling between the components it furthermore is possible to implement BHMs that mimic radiative effects [12][13][14][15] in waveguide QED, [14,[16][17][18][19][20][21] featuring fractional decay, bound states [22,23], and polaritons [15,24,25]. This also provides a connection to photon-based many-body physics [24,26,27] in the microwave domain [28][29][30][31][32].…”