1,793 research outputs found
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
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 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 sites of square, triangular and honeycomb
lattices. Strong evidence is given that the amplitudes of the ``analytical''
correction terms, , 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
Quantum Speedup by Quantum Annealing
We study the glued-trees problem of Childs et. al. in the adiabatic model of
quantum computing and provide an annealing schedule to solve an oracular
problem exponentially faster than classically possible. The Hamiltonians
involved in the quantum annealing do not suffer from the so-called sign
problem. Unlike the typical scenario, our schedule is efficient even though the
minimum energy gap of the Hamiltonians is exponentially small in the problem
size. We discuss generalizations based on initial-state randomization to avoid
some slowdowns in adiabatic quantum computing due to small gaps.Comment: 7 page
Numerical Studies of the Two Dimensional XY Model with Symmetry Breaking Fields
We present results of numerical studies of the two dimensional XY model with
four and eight fold symmetry breaking fields. This model has recently been
shown to describe hydrogen induced reconstruction on the W(100) surface. Based
on mean-field and renormalization group arguments,we first show how the
interplay between the anisotropy fields can give rise to different phase
transitions in the model. When the fields are compatible with each other there
is a continuous phase transition when the fourth order field is varied from
negative to positive values. This transition becomes discontinuous at low
temperatures. These two regimes are separated by a multicritical point. In the
case of competing four and eight fold fields, the first order transition at low
temperatures opens up into two Ising transitions. We then use numerical methods
to accurately locate the position of the multicritical point, and to verify the
nature of the transitions. The different techniques used include Monte Carlo
histogram methods combined with finite size scaling analysis, the real space
Monte Carlo Renormalization Group method, and the Monte Carlo Transfer Matrix
method. Our numerical results are in good agreement with the theoretical
arguments.Comment: 29 pages, HU-TFT-94-36, to appear in Phys. Rev. B, Vol 50, November
1, 1994. A LaTeX file with no figure
- âŠ