Pith. sign in

REVIEW 3 major objections 5 minor 65 references

Adaptive spectroscopy reveals qubit defects switching and drifting on second timescales.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

T0 review · deepseek-v4-flash

2026-08-04 15:21 UTC pith:QVSLYFIN

load-bearing objection Likely real observation of seconds-scale TLS dynamics, but the quantitative timescales are not secure without a static-TLS null control. the 3 major comments →

arxiv 2608.02086 v1 pith:QVSLYFIN submitted 2026-08-03 quant-ph cond-mat.mes-hall

Adaptive Spectroscopy of Fast Two-Level-System Dynamics in Superconducting Qubits

classification quant-ph cond-mat.mes-hall
keywords two-level systemsspectral diffusiontelegraphic switchingadaptive spectroscopysuperconducting qubitstransmonrandomized benchmarkingrelaxation rate
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper claims that in flux-tunable superconducting qubits, individual two-level-system defects change their transition frequency on timescales of seconds—roughly 300 times faster than conventional nonadaptive spectroscopy has been able to see. Using an FPGA-based adaptive Bayesian estimator of the qubit decay rate, the authors produce time-resolved maps of relaxation rate versus frequency and time, and extract a telegraphic switching time of about 2.2 seconds and a spectral diffusivity of about 0.93 MHz²/s. They also interleave randomized benchmarking and find that gate infidelity tracks the relaxation-rate fluctuations when the qubit is biased on a TLS resonance. If correct, these results mean that TLS-aware calibration and error mitigation in superconducting processors must operate on much shorter timescales than the hours-long intervals currently used.

Core claim

The central discovery is that TLS-induced energy-relaxation peaks in flux-tunable transmons are not static on the timescales of conventional spectroscopy: a single TLS switches between two transition frequencies with correlation time 2.2 seconds, and another TLS diffuses in frequency with diffusivity 0.93 MHz²/s. These dynamics become observable because each point in the relaxation-rate map is obtained in about 1.5 ms by adaptive Bayesian estimation, making a full sweep of 500 frequencies take about 0.78 seconds. The authors argue that the same fast dynamics appear in two independently fabricated device platforms, and that the fluctuation of the relaxation rate at a TLS resonance is correlat

What carries the argument

The key machinery is the real-time adaptive Bayesian estimator of the decay rate Γ1, implemented on an FPGA. The controller stores the posterior over Γ1 as a two-parameter distribution, chooses each probe wait time τ = 0.1/Γ̂ from the current estimate, updates after every single-shot readout, and reinitializes the prior at each flux point. Repeated sweeps of the flux-tunable qubit frequency turn these estimates into time-resolved TLS maps; a two-state hidden Markov model and an Ornstein-Uhlenbeck fit to the tracked trajectory then convert the maps into a switching time and a diffusivity.

Load-bearing premise

The results stand or fall on whether the adaptive estimator's Γ̂ values at N=30 probe cycles faithfully reflect the qubit's true decay rate at each frequency point; if estimator noise is mistaken for TLS motion, the inferred switching and diffusion are artifacts.

What would settle it

Repeat the same adaptive sweep while the qubit is parked on a TLS that is independently known to be static, for example verified by repeated conventional T1 measurements at much slower sampling; if the adaptive map still shows telegraphic switching or diffusion of that feature, or if the extracted timescales shift when N is raised from 30 to several hundred, the reported dynamics are tool artifacts rather than TLS motion.

Watch this falsifier. Get emailed when new claim-graph text bears on it.

If this is right

  • TLS-induced relaxation fluctuations in flux-tunable transmons occur on seconds timescales, so standard T1 tracking with minutes-to-hours sampling averages over the very dynamics that limit gate performance.
  • Spectroscopy with sub-second sampling can map these dynamics across a 500 MHz band in under a second per sweep, making two-dimensional TLS maps practical for routine qubit characterization.
  • The same seconds-scale dynamics appear in two independently fabricated, differently architected devices, suggesting a generic property of current transmon TLS environments rather than a single-device quirk.
  • Gate infidelity measured by randomized benchmarking tracks the relaxation-limited prediction when the qubit is on a TLS resonance, while a few MHz away the correlation vanishes, so fast TLS motion is a real gate-error contributor.
  • The low-latency online estimate of Γ1 is fast enough to be embedded in adaptive control loops for error mitigation, for example pausing operation when the estimated relaxation rate is high.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Inference: the same interleaving scheme could be extended to multi-qubit operations, where a TLS hitting one qubit during a two-qubit gate may produce correlated errors that the adaptive map would reveal in real time.
  • Inference: if confirmed on more devices, the measured diffusivity of about 0.9 MHz²/s and switching rates provide quantitative targets for microscopic models of TLS spectral diffusion, such as coupling to a fluctuating defect bath.
  • Inference: a direct stress test—tracking a known static TLS with the same N=30 estimator settings and with higher N—would separate physical TLS motion from estimator-induced fluctuations; the paper does not report such a control.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

