---
title: 'Routed Circuit Decomposition: Methods & Applications'
url: https://www.emergentmind.com/topics/routed-circuit-decomposition
type: topic
---

# Routed Circuit Decomposition: Methods & Applications

Searching arXiv for the cited papers and topic variants to ground the response.
In current research usage, routed circuit decomposition denotes a family of decomposition problems in which a circuit, unitary, process, or network is factored into simpler components while preserving, exploiting, or explicitly representing routing structure. In quantum compilation, this ranges from topology-agnostic synthesis passes whose outputs are later mapped and routed, to coupled decomposition-and-routing methods that optimize native-basis depth under sparse connectivity. In other literatures, closely related constructions describe route-constrained quantum processes built from direct sums and tensor products, pivot-based reduction of superconducting circuit topologies, uniquely routed path-length-matched photonic layouts, path decompositions obtained by cutting Eulerian routes, and sparse decompositions of circuit-switch schedules into permutations [2101.02993] [2308.03874] [2011.08120] [2412.06880] [1209.5144] [1910.02819] [2011.02752].

## 1. Terminological scope

Across these literatures, the common invariant is that decomposition is not performed on an abstract algebraic object alone. The object is decomposed together with a routing constraint: qubit connectivity, direct-sum sector flow, branch topology, waveguide trajectories, graph traversal order, or switch configurations. The resulting subcomponents are therefore not merely smaller; they are also route-compatible.

| Domain | Decomposition object | Route constraint |
|---|---|---|
| Quantum compilation | \(U \in U(2^n)\) or logical circuits | Hardware graph, native basis, SWAP overhead |
| Routed quantum processes | Linear maps, channels, supermaps | Relations between direct-sum sectors |
| Superconducting circuits | Lumped-element models | Node-loop or edge-network topology |
| Integrated photonics | 3D waveguide layouts | Unique trajectories, separation, path length |
| Graph theory | Eulerian circuits in quartic planar graphs | Local avoidance of short subcycles |
| Circuit switching | Traffic matrices | Permutation configurations over time |

A useful distinction in the quantum-compilation literature is between decomposition **before** routing and decomposition **with** routing. Exact unitary synthesis methods such as QSD and block-ZXZ produce structured gate lists over \(\{R_y,R_z,\mathrm{CNOT}\}\) or equivalent bases, but assume “perfect qubits” and leave topology to later compiler passes [2101.02993] [2403.13692]. By contrast, MIRAGE, TAP-plus-token-swapping routing, architecture-aware adaptive compression, and star-topology Cartan synthesis incorporate hardware structure directly into the decomposition objective or the decomposition algebra itself [2308.03874] [2206.01294] [2203.04426] [2506.16900].

This breadth also explains why the term is not tied to a single mathematical formalism. In some papers, routing is a graph-theoretic adjacency constraint; in others it is a relation on sector labels; in others it is literally the geometry of routed waveguides or the sequence of switch permutations. A plausible implication is that routed circuit decomposition is best understood as a structural perspective rather than a single algorithmic family.

## 2. Exact unitary synthesis before routing

A foundational strand treats routed decomposition as a downstream consequence of exact unitary synthesis. In the OpenQL-integrated QSD implementation, an arbitrary \(n\)-qubit unitary \(U \in U(2^n)\) is decomposed into single-qubit \(R_y(\theta)\), \(R_z(\theta)\), and two-qubit CNOT gates, intended as a mid-end synthesis pass before mapping, routing, and scheduling [2101.02993]. The core recursion applies Cosine–Sine Decomposition and demultiplexing, yielding at each level a circuit pattern consisting of a multiplexed \(R_z\), a multiplexed \(R_y\), another multiplexed \(R_z\), and four \((n-1)\)-qubit unitaries. Multi-controlled rotations are then realized by Gray-code ladders: a \(k\)-control rotation uses \(2^k\) CNOTs and \(2^k\) single-qubit \(R_a\) gates.

