Pith. sign in

REVIEW 4 major objections 5 minor 19 references

Resonant escape in Josephson tunnel junctions under millimeter-wave irradiation

T0 review · 4 major / 5 minor · reviewed 2026-08-12 · deepseek-v4-flash

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

Pith's one-line read Under 100–110 GHz irradiation, the switching-current distribution of Nb/Al-AlOx/Nb Josephson junctions splits into two peaks, a signature of resonant escape of the Josephson phase that the strong-driving model reproduces quantitatively.

desk verdict Real double-peak mm-wave data, but the quantitative fit to the strong-driving model has an unaddressed 7 GHz plasma-frequency discrepancy. read the letter →

arxiv 2411.15048 v1 pith:TSWTZ74R submitted 2024-11-22 cond-mat.supr-con

J. N. Kämmerer (1) , S. Masis (1) , K. Hambardzumyan (1) , P. Lenhard (1) , U. Strobel (1) , J. Lisenfeld (1) , H. Rotzinger (1 , 2)
show 5 more authors
A. V. Ustinov (1 2) ((1) Physikalisches Institut Karlsruher Institut für Technologie (2) Institut für Quantenmaterialien und Technologie Karlsruher Institut für Technologie)
This is my paper · ORCID
classification cond-mat.supr-con
keywords Josephsontunneljunctionsresonantescapeswitchingcurrentdistributionmillimeter-waveirradiationplasmafrequencyphasequbitsstrong-drivingmodelNb/Al-AlOx/Nb
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 tries to establish that the switching statistics of a current-biased Nb/Al-AlOx/Nb Josephson junction under millimeter-wave irradiation carry a clear fingerprint of resonant phase escape: as the radiation power increases, the single peak in the switching-current distribution splits into a primary peak and a resonant peak at lower current. The authors show that the power dependence of these two branches at 106–110 GHz is quantitatively captured by the strong-driving model of effective barrier suppression, with a single shared plasma frequency ωp/2π = 145.5 GHz and quality factors Q ≈ 90. This matters because junctions with plasma frequencies near 140 GHz have been difficult to fabricate and characterize, and operation of phase qubits near 100 GHz would relax cooling requirements well beyond the current millikelvin standard. The measurements are performed at 4.2 K in the thermally activated regime, so the work is a feasibility proof for mm-wave phase qubits rather than a demonstration of quantum-coherent operation.

What carries the argument

The central object is the Josephson phase difference δ, treated as the coordinate of a fictitious particle in the tilted-washboard potential of the resistively and capacitively shunted junction model. Escape over the current-reduced barrier is enhanced by radiation when the drive frequency matches the small-oscillation frequency ω0(I), and the analysis rests on the strong-driving limit (ω/ωp)^5 ≫ ℏωp/EJ (about 170 here), where the escape is governed by a radiation-induced suppression of the effective potential barrier. The load-bearing identity is Eq. (8), a transcendental equation for the average switching-current shift ⟨δIsw(P)⟩ whose multiple solutions produce the coexisting primary and resonant peaks; in the harmonic approximation the sum reduces to a single resonant term (Eq. (10)) with the phase matrix element of the harmonic oscillator.

What would settle it

An independent measurement of the plasma frequency that disagrees with the fitted 145.5 GHz by more than the stated uncertainty, for example direct microwave spectroscopy of the same junction as a qubit at base temperature, would falsify the identification; likewise, if only the primary peak appeared at higher drive power, or if fitting each frequency required a different plasma frequency rather than one common value, the strong-driving model's claim would fail.

Watch

Extended reading notes

Core claim

The central claim is that the double-peak structures observed in the switching-current distributions are not noise artifacts but resonant escape of the Josephson phase from a metastable well: the drive at 100–110 GHz enhances the escape rate when the small-oscillation frequency ω0(I) = ωp(1 − (I/Ic)²)^{1/4} matches the irradiation frequency, producing a second, power-dependent peak at lower current. Fitting the resonant-peak position versus frequency gives ωp/2π = 138.3(7) GHz, close to the 140 GHz design target. The full power dependence at three frequencies is fitted by the strong-driving model of Eq. (8), using Ic = 242 μA and T = 4.2 K, with a common plasma frequency of 145.5 GHz and effective quality factors of about 90. The authors take this agreement as evidence that low-loss Nb/Al-AlOx/Nb junctions with plasma frequencies around 140 GHz can be fabricated, and that the same switching-current measurement provides a rapid characterization route for mm-wave qubit development.

