Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

# A new class of bilayer kagome lattice compounds with Dirac nodal lines and pressure-induced superconductivity

## Abstract

Kagome lattice composed of transition-metal ions provides a great opportunity to explore the intertwining between geometry, electronic orders and band topology. The discovery of multiple competing orders that connect intimately with the underlying topological band structure in nonmagnetic kagome metals AV3Sb5 (A = K, Rb, Cs) further pushes this topic to the quantum frontier. Here we report a new class of vanadium-based compounds with kagome bilayers, namely AV6Sb6 (A = K, Rb, Cs) and V6Sb4, which, together with AV3Sb5, compose a series of kagome compounds with a generic chemical formula (Am-1Sb2m)(V3Sb)n (m = 1, 2; n = 1, 2). Theoretical calculations combined with angle-resolved photoemission measurements reveal that these compounds feature Dirac nodal lines in close vicinity to the Fermi level. Pressure-induced superconductivity in AV6Sb6 further suggests promising emergent phenomena in these materials. The establishment of a new family of layered kagome materials paves the way for designer of fascinating kagome systems with diverse topological nontrivialities and collective ground states.

## Introduction

The lattice geometry and crystalline symmetry are key factors determining the electronic properties of a crystal. One of the most prominent examples is the kagome lattice, a two-dimensional (2D) corner-sharing triangular network, which recently emerges as a rich frontier for exploring topological and correlated electronic phenomena1,2,3,4,5,6,7,8,9,10,11. Owing to its unique lattice geometry, the kagome lattice naturally incorporates linear band crossings hosting Dirac fermions4,5 as well as destructive interference-derived flat bands in its electronic structure. During the past decade, the topological aspects of kagome lattice have been intensively studied in transition-metal kagome magnets11,12,13,14,15,16,17,18; various topologically nontrivial electronic ground states associated with different forms of magnetism have been realized, such as (massive) Dirac fermions and flat bands in ferromagnetic Fe3Sn212,13 and antiferromagnetic FeSn14, Chern-gapped Dirac fermions in ferromagnetic TbMn6Sn615, Weyl fermions in the ferromagnet Co3Sn2S216,17 and the noncollinear antiferromagnet Mn3Sn18, and so on.

More interestingly, in the absence of magnetism, electron correlations could also provoke the emergence of unusual electronic states in kagome lattices7,8,9,10. A notable example is the recently discovered topological kagome metals AV3Sb5 (A = K, Rb, Cs)19,20. This family of materials not only carries a nontrivial Z2 topological index with band inversion20, but also hosts cascade of symmetry-breaking electronic orders including (potentially chiral) charge density wave (CDW)21, nematic/sematic order22 and superconductivity20. Subsequent studies have revealed involved intertwining between these electronic orders that gives rise to numerous exotic phenomena23, including intrinsic anomalous Hall effect24,25, unusual competition between CDW and superconductivity26,27, pair density wave order28, and possible Majorana zero modes inside the superconducting vortex core29. It has been suggested that AV3Sb5 resembles the high-Tc superconductors in view of the enriched low-temperature (T) orderings23. Whilst the ongoing investigations provide glimpse to the rich interplay between topology and correlations in kagome lattice, the experimental realization has been still limited, in great part owing to the rarity of kagome materials. Here, we report a new class of vanadium compounds which contains kagome bilayers in their crystal structures. The two members in this class, namely AV6Sb6 (A = K, Rb, Cs) and V6Sb4, host different types of Dirac nodal lines in their band structures. Moreover, superconductivity is realized in the former under pressure. These results identify the bilayer vanadium-based kagome compounds as ideal candidates for studying the topological nontriviality and its fingerprints on the electronic properties in kagome lattice.

## Results and discussion

### Crystal structures and physical properties for the bilayer kagome compounds

