pith. sign in

arxiv: 2605.17667 · v1 · pith:QLYYLM5Wnew · submitted 2026-05-17 · ❄️ cond-mat.supr-con · cond-mat.mes-hall

I-V characteristics of SNS junctions with a multivalley normal region

Liam Bonds , Shiang-Bin Chiu , Anton Andreev , Boris Spivak This is my paper

Pith reviewed 2026-05-19 21:55 UTC · model grok-4.3

classification ❄️ cond-mat.supr-con cond-mat.mes-hall
keywords SNS junctionmultivalleyI-V characteristicssuperconducting kineticsrelaxation timesLarkin-Ovchinnikovnonmonotonic current
0
0 comments X

The pith

SNS junctions with a multivalley normal region show nonmonotonic I-V characteristics with two current peaks.

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

In conductors with multiple valleys, the inter-valley and inelastic relaxation times can be much longer than the momentum relaxation time. The paper shows that this timescale separation dramatically alters the current-voltage relation in SNS junctions that use such a material as the normal bridge. Under voltage bias the current is nonmonotonic and displays two peaks whose locations are set by the long relaxation times; these peaks can be substantially larger than the critical current. Under current bias the effective resistance in the nonlinear regime drops by orders of magnitude below the normal-state value because the slow processes control the dynamics.

Core claim

Generalizing the Larkin-Ovchinnikov kinetic equations to multivalley superconductors yields a description of SNS junctions in which, at fixed voltage V, I(V) exhibits two peaks of comparable height at V1 approximately ħ/(e τ_in) and V2 approximately ħ/(e τ_v). These peaks may greatly exceed Ic(T). At fixed current in the window Ic(T) < I ≲ I_jump the nonlinear resistance is set by the long relaxation times and can be orders of magnitude smaller than the normal-state resistance.

What carries the argument

Generalized Larkin-Ovchinnikov equations that track separate intra-valley momentum relaxation, inelastic relaxation, and inter-valley relaxation in the normal region of an SNS junction.

If this is right

  • The current-voltage curve is nonmonotonic with peaks at voltages inversely proportional to the inelastic and inter-valley relaxation times.
  • Peak currents can greatly exceed the equilibrium critical current of the junction.
  • In the current-biased regime above Ic but below the jump current, the differential resistance is much smaller than the normal-state resistance.
  • The kinetics are dominated by the separation of the three relaxation time scales.

Where Pith is reading between the lines

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

  • Similar nonmonotonic features may appear in other hybrid superconducting structures built on multivalley semiconductors or semimetals.
  • Experimental tests could involve fabricating SNS junctions on silicon or transition-metal dichalcogenides and measuring the positions of the current peaks versus estimated relaxation times.
  • Accounting for valley relaxation might be necessary when interpreting excess-current data in superconducting point contacts or junctions.
  • The effect provides a new handle for engineering low-resistance states in superconducting circuits via valley degree of freedom.

Load-bearing premise

The inter-valley relaxation time and the inelastic relaxation time are significantly longer than the intra-valley momentum relaxation time.

What would settle it

A measured I-V curve that is monotonic or whose peak voltages do not scale as the inverse of independently measured τ_in and τ_v would contradict the predicted nonmonotonicity and peak locations.

Figures

Figures reproduced from arXiv: 2605.17667 by Anton Andreev, Boris Spivak, Liam Bonds, Shiang-Bin Chiu.

