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

Proceedings of XXIIIrd International Symposium on Lattice Field Theory — PoS(LAT2005)

QMC simulations of Heisenberg ferromagnet

L. Bogacz, Wolfhard Janke E

Published Dec 14, 2005• 0 citations
Publisher

Abstract

We use the Quantum Monte Carlo Stochastic Series Expansion (SSE) algorithm to investigate the properties of the one-dimensional Heisenberg ferromagnet.We consider the model for two spin values s = 1 2 and s = 1 and study the dependence of the specific heat, magnetisation, magnetic susceptibility and correlation functions on temperature and magnetic field.Obtained results are compared to those derived by analytical methods (Bethe ansatz, Green's function, Random Phase Approximation).In particular, the surprising prediction of a double-peak structure in the specific heat at low temperatures and small magnetic fields can be confirmed.

DOI
10.22323/1.020.0241

Related papers

Investigating the stochastic series expansion method for one-dimensional spin-1/2 antiferromagnetic Heisenberg model

Ali Akbar Hasan, Mohammad G. Merdan, Hayder M. Abdualjalil

We investigate the Stochastic Series Expansion (SSE) technique, guided by Sandvik's algorithm of the quantum Monte Carlo simulation for spin-1/2 antiferromagnetic Heisenberg linear chain. The impact of the Monte Carlo steps, lattice sizes, and various temperatures on the variability and fidelity of simulation outcomes are studied. The findings show that the Monte Carlo steps where measurements are made (Msteps), and lattice size are crucial for stabilizing the results. Our results show that for larger lattice sizes, such as L=1024, achieving reliable outcomes necessitates extensive Monte Carlo simulations, specifically. The convergence in the properties is achieved with a high number of Monte Carlo steps starting from 7×10^4. Strong fluctuations are calculated between the results for each magnetic property and the existence of negative and positive in the specific heat at low temperatures.

AIP conference proceedings•Jan 1, 2026•0 citations

Quantum Monte Carlo simulation of tow-dimensional antiferromagnetic Heisenberg model

Jinbin Li

The Antiferromagnetic Heisenberg Model with spin 1/2 on a square lattice was simulated by using the Stochastic Series Expansion(SSE) Quantum Monte Carlo Method.The temperature dependence of the internal energy,specific heat and uniform susceptibility was observed.The results showed that specific heat has a maximum at kT/J=0.6 and uniform susceptibility saturates at kT/J=1 for the isotropy case.And the size of lattice also affects the temperature dependence of these thermodynamic quantities.For the anisotropy case,energy decreases with anisotropic parameter g increasing.Susceptibility changes slowly with temperature for g1 and exponentially decay for g1 in the low temperature region and its behaviors approach consistent for different g in the high temperature region.

Journal of Guangxi University•Jan 1, 2011•0 citations

Thermodynamic properties of the two-dimensional quantum Heisenberg ferromagnet and the effects of bond dilution

E. F. Talbot

The energy, the specific heat, and the magnetic susceptibility of the bond-diluted S=(1/2) Heisenberg ferromagnet on a two-dimensional square lattice are studied using a Monte Carlo sampling of the terms contributing to the partition function. The full lattice is studied in detail for sizes up to 20\ifmmode\times\else\texttimes\fi{}20, and evidence is presented to suggest that the low-temperature correlation length behaves as \ensuremath{\xi}\ensuremath{\sim}${e}^{\ensuremath{\nu}\ensuremath{\beta}J}$ with \ensuremath{\nu}\ensuremath{\sim}1.91\ifmmode\pm\else\textpm\fi{}0.06. The temperature dependence of the susceptibility, however, is not quite a simple exponential in \ensuremath{\beta}J. Knowledge of the correlation length and finite-size scaling allows the determination of the prefactor. Our conclusion is that the low-temperature susceptibility behaves as ${T}^{\mathrm{\ensuremath{-}}1}$${e}^{2\ensuremath{\nu}\ensuremath{\beta}J}$. Square lattices of size ranging from 6\ifmmode\times\else\texttimes\fi{}6 to 12\ifmmode\times\else\texttimes\fi{}12 are studied for ${k}_{B}$T/J ranging from 0.1 to 3 and for several values of the bond probability p. Not surprisingly, the most drastic thermodynamic behavior near the critical bond-percolation probability is the sharp increase in the low-temperature susceptibility for p>${p}_{c}$. This corresponds to much larger correlation effects in the percolating lattice than below the percolation threshold. We also compare the thermodynamics at p=${p}_{c}$ to the scaling theories of Lubensky and Stanley et al. which suggest that the correlations at the percolation threshold are predominantly one dimensional.

Physical review. B, Condensed matter•Apr 1, 1986•5 citations

Green’s function study of one-dimensional spin-1/2 Ising-Heisenberg ferromagnets

W. J. Zhang, J. L. Shen, Jian Xu +1 more

The one-dimensional spin-1/2 Ising-Heisenberg ferromagnets are studied by the Green's function method at finite temperatures. The Green's function is decoupled in terms of the correlation functions of 〈${\mathit{S}}_{0}^{\mathit{z}}$${\mathit{S}}_{\mathit{n}}^{\mathit{z}}$〉 and 〈${\mathit{S}}_{0}^{+}$${\mathit{S}}_{\mathit{n}}^{\mathrm{\ensuremath{-}}}$〉 under the constraint of 〈${\mathit{S}}^{\mathit{z}}$〉=0, from which a set of self-consistent equations are derived and solved numerically. The thermodynamic properties such as internal energy, specific heat, and zero-field parallel susceptibility are calculated. Our results predict the correct behavior for the correlation functions and thermodynamic quantities in the entire temperature region, in agreement both with the exact numerical results on the finite chain and the experimental data.

Physical review. B, Condensed matter•Feb 1, 1995•7 citations

The low-temperature magnetic properties of the one-dimensional anisotropic Heisenberg ferromagnet in a magnetic field

Liangjie Fu

The spin-1/2 quantum Heisenberg ferromagnet with anisotropic spin exchange interactions in one dimension is investigated by means of the Green's function method and the spin-wave theory. In this paper we use the Holstein-Primakoff transformation to get the approximately Hamiltonian. We show that the magnetic field and anisotropy can have a rather drastic effect on the correlation functions and susceptibilities in the low-temperature region. Our results show that inverse correlation lengths kx and ky are independent of the temperature in the weak magnetic field. kx is less dependent on the anisotropy than ky. These results are consistent with the theoretical and experimental results.

Journal of Guangzhou University•Jan 1, 2009•0 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
We use the Quantum Monte Carlo Stochastic Series Expansion (SSE) algorithm to investigate the properties of the one-dimensional Heisenberg ferromagnet.We consider the model for two spin values s = 1 2 and s = 1 and study the dependence of the specific heat, magnetisation, magnetic susceptibility and correlation functions on temperature and magnetic field.Obtained results are compared to those derived by analytical methods (Bethe ansatz, Green's function, Random Phase Approximation).In particular, the surprising prediction of a double-peak structure in the specific heat at low temperatures and small magnetic fields can be confirmed.
We investigate the Stochastic Series Expansion (SSE) technique, guided by Sandvik's algorithm of the quantum Monte Carlo simulation for spin-1/2 antiferromagnetic Heisenberg linear chain. The impact of the Monte Carlo steps, lattice sizes, and various temperatures on the variability and fidelity of simulation outcomes are studied. The findings show that the Monte Carlo steps where measurements are made (Msteps), and lattice size are crucial for stabilizing the results. Our results show that for larger lattice sizes, such as L=1024, achieving reliable outcomes necessitates extensive Monte Carlo simulations, specifically. The convergence in the properties is achieved with a high number of Monte Carlo steps starting from 7×10^4. Strong fluctuations are calculated between the results for each magnetic property and the existence of negative and positive in the specific heat at low temperatures.
The Antiferromagnetic Heisenberg Model with spin 1/2 on a square lattice was simulated by using the Stochastic Series Expansion(SSE) Quantum Monte Carlo Method.The temperature dependence of the internal energy,specific heat and uniform susceptibility was observed.The results showed that specific heat has a maximum at kT/J=0.6 and uniform susceptibility saturates at kT/J=1 for the isotropy case.And the size of lattice also affects the temperature dependence of these thermodynamic quantities.For the anisotropy case,energy decreases with anisotropic parameter g increasing.Susceptibility changes slowly with temperature for g1 and exponentially decay for g1 in the low temperature region and its behaviors approach consistent for different g in the high temperature region.
The energy, the specific heat, and the magnetic susceptibility of the bond-diluted S=(1/2) Heisenberg ferromagnet on a two-dimensional square lattice are studied using a Monte Carlo sampling of the terms contributing to the partition function. The full lattice is studied in detail for sizes up to 20\ifmmode\times\else\texttimes\fi{}20, and evidence is presented to suggest that the low-temperature correlation length behaves as \ensuremath{\xi}\ensuremath{\sim}${e}^{\ensuremath{\nu}\ensuremath{\beta}J}$ with \ensuremath{\nu}\ensuremath{\sim}1.91\ifmmode\pm\else\textpm\fi{}0.06. The temperature dependence of the susceptibility, however, is not quite a simple exponential in \ensuremath{\beta}J. Knowledge of the correlation length and finite-size scaling allows the determination of the prefactor. Our conclusion is that the low-temperature susceptibility behaves as ${T}^{\mathrm{\ensuremath{-}}1}$${e}^{2\ensuremath{\nu}\ensuremath{\beta}J}$. Square lattices of size ranging from 6\ifmmode\times\else\texttimes\fi{}6 to 12\ifmmode\times\else\texttimes\fi{}12 are studied for kB{k}_{B}kB​T/J ranging from 0.1 to 3 and for several values of the bond probability p. Not surprisingly, the most drastic thermodynamic behavior near the critical bond-percolation probability is the sharp increase in the low-temperature susceptibility for p>pc{p}_{c}pc​. This corresponds to much larger correlation effects in the percolating lattice than below the percolation threshold. We also compare the thermodynamics at p=pc{p}_{c}pc​ to the scaling theories of Lubensky and Stanley et al. which suggest that the correlations at the percolation threshold are predominantly one dimensional.
The one-dimensional spin-1/2 Ising-Heisenberg ferromagnets are studied by the Green's function method at finite temperatures. The Green's function is decoupled in terms of the correlation functions of 〈S0z{\mathit{S}}_{0}^{\mathit{z}}S0z​Snz{\mathit{S}}_{\mathit{n}}^{\mathit{z}}Snz​〉 and 〈S0+{\mathit{S}}_{0}^{+}S0+​${\mathit{S}}_{\mathit{n}}^{\mathrm{\ensuremath{-}}}$〉 under the constraint of 〈Sz{\mathit{S}}^{\mathit{z}}Sz〉=0, from which a set of self-consistent equations are derived and solved numerically. The thermodynamic properties such as internal energy, specific heat, and zero-field parallel susceptibility are calculated. Our results predict the correct behavior for the correlation functions and thermodynamic quantities in the entire temperature region, in agreement both with the exact numerical results on the finite chain and the experimental data.
The spin-1/2 quantum Heisenberg ferromagnet with anisotropic spin exchange interactions in one dimension is investigated by means of the Green's function method and the spin-wave theory. In this paper we use the Holstein-Primakoff transformation to get the approximately Hamiltonian. We show that the magnetic field and anisotropy can have a rather drastic effect on the correlation functions and susceptibilities in the low-temperature region. Our results show that inverse correlation lengths kx and ky are independent of the temperature in the weak magnetic field. kx is less dependent on the anisotropy than ky. These results are consistent with the theoretical and experimental 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.