Insight into SRF cavity performance from simulations of Nb's surface oxide dissolution and diffusion

Rowan Becker; Ryan M. L. McFadden; Tobias Junginger

REVIEW 2 major objections 6 minor 74 references

Insight into SRF cavity performance from simulations of Nb's surface oxide dissolution and diffusion

T0 review · 2 major / 6 minor · reviewed 2026-08-14 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Oxygen released by vacuum annealing shifts niobium's Meissner screening current away from the surface.

desk verdict A careful, honest simulation paper whose headline current-redistribution result was partly in prior work and rests on a local-London assumption that is not well controlled near the surface. read the letter →

arxiv 2608.13540 v1 pith:6FXZDKOQ submitted 2026-08-13 cond-mat.supr-con cond-mat.mtrl-sciphysics.acc-ph

Ryan M. L. McFadden , Rowan Becker , Tobias Junginger This is my paper · ORCID

classification cond-mat.supr-concond-mat.mtrl-sciphysics.acc-ph

keywords niobiumSRFcavitiesoxygendopingMeissnerscreeningreaction-diffusionpenetrationdepthsupercurrentdensityvacuumheattreatment

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper builds a numerical bridge between a common industrial step—vacuum baking of niobium superconducting radio-frequency (SRF) cavities—and a microscopic superconducting quantity, the Meissner screening-current profile. It simulates the dissolution of Nb's native oxide and the inward diffusion of interstitial oxygen, then feeds the resulting depth-dependent oxygen profile into calculations of the electron mean-free path, the magnetic penetration depth, and the field profile through a generalized London equation. The central result is that oxygen doping changes where the screening current flows: the surface current is reduced and the current maximum moves several nanometres below the surface, into material whose depth-dependent critical current is larger. This matters because the near-surface screening current sets how high an RF field a cavity can sustain, so a quantitative relation between bake recipe and current profile offers a way to design heat treatments rather than adjust them empirically.

What carries the argument

The load-bearing machinery is a numerical chain with four links. First, Arrhenius rate constants for the stepwise oxide reduction Nb2O5 → NbO2 → NbO set the time-dependent oxygen source at the surface. Second, a Crank-Nicolson solution of the reaction-diffusion equation $\partial[O]/\partial t = D\,\partial^2[O]/\partial x^2 + q(x,t)$ produces the oxygen profile $[O](x)$. Third, that profile is converted into a depth-dependent London penetration depth $\lambda(T,x,B_0)$ through the residual resistivity, the electron mean-free path $\ell(x)$, the impurity-limit expression $\lambda_0(x)=\lambda_L\sqrt{1+\pi\xi_0/(2\ell(x))}$, and the nonlinear Meissner correction. Fourth, the generalized London equation $\lambda^2 B'' + 2\lambda\lambda' B' = B$ is solved numerically for the magnetic field, and the supercurrent follows from $J(x)=-(1/\mu_0)B'(x)$. The generalized London equation carries the central claim: with depth-dependent $\lambda$, the field profile is no longer a pure exponential, so the current maximum can move below the surface.

What would settle it

A depth-resolved measurement of the Meissner field profile of a niobium sample vacuum-annealed at 125 °C for 8 h would settle the claim: if the supercurrent density is largest at the surface and decays monotonically rather than peaking a few nanometres below it, the predicted redistribution is absent.

Watch

Extended reading notes

Core claim

The central claim is that oxygen doping redistributes the Meissner screening current in surface-treated niobium. For representative vacuum annealing, the supercurrent density $J(x)$ at the surface is lower than in clean Nb, its maximum is no longer at the surface but several nanometres inside the material, and both the surface value and the peak value of $J(x)$ are reduced relative to clean Nb. The normalized ratio $j(x)=J(x)/J_c(x)$, where $J_c(x)$ is the depth-dependent critical current density, stays below one for all depths in the simulated examples; the peak of $J(x)$ sits close to the depth where $J_c(x)$ is largest. The paper maps how the peak position and magnitude vary across temperatures from 50 °C to 200 °C and times from 0.5 h to 120 h, identifying a valley in peak-current reduction and a ridge in peak depth that locate favourable treatment conditions.

Load-bearing premise

The simulations assume that niobium's electromagnetic response stays local even where the oxygen profile makes the penetration depth vary over tens of nanometres, comparable to the coherence length; if nonlocal electrodynamics is required, the computed screening-current profile would change.

Editorial extensions

If this is right

If the claim is correct, a fixed bake dose does not merely change the magnitude of the screening current; it changes where that current flows, so surface currents and peak currents respond differently to temperature and time.
The parameter maps imply a low-temperature, moderate-time window in which the peak current is reduced and pushed toward the depth where the critical current is highest, and a hotter or longer window in which the profile homogenizes toward the dirty limit.
Because $j(x)=J(x)/J_c(x)<1$ everywhere in the simulated examples, the doped surface retains Meissner-state margin even at an applied field set equal to Nb's thermodynamic critical field.
The same numerical machinery is argued to extend to time-dependent and multi-step treatments and to temperatures up to about 400 °C, so the conclusions are not limited to the surveyed grid of recipes.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