3 major / 5 minor

Summary. The paper reports an FPGA-implemented adaptive Bayesian relaxation-rate estimation protocol applied to flux-tunable transmon qubits, producing time-resolved TLS loss maps with per-column sampling times near 0.78 s. From these maps it claims (i) telegraphic TLS switching with characteristic time τsw ≈ 2.2 s for one feature, (ii) spectral diffusion with D ≈ 0.93 MHz²/s for a second feature, stated as roughly 300× faster than previous nonadaptive spectroscopy, (iii) qualitatively similar dynamics on a second, independently fabricated device in a different laboratory, and (iv) a statistically significant correlation between tracked Γ1 fluctuations and randomized-benchmarking gate infidelity at the TLS resonance (ρ0 = 0.618) that is absent a few MHz away (ρ0 = −0.0507). The statistical treatment in the Supplemental Material includes an HMM with forward-filtering backward-sampling, credible intervals, block-bootstrap uncertainties, a permutation test, and a 3×3×3 robustness grid.

Significance. If the measurements are faithful, the paper opens a new observational regime: it would show that individual TLSs in flux-tunable transmons undergo telegraphic and diffusive spectral motion on seconds timescales, with direct implications for TLS-aware calibration, qubit screening, and error mitigation. The paper has real strengths: the adaptive estimator achieves a substantial speedup over conventional T1 measurements, the statistical analyses are unusually careful (posterior credible intervals, bootstrap, permutation tests, robustness grid), the central RB correlation is tested against an off-TLS control, and the effect is reproduced on two devices in different laboratories. The main open risk is estimator fidelity: whether the reported fast dynamics could be injected by the adaptive estimation procedure itself rather than by the TLS. This is a correctness question, not a stylistic one, and it is the primary reason the paper needs revision.

major comments (3)
  1. [Sec. II; Supp. Eqs. (S.18)–(S.22), (S.4)–(S.7)] The central claim—that the time-dependent features are physical TLS motion—requires a validation that the estimator does not generate spurious fast fluctuations. The protocol uses only N = 30 adaptive probe cycles per frequency point and reinitializes the prior at every flux-time point. No null control is reported. The problem is concrete: Eq. (S.18) subtracts the time median at each frequency, so a perfectly static TLS is removed by construction and only estimation noise remains; the Viterbi tracker of Eq. (S.19) will then produce a center-frequency trajectory even for a static TLS, and fitting its MSD to the OU form of Eq. (S.22) will return a nonzero D. Similarly, the two-component Gaussian HMM of Eqs. (S.4)–(S.7) can split a noise-only score sequence into two clusters and yield a nonzero switching rate. The off-TLS RB control in Supp. Sec. IV validates the correlation analysis, not t
  2. [Fig. 2; Eqs. (S.3), (S.14), (S.23)] The two headline numbers are each extracted from a single TLS feature, and the frequency bands B0 and B1 in Eq. (S.3) are chosen after inspecting the same data. The abstract and discussion state these as 'a characteristic timescale of a few seconds' and D ≈ 0.9 MHz²/s, which reads as a general property of TLS dynamics. The cross-device reproducibility in Fig. 3(a) is qualitative and is not analyzed with the same HMM/OU pipeline. A single feature is sufficient for an observation, but for a claimed characteristic timescale the authors should either apply the same analysis to all resolvable features across both devices and report the distribution, or explicitly temper the wording to 'one observed feature' for each quantitative value. The post hoc band choice should also be acknowledged as a selection effect; the current uncertainty intervals (S.14) and (S.23) do not include this selection.
  3. [Supp. Sec. II C] The Allan-deviation comparison is described in the text as an 'independent check' of the HMM rate, but it is not independent: the credible band in Fig. S7 is generated by simulating trajectories from the posterior draws of (k01, k10), so the agreement of the empirical curve with that band is a posterior predictive consistency check, not a validation against a null or against an independent estimator. This does not invalidate the HMM analysis, but the wording should be corrected, and the analysis should be supplemented by the null control requested above before the Allan-deviation curve can be used as supporting evidence.
