Pith. sign in

REVIEW 3 major objections 5 minor 24 references

Probing the minigap in topological insulator-based Josephson junctions under radio frequency irradiation

T0 review · 3 major / 5 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 The contact-resistance minigap probe in topological-insulator Josephson junctions survives rf irradiation.

desk verdict A useful, honest extension of the contact-resistance method to rf environments, with one hand-set coupling ratio that should be calibrated before the quantitative match is taken at face value. read the letter →

arxiv 1908.05010 v1 pith:PEAKDCY7 submitted 2019-08-14 cond-mat.mes-hall

Guang Yang , Zhaozheng Lyu , Xiang Zhang , Fanming Qu , Li Lu This is my paper · ORCID
classification cond-mat.mes-hall PACS 74.45.+c03.65.Vf85.25.Cp
keywords topologicalinsulatorJosephsonjunctionradiofrequencyirradiationcontact-resistancemeasurementminigapAndreevboundstatesShapirostepsBTKtheory
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

The paper extends a contact-resistance technique for measuring the minigap, and hence the Andreev bound states, in superconductor/topological-insulator/superconductor (S-TI-S) Josephson junctions into the regime of radio-frequency (rf) irradiation. It reports that as rf power increases, the measured minigap feature broadens and shifts to higher bias currents, mirroring the outer border of the Shapiro-step region in the same junction. The authors show that the rf-dependent contact resistance can be reproduced quantitatively by combining the resistively shunted Josephson junction (RSJ) model for the phase dynamics with the Blonder-Tinkham-Klapwijk (BTK) theory for the normal-metal contact. If correct, this establishes that the contact-resistance probe remains trustworthy under rf drive, a step toward using rf techniques for fast manipulation and readout in Majorana-related experiments.

What carries the argument

The load-bearing objects are the RSJ equation $I_J + I_{rf0}\cos(2\pi f_{rf} t) = \frac{\hbar}{2eR}\frac{d\varphi}{dt} + I_c\sin\varphi$ and the BTK scattering theory of a normal-metal/superconductor contact. The RSJ equation supplies the time-dependent phase $\varphi(t)$, which feeds the proximity minigap through the $4\pi$-periodic form $\Delta(t) = \Delta_0 \cos(\varphi(t)/2)$. The BTK theory supplies a precomputed resistance map giving the instantaneous differential contact resistance as a function of the instantaneous minigap and bias current; time-averaging over one rf cycle gives the predicted $dV_b/dI_b$. The same machinery reproduces the Shapiro map of the junction, tying the broadening of the minigap feature to the same mode-locking envelope traced on the Josephson side.

What would settle it

Measure the actual rf current amplitude across the Pd-Bi2Te3 interface directly, for example by calibrated on-chip rf detection or transport, over a range of frequencies and powers, and compare it with the $I'_{rf0}$ inferred from the RSJ-plus-BTK fit; if the inferred coupling ratio varies with frequency, power, or device geometry, or if the simulation fails at a power where the ratio is fixed by measurement, the claim that RSJ plus BTK fully accounts for the rf contact resistance is falsified.

Watch

Extended reading notes

Core claim

The central claim is that the differential contact resistance $dV_b/dI_b$ across a normal-metal Pd electrode touching the junction area of an Al-Bi2Te3-Al Josephson junction can still be read as a direct probe of the proximity-induced minigap when the junction is under rf irradiation. The rf field acts in two ways: it drives the superconducting phase difference $\varphi(t)$ through the RSJ circuit equation, making the surface-state minigap $\Delta(t) = \Delta_0 \cos(\varphi(t)/2)$ time-dependent, and it adds an alternating current across the Pd-Bi2Te3 interface. Using BTK-theory parameters fixed at zero rf power (minigap $\Delta_0 = 20.4\,\mu$eV, $N = 25.5$ channels, barrier strength $Z = 0.66$, temperature $T = 0.1$ K), the authors compute the instantaneous $dV_b/dI_b$ from the BTK resistance map, time-average it, and reproduce the measured rf-broadened dip-peak-dip structure. The claim is that this two-effect picture fully accounts for the data.

Load-bearing premise

The simulation assumes that the rf current reaching the Pd-Bi2Te3 interface is exactly $I'_{rf0} = I_{rf0}/14.4$, a factor chosen to make the model match the data with no independent calibration of the coupling between the antenna and the contact.

Editorial extensions

