Circuit Optimization
Short Definition
Section titled “Short Definition”Circuit optimization transforms an existing quantum or quantum–classical circuit into a circuit with the same declared behavior and lower cost under a stated objective. The transformation may cancel inverse gates, fuse rotations, commute operations to expose cancellations, resynthesize small windows, exploit Clifford or phase-polynomial structure, shorten a dependency graph, or search a larger equivalence class. It is not merely “removing gates”: an optimizer can deliberately add gates or ancillas to reduce depth, non-Clifford cost, duration, or another more important resource.
A complete contract has the schematic form
Here is the input circuit, defines behavioral equivalence, is the legal operation set, is the ancilla and classical-control contract, is the cost model, and contains limits on compiler time, memory, and approximation. The output comes with a certificate identifying the transformations, equivalence evidence, costs, assumptions, and tool versions.
This page is the canonical home for circuit rewriting, simplification, commutation, and abstract depth reduction. Gate Decomposition owns synthesis of an operation into a gate alphabet. Circuit Intermediate Representations owns the representation and capability contract. Later pages in this chapter own logical-to-physical placement, routing, calibration- and noise-aware choices, and pulse scheduling. A compiler may interleave those tasks, but their claims and evidence remain distinct.
Optimization Begins with Semantics
Section titled “Optimization Begins with Semantics”For a static unitary circuit on the same input and output wires, a common exact contract is
where a global phase is allowed only if the surrounding context makes it unobservable. If the block can later be controlled, used on one interferometer branch, or compared with a phase reference, literal equality or an explicitly propagated phase correction may be required.
Approximate optimization must name both metric and budget. One possible contract is
Near-equality is not an automatic license to delete a small rotation. It is a different compiler mode whose error must be composed with algorithmic, synthesis, numerical, and physical error budgets. Error Estimates develops the distinction between a residual, a bound, and a certified tolerance.
Unitary equality is insufficient for dynamic circuits. Measurements, resets, discarded systems, and classical branches define a channel or, when outcomes are retained, a quantum instrument. If is the trace-nonincreasing map associated with output record , exact instrument-level equivalence asks for
This requirement preserves both and the conditional output state. Preserving only the outcome distribution is weaker; preserving only the unconditioned channel is weaker still. Quantum Instruments is the canonical treatment of these distinctions.
Ancillas introduce another boundary. An exact clean-ancilla rewrite may need
for every data state , including restoration of all ancillas. Equality only on computational-basis inputs, or only after tracing out garbage, is a different contract and cannot be substituted silently.
Cost Is a Vector
Section titled “Cost Is a Vector”Useful circuit costs are rarely ordered by a single universal number. A target-independent optimizer might retain
where denotes counts, denotes abstract depths, and is compiler runtime. A present-day device may care most about entanglers and duration. A fault-tolerant architecture may prioritize magic-state demand, reaction depth, or spacetime volume. A simulator may prefer gates that admit efficient fusion even when the emitted circuit is not best for hardware.
There are three defensible ways to compare cost vectors:
- a lexicographic policy, such as minimizing count before Clifford count;
- a weighted score whose units and weights are recorded;
- a Pareto policy that retains every nondominated candidate.
The statement “ is optimized” is incomplete without one of these choices. Even “minimum depth” is ambiguous until gate durations, concurrency, resource conflicts, measurement latency, and connectivity are fixed.
Optimization is a translation-validation loop. Candidate rewrites are accepted only when they satisfy the declared equivalence and improve the selected cost policy. Rejected or dominated candidates return to pass selection; accepted circuits retain a pass trace and certificate.
Expose the Dependency Structure
Section titled “Expose the Dependency Structure”A textual gate list contains an arbitrary choice of order among operations that could run or move independently. Optimizers therefore construct a dependency directed acyclic graph
where vertices are operations and edges encode quantum-resource order, classical data flow, timing constraints, and semantic fences. Every legal schedule is a topological ordering that also satisfies the target capability contract.
Operations on disjoint subsystems commute. Operations sharing a subsystem may also commute, but this must follow from their semantics rather than their diagrammatic appearance. For unitary gates and ,
For Pauli strings, write a Pauli modulo phase as with support and support . Then
The strings commute exactly when . Their rotations then commute as well:
Gate-specific commutation rules are often cheaper to check. For ,
The analogous and moves are not generally valid. A rewrite engine should encode the positive rule and its side conditions, not infer a visual symmetry that the gate does not possess.
Barriers require a declared meaning. A semantic fence forbids movement because an external timing, calibration, or observation contract depends on the boundary. A display-only barrier may be ignored. Treating both spellings as the same object either blocks useful work or changes behavior.
Local Simplification and Canonicalization
Section titled “Local Simplification and Canonicalization”The cheapest passes remove algebraic redundancy:
This includes deleting identities, canceling adjacent self-inverse gates, fusing rotations about the same Pauli axis, normalizing angles, removing dead classical calculations, and canonicalizing equivalent gate spellings. A commutation pass can bring matching operations together before the local rule fires.
Angle normalization must respect the equivalence contract. With
one has
Reduction modulo is valid up to global phase for an isolated unitary, but modulo is needed for literal matrix equality. The sign can become a relative phase after adding a control. Canonicalization is therefore not semantics-free bookkeeping.
Local rules should also preserve parameter domains. The exact symbolic rule
holds for every allowed . Replacing a rotation because its angle is small at one sampled parameter value is approximate specialization, not symbolic optimization.
Worked example: merge two parity rotations
Section titled “Worked example: merge two parity rotations”Writing , the standard decomposition
can arise twice in succession. Expanding both blocks gives
The middle CNOTs cancel and the rotations fuse. The result reduces four entanglers to two, or to zero when is an allowed identity angle. This example also shows why pass order matters: expanding structured operations can expose a cancellation, while expanding too early can hide larger algebraic structure.
Peepholes, Templates, and Superoptimization
Section titled “Peepholes, Templates, and Superoptimization”A peephole optimizer selects a small subcircuit , computes or recognizes its semantics, and replaces it by a cheaper equivalent . Exact matrix comparison is practical only for small wire counts, but a window can also be classified by a Clifford tableau, binary linear map, phase polynomial, or another compact invariant.
Template methods begin with an identity
If commutation and pattern matching identify a costly product from this template inside the circuit, the relation implies
Replacing by is profitable when the chosen cost of is smaller and every side condition is satisfied. Templates can reach beyond adjacent cancellation while remaining interpretable and fast.
Window resynthesis searches more broadly. It may enumerate circuits up to a size bound, use meet-in-the-middle tables, invoke a canonical synthesizer, or apply numerical optimization to a parameterized ansatz. Superoptimization extends this idea by generating equivalence classes or candidate identities and searching them under a cost model. Because exhaustive search grows rapidly, practical systems restrict width, gate count, parameter form, or search time.
Rewrite direction needs a termination policy. If every accepted exact rule strictly decreases a well-founded lexicographic score, repeated rewriting terminates, though possibly at a poor local minimum. Bidirectional rules such as commutation can loop. Equality-saturation systems avoid committing to one direction by storing many equivalent forms and performing cost extraction later, but their equivalence graph can itself grow beyond practical limits.
Exploit Algebraic Structure
Section titled “Exploit Algebraic Structure”Local matrices are not the only useful semantics. Different circuit fragments admit different compact descriptions.
| Fragment | Compact semantics | Typical optimization |
|---|---|---|
| CNOT network | invertible binary matrix | block resynthesis and Gaussian-style elimination |
| CNOT plus diagonal phases | linear map plus phase polynomial | parity-term fusion and resynthesis |
| Clifford circuit | symplectic map and phase data | tableau reduction and Clifford resynthesis |
| Pauli rotations | labeled Pauli strings and angles | commuting sets, phase gadgets, shared parity networks |
| reversible classical circuit | basis-state permutation or Boolean map | templates and reversible-logic resynthesis |
| general bounded-width window | dense unitary or channel | exact or approximate peephole search |
| ZX diagram | typed open graph with rewrite semantics | graph simplification and circuit extraction |
The optimizer should move into a representation whose equality test and rewrite rules match the fragment, then return with a checked extraction procedure. No one representation makes every optimization easy.
Linear reversible blocks
Section titled “Linear reversible blocks”A CNOT-only circuit acts on computational-basis strings as
Each CNOT is an elementary row operation over . An optimizer can accumulate a maximal block into and resynthesize the map rather than preserve an inefficient historical sequence. On unrestricted connectivity, the worst-case CNOT count can be made , asymptotically improving ordinary one-row-at-a-time elimination. Connectivity-constrained synthesis belongs with mapping and routing; the algebraic target should be retained for that later pass.
Phase polynomials
Section titled “Phase polynomials”Let
A circuit over CNOT and these phase gates acts as
where and
Each is a parity function, embedded as or in the phase. A CNOT updates the live parity labels; a phase gate adds its angle to the label currently carried by its wire. Terms with the same can be combined, and a coefficient equal to zero modulo disappears.
This representation sees cancellations separated by long CNOT networks. It also distinguishes the final linear transformation from the phase function , so either part can be resynthesized. Since
converting between phase-gate and conventions requires tracking global phase whenever the context is phase-sensitive.
For Clifford+ circuits, coefficients are multiples of . Odd coefficients represent non-Clifford phase terms, while even coefficients are Clifford phases. Reducing the number or layering of odd terms can lower count or depth. For important diagonal families, -count reduction is related to decoding punctured Reed–Muller codes; more generally, common gate-count and depth optimization problems are NP-hard. Efficient heuristics and restricted exact algorithms are therefore central, not temporary substitutes for a known universal polynomial-time optimizer.
Clifford and Pauli structure
Section titled “Clifford and Pauli structure”A Clifford circuit maps Paulis to Paulis under conjugation, so its action can be stored as a binary symplectic transformation plus phase data. Large Clifford windows can then be compared and resynthesized without constructing a matrix. Stabilizer Formalism develops the tableau and symplectic representation.
Circuits dominated by Pauli rotations benefit from grouping commuting strings, combining identical rotations, and sharing parity-computation networks. A phase gadget represents a parity-dependent phase independently of one particular CNOT ladder. Keeping the gadget abstract allows a later extractor to trade entangler count against depth or connectivity.
ZX-calculus
Section titled “ZX-calculus”The ZX-calculus translates circuits into open diagrams whose sound rewrite rules can expose nonlocal Clifford and phase structure. Graph transformations such as local complementation and pivoting can simplify a diagram before a new circuit is extracted. This can reveal cancellations hidden from fixed-width peepholes.
Three claims must remain separate:
- every applied graphical rewrite is sound in the declared fragment;
- the chosen rewrite strategy reaches a useful reduced diagram;
- circuit extraction succeeds and improves the selected circuit cost.
Soundness does not imply that a heuristic finds a global optimum, and a smaller diagram need not extract to a cheaper circuit under every gate alphabet. Completeness is also fragment-dependent. The certificate should record the rule set, scalar convention, extraction algorithm, and post- extraction cost.
Depth Reduction and Scheduling
Section titled “Depth Reduction and Scheduling”For a dependency graph with operation duration , an abstract critical path is
Taking gives layer depth; assigning gate-class durations gives a weighted depth. An as-soon-as-possible schedule supplies an upper bound after resource conflicts are included, while an as-late-as-possible schedule reveals slack. Commutation can remove avoidable order edges, and resynthesis can change the graph itself.
Count and depth can move in opposite directions. Parallelizing parity computations may require more ancillas or CNOTs. Minimizing depth can leave count unchanged and increase Clifford work. A target-independent depth claim should therefore state:
- operation durations or the unit-layer convention;
- which operations conflict for resources;
- whether classical feedforward has zero or nonzero latency;
- whether connectivity and movement are included;
- whether ancillas may be introduced.
Physical placement, routing, crosstalk, calibrated durations, and controller timing can reverse an abstract schedule. Those refinements belong to their later canonical pages.
Pass Ordering and Fixed Points
Section titled “Pass Ordering and Fixed Points”Optimization passes generally do not commute:
Rotation fusion before finite-alphabet synthesis can approximate one angle instead of two. Decomposition can expose adjacent inverses. Local cleanup can make a larger resynthesis window recognizable. Routing can introduce SWAP or direction-correction gates that create new local cancellations. Conversely, premature decomposition can erase a high-level operation that a structure-aware pass would optimize better.
A robust pipeline commonly alternates:
- validation and canonical normalization;
- inexpensive local cleanup;
- dependency and commutation analysis;
- structure-specific block optimization;
- bounded global or window search;
- cleanup and cost extraction;
- translation validation.
Repeating passes to a fixed point is useful only with a stop rule. Record the iteration bound, score history, timeout, random seed, and whether the final candidate is deterministic. A fixed point under one pass set means only that those passes found no accepted rewrite; it is not a proof of global optimality.
Parameterized and Approximate Circuits
Section titled “Parameterized and Approximate Circuits”Parameterized circuits require identities valid on the declared parameter domain. An optimizer may:
- combine affine angle expressions exactly;
- propagate known periodicities and constants;
- factor shared parameter expressions;
- use commutation conditions independent of parameter values;
- specialize only after parameters are bound.
Numerical sampling of several parameter values is not proof of a symbolic identity. Branch cuts in Euler angles, floating-point comparisons near multiples of , and assumptions such as should be explicit. For variational workloads, an approximate rewrite can also alter derivatives and the objective landscape even when sampled unitaries are close.
Approximate block replacement should return a residual in a composable metric. If unitary replacements satisfy
then a telescoping argument gives the conservative bound
This bound may be loose, but it prevents an optimizer from spending the same global tolerance independently in every window. Gate Decomposition owns synthesis-error allocation and finite-alphabet approximation.
Measurements, Resets, and Classical Control
Section titled “Measurements, Resets, and Classical Control”Dynamic operations break many unitary intuitions:
- measurement is not invertible, so it has no cancellation partner;
- reset discards input information and is not an inverse of preparation;
- moving a gate through measurement can change the measured observable;
- branch-local rewrites must preserve classical labels and branch conditions;
- dead-branch elimination needs a proof that the condition is unreachable;
- moving work across a feedback boundary can change latency and noise exposure.
For example, applying before a computational-basis measurement flips the outcome distribution. Moving it after measurement is equivalent only if the classical bit is also negated and no retained post-measurement quantum state requires the original action. This is an instrument rewrite with a classical relabeling, not ordinary gate commutation.
Compiler IRs should expose quantum and classical dependencies, measurement result ownership, reset postconditions, and semantic fences. Optimizing only a drawn unitary skeleton can silently invalidate an adaptive program.
Verification and Certificates
Section titled “Verification and Certificates”Optimization is unusually well suited to translation validation: verify the particular input–output pair produced by a pass, even when the optimizer itself is not formally proved correct.
For static unitary circuits on modest width, form the miter
Exact phase-insensitive equivalence requires . Full matrices are exponential, so larger circuits need structure-aware methods: Clifford tableaux, phase polynomials, path sums, decision diagrams, tensor contraction, ZX reduction, symbolic algebra, or proof-assistant-verified rewrite rules. Each method has a scope and can return “unknown” without showing non-equivalence.
Use layered evidence:
- prove each primitive rewrite once, including side conditions;
- validate IR invariants before and after every pass;
- check compact semantic invariants for structured fragments;
- translation-validate the actual emitted circuit;
- use full matrices only for bounded width;
- run randomized state tests as diagnostics, not proofs;
- retain an independently checkable certificate where practical.
A useful optimization record includes:
| Field | Minimum content |
|---|---|
| input | circuit hash, IR version, gate definitions, parameter bindings |
| equivalence | unitary, phase-insensitive, subspace, channel, instrument, or approximate |
| resources | ancilla initialization and restoration, classical outputs, fences |
| objective | full cost vector, ordering policy, durations and connectivity assumptions |
| passes | ordered pass list, parameters, iterations, implementation versions |
| search | timeout, random seed, candidate budget, optimality claim if any |
| validation | checker, semantic representation, precision, residual, result |
| output | circuit hash, resource counts, depth model, provenance map |
Formal verification does not certify the physical device. It certifies a relation between modeled artifacts. Calibration drift, crosstalk, leakage, controller faults, and model mismatch require separate experimental evidence.
Common Mistakes
Section titled “Common Mistakes”- Reporting “gate-count reduction” without naming the gate alphabet and counted gate classes.
- Treating a lower count as a lower depth, duration, or logical cost.
- Canceling gates that are adjacent in text but separated by a classical, timing, or resource dependency.
- Commuting through a measurement or reset using a unitary identity alone.
- Reducing angles modulo when literal phase matters.
- Using a floating-point near-zero test as an exact symbolic rewrite.
- Dropping or borrowing ancillas without preserving their state contract.
- Assuming a local fixed point is globally optimal.
- Comparing optimizer benchmarks with different decomposition, routing, timeout, or verification settings.
- Running an optimization pass without updating provenance and approximation budgets.
- Treating successful random simulation as an equivalence proof.
- Calling a hardware-noise score “fidelity” without a calibrated model and validation data.
Exercises
Section titled “Exercises”1. Fuse rotations with the correct period
Section titled “1. Fuse rotations with the correct period”Simplify and state when an angle may be reduced modulo rather than .
Solution
The generators are identical, so their exponentials commute and combine:
Because , reduction modulo preserves the operation only up to global phase. It is valid when the equivalence contract quotients global phase and the block cannot later be controlled or placed in a phase-sensitive context. Literal matrix equality uses period .
2. Test Pauli commutation
Section titled “2. Test Pauli commutation”Determine whether commutes with , and decide whether their rotations may be reordered.
Solution
On qubit , and anticommute. On qubit , the two operators commute. On qubit , is identity, and on qubit , is identity. There is one anticommuting overlap, so
The strings do not commute, and generic rotations and cannot be reordered. Special angles can produce accidental identities, but they require a separate parameter-aware proof.
3. Optimize two ZZ gadgets
Section titled “3. Optimize two ZZ gadgets”Starting from the decomposed circuit for , count the CNOTs before and after local optimization.
Solution
Each gadget initially contains two CNOTs, so the concatenation has four. The two middle CNOTs cancel. The remaining rotations fuse:
The two outer CNOTs then become adjacent and cancel. The optimized circuit is the identity with zero CNOTs. The result is exact, including global phase, because .
4. Extract a phase polynomial
Section titled “4. Extract a phase polynomial”On basis input , perform in temporal order , , , and . Find and .
Solution
The first phase contributes . After the first CNOT, wire carries , so the second phase contributes . The final CNOT restores the input bits. Hence
Any later phase term on the same parity can be fused by adding its coefficient modulo .
5. Use an identity template
Section titled “5. Use an identity template”Suppose an exact template factors as . Prove the replacement and give one condition beyond gate count that can make it undesirable.
Solution
From , right multiplication by gives
The rewrite is therefore exact when the matched gates satisfy every template side condition. Even if uses fewer gates, it may have greater critical-path depth, more expensive gate classes, unsupported inverse gates, additional ancillas, or worse connectivity. Acceptance must use the declared cost and legality contract.
6. Compute weighted depth
Section titled “6. Compute weighted depth”A dependency graph has edges , , , and . Durations are
Find the critical-path depth.
Solution
The maximal path weights are
Thus . There are two critical paths. Removing only one of them does not necessarily reduce total depth.
7. Explain a pass-order advantage
Section titled “7. Explain a pass-order advantage”Two adjacent rotations and must be approximated over a finite gate alphabet. Why can fusion before synthesis outperform synthesis before fusion?
Solution
Fusion gives one exact symbolic target,
which needs one approximation and one allocated error budget. Synthesizing the two rotations separately produces two gate words, may spend two error allocations, and can hide the exact angle relation from a later local pass. The later pass might still cancel parts of the words, but it is not guaranteed to recover a near-optimal approximation of the combined angle. This is why high-level algebraic optimization should normally precede irreversible finite-alphabet lowering.
8. Preserve a measurement rewrite
Section titled “8. Preserve a measurement rewrite”A circuit applies , measures , and returns the classical result. Can the be moved after measurement while preserving the declared output?
Solution
Not without changing the classical processing. For input , applying before measurement exchanges the probabilities of outcomes and . Measurement first gives the unexchanged label. An equivalent rewrite may measure first and then return the complemented classical bit .
If the post-measurement quantum state is retained, the optimizer must also check its required convention. The rewrite is an equivalence of instruments with relabeled classical output, not a commutation of two unitary gates.
9. Design an optimization certificate
Section titled “9. Design an optimization certificate”What evidence is needed to support the claim “this pass reduced the circuit depth by without changing the computation”?
Solution
Record the exact input and output circuits or hashes; IR, operation, tensor, parameter, and phase conventions; the equivalence notion; ancilla, measurement, and output contracts; the pass sequence and versions; and the translation-validation method and result.
For the depth claim, record the old and new dependency graphs or reproducible schedules, gate durations, resource conflicts, connectivity assumptions, classical latency, and rounding convention. State whether refers to unit layers, weighted logical duration, or a target-native schedule. Without those details, neither the semantic nor the performance claim is auditable.
Research Status
Section titled “Research Status”Inverse cancellation, rotation fusion, Clifford/tableau methods, linear reversible synthesis, and phase-polynomial reasoning are established tools. Template, peephole, graph-rewrite, and verified-rewrite systems are also well-developed, but their achieved cost depends strongly on gate alphabet, benchmark family, pass order, and resource model.
Active work includes scalable equivalence checking, automatically discovered rewrite systems, equality saturation, topology-aware phase-polynomial extraction, measurement-assisted fault-tolerant optimization, machine-guided search, and multi-objective optimization across logical and physical layers. Recent learned systems have improved selected -count benchmarks, but this is evidence for a method on those benchmarks, not a proof of universal superiority or global optimality. Hardness results for common gate-count and depth objectives explain why restricted exact solvers and heuristic search coexist.
Further Connections
Section titled “Further Connections”- Quantum Software Stack places optimization in the full path from application intent to executable and evidence.
- Circuit Intermediate Representations defines the typed dependencies, semantic fences, gate profiles, and provenance that optimization must preserve.
- Gate Decomposition develops structure-aware synthesis, finite-alphabet approximation, and synthesis certificates.
- Resource Estimation Tools turns verified changes in count, depth, liveness, and temporal resource demand into layered architecture estimates and Pareto comparisons.
- Qubit Mapping and Routing places program qubits on constrained targets, legalizes nonlocal interactions, tracks changing identities, and verifies the final route.
- Error-Aware Compilation supplies dated, uncertainty-qualified physical objectives for choosing among semantically valid optimized circuits.
- Circuit Model fixes register order, dynamic-circuit semantics, and logical, compiled, and physical resource levels.
- Universal Gate Sets explains exact and approximate gate alphabets and why Clifford+ cost matters.
- Stabilizer Formalism develops Pauli, Clifford, symplectic, and tableau machinery used by structured optimizers.
- Control, Readout, and Calibration explains why calibrated duration and error information can change a circuit-level ranking.
- Quantum Gates Formula Card and Quantum Gates Reference Table provide compact convention checks for rewrite rules.
References
Section titled “References”- D. Maslov, G. W. Dueck, D. M. Miller, and C. Negrevergne, “Quantum circuit simplification and level compaction,” IEEE Transactions on Computer-Aided Design of Integrated Circuits and Systems 27, 436–444 (2008), doi:10.1109/TCAD.2007.911334.
- Y. Nam, N. J. Ross, Y. Su, A. M. Childs, and D. Maslov, “Automated optimization of large quantum circuits with continuous parameters,” npj Quantum Information 4, 23 (2018), doi:10.1038/s41534-018-0072-4.
- K. N. Patel, I. L. Markov, and J. P. Hayes, “Efficient synthesis of linear reversible circuits,” Quantum Information and Computation 8, 282–294 (2008), arXiv:quant-ph/0302002.
- M. Amy, D. Maslov, and M. Mosca, “Polynomial-time -depth optimization of Clifford+ circuits via matroid partitioning,” IEEE Transactions on Computer-Aided Design of Integrated Circuits and Systems 33, 1476–1489 (2014), doi:10.1109/TCAD.2014.2341953.
- M. Amy and M. Mosca, “-count optimization and Reed–Muller codes,” IEEE Transactions on Information Theory 65, 4771–4784 (2019), doi:10.1109/TIT.2019.2906374.
- L. E. Heyfron and E. T. Campbell, “An efficient quantum compiler that reduces count,” Quantum Science and Technology 4, 015004 (2018), doi:10.1088/2058-9565/aad604.
- R. Duncan, A. Kissinger, S. Perdrix, and J. van de Wetering, “Graph-theoretic simplification of quantum circuits with the ZX-calculus,” Quantum 4, 279 (2020), doi:10.22331/q-2020-06-04-279.
- A. Kissinger and J. van de Wetering, “Reducing -count with the ZX-calculus,” Physical Review A 102, 022406 (2020), doi:10.1103/PhysRevA.102.022406.
- A. Cowtan, S. Dilkes, R. Duncan, W. Simmons, and S. Sivarajah, “Phase gadget synthesis for shallow circuits,” in Proceedings of QPL 2019, EPTCS 318, 213–228 (2020), doi:10.4204/EPTCS.318.13.
- K. Hietala, R. Rand, S.-H. Hung, X. Wu, and M. Hicks, “A verified optimizer for Quantum circuits,” Proceedings of the ACM on Programming Languages 5, POPL, Article 37 (2021), doi:10.1145/3434318.
- M. Xu et al., “Quartz: Superoptimization of quantum circuits,” in Proceedings of PLDI 2022, 625–640 (2022), doi:10.1145/3519939.3523433.
- S. Yamashita and I. L. Markov, “Fast equivalence-checking for quantum circuits,” in Proceedings of NQCC 2010, 23–28 (2010), arXiv:0909.4119.
- T. Peham, L. Burgholzer, and R. Wille, “Equivalence checking of quantum circuits with the ZX-calculus,” IEEE Journal on Emerging and Selected Topics in Circuits and Systems 12, 662–675 (2022), doi:10.1109/JETCAS.2022.3202204.
- J. van de Wetering and M. Amy, “Optimising quantum circuits is generally hard,” arXiv:2310.05958, version 3 (2024), arXiv:2310.05958.
- F. J. R. Ruiz et al., “Quantum circuit optimization with AlphaTensor,” Nature Machine Intelligence 7, 374–385 (2025), doi:10.1038/s42256-025-01001-1.
- M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information, 10th anniversary ed., Cambridge University Press (2010), doi:10.1017/CBO9780511976667.