One consequence the authors leave implicit: if the mechanism is right, heat-treatment design becomes an inverse problem—choose a temperature–time trajectory that yields a target $J(x)$ profile instead of scanning recipes empirically.
The same impurity-profile-to-current-redistribution chain should apply to other interstitial dopants; nitrogen-doped or nitrogen-infused niobium may show analogous subsurface current peaks, a case the paper notes has not yet been modeled.
The local-electrodynamics assumption can be tested by comparing predicted screening profiles with depth-resolved measurements; disagreement would show exactly where the generalized London equation needs nonlocal corrections.
The shift of the current peak may also change the field at which vortices first enter, since the controlling quantity would be the depth-resolved ratio $j(x)$ rather than the surface current alone.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

The headline result — oxygen doping moves the Meissner screening current peak a few nanometres below the surface — is not new: the paper says the valley was identified in Ref. [9] and the deformed J(x) is consistent with that paper. What is new is the systematic scan over temperature and time, the maps of J~, J~0, and x~, and the j(x)<1 criterion. That is a useful extension of an existing framework, not a discovery.

The paper is careful. The reaction-diffusion setup is standard, the Crank-Nicolson solver is correctly described, and the parameter choices are documented with literature sources. They explicitly reject an outlier dataset with reasons. The self-cited LE-μSR work is directly relevant and not a circular fit: the output current redistribution does not feed back into the input parameters.

The soft spots. The load-bearing step is Eq. (30), the generalized London equation with a depth-dependent λ(x). The problem is that the system is not in the dirty limit: at the surface ℓ(x) is 30–50 nm while ξ0=39 nm, and λ(x) varies over 10–40 nm. Nonlocal effects are known to suppress the surface current and move the peak inward even in homogeneous clean Nb (Suter et al., PRB 72, 024506). The footnote in Sec. IIC2 cites their own Ref. [11] to wave this away, but that reference measures λ and ξ0; it does not validate a local treatment for an inhomogeneous profile with a gradient length comparable to ξ0. So the predicted peak shift and the 'valley' maps may be partly nonlocal boundary artifacts, not cleanly attributable to oxygen doping. This needs to be addressed — at minimum with a Pippard/BCS nonlocal calculation for a few representative profiles, or a controlled perturbative estimate.

Two smaller issues: the uncertainties in Table II are never propagated, so the recipe maps come without error bars; and there is no code or data deposit, which matters for a paper whose selling point is a fast simulation tool. Experimental validation of the predicted J(x) is absent, though the authors acknowledge that.

Who is this for? The SRF cavity community, as a quick screening tool for comparing bake recipes. It deserves a serious referee — the questions are relevant and the paper is honestly written — but the electrodynamics section needs major work before the headline numbers can be trusted. I would send it to peer review, but with a strong request for a nonlocal check and for the code/data. Not a desk reject.

Referee Report

2 major / 6 minor

Summary. The paper develops a one-dimensional reaction-diffusion model for dissolution of Nb's native oxide and diffusion of interstitial oxygen during vacuum annealing, and post-processes the computed oxygen profiles into depth-dependent electron mean free path, penetration depth, critical current density, and Meissner screening profiles via a generalized local London equation. The headline claim is that oxygen doping reduces the supercurrent at the surface and moves its maximum several nanometres below the surface, with a systematic scan of about 72,000 temperature-duration combinations mapping the resulting 'valley' and 'cliff' structure in metrics such as peak-current ratio and peak position. The authors conclude that low-temperature or long-time treatments can engineer current distributions favorable to SRF cavity operation.

Significance. If the local-electrodynamics calculation is justified, this is a useful contribution: it converts previously scattered kinetic data into falsifiable predictions for depth-resolved screening profiles, it documents the numerical scheme in enough detail to be reproduced (Appendix C), and it connects the simulated parameter space to empirical recipes such as 110 C/60 h and 145 C/3 h. A notable strength is that the input parameters are mostly independent literature values, so the central mechanism is not fitted to the output; the output current redistribution was not used to set any of the reaction, diffusion, or resistivity constants. The main risk is that the qualitative headline result depends on a local London assumption that is not established for the simulated mean-free-path range, and the paper does not quantify how this affects the claimed current redistribution.

major comments (2)

[Table II, Figs. 3-5, and Section III] The manuscript reports uncertainties for the reaction rates and diffusivity in Table II and Appendix B (for example, E_A,1 = 133.7(15) kJ/mol and D_0 = 0.59(9) x 10^-6 m^2/s), but all simulated profiles and the parameter-space maps in Figs. 3-5 use point estimates only. The locations of the 'valley' in Fig. 3 and the 'cliff' in Fig. 5 are set by the kinetic competition between oxide dissolution and oxygen diffusion, so a Monte Carlo or linearized propagation of the tabulated uncertainties is needed to show that the claimed optimal treatment regions are robust. Without such an analysis, the practical recommendations in Section IV are not quantitatively supported.
[Eq. (30) and Fig. 2] The comparison against the 'clean' and 'dirty' exponential baselines in Fig. 2 is not a sufficient control for the central claim. The clean-limit baseline is computed from the same local London equation, which in a homogeneous clean superconductor is known to overestimate the surface current because nonlocal effects are omitted; see Ref. [36]. The treated profile could therefore show a reduced surface current and an interior peak partly because the reference itself is the wrong local limit. This is a consequence of the same issue raised in the first major comment, and it should be addressed before the abstract's 'we find' claim is stated as established.

