We introduce a model of Dirac fermions in
2+1 dimensions with a semimetallic, a quantum spin-Hall insulating (QSHI), and an
s-wave superconducting (SSC) phase. The phase diagram features a multicritical point at which all three phases meet as well as a QSHI-SSC deconfined critical point. The QSHI and SSC orders correspond to mutually anticommuting mass terms of the Dirac Hamiltonian. Based on this algebraic property, SO(5) symmetric field theories have been put forward to describe both types of critical points. Using quantum Monte Carlo simulations, we directly study the operator that rotates between QSHI and SSC states. The results suggest that it commutes with the low-energy effective Hamiltonian at criticality but has a gap in the ordered phases. This implies an emergent SO(5) symmetry at both the multicritical and the deconfined critical points.
Quantum criticality is closely related to the existence of two phases with unrelated symmetry breaking. In this paper, we study Néel and Kekulé valence bond state (VBS) quantum criticality in Dirac semimetals with four-fermion interactions. Our results show that all possible dynamical masses yield the same critical coupling, which exhibits the phenomenon that all possible phases meet at a multicritical point (e.g., a tricritical point for the Néel, Kekulé-VBS and semimetallic phases). In terms of the well-established Wess–Zumino–Witten field theory, we investigate the typical criticality for the transition between Néel and Kekulé-VBS phases, and the compatible Néel–Kekulé-VBS mass matrices imply the existence of a non-Landau transition between the Néel and Kekulé-VBS phases. We show the existence of mutual duality in the defect-driven Néel–Kekulé-VBS transition near the non-Landau critical point and find that this mutual duality results from the presence of a mutual Chern–Simons term. We also study the mutual duality based on dual topological excitations.
We consider interacting (2+1)-dimensional Dirac fermions with competing symmetry-breaking electronic instabilities, as described by relativistic quantum field theories of the Gross-Neveu-Yukawa flavor with anti-commuting mass terms. We demonstrate, using a combination of nonperturbative field-theoretical analysis and an adapted quantum Monte Carlo approach, that such systems exhibit a strong-coupling quantum multicritical fixed point with an emerging enhanced symmetry. Moreover, an extended phase coexistence regime expands out from this high-symmetry point. Our results disagree with recent results on the presence of a deconfined quantum criticality in (2+1)-dimensional Dirac fermions for the particular case of O(3) Neel and Z(2) Kekule symmetry-breaking instabilities on the graphene lattice. The robustness of these phenomena with respect to the microscopic symmetries furthermore demonstrates their relevance for a wide range of Dirac materials of current interest, from both theory and ongoing experiments.
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.
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 toward 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 are 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.
According to the Landau criterion, a phase transition should be first order when cubic terms of order parameters are allowed in its effective Ginzburg-Landau free energy. Recently, it was shown by renormalization-group (RG) analysis that continuous transition can happen at putatively first-order
Z3 transitions in two-dimensional (2D) Dirac semimetals and such non-Landau phase transitions were dubbed ``fermion-induced quantum critical points'' (FIQCPs) [Z.-X. Li, Y.-F. Jiang, S.-K. Jian, and H. Yao, Nat. Commun. 8, 314 (2017)]. The RG analysis, controlled by the
1/N expansion with
N the number of flavors of four-component Dirac fermions, shows that the FIQCP occurs for $N\ensuremath{\ge}{N}_{c}$. Previous quantum Monte Carlo simulations of a microscopic model of
SU(N) fermions on the honeycomb lattice showed that the FIQCP occurs at the transition between Dirac semimetals and Kekule valence bond solids for $N\ensuremath{\ge}2$. However, the precise value of the lower bound
Nc has not been established. Especially, the case of
N=1 has not been explored by studying microscopic models so far. Here, by introducing a generalized
SU(N) fermion model with
N=1 (namely, spinless fermions on the honeycomb lattice), we perform large-scale sign-problem-free Majorana quantum Monte Carlo simulations and find convincing evidence of the FIQCP for
N=1. Consequently, our results suggest that the FIQCP can occur in 2D Dirac semimetals for all positive integers $N\ensuremath{\ge}1$.