If this is right

  • The contact-resistance minigap signature under rf broadens along the same envelope as the Shapiro-step region, indicating that both reflect the same rf-driven phase dynamics.
  • Without accounting for rf power, the apparent minigap would look larger and extend to higher energies; the broadening is an instrumental effect captured by RSJ plus BTK, not a change in the intrinsic gap.
  • The method can be used under rf to track the minigap, including the expected $4\pi$-periodic linear closing signature, provided the rf coupling is characterized.
  • In these devices the supercurrent is dominated by a $2\pi$-periodic component, and the absence of missing odd Shapiro steps is consistent with bulk-dominated transport; probing surface-state $4\pi$ signatures requires devices with a larger surface contribution.
  • The RSJ-plus-BTK simulation workflow can be applied to future S-TI-S devices under rf to extract the minigap from measured contact resistance.

Reading between the lines

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

  • An independent calibration of the rf coupling ratio, for example by directly measuring the rf current across the Pd-Bi2Te3 interface, would turn the current simulation into a parameter-free test; until then the quantitative agreement is consistency rather than confirmation.
  • In devices with a non-negligible $4\pi$-periodic surface supercurrent, the same two-effect picture predicts a qualitatively different rf dependence of the contact resistance, such as additional dips or half-integer features, which would be a testable distinction.
  • The rf-induced broadening of the minigap feature implies that rf readout schemes for Majorana qubits must either operate at low rf power or deconvolve the RSJ/BTK envelope, otherwise an apparent gap closing could be misinterpreted.
  • The simulation framework already produces instantaneous quantities, so a natural extension is pulsed-rf operation to time-resolve the gap dynamics rather than the time-averaged response reported here.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

3 major / 5 minor

Summary. The paper reports measurements of the differential contact resistance across a Pd-Bi2Te3 interface in an Al-Bi2Te3-Al Josephson junction under radio-frequency (rf) irradiation. The authors observe that the minigap-related dip-peak-dip structure broadens and extends to higher bias currents with increasing rf power, in a way they describe as similar to the rf-power dependence of the outer border of the Shapiro-step region. They simulate both the Shapiro map (using the RSJ model) and the rf-dependent contact resistance (using the RSJ model plus BTK theory) and conclude that the contact-resistance-measurement method remains valid for detecting the minigap under rf irradiation. The zero-rf BTK fit yields the minigap, channel number, barrier strength, and effective temperature; the rf simulation additionally introduces a ratio I'_rf0 = I_rf0/14.4 for the rf current amplitude across the Pd-Bi2Te3 interface.

Significance. If the rf generalization is quantitatively reliable, the work is useful for planned experiments that combine rf manipulation or readout with minigap detection in topological-insulator Josephson junctions, particularly for Majorana-related physics. The paper has real strengths: the zero-rf BTK fit in Fig. 3b is clean, the Shapiro-step simulation in Fig. 2c uses directly measured values of I_c and R, and the qualitative behavior of the rf-broadened contact resistance is physically plausible. However, the central quantitative claim rests on an uncalibrated coupling ratio introduced for the rf current at the Pd-Bi2Te3 interface, and the agreement between simulation and data is assessed visually rather than by a quantitative metric. These issues materially weaken the claim that the method is validated under rf.

major comments (3)
  1. [Supplemental Material, Sec. 3] The choice I'_rf0 = I_rf0/14.4 is introduced without derivation, independent calibration, or uncertainty estimate. Since I'_rf0 controls both the amplitude scale and the power-dependent broadening of the simulated dV_b/dI_b in Fig. 3(c,d), the agreement between the black simulations and red data in Fig. 3(d) is not a parameter-free prediction; it is a fit of one parameter to the same data it is used to validate. The manuscript should either calibrate this ratio independently (for example, by a direct rf-current measurement or a separate control experiment) or report a sensitivity analysis over I'_rf0 and state explicitly that the rf validation is conditional on this parameter.
  2. [Main text, Sec. 3 and Fig. 3(c)] The statement that the broadened minigap border follows a trace similar to the outer border of the Shapiro-step region is not an independent check, because the red dashed border I_e + I'_rf0 is constructed with I'_rf0 proportional to I_rf0, the same proportionality that sets the Shapiro border I_c + I_rf0 in Fig. 2(d). The similar rf-power dependence is therefore built in by construction rather than emerging as a separate prediction.
  3. [Main text, Sec. 3 (Figs. 2 and 3)] The agreement between experiment and simulation is assessed only visually; no goodness-of-fit measure, noise level, or uncertainty on the extracted parameters is reported. For the central claim that the contact-resistance method can still be well interpreted numerically under rf, the authors should provide a quantitative comparison for representative line cuts, including residuals or error bars, and should state how the uncertainty of the zero-rf BTK parameters propagates into the rf simulation.
