---
title: Heavy-Hex Graph Ising Models
url: https://www.emergentmind.com/topics/heavy-hex-graph-ising-models
type: topic
---

# Heavy-Hex Graph Ising Models

Heavy-hex graph Ising models are Ising-type spin systems whose interaction structure is either the heavy-hex graph itself or a lattice that is compiled onto heavy-hex hardware. In the literature, the topic therefore spans several closely related settings: classical quadratic and geometrically local cubic spin glasses defined directly on heavy-hex graphs, transverse-field and kicked Ising dynamics on the heavy-hex lattice, and line-graph Ising circuits that are deterministically routed onto heavy-hex processors. Across these settings, the heavy-hex graph functions simultaneously as a graph-theoretic object, a native superconducting-qubit coupling map, and a constrained platform for studying compilation overhead, dynamical phenomena, and classical-versus-quantum scaling [2306.05385] [2412.15572] [2402.01558].

## 1. Graph-theoretic structure of the heavy-hex lattice

The heavy-hex graph is a sparse, planar, bipartite, geometrically local graph with bounded degree. In the construction used for large-scale random Ising studies, every generated heavy-hex graph is bipartite with vertices partitioned as \(V = V_2 \sqcup V_3\) such that \(E \subset V_2 \times V_3\); vertices in \(V_2\) have maximum degree \(2\), vertices in \(V_3\) have maximum degree \(3\), and the number of edges scales linearly with the number of vertices \(n\), producing \(O(n)\) edges [2412.15572]. The transverse-field literature describes the same structure as a decorated hexagonal lattice: starting from a regular hexagonal lattice, every edge is decorated with one additional site, yielding two sublattices, \(A\) of coordination number \(2\) and \(B\) of coordination number \(3\), with no \(A\!-\!A\) or \(B\!-\!B\) nearest-neighbor bonds [2402.01558].

In the routing literature, heavy-hex appears as the heavy graph \(\mathrm{heavy}(G)\) of a hexagonal base lattice \(G\). One places a node on every edge of \(G\) and connects that new node to the two original vertices incident to the edge. The new nodes are the “heavy nodes,” while the original vertices serve as “light” mediator nodes. This representation is central for compiling line-graph interaction patterns onto heavy-hex hardware because algorithm qubits can be assigned to heavy nodes and routed through unique mediators on the light sublattice [2306.05385].

The architecture is also a hardware design choice. IBM’s superconducting processors use heavy-hex connectivity to mitigate frequency-collision and crosstalk errors while maintaining a manufacturable sparse topology. For the 127-qubit heavy-hex processor considered in magnetization-dynamics studies, the hardware graph has 144 two-qubit couplers and is 3-edge-colorable, a combinatorial property that enables parallel scheduling of all nearest-neighbor ZZ interactions in depth \(3\) [2311.01657].

## 2. Hamiltonian families defined on heavy-hex graphs

Several distinct Hamiltonian classes appear under the heading of heavy-hex graph Ising models. They differ in whether the model is classical or quantum, whether interactions are quadratic or higher order, and whether the graph is used directly or as a hardware target.

| Family | Hamiltonian | Structural rule |
|---|---|---|
| Quadratic classical Ising | \(C(z)=\sum_{v\in V} d_v z_v + \sum_{(i,j)\in E} d_{i,j} z_i z_j\) | Pairwise terms only on heavy-hex edges |
| Quadratic + local cubic Ising | \(C(z)=\sum_{v\in V} d_v z_v + \sum_{(i,j)\in E} d_{i,j} z_i z_j + \sum_{l\in W} d_{l,n_1(l),n_2(l)} z_l z_{n_1(l)} z_{n_2(l)}\) | Cubic terms centered on degree-2 vertices |
| Heavy-hex TFIM | \(H=-J\sum_{\langle ij\rangle}\sigma_i^z\sigma_j^z + h\sum_i \sigma_i^x\) | Uniform nearest-neighbor ferromagnetic couplings on heavy-hex |
| Kicked heavy-hex TFIM | \(H=-(\pi/4)\sum_{\langle i,j\rangle} Z_i Z_j + (\theta_h/2)\sum_i X_i\) | Stroboscopic IBM circuit model on heavy-hex |
| XX+Z transverse-field model | \(\hat H = B\sum_i \hat\sigma_i^Z + \sum_{\langle i,j\rangle\in E} J_{ij}\hat\sigma_i^X\hat\sigma_j^X\) | Native heavy-hex edge set used directly |