minor comments (6)

[Sections II-IIB, Appendix C] Typos and duplicated words should be corrected: 'impliment' in Section II, 'in agreement with with Ref. [9]' in Section IIB, 'where where' and 'matricies' in Appendix C, and 'aaclear' later in Appendix C.
[Table II, k1 row] The entry for k1 from Ref. [6] is garbled as '8.0x10 7.0(4) approximately 0.12(14)x10 9'; the exponent notation should be cleaned so the value and its uncertainty can be compared directly with the other rows.
[Table I] The row '[O]bulk x Delta x = 0.005 at.% nm' conflicts with the text's statement that [O]bulk = 0.005 at.%; clarify whether the bulk concentration is a volumetric concentration or an areal source term.
[Section III, after Fig. 3] The fitted expressions for the 'valley' and 'ridge' mix units and take logarithms of dimensioned quantities; define all symbols and state the unit conventions for T and t so that the fits are reproducible.
[Appendix D] The captions of Figs. 6-9 repeat the full methodology for each panel; consider abbreviating them and checking that the 'clean' and 'dirty' dotted lines remain legible in every panel.
[References] Ref. [50] is cited as an arXiv preprint; if a peer-reviewed version now exists, the citation should be updated.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the simulated current redistribution follows from independently parameterized reaction-diffusion and London-model inputs; no output quantity was used to set an input.

full rationale

The derivation chain is self-contained: oxygen profiles come from a reaction-diffusion equation (Eq. 17) with Arrhenius rate constants and diffusivities taken from the literature (XPS, SIMS, and tracer-diffusion studies), and the superconducting response is then computed from those profiles using standard local-limit relations (Eqs. 21-32). The central result, that J(x) is suppressed at the surface and peaked below it, is a numerical consequence of solving the inhomogeneous London equation (Eq. 30) with a depth-dependent lambda(x); it is not used to set any kinetic, diffusive, or electromagnetic parameter. The self-cited LE-muSR values for lambda_L and xi_0 (Ref. [11]) are independent experimental inputs, not outputs of this model's fit, and the paper does not claim to extract any parameter from the predicted current profile. The local-electrodynamics assumption behind Eq. (30) is a physical modeling assumption and a potential correctness risk, but it is not a circularity because the target result is not encoded in that assumption by construction. Figures 3-5 are descriptive metrics of the simulation output and are explicitly connected to the earlier independent calculation in Ref. [9]; presenting them as predictions of the model, rather than as fits to external data, does not make the derivation circular.

Assumptions & free parameters 16 free parameters · 9 assumptions · 0 invented entities

The simulation outcome is controlled by 16 literature-derived material and kinetic parameters, the plane-source dissolution model, the local generalized London equation, and implicit boundary conditions. The central qualitative result is robust to moderate parameter changes, but no uncertainty propagation or experimental validation is provided to test robustness quantitatively.

free parameters (16)

A0,1 = 9.75e8 s^-1
Arrhenius prefactor for k1, adopted as weighted average of Refs. [6,8,9].
EA,1 = 133.7 kJ/mol
Activation energy for k1 from the same weighted average; controls oxide dissolution rate at bake temperatures.
A0,2 = 1.5e11 s^-1
Arrhenius prefactor for k2, taken from Ref. [6] multi-point XPS method.
EA,2 = 177 kJ/mol
Activation energy for k2 from Ref. [6]; the slower NbO2 dissolution channel shapes the second valley feature.
D0 = 0.59e-6 m^2/s
Oxygen diffusivity prefactor via octahedral interstitial sites, from Refs. [21,22].
EA,D = 109.7 kJ/mol
Oxygen migration barrier in Nb from Refs. [21,22].
sigma_e = 3.7e-16 Ohm m^2
Empirical material constant linking electron mean free path to residual resistivity in Eq. (21), from Ref. [29].
a_O = 4.5e-12 Ohm m ppma^-1
Linear residual resistivity coefficient for oxygen in Nb in Eq. (22), from Ref. [30]; bridges chemistry to electrodynamics.
lambda_L = 29 nm
London penetration depth for pure Nb; literature average plus recent LE-muSR measurement in Ref. [11].
xi_0 = 39 nm
BCS coherence length for pure Nb; literature average plus recent LE-muSR measurement in Ref. [11].
[Nb2O5]0 * Delta x = 200 at.% nm
Initial surface oxygen reservoir in the native oxide, following Refs. [8,9].
[O]0 * Delta x = 5 at.% nm
Initial surface interstitial oxygen content, following Refs. [7,9,51].
[O]bulk = 0.005 at.%
Assumed homogeneous bulk oxygen concentration before annealing.
Bc = 199.3 mT
Thermodynamic critical field of Nb at 0 K from Ref. [37]; used in Jc and the NLME correction.
B0 = 199.3 mT
Applied field set equal to Bc in Eq. (25) and the screening calculation, chosen to represent high-field RF operation.
Tc = 9.25 K
Critical temperature of Nb, used in Eq. (23).

assumptions (9)

