Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Investigation of One- and Two-Dimensional Quantum Spin Systems by Monte Carlo Simulations

Luruper Chaussee

Published Jan 1, 2000• 0 citations
OpenAlex

Abstract

A method for simulating quantum spin systems with nearest neighbour interaction is reviewed. In the framework of a path-integral approach based on the Trotter formula, the quantum spin system is mapped into the limit of a sequence of classical spin systems with one extra dimension. Numerical questions connected to this limit are discussed. The role of particle and winding numbers is clarified. An overview of the problems investigated by our group is given. In particular, two new results are presented. A simulation of the one-dimensional spin S isotropic Heisenberg antiferromagnet, S=1/2, 1, 3/2, 2, supports Haldane's conjecture that the mass gap is zero for half­ integer S and nonzero for integer S. For the two-dimensional spin 1/2 xy model we discuss evidence for a Kosterlitz-Thouless transition. In particular, we show that the helicity modulus can be computed directly from our simulaton.

Related papers

Quantum Monte Carlo simulation of the one-dimensional spin S xxz model. I. The method

M Marcu, Andreas Wiesler

The authors develop a quantum Monte Carlo procedure for the xxz spin chain. The quantum model in one dimension is mapped by a path integral method based on the Trotter formula into the infinite anisotropy limit of a classical two-dimensional spin system. The expectation values of the energy density, susceptibility and static s z -s z structure function are computed in the classical model for finite but increasing anisotropy. The computation method is a Monte Carlo simulation of the classical system. The results are systematically extrapolated to the infinite anisotropy limit.

Journal of Physics A Mathematical and General•Sep 11, 1985•13 citations

MONTE CARLO STUDY OF TIIE TWO-DIMENSIONAL SPIN-I/2 XY MODEL ~

Hans A. De Raedt, Bart De Raedt

The generalized Trotter formula is used to map the two-dimensional spin-l/2 XY model onto a three-dimensional lsing model with complicated many-spin interactions. This relationship sused to construct an efficient Monte Carlo algorithm. Simulation data for the specific heat and vortex correlation functions give strong evidence for the existence of a phase transition. Although it is well known that there cannot exist long-range order in the two-dimensional (2D) spin-1/2 XY model at non-zero temperature [1], high-temper-ature series expansions [2] suggest he existence of critical behavior. Little is known about the quantum model mainly because most techniques which have been succesfully applied to the classical model (also referred to as planar rotator model) are very difficult to extend to the quantum case. In particular, quantum renormalization group analYsis of the 2D spin-1/2 XY model have given either contradictive or unreliable re-suits [3-8]. An alternative method to calculate the thermody-namic properties of this modelis to exploit the formal analogy between d-dimensional quantum spin systems and (d + 1)-dimensional Ising-spin models [9]. Using ~ ' This paper was first submitted to Physical Review Letters

•Jan 1, 1983•0 citations

Monte Carlo study of quantum spin chains

S. V. Meshkov

Results of a numerical experiment on one-dimensional quantum spin chains are reported. The Heisenberg model with a nearest-neighbor interaction is investigated for spins 1/2\ensuremath{\le}S\ensuremath{\le}2. The procedure of measurement in the conventional checkerboard world-line algorithm is modified to give lower systematic errors. In addition to spatial-correlation properties, the magnon spectral density is extracted from the magnon Green function via the fast maximum-entropy algorithm. In particular, the existence of the Haldane gap in S=1 and S=2 cases is explicitly demonstrated. Hidden antiferromagnetic order peculiar to the Haldane phase is measured for S=1 as a function of one-site anisotropy. The mapping onto the restricted solid-on-solid model is analyzed. The difference between systems with even and odd integer spins is discussed.

Physical review. B, Condensed matter•Sep 1, 1993•55 citations

Monte Carlo studies of two dimensional quantum spin systems

M. S. Makivić

Spin-1/2 nearest neighbor Heisenberg antiferromagnet and XY model on a square lattice are studied via large scale quantum Monte Carlo simulations using a fast and efficient multispin coding algorithm on the Caltech/JPL MarkIIIfp parallel supercomputer, based on the Suzuki-Trotter transformation. We performed simulations with very good statistics on lattices as large as 128x128 spins, in the temperature range from 0.1 to 2.5 in units of the effective exchange coupling J. We calculated energy, specific heat, magnetic susceptibilities and also spin correlation functions from which we deduce the correlation lengths. For the Heisenberg model, at temperatures higher than J the results are in excellent agreement with high-temperature series expansion. At low temperatures the long wavelength behavior is essentially classical. Our data show that the correlation length and staggered susceptibility are quantitatively well described by the renormalized classical picture at the 2-loop level of approximation. From the divergence of the correlation length, we deduce the value of the quantum renormalized spin stiffness, [rho][subscript s]/J = 0.199(2). We give evidence that the correlation function is of the Ornstein-Zernicke type. By comparing the largest measured correlation lengths with neutron scattering experiments on La2CuO4, we deduce the value of effective exchange coupling J = 1450±30 K. By measuring the imaginary time-dependent correlation functions, we show that the dynamics of the model can be well understood within a Bose liquid-type picture. The spin waves are rather sharp throughout most of the Brillouin zone and the damping is weakly dependent on the wave vector. In the case of the XY model, convincing numerical evidence is obtained on square lattices as large as 96x96 that the spin-1/2 XY model undergoes a Kosterlitz-Thouless (KT) phase transition at kTc/J = 0.350(4). The correlation length and in-plane susceptibility diverge at Tc precisely according to the form predicted by Kosterlitz and Thouless for the classical XY model. The specific heat increases very rapidly on heating near Tc and exhibits a peak around kT/J = 0.45. We also measure the spin stiffness and the correlation function exponent below the transition temperature. Within the statistical accuracy of the measurements, the results are well described by the square root singularity (with a nonuniversal amplitude) below Tc, and they have the universal values in agreement with KT theory at Tc.

