Monte Carlo determination of multiple extremal eigenpairs
Abstract
We present a Monte Carlo algorithm that allows the simultaneous determination of a few extremal eigenpairs of a very large matrix without the need to compute the inner product of two vectors or store all the components of any one vector. The algorithm, a Monte Carlo implementation of a deterministic one we recently benchmarked, is an extension of the power method. In the implementation presented, we used a basic Monte Carlo splitting and termination method called the comb, incorporated the weight cancellation method of Arnow , and exploited a sampling method, the sewing method, that does a large state space sampling as a succession of small state space samplings. We illustrate the effectiveness of the algorithm by its determination of the two largest eigenvalues of the transfer matrices for variously sized two-dimensional, zero-field Ising models. While very likely useful for other transfer-matrix problems, the algorithm is however quite general and should find application to a larger variety of problems requiring a few dominant eigenvalues of a matrix.
- Publication:
-
Physical Review E
- Pub Date:
- October 2009
- DOI:
- 10.1103/PhysRevE.80.046704
- arXiv:
- arXiv:0812.4854
- Bibcode:
- 2009PhRvE..80d6704B
- Keywords:
-
- 02.70.Tt;
- 05.10.Ln;
- Justifications or modifications of Monte Carlo methods;
- Monte Carlo methods;
- Physics - Computational Physics;
- Condensed Matter - Statistical Mechanics
- E-Print:
- 22 pages, no figures