key: cord-0901870-4xy41kz3
authors: Li, H. X.; Zhang, T. T.; Pai, Y.-Y.; Marvinney, C.; Said, A.; Yilmaz, T.; Yin, Q.; Gong, C.; Tu, Z.; Vescovo, E.; Moore, R. G.; Murakami, S.; Lei, H. C.; Lee, H. N.; Lawrie, B.; Miao, H.
title: Observation of Unconventional Charge Density Wave without Acoustic Phonon Anomaly in Kagome Superconductors AV3Sb5 (A=Rb,Cs)
date: 2021-03-17
journal: Physical Review X
DOI: 10.1103/physrevx.11.031050
sha: 57220009ed702d6e7e8602e857a6efa157cf5b87
doc_id: 901870
cord_uid: 4xy41kz3

The combination of non-trivial band topology and symmetry breaking phases gives rise to novel quantum states and phenomena such as topological superconductivity, quantum anomalous Hall effect and axion electrodynamics. Evidence of intertwined charge density wave (CDW) and superconducting order parameters has recently been observed in a novel kagome material AV3Sb5 (A=K,Rb,Cs) that features a Z2 topological invariant in the electronic structure. However, the origin of the CDW and its intricate interplay with topological state has yet to be determined. Here, using hard x-ray scattering, we demonstrate a three-dimensional (3D) CDW with 2*2*2 superstructure in (Rb,Cs)V3Sb5. Unexpectedly, we find that the CDW fails to induce acoustic phonon anomalies at the CDW wavevector but yields a novel Raman mode which quickly damps into a broad continuum below the CDW transition temperature. Our observations exclude strong electron-phonon coupling driven CDW in AV3Sb5 and point to an unconventional and electronic-driven mechanism that couples the CDW and the topological band structure.

The combination of non-trivial band topology and symmetry breaking phases gives rise to novel quantum states and phenomena such as topological superconductivity, quantum anomalous Hall effect and axion electrodynamics. Evidence of intertwined charge density wave (CDW) and superconducting order parameters has recently been observed in a novel kagome material AV3Sb5 (A=K,Rb,Cs) that features a Z2 topological invariant in the electronic structure. However, the origin of the CDW and its intricate interplay with topological state has yet to be determined. Here, using hard x-ray scattering, we demonstrate a three-dimensional (3D) CDW with 2 × 2 × 2 superstructure in (Rb,Cs)V3Sb5. Unexpectedly, we find that the CDW fails to induce acoustic phonon anomalies at the CDW wavevector but yields a novel Raman mode that quickly damps into a broad continuum below the CDW transition temperature. Our observations exclude strong electron-phonon coupling driven CDW in AV3Sb5 and point to an unconventional particle-hole condensation mechanism that couples CDW, superconductivity and topological band structure.

The charge density wave (CDW), a translational symmetry-breaking electronic fluid, plays a crucial role in unconventional superconductors and intertwined electronic orders [1] [2] [3] [4] [5] . While CDWs have been isolated from topological excitations, recently an experimental evidence of topological CDW with chiral flux is observed in a new kagome metal, AV 3 Sb 5 (A=K, Rb, Cs) [6] , whose crystal structure and 3D Brillouin zone are shown in Fig. 1(a) and (b), respectively. Density functional theory (DFT) calculations of the electronic structure find a Z 2 topological invariant and saddle points at the M point that may lead to a CDW instability at low temperature [6] [7] [8] [9] [10] . Interestingly, scanning tunneling microscopy (STM) studies of CsV 3 Sb 5 observed evidence of roton pair-density wave and Majorana zero mode, suggesting novel topological superconductivity intertwines with CDW [11, 12] . Despite the intimate correlations between the CDW, superconductivity (SC) and topological band structure [6] [7] [8] [9] [10] [11] [12] [13] [14] [15] [16] , the nature of the CDW, in particular the role of electron-phonon coupling (EPC), remains unresolved, hindering the proper interpretations of SC in AV 3 Sb 5 .

