Rotational symmetry breaking in Heisenberg model on triangular lattice

We study a finite-temperature phase transition in the two-dimensional classical Heisenberg model on a triangular lattice with a ferromagnetic nearest-neighbor interaction $J_1$ and an antiferromagnetic third-nearest-neighbor interaction $J_3$ using Monte Carlo simulation. Apart from a trivial degeneracy corresponding to O(3) spin rotations, the ground state for $J_3\ne 0$ has a threefold degeneracy corresponding to 120 degree lattice rotations. We find that this model exhibits a phase transition with breaking of the three-fold symmetry when $J3=J1/3$ and that the transition is of the first order.