We consider a model of Dirac fermions in 2+1 dimensions with dynamically generated, anticommuting SO(3) Néel and Z_{2} Kekulé mass terms that permits sign-free quantum Monte Carlo simulations. The phase diagram is obtained from finite-size scaling and includes a direct and continuous transition between the Néel and Kekulé phases. The fermions remain gapped across the transition, and our data support an emergent SO(4) symmetry unifying the two order parameters. While the bare symmetries of our model do not allow for spinon-carrying Z_{3} vortices in the Kekulé mass, the emergent SO(4) invariance permits an interpretation of the transition in terms of deconfined quantum criticality. The phase diagram also features a tricritical point at which the Néel, Kekulé, and semimetallic phases meet. The present sign-free approach can be generalized to a variety of other mass terms and thereby provides a new framework to study exotic critical phenomena.
We analyze emergent quantum multicriticality for strongly interacting, massless Dirac fermions in two spatial dimensions (
d=2) within the framework of Gross-Neveu-Yukawa models, by considering the competing order parameters that give rise to fully gapped (insulating or superconducting) ground states. We focus only on those competing orders which can be rotated into each other by generators of an exact or emergent chiral symmetry of massless Dirac fermions, and break
O(S1) and
O(S2) symmetries in the ordered phase. Performing a renormalization-group analysis by using the $\ensuremath{\epsilon}=(3\ensuremath{-}d)$ expansion scheme, we show that all the coupling constants in the critical hyperplane flow toward a new attractive fixed point, supporting an enlarged
O(S1+S2) chiral symmetry. Such a fixed point acts as an exotic quantum multicritical point (MCP), governing the continuous semimetal-insulator as well as insulator-insulator (for example, antiferromagnet to valence bond solid) quantum phase transitions. In comparison with the lower symmetric semimetal-insulator quantum critical points, possessing either
O(S1) or
O(S2) chiral symmetry, the MCP displays enhanced correlation length exponents, and anomalous scaling dimensions for both fermionic and bosonic fields. We discuss the scaling properties of the ratio of bosonic and fermionic masses, and the increased dc resistivity at the MCP. By computing the scaling dimensions of different local fermion bilinears in the particle-hole channel, we establish that most of the four fermion operators or generalized density-density correlation functions display faster power-law decays at the MCP compared to the free fermion and lower symmetric itinerant quantum critical points. Possible generalization of this scenario to higher-dimensional Dirac fermions is also outlined.
We resolve the nature of the quantum phase transition between antiferromagnetic and valence-bond solid ground states in square-lattice spin-1/2 magnets - the putative deconfined quantum-critical point. Using quantum Monte Carlo simulations of J-Q models, in which Heisenberg exchange J competes with multi-spin interactions Qn formed by n singlet projectors, we demonstrate first-order discontinuities increasing with n. For n=2 and n=3, the first-order signatures are very weak and we extract universal scaling dimensions of a proximate multicritical point. Their values coincide closely with known bounds for an SO(5) symmetric conformal field theory (CFT) containing a relevant SO(5) singlet operator that destabilizes the generic deconfined critical point. We characterize the emergent SO(5) symmetry by its scaling dimensions of the leading irrelevant perturbations. Using the J-Q2-Q6 model, we tune the transition from weak to stronger first-order and find scaling of the coexisting order parameters with a large exponent beta ~ 0.85, in agreement with the relevant scaling dimensions of the multicritical point. The exponent that had previously been identified with the correlation length does not have a corresponding CFT level, and we explain it as pseudo-critical length scale of bubbles of one phase inside the other that are compressed in size by an SO(5)-violating irrelevant repulsive interaction between the two order parameters. The bubble scale reflects a hierarchy of excitation gaps generated by the SO(5) perturbation in the neighborhood of the first-order transition, akin to a Higgs mechanism. We argue that the multicritical point is the tip of a gapless spin liquid phase in frustrated Heisenberg models into which the J-Q models can be deformed, but the spin liquid, and perhaps the multicritical point itself, exist outside the regime amenable to sign-free simulations.
We study the quantum multicritical point in a (2+1)-dimensional Dirac system between the semimetallic phase and two ordered phases that are characterized by anticommuting mass terms with
O(N1) and
O(N2) symmetries, respectively. Using $\ensuremath{\epsilon}$ expansion around the upper critical space-time dimension of four, we demonstrate the existence of a stable renormalization-group fixed point, enabling a direct and continuous transition between the two ordered phases directly at the multicritical point. This point is found to be characterized by an emergent
O(N1+N2) symmetry for arbitrary values of
N1 and
N2 and fermion flavor numbers
Nf as long as the corresponding representation of the Clifford algebra exists. Small
O(N)-breaking perturbations near the chiral
O(N) fixed point are therefore irrelevant. This result can be traced back to the presence of gapless Dirac degrees of freedom at criticality, and it is in clear contrast to the purely bosonic
O(N) fixed point, which is stable only when
N<3. As a by-product, we obtain predictions for the critical behavior of the chiral
O(N) universality classes for arbitrary
N and fermion flavor number
Nf. Implications for critical Weyl and Dirac systems in 3+1 dimensions are also briefly discussed.
Motivated by the physics of spin-orbital liquids, we study a model of interacting Dirac fermions on a bilayer honeycomb lattice at half filling, featuring an explicit global SO(3)
×U(1) symmetry. Using large-scale auxiliary- field quantum Monte Carlo (QMC) simulations, we locate two zero-temperature phase transitions as function of increasing interaction strength. First, we observe a continuous transition from the weakly-interacting semimetal to a different semimetallic phase in which the SO(3) symmetry is spontaneously broken and where two out of three Dirac cones acquire a mass gap. The associated quantum critical point can be understood in terms of a Gross-Neveu-SO(3) theory. Second, we subsequently observe a transition towards an insulating phase in which the SO(3) symmetry is restored and the U(1) symmetry is spontaneously broken. While strongly first order at the mean-field level, the QMC data is consistent with a direct and continuous transition. It is thus a candidate for a new type of deconfined quantum critical point that features gapless fermionic degrees of freedom.
Deconfined quantum criticality with emergent SO(5) symmetry in correlated systems remains elusive. Here, by performing numerically-exact state-of-the-art quantum Monte Carlo (QMC) simulations, we show convincing evidences of deconfined quantum critical points (DQCP) between antiferromagnetic and valence-bond-solid phases in the extended Hubbard model of fermions on the honeycomb lattice with large system sizes. We further demonstrate evidences of the SO(5) symmetry at the DQCP. It is important to note that the critical exponents obtained by finite-size scaling at the DQCP here are consistent with the rigourous conformal bounds. Consequently, we established a promising arena of DQCP with emergent SO(5) symmetry in interacting systems of fermions. Its possible experimental relevances in correlated systems of Dirac fermions will be discussed briefly.