We present new numerical tools to analyze symmetry-broken phases in the context of $SU(2)$-symmetric translation-invariant matrix product states (MPS) and density-matrix renormalization-group (DMRG) methods for infinite cylinders, and determine the phase diagram of the geometrically-frustrated triangular Heisenberg model with nearest and next-nearest neighbor (NN and NNN) interactions. The appearance of Nambu-Goldstone modes in the excitation spectrum is characterized by "tower of states" levels in the momentum-resolved entanglement spectrum. Symmetry-breaking phase transitions are detected by a combination of the correlation lengths and second and fourth cumulants of the magnetic order parameters (which we call the Binder ratio), even though symmetry implies that the order parameter itself is strictly zero. Using this approach, we have identified $120^{\circ}$ order, a columnar order, and an algebraic spin liquid (specific to width-6 systems), alongside the previously studied topological spin liquid phase. For the latter, we also demonstrate robustness against chiral perturbations.
Read full abstract