domain assumption First-order oxide dissolution kinetics with k3 = 0
Section IIA, Eqs. (4)-(7); NbO is assumed not to dissolve up to at least 1000 K based on XPS in Ref. [6].
domain assumption Oxygen source is a plane source at the surface (delta function)
Eq. (18), justified by native oxide thickness below about 5 nm (Ref. [10]); used to inject dissolved oxygen into the diffusion domain.
domain assumption Oxygen diffusivity in the surface oxide equals diffusivity in bulk Nb
Section IID assumptions; authors note this is questionable for thick anodic oxides but reasonable for thin native oxides, following Refs. [7,9].
domain assumption Interstitial oxygen migrates through octahedral sites in bcc Nb
Appendix B; octahedral D0 and EA,D are adopted, with the tetrahedral path giving similar values.
domain assumption Local electrodynamics are adequate for Nb screening
Section IIC2 footnote; nonlocal effects are assumed small even in clean Nb, which underlies Eq. (30).
domain assumption Generalized London equation (30) governs B(x) for depth-dependent lambda
Section IIC4, Eq. (30), cited to Refs. [46-48]; the entire supercurrent profile is derived from this equation.
domain assumption Neumann boundary conditions confine oxygen inside the simulation domain
Appendix C Eqs. (C8)-(C9); no oxygen flux leaves through x=0 or x=L, so all dissolved oxygen stays in Nb.
domain assumption Residual resistivity is a linear sum of impurity contributions
Section IIC1, Eq. (22); uses dilute-limit proportionality a_O for oxygen, applied at near-surface concentrations around 0.2-0.4 at.% in the simulations.
standard math BCS impure-superconductor formula for penetration depth
Eq. (24), from Refs. [34,39]; the factor pi/2 is taken from BCS theory for impure superconductors.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Insight into SRF cavity performance from simulations of Nb's surface oxide dissolution and diffusion." pith.science (2026). https://pith.science/paper/6FXZDKOQ

