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

Physical review. B, Condensed matter

Monte Carlo renormalization-group analysis of the classical Heisenberg model in two dimensions

Stephen H. Shenker, Jan Tobochnik

Published Nov 1, 1980• 259 citations
Publisher

Abstract

Monte Carlo renormalization-group methods were applied to the classical three-component Heisenberg model on a two-dimensional lattice. Expectation values of local correlations of spins and various sized block spins were computed using traditional Monte Carlo methods. By matching quantities at different length scales generated by different Monte Carlo Hamiltonians we directly determined the renormalization of the nearest-neighbor coupling. Using these data and the results of high-temperature expansions and low-temperature renormalization-group calculations we have determined that this model does not have a phase transition. We have also obtained the amplitude for the low-temperature divergence of the susceptibility and the correlation length.

DOI
10.1103/physrevb.22.4462

Related papers

Low-temperature behavior of the correlation length and the susceptibility of a quantum Heisenberg ferromagnet in two dimensions

Peter Kopietz, Sudip Chakravarty

The low-temperature behavior of quantum Heisenberg models can be obtained by combining a renormalization-group analysis of the quantum models with a Monte Carlo simulation of the corresponding classical model, where the ``spins'' are represented by unit vectors defined on a lattice. In this approach renormalization-group equations for the quantum models lead to an effective classical O(3) nonlinear \ensuremath{\sigma} model in appropriate regimes of the parameter space. In the present work we apply this technique, which was recently used by Chakravarty, Halperin, and Nelson to discuss the low-temperature behavior of quantum antiferromagnets, to ferromagnets and calculate the correlation length and the susceptibility of a two-dimensional quantum ferromagnet at low temperatures. The results obtained with this method are then compared with those obtained from a direct simulation of the spin S=1/2 quantum ferromagnet. Excellent agreement between the two approaches leads us to believe that the analytic expressions obtained in the renormalization-group method are accurate at low temperatures. This strengthens the validity of the procedure adopted for antiferromagnets as well, where no such simulation results are yet known. The present paper also further substantiates recent work on the properties of solid $^{3}\mathrm{He}$ in two dimensions.

Physical review. B, Condensed matter•Sep 1, 1989•63 citations

Low-temperature behavior of the correlation length and the susceptibility of the ferromagnetic quantum Heisenberg chain

Peter Kopietz

The renormalization-group recursion relations for the spin-S nearest-neighbor quantum Heisenberg ferromagnet recently derived by Chakravarty and the present author are analyzed in one dimension. It is found that the susecptibility \ensuremath{\chi} is given at low temperatures T by T\ensuremath{\chi}=${C}_{\ensuremath{\chi}}^{(0)}$(${\mathrm{JS}}^{2}$/T) [1+${C}_{\ensuremath{\chi}}^{(1)}$\ensuremath{\eta}+O(${\ensuremath{\eta}}^{2}$)], where J is the exchange coupling and the dimensionless expansion parameter \ensuremath{\eta} is given by \ensuremath{\eta}=${\ensuremath{\pi}}^{\mathrm{\ensuremath{-}}1}$(T/${\mathrm{JS}}^{3}$${)}^{1/2}$. For the correlation length \ensuremath{\xi} we find (\ensuremath{\xi}/a)=${C}_{\ensuremath{\xi}}^{(0)}$(${\mathrm{JS}}^{2}$/T) [1+${C}_{\ensuremath{\xi}}^{(1)}$\ensuremath{\eta}+O(${\ensuremath{\eta}}^{2}$)], where a is the lattice spacing. The values of the coefficients ${C}_{\ensuremath{\chi}}^{(j)}$ and ${C}_{\ensuremath{\xi}}^{(j)}$ (j=0,1), have been determined by fitting the above expressions with Monte Carlo data obtained for the S=1/2 Heisenberg chain at temperatures as low as T/J=0.01.

Physical review. B, Condensed matter•Sep 1, 1989•17 citations

Block-spin renormalization group in the large- n limit

J. E. Hirsch, Stephen H. Shenker

We present a real-space renormalization-group (RG) study of the $n$-component classical Heisenberg model in the large-$n$ limit. We obtain exact expressions for correlation functions of block spins, given those of the original spins. We show how to extract the critical behavior of the model from such information. The model provides a good testing ground for evaluating the effectiveness of Monte Carlo RG techniques. We study numerically the implementation of the real-space RG transformation on finite lattices. In two dimensions, we obtain the $\ensuremath{\beta}$ function as a function of temperature. In three dimensions, we calculate the critical temperature and the thermal exponent. We discuss the effects of the use of finite lattices and of the use of different forms for the initial action on the accuracy of the results.

