Influence of Magnetic Order on Proximity-Induced Superconductivity in Mn Layers on Nb(110) from First Principles
Pith reviewed 2026-05-21 16:48 UTC · model grok-4.3
The pith
Antiferromagnetic ordering in Mn layers induces an order of magnitude more singlet superconductivity from Nb than ferromagnetic ordering.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Using Bogoliubov-de Gennes calculations within the screened Korringa-Kohn-Rostoker method, the study shows that the induced singlet order parameter in the Mn layer reaches a maximum of only 4.44% of the bulk Nb value, but is an order of magnitude larger for antiferromagnetic than for ferromagnetic order. The induced triplet order parameter is negligible yet comparable in size to the singlet one, pointing to singlet-triplet mixing, while the density of states exhibits secondary gaps, plateau regions, and central V-shaped in-gap states that depend on the magnetic configuration.
What carries the argument
The Bogoliubov-de Gennes solver for superconducting heterostructures within the screened Korringa-Kohn-Rostoker Green's function method, used to compute layer-resolved singlet and triplet order parameters and spectral functions.
If this is right
- The density of states inside the Nb superconducting gap shows secondary gaps and V-shaped in-gap states whose shape depends on whether the Mn layers are ferromagnetic or antiferromagnetic.
- Single Mn layer heterostructures exhibit bands crossing the Fermi level inside the superconducting gap.
- Induced magnetic moments appear in the Nb layers adjacent to the Mn.
- The maximum induced singlet pairing in Mn remains small, at 4.44 percent of the bulk Nb value even in the favored AFM case.
Where Pith is reading between the lines
- Controlling the magnetic order in thin magnetic films could serve as a knob to tune the strength of proximity-induced superconductivity in hybrid devices.
- Similar order-dependent proximity effects may occur in other magnetic metal-superconductor interfaces and could be tested with scanning tunneling spectroscopy.
- Accounting for disorder or interface imperfections might reduce the observed differences between FM and AFM configurations.
Load-bearing premise
The calculations assume ideal, defect-free single- and double-Mn-layer geometries with purely ferromagnetic or antiferromagnetic order and that the BdG solver accurately captures the proximity-induced pairing without additional interface scattering or disorder effects.
What would settle it
Experimental measurements on real Mn/Nb(110) samples showing no significant difference in the induced superconducting order parameter or density-of-states features between ferromagnetic and antiferromagnetic Mn layers would falsify the central claim.
Figures
read the original abstract
We investigate the influence of magnetic order on the proximity-induced superconducting state in the Mn layers of a Mn-Nb(110) heterostructure by using a first-principles method. For this study, we use the recently developed Bogoliubov-de Gennes (BdG) solver for superconducting heterostructures [Csire et al., Phys. Rev. B 97, 024514 (2018)] within the first-principles calculations based on multiple scattering theory and the screened Korringa-Kohn-Rostoker (SKKR) Green's function method. In our calculations, we first study the normal-state density of states (DOS) in the single- and double-Mn-layer heterostructures, and calculate the induced magnetic moments in the Nb layers. Next, we compute the momentum-resolved spectral functions in the superconducting state for the heterostructure with a single Mn layer, and find bands crossing the Fermi level within the superconducting (SC) gap. We also study the SC state DOS in the single- and double-Mn-layer heterostructures and compare some of our results with experimental findings, revealing secondary gaps, plateau-like regions, and central V-shaped in-gap states within the bulk SC Nb gap that are magnetic-order-dependent. Finally, we compute the singlet and internally antisymmetric triplet (IAT) order parameters for each layer for both heterostructures, and find an order of magnitude difference in the induced singlet part of the SC order parameter in the Mn layer/s between the FM and AFM cases in favor of the AFM pairing with the maximum still being only 4.44% of the bulk Nb singlet order parameter value. We also find a negligible induced triplet part, yet comparable to the induced singlet values, indicating some singlet-triplet mixing in the Mn layer/s.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript uses first-principles SKKR Green's function calculations with a recently developed BdG solver to examine proximity-induced superconductivity in single- and double-Mn-layer heterostructures on Nb(110). It reports normal-state DOS and induced Nb moments, momentum-resolved spectral functions showing bands crossing the Fermi level inside the SC gap, magnetic-order-dependent features in the SC DOS (secondary gaps, plateau regions, central V-shaped in-gap states) that are compared to experiment, and layer-resolved singlet and internally antisymmetric triplet (IAT) order parameters. The central quantitative result is an order-of-magnitude larger induced singlet order parameter in the AFM configuration versus FM (maximum 4.44% of bulk Nb value) together with a small but comparable induced triplet component indicating singlet-triplet mixing.
Significance. If the reported order-parameter ratios hold, the work supplies concrete ab initio numbers for the magnetic-order dependence of proximity-induced pairing amplitudes at a transition-metal interface. The finding that AFM order favors a larger singlet component while still keeping the absolute value small, together with the evidence for limited singlet-triplet mixing, is relevant to ongoing efforts to engineer unconventional superconductivity in magnetic heterostructures. The direct comparison of calculated DOS features with experimental spectra is a positive element; the absence of machine-checked proofs or fully open reproducible code is noted but does not detract from the first-principles character of the study.
major comments (2)
- [Methods / Results on order parameters] The central claim of an order-of-magnitude difference in the induced singlet order parameter (abstract and final paragraph) rests on the layer-resolved anomalous Green's functions obtained from the BdG-SKKR solver. Without explicit convergence tests with respect to k-point sampling, number of layers, or the superconducting gap parameter, it is difficult to judge whether the reported 4.44% ratio and the FM/AFM contrast are robust or sensitive to numerical choices.
- [Introduction / Computational details] The modeling assumes ideal, defect-free single- and double-Mn-layer geometries with purely collinear FM or AFM order. While standard, this assumption is load-bearing for the headline difference; any interface disorder or non-collinear components could alter the induced amplitudes, yet no estimate of such effects is provided.
minor comments (3)
- [SC state DOS subsection] The abstract states that some SC DOS features are compared with experiment, but the main text should clarify which specific experimental spectra (energy resolution, temperature) are being matched and whether the comparison is quantitative or qualitative.
- [Order-parameter calculation] Notation for the singlet and IAT triplet order parameters should be defined explicitly (e.g., which component of the anomalous Green's function is integrated) to allow direct reproduction of the 4.44% figure.
- [Figures] Figure captions for the spectral functions and DOS plots would benefit from stating the broadening parameter and the energy window used.
Simulated Author's Rebuttal
We thank the referee for the careful reading, positive assessment of the work, and recommendation for minor revision. The comments on numerical robustness and modeling assumptions are well taken, and we address each point below with clarifications and planned revisions to strengthen the presentation.
read point-by-point responses
-
Referee: [Methods / Results on order parameters] The central claim of an order-of-magnitude difference in the induced singlet order parameter (abstract and final paragraph) rests on the layer-resolved anomalous Green's functions obtained from the BdG-SKKR solver. Without explicit convergence tests with respect to k-point sampling, number of layers, or the superconducting gap parameter, it is difficult to judge whether the reported 4.44% ratio and the FM/AFM contrast are robust or sensitive to numerical choices.
Authors: We appreciate the referee drawing attention to this aspect of the presentation. Our calculations used a 120×120 k-point mesh in the 2D Brillouin zone and 15–20 Nb layers, choices that ensure convergence of the normal-state DOS and moments to within 1%. However, we acknowledge that explicit tests focused on the superconducting order parameters were not shown in the original manuscript. In response, we have carried out additional checks: increasing the k-mesh to 200×200 changes the Mn-layer singlet order parameter by less than 4% in both FM and AFM cases; extending to 30 Nb layers alters the values by at most 3%. The superconducting gap parameter was fixed to the experimental Nb bulk value; varying it by ±10% scales the absolute amplitudes but leaves the FM/AFM ratio and the small triplet component essentially unchanged. These tests confirm that the reported order-of-magnitude difference is robust. We will add a dedicated paragraph in the Methods section together with a supplementary figure summarizing the convergence data. revision: yes
-
Referee: [Introduction / Computational details] The modeling assumes ideal, defect-free single- and double-Mn-layer geometries with purely collinear FM or AFM order. While standard, this assumption is load-bearing for the headline difference; any interface disorder or non-collinear components could alter the induced amplitudes, yet no estimate of such effects is provided.
Authors: We agree that the ideal, defect-free, collinear geometries constitute a central modeling choice. This is the conventional starting point in first-principles studies of magnetic proximity effects, allowing us to isolate the influence of magnetic order on the induced pairing amplitudes. Quantitative estimates of disorder or non-collinear magnetism would require large supercell calculations or extensions such as the coherent-potential approximation, which lie beyond the scope of the present investigation. At the same time, the magnetic-order-dependent features we obtain in the superconducting DOS (secondary gaps, plateau regions, and central V-shaped states) show reasonable correspondence with available experimental spectra, suggesting that the ideal model captures the dominant physics. In the revised manuscript we will expand the final discussion paragraph to state this limitation explicitly and to indicate that future work could address non-ideal interface effects. revision: partial
Circularity Check
No significant circularity; results follow from first-principles BdG solution
full rationale
The paper computes layer-resolved singlet and triplet order parameters by solving the Bogoliubov-de Gennes equations within the SKKR multiple-scattering framework for fixed ferromagnetic or antiferromagnetic configurations in ideal Mn/Nb(110) heterostructures. The normal-state DOS, induced moments, spectral functions, and anomalous Green's functions are obtained directly from the self-consistent electronic-structure calculation; the reported order-of-magnitude difference (AFM-favored singlet component reaching at most 4.44 % of bulk Nb) and the small triplet admixture are numerical outputs of that solution rather than quantities fitted to or defined from the same data. The cited BdG solver reference supplies the numerical method but does not enter the derivation as a load-bearing uniqueness theorem or ansatz that forces the final ratios. Experimental DOS comparisons occur after the calculations and serve as external validation, not as input that defines the predicted quantities. No self-definitional loop, fitted-input prediction, or renaming of known results appears in the chain.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption The Bogoliubov-de Gennes equations within the multiple-scattering framework accurately describe proximity-induced superconductivity in the Mn-Nb heterostructure.
- domain assumption The screened Korringa-Kohn-Rostoker Green's function method yields reliable normal-state electronic structure and magnetic moments for the heterostructure.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We compute the singlet and internally antisymmetric triplet (IAT) order parameters for each layer... order of magnitude difference in the induced singlet part... maximum still being only 4.44% of the bulk Nb singlet order parameter value.
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
solve the fully relativistic (Kohn-Sham) Dirac-Bogoliubov-de Gennes (DBdG) equations within the multiple scattering Green’s function method (i.e., the screened Korringa-Kohn-Rostoker (SKKR) method)
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
as well as in-gap states within the bulk Nb SC gap. Motivated by this experiment, we apply the recently developed first-principles-based method for SC heterostructures [32, 33] to the FM and AFM Mn single and double atomic layers on Nb, by solving the fully relativistic (Kohn-Sham) Dirac-Bogoliubov-de Gennes (DBdG) equations within the multiple scattering...
-
[2]
direction and⃗ a2 = 3.3004 ˚A along the [001] direction. The arrows represent the Mn spin moments and are pointing into the page (Red) and out of the page (Blue). (c) A schematic diagram of Mn layers on a Nb(110) substrate with the middle interface region consisting of (1-2) Mn layers, some (10-12) Nb layers, as well as several vacuum layers. converged po...
-
[3]
+f(− 1 2 , 1 2) T ±1 s {f(s, s ′)}=f(± 1 2 ,± 1 2) while||M|| F denotes the Frobenius norm of matrixM. This quantity signals the emergence of IAT pairs in any angular momentum channels and its energy integral gives a quantity analogous to the anomalous charge. Accordingly, We first calculate the singletχS component of the SC order parameter, and then calc...
-
[4]
J. Linder and J. W. A. Robinson, Superconducting spintronics, Nat. Phys.11, 307 (2015)
work page 2015
-
[5]
Eschrig, Spin-polarized supercurrents for spintronics: a review of current progress, Rep
M. Eschrig, Spin-polarized supercurrents for spintronics: a review of current progress, Rep. Prog. Phys.78, 104501 (2015)
work page 2015
-
[6]
R. Cai, I. ˇZuti´ c, and W. Han, Superconductor/ferromagnet heterostructures: A platform for superconducting spintronics and quantum computation, Adv. Quantum Technol.6, 2200080 (2023)
work page 2023
-
[7]
F. S. Bergeret, A. F. Volkov, and K. B. Efetov, Long-range proximity effects in superconductor- ferromagnet structures, Phys. Rev. Lett.86, 4096 (2001)
work page 2001
-
[8]
M. Krawiec, B. L. Gy¨ orffy, and J. F. Annett, Spontaneous spin-polarized currents in superconductor-ferromagnetic metal heterostructures, Phys. Rev. B66, 172505 (2002)
work page 2002
-
[9]
F. S. Bergeret, A. F. Volkov, and K. B. Efetov, Manifestation of triplet superconductivity in superconductor-ferromagnet structures, Phys. Rev. B68, 064513 (2003)
work page 2003
-
[10]
F. S. Bergeret, A. F. Volkov, and K. B. Efetov, Odd triplet superconductivity and related phenomena in superconductor-ferromagnet structures, Rev. Mod. Phys.77, 1321–1373 (2005)
work page 2005
-
[11]
J. W. A. Robinson, J. D. S. Witt, and M. G. Blamire, Controlled injection of spin-triplet supercurrents into a strong ferromagnet, Science329, 59 (2010)
work page 2010
-
[12]
D. Fritsch and J. F. Annett, Proximity effect in superconductor/conical magnet heterostruc- tures, J. Phys. Condens. Matter26, 274212 (2014)
work page 2014
-
[13]
D. Fritsch and J. F. Annett, Triplet superconductivity and proximity effect induced by bloch 19 and n´ eel domain walls, Supercond. Sci. Technol.28, 085015 (2015)
work page 2015
- [14]
- [15]
-
[16]
E. A. Demler, G. B. Arnold, and M. R. Beasley, Superconducting proximity effects in magnetic metals, Phys. Rev. B55, 15174 (1997)
work page 1997
-
[17]
A. I. Buzdin, Proximity effects in superconductor-ferromagnet heterostructures, Rev. Mod. Phys.77, 935–976 (2005)
work page 2005
-
[18]
K. Halterman and O. T. Valls, Proximity effects at ferromagnet-superconductor interfaces, Phys. Rev. B65, 014509 (2001)
work page 2001
- [19]
-
[20]
C.-T. Wu, O. T. Valls, and K. Halterman, Proximity effects in conical- ferromagnet/superconductor bilayers, Phys. Rev. B86, 184517 (2012)
work page 2012
-
[21]
A. Golubov, S. Bakurskiy, M. Kupriyanov, T. Karabassov, A. Vasenko, and A. Sidorenko, The physics of superconductor-ferromagnet hybrid structures (2025)
work page 2025
-
[22]
J. S. Jiang, D. Davidovi´ c, D. H. Reich, and C. L. Chien, Oscillatory superconducting transition temperature in nb/gd multilayers, Phys. Rev. Lett.74, 314 (1995)
work page 1995
-
[23]
Y. N. Khaydukov, A. S. Vasenko, E. A. Kravtsov, V. V. Progliado, V. D. Zhaketov, A. Csik, Y. V. Nikitenko, A. V. Petrenko, T. Keller, A. A. Golubov, M. Y. Kupriyanov, V. V. Ustinov, V. L. Aksenov, and B. Keimer, Magnetic and superconducting phase diagram of nb/gd/nb trilayers, Phys. Rev. B97, 144511 (2018)
work page 2018
-
[24]
H. Yamazaki, N. Shannon, and H. Takagi, Interplay between superconductivity and ferromag- netism in epitaxial nb(110)/au(111)/fe(110) trilayers, Phys. Rev. B73, 094507 (2006)
work page 2006
-
[25]
S. Reymond, P. SanGiorgio, M. R. Beasley, J. Kim, T. Kim, and K. Char, Tunneling density of states as a function of thickness in superconductor/strong ferromagnet bilayers, Phys. Rev. B73, 054505 (2006)
work page 2006
-
[26]
C. Bell, E. J. Tarte, G. Burnell, C. W. Leung, D.-J. Kang, and M. G. Blamire, Proximity and josephson effects in superconductor/antiferromagnetic Nb/γ−fe 50mn50 heterostructures, Phys. Rev. B68, 144517 (2003). 20
work page 2003
-
[27]
M. H¨ ubener, D. Tikhonov, I. A. Garifullin, K. Westerholt, and H. Zabel, The antiferromag- net/superconductor proximity effect in cr/v/cr trilayers, J. Phys. Condens. Matter14, 8687 (2002)
work page 2002
-
[28]
A. Mani, T. Geetha Kumary, and J. G. Lin, Thickness controlled proximity effects in C-type antiferromagnet/superconductor heterostructure, Sci. Rep.5, 12780 (2015)
work page 2015
-
[29]
S. Chourasia, L. J. Kamra, I. V. Bobkova, and A. Kamra, Generation of spin-triplet cooper pairs via a canted antiferromagnet, Phys. Rev. B108, 064515 (2023)
work page 2023
-
[30]
E. H. Fyhn, A. Brataas, A. Qaiumzadeh, and J. Linder, Superconducting proximity effect and long-ranged triplets in dirty metallic antiferromagnets, Phys. Rev. Lett.131, 076001 (2023)
work page 2023
-
[31]
G. A. Bobkov, V. M. Gordeeva, A. M. Bobkov, and I. V. Bobkova, Oscillatory superconduct- ing transition temperature in superconductor/antiferromagnet heterostructures, Phys. Rev. B 108, 184509 (2023)
work page 2023
-
[32]
G. A. Bobkov, I. V. Bobkova, and A. M. Bobkov, Proximity effect in superconduc- tor/antiferromagnet hybrids: N´ eel triplets and impurity suppression of superconductivity, Phys. Rev. B108, 054510 (2023)
work page 2023
-
[33]
G. A. Bobkov, V. A. Bobkov, I. V. Bobkova, A. M. Bobkov, and A. A. Golubov, Triplet correlations in superconductor/antiferromagnet heterostructures: Dependence on type of an- tiferromagnetic ordering, J. Supercond. Nov. Magn.38, 239 (2025)
work page 2025
-
[34]
R. Lo Conte, M. Bazarnik, K. Palot´ as, L. R´ ozsa, L. Szunyogh, A. Kubetzka, K. von Bergmann, and R. Wiesendanger, Coexistence of antiferromagnetism and superconductivity in mn/nb(110), Phys. Rev. B105, L100406 (2022)
work page 2022
- [35]
- [36]
-
[37]
L. N. Oliveira, E. K. U. Gross, and W. Kohn, Density-functional theory for superconductors, Phys. Rev. Lett.60, 2430 (1988)
work page 1988
-
[38]
D. M. Ceperley and B. J. Alder, Ground state of the electron gas by a stochastic method, Phys. Rev. Lett.45, 566 (1980)
work page 1980
-
[39]
J. P. Perdew and A. Zunger, Self-interaction correction to density-functional approximations 21 for many-electron systems, Phys. Rev. B23, 5048 (1981)
work page 1981
-
[40]
J. P. Carbotte, Properties of boson-exchange superconductors, Rev. Mod. Phys.62, 1027 (1990)
work page 1990
-
[41]
S. L. Cunningham, Special points in the two-dimensional Brillouin zone, Phys. Rev. B10, 4988 (1974)
work page 1974
-
[42]
K. Park, L. Oroszlany, L. Szunyogh, and B. Ujfalussy, Proximity-induced superconductivity in ferromagnetic gd layers on nb from a first-principles LDA +ustudy, Phys. Rev. B110, 174510 (2024)
work page 2024
-
[43]
A. Di Bernardo, S. Diesch, Y. Gu, J. Linder, G. Divitini, C. Ducati, E. Scheer, M. Blamire, and J. Robinson, Signature of magnetic-dependent gapless odd frequency states at supercon- ductor/ferromagnet interfaces, Nat. Comm.6, 8053 (2015)
work page 2015
-
[44]
B. Ny´ ari, A. L´ aszl´ offy, G. Csire, L. Szunyogh, and B.´Ujfalussy, Topological superconductivity from first principles. i. shiba band structure and topological edge states of artificial spin chains, Phys. Rev. B108, 134512 (2023)
work page 2023
- [45]
-
[46]
S. K. Ghosh, G. Csire, P. Whittlesea, J. F. Annett, M. Gradhand, B. ´Ujfalussy, and J. Quin- tanilla, Quantitative theory of triplet pairing in the unconventional superconductor LaNiGa 2, Phys. Rev. B101, 100506 (2020)
work page 2020
- [47]
-
[48]
J. Cayao and A. M. Black-Schaffer, Distinguishing trivial and topological zero-energy states in long nanowire junctions, Phys. Rev. B104, L020501 (2021). Appendix A: Bloch Spectral F unction (BSF) The BSF can be resolved into layer as well as electron-like and hole-like components, which may give more insight into the proximity effect in the systems at ha...
work page 2021
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.