  We consider the two-dimensional classical XY model on a square lattice in the
thermodynamic limit using tensor renormalization group and precisely determine
the critical temperature corresponding to the Berezinskii-Kosterlitz-Thouless
(BKT) phase transition to be 0.89290(5) which is an improvement compared to
earlier studies using tensor network methods.
