We consider the Bose-Hubbard model in two and three spatial dimensions numerically compute quantum circuit complexity of ground state Mott insulator superfluid phases using a mean field approximation with additional quadratic fluctuations. After mapping to qubit system, result is given by associated Bogoliubov transformation applied reference taken be state. In particular, has peaks at $O(2)$ c...