minor comments (5)
  1. [Sec. III A; Supp. Sec. III] The 300× comparison with previous diffusivities should be made reproducible: please give the exact previous D values used, the conversion to the same short-time MSD convention, and the references for each number. As written, this claim is difficult to verify.
  2. [Supp. Eq. (S.4) and Fig. S3] The raw band-score correlation ρ = −0.36 and integrated-weight correlation ρ = −0.52 are not exact conservation checks. The text should avoid implying quantitative spectral-weight transfer between the two TLS states unless a conservation test with uncertainties is provided.
  3. [Fig. 3(b)] The moving-average correlation is computed after smoothing with a 5 s window. This is appropriate for trend-level comparison, but the p-value should be interpreted as testing correlation of the smoothed traces, not of the raw single-repetition values. A sentence clarifying this would help.
  4. [Abstract and Sec. IV] The phrase 'redefine the timescales relevant to TLS-aware characterization' is stronger than the single-feature evidence supports. I suggest a more measured formulation such as 'reveal a previously inaccessible regime.'
  5. [General] A data-availability statement and, if possible, release of the analysis code would strengthen reproducibility, given that the full pipeline (smoothing, median subtraction, Viterbi tracking, HMM, bootstrap) is described only in prose.

Circularity Check

0 steps flagged

No significant circularity: empirical measurement with independent cross-checks; self-cited estimator is an external tool.

full rationale

The paper's central claims are empirical observations of TLS dynamics, not derivations from assumptions. The adaptive estimator is taken from the authors' prior work (Ref. [27]), but that prior work is a published, independently validated measurement method; using it does not make the observed dynamics circular, because the TLS maps are new data and are not constructed to equal the estimator's assumptions. The spectral-diffusion analysis fits an OU model to the MSD of a Viterbi-tracked trajectory, and the telegraphic analysis fits a two-state HMM; these are standard descriptive model fits, not predictions forced by the choice of model, and the paper includes robustness checks and an Allan-deviation cross-check. The RB correlation compares measured gate infidelity to an independently estimated relaxation-limited contribution from interleaved T1 data, with an off-TLS null control; it is not a fit to the same data. The absence of a static-TLS null control is a potential correctness risk (estimator noise could masquerade as dynamics), but it is not a circularity: no equation in the paper reduces a claimed prediction to an input by definition.

Axiom & Free-Parameter Ledger

5 free parameters · 6 axioms · 0 invented entities

The paper introduces no new physical entities. Its accounting rests on a set of analysis choices (OU model, Viterbi penalty, HMM bands) and on the fidelity of the prior adaptive method. The central measured quantities D and tau_sw are outputs of fits, not free parameters, but the analysis hyperparameters are listed above.

free parameters (5)
  • OU correlation time tau_c = 32 s (68% CI 21-45 s; 95% CI 13-61 s)
    Fitted to the MSD of the tracked TLS trajectory in the OU model; a model output, not a central claim.
  • Viterbi jump penalty lambda = 10
    Chosen by hand for trajectory extraction; robustness grid (5, 10, 20) shows modest effect, so it is not critical.
  • Bayesian wait-time constant c = 0.1
    From Ref [27]; sets the adaptive wait time tau=c/hat_Gamma1. Not fitted here, but load-bearing for the estimator's behavior.
  • Telegraphic band edges B0, B1 = [-335,-285] MHz and [-230,-185] MHz
    Chosen post hoc from the data to define the two TLS states; HMM results depend on this choice, and no band sensitivity analysis is provided.
  • HMM emission and transition parameters = Fitted by Baum-Welch (500 iterations)
    Gaussian means, covariances, and transition matrix fitted to the band-score sequence; part of the statistical analysis, not the physics claim per se.
axioms (6)
  • domain assumption Exponential-decay likelihood with misclassification probabilities alpha and beta
    The Bayesian update assumes a single exponential decay during each wait time; if TLS switching occurs within a probe cycle, estimates would be biased.
  • domain assumption Quasistatic approximation for TLS during each probe cycle
    The likelihood treats Gamma1 as constant over the ~10-100 us wait time; the authors assume this holds.
  • ad hoc to paper Two-parameter posterior approximation of Ref [27] is faithful
    The FPGA stores only two parameters for the posterior; this is a simplification from the prior work and is not independently validated here.
  • domain assumption Ornstein-Uhlenbeck model for spectral diffusion
    The MSD is fit to the OU form MSD=2D/gamma (1 - exp(-gamma tau)); deviations could bias D.
  • domain assumption Viterbi tracker with penalty lambda recovers the true TLS trajectory
    The extracted frequency trajectory is the output of a regularized optimization; robustness grid suggests stability.
  • domain assumption Observed fluctuations are TLS dynamics rather than qubit-frequency noise
    Interleaved Ramsey checks keep frequency drift within +/-1 MHz, supporting the interpretation but not eliminating other noise sources.

