Parity-dependent coupling of molecular spin chains to a superconductor

Data availability

All data that support the plots in this paper and other findings of this study are available from the corresponding author. Data analysis was done with WsXM54 and SpectraFox55.

References

  1. Haldane, F. D. M. Nonlinear field theory of large-spin Heisenberg antiferromagnets: semiclassically quantized solitons of the one-dimensional easy-axis Néel state. Phys. Rev. Lett. 50, 1153–1156 https://doi.org/10.1103/PhysRevLett.50.1153 (1983).

    Article  ADS  MathSciNet  Google Scholar 

  2. Affleck, I., Kennedy, T., Lieb, E. H. & Tasaki, H. Rigorous results on valence-bond ground states in antiferromagnets. Phys. Rev. Lett. 59, 799 https://doi.org/10.1103/PhysRevLett.59.799 (1987).

    Article  ADS  CAS  PubMed  Google Scholar 

  3. Su, W. P., Schrieffer, J. R. & Heeger, A. J. Solitons in polyacetylene. Phys. Rev. Lett. 42, 1698 https://doi.org/10.1103/PhysRevLett.42.1698 (1979).

    Article  ADS  CAS  Google Scholar 

  4. Jackiw, R. & Rebbi, C. Solitons with fermion number 1/2. Phys. Rev. D 13, 3398 https://doi.org/10.1103/PhysRevD.13.3398 (1976).

    Article  ADS  MathSciNet  CAS  Google Scholar 

  5. Kitaev, A. Y. Unpaired Majorana fermions in quantum wires. Phys.-Uspekhi 44, 131–136 https://doi.org/10.1070/1063-7869/44/10S/S29 (2001).

    Article  ADS  Google Scholar 

  6. Tasaki, H. Ground state of the S = 1 antiferromagnetic Heisenberg chain is topologically nontrivial if gapped. Phys. Rev. Lett. 134, 076602 https://doi.org/10.1103/PhysRevLett.134.076602 (2025).

    Article  ADS  MathSciNet  CAS  PubMed  Google Scholar 

  7. Miyake, A. Quantum computation on the edge of a symmetry-protected topological order. Phys. Rev. Lett. 105, 040501 https://doi.org/10.1103/PhysRevLett.105.040501 (2010).

    Article  ADS  CAS  PubMed  Google Scholar 

  8. Jaworowski, B., Rogers, N., Grabowski, M. & Hawrylak, P. Macroscopic singlet-triplet qubit in synthetic spin-one chain in semiconductor nanowires. Sci. Rep. 7, 5529 https://doi.org/10.1038/s41598-017-05655-9 (2017).

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  9. Brennen, G. K. & Miyake, A. Measurement-based quantum computer in the gapped ground state of a two-body Hamiltonian. Phys. Rev. Lett. 101, 010502 https://doi.org/10.1103/PhysRevLett.101.010502 (2008).

    Article  ADS  CAS  PubMed  Google Scholar 

  10. del Castillo, Y., Ferrón, A. & Fernández-Rossier, J. Remote spin control in Haldane spin chains. SciPost Phys. Core 9, 024 https://doi.org/10.21468/SciPostPhysCore.9.2.024 (2026).

  11. Huang, C.-H., Ortuzar, J. & Cazalilla, M. A. Superconducting spin-singlet qubit in a triangulene spin chain. Phys. Rev. B 113, 064503 https://doi.org/10.1103/58rh-rjml (2026).

    Article  ADS  CAS  Google Scholar 

  12. Baran, V. V. & Paaske, J. Spin-1 Haldane chains of superconductor-semiconductor hybrids. Phys. Rev. B 110, 064503 https://doi.org/10.1103/PhysRevB.110.064503 (2024).

    Article  ADS  CAS  Google Scholar 

  13. ten Haaf, S. L. D. et al. Observation of edge and bulk states in a three-site Kitaev chain. Nature 641, 890–895 https://doi.org/10.1038/s41586-025-08892-5 (2025).

    Article  CAS  PubMed  Google Scholar 

  14. Clair, S. & de Oteyza, D. G. Controlling a chemical coupling reaction on a surface: tools and strategies for on-surface synthesis. Chem. Rev. 119, 4717–4776 https://doi.org/10.1021/acs.chemrev.8b00601 (2019).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Cao, T., Zhao, F. & Louie, S. G. Topological phases in graphene nanoribbons: junction states, spin centers, and quantum spin chains. Phys. Rev. Lett. 119, 076401 https://doi.org/10.1103/PhysRevLett.119.076401 (2017).

    Article  ADS  PubMed  Google Scholar 

  16. Rizzo, D. J. et al. Topological band engineering of graphene nanoribbons. Nature 560, 204–208 https://doi.org/10.1038/s41586-018-0376-8 (2018).

    Article  ADS  CAS  PubMed  Google Scholar 

  17. Gröning, O. et al. Engineering of robust topological quantum phases in graphene nanoribbons. Nature 560, 209–213 https://doi.org/10.1038/s41586-018-0375-9 (2018).

    Article  ADS  CAS  PubMed  Google Scholar 

  18. Cirera, B. et al. Tailoring topological order and π-conjugation to engineer quasi-metallic polymers. Nat. Nanotechnol. 15, 437–443 https://doi.org/10.1038/s41565-020-0668-7 (2020).

    Article  ADS  CAS  PubMed  Google Scholar 

  19. Li, J. et al. Topological phase transition in chiral graphene nanoribbons: from edge bands to end states. Nat. Commun. 12, 5538 https://doi.org/10.1038/s41467-021-25688-z (2021).

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  20. Mishra, S. et al. Observation of fractional edge excitations in nanographene spin chains. Nature 598, 287–292 https://doi.org/10.1038/s41586-021-03842-3 (2021).

    Article  ADS  CAS  PubMed  Google Scholar 

  21. Zhao, C. et al. Tunable topological phases in nanographene-based spin-1/2 alternating-exchange Heisenberg chains. Nat. Nanotechnol. 19, 1789–1795 https://doi.org/10.1038/s41565-024-01805-z (2024).

    Article  ADS  CAS  PubMed  Google Scholar 

  22. Zhao, C. et al. Spin excitations in nanographene-based antiferromagnetic spin-1/2 Heisenberg chains. Nat. Mater. 24, 722–727 https://doi.org/10.1038/s41563-025-02166-1 (2025).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Sun, K. et al. On-surface synthesis of Heisenberg spin-1/2 antiferromagnetic molecular chains. Sci. Adv. 11, eads1641 https://doi.org/10.1126/sciadv.ads1641 (2025).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Su, X. et al. Fabrication of spin-1/2 Heisenberg antiferromagnetic chains via combined on-surface synthesis and reduction for spinon detection. Nat. Synth. 4, 694–701 https://doi.org/10.1038/s44160-025-00744-4 (2025).

    Article  CAS  Google Scholar 

  25. Yuan, Z. et al. Fractional spinon quasiparticles in open-shell triangulene spin-1/2 chains. J. Am. Chem. Soc. 147, 5004–5013 https://doi.org/10.1021/jacs.4c14712 (2025).

    Article  ADS  CAS  PubMed  Google Scholar 

  26. Heinrich, B. W., Braun, L., Pascual, J. I. & Franke, K. J. Protection of excited spin states by a superconducting energy gap. Nat. Phys. 9, 765–768 https://doi.org/10.1038/nphys2794 (2013).

    Article  CAS  Google Scholar 

  27. Vaxevani, K. et al. Extending the spin excitation lifetime of a magnetic molecule on a proximitized superconductor. Nano Lett. 22, 6075–6082 https://doi.org/10.1021/acs.nanolett.2c00924 (2022).

    Article  ADS  CAS  PubMed  Google Scholar 

  28. Yu, L. Bound state in superconductors with paramagnetic impurities. Acta Phys. Sin. 21, 75–91 https://doi.org/10.7498/aps.21.75 (1965).

    Article  CAS  Google Scholar 

  29. Shiba, H. Classical spins in superconductors. Prog. Theor. Phys. 40, 435–451 https://doi.org/10.1143/PTP.40.435 (1968).

    Article  ADS  Google Scholar 

  30. Rusinov, A. Theory of gapless superconductivity in alloys containing paramagnetic impurities. Sov. Phys. JETP 29, 1101–1106 (1969).

    ADS  Google Scholar 

  31. Heinrich, B. W., Pascual, J. I. & Franke, K. J. Single magnetic adsorbates on s-wave superconductors. Prog. Surf. Sci. 93, 1–19 https://doi.org/10.1016/j.progsurf.2018.01.001 (2018).

    Article  ADS  CAS  Google Scholar 

  32. Liu, J.-C. et al. Proximity-induced superconductivity in atomically precise nanographene on Ag/Nb(110). ACS Mater. Lett. 5, 1083–1090 https://doi.org/10.1021/acsmaterialslett.2c00955 (2023).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Wang, D. et al. Epitaxial growth of gold films on the elemental superconductors V(100), Nb(100), and Nb(110). Phys. Rev. Mater. 9, 066201 https://doi.org/10.1103/PhysRevMaterials.9.066201 (2025).

    Article  CAS  Google Scholar 

  34. Delgado, F., Batista, C. D. & Fernández-Rossier, J. Local probe of fractional edge states of S = 1 Heisenberg spin chains. Phys. Rev. Lett. 111, 167201 https://doi.org/10.1103/PhysRevLett.111.167201 (2013).

    Article  ADS  CAS  PubMed  Google Scholar 

  35. Wei, P., Manna, S., Eich, M., Lee, P. & Moodera, J. Superconductivity in the surface state of noble metal gold and its Fermi level tuning by EuS dielectric. Phys. Rev. Lett. 122, 247002 https://doi.org/10.1103/PhysRevLett.122.247002 (2019).

    Article  ADS  CAS  PubMed  Google Scholar 

  36. Vaxevani, K. et al. Inducing superconductivity in atomically thin Au(111) films via interface engineering. Zenodo https://doi.org/10.5281/zenodo.22815538 (2026).

  37. de Gennes, P. & Saint-James, D. Elementary excitations in the vicinity of a normal metal-superconducting metal contact. Phys. Lett. 4, 151–152 https://doi.org/10.1016/0031-9163(63)90148-3 (1963).

    Article  ADS  Google Scholar 

  38. Schneider, L. et al. Proximity superconductivity in atom-by-atom crafted quantum dots. Nature 621, 60–65 https://doi.org/10.1038/s41586-023-06312-0 (2023).

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  39. Ortuzar, J., Pascual, J. I., Bergeret, F. S. & Cazalilla, M. A. Theory of a single magnetic impurity on a thin metal film in proximity to a superconductor. Phys. Rev. B 108, 024511 https://doi.org/10.1103/PhysRevB.108.024511 (2023).

    Article  ADS  CAS  Google Scholar 

  40. Trivini, S. et al. Diluted Yu-Shiba-Rusinov arrays on the anisotropic superconductor β− Bi2Pd. Phys. Rev. B 110, 235405 https://doi.org/10.1103/PhysRevB.110.235405 (2024).

    Article  ADS  CAS  Google Scholar 

  41. Gross, L. et al. High-resolution molecular orbital imaging using a p-wave STM tip. Phys. Rev. Lett. 107, 086101 https://doi.org/10.1103/PhysRevLett.107.086101 (2011).

    Article  ADS  CAS  PubMed  Google Scholar 

  42. Li, J. et al. Single spin localization and manipulation in graphene open-shell nanostructures. Nat. Commun. 10, 200 https://doi.org/10.1038/s41467-018-08060-6 (2019).

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  43. Aligia, A. A., Lobos, A. M., Peralta Gavensky, L. & Gazza, C. J. Probing boundary spins in the Su-Schrieffer-Heeger-Hubbard model. Phys. Rev. B 112, 235156 https://doi.org/10.1103/vhfn-6f3d (2025).

    Article  ADS  CAS  Google Scholar 

  44. Franke, K. J., Schulze, G. & Pascual, J. I. Competition of superconducting phenomena and Kondo screening at the nanoscale. Science 332, 940–944 https://doi.org/10.1126/science.1202204 (2011).

    Article  ADS  CAS  PubMed  Google Scholar 

  45. Turco, E. et al. Demonstrating Kondo behavior by temperature-dependent scanning tunneling spectroscopy. Phys. Rev. Res. 6, L022061 https://doi.org/10.1103/PhysRevResearch.6.L022061 (2024).

    Article  CAS  Google Scholar 

  46. von Oppen, F. & Franke, K. J. Yu-Shiba-Rusinov states in real metals. Phys. Rev. B 103, 205424 https://doi.org/10.1103/PhysRevB.103.205424 (2021).

    Article  ADS  Google Scholar 

  47. Trivini, S. et al. Cooper pair excitation mediated by a molecular quantum spin on a superconducting proximitized gold film. Phys. Rev. Lett. 130, 136004 https://doi.org/10.1103/PhysRevLett.130.136004 (2023).

    Article  ADS  CAS  PubMed  Google Scholar 

  48. Zhao, C. et al. Tailoring magnetism of graphene nanoflakes via tip-controlled dehydrogenation. Phys. Rev. Lett. 132, 046201 https://doi.org/10.1103/PhysRevLett.132.046201 (2024).

    Article  ADS  CAS  PubMed  Google Scholar 

  49. Spinelli, A., Bryant, B., Delgado, F., Fernández-Rossier, J. & Otte, A. F. Imaging of spin waves in atomically designed nanomagnets. Nat. Mater. 13, 782–785 https://doi.org/10.1038/nmat4018 (2014).

    Article  CAS  PubMed  Google Scholar 

  50. Drost, R., Kezilebieke, S., Lado, J. L. & Liljeroth, P. Real-space imaging of triplon excitations in engineered quantum magnets. Phys. Rev. Lett. 131, 086701 https://doi.org/10.1103/PhysRevLett.131.086701 (2023).

    Article  ADS  CAS  PubMed  Google Scholar 

  51. Machens, A., Konstantinidis, N. P., Waldmann, O., Schneider, I. & Eggert, S. Even-odd effect in short antiferromagnetic Heisenberg chains. Phys. Rev. B 87, 144409 https://doi.org/10.1103/PhysRevB.87.144409 (2013).

    Article  ADS  CAS  Google Scholar 

  52. Landé, A. Termstruktur und Zeemaneffekt der Multipletts. Z. Phys. 15, 189–205 https://doi.org/10.1007/BF01330473 (1923).

    Article  ADS  Google Scholar 

  53. Haldane, F. D. M. Continuum dynamics of the 1-D Heisenberg antiferromagnet: identification with the O(3) nonlinear sigma model. Phys. Lett. A 93, 464–468 https://doi.org/10.1016/0375-9601(83)90631-X (1983).

    Article  ADS  MathSciNet  Google Scholar 

  54. Horcas, I. et al. WSXM: a software for scanning probe microscopy and a tool for nanotechnology. Rev. Sci. Instrum. 78, 013705 https://doi.org/10.1063/1.2432410 (2007).

  55. Ruby, M. SpectraFox: a free open-source data management and analysis tool for scanning probe microscopy and spectroscopy. SoftwareX 5, 31–36 https://doi.org/10.1016/j.softx.2016.04.001 (2016).

