pith. sign in

arxiv: 2606.25862 · v1 · pith:ZSXLWZO2new · submitted 2026-06-24 · ❄️ cond-mat.supr-con · physics.comp-ph

First-Principles Quantum-Spectral framework for Elementary Vortex Pinning in superconductors

Haozhe Shi , Yuncheng Xie , Tong Zhang , Weibin Chu , Xin-Gao Gong This is my paper

Pith reviewed 2026-06-25 19:31 UTC · model grok-4.3

classification ❄️ cond-mat.supr-con physics.comp-ph
keywords vortex pinningBogoliubov-de Gennesfirst-principlesCdGM statespoint defectsFeSesuperconductivitycritical current
0
0 comments X

The pith

A defect reorganizes vortex-core states to generate the elementary pinning force in type-II superconductors.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper builds a computational method that turns the measured reorganization of vortex-core electronic states into a calculated pinning energy and force. It embeds atomically detailed defect structures from density-functional theory into a finite-box Bogoliubov-de Gennes free-energy calculation and isolates the pinning contribution through a four-configuration subtraction. Applied to FeSe, the method recovers both the pinning strength and the spectral pattern seen in scanning tunneling microscopy for the iron-site vacancy. The same approach shows attractive pinning for all five point defects examined in FeSe and FeTe, with the tellurium vacancy in FeTe strongest. If the framework is correct, elementary pinning becomes a quantity that can be computed directly from electronic structure rather than fitted after measurement.

Core claim

The defect-induced reorganization of the vortex-core spectrum is the microscopic origin of the elementary pinning force; the framework reproduces the microscopic STM value for the FeSe Fe-site vacancy together with the measured spectral reorganization. The force is evaluated as a finite-box vortex-insertion free energy whose four-configuration subtraction isolates the meV-scale interaction from much larger backgrounds. With the superconducting gap scale and vortex-core profile fixed from experiments, the FeSe Fe-site vacancy reproduces the microscopic STM value together with the measured spectral reorganization. All five point defects in FeSe and FeTe pin attractively, with FeTe Te-site vaca

What carries the argument

finite-box projected Bogoliubov-de Gennes free-energy formalism with four-configuration subtraction that converts defect-induced quasiparticle spectral reorganization into pinning energies

If this is right

  • The framework reproduces both the pinning energy and the spectral changes measured by STM for the Fe-site vacancy in FeSe.
  • All five point defects studied in FeSe and FeTe produce attractive pinning forces.
  • The Te-site vacancy in FeTe yields the largest pinning energy among the defects examined.
  • Elementary vortex pinning can be treated as a first-principles electronic-structure quantity.
  • Point-defect screening for higher critical currents becomes feasible by direct computation.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same embedding procedure could be used to rank candidate defects in other type-II materials before synthesis.
  • Connecting the computed elementary forces to macroscopic critical-current models would test whether single-defect pinning dominates in real samples.
  • The framework might be extended to calculate pinning by line defects or grain boundaries by enlarging the simulation cell while keeping the spectral subtraction.

Load-bearing premise

The superconducting gap scale and vortex-core profile are taken as fixed inputs from experiment, and the four-configuration free-energy subtraction cleanly isolates the meV-scale pinning interaction from much larger background energies without residual contamination.

What would settle it

A measured pinning energy or force for the Te-site vacancy in FeTe that differs substantially from the value computed by the framework would show that spectral reorganization does not account for the pinning force.

read the original abstract

The critical current of a type-II superconductor is controlled by vortex pinning, whose microscopic input is the elementary pinning force. Scanning tunneling spectroscopy has shown that a defect pins a vortex by reorganizing the Caroli-de Gennes-Matricon (CdGM) states in its core, but why this spectral reorganization amounts to a pinning force has lacked a quantum-mechanical, first-principles account. Here we establish a transferable first-principles computational framework for elementary vortex pinning, in which defect-resolved DFT/Wannier electronic structures are embedded into a finite-box projected Bogoliubov-de Gennes free-energy formalism to convert quasiparticle spectral reorganization into vortex-pinning energies and forces. Using this framework, we confirm that the defect-induced reorganization of the vortex-core spectrum is the microscopic origin of the elementary pinning force. The force is evaluated as a finite-box vortex-insertion free energy whose four-configuration subtraction isolates the meV-scale interaction from much larger backgrounds. With the superconducting gap scale and vortex-core profile fixed from experiments, the FeSe Fe-site vacancy reproduces the microscopic STM value together with the measured spectral reorganization. All five point defects in FeSe and FeTe pin attractively, with FeTe Te-site vacancy strongest. Elementary vortex pinning thereby becomes a computable electronic-structure quantity, opening the first-principles screening of point defects toward higher critical currents.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 3 minor

