{"id":"5ab69cc1-a87c-4ff4-af62-43b7d3b9d9c2","arxiv_id":"2507.23169","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"In a gap-engineered transmon, quasiparticles cause a resonant enhancement of frequency shift and relaxation when the qubit frequency matches the gap difference, giving a new spectroscopic probe.","lead":"This paper calculates how quasiparticles change the frequency and relaxation rate of a superconducting qubit whose two junction electrodes have different energy gaps. It finds a resonance when the qubit frequency matches the gap difference, which could be used to probe the energy distribution of quasiparticles.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Resonant enhancement is robust, but the quantitative lineshapes and inferred 'temperature' are conditional on the quasi-thermal single-distribution ansatz (Eqs. 5, 24).","rationale":"The derivation is careful and internally consistent; the results reduce to the symmetric-junction limit, and the Kubo-formula appendix strengthens the admittance calculation. The central resonant feature at ω10=δΔ is driven by the coincidence of the two BCS density-of-states singularities and should survive for any distribution with finite occupation at the lower gap edge. What is not robust is the specific temperature dependence: Eqs. (48), (58), and (60) insert the Boltzmann form Eq. (24) inside the integrals, and f0 is set by the fixed-density condition Eq. (44). A non-thermal low-energy tail—or different distributions in the two leads—would change the shape functions F and G, the collapse, and the amplitude of the predicted effects. The paper does provide integral representations for arbitrary f, so the spectroscopic idea is not invalidated, but the quantitative predictions and the meaning of the extracted effective temperature are conditional on the quasi-thermal model. This is exactly the reader's weakest assumption, and the conditional verdict remains appropriate. The abstract's claim of an 'anomalous positive frequency shift' is a separate presentation error, contradicted by the body's statement in Sec. V B that δω10^Y is negative at all parameter values; it should be corrected before publication.","tokens_in":22225,"tokens_out":19173,"duration_ms":234436,"concrete_test":"Numerically evaluate Eqs. (23), (29), and (56) with a non-thermal but physically motivated distribution, e.g., f(Δ_L+ε)=f0 (ε/Δ)^α e^{-ε/T_eff} for α ∈ {0, 0.5, 1}, renormalizing f0 at each temperature so that Eqs. (42)–(43) keep xqp fixed, and also with an independent right-lead distribution (including the f_R=0 limit). Recompute the curves of Figs. 2 and 4 in the regime |δΔ−ω10|≪δΔ. If the resonance position, weight, and universal collapse onto F(T/T*) and G(T/T*) persist for all α and for f_R=0, the quasi-thermal assumption is not load-bearing; otherwise, Eqs. (48), (58), and (60) should be presented as model-specific and the inversion proposal needs a separate demonstration.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The quantitative core—the T* scale, the universal functions F(T/T*) and G(T/T*), the data collapse in Figs. 2 and 4, and the amplitude of the resonant frequency shift and relaxation rate—is obtained by inserting the quasi-thermal ansatz of Eq. (24), f_{ΔL+ε}=f0 e^{-ε/T} and f_{ΔR+ε}=f0 e^{-δΔ/T} e^{-ε/T}, with f0 fixed by the total-density condition, Eqs. (42)–(44). Equations (45), (48), (58), and (60) are closed forms of this one-parameter model. If the real low-energy quasiparticle distribution is non-thermal (a point the paper's own Refs. [35,36] show is unresolved) or differs between the two leads, then: (i) the logarithmic resonance at ω10=δΔ may survive, but its weight and its temperature shape are controlled by f(ε) near the gap edge, not by a single T; (ii) the collapse onto F(T/T*) and G(T/T*) is a consequence of the Boltzmann ansatz, not a generic prediction; and (iii) the proposal to extract the quasiparticle distribution from measurements is not demonstrated—the integral representations (23) and (29) are stated, but no inversion, uniqueness, or sensitivity analysis is given. The same-distribution-in-both-leads assumption (Eq. 5) is also physically nontrivial for weakly coupled leads with different gaps. Thus the central quantitative predictions are conditional on quasiparticles being quasi-thermal with a common effective temperature.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper develops a theory of quasiparticle-induced frequency shifts and relaxation rates in a gap-engineered transmon, i.e., a split transmon whose Josephson junctions connect superconductors with unequal gaps ΔL and ΔR. Using the tunnel Hamiltonian, Fermi's golden rule for the dissipative response, and Kramers-Kronig plus inductive matching for the reactive response, the authors derive closed-form expressions for the finite-frequency admittance (Eqs. (20)–(30)), the transmon frequency shift (Eqs. (45), (48), (52)), the downward relaxation rate (Eqs. (58), (60)), and the parity-switching rate (Eqs. (64), (67)). The central predictions are: (i) a resonant enhancement of the quasiparticle-induced relaxation rate when the qubit frequency ω10 matches the gap difference δΔ, with a logarithmic divergence and a temperature scale T* = |δΔ − ω10|; (ii) a correspondingly narrow and deep dip in the frequency shift on the δΔ > ω10 side; and (iii) universal scaling functions F(T/T*) and G(T/T*) in the quasi-thermal model. The paper also proposes that these effects could be used to probe the low-energy quasiparticle distribution. The derivation is analytic and detailed, with the full Kubo-formalism expression relegated to Appendix A.","tokens_in":22477,"tokens_out":15660,"duration_ms":166916,"significance":"If the quantitative predictions hold, this work provides a useful theoretical foundation for interpreting experiments on gap-engineered transmons and for designing spectroscopic probes of nonequilibrium quasiparticles. The paper's strengths include its transparent microscopic starting point, absence of fitted parameters in the final predictions (xqp and T are inputs), closed-form analytic results, reproduction of the symmetric-junction limit δΔ → 0, and explicit falsifiable predictions such as the resonance at ω10 = δΔ and the collapse of curves onto F(T/T*) and G(T/T*). The main caveat is that the quantitative lineshapes and the proposed distribution-extraction protocol are derived under a quasi-thermal Boltzmann distribution and a common distribution in both leads; this is acknowledged as an assumption but its consequences for the proposed spectroscopy are not fully explored.","major_comments":[{"comment":"The prefactor in Eq. (60) appears to be inconsistent with the low-temperature limit of Eq. (58). For T ≪ T* ≡ |δΔ − ω10|, Eq. (58) gives Γ1→0 ≈ (xqpζ/2π) F [sqrt(Δ/(2T*))] G(T/T*) with F = (ω10(0)^2 + ω10(Φ)^2)/ω10(Φ), i.e., the amplitude contains sqrt(Δ/(2T*)) = sqrt(πΔ/(2πT*)). Eq. (60) instead displays sqrt(Δ/(2πT*)), which is smaller by a factor sqrt(π). This affects the quantitative amplitude and the data collapse in Fig. 4. Please correct the prefactor or clarify the definition of G.","section":"§VI A, Eq. (60)"},{"comment":"The resonant first term of Eq. (48) is a factor of 2 larger than the corresponding T ≪ T* limit of Eq. (45). Setting Φ = 0, so that ω10(0) = ω10(Φ) = ω, and using h(x) ≈ 1/sqrt(πx) for x ≫ 1, Eq. (45) yields δωY10 ≈ −(xqpζω/2π) sqrt(Δ/(2T*)), whereas Eq. (48) yields −(xqpζω/π) sqrt(Δ/(2T*)). The discrepancy apparently arises from the missing factor 1/2 in the flux factor (ω10(0)^2 + ω10(Φ)^2)/(2ω10(Φ)^2) when passing from Eq. (45) to Eq. (48). This changes the claimed resonant amplitude and the inset in Fig. 2; please verify the other regular terms in Eq. (48) as well.","section":"§V B, Eq. (48)"},{"comment":"The paper states that the integral representations, Eqs. (23) and (29), allow one to access the quasiparticle energy distribution function, but no inversion, uniqueness, or sensitivity analysis is provided. The extraction of an 'effective temperature' and the universal collapse onto F(T/T*) and G(T/T*) are consequences of the Boltzmann ansatz Eq. (24), not of a model-independent inversion. As written, the manuscript demonstrates a quasi-thermal model and its observable signatures, but it does not demonstrate that measurements of δω10 and Γ1→0 can reconstruct a general non-thermal distribution. Please either add such an analysis or explicitly limit the proposal to the quasi-thermal case.","section":"§VII and Eqs. (23), (29)"},{"comment":"The assumption that quasiparticles in both leads share the same distribution function, Eq. (5), together with the fixed-total-density condition Eq. (42), implies that the two leads are equilibrated with a common effective temperature and chemical potential on the relevant time scales. For weakly coupled leads with different gaps this is physically nontrivial; if the two leads instead support independent quasiparticle populations, the volume factor ζ(T) and the temperature dependencies of Eqs. (45), (58), and (64) would change. Please discuss the parameter regime in which this equilibration assumption is expected to hold, or present the two-population generalization.","section":"§II, Eq. (5) and §V B, Eqs. (42)–(44)"}],"minor_comments":[{"comment":"The text refers to 'δωY01' in two places ('The temperature dependence of δωY01 in this case' and 'The resonant enhancement of the frequency shift δωY01'); both should read δωY10.","section":"§V B"},{"comment":"There is a typo in the Introduction: 'gap-engieered devices' should be 'gap-engineered devices'.","section":"§I"},{"comment":"The notation in Eq. (29) is somewhat dense: the integration limits 'Z |ω−δ∆| 0' should be written as ∫_0^{|ω−δΔ|} for clarity, and the step-function arguments should be explicitly parenthesized.","section":"§IV B, Eq. (29)"},{"comment":"The figure captions state that ζ(T) = 1 in the limit VL ≫ VR; for Fig. 3, where VL/VR is varied, it would help to state explicitly that the plotted curves use the full ζ(T) of Eq. (44).","section":"§V B, Figs. 2 and 3"},{"comment":"The definition T* = |δΔ − ω10(Φ)| in Eq. (60) is said to 'extend' the earlier definition; for consistency with Eq. (48), where T* = δΔ − ω10, please state clearly that the absolute value is taken in Eq. (60).","section":"§VI A, Eq. (60)"}],"recommendation":"major_revision","confidential_remarks":"The two factor discrepancies in Eqs. (48) and (60) are likely typographical rather than conceptual, since the full expressions Eqs. (45) and (58) appear internally consistent and the Kubo derivation reproduces known limits. Nevertheless, because these equations are the quantitative core of the paper and are used in the main figures, they must be corrected before publication. I would also ask the authors to temper or substantiate the distribution-retrieval claim in Section VII. Once these points are addressed, the paper would be suitable for the journal."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper does something real: it works out the finite-frequency admittance of a Josephson junction with unequal gaps, and from it the quasiparticle-induced frequency shift and relaxation rate of a gap-engineered transmon. The new physics is a resonance when the qubit frequency ω10 matches the gap difference δΔ, controlled by a new temperature scale T*=|δΔ−ω10|. Near the resonance the relaxation rate is enhanced by sqrt(Δ/T*) and the frequency shift develops a dip. These predictions are concrete enough to serve as a spectroscopic probe of the low-energy quasiparticle distribution, which is genuinely useful.\n\nThe derivation is solid. The dissipative part comes from Fermi's golden rule, the reactive part from Kramers-Kronig plus matching to the dc Josephson inductance, and the symmetric-gap limit reduces to Catelani et al. The high-temperature limit also checks out. No fitted parameters; xqp and T are inputs. The integral representations for arbitrary distribution functions are given, which is good practice.\n\nTwo soft spots. First, the abstract says the gap difference can induce an 'anomalous positive frequency shift.' The body states that δωY10 is negative at all parameter values. What they mean is that the shift becomes less negative (has a positive slope) with temperature, but the wording is wrong and will mislead readers. Second, the quantitative closed forms, the universal collapse in Figs. 2 and 4, and the inversion proposal all rest on the quasi-thermal Boltzmann ansatz, Eq. (24), and on the same distribution in both leads, Eq. (5). If the real distribution is non-thermal or differs between leads, the curves change. The paper acknowledges this and provides the integral representations, but it does not demonstrate that a measurement can actually invert those integrals to extract f(ε). So the spectroscopy proposal is a promising idea, not a demonstrated method.\n\nNone of this undermines the central physics. The resonance is robust because it comes from the divergent density of states, not from the thermal ansatz; the ansatz controls the lineshape and weight. This paper deserves a serious referee. I would send it to PRB or similar. I'd also ask the authors to fix the abstract and add a short paragraph on the inversion question.\n\nFor us, it's worth a reading-group slot and I'd cite the admittance result.","headline":"Genuinely new resonance physics in gap-engineered transmons, solidly derived, with an abstract overclaim and an unproven inversion step.","tokens_in":23119,"tokens_out":2840,"would_cite":true,"duration_ms":30431,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"In a gap-engineered transmon, quasiparticles produce a resonant relaxation peak and a frequency-shift dip when the qubit frequency equals the gap difference, giving a spectroscopic window on the quasiparticle energy distribution.","keywords":["quasiparticles","transmon","gap-engineered Josephson junction","qubit relaxation","frequency shift","junction admittance","superconducting qubit spectroscopy","density of states singularity"],"falsifier":"In a split transmon with a known gap difference $\\delta\\Delta$, sweep $\\omega_{10}$ through $\\delta\\Delta$ by magnetic flux at fixed quasiparticle density and measure $\\Gamma_{1\\to 0}(T)$ and $\\delta\\omega_{10}(T)$. The paper's claim is a logarithmic divergence of the rate at $\\omega_{10} = \\delta\\Delta$, an additional $e^{-T_*/T}$ suppression on the $\\delta\\Delta > \\omega_{10}$ side, and a frequency-shift dip of width $T_*$ on that same side; observing none of these under quasi-thermal conditions would falsify the central resonance picture.","tokens_in":21929,"feed_emoji":"⚛️","tokens_out":10862,"duration_ms":124185,"temperature":0.7,"pith_summary":"This paper predicts that in a transmon—a superconducting qubit whose Josephson junction connects two superconductors with different gaps—quasiparticles create sharp, measurable features in the qubit's properties. When the qubit frequency $\\omega_{10}$ matches the gap difference $\\delta\\Delta$, the quasiparticle-induced relaxation rate develops a logarithmic singularity and the frequency shift develops a narrow, temperature-dependent dip whose width is set by $T_* = |\\delta\\Delta - \\omega_{10}|$. Near this resonance, both observables collapse onto universal functions of $T/T_*$. Because a split transmon's frequency is flux-tunable, the resonance can be swept experimentally and used to probe the low-energy part of the quasiparticle distribution, a quantity that is otherwise difficult to access.","feed_headline":"Transmon resonance maps quasiparticle energies","feed_subtitle":"Tuning the qubit across the gap difference yields a resonant rate and a frequency dip that encode the quasiparticle distribution.","key_machinery":"The machinery is the finite-frequency admittance $Y(\\omega)$ of an asymmetric Josephson junction, obtained from a tunnel Hamiltonian with BCS quasiparticles and evaluated through the Kubo formalism and Fermi's golden rule. The load-bearing identity is the spectral density in Eq. (25), built from two square-root density-of-states factors $1/\\sqrt{\\varepsilon}$ at the gap edges; it produces a logarithmic divergence in $\\mathrm{Re}\\,Y$ and a step in $\\mathrm{Im}\\,Y$ at $\\omega = \\delta\\Delta$. The paper then evaluates the split-transmon frequency shift and relaxation rate from this admittance, assuming a quasi-thermal Boltzmann distribution $f_{\\Delta+\\varepsilon} = f_0 e^{-\\varepsilon/T}$ with fixed total density $x_{qp}$, giving closed forms in terms of $H(x) = e^{x/2}K_0(x/2)$ and $h(x) = e^{-x/2}I_0(x/2)$.","core_discovery":"The central claim is that a finite gap difference $\\delta\\Delta$ between the superconductors of a Josephson junction qualitatively changes the quasiparticle back-action on a transmon. The singular density of states at the two gap edges makes the dissipative part of the junction admittance diverge logarithmically at $\\omega = \\delta\\Delta$ while the reactive part jumps, and this singularity propagates into the qubit observables: the relaxation rate $\\Gamma_{1\\to 0}$ is resonantly enhanced by a factor $\\sqrt{\\Delta/|\\delta\\Delta - \\omega_{10}|}$ and the frequency shift develops a dip that is deeper and narrower than in a symmetric junction. The paper gives closed-form expressions for both effects under a quasi-thermal quasiparticle distribution—Eq. (60) for the rate and Eq. (48) for the shift—and shows that near resonance each has a universal temperature dependence, $G(T/T_*)$ for the rate and $F(T/T_*)$ for the shift. It argues that measuring these features, together with the quasi-elastic parity-switching rate, provides a practical spectroscopic method to probe the low-energy quasiparticle distribution.","pith_inferences":["The integral representations in Eqs. (23) and (29) relate the measured shift and rate directly to the distribution function $f_\\varepsilon$; inverting them at several flux settings could recover a non-thermal $f_\\varepsilon$ rather than merely fitting an effective temperature, extending the proposed spectroscopy beyond the quasi-thermal assumption.","A device with a small low-gap volume $V_L \\ll V_R$ concentrates quasiparticles in the low-gap lead and amplifies the resonant signal by a factor $(V_L+V_R)/V_L$ up to a Saha-type ionization temperature $T_{\\mathrm{Saha}}$; this trade-off suggests a design principle for quasiparticle sensors even though it worsens qubit coherence.","If the resident quasiparticle distribution is strongly non-thermal, the predicted collapse of the curves in the paper's Figs. 2 and 4 would fail, and the shape of that failure could itself diagnose the energy dependence of quasiparticle relaxation, a debated quantity the paper discusses.","A wording discrepancy: the abstract's phrase \"anomalous positive frequency shift\" does not match the main text, which states that the admittance-induced frequency shift is negative at all parameter values; the non-monotonic feature shown in the figures is a dip toward more negative values."],"forward_implications":["Sweeping a split transmon's frequency through $\\delta\\Delta$ should reveal a resonant enhancement in $\\Gamma_{1\\to 0}$ whose height grows as $\\sqrt{\\Delta/|\\delta\\Delta - \\omega_{10}|}$, making weak quasiparticle populations measurable.","Near resonance the temperature dependence of the relaxation rate has the universal form $G(T/T_*)$ on the $\\omega_{10} > \\delta\\Delta$ side and $G(T/T_*)e^{-T_*/T}$ on the $\\delta\\Delta > \\omega_{10}$ side, so a single experiment can separate quasiparticle density from spectral shape.","The frequency shift near $\\delta\\Delta > \\omega_{10}$ develops a dip at $T \\sim T_*$ with amplitude $x_{qp}\\sqrt{\\Delta/T_*}$, giving a second, independent observable for the same quasiparticle parameters.","At $T \\ll \\delta\\Delta$ and $\\omega_{10} \\ll \\delta\\Delta$, dissipation is exponentially suppressed and the frequency shift varies on the temperature scale $\\delta\\Delta$ rather than $\\omega_{10}$, confirming the gap-engineering strategy for protecting qubits from quasiparticles.","Parity-switching rates are exponentially suppressed at low temperature and show no resonance at $\\omega_{10} = \\delta\\Delta$, so they can serve as a calibration reference for the resonant rates."],"supporting_citations":[{"why":"Supplies the Fermi golden rule treatment of quasiparticle dissipation in Josephson junctions and the general theory of quasiparticle effects in superconducting qubits.","marker":"[1]"},{"why":"Gives the baseline quasiparticle-induced frequency shift and relaxation rate for a symmetric-junction transmon that the gap-engineered results extend.","marker":"[33]"},{"why":"Describes the gap-engineered transmon experiment this theory is built to explain and reproduces its limiting expression for the gap-suppression frequency shift.","marker":"[26]"},{"why":"Provides the split-transmon architecture with flux-tunable frequency used to sweep $\\omega_{10}$ across $\\delta\\Delta$.","marker":"[34]"},{"why":"Supplies the zero-temperature Ambegaokar-Baratoff formula underlying the Josephson energy and kinetic inductance definitions.","marker":"[39]"},{"why":"Gives the quasiparticle-induced suppression of the superconducting gap used for the $\\delta\\omega^{\\Delta}_{10}$ part of the shift.","marker":"[37]"},{"why":"Provides the Kramers-Kronig relation and the matching of the purely inductive term in the junction admittance.","marker":"[40]"},{"why":"Supplies the demonstration of transmon-based quasiparticle spectroscopy and the use of quasiparticle temperature as a proxy for their typical energy.","marker":"[31]"},{"why":"Gives the symmetric-junction parity-switching rate used as the high-temperature limit and reference point for the gap-engineered rates.","marker":"[47]"}],"fun_headline_variants":["Transmon gap mismatch maps quasiparticle energies","Qubit resonance reads quasiparticle spectrum","Gap-engineered transmon probes quasiparticle temperature","Frequency dip and rate spike reveal quasiparticles","Anomalous shift marks quasiparticle gap difference"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The quantitative formulas assume quasiparticles in both leads share a single quasi-thermal Boltzmann distribution with one effective temperature and that the total quasiparticle density stays fixed as temperature changes; if the real distribution is non-thermal or differs between the two leads, the predicted temperature curves and universal collapse change.","fun_headline_variants_meta":{"raw":{"variants":["Transmon gap mismatch maps quasiparticle energies","Qubit resonance reads quasiparticle spectrum","Gap-engineered transmon probes quasiparticle temperature","Frequency dip and rate spike reveal quasiparticles","Anomalous shift marks quasiparticle gap difference"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.0002,"raw_usage":{"total_tokens":1357,"prompt_tokens":910,"completion_tokens":447,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":526,"completion_tokens_details":{"reasoning_tokens":373}},"tokens_in":526,"tokens_out":447,"duration_ms":5491,"temperature":1.0,"reasoning_tokens":373,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T10:59:52.948853+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"In a split transmon with a known gap difference $\\delta\\Delta$, sweep $\\omega_{10}$ through $\\delta\\Delta$ by magnetic flux at fixed quasiparticle density and measure $\\Gamma_{1\\to 0}(T)$ and $\\delta\\omega_{10}(T)$. The paper's claim is a logarithmic divergence of the rate at $\\omega_{10} = \\delta\\Delta$, an additional $e^{-T_*/T}$ suppression on the $\\delta\\Delta > \\omega_{10}$ side, and a frequency-shift dip of width $T_*$ on that same side; observing none of these under quasi-thermal conditions would falsify the central resonance picture.","supporting_citations":[{"cited_title":"Catelani , author R","cited_arxiv_id":null,"evidence_quote":"Gives the baseline quasiparticle-induced frequency shift and relaxation rate for a symmetric-junction transmon that the gap-engineered results extend."},{"cited_title":"Diamond , author V","cited_arxiv_id":null,"evidence_quote":"Provides the split-transmon architecture with flux-tunable frequency used to sweep $\\omega_{10}$ across $\\delta\\Delta$."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the quasiparticle-induced suppression of the superconducting gap used for the $\\delta\\omega^{\\Delta}_{10}$ part of the shift."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the Kramers-Kronig relation and the matching of the purely inductive term in the junction admittance."},{"cited_title":"Connolly , author P","cited_arxiv_id":null,"evidence_quote":"Supplies the demonstration of transmon-based quasiparticle spectroscopy and the use of quasiparticle temperature as a proxy for their typical energy."},{"cited_title":"Catelani ,\\ 10.1103/PhysRevB.89.094522 journal journal Phys","cited_arxiv_id":null,"evidence_quote":"Gives the symmetric-junction parity-switching rate used as the high-temperature limit and reference point for the gap-engineered rates."}],"review_version":1}