Download references

Acknowledgements

We thank L. Edens for technical support and J. Wolf for help with the substrate growth. We also thank the Quantum-Probe Microscopy group for helpful discussions.

Funding

This research was financed by MCIN/AEI/10.13039/501100011033 and the European Regional Development Fund (ERDF) through grant nos. PID2022-140845OB-C61, PID2022-140845OB-C62 and CEX2020-001038-M, by Xunta de Galicia (Centro de Investigación do Sistema Universitario de Galicia, 2023–2027, grant number ED431G 2023/03) and by the European Union Horizon 2020 and HE programmes through FET-Open project SPRING (no. 863098), the ERC Synergy Grant MolDAM (no. 951519) and the ERC-AdG CONSPIRA (no. 101097693). K.V. acknowledges support from MCIN/AEI/10.13039/501100011033 through PhD scholarship PRE2022-102887. J.O. acknowledges support from the Basque Government through PhD scholarship PRE_2021_1_0350. D.W. acknowledges support from ‘Juan de la Cierva’ grant FJC2020-043831-I financed by MCIN/AEI/10.13039/501100011033 and by the European Union NextGenerationEU/PRTR, as well as from NSFC project no. 12474474.

Author information

Author notes

  1. These authors contributed equally: Katerina Vaxevani, Jon Ortuzar, Stefano Trivini