pith-pipeline@v1.3.0-daily-deepseek · 19037 in / 13973 out tokens · 94986 ms · 2026-08-04T15:21:25.867723+00:00 · methodology

0 comments
read the original abstract

Parasitic two-level-system (TLS) defects are a major source of energy relaxation and temporal instability in superconducting quantum processors. Our sub-second adaptive spectroscopy reveals telegraphic switching of TLSs with a characteristic timescale of a few seconds and spectral diffusion with diffusivity $D \approx 0.9~\mathrm{MHz}^2/\mathrm{s}$. These timescales are about $3 \times 10^2$ times faster than what is observed in conventional nonadaptive spectroscopy, which typically requires hours of measurement time. We resolve such fast dynamics on a field-programmable gate array (FPGA)-based controller that enables measurement of frequency- and time-resolved relaxations with sub-second temporal resolution in flux-tunable superconducting qubits. We observe similar defect dynamics across multiple qubits in independently fabricated devices measured in different laboratories. We correlate TLS-induced fluctuations with gate-level errors using randomized benchmarking. Our results reveal a previously inaccessible regime of frequency-resolved TLS dynamics and redefine the timescales relevant to TLS-aware characterization and calibration of superconducting quantum processors.

Figures

Figures reproduced from arXiv: 2608.02086 by Bethany M. Niedzielski, David Pahl, Fabrizio Berritta, Gabriel Cutter, Jan A. Krzywda, Jeffrey A. Grover, Kyle Serniak, Lukas Pahl, Max Hays, Michael Gingras, Mollie E. Schwartz, Paul Buttles, Robert McDermott, Shravan Patel, Spencer Weeden, Stanislav Eilhart, William D. Oliver, William P. Banner.

