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.
| Path | Content |
|---|---|
*.py | analysis scripts (one per claim, see the map below) |
barrido.csv, pares.csv, serie.csv, mejora_por_cubit_salida.txt | archived 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.py | the figures of the paper |
ficha_ideas.toml | the hypotheses H1-H3, written before the sweep |
LICENSE, NOTICE | Apache-2.0 for our code; origin and licence of the Qiskit data |
cz_best_*.npy | the 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.txt | archived logs of the optimal-control runs |
pulso_x_r*.npy, truncacion_L6_salida.txt, rfull_salida.txt | saved single-qubit X pulses and archived outputs of truncacion_1q.py (not re-run when packaging) |
*_salida.txt of the quick scripts | regenerated 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.txt | higher-order cross-resonance check (only our script; the third-party papers are not redistributed) |
| Paper section / claim | Script | Output |
|---|---|---|
| Drive convention c = 1: calibrated π pulse gives 1.148π (c = 2 gives 2.296π), 32 backends | pulsos_x.py | printed table |
| omegad · amp identical to 9 digits for equal σ; no closed form for 1.148 | identidad_omegad.py | printed |
| Two-level controls of the convention (π/omegad, doubling) | convencion.py | printed |
| Single-qubit gap: 41 backends, 1646 qubits, median 11.30 [7.19, 15.54]; falsifier and exclusions | barrido.py | barrido.csv |
| Two-qubit pairs: 1758 pairs, 43 backends; median ratio 4.22 to π/(2 jq) | pares.py | pares.csv |
| Cross resonance vs first-order formula: 1212 CX, median θ/(π/2) 1.11; Ω_peak/\|δ\| 0.0907; Ω_CR/Δ 0.451 | cruce_formulas.py | printed |
Per-qubit durations: ×1.09; cairo above the cap | mejora_por_cubit.py | mejora_por_cubit_salida.txt |
| Time series 2021–2025 (13 Qiskit releases) | serie_temporal.py, serie_porque.py | serie.csv, printed |
| Lengthening does not lower the error; stratified controls (retraction of a pooled effect) | duracion_error.py, degradacion.py, frecuencia_degradados.py | printed |
| Certified two-qubit bound: numerical controls (2.6 % from π/(4J); Duhamel; random local control) | control_dos_cubits.py | printed |
| m_d for CNOT with extra levels (1/√2 for d = 2, 3) | cnot_qudit.py, cnot_m3.py | printed |
| Zieliński–Zych inequality scanned for α ∈ (0, 10] | frontera_zz.py, frontera_zz2.py | printed |
| Optimal control: single qubit law T_min ≈ max(π/Ω, ~1.1·2π/\|δ\|) | grape_1q.py | printed |
| Optimal control: cross resonance, qutrit-assisted CZ, Bell with ancillas (scipy versions) | grape_cr.py, grape_cz.py, grape_bell.py | printed |
| Convention from calibrated π pulses (θ/π range over pulse shapes, K(σ)) | identidad_calibracion.py | identidad_calibracion_salida.txt |
| Backend-weighted median 11.37, clustered-bootstrap interval [11.02, 11.65] | mediana_por_backend.py | mediana_por_backend_salida.txt |
| Toronto 2021 sensitivity (2.78 with, 2.74 without) | toronto_sensibilidad.py | toronto_sensibilidad_salida.txt |
| CR with 0.7 drive fraction: monotonicity of the feasible time | grape_cr_continuacion.py | cr_continuacion_salida.txt |
| Qutrit-assisted CZ, amplitude-cap scan (PyTorch; GPU or CPU, long) | grape_cz_torch.py | cz_tope_salida.txt, cz_tope_*.txt, cz_best_*.npy |
| Independent re-evaluation of every saved CZ pulse (infidelity, leakage, phases, \|\|A_k\|\|) | verificar_cz.py | verificar_cz_salida.txt |
| Robust epsilon for the entangling-time bound: eps(\|++>) and (π/4 − arcsin √eps)/(2J) | epsilon_cz.py | epsilon_cz_salida.txt |
| Population outside the qubit subspace during the CZ (up to 60 %) | poblacion_cz.py | poblacion_cz_salida.txt |
| Bell-state preparation with ancillas, width scan | grape_bell_torch.py, grape_anchura.py | bell_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.15652 | cruces_area_cz.py | cruces_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 table | cruces_area_cz.py | cruces_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.py | cota_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.351 | tercera_revision.py | tercera_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.py | truncacion_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.py | truncacion_reevaluar_salida.txt (regenerated) |
| Higher-order (Magesan–Gambetta 2020) cross-resonance check | fuentes/cruce_mg2020.py | fuentes/cruce_mg2020_salida.txt |
| Claim | Script | Data / 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.py | printed |
| IBM's coherence model reproduced (0.625 % against 0.63 % on Q240–Q247) | relevo/diag_relevo.py modelo datos/ibm_2610.29986.toml Q240–Q247 | printed |
| Coherence error of the optimal protocol 0.13–0.26 % (master equation) | relevo/lindblad_relevo.py, relevo/suelo_lindblad.py datos/ibm_2610.29986.toml | printed |
| Independent check of the master equation (RK4, other fidelity estimator) | relevo/lindblad_rk4.py | printed (0.1289 %) |
| Additive model exact for relaxation, not for dephasing | relevo/comparar_modelos.py | printed |
| A slower protocol can lose less (not a floor); first-order dephasing formula | relevo/falsar_v198.py | printed |
| Numerical minimal time π/√2 (real and complex controls) | relevo/tiempo_minimo.py | printed |
| Validation against the published Lindblad simulation of Heya et al. (2e-3, 6.2e-3) | relevo/heya_satd.py cpw, relevo/heya_satd.py todo | printed (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-checked | isabelle/Quantum_Speed_Limits/ (Isabelle2026 + HOL-Analysis; isabelle build -D .) | — |
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