Niobium's intrinsic coherence length and penetration depth revisited using low-energy muon spin spectroscopy and secondary-ion mass spectrometry
Source abstract
Abstract
We report direct, simultaneous measurements of the London penetration depth ($λ_L$) and Bardeen-Cooper-Schrieffer (BCS) coherence length ($ξ_0$) in oxygen-doped niobium, with impurity concentrations spanning the "clean" to "dirty" limits. Two depth-resolved techniques - low-energy muon spin spectroscopy (LE-$μ$SR) and secondary-ion mass spectrometry (SIMS) - were used to quantify the element's Meissner screening profiles, analyzed within a framework that accounts for nonlocal electrodynamics. The analysis indicates intrinsic length scales of $λ_L = 29.1(10)$ nm and $ξ_0 = 39.9(25)$ nm, corresponding to a Ginzburg-Landau (GL) parameter of $κ= 0.70(5)$. The obtained $λ_L$ and $κ$ values, accurately quantified at the nanoscale, are smaller than values commonly used in applications and modeling, and indicate that clean niobium lies at the boundary between type-I and type-II superconductivity, supporting the contemporary view that its intrinsic state may be type-I.
Source: openalex_inverted_index · Quality: baseline
Research context
OpenAlex topics: Muon and positron interactions and applications · High-Energy Particle Collisions Research · Neutrino Physics Research
SciVX discovery tags are terminology-based. OpenAlex topics are source-provided classifications.
Affiliations
TRIUMF
V6T 2A3
University of Victoria
V8P 5C2
Virginia Polytechnic Institute and State University
Pacific Northwest National Laboratory
Thomas Jefferson National Accelerator Facility
Paul Scherrer Institute
Funding
Natural Sciences and Engineering Research Council of Canada · Nuclear Physics · U.S. Department of Energy · Swiss Muon Source, Paul Scherrer Institute, Villigen, Switzerland
References (93)
Show source-provided references
1. 10.1098/rspa.1935.004810.1098/rspa.1935.0048 ↗
2. 10.1098/rspa.1953.004010.1098/rspa.1953.0040 ↗
3. 10.1103/PhysRev.108.117510.1103/PhysRev.108.1175 ↗
4. 10.1007/BF0139485110.1007/BF01394851 ↗
5. 10.1080/0010751810823154110.1080/00107518108231541 ↗
6. 10.1155/2013/10437910.1155/2013/104379 ↗
7. 10.1016/j.ssnmr.2015.02.00410.1016/j.ssnmr.2015.02.004 ↗
8. 10.1515/zpch-2021-315410.1515/zpch-2021-3154 ↗
9. 10.1080/0010751041000167680310.1080/00107510410001676803 ↗
10. 10.1088/0953-8984/16/40/01010.1088/0953-8984/16/40/010 ↗
11. 10.1103/PhysRevLett.92.08700110.1103/PhysRevLett.92.087001 ↗
12. 10.1103/PhysRevB.72.02450610.1103/PhysRevB.72.024506 ↗
13. 10.1016/j.physb.2005.11.06510.1016/j.physb.2005.11.065 ↗
14. 10.1103/PhysRevB.85.22050110.1103/PhysRevB.85.220501 ↗
15. 10.1103/PhysRevB.87.10450810.1103/PhysRevB.87.104508 ↗
16. 10.1103/PhysRevLett.84.495810.1103/PhysRevLett.84.4958 ↗
17. 10.1103/PhysRevLett.92.05760210.1103/PhysRevLett.92.057602 ↗
18. 10.1103/PhysRevB.81.18050210.1103/PhysRevB.81.180502 ↗
19. 10.1103/PhysRevB.88.18050110.1103/PhysRevB.88.180501 ↗
20. 10.1103/PhysRevB.85.06050610.1103/PhysRevB.85.060506 ↗
21. 10.1088/1361-6668/ab119e10.1088/1361-6668/ab119e ↗
22. 10.1103/PhysRevB.89.18050810.1103/PhysRevB.89.180508 ↗
23. 10.1103/PhysRevB.79.14451810.1103/PhysRevB.79.144518 ↗
24. For example, Nb has the highest critical temperature Tc≈9.25K [76] of the elements and the largest lower critical field Bc1≈170mT [76] of all type-II superconductors.
25. H. Padamsee · 202310.1002/9783527836314 ↗
26. 10.1063/5.017553210.1063/5.0175532 ↗
27. 10.1063/1.486601310.1063/1.4866013 ↗
28. 10.1063/5.015700010.1063/5.0157000 ↗
29. 10.1088/1361-6668/aa892610.1088/1361-6668/aa8926 ↗
30. 10.1103/PhysRevApplied.19.04401810.1103/PhysRevApplied.19.044018 ↗
31. 10.1103/PhysRevX.5.04102110.1103/PhysRevX.5.041021 ↗
32. 10.1103/PhysRevLett.120.24700110.1103/PhysRevLett.120.247001 ↗
33. 10.1103/PhysRevLett.125.02680210.1103/PhysRevLett.125.026802 ↗
34. 10.1088/1361-6668/ad146210.1088/1361-6668/ad1462 ↗
35. 10.1103/PhysRev.139.A151510.1103/PhysRev.139.A1515 ↗
36. 10.1016/j.nimb.2025.16595410.1016/j.nimb.2025.165954 ↗
37. 10.1063/5.022703910.1063/5.0227039 ↗
38. 10.1116/6.000084810.1116/6.0000848 ↗
39. 10.1116/6.000174110.1116/6.0001741 ↗
40. 10.1103/PhysRevApplied.13.01402410.1103/PhysRevApplied.13.014024 ↗
41. 10.1093/ptep/ptab05610.1093/ptep/ptab056 ↗
42. 10.1063/5.005946410.1063/5.0059464 ↗
43. 10.1063/5.019123410.1063/5.0191234 ↗
44. In a material's Meissner response, nonlocal effects manifest predominantly as (1) a reduced initial slope in the spatial decay of B(z), and (2) as a negative curvature in B(z), culminating in sign reversal when z≫λL. Qualitatively, these can be understood as a consequence of the finite spatial extent of Cooper pairs, causing each “partner” to experience different local fields, leading to a diminished screening capacity that is eventually “overcompensated” deeper below the surface [12].
45. 10.3389/fphy.2023.126987210.3389/fphy.2023.1269872 ↗
46. 10.1016/S0921-4534(99)00207-510.1016/S0921-4534(99)00207-5 ↗
47. See, for example, the literature average reported in Refs. [26, 30].
48. 10.1103/PhysRevAccelBeams.26.10310110.1103/PhysRevAccelBeams.26.103101 ↗
49. 10.1007/BF0061520110.1007/BF00615201 ↗
50. 10.1016/0040-6090(90)90079-S10.1016/0040-6090(90)90079-S ↗
51. 10.1515/ijmr-1976-67110510.1515/ijmr-1976-671105 ↗
52. 10.1051/jphys:01968002902-302400010.1051/jphys:01968002902-3024000 ↗
53. See Supplemental Material at http://link.aps.org/supplemental/10.1103/2nsw-n8gf for additional information on the sample preparation steps, the LE−µSR and SIMS measurement details, the analysis of the raw data, and tables summarizing all fitted values.
54. 10.1016/j.nima.2008.07.08110.1016/j.nima.2008.07.081 ↗
55. 10.1103/PhysRevLett.72.279310.1103/PhysRevLett.72.2793 ↗
56. 10.1016/S0169-4332(00)00857-610.1016/S0169-4332(00)00857-6 ↗
57. 10.1016/S0921-4526(00)00303-310.1016/S0921-4526(00)00303-3 ↗
58. 10.1016/j.phpro.2012.04.03910.1016/j.phpro.2012.04.039 ↗
59. 10.1016/S0168-583X(01)01166-110.1016/S0168-583X(01)01166-1 ↗
60. W. Eckstein · 199110.1007/978-3-642-73513-4 ↗
61. Z. Salman, TRIM.SP, https://bitbucket.org/zaher-salman/trimsp/ (2025).
62. A. Amato · 202410.1007/978-3-031-44959-8 ↗
63. 10.1103/RevModPhys.93.02501010.1103/RevModPhys.93.025010 ↗
64. 10.1093/ptep/ptac09710.1093/ptep/ptac097 ↗
65. 10.1016/j.phpro.2012.04.04210.1016/j.phpro.2012.04.042 ↗
66. As the ρ(z,E)s are obtained from Monte Carlo simulations performed at discrete E, here we adopt an empirical “interpolation” scheme (see Refs. [28, 30, 34] and the Supplemental Material [53]) to facilitate evaluation of Eq. (3) at arbitrary E.
67. The assumption of a mirrorlike reflection of electrons at a superconductor's surface can be contrasted with diffusive scattering. In reality, most materials likely exhibit behavior between these two limits. Pragmatically, either choice is permissible as both produce nearly identical screening lengths [68].
68. M. Tinkham · 1996
69. While the “dead layer” d is a sample-dependent quantity, nonzero values are common for Nb [12, 26, 27, 29, 30], likely due to its natural surface oxidation [49].
70. Even though Bapplied is parallel the sample surface in the measurements, the finite thickness of our samples (i.e., a few millimeters) is enough to ensure that Ñ>0 (cf. thin-film samples where Ñ→0 [12]).
71. Nonanalytic expressions equivalent to Eq. (6) from Pippard theory [2] can be derived from BCS theory [3, 4]. Importantly, quantitative differences between the two models are negligible [12].
72. Note that the “two-fluid” temperature dependence in Eq. (7) is known to deviate from BCS theory [3] when T≪Tc; however, it is known to accurately reproduce Nb's λ(T) (see, e.g., Ref. [12]).
73. Note that the factor of (π/2) included in Eq. (8) ensures equivalent length scales with BCS theory [4].
74. 10.1103/PhysRev.149.36810.1103/PhysRev.149.368 ↗
75. From fitting, we find that the effective demagnetization factor Ñ ranges from 0.03 to 0.06, with its exact value depending on the sample and analysis model. The staged analysis generally finds Ñ≳0.05, while Ñ≲0.04 is found in the direct approach. Exact values are given in Tables S4–S7 in the Supplemental Material [53].
76. 10.1103/PhysRev.149.23110.1103/PhysRev.149.231 ↗
77. 10.1103/PhysRev.128.59110.1103/PhysRev.128.591 ↗
78. 10.1103/PhysRev.113.120910.1103/PhysRev.113.1209 ↗
79. Explicitly, the weighed average of a set of n values vi with uncertainties δvi is given by (∑inwivi)/(∑inwi), where the weights wi≡1/δvi2. Similarly, the uncertainty in the weighted average is given by: 1/(∑inwi).
80. 10.1103/PhysRevB.106.L18050510.1103/PhysRevB.106.L180505 ↗
81. 10.3389/femat.2023.125940110.3389/femat.2023.1259401 ↗
82. 10.1103/PhysRevB.19.454510.1103/PhysRevB.19.4545 ↗
83. N. E. Alekseevskiy · Fiz. Metal. Metalloved. · 1974 · 37 · 63
84. 10.1103/PhysRevB.7.13610.1103/PhysRevB.7.136 ↗
85. 10.1007/BF0068162110.1007/BF00681621 ↗
86. 10.1103/PhysRevB.96.14450610.1103/PhysRevB.96.144506 ↗
87. 10.1103/PhysRevB.100.06450310.1103/PhysRevB.100.064503 ↗
88. 10.7567/JJAPS.26S3.91710.7567/JJAPS.26S3.917 ↗
89. 10.1016/0921-4534(89)90337-710.1016/0921-4534(89)90337-7 ↗
90. 10.1103/PhysRevB.104.06450410.1103/PhysRevB.104.064504 ↗
91. 10.1103/PhysRevB.111.09451910.1103/PhysRevB.111.094519 ↗
92. 10.1088/1361-6668/30/3/03400410.1088/1361-6668/30/3/034004 ↗
93. All LE−µSR data generated at PSI are available for download from http://musruser.psi.ch/ (using the SearchDB tool).