This structure is directly relevant to routed decomposition even though topology is ignored during synthesis. The output has a regular interaction pattern in which one pivot qubit repeatedly couples to many controls. The baseline QSD CNOT count is
\[
N_{\mathrm{CNOT}}(n)=\frac{3}{4}4^n-\frac{3}{2}2^n,
\]
while the optimized variant that stops at two-qubit gates gives
\[
N_{\mathrm{CNOT}}^{\mathrm{opt}}(n)=\frac{23}{48}4^n-\frac{3}{2}2^n+\frac{4}{3}.
\]
The same implementation reports that for \(n=10\), OpenQL produces about \(784{,}896\) CNOTs and about \(2.36\) million total gates, compared with about \(5.24\) million CNOTs and about \(6.81\) million total gates for Qubiter’s CSD-based flow; decomposition time is also \(10\)–\(100\times\) faster up to \(10\) qubits [2101.02993]. These numbers matter for routed decomposition because CNOT count is the main cost driver for later mapping and SWAP insertion.

The block-ZXZ program refines the same topology-agnostic paradigm. Its top-level factorization writes an \(n\)-qubit unitary in terms of block-diagonal \(A_1\oplus A_2\), a central block-ZXZ mixing term, and another block-diagonal factor, then demultiplexes each block into smaller unitaries plus uniformly controlled \(R_z\) gates [2403.13692]. The recursion remains QSD-like—four \((n-1)\)-qubit unitaries and three multiplexors per level—but replaces the central uniformly controlled \(R_y\) by a block-ZXZ central multiplexor that absorbs more two-qubit gates. The resulting CNOT count is
\[
N_{\mathrm{CNOT}}(n)=\frac{22}{48}4^n-\frac{3}{2}2^n+\frac{5}{3},
\]
which improves on optimized QSD by
\[
\frac{4^{n-2}-1}{3}
\]
CNOT gates, and gives a generic 3-qubit exact decomposition with \(19\) CNOT gates rather than \(20\) [2403.13692].

Both methods are therefore best viewed as topology-agnostic front ends to routed circuit decomposition. Their contribution is not hardware feasibility per se, but a recursively regular interaction skeleton—multiplexed rotations, Gray-code control ladders, and small exact leaf circuits—that later routing passes can exploit or, if necessary, must undo.

## 3. Coupled decomposition and routing on constrained quantum hardware