In Fig. 1a, b we illustrate the comparison of the crystal structures for the single-layer kagome compounds AV3Sb5 (A = K, Rb, Cs) and the bilayer kagome compounds: AV6Sb6 (A = K, Rb, Cs) and V6Sb4. All these materials consist of layered structural units stacking along the crystallographic c axis. A unit cell of AV3Sb5 hosts a single V3Sb layer in which the vanadium sublattice forms a perfect 2D kagome net and the Sb atoms locate at the center of all the kagome hexagons19 (upper panel in Fig. 1a); such V3Sb layers are separated by graphite-like Sb2 net and a triangular sublattice of alkaline cations A (left panel in Fig. 1b). The resulted hexagonal structure corresponds to the space group P6/mmm. By contrast, the bilayer compounds AV6Sb6 and V6Sb4 crystalize in the rhombohedral space group R$$\bar{3}$$m (No. 166); in their unit cells two adjacent V3Sb kagome slabs with in-plane offsets of V sites30 form a structural unit (lower panel in Fig. 1a). In AV6Sb6 these (V3Sb)2 bilayers are sandwiched by Sb2 sheets and triangular A sublattices that are the same as those in AV3Sb5 (Fig. 1b), whilst in V6Sb4 the intercalated A cations are absent and there is a single Sb2 net between two (V3Sb)2 bilayers. We note that the structure of V6Sb4 is identical to the ferromagnetic kagome compound Fe3Sn212,30, whereas the AV6Sb6 series adopts an unprecedented crystal structure of kagome compounds (for detailed structural parameters, see Supplementary Tables 1 and 2). Although the single-layer and bilayer compounds have different crystalline symmetries, we stress that they can be represented by a generic chemical formula (Am-1Sb2m)(V3Sb)n (m = 1, 2; n = 1, 2) with the values of (m, n) for AV3Sb5, AV6Sb6 and V6Sb4 being (2, 1), (2, 2) and (1, 2), respectively. Hence, the crystal structure of the entire series (Am-1Sb2m)(V3Sb)n can be viewed as an alternate stacking of the (Am-1Sb2m) blocks and the (V3Sb)n blocks.

Single crystals of AV6Sb6 (A = K, Rb, Cs) and V6Sb4 were synthesized using a self-flux growth technique (Methods). The high quality of these single crystals is confirmed by our X-ray diffraction (XRD) measurements: the reconstructed (hk0) plane determined from single-crystal XRD displays sharp spots with sixfold rotational symmetry (Fig. 1d, e) and a 2θ-ω scan yields only a series of sharp (00l) Bragg peaks (Fig. 1c). We summarize the results of the single-crystal XRD analysis for V6Sb4, RbV6Sb6, and CsV6Sb6 in Supplementary Tables 1 and 2. Electrical transport and magnetization measurements identify all the bilayer compounds as nonmagnetic metals without indications of additional orderings (Fig. 1f and Supplementary Fig. S1). We note that no CDW ordering is observed in our bilayer kagome samples, contrary to the case in the single-layer AV3Sb5 (A = K, Rb, Cs) where a CDW transition occurs at 80–100 K19. Hall resistivity measured in CsV6Sb6 (inset of Fig. 1f) indicates dominant electron-type carriers with the density ne = 2.1 × 1021 cm−3.

### First-principles calculations on CsV6Sb6 and V6Sb4

To further investigate the electronic structures of those bilayer kagome compounds, we performed first-principles calculations based on the density functional theory (DFT). The calculation results approve that AV6Sb6 (A = K, Rb, Cs) and V6Sb4 lack local magnetic interactions and possess nonmagnetic ground states, consistent with the experimental observations (Supplementary Fig. 1). Here, we mainly focus on the compounds CsV6Sb6 and V6Sb4. The calculated band structures for the other two members in the isostructural series AV6Sb6 with A = K, Rb are presented in Supplementary Fig. 2. In Fig. 2a, b we show the band structures as well as the projected density of states (DOS) for CsV6Sb6 and V6Sb4, respectively. The corresponding high-symmetry paths of the BZ of a rhombohedral lattice are depicted in Fig. 2c. The DOS in the vicinity of the Fermi level are dominated by V-3d orbitals for both CsV6Sb6 and V6Sb4, whereas the Sb-p orbitals have rather weak contributions (Fig. 2a, b). Remarkably, the calculated band structures of both compounds exhibit linear band crossings (i.e, Dirac points) that are close to the Fermi level along the specific high-symmetry paths, Г-Z-F-Г-L-Z-P (Fig. 2a, b). To be mentioned, the electronic structure of the bilayer kagome compounds is distinct from that of the single-layer AV3Sb5 which are hallmarked by multiple Dirac crossings and saddle points near the Fermi level20,31.