For random combinatorial optimization instances, the coefficients \(d_v\), \(d_{i,j}\), and, when present, \(d_{l,n_1(l),n_2(l)}\) are independently drawn from \(\{+1,-1\}\) with probability \(0.5\) each. The cubic extension is geometrically local by construction: each 3-local term involves a degree-2 vertex and its two neighbors, so every cubic monomial occupies a local heavy-hex motif rather than an arbitrary hyperedge [2412.15572].

For heavy-hex-compatible QAOA studies, the cost Hamiltonian is written as
\[
H_C = \sum_{i \in V} h_i Z_i \;+\; \sum_{(i,j)\in E} J_{ij}\, Z_i Z_j \;+\; \sum_{(l,n_1(l),n_2(l))\in T} K_{l n_1(l) n_2(l)}\, Z_l Z_{n_1(l)} Z_{n_2(l)} ,
\]
with \(E \subset V_2 \times V_3\) and \(T\) consisting of triples centered on degree-2 vertices. This model class is “heavy-hex compatible” because all terms are supported on native couplers or native two-edge motifs of the hardware graph [2312.00997].

A broader equilibrium formulation also exists. For ferromagnetic Ising models with general external field on a finite graph \(G=(V,E)\), the Hamiltonian is
\[
H(\sigma)= -\!\!\sum_{e=(i,j)\in E}J_e\,\sigma_i\sigma_j - \sum_{i\in V} h_i\,\sigma_i,
\]
with \(J_{ij}\ge 0\). Because the random-cluster and Edwards–Sokal formalism applies on any finite graph and the heavy-hex graph is planar, bounded-degree, and amenable, this formulation specializes directly to finite heavy-hex patches [1506.06818].

## 3. Native compilation, mediated routing, and hardware realization

When the Ising model is already defined on the heavy-hex graph, heavy-hex processors admit direct spin-to-qubit mapping with no routing SWAPs. This is the regime targeted in heavy-hex QAOA and short-time evolution studies. A 3-edge-coloring of the heavy-hex coupling map partitions couplers into three disjoint matchings, so all pairwise phase-separation gates can be scheduled in parallel within each color class. For a quadratic term \(e^{-i \gamma J_{ij} Z_i Z_j}\), the compilation uses a parity-compute \(R_Z\)-uncompute pattern,
\[
\mathrm{CNOT}(i\!\to\!j),\ R_Z(2\gamma J_{ij})\ \text{on }j,\ \mathrm{CNOT}(i\!\to\!j).
\]
For a cubic term on a degree-2 center \(l\) with neighbors \(n_1,n_2\), the same parity chains are reused:
\[
\mathrm{CNOT}(n_1\!\to\!l),\ \mathrm{CNOT}(n_2\!\to\!l),\ R_Z(2\gamma K_{l n_1 n_2})\ \text{on }l,\ \mathrm{CNOT}(n_2\!\to\!l),\ \mathrm{CNOT}(n_1\!\to\!l).
\]
In the layered schedule, cubic terms therefore add only one additional \(R_Z\) per triple and do not increase the CNOT depth, which remains \(6\) per QAOA layer independent of instance size [2312.00997]. The same strict two-qubit depth \(6\) per layer is retained in later high-round parameter-transfer experiments on 127-, 133-, and 156-qubit heavy-hex devices, where a \(156\)-qubit circuit at \(p=49\) contains \(17{,}248\) CZ gates [2509.13528].

A different use of heavy-hex arises when the algorithmic interaction graph is not heavy-hex but a line graph \(L(G)\). In that setting, the hardware graph is \(\mathrm{heavy}(G)\), and every logical two-qubit ZZ interaction on \(L(G)\) is mediated through a unique heavy-hex light node. If
\[
U_{ij}(\theta)=e^{-i\theta Z_i Z_j},
\]
then line-graph routing replaces it by
\[
M\!U^{(\theta)}_{i m j}:=\mathrm{SWAP}_{m i}\,\Big(e^{-i \theta Z_m Z_j}\Big)\,\mathrm{SWAP}_{i m},
\]
where \(m\) is the unique mediator with \((i,m)\) and \((m,j)\) hardware edges. The routing algorithm first tests whether the coupling graph is a line graph via \(L^{-1}\), builds \(H=\mathrm{heavy}(G)\), replaces every ZZ gate by the mediated form, removes superfluous SWAPs, and optionally embeds \(H\) into the labeled device graph using VF2 [2306.05385].

