---
title: Circuit Cutting in Quantum Circuits
url: https://www.emergentmind.com/topics/circuit-cutting
type: topic
---

# Circuit Cutting in Quantum Circuits

Circuit cutting is a hybrid quantum–classical technique in which a circuit that is too large, too deep, or too nonlocal for the available hardware is partitioned into smaller fragments, the fragments are executed independently, and the target quantity for the original circuit is recovered by classical post-processing [2411.17756]. It was originally designed for expectation-value estimation, but the literature now spans spatial gate cuts for distributed execution, temporal wire cuts for sequential execution, full-output reconstruction, classical-shadow prediction, error mitigation, near-Clifford simulation, and direct sampling from trained variational circuits [2504.09334][2507.06649].

## 1. Foundations and cut models

Circuit cutting is motivated by the standard NISQ constraints: limited qubit count, limited coherence and depth, imperfect gate fidelity, and, in distributed settings, expensive inter-device communication [2504.09334]. In the most common taxonomy, **spatial circuit cutting** is the gate-cut approach that reduces circuit width by partitioning qubits across multiple QPUs and cutting inter-partition gates, whereas **temporal circuit cutting** is the qubit-wire-cut approach that reduces circuit depth by partitioning operations into sequential segments while keeping the same width [2504.09334]. The literature also distinguishes ancilla-free and ancilla-assisted wire-cutting models, and with or without classical communication between fragments [2410.20313][2306.12084].

The object being reconstructed depends on the application. In the canonical setting, cutting estimates an observable of the original circuit. In more recent work, the target may instead be the full output distribution, a family of observables reused through classical shadows, or the bitstring distribution of a trained QAOA ansatz used for optimization [2005.12702][2212.00761][2507.06649]. That distinction is substantive: recovering an expectation value is a linear task, whereas recovering a usable sample distribution is sensitive to distributional distortion, tail behavior, and finite-shot regularization.

The cut itself may be imposed on a wire, on a gate, or on a structured interface derived from a graph or hypergraph model of the circuit. Hypergraph-based formulations make this explicit by using a primal hypergraph for spatial cutting and a dual hypergraph for temporal cutting, so that the same partitioning machinery targets different physical costs: communication across QPUs in the spatial case, and qubit initialization or reset overhead in the temporal case [2504.09334].

## 2. Mathematical formalisms and reconstruction overhead

A standard quasiprobabilistic formulation expresses a nonlocal channel as a signed combination of local or communication-assisted channels,
\[
\mathcal U_{AB} = \sum_{i=1}^m a_i \mathcal F_i,
\]
and measures cost through the quasiprobability extent
\[
\gamma_S(\mathcal E) \coloneqq \inf\Big\{ \sum_{i=1}^m |a_i| : \mathcal E = \sum_{i=1}^m a_i \mathcal F_i,\; \mathcal F_i\in S,\; a_i\in\mathbb R \Big\},
\]
with single-cut sampling overhead factor \(\gamma_S(\mathcal E)^2\) [2503.22384]. In the corresponding Monte Carlo view, the number of shots scales as
\[
\sim \frac{2 a_1^2}{\epsilon^2}\log\!\left(\frac{2}{\delta}\right),
\]
so the \(\ell_1\)-norm of the reconstruction coefficients is the central cost parameter [2503.22384].

For wire cutting, a widely used Pauli-basis identity is
\[
\rho \simeq \frac{1}{2}\sum_{M\in\{I,X,Y,Z\}} M \otimes \operatorname{tr}_n(M_{(n)}\rho),
\]
which underlies the separation of upstream measurement and downstream preparation across a cut wire [2212.07335]. CutQC-style exact wire cutting generalizes this to \(k\) cuts through
\[
\rho=\frac{1}{2^k}\sum_{M\in\mathcal B^k}\bigotimes_{j=1}^{N_f}\rho_{f_j}(M(f_j)),
\]
but the number of fragment configurations grows exponentially with the number of cut interfaces [2410.20313].