Our DFT calculations reveal that CsV6Sb6 is a Dirac nodal line semimetal. As shown in Fig. 2a, the band crossings occur along the high-symmetry paths from Γ/Z to the high-symmetry points at the boundaries of the Brillouin zone (BZ) such as the P, F, and L. All these high-symmetry paths lie in a middle plane of the BZ (highlighted in yellow in Fig. 2c); on this plane, the band crossings form two type-II Dirac nodal lines (shown as the red lines in Fig. 2c) characterized by tilted Dirac cones32. By further checking the little group of the BZ, we find that these type-II Dirac nodal lines are symmetry protected because the two intersecting bands along the above arbitrary high-symmetry paths belong to different irreducible representations (IR) Γ1 and Γ2 of the mirror symmetry Cs (see the insets in Fig. 2a). Due to the three-fold rotational symmetry, there are three equivalent middle planes, and thus six type-II Dirac nodal lines that are symmetrically distributed in the BZ of CsV6Sb6. On the other hand, the crossings of valence and conduction bands are absent along Γ/Z-L in the band structure of V6Sb4 (Fig. 2b). However, the preserved band crossings along Γ/Z-F can still form symmetry-protected nodal lines on three equivalent middle planes. Contrary to the type-II nodal lines in CsV6Sb6, these nodal lines in V6Sb4 are type-I and feature closed loops around the band inverted point F (Supplementary Fig. 2).

In the presence of the spin-orbit coupling (SOC), the spin-rotation symmetry is broken, subsequently, the nodal lines in a system with the coexistence of spatial inversion and time-reversal symmetries are always destroyed33. As shown in the insets of Fig. 2a, b, when the SOC effect is included, the two intersecting bands without SOC now belong to the same IR Γ4 of the mirror symmetry Cs. Thus the band crossings are avoided. Nonetheless, the SOC effect is rather weak in the vanadium bilayer kagome compounds: the gap opened at the Dirac band crossings is almost negligible; in particular, for CsV6Sb6 the gap width is less than 1 meV. Therefore, its type-II Dirac nodal lines are nearly intact. With the existence of the SOC gaps, we can use parity products of occupied bands at the time-reversal invariant momenta (TRIM) points to reveal the topological nontriviality34. For a nonmagnetic compound crystalizing in the R$$\bar{3}$$m space group, the three F (and the three L) points are equivalent, thus we only calculate the parity eigenvalues at four TRIM points, i.e., the Γ, Z, F, and L points (for detailed calculation results, see Supplementary Table 3). As illustrated in Fig. 2d, the parity analysis indicates that bands 115 and 117 (95 and 99) are topologically nontrivial in CsV6Sb6 (V6Sb4). These nontrivial bands are close to the Fermi level and their cooperation with nodal fermions would be expected to generate rich exotic quantum phenomena. It should be noticed that the inclusion of Hubbard U in the DFT calculations does not change the band topology (Supplementary Fig. 3), and the Dirac nodal lines remain intact with U = 2.0 eV.

### ARPES measurements on CsV6Sb6