Load-bearing premise

The load-bearing premise is that radiation suppresses the escape barrier in the way the paper's strong-driving model describes, with a known plasma frequency and damping; if that simplified model is wrong, the inferred plasma frequency and quality factor do not follow from the data, especially since the microwave power axis is in arbitrary units and the coupling is fit rather than independently measured.

Editorial extensions

If this is right

  • If the strong-driving model is right, the position of the resonant peak as a function of irradiation frequency is a direct measure of the junction plasma frequency, so switching-statistics measurements can be used as a fast, non-destructive diagnostic for high-current-density junction fabrication.
  • Junctions with plasma frequencies around 140 GHz are within reach of the standard Nb/Al-AlOx/Nb trilayer process, which removes a material barrier to phase qubits operating near 100 GHz.
  • Operating a phase qubit near 100 GHz raises the temperature scale T0 = hf/kB to about 4.8 K, suggesting that qubit operation at 1 K or above is not ruled out by level-spacing requirements.
  • The agreement of the fitted quality factors Q ≈ 90 at 4.2 K with quasiparticle-limited damping implies the same junctions should be substantially less dissipative at millikelvin temperatures, where coherent mm-wave qubit operation would be attempted.
  • The multivalued solutions of Eq. (8) explain the coexistence of primary and resonant peaks, giving a predictive handle on the microwave power levels needed to resonantly switch a junction in applications.

Reading between the lines

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

  • One extension the authors leave implicit is that the harmonic approximation keeps only the 0→1 transition, so the fit cannot test the anharmonic part of the well; repeating the measurement at lower frequencies or higher power could expose deviation from Eq. (10) and thus map the level structure directly.
  • The arbitrary power axis means the extracted coupling coefficient k is not an absolute microwave amplitude; an independent calibration of the on-chip field would turn Eq. (8) into a quantitative mm-wave power meter.
  • If the quality factor at 4.2 K is indeed quasiparticle-limited at Q ≈ 90, cooling the same junctions below 1 K should raise Q substantially; a straightforward test is to repeat the double-peak spectroscopy at 100–200 mK and check whether the resonant-peak width narrows as expected.
  • The technique could be transferred to other high-gap materials such as niobium nitride, where even larger plasma frequencies would push the qubit operating frequency further and potentially enable operation at liquid-nitrogen temperatures, though this goes well beyond what the data demonstrate.
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.

Referee Report

4 major / 5 minor

Summary. The paper reports switching-current measurements of Nb/Al-AlOx/Nb Josephson junctions under millimeter-wave irradiation at 100–110 GHz. At these frequencies the measured switching-current histograms show double-peak structures, which the authors interpret as resonant escape of the Josephson phase from a metastable well. A spectroscopic scan of the resonant-peak position versus frequency is used to extract a plasma frequency fp=138.3(7) GHz, and power-dependent switching-current branches at three frequencies are fitted with the strong-driving model of Eqs. (8)–(10), using a common plasma frequency ωp/2π=145.5 GHz and quality factors Q≈90. The paper concludes that junctions with plasma frequencies around 140 GHz have been fabricated and that phase qubits operating near 100 GHz are feasible.

Significance. The qualitative observation of clear double-peak switching-current structures at millimeter-wave frequencies is a credible and potentially important experimental result, and it extends earlier microwave-frequency resonant-escape studies to a regime relevant for higher-temperature phase qubits. The use of the published strong-driving model (Fistul et al., ref. [6]) is appropriate, and the data set is systematic: three irradiation frequencies, power sweeps, and a frequency scan of the resonant peak. If the quantitative parameter inconsistencies are resolved, the work would provide a useful rapid characterization method for high-plasma-frequency junctions. The paper does not provide code or machine-checked proofs, but the experimental protocol is clearly described. The main weakness is that the central quantitative claim rests on fits with several free parameters and on a common plasma frequency that is not reconciled with the independently extracted spectroscopic value.

