Computation of Dominant Eigenvalues and Eigenvectors: A Comparative Study of Algorithms

We investigate two widely used recursive algorithms for the computation of eigenvectors with extreme eigenvalues of large symmetric matrices -- the modified Lanczös method and the conjugate-gradient method. The goal is to establish a connection between their underlying principles and to evaluate their performance in applications to Hamiltonian and transfer matrices of selected model systems of interest in condensed matter physics and statistical mechanics. The conjugate-gradient method is found to converge more rapidly for understandable reasons, while storage requirements are the same for both methods

Improved Phenomenological Renormalization Schemes

An analysis is made of various methods of phenomenological renormalization based on finite-size scaling equations for inverse correlation lengths, the singular part of the free energy density, and their derivatives. The analysis is made using two-dimensional Ising and Potts lattices and the three-dimensional Ising model. Variants of equations for the phenomenological renormalization group are obtained which ensure more rapid convergence than the conventionally used Nightingale phenomenological renormalization scheme. An estimate is obtained for the critical finite-size scaling amplitude of the internal energy in the three-dimensional Ising model. It is shown that the two-dimensional Ising and Potts models contain no finite-size corrections to the internal energy so that the positions of the critical points for these models can be determined exactly from solutions for strips of finite width. It is also found that for the two-dimensional Ising model the scaling finite-size equation for the derivative of the inverse correlation length with respect to temperature gives the exact value of the thermal critical exponent.Comment: 14 pages with 1 figure in late

The Dynamic Exponent of the Two-Dimensional Ising Model and Monte Carlo Computation of the Sub-Dominant Eigenvalue of the Stochastic Matrix

We introduce a novel variance-reducing Monte Carlo algorithm for accurate determination of autocorrelation times. We apply this method to two-dimensional Ising systems with sizes up to $15 \times 15$, using single-spin flip dynamics, random site selection and transition probabilities according to the heat-bath method. From a finite-size scaling analysis of these autocorrelation times, the dynamical critical exponent $z$ is determined as $z=2.1665$ (12)

Critical temperature of a fully anisotropic three-dimensional Ising model

The critical temperature of a three-dimensional Ising model on a simple cubic lattice with different coupling strengths along all three spatial directions is calculated via the transfer matrix method and a finite size scaling for L x L oo clusters (L=2 and 3). The results obtained are compared with available calculations. An exact analytical solution is found for the 2 x 2 oo Ising chain with fully anisotropic interactions (arbitrary J_x, J_y and J_z).Comment: 17 pages in tex using preprint.sty for IOP journals, no figure

Comment on "Two Phase Transitions in the Fully frustrated XY Model"

The conclusions of a recent paper by Olsson (Phys. Rev. Lett. 75, 2758 (1995), cond-mat/9506082) about the fully frustrated XY model in two dimensions are questioned. In particular, the evidence presented for having two separate chiral and U(1) phase transitions are critically considered.Comment: One page one table, to Appear in Physical Review Letter

High precision Monte Carlo study of the 3D XY-universality class

We present a Monte Carlo study of the two-component $\phi^4$ model on the simple cubic lattice in three dimensions. By suitable tuning of the coupling constant $\lambda$ we eliminate leading order corrections to scaling. High statistics simulations using finite size scaling techniques yield $\nu=0.6723(3)[8]$ and $\eta=0.0381(2)[2]$, where the statistical and systematical errors are given in the first and second bracket, respectively. These results are more precise than any previous theoretical estimate of the critical exponents for the 3D XY universality class.Comment: 13 page

The phase diagram of the anisotropic Spin-1 Heisenberg Chain

We applied the Density Matrix Renormalization Group to the XXZ spin-1 quantum chain. In studing this model we aim to clarify controversials about the point where the massive Haldane phase appears.Comment: 2 pages (standart LaTex), 1 figure (PostScript) uuencode

High-precision estimate of g4 in the 2D Ising model

We compute the renormalized four-point coupling in the 2d Ising model using transfer-matrix techniques. We greatly reduce the systematic uncertainties which usually affect this type of calculations by using the exact knowledge of several terms in the scaling function of the free energy. Our final result is g4=14.69735(3).Comment: 17 pages, revised version with minor changes, accepted for publication in Journal of Physics

Random walks near Rokhsar-Kivelson points

There is a class of quantum Hamiltonians known as Rokhsar-Kivelson(RK)-Hamiltonians for which static ground state properties can be obtained by evaluating thermal expectation values for classical models. The ground state of an RK-Hamiltonian is known explicitly, and its dynamical properties can be obtained by performing a classical Monte Carlo simulation. We discuss the details of a Diffusion Monte Carlo method that is a good tool for studying statics and dynamics of perturbed RK-Hamiltonians without time discretization errors. As a general result we point out that the relation between the quantum dynamics and classical Monte Carlo simulations for RK-Hamiltonians follows from the known fact that the imaginary-time evolution operator that describes optimal importance sampling, in which the exact ground state is used as guiding function, is Markovian. Thus quantum dynamics can be studied by a classical Monte Carlo simulation for any Hamiltonian that is free of the sign problem provided its ground state is known explicitly.Comment: 12 pages, 9 figures, RevTe

Finite-size scaling corrections in two-dimensional Ising and Potts ferromagnets

Finite-size corrections to scaling of critical correlation lengths and free energies of Ising and three-state Potts ferromagnets are analysed by numerical methods, on strips of width $N$ sites of square, triangular and honeycomb lattices. Strong evidence is given that the amplitudes of the ``analytical'' correction terms, $N^{-2}$, are identically zero for triangular-- and honeycomb Ising systems. For Potts spins, our results are broadly consistent with this lattice-dependent pattern of cancellations, though for correlation lengths non-vanishing (albeit rather small) amplitudes cannot be entirely ruled out.Comment: 11 pages, LaTeX with Institute of Physics macros, 2 EPS figures; to appear in Journal of Physics