The band structure calculation results for CsV6Sb6 are supported by our angle-resolved photoemission spectroscopy (ARPES) measurements. As shown in Fig. 3a–f, the constant energy contours at binding energies Eb = 0, 200 and 400 meV, as measured with 60 eV photons, are similar to those of the DFT+ U calculations (U = 2eV). At the Fermi energy, we observe relatively high intensities close to the boundary of the projected 2D BZ (marked by yellow lines in Fig. 3b) and a hexagonal contour around the $$\bar{\Gamma }$$ point; both of which are consistent with the DFT results (Fig. 3a). At Eb = 200 meV, the most prominent features of the calculated constant energy contours are the hexagonal pocket around the $$\bar{\Gamma }$$ point and rounded-triangular pockets centered at the $$\bar{{{{{{\rm{K}}}}}}}$$ points (Fig. 3c), whereas a gear-shaped pocket centered at the $$\bar{\Gamma }$$ develops at higher binding energies (Fig. 3e). These features are well reproduced in the photoemission intensity maps (Fig. 3d, f). Figure 3g shows the calculated bulk band dispersion along the $$\bar{\Gamma }$$-$$\bar{{{{{{\rm{M}}}}}}}$$ direction (red dashed line in Fig. 3b) with different kz ranging from 0 to 2π/c, and one can find that the kz dependence of the bulk bands in this material is very weak, which is also confirmed by the photon-energy-dependent ARPES measurements (Supplementary Figs. 4, 5). For a direct comparison, in Fig. 3h we display the corresponding ARPES-intensity plot. An overall agreement between the ARPES data and the calculated bulk bands (Fig. 3g). In particular, the two crossing bands (bands 117 and 119 in Fig. 2a) contributing to the type-II Dirac nodal lines in the $$\bar{\Gamma }$$-$$\bar{{{{{{\rm{M}}}}}}}$$ direction are resolved, with the Dirac crossings located in close vicinity to the EF. We note that the band structure of CsV6Sb6 is significantly different from its single-layer counterpart CsV3Sb513,31, indicating the importance of inter-(V3Sb)-layer coupling in these kagome lattice compounds.

### Pressure-induced superconductivity in AV3Sb5

Despite that the bilayer vanadium-based kagome compounds, unlike AV3Sb5 (A = K, Rb, Cs)20,26, do not exhibit superconductivity at ambient pressure, we realize superconductivity in AV6Sb6 (A = K, Rb, Cs) by applying quasi-hydrostatic pressures. The results of high-pressure resistance measurements on two CsV6Sb6 samples are presented in Fig. 4a. The T-dependences of the normalized resistance R/R300K show metallic behavior in the entire pressure range with the residual resistivity ratio (RRR) gradually decreases below 15 GPa. At 15.8 GPa it drops to ~1.2. With pressure further increasing, superconducting transition emerges at 21.2 GPa manifested by a pronounced resistance drop (Fig. 4b). Magnetic fields can gradually suppress the transition temperature, which confirms that the resistance drop is due to a superconducting transition (Supplementary Fig. 6). The transition temperature determined by $${T}_{c}^{90 \% }$$ (i.e., where the resistance drop to 90% of the normal state value) is 1.04 K. The evolution of $${T}_{c}^{90 \% }$$ under applied pressure is nonmonotonic: it first increases rapidly and reaches a maximum of 1.48 K at 33.0 GPa, then starts to decrease slowly; superconductivity persists up to the highest pressure (79.5 GPa) achieved in this measurement, where $${T}_{c}^{90 \% }$$ is ~1 K (Fig. 4b). Such nonmonotonic behavior gives rise to a broad dome-shaped superconducting regime in the high-pressure phase diagram (Fig. 4c), suggesting a complex interplay of the pressure-dependent DOS and structural instabilities35,36. Similar superconducting dome is also observed in RbV6Sb6 and KV6Sb6 under high pressure, yet with lower maximum Tc (Supplementary Fig. 7). Intriguingly, superconductivity in all three materials appears in the vicinity of the minimum of RRR (Fig. 4c and Supplementary Fig. 7). The correlation between Tc and RRR indicates that the emergence of superconductivity is associated with electronic structure modifications.

We note that such modification stems from a structural phase transition that occurs above ~20 GPa, where the structure changes from rhombohedral to monoclinic as revealed by our high-pressure XRD measurements (Supplementary Figs. 8, 9). It is most likely that the monoclinic phase hosts superconductivity. Refinements of the XRD data provide further details of the structural transition: the ambient-pressure rhombohedral phase evolves into the monoclinic phase via a lattice distortion at which the in-plane lattice parameter a develops into two unequal values a and c, whilst the angle β changes from 120° to ~110° (Supplementary Fig. 8). The pressure-induced superconductivity accompanied by a structural transition in AV6Sb6 resembles that observed in numerous topological semimetals37,38,39,40,41. Most notably, in some of these materials, the superconductivity is proposed to emerge from topological electronic bands37,39,40, offering a good opportunity to probe the feasible realization of topological superconductivity. Future studies are needed to determine the topological properties of the high-pressure monoclinic phase and to clarify how superconductivity develops in this phase with lower symmetry.