Authors and Affiliations

  1. CIC nanoGUNE - BRTA, Donostia-San Sebastián, Spain

    Katerina Vaxevani, Jon Ortuzar, Stefano Trivini, Georg Monninger, Dongfei Wang & Jose Ignacio Pascual

  2. Quantronics group, Service de Physique de l’État Condensé (CNRS, UMR 3680), IRAMIS, CEA-Saclay, Université Paris-Saclay, Gif-sur-Yvette, France

    Jon Ortuzar

  3. School of Physics, Beijing Institute of Technology, Beijing, China

    Dongfei Wang

  4. Centro Singular de Investigación en Química Biolóxica e Materiais Moleculares (CiQUS) and Departamento de Química Orgánica, Universidad de Santiago de Compostela, Santiago de Compostela, Spain

    Manuel Vilas-Varela, Lucía Gómez-Rodrigo & Diego Peña

  5. Oportunius, Galician Innovation Agency (GAIN), Santiago de Compostela, Spain

    Diego Peña

  6. Ikerbasque, Basque Foundation for Science, Bilbao, Spain

    Jose Ignacio Pascual

Authors

  1. Katerina Vaxevani
  2. Jon Ortuzar
  3. Stefano Trivini
  4. Georg Monninger
  5. Dongfei Wang
  6. Manuel Vilas-Varela
  7. Lucía Gómez-Rodrigo
  8. Diego Peña
  9. Jose Ignacio Pascual

