Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2003.10371 · Physical review. B./Physical review. B

Nonperturbative Dyson-Schwinger equation approach to strongly interacting Dirac fermion systems

Xiao‐Yin Pan, Zhao-Kun Yang, Xin Li, Guo-Zhu Liu

Published Aug 23, 2021• 10 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physrevb.104.085141

Related papers

Fermion dispersion renormalization in a two-dimensional semi-Dirac semimetal

Hao-Fu Zhu, Xiao‐Yin Pan, Guo-Zhu Liu

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.

arXiv:2112.06440•Physical review. B./Physical review. B•Feb 22, 2022•3 citations

Field Theories of Interacting and Disordered Dirac Fermions in Graphene

Christian Fräßdorf

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.

•Jan 1, 2019•0 citations

Bosonization of the Q=0 continuum of Dirac fermions

S. F. Mantilla, Inti Sodemann

We develop a bosonization formalism that captures nonperturbatively the interaction effects on the $\mathbf{Q}=0$ continuum of excitations of nodal fermions above one dimension. Our approach is a natural extension of the classic bosonization scheme for higher dimensional Fermi surfaces to include the $\mathbf{Q}=0$ neutral excitations that would be absent in a single-band system. The problem is reduced to solving a boson bilinear Hamiltonian. We establish a rigorous microscopic footing for this approach by showing that the solution of such boson bilinear Hamiltonian is exactly equivalent to performing the infinite sum of Feynman diagrams associated with the Kadanoff-Baym particle-hole propagator that arises from the self-consistent Hartree-Fock approximation to the single-particle Green's function. We apply this machinery to compute the interaction corrections to the optical conductivity of two-dimensional Dirac fermions with Coulomb interactions reproducing the results of perturbative renormalization group at weak coupling and extending them to the strong-coupling regime.

arXiv:2002.05732•Physical review. B./Physical review. B•Sep 8, 2020•1 citation

Correlated interaction effects in an anisotropic flat band fermion system

Jing-Rong Wang, Changjin Zhang

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.

arXiv:2412.13980•Dec 18, 2024•0 citations

Stability of Dirac Liquids with Strong Coulomb Interaction

I. S. Tupitsyn, Nikolay V. Prokof'ev

We develop and apply the diagrammatic Monte Carlo technique to address the problem of the stability of the Dirac liquid state (in a graphene-type system) against the strong long-range part of the Coulomb interaction. So far, all attempts to deal with this problem in the field-theoretical framework were limited either to perturbative or random phase approximation and functional renormalization group treatments, with diametrically opposite conclusions. Our calculations aim at the approximation-free solution with controlled accuracy by computing vertex corrections from higher-order skeleton diagrams and establishing the renormalization group flow of the effective Coulomb coupling constant. We unambiguously show that with increasing the system size L (up to ln(L)∼40), the coupling constant always flows towards zero; i.e., the two-dimensional Dirac liquid is an asymptotically free T=0 state with divergent Fermi velocity.

arXiv:1608.00133•Physical Review Letters•Jan 12, 2017•43 citations

Quantum critical behavior in a graphenelike model

Simon Hands, Costas Strouthos

We present the results of numerical simulations of a $(2+1)$-dimensional fermion field theory based on a recent proposal for a model of graphene consisting of ${N}_{f}$ four-component Dirac fermions moving in the plane and interacting via an instantaneous Coulomb interaction. In the strong-coupling limit we identify a critical number of flavors ${N}_{fc}=4.8(2)$ separating an insulating from a conducting phase. This transition corresponds to the location of a quantum critical point, and we use a fit to the equation of state for the chiral order parameter to estimate the critical exponents. Next we simulate ${N}_{f}=2$ corresponding to real graphene and approximately locate a transition from strong- to weak-coupling behavior. Strong correlations are evident in the weak-coupling regime.

arXiv:0806.4877•Physical Review B•Oct 24, 2008•140 citations
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/N1/N1/N expansion with NNN being the fermion flavor is also in doubt since NNN equals to 2 or 4 in realistic systems. Here, we investigate the interaction between (1+2)(1+2)(1+2)- and (1+3)(1+3)(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.
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.
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.
We develop a bosonization formalism that captures nonperturbatively the interaction effects on the Q=0\mathbf{Q}=0Q=0 continuum of excitations of nodal fermions above one dimension. Our approach is a natural extension of the classic bosonization scheme for higher dimensional Fermi surfaces to include the Q=0\mathbf{Q}=0Q=0 neutral excitations that would be absent in a single-band system. The problem is reduced to solving a boson bilinear Hamiltonian. We establish a rigorous microscopic footing for this approach by showing that the solution of such boson bilinear Hamiltonian is exactly equivalent to performing the infinite sum of Feynman diagrams associated with the Kadanoff-Baym particle-hole propagator that arises from the self-consistent Hartree-Fock approximation to the single-particle Green's function. We apply this machinery to compute the interaction corrections to the optical conductivity of two-dimensional Dirac fermions with Coulomb interactions reproducing the results of perturbative renormalization group at weak coupling and extending them to the strong-coupling regime.
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.
We develop and apply the diagrammatic Monte Carlo technique to address the problem of the stability of the Dirac liquid state (in a graphene-type system) against the strong long-range part of the Coulomb interaction. So far, all attempts to deal with this problem in the field-theoretical framework were limited either to perturbative or random phase approximation and functional renormalization group treatments, with diametrically opposite conclusions. Our calculations aim at the approximation-free solution with controlled accuracy by computing vertex corrections from higher-order skeleton diagrams and establishing the renormalization group flow of the effective Coulomb coupling constant. We unambiguously show that with increasing the system size L (up to ln(L)∼40), the coupling constant always flows towards zero; i.e., the two-dimensional Dirac liquid is an asymptotically free T=0 state with divergent Fermi velocity.
We present the results of numerical simulations of a (2+1)(2+1)(2+1)-dimensional fermion field theory based on a recent proposal for a model of graphene consisting of Nf{N}_{f}Nf​ four-component Dirac fermions moving in the plane and interacting via an instantaneous Coulomb interaction. In the strong-coupling limit we identify a critical number of flavors Nfc=4.8(2){N}_{fc}=4.8(2)Nfc​=4.8(2) separating an insulating from a conducting phase. This transition corresponds to the location of a quantum critical point, and we use a fit to the equation of state for the chiral order parameter to estimate the critical exponents. Next we simulate Nf=2{N}_{f}=2Nf​=2 corresponding to real graphene and approximately locate a transition from strong- to weak-coupling behavior. Strong correlations are evident in the weak-coupling regime.