We use the half-filled zeroth Landau level in graphene as a regularization scheme to study the physics of the SO(5) nonlinear sigma model subject to a Wess-Zumino-Witten topological term in 2+1 dimensions. As shown by Ippoliti et al. [Phys. Rev. B 98, 235108 (2019)PRBMDO2469-995010.1103/PhysRevB.98.235108], this approach allows for negative sign free auxiliary field quantum Monte Carlo simulations. The model has a single free parameter U_{0} that monitors the stiffness. Within the parameter range accessible to negative sign free simulations, we observe an ordered phase in the large U_{0} or stiff limit. Remarkably, upon reducing U_{0} the magnetization drops substantially, and the correlation length exceeds our biggest system sizes, accommodating 100 flux quanta. The implications of our results for deconfined quantum phase transitions between valence bond solids and antiferromagnets are discussed.
We use the half-filled zeroth Landau level in graphene as a regularization scheme to study the physics of the SO(5) nonlinear sigma model subject to a Wess-Zumino-Witten topological term in 2+1 dimensions. As shown by Ippoliti et al. [Phys. Rev. B 98, 235108 (2019)PRBMDO2469-995010.1103/PhysRevB.98.235108], this approach allows for negative sign free auxiliary field quantum Monte Carlo simulations. The model has a single free parameter U_{0} that monitors the stiffness. Within the parameter range accessible to negative sign free simulations, we observe an ordered phase in the large U_{0} or stiff limit. Remarkably, upon reducing U_{0} the magnetization drops substantially, and the correlation length exceeds our biggest system sizes, accommodating 100 flux quanta. The implications of our results for deconfined quantum phase transitions between valence bond solids and antiferromagnets are discussed.
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(βNqNϕ2), yielding a speedup of a factor of Nϕ (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ϕ=140 on the torus and Nϕ=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 U(1)×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.
We explore a method for regulating 2+1D quantum critical points in which the ultraviolet cutoff is provided by the finite density of states of particles in a magnetic field rather than by a lattice. Such Landau-level quantization allows for numerical computations on arbitrary manifolds, like spheres, without introducing lattice defects. In particular, when half-filling a Landau level with N=4 electron flavors, with appropriate interaction anisotropies in flavor space, we obtain a fully continuum regularization of the O(5) nonlinear sigma model with a topological term, which has been conjectured to flow to a deconfined quantum critical point. We demonstrate that this model can be solved by both infinite density-matrix renormalization group (DMRG) calculations and sign-free determinantal quantum Monte Carlo. DMRG calculations estimate the scaling dimension of the O(5) vector operator to be in the range ${\mathrm{\ensuremath{\Delta}}}_{V}\ensuremath{\sim}0.55--0.7$, depending on the stiffness of the nonlinear sigma model. Future Monte Carlo simulations will be required to determine whether this dependence is a finite-size effect or further evidence for a weak first-order transition.
Critical phenomena beyond the Landau-Ginzburg-Wilson paradigm have been long sought after. Among many candidate scenarios, the deconfined quantum critical point (DQCP) constitutes the most fascinating one, and its lattice model realization has been debated over the past two decades. Following pioneering works with the fuzzy sphere methods [Zhu et al., Phys. Rev. X 13, 021009 (2023); Hu et al., Phys. Rev. Lett. 131, 031601 (2023); Han et al., arXiv:2312.04047; Zhou et al., Phys. Rev. X 14, 021044 (2024)], we apply the spherical Landau level regularization to study the effective (2+1)D SO(5) nonlinear sigma model with a topological term and the potential DQCP therein. Utilizing the state-of-the-art density matrix renormalization group method with explicit ${\text{SU(2)}}_{\text{spin}}\ifmmode\times\else\texttimes\fi{}{\text{U(1)}}_{\text{charge}}\ifmmode\times\else\texttimes\fi{}{\text{U(1)}}_{\text{angular-momentum}}$ symmetry as well as exact diagonalization simulations, we provide a comprehensive phase diagram for the model with a SO(5) continuous transition line---extension of the previous identified SO(5) multicritical point [Chen et al., Phys. Rev. Lett. 132, 246503 (2024)]---while tuning interaction length. The state-operator correspondence with the conformal tower structure is used to identify the emergent conformal symmetry with the best scaling dimension of relevant primary fields, and they match well with the critical exponents obtained from the crossing point analysis of the correlation ratio. Our results thus further support the rich structure of the phase diagram of the SO(5) model.
Novel critical phenomena beyond the Landau-Ginzburg-Wilson paradigm have been long sought after. Among many candidate scenarios, the deconfined quantum critical point (DQCP) constitutes the most fascinating one, and its lattice model realization has been debated over the past two decades. Here we apply the spherical Landau level regularization upon the exact (2+1)D SO(5) nonlinear sigma model with a topological term to study the potential DQCP therein. We perform a density matrix renormalization group (DMRG) simulation with SU(2)_{spin}×U(1)_{charge}×U(1)_{angular-momentum} symmetries explicitly implemented. Using crossing point analysis for the critical properties of the DMRG data, accompanied by quantum Monte Carlo simulations, we accurately obtain the comprehensive phase diagram of the model and find various novel quantum phases, including Néel, ferromagnet (FM), valence bond solid (VBS), valley polarized (VP) states and a gapless quantum disordered phase occupying an extended area of the phase diagram. The VBS-disorder and Néel-disorder transitions are continuous with non-Wilson-Fisher exponents. Our results show the VBS and Néel states are separated by either a weakly first-order transition or the disordered region with a multicritical point in between, thus opening up more interesting questions on the two-decade long debate on the nature of the DQCP.
It has been proposed that the deconfined criticality in (2+1)d---the quantum phase transition between a N\'eel antiferromagnet and a valence-bond solid (VBS)---may actually be pseudocritical, in the sense that it is a weakly first-order transition with a generically long correlation length. The underlying field theory of the transition would be a slightly complex (nonunitary) fixed point as a result of fixed points annihilation. This proposal was motivated by existing numerical results from large scale Monte Carlo simulations as well as a conformal bootstrap. However, an actual theory of such a complex fixed point, incorporating key features of the transition such as the emergent SO(5) symmetry, is so far absent. Here we propose a Wess-Zumino-Witten (WZW) nonlinear sigma model with level k=1, defined in $2+\ensuremath{\epsilon}$ dimensions, with target space ${S}^{3+\ensuremath{\epsilon}}$ and global symmetry $\text{SO}(4+\ensuremath{\epsilon})$. This gives a formal interpolation between the deconfined criticality at d=3 and the SU(2)1 WZW theory at d=2 describing the spin-1/2 Heisenberg chain. The theory can be formally controlled, at least to leading order, in terms of the inverse of the WZW level 1/k. We show that at leading order there is a fixed point annihilation at ${d}^{*}\ensuremath{\approx}2.77$, with complex fixed points above this dimension including the physical d=3 case. The pseudocritical properties such as correlation length, scaling dimensions, and the drifts of scaling dimensions as the system size increases, calculated crudely to leading order, are qualitatively consistent with existing numerics.
Dirac fermions in (2 + 1) dimensions with dynamically generated anticommuting SO(3) antiferromagnetic and Z2 Kekul\'e valence-bond solid (KVBS) masses map onto a field theory with a topological $\ensuremath{\theta}$ term. This term provides a mechanism for continuous phase transitions between different symmetry-broken states: topological defects of one phase carry the charge of the other and proliferate at the transition. The $\ensuremath{\theta}$ term implies that a domain wall of the Z2 KVBS order parameter harbors a spin-1/2 Heisenberg chain, as described by a (1 + 1)-dimensional SO(3) nonlinear sigma model with $\ensuremath{\theta}$ term at $\ensuremath{\theta}=\ensuremath{\pi}$. Using pinning fields to stabilize the domain wall, we show that our auxiliary-field quantum Monte Carlo simulations indeed support the emergence of a spin-1/2 chain at the Z2 topological defect. Surprisingly, the consequences of the topological term are seen far from the critical point such that the physics of apparently unrelated model systems are naturally understood by invoking them. This concept can be generalized to higher dimensions where (2 + 1)-dimensional SO(4) or SO(5) theories with topological terms are realized at a domain wall.