The formal guarantee is explicit: every circuit with coupling graph \(L(G)\) can be performed on hardware with coupling graph \(\mathrm{heavy}(G)\) with
\[
\text{SWAP overhead} \le 2 \times |G_{2q}|,
\]
and overall classical routing cost
\[
T_{\text{route}} = O(N^2),
\]
with \(N\) the total gate count. For the kagome-to-heavy-hex benchmark, the routed circuit has depth \(226\), SWAP count \(7{,}968\), active qubits \(300\), and routing time \(43\) s; compared to SABRE at its highest optimization level, line-graph routing achieved \(\sim 3.5\times\) smaller depth while also using fewer SWAPs [2306.05385].

A common misconception is that heavy-hex itself is a line graph. The routing construction states the opposite: heavy-hex is the hardware target \(\mathrm{heavy}(G)\) when the algorithmic coupling graph is the line graph \(L(G)\), as in kagome \(=L(\text{hexagonal})\) [2306.05385].

## 4. Dynamical phenomena on heavy-hex Ising lattices

The heavy-hex transverse-field Ising model exhibits a dynamical regime in which the lattice geometry itself induces confinement. For
\[
H = -J \sum_{\langle i j \rangle} \sigma_i^z \sigma_j^z + h \sum_i \sigma_i^x,
\]
defined on the decorated hexagonal heavy-hex lattice, a quench from the broken-symmetry state \(|Z^+\rangle\) at small \(h/J\) produces persistent oscillations, suppressed correlation spreading, and saturation of the entanglement entropy. The mechanism is geometric: flipping a single \(A\) site creates two domain walls and costs \(4J\), whereas flipping a single \(B\) site together with any subset of its three neighboring \(A\) sites creates three domain walls and costs \(6J\). The effective confining potential is summarized heuristically as
\[
V_{\text{eff}} \approx 2J \times N_{\text{dw}},
\]
and the low-energy zero-momentum sector is captured by a \(5\times 5\) effective Hamiltonian whose eigenvalues define five hadronic quasiparticle masses. For quenches from \(|Z^+\rangle\), the regime \(h/J=0.1,0.2,0.3\) shows long-lived oscillations and bounded entanglement, whereas \(h/J=1.0\) shows linear growth of entanglement entropy density, spreading correlations, and relaxation consistent with ETH [2402.01558].

A closely related heavy-hex dynamics problem was implemented as a kicked ferromagnetic TFIM on IBM’s 127-qubit heavy-hex architecture. The Hamiltonian is parameterized as
\[
H = - \frac{\pi}{4} \sum_{\langle i,j\rangle} Z_i Z_j + \frac{\theta_h}{2} \sum_i X_i,
\]
with first-order Trotterization
\[
e^{-i N \delta_t H} = \prod_{p=1}^N \left[ \left(\prod_{\langle i,j\rangle} R_{Z_i Z_j}(-\pi/2)\right)\left(\prod_i R_x(\theta_h)\right) \right].
\]
Because the heavy-hex graph is 3-edge-colorable, the ZZ layer fits in ECR/CNOT depth \(3\). This model was simulated on D-Wave Pegasus quantum annealers by matching the ratio of transverse-field and Ising energy scales through an anneal pause \(s_p\) satisfying
\[
\frac{A(s_p)}{B(s_p)J} = - \frac{2\theta_h}{\pi},
\]
together with pause-time matching equations for the ZZ and \(X\) actions. Using reverse quantum annealing and h-gain state encoding, the study reported equivalent magnetization dynamics for \(N=20\), \(50\), and up to \(10{,}000\) time steps, consistent with exact classical \(27\)-qubit heavy-hex Trotterized circuit dynamics for whole-lattice magnetization [2311.01657].

