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

AIP conference proceedings

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

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

Published Jan 1, 2026• 0 citations
Publisher

Abstract

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.

DOI
10.1063/5.0319339

Related papers

Impact of Equilibration Steps on Quantum Monte Carlo for the Spin-1/2 Heisenberg Model

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

This study delves into the pivotal role of equilibration Monte Carlo steps (𝑖 steps ) in the Stochastic Series Expansion method-a cutting-edge quantum Monte Carlo technique known for its precision in simulating low-dimensional quantum systems.The research focuses on how varying equilibration steps impact the accuracy and stability of ground state energy, sublattice magnetization, specific heat, and susceptibility results, by analyzing the effect of equilibration steps on these physical properties across different temperatures and lattice sizes.Our findings demonstrate that while ground state energy and sublattice magnetization achieve remarkable stability, susceptibility and specific heat present variabilities, especially at low temperatures.Ground state energy stabilizes after approximately 50 equilibration steps.However, excessive equilibration steps lead to unnecessary computational costs without further accuracy improvements.Furthermore, increasing the measurement steps significantly reduces error bars, enhancing reliability, as demonstrated by a 30% reduction in standard deviation with higher counts.Low temperatures exhibit greater variability in specific heat and susceptibility due to complex quantum effects.

Mathematical Modelling and Engineering Problems•Aug 31, 2025•0 citations

QMC simulations of Heisenberg ferromagnet

L. Bogacz, Wolfhard Janke E

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.

Proceedings of XXIIIrd International Symposium on Lattice Field Theory — PoS(LAT2005)•Dec 14, 2005•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

Deterministic Loop Stochastic Series Expansion Algorithm for Quantum Spin Models in Magnetic Fields

Liuyun Dao, Yan-Cheng Wang, Hui Shao

The stochastic series expansion (SSE) algorithm is one of the most powerful quantum Monte Carlo methods and has been extensively applied to the study of quantum many body systems. Its efficiency is particularly enhanced with a deterministic loop update scheme in the study of the S=1/2 quantum spin systems that preserve SU(2) spin rotational symmetry. Once the symmetry is broken, such as by an external field, a directed loop method is typically required, resulting in a significant reduction in efficiency. Inspired by the SSE approach developed for the quantum Ising model, we introduce a deterministic loop SSE method that is particularly suited for antiferromagnetic systems under a staggered magnetic field. This method enables separate investigations of longitudinal and transverse modes in magnetically ordered phases arising from spontaneous symmetry breaking. We benchmark the performance of our algorithm against the standard directed loop approach applied to the antiferromagnetic Heisenberg chain and demonstrate that our method substantially reduces CPU time per Monte Carlo step, thereby can outperform the directed loop algorithm in efficiency.

arXiv:2604.04635•Apr 6, 2026•0 citations

Quantum Monte Carlo simulation for the conductance of one-dimensional quantum spin systems

Kim Louis, Claudius Gros

Recently, the stochastic series expansion (SSE) has been proposed as a powerful Monte Carlo method, which allows simulations at low T for quantum-spin systems. We show that the SSE allows us to compute the magnetic conductance for various one-dimensional spin systems without further approximations. We consider various modifications of the anisotropic Heisenberg chain. We recover the Kane-Fisher scaling for one impurity in a Luttinger liquid and study the influence of noninteracting leads for the conductance of an interacting system.

arXiv:cond-mat/0309376•Physical review. B, Condensed matter•Nov 26, 2003•17 citations

Loop ensembles in Stochastic Series Expansion of Two-Dimensional Heisenberg Antiferromagnets

Vedant Motamarri