major comments (4)
  1. [Sec. IV, Fig. 2 and Fig. 3] The power-branch fits in Fig. 4 yield a common plasma frequency ωp/2π=145.5 GHz, whereas the independent spectroscopic fit in Fig. 3 gives ωp/2π=138.3(7) GHz. This ≈7 GHz (5%) discrepancy is not discussed, yet the abstract and summary claim that the strong-driving model 'well explains' the data and that 'these measurements yield similar results.' Because the common ωp is a central output of the fits and the spectroscopic value is the physical plasma frequency, the authors should either reconcile the difference through a quantitative uncertainty analysis (e.g., fit covariance, systematic shifts from power-dependent peak positions, or model approximations) or soften the quantitative claim. As written, the specific assertion of agreement with Eqs. (8)–(10) is not fully supported.
  2. [Sec. IV, Fig. 2 and Fig. 3] The undriven escape-rate fit gives Ic=242(1) μA (Fig. 2), and this value is used as a fixed parameter for the Fig. 4 fits, while the spectroscopic fit of Eq. (3) gives Ic=258(3) μA (Fig. 3). The 16 μA difference between the two determinations is not mentioned. Since the small-oscillation frequency ω0(I) in Eq. (3) and the normalized shift ⟨δIsw⟩ depend on Ic, using Ic=242 μA may bias the extracted ωp and Q. The authors should justify the choice of one Ic value, include the uncertainty in Ic in the fit propagation, or fit Ic as a constrained parameter.
  3. [Sec. IV, Fig. 4] The Fig. 4 fits use a common ωp but allow Q and k^{-1}C^{-2} to be independent for each frequency, and the mm-wave power axis is in arbitrary units. With these free parameters, the visual agreement shown in Fig. 4 does not by itself establish that the common value ωp/2π=145.5 GHz is uniquely determined; model error or the Ic inconsistency could be absorbed by the per-frequency parameters. The reported values k^{-1}C^{-2}=0.82, 0.115, and 0.6 (arbitrary units) differ by almost an order of magnitude across the three frequencies, which is unexplained if k and C are junction properties. Please report residuals, parameter uncertainties, and a sensitivity analysis (e.g., fixing Q to a common value or using the spectroscopic ωp) to demonstrate the robustness of the common-ωp claim.
  4. [Sec. IV, Eqs. (8) and (10)] The text states that the red curves in Fig. 4 are fits to 'solutions of Eq. (8)', but the reported fit parameter k^{-1}C^{-2} appears explicitly only in the harmonic-approximation form Eq. (10), not in the multilevel Eq. (8). This ambiguity is important for reproducibility: if Eq. (10) was actually solved, the text and figure caption should say so; if the full Eq. (8) was used, the definition of C^{-2}, the summation range over n,m, and the truncation of the level sum should be provided.
minor comments (5)
  1. [Sec. IV, Fig. 4] The red fits in Fig. 4 are shown without residuals, error bars, or a goodness-of-fit measure, so the reader cannot independently assess how well the model reproduces the data.
  2. [Sec. IV, Fig. 3] The double-Lorentzian fitting procedure used to extract the peak positions is not described; please state the fit function, any shared parameters, and the uncertainties on the peak positions, since those uncertainties propagate into the spectroscopic plasma frequency.
  3. [Sec. II, Eq. (8)] Please define all symbols in Eq. (8), including the summation range over n,m and the normalization of the matrix elements f_nm; currently k^{-1}C^{-2} appears only in Eq. (10).
  4. [Sec. III, Fig. 1] The schematic and text omit the current-divider ratio and the cut-off frequencies of the low-pass filters; stating these values would improve reproducibility of the noise environment.
  5. [Throughout] There are minor typographical issues, including 'T op' in the Fig. 2 caption and 'V A100' in Sec. III A, that should be corrected.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the strong-driving model is an externally published theory, and the Fig. 4 comparisons are fits rather than parameter-free predictions, which is a fitting-strength issue, not a circular definition.

