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.
Spinless fermions on the honeycomb lattice with repulsive nearest-neighbor interactions are known to harbour a quantum critical point at half-filling, with critical behavior in the Gross-Neveu (chiral Ising) universality class. The critical interaction strength separates a weak-coupling semimetallic regime from a commensurate charge-density-wave phase. The phase diagram of this basic model of correlated fermions on the honeycomb lattice beyond half-filling is, however, less well established. Here, we perform an analysis of its many-body instabilities using the functional renormalization group method with a basic Fermi surface patching scheme, which allows us to treat instabilities in competing channels on equal footing also away from half-filling. Between half-filling and the Van Hove filling, the free Fermi surface is holelike and we again find a charge-density wave instability to be dominant at large interactions. Moreover, its characteristics are those of the half-filled case. Directly at the Van Hove filling, the nesting property of the free Fermi surface stabilizes a dimerized bond-order phase. At lower filling, the free Fermi surface becomes electronlike and a superconducting instability with
f-wave symmetry is found to emerge from the interplay of intra-unit-cell repulsion and collective fluctuations in the proximity to the charge-density wave instability. We estimate the extent of the various phases and extract the corresponding order parameters from the effective low-energy Hamiltonians.
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$.
We study the phase diagram of interacting spinless fermions on the honeycomb bilayer at charge neutrality using large-scale quantum Monte Carlo simulations. In the noninteracting limit, the low-energy spectrum features quadratically dispersing bands that touch at the corners of the hexagonal Brillouin zone. Weak to intermediate interactions induce a splitting of each of the quadratic band touching points into four Dirac points, located along high-symmetry directions of the reciprocal lattice. Strong interactions lead to the formation of a layer-polarized charge density wave, which spontaneously breaks the
Z2 layer inversion symmetry and opens an insulating gap in the spectrum. We show that the semimetal-to-insulator quantum phase transition as a function of interaction is continuous and characterized by emergent relativistic symmetry. Our results for the values of the correlation-length exponent
ν, the order-parameter anomalous dimension
ηϕ, and the fermion anomalous dimension
ηψ agree with those of the theoretically predicted 2+1D Gross-Neveu-Ising universality class with eight two-component Dirac fermions within less than 5\%\ deviation. We also determine the crossover scale as a function of interaction strength between the nonrelativistic semimetal state at high temperatures, characterized by dynamical critical exponent
z=2, and the Dirac semimetal state at intermediates temperatures, characterized by
z=1. Further reducing the temperature below the crossover scale at a fixed value of the interaction strength above the quantum critical point results in a classical ordering transition in the 2D Ising universality class.
This Thesis explores the formation of novel phases in Dirac systems that are driven by strong many-body interactions and collective quantum fluctuations. The central focus is the half-filled honeycomb lattice, which is the prototypical Dirac semimetal with low-energy excitations described by Dirac fermions. Of particular interest are exotic phases of matter that are not favoured at mean-field but are stabilised by fluctuations. Self-consistent path integral techniques are used to demonstrate that extended repulsive interactions favour charge order that breaks translational and other lattice symmetries. This charge order outcompetes interaction induced topological insulating phases, which are dominant at the mean-field level. The interacting quantum critical point was found to be near the critical point of a topological Lifshitz transition. Interacting Dirac fermions exhibit perhaps the simplest example of fermionic quantum criticality. Gross-Neveu-Yukawa field theories describe the universality class of semimetal-insulator transitions, which are driven by on-site and nearest neighbour repulsive interactions on the honeycomb lattice. The Dirac fermions couple to the dynamical order parameter fields and play a crucial role in determining the universal behaviour. Naturally, this variety of criticality is outside of the Ginzburg-Landau-Wilson paradigm. This idea is extended to include the translational symmetry breaking charge order. Here it is found that the breaking of crystal symmetries corresponds to Lorentz violating effective field theories. The criticality is analysed with the renormalisation group. The low energy excitations at the topological transition are gapless semi-Dirac fermionic quasiparticles that disperse linearly in one direction and quadratically in the other. The instabilities induced by generic short-range repulsions are investigated. The criticality of this anisotropic system is analysed with the renormalisation group in two spatial dimensions, and it is found to exhibit anisotropic order parameter correlations.
The fluctuations-driven continuous quantum criticality has sparked tremendous interest in condensed matter physics. It has been verified that the fluctuations of gapless fermions can change the nature of phase transition at criticality. In this paper, we study the fermionic quantum criticality with enlarged Ising
×Ising fluctuations in honeycomb lattice materials. The perturbative renormalization group approach is employed to investigated the flow equations and stable fixed-point structure. By including cubic terms of order parameter, we found that the phase transition is always continuous for any finite value of the number of Dirac fermions
Nf. Further, we also compute the critical exponents and predict the critical scaling behavior for
Nf=2 case which corresponds to spin-1/2 fermions on graphene-like materials.