From trivial to topological paramagnets: The case of mathbb{Z}₂ and mathbb{Z}₂³ symmetries in two dimensions

Using quantum Monte Carlo simulations, we map out the phase diagram of Hamiltonians interpolating between trivial and non-trivial bosonic symmetry-protected topological phases, protected by $\mathbb{Z}_2$ and $\mathbb{Z}_2^3$ symmetries, in two dimensions. In all cases, we find that the trivial and the topological phases are separated by an intermediate phase in which the protecting symmetry is spontaneously broken. Depending on the model, we identify a variety of magnetic orders on the triangular lattice, including ferromagnetism, $\sqrt{3}\times\sqrt{3}$ order, and stripe orders (both commensurate and incommensurate). Critical properties are determined through a finite-size scaling analysis. Possible scenarios regarding the nature of the phase transitions are discussed.