Heavy-hex graph Ising models have also been used for short-time quantum-evolution-based ground-state estimation in the XX+Z transverse-field model
\[
\hat H = B \sum_i \hat\sigma_i^Z + \sum_{\langle i,j\rangle\in E} J_{ij} \hat\sigma_i^X \hat\sigma_j^X.
\]
Here the heavy-hex edge set is used directly, so no minor embedding is required. The computation combines shallow Trotterized sampling with a classical subspace diagonalization step, and success is diagnosed by the information-theoretic ratio
\[
R = \log_2 \left[ \frac{I_{\mathcal B}}{K I_\Psi} \right].
\]
On ibm_torino, with median CZ error \(2.5\times 10^{-3}\) and median readout assignment error \(3.0\times 10^{-2}\), reliable ground-state energy approximations were reported up to \(63\) heavy-hex qubits when \(|J|/|B|\) is not too large and \(R>0\) [2604.25715].

These dynamical studies collectively show that heavy-hex geometry is not merely a hardware constraint. It also organizes the excitation spectrum, limits or enhances entanglement growth depending on the quench protocol, and strongly affects which observables remain accessible on current devices.

## 5. Classical complexity, equilibrium structure, and exact graphical representations

For random quadratic and geometrically local cubic Ising models defined directly on 2D heavy-hex graphs, the dominant structural fact is sparsity. The graph has maximum degree \(3\), is planar and bipartite, and supports only local pairwise and local cubic interactions. In an empirical scaling study covering \(10\) random instances per size for \(n\) from \(100\) up to \(160{,}000\), simulated annealing with \(100\) or \(1000\) sweeps showed evidence of exponential Time-to-Solution scaling on quadratic-only instances, whereas Gurobi solved the same sparse instances efficiently on average. For quadratic-only instances, the best fit for Gurobi IQP runtime was linear with \(a \approx 10^{-4}\), \(b \approx 10^{-4}\), and RMSE \(=0.683\). With cubic terms, ILP lifting showed linear scaling with \(a = 0.00791\), \(b = 10^{-4}\), and RMSE \(=27.58\), while IQP order reduction showed weakly quadratic scaling with \(a \approx 10^{-6}\), \(b = 0.0802\), \(c \approx 10^{-6}\), and RMSE \(=192.6\). The central conclusion was that these sparse heavy-hex instances are not classically hard on average for strong exact solvers, even though a simple heuristic such as simulated annealing can display weak exponential scaling [2412.15572].

This empirical result is important because heavy-hex-native instances are often chosen precisely to incur little or no embedding overhead on current quantum hardware. A second common misconception is therefore that hardware-native automatically implies classically challenging. The large-instance Gurobi study shows the opposite for the random heavy-hex ensembles tested: sparse locality can make exact classical optimization efficient on average [2412.15572].

The equilibrium statistical mechanics of ferromagnetic heavy-hex Ising models can be analyzed exactly through random-cluster and Edwards–Sokal representations. For a finite heavy-hex patch \(G=(V,E)\), with \(p_{ij}=1-e^{-2\beta J_{ij}}\), the random-cluster weight is
\[
\phi_{\boldsymbol{p},\boldsymbol{h},G}(\omega)
=
\frac{1}{\mathscr{Z}^{\mathrm{RC}}_{\boldsymbol{p},\boldsymbol{h},G}}
\,B_{\boldsymbol{J}}(\omega)\,
\prod_{\alpha=1}^{k(\omega,G)} 2\cosh\big(\pmb{h}(K_\alpha)\big),
\]
where \(\pmb{h}(K)=\beta\sum_{i\in K} h_i\). The exact single-site magnetization is
\[
m_x = \phi_{\boldsymbol{p},\boldsymbol{h},G}\big(\tanh(\pmb{h}(K_t))\big),
\]
and the exact spin-spin correlation is
\[
\lambda_{\beta,\boldsymbol{h},V}(\sigma_x\sigma_y)
=
\phi_{\boldsymbol{p},\boldsymbol{h},G}(x\leftrightarrow y)
+
\phi_{\boldsymbol{p},\boldsymbol{h},G}\!\big(
\mathds{1}_{\{x\not\leftrightarrow y\}}\,
\tanh(\pmb{h}(K_t))\,\tanh(\pmb{h}(K_u))
\big).
\]
Because heavy-hex graphs are planar, periodic, bounded-degree, and amenable, the associated GRC measures inherit strong FKG, uniqueness of the infinite connected component under free and max-wired limits, and quasilocality [1506.06818].

