We study the connection between the exponent of the order parameter of the Mott insulator-tosuperfluid transition occurring in the two-dimensional Bose-Hubbard model, and the divergence exponents of its one-and two-particle correlation functions. We find that at the multicritical points all divergence exponents are related to each other, allowing us to express the critical exponent in terms of one single divergence exponent. This approach correctly reproduces the critical exponent of the threedimensional XY universality class. Because divergence exponents can be computed in an efficient manner by hypergeometric analytic continuation, our strategy is applicable to a wide class of systems.

We explore in detail how analytic continuation of divergent perturbation series by generalized hypergeometric functions is achieved in practice. Using the example of strong-coupling perturbation series provided by the two-dimensional Bose-Hubbard model, we compare hypergeometric continuation to Shanks and Padé techniques, and demonstrate that the former yields a powerful, efficient and reliable alternative for computing the phase diagram of the Mott insulator-to-superfluid transition. In contrast to Shanks transformations and Padé approximations, hypergeometric continuation also allows us to determine the exponents which characterize the divergence of correlation functions at the transition points. Therefore, hypergeometric continuation constitutes a promising tool for the study of quantum phase transitions.

We develop a scheme for analytic continuation of the strong-coupling perturbation series of the pure Bose-Hubbard model beyond the Mott insulator-to-superfluid transition at zero temperature, based on hypergeometric functions and their generalizations. We then apply this scheme for computing the critical exponent of the order parameter of this quantum phase transition for the twodimensional case, which falls into the universality class of the three-dimensional XY model. This leads to a nontrivial test of the universality hypothesis.

The Mott insulator-to-superfluid transition exhibited by the Bose-Hubbard model on a two-dimensional square lattice occurs for any value of the chemical potential, but becomes critical at the tips of the so-called Mott lobes only. Employing a numerical approach based on a combination of high-order perturbation theory and hypergeometric analytic continuation we investigate how quantum critical properties manifest themselves in computational practice. We consider two-dimensional triangular lattices and three-dimensional cubic lattices for comparison, providing accurate parametrizations of the phase boundaries at the tips of the respective first lobes. In particular, we lend strong support to a recently suggested inequality which bounds the divergence exponent of the one-particle correlation function in terms of that of the two-particle correlation function, and which sharpens to an equality if and only if a system becomes critical. arXiv:1912.00679v1 [cond-mat.stat-mech]