Our theoretical and experimental investigations reveal unusual topological metal phases in the bilayer vanadium-based kagome compounds. Distinct from the typical kagome physics concerning Dirac points and van Hove singularity in single-layered AV3Sb520,23, the most important characteristics of the bilayer compounds are denoted by the Dirac nodal lines near the Fermi level. The pressure-induced superconductivity discovered in AV6Sb6 further suggests promising emergent phenomena in the bilayer kagome materials. All these results provide inspiring perspectives for future explorations of the enriched topological physics in kagome lattice.

## Methods

### Single crystal growth and characterization

Single crystals of AV6Sb6 (A = K, Rb, Cs) and V6Sb4 were grown using a self-flux method. For AV6Sb6, K/Rb/Cs (ingot or liquid 99.5%), V (powder, 99.5%), Sb (shot, 99.999%) were loaded into an alumina container with the molar ratio of 2:3:6 and then sealed into a double-wall silica tube. Mixture was subsequently heated to 1150 °C and kept for 24 h, then slowly cooled to 1050 °C for 72 h. The excess flux was removed by centrifuging at 1050 °C. The products were single crystals with sizes of a few millimeters in ab plane and less than 100 μm in thickness. For V6Sb4, K (ingot, 99.5%), V (powder, 99.5%), Sb (shot, 99.999%) were mixed with the molar ratio of 2:3:6 and sealed in the same double-wall silica tube. The tube was heated to 1150 °C and kept for 24 h, then slowly cooled to 1050 °C in 72 h. After soaked at 1050 °C for 5 h, the temperature was increased to 1100 °C in 20 h. Large crystals of V6Sb4 with sizes of several millimeters were obtained from the flux by centrifuging at 1100 °C.

Single-crystal X-ray diffraction measurements were carried out on a XtaLAB AFC12 (RINC): Kappa single diffractometer (Rigaku, Japan) with a charge-coupled device detector and Cu source in Core Facility Center for Life Sciences, USTC. The data was processed and reduced using CrysAlisPro42. Using Olex-243, the structure was solved with the ShelXT structure solution program44 via direct methods and refined with the ShelXL refinement package45. Magnetization measurements were performed on a Quantum Design Magnetic Properties Measurement Systems (MPMS-5). Plate-shaped single crystals were attached to a quartz rod with the magnetic field applied parallel to and perpendicular to the c-axis. The transport properties were measured on the Quantum Design Physical Properties Measurement System (PPMS-9) using a standard six-probe configuration.

### First-principles calculations

To depict electronic properties of these novel kagome compounds of AV6Sb6 and V6Sb4, we carried out first-principles calculations based on the density functional theory46 as implemented in the Vienna ab initio simulation package47. The exchange-correlation functional was described by generalized gradient approximation with Perdew–Burke–Ernzerhof formalism48. The core-valence interactions were treated by projector augmented-wave potentials49 with a plane-wave-basis cutoff of 450 eV. The Brillouin zone (BZ) was sampled by a 12 × 12 × 12 Monkhorst-Pack grid50 to simulate the rhombohedral structure. The crystal structures were fully relaxed by minimizing the forces on each atom smaller than 1.0 × 10−3 eV/Å, and the van der Waals interactions along the c-layer stacking direction were considered by the Crimme (DFT-D3) method51. The topological class was characterized by the Z2 invariants, which are calculated from the parity eigenvalues at TRIM points using IRVSP package52. We also employed the DFT + U method to calculate the band structure of CsV6Sb6, in which the Hubbard U correction represents the on-site Coulomb interactions on the d-orbital of vanadium; the band topology persists upon varying U values (Supplementary Information Fig. S3).

### ARPES measurements

ARPES measurements were performed at the beamline 13 U of the National Synchrotron Radiation Laboratory at University of Science and Technology of China (photon energy 35 eV), and the beamline 03U of the Shanghai Synchrotron Radiation Facility (SSRF) (photon energy 42–60 eV). The samples were cleaved in situ with a base pressure less than 6 × 10−11 torr. We note that the terraced surfaces of the cleaved samples usually hamper a clear observation of the fine electronic features. Further ARPES investigations with improved resolution are required to reveal the possible topological surface states53 in the AV6Sb6 materials.

