Monte Carlo study of degenerate ground states and residual entropy in a frustrated honeycomb lattice Ising model

Shawn Andrews, Hans De Sterck, Stephen Inglis, Roger G. Melko

Research output: Contribution to journalArticleResearchpeer-review

9 Citations (Scopus)


We study a classical fully frustrated honeycomb lattice Ising model using Markov-chain Monte Carlo methods and exact calculations. The Hamiltonian realizes a degenerate ground-state manifold of equal-energy states, where each hexagonal plaquette of the lattice has one and only one unsatisfied bond, with an extensive residual entropy that grows as the number of spins N. Traditional single-spin-flip Monte Carlo methods fail to sample all possible spin configurations in this ground state efficiently, due to their separation by large energy barriers. We develop a nonlocal "chain-flip" algorithm that solves this problem, and demonstrate its effectiveness on the Ising Hamiltonian with and without perturbative interactions. The two perturbations considered are a slightly weakened bond and an external magnetic field h. For some cases, the chain-flip move is necessary for the simulation to find an ordered ground state. In the case of the magnetic field, two magnetized ground states with nonextensive entropy are found, and two special values of h exist where the residual entropy again becomes extensive, scaling proportionally to Nln, where is the golden ratio.

Original languageEnglish
Article number041127
JournalPhysical Review E
Issue number4
Publication statusPublished - 1 Apr 2009
Externally publishedYes

Cite this