minor comments (5)
  1. [Supplemental Material, Sec. 3] There are several typographical errors: 'Blonder-Tinkham-Klapeijk' should be 'Blonder-Tinkham-Klapwijk', 'trsnsport' should be 'transport', and 'dimentional' should be 'dimensional'.
  2. [Main text, Sec. 3, Eq. (1)] Equation (1) is typeset with garbled notation, including 'd () 1sin(2π )' and 'd4 π r'; please rewrite the equation in standard form and define every symbol at first use.
  3. [Main text, Sec. 3 (first paragraph)] The symbol JI is used both as a general dc current and as the Josephson current; this notation should be clarified in the text and figure captions to avoid confusion.
  4. [Main text, Sec. 3 and Fig. 3(c)] The characteristic current I_e is not defined in the text; please provide a formal definition when it first appears, rather than only describing it in the figure caption.
  5. [Supplemental Material, Sec. 2] The conversion of rf power to I_rf0 uses a 50-ohm reference impedance without discussing cable attenuation, impedance mismatch, or antenna coupling efficiency; please state these assumptions or quantify their effect.

Circularity Check

3 steps flagged · score 5.0 of 10

The rf simulation's agreement is partly fitted: I'_rf0=I_rf0/14.4 is chosen by hand, and the 'similar' minigap/Shapiro border is built into the proportional construction.

  1. fitted input called prediction [Supplemental Material, Sec. 3; Figs. 3c and 3d]
    "Here, we choose I'_rf0 = I_rf0/14.4 to simulate the experimental data."

    This is the load-bearing parameter of the rf contact-resistance simulation. I'_rf0 controls both the overall scale and the power-dependent broadening of the simulated dV_b/dI_b curves, and it is introduced as a free choice ('we choose') with no independent calibration, error estimate, or sensitivity analysis. The subsequent agreement between the simulated black curves and the measured red curves in Fig. 3d is therefore partly enforced by this fitted value rather than being a parameter-free prediction. The zero-rf BTK parameters are fixed from the no-rf fit, so the circularity is partial, but the claim that RSJ+BTK 'can still reproduce the measured data' under rf is weakened by construction.

  2. self definitional [Main text Sec. 3; Fig. 2d and Fig. 3c captions]
    "The red dashed line in Fig. 3c represents the border of I_b ≤ I_e + I'_rf0, where I_e is a characteristic current beyond which there is no signature of superconductivity without rf irradiation. ... The black line represents the critical supercurrent I_c, and the red curve represents the sum of the critical supercurrent I_c and the amplitude of rf-driving current I_rf0."

    Because I'_rf0 is set proportional to I_rf0 (with the 1/14.4 ratio adjusted in the Supplemental Material), the plotted minigap border I_e + I'_rf0 has exactly the same rf-power scaling as the Shapiro-step border I_c + I_rf0. The paper's observation that the minigap broadening follows 'a similar trace' as the Shapiro-step region is therefore not an independent consequence of the RSJ+BTK dynamics; it is built into the proportional construction of the two guide lines. The comparison only tests the proportionality constant, which is the very parameter chosen to make the simulation match.

1 more flagged steps
  1. ansatz smuggled in via citation [Main text Sec. 3; Supplemental Sec. 3; refs. [10], [2-4]]
    "In the presence of rf irradiation, the time-dependent minigap of the surface state should follow the 4π-period form[10] and can be expressed as: Δ(t)=Δ0 cos(φ(t)/2)."

    The rf simulation assumes the very quantity the contact-resistance method is supposed to detect—the minigap—obeys the functional form Δ(t)=Δ0 cos(φ(t)/2), cited to prior work from the same group (ref. [10]; Supplemental refs. [2-4]). The simulated rf-induced broadening of the minigap then follows directly from this assumed form, and the match with measured dV_b/dI_b is offered as confirmation that the method still works under rf. If the cited form itself was adopted as an ansatz rather than independently established in this paper, the simulation does not independently validate the minigap's rf dependence; it reconstructs the data from the assumed model. This is a load-bearing self-citation rather than an external check.

full rationale

