Monte-Carlo study of Dirac semimetals phase diagram

In this paper the phase diagram of Dirac semimetals is studied within lattice Monte-Carlo simulation. In particular, we concentrate on the dynamical chiral symmetry breaking which results in semimetal/insulator transition. Using numerical simulation we determined the values of the critical coupling constant of the semimetal/insulator transition for different values of the anisotropy of the Fermi velocity. This measurement allowed us to draw tentative phase diagram for Dirac semimetals. It turns out that within the Dirac model with Coulomb interaction both Na$_3$Bi and Cd$_3$As$_2$ known experimentally to be Dirac semimetals would lie deeply in the insulating region of the phase diagram. It probably shows a decisive role of screening of the interelectron interaction in real materials, similar to the situation in graphene.