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

Abstract

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 frustrated 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 non-local "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 non-extensive entropy are found, and two special values of h exist where the residual entropy again becomes extensive, scaling proportional to N ln phi, where phi is the golden ratio.

arXiv comment12+ pages, 11 figures: to appear in PRE

References in corpus (17)

• Spin Ice State in Frustrated Magnetic Pyrochlore Materials
• An RVB phase in the triangular lattice quantum dimer model
• Dipolar Interactions and Origin of Spin Ice in Ising Pyrochlore Magnets
• Ising models of quantum frustration
• Supersolid hardcore bosons on the triangular lattice
• Supersolid Order from Disorder: Hard-Core Bosons on the Triangular Lattice
• Long Range Order at Low Temperature in Dipolar Spin Ice
• Quantum dimer model on the kagome lattice: solvable dimer liquid and Ising gauge theory
• Persistent supersolid phase of hard-core bosons on the triangular lattice
• Why spin ice obeys the ice rules
• Topologically protected quantum bits from Josephson junction arrays
• Thermodynamic Properties of the Dipolar Spin Ice Model
• Classical dimers on the triangular lattice
• Superfluids and Supersolids on Frustrated 2D Lattices
• Monte Carlo simulation of ice models
• Magnetization process of spin ice in a [111] magnetic field
• Magnetic properties of a helical spin chain with alternating isotropic and anisotropic spins: magnetization plateaus and finite entropy

Cited by in corpus (7)

• Magnetic ordered structure dependence of magnetic refrigeration efficiency
• Thermodynamic description of the Ising antiferromagnet on triangular lattice with selective dilution by modified Pair Approximation method
• A generalized Monte Carlo loop algorithm for frustrated Ising models
• Spin disorder in an Ising honeycomb chain cobaltate
• Residual Entropy of Ice: A Study Based on Transfer Matrices
• Mott Insulators in a Fully-Frustrated Bose Hubbard Model on the Honeycomb Lattice
• Boltzmann Sampling of Frustrated J1 - J2 Ising Models with Programmable Quantum Annealers