Arxiv:1803.11480V3 [Cond-Mat.Str-El] 23 May 2020 Submitted To: J
Total Page:16
File Type:pdf, Size:1020Kb
Interacting Dirac Materials S. Banerjee1;2;3, D.S.L. Abergel2, H. Ågren3;4, G. Aeppli5;6;7, A.V. Balatsky2;8 1 Theoretical Physics III, Center for Electronic Correlations and Magnetism, Institute of Physics, University of Augsburg, 86135 Augsburg Germany 2Nordita, KTH Royal Institute of Technology and Stockholm University, Roslagstullsbacken 23, 10691 Stockholm, Sweden 3Division of Theoretical Chemistry and Biology, Royal Institute of Technology, SE-10691 Stockholm, Sweden 4 Department of Physics and Astronomy, Uppsala University, Box 516, SE-751 20 Uppsala, Sweden 5 Paul Scherrer Institute, CH-5232 Villigen PSI, Switzerland 6 Laboratory for Solid State Physics, ETH Zurich, Zurich, CH-8093, Switzerland 7 Institut de Physique, EPF Lausanne, Lausanne, CH-1015, Switzerland 8 Department of Physics, University of Connecticut, Storrs, Connecticut 06269, USA E-mail: [email protected] 26 May 2020 Abstract. We investigate the extent to which the class of Dirac materials in two-dimensions provides general statements about the behavior of both fermionic and bosonic Dirac quasiparticles in the interacting regime. For both quasiparticle types, we find common features for the interaction induced renormalization of the conical Dirac spectrum. We perform the perturbative renormalization analysis and compute the self-energy for both quasiparticle types with different interactions and collate previous results from the literature whenever necessary. Guided by the systematic presentation of our results in Table 1, we conclude that long- range interactions generically lead to an increase of the slope of the single-particle Dirac cone, whereas short-range interactions lead to a decrease. The quasiparticle statistics does not qualitatively impact the self-energy correction for long-range repulsion but does affect the behavior of short-range coupled systems, giving rise to different thermal power-law contributions. The possibility of a universal description of the Dirac materials based on these features is also mentioned. Keywords: Fermi liquid, Dirac bosons, perturbation theory, quasiparticle arXiv:1803.11480v3 [cond-mat.str-el] 23 May 2020 Submitted to: J. Phys.: Condens. Matter Interacting Dirac Materials 2 1. Introduction materials bear interesting similarities: both have symmetry protected band crossing points around which The discovery of graphene [1, 2, 3] and other materials the excitations have linear dispersion parameterized by with Dirac nodes has led to a strong interest in a velocity. However, for fermions, these are typically materials with Dirac dispersion [4, 5, 6] in condensed low-energy excitations around the Fermi energy at the matter physics. In a typical d-dimensional metal, the Dirac nodes; whereas the bosonic Dirac excitations lie Fermi surface is a (d − 1)-dimensional manifold. In at a higher energy from the ground state (see Fig. 1). Dirac materials, on the other hand, the locus of points Several important questions about this new class with zero quasiparticle energy shrinks to d − 2 or of materials remain unanswered: is there a universal d − 3 dimensions due to some additional symmetries behavior for bosonic and fermionic Dirac excitations present in the system [4, 6]. Specifically, for two- when interactions are taken into account? What dimensional (2D) Dirac materials, Dirac nodes have role does the form of the interaction play? And to be point objects. The low-energy quasiparticles how are these results modified at finite temperature? around this Fermi point have linear dispersion and are The first question is particularly important because described by a massless Dirac-like equation. Because the quasiparticle statistics will have a strong impact of the linear dispersion, the quasiparticles behave like on interacting behavior. Our results in Sec. 3, relativistic particles with an effective speed of light taken together with the existing literature in the replaced by the Fermi velocity. Many realizations of field of Dirac materials, allows us to extract certain Dirac fermions have been discussed in, for example, striking similarities between the renormalization of superfluid phases of 3He [7], the surface states of fermionic and bosonic Dirac materials. On a topological insulators [8, 9, 10], high-temperature qualitative level, long-range interactions (Coulomb d-wave superconductors [11], graphene [12] and type) lead to an increase in the slope of the several artificial electronic systems [13, 14]. The Dirac cone, whereas short-range interactions decrease resulting Dirac cone in the spectrum is protected by this slope. Furthermore, both the fermionic and symmetries [15] and is robust under disorder and bosonic excitations acquire logarithmic temperature- interaction. For example, Dirac nodes are protected dependent self-energy corrections (see Table 1). Such by parity and time-reversal symmetries in graphene [3], temperature dependence arises from Dirac dispersions by time-reversal and mirror symmetries in topological in the fermionic illustration, whereas for bosons it insulators [8] and by a combination of the C4 rotation is associated with the typical parabolic dispersions. and gauge transformation in the case of d-wave On the other hand, different power-law corrections superconductors in two-dimensional square lattices [4]. in temperature and momentum are obtained for the The Dirac nodes lead to commonly observed universal renormalized dispersions for short-range interactions. features of excitations such as the scalingp of the Table 1 summarizes the trends in renormalization Landau levels with the magnetic field as B [3], of Dirac node velocity. These trends indicate common universal metallic and thermal conductivities related features in velocity dependence in short and long range to the nodal structure of carriers, suppressed back- interacting cases. One indeed might suspect there scattering [4], and Klein tunneling [16]. The common are universal trends in renormalzations of Dirac nodes properties shared by these systems bring about the depending on dimensionality and nature of interaction. unified concept of Dirac materials as a useful category To determine if indeed these trends are universal one for explaining their behavior. would need to investigate a broader class of models in Recently, we are witnessing an upsurge of D = 1,2,3. This investigation of broader class of models theoretical and experimental investigations in bosonic is left for the future work. systems possessing Dirac-like spectra. Plausible In this paper, we analyze the renormalized physical realizations of such bosonic Dirac materials dispersions and the induced quasiparticle lifetimes range from artificial honeycomb lattices made out for the weakly interacting two-dimenisonal Dirac of superconducting grains [17] to magnets [18], excitations. We therefore set up our models for : photonic and acoustic crystals [19, 20], and plasmonic (a) Dirac fermions interacting via an onsite Hubbard devices [21]. At the level of the single-particle physics, repulsion and (b) Dirac bosons subject to a long-range the band structures of bosonic and fermionic Dirac Coulomb interaction. We purposefully recollect the Interacting Dirac Materials 3 (a) (b) various fermionic Dirac excitations; whereas various Refs. [17, 18, 19, 20, 21] (and more) investigated the bosonic counterparts. Band-structures like the ones shown in Fig. 1 containing Dirac-like nodes typically arise from a tight-binding model defined on the non- symmorphic lattice such as honeycomb, kagome, etc. A generic Hamiltonian, in such a case, can be written as X y y X y y Figure 1. Comparison of the band-structure for non-interacting H = −t ci cj(bi bj) + h:c: ∓ "F ("B) ci ci(bi bi); Dirac bosons (left panel) vs fermions (right panel). (a) The hiji i energy spectrum for a single magnon excitation induced by a ferromagnetic Heisenberg model defined on a honeycomb lattice (1) y y (see Eq. 15 and Refs. [18, 22]). The magnons can occupy any where ci (bi ) is the creation operator for the fermionic state belonging to the two bands which cross each other at (bosonic) quasiparticles, like electrons (magnons, the Dirac points (marked in red and black points). The near phonons, plasmons, etc.), at site i within the unit cell vicinity area of the spectrum around these crossing points can be effectively described by a Dirac-like Hamiltonian. (b) In of the lattice. t- is the amplitude for the quasiparticle contrast, the Dirac fermionic excitations are described by the hopping between the nearest-neighbor sites hiji. The low-lying states from the Fermi energy " at the Dirac node F lattice constant is denoted by a. "F ("B) denotes within a bandwidth ±W . the Fermi (an onsite) energy for fermionic (bosonic) quasiparticles. relevant results for the interacting Dirac magnons (see For example, in a honeycomb lattice, the Ref. [22]) as a prototype for Dirac bosons coupled Hamiltonian in Eq. 1 gives rise to band-structures as via short-range interaction. In a similar spirit, any shown in Fig. 1. In particular, for magnons (bosonic explicit analysis for the renormalization of the Dirac excitations, t = JS), one obtains a band-structure fermions subject to long-range interactions is dropped. as in Fig. 1a (see Ref. [22] for the details). The Instead, we rely on the well-known results obtained onsite energy "B(= 3JS) in this case, accounts for the in Refs. [3, 23]. We do not consider the role of positivity of the magnonic dispersions. In contrast, spin-orbit coupling or similar interactions that might for fermionic