REVIEW 3 major objections 3 minor 44 references
Magnetic Field Dependence of the Spin Susceptibility on Conventional s-wave Superconductor LaRu$_4$P$_{12}$ Revealed by $^{31}$P-NMR and $^{139}$La-NMR
T0 review · 3 major / 3 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read In the conventional s-wave superconductor LaRu4P12, the superconducting-state spin susceptibility is proportional to the applied magnetic field and connects smoothly to the normal-state value at the upper critical field.
desk verdict Direct NMR confirmation of H-linear spin susceptibility in an orbital-limited s-wave superconductor; clean two-site subtraction, with minor caveats about the K_dia assumption and fit range. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the two-site Knight-shift difference. Because the diamagnetic shift $K_{\rm dia}$ is assumed to be identical at the $^{31}$P and $^{139}$La sites, the difference $\Delta K^{\rm La} - \Delta K^{\rm P}$ cancels it, leaving $(A_{\rm hf}^{\rm La} - A_{\rm hf}^{\rm P})\Delta\chi_{\rm spin}$. The hyperfine ratio $A_{\rm hf}^{\rm P}/A_{\rm hf}^{\rm La} = 0.60$ is read off the linear $K^{\rm P}$-versus-$K^{\rm La}$ plot from 10 K to 220 K, converting the two-site difference into the spin shift at each site. This subtraction is what lets the authors extract a field-linear spin signal from spectra that otherwise contain vortex-lattice broadening (the Redfield pattern, an asymmetric line shape from the vortex lattice) and a large diamagnetic shift at low fields.
What would settle it
A direct calculation of the local magnetic field at the phosphorus and lanthanum sites in the vortex lattice, using the reported $\kappa = 24.7$ and $H_{c2} = 3.2$ T, would show whether the diamagnetic Knight shift is site independent; if it is not, the two-site subtraction does not cancel it and the reported linear spin susceptibility would be an artifact. Comparing the extraction with a third nuclear site would test this directly.
Extended reading notes
Core claim
The central claim is that the spin susceptibility in the superconducting state of LaRu$_4$P$_{12}$ obeys $\chi_{\rm spin} = \alpha(H - H_{c2}) + \chi_{\rm normal}$, with a linear recovery in field that joins the normal-state susceptibility at $H_{c2}$. The evidence comes from the $^{31}$P and $^{139}$La Knight shifts at 1.4 K: measured against the 10 K normal-state values, the two-site difference removes the common diamagnetic shift, and the extracted $\Delta K_{\rm spin}$ increases linearly with $H$ up to $\mu_0 H_{c2} = 3.2$ T. The authors interpret this linear field dependence as the expected consequence of the quasiparticle density of states being proportional to $H/H_{c2}$ in the mixed state, realized when the upper critical field is governed by orbital pair-breaking rather than Pauli limiting.
Load-bearing premise
The analysis assumes that the magnetic-field shift caused by superconducting screening currents is exactly the same at the 31P and 139La sites, so subtracting the two-site difference removes it completely; if the vortex lattice shifts one site more than the other, the extracted spin susceptibility is systematically wrong.
Editorial extensions
If this is right
- In LaRu$_4$P$_{12}$, the superconducting-state spin susceptibility follows $\chi_{\rm spin} = \alpha(H - H_{c2}) + \chi_{\rm normal}$, so the spin response recovers linearly with field and matches the normal-state value at $H_{c2}$.
- The smooth recovery means no Pauli-paramagnetic anomaly or Fulde-Ferrell-Larkin-Ovchinnikov-type spin response appears near the upper critical field, consistent with orbital pair-breaking dominance.
- The two-site subtraction recipe can be applied to other multi-site conventional superconductors to isolate the spin susceptibility from the superconducting diamagnetic background.
- The fitted Ginzburg-Landau parameter $\kappa \approx 24.7$ and estimated lower critical field $\mu_0 H_{c1} \approx 10$ mT place LaRu$_4$P$_{12}$ in the strong type-II regime.
- The result is consistent with the earlier $1/T_1T \propto H^2$ observation, reinforcing the orbital-pair-breaking picture from a different observable.
Reading between the lines
- If the two-site cancellation of the diamagnetic shift is exact, the same linear $\chi_{\rm spin}(H)$ should be seen in other orbital-limited s-wave superconductors that have two NMR-accessible sites; the two-site subtraction could become a quick diagnostic of orbital pair-breaking.
- The data deviate from linearity below 1 T where the diamagnetic term dominates, so the true low-field behavior remains unresolved; measurements closer to the lower critical field with oriented crystals could reveal a vortex-lattice correction hidden by the current error bars.
- The assumed proportionality between $\chi_{\rm spin}$ and the quasiparticle density of states could be tested by comparing the NMR-derived slope $\alpha$ with the specific-heat-derived $D(E_F)$ on the same crystals; a mismatch would indicate corrections beyond the simple single-particle picture.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports 31P and 139La NMR Knight-shift measurements on the conventional s-wave superconductor LaRu4P12 in the mixed state at 1.4 K. By taking the difference of the Knight shifts at the two nuclear sites, the authors argue that the superconducting diamagnetic shift Kdia cancels, allowing them to extract the spin part of the shift ΔKspin as a function of magnetic field. The extracted ΔKspin is reported to increase linearly with H and to connect smoothly to the normal-state value at Hc2, leading to the conclusion that the superconducting-state spin susceptibility follows χspin = α(H - Hc2) + χnormal, as expected when the upper critical field is governed by orbital pair breaking. The authors also estimate Kdia and fit it with a Ginzburg-Landau expression to obtain κ ≈ 24.7. The central claim is that this is a textbook demonstration of orbital-limited superconductivity with a linear field recovery of the spin susceptibility.
Significance. If the central claim holds, the paper provides a direct, two-site NMR measurement of the field-dependent spin susceptibility in a conventional s-wave superconductor, complementing earlier 1/T1T data on the same material and giving a clean experimental confirmation of the orbital pair-breaking scenario. The two-site subtraction technique is a sound and potentially reusable approach, and the hyperfine-coupling ratio used in the analysis is calibrated from independent 10–220 K data rather than fitted to the superconducting-state data. Error bars are provided for the reported Knight shifts. The main weaknesses are the reliance on an unquantified assumption that Kdia is identical at the 31P and 139La sites, the exclusion of the two lowest-field data points from the linear fit, and the absence of a finite-temperature discussion for the comparison with the T = 0 theory. These issues are addressable but require additional analysis and clarification.
major comments (3)
- [Eq. (4)] The central subtraction in Eq. (4) cancels Kdia only if the diamagnetic Knight shift is strictly identical at the 31P and 139La sites. The manuscript justifies this by the statement 'Since Kdia is a bulk effect, it is reasonable to assume that Kdia works at the 31P and 139La sites in the same manner' (p. 3). This assumption is load-bearing: a site-dependent residual Kdia_La - Kdia_P that varies with H would directly bias the extracted Δχspin and could produce or distort the claimed linear H dependence. No quantitative argument or experimental check is provided. In particular, the Kdia values in Fig. 6 are obtained from Eq. (2) using the same two-site subtraction, so the fit yielding κ = 24.7 is not an independent validation of the cancellation. Please add a quantitative justification for site-independent Kdia (for example, that the vortex-lattice field varies on scales λ and ξ which are orders of magnitude larger than the unit cell, so the macroscopic diamagnetic field is the same at both sites) or estimate the maximum possible field-dependent residual from the data and propagate it into the systematic uncertainty of ΔKspin.
- [Fig. 5] The linear fit in Fig. 5 excludes the µ0H = 0.5 T and 1.1 T points, which the text says 'deviate from this behavior with large error bar.' Because Kdia is largest at low fields (Fig. 6), this is exactly the field region where a failure of the Kdia-cancellation assumption would be most visible; excluding these points prevents the data from falsifying the cancellation. Please report the fit details (number of points, slope, intercept, reduced χ²), perform the fit with all points weighted by their errors, and show the residuals. If the low-field points are excluded, provide a quantitative justification beyond the large error bars, and state explicitly that the claimed proportionality is restricted to the fitted field range rather than the full 0 < H < Hc2 interval.
- [Conclusions] The linear relation χspin = α(H - Hc2) + χnormal is a zero-temperature expectation: the quasiparticle density of states in the mixed state of a type-II superconductor is proportional to H/Hc2 at T = 0. The measurements are performed at 1.4 K, and for fields close to Hc2 the superconducting transition temperature at that field is only slightly above 1.4 K (e.g., from Fig. 1(c), Tc(3 T) is only about 1.8 K), so thermal quasiparticle contributions are not negligible. Please either quantify the finite-temperature correction to the linear form in the field range used for the fit, or soften the conclusion to state that the data are consistent with the T = 0 linear form within the present experimental accuracy.
minor comments (3)
- [p. 2] In the sentence 'Magnetic field calibration was carried out using 63Cu and 65Cu signal arising from the NMR coil', the word 'signal' should be plural ('signals').
- [Fig. 5] The vertical axis label in Fig. 5 is '-ΔKspin (%)' while the text and caption refer to 'ΔKspin'. Please define the sign convention explicitly, because the description 'the value increases linearly with increasing H' is confusing if the plotted quantity is negative.
- [p. 3] The sentence 'If one extrapolates ΔKspin to zero magnetic field, the resulting value is nearly equal to K (10 K)' is not clear: K(10 K) is the absolute Knight shift at one site, while ΔKspin is a difference. Please restate this as an extrapolation of Kspin(1.4 K, H → 0) to zero or to a value consistent with the spin susceptibility vanishing at H = 0.
Circularity Check
No significant circularity: the linear field recovery of the spin susceptibility is an empirical result inferred from independently calibrated Knight-shift differences, not a consequence of the fitting or self-citation chain.
full rationale
The paper's central claim, that the superconducting-state spin susceptibility recovers linearly with magnetic field and connects smoothly to the normal-state value, is obtained from measured 31P and 139La NMR Knight shifts. The crucial two-site subtraction in Eq. (4) assumes, rather than derives, that the diamagnetic Knight shift Kdia is identical at both sites; this is an unverified assumption and a possible systematic error, but it is not circular because it does not encode the target linear H dependence. The hyperfine ratio A_P/A_La = 0.60 used in Eq. (6) is calibrated independently from the normal-state K_P versus K_La relation between 10 K and 220 K (Fig. 4), outside the superconducting state and outside the field range of the claimed effect. The linear dependence in Fig. 5 is read from the data and compared with the theoretical expectation D(E_F) ∝ H/H_c2 and with the prior 1/T1T result; it is not manufactured by the fitting parameters. The Ginzburg-Landau parameter κ = 24.7 obtained from fitting Eq. (7) to Kdia is a consistency check and is not fed back into the extraction of ΔKspin. Self-citations to refs. 34-36 and 39 supply prior spectra, Hc2 data, and a related compound comparison, but the central conclusion does not reduce to these citations. No derivation step is equivalent by construction to its own input, so there is no circularity.
Assumptions & free parameters
free parameters (2)
- Hyperfine coupling ratio A_P_hf / A_La_hf =
0.60
- Ginzburg-Landau parameter κ =
24.7
assumptions (4)
- domain assumption Normal-state Knight shift is independent of magnetic field.
- domain assumption The diamagnetic Knight shift Kdia is identical at the 31P and 139La sites.
- domain assumption The hyperfine coupling ratio measured between 10 K and 220 K remains valid at 1.4 K.
- domain assumption The peak position of the broadened vortex-state NMR spectrum is a valid measure of the average Knight shift needed for the analysis.
Cite this review
Pith. "Pith review of Magnetic Field Dependence of the Spin Susceptibility on Conventional s-wave Superconductor LaRu$_4$P$_{12}$ Revealed by $^{31}$P-NMR and $^{139}$La-NMR." pith.science (2026). https://pith.science/paper/PVFL5YJD
@misc{pith2026250520702,
author = {Pith},
title = {Pith review of: Magnetic Field Dependence of the Spin Susceptibility on Conventional s-wave Superconductor LaRu$_4$P$_12$ Revealed by $^31$P-NMR and $^139$La-NMR},
year = {2026},
howpublished = {\url{https://pith.science/paper/PVFL5YJD}},
note = {Machine review of arXiv:2505.20702}
}
abstract
The magnetic field dependence of the spin part of Knight shift, which is proportional to the superconducting-state spin susceptibility, was investigated at two nuclear sites, $^{31}$P and $^{139}$La in a conventional s-wave superconductor LaRu$_4$P$_{12}$. After the analyses, we confirmed that the superconducting-state spin susceptibility is proportional to magnetic field, and connects to the normal-state spin susceptibility smoothly. This is a textbook example, when the superconductivity is broken with the orbital pair-breaking effect.
Figures
Figures from the paper (2 more)
Reference graph
Works this paper leans on
-
[1]
Bardeen, L
J. Bardeen, L. N. Cooper, and J. R. Schrieffer, Phys. Rev.108, 1175 (1957)
1957
- [2]
-
[3]
H. Fine, M. Lipsicas, and M. Strongin, Physics Letters A 29, 366–367 (1969)
work page 1969
-
[4]
S. Kitagawa, T. Higuchi, M. Manago, T. Yamanaka, K. Ishida, H. S. Jeevan, and C. Geibel, Phys. Rev. B 96, 134506 (2017)
work page 2017
- [5]
-
[6]
H. Tou, Y. Kitaoka, K. Asayama, N. Kimura, Y. ¯Onuki, E. Ya- mamoto, and K. Maezawa, Phys. Rev. Lett. 77, 1374 (1996)
work page 1996
-
[7]
H. Tou, Y. Kitaoka, K. Ishida, K. Asayama, N. Kimura, Y. ¯Onuki, E. Yamamoto, Y. Haga, and K. Maezawa, Phys. Rev. Lett. 80, 3129 (1998)
work page 1998
- [8]
Show all 44 references
-
[9]
A. J. Leggett, Rev. Mod. Phys. 47, 331 (1975)
1975
-
[10]
Matsumura, H
H. Matsumura, H. Fujibayashi, K. Kinjo, S. Kitagawa, K. Ishida, Y. Tokunaga, H. Sakai, S. Kambe, A. Nakamura, Y. Shimizu, Y. Homma, D. Li, F. Honda, and D. Aoki, J. Phys. Soc. Jpn. 92, 063701 (2023)
2023
-
[11]
J. Yang, J. Luo, C. Yi, Y. Shi, Y. Zhou, and G.-q. Zheng, Sci. Adv. 7, eabl4432 (2021)
2021
-
[12]
A. M. Clogston, Phys. Rev. Lett. 9, 266 (1962)
1962
-
[13]
B. S. Chandrasekhar, App. Phys. Lett. 1, 7 (1962)
1962
-
[14]
Maki and T
K. Maki and T. Tsuneto, Prog. Theor. Phys. 31, 945 (1964)
1964
-
[15]
N. R. Werthamer, E. Helfand, and P. C. Hohenberg, Phys. Rev. 147, 295 (1966)
1966
-
[16]
Tayama, A
T. Tayama, A. Harita, T. Sakakibara, Y. Haga, H. Shishido, R. Settai, and Y. Onuki, Phys. Rev. B 65, 180504 (2002)
2002
-
[17]
Maeno, S
Y. Maeno, S. Yonezawa, and A. Ramires, J. Phys. Soc. Jpn. 93, 062001 (2024)
2024
-
[18]
S. Uji, T. Konoike, Y. Hattori, T. Terashima, T. Oguchi, T. Yamada, D. Hirai, T. Ikenobe, and Z. Hiroi, J. Phys. Soc. Jpn. 92, 074703 (2023)
2023
-
[19]
Hagawa, M
M. Hagawa, M. Matsukawa, K. Niinuma, R. Kudo, Y. Mizushima, N. Kawarada, H. Yamamoto, K. Sano, Y. Ono, and T. Sasaki, J. Phys. Soc. Jpn. 93, 044705 (2024)
2024
-
[20]
Young, R
B.-L. Young, R. R. Urbano, N. J. Curro, J. D. Thompson, J. L. Sarrao, A. B. Vorontsov, and M. J. Graf, Phys. Rev. Lett. 98, 036402 (2007)
2007
-
[21]
Kitagawa, R
S. Kitagawa, R. Takaki, M. Manago, K. Ishida, and N. K. Sato, J. Phys. Soc. Jpn. 87, 013701 (2018)
2018
-
[22]
Kitagawa, G
S. Kitagawa, G. Nakamine, K. Ishida, H. S. Jeevan, C. Geibel, and F. Steglich, Phys. Rev. Lett. 121, 157004 (2018)
2018
-
[23]
Chronister, A
A. Chronister, A. Pustogow, N. Kikugawa, D. A. Sokolov, F. Jerzembeck, C. W. Hicks, A. P. Mackenzie, E. D. Bauer, and S. E. Brown, Proceedings of the National Academy of Sciences 118, e2025313118 (2021)
2021
-
[24]
Kinjo, M
K. Kinjo, M. Manago, S. Kitagawa, Z. Q. Mao, S. Yonezawa, Y. Maeno, and K. Ishida, Science 376, 397 (2022)
2022
-
[25]
Nakamine, K
G. Nakamine, K. Kinjo, S. Kitagawa, K. Ishida, Y. Tokunaga, H. Sakai, S. Kambe, A. Nakamura, Y. Shimizu, Y. Homma, D. Li, F. Honda, and D. Aoki, J. Phys. Soc. Jpn. 90, 064709 (2021)
2021
-
[26]
Nakamine, K
G. Nakamine, K. Kinjo, S. Kitagawa, K. Ishida, Y. Tokunaga, H. Sakai, S. Kambe, A. Nakamura, Y. Shimizu, Y. Homma, D. Li, F. Honda, and D. Aoki, Phys. Rev. B 103, L100503 (2021)
2021
-
[27]
Kinjo, H
K. Kinjo, H. Fujibayashi, S. Kitagawa, K. Ishida, Y. Toku- naga, H. Sakai, S. Kambe, A. Nakamura, Y. Shimizu, Y. Homma, D. X. Li, F. Honda, D. Aoki, K. Hiraki, M. Ki- mata, and T. Sasaki, Phys. Rev. B 107, L060502 (2023)
2023
-
[28]
Kitagawa, K
S. Kitagawa, K. Nakanishi, H. Matsumura, Y. Takahashi, K. Ishida, Y. Tokunaga, H. Sakai, S. Kambe, A. Nakamura, Y. Shimizu, Y. Homma, D. Li, F. Honda, A. Miyake, and D. Aoki, J. Phys. Soc. Jpn. 93, 123701 (2024)
2024
-
[29]
Nohara, M
M. Nohara, M. Isshiki, F. Sakai, and H. Takagi, J. Phys. Soc. Jpn. 68, 1078 (1999)
1999
-
[30]
Hanaguri, A
T. Hanaguri, A. Koizumi, K. Takaki, M. Nohara, H. Tak- agi, and K. Kitazawa, Physica B: Cond. Mat. 329–333, 1355 (2003)
2003
-
[31]
Meisner, Physica B+C 108, 763 (1981)
G. Meisner, Physica B+C 108, 763 (1981)
1981
-
[32]
Shirotani, T
I. Shirotani, T. Adachi, K. Tachi, S. Todo, K. Nozawa, T. Yagi, and M. Kinoshita, J. Phys. Chem. Solids 57, 211 (1996)
1996
-
[33]
Tsuda, T
S. Tsuda, T. Yokoya, T. Kiss, T. Shimojima, S. Shin, T. To- gasi, S. Watanabe, C. Q. Zhang, C. T. Chen, H. Sugawara, H. Sato, and H. Harima, J. Phys. Soc. Jpn. 75, 064711 (2006)
2006
-
[34]
Nakai, Y
Y. Nakai, Y. Hayashi, K. Ishida, H. Sugawara, D. Kikuchi, and H. Sato, Physica B: Cond. Matt. 403, 1109 (2008)
2008
-
[35]
Nakai, Y
Y. Nakai, Y. Hayashi, K. Kitagawa, K. Ishida, H. Sugawara, D. Kikuchi, and H. Sato, J. Phys. Soc. Jpn. 77, 333 (2008)
2008
-
[36]
Kinjo, S
K. Kinjo, S. Kitagawa, Y. Nakai, K. Ishida, H. Sugawara, and H. Sato, J. Phys. Soc. Jpn. 88, 065002 (2019)
2019
-
[37]
Jeitschko and D
W. Jeitschko and D. Braun, Acta Cryst. 33, 3401 (1977)
1977
-
[38]
V. A. Stenger, C. H. Pennington, D. R. Buffinger, and R. P. Ziebarth, Phys. Rev. Lett. 74, 1649 (1995)
1995
-
[39]
Nakai, K
Y. Nakai, K. Ishida, D. Kikuchi, H. Sugawara, and H. Sato, J. Phys. Soc. Jpn. 74, 3370 (2005)
2005
-
[40]
Momma and F
K. Momma and F. Izumi, J. Appl. Crystallogr.44, 1272 (2011)
2011
-
[41]
Uchiumi, I
T. Uchiumi, I. Shirotani, C. Sekine, S. Todo, T. Yagi, Y. Nakazawa, and K. Kanoda, J. Phys. Chem. Solids 60, 689 (1999)
1999
-
[42]
Fite and A
W. Fite and A. G. Redfield, Phys. Rev. Lett. 17, 381 (1966)
1966
-
[43]
E. H. Brandt, Phys. Rev. B 68, 054506 (2003)
2003
-
[44]
J. Li, B. L. Kang, D. Zhao, B. Lei, Y. B. Zhou, D. W. Song, S. J. Li, L. X. Zheng, L. P. Nie, T. Wu, and X. H. Chen, Phys. Rev. B 105, 054514 (2022). 5
2022
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.