Contributions

J.I.P. conceived of and supervised the project. M.V.-V. and L.G.-R. synthesized the molecular precursors in solution under the supervision of D.P. K.V., G.M. and D.W. developed the superconducting substrate and grew the triangulene chains. J.O., S.T. and G.M. performed the STM/STS measurements. K.V., J.O., S.T. and G.M. analysed the data. J.O. developed the theoretical model and simulations. J.I.P. wrote the manuscript, with input from K.V., J.O. and S.T. All authors discussed the results and contributed to the final manuscript.

Corresponding author

Correspondence to Jose Ignacio Pascual.

Ethics declarations

Competing interests

The authors declare no competing interests.

Peer review

Peer review information

Nature thanks Ping Yu and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.

Additional information

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

Extended data figures and tables

Extended Data Fig. 1 Results for N = 3 and N = 4 [3]-triangulene chains.

a, Bond-resolved current image of triangulene chain with three units (V = 5 mV, constant-height at tunnel resistance RT ≈ 16 MΩ). b, Low-energy dI/dV spectra over the termini. Dashed grey line is a reference spectrum on the Au(111) substrate. Spectra on the chain termini show YSR peaks developing inside the dGSJ excitations. c, Bond-resolved current image of triangulene chain with four units (V = 8.5 mV, constant-height at RT ≈ 16 MΩ). d, Low-energy dI/dV spectra over the termini. Dashed grey line is a reference spectrum on the Au(111) substrate. Spectra on the chain termini show peaks at about 8.5 meV, corresponding to ST excitations with energy EST = 6.3 meV. e, Differential conductance map at the peak energy, measured simultaneously to the image in c, showing larger conductance over the edge triangulene units.

