Mott-insulator-superfluid-liquid transition in a one-dimensional bosonic Hubbard model: Quantum Monte Carlo method

The values of the insulator gap Δ in one-dimensional systems of interacting bosons described by the Hubbard Hamiltonian are calculated at low temperatures by the quantum world-line Monte Carlo algorithm. The dependence of Δ on the size of the system, the temperature, and the parameters of the model is investigated. It is shown that a chain with N _a=50 sites is already sufficient to estimate the thermodynamic value of the critical quantity (t/U)_c for which a transition from the insulator into the superfluid state occurs in a commensurate system. To within the computational error, this value, (t/U)_c=0.300±0.005, agrees with the value (t/U)_c=0.304±0.002 obtained previously by the combined “exact diagonalization + renormalization-group analysis” method. The characteristic Kosterlitz-Thouless behavior of the insulator gap is demonstrated near the critical region: Δ∼exp[−b(1−t/t _c)^−1/2]. Pis’ma Zh. éksp. Teor. Fiz. 64, No. 2, 92–96 (25 July 1996)