{"id":"6fe9d5d7-8c5f-42b3-9808-bb6099d01b94","arxiv_id":"2505.09322","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Including the superconductor's frequency-dependent surface impedance in circuit QED quantization shows that resonator losses, not just junction capacitance, set the effective light-matter coupling cutoff near the superconducting gap.","lead":"This paper builds a quantum model of superconducting microwave resonators that includes how the metal's surface impedance changes with frequency. It shows that these material effects set the natural cutoff for how strongly a qubit couples to high-frequency resonator modes.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The spectral density J(ω) in Eq. (18) is asserted, and the complex-pole Lamb shift Eq. (19) appears not to reduce to the known dispersionless result as g→0, suggesting a normalization error.","rationale":"The paper contains substantial constructive work: the Lagrangian construction in the Supplement, the Kramers-Kronig reduction, and the Green identity (S52) are nontrivial and internally coherent as far as I can verify. I do not see an outright contradiction that would justify rejection. However, the central quantitative claim depends entirely on the asserted spectral density J(ω) and the contour-integral reduction in Eq. (19). The reader's weakest assumption correctly identifies the missing derivation of J(ω); my read adds a concrete and sharper symptom: the framework does not visibly reduce to the known dispersionless CC result in the small-g limit, and Table I shows a large gap between Δ_dispersion and Δ_nodispersion precisely for the largest width, where g is smallest. This is a red flag rather than a proven error, so I would keep the verdict CONDITIONAL: the paper should either derive Eq. (18) from the quantized fields and confirm the lossless limit, or the central cutoff claim should be revised. I therefore set verdict_should_be to UNCHANGED, meaning the reader's CONDITIONAL verdict remains appropriate.","tokens_in":30364,"tokens_out":24614,"duration_ms":277730,"concrete_test":"Repeat the Table I Lamb-shift calculation for one fixed geometry and material (e.g., Aluminum, s = 20 μm) with g scaled by 10^{-1}, 10^{-2}, 10^{-3}, and 10^{-4}, keeping all other parameters fixed, and verify that Δ_dispersion from Eq. (19) approaches Δ_nodispersion monotonically as g→0. If it instead saturates near −45 MHz or fails to converge to the dispersionless value, the J(ω) weights or the residue normalization in Eqs. (18)–(19) are incorrect. As a complementary check, derive J(ω) explicitly by squaring the reservoir coupling amplitudes in Eq. (12) and integrating over the continuum label x'.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Eqs. (18)–(19) are the only bridge from the quantized fields to the central claim, but J(ω) is introduced with the phrase 'which gives way to' and is not derived from Eq. (12) or the Hamiltonian (4). The residue sum in Eq. (19) is never checked against the one limit where the answer is known. As the geometry factor g→0 (or Zs→0), Eq. (10) reduces to the lossless dispersion relation and Eq. (17) reduces to the current-conserving coupling gn ∝ √ωn of Ref. [3]; hence Δ_dispersion should tend to Δ_nodispersion. Table I does not show this: the smallest-g entry, s = 20 μm, has Δ_dispersion = −47.62 MHz while Δ_nodispersion = −276.16 MHz, a factor ~5.8 discrepancy, and the trend from s = 10 to s = 20 (−51.1 to −47.6 MHz) suggests the g→0 limit lies near −45 MHz, not −276 MHz. If the residue weights or biorthogonal-mode normalization in Eq. (19) are off, the complex-pole machinery can produce a plausible-looking geometry dependence while giving quantitatively wrong Lamb shifts and an incorrect s ≳ 10 μm threshold.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a quantization of a coplanar-waveguide resonator coupled to a charge qubit that includes the frequency-dependent surface impedance of the superconductor. The central object is a resonator-reservoir Hamiltonian whose classical limit recovers a non-Hermitian wave equation with complex refractive index. The authors solve for complex eigenfrequencies as fixed points of the dispersion relation, extend the Mattis-Bardeen conductivity to complex frequencies, and use the resulting spectral density to compute a Lamb shift. Their main claim is that superconducting loss effectively decouples light and matter at frequencies beyond the gap, so that the Lamb-shift partial sums converge more rapidly for narrow center conductors, with the contribution of beyond-gap modes negligible only for s ≳ 10 µm.","tokens_in":30639,"tokens_out":4828,"duration_ms":51132,"significance":"If the quantitative claims survive scrutiny, the result is significant: it identifies a material- and geometry-dependent cutoff for multimode circuit-QED interactions, replacing the common purely capacitive cutoff estimate. The paper has notable strengths: no free parameters are fitted to the computed Lamb shifts; the complex poles are determined as fixed points of the dispersion relation; and the analytic continuation of the Mattis-Bardeen conductivity is worked out in explicit form. The prediction of a geometry-dependent convergence threshold is falsifiable and practically relevant. However, the central quantitative bridge from the quantized fields to the Lamb shift, namely the spectral density in Eq. (18) and the residue sum in Eq. (19), is asserted rather than derived, and the g→0 limit in Table I shows a large unexplained discrepancy against the no-dispersion baseline. These issues prevent me from endorsing the paper's quantitative conclusions in their present form.","major_comments":[{"comment":"The spectral density J(ω) is asserted rather than derived. The text states that Eq. (17) for gn 'gives way to' Eq. (18), but no derivation of the specific combination ω²(gRs/(ωℓm))Re{ϵ}Re{G} + ω²|ϵ|²Im{G} is provided from the field expansions in Eqs. (11) and (12), from Eq. (4), or in the Supplemental Material. The complex poles of the Green's function determine mode frequencies and decay rates, but the residue weights entering the Lamb shift in Eq. (19) are governed by the numerator of J(ω). A different normalization of the dipole moment, the biorthogonal modes, or the Fano diagonalization would change the weight of each pole and hence the values in Figs. 3 and 4. Because the geometry-dependent cutoff claim rests entirely on this bridge, the derivation of Eq. (18) is load-bearing and must be supplied.","section":"Main text, Eq. (18)"},{"comment":"The manuscript never checks the Lamb shift against the known lossless limit, and the table appears inconsistent with it. As g→0, the surface-impedance term in Eq. (10) vanishes, and Eq. (17) reduces to the current-conserving coupling gn ∝ √ωn of Ref. [3]; hence Δ_dispersion should tend to Δ_nodispersion. Table I shows that the smallest-g entry, s = 20 µm, has Δ_dispersion = −47.62 MHz while Δ_nodispersion = −276.164 MHz, a factor of about 5.8, and the trend from s = 10 µm (−51.07 MHz) to s = 20 µm (−47.62 MHz) suggests that the g→0 limit lies near −45 MHz, not −276 MHz. This discrepancy indicates a normalization or residue-weight error in Eq. (19), or an incorrect proportionality factor in Eq. (18). This issue must be resolved before the central claim of a geometry-dependent cutoff can be accepted.","section":"Table I and the g→0 limit"},{"comment":"The analytic continuation of the conductivity drops two integrals on the basis of an unquantified numerical statement: 'Numerical integration showed that the contribution of the last two integrals in (S79) is negligible, hence we have dropped them from further analysis.' No parameter range, error estimate, or supporting plot is given. Since Eqs. (15) and (16) are used to compute the complex fixed points {ω_n, γ_n} that feed Eq. (19), this truncation is load-bearing. The authors should provide a quantitative bound or a numerical comparison showing that the dropped terms are small over the full range of frequencies, widths s, and damping rates used in Figs. 3 and 4.","section":"Supplement, 'Extreme Anomalous Limit', Eq. (S81)"},{"comment":"The contour-integral evaluation leading to Eq. (19) is not shown. The result is stated as a sum over complex poles of the Green's function plus mirrored terms, but the derivation should also justify the neglect of possible branch-cut contributions from the analytically continued conductivity and from ϵ(ω). In addition, the expression should be validated against the known result in the limit g→0, which would also resolve the discrepancy noted above. As written, Eq. (19) is the only quantitative output of the paper, and its derivation cannot be left as an assertion.","section":"Main text, Eq. (19)"}],"minor_comments":[{"comment":"The sentence 'leading to É ∼ O(10^3) GHz' should read 'leading to a cutoff frequency ω ∼ 10^3 GHz'; also, the following sentence refers to 'Equation (1) is non-hermitian', but the non-Hermitian object is the wave equation (2), not Eq. (1).","section":"Main text, Eq. (1) and surrounding text"},{"comment":"There are several typos and garbled symbols: 'josephon junction' should be 'Josephson junction'; 'Kirchoff's laws' should be 'Kirchhoff's laws'; 'Kramer-Kronig' should be 'Kramers-Kronig'; 'off-set' should be 'offset'; 'de-tuning' should be 'detuning'; and 'Én f 2∆/ℏ' and 's g 10µm' should be 'ω_n < 2Δ/ℏ' and 's ≳ 10 µm'.","section":"Throughout"},{"comment":"The dashed horizontal line is said to denote the 70% convergence, but it is not specified in the caption which quantity the partial sums are normalized by and why the normalization denominator is the 2500-term sum; this should be stated explicitly.","section":"Fig. 4 and accompanying text"},{"comment":"References [11] and [33] are the same dissertation by J. Gao; citing the same work twice with different numbers is confusing and should be consolidated.","section":"References"},{"comment":"The term 'di-magnetic term' should be 'diamagnetic term', and the notation should be defined where first introduced to avoid confusion with the magnetic inductance ℓm.","section":"Supplement, Eq. (S5)"}],"recommendation":"major_revision","confidential_remarks":"The strongest concern is the g→0 limit in Table I: a factor-of-5.8 discrepancy between Δ_dispersion and Δ_nodispersion at the smallest geometric factor suggests a normalization or residue-weight error rather than a physical effect. If the authors can supply the missing derivation of Eq. (18), a proper contour-integration derivation of Eq. (19), and a validated g→0 limit, the paper could become publishable. I would ask the editor to require those additions before sending the manuscript back for another round, rather than accepting the current form."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know up front. This is a genuine extension of cQED quantization, not a repackaging: the paper puts the complex, frequency-dependent Mattis-Bardeen surface impedance into the current-conserving Hamiltonian and continues the conductivity into the complex plane to find the complex poles that set the light-matter cutoff. The supplement does real work—the distributed-element Lagrangian, the Fano diagonalization, the Green's-identity commutator, and a long elliptic-integral evaluation. No free parameters are fitted, and the physical idea, that quasiparticle loss above 2Δ acts as a natural cutoff, is sensible. But the quantitative claim—that the cutoff is geometry-dependent with a threshold near s ≳ 10 µm—rests on an asserted spectral density and a limit that is never checked. I would not take the numbers at face value until that is resolved.\n\nSoft spots, in proportion. First, Eq. (18), the spectral density J(ω) feeding the Lamb shift, is introduced with \"which gives way to\" rather than derived from the Hamiltonian. The residue sum in Eq. (19) inherits its normalization, and the supplement derives the mode expansions but stops short of producing J(ω). A referee is entitled to see that step.\n\nSecond, the stress-test concern about the lossless limit lands, in my reading. As g→0 (or Zs→0), the framework should reduce to the no-dispersion result of Ref. [3]. Table I suggests otherwise: the smallest-g row, s=20 µm, gives Δ_dispersion = −47.62 MHz against Δ_nodispersion = −276 MHz, a factor of six. One can construct a defense: the beyond-gap modes are a damped continuum whose coupling vanishes with Rs, while the no-dispersion model retains discrete modes with a capacitive cutoff, so the two need not match. But the paper never makes that argument, and the trend in s points to a small-g limit near −45 MHz, far from −276 MHz. If the residue weights in Eq. (19) are normalized incorrectly, the width dependence—and the s ≳ 10 µm threshold—loses its support.\n\nMinor: the supplement drops two integrals after Eq. (S81) as \"negligible\" without numbers, and there are typos (\"josephon\", \"di-magnetic\", Eq. 1 vs Eq. 2).\n\nWho this is for: cQED theorists working on multimode effects, Lamb shifts, and wide-band devices. The framework deserves engagement; the specific predictions do not yet. I would send it to peer review—the novelty and the technical substance justify referee time—but the referee should push hard on the derivation of J(ω) and on the lossless limit. Expect a major revision, not a desk reject.","headline":"A genuine extension of cQED quantization with a sensible physical claim; the central Lamb-shift calculation rests on an asserted spectral density and an unchecked lossless limit.","tokens_in":31155,"tokens_out":12897,"would_cite":false,"duration_ms":119414,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Circuit QED's light-matter coupling cutoff is set by superconducting dispersion, not by the junction capacitance, with loss above the gap effectively decoupling qubit and resonator.","keywords":["circuit quantum electrodynamics","surface impedance","superconducting gap","Lamb shift","light-matter coupling cutoff","charge qubits","coplanar waveguide resonators","BCS conductivity"],"falsifier":"Re-derive $J(\\omega)$ by Fano-diagonalizing the Hamiltonian in Eq. (4) and compare the complex-pole weights with Eq. (18); if the weights differ, the width dependence in Figs. 3 and 4 is not secured. Experimentally, one could measure the qubit Lamb shift at fixed qubit frequency for Aluminum resonators with center-conductor widths 1.2 µm and 20 µm, where the paper predicts roughly $-95$ MHz and $-48$ MHz, so a substantially different ratio would refute the central claim.","tokens_in":30168,"feed_emoji":"⚛️","tokens_out":7282,"duration_ms":72408,"temperature":0.7,"pith_summary":"This paper sets out to show that circuit-QED quantization must include the superconductor's frequency-dependent surface impedance, because that impedance controls how strongly a qubit couples to high-frequency resonator modes. The central claim is that above the superconducting gap, the loss that breaks Cooper pairs also makes the qubit and resonator effectively decouple, so the light-matter interaction has a material- and geometry-dependent cutoff rather than one fixed by the junction capacitance. A sympathetic reader would care because off-resonant modes contribute to Lamb shifts and multi-qubit couplings, and current models treat the resonator as ideal over the whole frequency range. The authors demonstrate the effect by computing Lamb shifts for Aluminum and Niobium resonators of varying center-conductor width, finding a roughly fourfold drop in total Lamb shift for Aluminum between 0.6 µm and 20 µm widths.","feed_headline":"Beyond the gap, superconducting loss decouples qubit and resonator","feed_subtitle":"Dispersion, not junction capacitance, sets the cutoff; narrow resonators show much larger Lamb shifts.","key_machinery":"The central object is the frequency-dependent surface impedance $Z_s(\\omega)=R_s(\\omega)+iX_s(\\omega)$, which enters the transmission-line description through the complex refractive index $\\epsilon(\\omega)=1+\\frac{g Z_s(\\omega)}{i\\omega\\ell_m}$. The paper analytically continues the BCS conductivity to complex frequencies so that the dispersion relation has complex fixed points $\\omega_n=\\Omega_n+i\\Gamma_n$, and it constructs the Green's function of the resulting non-Hermitian wave equation from bi-orthogonal eigenmodes. That Green's function feeds the spectral density $J(\\omega)$ used in the Lamb-shift integral, converting the lossy broadened modes into a sum over complex poles.","core_discovery":"The paper claims that when the frequency-dependent surface impedance of the superconductor is included in the canonical quantization of a coplanar-waveguide resonator, resonator modes above the superconducting gap acquire finite lifetimes, and this loss effectively decouples the light and matter degrees of freedom at high frequencies. Consequently, the effective cutoff of the qubit-resonator coupling is set by material dispersion and resonator geometry, not by the series capacitance of the qubit junction. The authors further claim that the contribution from beyond-gap modes to the Lamb shift is negligible only for center-conductor widths of roughly 10 µm or larger, while narrow resonators show much larger total Lamb shifts and faster convergence of the partial-sum series.","pith_inferences":["If the spectral density in Eq. (18) is correct, the same geometry-dependent decoupling should appear in other observables that sum over high-frequency modes, such as dispersive qubit shifts and inter-qubit couplings; the paper does not compute these.","The predicted red-shift of high modes and the appearance of extra modes near the gap could be probed directly in broadband transmission measurements on Niobium resonators, a signature the paper mentions but does not calculate.","For high-kinetic-inductance materials such as granular aluminum, where the surface impedance is much larger, the effective cutoff should move to lower frequencies and dispersion effects should be stronger; this is an extension beyond the paper's Aluminum and Niobium examples.","The results suggest resonator center-conductor width can be used as a design knob to suppress beyond-gap contributions without changing the qubit itself."],"forward_implications":["The effective qubit-resonator coupling cutoff is determined by material dispersion and resonator geometry, not by the junction capacitance that previous models used.","Beyond-gap modes can be safely ignored only for wide center conductors ($s \\gtrsim 10\\,\\mu\\mathrm{m}$); for narrow resonators they contribute substantially to the Lamb shift.","Narrower resonators exhibit larger total Lamb shifts, for example about $-212.5$ MHz at 0.6 µm versus about $-47.6$ MHz at 20 µm for the Aluminum case considered.","The convergence of Lamb-shift partial sums is faster for narrow resonators, so truncating the mode sum is safer there despite the larger total shift.","Material-specific behavior differs between the extreme-anomalous limit relevant to Aluminum and the dirty limit relevant to Niobium, including extra modes appearing near the gap."],"supporting_citations":[{"why":"Supplies the current-conserving Hamiltonian and the coupling scaling $g_n \\propto \\sqrt{\\omega_n}$ that the paper generalizes to include surface impedance.","marker":"[3]"},{"why":"Provides the quasi-TEM field-penetration picture and the geometric factor relating CPW geometry to the surface-impedance contribution.","marker":"[33]"},{"why":"Supplies the BCS conductivity model whose analytic continuation to complex frequencies underlies the lossy-mode calculation.","marker":"[35]"},{"why":"The supplemental material derives the resonator-reservoir Hamiltonian, the Green's identity, the equal-time commutation relations, and the complex-plane conductivity extension.","marker":"[39]"},{"why":"Provides the Green's-function treatment of local resonator perturbations that the paper adapts to the lossy, dispersive wave equation.","marker":"[40]"},{"why":"Supplies the spectral-density form used to express the Lamb shift in circuit QED.","marker":"[43]"},{"why":"Gives the broadband Lamb-shift integral against the spectral density used in Eq. (19).","marker":"[44]"}],"fun_headline_variants":["Dispersion, not capacitor, dictates circuit QED cutoff","Superconducting loss decouples qubit from high-frequency modes","Narrow resonators amplify Lamb shift, loss sets interaction cutoff","Material loss, not junction capacitance, ends qubit-photon coupling","Beyond-gap modes fade, qubit coupling cut by material loss"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the spectral density $J(\\omega)$ stated in Eq. (18), introduced with the phrase \"which gives way to,\" is the correct coupling-weighted mode density for a lossy, dispersive resonator; the paper does not derive it from the resonator-reservoir Hamiltonian.","fun_headline_variants_meta":{"raw":{"variants":["Dispersion, not capacitor, dictates circuit QED cutoff","Superconducting loss decouples qubit from high-frequency modes","Narrow resonators amplify Lamb shift, loss sets interaction cutoff","Material loss, not junction capacitance, ends qubit-photon coupling","Beyond-gap modes fade, qubit coupling cut by material loss"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000247,"raw_usage":{"total_tokens":1459,"prompt_tokens":777,"completion_tokens":682,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":393,"completion_tokens_details":{"reasoning_tokens":593}},"tokens_in":393,"tokens_out":682,"duration_ms":6744,"temperature":1.0,"reasoning_tokens":593,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T21:35:10.926699+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-derive $J(\\omega)$ by Fano-diagonalizing the Hamiltonian in Eq. (4) and compare the complex-pole weights with Eq. (18); if the weights differ, the width dependence in Figs. 3 and 4 is not secured. Experimentally, one could measure the qubit Lamb shift at fixed qubit frequency for Aluminum resonators with center-conductor widths 1.2 µm and 20 µm, where the paper predicts roughly $-95$ MHz and $-48$ MHz, so a substantially different ratio would refute the central claim.","supporting_citations":[{"cited_title":"Silveri, S","cited_arxiv_id":null,"evidence_quote":"Gives the broadband Lamb-shift integral against the spectral density used in Eq. (19)."},{"cited_title":"Malekakhlagh and H","cited_arxiv_id":null,"evidence_quote":"Supplies the current-conserving Hamiltonian and the coupling scaling $g_n \\propto \\sqrt{\\omega_n}$ that the paper generalizes to include surface impedance."},{"cited_title":"Malekakhlagh, A","cited_arxiv_id":null,"evidence_quote":"Provides the Green's-function treatment of local resonator perturbations that the paper adapts to the lossy, dispersive wave equation."},{"cited_title":"H¨ ummer, F","cited_arxiv_id":null,"evidence_quote":"Supplies the spectral-density form used to express the Lamb shift in circuit QED."}],"review_version":1}