The zero-field BTK fit (Δ0, N, Z, T) and the RSJ Shapiro-map simulation use independently measured I_c and R and are not circular. The circularity is confined to the rf contact-resistance branch: the simulation's quantitative agreement depends on the hand-set, uncalibrated ratio I'_rf0=I_rf0/14.4, and the 'similar' power dependence of the minigap and Shapiro borders is a consequence of assuming proportionality between I'_rf0 and I_rf0. The minigap's time-dependent form is also imported from same-group prior work. Because these choices are transparent and the core zero-rf modeling is independent, the paper is not wholly circular, but the rf validation claim is substantially weaker than a parameter-free prediction. Score 5 reflects partial circularity of the central rf claim.

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

The central claim rests on the zero-rf BTK fit for the intrinsic gap parameters, RSJ dynamics for the phase, a hand-set interface rf coupling ratio, and the assumption of adiabatic instant response. No new physical entities are introduced. The parameters Delta0, N, Z, and T are fitted to the paper's own zero-rf data, while the ratio 1/14.4 is fitted to the rf data.

free parameters (5)
  • Minigap Delta0 = 20.4 micro-eV
    Fitted to zero-rf contact-resistance curve via BTK theory in Sec. 3.
  • Effective number of channels N = 25.5
    Fitted with BTK theory to zero-rf data in Sec. 3.
  • Barrier strength Z = 0.66
    Fitted with BTK theory to zero-rf data in Sec. 3.
  • Effective electron temperature T = 0.1 K
    Fitted with BTK theory to zero-rf data in Sec. 3.
  • Pd-Bi2Te3 interface rf-current ratio I'_rf0 / I_rf0 = 1/14.4
    Chosen by hand in Supplemental Sec. 3 to make the simulated contact resistance match the experimental rf data.
assumptions (6)
  • domain assumption The RSJ model with a sinusoidal 2pi-periodic current-phase relation governs junction phase dynamics under rf irradiation.
    Used to solve phi(t) via Eq. (1) in Sec. 3 and to simulate the Shapiro map in Fig. 2c.
  • domain assumption BTK theory describes the differential conductance across the Pd-Bi2Te3 interface with proximity-induced superconductivity.
    The zero-rf line shape in Fig. 3b is fitted with BTK, and the R map is generated from the BTK parameters.
  • domain assumption The instantaneous minigap follows Delta(t) = Delta0 cos(phi(t)/2), inherited from prior theory for 4pi-periodic Andreev bound states.
    Used in Sec. 3 for the contact-resistance simulation; assumes the Andreev spectrum adiabatically follows the phase.
  • domain assumption The measured differential contact resistance equals the time average of the instantaneous BTK differential resistance over an rf cycle.
    The simulation time-averages dV_b/dI_b(t) to compare with the dc measurement, as described in Supplemental Sec. 3.
  • ad hoc to paper The rf-induced current across the Pd-Bi2Te3 interface is sinusoidal with amplitude I'_rf0 = I_rf0 / 14.4.
    This ratio is chosen to match the data; no independent calibration is provided, so it is both an assumption and a free parameter.
  • domain assumption The rf power conversion formula I_rf0 = sqrt(10^(P_dBm/10) * 1000 / 50) assumes a 50 ohm termination and no significant attenuation.
    Used in Supplemental Sec. 2 to relate rf power to rf current; any uncorrected attenuation would change the absolute scale.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Probing the minigap in topological insulator-based Josephson junctions under radio frequency irradiation." pith.science (2026). https://pith.science/paper/PEAKDCY7

@misc{pith2026190805010,
  author       = {Pith},
  title        = {Pith review of: Probing the minigap in topological insulator-based Josephson junctions under radio frequency irradiation},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/PEAKDCY7}},
  note         = {Machine review of arXiv:1908.05010}
}
read the original abstract

Recently, a contact-resistance-measurement method was developed to detect the minigap, hence the Andreev bound states (ABSs), in Josephson junctions constructed on the surface of three-dimensional topological insulators (3D TIs). In this work, we further generalize that method to the circumstance with radio frequency (rf) irradiation. We find that with the increase of rf power, the measured minigap becomes broadened and extends to higher energies, in a way similar to the rf power dependence of the outer border of the Shapiro step region. We show that the corresponding data of contact resistance under rf irradiation can be well interpreted by using the resistively shunted Josephson junction model (RSJ model) and the Blonder-Tinkham-Klapwijk (BTK) theory. Our findings could be useful when using the contact-resistance-measurement method to study the Majorana-related physics in topological insulator-based Josephson junctions under rf irradiation.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