•Jul 20, 2007•0 citations

Monte Carlo simulation of the one-dimensional t-J-J' model

Matthias Troyer

I present a new and efficient disconnected cluster decomposition to perform Monte Carlo Simulations based on the Trotter-Suzuki method for 1-d quantum systems. As opposed to the normal checkerboard world line algorithm this mapping allows one to simulate models with interactions up to next nearest neighbours or even one-dimensional double layers. I have used this cluster decomposition to simulate the one-dimensional t-J-J' model at α=J'/J=0.5 and have investigated the sign problem in this model. At half band filling the system is dimerised into singlet spin pairs when α=0.5 due to frustration. I believe that the spin gap is no longer present around quarter band filling

•Jul 5, 1991•0 citations

Meron-cluster simulation of the quantum antiferromagnetic Heisenberg model in a magnetic field in one- and two-dimensions

Palma, A. Riveros

Motivated by the numerical simulation of systems which display quantum phase transitions, we present a novel application of the meron-cluster algorithm to simulate the quantum antiferromagnetic Heisenberg model coupled to an external uniform magnetic field both in one and in two dimensions.In the infinite volume limit and at zero temperature we found numerical evidence that supports a quantum phase transition very close to the critical values B c = 2 and B c = 4 for the system in one and two dimensions, respectively.For the one dimensional system, we have compared the numerical data obtained with analytical predictions for the magnetization density as a function of the external field obtained by scaling-behaviour analysis and Bethe Ansatz techniques.Since there is no analytical solution for the two dimensional case, we have compared our results with the magnetization density obtained by scaling relations for small lattice sizes and with the approximated thermodynamical limit at zero temperature guessed by scaling relations.Moreover, we have compared the numerical data with other numerical simulations performed by using different algorithms in one and two dimensions, like the directed loop method.The numerical data obtained are in perfect agreement with all these previous results, which confirms that the meron-algorithm is reliable for quantum Monte Carlo simulations and applicable both in one and two dimensions.Finally, we have computed the integrated autocorrelation time to measure the efficiency of the meron algorithm in one dimension.