Figure 1
Figure 1. Figure 1: Sketch of a multivalley SNS junction with an [PITH_FULL_IMAGE:figures/full_fig_p001_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Plot of average current I versus bias voltage V in an SNS junction with a multivalley normal metal. At bias voltages of the order eV1 ∼ τ −1 in and eV2 ∼ τ −1 v there exist two maxima of the average current; Imax,1, Imax,2. These maxima are of the same order and may greatly exceed the temperature-dependent critical current Ic(T). tegrating Eq. (10) over space and dividing it by ˜ν α(ϵ, χ). This yields (cf.… view at source ↗
Figure 3
Figure 3. Figure 3: I-V characteristics of current-biased SNS junctions obtained by solving the model in Appendix [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
read the original abstract

In multivalley conductors the inter-valley relaxation time $\tau_v$ and the inelastic relaxation time $\tau_{in}$ may be significantly longer than the intra-valley momentum relaxation time $\tau$. We show that this separation of time scales has dramatic effects on the I-V characteristics of SNS junctions with a multivalley normal region. We generalize the Larkin-Ovchinnikov equations describing superconducting kinetics to the case of multivalley superconductors. We use this generalization to obtain a kinetic description of multivalley SNS junctions. We find that at constant voltage bias $V$, the current $I(V)$ is nonmonotonic; it exhibits two peaks of similar magnitude $I_\text{max,1} \sim I_\text{max,2}$ at $V_1 \sim \hbar(e\tau_{in})^{-1}$ and $V_2\sim \hbar(e\tau_v)^{-1}$, which may greatly exceed the critical current $I_c(T)$. At constant current bias $I$ we find that in a wide interval, $I_c(T) < I \lesssim I_{\text{jump}}$, the nonlinear resistance of the junction is controlled by the long relaxation times and may be several orders of magnitude smaller than the normal state resistance.

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 / 2 minor

Summary. The manuscript generalizes the Larkin-Ovchinnikov kinetic equations to multivalley normal regions in SNS junctions. Under the assumption that inter-valley relaxation time τ_v and inelastic time τ_in are much longer than intra-valley momentum relaxation τ, the authors derive non-monotonic I(V) curves at constant voltage bias featuring two peaks of comparable height at V ~ ħ/(e τ_in) and V ~ ħ/(e τ_v), both potentially exceeding I_c(T). At constant current bias they report a wide interval where the nonlinear resistance is controlled by the long times and can be orders of magnitude below the normal-state value.

Significance. If the central derivation holds, the work identifies a concrete mechanism by which multivalley structure produces distinctive non-monotonic transport in superconducting weak links, with possible relevance to materials such as TMDs or multivalley semiconductors. The explicit generalization of the LO collision integrals supplies a reusable kinetic framework; the prediction of two distinct voltage scales tied to independent relaxation channels is falsifiable and could guide experiments.

major comments (2)
  1. [§2.3] §2.3, after Eq. (11): the collision integrals for inter-valley and inelastic scattering are written with independent phenomenological rates τ_v and τ_in; the manuscript does not derive the microscopic conditions (e.g., valley-orbit coupling strength or phonon matrix elements) under which τ_v, τ_in ≫ τ remains valid in a multivalley metal, leaving the separation of timescales as an assumption rather than a controlled limit.
  2. [§4] §4, Eq. (25) and the paragraph following Fig. 3: the claim that both peaks reach I_max,1 ≈ I_max,2 ≫ I_c(T) is obtained from the steady-state solution of the generalized kinetic equation; however, the ratio I_max/I_c is shown only for a narrow range of τ_v/τ_in values, and no error estimate or sensitivity analysis is provided when the separation is only moderate (τ_v/τ ≈ 10–30), which weakens the quantitative assertion that the peaks “may greatly exceed” I_c.
minor comments (2)
  1. Notation: the symbol τ is used both for the intra-valley momentum time and, in a few places, for a generic relaxation time; a consistent subscript (e.g., τ_m) would remove ambiguity.
  2. Figure 4 caption: the curves for different bias regimes are not labeled with the corresponding values of τ_in/τ_v; adding this information would make the comparison to the analytic limits clearer.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the positive evaluation of the significance of our work and for the detailed and constructive major comments. We respond to each point below, indicating the revisions we will make to the manuscript.

read point-by-point responses

  1. Referee: [§2.3] §2.3, after Eq. (11): the collision integrals for inter-valley and inelastic scattering are written with independent phenomenological rates τ_v and τ_in; the manuscript does not derive the microscopic conditions (e.g., valley-orbit coupling strength or phonon matrix elements) under which τ_v, τ_in ≫ τ remains valid in a multivalley metal, leaving the separation of timescales as an assumption rather than a controlled limit.

    Authors: We agree that the separation τ_v, τ_in ≫ τ is introduced as a phenomenological assumption rather than derived from microscopic parameters such as valley-orbit coupling or phonon matrix elements. This is consistent with the standard approach in kinetic-equation treatments of nonequilibrium superconductivity, where the focus is on the consequences of timescale separation once it is established. The assumption is motivated by experimental observations in multivalley materials (e.g., TMDs), where momentum mismatch suppresses inter-valley scattering. In the revised manuscript we will expand the discussion immediately after Eq. (11) to cite representative literature on measured relaxation times and to state explicitly the regime of validity we assume. A first-principles calculation of the microscopic rates lies beyond the scope of this kinetic-theory paper. revision: partial

  2. Referee: [§4] §4, Eq. (25) and the paragraph following Fig. 3: the claim that both peaks reach I_max,1 ≈ I_max,2 ≫ I_c(T) is obtained from the steady-state solution of the generalized kinetic equation; however, the ratio I_max/I_c is shown only for a narrow range of τ_v/τ_in values, and no error estimate or sensitivity analysis is provided when the separation is only moderate (τ_v/τ ≈ 10–30), which weakens the quantitative assertion that the peaks “may greatly exceed” I_c.

    Authors: The referee is correct that the presented numerical results for the peak currents are limited to a specific parameter set. We will revise the paragraph following Fig. 3 and add a new inset or supplementary figure that displays I_max,1/I_c and I_max,2/I_c as functions of τ_v/τ and τ_in/τ over a wider interval, explicitly including moderate separations (τ_v/τ ≈ 10–30). A brief sensitivity discussion will be included to quantify how the excess over I_c(T) degrades when the timescale separation is reduced. These additions will be incorporated in the revised version. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation follows from generalized kinetic equations under stated timescale separation

full rationale

The paper takes the separation τ_v, τ_in ≫ τ as a physical premise for multivalley metals and generalizes the standard Larkin-Ovchinnikov kinetic equations to this multivalley setting. The non-monotonic I(V) with peaks at voltages ħ/(e τ_in) and ħ/(e τ_v) is obtained by solving the resulting collision integrals and kinetic equations under constant voltage or current bias. No equation or result is shown to reduce to a fitted parameter, a self-referential definition, or a load-bearing self-citation whose validity itself depends on the target I-V features. The central claims therefore remain independent of the inputs by construction.

Axiom & Free-Parameter Ledger

2 free parameters · 1 axioms · 0 invented entities

The central claim rests on the validity of generalizing the Larkin-Ovchinnikov kinetic equations to multivalley systems and on the physical assumption that inter-valley and inelastic times are much longer than the momentum relaxation time. No explicit free parameters are fitted in the abstract description.

free parameters (2)
  • τ_v
    Inter-valley relaxation time introduced as significantly longer than intra-valley τ; sets one of the voltage scales for the peaks.
  • τ_in
    Inelastic relaxation time introduced as significantly longer than intra-valley τ; sets the other voltage scale for the peaks.
axioms (1)
  • domain assumption Larkin-Ovchinnikov equations can be generalized to multivalley superconductors while preserving their kinetic structure.
    Invoked to obtain the kinetic description of the multivalley SNS junction.

pith-pipeline@v0.9.0 · 5766 in / 1540 out tokens · 77831 ms · 2026-05-19T21:55:55.638764+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

45 extracted references · 45 canonical work pages

  1. [1]

    Voltage bias It is shown below that the I-V characteristic of voltage- biased junctions turns out to be nonmonotonic, with two pronounced peaks of similar heightI max,1 ∼I max,2, as illustrated in Fig. 2. The first peak is reached ateV∼ τ −1 in , and the second ateV∼τ −1 v . We will show that the peak currentsI max,1 ∼I max,2 can be significantly larger t...

    work page

  2. [2]

    (28) that only the quasiparticles in energy levels within a narrow windowϵ∼E T ≪T contribute to the current (below, we refer to them as “active” levels)

    Current bias It follows from Eq. (28) that only the quasiparticles in energy levels within a narrow windowϵ∼E T ≪T contribute to the current (below, we refer to them as “active” levels). This observation enables us to introduce a simple model of the junction, which describes the I- V characteristics of the junction with accuracy of order unity. In this mo...

    work page

  3. [3]

    High bias regime It is clear from Fig. 3 that even in the high bias regime, where the inelastic and inter-valley relaxation become ineffective, the nonequilibrium population of active quasiparticle levels significantly affects the average voltage across the current-biased junction. In the high bias regime,I≫G N /eτin, the solution of the model simplifies ...

    work page

  4. [4]

    E. I. Blount, Phys. Rev.114, 418 (1959)

    work page 1959

  5. [5]

    Adler, IBM J

    E. Adler, IBM J. Res. Dev.8, 430 (1964)

    work page 1964

  6. [6]

    S. V. Gantsevich and V. L. Gurevich, Phys. Rev.161, 736 (1967)

    work page 1967

  7. [7]

    E. I. Rashba, Z. S. Gribnikov, and V. Y. Kravchenko, Phys. Usp.19, 361 (1976)

    work page 1976

  8. [8]

    Takayanagi and T

    H. Takayanagi and T. Kawakami, Phys. Rev. Lett.54, 2449 (1985)

    work page 1985

  9. [9]

    A. W. Kleinsasser, T. N. Jackson, D. McInturff, F. Rammo, G. D. Pettit, and J. M. Woodall, Applied Physics Letters55, 1909 (1989)

    work page 1909

  10. [10]

    Akazaki, H

    T. Akazaki, H. Takayanagi, J. Nitta, and T. Enoki, Ap- plied Physics Letters68, 418 (1996)

    work page 1996

  11. [11]

    Giazotto, K

    F. Giazotto, K. Grove-Rasmussen, R. Fazio, F. Beltram, E. H. Linfield, and D. A. Ritchie, J. Supercond.17, 317 (2004)

    work page 2004

  12. [12]

    Shabani, M

    J. Shabani, M. Kjaergaard, H. J. Suominen, Y. Kim, F. Nichele, K. Pakrouski, T. Stankevic, R. M. Lutchyn, P. Krogstrup, R. Feidenhans’l, S. Kraemer, C. Nayak, M. Troyer, C. M. Marcus, and C. J. Palmstrøm, Phys. Rev. B93, 155402 (2016)

    work page 2016

  13. [13]

    Kjærgaard, F

    M. Kjærgaard, F. Nichele, H. J. Suominen, M. P. Nowak, M. Wimmer, A. R. Akhmerov, J. A. Folk, K. Flensberg, J. Shabani, C. J. Palmstrøm, and C. M. Marcus, Nat. Commun.7, 12841 (2016)

    work page 2016

  14. [14]

    Kjaergaard, H

    M. Kjaergaard, H. J. Suominen, M. P. Nowak, A. R. Akhmerov, J. Shabani, C. J. Palmstrøm, F. Nichele, and C. M. Marcus, Phys. Rev. Appl.7, 034029 (2017)

    work page 2017

  15. [15]

    C. G. L. Bøttcher, F. Nichele, M. Kjaergaard, H. J. Suominen, J. Shabani, C. J. Palmstrøm, and C. M. Mar- cus, Nat. Phys.14, 1138 (2018)

    work page 2018

  16. [16]

    Nishino, E

    T. Nishino, E. Yamada, and U. Kawabe, Phys. Rev. B 33, 2042(R) (1986)

    work page 2042

  17. [17]

    Hatano, T

    M. Hatano, T. Nishino, F. Murai, and U. Kawabe, Ap- plied Physics Letters53, 409 (1988)

    work page 1988

  18. [18]

    Chiodi, J.-E

    F. Chiodi, J.-E. Duvauchelle, C. Marcenat, D. D´ ebarre, and F. Lefloch, Phys. Rev. B96, 024503 (2017)

    work page 2017

  19. [19]

    N. W. Hendrickx, M. L. V. Tagliaferri, M. Kouwenhoven, R. Li, D. P. Franke, A. Sammak, A. Brinkman, G. Scap- pucci, and M. Veldhorst, Phys. Rev. B99, 075435 (2019)

    work page 2019

  20. [20]

    Aggarwal, A

    K. Aggarwal, A. Hofmann, D. Jirovec, I. Prieto, A. Sam- mak, M. Botifoll, S. Mart´ ı-S´ anchez, M. Veldhorst, J. Ar- biol, G. Scappucci, J. Danon, and G. Katsaros, Phys. Rev. Res.3, L022005 (2021)

    work page 2021

  21. [21]

    Tosato, V

    A. Tosato, V. Levajac, J.-Y. Wang, C. J. Boor, F. Borsoi, M. Botifoll, C. N. Borja, S. Mart´ ı-S´ anchez, J. Arbiol, A. Sammak, M. Veldhorst, and G. Scappucci, Commun. Mater.4, 23 (2023)

    work page 2023

  22. [22]

    Valentini, O

    M. Valentini, O. Sagi, L. Baghumyan, T. de Gi- jsel, J. Jung, S. Calcaterra, A. Ballabio, J. Aguil- era Servin, K. Aggarwal, M. Janik, T. Adletzberger, R. Seoane Souto, M. Leijnse, J. Danon, C. Schrade, E. Bakkers, D. Chrastina, G. Isella, and G. Katsaros, Nat. Commun.15, 169 (2024)

    work page 2024

  23. [23]

    A. I. Larkin and Y. N. Ovchinnikov, Sov. Phys. JETP 24, 1035 (1967)

    work page 1967

  24. [24]

    Tinkham,Introduction to superconductivity(Courier Corporation, 2004)

    M. Tinkham,Introduction to superconductivity(Courier Corporation, 2004)

    work page 2004

  25. [25]

    S. N. Artemenko, A. F. Volkov, and A. V. Zaitsev, Zh. Eks. Teor. Fiz.76, 1816 (1979)

    work page 1979

  26. [26]

    Averin and A

    D. Averin and A. Bardas, Phys. Rev. B53, R1705(R) (1996)

    work page 1996

  27. [27]

    T. Liu, A. V. Andreev, and B. Z. Spivak, Ann. Phys.456 (2023)

    work page 2023

  28. [28]

    T. Liu, M. Smith, A. V. Andreev, and B. Z. Spivak, Phys. Rev. B109, L020501 (2024)

    work page 2024

  29. [29]

    A. I. Larkin and Y. N. Ovchinnikov, Zh. Eksp. Teor. Fiz. 73, 7 (1977)

    work page 1977

  30. [30]

    J. B. Nielsen and H. Smith, Phys. Rev. B31, 2831 (1985)

    work page 1985

  31. [31]

    A. L. Shelankov, J. Low Temp. Phys.60, 29 (1985)

    work page 1985

  32. [32]

    In this case, one can use the Boltzmann equation for the quasiparticle dis- tribution functionn α(p) [42]

    The physical meaning of the two distribution functions may be understood by considering relatively clean su- perconductors, with ∆τ≫1. In this case, one can use the Boltzmann equation for the quasiparticle dis- tribution functionn α(p) [42]. In the diffusive regime, the distribution function becomes isotropicn α(p)→ nα(ξp). The quasiparticle energy isϵ= p...

    work page

  33. [33]

    A. I. Larkin and Y. N. Ovchinnikov, inNonequilibrium Superconductivity, edited by D. N. Langenberg and A. I. Larkin (Elsevier, Amsterdam, 1986) Chap. 11

    work page 1986

  34. [34]

    (17) reproduces the expression for equilibrium supercurrent

    Forf α 1 = 0 andf α 0 (ϵ) =f F (ϵ) Eq. (17) reproduces the expression for equilibrium supercurrent. We note that in a normal metalK 1 = 4, so the expression for the cur- 14 rent isJ α =−eν α N D R ∞ 0 dϵ∇f α 1 =eν α N D R ∞ −∞ dϵ∇αδn. This is consistent with the correspondence off α 1 to the particle-hole branch imbalance, see the footnote above Eq. (10)

    work page

  35. [35]

    (19) does not preclude deviations of the elec- tron densityn α in a given valley from equilibrium

    In multivalley conductors, the charge neutrality condi- tion Eq. (19) does not preclude deviations of the elec- tron densityn α in a given valley from equilibrium. It follows from Eqs. (18) and (19) that in the normal state, the electric field is related tof α 1 (ϵ) byeE= 1P α να N P α R ∞ 0 να N ∇f α 1 (ϵ)dϵ

    work page

  36. [36]

    H. Suhl, B. T. Matthias, and L. R. Walker, Phys. Rev. Lett.3, 552 (1959)

    work page 1959

  37. [37]

    F. Zhou, P. Charlat, B. Spivak, and B. Pannetier, J. Low Temp. Phys.110, 841 (1998)

    work page 1998

  38. [38]

    Smith, A

    M. Smith, A. V. Andreev, and B. Z. Spivak, Phys. Rev. B101, 134508 (2020)

    work page 2020

  39. [39]

    Smith, A

    M. Smith, A. V. Andreev, and B. Z. Spivak, Annals of Physics417, 168105 (2020)

    work page 2020

  40. [40]

    The upper limit in voltage is determined by the condition G(V)/G N ∼1

    work page

  41. [41]

    J. E. Moyal, inMathematical Proceedings of the Cam- bridge Philosophical Society, Vol. 45 (Cambridge Univer- sity Press, 1949) pp. 99–124

    work page 1949

  42. [42]

    H. J. Groenewold, inOn the principles of elementary quantum mechanics(Springer, 1946) pp. 1–56

    work page 1946

  43. [43]

    G. A. Baker, Phys. Rev.109, 2198 (1958)

    work page 1958

  44. [44]

    N. B. Kopnin,Theory of nonequilibrium superconductiv- ity, Vol. 110 (Oxford University Press, 2001)

    work page 2001

  45. [45]

    A. G. Aronov, Y. M. Gal’perin, V. L. Gurevich, and V. I. Kozub, Adv. Phys.30, 539 (1981)

    work page 1981