A new Monte Carlo power method for the eigenvalue problem of transfer matrices |
| |
Authors: | Tohru Koma |
| |
Institution: | (1) Department of Physics, Gakushuin University, Mejiro, Toshima-ku, 171 Tokyo, Japan |
| |
Abstract: | We propose a new Monte Carlo method for calculating eigenvalues of transfer matrices leading to free energies and to correlation lengths of classical and quantum many-body systems. Generally, this method can be applied to the calculation of the maximum eigenvalue of a nonnegative matrix  such that all the matrix elements of Âk are strictly positive for an integerk. This method is based on a new representation of the maximum eigenvalue of the matrix  as the thermal average of a certain observable of a many-body system. Therefore one can easily calculate the maximum eigenvalue of a transfer matrix leading to the free energy in the standard Monte Carlo simulations, such as the Metropolis algorithm. As test cases, we calculate the free energies of the square-lattice Ising model and of the spin-1/2XY Heisenberg chain. We also prove two useful theorems on the ergodicity in quantum Monte Carlo algorithms, or more generally, on the ergodicity of Monte Carlo algorithms using our new representation of the maximum eigenvalue of the matrixÂ. |
| |
Keywords: | Monte Carlo simulations power method eigenvalue problems transfer matrices free energy calculation |
本文献已被 SpringerLink 等数据库收录! |
|