Theoretically, the CDW mechanisms are broadly separated into two categories: weak coupling scenarios based on Fermi-surface instabilities [17] [18] [19] and strong coupling theories derived from local electron-electron or electronphonon interactions [20] [21] [22] [23] . As we show in Fig. 1(c) and (d), due to the finite EPC, the formation of CDWs distorts the underlying lattice and typically results in acoustic phonon anomalies. For a weak coupling mechanism, the phonon softening is sharply confined near the CDW wavevector, Q CDW , and is known as Kohn anomaly [19] [illustrated in Fig. 1(e) ]. In contrast, strong coupling mechanisms usually give rise to phonon anomalies in an extended Q-range near Q CDW [24] [25] [26] [27] , as shown in Fig. 1 (f). Here we combine inelastic x-ray scattering (IXS), Raman spectroscopy and angle-resolved photoemission spectroscopy (ARPES) to uncover a novel particle-hole condensation in AV 3 Sb 5 (A=Cs, Rb). Unlike conventional CDW materials, we demonstrate that neither longitudinal (LA) nor transverse (TA) acoustic phonon modes show CDW induced phonon anomalies in AV 3 Sb 5 (A=Cs, Rb). Instead, we discover a novel Raman mode that emerges at the CDW temperature, T CDW , and is quickly damped below T CDW into a continuum that broadly peaked at 20 meV. Intriguingly, the three-dimensional (3D) 2 × 2 × 2 superstructure, determined by hard x-ray diffraction (XRD), connects particle-hole scattering channels in the electronic band structure and yields a band-inversion near the Fermi level, E F , below T CDW . Our results uncover the intricate interplay between non-trivial band topology and intertwined symmetry breaking phases and shed light on the unusual SC in CsV 3 Sb 5 .

We first determine the bulk CDW superlattice peaks using meV resolution XRD. (r.l.u.). The onset temperature at T = 103 K is consistent with the CDW transition temperature in RbV 3 Sb 5 (see Appendix A). The absence of hysteresis, as shown in Fig. 2 (c) and (d), suggests that the CDW order is a second order phase transition. The observation of a CDW peak at half integer L in Fig. 2 (b) proves that the CDW in RbV 3 Sb 5 forms a 3D 2 × 2 × 2 superstructure rather than a 2D 2 × 2 superstructure as suggested by previous STM and XRD measurements [6, 8, 11, 12, 14] . Indeed, the CDW correlation length, determined by the half-widthat-half-maximum (HWHM), is larger than 1000Å along the L-direction, establishing a long range ordered 3D-CDW. As shown in Fig. 2 (e) and (f), the same 3D-CDW superlattice structure is also observed in CsV 3 Sb 5 , suggesting a ubiquitous CDW mechanism for all AV 3 Sb 5 . Interestingly, we find that while the CDW is long-range ordered, the integrated CDW intensity that is proportional to the CDW order parameter is extremely small. Comparing with fundamental Bragg peaks, the CDW peak intensity is 3 ∼ 5 orders of magnitude weaker, demonstrating small lattice distortions.