Physical review. B, Condensed matter•Feb 1, 1983•36 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

Critical exponents of the classical three-dimensional Heisenberg model: A single-cluster Monte Carlo study

Christian Holm, Wolfhard Janke E

We have simulated the three-dimensional Heisenberg model on simple cubic lattices, using the single-cluster Monte Carlo update algorithm. The expected pronounced reduction of critical slowing down at the phase transition is verified. This allows simulations on significantly larger lattices than in previous studies and consequently a better control over systematic errors. In one set of simulations we employ the usual finite-size scaling methods to compute the critical exponents \ensuremath{\nu},\ensuremath{\alpha},\ensuremath{\beta},\ensuremath{\gamma},\ensuremath{\eta} from a few measurements in the vicinity of the critical point, making extensive use of histogram reweighting and optimization techniques. In another set of simulations we report measurements of improved estimators for the spatial correlation length and the susceptibility in the high-temperature phase, obtained on lattices with up to ${100}^{3}$ spins. This enables us to compute independent estimates of \ensuremath{\nu} and \ensuremath{\gamma} from power-law fits of their critical divergencies.

arXiv:hep-lat/9301002•Physical review. B, Condensed matter•Jul 1, 1993•232 citations

Monte Carlo renormalization-group calculations of critical behavior in the simple-cubic Ising model

Godfrey Stuart Pawley, Robert H. Swendsen, D. J. Wallace +1 more

The Monte Carlo renormalization group is applied to the three-dimensional Ising model on simple cubic lattices with ${8}^{3}$, ${16}^{3}$, ${32}^{3}$, and ${64}^{3}$ sites. The comparison of block-spin correlation functions from the largest lattices yields the nearest-neighbor critical coupling ${K}_{1}^{c}=0.221654(6)$. After allowing for (i) interpolation to this best estimate for ${K}_{1}^{c}$, (ii) an apparent finite-size effect in the renormalization-group transformation due to the measurement of correlation functions of too few (seven) operators, and (iii) the extrapolation for the effect of a slow transient towards the fixed point, the values $\ensuremath{\nu}=0.629(4)$ and $\ensuremath{\eta}=0.031(5)$ are obtained for the thermal and magnetic exponents. The correction-to-scaling exponent $\ensuremath{\omega}$ is estimated to be around 1; to obtain an accuracy competitive with other methods requires measurements with more than seven operators. We briefly review the problem of redundant operators and indicate the future prospects for this kind of calculation.

