{"id":"17d276c8-2d60-4d27-a11d-60dda18b8a21","arxiv_id":"1908.05010","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Contact-resistance detection of the minigap in Al-Bi2Te3-Al Josephson junctions remains valid under rf irradiation, with the rf-broadened gap feature reproduced by an RSJ-plus-BTK simulation.","lead":"This paper tests whether a contact-resistance method for measuring the minigap in a topological insulator Josephson junction remains reliable under radio frequency irradiation. It finds the gap feature broadens with rf power and can be reproduced by an RSJ-plus-BTK simulation.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The rf simulation's quantitative agreement is not parameter-free: I'_rf0=I_rf0/14.4 is a hand-set coupling ratio (Supplemental Sec. 3), so the claimed validation of the contact-resistance method under rf rests on an uncalibrated fit parameter.","rationale":"I read the paper as making a modest but useful claim: the contact-resistance method, previously demonstrated at zero rf, still yields data under rf that can be interpreted with the same physical ingredients (RSJ for the Josephson phase dynamics and BTK for the normal-TI interface), provided one allows the rf current at the interface to differ from the rf current through the junction. The reader correctly identified the I'_rf0 = I_rf0/14.4 ratio in Supplemental Sec. 3 as the least secure link. My stress-test agrees and sharpens the point: this ratio is not merely a missing error bar; it is a free parameter that controls the simulated line shapes and the power dependence, so the quoted agreement is a fit rather than a parameter-free prediction. That said, the concern does not overturn the paper. The zero-rf parameters are fixed before the rf comparison, the qualitative evolution with rf power is reproduced, and the authors do not overclaim a new theory. The appropriate verdict remains CONDITIONAL as the reader stated, so no verdict change is needed. The concrete test I propose would settle the issue by removing the free parameter through an independent normal-state calibration of the rf current at the interface; if that is experimentally impractical, a sensitivity analysis varying the ratio around 14.4 and reporting residuals would be a minimum alternative.","tokens_in":6872,"tokens_out":7632,"duration_ms":89492,"concrete_test":"Directly calibrate the rf coupling ratio by measuring the rf-induced current amplitude at the Pd-Bi2Te3 interface with the junction driven into the normal state (e.g., T > Tc or B > Bc), using the same antenna, frequency, and power range as in Fig. 3. Then rerun the Fig. 3d simulation with I'_rf0 fixed to this measured value (no free parameter) and compare to the experimental curves. If the simulated curves deviate from the data by more than the experimental noise, the agreement reported in the paper is not an independent confirmation of the method; if they match, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that the rf-dependent contact resistance can be quantitatively interpreted by RSJ plus BTK using the zero-rf parameters. The load-bearing link is the rf current amplitude at the Pd-Bi2Te3 interface, set in Supplemental Material Sec. 3 to I'_rf0 = I_rf0/14.4. This ratio is introduced without derivation, independent calibration, or error estimate. Because I'_rf0 controls both the overall scaling of the simulated dV/dI versus bias current and the power-dependent broadening, the agreement in Fig. 3d is partly built into the choice of this parameter. Moreover, the similarity between the minigap border (I_e + I'_rf0) and the Shapiro-step border (I_c + I_rf0) is guaranteed by the proportionality I'_rf0 ∝ I_rf0, so that specific 'similar rf power dependence' is not an independent prediction. The zero-rf BTK fit to Fig. 3b is legitimate, and the qualitative rf broadening is physically reasonable, but the one-free-parameter nature of the rf simulation, with no reported goodness-of-fit or sensitivity analysis, substantially weakens the strength of the claim that the method is validated under rf.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":7174,"tokens_out":3262,"duration_ms":33419,"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":[{"comment":"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.","section":"Supplemental Material, Sec. 3"},{"comment":"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.","section":"Main text, Sec. 3 and Fig. 3(c)"},{"comment":"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.","section":"Main text, Sec. 3 (Figs. 2 and 3)"}],"minor_comments":[{"comment":"There are several typographical errors: 'Blonder-Tinkham-Klapeijk' should be 'Blonder-Tinkham-Klapwijk', 'trsnsport' should be 'transport', and 'dimentional' should be 'dimensional'.","section":"Supplemental Material, Sec. 3"},{"comment":"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.","section":"Main text, Sec. 3, Eq. (1)"},{"comment":"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.","section":"Main text, Sec. 3 (first paragraph)"},{"comment":"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.","section":"Main text, Sec. 3 and Fig. 3(c)"},{"comment":"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.","section":"Supplemental Material, Sec. 2"}],"recommendation":"major_revision","confidential_remarks":"The zero-rf BTK analysis and the qualitative physics are sound, but the paper's main quantitative conclusion depends on an uncalibrated ratio I'_rf0/I_rf0 = 1/14.4. If the authors cannot provide an independent calibration or a sensitivity analysis, the conclusion should be softened to a qualitative demonstration. The manuscript is otherwise within the journal's scope and the experimental data appear to be of good quality."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper does what it says: it takes an already published contact-resistance method for detecting the minigap in TI-based Josephson junctions and shows it still works under rf irradiation. The zero-rf BTK fit is clean, and the Shapiro-step simulation uses measured I_c and R rather than fitted parameters. That is real content, not a repackaging.\n\nWhat is new is the rf-irradiated dataset and the time-averaging simulation that connects the rf-driven phase dynamics (RSJ) to the BTK contact-resistance calculation. The authors are candid that their devices are bulk-dominated 2π junctions, so this is not a Majorana claim; it is a methodological extension. That is useful for experimentalists who want to use contact-resistance readout in rf environments, and the paper is clearly written for that audience.\n\nThe soft spot is the rf contact-resistance simulation. It relies on I'_rf0 = I_rf0/14.4, set in the Supplemental Material to match the data, with no independent calibration, no error bar, and no sensitivity analysis. Since I'_rf0 controls both the broadening and the location of the minigap border, the visual agreement in Fig. 3d is partly built into that choice. Also, the claimed similarity between the minigap border and the Shapiro-step border is less independent than it appears, because one is I_e + I'_rf0 and the other is I_c + I_rf0, with I'_rf0 proportional to I_rf0. That said, this is a one-free-parameter fit, not a circular validation: the zero-rf parameters are fixed before the rf comparison, and the qualitative rounding and extension of the dip-peak-dip structure is physically reasonable.\n\nI would not reject the paper over this. I would ask the authors to calibrate the coupling ratio independently, or at least show a sensitivity scan over I'_rf0/I_rf0 and report a quantitative goodness-of-fit. The self-citation pattern is appropriate here because the earlier work is theirs and directly relevant.\n\nIf I were refereeing, I would engage seriously. The paper deserves a serious referee session; it is honest, specialized, and methodologically useful.","headline":"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.","tokens_in":7720,"tokens_out":1956,"would_cite":true,"duration_ms":20663,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["74.45.+c","03.65.Vf","85.25.Cp"],"model":"deepseek-v4-flash","headline":"The contact-resistance minigap probe in topological-insulator Josephson junctions survives rf irradiation.","keywords":["topological insulator","Josephson junction","radio frequency irradiation","contact-resistance measurement","minigap","Andreev bound states","Shapiro steps","BTK theory"],"falsifier":"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.","tokens_in":6648,"feed_emoji":"⚡","tokens_out":6189,"duration_ms":56814,"temperature":0.7,"pith_summary":"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.","feed_headline":"Minigap probe survives radio-frequency irradiation","feed_subtitle":"Contact-resistance detection of Andreev bound states stays valid in rf-driven Josephson junctions.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Sets the theoretical target: proximity-induced p-wave superconductivity on a 3D topological insulator surface, $4\\pi$-periodic Andreev levels, and Majorana bound states that the minigap probe is designed to detect.","marker":"[1]"},{"why":"Earlier demonstration of the contact-resistance measurement approach for reading out Andreev bound states in S-TI-S junctions; supplies the method being generalized.","marker":"[9]"},{"why":"Establishes the zero-rf baseline: BTK fitting of the measured contact resistance gives $\\Delta_0$, $N$, $Z$, and $T$ used in all rf simulations.","marker":"[10]"},{"why":"Extends the contact-resistance method to trijunctions and to the complete closing of the minigap, motivating the need to keep the method valid under rf.","marker":"[11]"},{"why":"Supplies the resistively shunted junction model and the current-continuity equation used to simulate both the Shapiro map and the rf-modulated phase $\\varphi(t)$.","marker":"[12]"},{"why":"Supplies the BTK theory connecting the instantaneous minigap and bias current to the differential contact resistance across the Pd-Bi2Te3 interface.","marker":"[13]"},{"why":"Gives the mode-locking criterion for Shapiro steps, used to interpret the outer border of the Shapiro region as the envelope of the broadened minigap feature.","marker":"[14]"}],"fun_headline_variants":["Contact resistance reveals minigap under rf drive","Minigap detection works under rf irradiation","rf-broadened minigap fits RSJ-BTK model","Probing Andreev states in rf-driven junctions","Minigap probe extends to rf regime"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Contact resistance reveals minigap under rf drive","Minigap detection works under rf irradiation","rf-broadened minigap fits RSJ-BTK model","Probing Andreev states in rf-driven junctions","Minigap probe extends to rf regime"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000455,"raw_usage":{"total_tokens":2286,"prompt_tokens":944,"completion_tokens":1342,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":560,"completion_tokens_details":{"reasoning_tokens":1268}},"tokens_in":560,"tokens_out":1342,"duration_ms":9496,"temperature":1.0,"reasoning_tokens":1268,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:25:38.188210+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Sets the theoretical target: proximity-induced p-wave superconductivity on a 3D topological insulator surface, $4\\pi$-periodic Andreev levels, and Majorana bound states that the minigap probe is designed to detect."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Earlier demonstration of the contact-resistance measurement approach for reading out Andreev bound states in S-TI-S junctions; supplies the method being generalized."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the zero-rf baseline: BTK fitting of the measured contact resistance gives $\\Delta_0$, $N$, $Z$, and $T$ used in all rf simulations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Extends the contact-resistance method to trijunctions and to the complete closing of the minigap, motivating the need to keep the method valid under rf."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the resistively shunted junction model and the current-continuity equation used to simulate both the Shapiro map and the rf-modulated phase $\\varphi(t)$."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the BTK theory connecting the instantaneous minigap and bias current to the differential contact resistance across the Pd-Bi2Te3 interface."},{"cited_title":"Probing the minigap in topological insulator-based Josephson junctions under radio frequency irradiation","cited_arxiv_id":null,"evidence_quote":"Gives the mode-locking criterion for Shapiro steps, used to interpret the outer border of the Shapiro region as the envelope of the broadened minigap feature."}],"review_version":1}