24 extracted references · 24 canonical work pages

  1. [1]

    Fu L and Kane C L 2008 Phys Rev Lett 100 096407

    work page 2008

  2. [2]

    Sochnikov I, Maier L, Watson C A, Kirtley J R, Gould C, Tkachov G, Hankiewicz E M, Brüne C, Buhmann H, Molenkamp L W and Moler K A 2015 Phys Rev Lett 114 066801

    work page 2015

  3. [3]

    Kurter C, Finck A D K, Hor Y S and Van Harlingen D J 2015 Nature communications 6 7130

    work page 2015

  4. [4]

    Kayyalha M, Kazakov A, Miotkowski I, Khlebnikov S, Rokhinson L P and Chen Y P

    work page

  5. [5]

    Rokhinson L P, Liu X and Furdyna J K 2012 Nat Phys 8 795

    work page 2012

  6. [6]

    Wiedenmann J, Bocquillon E, Deacon R S, Hartinger S, Herrmann O, Klapwijk T M, Maier L, Ames C, Brüne C, Gould C, Oiwa A, Ishibashi K, Tarucha S, Buhmann H and Molenkamp L W 2016 Nature communications 7 10303

    work page 2016

  7. [7]

    Bocquillon E, Deacon R S, Wiedenmann J, Leubner P, Klapwijk T M, Brune C, Ishibashi K, Buhmann H and Molenkamp L W 2017 Nature nanotechnology 12 137

    work page 2017

  8. [8]

    Le Calvez K, Veyrat L, Gay F, Plaindoux P, Winkelmann C B, Courtois H and Sacépé B 2019 Communications Physics 2

    work page 2019

Show all 24 references
  1. [9]

    Pang Y , Wang J, Lyu Z, Yang G, Fan J, Liu G, Ji Z, Jing X, Yang C and Lu L 2016 Chinese Physics B 25

    2016

  2. [10]

    Lyu Z, Pang Y , Wang J, Yang G, Fan J, Liu G, Ji Z, Jing X, Yang C, Qu F and Lu L 2018 Physical Review B 98

    2018

  3. [11]

    Yang G, Lyu Z, Wang J, Ying J, Zhang X, Shen J, Liu G, Fan J, Ji Z, Jing X, Qu F and Lu L

  4. [12]

    Barone A and Paterno G 1982 Physics and Application of the Josephson Effect

    1982

  5. [13]

    Blonder G E, Tinkham M and Klapwijk T M 1982 Physical Review B 25 4515

    1982

  6. [14]

    Probing the minigap in topological insulator-based Josephson junctions under radio frequency irradiation

    Domínguez F, Kashuba O, Bocquillon E, Wiedenmann J, Deacon R S, Klapwijk T M, Platero G, Molenkamp L W, Trauzettel B and Hankiewicz E M 2017 Physical Review B 95 195430 Supplemental Material for “Probing the minigap in topological insulator-based Josephson junctions under radi...

    2017

  7. [15]

    The resistively shunted junction model

  8. [16]

    Simulating the Shapiro map

  9. [17]

    Simulation of the differential contact resistance

  10. [18]

    The resistively shunted junction (RSJ) model A Josephson junction can be treated as an equivalent circuit with a pure Josephson junction J in parallel with a shunted resistor R , as shown below.[1] The current continuity equation for the circuit is: J Jc () sin( ( ))VtI It R φ...

  11. [19]

    Simulating the Shapiro map When the Josephson junction is exposed to r adio frequency (rf) irradiation, an ac current rfI will be applied to the Josepshon junciton. The equivalent circuit now becomes: The current continuity equation now becomes: J rf0 rf c d ()cos(2 ) sin( ( )...

  12. [20]

    Simulation of the differential contact resistance We have proven that the Blonder -Tinkham-Klapeijk (BTK) theory can be used to describe the electron trsnsport processes across the interface between a normal metal and the surface of a three dimentional topological insulators (...

  13. [21]

    Barone A and Paterno G 1982 Physics and Application of the Josephson Effect John Wiley and Sons

    1982

  14. [22]

    Lyu Z, Pang Y , Wang J, Yang G, Fan J, Liu G, Ji Z, Jing X, Yang C, Qu F and Lu L 2018 Physical Review B 98 155403

    2018

  15. [23]

    Yang G, Lyu Z, Wang J, Ying J, Zhang X , Shen J, Liu G, Fan J, Ji Z, Jing X, Qu F and Lu L 2019 arXiv: 1904.02677 [cond-mat.mes-hall]

    2019 arXiv

  16. [24]

    Pang Y , Wang J, Lyu Z, Yang G, Fan J, Liu G, Ji Z, Jing X, Yang C and Lu L 2016 Chinese Physics B 25 117402

    2016

Pith tools

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