Extended Data Fig. 2 Results for a N = 5 [3]-triangulene chain.

a, Bond-resolved current image of triangulene chain with five units (V = 5 mV, constant-height at RT ≈ 16 MΩ). b, Low-energy dI/dV spectra over the termini. Dashed grey line is a reference spectrum on the Au(111) substrate. Spectra on the chain termini show YSR states developing from the dGSJ peaks and TS excitation peaks, shown in c. c, Long-range dI/dV spectra averaged at several positions over the edge and the bulk of the triangulene chain in a. d, Differential conductance maps of the triangulene chain in a. The maps are chosen at bias values of the YSR peak, the negative differential resistance (NDR) dip after the dGSJ peaks, the TS excitation and two CS excitation steps, respectively. Both c and d were obtained from a spectral grid over the whole triangulene chain.

Extended Data Fig. 3 Results for N = 6 and N = 7 [3]-triangulene chains.

a, Bond-resolved current image of triangulene chain with six units (V = 5 mV, constant-height at RT ≈ 16 MΩ). b, Low-energy dI/dV spectra over the termini. Dashed grey line is a reference spectrum on the Au(111) substrate. Spectra on the chain termini show ST excitation peaks with EST = 3.1 meV. c, Bond-resolved current image of triangulene chain with seven units (V = 5 mV, constant-height at RT ≈ 16 MΩ). d, Low-energy dI/dV spectra over the termini. Dashed grey line is a reference spectrum on the Au(111) substrate. Spectra on the chain termini show YSR excitation peaks developing inside the gap between dGSJ peaks. Weak TS excitation peaks are observed outside the gap, with ETS = 2.7 meV.

