“…The phase diagram of the disordered Bose-Hubbard model has been studied by a number of methods including the quantum Monte Carlo [59][60][61][62][63], renormalization group [64,65], density-matrix renormalization group techniques [66][67][68], tensor networks-based algorithms, or various mean-field approaches [6,52,53,69]. In this work we propose an extension of the local mean-field method, thus let us first briefly review the mean-field approaches used earlier.…”