Symmetry breaking and competing valence bond states in the star lattice Heisenberg antiferromagnet

TL;DR

This study explores the complex phase diagram of the star lattice Heisenberg antiferromagnet, revealing a first-order transition and a subtle competition between different valence bond crystal states using advanced numerical methods.

Contribution

It uncovers a low critical ratio for phase transition and demonstrates the competition between two valence bond crystal states, supported by both iPEPS simulations and high-order series expansions.

Findings

First-order phase transition at J_d/J_t ≈ 0.18

Close competition between columnar and √3×√3 VBC states

√3×√3 VBC energetically favored at sixth order in perturbation theory

Abstract

We investigate the ground state phase diagram of the spin-$1/2$ antiferromagnetic Heisenberg model on the star lattice using infinite projected entangled pair states (iPEPS) and high-order series expansions. The model includes two distinct couplings: $J_d$ on the dimer bonds and $J_t$ on the trimer bonds. While it is established that the system hosts a valence bond solid (VBS) phase for $J_d \ge J_t$, the ground state phase diagram for $J_d < J_t$ has remained unsettled. Our iPEPS simulations uncover a first-order phase transition at $J_d/J_t \approx 0.18$, significantly lower than previously reported estimates. Beyond this transition, we identify a close competition between two valence bond crystal (VBC) states: a columnar VBC and a $\sqrt{3} \times \sqrt{3}$ VBC, with the latter consistently exhibiting lower energy across all finite bond dimensions. The high-order series expansion supports this by finding that the $\sqrt{3} \times \sqrt{3}$ VBC state indeed becomes energetically favorable, but only at sixth order in perturbation theory, revealing the subtle nature of the competition between candidate states.