Several later papers sharpen the overhead theory. Randomized-measurement wire cutting introduces the identity-channel decomposition
\[
\mathrm{id} = (2d+1)\,\mathbb E_z\!\left[(-1)^z \Psi_z\right]
\]
for a \(d=2^k\) dimensional cut subsystem, yielding \(\widetilde O(4^k/\varepsilon^2)\) sample complexity and an information-theoretic lower bound of \(\Omega(2^k/\varepsilon^2)\) for comparable wire-cutting procedures [2207.14734]. For gate cutting, the optimal overhead for an arbitrary two-qubit unitary with KAK coefficients \(u_i\) is
\[
\gamma_{\mathrm{LO}}(U)=\gamma_{\mathrm{LOCC}}(U)=1+2\sum_{i\ne j}|u_i||u_j|,
\]
and for \(n\) identical copies,
\[
\gamma(U^{\otimes n})=2(1+\Delta_U)^n-1,
\]
which is strictly cheaper than cutting the gates independently; the same work shows no advantage from classical communication for arbitrary two-qubit unitary cuts in this setting [2312.11638].

The role of side information has also been formalized through quantum instruments. In that framework, side information enlarges the operational decomposition set from channels to signed combinations of instrument branches, but its value depends strongly on the free-operation class: it is essential for feasibility in the no-communication LO setting, gives no extent reduction for states, and gives no advantage in SEP or PPT for Hermitian-preserving trace-preserving targets; whether \(\gamma_{\mathrm{LOCC}^\star}(\mathcal E) < \gamma_{\mathrm{LOCC}}(\mathcal E)\) can occur remains open [2503.22384].

## 3. Partitioning, placement, and decomposition-aware optimization

Because overhead grows rapidly with the cut interface, cut placement is itself an optimization problem. A general architecture-aware strategy represents the circuit as a hypergraph \(H=(V,E)\), with qubits as vertices and multi-qubit gates as hyperedges in the primal representation, then uses the dual hypergraph \(H^*=(E,X)\) when the goal is temporal wire cutting [2504.09334]. On top of this representation, Stoer–Wagner, Kernighan–Lin, and especially Fiduccia–Mattheyses heuristics are used to seek balanced bipartitions with few crossed hyperedges; the same work introduces the coupling base
\[
C_b=C(n,2)=\frac{n(n-1)}{2}
\]
and coupling ratio
\[
C_r=\frac{\text{total number of n-qubits gates}}{C_b}
\]
as structural indicators of cut difficulty [2504.09334].

Problem-graph-aware placement is particularly prominent for QAOA. In separator-based QAOA cutting for MaxCut, the minimum number of wire cuts per layer is tied to the cardinality of a minimum vertex separator of the problem graph, so the preprocessing problem becomes finding a low-cardinality balanced vertex separator, then shrinking it before any quantum decomposition is performed [2507.06649]. For \(p=2\) QAOA, the reported workflow shrinks the separator to a single vertex, producing exactly two wire cuts, one per layer, and this cut minimization is crucial because reconstruction overhead “increases exponentially with the number of cuts and quickly becomes prohibitive if circuit cutting is performed naively” [2507.06649]. In a related randomized-measurement analysis, a balanced vertex separator of size \(\kappa\) yields QAOA simulation overhead roughly \(2^{O(p\kappa)}\) [2207.14734].

Distributed systems place additional pressure on placement quality. FitCut converts the circuit DAG into a weighted graph whose vertices are two-qubit gates and whose edge weights count shared qubits, then applies constrained community detection and a capacity-aware merge heuristic [2405.04514]. Its primary objective proxies cut count through
\[
Obj_{nc} = \sum_{P_i \in P} QC_{P_i} - QC_{input},
\]
while a secondary objective maximizes worker utilization through idle-qubit penalties. In the reported experiments, FitCut reduces cut-search time by factors ranging from \(3\) to \(2000\) relative to Qiskit CKT and improves worker-side utilization by up to \(3.88\) times, with a system-wide improvement of \(2.86\) times [2405.04514].

A separate line of work observes that decomposition and cut placement should not be treated independently. For MCX and CCCX gates, partition-aware decompositions using a small number of ancillas can compress control information crossing the partition so that a multi-control gate is replaced by one or two cut CNOTs rather than many cross-boundary two-qubit gates [2603.26278]. The strongest reported variant, dec2Ad, lowers the asymptotic sampling overhead from \(O(6^{2n})\) in prior direct multi-control cutting to \(O(3^{2n})\) without classical communication and \(O(2^{2n})\) with classical communication, at the price of two ancillas and dirty-ancilla semantics [2603.26278].

## 4. Execution-level optimization and measurement reduction