The Stochastic Series Expansion (SSE) method along with resummation over the spin or flavor values maps the partition function of a quantum antiferromagnet to a closely-packed loop gas model in one higher dimension. Earlier work by Nahum et al. has shown that certain closely-packed three-dimensional loop gas models exhibit phases dominated by macroscopic loops, wherein the corresponding joint distribution of loop lengths is Poisson-Dirichlet. On grounds of universality, the same is expected of the ensemble of loops obtained in (2+1)-dimensional quantum antiferromagnets, albeit the loops emerge from a different microscopic origin. We sample the SSE loop ensemble for SU(N) antiferromagnets on a square lattice using Monte Carlo and study how the joint distribution varies with the degree of representation N and inverse temperature $β$. We observe that, for low temperatures and small N($\leq$ 4), the distribution indeed shows characteristics of Poisson-Dirichlet behaviour when antiferromagnetic correlations dominate the system.

arXiv:2306.09366•Jun 14, 2023•0 citations
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.
This study delves into the pivotal role of equilibration Monte Carlo steps (𝑖 steps ) in the Stochastic Series Expansion method-a cutting-edge quantum Monte Carlo technique known for its precision in simulating low-dimensional quantum systems.The research focuses on how varying equilibration steps impact the accuracy and stability of ground state energy, sublattice magnetization, specific heat, and susceptibility results, by analyzing the effect of equilibration steps on these physical properties across different temperatures and lattice sizes.Our findings demonstrate that while ground state energy and sublattice magnetization achieve remarkable stability, susceptibility and specific heat present variabilities, especially at low temperatures.Ground state energy stabilizes after approximately 50 equilibration steps.However, excessive equilibration steps lead to unnecessary computational costs without further accuracy improvements.Furthermore, increasing the measurement steps significantly reduces error bars, enhancing reliability, as demonstrated by a 30% reduction in standard deviation with higher counts.Low temperatures exhibit greater variability in specific heat and susceptibility due to complex quantum effects.
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.
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 stochastic series expansion (SSE) algorithm is one of the most powerful quantum Monte Carlo methods and has been extensively applied to the study of quantum many body systems. Its efficiency is particularly enhanced with a deterministic loop update scheme in the study of the S=1/2 quantum spin systems that preserve SU(2) spin rotational symmetry. Once the symmetry is broken, such as by an external field, a directed loop method is typically required, resulting in a significant reduction in efficiency. Inspired by the SSE approach developed for the quantum Ising model, we introduce a deterministic loop SSE method that is particularly suited for antiferromagnetic systems under a staggered magnetic field. This method enables separate investigations of longitudinal and transverse modes in magnetically ordered phases arising from spontaneous symmetry breaking. We benchmark the performance of our algorithm against the standard directed loop approach applied to the antiferromagnetic Heisenberg chain and demonstrate that our method substantially reduces CPU time per Monte Carlo step, thereby can outperform the directed loop algorithm in efficiency.
Recently, the stochastic series expansion (SSE) has been proposed as a powerful Monte Carlo method, which allows simulations at low T for quantum-spin systems. We show that the SSE allows us to compute the magnetic conductance for various one-dimensional spin systems without further approximations. We consider various modifications of the anisotropic Heisenberg chain. We recover the Kane-Fisher scaling for one impurity in a Luttinger liquid and study the influence of noninteracting leads for the conductance of an interacting system.
The Stochastic Series Expansion (SSE) method along with resummation over the spin or flavor values maps the partition function of a quantum antiferromagnet to a closely-packed loop gas model in one higher dimension. Earlier work by Nahum et al. has shown that certain closely-packed three-dimensional loop gas models exhibit phases dominated by macroscopic loops, wherein the corresponding joint distribution of loop lengths is Poisson-Dirichlet. On grounds of universality, the same is expected of the ensemble of loops obtained in (2+1)-dimensional quantum antiferromagnets, albeit the loops emerge from a different microscopic origin. We sample the SSE loop ensemble for SU(N) antiferromagnets on a square lattice using Monte Carlo and study how the joint distribution varies with the degree of representation N and inverse temperature βββ. We observe that, for low temperatures and small N(≤\leq≤ 4), the distribution indeed shows characteristics of Poisson-Dirichlet behaviour when antiferromagnetic correlations dominate the system.