---
title: GHZ-Preserving Operations
url: https://www.emergentmind.com/topics/ghz-preserving-operations
type: topic
---

# GHZ-Preserving Operations

A Greenberger–Horne–Zeilinger (GHZ)-preserving operation is a unitary or channel on $n$-qubit systems that permutes the GHZ basis states (up to a global phase), ensuring that evolution never leaves the subspace of states locally Pauli-equivalent to the canonical GHZ state. This property enables efficient simulation and optimization of entanglement distillation protocols for quantum networks. Recent advances provide a complete algebraic classification of such operations, a constant-time simulation method, and circuit optimization frameworks delivering superior performance under realistic noise, with extensions to graph states local-Clifford-equivalent to GHZ states [2510.25854].

## 1. Definition: GHZ States and GHZ-Preserving Operations

The $n$-qubit GHZ state is
\[
|\mathrm{GHZ}_n\rangle = \frac{1}{\sqrt{2}}(|0\rangle^{\otimes n} + |1\rangle^{\otimes n})\,.
\]
The GHZ basis comprises all states of the form
\[
\{ P_1\otimes P_2\otimes\cdots\otimes P_n\,|\,P_i\in\{I,X,Y,Z\}\}\,|\mathrm{GHZ}_n\rangle\,,
\]
i.e., local Pauli operations on each qubit applied to $|\mathrm{GHZ}_n\rangle$.

A unitary $U$ (or quantum channel $\mathcal{E}$) is GHZ-preserving if, for any basis vector $|\mathrm{GHZ}_i\rangle$, 
\[
U\,|\mathrm{GHZ}_i\rangle = e^{i\phi}\,|\mathrm{GHZ}_j\rangle,
\]
for some $j$ and phase $\phi$. In the density-matrix formalism, this reads
\[
\mathcal E(\,|\mathrm{GHZ}_n\rangle\langle\mathrm{GHZ}_n|\,) = |\mathrm{GHZ}_n\rangle\langle\mathrm{GHZ}_n| 
\quad \text{(up to phase for pure states).}
\]
Intuitively, GHZ-preserving gates act exclusively within the GHZ stabilizer structure, only relabeling basis states.

## 2. Algebraic Structure and Enumeration of GHZ-Preserving Gates

GHZ-preserving unitaries, up to phase and local Paulis, can be decomposed into three distinct subgroups [2510.25854, Thm. II.3]:
- **Local Pauli gates:** Arbitrary single-qubit $P_j$ on each of $2n$ qubits (for multiple GHZ copies), absorbing global phases.
- **B group (bilocal):** Eight elements for each pair of nodes, generated by taking the same two-qubit gate $g(f_1\otimes f_2)$ on corresponding pairs, with $g\in\{I, CZ\}$ and $f_1,f_2\in\{I, S\}$.
- **H group (homogeneous):** Six elements generated by applying the identical two-qubit gate at every node, drawn from $\{\mathrm{SWAP},\ \mathit{CNOT}_{12},\ \mathit{CNOT}_{21},\ \mathit{DCX}_{12},\ \mathit{DCX}_{21},\ I\}$.

Any GHZ-preserving unitary may be written as
\[
U = \left(\prod_{j=1}^{2n-2} P_j\right)\left(\prod_{i=1}^{n-1} b_{i,i+1}\right) h,
\]
where $P_j$ are local Paulis, $b_{i,i+1}$ are B-group bilocals, and $h$ is an H-group homogeneous gate.

The total number of distinct (phaseless) GHZ-preserving unitaries is $6\times8^{n-1}$.

| Subgroup | Generator Elements            | Algebraic Structure      |
|----------|-------------------------------|-------------------------|
| B group  | $\{I\otimes I$, $S\otimes I$, $I\otimes S$, $CZ$, etc.$\}$ | $\mathbb{Z}_2^3$         |
| H group  | $\{I$, SWAP, $\mathit{CNOT}_{12}$, $\mathit{CNOT}_{21}$, $\mathit{DCX}_{12}$, ...$\}$ | $D_3$ (dihedral), order $6$ |