full rationale

The derivation chain is largely non-circular. The central model, Eq. (8), is taken from ref. [6] (Fistul, Wallraff, and Ustinov), but that is a published physical theory with independent mathematical content that has been applied and tested elsewhere, including in later work by other groups (e.g., ref. [16]). The present paper does not assume the conclusion it claims to observe: the double-peak switching-current distributions are raw experimental data, and the statement that they 'result from resonant escape from a stationary state [6]' is an interpretation of those data using an externally stated model, not a definitional equivalence. The power-dependence fits in Fig. 4 use Eq. (8)/(10) with k, Q, and a common plasma frequency omega_p as fit parameters, with the power axis in arbitrary units; the resulting agreement is therefore a fit rather than an independent parameter-free prediction. That reduces the evidential force of the quantitative claim, but it is a limitation in inference, not circularity. The internal discrepancies between the spectroscopic plasma frequency from Fig. 3 (omega_p/2pi = 138.3 GHz) and the fitted common value from Fig. 4 (omega_p/2pi = 145.5 GHz), and the different critical currents used (242 microA vs 258 microA), are important consistency and model-validity concerns, but they do not mean the model is defined in terms of the data it purports to explain. The only empirically load-bearing self-citation, Eq. (8) from ref. [6], is a falsifiable theory that was published before and independently of the present data; the fabrication citation [18] merely documents the sample process and does not bear on the physical reasoning. No step reduces by construction to its inputs, and no fitted parameter is renamed as a prediction.

Assumptions & free parameters 8 free parameters · 5 assumptions · 0 invented entities

The central claim rests on established RCSJ and Kramers theories, plus the strong-driving escape model of Fistul et al. (ref. 6), which is the authors' own framework. The quantitative fits require eight fitted or calibration parameters, including an arbitrary power scale, so the model agreement is partly a fitting exercise rather than a parameter-free prediction.

free parameters (8)
  • Critical current (no irradiation fit) = Ic = 242(1) uA
    Fitted to the switching current histogram without irradiation via Eq. (7) (Fig. 2).
  • Escape temperature = Tesc = 6.3(4) K
    Fitted parameter in Eq. (7) fit to the no-irradiation histogram (Fig. 2); higher than bath temperature 4.2 K, not discussed.
  • Plasma frequency (no irradiation fit) = omega_p/(2 pi) = 0.1(2) THz
    Fit parameter in Eq. (7) fit to the no-irradiation histogram; the paper notes it is not precise.
  • Plasma frequency (spectroscopic) = omega_p/(2 pi) = 138.3(7) GHz
    Fit to resonant peak position vs frequency using f0(I) = fp (1 - (I/Ic)^2)^(1/4) (Fig. 3).
  • Critical current (spectroscopic) = Ic = 258(3) uA
    Fit parameter in the same spectroscopic fit (Fig. 3); differs from the 242 uA used in the power fits.
  • Plasma frequency (power fits) = omega_p/(2 pi) = 145.5 GHz
    Common fit parameter for the three power-dependent fits in Fig. 4; differs from the spectroscopic value of 138.3 GHz.
  • Quality factor Q = Q about 90-97 (frequency dependent)
    Fitted per frequency in Fig. 4; reported as about 97 at 109.8 GHz, about 91 at 108.0 GHz, about 90 at 106.2 GHz, with no error bars.
  • Amplitude parameter k^-1 C^-2 = 0.82, 0.115, 0.6 x 10^50 (arb. unit)
    Fitted per frequency in Fig. 4; absorbs the unknown mm-wave coupling k and capacitance C. The power axis is in arbitrary units.
