← Blueprint · Programs and data of «Certified Speed-Limit Diagnostics for Deployed Superconducting Quantum Processors»

Draft package, local preview. Everything needed to reproduce every number and figure of the paper: the analysis scripts, their outputs, the saved optimal CZ pulses (re-evaluated deterministically by verificar_cz.py and epsilon_cz.py), and the public IBM snapshots they read (Apache 2.0, Qiskit). Our code is under Apache-2.0 (LICENSE); origin of the data in NOTICE.

Download the whole package (zip, 5.2 MB) · 555 files · sha256 e46fa1f156cd883359c2b0e61dc94303a4ae45d8018c3202ef4f28a0cba011a7 · README.md

The blueprint’s own checks (recompute the headline numbers from the archived CSV with controls that must fail): comprueba_numeros.py, prueba_labs.py, qsl_diagnostics_lab.py.

Claim → script → output

PathContent
*.pyanalysis scripts (one per claim, see the map below)
barrido.csv, pares.csv, serie.csv, mejora_por_cubit_salida.txtarchived outputs
datos/qiskit-ibm-runtime 0.51.0 fake-backend snapshots (props_*.json, conf_*.json)
datos_pulsos/qiskit-terra 0.24.1 snapshots with calibrated pulses (defs_*.json, conf_*.json)
datos_serie/snapshots from 13 Qiskit releases, one folder per (backend, date)
paper/figuras.py, paper/figuras_3d.pythe figures of the paper
ficha_ideas.tomlthe hypotheses H1-H3, written before the sweep
LICENSE, NOTICEApache-2.0 for our code; origin and licence of the Qiskit data
cz_best_*.npythe saved optimal CZ pulses (witnesses): re-evaluated deterministically by verificar_cz.py, epsilon_cz.py, poblacion_cz.py
cz_tope_*.txt, cz_tope_salida.txt, anchura_u4d_salida.txt, bell_anchura_salida.txt, cr_continuacion_salida.txtarchived logs of the optimal-control runs
pulso_x_r*.npy, truncacion_L6_salida.txt, rfull_salida.txtsaved single-qubit X pulses and archived outputs of truncacion_1q.py (not re-run when packaging)
*_salida.txt of the quick scriptsregenerated when the package is built (cruces_area_cz, tercera_revision, mediana_por_backend, toronto_sensibilidad, identidad_calibracion, cota_cz_refinada, verificar_cz, epsilon_cz, poblacion_cz)
fuentes/cruce_mg2020.py, fuentes/cruce_mg2020_salida.txthigher-order cross-resonance check (only our script; the third-party papers are not redistributed)
Paper section / claimScriptOutput
Drive convention c = 1: calibrated π pulse gives 1.148π (c = 2 gives 2.296π), 32 backendspulsos_x.pyprinted table
omegad · amp identical to 9 digits for equal σ; no closed form for 1.148identidad_omegad.pyprinted
Two-level controls of the convention (π/omegad, doubling)convencion.pyprinted
Single-qubit gap: 41 backends, 1646 qubits, median 11.30 [7.19, 15.54]; falsifier and exclusionsbarrido.pybarrido.csv
Two-qubit pairs: 1758 pairs, 43 backends; median ratio 4.22 to π/(2 jq)pares.pypares.csv
Cross resonance vs first-order formula: 1212 CX, median θ/(π/2) 1.11; Ω_peak/\|δ\| 0.0907; Ω_CR/Δ 0.451cruce_formulas.pyprinted
Per-qubit durations: ×1.09; cairo above the capmejora_por_cubit.pymejora_por_cubit_salida.txt
Time series 2021–2025 (13 Qiskit releases)serie_temporal.py, serie_porque.pyserie.csv, printed
Lengthening does not lower the error; stratified controls (retraction of a pooled effect)duracion_error.py, degradacion.py, frecuencia_degradados.pyprinted
Certified two-qubit bound: numerical controls (2.6 % from π/(4J); Duhamel; random local control)control_dos_cubits.pyprinted
m_d for CNOT with extra levels (1/√2 for d = 2, 3)cnot_qudit.py, cnot_m3.pyprinted
Zieliński–Zych inequality scanned for α ∈ (0, 10]frontera_zz.py, frontera_zz2.pyprinted
Optimal control: single qubit law T_min ≈ max(π/Ω, ~1.1·2π/\|δ\|)grape_1q.pyprinted
Optimal control: cross resonance, qutrit-assisted CZ, Bell with ancillas (scipy versions)grape_cr.py, grape_cz.py, grape_bell.pyprinted
Convention from calibrated π pulses (θ/π range over pulse shapes, K(σ))identidad_calibracion.pyidentidad_calibracion_salida.txt
Backend-weighted median 11.37, clustered-bootstrap interval [11.02, 11.65]mediana_por_backend.pymediana_por_backend_salida.txt
Toronto 2021 sensitivity (2.78 with, 2.74 without)toronto_sensibilidad.pytoronto_sensibilidad_salida.txt
CR with 0.7 drive fraction: monotonicity of the feasible timegrape_cr_continuacion.pycr_continuacion_salida.txt
Qutrit-assisted CZ, amplitude-cap scan (PyTorch; GPU or CPU, long)grape_cz_torch.pycz_tope_salida.txt, cz_tope_*.txt, cz_best_*.npy
Independent re-evaluation of every saved CZ pulse (infidelity, leakage, phases, \|\|A_k\|\|)verificar_cz.pyverificar_cz_salida.txt
Robust epsilon for the entangling-time bound: eps(\|++>) and (π/4 − arcsin √eps)/(2J)epsilon_cz.pyepsilon_cz_salida.txt
Population outside the qubit subspace during the CZ (up to 60 %)poblacion_cz.pypoblacion_cz_salida.txt
Bell-state preparation with ancillas, width scangrape_bell_torch.py, grape_anchura.pybell_anchura_salida.txt, anchura_u4d_salida.txt
Pulse-area identity omegad·\|amp\|·A_rec = π on 777 qubits (max deviation 1.18e-7); residual 1/(1−e^−2) = 1.15652cruces_area_cz.pycruces_area_cz_salida.txt (regenerated)
CZ construction in π/(6J) (twirl by W, code in the extreme eigenvectors of Q; code state 5/9 outside the qubit subspace); d-dimensional bounds tablecruces_area_cz.pycruces_area_cz_salida.txt
Refined CZ bound T_* ≥ 0.4511/J (lemma Tr[(P_A⊗P_B)H_c²]/4 ≤ 97/32 J², equality case, resulting time bound, T_*(2e-8) ≥ 0.4510/J, the e ≥ 1.9e-6 exclusion at T = 0.45/J); the Lean lemma is QSL.Exchange.refined_cz_lemma (the time argument is not formalized)cota_cz_refinada.pycota_cz_refinada_salida.txt (regenerated)
Exact Bell preparation by \|12⟩↔\|21⟩ (T = 0.3926/J); Choi-based gate bound 0.392599/J for e = 2e-8; X bound in the trace metric (6 % margin at r ≤ 0.2; IBM cap 2.23–2.37×); r_full = omegad/\|δ\| median 0.351tercera_revision.pytercera_revision_salida.txt (regenerated)
Truncation of the fast X pulse: 4-level pulse at r = 1 degrades with ≥ 5 levels, 6-level pulse converged to 12 levels (long optimization; archived)truncacion_1q.pytruncacion_L6_salida.txt, rfull_salida.txt, pulso_x_r*.npy (saved pulses)
Re-evaluation of the saved X pulses with 4, 5, 6, 8, 12 levels and the Duhamel certificate η₄ (4-level pulse at r = 1: 7e-9 in 4 levels, ≈ 7e-4 with more; η₄ = 0.551)truncacion_reevaluar.pytruncacion_reevaluar_salida.txt (regenerated)
Higher-order (Magesan–Gambetta 2020) cross-resonance checkfuentes/cruce_mg2020.pyfuentes/cruce_mg2020_salida.txt
ClaimScriptData / output
Published values of the four links (each with its table of arXiv:2610.29986)—relevo/datos/ibm_2610.29986.toml
Ratios T_g/T_min = 5.8–7.9 (single-mode optimum π/(√2 g), g/2π = 8–9 MHz)relevo/cocientes.pyprinted
IBM's coherence model reproduced (0.625 % against 0.63 % on Q240–Q247)relevo/diag_relevo.py modelo datos/ibm_2610.29986.toml Q240–Q247printed
Coherence error of the optimal protocol 0.13–0.26 % (master equation)relevo/lindblad_relevo.py, relevo/suelo_lindblad.py datos/ibm_2610.29986.tomlprinted
Independent check of the master equation (RK4, other fidelity estimator)relevo/lindblad_rk4.pyprinted (0.1289 %)
Additive model exact for relaxation, not for dephasingrelevo/comparar_modelos.pyprinted
A slower protocol can lose less (not a floor); first-order dephasing formularelevo/falsar_v198.pyprinted
Numerical minimal time π/√2 (real and complex controls)relevo/tiempo_minimo.pyprinted
Validation against the published Lindblad simulation of Heya et al. (2e-3, 6.2e-3)relevo/heya_satd.py cpw, relevo/heya_satd.py todoprinted (2.042e-3, 6.037e-3)
The whole chain with controls and replica (excovery.lanzadera missions)relevo/mision_heya_satd.toml, relevo/mision_relevo_ibm.toml—
Minimal-time lower bound gT ≥ π/√2 and the other relay theorems, machine-checkedisabelle/Quantum_Speed_Limits/ (Isabelle2026 + HOL-Analysis; isabelle build -D .)—

