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
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
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
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
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
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
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
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
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
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
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).
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
Rusinov, A. Theory of gapless superconductivity in alloys containing paramagnetic impurities. Sov. Phys. JETP 29, 1101–1106 (1969).
ADS Google Scholar
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
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
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
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
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
Vaxevani, K. et al. Inducing superconductivity in atomically thin Au(111) films via interface engineering. Zenodo https://doi.org/10.5281/zenodo.22815538 (2026).
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
Landé, A. Termstruktur und Zeemaneffekt der Multipletts. Z. Phys. 15, 189–205 https://doi.org/10.1007/BF01330473 (1923).
Article ADS Google Scholar
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
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).
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.
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 )
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.
Reprints and permissions
About this article
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