“…This method works very well for J/U below the critical point (J/U ) c , which is less than one in all dimensions, where all physical quantities are analytical functions. It was used in the studies of spinless bosons with local [51][52][53][54][250][251][252][253][254][255][256][257][258][259] and nearest-neighbor interactions [260,261], two-species bosons [262], and spin-1 bosons [263,264] in different types of regular lattices like hypercubic isotropic [51,52,54,255,[260][261][262][263]265] and anisotropic [253], superlattices [254,255], two-dimensional triangular and kagome lattices [256][257][258], and in the presence of artificial gauge fields [259,265], as well as in disordered lattices [52,67,[266][267][268]. The strong-coupling expansion allows to avoid finite-size effects and shows excellent agreement with exact numerical data.…”