We investigate nematic quantum phase transitions in two different Dirac fermion models. The models feature twofold and fourfold, respectively, lattice rotational symmetries that are spontaneously broken in the ordered phase. Using negative-sign-free quantum Monte Carlo simulations and an ε-expansion renormalization group analysis, we show that both models exhibit continuous phase transitions. In contrast to generic Gross-Neveu dynamical mass generation, the quantum critical regime is characterized by large velocity anisotropies, with fixed-point values being approached very slowly. Both experimental and numerical investigations will not be representative of the infrared fixed point, but of a quasiuniversal regime where the drift of the exponents tracks the velocity anisotropy.
Abstract A central concept in the theory of phase transitions beyond the Landau–Ginzburg–Wilson paradigm is fractionalization: the formation of new quasiparticles that interact via emergent gauge fields. This concept has been extensively explored in the context of continuous quantum phase transitions between distinct orders that break different symmetries. We propose a mechanism for continuous order-to-order quantum phase transitions that operates independently of fractionalization. This mechanism is based on the collision and annihilation of two renormalization group fixed points: a quantum critical fixed point and an infrared stable fixed point. The annihilation of these fixed points rearranges the flow topology, eliminating the disordered phase associated with the infrared stable fixed point and promoting a second critical fixed point, unaffected by the collision, to a quantum critical point between distinct orders. We argue that this mechanism is relevant to a broad spectrum of physical systems. In particular, it can manifest in Luttinger fermion systems in three spatial dimensions, leading to a continuous quantum phase transition between an antiferromagnetic Weyl semimetal state, which breaks time-reversal symmetry, and a nematic topological insulator, characterized by broken lattice rotational symmetry. This continuous antiferromagnetic-Weyl-to-nematic-insulator transition might be observed in rare-earth pyrochlore iridates R 2 Ir2O7. Other possible realizations include kagome quantum magnets, quantum impurity models, and quantum chromodynamics with supplemental four-fermion interactions.
The canonical Gross-Neveu model for N two-component Dirac fermions in 2+1 dimensions suffers a continuous phase transition at a critical interaction gc1∼1/N at large N , at which its continuous symmetry SO(2N) is preserved and a discrete (Ising) symmetry becomes spontaneously broken. A recent mean-field calculation, however, points to an additional transition at a different critical gc2∼−Ngc1 , at which SO(2N)→SO(N)×SO(N) . To study the latter phase transition we rewrite the Gross-Neveu interaction g(ψ¯ψ)2 in terms of three different quartic terms for the single ( L=1 ) 4N -component real (Majorana) fermion, and then extend the theory to L>1 . This allows us to track the evolution of the fixed points of the renormalization group transformation starting from L≫1 , where one can discern three distinct critical points which correspond to continuous phase transitions into (1) SO(2N) -singlet mass-order-parameter, (2) SO(2N) -symmetric-tensor mass-order-parameters, and (3) SO(2N) -adjoint nematic-order-parameters, down to L=1 value that is relevant to the standard Gross-Neveu model. Below the critical value of Lc(N)≈0.35N for N≫1 only the Gross-Neveu critical point (1) still implies a diverging susceptibility for its corresponding ( SO(2N) -singlet) order parameter, whereas the two new critical points that existed at large L ultimately become equivalent to the Gaussian fixed point at L=1 . We interpret this metamorphosis of the SO(2N) -symmetric-tensor fixed point from critical to spurious as an indication that the transition at gc2 in the original Gross-Neveu model is turned first-order by fluctuations. Published by the American Physical Society 2024
Quantum critical phenomena may be qualitatively different when massless Dirac fermions are present at criticality. Using our recently-discovered fermion-sign-free Majorana quantum Monte Carlo method introduced by us in (Li et al 2015 Phys. Rev. B 91 241117 ), we investigate the quantum critical phenomena of spinless Dirac fermions at their charge-density-wave phase transitions on the honeycomb lattice having sites with largest L = 24. By finite-size scaling, we accurately obtain critical exponents of this so-called Gross–Neveu chiral-Ising universality class of two (two-component) Dirac fermions in 2+1D: , , and , which are qualitatively different from the mean-field results but are reasonably close to the ones obtained from renormalization group calculations.
We present quantum Monte Carlo simulations for the chiral Heisenberg Gross-Neveu-Yukawa quantum phase transition of relativistic fermions with N=4 Dirac spinor components subject to a repulsive, local four fermion interaction in (2+1)D. Here we employ a two-dimensional lattice Hamiltonian with a single, spin-degenerate Dirac cone, which exactly reproduces a linear energy-momentum relation for all finite size lattice momenta in the absence of interactions. This allows us to significantly reduce finite size corrections compared to the widely studied honeycomb and π-flux lattices. A Hubbard term dynamically generates a mass beyond a critical coupling of U_{c}=6.76(1) as the system acquires antiferromagnetic order and SU(2) spin rotational symmetry is spontaneously broken. At the quantum phase transition, we extract a self-consistent set of critical exponents ν=0.98(1), η_{ϕ}=0.53(1), η_{ψ}=0.18(1), and β=0.75(1). We provide evidence for the continuous degradation of the quasiparticle weight of the fermionic excitations as the critical point is approached from the semimetallic phase. Finally, we study the effective "speed of light" of the low-energy relativistic description, which depends on the interaction U, but is expected to be regular across the quantum phase transition. We illustrate that the strongly coupled bosonic and fermionic excitations share a common velocity at the critical point.
Evidence for relativistic quantum criticality of antiferromagnetism and superconductivity in two-dimensional Dirac fermion systems has been found in large-scale quantum Monte Carlo simulations. However, the corresponding (
2+1)-dimensional Gross--Neveu--Yukawa field theory with
Nf=2 four-component Dirac fermions coupled to two triplets of order parameters does not exhibit a renormalization group fixed point at one-loop order. Instead, the theory only features a critical point for a large or very small fractional number of fermion flavors
Nf, which disappears for a broad range of flavor numbers around the physical case,
Nf=2, due to fixed-point annihilation. This raises the question on how to explain the observed scaling collapse in the quantum Monte Carlo data. Here, we extend previous renormalization-group analyses by studying a generalized model at two-loop order in $4\ensuremath{-}\ensuremath{\epsilon}$ spacetime dimensions. We determine the $\ensuremath{\epsilon}$ correction to the upper and lower critical flavor numbers for the fixed-point annihilation and find that they both go towards the physical case
Nf=2. However, this only happens very slowly, such that an extrapolation to $\ensuremath{\epsilon}=1$ still suggests the absence of criticality in
2+1 dimensions. Thereby, we consolidate the finding that the continuum field theory does not feature a stable renormalization-group fixed point and no true quantum criticality would be expected for the considered system. We briefly discuss a possible reconciliation in terms of a complex conformal field theory. Further, we also explore the fixed-point structure in an enlarged theory space and identify a candidate stable fixed-point solution.
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.