### High-pressure transport measurements

Diamond anvils with various culets (200–300 µm) were used for high-pressure transport measurements. NaCl was used as a pressure transmitting medium and the pressure was calibrated by using the shift of ruby fluorescence and diamond anvil Raman at room temperature. For each measurement cycle, the pressure was applied at room temperature using the miniature diamond anvil cell. The transport measurements were performed in a dilution refrigerator (Kelvinox JT, Oxford Instruments) or a 3He cryostat (HelioxVT, Oxford Instruments). Single-crystalline samples of AV6Sb6 and V6Sb4 were cut into typical dimensions of 50 × 50 × 10 μm3. The resistivity was measured using a four-probe configuration. V6Sb4 does not show indication of superconducting transition up to the highest pressure we achieve, i.e., P ≈ 80 GPa (Supplementary Fig. 7).

### High-pressure X-ray diffraction

The high-pressure synchrotron XRD measurements were performed at room temperature at the beamline BL15U1 of SSRF with a wavelength of λ = 0.6199 Å. A symmetric diamond anvil cell with a pair of 200 µm culet size anvils was used to generate pressure. 70 µm sample chamber is drilled from the Re gasket and Daphne 7373 oil was loaded as a pressure transmitting medium.

## Data availability

All data supporting the findings of this study are available from the corresponding authors upon reasonable request.

## References

1. Balents, L. Spin liquids in frustrated magnets. Nature 464, 199 (2010).

2. Zhou, Y., Kanoda, K. & Ng, T.-K. Quantum spin liquid states. Rev. Mod. Phys. 89, 025003 (2017).

3. Broholm, C. et al. Quantum spin liquids. Science 367, 263 (2020).

4. Guo, H. M. & Franz, M. Topological insulator on the kagome lattice. Phys. Rev. B 80, 113102 (2009).

5. Mazin, I. I. et al. Theoretical prediction of a strongly correlated Dirac metal. Nat. Commun. 5, 4261 (2014).

6. Xu, G., Lian, B. & Zhang, S.-C. Intrinsic quantum anomalous Hall effect in the kagome lattice Cs2LiMn3F12. Phys. Rev. Lett. 115, 186802 (2015).

7. Ko, W.-H., Lee, P. A. & Wen, X.-G. Doped kagome system as exotic superconductor. Phys. Rev. B 79, 214502 (2009).

8. Yu, S.-L. & Li, J.-X. Chiral superconducting phase and chiral spin-density-wave phase in a Hubbard model on the kagome lattice. Phys. Rev. B 85, 144402 (2012).

9. Kiesel, M. L., Platt, C. & Thomale, R. Unconventional Fermi surface instabilities in the Kagome Hubbard model. Phys. Rev. Lett. 110, 126405 (2013).

10. Wang, W.-S., Li, Z.-Z., Xiang, Y.-Y. & Wang, Q.-H. Competing electronic orders on kagome lattices at van Hove filling. Phys. Rev. B 87, 115135 (2013).

11. Yin, J.-X., Pan, S. H. & Hasan, M. Z. Probing topological quantum matter with scanning tunneling microscopy. Nat. Rev. Phys. 3, 249 (2021).

12. Ye, L. et al. Massive Dirac fermions in a ferromagnetic kagome metal. Nature 555, 638 (2018).

13. Lin, Z. et al. Flatbands and emergent ferromagnetic ordering in Fe3Sn2 Kagome lattices. Phys. Rev. Lett. 121, 096401 (2018).

14. Kang, M. et al. Dirac fermions and flat bands in the ideal kagome metal FeSn. Nat. Mater. 19, 163–169 (2020).

15. Yin, J.-X., et al. Quantum limit Chern topological magnetism in TbMn6Sn6, Nature 583–533 (2020).

16. Liu, E. et al. Giant anomalous Hall effect in a ferromagnetic kagome-lattice semimetal. Nat. Phys. 14, 1125–1131 (2018).

17. Liu, D. F. et al. Magnetic Weyl semimetal phase in a Kagome crystal. Science 365, 1282 (2019).

