pith. machine review for the scientific record. sign in

arxiv: 2604.21458 · v1 · submitted 2026-04-23 · 🪐 quant-ph

Recognition: unknown

HEOM-in-Calibration-Loop: Exposing Non-Markovian Bath Signatures That Markovian Calibration Elides in Superconducting-Qubit Tune-Up

Jun Ye
Authors on Pith no claims yet

Pith reviewed 2026-05-09 22:36 UTC · model grok-4.3

classification 🪐 quant-ph
keywords superconducting qubitsnon-Markovian dynamicsHEOMqubit calibrationRamsey decay1/f noisecoherence time T2*
0
0 comments X

The pith

Embedding a non-Markovian HEOM solver in superconducting-qubit calibration recovers physical revival envelopes that Markovian fits suppress, producing T2* estimates 13 to 72 times longer.

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

The paper integrates a hierarchical-equations-of-motion solver driven by a 1/f noise bath directly into a closed-loop calibration protocol chain for superconducting qubits and compares it against standard Markovian solvers in simulation. In the Ramsey channel the non-Markovian treatment yields a coherence revival whose extracted primary T2* differs from the Markovian reference by at least 13 times at 95 percent bootstrap , with point estimates reaching 28 times and 72 times on denser grids. Rabi contrast shows a small systematic drop while T1 shape remains similar except for a reduced initial occupation, indicating that bath structure appears as an explicit residual rather than being absorbed into fit errors. The added scheduling overhead stays under 10 microseconds per protocol, showing the approach fits inside existing calibration budgets.

Core claim

The Markovian fit is censored by its exponential-family numerical ceiling, while HEOM recovers a physical revival envelope whose primary T2* separates from the Markovian reference by at least 13x at 95% independent-bootstrap within the HEOM-feasible budget; the point-estimate ratio reaches >=28x on the 50-point primary-t1 grid and ~72x on the 30-point biexp-family tau_aw pivot at L=5.

What carries the argument

A Tier-1 1/f Burkard bath inside a QuTiP hierarchical-equations-of-motion solver placed inside a multi-protocol calibration DAG (Rabi followed by Ramsey or T1) that runs against pulse-level simulation backends.

If this is right

  • Calibration output can now treat bath structure as a quantifiable diagnostic rather than a hidden residual.
  • The Ramsey channel provides the strongest statistical separation while Rabi and T1 channels supply corroborative or null results.
  • The HEOM-in-loop DAG adds only 9.62 microseconds of average scheduling overhead with no meaningful latency penalty at protocol scale.
  • T1 decay shape remains identical across backends but initial occupation drops to 0.879 under the non-Markovian model and stays stable under grid densification.

Where Pith is reading between the lines

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

  • Hardware validation on real devices would test whether the simulated revival directly improves coherence estimates used in error budgeting.
  • The same loop could be applied to other noise spectra or qubit modalities to expose platform-specific non-Markovian features.
  • Explicit reporting of bath-induced residuals might allow calibration routines to flag when Markovian assumptions break down before they affect gate fidelity.

Load-bearing premise

The Tier-1 1/f Burkard bath model and the pulse-level simulator faithfully represent the physics of real superconducting qubits so that differences between solvers indicate genuine non-Markovian signatures.

What would settle it

Running the identical Ramsey protocol on physical superconducting qubits and checking whether the measured coherence envelope exhibits the predicted revival peak and T2* ratio that the HEOM simulation produces but the Markovian simulation does not.

Figures

Figures reproduced from arXiv: 2604.21458 by Jun Ye.

