Arxiv:1901.05432V3 [Cond-Mat.Str-El] 11 May 2019
Total Page:16
File Type:pdf, Size:1020Kb
Magic of high-order van Hove singularity Noah F. Q. Yuan, Hiroki Isobe, and Liang Fu Department of Physics, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA We introduce a new type of van Hove singularity in two dimensions, where a saddle point in momentum space is changed from second-order to high-order. Correspondingly, the density of states near such \high-order van Hove singularity" is significantly enhanced from logarithmic to power-law divergence, which promises stronger electron correlation effects. High-order van Hove singularity can be generally achieved by tuning the band structure with a single parameter in moir´esuperlattices, such as twisted bilayer graphene by tuning twist angle or applying pressure, and trilayer graphene by applying vertical electric field. The recent discovery of unconventional insulating (a) 50 states and superconductivity in twisted bilayer graphene (b) 45 95 (TBG) [1, 2] has attracted enormous interest. At am- bient pressure, these intriguing phenomena of correlated 40 90 electrons only occur near a specific twist angle θ ≈ 1:1◦, 35 85 30 widely referred to as the magic angle. The existence of 80 (nS) such a magic angle was first predicted from band struc- 25 c G(nS) 75 ture calculations. Based on a continuum model [3, 4], 20 G+G Bistritzer and MacDonald showed [5] that the lowest 15 70 moir´ebands in TBG become exceptionally flat at this 10 =1.32 2 65 magic angle, and are hence expected to exhibit strong 5 electron correlation. This pioneering work inspired a 0 -60 -40 -20 0 20 40 60 1 2 4 6 8 10 large body of experimental works in recent years, which E(meV) |E-E |(meV) v culminated in the discovery reported in Ref. [1, 2]. The origin of magic-angle phenomena, i.e., the emergence of superconductivity and correlated insulator, is now under FIG. 1: Theoretical fit to the tunneling conductance mea- intensive study [6{25]. surement [26]. (a) Open circles are tunneling conductance G of twisted bilayer graphene at twist angle 1:10◦, and solid While the reduction of bandwidth is undoubtedly im- lines are fitting according to Eq. (5) with details given in the portant, it is evident that the phenomenology of magic- Supplemental Material. Dashed lines denote singularity ener- angle TBG cannot be ascribed to a completely flat gies, indicating the asymmetry of conductance peaks. (b) The band. Quantum oscillation at low magnetic fields re- peak at Ev = 16:72 meV plotted in logarithmic scales. Elec- veals doping-dependent Fermi surfaces within the narrow tron and hole sides of the peak fall into two parallel lines withp moir´eband. Detailed scanning tunneling spectroscopy the same slope −1=4 and asymmetry ratio η = 1:32 ≈ 2. (STS) studies [26{29] and compressibility measurements The only parameter in (b) is the background offset Gc = 57:6 nS. [30] show that the moir´ebandwidth at a magic angle is a few tens of meV, larger than what previous calcula- tions found [5, 31{36]. Moreover, there is no direct evi- recent STS measurements on gate-tunable TBG around dence from existing STS measurements at various twist magic angle [26, 27] find that when the VHS shifts to angles that the moir´eband is exceptionally flat right at the Fermi level under gating (the corresponding density the magic angle (see also Ref. [29]). These experimen- is within 10% of half filling), the VHS peak in tunneling tal results motivated us to consider additional feature of density of states splits into two new peaks (see also Ref. electronic structure which may create favorable condition [46]). These findings clearly demonstrate the prominent for correlated electron phenomena in magic-angle TBG. role of VHS in magic-angle TBG. In this work, we propose a new perspective that re- On the other hand, it is also known from general con- arXiv:1901.05432v3 [cond-mat.str-el] 11 May 2019 lates magic angle to a new type of van Hove singularity sideration and previous STS measurements [47{49] that (VHS) in the single-particle energy spectrum of TBG. VHS are present at all twist angles. It is therefore unclear The importance of VHS has already been recognized in whether VHS at magic angle is anything special. our theory of correlated TBG [37, 38] and related stud- We propose the following scenario. As the twist an- ies [39{45] using the weak coupling approach. Generally gle decreases below a critical value θc, the van Hove speaking, VHS in two-dimensional systems are associated saddle point|which marks the change of topology in with saddle points of energy dispersion in k space. When Fermi surface (Lifshitz transition)|undergoes a topolog- a VHS is close to Fermi energy, the increased density of ical transition whereby a single saddle point splits into states (DOS) amplifies electron correlation, resulting in two new ones. Right at θc, the saddle point changes various ordering instabilities such as density wave and su- from second-order to higher-order; the DOS at VHS is perconductivity at low temperature [37]. Indeed, the two significantly enhanced from logarithmic to power-law di- 2 vergence, which promises stronger electron correlation. of two Fermi surfaces, but they touch tangentially (gen- We propose that proximity to such \high-order van Hove erally speaking) rather than intersect at a finite angle as singularity", which requires tuning to the critical twist in the case of ordinary VHS. angle or pressure, is an important factor responsible for VHS in TBG. We now turn to VHS in TBG, whose correlated electron phenomena in TBG near half filling. low-energy moir´ebands arise from inter-layer coupling of We demonstrate by topological argument that high- Dirac states on the two graphene layers. Since single- order VHS can be generally achieved by tuning the band particle scattering between K and K0 points requires structure with just a single control parameter. For TBG, large momentum transfer, it is suppressed at small twist besides twist angle, pressure can also induce supercon- angles due to the long wavelength of moir´epotential. In ductivity at a twist angle larger than 1:1◦ [50]. By track- the absence of valley hybridization, the moir´ebands from ing the evolution of the moir´eband with twist angle and K and K0 valleys are decoupled [3{5, 56], with energy pressure, we locate the high-order VHS in experimentally dispersions denoted by E+(k) and E−(k) respectively. relevant parameter regime and predict its key feature, a Time-reversal symmetry implies E−(k) = E+(−k), so it distinctive asymmetric peak in DOS. This feature com- suffices to consider E+(k) only in the following. pares well with the experimentally observed VHS peak The number and location of VHS points in TBG de- in magic-angle TBG, as shown in Fig. 1. We also discuss pend on the twist angle. For example, the band struc- the splitting of VHS peak by nematic or density wave ture calculation for θ = 2◦ [57] reveals the existence of order. three symmetry-related van Hove saddle points on ΓM Ordinary and high-order VHS. In two-dimensional lines in the mini-Brillouin zone (MBZ), shown in Fig. (2D) electron systems with energy dispersion E(k), an 2a. Across the VHS energy, the Fermi contour changes ordinary VHS with logarithmically diverging DOS occurs from two disjoint Dirac pockets around MBZ corners to at a saddle point ks, determined by a single pocket enclosing the MBZ center. This leads to the conversion between electron and hole charge carriers, rkE(ks) = 0 and det D < 0; (1) as evidenced by the sign change of Hall coefficient [57]. On the other hand, at θ = 1:05◦, on each ΓM line there 1 where Dij ≡ 2 @i@jE is the 2 × 2 Hessian matrix of E at is a local energy maximum|instead of saddle point|in ks. Since D is symmetric by definition, we can rotate the the moir´evalence band. As doping increases, an addi- axes to diagonalize D and Taylor expand the energy dis- tional hole pocket (not present at θ = 2◦) emerges out 2 2 persion near ks as E −Ev = −αpx +βpy, where Ev is the of each energy maximum and eventually intersects two VHS energy, the momentum p = k−ks is measured from Dirac pockets at two new saddle points on opposite sides the saddle point, and the coefficients −α, β are the two of ΓM [33]. In such band structure there are a total of 6 eigenvalues of D with −αβ = det D < 0. This disper- VHS points, shown in Fig. 2c. sion describes two pieces of Fermi contours that intersect By continuity, we deduce from this change in the num- at the saddle point ks. It is known from a topological ber of VHS points that a topological transition of saddle consideration that due to the periodicity of E(k) on the points must occur, at a hitherto unknown critical twist Brillouin zone (a torus), van Hove saddle points appear angle, in such a way that a saddle point on ΓM (denoted quite generally in energy bands of two-dimensional ma- by Λ0) splits into a pair of new ones off ΓM (denoted terials [51]. by Λ±). Importantly, the behavior of these VHS points We now introduce the notion of high-order VHS by the in the vicinity of this transition is solely governed by the following conditions at a saddle point local energy dispersion near Λ0. We now expand E(k) near Λ0 to higher orders: rkE(ks) = 0 and det D = 0; (2) 2 2 2 4 E − Ev = −αpx + βpy + γpxpy + κpy + :::; (3) which imply αβ = 0.