With the 3D-CDW established, we explore the origin of the CDW in AV 3 Sb 5 using IXS and Raman spectroscopy. As we have shown in Fig. 1 , the formation of CDWs is always accompanied by acoustic phonon anomalies because of EPC. Indeed, even in the cuprates, where CDWs are short ranged and possibly dynamical, a 15% LA phonon softening has been observed [24, 25] , demonstrating acoustic phonon is a sensitive probe of CDW correlations. Recent DFT calculations of the phonon spectrum of AV 3 Sb 5 show large negative frequency of the LA near the M and L point, suggesting EPC induced CDW [13] (see supplementary materials in ref.

[28]). To test this scenario, we performed IXS measurement along (3, 0, 0) − (3.5, 0, 0) and (3, 3, 0) − (3.5, 2.5, 0) to selectively enhance LA and TA phonon modes in our experimental geometry (see [28] . Unexpectedly, as we show in Fig. 3 (a)-(e), CDWs in AV 3 Sb 5 failed to induce phonon anomalies on both TA and LA modes from T = 300 K > T CDW to 50 K < T CDW . The absence of phonon softening excludes strong EPC as the driven force of the CDW in AV 3 Sb 5 . Instead, our inelastic and elastic results indicate an electronic-driven CDW mechanism with extremely weak electron acoustic-phonon interactions.

To confirm this scenario, we performed temperature dependent Raman spectroscopy on CsV 3 Sb 5 with T CDW =94 K. The results are shown in Fig. 3 (f) and (g). According to DFT calculations (see supplementary materials in [28]), four Raman active modes at 8.6, 15.1, 16.2 and 17.4 meV, corresponding to E 1g , E 2g , A 1g and B 2g modes, respectively, are allowed for the normal state structure. At 150 K, we observe three optical phonon modes at 11, 15 and 17 meV, which we attribute to E 1g , E 2g and B 2g modes (The 11 meV peak may actually be at lower energy but appeared to be at 11 meV due to the cutoff of the longpass filter at 11 meV). Interestingly, we observe a continuum that is broadly cen- tered around 19.5 meV. As we cool down towards the T CDW , the center of this continuum gradually shifts to lower energy and eventually develops into a strong, sharp and asymmetric peak at 19 meV, demonstrating its intimate correlation with the CDW. Furthermore, 19 meV is nearly two times of the STM determined single-particle gap at E F [6, 11, 12, 14] , consistent with the two-particle process of Raman scattering. Unexpectedly, moving to lower temperature, this peak is suppressed and hardened to 20 meV at 30 K. The novel temperature dependence of the broad Raman mode strongly supports an electronic origin of the CDW. Remarkably, although Raman spectroscopy measures excitations at the Brillouin zone center, this unconventional mode shows clear softening above T CDW and hardening below T CDW , reminiscent of the CDW phason and amplitudon excitations near the Q CDW [26, 27] . A similar effect has been observed in CDW material TaSe 2 [27] , which also has a 3D 2 × 2 × 2 superstructure. We suspect that CDWs in AV 3 Sb 5 might be a novel realization of particle-hole condensation in a metallic system, where the particle-hole excitation is highly damped by itinerant electrons. In this scenario, the absence of phonon softening might be a signature of large CDW phason gap that is induced by commensurability effect [29] . Indeed, As we discuss in more details in Appendix B, the combination of CDW commensurate period and the unique electronic structure of AV 3 Sb 5 allows a large gap in the CDW phase mode. We note, below 30 K, we observe additional Raman active modes at 25.4 and 27.5 meV, demonstrating additional symmetry breaking phase at low temperature. Comparing with previous STM study [14] , this new low-temperature phase is possibly the unidirectional CDW. Interestingly, this new symmetry breaking phase correlates strongly with the CDW-related continuum, where the intensity of the continuum significantly suppressed when the new Raman modes emerge.

To disclose the particle-hole scattering channel that is connected by the 3D-CDW, in Fig. 4 , we show the lowtemperature (T = 15 K) electronic structure of RbV 3 Sb 5 determined by ARPES measurements and DFT calculations. As shown in Fig. 4(a) and (b) , the overall Fermi surface topology is similar to previous ARPES studies of KV 3 Sb 5 and CsV 3 Sb 5 [8, 34] . Photon energies of Raman spectra on CsV3Sb5 from 150 K to 3 K. Beside the sharp optical phonon peak at 11, 15 and 17 meV, there is a continuum that is broadly centered near 19 meV (The 11 meV peak may actually be at lower energy but appeared to be at 11 meV due to the cutoff of the longpass filter at 90 cm −1 ). The black dots in (f) mark the peak positions of the continuum near 19 meV. When cooling towards TCDW, the center of this continuum gradually shifts from 19.5 to 19 meV and starts to build a strong, sharp and asymmetric peak at TCDW = 94 K. Moving to lower temperature, this peak is quickly damped and hardened to 20 meV at 30 K. Below 30 K, two new phonon modes at 25.4 and 27.5 meV ared observed, which may correspond to the stripe phase observed by a recent STM study [14] . The error bars in panel (d), (e) represent one standard deviation assuming Poisson counting statistics. The error bars in panel (b) denote the 2σ returned from the fitting that extract the spectral peak position. hν = 100 and 90 eV are corresponding to k z = π and 0.15π around the M point, respectively (see supplementary material [28] for the k z map). The most important observations are shown in Fig. 4 (c) and (d), corresponding to the blue-cut in Fig. 4(a) and (b), respectively. At k z = π, we observe a shallow electron-band with a band bottom of ∼ 50 meV, while at k z = 0.15π, we observe a steep hole-like band dispersing towards the Fermi level. The shallow electron and hole band can be further revealed by the stacking momentum distribution curves (MDCs) displayed in Fig. 4 (e) and (f). Similar electronic structure is also observed in CsV 3 Sb 5 (see supplementary materials [28] ). Comparing with the DFT calculations shown in Fig. 4(g) , we find that the electron and hole bands correspond to band #49,50 at the L point [ Fig. 4(g) ] and band #47,48 at the M point, respectively. The parity distribution confirms the opposite parity between these two bands, which is reported to hold a nontrivial Z 2 topology [8] . We note that at k z = 0.15π, a weak electron-like spectral weight is observed near the Fermi level. It likely comes from the bulk projected sur-face state as suggested in a previous study [8] . Our observations establish two particle-hole channels that are relevant to the CDW: (i) the quasi-nested electron and hole pockets at the M and L point [see Fig. 1 

the saddle points at the M point [4, 5, 9, 10] . The former is likely responsible for the L-component of the 3D-CDW , while the latter may play crucial role for the in-plane component as suggested by previous studies [5, 10, 35] . Interestingly, in the CDW phase, the shallow electron and hole-like bands will fold to the Brillouin zone center and open a band gap between bands with opposite parity. This CDW induced band inversion near the Fermi level may possibly responsible for the novel SC at low temperature [11, 12, 14, 16] .

In summary, we demonstrated a novel electronic driven 3D-CDW in the Z 2 kagome superconductor AV 3 Sb 5 , where the formation of CDW failed to induce phonon anomalies near Q CDW . Our observations exclude strong EPC driven CDW in AV 3 Sb 5 and point to an unconventional particle-hole condensation mechanism that couples the CDW and the topological band structure. 

Sample: Single crystals of AV 3 Sb 5 (A = Rb, Cs) were grown from Rb ingot (purity 99.9%), Cs ingot (purity 99.9%), V powder (purity 99.9%) and Sb grains (purity 99.999%) using the self-flux method, similar to the growth of RbV 3 Sb 5 [13] . The mixture was put into an alumina crucible and sealed in a quartz ampoule under partial argon atmosphere. The sealed quartz ampoule was heated to 1273 K for 12 h and soaked there for 24 h. Then it was cooled down to 1173 K at 50 K/h and further to 923 K at a slowly rate. Finally, the ampoule was taken out from the furnace and decanted with a centrifuge to separate CsV 3 Sb 5 single crystals from the flux. The obtained crystals have a typical size of 2 × 2 × 0.02 mm 3 . CDW transition is clearly observed in transport and specific heat measurement shown in Fig. A1(a) and (b), respectively. Inelastic x-ray scattering: The experiments were conducted at beam line 30-ID-C (HERIX) at the Advanced Photon Source (APS). The highly monochromatic x-ray beam of incident energy E i = 23.7 keV (l = 0.5226 A) was focused on the sample with a beam cross section of 35×15 µm 2 (horizontal × vertical). The total energy resolution of the monochromatic x-ray beam and analyzer crystals was ∆E ∼1.5 meV (full width at half maximum). The measurements were performed in transmission geometry. Typical counting times were in the range of 30 to 360 seconds per point in the energy scans at constant momentum transfer Q. H, K, L are defined in the hexagonal structure with a=b=5.472Å, c=9.073Å at the room temperature for RbV 3 Sb 5 , and a=b=5.495Å, c=9.309Å at the room temperature for CsV 3 Sb 5 . ARPES experiment: The ARPES experiments were performed on single crystals of RbV 3 Sb 5 and CsV 3 Sb 5 . The samples were cleaved in-situ in a vaccum better than 5×10 −11 torr. The experiment was performed at beamline 21-ID-1 at National Synchrotron Light Source II, Brookhaven. The measurements were taken with synchrotron light source and a Scienta-Omicron DA30 electron analyzer. The total energy resolution of the ARPES measurement is ∼15 meV. The sample stage was maintained at low temperature (T =15 K) throughout the experiment. Raman experiment: Raman spectroscopy was per-formed in a Montana Instruments closed cycle Cryostation s100 and utilized a Hubner Photonics 532 nm diode pumped laser excitation and an Isoplane SCT-320 imaging spectrograph with a 400B-eXcelon CCD camera and a 2400 groove/mm visible-holographic grating. Semrock dichroic and long pass filters were integrated in the optics train. For all Raman spectra reported here, a 2 mW laser excitation power was used, and Raman spectra were acquired for 300 s each. The sample was initially cooled to 3 K then heated to 150 K. All of the reported Raman spectra were then acquired while monotonically recooling the sample from 150 K to 3 K. The optics train was refocused at each temperature after waiting for the sample temperature to stabilize in order to correct for changes due to thermal expansion and to ensure that the sample was in thermal equilibrium before acquiring a spectrum. DFT calculations: The electronic band structure of RbV 3 Sb 5 calculated by VASP [30] [31] [32] after a full relaxation on atomic positions with an atom pair-wise correction method (DFT-D3) [33] . The phonon band structure of RbV 3 Sb 5 is calculated using VASP within the Perdew-Burke-Ernzerhof exchange correlation based on density functional perturbation theory. An equivalent k-point mesh of 9 × 9 × 6 is used in the self-consistent calculation, and the cutoff energy for the plane-wave basis is 400 eV. Prior to the phonon spectra calculation, crystal structure is relaxed with the residual force on each atom is less than 0.001 eV/Å.

For incommensurate CDW (IC-CDW), the condensation energy, E cond is independent of the phase of CDW, φ. For commensurate CDW (C-CDW), however, the condensation energy becomes phase dependent:

where λ is the dimensionless EPC constant, n( F ) and ∆ are the density-of-state at E F and the CDW gap, respectively. M is the CDW period relative to the lattice constant a 0 . D is the electronic band width. This phase dependent E cond (φ) means that gliding the CDW requires finite energy, corresponding to a finite gap in the CDW phase mode:

where ω 2 F = 4M 2 λ m∆ 2 m * ( ∆ D ) M −2 . m * is the effective mass. At q = q CDW , Eq. B2 gives:

Usually D ∆ and therefore, ω φ ∼ 0. However, for M = 2, which is the case of AV 3 Sb 5 , ω φ = 4M λ m m * ∆.

Most interestingly, a simple parabolic fitting of the ARPES data shown in Fig. 4 gives m m * ∼3 for the shallow electron and hole pockets, suggesting large ω φ in AV 3 Sb 5 . As we shown in Fig. A2c , due to the large phason gap, the particle-hole condensation will not interacting with acoustic phonon modes, as these modes are laying inside the phason gap.

Colloquium : Theory of intertwined orders in high temperature superconductors

Amperean pairing and the pseudogap phase of cuprate superconductors

Formation of incommensurate charge density waves in cuprates

Unconventional fermi surface instabilities in the kagome hubbard model

Competing electronic orders on kagome lattices at van hove filling

Discovery of topological charge order in kagome superconductor kv3sb5 (2020)

Toberer, New kagome prototype materials: discovery of kv3sb5, rbv3sb5, and csv3sb5

Csv3sb5: A 2 topological kagome metal with a superconducting ground state

Charge density waves and electronic properties of superconducting kagome metals

Chiral flux phase in the kagome superconductor av3sb5 (2021)

Three-dimensional charge density wave and robust zero-bias conductance peak inside the superconducting vortex core of a kagome superconductor csv3sb5

Roton pair density wave and unconventional strong-coupling superconductivity in a topological kagome metal (2021)

Superconductivity and normal-state properties of kagome metal rbv3sb5 single crystals

Wilson, and I. Zeljkovic, Cascade of correlated electron states in a kagome superconductor csv3sb5

Concurrence of anomalous hall effect and charge density wave in a superconducting topological kagome metal (2021)

Double superconducting dome and triple enhancement of tc in the kagome superconductor csv3sb5 under high pressure

Fermi surface nesting and the origin of charge density waves in metals

New mechanism for a chargedensity-wave instability

Giant kohn anomaly and the phase transition in charge density wave zrte3

Charged magnetic domain lines and the magnetism of high-Tc oxides

Magnetism in la2cuo4 based compounds

Charged solitons in the hartree-fock approximation to the large-u hubbard model

Strong-coupling theory of charge-density-wave transitions

Incommensurate phonon anomaly and the nature of charge density waves in cuprates

Inelastic X-ray scattering in YBa2Cu3O6.6 reveals giant phonon anomalies and elastic central peak due to charge-density-wave formation

(a) schematically shows commensurate (yellow) and incommensurate (cyan) CDW periods with respect to the lattice (red). (b) Yellow and cyan curves are the calculated phase modes for C-CDW and IC-CDW, respectively. (c) shows the CDW mechanism without electron-acoustic phonon interaction. The temperature dependent particle-hole excitations are shown as red, yellow and blue curves. The dashed line show LA mode, whose energy at qCDW is less than ω φ . Shaded pink area, indicates the momentum position of Raman measurement. of the charge-density wave in 2h−nbse2

Signatures of exciton condensation in a transition metal dichalcogenide

Conductivity from charge or spin density waves, Solid State Communications

Ab initio molecular-dynamics simulation of the liquid-metal-amorphous-semiconductor transition in germanium

Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set

Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set

Effect of the damping function in dispersion corrected density functional theory

Giant, unconventional anomalous hall effect in the metallic frustrated magnet candidate, kv3sb5

Comparison of charge modulations in La1.875Ba0.125CuO4 and YBa2Cu3O6