Extended Data Fig. 4 Results for a N = 8 [3]-triangulene chain.

a, Constant-current image of triangulene chain with eight units (V = 5 mV). b, Low-energy dI/dV spectra over the termini. Dashed grey line is a reference spectrum on the Au(111) substrate. Spectra on the chain termini show ST excitation peaks with EST = 1.9 meV. c, Comparison of spin excitation spectra for two distinct N = 8 chains. Both host the many-body S = 1 spin excitations in the same frequencies. The peak at about 35 meV is related to the vibrations of the CO at the scanning tunnelling microscope tip.

Extended Data Fig. 5 Results for N = 9 and N = 10 [3]-triangulene chains.

a, Bond-resolved current image of triangulene chain with nine units (V = 5 mV, constant-height at RT ≈ 16 MΩ). b, Low-energy dI/dV spectra over the termini. Dashed grey line is a reference spectrum on the Au(111) substrate. Spectra on the chain termini show YSR excitation peaks developing inside the dGSJ excitations. Weak TS excitation peaks are observed outside the gap, with ETS = 1.6 meV. c, Bond-resolved current image of triangulene chain with 10 units (V = 5 mV, constant-height at RT ≈ 16 MΩ). d, Low-energy dI/dV spectra over the termini. Dashed grey line is a reference spectrum on the Au(111) substrate. Spectra on the chain termini show ST excitation peaks with EST = 1.0 meV.

Extended Data Fig. 6 Results for a N = 11 [3]-triangulene chain.

a, Low-energy dI/dV spectra over both termini of the N = 11 triangulene chain in Fig. 2d. b, Low-energy dI/dV spectra over both termini of another N = 11 triangulene chain with less YSR bound energy δEYSR. In both plots, the insets are bond-resolved images and dashed grey lines are reference spectra on the Au(111) substrate. Both chains show weak TS excitation shoulders as marked, with ETS ≈ 0.8 meV, the one in b with larger peak amplitude. c, Spatial distribution of YSR peak amplitude and TS onset over the top-end triangulene in a, obtained from a 20 × 20 spectral grid. Both signals are distributed around the zigzag periphery, with similar shape at both polarities.

Extended Data Fig. 7 Results for N = 12 and N = 13 [3]-triangulene chains.

