Numerical evidence of a critical point in the (2+1)D SO(5) nonlinear sigma model with Wess-Zumino-Witten term

Abstract

We develop an optimized continuous-field quantum Monte Carlo (QMC) algorithm to investigate the projected SO(5) nonlinear sigma model with a Wess-Zumino-Witten term, which describes half-filled Dirac fermions in 2+1 space-time dimensions akin to graphene and Yukawa coupled to a quintuplet of compatible mass terms. Our algorithm reduces the computational complexity to $O(\beta N_{\mathbf{q}}{N}_{\varphi }^2)$, yielding a speedup of a factor of $N_{\varphi }$ (the number of magnetic fluxes, i.e., system size) relative to prior works [1-4]. This advance enables us to simulate system sizes up to $N_{\varphi }=140$ on the torus and $N_{\varphi }=59$ on the sphere, far exceeding the maximum sizes previously accessed, and to map out the universal phase diagram of the model on both geometries. Most notably, we identify and characterize a critical point that separates an SO(5)-broken ordered phase at small coupling from an SO(5)-symmetric disordered phase at large coupling. The critical point becomes multicritical upon the inclusion of terms that break the SO(5) symmetry down to $\mathrm{U}(1)\times \mathrm{SU}(2)$, relevant for the deconfined phase transition between Néel antiferromagnetic and valence-bond-solid orders in quantum magnets. Our finding of a multicritical point in the phase diagram of the SO(5) nonlinear sigma model with Wess-Zumino-Witten term resolves the long-standing open question of its global structure, and our QMC algorithm opens a new avenue for systematic studies of projected Hamiltonians, ranging from correlated flat bands to fractional quantum (anomalous) Hall systems.

16+5 pages, 5+4 figures