For classical sampling beyond cluster algorithms, the Hybrid Monte Carlo formalism offers a graph-agnostic continuous-variable representation of the Ising model. After a Hubbard–Stratonovich transformation with a stabilizing shift \(K̃ = K + C I\), HMC targets
\[
S(\psi) = \frac{1}{2}\psi^\top\psi - \sqrt{J'}\, \psi\!\cdot\! h' - \sum_{i=1}^N \log\cosh(\sqrt{J'}\,\psi_i),
\]
with leapfrog dynamics driven by
\[
\dot \psi = p,\qquad
\dot p = -\psi + \sqrt{J'}\, h' + \sqrt{J'}\, \tanh(\sqrt{J'}\,\psi).
\]
In the reported benchmarks, acceptance near criticality was about \(70\)–\(80\%\), and the integrated autocorrelation time of \(|m|\) scaled with dynamic exponent \(z \approx 2\) in \(d=2\) and \(d=3\). Since the method is explicitly formulated for arbitrary lattices and couplings, it carries over directly to heavy-hex graphs by building the symmetric coupling matrix from the heavy-hex edge set [1912.03278].

## 6. QAOA, parameter transfer, and the current research frontier

Heavy-hex-compatible higher-order Ising spin glasses have become a standard benchmark for whole-chip QAOA because the graph structure allows constant-depth compilation. In noiseless simulations on random \(\pm 1\) linear, quadratic, and cubic heavy-hex instances, a single set of angles trained on one 16-qubit instance transferred effectively to ensembles at 16, 27, and 127 qubits for \(p=1\) through \(p=5\). For 127 qubits, the average noiseless expectations across the ensemble were reported as roughly \(-85\), \(-110\), \(-130\), \(-140\), and \(-145\) for \(p=1,\dots,5\). The \(p=1\) mean energy landscapes were nearly identical across 27-, 127-, and 414-qubit devices, with the best angles concentrating near \(\beta \approx 0.4\)–\(0.47\) and \(\gamma \approx 2.83\)–\(2.86\) modulo periodicity. On hardware, the best 27-qubit processors improved mean energies up to \(p=3\), the best 127-qubit processors up to \(p=2\), and noisy means remained about a factor of two worse in magnitude than the noiseless expectations [2312.00997].

Later work pushed the same heavy-hex-compatible program to much larger QAOA depths. Angles trained classically on 16- and 27-qubit donor instances were transferred to 127-, 133-, and 156-qubit targets, with classical validation by statevector, PEPS, MPS, and LOWESA. The central result is qualified rather than monotone: transferred schedules do not in general improve strictly with \(p\), and temporary performance decreases occur, but the overall trend is toward improved expectation value with increasing depth, in many cases converging close to the true ground-state energy of the \(100+\)-qubit instances. On hardware, continuous improvement of solution quality was observed up to \(p=5\) on ibm_fez, \(p=9\) on ibm_torino, and \(p=10\) on ibm_pittsburgh, with no dynamical decoupling or error mitigation [2509.13528].

Several boundaries of the field are now clear. Heavy-hex-native compilation is exceptionally favorable: it avoids routing, fixes two-qubit depth at \(6\) per QAOA layer, and accommodates geometrically local cubic terms without extra entangling depth [2312.00997] [2509.13528]. At the same time, the random heavy-hex optimization problems most compatible with current hardware are structurally sparse, and sparse heavy-hex instances can be easy for strong classical exact solvers on average [2412.15572]. In dynamical settings, heavy-hex geometry can either suppress entanglement through geometric confinement or support long-time analog emulation through native embedding on annealers, depending on the Hamiltonian and observable [2402.01558] [2311.01657].

Open directions reported in the literature include denser and/or nonlocal heavy-hex-compatible benchmarks, alternative classical heuristics such as cluster updates and parallel tempering, measurement of structural quantities such as treewidth and cutwidth, dynamic endpoint selection in line-graph routing to enhance SWAP cancellation, angle extrapolation beyond the trained QAOA depth, and lightweight online refinement of transferred schedules [2412.15572] [2306.05385] [2509.13528]. Taken together, these directions suggest that heavy-hex graph Ising models are best understood not as a single model, but as a tightly connected family of graph-native statistical-mechanical and quantum-circuit problems whose defining feature is the persistent interaction between lattice geometry, hardware constraints, and algorithm design.

Source: https://www.emergentmind.com/topics/heavy-hex-graph-ising-models