B gates enable bilocal error correlation; H gates homogeneously shuffle errors; Paulis absorb phases.

## 3. Simulation Framework: Constant-Time GHZ Stabilizer Updating

Traditional simulation methods require $O(2^n)$ (full state vector) or $O(n)$ (tableau) time per gate. For GHZ-preserving operations, all action is limited to permutations of the "phase bits" in a block-diagonal stabilizer tableau; there is no mixing of basis elements.

The entire configuration of $m$ GHZ basis states is encoded in an $n\times m$ bitstring, tracking the $\pm$ sign of each generator. For two GHZ states ($m=2$), there are $2^{2n}$ possible basis elements.

For each elementary B or H gate $g$, precompute a permutation $\Pi_g:\{0,\ldots,2^{2n}-1\}\to\{0,\ldots,2^{2n}-1\}$. Application of any GHZ-preserving gate $g$ to a state with index $\text{stateIndex}$ then reduces to:
\[
\text{stateIndex} \leftarrow \Pi_g[\text{stateIndex}]
\]
which is $O(1)$ time.

All permutations are computed once at initialization and reused for subsequent simulations, eliminating further bookkeeping. Since gate application never exits the GHZ basis, this method maintains high efficiency irrespective of $n$ [2510.25854].

## 4. Circuit-Level Optimization for GHZ Distillation

Leveraging the $O(1)$ simulation of GHZ-preserving gates and the drastically shrunken gate set ($6\times8^{n-1}$ per layer), a genetic-algorithm optimizer with simulated annealing can efficiently search for short distillation circuits maximizing output fidelity under realistic noise (independent depolarizing noise $p$ per gate, measurement error $\eta$, input fidelity $f_\text{in}$). Hardware constraints include register size $R$ and total raw states $N$.

In a demonstrative scenario (three nodes, $5\!\to\!1$ distillation, $R=3$, $p=0.01$, $\eta=0.01$, $f_\text{in}=0.90$), circuits constructed from H-group gates achieved:
- Higher output fidelity at fixed success probability than standard recurrence methods (e.g., BBPSSW-style).
- Fewer raw states to reach a desired fidelity versus nested protocols.
- Uniform gains across $f_\text{in}$ up to $p,\eta\sim1\%$.

For $P_\mathrm{succ}\ge 0.5$ and $f_\mathrm{in}=0.9$, output fidelity improved by several percent absolute over previous protocols, saving on average one raw input copy [2510.25854].

## 5. GHZ-Preserving Operations for Local-Clifford-Equivalent Graph States

GHZ states correspond to special two-colorable graph states (star graphs with an additional Hadamard on the central node). Any graph state that is local-Clifford (LC)-equivalent to a GHZ state inherits a preserving-gate algebra via conjugation: the B and H generators are mapped through the LC that interconverts the graph and GHZ states.

For example, the 3-qubit complete graph (triangle) with standard generators $\{XZZ, ZXZ, ZZX\}$ admits a GHZ-matching generator set under local basis change ($\{XXX, YYI, IYY\}$). The same enumeration of preserving operations and $O(1)$ simulation extend, with minor relabeling of state indexing.

Full extension to arbitrary graph state topologies is ongoing work [2510.25854].

## 6. Impact on Quantum Networking and Open Resources

By restricting attention to the GHZ basis and its stabilizer group, the search space for circuit synthesis shrinks from superexponential (over the full Clifford group) to $6\times8^{n-1}$ elements per layer. Simulation complexity drops from $O(2^n)$ or $O(n)$ per gate to a constant $O(1)$. These advances enable rapid circuit-level optimization, discovering practical distillation circuits outperforming all known protocols under realistic noise regimes.

These methods are directly relevant for near-term quantum networks conducting multipartite entanglement distillation under resource and noise constraints. An open-source Julia implementation is publicly available at https://github.com/umass-qc/GHZpreserving.jl [2510.25854].

Source: https://www.emergentmind.com/topics/ghz-preserving-operations