Rust acceleration, quantum vs classical scaling, and IBM hardware validation results.
PyO3 Rust extension with Rayon parallelisation. Benchmarked against pure Qiskit on identical workloads.
| Operation | Qiskit (Python) | Rust Engine | Speedup |
|---|---|---|---|
| Hamiltonian construction (8q) | 2,701 ms | 0.5 ms | 5,401× |
| ICI three-level evolution (2,000 pts) | 68.30 ms | 0.04 ms | 1,665× |
| Pauli string operations | 95 ms | 0.12 ms | 792× |
| Knm matrix utilities | 180 ms | 0.4 ms | 450× |
| DLA dimension computation | 340 ms | 1.2 ms | 283× |
| OTOC evaluation | 520 ms | 2.1 ms | 248× |
| Krylov complexity | 410 ms | 1.8 ms | 228× |
| Lindblad operators | 290 ms | 1.5 ms | 193× |
| Kuramoto order parameter | 45 ms | 0.3 ms | 150× |
| (α,β)-hypergeometric envelope (2F1) | scipy loop | custom series | 44× |
37 Rust functions exported via PyO3 across 21 source files (3,983 LOC). ndarray + num-complex + rayon. Edition 2021. Rust vs Python parity verified at 4.97×10−14 for the ICI integrator.
| Experiment | Qubits | Depth | Hardware | Exact | Error |
|---|---|---|---|---|---|
| VQE ground state | 4 | 12 CZ | −6.2998 | −6.3030 | 0.05% |
| Kuramoto XY (1 rep) | 4 | 85 | R=0.743 | R=0.802 | 7.3% |
| Qubit scaling | 6 | 147 | R=0.482 | R=0.532 | 9.3% |
| UPDE-16 snapshot | 16 | 770 | R=0.332 | R=0.615 | 46% |
IBM ibm_fez (Heron r2, 156 qubits). February 2026. 33 jobs, 176,000+ shots.
| Depth | Reps | Even sector leakage | Odd sector leakage | Asymmetry | Welch p |
|---|---|---|---|---|---|
| 2 | 12 | 8.06% | 8.27% | −2.5% | 0.446 |
| 4 | 21 | 9.82% | 8.62% | +13.98% | 1.45×10−6 |
| 6 | 21 | 12.91% | 10.99% | +17.48% | 6.61×10−6 |
| 8 | 21 | 14.43% | 12.84% | +12.41% | 8.89×10−5 |
| 10 | 21 | 16.58% | 14.95% | +10.91% | 6.67×10−6 |
| 14 | 21 | 18.98% | 17.97% | +5.58% | 0.00995 |
| 20 | 12 | 22.95% | 21.14% | +8.55% | 0.00666 |
| 30 | 12 | 27.71% | 25.76% | +7.58% | 0.00955 |
342 circuits, 4 qubits, 8 Trotter depth points. Fisher's combined χ2 = 123.40 (df = 16), combined p ≪ 10−16. 7 of 8 depths significant at p < 0.05. Readout baseline mean error 1.67%. First hardware confirmation of the SCPN DLA parity prediction: the odd (“feedback”) sector of the XY Hamiltonian dynamical Lie algebra is robustly more decoherence-resistant than the even (“projection”) sector.
Beyond ~250 CZ gates, decoherence dominates. The 16-qubit UPDE experiment (depth 770) shows 46% error — consistent with known T1/T2 limits. Error mitigation (ZNE, PEC, DD) reduces but does not eliminate this. Future fault-tolerant hardware will lift this ceiling.
Classical Kuramoto simulation scales as O(N²) per timestep. Quantum XY evolution on ideal hardware scales as O(N) circuit depth with O(log N) measurement overhead. The quantum advantage crossover is estimated at N ≈ 200–500 oscillators on fault-tolerant hardware with logical error rates < 10−6.
Matrix Product State simulation provides efficient classical baseline for low-entanglement regimes. Included in benchmarks subpackage for honest comparison. MPS wins for weakly-coupled systems; quantum hardware wins for strongly-entangled synchronisation.