arXiv:1312.4915•Condensed Matter Physics•May 25, 2015•8 citations
A method for simulating quantum spin systems with nearest neighbour interaction is reviewed. In the framework of a path-integral approach based on the Trotter formula, the quantum spin system is mapped into the limit of a sequence of classical spin systems with one extra dimension. Numerical questions connected to this limit are discussed. The role of particle and winding numbers is clarified. An overview of the problems investigated by our group is given. In particular, two new results are presented. A simulation of the one-dimensional spin S isotropic Heisenberg antiferromagnet, S=1/2, 1, 3/2, 2, supports Haldane's conjecture that the mass gap is zero for half­ integer S and nonzero for integer S. For the two-dimensional spin 1/2 xy model we discuss evidence for a Kosterlitz-Thouless transition. In particular, we show that the helicity modulus can be computed directly from our simulaton.
The authors develop a quantum Monte Carlo procedure for the xxz spin chain. The quantum model in one dimension is mapped by a path integral method based on the Trotter formula into the infinite anisotropy limit of a classical two-dimensional spin system. The expectation values of the energy density, susceptibility and static s z -s z structure function are computed in the classical model for finite but increasing anisotropy. The computation method is a Monte Carlo simulation of the classical system. The results are systematically extrapolated to the infinite anisotropy limit.
The generalized Trotter formula is used to map the two-dimensional spin-l/2 XY model onto a three-dimensional lsing model with complicated many-spin interactions. This relationship sused to construct an efficient Monte Carlo algorithm. Simulation data for the specific heat and vortex correlation functions give strong evidence for the existence of a phase transition. Although it is well known that there cannot exist long-range order in the two-dimensional (2D) spin-1/2 XY model at non-zero temperature [1], high-temper-ature series expansions [2] suggest he existence of critical behavior. Little is known about the quantum model mainly because most techniques which have been succesfully applied to the classical model (also referred to as planar rotator model) are very difficult to extend to the quantum case. In particular, quantum renormalization group analYsis of the 2D spin-1/2 XY model have given either contradictive or unreliable re-suits [3-8]. An alternative method to calculate the thermody-namic properties of this modelis to exploit the formal analogy between d-dimensional quantum spin systems and (d + 1)-dimensional Ising-spin models [9]. Using ~ ' This paper was first submitted to Physical Review Letters
Results of a numerical experiment on one-dimensional quantum spin chains are reported. The Heisenberg model with a nearest-neighbor interaction is investigated for spins 1/2\ensuremath{\le}S\ensuremath{\le}2. The procedure of measurement in the conventional checkerboard world-line algorithm is modified to give lower systematic errors. In addition to spatial-correlation properties, the magnon spectral density is extracted from the magnon Green function via the fast maximum-entropy algorithm. In particular, the existence of the Haldane gap in S=1 and S=2 cases is explicitly demonstrated. Hidden antiferromagnetic order peculiar to the Haldane phase is measured for S=1 as a function of one-site anisotropy. The mapping onto the restricted solid-on-solid model is analyzed. The difference between systems with even and odd integer spins is discussed.
Spin-1/2 nearest neighbor Heisenberg antiferromagnet and XY model on a square lattice are studied via large scale quantum Monte Carlo simulations using a fast and efficient multispin coding algorithm on the Caltech/JPL MarkIIIfp parallel supercomputer, based on the Suzuki-Trotter transformation. We performed simulations with very good statistics on lattices as large as 128x128 spins, in the temperature range from 0.1 to 2.5 in units of the effective exchange coupling J. We calculated energy, specific heat, magnetic susceptibilities and also spin correlation functions from which we deduce the correlation lengths. For the Heisenberg model, at temperatures higher than J the results are in excellent agreement with high-temperature series expansion. At low temperatures the long wavelength behavior is essentially classical. Our data show that the correlation length and staggered susceptibility are quantitatively well described by the renormalized classical picture at the 2-loop level of approximation. From the divergence of the correlation length, we deduce the value of the quantum renormalized spin stiffness, [rho][subscript s]/J = 0.199(2). We give evidence that the correlation function is of the Ornstein-Zernicke type. By comparing the largest measured correlation lengths with neutron scattering experiments on La2CuO4, we deduce the value of effective exchange coupling J = 1450±30 K. By measuring the imaginary time-dependent correlation functions, we show that the dynamics of the model can be well understood within a Bose liquid-type picture. The spin waves are rather sharp throughout most of the Brillouin zone and the damping is weakly dependent on the wave vector. In the case of the XY model, convincing numerical evidence is obtained on square lattices as large as 96x96 that the spin-1/2 XY model undergoes a Kosterlitz-Thouless (KT) phase transition at kTc/J = 0.350(4). The correlation length and in-plane susceptibility diverge at Tc precisely according to the form predicted by Kosterlitz and Thouless for the classical XY model. The specific heat increases very rapidly on heating near Tc and exhibits a peak around kT/J = 0.45. We also measure the spin stiffness and the correlation function exponent below the transition temperature. Within the statistical accuracy of the measurements, the results are well described by the square root singularity (with a nonuniversal amplitude) below Tc, and they have the universal values in agreement with KT theory at Tc.
I present a new and efficient disconnected cluster decomposition to perform Monte Carlo Simulations based on the Trotter-Suzuki method for 1-d quantum systems. As opposed to the normal checkerboard world line algorithm this mapping allows one to simulate models with interactions up to next nearest neighbours or even one-dimensional double layers. I have used this cluster decomposition to simulate the one-dimensional t-J-J' model at α=J'/J=0.5 and have investigated the sign problem in this model. At half band filling the system is dimerised into singlet spin pairs when α=0.5 due to frustration. I believe that the spin gap is no longer present around quarter band filling
Motivated by the numerical simulation of systems which display quantum phase transitions, we present a novel application of the meron-cluster algorithm to simulate the quantum antiferromagnetic Heisenberg model coupled to an external uniform magnetic field both in one and in two dimensions.In the infinite volume limit and at zero temperature we found numerical evidence that supports a quantum phase transition very close to the critical values B c = 2 and B c = 4 for the system in one and two dimensions, respectively.For the one dimensional system, we have compared the numerical data obtained with analytical predictions for the magnetization density as a function of the external field obtained by scaling-behaviour analysis and Bethe Ansatz techniques.Since there is no analytical solution for the two dimensional case, we have compared our results with the magnetization density obtained by scaling relations for small lattice sizes and with the approximated thermodynamical limit at zero temperature guessed by scaling relations.Moreover, we have compared the numerical data with other numerical simulations performed by using different algorithms in one and two dimensions, like the directed loop method.The numerical data obtained are in perfect agreement with all these previous results, which confirms that the meron-algorithm is reliable for quantum Monte Carlo simulations and applicable both in one and two dimensions.Finally, we have computed the integrated autocorrelation time to measure the efficiency of the meron algorithm in one dimension.