We investigate the interplay between charge order and superconductivity near an antiferromagnetic quantum critical point using sign-problem-free Quantum Monte Carlo simulations. We establish that, when the electronic dispersion is particle-hole symmetric, the system has an emergent SU(2) symmetry that implies a degeneracy between d-wave superconductivity and charge order with d-wave form factor. Deviations from particle-hole symmetry, however, rapidly lift this degeneracy, despite the fact that the SU(2) symmetry is preserved at low energies. As a result, we find a strong suppression of charge order caused by the competing, leading superconducting instability. Across the antiferromagnetic phase transition, we also observe a shift in the charge order wave vector from diagonal to axial. We discuss the implications of our results to the universal phase diagram of antiferromagnetic quantum-critical metals and to the elucidation of the charge order experimentally observed in the cuprates.
Bibliographical noteFunding Information:
Wang Xiaoyu 1,2 Wang Yuxuan 3 Schattner Yoni 4,5,6 Berg Erez 2 Fernandes Rafael M. 1 School of Physics and Astronomy, 1 University of Minnesota , Minneapolis, Minnesota 55455, USA Department of Physics, 2 University of Chicago , Chicago, Illinois 60637, USA Institute for Condensed Matter Theory and Department of Physics, 3 University of Illinois , Urbana-Champaign, Illinois 61801, USA 4 Weizmann Institute of Science , Rehovot 761000, Israel Department of Physics, 5 Stanford University , Stanford, California 94305, USA Stanford Institute for Materials and Energy Sciences, 6 SLAC National Accelerator Laboratory and Stanford University , Menlo Park, California 94025, USA 15 June 2018 15 June 2018 120 24 247002 5 October 2017 © 2018 American Physical Society 2018 American Physical Society We investigate the interplay between charge order and superconductivity near an antiferromagnetic quantum critical point using sign-problem-free Quantum Monte Carlo simulations. We establish that, when the electronic dispersion is particle-hole symmetric, the system has an emergent SU(2) symmetry that implies a degeneracy between d -wave superconductivity and charge order with d -wave form factor. Deviations from particle-hole symmetry, however, rapidly lift this degeneracy, despite the fact that the SU(2) symmetry is preserved at low energies. As a result, we find a strong suppression of charge order caused by the competing, leading superconducting instability. Across the antiferromagnetic phase transition, we also observe a shift in the charge order wave vector from diagonal to axial. We discuss the implications of our results to the universal phase diagram of antiferromagnetic quantum-critical metals and to the elucidation of the charge order experimentally observed in the cuprates. U.S. Department of Energy 10.13039/100000015 Office of Science 10.13039/100006132 Basic Energy Sciences 10.13039/100006151 DE-SC0012336 Research Corporation for Science Advancement 10.13039/100001309 University of Minnesota 10.13039/100007249 Gordon and Betty Moore Foundation 10.13039/100000936 GBMF4305 University of Illinois 10.13039/100010443 Division of Materials Sciences and Engineering 10.13039/100013055 DE-AC02-76SF00515 The phase diagrams of a number of strongly correlated materials display putative quantum critical points (QCP), in which the transition temperature of an electronically ordered state is suppressed to zero. In systems displaying antiferromagnetic (AFM) order, such as heavy fermions, cuprates, and iron pnictides, unconventional superconductivity (SC) is found to emerge near the QCP  . Although it is well established that the interactions mediated by fluctuations near an AFM QCP favor a sign-changing SC gap, the extent to which this physics describes the actual materials remains widely debated. In this context, analytical investigations of metallic AFM QCP in two dimensions revealed a surprising result: the same electronic interaction that promotes sign-changing SC also promotes an unusual sign-changing bond charge order (CO) [2–6] . This magnetic mechanism for CO is sharply distinct from the usual mechanisms involving phonons and Fermi surface nesting. Taken at face value, this result would suggest that CO should emerge generically in the phase diagrams of AFM systems. These theoretical results were brought to the spotlight by the experimental observation of the sign-changing bond CO in cuprate high- T c superconductors [7–19] , spurring many ideas on the interplay among AFM QCP, SC, and CO [5,6,20–26] . It has been proposed, for instance, that pseudogap physics is a manifestation of a more fundamental symmetry between SC and CO near an AFM QCP  . However, most of these theoretical works have relied on certain uncontrolled approximations, which are required for an analytical treatment of an AFM QCP in a metal. The fundamental question about the universality of CO near an AFM QCP, and the more specific question about the relevance of this result to explain charge order in cuprates, begs for unbiased methods to probe this phenomenon. In this Letter, we employ the determinantal Quantum Monte Carlo (QMC) method to address these questions. We consider the two-band version of the spin-fermion model, for which the QMC method does not suffer from the fermionic sign problem  . The model consists of free electrons coupled to an AFM order parameter tuned to its QCP  . Analytical and sign-problem-free QMC calculations have established the existence of a SC dome surrounding the QCP [28,29] . As for CO, an emergent SU(2) symmetry between CO and SC was found analytically within an approximation that considers only the vicinity of the AFM hot spots—the points on the Fermi surface separated by the AFM wave vector Q AFM = ( π , π ) [3,4] . The resulting CO wave vector lies along the diagonal of the Brillouin zone Q CO = ( Q 0 , Q 0 ) with 2 Q 0 being the distance between hot spots. CO with axial wave vectors ( Q 0 , 0 ) and/or ( 0 , Q 0 ) , which are those experimentally observed in cuprates, has also been proposed within the spin-fermion model [5,30,31] . Although QMC investigations have not yet found CO, they have focused on very narrow parameter regimes  or were performed in the superconducting phase  . Here, we report QMC results showing the existence of a SU(2) symmetry between CO and SC near the AFM QCP when the noninteracting band structure has particle-hole symmetry. This SU(2) symmetry implies a degeneracy between SC and CO, manifested by a sharp enhancement of both susceptibilities as the QCP is approached. This result demonstrates the nontrivial mechanism of magnetically mediated CO, and it establishes that the same interaction that promotes SC in the spin-fermion model is also capable of promoting CO. As the particle-hole symmetry is broken, however, we find that while the enhancement of the SC susceptibility is preserved, the CO susceptibility shows a very weak enhancement near the QCP. Furthermore, near the onset of SC, the CO susceptibility is even suppressed with respect to its noninteracting value, signaling a strong competition between these two states already in the fluctuating regime. This happens even though the SU(2) symmetry is preserved locally at the hot spots. The fragility of the CO-SC degeneracy implies that CO near an AFM QCP is not a universal phenomenon, but instead it requires a fine-tuned band structure that goes beyond just hot-spot properties. We also investigate the wave vector for which the CO susceptibility is maximal. When CO and SC are degenerate, the wave vector is diagonal, in agreement with the analytical approximations. However, once CO and SC are no longer degenerate, the wave vector tends to change from diagonal to axial as the AFM phase is approached. This is consistent with theoretical proposals that axial CO is favored over the diagonal one if the antinodal region of the Brillouin zone is gapped [33,34] . Finally, we discuss the implications of our results to materials that display putative AFM QCPs and their relevance to understand CO in cuprates. The spin-fermion model is a low-energy model describing electrons interacting via the exchange of AFM fluctuations. In its two-band version (whose physics has been argued to be similar to the one-band version  ), the model is described by the following action S = S ψ + S ϕ + S λ defined on a two-dimensional square lattice: S ψ = ∫ τ , r r ′ ∑ i = c , d [ ( ∂ τ - μ ) δ r r ′ - t i , r r ′ ] ψ i , r α † ψ i , r ′ α , S ϕ = 1 2 ∫ τ , r ( 1 v s 2 ( ∂ τ ϕ ) 2 + ( ∇ ϕ ) 2 + r 0 ϕ 2 + u 2 ( ϕ 2 ) 2 ) , S λ = λ ∫ τ , r e i Q AFM · r ϕ · ( ψ c , r α † σ α β ψ d , r β + H . c . ) . (1) Here, ∫ τ , r is shorthand for ∫ d τ ∑ r , τ ∈ [ 0 , β ) is the imaginary time, and β = 1 / T is the inverse temperature. The action S ψ describes the fermionic degrees of freedom, with the operator ψ i , r α annihilating an electron of spin α at site r and band i . Summation over α , β is implied. The two electron bands are labeled c and d . The band dispersion is parametrized by the chemical potential μ and the hopping amplitudes t i , r r ′ . Here, we consider only nearest-neighbor hopping and set t c , x = t d , y ≡ t x and t c , y = t d , x ≡ t y to enforce the system to remain invariant under a 90° rotation followed by a c ↔ d exchange. The action S ϕ describes the spin degrees of freedom, with the bosonic field ϕ denoting the AFM order parameter with the ordering wave vector Q AFM = ( π , π ) , and σ denoting Pauli matrices. The lattice derivative ( ∇ ϕ ) 2 is short for ∑ δ = x ^ , y ^ ( ϕ r - ϕ r + δ ) 2 . The parameter r 0 tunes the AFM transition, whereas v s and u describe the stiffness of the temporal and amplitude fluctuations, respectively. To save computational time, we follow previous works and consider easy-plane antiferromagnetism; i.e., ϕ = ( ϕ x , ϕ y ) [29,35,36] . The action S λ couples spins and fermions via the parameter λ . The two-band structure of the model ensures the absence of the sign problem in our simulations  . The fermionic, magnetic, and superconducting properties of this model have been thoroughly studied recently, revealing a SC dome surrounding the QCP [29,36] . In particular, the SC order parameter Δ was found to have a “ d -wave” symmetry, i.e., to change its sign between the two bands: Δ = ∫ τ , r i σ α β y ( ψ c , r α ψ c , r β - ψ d , r α ψ d , r β ) . The CO order parameter ρ investigated here also has opposite signs in the two bands (and is thus analogous to the d -wave bond CO in the one-band version of the model): ρ = ∫ τ , r e i Q CO · r σ α β 0 ( ψ c , r α † ψ c , r α - ψ d , r α † ψ d , r β ) , where Q CO is the CO wave vector. Analytical studies found a symmetry relating the SC and CO order parameters under an approximation that focuses on the vicinity of the hot spots [3–5,37] . In the two-band version of the model, each hot spot of a given pair ( c i , d i ) is located on a different band, as shown in Fig. 1(a) . The hot-spot model with linearized dispersions has an emergent SU(2) symmetry that rotates the SC order parameter Δ HS = i σ α β y ( ψ c 1 , α ψ c 2 , β - ψ d 1 , α ψ d 2 , β ) onto the CO order parameter ρ HS = σ α β 0 ( ψ c 1 , α † ψ c 2 , β - ψ d 1 , α † ψ d 2 , β ) . This CO has a diagonal wave vector Q CO ≡ ( Q 0 , Q 0 ) , which separates two hot spots belonging to different pairs but to the same band [see Fig. 1(a) ]. Our goal here is to investigate (i) to what extent does this symmetry play a role in the vicinity of an AFM QCP and (ii) more broadly, if CO is a generic feature near such a QCP. To this end, we perform a systematic investigation of the SC and CO susceptibilities in the two-band spin-fermion model. 1 10.1103/PhysRevLett.120.247002.f1 FIG. 1. (a) Schematic Fermi surface of the spin-fermion model with two bands ( c , dashed line and d , solid line). Hot spots are marked by solid symbols. Two pairs of hot spots ( c 1 , d 1 ) and ( c 2 , d 2 ) are highlighted to illustrate the relationship between the AFM wave vector Q AFM and the CO wave vector Q CO . The band dispersions used in our QMC calculations are shown in (b) (particle-hole symmetric dispersion, μ = 0 , with t y = t x / 2 ) and (c) (particle-hole asymmetric dispersion, μ / t x = - 3 , with t y = 0 ). Changing μ tunes the CO wave vector Q CO = ( Q 0 , Q 0 ) since Q 0 = 2 arccos ( - μ / 2 t x ) . We choose as our starting point the parameters for which the approximate SU(2) symmetry is promoted to an exact lattice symmetry. This corresponds to the case in which the c and d bands are particle-hole symmetric; i.e., μ = 0 . This allows us to systematically study the effect of breaking the particle-hole symmetry at the lattice level. For μ = 0 , the electronic action for a given AFM field configuration—corresponding to the S ψ and S λ terms of the action in Eq. (1) —is invariant under a rotation in particle-hole space, ψ i r α → e i Q AFM · r ( i σ α β y ) ψ i , r β † . This invariance can be seen by constructing a four-dimensional spinor that combines rotated and nonrotated operators at each band, Ψ i , r ≡ ( ψ i , r ↑ , ψ i , r ↓ , e i Q AFM · r ψ i , r ↓ † , - e i Q AFM · r ψ i , r ↑ † ) T . In this representation, when μ = 0 , the Hamiltonian commutes with all SU(2) generators τ in the particle-hole space. Importantly, the SC and CO order parameters form a three-component vector Φ ≡ ( Re Δ , Im Δ , ρ ) and couples to the electrons as ∑ r e i Q AFM · r Φ · ( σ 0 ⊗ τ ) ( Ψ c , r † Ψ c , r - Ψ d , r † Ψ d , r ) . Note that Q CO = Q AFM , enforcing ρ to be real. As a result, an enhancement of the SC susceptibility implies an equally strong enhancement in the CO channel. This symmetry is analogous to that observed in the half-filled negative- U Hubbard model  . Here, however, both the SC and CO have a d -wave symmetry. To demonstrate the existence of this SU(2) symmetry for μ = 0 , we perform QMC simulations on a square lattice of size L = 12 . All energies are expressed in terms of the hopping t x ≡ t , and the parameters are set to v s = 2 t , u = t - 1 , λ 2 = 4 t , and t y = t / 2 , resulting in the Fermi surface shown in Fig. 1(b) (other particle-hole symmetric dispersions are presented in the Supplemental Material  ). Figure 2(a) shows the SC susceptibility χ SC ; the CO susceptibility χ CO diag with diagonal wave vector Q CO = ( Q 0 , Q 0 ) , where Q 0 = π ; and the CO susceptibility χ CO axial with axial wave vector Q CO = ( Q 0 , 0 ) / ( 0 , Q 0 ) as a function of the distance to the AFM QCP for β t = 12 . The position r c of the AFM QCP is determined via the AFM susceptibility  . The degeneracy between diagonal CO and SC is evident, as well as the enhancement of both susceptibilities at the AFM QCP. The fact that χ SC = 2 χ CO diag is because the complex SC order parameter has two components whereas the real CO order parameter has one. In contrast, the axial CO susceptibility remains small and nearly unaffected by the proximity to the QCP. Figure 2(b) , which shows the behavior at the QCP, confirms that the degeneracy is present at all temperatures. 2 10.1103/PhysRevLett.120.247002.f2 FIG. 2. SC susceptibility χ SC (circles) and CO susceptibilities for diagonal wave vector Q CO = ( Q 0 , Q 0 ) , χ CO diag (triangles) and axial wave vector Q CO = ( Q 0 , 0 ) / ( 0 , Q 0 ) , χ CO axial (inverted triangles), as a function of (a) the distance r 0 - r c to the AFM QCP (fixed temperature β t = 12 ) and (b) temperature T / t (fixed r 0 = r c at the AFM QCP). The particle-hole symmetric dispersion used here is that of Fig. 1(b) . We now proceed to investigate whether there is a remnant near-degeneracy between SC and CO when the particle-hole symmetry is broken ( μ ≠ 0 ). In this case, although there is no lattice SU(2) symmetry, an approximate SU(2) symmetry of the low-energy theory near the hot spots is preserved [37,40] . To favor the CO state, we consider one-dimensional dispersions ( t y = 0 ), as shown in Fig. 1(c) , although the results are similar for finite t y (see the Supplemental Material  ). To be able to assess the relevant CO wave vectors Q CO = ( Q 0 , Q 0 ) in the finite-size QMC simulations, we choose μ values that yield commensurate Q 0 ≡ 2 arccos ( - μ / 2 t ) = ( 2 π n / L ) , with n = 2 , 3, 4, 5. Figure 3(a) displays the behavior of χ SC and χ CO diag , normalized by their noninteracting ( λ = 0 ) value, as a function of the distance to the QCP for different values of μ . While the sharp enhancement of χ SC at r 0 = r c is preserved, the enhancement of χ CO diag is small for r 0 > r c . This enhancement of χ CO diag is larger the closer μ is to zero, i.e., the closer the global lattice SU(2) symmetry is to be restored. The CO-SC degeneracy observed for μ = 0 is absent, with SC clearly winning over CO. The competition between the two orders is highlighted in Fig. 3(b) , where the T dependences of 1 / χ SC and 1 / χ CO diag are plotted at the QCP, r 0 = r c . Interestingly, right above the Berezinskii-Kosterlitz-Thouless superconducting transition temperature T c (extracted from the superfluid density; see Ref.  ), χ CO diag reverses its trend and starts to decrease upon lowering the temperature. This provides evidence that the competition between SC and CO takes place already in the fluctuating regime. 3 10.1103/PhysRevLett.120.247002.f3 FIG. 3. (a) SC (circles) and diagonal CO (triangles) susceptibilities, normalized by their noninteracting values, as a function of the distance to the QCP r 0 - r c and for a fixed temperature β t = 10 . The dispersion is represented in Fig. 1(c) , with different chemical potentials parametrized by Q 0 = 2 arccos ( - μ / 2 t ) . Panel (b) shows the temperature dependence of the inverse susceptibilities at the AFM QCP ( r 0 = r c ) for μ / t = - 2 ( Q 0 = π / 2 ). Another important result from our QMC simulations is that, when the AFM hot spots are near the antinodal region of the Brillouin zone [i.e., ( π , 0 ) / ( 0 , π ) ], the CO wave vector tends to shift from diagonal to axial inside the AFM phase. To illustrate this, in Fig. 4 we plot χ CO ( q ) for the system with μ / t = - 2 ( Q 0 = π / 2 ). Results for other fillings are discussed in the Supplemental Material  . The tendency of shifting Q CO from diagonal (paramagnetic phase) to axial (AFM phase) is evident. To quantify this behavior, we plot in Fig. 4(d) the ratio between the maxima of χ CO along the diagonal and axial directions as a function of r 0 for β t = 14 . The maximum is along the diagonal direction in the paramagnetic phase. Inside the AFM phase, however, the maximum quickly shifts to the axial direction. This effect is consistent with theoretical proposals that axial CO is favored over diagonal CO if the antinodal regions of the Brillouin zone are gapped either by AFM order discussed here or by a more exotic pseudogap state [33,34] . 4 10.1103/PhysRevLett.120.247002.f4 FIG. 4. Panels (a)–(c) show the momentum dependence of the CO susceptibility χ CO ( q ) in the paramagnetic phase r 0 = 2.2 r c (a), at the QCP (b), and in the AFM phase r 0 = 0.44 r c (c). The dispersion is represented in Fig. 1(c) with μ / t = - 2 . (d) Ratio η of the maximum values of χ CO ( q ) along the diagonal direction q = ( q , q ) and along the axial direction q = ( q , 0 ) for different inverse temperatures β (inset), as a function of the distance to the QCP at r 0 = r c . The dashed line at η = 1 marks the location where the maximum of the CO susceptibility changes between the diagonal and axial directions. The results are obtained for β t = 14 and L = 16 . In summary, we showed that the spin-fermion model with particle-hole symmetric bands has an exact SU(2) symmetry that relates d -wave SC and d -wave CO. The breaking of this symmetry at the lattice level strongly suppresses the CO susceptibility, even though it remains a low-energy symmetry near the hot spots. Compared with previous QMC investigations of the spin-fermion model, which showed that the SC instability is governed by the hot spots  , our results indicate that the CO instability is instead governed by the full electronic dispersion. Such an asymmetry between CO and SC implies that CO is not a universal phenomenon associated with AFM quantum criticality, in contrast to SC. The fragility of the SC-CO symmetry can be understood in renormalization group language. At the fixed point with a high emergent symmetry, small symmetry-breaking perturbations are strongly relevant  and they quickly drive the system away from the high-symmetry regime. The high symmetry may yet emerge at very low energies in an ideal, weak-coupling critical point. At a finite value of the coupling constant, however, interactions between degrees of freedom at different energy scales leads to a loss of the emergent symmetry. The SC channel wins because it is enhanced by both high- and low-energy degrees of freedom, in contrast to density-wave-like instabilities. The applicability of these results to specific materials—and particularly the cuprates—remains an open question. On the one hand, the CO observed in most cuprates only acquires a substantially long correlation length once SC is fully suppressed, and CO fluctuations are found to be suppressed by the onset of SC [8,18,19] . Furthermore, in the pseudogap state in which CO is experimentally observed, the CO wave vector is axial, not diagonal. All of these observations seem at least qualitatively consistent with our results for systems without particle-hole symmetric band dispersions. On the other hand, in hole-doped cuprates, AFM fluctuations become weaker as the system approaches optimal doping and CO is observed. The fact that CO is strongest near a specific doping close to 1 / 8 , where AFM fluctuations are not particularly enhanced, suggests that lattice commensuration effects, rather than AFM criticality, may play an important role in these systems. We thank A. Chubukov for fruitful discussions. X. W. and R. M. F. were supported by the U.S. Department of Energy, Office of Science, Basic Energy Sciences, under Award No. DE-SC0012336. R. M. F. also acknowledges partial support from the Research Corporation for Science Advancement via the Cottrell Scholar Award, and X. W. acknowledges support from the Doctoral Dissertation Fellowship offered by the University of Minnesota. Y. W. is supported by the Gordon and Betty Moore Foundation’s EPiQS Initiative through Grant No. GBMF4305 at the University of Illinois. Y. S. is supported by the U.S. Department of Energy, Office of Science, Basic Energy Sciences, Materials Sciences and Engineering Division, under Contract No. DE-AC02-76SF00515. R. M. F. and X. W. thank the Minnesota Supercomputing Institute (MSI) at the University of Minnesota, where part of the numerical computations was performed.  1 D. J. Scalapino , Rev. Mod. Phys. 84 , 1383 ( 2012 ). RMPHAT 0034-6861 10.1103/RevModPhys.84.1383  2 M. A. Metlitski and S. Sachdev , New J. Phys. 12 , 105007 ( 2010 ). NJOPFM 1367-2630 10.1088/1367-2630/12/10/105007  3 M. A. Metlitski and S. Sachdev , Phys. Rev. B 82 , 075128 ( 2010 ). PRBMDO 1098-0121 10.1103/PhysRevB.82.075128  4 K. B. Efetov , H. Meier , and C. Pepin , Nat. Phys. 9 , 442 ( 2013 ). NPAHAX 1745-2473 10.1038/nphys2641  5 Y. Wang and A. Chubukov , Phys. Rev. B 90 , 035149 ( 2014 ). PRBMDO 1098-0121 10.1103/PhysRevB.90.035149  6 A. Allais , J. Bauer , and S. Sachdev , Phys. Rev. B 90 , 155114 ( 2014 ). PRBMDO 1098-0121 10.1103/PhysRevB.90.155114  7 T. Wu , H. Mayaffre , S. Kramer , M. Horvatic , C. Berthier , W. N. Hardy , R. Liang , D. A. Bonn , and M.-H. Julien , Nature (London) 477 , 191 ( 2011 ). NATUAS 0028-0836 10.1038/nature10345  8 J. Chang , E. Blackburn , A. T. Holmes , N. B. Christensen , J. Larsen , J. Mesot , R. Liang , D. A. Bonn , W. N. Hardy , A. Watenphul , M. v. Zimmermann , E. M. Forgan , and S. M. Hayden , Nat. Phys. 8 , 871 ( 2012 ). NPAHAX 1745-2473 10.1038/nphys2456  9 A. J. Achkar , R. Sutarto , X. Mao , F. He , A. Frano , S. Blanco-Canosa , M. Le Tacon , G. Ghiringhelli , L. Braicovich , M. Minola , M. Moretti Sala , C. Mazzoli , R. Liang , D. A. Bonn , W. N. Hardy , B. Keimer , G. A. Sawatzky , and D. G. Hawthorn , Phys. Rev. Lett. 109 , 167001 ( 2012 ). PRLTAO 0031-9007 10.1103/PhysRevLett.109.167001  10 G. Ghiringhelli , M. Le Tacon , M. Minola , S. Blanco-Canosa , C. Mazzoli , N. B. Brookes , G. M. De Luca , A. Frano , D. G. Hawthorn , F. He , T. Loew , M. M. Sala , D. C. Peets , M. Salluzzo , E. Schierle , R. Sutarto , G. A. Sawatzky , E. Weschke , B. Keimer , and L. Braicovich , Science 337 , 821 ( 2012 ). SCIEAS 0036-8075 10.1126/science.1223532  11 E. Blackburn , J. Chang , M. Hücker , A. T. Holmes , N. B. Christensen , R. Liang , D. A. Bonn , W. N. Hardy , U. Rütt , O. Gutowski , M. v. Zimmermann , E. M. Forgan , and S. M. Hayden , Phys. Rev. Lett. 110 , 137004 ( 2013 ). PRLTAO 0031-9007 10.1103/PhysRevLett.110.137004  12 N. Doiron-Leyraud , S. Lepault , O. Cyr-Choinière , B. Vignolle , G. Grissonnanche , F. Laliberté , J. Chang , N. Barišić , M. K. Chan , L. Ji , X. Zhao , Y. Li , M. Greven , C. Proust , and L. Taillefer , Phys. Rev. X 3 , 021019 ( 2013 ). PRXHAE 2160-3308 10.1103/PhysRevX.3.021019  13 D. LeBoeuf , S. Kramer , W. N. Hardy , R. Liang , D. A. Bonn , and C. Proust , Nat. Phys. 9 , 79 ( 2013 ). NPAHAX 1745-2473 10.1038/nphys2502  14 R. Comin , A. Frano , M. M. Yee , Y. Yoshida , H. Eisaki , E. Schierle , E. Weschke , R. Sutarto , F. He , A. Soumyanarayanan , Y. He , M. Le Tacon , I. S. Elfimov , J. E. Hoffman , G. A. Sawatzky , B. Keimer , and A. Damascelli , Science 343 , 390 ( 2014 ). SCIEAS 0036-8075 10.1126/science.1242996  15 K. Fujita , C. K. Kim , I. Lee , J. Lee , M. H. Hamidian , I. A. Firmo , S. Mukhopadhyay , H. Eisaki , S. Uchida , M. J. Lawler , E.-A. Kim , and J. C. Davis , Science 344 , 612 ( 2014 ). SCIEAS 0036-8075 10.1126/science.1248783  16 E. H. da Silva Neto , P. Aynajian , A. Frano , R. Comin , E. Schierle , E. Weschke , A. Gyenis , J. Wen , J. Schneeloch , Z. Xu , S. Ono , G. Gu , M. Le Tacon , and A. Yazdani , Science 343 , 393 ( 2014 ). SCIEAS 0036-8075 10.1126/science.1243479  17 A. Mesaros , K. Fujita , S. D. Edkins , M. H. Hamidian , H. Eisaki , S.-I. Uchida , J. C. S. Davis , M. J. Lawler , and E.-A. Kim , Proc. Natl. Acad. Sci. U.S.A. 113 , 12661 ( 2016 ). PNASA6 0027-8424 10.1073/pnas.1614247113  18 H. Jang , W.-S. Lee , H. Nojiri , S. Matsuzawa , H. Yasumura , L. Nie , A. V. Maharaj , S. Gerber , Y.-J. Liu , A. Mehta , D. A. Bonn , R. Liang , W. N. Hardy , C. A. Burns , Z. Islam , S. Song , J. Hastings , T. P. Devereaux , Z.-X. Shen , S. A. Kivelson , C.-C. Kao , D. Zhu , and J.-S. Lee , Proc. Natl. Acad. Sci. U.S.A. 113 , 14645 ( 2016 ). PNASA6 0027-8424 10.1073/pnas.1612849113  19 J. Chang , E. Blackburn , O. Ivashko , A. T. Holmes , N. B. Christensen , M. Hücker , R. Liang , D. A. Bonn , W. N. Hardy , U. Rütt , M. v. Zimmermann , E. M. Forgan , and S. M. Hayden , Nat. Commun. 7 , 11494 ( 2016 ). NCAOBW 2041-1723 10.1038/ncomms11494  20 S. Bulut , W. A. Atkinson , and A. P. Kampf , Phys. Rev. B 88 , 155132 ( 2013 ). PRBMDO 1098-0121 10.1103/PhysRevB.88.155132  21 S. Sachdev and R. La Placa , Phys. Rev. Lett. 111 , 027202 ( 2013 ). PRLTAO 0031-9007 10.1103/PhysRevLett.111.027202  22 L. E. Hayward , D. G. Hawthorn , R. G. Melko , and S. Sachdev , Science 343 , 1336 ( 2014 ). SCIEAS 0036-8075 10.1126/science.1246310  23 S. Caprara , C. Di Castro , G. Seibold , and M. Grilli , Phys. Rev. B 95 , 224511 ( 2017 ). PRBMDO 2469-9950 10.1103/PhysRevB.95.224511  24 K. Kawaguchi , Y. Yamakawa , M. Tsuchiizu , and H. Kontani , J. Phys. Soc. Jpn. 86 , 063707 ( 2017 ). JUPSAU 0031-9015 10.7566/JPSJ.86.063707  25 M. Tsuchiizu , K. Kawaguchi , Y. Yamakawa , and H. Kontani , Phys. Rev. B 97 , 165131 ( 2018 ). PRBMDO 2469-9950 10.1103/PhysRevB.97.165131  26 W. Atkinson , S. Ufkes , and A. P. Kampf , Phys. Rev. B 97 , 125147 ( 2018 ). PRBMDO 2469-9950 10.1103/PhysRevB.97.125147  27 E. Berg , M. A. Metlitski , and S. Sachdev , Science 338 , 1606 ( 2012 ). SCIEAS 0036-8075 10.1126/science.1227769  28 A. Abanov , A. V. Chubukov , and J. Schmalian , Adv. Phys. 52 , 119 ( 2003 ). ADPHAH 0001-8732 10.1080/0001873021000057123  29 Y. Schattner , M. H. Gerlach , S. Trebst , and E. Berg , Phys. Rev. Lett. 117 , 097002 ( 2016 ). PRLTAO 0031-9007 10.1103/PhysRevLett.117.097002  30 H. Meier , C. Pépin , M. Einenkel , and K. B. Efetov , Phys. Rev. B 89 , 195115 ( 2014 ). PRBMDO 1098-0121 10.1103/PhysRevB.89.195115  31 D. Chowdhury and S. Sachdev , Phys. Rev. B 90 , 134516 ( 2014 ). PRBMDO 1098-0121 10.1103/PhysRevB.90.134516  32 Z.-X. Li , F. Wang , H. Yao , and D.-H. Lee , Phys. Rev. B 95 , 214505 ( 2017 ). PRBMDO 2469-9950 10.1103/PhysRevB.95.214505  33 D. Chowdhury and S. Sachdev , Phys. Rev. B 90 , 245136 ( 2014 ). PRBMDO 1098-0121 10.1103/PhysRevB.90.245136  34 W. A. Atkinson , A. P. Kampf , and S. Bulut , New J. Phys. 17 , 013025 ( 2015 ). NJOPFM 1367-2630 10.1088/1367-2630/17/1/013025  35 X. Wang , Y. Schattner , E. Berg , and R. M. Fernandes , Phys. Rev. B 95 , 174520 ( 2017 ). PRBMDO 2469-9950 10.1103/PhysRevB.95.174520  36 M. H. Gerlach , Y. Schattner , E. Berg , and S. Trebst , Phys. Rev. B 95 , 035124 ( 2017 ). PRBMDO 2469-9950 10.1103/PhysRevB.95.035124  37 Y. Wang , D. F. Agterberg , and A. Chubukov , Phys. Rev. B 91 , 115103 ( 2015 ). PRBMDO 1098-0121 10.1103/PhysRevB.91.115103  38 A. Moreo and D. J. Scalapino , Phys. Rev. Lett. 66 , 946 ( 1991 ). PRLTAO 0031-9007 10.1103/PhysRevLett.66.946  39 See Supplemental Material at http://link.aps.org/supplemental/10.1103/PhysRevLett.120.247002 for additional numerical studies of band parameters, evolution of the charge ordering wave vector, and finite size effects.  Technically, the hot-spot symmetry is SU ( 2 ) × SU ( 2 ) ∼ SO ( 4 ) and not SU(2) because for μ ≠ 0 the CO order parameter is complex, giving rise to a four-component supervector.  41 E. Fradkin , S. A. Kivelson , and J. M. Tranquada , Rev. Mod. Phys. 87 , 457 ( 2015 ). RMPHAT 0034-6861 10.1103/RevModPhys.87.457
© 2018 American Physical Society.