Summary. The manuscript develops a first-principles framework that embeds defect-resolved DFT/Wannier electronic structures into a finite-box projected Bogoliubov-de Gennes free-energy formalism. It converts defect-induced reorganization of Caroli-de Gennes-Matricon states into elementary vortex-pinning energies via a four-configuration subtraction (vortex+defect, vortex-only, defect-only, pristine). With the gap magnitude and vortex-core profile fixed from experiment, the method is reported to reproduce the STM pinning energy and spectral changes for the Fe-site vacancy in FeSe, and to find attractive pinning for all five point defects examined in FeSe and FeTe (strongest for Te-site vacancy in FeTe).

Significance. If the numerical isolation of the pinning energy is robust, the work supplies a transferable route from electronic-structure calculations to the microscopic pinning force that limits critical current in type-II superconductors. It directly links quasiparticle spectral reorganization to pinning energetics and enables computational screening of point defects. The approach is notable for its attempt to keep the framework first-principles once the experimental gap and core profile are supplied as inputs.

major comments (2)
  1. [methods / four-configuration subtraction] Four-configuration subtraction (methods section describing the finite-box BdG free-energy difference): the central numerical claim that this procedure isolates a clean meV-scale pinning interaction from orders-of-magnitude larger background terms is load-bearing. The manuscript does not report explicit null-tests (e.g., subtraction on the pristine system) or systematic box-size convergence to bound residual contamination at the meV level; without such quantification the reported reproduction of the STM value cannot be verified as free of cancellation artifacts.
  2. [results for FeSe vacancy] § on numerical results for FeSe Fe-site vacancy: the reproduction of the STM pinning value is presented as confirmation that CdGM reorganization is the microscopic origin, yet the sensitivity of the extracted energy to the fixed experimental inputs (gap scale and core profile) is not quantified. A modest variation in either input could shift the difference by an amount comparable to the reported meV pinning energy, weakening the claim that the result is determined by the defect-induced spectral change alone.
minor comments (3)
  1. [methods] Notation for the projected BdG Hamiltonian: an explicit equation defining the projection of the DFT/Wannier bands onto the finite box would improve reproducibility.
  2. [figures] Figure captions (spectral plots): the energy scale and broadening used for the density of states should be stated explicitly so that the reported reorganization can be compared directly with STM data.
  3. [introduction] Reference list: prior STM work on CdGM states in FeSe is cited, but an additional reference to earlier theoretical treatments of vortex pinning energetics would help situate the four-configuration approach.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive report and the recognition of the framework's potential. We address the two major comments below and will revise the manuscript accordingly to strengthen the numerical validation.

read point-by-point responses

  1. Referee: [methods / four-configuration subtraction] Four-configuration subtraction (methods section describing the finite-box BdG free-energy difference): the central numerical claim that this procedure isolates a clean meV-scale pinning interaction from orders-of-magnitude larger background terms is load-bearing. The manuscript does not report explicit null-tests (e.g., subtraction on the pristine system) or systematic box-size convergence to bound residual contamination at the meV level; without such quantification the reported reproduction of the STM value cannot be verified as free of cancellation artifacts.

    Authors: We agree that explicit null-tests and box-size convergence tests are necessary to rigorously bound possible cancellation artifacts. In the revised manuscript we will add (i) a null-test subtraction performed on the pristine supercell (expected to return zero within numerical tolerance) and (ii) a systematic study of pinning energy versus box linear size (up to the largest computationally feasible supercells), demonstrating that the extracted meV-scale value stabilizes to within <0.1 meV. These additions will be placed in the Methods section and a new supplementary figure. revision: yes

  2. Referee: [results for FeSe vacancy] § on numerical results for FeSe Fe-site vacancy: the reproduction of the STM pinning value is presented as confirmation that CdGM reorganization is the microscopic origin, yet the sensitivity of the extracted energy to the fixed experimental inputs (gap scale and core profile) is not quantified. A modest variation in either input could shift the difference by an amount comparable to the reported meV pinning energy, weakening the claim that the result is determined by the defect-induced spectral change alone.

    Authors: We acknowledge that a quantitative sensitivity analysis to the experimental inputs was not provided. In revision we will add a supplementary section that varies the gap magnitude by ±10 % around the experimental value and perturbs the core-profile parameters within the range consistent with STM data. The resulting spread in the computed pinning energy will be shown to remain well below the reported STM value, thereby confirming that the dominant contribution indeed originates from the defect-induced CdGM reorganization rather than from the fixed inputs. revision: yes

Circularity Check

0 steps flagged

No significant circularity in the derivation chain

full rationale