Figure 1
Figure 1. Figure 1: (c), which is a zoom-in of the data in [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: (a) shows the controller estimates Γˆ 1 in Q1. The controller estimates Γ1 using N = 30 probe cy￾cles [47] at each value of Φ, sweeping the flux so that fQ1 spans a 500 MHz detuning range from the upper degen￾eracy point in 1 MHz steps. Each Γ1 estimate requires only about 1.5 ms of laboratory time. Every frequency sweep is interleaved with a Ramsey experiment at Φ = 0 to verify afterwards that the qubit f… view at source ↗
Figure 3
Figure 3. Figure 3: (a), is qualitatively similar in terms of timescales to [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

65 extracted references · 4 canonical work pages · 3 internal anchors

  1. [1]

    Müller, J

    C. Müller, J. H. Cole, and J. Lisenfeld, Towards under- standing two-level-systems in amorphous solids: insights from quantum circuits, Reports on Progress in Physics 82, 124501 (2019)

  2. [2]

    Siddiqi, Engineering high-coherence superconducting qubits, Nature Reviews Materials6, 875 (2021)

    I. Siddiqi, Engineering high-coherence superconducting qubits, Nature Reviews Materials6, 875 (2021)

  3. [3]

    C. E. Murray, Material matters in superconducting qubits, Materials Science and Engineering: R: Reports 146, 100646 (2021)

  4. [4]

    J. S. Rojas-Arias, A. Noiri, P. Stano, T. Nakajima, J. Yoneda, K. Takeda, T. Kobayashi, A. Sammak, G. Scappucci, D. Loss,et al., Spatial noise correlations beyond nearest neighbors in 28 Si/Si-Ge spin qubits, Physical Review Applied20, 054024 (2023)

  5. [5]

    F. Ye, A. Ellaboudy, D. Albrecht, R. Vudatha, N. T. Jacobson, and J. M. Nichol, Characterization of individ- ual charge fluctuators in Si/SiGe quantum dots, Physical Review B110, 235305 (2024)

  6. [6]

    Donnelly, J

    M. Donnelly, J. Rowlands, L. Kranz, Y. Hsueh, Y. Chung, A. Timofeev, H. Geng, P. Singh-Gregory, S. Gorman, J. Keizer,et al., Noise correlations in an atom-based quantum dot array, Physical Review Applied 23, 064058 (2025)

  7. [7]

    J. S. Rojas-Arias, A. Noiri, J. Yoneda, P. Stano, T. Nakajima, K. Takeda, T. Kobayashi, G. Scappucci, S. Tarucha, and D. Loss, Inferring charge-noise source locations from correlations in spin qubits, Physical Re- view Letters136, 027001 (2026)

  8. [8]

    Lisenfeld, A

    J. Lisenfeld, A. Bilmes, S. Matityahu, S. Zanker, M. Marthaler, M. Schechter, G. Schön, A. Shnirman, G. Weiss, and A. V. Ustinov, Decoherence spectroscopy with individual two-level tunneling defects, Scientific re- ports6, 23786 (2016)

  9. [9]

    Schlör, J

    S. Schlör, J. Lisenfeld, C. Müller, A. Bilmes, A. Schnei- der, D. P. Pappas, A. V. Ustinov, and M. Weides, Cor- relating decoherence in transmon qubits: Low frequency noise by single fluctuators, Physical Review Letters123, 190502 (2019)

  10. [10]

    Y. Gao, Y. Zhang, H. Xu, P. Shi, F. Li, Y. Feng, W. Sun, J. Ding, Y. Liu, H. Wang,et al., Non- local and non-Markovian effects of a microscopic two- level defect in superconducting quantum circuits (2026), 10.48550/arXiv.2605.23385

  11. [11]

    J. F. Kam, S. Gicev, K. Modi, A. Southwell, and M. Us- man, Detrimental non-Markovian errors for surface code memory (2024), 10.48550/arXiv.2410.23779

  12. [12]

    Hakoshima, Y

    H. Hakoshima, Y. Matsuzaki, and S. Endo, Relationship between costs for quantum error mitigation and non- Markovian measures, Physical Review A103, 012611 (2021)

  13. [13]

    Shalibo, Y

    Y. Shalibo, Y. Rofe, D. Shwa, F. Zeides, M. Neeley, J. M. Martinis, and N. Katz, Lifetime and coherence of two- level defects in a Josephson junction, Physical Review Letters105, 177001 (2010)

  14. [14]

    Colao Zanuz, Q

    D. Colao Zanuz, Q. Ficheux, L. Michaud, A. Orekhov, K. Hanke, A. Flasby, M. Bahrami Panah, G. J. Nor- ris, M. Kerschbaum, A. Remm, F. Swiadek, C. Hellings, S. Lazăr, C. Scarato, N. Lacroix,et al., Mitigating losses of superconducting qubits strongly coupled to defect modes, Physical Review Applied23, 044054 (2025)

  15. [15]

    Lisenfeld, A

    J. Lisenfeld, A. Bilmes, A. Megrant, R. Barends, J. Kelly, P. Klimov, G. Weiss, J. M. Martinis, and A. V. Ustinov, Electric field spectroscopy of material defects in trans- mon qubits, npj Quantum Information5, 105 (2019)

  16. [16]

    C. R. H. McRae, G. M. Stiehl, H. Wang, S.-X. Lin, S. A. Caldwell, D. P. Pappas, J. Mutus, and J. Combes, Repro- ducible coherence characterization of superconducting quantum devices, Applied Physics Letters119, 100501 (2021)

  17. [17]

    Weeden, D

    S. Weeden, D. Harrison, S. Patel, M. Snyder, E. Black- well, G. Spahn, S. Abdullah, Y. Takeda, B. Plourde, J. Martinis, and R. McDermott, Statistics of strongly coupled defects in superconducting qubits, Physical Re- view Applied25, 044050 (2026)

  18. [18]

    O. F. Wolff, H. Mantry, R. Raja, W.-H. Peng, K. Sin- girikonda, S. Lee, S. Sudhaman, R. Goncalves, P. Y. Huang, A. Kou,et al., Structural control of two- level defect density revealed by high-throughput cor- relative measurements of Josephson junctions (2026), 10.48550/arXiv.2602.11469

  19. [19]

    Degnan, C.-C

    Z. Degnan, C.-C. Chiu, Y.-H. Chen, D. Som- mers, L. Abdurakhimov, L. Zhu, A. Fedorov, and P. Jacobson, Reducing TLS loss in tantalum CPW resonators using titanium sacrificial layers (2026), 10.48550/arXiv.2601.16369

  20. [20]

    G. J. Grabovskij, T. Peichl, J. Lisenfeld, G. Weiss, and A. V. Ustinov, Strain tuning of individual atomic tunnel- ing systems detected by a superconducting qubit, Science 338, 232 (2012)

  21. [21]

    Bilmes, A

    A. Bilmes, A. Megrant, P. Klimov, G. Weiss, J. M. Mar- tinis, A. V. Ustinov, and J. Lisenfeld, Resolving the po- sitions of defects in superconducting quantum bits, Sci- entific reports10, 3090 (2020)

  22. [22]

    Y. Kim, L. C. Govia, A. Dane, E. v. d. Berg, D. M. Zajac, B. Mitchell, Y. Liu, K. Balakrishnan, G. Keefe, A. Stabile,et al., Error mitigation with stabilized noise in superconducting quantum processors (2024), 10.48550/arXiv.2407.02467

  23. [23]

    Chen, K.-H

    L. Chen, K.-H. Lee, C.-H. Liu, B. Marinelli, R. K. Naik, Z. Kang, N. Goss, H. Kim, D. I. Santiago, and I. Siddiqi, Scalable and site-specific frequency tuning of two-level system defects in superconducting qubit arrays (2025), 10.48550/arXiv.2503.04702

  24. [24]

    A. Dane, K. Balakrishnan, B. Wacaser, L.-W. Hung, H. J. Mamin, D. Rugar, R. M. Shelby, C. Mur- ray, K. Rodbell, and J. Sleight, Performance stabiliza- tion of high-coherence superconducting qubits (2025), 10.48550/arXiv.2503.12514. 7

  25. [25]

    F. Ye, A. Ellaboudy, and J. M. Nichol, Stabilizing an individual charge fluctuator in a Si/Si-Ge quantum dot, Physical Review Applied23, 044063 (2025)

  26. [26]

    Operating a bistable qubit

    F. Berritta, J. A. Krzywda, T. Dvir, P. Buttles, S. Eil- hart, J. Danon, and F. Kuemmeth, Operating a bistable qubit (2026), 10.48550/arXiv.2605.03187

  27. [27]

    Berritta, J

    F. Berritta, J. Benestad, J. A. Krzywda, O. Krause, M. A. Marciniak, S. Krøjer, C. W. Warren, E. Hogedal, A. Nylander, I. Ahmad, A. Osman, J. Biznárová, M. Rommel, A. F. Roudsari, J. Bylander, G. Tancredi, J. Danon, J. Hastrup, F. Kuemmeth, and M. Kjaer- gaard, Real-time adaptive tracking of fluctuating relax- ation rates in superconducting qubits, Physic...

  28. [28]

    P. V. Klimov, J. Kelly, Z. Chen, M. Neeley, A. Megrant, B. Burkett, R. Barends, K. Arya, B. Chiaro, Y. Chen, A. Dunsworth, A. Fowler, B. Foxen, C. Gidney, M. Giustina, R. Graff,et al., Fluctuations of energy- relaxation times in superconducting qubits, Physical Re- view Letters121, 090502 (2018)

  29. [29]

    Béjanin, C

    J. Béjanin, C. Earnest, A. Sharafeldin, and M. Mariantoni, Interacting defects generate stochastic fluctuations in superconducting qubits, Physical Review B104, 094106 (2021)

  30. [30]

    Carroll, S

    M. Carroll, S. Rosenblatt, P. Jurcevic, I. Lauer, and A. Kandala, Dynamics of superconducting qubit relax- ation times, npj Quantum Information8, 132 (2022)

  31. [31]

    T. Roy, X. You, D. van Zanten, F. Crisa, S. Garat- toni, S. Zhu, A. Grassellino, and A. Romanenko, Two-level system spectroscopy from correlated mul- tilevel relaxation in superconducting qubits (2026), 10.48550/arXiv.2602.11127

  32. [32]

    Leroux, S

    C. Leroux, S. F. Lin, P. Bienias, K. R. Sankar, A. Ben- hemou, A. Kubica, and J. K. Iverson, Snakes and ladders: Adapting the surface code to defects, PRX Quantum6, 040302 (2025)

  33. [33]

    J. M. Drouet, X. C. Kolesnikow, C. K. McLauchlan, G. M. Nixon, S.-H. Lee, D. J. Williamson, S. D. Bartlett, B. J. Brown, and R. Harper, Error correction on an ar- ray of superconducting qubits with defective components (2026), 10.48550/arXiv.2607.12118

  34. [34]

    P. V. Klimov, A. Bengtsson, C. Quintana, A. Bourassa, S. Hong, A. Dunsworth, K. J. Satzinger, W. P. Liv- ingston, V. Sivak, M. Y. Niu,et al., Optimizing quantum gates towards the scale of logical qubits, Nature Commu- nications15, 2442 (2024)

  35. [35]

    Krantz, M

    P. Krantz, M. Kjaergaard, F. Yan, T. P. Orlando, S. Gus- tavsson, and W. D. Oliver, A quantum engineer’s guide to superconducting qubits, Applied Physics Reviews6, 021318 (2019)

  36. [36]

    Blais, A

    A. Blais, A. L. Grimsmo, S. M. Girvin, and A. Wall- raff, Circuit quantum electrodynamics, Reviews of Mod- ern Physics93, 025005 (2021)

  37. [37]

    [53– 57]

    See the Supplemental Material, which includes Ref. [53– 57]

  38. [38]

    Gebhart, R

    V. Gebhart, R. Santagati, A. A. Gentile, E. M. Gauger, D. Craig, N. Ares, L. Banchi, F. Marquardt, L. Pezzè, and C. Bonato, Learning quantum systems, Nature Re- views Physics5, 141 (2023)

  39. [39]

    R. C. Kurchin, Using Bayesian parameter estimation to learn more from data without black boxes, Nature Re- views Physics6, 152 (2024)

  40. [40]

    M. J. Arshad, C. Bekker, B. Haylock, K. Skrzypczak, D. White, B. Griffiths, J. Gore, G. W. Morley, P. Salter, J. Smith, I. Zohar, A. Finkler, Y. Altmann, E. M. Gauger, and C. Bonato, Real-time adaptive estimation of decoherence timescales for a single qubit, Physical Re- view Applied21, 024026 (2024)

  41. [41]

    Berritta, T

    F. Berritta, T. Rasmussen, J. A. Krzywda, J. van der Heijden, F. Fedele, S. Fallahi, G. C. Gardner, M. J. Man- fra, E. van Nieuwenburg, J. Danon, A. Chatterjee, and F. Kuemmeth, Real-time two-axis control of a spin qubit, Nature Communications15, 1676 (2024)

  42. [42]

    Berritta, J

    F. Berritta, J. A. Krzywda, J. Benestad, J. van der Hei- jden, F. Fedele, S. Fallahi, G. C. Gardner, M. J. Man- fra, E. van Nieuwenburg, J. Danon, A. Chatterjee, and F. Kuemmeth, Physics-informed tracking of qubit fluctu- ations, Physical Review Applied22, 014033 (2024)

  43. [43]

    J. Park, H. Jang, H. Sohn, J. Yun, Y. Song, B. Kang, L. E. A. Stehouwer, D. D. Esposti, G. Scappucci, and D. Kim, Passive and active suppression of transduced noise in silicon spin qubits, Nature Communications16, 78 (2025)

  44. [44]

    Berritta, J

    F. Berritta, J. Benestad, L. Pahl, M. Mathews, J. A. Krzywda, R. Assouly, Y. Sung, D. K. Kim, B. M. Niedzielski, K. Serniak, M. E. Schwartz, J. L. Yoder, A. Chatterjee, J. A. Grover, J. Danon, W. D. Oliver, and F. Kuemmeth, Efficient qubit calibration by binary- search Hamiltonian tracking, PRX Quantum6, 030335 (2025)

  45. [45]

    Belliardo, E

    F. Belliardo, E. M. Gauger, M. H. Abobeih, T. H. Taminiau, Y.Altmann,andC.Bonato,Multidimensional quantumestimationandmodellearningframeworkbased on variational Bayesian inference, PRX Quantum7, 020360 (2026)

  46. [46]

    The likelihood function isP(m|Γ 1, τ) = 1−m− (−1)m[β+ (1−α−β)e −Γ1τ ],whereαandβare the misclassification probabilities for measuring|0⟩when the true state at the beginning of the measurement is|1⟩and measuring|1⟩when the true state is|0⟩, respectively

  47. [47]

    The two-parameter posterior approximation keeps the Bayesian update time to about2.2µs

    A probe cycle consists of a1.5µsreadout, an approxi- mately3µsresonator-depletion wait, and the adaptive wait time. The two-parameter posterior approximation keeps the Bayesian update time to about2.2µs

  48. [48]

    Knill, D

    E. Knill, D. Leibfried, R. Reichle, J. Britton, R. B. Blakestad, J. D. Jost, C. Langer, R. Ozeri, S. Seidelin, and D. J. Wineland, Randomized benchmarking of quan- tum gates, Physical Review A77, 012307 (2008)

  49. [50]

    Lisenfeld, A

    J. Lisenfeld, A. K. Händel, E. Daum, B. Berlitz, A. Bilmes, and A. V. Ustinov, Mapping the positions of two-level-systems on the surface of a superconducting transmon qubit (2026), 10.48550/arXiv.2511.05365

  50. [51]

    Thorbeck, A

    T. Thorbeck, A. Eddins, I. Lauer, D. T. McClure, and M. Carroll, Two-level-system dynamics in a supercon- ducting qubit due to background ionizing radiation, PRX Quantum4, 020356 (2023)

  51. [52]

    M. A. Marciniak, R. T. Birke, J. B. Severin, F. Berritta, D. Kjær, F. Nilsson, S. N. Themadath, S. Kallatt, J. L. Webb, K. Bentsen, T. Madsen, Z. Sun, S. Krøjer, C. W. Warren, J. Hastrup, and M. Kjaergaard, Millisecond- scale calibration and benchmarking of superconducting qubits (2026), 10.48550/arXiv.2602.11912. 8

  52. [53]

    M. A. Gingras, B. M. Niedzielski, K. A. Grossklaus, D. Miller, F. Contipelli, K. Azar, L. D. Burkhart, G. Calusine, D. Davis, R. D. Piñero, J. M. Gertler, T. M. Hazard, C. F. Hirjibehedin, D. K. Kim, J. M. Knecht, A. J. Melville, C. O’Connell, R. A. Rood, A. Sabbah, H. Stickler, J. L. Yoder, W. D. Oliver, M. E. Schwartz, and K. Serniak, Improving the perf...

  53. [54]

    S.Frühwirth-Schnatter,Dataaugmentationanddynamic linear models, Journal of Time Series Analysis15, 183 (1994)

  54. [56]

    L. E. Baum, T. Petrie, G. Soules, and N. Weiss, A max- imization technique occurring in the statistical analysis of probabilistic functions of Markov chains, The Annals of Mathematical Statistics41, 164 (1970)

  55. [57]

    Adaptive Spectroscopy of Fast Two-level-system Dynamics in Superconducting Qubits

    L. R. Rabiner, A tutorial on hidden Markov models and selected applications in speech recognition, Proceedings of the IEEE77, 257 (1989). Supplemental Material for “Adaptive Spectroscopy of Fast Two-level-system Dynamics in Superconducting Qubits” Fabrizio Berritta ∗ Research Laboratory of Electronics, Massachusetts Institute of Technology, Cambridge, MA ...

  56. [58]

    Forward-filtering backward-sampling 5

  57. [59]

    Allan deviation 7 III

    Counting uncertainty 6 C. Allan deviation 7 III. Spectral diffusion analysis 9 IV. On- and off-TLS randomized benchmarking comparison 11 V. Wide frequency scans 12 References 12 I. EXPERIMENT AL SETUPS Device A is operated at the Israeli Quantum Computing Center (IQCC), while Device B is operated in an academic laboratory at MIT. Figure S1(a,b) shows the ...

  58. [60]

    Forward-filtering backward-sampling To propagate the state-assignment uncertainty we drawN p = 400 full state paths from the exact posteriorP(s 1:T | x1:T ) by forward-filtering backward-sampling [4, 5]. The forward pass computes the filtered joint likelihood, and the 6 0 200 400 600 laboratory time (s) 325 300 275 250 225 200 frequency detuning (MHz) ban...

  59. [61]

    Counting uncertainty Given a sampled path the two states are fixed, but the rates are still only finitely sampled: we observe a limited number of switches, and this counting noise must be propagated. Conditional on a path, each time step in state 0 independently transitions to state 1 with probabilityp 01, so the number of such transitions among then 0 ti...

  60. [62]

    M. A. Gingras, B. M. Niedzielski, K. A. Grossklaus, D. Miller, F. Contipelli, K. Azar, L. D. Burkhart, G. Calusine, D. Davis, R. D. Pi˜ nero, J. M. Gertler, T. M. Hazard, C. F. Hirjibehedin, D. K. Kim, J. M. Knecht, A. J. Melville, C. O’Connell, R. A. Rood, A. Sabbah, H. Stickler, J. L. Yoder, W. D. Oliver, M. E. Schwartz, and K. Serniak, Improving the pe...

  61. [63]

    L. R. Rabiner, A tutorial on hidden Markov models and selected applications in speech recognition, Proceedings of the IEEE77, 257 (1989)

  62. [64]

    L. E. Baum, T. Petrie, G. Soules, and N. Weiss, A maximization technique occurring in the statistical analysis of probabilistic functions of Markov chains, The Annals of Mathematical Statistics41, 164 (1970)

  63. [65]

    C. K. Carter and R. Kohn, On Gibbs sampling for state space models, Biometrika81, 541 (1994)

  64. [66]

    Fr¨ uhwirth-Schnatter, Data augmentation and dynamic linear models, Journal of Time Series Analysis15, 183 (1994)

    S. Fr¨ uhwirth-Schnatter, Data augmentation and dynamic linear models, Journal of Time Series Analysis15, 183 (1994)

  65. [67]

    P. J. J. O’Malley, J. Kelly, R. Barends, B. Campbell, Y. Chen, Z. Chen, B. Chiaro, A. Dunsworth, A. G. Fowler, I.-C. Hoi, E. Jeffrey, A. Megrant, J. Mutus, C. Neill, C. Quintana,et al., Qubit metrology of ultralow phase noise using randomized benchmarking, Physical Review Applied3, 044009 (2015)