Once a cut plan is fixed, a second layer of optimization targets the fragment execution burden itself. One example is ancilla-assisted, commuting-group wire cutting: quantum input initializations are converted into output measurements using ancillas, and the resulting Pauli observables are grouped into commuting families via MUBs-based grouping [2410.20313]. For a fragment with \(n_j^{qi}\) quantum inputs and \(n_j^{qo}\) quantum outputs, this reduces the number of required fragment executions from
\[
4^{n_j^{qi}}3^{n_j^{qo}}
\]
to
\[
2^{n_j^{qi}+n_j^{qo}+1}
\]
when \(n_j^{qo}>0\), while preserving exact reconstruction in the CutQC-style wire-cutting model [2410.20313].

Another strategy is to optimize the basis in which the cut state is represented. RICCO introduces a parameterized unitary at the cut so that, after optimization, many boundary observables vanish and the upstream fragment can be reconstructed from fewer measurement settings [2211.07358]. In the reported benchmarks, one cut required \(7\) circuit executions for QCUT and \(3\) for RICCO, while two cuts required \(25\) for QCUT and \(5\) for RICCO; however, the optimization overhead itself was large, about \(400\) circuit executions on average for the two-qubit-cut random-circuit case, and at least \(8\times\) the QCUT execution count in the VQE experiment [2211.07358].

A simpler exact pruning rule is to omit basis elements whose total contribution vanishes. “Golden cutting points” are cuts at which one Pauli basis element passes “no information” through the cut, so its downstream fragment executions can be skipped without changing the reconstructed observable [2304.04093]. For \(K_g\) golden cuts and \(K_r\) regular cuts, the paper reports reconstruction complexity \(O(4^{K_r}3^{K_g})\) and circuit-evaluation cost \(O(6^{K_r}4^{K_g})\), with up to \(33\%\) wall-time reduction on IBM hardware and no loss of accuracy in the tested instances [2304.04093].

State-dependent optimization extends this logic further. Because downstream cut fragments begin in known states and upstream fragments are measured in selected observables, state-dependent simplifications become broadly applicable after cutting [2506.05722]. The proposed SDO framework combines initial-state dependent optimization, measure-state dependent optimization, and biased observable selection; in noisy simulations of QAOA, QFT, and Bernstein–Vazirani, the method consistently mitigates noise and improves overall circuit performance [2506.05722].

## 5. Applications and specialized variants

A major branch of the literature uses circuit cutting to reconstruct **full output distributions** rather than single expectation values. Maximum-likelihood fragment tomography recasts each fragment as a Choi-like tomographic object, projects noisy fragment estimates back to the physical positive semidefinite cone, and then contracts the fragment models into a valid global probability distribution [2005.12702]. In numerical experiments on clustered random unitary circuits, MLFT consistently improves over direct fragment recombination and can estimate the output of a fragmented circuit with higher fidelity than full execution under the same shot budget [2005.12702].

Classical-shadow variants replace full fragment tomography by reusable randomized measurements on fragment Choi states. For observables that factorize across fragments, the divide-and-conquer shadow method gives per-fragment sample complexity
\[
K_f = 2 \log \frac{ 2 |F| \cdot 4^{\qdeg(f)} }{\delta}, \qquad
N_f = \frac{34 |F|}{\epsilon^2} 4^{\deg(f)+2|E|} \|O\|_{S}^4,
\]
and the numerical results show an advantage over uncut shadow tomography for high-weight observables [2212.00761]. The underlying mechanism is that the exponential dependence on observable weight can be traded for smaller per-fragment operator size, provided the cut overhead remains manageable.

Circuit cutting has also been repurposed for **error mitigation**. Simulated Quantum Error Mitigation cuts away the Pauli Check Sandwiching mitigation gadget, executes the main ansatz on hardware, simulates the small correction fragment classically, and recombines qubitwise-mitigated distributions into a full output distribution [2212.07335]. In hardware experiments on IBM Auckland with a two-local SU(2) VQE ansatz, the Bayesian-recombined SQEM energies were \(-7.8625\) vs \(-7.8634\) for LiH, \(-98.5653\) vs \(-98.5684\) for HF, and \(-15.4678\) vs \(-15.4706\) for BeH\(_2\), all relative to the corresponding noiseless ansatz values [2212.07335].