The paper computes the pinning energy explicitly as a four-configuration finite-box free-energy difference (vortex+defect minus vortex minus defect minus pristine) after embedding DFT/Wannier bands into a projected BdG Hamiltonian. The gap magnitude and vortex-core profile are fixed external experimental inputs, and the central numerical claim is a direct match to an independent STM measurement for the FeSe Fe-site vacancy. No step reduces by construction to a fitted parameter renamed as prediction, no self-citation chain justifies a uniqueness theorem or ansatz, and the spectral-reorganization origin is shown via the explicit energy subtraction rather than by definition. The framework is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

Only the abstract is available; the ledger is therefore limited to elements explicitly named. The gap scale and vortex-core profile are taken from experiment and function as external inputs. The BdG formalism and Wannier projection are standard methods.

free parameters (2)
  • superconducting gap scale
    Fixed from experiments as stated in the abstract; used as input to the BdG calculation.
  • vortex-core profile
    Fixed from experiments as stated in the abstract; used to set the spatial form inside the finite-box model.
axioms (2)
  • standard math Standard Bogoliubov-de Gennes mean-field equations govern the quasiparticle spectrum.
    Invoked when the DFT/Wannier structures are embedded into the projected BdG free-energy formalism.
  • domain assumption Wannier projection accurately represents the defect-resolved electronic structure for subsequent BdG embedding.
    Required for the transferability claim of the framework.

pith-pipeline@v0.9.1-grok · 5786 in / 1480 out tokens · 35375 ms · 2026-06-25T19:31:03.295873+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

