Abstract
Localized or propagating Majorana boundary modes are the key feature of topological superconductors. They are rare in naturally occurring compounds, but the tailored manipulation of quantum matter offers opportunities for their realization. Specifically, lattices of Yu–Shiba–Rusinov bound states—Shiba lattices—that arise when magnetic adatoms are placed on the surface of a conventional superconductor can be used to create topological bands within the superconducting gap of the substrate. Here we reveal two signatures consistent with the realization of two types of mirror-symmetry-protected topological superconductor using scanning tunnelling microscopy to create and probe adatom lattices with single-atom precision. The first has edge modes as well as higher-order corner states, and the second has symmetry-protected bulk nodal points. In principle, their topological character and boundary modes should be protected by the spatial symmetries of the adatom lattice. Our results highlight the potential of Shiba lattices as a platform to design the topology and sample geometry of two-dimensional superconductors.
This is a preview of subscription content, access via your institution
Access options
Access Nature and 54 other Nature Portfolio journals
Get Nature+, our best-value online-access subscription
$29.99 / 30 days
cancel any time
Subscribe to this journal
Receive 12 print issues and online access
$209.00 per year
only $17.42 per issue
Buy this article
- Purchase on Springer Link
- Instant access to full article PDF
Prices may be subject to local taxes which are calculated during checkout
Similar content being viewed by others
Data availability
The data underlying Figs. 1, 2 and 4 are available at https://doi.org/10.6084/m9.figshare.22794413. Source data are provided with this paper.
References
Meissner, W. & Ochsenfeld, R. Ein neuer Effekt bei Eintritt der Supraleitfähigkeit. Naturwissenschaften 21, 787–788 (1933).
Yu, L. Bound state in superconductors with paramagnetic impurities. Acta Phys. Sin. 21, 75–91 (1965).
Shiba, H. Classical spins in superconductors. Prog. Theor. Phys. 40, 435–451 (1968).
Rusinov, A. I. Superconductivity near a paramagnetic impurity. Zh. Eksp. Teor. Fiz. Pisma Red. 9, 146 (1969).
Küster, F. et al. Correlating Josephson supercurrents and Shiba states in quantum spins unconventionally coupled to superconductors. Nat. Commun. 12, 1108 (2021).
Schneider, L. et al. Magnetism and in-gap states of 3d transition metal atoms on superconducting Re. npj Quantum Mater. 4, 42 (2019).
Nadj-Perge, S., Drozdov, I. K., Bernevig, B. A. & Yazdani, A. Proposal for realizing Majorana fermions in chains of magnetic atoms on a superconductor. Phys. Rev. B 88, 020407 (2013).
Pientka, F., Peng, Y., Glazman, L. & von Oppen, F. Topological superconducting phase and Majorana bound states in Shiba chains. Phys. Scr. T164, 014008 (2015).
Brydon, P. M. R., Das Sarma, S., Hui, H.-Y. & Sau, J. D. Topological Yu-Shiba-Rusinov chain from spin-orbit coupling. Phys. Rev. B 91, 064505 (2015).
Nadj-Perge, S. et al. Observation of Majorana fermions in ferromagnetic atomic chains on a superconductor. Science 346, 602–607 (2014).
Ruby, M. et al. End states and subgap structure in proximity-coupled chains of magnetic adatoms. Phys. Rev. Lett. 115, 197204 (2015).
Ruby, M., Heinrich, B. W., Peng, Y., von Oppen, F. & Franke, K. J. Exploring a proximity-coupled Co chain on Pb(110) as a possible Majorana platform. Nano Lett. 17, 4473–4477 (2017).
Schneider, L. et al. Controlling in-gap end states by linking nonmagnetic atoms and artificially-constructed spin chains on superconductors. Nat. Commun. 11, 4707 (2020).
Mier, C. et al. Atomic manipulation of in-gap states in the β−Bi2Pd superconductor. Phys. Rev. B 104, 045406 (2021).
Schneider, L. et al. Precursors of Majorana modes and their length-dependent energy oscillations probed at both ends of atomic Shiba chains. Nat. Nanotechnol. 17, 384–389 (2022).
Kim, H. et al. Toward tailoring Majorana bound states in artificially constructed magnetic atom chains on elemental superconductors. Sci. Adv. 4, eaar5251 (2018).
Kitaev, A. Y. Unpaired Majorana fermions in quantum wires. Phys. Usp. 44, 131–136 (2001).
Li, J. et al. Two-dimensional chiral topological superconductivity in Shiba lattices. Nat. Commun. 7, 12297 (2016).
Ménard, G. C. et al. Two-dimensional topological superconductivity in Pb/Co/Si(111). Nat. Commun. 8, 2040 (2017).
Palacio-Morales, A. et al. Atomic-scale interface engineering of Majorana edge modes in a 2D magnet-superconductor hybrid system. Sci. Adv. 5, eaav6600 (2019).
Kezilebieke, S. et al. Topological superconductivity in a van der Waals heterostructure. Nature 588, 424–428 (2020).
Kitaev, A. Periodic table for topological insulators and superconductors. AIP Conf. Proc. 1134, 22–30 (2009).
Ryu, S., Schnyder, A. P., Furusaki, A. & Ludwig, A. W. W. Topological insulators and superconductors: tenfold way and dimensional hierarchy. New J. Phys. 12, 065010 (2010).
Ono, S., Po, H. C. & Watanabe, H. Refined symmetry indicators for topological superconductors in all space groups. Sci. Adv. 6, eaaz8367 (2020).
Shiozaki, K. & Sato, M. Topology of crystalline insulators and superconductors. Phys. Rev. B 90, 165114 (2014).
Wang, Q.-Z. & Liu, C.-X. Topological nonsymmorphic crystalline superconductors. Phys. Rev. B 93, 020505 (2016).
Wu, X. et al. Boundary-obstructed topological high-Tc superconductivity in iron pnictides. Phys. Rev. X 10, 041014 (2020).
Benalcazar, W. A., Bernevig, B. A. & Hughes, T. L. Quantized electric multipole insulators. Science 357, 61–66 (2017).
Schindler, F. et al. Higher-order topological insulators. Sci. Adv. 4, eaat0346 (2018).
Geier, M., Trifunovic, L., Hoskam, M. & Brouwer, P. W. Second-order topological insulators and superconductors with an order-two crystalline symmetry. Phys. Rev. B 97, 205135 (2018).
Khalaf, E. Higher-order topological insulators and superconductors protected by inversion symmetry. Phys. Rev. B 97, 205136 (2018).
Schnyder, A. P. & Brydon, P. M. R. Topological surface states in nodal superconductors. J. Phys. Condens. Matter 27, 243201 (2015).
Schneider, L. et al. Topological Shiba bands in artificial spin chains on superconductors. Nat. Phys. 17, 943–948 (2021).
Friedrich, F., Boshuis, R., Bode, M. & Odobesko, A. Coupling of Yu–Shiba–Rusinov states in one-dimensional chains of Fe atoms on Nb(110). Phys. Rev. B 103, 235437 (2021).
Küster, F. et al. Non-Majorana modes in diluted spin chains proximitized to a superconductor. Proc. Natl Acad. Sci. USA 119, e2210589119 (2022).
Odobesko, A. B. et al. Preparation and electronic properties of clean superconducting Nb(110) surfaces. Phys. Rev. B 99, 115437 (2019).
Odobesko, A. et al. Observation of tunable single-atom Yu–Shiba–Rusinov states. Phys. Rev. B 102, 174504 (2020).
Küster, F., Brinker, S., Lounis, S., Parkin, S. S. P. & Sessi, P. Long range and highly tunable interaction between local spins coupled to a superconducting condensate. Nat. Commun. 12, 6722 (2021).
Stroscio, J. A. & Eigler, D. M. Atomic and molecular manipulation with the scanning tunneling microscope. Science 254, 1319–1326 (1991).
Franke, K. J., Schulze, G. & Pascual, J. I. Competition of superconducting phenomena and Kondo screening at the nanoscale. Science 332, 940–944 (2011).
Yazdani, A., Jones, B. A., Lutz, C. P., Crommie, M. F. & Eigler, D. M. Probing the local effects of magnetic impurities on superconductivity. Science 275, 1767–1770 (1997).
Ruby, M., Peng, Y., von Oppen, F., Heinrich, B. W. & Franke, K. J. Orbital picture of Yu–Shiba–Rusinov multiplets. Phys. Rev. Lett. 117, 186801 (2016).
Choi, D.-J. et al. Mapping the orbital structure of impurity bound states in a superconductor. Nat. Commun. 8, 15175 (2017).
Schneider, L., Beck, P., Wiebe, J. & Wiesendanger, R. Atomic-scale spin-polarization maps using functionalized superconducting probes. Sci. Adv. 7, eabd7302 (2021).
Ruby, M., Heinrich, B. W., Peng, Y., von Oppen, F. & Franke, K. J. Wave-function hybridization in Yu–Shiba–Rusinov dimers. Phys. Rev. Lett. 120, 156803 (2018).
Beck, P. et al. Spin-orbit coupling induced splitting of Yu–Shiba–Rusinov states in antiferromagnetic dimers. Nat. Commun. 12, 2040 (2021).
Ding, H. et al. Tuning interactions between spins in a superconductor. Proc. Natl Acad. Sci. USA 118, e2024837118 (2021).
Kezilebieke, S., Dvorak, M., Ojanen, T. & Liljeroth, P. Coupled Yu–Shiba–Rusinov states in molecular dimers on NbSe2. Nano Lett. 18, 2311–2315 (2018).
Choi, D.-J. et al. Influence of magnetic ordering between Cr adatoms on the Yu–Shiba–Rusinov states of the β−Bi2Pd superconductor. Phys. Rev. Lett. 120, 167001 (2018).
Liebhaber, E. et al. Quantum spins and hybridization in artificially-constructed chains of magnetic adatoms on a superconductor. Nat. Commun. 13, 2160 (2022).
Huang, H. et al. Tunnelling dynamics between superconducting bound states at the atomic limit. Nat. Phys. 16, 1227–1231 (2020).
Plotnik, Y. et al. Observation of unconventional edge states in ‘photonic graphene’. Nat. Mater. 13, 57–62 (2014).
Imhof, S. et al. Topolectrical-circuit realization of topological corner modes. Nat. Phys. 14, 925–929 (2018).
Freeney, S. E., van den Broeke, J. J., Harsveld van der Veen, A. J. J., Swart, I. & Morais Smith, C. Edge-dependent topology in Kekulé lattices. Phys. Rev. Lett. 124, 236404 (2020).
Papanikolaou, N., Zeller, R. & Dederichs, P. H. Conceptual improvements of the KKR method. J. Phys. Condens. Matter 14, 2799 (2002).
Bauer, D. S. G. Development of a Relativistic Full-Potential First-Principles Multiple Scattering Green Function Method Applied to Complex Magnetic Textures of Nanostructures at Surfaces. PhD thesis, RWTH Aachen (2014).
Vosko, S. H., Wilk, L. & Nusair, M. Accurate spin-dependent electron liquid correlation energies for local spin density calculations: a critical analysis. Can. J. Phys. 58, 1200–1211 (1980).
Liechtenstein, A. I., Katsnelson, M. I., Antropov, V. P. & Gubanov, V. A. Local spin density functional approach to the theory of exchange interactions in ferromagnetic metals and alloys. J. Magn. Magn. Mater. 67, 65–74 (1987).
Ebert, H. & Mankovsky, S. Anisotropic exchange coupling in diluted magnetic semiconductors: ab initio spin-density functional theory. Phys. Rev. B 79, 045209 (2009).
Acknowledgements
A.A. and S.L. thank M. dos Santos Dias and S. Brinker for fruitful discussions. G.W. acknowledges NCCR MARVEL funding from the Swiss National Science Foundation. T.N. acknowledges funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (ERC-StG-Neupert-757867-PARATOP). M.O.S. acknowledges funding from the Swiss National Science Foundation (Project 200021E_198011) as part of the FOR 5249 (QUAST) lead by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation). A.A. is funded by the the Palestinian-German Science Bridge (BMBF grant no. 01DH16027) and her DFT simulations were made on the supercomputer JURECA at Forschungszentrum Jülich with computing time granted through JARA. R.T. acknowledges the Wuerzburg-Dresden Cluster of Excellence ct.qmat.
Author information
Authors and Affiliations
Contributions
F.K., S.D. and P.S. conceived and performed the experiment. M.O.S., G.W., R.T. and T.N. carried out the theoretical analysis. A.A. and S.L. performed the first-principles calculations. M.O.S., G.W., F.K., P.S., T.N. and S.S.P.P. wrote the paper with input from all the authors.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Physics thanks the anonymous reviewers for their contribution to the peer review of this work
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Extended data
Extended Data Fig. 1 Scanning tunneling spectroscopy measured with a superconducting Nb tip on clean Nb(110).
A bulk superconducting Nb cluster at the tip was obtained by deep indentation into the sample. a Tunneling spectroscopy shows the convolution of tip and sample density of states resulting in a large gap with size 2(Δtip + Δsample)/e. From the measured convoluted gap size of 3.05 mV and the database value for bulk Nb (ΔNb = 1.52 meV), we obtain a tip gap approximately equal to the bulk value, that is Δtip ≊ ΔNb, which is indicated by the gray area. Measurement parameters: stabilized at sample bias -5 mV b, tunneling current 500 pA, bias AC modulation amplitude 40 μV, temperature 500 mK. b Zero bias differential conductance peak. The measured FWHM of 135 μV provides a reference for our experimental energy resolution. Measurement parameters: stabilized at sample bias -5 mV, tunneling current 30 nA, bias AC modulation amplitude 10 μV, temperature 500 mK.
Extended Data Fig. 2 Spectroscopy with spin-sensitive tip used for Fig. 1e,g of the main text.
Both measurements are acquired on the clean Nb substrate with a superconducting tip featuring single Cr atoms at the apex that induce Shiba states indicated by green arrows for the zero-field measurement (black curve). The applied external magnetic field of 0.7 T in-plane overcomes the second critical field for the Nb sample while superconductivity is maintained in the Nb cluster at the tip (blue curve) where Shiba states are still visible as shoulders on the flank of the gap. The observed spin contrast in constant-height dI/dU maps was strongest for the indicated bias at 0.33 mV.
Extended Data Fig. 3 Deconvoluted tunneling spectroscopy data.
Tunneling spectroscopy data corresponding to Fig. 2 of the main article after performing numerical deconvolution, and the panels correspond to the terminations A1 (a), A2 (b), A3 (c), Β1 (d) and Β2 (e). Respective lattice positions are indicated by colored circles corresponding to STS lines. Numerical deconvolution was carried out following the process described in Ref. 15.
Extended Data Fig. 4 Formation of two dimensional Shiba bands inside the bulk of artificial spin structures.
a, Topographic image from the full spectroscopic measurement grid on structure A1. b–d, Two representative dI/dU maps for bulk Shiba bands at sample biases -2.37 mV (b), -2.21 mV (c) and -2.09 mV (d), demonstrating their 2D character by intensity modulations in both orthogonal directions ([001] and [\(1\overline{1}0\)]).
Extended Data Fig. 5 Analysis of Shiba orbital characters.
a, Topographic image of the A1 lattice analyzed in the main text. b–c, dI/dU maps acquired at biases corresponding to the two most prominent peaks in the dI/dU spectrum reported in d. The spectrum has been obtained by averaging over the region identified by a dashed black box in a. Two peaks are dominating the scene inside the superconducting gap, centered at U = −2.24mV and U = −2.24mV (see black lines). Their strong intensity allow assigning them to d\({}_{{{{{\rm{z}}}}}^{2}}\) -derived states which can be effectively measured by STS because of their large extension into the vacuum. This orbital assignment is corroborated by the respective dI/dU maps, which show two distinct patterns, both characterized by intensity maxima centered around the position of the adatoms, as expected for d\({}_{{{{{\rm{z}}}}}^{2}}\)-derived states.
Extended Data Fig. 6 Averaged bulk signals of structures A and B.
a,b, Topography images of structures A and B with the respective bulk area indicated by a dashed rectangle used for creating a bulk STS signal which is shown in b. c, Plots each signal after deconvolution by a model tip DOS corresponding to a Nb superconducting gap in order to get a good approximation of the sample DOS. Structure B has clearly no gap in the bulk while A features a dip at zero energy which is very close to zero intensity.
Extended Data Fig. 7 Tunneling spectroscopy data acquired in several points along straight lines over the A3 termination.
a, Topography image of structure A3 with two dashed lines indicating STS line paths A (through bulk) and B (along edge) as well as the corner C. b,c Corresponding dI/dU spectra along those lines. At the corner C, indicated by a black solid arrow, we observe the strongest residual intensity at -Δtip (dashed black line) corresponding to zero energy.
Extended Data Fig. 8 Shiba bands along [001] in momentum space.
Evolution of Shiba bands from 1D to 2D comparing a long 1D wire of Cr adatoms along the [001] direction and 0.66 nm inter-atomic distance to a two dimensional arrangement according to structure A. a,b, Topography images. c,d, FFT calculated from the STS data taken along the [001] direction. For better visibility, the plots are separated to show the energy range −4mV… − Δtip and Δtip…4mV with adjusted dI/dU colorscale contrast. The momentum space plots suggest dispersing bands. The signal with highest intensity around ± 2mV is assigned to the dz2 atomic orbital (white dashed lines). Other bands arise from different single-atomic orbitals (red dashed line). From 1D to 2D, a change of sign in the dispersion relation is observed. This demonstrates the potential to engineer the band structure by advancement into two dimensions.
Supplementary information
Supplementary Information
Supplementary Sections I and II.
Source data
Source Data Fig. 2
Spectroscopic data for lattice terminations A1, A2, A3, B1 and B2. Zeroth column: sample bias (meV). Columns 1–4: dI/dU (nS) data for termination A1 (Fig. 2a). The four columns correspond to corner lower left, edge bottom, edge left and bulk. Columns 5–8: dI/dU (nS) data for termination A2 (Fig. 2b). The four columns correspond to corner lower, bulk, corner left and edge lower left. Columns 9–12: dI/dU (nS) data for termination A3 (Fig. 2c). The four columns correspond to edge lower left, bulk, corner left and corner lower. Columns 13–16: dI/dU (nS) data for termination B1 (Fig. 2d). The four columns correspond to bulk, lower corner, lower left edge and left corner. Columns 17–20: dI/dU (nS) data for termination B2 (Fig. 2e). The four columns correspond to lower edge, lower left corner, left edge and bulk.
Source Data Fig. 3
Magnetic coupling extracted from ab initio calculations. Column 0: distance from the reference site. Column 1: distance along the x direction from the reference site. Column 2: distance along the y direction from the reference site. Column 3: ab initio Heisenberg exchange interaction J (meV). Rows 0–21: positions compatible with the rectangular structure (A). Rows 22–39: positions compatible with the rhombic structure (B).
Source Data Fig. 4
Source data for Fig. 4b–d,j. It comprises the following files: data_ribbons_Fig4.txt (for Fig. 4b,c,j) and Spectrum_A3_Fig4d_lower.txt and Spectrum_A3_Fig4d_upper.txt (for Fig. 4d). File data_ribbons_Fig4.txt: ribbon spectra for structures A1 (Fig. 4b), A2 (Fig. 4c) and B1 (Fig. 4j). Ribbon spectra shown in Fig. 4b,c,j. The rows correspond to different values of the momentum. The nth row: k = –π + 2π × n/200. Columns 0–959: eigenvalues of the ribbon structure shown in Fig. 4b. Columns 960–1,439: eigenvalues of the ribbon structure shown in Fig. 4c. Columns 1,440–1,915: eigenvalues of the ribbon structure shown in Fig. 4j. File Spectrum_A3_Fig4d_lower.txt: energy spectrum for A3 termination with n = 101 rows. File Spectrum_A3_Fig4d_upper.txt: energy spectrum for A3 termination with n = 17 rows.
Source Data Extended Data Fig. 1
Source data for Extended Data Fig. 1a,b. The zip files contains ED1_a.txt (Extended Data Fig. 1a) and ED1_b.txt (Extended Data Fig. 1b). File ED1_a.txt: scanning tunnelling spectroscopy measured with a superconducting Nb tip on clean Nb(110). Zeroth column: sample bias (mV). First column: dI/dU (µS). File ED1_b.txt: scanning tunnelling spectroscopy measured with a superconducting Nb tip on clean Nb(110). Zeroth column: sample bias (mV). First column: dI/dU (µS).
Source Data Extended Data Fig. 2
Spectroscopy with spin-sensitive tip used for Extended Data Fig. 1e,g. Zeroth column: sample bias (mV). First column: dI/dU (µS) with 0 T external magnetic field. Second column: dI/dU (µS) with 0.7 T external magnetic field applied in plane.
Source Data Extended Data Fig. 3
Deconvoluted spectroscopic data for lattice terminations A1, A2, A3, B1 and B2. Zeroth column: sample bias (meV). Columns 1–4: density of states (DOS) (arb. units) data for termination A1 (Extended Data Fig. 3a). The four columns correspond to corner lower left, edge bottom, edge left and bulk. Columns 5–8: DOS (arb. units) data for termination A2 (Extended Data Fig. 3b). The four columns correspond to corner lower, bulk, corner left and edge lower left. Columns 9–12: DOS (arb. units) data for termination A3 (Extended Data Fig. 3c). The four columns correspond to edge lower left, bulk, corner left and corner lower. Columns 13–16: DOS (arb. units) data for termination B1 (Extended Data Fig. 3d). The four columns correspond to bulk, lower corner, lower left edge and left corner. Columns 17–20: DOS (arb. units) data for termination B2 (Extended Data Fig. 3e). The four columns correspond to lower edge, lower left corner, left edge and bulk.
Source Data Extended Data Fig. 5
Analysis of Shiba orbital characters. Zeroth column: sample bias (mV). First column: dI/dU (nS).
Source Data Extended Data Fig. 6
Two files with source data for Extended Data Fig. 6b,c, namely, files ED_6b.txt and ED_6c.txt, respectively. File ED_6b.txt: averaged bulk signals of structures A and B. Zeroth column: sample bias (mV). First column: dI/dU (nS), structure A. Second column: dI/dU (nS), structure B. File ED_6c.txt: averaged bulk signals of structures A and B deconvoluted. Zeroth column: energy with respect to the Fermi energy (meV). First column: sample LDOS (arb. units), structure A. Second column: sample LDOS (arb. units), structure B.
Rights and permissions
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Soldini, M.O., Küster, F., Wagner, G. et al. Two-dimensional Shiba lattices as a possible platform for crystalline topological superconductivity. Nat. Phys. 19, 1848–1854 (2023). https://doi.org/10.1038/s41567-023-02104-5
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/s41567-023-02104-5
