REVIEW 3 major objections 4 minor 34 references
Switchable Josephson current in junctions with spin-orbit coupling
T0 review · 3 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read Spin-orbit coupling switches Josephson current from 0 to pi.
desk verdict A solid Usadel-based paper with one clean analytic result and an attractive but numerically supported gate-switching claim, held back by an unquantified thin-layer reduction. 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 linearized Usadel equation for the anomalous Green's function, generalized to linear-in-momentum spin-orbit coupling via the SU(2) covariant derivative $\tilde\nabla_k=\partial_k-i[\hat A_k,\cdot]$. The spin-orbit coupling enters through two averaged tensors obtained after integrating the thin layers along the transverse direction: the spin-precession tensor $\bar C^{ab}_k=\varepsilon^{acb} A^c_k d/(W+d)$ and the Dyakonov-Perell spin-relaxation tensor $\bar\Gamma^{ab}$. The argument proceeds by expanding the singlet and triplet components in powers of the spin-orbit coupling and showing that the precession term generates the long-range component $\tilde f^z_t$ at first order, while the relaxation term generates $\tilde f^y_t$ at second order; their competition fixes the sign of the critical current in Eq. (31).
What would settle it
Measure the critical current of a symmetric type 1 junction with pure Rashba spin-orbit coupling as a function of the exchange-field angle $\vartheta$; the theory predicts $j_c>0$ for all $\vartheta$, so any observed sign change would falsify the central claim that both Rashba and Dresselhaus coupling are required for reversal in such junctions.
Extended reading notes
Core claim
The central claim is that the long-range triplet Josephson current in long lateral diffusive junctions is switchable in sign, and the switch is controlled by the direction of the exchange field and by the Rashba/Dresselhaus composition and strength of the spin-orbit coupling. The paper derives, to second order in the spin-orbit coupling for a symmetric type 1 junction, the critical current $j_c \propto \left[(\bar\alpha\cos\vartheta+\bar\beta\sin\vartheta)^2/(2\kappa_\omega) - 8\bar\alpha^2\bar\beta^2\cos^2 2\vartheta/\kappa_\omega^3\right]$, where the first term comes from spin precession (positive, favoring the 0 state) and the second from anisotropic spin relaxation (negative, favoring the pi state). The reversal therefore requires both Rashba and Dresselhaus coupling in a symmetric type 1 junction, and it disappears when either coupling is zero. In type 2 junctions, numerical solution of the same z-averaged Usadel equations shows that a 0-pi transition already occurs for a pure Rashba or Dresselhaus coupling as its strength is increased, and the range of parameters supporting the pi state is wider.
Load-bearing premise
The one-dimensional description of the junction relies on the ferromagnetic bridge and the spin-orbit-coupled interlayers being thin enough that the superconducting correlations are nearly constant across the thickness, so the transverse integration in Eq. (11) captures the physics; if this film-thickness assumption is violated, the quantitative current-angle curves and the 0-pi boundaries would shift even if the qualitative mechanism survived.
Editorial extensions
If this is right
- In symmetric type 1 junctions, rotating an in-plane exchange field switches the current sign only when both Rashba and Dresselhaus spin-orbit coupling are present; the current stays non-negative for pure Rashba or pure Dresselhaus coupling.
- In type 2 junctions, a gate voltage changing the Rashba strength alone can drive the junction from a 0 to a pi ground state, making a voltage-controlled supercurrent switch feasible in semiconducting bridges.
- Asymmetric type 1 junctions with different spin-orbit coupling in the two lead interfaces can show current reversal from spin precession alone, for example when one interface has only Rashba and the other only Dresselhaus coupling.
- A single lateral junction can serve as a detector of long-range triplet correlations: the sign of the critical current as the field direction is rotated provides a way to identify the presence of such correlations in a device.
Reading between the lines
- The same precession-versus-relaxation competition likely appears in other spin-orbit-coupled hybrid junctions, suggesting that 0-pi switching could serve as a direct probe of the relative strength of spin precession and Dyakonov-Perell relaxation in a material, beyond the two geometries studied here.
- If the thin-layer z-averaging holds, the type 2 prediction implies a quantitative threshold in the gate-tunable Rashba parameter; measuring the current-angle curve below and above that threshold would give a clean experimental test of the theory.
- The framework could be extended to time-dependent fields: because the sign reversal is controlled by orientation, a rotating exchange field may pump a switching supercurrent, a route the paper does not explore.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies two lateral diffusive Josephson junctions with Rashba and Dresselhaus spin-orbit coupling (SOC) in the presence of an in-plane exchange field. In type-1 junctions the SOC is confined to heavy-metal interlayers under the superconducting electrodes; in type-2 junctions the SOC is present in the ferromagnetic bridge. In both cases the junction length is assumed to exceed the magnetic decay length, so that the Josephson current is carried by the long-range triplet component. The authors derive linearized Usadel equations after averaging over the layer thickness, obtain an analytic perturbative expression for the critical current in type-1 junctions to leading and next-to-leading order in the SOC strength, and present numerical results for both junction types. The central claims are that the critical current can be controlled by rotating the exchange field, and that 0–π transitions can be induced by the competition between spin precession and spin relaxation; in type-2 junctions, a 0–π transition can be induced by a pure Rashba or Dresselhaus SOC when its strength is increased.
Significance. If the central claims hold, the results are of clear experimental relevance: they predict a voltage-gate-controlled 0–π transition in a diffusive lateral junction, and they provide a single-junction diagnostic for long-range triplet correlations. The paper has genuine strengths. The analytic perturbative solution in Section III is internally consistent and is explicitly compared with the small-SOC numerical curves, with no fitting parameters. The distinction between spin-precession and spin-relaxation mechanisms is clearly developed, and the predicted dependence on the exchange-field orientation is specific and falsifiable. The main limitations concern the validation of the one-dimensional reduction that underlies all quantitative results, and the type-2 large-SOC regime, which is supported only by numerical solutions of the same reduced model.
major comments (3)
- [Section II.A and Appendix A, Eqs. (12)–(13) and (A5)–(A6)] The z-averaged exchange field hbar = hW/(W+d) is used in the bridge region for the type-1 junction, but the SOC interlayer of thickness d exists only under the superconducting electrodes; the true z-average in the bridge, |x|<L/2, should be h, not hbar. With d/W = 1 in Fig. 2, this is not a small correction: the short-range component and the boundary conditions entering the analytic current, Eq. (31), are computed with an exchange field in the bridge that is too small by the factor W/(W+d). Please either introduce an x-dependent hbar(x) (hbar in the electrodes, h in the bridge) or state and justify the d≪W limit in which the discrepancy is negligible.
- [Section II, Eq. (11), and Figs. 2–5] The reduction of the two-dimensional Usadel problem to the averaged one-dimensional equations is justified by the assumption that d and W are small compared with the scale over which the anomalous Green function varies, but no quantitative check is provided. The numerical calculations use d/W = 1 and L = 5ξ0 without stating the ratio (W+d)/ξ0, so the reader cannot tell whether the averaging condition is satisfied in the parameter regime where the 0–π boundaries are determined. A full two-dimensional numerical solution for representative parameters, or an explicit estimate of the neglected O((W+d)^2) terms, would establish that the computed sign-change boundaries are not an artifact of the one-dimensional reduction.
- [Section IV, Figs. 3(d)–(f) and 5] The most novel claim, that type-2 junctions allow a 0–π transition for pure Rashba or Dresselhaus SOC at large SOC strength, is supported only by numerical integration of the same reduced one-dimensional equations, without an analytic large-SOC limit or a convergence study. The analytic result in Eq. (31) is derived for small SOC and cannot validate this regime. Please add either an analytic large-SOC estimate or additional numerical evidence, such as agreement with the one-dimensional pure-Rashba result of Ref. 27 or a check of numerical convergence, to make this claim robust.
minor comments (4)
- [Abstract] There is a stray comma before the period in 'tuning the strength of the spin-orbit coupling in type 2 junctions., and also discuss'.
- [Section III heading area] The text begins 'Her we focus' in the sentence before Eq. (29); this should be 'Here'.
- [Fig. 3 caption] The caption reads 'for in junction of type 2'; this should be 'for a junction of type 2'.
- [Section IV, paragraph after Fig. 2] The phrase 'The case whenα,β ⁄= 0' lacks spacing; it should read 'The case when α,β ≠ 0'.
Circularity Check
No fitted-input circularity: the predicted 0-pi transitions are outputs of the Usadel equations, not inputs.
full rationale
The derivation chain is self-contained from the stated model. Equation (1) and the boundary condition (3) are taken from Refs. 10 and 11, which include one coauthor of the present paper, but these are the standard linearized Usadel equations generalized to linear-in-momentum spin-orbit coupling; they do not already contain the direction-or-strength dependence of the Josephson current sign that is claimed. The z-averaging reduction in Eq. (11) is performed explicitly under a stated thin-layer assumption rather than imported as an unexplained ansatz, and it is this averaged set of equations, (12)-(13) and (24)-(25), that is solved. The analytical current in Eq. (31) is obtained perturbatively from those equations, with the positive spin-precession and negative spin-relaxation terms both derived rather than imposed. The numerical calculations solve the same equations with fixed model parameters (e.g. h = 10 Delta, L = 5 xi0, T = 0.01 Delta) and no fitting, and the 0-pi boundaries as functions of theta or SOC strength are outputs. Comparisons with Refs. 27 and 28 occur after the results are obtained and are not used to set parameters or select the final expressions. The only mild self-citation is the underlying SU(2)/Usadel formalism, which is load-bearing but does not encode the target sign-reversal result; the thin-layer validity concern is a modeling and correctness risk, not a circularity.
Assumptions & free parameters
free parameters (7)
- Rashba SOC strength α =
scanned from 0.01 to 0.5 in units of 1/ξ0 in Figs. 2, 3 and 5
- Dresselhaus SOC strength β =
0, 0.01 to 0.5 in units of 1/ξ0 in figures
- Exchange field orientation ϑ =
scanned over [0, π] in figures
- Exchange field magnitude h =
10 Δ in the numerical figures
- Junction length L =
5 ξ0 in the numerical figures
- Temperature T =
0.01 Δ in the numerical figures
- SOC-layer to bridge thickness ratio d/W =
1 in the type 1 numerical figures
assumptions (4)
- domain assumption The diffusive limit holds and the linearized Usadel equation, Eq. (1), generalized to linear-in-momentum SOC, describes the induced correlations in the non-superconducting region.
- domain assumption The SOC interlayer thickness d and bridge thickness W are small compared with the scale over which the anomalous Green function varies, allowing z-integration to reduce the 2D problem to a 1D problem.
- domain assumption The junction length satisfies L >> ξ_h, so short-range triplet components decay in the bridge and the Josephson current is carried only by the long-range triplet components.
- domain assumption For the analytic type 1 solution, the exchange interaction is the dominant energy scale, h >> Dα², Dβ², Dαβ, T, and the current is computed up to second order in the SOC fields.
Cite this review
Pith. "Pith review of Switchable Josephson current in junctions with spin-orbit coupling." pith.science (2026). https://pith.science/paper/OJIARHV5
@misc{pith2026190900401,
author = {Pith},
title = {Pith review of: Switchable Josephson current in junctions with spin-orbit coupling},
year = {2026},
howpublished = {\url{https://pith.science/paper/OJIARHV5}},
note = {Machine review of arXiv:1909.00401}
}
abstract
We study the Josephson current in two types of lateral junctions with spin-orbit coupling and an exchange field. The first system (type 1 junction) consists of superconductors with heavy metal interlayers linked by a ferromagnetic bridge, such that the spin-orbit coupling is finite only at the superconductor/heavy metal interface. In the second type (type 2) of system we assume that the spin orbit coupling is finite in the bridge region. The length of both junctions is larger than the magnetic decay length such that the Josephson current is carried uniquely by the long-range triplet component of the condensate. The latter is generated by the spin-orbit coupling via two mechanisms, spin precession and inhomogeneous spin-relaxation. We show that the current can be controlled by rotating the magnetization of the bridge or by tuning the strength of the spin-orbit coupling in type 2 junctions., and also discuss how the ground-state of the junction can be tuned from a $0$ to a $\pi$ phase difference between the superconducting electrodes. In leading order in the spin-orbit coupling, the spin precession dominates the behavior of the triplet component and both junctions behave similarly. However, when spin relaxation effects are included junction of type 2 offers a wider parameter range in which $0$-$\pi$ transitions take place.
Figures
Reference graph
Works this paper leans on
-
[1]
write newline
" write newline "" before.all 'output.state := FUNCTION fin.entry add.period write newline FUNCTION new.block output.state before.all = 'skip after.block 'output.state := if FUNCTION new.sentence skip FUNCTION not #0 #1 if FUNCTION and 'skip pop #0 if FUNCTION or pop #1 'skip if FUNCTION new.block.checka empty 'skip 'new.block if FUNCTION new.block.checkb...
-
[2]
F. S. Bergeret, A. F. Volkov, and K. B. Efetov, Phys. Rev. Lett. 86 , 4096 (2001)
work page 2001
-
[3]
F. S. Bergeret, A. F. Volkov, and K. B. Efetov, Rev. Mod. Phys. 77 , 1321 (2005), 0506047
work page 2005
-
[4]
Spin-polarized supercurrents for spintronics: a review of current progress
M. Eschrig, Phys. Today 64 , 43 (2011), 1509.02242
work page Pith review arXiv 2011
- [5]
- [6]
-
[7]
S. V. Panyukov, L. N. Bulaevskil, and A. I. Buzdin, Sov. Phys. JETP 35 , 178 (1982)
work page 1982
-
[8]
V. V. Ryazanov et al. , Phys. Rev. Lett. 86 , 2427 (2001), 0008364
work page 2001
Show all 34 references
-
[9]
A. F. Volkov, K. B. Efetov, and F. S. Bergeret, Phys. Rev. B 64 , 1 (2001)
2001
-
[10]
Kontos et al
T. Kontos et al. , Physical Review Letters 89 , 1 (2002)
2002
-
[11]
F. S. Bergeret and I. V. Tokatly, Phys. Rev. Lett. 110 , 1 (2013), 1211.3084
2013 arXiv
-
[12]
F. S. Bergeret and I. V. Tokatly, Phys. Rev. B 89 , 134517 (2014), 1402.1025
2014 arXiv
-
[13]
Robinson, J
J. Robinson, J. Witt, and M. Blamire, Science 329 , 59 (2010)
2010
-
[14]
M. S. Anwar, F. Czeschka, M. Hesselberth, M. Porcu, and J. Aarts, Phys. Rev. B 82 , 2 (2010)
2010
-
[15]
X. L. Wang et al. , Phys. Rev. B 89 , 3 (2014)
2014
-
[16]
E. C. Gingrich et al. , Phys. Rev. B 86 , 1 (2012)
2012
-
[17]
J. W. A. Robinson, F. Chiodi, M. Egilmez, G. B. Hal \' a sz, and M. G. Blamire, Sci. Rep. 2 , 699 (2012)
2012
-
[18]
Chiodi et al
F. Chiodi et al. , EPL 101 , 37002 (2013)
2013
-
[19]
A. Pal, Z. Barber, J. Robinson, and M. Blamire, Nat. Commun. 5 , 3340 (2014)
2014
-
[20]
J. W. A. Robinson, N. Banerjee, and M. G. Blamire, Phys. Rev. B 89 , 104505 (2014)
2014
-
[21]
Kalcheim, O
Y. Kalcheim, O. Millo, M. Egilmez, J. W. A. Robinson, and M. G. Blamire, Phys. Rev. B 85 , 1 (2012)
2012
-
[22]
Banerjee et al
N. Banerjee et al. , Nat. Commun. 5 , 3048 (2014)
2014
-
[23]
T. S. Khaire, M. A. Khasawneh, W. P. Pratt, and N. O. Birge, Phys. Rev. Lett. 104 , 137002 (2010)
2010
-
[24]
Satchell and N
N. Satchell and N. O. Birge, Phys. Rev. B 97 , 1 (2018)
2018
-
[25]
Banerjee et al
N. Banerjee et al. , Phys. Rev. B 97 , 1 (2018)
2018
-
[26]
Satchell, R
N. Satchell, R. Loloee, and N. O. Birge, Phys. Rev. B 99 , 1 (2019), arXiv:1904.08798v2
2019 arXiv
-
[27]
X. Liu, J. K. Jain, and C. X. Liu, Phys. Rev. Lett. 113 , 1 (2014)
2014
-
[28]
Arjoranta and T
J. Arjoranta and T. T. Heikkil \" a , Phys. Rev. B 93 , 1 (2016)
2016
-
[29]
J. R. Eskilt, M. Amundsen, N. Banerjee, and J. Linder, arxiv (2019), 1906.07725
2019 arXiv
-
[30]
Usadel, Phys
K. Usadel, Phys. Rev. Lett. 25 , 507 (1970)
1970
-
[31]
Eq. ( Usadel_basic ) is written in the strict diffusive limit and do not take into account charge-spin conversion terms which are higher order in the momentum relaxation rate Bergeret2015,Konschelle2015
-
[32]
M. Y. Kupriyanov and V. F. Lukichev, Sov. Phys. JETP 67 , 1163 (1988)
1988
-
[33]
Konschelle, I
F. Konschelle, I. V. Tokatly, and F. S. Bergeret, Phys. Rev. B 92 (2015), 1506.02977
2015 arXiv
-
[34]
F. S. Bergeret and I. V. Tokatly, Epl 110 (2015)
2015
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.