a, Bond-resolved current image of triangulene chain with 12 units (V = 5 mV, constant-height at RT ≈ 16 MΩ). b, Low-energy dI/dV spectra over the termini. Dashed grey line is a reference spectrum on the Au(111) substrate. Spectra on the chain termini show faint ST excitation peaks with EST ≈ 0.5 meV. The dGSJ resonances appear with particle–hole asymmetry characteristic of YSR states with finite Coulomb potential scattering. Here we attribute it to the incipient emergence of YSR states, still overlapping with the quasiparticle peaks of the film. c, Bond-resolved current image of triangulene chain with 13 units (V = 5 mV, constant-height at RT ≈ 16 MΩ). d, Low-energy dI/dV spectra over the termini. Dashed grey line is a reference spectrum on the Au(111) substrate. Spectra on the chain termini show YSR excitation peaks developing inside the dGSJ excitations. Spectra on the chain termini show faint TS excitation peaks with ETS ≈ 0.5 meV.

Extended Data Fig. 8 Results for a N = 14 [3]-triangulene chain.

a, Bond-resolved current image of triangulene chain with 14 units (V = 5 mV, constant-height at RT ≈ 16 MΩ). b, Low-energy dI/dV spectra over the termini. Dashed grey line is a reference spectrum on the Au(111) substrate. In the inset, we show the scaled components of the dGSJ states for each termini, to emphasize the peak components owing to YSR state.

Extended Data Fig. 9 Edge state passivation by dehydrogenation of a terminal triangulene, N = 14.

a, Bond-resolved current STM image of the N = 14 chain. The insets show a zoomed image of the end triangulene unit n = 14 before and after applying a bias pulse over the marked C site. The pulse produces dehydrogenation of the targeted C site, after which the triangulene appears with characteristic orbital shape in the images (V = 5 mV). b, Differential conductance spectra on the end unit n = 1 before and after the dehydrogenation in a. The YSR peak of Extended Data Fig. 8 over the n = 1 triangulene shifts to a binding energy δEYSR = 0.36 meV when the opposite end is dehydrogenated (V = 5 mV, I = 500 pA). c, Differential conductance spectral map versus tip–sample resistance over the n = 1 triangulene unit. The YSR state shifts towards the gap edges and merges with dGSJ resonances as the tip approaches, indicating that the edge spin lies in the weak interaction regime (JK ≪ Δ; Supplementary Information section 5).

Extended Data Fig. 10 Edge state passivation by dehydrogenation of a terminal triangulene, N = 8.

a, Constant-current STM image of a N = 8 chain before and after dehydrogenation of the lower triangulene end. b, Comparison of spectra on the top triangulene before and after the dehydrogenation in a. The passivation of the lower edge state leads to disappearance of ST features and the emergence of an in-gap YSR state. As shown in Supplementary Information section 4, the H atom removal removes intra-chain exchange correlations and results in a free-like edge spin on the opposite end.

Extended Data Fig. 11 STM image and dI/dV spectra of the surface of a 25-monolayer Au(111) film.

a, STM image of the surface of a 25-monolayer Au(111) film, exhibiting the herringbone reconstruction. b, dI/dV spectra measured at 1.2 K, showing the characteristic onset of its surface state close to −0.5 V sample bias. This value is close to the surface state onset on bulk Au(111) substrates.

Supplementary information

Supplementary Information (download PDF )

This file contains five supplementary figures, further references and the following sections: 1. Deconvolution of the superconducting tip gap and YSR bound state energies; 2. Spin model of the S = 1 chain; 3. Two-site model of YSR excitations in fractional spin end states; 4. Theoretical analysis of the tip-induced passivation of fractional spin end states; and 5. Determination of interaction regime with the superconductor.

Peer Review File (download PDF )

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Vaxevani, K., Ortuzar, J., Trivini, S. et al. Parity-dependent coupling of molecular spin chains to a superconductor. Nature 658, 617–622 (2026). https://doi.org/10.1038/s41586-026-11117-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • Version of record:

  • Issue date:

  • DOI: https://doi.org/10.1038/s41586-026-11117-y