36 extracted references

  1. [1]

    Larbalestier, D., Gurevich, A., Feldmann, D. M. & Polyanskii, A. High-T c superconducting materials for electric power applications.Nature414, 368–377 (2001)

    2001

  2. [2]

    Abrikosov, A. A. On the magnetic properties of superconductors of the second group.Sov. Phys. JETP5, 1174–1182 (1957). URL https://jetp.ras.ru/cgi-bin/e/ index/e/5/6/p1174?a=list. Russian original: Zh. Eksp. Teor. Fiz. 32, 1442 (1957)

    1957

  3. [3]

    Bean, C. P. Magnetization of high-field superconductors.Rev. Mod. Phys.36, 31–39 (1964)

    1964

  4. [4]

    Anderson, P. W. & Kim, Y. B. Hard superconductivity: Theory of the motion of abrikosov flux lines.Rev. Mod. Phys.36, 39–43 (1964)

    1964

  5. [5]

    Flux pinning mechanisms in type ii superconductors.Philos

    Dew-Hughes, D. Flux pinning mechanisms in type ii superconductors.Philos. Mag.30, 293–305 (1974)

    1974

  6. [6]

    Larkin, A. I. & Ovchinnikov, Y. N. Pinning in type ii superconductors.J. Low Temp. Phys.34, 409–428 (1979)

    1979

  7. [7]

    V., Geshkenbein, V

    Blatter, G., Feigel’man, M. V., Geshkenbein, V. B., Larkin, A. I. & Vinokur, V. M. Vortices in high-temperature superconductors.Rev. Mod. Phys.66, 1125–1388 (1994)

    1994

  8. [8]

    V., Kurkij¨ arvi, J

    Thuneberg, E. V., Kurkij¨ arvi, J. & Rainer, D. Elementary-flux-pinning potential in type-ii superconductors.Phys. Rev. B29, 3913–3923 (1984)

    1984

  9. [9]

    B., Finnemore, D

    Hyun, O. B., Finnemore, D. K., Schwartzkopf, L. & Clem, J. R. Elementary pinning force for a superconducting vortex.Phys. Rev. Lett.58, 599–601 (1987)

    1987

  10. [10]

    Caroli, C., de Gennes, P. G. & Matricon, J. Bound fermion states on a vortex line in a type ii superconductor.Phys. Lett.9, 307–309 (1964)

    1964

  11. [11]

    F., Robinson, R

    Hess, H. F., Robinson, R. B., Dynes, R. C., Valles, J. M., Jr. & Waszczak, J. V. Scanning-tunneling-microscope observation of the abrikosov flux lattice and the density of states near and inside a fluxoid.Phys. Rev. Lett.62, 214–216 (1989). 15

    1989

  12. [12]

    Chen, C.et al.Observation of discrete conventional caroli–de gennes–matricon states in the vortex core of single-layer FeSe/SrTiO3.Phys. Rev. Lett.124, 097001 (2020)

    2020

  13. [13]

    Z.et al.Observation of distinct spatial distributions of the zero and nonzero energy vortex modes in (Li 0.84Fe0.16)OHFeSe.Phys

    Zhang, T. Z.et al.Observation of distinct spatial distributions of the zero and nonzero energy vortex modes in (Li 0.84Fe0.16)OHFeSe.Phys. Rev. Lett.126, 127001 (2021)

    2021

  14. [14]

    Z.et al.Phase shift and magnetic anisotropy induced field splitting of impurity states in (Li 1−xFex)OHFeSe superconductor.Phys

    Zhang, T. Z.et al.Phase shift and magnetic anisotropy induced field splitting of impurity states in (Li 1−xFex)OHFeSe superconductor.Phys. Rev. Lett.130, 206001 (2023)

    2023

  15. [15]

    Chen, C.et al.Revealing the microscopic mechanism of elementary vortex pinning in superconductors.Phys. Rev. X14, 041039 (2024)

    2024

  16. [16]

    G.Superconductivity of Metals and Alloys(W

    de Gennes, P. G.Superconductivity of Metals and Alloys(W. A. Benjamin, New York, 1966)

    1966

  17. [17]

    & Schl¨ uter, M

    Gygi, F. & Schl¨ uter, M. Self-consistent electronic structure of a vortex line in a type-ii superconductor.Phys. Rev. B43, 7609–7621 (1991)

    1991

  18. [18]

    & Hosono, H

    Kamihara, Y., Watanabe, T., Hirano, M. & Hosono, H. Iron-based layered super- conductor La[O1−xFx]FeAs (x= 0.05–0.12) withT c = 26 K.J. Am. Chem. Soc. 130, 3296–3297 (2008)

    2008

  19. [19]

    Huang, Y.et al.Superconducting (Li,Fe)OHFeSe film of high quality and high critical parameters.Chin. Phys. Lett.34, 077404 (2017)

    2017

  20. [20]

    Nano Lett.21, 1327–1334 (2021)

    Qin, H.et al.Moir´ e superlattice-induced superconductivity in one-unit-cell FeTe. Nano Lett.21, 1327–1334 (2021)

    2021

  21. [21]

    Ren, W.et al.Oxygen adsorption induced superconductivity in ultrathin FeTe film on SrTiO3(001).Materials14, 4584 (2021)

    2021

  22. [22]

    Yan, Z.-J.et al.Stoichiometric FeTe is a superconductor.Nature652, 342–348 (2026)

    2026

  23. [23]

    Bardeen, J., Cooper, L. N. & Schrieffer, J. R. Theory of superconductivity.Phys. Rev.108, 1175–1204 (1957)

    1957

  24. [24]

    Gor’kov, L. P. Microscopic derivation of the Ginzburg–Landau equations in the theory of superconductivity.Sov. Phys. JETP9, 1364–1367 (1959). URL http: //jetp.ras.ru/cgi-bin/e/index/e/9/6/p1364?a=list

    1959

  25. [25]

    & Sham, L

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

    1965

  26. [26]

    P., Burke, K

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

    1996

  27. [27]

    & Furthm¨ uller, J

    Kresse, G. & Furthm¨ uller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set.Phys. Rev. B54, 11169–11186 (1996)

    1996

  28. [28]

    & Vanderbilt, D

    Marzari, N. & Vanderbilt, D. Maximally localized generalized Wannier functions for composite energy bands.Phys. Rev. B56, 12847–12865 (1997)

    1997

  29. [29]

    A.et al.wannier90: A tool for obtaining maximally-localised Wannier functions.Comput

    Mostofi, A. A.et al.wannier90: A tool for obtaining maximally-localised Wannier functions.Comput. Phys. Commun.178, 685–699 (2008)

    2008

  30. [30]

    A., Yates, J

    Marzari, N., Mostofi, A. A., Yates, J. R., Souza, I. & Vanderbilt, D. Maximally localized Wannier functions: Theory and applications.Rev. Mod. Phys.84, 1419– 1475 (2012)

    2012

  31. [31]

    Sch¨ onemann, P. H. A generalized solution of the orthogonal Procrustes problem. Psychometrika31, 1–10 (1966)

    1966

  32. [32]

    Kuhn, H. W. The Hungarian method for the assignment problem.Nav. Res. Logist. Q.2, 83–97 (1955)

    1955

  33. [33]

    Bl¨ ochl, P. E. Projector augmented-wave method.Phys. Rev. B50, 17953–17979 (1994)

    1994

  34. [34]

    & Joubert, D

    Kresse, G. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method.Phys. Rev. B59, 1758–1775 (1999)

    1999

  35. [35]

    Clem, J. R. Simple model for the vortex core in a type II superconductor.J. Low Temp. Phys.18, 427–434 (1975)

    1975

  36. [36]

    N., Gross, E

    Oliveira, L. N., Gross, E. K. U. & Kohn, W. Density-functional theory for superconductors.Phys. Rev. Lett.60, 2430–2433 (1988). 17

    1988