@misc{pith2026260813540,
  author       = {Pith},
  title        = {Pith review of: Insight into SRF cavity performance from simulations of Nb's surface oxide dissolution and diffusion},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/6FXZDKOQ}},
  note         = {Machine review of arXiv:2608.13540}
}

abstract

We report simulations of the dissolution and diffusion of Nb's surface oxide layer in vacuum. While this chemical doping process is important for the surface preparation of Nb superconducting radio frequency (SRF) cavities - common components of particle accelerators - quantitatively linking the resulting oxygen distributions to superconducting performance remains challenging. In this work, we simulate the reaction-diffusion process numerically for treatment temperatures $T = 50^{\circ}$C to $200^{\circ}$C and times $t = 0.5$ h to $120$ h, and calculate the effect of the spatially inhomogeneous oxygen doping on Nb's superconducting properties. We find that oxygen doping redistributes the Meissner screening current, reducing its value at the surface and shifting its maximum several nanometres into the material. These results provide a microscopic link between oxygen diffusion profiles and the electromagnetic response of Nb relevant for SRF cavity operation. This work provides a quantitative framework linking oxygen diffusion profiles to superconducting performance and establishes a foundation for future studies involving time-dependent and multi-step heat treatment protocols.

Figures

Figures reproduced from arXiv: 2608.13540 by the authors.

**Figure 2.** FIG. 2. Supercurrent density [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗

**Figure 3.** FIG. 3. Ratio of the maximum supercurrent density in treated Nb [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗

Figures from the paper (5 more)

**Figure 5.** Figure 5: FIG. 5. Position of the maximum supercurrent density in treated [PITH_FULL_IMAGE:figures/full_fig_p009_5.png]

**Figure 6.** Figure 6: FIG. 6. Simulated oxygen defect profile [PITH_FULL_IMAGE:figures/full_fig_p014_6.png]

**Figure 7.** Figure 7: FIG. 7. Simulated oxygen defect profile [PITH_FULL_IMAGE:figures/full_fig_p014_7.png]

**Figure 8.** Figure 8: FIG. 8. Simulated oxygen defect profile [PITH_FULL_IMAGE:figures/full_fig_p015_8.png]

**Figure 9.** Figure 9: FIG. 9. Simulated oxygen defect profile [PITH_FULL_IMAGE:figures/full_fig_p015_9.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

74 extracted references · 72 canonical work pages

[9]

Kubo, How high a field has been and can be achieved in superconducting bulk niobium cavities: the role of RRR, Jpn

T. Kubo, How high a field has been and can be achieved in superconducting bulk niobium cavities: the role of RRR, Jpn. J. Appl. Phys.64, 018002 (2025)

work page 2025
[11]

Posen, A

S. Posen, A. Romanenko, A. Grassellino, O. Melnychuk, and D. Sergatskov, Ultralow surface resistance via vacuum heat treatment of superconducting radio-frequency cavities, Phys. Rev. Appl.13, 014024 (2020)

work page 2020
[36]

B. B. Goodman and G. Kuhn, Influence des défauts étendus sur les propriétés supraconductrices du niobium, J. Phys. France29, 240 (1968)

work page 1968
[1]

scattering

Electron Mean-Free-Path The electron mean-free-pathℓ defines the average distance thechargecarriertravelsinamaterialbeforeitsenergy/direction is(substantially)alteredbya“scattering”event(e.g.,“colliding” with an impurity atom or phonon). In this sense,ℓ is directly relatedtotheabundanceoflatticedefects(e.g.,impurityatoms, vacancies, dislocations, grain bo...

work page
[2]

effective

Magnetic Penetration Depth In brief, for a superconductor in its Meissner state, the magneticpenetrationdepth 𝜆describesthecharacteristiclength scale(typically10sto100sofnanometers)overwhichexternal magnetic flux may “creep” below the surface before being completely screened deeper in its interior. In the limit that the superconductor’selectrodynamicsarep...

work page
[3]

correction

Non-linear Meissner Effect Empirically, it is well-known that𝜆(𝑇,𝑥) is affected by an applied magnetic field𝐵0, a result that was first understood theoreticallyusingGinzburg-Landau(GL)theory[ 40,41],with later developments extending it to unconventional supercon- ductors (see, e.g., [42]) Qualitatively, the effect accounts for 𝜆(𝑇,𝑥) ’s quadratic increase...

work page
[4]

Meissner Screening Profile Keeping with our assumption of local electrodynamics (see Section IIC2), the “classic” London result [31] describing the spatial dependence of magnetic-flux-penetration below a Meissner state superconductor’s surface can be written in a generalized form so as to account for spatial inhomogeneities (e.g., from a depth-dependent i...

work page
[5]

Supercurrent Density In 1D, the supercurrent density𝐽(𝑥) is related to𝐵(𝑥) by its first spatial derivative: 𝐽(𝑥)≡− 1 𝜇0 d d𝑥𝐵(𝑥),(31) where𝜇0 =1.256 637 061 27(20)N A −2 is the vacuum mag- netic permeability [19]. Ultimately, this is the most interesting quantity for SRF applications, where there is a strong desire to: minimize its value right at the surf...

work page

Show all 74 references

[6]

clean” and “dirty

Critical Supercurrent Density The critical supercurrent density𝐽𝑐 defines the maximum supercurrent that the material can withstand before magnetic vortices penetrate its surface (see, e.g., [49]). Within London theory, which assumes a local response, it can be written as: 𝐽𝑐(𝑇...
[7]

Goldilocks

that a well-defined supercurrent peak exceeding the value at the surface is most favourable at low temperatures and modest baking times; and 2) there is a “Goldilocks” amount of baking thatmaximizesthisquantity(i.e.,toolittleandthepeakremains close to the surface, too much and...
[8]

Padamsee,Superconducting Radiofrequency Technology for Accelerators: State of the Art and Emerging Trends(Wiley, Weinheim, 2023)

H. Padamsee,Superconducting Radiofrequency Technology for Accelerators: State of the Art and Emerging Trends(Wiley, Weinheim, 2023)

2023
[10]

Ciovati, Effect of low-temperature baking on the radio- frequency properties of niobium superconducting cavities for particle accelerators, J

G. Ciovati, Effect of low-temperature baking on the radio- frequency properties of niobium superconducting cavities for particle accelerators, J. Appl. Phys.96, 1591 (2004)

2004
[12]

H. Ito, H. Araki, K. Takahashi, and K. Umemori, Influence of furnacebakingonQ-Ebehaviorofsuperconductingaccelerating cavities, Prog. Theor. Exp. Phys.2021, 071G01 (2021)

2021
[13]

B.R.King,H.C.Patel,D.A.Gulino,andB.J.Tatarchuk,Kinetic measurements of oxygen dissolution into niobium substrates: In situ X-ray photoelectron spectroscopy studies, Thin Solid Films 192, 351 (1990)

1990
[14]

Ciovati, Improved oxygen diffusion model to explain the effect of low-temperature baking on high field losses in niobium superconducting cavities, Appl

G. Ciovati, Improved oxygen diffusion model to explain the effect of low-temperature baking on high field losses in niobium superconducting cavities, Appl. Phys. Lett.89, 022507 (2006)

2006
[15]

E. M. Lechner, J. W. Angle, F. A. Stevie, M. J. Kelley, C. E. Reece, and A. D. Palczewski, RF surface resistance tuning of superconducting niobium via thermal diffusion of native oxide, Appl. Phys. Lett.119, 082601 (2021)

2021
[16]

E. M. Lechner, J. W. Angle, A. D. Palczewski, F. A. Stevie, M. J. Kelley, and C. E. Reece, Oxide dissolution and oxygen diffusion scenarios in niobium and implications on the Bean-Livingston barrier in superconducting cavities, J. Appl. Phys.135, 133902 (2024)

2024
[17]

Halbritter, On the oxidation and on the superconductivity of niobium, Appl

J. Halbritter, On the oxidation and on the superconductivity of niobium, Appl. Phys. A43, 1 (1987)

1987
[18]

R. M. L. McFadden, J. W. Angle, E. M. Lechner, M. J. Kelley, C. E. Reece, M. A. Coble, T. Prokscha, Z. Salman, A. Suter, and T. Junginger, Niobium’s intrinsic coherence length and penetration depth revisited using low-energy muon spin spec- troscopy and secondary-ion mass spec...

2026
[19]

R. M. L. McFadden, M. Asaduzzaman, T. Prokscha, Z. Salman, A. Suter, and T. Junginger, Depth-resolved measurements of the Meissner screening profile in surface-treated Nb, Phys. Rev. Appl.19, 044018 (2023)

2023
[20]

Strong Meissner screening change in superconducting radio frequency cavities due to mild baking

R. M. L. McFadden, M. Asaduzzaman, and T. Junginger, Com- ment on “Strong Meissner screening change in superconducting radio frequency cavities due to mild baking” [Appl. Phys. Lett. 104, 072601 (2014)], Appl. Phys. Lett.124, 086101 (2024)

2014
[21]

R.M.L.McFaddenandT.Junginger,Searchforinhomogeneous Meissner screening in Nb induced by low-temperature surface treatments, AIP Adv.14, 095320 (2024)

2024
[22]

Wiencek, and J

C.Bate, K.Kasprzak, D.Reschke, L.Steder, L.Trelle, H.Weise, M. Wiencek, and J. Wolff, Correlation of srf performance to oxygen diffusion length of medium temperature heat treated cavities, Supercond. Sci. Technol.38, 025003 (2025)

2025
[23]

Prudnikava, Y

A. Prudnikava, Y. Tamashevich, A. Makarova, D. Smirnov, and J. Knobloch, In-situ synchrotron x-ray photoelectron spec- troscopy study of medium-temperature baking of niobium for srf application, Supercond. Sci. Technol.37, 075007 (2024)

2024
[24]

Tamashevich, A

Y. Tamashevich, A. Prudnikava, A. Matveenko, A. Neumann, O.Kugeler,andJ.Knobloch,Improvedrfperformanceofniobium cavities via in-situ vacuum heat treatment technique, Supercond. Sci. Technol.38, 045006 (2025)

2025
[25]

Checchin and A

M. Checchin and A. Grassellino, High-field Q-slope mitiga- tion due to impurity profile in superconducting radio-frequency cavities, Appl. Phys. Lett.117, 032601 (2020)

2020
[26]

E.Tiesinga,P.J.Mohr,D.B.Newell,andB.N.Taylor,CODATA recommended values of the fundamental physical constants: 2018, Rev. Mod. Phys.93, 025010 (2021)

2021
[27]

A. D. LeClaire, Diffusion of C, N, and O in metals, inLandolt- Börnstein — Group III Condensed Matter, Vol. 26, edited by H. Mehrer (Springer-Verlag, Berlin, 1990) Chap. 8, pp. 471–503

1990
[28]

F. J. M. Boratto and R. E. Reed-Hill, Computer optimization of elastic after-effect curves: Niobium relaxation from room temperature to114◦C, Scr. Metall.11, 709 (1977)

1977
[29]

F. J. M. Boratto and R. E. Reed-Hill, Statistical calculations on Arrhenius lines for the diffusion of oxygen and nitrogen in group V refractory metals, Mater. Sci. Eng.43, 97 (1980)

1980
[30]

D. N. Beshers, An investigation of interstitial sites in the bcc lattice, J. Appl. Phys.36, 290 (1965)

1965
[31]

Crank,The Mathematics of Diffusion, 2nd ed

J. Crank,The Mathematics of Diffusion, 2nd ed. (Oxford Univer- sity Press, London, 1975)

1975
[32]

Hundsdorfer and J

W. Hundsdorfer and J. Verwer,Numerical Solution of Time- Dependent Advection-Diffusion-Reaction Equations, Springer Series in Computational Mathematics, Vol. 33 (Springer, Berlin, 2003)

2003
[33]

Crank and P

J. Crank and P. Nicolson, A practical method for numerical evaluation of solutions of partial differential equations of the heat-conduction type, Math. Proc. Cambridge Philos. Soc.43, 50 (1947)

1947
[34]

K. K. Schulze, Preparation and characterization of ultra-high- purity niobium, J. Met.33, 33 (1981)

1981
[35]

Koethe and J

A. Koethe and J. I. Moench, Preparation of ultra high purity niobium, Mater. Trans., JIM41, 7 (2000)

2000
[37]

Schulze, J

K. Schulze, J. Fuß, H. Schultz, and S. Hofmann, Einfluß intersti- tieller frremdatome auf den restwiderstand von reinem niob, Z. Metallkd.67, 737 (1976)

1976
[38]

London and H

F. London and H. London, The electromagnetic equations of the supraconductor, Proc. R. Soc. London A149, 71 (1935)

1935
[39]

A. B. Pippard, An experimental and theoretical study of the relation between magnetic field and current in a superconductor, Proc. R. Soc. London A216, 547 (1953)

1953
[40]

Bardeen, L

J. Bardeen, L. N. Cooper, and J. R. Schrieffer, Theory of super- conductivity, Phys. Rev.108, 1175 (1957)

1957
[41]

Halbritter, On the penetration of the magnetic field into a superconductor, Z

J. Halbritter, On the penetration of the magnetic field into a superconductor, Z. Phys.243, 201 (1971)

1971
[42]

Tinkham,Introduction to Superconductivity, 2nd ed., Inter- national Series in Pure and Applied Physics (McGraw-Hill, New York, 1996)

M. Tinkham,Introduction to Superconductivity, 2nd ed., Inter- national Series in Pure and Applied Physics (McGraw-Hill, New York, 1996)

1996
[43]

Suter, E

A. Suter, E. Morenzoni, N. Garifianov, R. Khasanov, E. Kirk, H. Luetkens, T. Prokscha, and M. Horisberger, Observation of nonexponential magnetic penetration profiles in the Meissner state: A manifestation of nonlocal effects in superconductors, Phys. Rev. B72, 024506 (2005)

2005
[44]

D. K. Finnemore, T. F. Stromberg, and C. A. Swenson, Super- conducting properties of high-purity niobium, Phys. Rev.149, 231 (1966)

1966
[45]

D. A. Turner, G. Burt, and T. Junginger, No interface energy barrier and increased surface pinning in low temperature baked 17 niobium, Sci. Rep.12, 5522 (2022)

2022
[46]

P. B. Miller, Penetration depth in impure superconductors, Phys. Rev.113, 1209 (1959)

1959
[47]

V. L. Ginzburg and L. D. Landau, K teorii sverkhprovodimosti, Zh. Eksp. Teor. Fiz.20, 1064 (1950)

1950
[48]

V. L. Ginzburg and L. D. Landau, On the theory of superconduc- tivity,inCollectedPapersofL.D.Landau,editedbyD.TerHaar (Pergamon Press, Oxford, 1965) Chap. 73, pp. 546–568

1965
[49]

D. Xu, S. K. Yip, and J. A. Sauls, Nonlinear Meissner effect in unconventional superconductors, Phys. Rev. B51, 16233 (1995)

1995
[50]

Makita, C

J. Makita, C. Sundahl, G. Ciovati, C. B. Eom, and A. Gurevich, Nonlinear Meissner effect in Nb3Sn coplanar resonators, Phys. Rev. Res.4, 013156 (2022)

2022
[51]

L. P. Gor’kov, Microscopic derivation of the Ginzburg-Landau equations in the theory of superconductivity, Sov. Phys. JETP9, 1364 (1959)

1959
[52]

L. P. Gor’kov, Theory of superconducting alloys in a strong magnetic field near the critical temperature, Sov. Phys. JETP10, 998 (1960)

1960
[53]

R. W. Simon and P. M. Chaikin, Josephson tunneling studies of magnetic screening in proximity-superconducting silver, Phys. Rev. B23, 4463 (1981)

1981
[54]

J. R. Cave and J. E. Evetts, Critical temperature profile determi- nation using a modified London equation for inhomogeneous superconductors, J. Low Temp. Phys.63, 35 (1986)

1986
[55]

M. S. Pambianchi, J. Mao, and S. M. Anlage, Magnetic screen- ing in proximity-coupled superconductor/normal-metal bilayers, Phys. Rev. B50, 13659 (1994)

1994
[56]

Kubo, Multilayer coating for higher accelerating fields in superconductingradio-frequencycavities: areviewoftheoretical aspects, Supercond

T. Kubo, Multilayer coating for higher accelerating fields in superconductingradio-frequencycavities: areviewoftheoretical aspects, Supercond. Sci. Technol.30, 023001 (2017)

2017
[57]

Grassellino, A

A. Grassellino, A. Romanenko, D. Bice, O. Melnychuk, A. C. Crawford, S. Chandrasekaran, Z. Sung, D. A. Sergatskov, M. Checchin, S. Posen, M. Martinello, and G. Wu, Acceler- ating fields up to49 MV m−1 in TESLA-shape superconducting RF niobium cavities via75◦Cvacuum bake, arXiv...

arXiv
[58]

Arfaoui, C

I. Arfaoui, C. Guillot, J. Cousty, and C. Antoine, Evidence for a large enrichment of interstitial oxygen atoms in the nanometer- thick metal layer at the NbO/Nb (110) interface, J. App. Phys. 91, 9319 (2002)

2002
[59]

Oechsner, J

H. Oechsner, J. Giber, H. J. Füßer, and A. Darlinski, Phase transition and oxide dissolution processes in vacuum-annealed anodic Nb2O5/Nb systems, Thin Solid Films124, 199 (1985)

1985
[60]

R. M. L. McFadden, M. Asaduzzaman, T. J. Buck, D. L. Cortie, M. H. Dehn, S. R. Dunsiger, R. F. Kiefl, R. E. Laxdal, C. D. P. Levy,W.A.MacFarlane,G.D.Morris,M.R.Pearson,E.Thoeng, and T. Junginger, Depth-resolved measurement of the Meissner screening profile in a niobium thin fi...

2023
[61]

Visentin, Y

B. Visentin, Y. Gasser, and J. P. Charrier, First results on fast baking, Physica C441, 66 (2006)

2006
[62]

Visentin, Y

B. Visentin, Y. Gasser, M. Brunchon, F. Eozenou, and J. P. Charrier, Optimization of baking parameters for electropolished niobium superconducting cavities, inProceedings of SRF2007, InternationalWorkshoponRFSuperconductivityNo.13,Peking University (JACoW Publishing, Geneva, S...

2007
[63]

Dhakal, Nitrogen doping and infusion in SRF cavities: A review, Phys

P. Dhakal, Nitrogen doping and infusion in SRF cavities: A review, Phys. Open5, 100034 (2020)

2020
[64]

Y. Ye, Z. Yang, S. Wu, J. Li, T. Liu, G. Jiang, H. Li, Z. Jiang, J.Liu, Y.He,andW.Zhan,Calculationandunderstandingofthe nitrogen doping process in niobium for superconducting radio frequency accelerator, Physica C628, 1354616 (2025)

2025
[65]

R. W. Powers and M. V. Doyle, Diffusion of interstitial solutes in the group V transition metals, J. Appl. Phys.30, 514 (1959)

1959
[66]

Gibala and C

R. Gibala and C. A. Wert, The clustering of oxygen in solid solution in niobium — I. experimental, Acta Metall.14, 1095 (1966)

1966
[67]

Kirchheim, J

R. Kirchheim, J. Mathuni, and E. Fromm, Effektive ladung und diffusionskonstante von sauerstoff in niob, Z. Metallkde.68, 97 (1977)

1977
[68]

R. A. Perkins and R. A. Padgett Jr., Oxygen diffusion in niobium and Nb Zr alloys, Acta Metall.25, 1221 (1977)

1977
[69]

F. J. M. Boratto and R. E. Reed-Hill, On the calculations of the diffusion coefficients of oxygen and nitrogen in niobium, Metall. Trans. A8, 1233 (1977)

1977
[70]

R.Kirchheim,Metalsassinksandbarriersforinterstitialdiffusion with examples for oxygen diffusion in copper, niobium and tantalum, Acta Metall.27, 869 (1979)

1979
[71]

R. J. Lauf and C. J. Altstetter, Diffusion and trapping of oxygen in refractory metal alloys, Acta Metall.27, 1157 (1979)

1979
[72]

R.FarraroandR.B.McLellan,Thediffusionofheavyinterstitial solute atoms in body-centered cubic metals, Mater. Sci. Eng.39, 47 (1979)

1979
[73]

Michel, O

B. Michel, O. Kanert, and B. Günther, Motion of oxygen in vanadium and niobium studied by nuclear magnetic resonance, Acta Metall. Mater.42, 3409 (1994)

1994
[74]

J. G. Charney, R. Fjörtoft, and J. von Neumann, Numerical integration of the barotropic vorticity equation, Tellus2, 237 (1950)

1950

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.

[1] [9]

Kubo, How high a field has been and can be achieved in superconducting bulk niobium cavities: the role of RRR, Jpn

T. Kubo, How high a field has been and can be achieved in superconducting bulk niobium cavities: the role of RRR, Jpn. J. Appl. Phys.64, 018002 (2025)

work page 2025

[2] [11]

Posen, A

S. Posen, A. Romanenko, A. Grassellino, O. Melnychuk, and D. Sergatskov, Ultralow surface resistance via vacuum heat treatment of superconducting radio-frequency cavities, Phys. Rev. Appl.13, 014024 (2020)

work page 2020

[3] [36]

B. B. Goodman and G. Kuhn, Influence des défauts étendus sur les propriétés supraconductrices du niobium, J. Phys. France29, 240 (1968)

work page 1968

[4] [1]

scattering

Electron Mean-Free-Path The electron mean-free-pathℓ defines the average distance thechargecarriertravelsinamaterialbeforeitsenergy/direction is(substantially)alteredbya“scattering”event(e.g.,“colliding” with an impurity atom or phonon). In this sense,ℓ is directly relatedtotheabundanceoflatticedefects(e.g.,impurityatoms, vacancies, dislocations, grain bo...

work page

[5] [2]

effective

Magnetic Penetration Depth In brief, for a superconductor in its Meissner state, the magneticpenetrationdepth 𝜆describesthecharacteristiclength scale(typically10sto100sofnanometers)overwhichexternal magnetic flux may “creep” below the surface before being completely screened deeper in its interior. In the limit that the superconductor’selectrodynamicsarep...

work page

[6] [3]

correction

Non-linear Meissner Effect Empirically, it is well-known that𝜆(𝑇,𝑥) is affected by an applied magnetic field𝐵0, a result that was first understood theoreticallyusingGinzburg-Landau(GL)theory[ 40,41],with later developments extending it to unconventional supercon- ductors (see, e.g., [42]) Qualitatively, the effect accounts for 𝜆(𝑇,𝑥) ’s quadratic increase...

work page

[7] [4]

Meissner Screening Profile Keeping with our assumption of local electrodynamics (see Section IIC2), the “classic” London result [31] describing the spatial dependence of magnetic-flux-penetration below a Meissner state superconductor’s surface can be written in a generalized form so as to account for spatial inhomogeneities (e.g., from a depth-dependent i...

work page

[8] [5]

Supercurrent Density In 1D, the supercurrent density𝐽(𝑥) is related to𝐵(𝑥) by its first spatial derivative: 𝐽(𝑥)≡− 1 𝜇0 d d𝑥𝐵(𝑥),(31) where𝜇0 =1.256 637 061 27(20)N A −2 is the vacuum mag- netic permeability [19]. Ultimately, this is the most interesting quantity for SRF applications, where there is a strong desire to: minimize its value right at the surf...

work page