We present a nonperturbative study of the quantum many-body effects caused by the long-range Coulomb interaction in a two-dimensional semi-Dirac semimetal. This kind of semimetal may be realized in deformed graphene and a class of other realistic materials. In the noninteracting limit, the dispersion of semi-Dirac fermion is linear in one direction and quadratic in the other direction. When the impact of Coulomb interaction is taken into account, such a dispersion can be significantly modified. To reveal the correlation effects, we first obtain the exact self-consistent Dyson-Schwinger equation of the full fermion propagator and then extract the momentum dependence of the renormalized fermion dispersion from the numerical solutions. Our results show that the fermion dispersion becomes linear in two directions. These results are compared to previous theoretical works on semi-Dirac semimetals.
We present a nonperturbative study of the quantum many-body effects caused by the long-range Coulomb interaction in a two-dimensional semi-Dirac semimetal. This kind of semimetal may be realized in deformed graphene and a class of other realistic materials. In the noninteracting limit, the dispersion of semi-Dirac fermion is linear in one direction and quadratic in the other direction. When the impact of Coulomb interaction is taken into account, such a dispersion can be significantly modified. To reveal the correlation effects, we first obtain the exact self-consistent Dyson-Schwinger equation of the full fermion propagator and then extract the momentum dependence of the renormalized fermion dispersion from the numerical solutions. Our results show that the fermion dispersion becomes linear in two directions. These results are compared to previous theoretical works on semi-Dirac semimetals.
Studying the strong correlation effects in interacting Dirac fermion systems is one of the most challenging problems in modern condensed matter physics. The long-range Coulomb interaction and the fermion-phonon interaction can lead to a variety of intriguing properties. In the strong-coupling regime, weak-coupling perturbation theory breaks down. The validity of 1/N expansion with N being the fermion flavor is also in doubt since N equals to 2 or 4 in realistic systems. Here, we investigate the interaction between (1+2)- and (1+3)-dimensional massless Dirac fermions and a generic scalar boson, and develop an efficient nonperturbative approach to access the strong-coupling regime. We first derive a number of self-consistently coupled Ward-Takahashi identities based on a careful symmetry analysis and then use these identities to show that the full fermion-boson vertex function is solely determined by the full fermion propagator. Making use of this result, we rigorously prove that the full fermion propagator satisfies an exact and self-closed Dyson-Schwinger integral equation, which can be solved by employing numerical methods. A major advantage of our nonperturbative approach is that there is no need to employ any small expansion parameter. Our approach provides a unified theoretical framework for studying strong Coulomb and fermion-phonon interactions. It may also be used to approximately handle the Yukawa coupling between fermions and order-parameter fluctuations around continuous quantum critical points. Our approach is applied to treat the Coulomb interaction in undoped graphene. We find that the renormalized fermion velocity exhibits a logarithmic momentum dependence but is nearly energy independent, and that no excitonic gap is generated by the Coulomb interaction. These theoretical results are consistent with experiments in graphene.
Recent experimental [I. Jo et al ., Phys. Rev. Lett. 119, 016402 (2017)] and numerical [M. Ippoliti, S. D. Geraedts, R. N. Bhatt, Phys. Rev. B 95, 201104 (2017)] evidence suggests an intriguing universal relationship between the Fermi surface anisotropy of the noninteracting parent 2-dimensional (2D) electron gas and the strongly correlated composite Fermi liquid formed in a strong magnetic field close to half-filling. Inspired by these observations, we explore more generally the question of anisotropy renormalization in interacting 2D Fermi systems. Using a recently developed [H. -K. Tang et al ., Science 361, 570 (2018)] nonperturbative and numerically exact projective quantum Monte Carlo simulation as well as other numerical and analytic techniques, only for Dirac fermions with long-range Coulomb interactions do we find a universal square-root decrease of the Fermi-surface anisotropy. For the ν = 1 / 2 composite Fermi liquid, this result is surprising since a Dirac fermion ground state was only recently proposed as an alternative to the usual Halperin–Lee–Read state. Our proposed universality can be tested in several anisotropic Dirac materials including graphene, topological insulators, organic conductors, and magic-angle twisted bilayer graphene.
Two-dimensional semi-Dirac fermions are quasiparticles that disperse linearly in one direction and quadratically in the other. We investigate instabilities of semi-Dirac fermions toward charge and spin density wave and superconducting orders, driven by short-range interactions. We analyze the critical behavior of the Yukawa theories for the different order parameters using Wilson momentum shell renormalization group. We generalize to a large number Nf of fermion flavors to achieve analytic control in 2+1 dimensions and calculate critical exponents at one-loop order, systematically including 1/Nf corrections. The latter depend on the specific form of the bosonic infrared propagator in 2+1 dimensions, which needs to be included to regularize divergencies. The 1/Nf corrections are surprisingly small, suggesting that the expansion is well controlled in the physical dimension. The order parameter correlations inherit the electronic anisotropy of the semi-Dirac fermions, leading to correlation lengths that diverge along the spatial directions with distinct exponents, even at the mean-field level. We conjecture that the proximity to the critical point may stabilize novel modulated order phases.
Topological quantum phase transitions intrinsically intertwine self-similarity and topology of many-electron wave-functions, and divining them is one of the most significant ways to advance understanding in condensed matter physics. Our focus is to investigate an unconventional class of the transitions between insulators and Dirac semimetals whose description is beyond conventional pseudo relativistic Dirac Hamiltonian. At the transition without the long-range Coulomb interaction, the electronic energy dispersion along one direction behaves like a relativistic particle, linear in momentum, but along the other direction it behaves like a non-relativistic particle, quadratic in momentum. Various physical systems ranging from TiO2-VO2 heterostructure to organic material α-(BEDT-TTF)2I3 under pressure have been proposed to have such anisotropic dispersion relation. Here, we discover a novel quantum criticality at the phase transition by incorporating the long range Coulomb interaction. Unique interplay between the Coulomb interaction and electronic critical modes enforces not only the anisotropic renormalization of the Coulomb interaction but also marginally modified electronic excitation. In connection with experiments, we investigate several striking effects in physical observables of our novel criticality.
An anisotropic flat band fermion system with a novel dispersion that is linear along one direction and cubic along another is proposed in Phys. Rev. X. 13, 021012 (2023). We study the effects of Coulomb interaction in this fermion system by renormalization group theory and Dyson-Schwinger gap equation. We perform renormalizaton group analysis and find that fermion velocity is always restored along the direction that the fermions take cubic dispersion originally. Accordingly, the system takes the similar behaviors to the two-dimensional Dirac fermion system with Coulomb interaction in the low energy regime. Based on Dyson-Schwinger gap equation method, we find that an excitonic gap is generated if the Coulomb strength is large enough, and the system becomes a novel excitonic Chern insulator with quantized anomalous Hall conductivity. Observable quantities of this system in free case, under weak and strong enough Coulomb interaction are all analyzed.
After its first isolation by Geim and Novoselov in 2004, graphene has grown quickly into an independent subfield of contemporary condensed matter physics, counting thousands of publications. This material has attracted so much attention not only because it is truly two-dimensional -- a feature long thought to be impossible to be realized in nature -- but also because it is the prime representative of a new class of condensed matter systems, the Dirac materials, whose electronic spectrum disperses linearly around isolated points in the Brillouin zone. An outstanding consequence of this linear energy-momentum dispersion is the collapse of the Fermi surface to isolated points at charge neutrality, accompanied by a vanishing single-particle density of states. While this feature of the noninteracting model is what separates graphene from most other, more conventional condensed matter systems, giving rise to many of its astounding physical properties, it leads to severe computational complications when extending the model to contain two-body interactions and/or disorder. In this thesis we consider both extensions separately. After rigorously deriving the low-energy quantum field theory from a tight-binding model and establishing the one-particle irreducible (1PI) vertex functions as the fundamental building blocks of correlation functions, we construct a nonperturbative formalism for the calculation of these vertex functions using the well-established functional renormalization group (fRG). By combining the fRG with the Keldysh formalism, accounting for possible external electromagnetic fields and nontrivial initial correlations, this framework is capable of handling thermal equilibrium and true nonequilibrium alike. To circumvent some of the technical problems of the standard fRG formalism when treating fermionic systems at finite density, we also explore a variant of the fRG, where the chemical potential is interpreted as a flow parameter. We obtain hierarchical sets of flow equations, which describe the change of the 1PI vertex functions upon varying an artificial cutoff scale or the chemical potential, respectively. These nonperturbative formalisms have been used to calculate the renormalization of the Fermi velocity and the static dielectric function at finite temperature and density in the strong coupling regime as an explicit demonstration of their capabilities. Since the Keldysh fRG is an exact formulation of quantum field theory, in particular when external magnetic fields are present, it should be possible to describe the fractional quantum Hall effect in this framework. However, the typical truncation schemes that are currently available to approximately solve the fRG flow equation fail to account for the nontrivial correlation effects necessary to access the fractional quantum Hall regime. To circumvent this problem we consider a modified field theory, where these correlations are implemented ``by hand'' via a Chern-Simons gauge field, coupling to the fermions in addition to the Coulomb interaction term. This modified field theory is analyzed in a stationary phase approximation including Gaussian fluctuations. The electromagnetic response tensor is calculated in the random phase approximation, yielding the Hall conductivities. In the second part of this thesis we consider disordered Dirac fermions in the absence of two-particle interactions. After sketching how the most general disorder potential for Dirac fermions arises from a tight-binding model, we first consider an explicit physical scenario, a disordered graphene pn junction in the presence of a quantizing magnetic field perpendicular to the graphene sheet. We derive an effective one-dimensional theory for the chiral states that propagate along the junction interface and calculate the full conductance distribution in the crossover between the clean and the strong-disorder limit, via an exact solution of a Fokker-Planck equation. Lastly, we develop a nonperturbative approach to calculate the disorder induced self-energy that is based on exact Schwinger-Dyson equations and Ward identities, not only for graphene, but also for other paradigmatic semimetals with a nodal point dispersion.