Loop algorithm for classical antiferromagnetic Heisenberg models with biquadratic interactions

Hiroshi Shinaoka, Yusuke Tomita, Yukitoshi Motome

Research output: Contribution to journalArticle

2 Citations (Scopus)

Abstract

Monte Carlo simulation using the standard single-spin flip algorithm often fails to sample over the entire configuration space at low temperatures for frustrated spin systems. A typical example is a class of spin-ice type Ising models. In this case, the difficulty can be avoided by introducing a global-flip algorithm, the loop algorithm. Similar difficulty is encountered in O(3) Heisenberg models in the presence of biquadratic interaction. The loop algorithm, however, is not straightforwardly applied to this case, since the system does not have a priori spin-anisotropy axis for constructing the loops. We propose an extension of the loop algorithm to the bilinear-biquadratic models. The efficiency is tested for three different ways to flip spins on a loop in Monte Carlo simulation. We show that the most efficient method depends on the strength of the biquadratic interaction.

Original languageEnglish
Article number012009
JournalJournal of Physics: Conference Series
Volume320
DOIs
Publication statusPublished - 2011
Externally publishedYes

Fingerprint

interactions
Ising model
ice
simulation
anisotropy
configurations

ASJC Scopus subject areas

  • Physics and Astronomy(all)

Cite this

Loop algorithm for classical antiferromagnetic Heisenberg models with biquadratic interactions. / Shinaoka, Hiroshi; Tomita, Yusuke; Motome, Yukitoshi.

In: Journal of Physics: Conference Series, Vol. 320, 012009, 2011.

Research output: Contribution to journalArticle

@article{f1923ddef0f84eb888cc200a6c393bc1,
title = "Loop algorithm for classical antiferromagnetic Heisenberg models with biquadratic interactions",
abstract = "Monte Carlo simulation using the standard single-spin flip algorithm often fails to sample over the entire configuration space at low temperatures for frustrated spin systems. A typical example is a class of spin-ice type Ising models. In this case, the difficulty can be avoided by introducing a global-flip algorithm, the loop algorithm. Similar difficulty is encountered in O(3) Heisenberg models in the presence of biquadratic interaction. The loop algorithm, however, is not straightforwardly applied to this case, since the system does not have a priori spin-anisotropy axis for constructing the loops. We propose an extension of the loop algorithm to the bilinear-biquadratic models. The efficiency is tested for three different ways to flip spins on a loop in Monte Carlo simulation. We show that the most efficient method depends on the strength of the biquadratic interaction.",
author = "Hiroshi Shinaoka and Yusuke Tomita and Yukitoshi Motome",
year = "2011",
doi = "10.1088/1742-6596/320/1/012009",
language = "English",
volume = "320",
journal = "Journal of Physics: Conference Series",
issn = "1742-6588",
publisher = "IOP Publishing Ltd.",

}

TY - JOUR

T1 - Loop algorithm for classical antiferromagnetic Heisenberg models with biquadratic interactions

AU - Shinaoka, Hiroshi

AU - Tomita, Yusuke

AU - Motome, Yukitoshi

PY - 2011

Y1 - 2011

N2 - Monte Carlo simulation using the standard single-spin flip algorithm often fails to sample over the entire configuration space at low temperatures for frustrated spin systems. A typical example is a class of spin-ice type Ising models. In this case, the difficulty can be avoided by introducing a global-flip algorithm, the loop algorithm. Similar difficulty is encountered in O(3) Heisenberg models in the presence of biquadratic interaction. The loop algorithm, however, is not straightforwardly applied to this case, since the system does not have a priori spin-anisotropy axis for constructing the loops. We propose an extension of the loop algorithm to the bilinear-biquadratic models. The efficiency is tested for three different ways to flip spins on a loop in Monte Carlo simulation. We show that the most efficient method depends on the strength of the biquadratic interaction.

AB - Monte Carlo simulation using the standard single-spin flip algorithm often fails to sample over the entire configuration space at low temperatures for frustrated spin systems. A typical example is a class of spin-ice type Ising models. In this case, the difficulty can be avoided by introducing a global-flip algorithm, the loop algorithm. Similar difficulty is encountered in O(3) Heisenberg models in the presence of biquadratic interaction. The loop algorithm, however, is not straightforwardly applied to this case, since the system does not have a priori spin-anisotropy axis for constructing the loops. We propose an extension of the loop algorithm to the bilinear-biquadratic models. The efficiency is tested for three different ways to flip spins on a loop in Monte Carlo simulation. We show that the most efficient method depends on the strength of the biquadratic interaction.

UR - http://www.scopus.com/inward/record.url?scp=81055147521&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=81055147521&partnerID=8YFLogxK

U2 - 10.1088/1742-6596/320/1/012009

DO - 10.1088/1742-6596/320/1/012009

M3 - Article

AN - SCOPUS:81055147521

VL - 320

JO - Journal of Physics: Conference Series

JF - Journal of Physics: Conference Series

SN - 1742-6588

M1 - 012009

ER -