All files

Top level: anchura_u4d_salida.txt · barrido.csv · barrido.py · bell_anchura_salida.txt · cnot_m3.py · cnot_qudit.py · control_dos_cubits.py · convencion.py · cota_cz_refinada.py · cota_cz_refinada_salida.txt · cr_continuacion_salida.txt · cruce_formulas.py · cruces_area_cz.py · cruces_area_cz_salida.txt · cz_best_d2_M8_T1.4000_aNone.npy · cz_best_d2_M8_T1.5710_aNone.npy · cz_best_d2_M8_T1.6500_aNone.npy · cz_best_d3_M16_T0.4800_aNone.npy · cz_best_d3_M16_T0.5000_aNone.npy · cz_best_d3_M16_T0.5100_aNone.npy · cz_best_d3_M16_T0.5150_aNone.npy · cz_best_d3_M16_T0.5200_aNone.npy · cz_best_d3_M16_T0.5236_aNone.npy · cz_best_d3_M32_T0.5200_aNone.npy · cz_best_d3_M8_T0.4000_aNone.npy · cz_best_d3_M8_T0.4500_a20.0.npy · cz_best_d3_M8_T0.4500_a40.0.npy · cz_best_d3_M8_T0.4500_aNone.npy · cz_best_d3_M8_T0.4800_aNone.npy · cz_best_d3_M8_T0.5000_a10.0.npy · cz_best_d3_M8_T0.5000_a20.0.npy · cz_best_d3_M8_T0.5000_a40.0.npy · cz_best_d3_M8_T0.5000_aNone.npy · cz_best_d3_M8_T0.5100_aNone.npy · cz_best_d3_M8_T0.5150_aNone.npy · cz_best_d3_M8_T0.5200_aNone.npy · cz_best_d3_M8_T0.5236_aNone.npy · cz_best_d3_M8_T0.5500_a10.0.npy · cz_best_d3_M8_T0.5500_a20.0.npy · cz_best_d3_M8_T0.5500_a40.0.npy · cz_best_d3_M8_T0.5500_aNone.npy · cz_best_d3_M8_T0.6000_a10.0.npy · cz_best_d3_M8_T0.6000_a20.0.npy · cz_best_d3_M8_T0.6000_a40.0.npy · cz_best_d3_M8_T0.6000_aNone.npy · cz_tope_10_0.55.txt · cz_tope_20_0.45.txt · cz_tope_20_0.5.txt · cz_tope_20_0.55.txt · cz_tope_20_0.6.txt · cz_tope_40_0.45.txt · cz_tope_40_0.5.txt · cz_tope_40_0.55.txt · cz_tope_40_0.6.txt · cz_tope_salida.txt · cz_v3_certificado_salida.txt · cz_v3_norma_K.py · cz_v3_norma_K_salida.txt · cz_v3_revival.py · cz_v3_revival_salida.txt · degradacion.py · duracion_error.py · epsilon_cz.py · epsilon_cz_salida.txt · ficha_ideas.toml · frecuencia_degradados.py · frontera_zz.py · frontera_zz2.py · grape_1q.py · grape_anchura.py · grape_bell.py · grape_bell_torch.py · grape_cr.py · grape_cr_continuacion.py · grape_cz.py · grape_cz_torch.py · identidad_calibracion.py · identidad_calibracion_salida.txt · identidad_omegad.py · LICENSE · mediana_por_backend.py · mediana_por_backend_salida.txt · mejora_por_cubit.py · mejora_por_cubit_salida.txt · NOTICE · pares.csv · pares.py · poblacion_cz.py · poblacion_cz_salida.txt · pulso_x_r0.1_T33.0_L4.npy · pulso_x_r0.351_T10.0_L6.npy · pulso_x_r0.351_T10.8_L6.npy · pulso_x_r0.351_T11.6_L6.npy · pulso_x_r0.351_T9.4_L6.npy · pulso_x_r1.0_T10.0_L6.npy · pulso_x_r1.0_T12.0_L6.npy · pulso_x_r1.0_T7.0_L4.npy · pulso_x_r1.0_T7.0_L6.npy · pulso_x_r1.0_T8.2_L6.npy · pulsos_x.py · README.md · rfull_salida.txt · serie.csv · serie_porque.py · serie_temporal.py · tercera_revision.py · tercera_revision_salida.txt · toronto_sensibilidad.py · toronto_sensibilidad_salida.txt · truncacion_1q.py · truncacion_L6_salida.txt · truncacion_reevaluar.py · truncacion_reevaluar_salida.txt · verificar_cz.py · verificar_cz_salida.txt

datos/ (137 files)
datos_pulsos/ (67 files)
datos_serie/ (203 files)
fuentes/ (2 files)
isabelle/ (15 files)
paper/ (2 files)
relevo/ (14 files)