assumptions (5)
  • domain assumption RCSJ model describes the junction dynamics (Eq. 1).
    Standard model of current-biased Josephson junctions; used throughout the theory section.
  • standard math Kramers thermal escape rate formula (Eq. 4).
    Standard result for thermally activated escape; used to fit the no-irradiation distribution.
  • domain assumption Strong-driving model of Fistul et al. (Eq. 8) applies.
    The paper assumes the strong-driving limit (omega/omega_p)^5 much greater than hbar omega_p / E_J (computed about 170) and uses the effective barrier suppression model from the authors' own prior work (ref. 6).
  • domain assumption Harmonic approximation reduces the sum in Eq. (8) to the n=0, m=1 term (Eq. 10).
    Used to compute the power dependence fits; ignores multi-level contributions beyond the lowest transition, which is not explicitly justified for the data.
  • standard math The escape temperature Tesc is approximately constant and given by Eq. (6).
    Used to derive Eq. (7) for the escape rate fit.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Resonant escape in Josephson tunnel junctions under millimeter-wave irradiation." pith.science (2026). https://pith.science/paper/TSWTZ74R

@misc{pith2026241115048,
  author       = {Pith},
  title        = {Pith review of: Resonant escape in Josephson tunnel junctions under millimeter-wave irradiation},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/TSWTZ74R}},
  note         = {Machine review of arXiv:2411.15048}
}
read the original abstract

The microwave-driven dynamics of the superconducting phase difference across a Josephson junction is now widely employed in superconducting qubits and quantum circuits. With the typical energy level separation frequency of several GHz, cooling these quantum devices to the ground state requires temperatures below 100 mK. Pushing the operation frequency of superconducting qubits up may allow for operation of superconducting qubits at 1 K and even higher temperatures. Here we present measurements of the switching currents of niobium/aluminum-aluminum oxide/niobium Josephson junctions in the presence of millimeter-wave radiation at frequencies above 100 GHz. The observed switching current distributions display clear double-peak structures, which result from the resonant escape of the Josephson phase from a stationary state. We show that the data can be well explained by the strong-driving model including the irradiation-induced suppression of the potential barrier. While still being measured in the quasi-classical regime, our results point towards a feasibility of operating phase qubits around 100 GHz.

Figures

Figures reproduced from arXiv: 2411.15048 by the authors.

