We study the spin-1/2 Ising-XXZ model on a decorated honeycomb lattice composed of five spins per unit cell, one Ising spin, and four Heisenberg spins. This model involving the Heisenberg exchange interaction is one of the few models that can be exactly solvable through the generalized star-triangle transformation. The significance of this model is its close relationship to the fully decorated quantum Heisenberg honeycomb lattice since 4/5 of the particles are Heisenberg spins. We investigate the phase diagram at zero temperature and identify a relevant quantum spin frustrated phase resulting from the contribution of quantum Heisenberg exchange interaction. We obtain an exact residual entropy for the quantum spin frustrated phase, which coincides with the residual entropy of the antiferromagnetic spin-1/2 Ising model on a triangular lattice. We also thoroughly explore its thermodynamic properties, focusing mainly on the frustrated region such as entropy, specific heat, spontaneous magnetization, and critical temperature under several conditions.
The zero-temperature quantum phase diagram of the spin-
21 J1--
J2--
J1⊥ model on an
AA-stacked bilayer honeycomb lattice is investigated using the coupled cluster method (CCM). The model comprises two monolayers in each of which the spins, residing on honeycomb-lattice sites, interact via both nearest-neighbor (NN) and frustrating next-nearest-neighbor isotropic antiferromagnetic (AFM) Heisenberg exchange iteractions, with respective strengths
J1>0 and
J2≡κJ1>0. The two layers are coupled via a comparable Heisenberg exchange interaction between NN interlayer pairs, with a strength
J1⊥≡δJ1. The complete phase boundaries of two quasiclassical collinear AFM phases, namely the Neel and Neel-II phases, are calculated in the
κδ half-plane with
κ>0. Whereas on each monolayer in the Neel state all NN pairs of spins are antiparallel, in the Neel-II state NN pairs of spins on zigzag chains along one of the three equivalent honeycomb-lattice directions are antiparallel, while NN interchain spins are parallel. We calculate directly in the thermodynamic (infinite-lattice) limit both the magnetic order parameter
M and the excitation energy
Δ from the
sTz=0 ground state to the lowest-lying
∣sTz∣=1 excited state (where
sTz is the total
z component of spin for the system as a whole, and where the collinear ordering lies along the
z direction) for both quasiclassical states used (separately) as the CCM model state, on top of which the multispin quantum correlations are then calculated to high orders (
n≤10) in a systematic series of approximations involving
n-spin clusters. The sole approximation made is then to extrapolate the sequences of
nth-order results for
M and
Δ to the exact limit,
n→∞.