The two-dimensional ($2$D) spin-$S=1$ $XY$ model was investigated numerically as a realization of the $(2+1)$D superfluid-Mott-insulator (SF-MI) transition. The interaction parameters are extended so as to suppress corrections to finite-size scaling. Thereby, the external field of a unit flux quantum ($\Phi=2\pi$) is applied to the 2D cluster by incorporating the phase factor $e^{ i\phi_{ ij}}$ ($\phi_{ ij}$: gauge angle between the $i$ and $j$ sites) into the hopping amplitudes. Taking the advantage in that the exact-diagonalization method allows us to treat such a complex-valued matrix element, we evaluated the excess energy cost $\Delta E(2\pi)$ due to the magnetic flux $\Phi=2\pi$ explicitly in the SF ($XY$) phase. As a result, we found that the amplitude ratio $\rho_s / \Delta E(2\pi)$ ($\rho_s$: spin stiffness) makes sense in proximity to the critical point, exhibiting a notable plateau in the SF-phase side. The plateau height is estimated, and compared to the related studies.
Categories: cond-mat.stat-mech