Figure 1
Figure 1. A current ramp with constant ramp rate ˙I is generated by a sawtooth generator attached to a voltage￾controlled current source at room temperature. A pulse generator connected to the sawtooth generator restarts the ramp in an adjustable time interval ∆tpuls and the zero-crossing of the current ramp triggers a start signal. The voltage across the junction is amplified and a stop signal is triggered as soon as a thres… view at source ↗
Figure 3
Figure 3. FIG. 3. Results of the spectroscopic determination of the [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 2
Figure 2. FIG. 2 [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figures from the paper (1 more)
Figure 4
Figure 4. Figure 4: FIG. 4. Measurement of the switching current distribution in dependence of the applied mm-wave power for three different [PITH_FULL_IMAGE:figures/full_fig_p005_4.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

19 extracted references · 17 canonical work pages

  1. [6]

    M. V. Fistul, A. Wallraff, and A. V. Ustinov, Quantum escape of the phase in a strongly driven Josephson junc- tion, Physical Review B 68, 060504 (2003)

    work page 2003

  2. [1]

    Kjaergaard, M

    M. Kjaergaard, M. E. Schwartz, J. Braum¨ uller, P. Krantz, J. I.-J. Wang, S. Gustavsson, and W. D. Oliver, Superconducting Qubits: Current State of Play, Annual Review of Condensed Matter Physics 11, 369 (2020)

    2020

  3. [2]

    B. W. Roberts, Survey of superconductive materials and critical evaluation of selected properties, Journal of Phys- ical and Chemical Reference Data 5, 581 (1976)

    work page 1976

  4. [3]

    J. M. Martinis, M. H. Devoret, and J. Clarke, Energy- Level Quantization in the Zero-Voltage State of a Current-Biased Josephson Junction, Physical Review Letters 55, 1543 (1985)

    work page 1985

  5. [4]

    J. M. Martinis, M. H. Devoret, and J. Clarke, Experi- mental tests for the quantum behavior of a macroscopic degree of freedom: The phase difference across a Joseph- son junction, Physical Review B 35, 4682 (1987)

    work page 1987

  6. [5]

    Wallraff, T

    A. Wallraff, T. Duty, A. Lukashenko, and A. V. Usti- nov, Multiphoton Transitions between Energy Levels in a Current-Biased Josephson Tunnel Junction, Physical Review Letters 90, 037003 (2003)

    work page 2003

  7. [7]

    Tinkham, Introduction to superconductivity, 2nd ed., Dover books on physics (Dover Publ., Mineola, NY, 2004)

    M. Tinkham, Introduction to superconductivity, 2nd ed., Dover books on physics (Dover Publ., Mineola, NY, 2004)

    work page 2004

  8. [8]

    Ambegaokar and B

    V. Ambegaokar and B. I. Halperin, Voltage Due to Ther- mal Noise in the dc Josephson Effect, Physical Review Letters 22, 1364 (1969)

    work page 1969

Show all 19 references
  1. [9]

    T. A. Fulton and L. N. Dunkleberger, Lifetime of the zero-voltage state in Josephson tunnel junctions, Physical Review B 9, 4760 (1974)

    1974

  2. [10]

    H. A. Kramers, Brownian motion in a field of force and the diffusion model of chemical reactions, Physica 7, 284 (1940)

    1940

  3. [11]

    Weiss, Quantum dissipative systems, 2nd ed., Series in modern condensed matter physics No

    U. Weiss, Quantum dissipative systems, 2nd ed., Series in modern condensed matter physics No. 10 (World Scien- tific, Singapore ; River Edge, NJ, 1999)

    1999

  4. [12]

    Grabert, P

    H. Grabert, P. Olschowski, and U. Weiss, Quantum de- cay rates for dissipative systems at finite temperatures, Physical Review B 36, 1931 (1987)

    1987

  5. [13]

    B¨ uttiker, E

    M. B¨ uttiker, E. P. Harris, and R. Landauer, Thermal ac- tivation in extremely underdamped Josephson-junction circuits, Physical Review B 28, 1268 (1983)

    1983

  6. [14]

    Grabert and U

    H. Grabert and U. Weiss, Crossover from Thermal Hop- ping to Quantum Tunneling, Physical Review Letters53, 1787 (1984)

    1984

  7. [15]

    Grønbech-Jensen, M

    N. Grønbech-Jensen, M. G. Castellano, F. Chiarello, M. Cirillo, C. Cosmelli, L. V. Filippenko, R. Russo, and G. Torrioli, Microwave-Induced Thermal Escape in Josephson Junctions, Physical Review Letters93, 107002 (2004)

    2004

  8. [16]

    H. F. Yu, X. B. Zhu, J. K. Ren, Z. H. Peng, D. J. Cui, H. Deng, W. H. Cao, Y. Tian, G. H. Chen, D. N. Zheng, X. N. Jing, L. Lu, and S. P. Zhao, Resonant phase escape in Bi 2Sr2CaCu2O8+δ surface intrinsic Josephson junc- tions, New Journal of Physics 15, 095006 (2013)

    2013

  9. [17]

    J. M. Martinis, S. Nam, J. Aumentado, K. M. Lang, and C. Urbina, Decoherence of a superconducting qubit due to bias noise, Physical Review B 67, 094510 (2003)

    2003

  10. [18]

    Strobel, L

    U. Strobel, L. Radtke, L. Kamps, J. Voss, J. Lisenfeld, H. Rotzinger, and A. V. Ustinov, (in preparation)

  11. [19]

    Wallraff, Fluxon dynamics in annular Josephson junctions: from relativistic strings to quantum par- ticles, Mesoscopic Structures and Superconductivity No

    A. Wallraff, Fluxon dynamics in annular Josephson junctions: from relativistic strings to quantum par- ticles, Mesoscopic Structures and Superconductivity No. 1 (Lehrstuhl f¨ ur Mikrocharakterisierung, Friedrich- Alexander-Univ., Erlangen, 2001)

    2001

Pith tools

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