18. Kuroda, K. et al. Evidence for magnetic Weyl fermions in a correlated metal. Nat. Mater. 16, 1090–1095 (2017).

19. Ortiz, B. R. et al. New kagome prototype materials: discovery of KV3Sb5, RbV3Sb5, and CsV3Sb5. Phys. Rev. Mater. 3, 094407 (2019).

20. Ortiz, B. R. et al. CsV3Sb5: a Z2 topological kagome metal with a superconducting ground state. Phys. Rev. Lett. 125, 247002 (2020).

21. Jiang, Y.-X. et al. Unconventional chiral charge order in kagome superconductor KV3Sb5. Nat. Mater. https://doi.org/10.1038/s41563-021-01034-y (2021).

22. Zhao, H. et al. Cascade of correlated electron states in a kagome superconductor CsV3Sb5. Nature https://doi.org/10.1038/s41586-021-03946-w (2021).

23. Jiang, K. et al. Kagome superconductors AV3Sb5 (A = K, Rb, Cs), https://arxiv.org/abs/2109.10809 (2021).

24. Yang, S.-Y. et al. Unconventional anomalous hall effect in the metallic frustrated magnet candidate, KV3Sb5. Sci. Adv. 6, eabb6003 (2020).

25. Yu, F. H. et al. Concurrence of anomalous hall effect and charge density wave in a superconducting topological Kagome metal. Phys. Rev. B 104, L041103 (2021).

26. Yu, F. H. et al. Unusual competition of superconductivity and charge-density-wave state in a compressed topological kagome metal. Nat. Commun. 12, 3645 (2021).

27. Chen, K. Y. et al. Double superconducting dome and triple enhancement of Tc in the Kagome superconductor CsV3Sb5 under high pressure. Phys. Rev. Lett. 126, 247001 (2021).

28. Chen, H. et al. Roton pair density wave in a strong-coupling kagome superconductor. Nature https://doi.org/10.1038/s41586-021-03983-5 (2021).

29. Liang, Z. et al. Three-dimensional charge density wave and surface-dependent vortex-core states in a Kagome superconductor CsV3Sb5. Phys. Rev. X 11, 031026 (2021).

30. Fenner, L. A., Dee, A. A. & Wills, A. S. Non-collinearity and spin frustration in the itinerant kagome ferromagnet Fe3Sn2. J. Phys.: Condens. Matter 21, 452202 (2009).

31. Kang, M. et al. Twofold van Hove singularity and origin of charge order in topological kagome superconductor CsV3Sb5. Nat. Phys. 18, 301–308 (2022).

32. Soluyanov, A. A. et al. Type-II Weyl semimetals. Nature 527, 495 (2015).

33. Kim, Y., Wieder, B. J., Kane, C. L. & Rappe, A. M. Dirac line nodes in inversion-symmetric crystals. Phys. Rev. Lett. 115, 036806 (2015).

34. Fu, L. & Kane, C. L. Topological insulators with inversion symmetry. Phys. Rev. B 76, 045302 (2007).

35. Pan, X.-C. et al. Pressure-driven dome-shaped superconductivity and electronic structural evolution in tungsten ditelluride. Nat. Commun. 6, 7805 (2015).

36. Kang, D. et al. Superconductivity emerging from a suppressed large magnetoresistant state in tungsten ditelluride. Nat. Commun. 6, 7804 (2015).

37. Zhou, Y. et al. Pressure-induced superconductivity in a three-dimensional topological material ZrTe5. Proc. Natl Acad. Sci. USA 113, 2904–2909 (2016).

38. Li, X. et al. Pressure-induced phase transitions and superconductivity in a quasi–1-dimensional topological crystalline insulator α-Bi4Br4. Proc. Natl Acad. Sci. USA 116, 17696–17700 (2019).

39. ElGhazali, M. A. et al. Pressure-induced superconductivity up to 13.1 K in the pyrite phase of palladium diselenide PdSe2. Phys. Rev. B 96, 060509(R) (2017).

40. Cheng, E. et al. Pressure-induced superconductivity and topological phase transitions in the topological nodal-line semimetal SrAs3. NPJ Quantum Mater. 5, 38 (2020).

