No Arabic abstract
We use determinant quantum Monte Carlo (DQMC) simulations to study the role of electron-electron interactions on three-dimensional (3D) Dirac fermions based on the $pi$-flux model on a cubic lattice. We show that the Hubbard interaction drives the 3D Dirac semimetal to an antiferromagnetic (AF) insulator only above a finite critical interaction strength and the long-ranged AF order persists up to a finite temperature. We evaluate the critical interaction strength and temperatures using finite-size scaling of the spin structure factor. The critical behaviors are consistent with the (3+1)d Gross-Neveu universality class for the quantum critical point and 3D Heisenberg universality class for the thermal phase transitions. We further investigate correlation effects in birefringent Dirac fermion system. It is found that the critical interaction strength $U_c$ is decreased by reducing the velocity of the Dirac cone, quantifying the effect of velocity on the critical interaction strength in 3D Dirac fermion systems. Our findings unambiguously uncover correlation effects in 3D Dirac fermions, and may be observed using ultracold atoms in an optical lattice.
The analogues of elementary particles have been extensively searched for in condensed matter systems because of both scientific interests and technological applications. Recently massless Dirac fermions were found to emerge as low energy excitations in the materials named Dirac semimetals. All the currently known Dirac semimetals are nonmagnetic with both time-reversal symmetry $mathcal{T}$ and inversion symmetry $mathcal{P}$. Here we show that Dirac fermions can exist in one type of antiferromagnetic systems, where $mathcal{T}$ and $mathcal{P}$ are broken but their combination $mathcal{PT}$ is respected. We propose orthorhombic antiferromagnet CuMnAs as a candidate, analyze the robustness of the Dirac points with symmetry protections, and demonstrate its distinctive bulk dispersions as well as the corresponding surface states by emph{ab initio} calculations. Our results give a new route towards the realization of Dirac materials, and provide a possible platform to study the interplay of Dirac fermion physics and magnetism.
Two-dimensional Dirac fermions are subjected to two types of interactions, namely the long-range Coulomb interaction and the short-range on-site interaction. The former induces excitonic pairing if its strength $alpha$ is larger than some critical value $alpha_c$, whereas the latter drives an antiferromagnetic Mott transition when its strength $U$ exceeds a threshold $U_c$. Here, we study the impacts of the interplay of these two interactions on excitonic pairing with the Dyson-Schwinger equation approach. We find that the critical value $alpha_c$ is increased by weak short-range interaction. As $U$ increases to approach $U_c$, the quantum fluctuation of antiferromagnetic order parameter becomes important and interacts with the Dirac fermions via the Yukawa coupling. After treating the Coulomb interaction and Yukawa coupling interaction on an equal footing, we show that $alpha_c$ is substantially increased as $U rightarrow U_c$. Thus, the excitonic pairing is strongly suppressed near the antiferromagnetic quantum critical point. We obtain a global phase diagram on the $U$-$alpha$ plane, and illustrate that the excitonic insulating and antiferromagnetic phases are separated by an intermediate semimetal phase. These results provide a possible explanation of the discrepancy between recent theoretical progress on excitonic gap generation and existing experiments in suspended graphene.
We study a class of Dirac semimetals that feature an eightfold-degenerate double Dirac point. We show that 7 of the 230 space groups can host such Dirac points and argue that they all generically display linear dispersion. We introduce an explicit tight-binding model for space groups 130 and 135, showing that 135 can host an intrinsic double Dirac semimetal -- one with no additional degeneracies at the Fermi energy. We consider symmetry-lowering perturbations and show that uniaxial compressive strain in different directions leads to topologically distinct insulating phases. In addition, the double Dirac semimetal can accommodate topological line defects that bind helical modes. Potential materials realizations are discussed.
The spontaneous generation of charge-density-wave order in a Dirac fermion system via the natural mechanism of electron-phonon coupling is studied in the framework of the Holstein model on the honeycomb lattice. Using two independent and unbiased quantum Monte Carlo methods, the phase diagram as a function of temperature and coupling strength is determined. It features a quantum critical point as well as a line of thermal critical points. Finite-size scaling appears consistent with fermionic Gross-Neveu-Ising universality for the quantum phase transition, and bosonic Ising universality for the thermal phase transition. The critical temperature has a maximum at intermediate couplings. Our findings motivate experimental efforts to identify or engineer Dirac systems with sufficiently strong and tunable electron-phonon coupling.
We study the role of long-range electron-electron interactions in a system of two-dimensional anisotropic Dirac fermions, which naturally appear in uniaxially strained graphene, graphene in external potentials, some strongly anisotropic topological insulators, and engineered anisotropic graphene structures. We find that while for small interactions and anisotropy the system restores the conventional isotropic Dirac liquid behavior, strong enough anisotropy can lead to the formation of a quasi-one dimensional electronic phase with dominant charge order (anisotropic excitonic insulator).