A second strand removes the separation between synthesis and routing. MIRAGE is explicit on this point: it treats transpilation as a collaborative design problem in which decomposition into native basis gates and connectivity-constrained routing are optimized together for iSWAP-family hardware [2308.03874]. Its key abstraction is the mirror gate \(U'=\mathrm{SWAP}\cdot U\), together with monodromy-polytope estimates of the minimum number of basis gates required for a two-qubit unitary. In the \(\sqrt{\mathrm{iSWAP}}\) basis, the paper emphasizes that CNOT and CNS \(=\) CNOT\(+\)SWAP have the same circuit-depth cost; more generally, mirror gates can absorb SWAPs that would otherwise be explicit routing overhead. For \(\sqrt{\mathrm{iSWAP}}\) with \(k=2\), Haar-weighted Weyl-chamber coverage rises from \(79\%\) to \(94.4\%\) when mirrors are allowed. On heavy-hex topologies, MIRAGE reports average depth reduction \(31.19\%\), weighted depth reduction \(33.41\%\), average SWAP reduction \(56.19\%\), and weighted SWAP reduction \(75.82\%\); on a \(6\times6\) square lattice, the corresponding figures are \(29.58\%\), \(32.09\%\), \(59.86\%\), and \(77.61\%\) [2308.03874].

A related but distinct decomposition treats routing itself as a two-level problem. The TAP-plus-token-swapping framework decomposes qubit routing by SWAP insertion into a global allocation subproblem and a family of token-swapping subproblems [2206.01294]. The token allocation problem is solved as a binary program whose objective is a lower bound on the number of SWAPs, strengthened by subgraph-isomorphism inequalities; the transition between consecutive allocations is then realized by exact or approximate token swapping. On real IBM hardware for \(p=1\) QAOA MaxCut on an 8-node ring, the routed circuit produced by this method uses \(25\) two-qubit gates and depth \(21\), compared with \(71\) and \(53\) for TKET and \(80\) and \(65\) for StochasticSwap, with measured approximation ratio \(0.64\) versus \(0.52\) and \(0.51\) [2206.01294]. Here routed circuit decomposition is literally a decomposition of a compiler optimization problem into route planning and route realization.

Adaptive circuit compression provides a third architecture-aware route. SQUANDER replaces the discrete search over gate placements by continuous optimization over circuits built from controlled-\(R_y(\theta)\) blocks and then sequentially removes two-qubit blocks while re-optimizing the remainder [2203.04426]. Connectivity enters by restricting which qubit pairs may host two-qubit building blocks. This makes the method usable either as direct architecture-aware synthesis or as post-routing compression of an already routed circuit. On a 5-qubit linear chain, the example `4gt13_91` is reported at \(26\) CNOTs for SQUANDER, versus \(35\) for QSearch LEAP, \(55\) for QFAST, and \(1042\) for Qiskit [2203.04426]. The route constraint is therefore not enforced through SWAP templates alone, but through the admissible parameterized ansatz.

Hardware-specific routed decomposition can also be algebraic rather than graph-search-based. For a star-shaped NV-center-plus-nuclear-spin processor, Cartan decomposition is tailored to the Lie algebra generated by the hardware Hamiltonian, so that arbitrary logical unitaries are factored into nuclear single-qubit rotations and diagonal multi-qubit phase gates implemented through the central electron spin [2506.16900]. In that architecture, there is no direct nuclear–nuclear coupling, and routed decomposition avoids explicit SWAP networks by using global diagonal control as a mediated interaction resource. The paper reports a \(3\,\mu\mathrm{s}\) optimized \(ZZ\) gate with fidelity \(0.9997\), a CNOT between two nuclear spins with total time \(11\,\mu\mathrm{s}\) and fidelity \(0.99992\), and random-CNOT experiments in which the decomposed circuit depth is roughly independent of the logical circuit depth and scales mainly with circuit width [2506.16900].

Taken together, these approaches show that routed circuit decomposition can mean at least three different things in constrained quantum hardware: absorbing routing into basis selection, decomposing routing into optimization subproblems, or deriving a hardware-native factorization from the device’s control algebra.

## 4. Routed quantum processes and causal structure

A more formal meaning of routed circuit decomposition appears in the theory of routed linear maps and routed quantum circuits. Standard circuit formalisms handle subsystem composition through tensor products, but several process-theoretic settings require nontrivial blends of direct sums and tensor products, such as one-particle subspaces or causal decompositions of unitary channels [2011.08120]. The canonical example is
\[
\widetilde{\mathcal H}_{AB}
=
\mathcal H_A^1\otimes\mathcal H_B^0
\;\oplus\;
\mathcal H_A^0\otimes\mathcal H_B^1,
\]
which is the one-particle sector of two lines but not the full \(\mathcal H_A\otimes\mathcal H_B\). Routed linear maps attach a relation \(\lambda: I_A\to I_B\) to a linear map \(f:\mathcal H_A\to\mathcal H_B\), restricting which input sectors may connect to which output sectors. This yields practical domains and codomains, practical isometries and unitaries, and a consistent dagger symmetric monoidal category of routed maps. An important notion is the accessible space of a slice of a routed circuit: it is not simply the tensor product of the visible wires, but the subspace selected by all routes in the full diagram [2011.08120].

Within this framework, routed circuit decomposition becomes a factorization of a process into routed components whose sector relations make direct-sum constraints explicit. The formalism accommodates superpositions of trajectories, superpositions of channels, and causal decompositions of unitaries that cannot be expressed as ordinary standard-circuit decompositions without enlarging Hilbert spaces or losing unitarity [2011.08120]. Extended circuit diagrams of Lorenz and Barrett are recovered as a special case in which routes are index matchings.

The 2025 QC-QC result turns this into a constructive existence theorem for indefinite causal order. For every \(N\), a single routed graph \(G_{\mathrm{QC\text{-}QC}(N)}\) suffices to obtain a routed circuit decomposition for any \(N\)-party quantum circuit with quantum control of causal order [2507.08781]. The routed graph is built from agent nodes \(\mathbf A_k\), internal nodes \(\mathbf V_{n+1}\), and index sets \(S_n^k\) that encode which parties have already acted and which party acts next. Validity is established through bi-univocality and branch-graph analysis, and the associated skeletal supermap is then fleshed out by routed isometries corresponding to the QC-QC internal maps \(\tilde V_{n+1}\) and adapter supermaps at the agent nodes. The result is that every QC-QC, including the quantum switch and higher-party generalizations, admits a routed circuit decomposition on a universal template graph [2507.08781].

This process-theoretic literature therefore uses “routing” in a stronger sense than hardware mapping. It refers to fine-grained information flow through direct-sum sectors, branch choices, and causal-control registers. A plausible implication is that routed circuit decomposition unifies hardware routing and causal routing only at a high level; technically, the two use different state spaces, different composition laws, and different notions of feasibility.

## 5. Physical-layout and network formulations

Outside quantum compilation proper, routed circuit decomposition often refers to extracting or designing physical topology under explicit routing constraints. In superconducting lumped-element models, the central object is the node–loop network matrix \(\mathbf{\Omega}\), whose entries encode how inductive loops are routed between capacitive nodes [2412.06880]. By passing to a tree–cotree basis and the edge network matrix
\[
\mathbf{\Omega}_E
=
\mathbf{A}_{C_\mathcal T}^{-1}\mathbf{\Omega}\mathbf{B}_{L_\mathcal T}^{-T},
\]
the decomposition becomes a sequence of row and column pivoting operations. Row pivots correspond graphically to removing a capacitive tree edge from the tree and placing it in parallel with an inductive edge; column pivots contract an inductive cotree edge and place it in series with a capacitive edge. The resulting “fundamental form” isolates harmonic LC modes, removes free modes by Schur complements, and separates the nonlinear Josephson-junction and phase-slip sectors from auxiliary linear modes [2412.06880]. The same network-matrix machinery is then used for model extraction from full-wave electromagnetic simulation, because the DC block of the hybrid response matrix directly yields \(\mathbf{\Omega}_E\).

In three-dimensional integrated photonics, routed decomposition is geometric rather than algebraic. The path-length-matched 8-waveguide pupil-remapping beam combiner is decomposed into uniquely routed waveguides connecting a hexagonal input arrangement to a linear output array, subject to minimum separation, fabrication depth, and non-cross-through writing constraints [1209.5144]. The design first used cubic splines, then piecewise circular arcs, and matched physical path lengths to within \(0.1\,\mu\mathrm{m}\) in the design trajectories. Interferometric measurements on the fabricated device showed optical path-length matching to within \(45\,\mu\mathrm{m}\). The optimized design used radii of curvature greater than \(23\,\mathrm{mm}\), a minimum waveguide separation of \(30\,\mu\mathrm{m}\), and a \(30\,\mathrm{mm}\) chip length [1209.5144]. Here routed decomposition means that the device-level circuit is decomposed into individual trajectories whose geometry simultaneously satisfies routing, loss, and phase constraints.

A network-scheduling variant appears in high-speed circuit switching. A normalized traffic demand matrix \(X^\star\) is decomposed into a sparse sum of permutation matrices,
\[
X_k=\sum_{i=1}^k \theta_i P_i,
\]
where each permutation is a switch configuration and each \(\theta_i\) is its duration [2011.02752]. Classical Birkhoff decomposition gives exact realizability but poor sparsity; the revisited theory establishes that with an admissible subset of permutations, Birkhoff-type algorithms achieve an \(\epsilon\)-approximate decomposition with \(O(\log(1/\epsilon))\) permutations, and Birkhoff+ combines Frank–Wolfe direction selection with Birkhoff step sizing to obtain sparse decompositions using one LP per configuration [2011.02752]. In circuit-switch performance experiments, Birkhoff+(10) achieves about \(7\%\) higher throughput than Eclipse at \(\delta/W=10^{-2}\), while Birkhoff+ and Birkhoff+(10) are about an order of magnitude faster than Solstice and Eclipse for small \(\epsilon\) [2011.02752]. The circuit here is not a quantum circuit, but the decomposition principle is the same: global demand is represented as a sparse sequence of routed configurations.

These physical-layout literatures show that routed circuit decomposition is not inherently computational. It can equally describe topological reduction, geometric embedding, or schedule sparsification, provided route compatibility is part of the decomposition itself.

## 6. Abstract graph models, complexity, and open problems

At a combinatorial extreme, routed circuit decomposition becomes decomposition of a route into prescribed segments. For connected quartic planar graphs, a 4-locally self-avoiding Eulerian circuit exists if and only if the graph does not contain \(F_6\) as a subgraph [1910.02819]. This route can then be cut into paths of prescribed lengths: a quartic planar graph of order \(n\) admits a decomposition into \(k_1+k_2+k_3+k_4\) paths with \(k_i\) copies of \(P_{i+1}\) if and only if
\[
k_1+2k_2+3k_3+4k_4=2n.
\]
In particular, every connected quartic planar graph of even order admits a \(P_5\)-decomposition [1910.02819]. The same “route first, cut later” logic reappears in more concrete routed decomposition problems, although there the route is usually a sequence of hardware interactions, sector transitions, or switch configurations rather than an Eulerian traversal.

Across the surveyed literatures, the dominant limitations are different but structurally similar. Exact generic quantum synthesis scales exponentially in circuit width: QSD has gate counts \(O(4^n)\) and decomposition time roughly \(O(6^n)\), while block-ZXZ, though more CNOT-efficient, remains above the lower bound \(\frac14(4^n-3n-1)\) and does not optimize depth or topology directly [2101.02993] [2403.13692]. MIRAGE’s monodromy analysis is limited to two-qubit blocks, and its routing-plus-decomposition decisions remain heuristic even though they are physically informed [2308.03874]. SQUANDER’s continuous optimization is powerful on 3–5 qubits but expensive for larger blocks [2203.04426]. The NV-center Cartan approach is exponential in \(N\) in the generic case and is therefore suited primarily to few-qubit central-spin registers [2506.16900].

In the process-theoretic domain, routed circuits now cover all QC-QCs, but whether all purifiable process matrices admit routed circuit decompositions remains unresolved [2507.08781] [2011.08120]. In superconducting model extraction, the methods assume reciprocal, lossless linear parts and no ideal transformers [2412.06880]. In photonics, path-length matching remained highly manual and optical mismatches were still tens of microns despite sub-micron geometric matching [1209.5144]. In graph theory, extensions beyond quartic planar graphs or beyond path lengths up to \(4\) edges remain open [1910.02819].

These recurrent limitations indicate that routed circuit decomposition is usually successful when route structure is regular enough to be represented explicitly—Gray-code ladders, mirror-gate choices, branch graphs, pivotable network matrices, uniquely parameterized waveguide paths, or Eulerian tours. When the route space is high-dimensional, noisy, or only weakly structured, decomposition remains possible but typically loses optimality guarantees, scalability, or both.

Source: https://www.emergentmind.com/topics/routed-circuit-decomposition