41. Pei, C. et al. Pressure-induced superconductivity and structure phase transition in Pt2HgSe3. NPJ Quantum Mater. 6, 98 (2021).

42. CrysAlisPro Software system. Version 1.171.37.35. (Agilent Technologies Ltd, Yarnton, Oxfordshire, England, 2014).

43. Bourhis, L. J., Dolomanov, O. V., Gildea, R. J., Howard, J. A. K. & Puschmann, H. Olex2: a complete structure solution, refinement and anlysis program. J. Appl. Crystallogr. 42, 339–341 (2009).

44. Sheldrick, G. M. SHELX—Integrated space-group and crystal-structure determination. Acta Crystallogr. Sect. A 71, 3–8 (2015).

45. Sheldrick, G. M. A short history of SHELX. Acta Crystallogr. A 64, 112–122 (2008).

46. Kohn, W. & Sham, L. J. Self-consistent equations including exchange and correlation effects. Phys. Rev. 140, A1133 (1965).

47. Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 54, 11169 (1996).

48. Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 77, 3865 (1996).

49. Kresse, G. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. Rev. B 59, 1758 (1999).

50. Monkhorst, H. J. & Pack, J. D. Special points for Brillouin-zone integrations. Phys. Rev. B 13, 5188 (1976).

51. Grimme, S., Ehrlich, S. & Goerigk, L. Effect of the damping function in dispersion corrected density functional theory. J. Comput. Chem. 32, 1456 (2011).

52. Gao, J., Wu, Q., Persson, C. & Wang, Z. Irvsp: to obtain irreducible representations of electronic states in the VASP. Comput. Phys. Commun. 261, 107760 (2021).

53. Yang, Y. et al. Type-II nodal line fermions in New Z2 topological semimetals AV6Sb6 (A = K, Rb, and Cs) with kagome bilayer. Phys. Rev. B 104, 245128 (2021).

## Acknowledgements

We thank Shengtao Cui and Soohyun Cho for their assistance in synchrotron ARPES measurements. This work is supported by the National Key Research and Development Program of the Ministry of Science and Technology of China (2017YFA0303001, 2019YFA0704901, and 2016YFA0300201), the Anhui Initiative in Quantum Information Technologies (AHY160000), the Strategic Priority Research Program of the Chinese Academy of Sciences (XDB25000000), the National Natural Science Foundation of China (NSFC, Grants No. 11888101, 11974062, and U20322), the Science Challenge Project of China (TZ2016004) and the Key Research Program of Frontier Sciences, CAS, China (QYZDY-SSW-SLH021). The DFT calculations in this work are supported by the Supercomputing Center of University of Science and Technology of China. High-pressure synchrotron XRD work was performed at the BL15U1 beamline, SSRF in China. Part of this research used Beamline 03U of the Shanghai Synchrotron Radiation Facility, which is supported by ME2 project under contract No. 11227902 from the National Natural Science Foundation of China.

## Author information

Authors

### Contributions

X.H.C. conceived the project and supervised the overall research. M.S., F.Y., and J.Y. grew the single crystal samples. M. S. performed the XRD, electrical transport, and magnetization measurements with the help of B.L., and analyzed the data with T.W., Z.W., and Z.X. F.M., Y.L., Z.S., Z.J., Z.L., D.S., and J.H. performed the ARPES experiments and analyzed the resultant data. F. Y and J. Y. performed the high-pressure measurements. Y.Y. and R.W. performed the DFT calculations. Z.W., Z.X., and X.H.C. wrote the manuscript with input from all authors.

### Corresponding authors

Correspondence to Jianjun Ying or Xianhui Chen.

## Ethics declarations

### Competing interests

The authors declare no competing interests.

## Peer review

### Peer review information

Nature Communications thanks the anonymous reviewers for their contribution to the peer review of this work. Peer reviewer reports are available.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Rights and permissions

Reprints and Permissions

Shi, M., Yu, F., Yang, Y. et al. A new class of bilayer kagome lattice compounds with Dirac nodal lines and pressure-induced superconductivity. Nat Commun 13, 2773 (2022). https://doi.org/10.1038/s41467-022-30442-0

• Accepted:

• Published:

• DOI: https://doi.org/10.1038/s41467-022-30442-0