High-precision Monte Carlo study of several models in the three-dimensional U(1) universality class

Abstract

We present a worm-type Monte Carlo study of several typical models in the three-dimensional (3D) U(1) universality class, which include the classical 3D XY model in the directed flow representation and its Villain version, as well as the 2D quantum Bose-Hubbard (BH) model with unitary filling in the imaginary-time world-line representation. From the topology of the configurations on a torus, we sample the superfluid stiffness ${\rho }_s$ and the dimensionless wrapping probability $R$. From the finite-size scaling analyses of ${\rho }_s$ and of $R$, we determine the critical points as $T_c(\mathrm{XY})=2.2018441(5)$ and $T_c(\mathrm{Villain})=0.33306704(7)$ and $(t/U{)}_c(\mathrm{BH})=0.0597291(8)$, where $T$ is the temperature for the classical models, and $t$ and $U$ are respectively the hopping and on-site interaction strength for the BH model. The precision of our estimates improves significantly over that of the existing results. Moreover, it is observed that at criticality, the derivative of a wrapping probability with respect to $T$ suffers from negligible leading corrections and enables a precise determination of the correlation length critical exponent as $\nu =0.67183(18)$. In addition, the critical exponent $\eta$ is estimated as $\eta =0.03853(48)$ by analyzing a susceptibility-like quantity. We believe that these numerical results would provide a solid reference in the study of classical and quantum phase transitions in the 3D U(1) universality, including the recent development of the conformal bootstrap method.

Comments: 16 pages, 8 figures, 12 tables