When you tile a quantum processor with 6-qubit cells and evolve them under a Hamiltonian that lives on the hypercube Q₆, something useful happens: three ZZ correlator measurements inside each cell are forced to be exactly equal by the Hamiltonian's symmetry. On noisy hardware, they disagree — and the disagreement tells you which cells are corrupted.
This gives you a spatially resolved error map at zero qubit overhead. No surface codes, no additional circuits, no classical post-processing beyond comparing three numbers.
We validated this across 230+ jobs on IBM ibm_fez (48–156 qubits), spanning four unrelated application domains. Here's the mathematical structure behind it, and what we found.
Q₆: The Arena
The 6-dimensional hypercube Q₆ has 64 vertices and 192 edges. Each vertex is a 6-bit string; two vertices are adjacent if they differ in exactly one bit.
(* Construct Q₆ *)
Q6 = HypercubeGraph[6];
(* Basic properties *)
{VertexCount[Q6], EdgeCount[Q6]}
(* {64, 192} *)
(* Visualize — the 6D hypercube projected to 2D *)
Graph[Q6,
GraphLayout -> "SpringElectricalEmbedding",
VertexSize -> Small,
VertexStyle -> Directive[Opacity[0.6], RGBColor[0.2, 0.4, 0.8]],
EdgeStyle -> Directive[Opacity[0.15], Gray],
ImageSize -> 500,
PlotLabel -> "Q₆: 64 vertices, 192 edges"]
Each vertex naturally splits into two 3-bit halves — the "upper" and "lower" trigrams:
(* Every vertex s ∈ {0,...,63} has upper trigram (bits 3-5) and lower trigram (bits 0-2) *)
upperTrigram[s_] := BitShiftRight[s, 3]
lowerTrigram[s_] := BitAnd[s, 7]
(* The three natural qubit pairings within a 6-qubit cell: (0,3), (1,4), (2,5) *)
(* These pair bit k of the lower trigram with bit k of the upper trigram *)
gaugePairs = {{0, 3}, {1, 4}, {2, 5}};
This trigram decomposition is the origin of the gauge symmetry.
The Gauge Symmetry
Consider a Hamiltonian on 6 qubits whose diagonal potential depends only on the relationship between the upper and lower trigrams — not on either trigram individually. Such a Hamiltonian has an exact symmetry under simultaneous permutation of the three bit-pair axes.
The consequence: the three ZZ correlators
$$\langle Z_0 Z_3 \rangle = \langle Z_1 Z_4 \rangle = \langle Z_2 Z_5 \rangle$$
are exactly equal in any eigenstate (and throughout unitary evolution).
(* The three ZZ correlators on a 6-qubit cell *)
(* For a state |ψ⟩, ZZ_{i,j} = ⟨ψ| Z_i Z_j |ψ⟩ *)
(* On Q₆ vertices (computational basis), Z_k |s⟩ = (-1)^{s_k} |s⟩ *)
zzCorrelator[s_, {i_, j_}] :=
(-1)^BitGet[s, i] * (-1)^BitGet[s, j]
(* For any state that respects the trigram symmetry,
the three ZZ values must agree *)
(* Example: check that the mismatch function mm(s) = popcount(upper XOR lower)
treats all three bit-pair axes equivalently *)
mm[s_] := DigitCount[BitXor[upperTrigram[s], lowerTrigram[s]], 2, 1]
(* mm values for all 64 states *)
mmValues = Table[mm[s], {s, 0, 63}];
Tally[mmValues]
(* {{0, 8}, {1, 24}, {2, 24}, {3, 8}} — symmetric distribution *)
(* The mm function is invariant under permutation of the three axes *)
(* This is what forces the three ZZ correlators to agree *)
On Noisy Hardware: The Gauge Spread
On a real quantum processor, noise breaks this symmetry. We define the gauge spread for each cell:
$$\sigma_i = \max_p \langle ZZ_{i,p} \rangle - \min_p \langle ZZ_{i,p} \rangle, \quad p \in \{0, 1, 2\}$$
In an exact computation, σᵢ = 0. On hardware, σᵢ > 0 tells you how badly noise has corrupted that cell.
(* Simulated gauge spread distribution from 230+ IBM jobs *)
(* Data from actual ibm_fez measurements *)
gaugeData = {
(* {topology, qubits, cells, gaugePass, meanSpread} *)
{"Cubic 2×2×2", 48, 8, "8/8", 0.038},
{"Slab 3×3×2", 108, 18, "15/18", 0.069},
{"Molecular graph", 108, 18, "16-18/18", 0.045},
{"Biological network", 108, 18, "12-15/18", 0.082},
{"Protein surface", 156, 26, "21-23/26", 0.061}
};
BarChart[
gaugeData[[All, 5]],
ChartLabels -> gaugeData[[All, 1]],
PlotLabel -> "Mean gauge spread by topology (ibm_fez)",
AxesLabel -> {None, "Mean σ"},
GridLines -> {None, {0.10}},
GridLinesStyle -> Directive[Red, Dashed, Thick],
ChartStyle -> {
RGBColor[0.2, 0.6, 0.3],
RGBColor[0.3, 0.5, 0.7],
RGBColor[0.2, 0.7, 0.4],
RGBColor[0.8, 0.3, 0.3],
RGBColor[0.9, 0.6, 0.2]
},
ImageSize -> 600,
Epilog -> {Red, Dashed,
Text["threshold σ = 0.10", Scaled[{0.85, 0.75}]]}
]
A cell is classified as reliable when σᵢ < 0.10 and unreliable when σᵢ ≥ 0.10. Across 230+ jobs, 86% of all cells pass gauge.
The Key Finding: Topology, Not Size
The most striking result is that gauge quality depends on the structural regularity of the input, not on the number of qubits:
(* Gauge pass rate vs qubit count — topology matters more than scale *)
points = {
{48, 100.0, "Cubic"}, (* 8/8 *)
{108, 83.3, "Slab"}, (* 15/18 *)
{108, 94.4, "Molecular"}, (* 17/18 avg *)
{108, 75.0, "Bio network"}, (* ~13.5/18 avg *)
{108, 83.3, "Cardiac"}, (* 15/18 *)
{156, 84.6, "Protein"} (* 22/26 avg *)
};
ListPlot[
{#[[1]], #[[2]]} & /@ points,
PlotStyle -> PointSize[Large],
PlotLabels -> points[[All, 3]],
AxesLabel -> {"Qubits", "Gauge pass rate (%)"},
PlotRange -> {{40, 170}, {60, 105}},
PlotLabel -> "Gauge quality: topology > qubit count",
GridLines -> Automatic,
ImageSize -> 600
]
The 108-qubit molecular graph (all iron atoms in equivalent oxidation states — maximally symmetric) achieves perfect 18/18 gauge. The 108-qubit biological network (irregular connectivity, heterogeneous perturbation) drops to 12–15/18. And the 156-qubit protein surface (26 cells!) achieves 21–23/26 — outperforming several 108-qubit configurations despite being 45% larger.
Symmetric inputs produce tighter gauge. This is the central empirical observation.
Cross-Domain Universality
The same architecture — identical circuit construction, identical gauge diagnostic, identical energy extraction — was applied to four completely unrelated domains. Only the perturbation values and bond topology changed:
(* Cross-domain results from ibm_fez *)
domainResults = <|
"Pharmacogenomics\n(108q, 100+ jobs)" -> <|
"Key result" -> "3 known drug interactions detected",
"Gauge" -> "14-16/18",
"Notable" -> "Zero training data"
|>,
"Catalyst selectivity\n(108q, 10 jobs)" -> <|
"Key result" -> "N₂/H₂O selectivity = +36.75",
"Gauge" -> "16-18/18",
"Notable" -> "Perfect gauge on symmetric Fe cluster"
|>,
"Cardiac stratification\n(108q, 20 jobs)" -> <|
"Key result" -> "48× disruption: disease vs healthy",
"Gauge" -> "11-15/18",
"Notable" -> "Disease degrades gauge (expected)"
|>,
"Protein formulation\n(156q, 2 jobs)" -> <|
"Key result" -> "ΔE_epistasis = −79.30",
"Gauge" -> "21-23/26",
"Notable" -> "First quantum 3-body excipient computation"
|>
|>;
Grid[
Prepend[
KeyValueMap[{#1, #2["Key result"], #2["Gauge"], #2["Notable"]} &,
domainResults],
Style[#, Bold] & /@ {"Domain", "Key Result", "Gauge", "Note"}
],
Frame -> All,
Spacings -> {2, 1},
Background -> {None, {LightBlue, {White, GrayLevel[0.95]}}},
ItemSize -> {{18, 22, 8, 22}},
Alignment -> Left
]
The gauge diagnostic operates identically in all four domains. This confirms it's a property of the Hamiltonian geometry, not an artifact of any specific application.
Non-Additive Interactions: What Entanglement Buys You
The lattice naturally computes epistatic (non-additive) interactions between system components:
$$\Delta E(A, B) = E(AB) - E(A) - E(B) + E(\text{WT})$$
(* The epistasis formula — domain-agnostic *)
epistasis[eAB_, eA_, eB_, eWT_] := eAB - eA - eB + eWT
(* Protein surface example (156 qubits, adalimumab) *)
eBare = 903.03;
eTrehalose = 951.37;
ePS80 = 986.67;
eBoth = 955.71;
expected = eTrehalose + ePS80 - eBare;
actual = eBoth;
dE = epistasis[eBoth, eTrehalose, ePS80, eBare];
Print["Expected (additive): ", expected] (* 1035.01 *)
Print["Actual (measured): ", actual] (* 955.71 *)
Print["ΔE_epistasis: ", dE] (* -79.30 *)
Print["Interpretation: synergistic competition"]
The two excipients interfere with each other on the protein surface — their combined effect is 79 energy units less than the sum of their individual effects. This is invisible to single-excipient screening and requires the inter-cell entanglement generated by the lattice.
What Makes This Different
| Property | VQE / QAOA | This architecture | |----------|-----------|------------------| | Parameters | Variational (many) | Zero (from N=6) | | Error handling | External (ZNE, PEC) | Built-in gauge | | Spatial resolution | None | Per-cell | | Domain change | Redesign circuit | Change dose + bonds | | Overhead | Additional circuits | Zero |
The Hamiltonian has zero free parameters — all coupling constants are determined by the cell dimension N = 6. Domain-specific information enters only through the perturbation values and bond topology.
Try It Yourself
The Q₆ hypercube and its gauge structure can be explored entirely within Wolfram Language:
(* Verify: the three qubit-pair axes of Q₆ are symmetric *)
(* Swap axes 0↔1 (permute bits {0,1,2,3,4,5} → {1,0,2,4,3,5}) *)
swapAxes01[s_] := Module[{bits = IntegerDigits[s, 2, 6]},
FromDigits[bits[[{2, 1, 3, 5, 4, 6}]], 2]
]
(* Check that mm is invariant under this swap *)
And @@ Table[mm[s] == mm[swapAxes01[s]], {s, 0, 63}]
(* True — the potential landscape has the claimed 3-fold symmetry *)
(* This symmetry is what forces the three ZZ correlators to agree *)
(* On real hardware, noise breaks it — and the breaking tells you where *)
The full preprint with hardware data (230+ IBM jobs, four domains) is available at: academia.edu/s/ddbebfc5c6
Patent pending: US Provisional 64/027,290 (April 3, 2026).
All quantum computations performed on IBM ibm_fez (Eagle r3, 156 qubits) via IBM Quantum Platform. The Hamiltonian structure is proprietary; the gauge symmetry it produces is fully characterized in the preprint.
Attachments:
