Tensor network based finite-size scaling for two-dimensional Ising model

We propose a scheme to perform tensor network based finite-size scaling analysis for two-dimensional classical models. In the tensor network representation of the partition function, we use higher-order tensor renormalization group (HOTRG) method to coarse grain the weight tensor. The renormalized tensor is then used to construct the approximated transfer matrix of an infinite strip of finite width. By diagonalizing the transfer matrix we obtain the correlation length, the magnetization, and the energy density, which are used in finite-size scaling analysis to determine the critical temperature and the critical exponents. As a benchmark we study the two-dimensional classical Ising model. We show that the critical temperature and the critical exponents can be accurately determined. With HOTRG bond dimension $D=70$, the absolute errors of the critical temperature $T_c$ and the critical exponent $\nu$, $\beta$ are at the order of ${10}^{-7},{10}^{-5}$, ${10}^{-4}$, respectively. Furthermore, the results can be systematically improved by increasing the bond dimension of the HOTRG method. Finally, we study the length scale induced by the finite cutoff in bond dimension and elucidate its physical meaning in this context.