Figure 1
Figure 1. Figure 1: DAG Gantt of the RABI→ {RAMSEY ∥ T1} chain. Parallel execution saves 36.9 s (43 %) over serial; average scheduling latency is 9.62 µs. Ryzen 9 7900X, 12 cores / 24 threads, 30 GiB RAM). On this machine, the headline three-backend pipeline completes in under 10 minutes of wall time.3 III. RESULTS Running the same calibration DAG once per backend yields an asymmetric three-channel fingerprint across the 3×3 … view at source ↗
Figure 2
Figure 2. Figure 2: Ramsey coherent-channel gap. (a) Three-backend 30-point com [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
read the original abstract

Closed-loop superconducting-qubit calibration has matured into DAG-orchestrated protocol chains, yet published frameworks treat the bath via a Markovian master equation or a phenomenological likelihood, absorbing bath structure into fit residuals instead of reporting it as a diagnostic. We integrate a QuTiP 5.x hierarchical-equations-of-motion (HEOM) solver driven by a Tier-1 1/f Burkard bath into a multi-protocol calibration DAG (Rabi -> {Ramsey || T1}) and benchmark it against sesolve and mesolve on a frozen platform in a pulse-level simulator (no hardware validation). The Ramsey channel carries the headline: the Markovian fit is censored by its exponential-family numerical ceiling, while HEOM recovers a physical revival envelope whose primary T2* separates from the Markovian reference by at least 13x at 95% independent-bootstrap confidence within the HEOM-feasible budget; the point-estimate ratio reaches >=28x on the 50-point primary-t1 grid and ~72x on the 30-point biexp-family tau_aw pivot at L=5. Rabi contrast falls 2.17% below mesolve on a noise-limited 30-point grid; the paired-bootstrap CI crosses zero, so this channel corroborates rather than independently establishes the non-Markovian signature. T1 decay shape matches across backends (beta=1.000), yet HEOM's initial occupation drops from 1.000 to 0.879 -- a bath-dressed contamination stable under a 16-point densification. The DAG adds 9.62 us average per-protocol scheduling overhead, no meaningful latency penalty at protocol granularity. HEOM-in-loop thereby changes what calibration reports: bath structure appears as a quantifiable residual rather than a hidden confound.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript integrates a QuTiP 5.x HEOM solver driven by a Tier-1 1/f Burkard bath into a multi-protocol calibration DAG (Rabi to Ramsey/T1) and benchmarks it against sesolve and mesolve in a pulse-level simulator with no hardware validation. It claims that the Markovian fit is limited by its exponential-family ceiling while HEOM recovers a physical revival envelope, with the primary T2* separating from the Markovian reference by at least 13x at 95% independent-bootstrap CI (point estimates >=28x on the 50-point primary-t1 grid and ~72x on the 30-point biexp-family tau_aw pivot at L=5); Rabi contrast is 2.17% lower but CI crosses zero, T1 shape matches but initial occupation drops to 0.879, and DAG overhead is negligible.

Significance. If the simulator and bath model faithfully capture real-device physics, the work would show that non-Markovian effects produce measurable, calibration-relevant signatures (especially T2* revival envelopes) that standard Markovian protocols absorb into residuals. The bootstrap CIs on the simulated differences and the explicit comparison of solvers on identical bath realizations are strengths that provide reproducible numerical evidence. However, the complete absence of hardware validation means the transfer of these T2* separations and envelopes to actual superconducting-qubit tune-up remains untested, limiting significance to a methodological demonstration within simulation.

major comments (2)
  1. Abstract: the central claim that HEOM 'exposes real non-Markovian bath signatures' and thereby 'changes what calibration reports' in physical tune-up rests on simulation differences alone; the manuscript explicitly states no hardware validation was performed, so any mismatch between the Burkard spectral density or QuTiP HEOM implementation and real-device channels (TLS, quasiparticles, readout) would render the reported 13x–72x T2* separations non-diagnostic for actual calibration loops.
  2. Ramsey-channel results (50-point and 30-point grids, L=5): the headline T2* separation ratios depend on post-hoc grid density and biexp-family pivot choices; the manuscript should report sensitivity of the >=13x 95% CI and point estimates to variations in these analysis parameters to establish that the separation is not an artifact of the chosen discretization.
minor comments (2)
  1. Abstract: the phrase 'frozen platform' is undefined and should be clarified in the context of the pulse-level simulator.
  2. T1 results: the reported drop in initial occupation to 0.879 under HEOM is presented as bath-dressed contamination; a brief explanation of how this quantity is extracted from the HEOM hierarchy would improve reproducibility.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The central claim rests on the accuracy of the chosen 1/f bath model and the simulator's fidelity; these are domain assumptions rather than derived results.

free parameters (1)
  • Burkard bath parameters (alpha, cutoff, etc.)
    Specific values for the Tier-1 1/f bath are required to run HEOM but are not stated as derived from first principles or external data.
axioms (2)
  • domain assumption The 1/f Burkard bath model is an appropriate representation of the noise environment in superconducting qubits
    Invoked when driving the HEOM solver inside the calibration loop.
  • domain assumption The pulse-level simulator faithfully reproduces the relevant open-system dynamics without hardware-specific discrepancies
    Required for the claim that simulation differences indicate real non-Markovian signatures.

pith-pipeline@v0.9.0 · 5633 in / 1610 out tokens · 28230 ms · 2026-05-09T22:36:59.171437+00:00 · methodology

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

27 extracted references · 7 canonical work pages

  1. [1]

    Physical qubit calibration on a directed acyclic graph,

    J. Kelly, P. O’Malley, M. Neeley, H. Neven, and J. M. Martinis, “Physical qubit calibration on a directed acyclic graph,” 2018, arXiv:1803.03226 [quant-ph]

    work page arXiv 2018

  2. [2]

    Qibocal: an open- source framework for calibration of self-hosted quantum devices,

    A. Pasquale, E. Pedicillo, J. Cereijo, S. Ramos-Calderer, A. Candido, G. Palazzo, R. Carobene, M. Gobbo, S. Efthymiou, Y . P. Tan, I. Roth, M. Robbiati, J. Wilkens, A. Orgaz-Fuertes, D. Fuentes-Ruiz, A. Gi- achero, F. Brito, J. I. Latorre, and S. Carrazza, “Qibocal: an open- source framework for calibration of self-hosted quantum devices,” 2024, arXiv:241...

    work page arXiv 2024

  3. [3]

    Qiskit Experiments: A Python package to characterize and calibrate quantum computers,

    N. Kanazawa, D. J. Egger, Y . Ben-Haim, H. Zhang, W. E. Shanks, G. Aleksandrowicz, and C. J. Wood, “Qiskit Experiments: A Python package to characterize and calibrate quantum computers,”Journal of Open Source Software, vol. 8, no. 84, p. 5329, 2023

    2023

  4. [4]

    Leakage reduction in fast superconducting qubit gates via optimal control,

    M. Werninghaus, D. J. Egger, F. Roy, S. Machnes, F. K. Wilhelm, and S. Filipp, “Leakage reduction in fast superconducting qubit gates via optimal control,”npj Quantum Information, vol. 7, p. 14, 2021

    2021

  5. [5]

    Integrated tool set for control, calibration, and characterization of quantum devices applied to superconducting qubits,

    N. Wittler, F. Roy, K. Pack, M. Werninghaus, A. S. Roy, D. J. Egger, S. Filipp, F. K. Wilhelm, and S. Machnes, “Integrated tool set for control, calibration, and characterization of quantum devices applied to superconducting qubits,”Physical Review Applied, vol. 15, p. 034080, 2021

    2021

  6. [6]

    Non-Markovian qubit dynamics in the presence of1/f noise,

    G. Burkard, “Non-Markovian qubit dynamics in the presence of1/f noise,”Physical Review B, vol. 79, p. 125317, 2009

    2009

  7. [7]

    Noise spectroscopy through dynamical decoupling with a superconducting flux qubit,

    J. Bylander, S. Gustavsson, F. Yan, F. Yoshihara, K. Harrabi, G. Fitch, D. G. Cory, Y . Nakamura, J.-S. Tsai, and W. D. Oliver, “Noise spectroscopy through dynamical decoupling with a superconducting flux qubit,”Nature Physics, vol. 7, pp. 565–570, 2011

    2011

  8. [8]

    1/fnoise: Implications for solid-state quantum information,

    E. Paladino, Y . M. Galperin, G. Falci, and B. L. Altshuler, “1/fnoise: Implications for solid-state quantum information,”Reviews of Modern Physics, vol. 86, pp. 361–418, 2014

    2014

  9. [9]

    Simulation of one and two qubit superconducting quantum gates under the non-Markovian1/fnoise,

    Y . Chen, S. Zhang, and Q. Shi, “Simulation of one and two qubit superconducting quantum gates under the non-Markovian1/fnoise,” Physical Review B, vol. 113, p. 094302, 2026

    2026

  10. [10]

    Gate operations for superconducting qubits and non-Markovianity,

    K. Nakamura and J. Ankerhold, “Gate operations for superconducting qubits and non-Markovianity,”Physical Review Research, vol. 6, p. 033215, 2024

    2024

  11. [11]

    Impact of time-retarded noise on dynamical decoupling schemes for qubits,

    K. Nakamura and J. Ankerhold, “Impact of time-retarded noise on dynamical decoupling schemes for qubits,”Physical Review B, vol. 111, no. 6, p. 064503, Feb. 2025

    2025

  12. [12]

    Non-perturbative cpmg scaling and qutrit-driven breakdown under compiled superconducting-qubit control: a single-qubit study,

    J. Ye, “Non-perturbative cpmg scaling and qutrit-driven breakdown under compiled superconducting-qubit control: a single-qubit study,” 2026, arXiv:2603.29525 [quant-ph]

    work page arXiv 2026

  13. [13]

    Non-Markovian relaxation spectroscopy of fluxonium qubits,

    Z.-T. Zhuang, D. Rosenstock, B.-J. Liu, A. Somoroff, V . E. Manucharyan, and C. Wang, “Non-Markovian relaxation spectroscopy of fluxonium qubits,”Nature Communications, vol. 17, p. 1212, 2026

    2026

  14. [14]

    Fast quantum calibration using Bayesian optimization with state parameter estimator for non-Markovian environment,

    P. Qian, S. Qamar, X. Xiao, Y . Gu, X. Chai, Z. Zhao, N. Forcellini, and D. E. Liu, “Fast quantum calibration using Bayesian optimization with state parameter estimator for non-Markovian environment,” 2022, arXiv:2205.12929 [quant-ph]

    work page arXiv 2022

  15. [15]

    Calibration of quantum devices via robust statistical methods,

    A. Ramôa, R. Santagati, and N. Wiebe, “Calibration of quantum devices via robust statistical methods,” 2025, arXiv:2507.06941 [quant-ph]

    work page arXiv 2025

  16. [16]

    Efficient qubit calibration by binary-search Hamiltonian tracking,

    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,” 2025, arXiv:2501.05386 [quant-ph]

    work page arXiv 2025

  17. [17]

    Experimental graybox quantum system identification and control,

    A. Youssry, Y . Yang, R. J. Chapman, B. Haylock, F. Lenzini, M. Lobino, and A. Peruzzo, “Experimental graybox quantum system identification and control,”npj Quantum Information, vol. 10, p. 9, 2024

    2024

  18. [18]

    QuTiP 5: The quantum toolbox in Python,

    N. Lambert, E. Giguère, P. Menczel, B. Li, P. Hopf, G. Suárez, M. Gali, J. Lishman, R. Gadhvi, R. Agarwal, A. Galicia, N. Shammah, P. Nation, J. R. Johansson, S. Ahmed, S. Cross, A. Pitchford, and F. Nori, “QuTiP 5: The quantum toolbox in Python,”Physics Reports, vol. 1153, pp. 1–62, 2026, hEOM solver documented in Sec. V

    2026

  19. [19]

    HierarchicalEOM.jl: An efficient Julia frame- work for hierarchical equations of motion in open quantum systems,

    Y .-T. Huang, P.-C. Kuo, N. Lambert, M. Cirio, S. Cross, S.-L. Yang, F. Nori, and Y .-N. Chen, “HierarchicalEOM.jl: An efficient Julia frame- work for hierarchical equations of motion in open quantum systems,” Communications Physics, vol. 6, p. 313, 2023

    2023

  20. [20]

    Efficient non-Markovian quantum dynamics using time-evolving matrix product operators,

    A. Strathearn, P. Kirton, D. Kilda, J. Keeling, and B. W. Lovett, “Efficient non-Markovian quantum dynamics using time-evolving matrix product operators,”Nature Communications, vol. 9, p. 3322, 2018

    2018

  21. [21]

    Markovian embeddings of non-Markovian open system dynamics,

    M. Xu, J. T. Stockburger, and J. Ankerhold, “Markovian embeddings of non-Markovian open system dynamics,” 2026, arXiv:2602.21430 [quant- ph]

    work page arXiv 2026

  22. [22]

    Universal graph-based scheduling for quantum systems,

    L. Riesebos, B. Bondurant, and K. R. Brown, “Universal graph-based scheduling for quantum systems,”IEEE Micro, vol. 41, no. 5, pp. 57–65, 2021

    2021

  23. [23]

    How to enhance dephasing time in superconducting qubits,

    Ł. Cywi ´nski, R. M. Lutchyn, C. P. Nave, and S. Das Sarma, “How to enhance dephasing time in superconducting qubits,”Physical Review B, vol. 77, p. 174509, 2008

    2008

  24. [24]

    Noise-specific beating in the higher-level Ramsey curves of a transmon qubit,

    L. A. Martinez, Z. Peng, D. Appelö, D. M. Tennant, N. A. Petersson, J. L. DuBois, and Y . J. Rosen, “Noise-specific beating in the higher-level Ramsey curves of a transmon qubit,”Applied Physics Letters, vol. 122, no. 11, p. 114002, 2023

    2023

  25. [25]

    Environmental dynamics, correlations, and the emergence of noncanonical equilibrium states in open quantum systems,

    J. Iles-Smith, N. Lambert, and A. Nazir, “Environmental dynamics, correlations, and the emergence of noncanonical equilibrium states in open quantum systems,”Physical Review A, vol. 90, p. 032114, 2014

    2014

  26. [26]

    Strong coupling effects in quantum thermal transport with the reaction coordinate method,

    N. Anto-Sztrikacs and D. Segal, “Strong coupling effects in quantum thermal transport with the reaction coordinate method,”New Journal of Physics, vol. 23, p. 063036, 2021

    2021

  27. [27]

    Nonequilibrium thermodynamics in the strong coupling and non-Markovian regime based on a reaction coordinate mapping,

    P. Strasberg, G. Schaller, N. Lambert, and T. Brandes, “Nonequilibrium thermodynamics in the strong coupling and non-Markovian regime based on a reaction coordinate mapping,”New Journal of Physics, vol. 18, p. 073007, 2016. 7

    2016