Physical review. B, Condensed matter•Apr 1, 1984•318 citations
Monte Carlo renormalization-group methods were applied to the classical three-component Heisenberg model on a two-dimensional lattice. Expectation values of local correlations of spins and various sized block spins were computed using traditional Monte Carlo methods. By matching quantities at different length scales generated by different Monte Carlo Hamiltonians we directly determined the renormalization of the nearest-neighbor coupling. Using these data and the results of high-temperature expansions and low-temperature renormalization-group calculations we have determined that this model does not have a phase transition. We have also obtained the amplitude for the low-temperature divergence of the susceptibility and the correlation length.
The low-temperature behavior of quantum Heisenberg models can be obtained by combining a renormalization-group analysis of the quantum models with a Monte Carlo simulation of the corresponding classical model, where the ``spins'' are represented by unit vectors defined on a lattice. In this approach renormalization-group equations for the quantum models lead to an effective classical O(3) nonlinear \ensuremath{\sigma} model in appropriate regimes of the parameter space. In the present work we apply this technique, which was recently used by Chakravarty, Halperin, and Nelson to discuss the low-temperature behavior of quantum antiferromagnets, to ferromagnets and calculate the correlation length and the susceptibility of a two-dimensional quantum ferromagnet at low temperatures. The results obtained with this method are then compared with those obtained from a direct simulation of the spin S=1/2 quantum ferromagnet. Excellent agreement between the two approaches leads us to believe that the analytic expressions obtained in the renormalization-group method are accurate at low temperatures. This strengthens the validity of the procedure adopted for antiferromagnets as well, where no such simulation results are yet known. The present paper also further substantiates recent work on the properties of solid 3He^{3}\mathrm{He}3He in two dimensions.
The renormalization-group recursion relations for the spin-S nearest-neighbor quantum Heisenberg ferromagnet recently derived by Chakravarty and the present author are analyzed in one dimension. It is found that the susecptibility \ensuremath{\chi} is given at low temperatures T by T\ensuremath{\chi}=${C}_{\ensuremath{\chi}}^{(0)}$(JS2{\mathrm{JS}}^{2}JS2/T) [1+${C}_{\ensuremath{\chi}}^{(1)}$\ensuremath{\eta}+O(${\ensuremath{\eta}}^{2}$)], where J is the exchange coupling and the dimensionless expansion parameter \ensuremath{\eta} is given by \ensuremath{\eta}=${\ensuremath{\pi}}^{\mathrm{\ensuremath{-}}1}$(T/JS3{\mathrm{JS}}^{3}JS3)1/2{)}^{1/2})1/2. For the correlation length \ensuremath{\xi} we find (\ensuremath{\xi}/a)=${C}_{\ensuremath{\xi}}^{(0)}$(JS2{\mathrm{JS}}^{2}JS2/T) [1+${C}_{\ensuremath{\xi}}^{(1)}$\ensuremath{\eta}+O(${\ensuremath{\eta}}^{2}$)], where a is the lattice spacing. The values of the coefficients ${C}_{\ensuremath{\chi}}^{(j)}$ and ${C}_{\ensuremath{\xi}}^{(j)}$ (j=0,1), have been determined by fitting the above expressions with Monte Carlo data obtained for the S=1/2 Heisenberg chain at temperatures as low as T/J=0.01.
We present a real-space renormalization-group (RG) study of the nnn-component classical Heisenberg model in the large-nnn limit. We obtain exact expressions for correlation functions of block spins, given those of the original spins. We show how to extract the critical behavior of the model from such information. The model provides a good testing ground for evaluating the effectiveness of Monte Carlo RG techniques. We study numerically the implementation of the real-space RG transformation on finite lattices. In two dimensions, we obtain the $\ensuremath{\beta}$ function as a function of temperature. In three dimensions, we calculate the critical temperature and the thermal exponent. We discuss the effects of the use of finite lattices and of the use of different forms for the initial action on the accuracy of the results.
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.
We have simulated the three-dimensional Heisenberg model on simple cubic lattices, using the single-cluster Monte Carlo update algorithm. The expected pronounced reduction of critical slowing down at the phase transition is verified. This allows simulations on significantly larger lattices than in previous studies and consequently a better control over systematic errors. In one set of simulations we employ the usual finite-size scaling methods to compute the critical exponents \ensuremath{\nu},\ensuremath{\alpha},\ensuremath{\beta},\ensuremath{\gamma},\ensuremath{\eta} from a few measurements in the vicinity of the critical point, making extensive use of histogram reweighting and optimization techniques. In another set of simulations we report measurements of improved estimators for the spatial correlation length and the susceptibility in the high-temperature phase, obtained on lattices with up to 1003{100}^{3}1003 spins. This enables us to compute independent estimates of \ensuremath{\nu} and \ensuremath{\gamma} from power-law fits of their critical divergencies.
The Monte Carlo renormalization group is applied to the three-dimensional Ising model on simple cubic lattices with 83{8}^{3}83, 163{16}^{3}163, 323{32}^{3}323, and 643{64}^{3}643 sites. The comparison of block-spin correlation functions from the largest lattices yields the nearest-neighbor critical coupling K1c=0.221654(6){K}_{1}^{c}=0.221654(6)K1c​=0.221654(6). After allowing for (i) interpolation to this best estimate for K1c{K}_{1}^{c}K1c​, (ii) an apparent finite-size effect in the renormalization-group transformation due to the measurement of correlation functions of too few (seven) operators, and (iii) the extrapolation for the effect of a slow transient towards the fixed point, the values $\ensuremath{\nu}=0.629(4)$ and $\ensuremath{\eta}=0.031(5)$ are obtained for the thermal and magnetic exponents. The correction-to-scaling exponent $\ensuremath{\omega}$ is estimated to be around 1; to obtain an accuracy competitive with other methods requires measurements with more than seven operators. We briefly review the problem of redundant operators and indicate the future prospects for this kind of calculation.