For **classical simulation**, circuit cutting can be combined with other structure. SuperSim isolates non-Clifford operations using cuts so that the large resulting fragments are mostly Clifford and can be simulated with Stim, while the residual non-Clifford fragments are handled exactly [2303.10788]. In the reported HWEA benchmarks, a runtime crossover appears around the mid-20-qubit range, and the framework scales to hundreds of qubits with modest runtimes on suitable near-Clifford instances [2303.10788].

The most conceptually distinct application is **pure sampling** from trained variational circuits. For depth-\(p=2\) QAOA on MaxCut, circuit cutting was applied not to expectation-value evaluation but to final bitstring sampling [2507.06649]. In noiseless simulation, the reconstructed distribution is broadened and shifted toward suboptimal objective values; the paper gives a per-bitstring suppression guarantee of at most \(1/\kappa=1/12\), so \(\kappa=12\) additional samples compensate the worst-case probability loss [2507.06649]. On IBM hardware, however, the 25-node instance shows the opposite effect: the reduced fragment width and depth lower hardware noise enough that the cut distribution is shifted closer to the optimum than the uncut distribution, making the relevant trade-off “distribution distortion versus noise mitigation” rather than the older “exactness versus overhead” framing [2507.06649].

## 6. Scalability, systems perspective, and open questions

Resource-estimation studies consistently show that circuit cutting is not a universally scalable substitute for larger hardware. A broad fault-tolerant benchmarking study using Qiskit cutting and Azure Quantum Resource Estimator reports that cutting can reduce physical-qubit requirements by about \(30\%\) on average, but for realistic Hamiltonian simulation, QFT, QPE, QAOA, block encoding, and random circuits, the exponential growth of quantum runtime and classical recombination usually dominates [2411.17756]. Representative figures include QFT on 10 qubits, where the baseline runtime is \(0.07\) s and the cut runtime is \(2070\) s, and Ising \(6\times6\), where the baseline is \(103.70\) s and the cut runtime is \(1.6\times 10^{11}\) s, despite substantial qubit reductions [2411.17756].

A systems view sharpens that conclusion. In distributed QNN training, a cut-aware estimator call decomposes into partitioning, subexperiment generation, execution, and reconstruction,
\[
T_{\mathrm{total}} = T_{\mathrm{part}} + T_{\mathrm{gen}} + T_{\mathrm{exec}} + T_{\mathrm{rec}},
\]
and the observed bottleneck is often reconstruction rather than fragment execution [2602.16233]. The reported median reconstruction shares are \(0.430\), \(0.489\), and \(0.530\) for one, two, and three cuts, respectively, with poor parallel scaling and speed-up often near \(1\) or below it; nevertheless, test accuracy and robustness are preserved in the measured regimes [2602.16233]. This suggests that asymptotic cut-count analysis alone is insufficient for estimator-driven workloads.

MAESTROCUT pushes this systems interpretation further by treating circuit cutting as an online control problem over partitioning, shot allocation, estimator choice, and secure dispatch [2509.00811]. Its partitioner optimizes
\[
J(\Pi)=\alpha\,\bar c(\Pi)+\beta\,\bar e(\Pi)+\gamma\,\bar q(\Pi),
\]
re-partitions under a CUSUM drift trigger, and allocates shots by the closed-form rule
\[
s_i^\star \;=\; \frac{S\,u_i^{2/3}}{\sum_{j=1}^n u_j^{2/3}},
\]
derived from a topology-aware variance surrogate [2509.00811]. Tier-2 emulation reports stable jitter, high reliability, and about \(1\%\) software overhead for the confidentiality layer under stress scenarios [2509.00811]. The underlying implication is that practical cutting on cloud-era NISQ hardware may depend as much on runtime orchestration as on the reconstruction identity itself.

Several open questions remain structural rather than implementational. The power of side information is still unresolved in the physically important LOCC setting [2503.22384]. For QAOA, simple “classical cutting” by enumerating separator assignments can be cheaper than layer-by-layer quantum cutting, especially as \(p\) grows [2507.06649]. More generally, the literature repeatedly returns to the same barrier: the number of cuts, the size of the cut interface, and the reconstruction layer jointly determine whether cutting functions as a practical scaling tool, a noise-mitigation tool, a classical-simulation aid, or merely an expensive formal reduction. This suggests that the modern subject is no longer just circuit decomposition, but the joint optimization of cut placement, fragment structure, estimator family, shot allocation, and runtime policy.

Source: https://www.emergentmind.com/topics/circuit-cutting