---
title: Gaussian Stabilizer Group in Quantum Theory
url: https://www.emergentmind.com/topics/gaussian-stabilizer-group
type: topic
---

# Gaussian Stabilizer Group in Quantum Theory

The expression **Gaussian Stabilizer Group** appears in two closely related but technically distinct settings. In discrete-variable quantum information, the phrase denotes the role played by the **Clifford group** as the discrete analogue of the continuous-variable Gaussian unitary group: Clifford unitaries normalize the Pauli/Weyl group, implement symplectic linear maps on finite phase space, preserve stabilizer states, and underlie a discrete-variable convolution theory in which stabilizer states behave as “discrete quantum Gaussians” [2302.08423]. In the square-GKP setting, the same expression refers to the subgroup of **Gaussian unitaries** that act as encoded identities on the logical code space; this group is characterized by affine symplectic automorphisms of the dual lattice and is used to search over logically equivalent physical implementations of Clifford circuits under noise [2509.12502]. In both usages, the core structure is symplectic phase-space dynamics constrained by stabilizer preservation, but the ambient Hilbert spaces, group definitions, and operational roles differ.

## 1. Terminological scope and conceptual setting

In the discrete-variable formulation, the guiding claim is that **stabilizer states play a role in DV quantum systems similar to the role Gaussian states play in continuous-variable systems**, which motivates the designation “discrete quantum Gaussians” for stabilizer states [2302.08423]. From this viewpoint, the natural group of “Gaussian operations” is the Clifford group, because it preserves stabilizer states and acts symplectically on finite phase space.

In the GKP formulation, the starting point is different. For grid codes, the unitary stabilizer group consists of those unitaries that act as encoded identities on the logical subspace, and the **Gaussian stabilizer group** is obtained by restricting this unitary stabilizer group to Gaussian unitaries [2509.12502]. Here the question is not whether Gaussian states are replaced by stabilizer states, but which Gaussian transformations leave the logical action trivial while possibly changing the physical implementation.

These usages are structurally aligned because both are formulated in terms of Heisenberg–Weyl covariance and symplectic transformations. They are nevertheless not identical. In the DV setting, “Gaussian stabilizer group” is effectively a group-theoretic interpretation of the Clifford group as the discrete Gaussian group. In the GKP setting, it is a subgroup of Gaussian unitaries defined by trivial encoded action on a bosonic stabilizer code.

## 2. Finite-dimensional formulation: Weyl systems, stabilizers, and the Clifford group

For prime local dimension \(d\) and \(n\) qudits, the finite phase space is
\[
V^n = \mathbb{Z}_d^n \times \mathbb{Z}_d^n.
\]
The single-qudit Pauli operators satisfy
\[
X|k\rangle = |k+1\rangle,\qquad Z|k\rangle = \chi(k)|k\rangle,\qquad \chi(k)=\exp(2\pi i k/d).
\]
The Weyl operators are
\[
w(p,q)=\chi(-2^{-1}pq)Z^pX^q \quad \text{for odd } d,
\]
and
\[
w(p,q)=i^{-pq}Z^pX^q \quad \text{for } d=2,
\]
with symplectic inner product
\[
\{(p,q),(p',q')\}_s = pq' - qp'.
\]
For \(n\) qudits,
\[
w(\vec p,\vec q)=\prod_k w(p_k,q_k).
\]

A pure stabilizer state is a common eigenvector of an abelian subgroup \(S\) of the Weyl group of size \(d^n\), with \(n\) commuting generators. General stabilizer states are convex mixtures of pure stabilizers. Equivalently, the stabilizer algebra \(C^*(S)\) is abelian and has minimal projections; normalized minimal stabilizer projections (MSPS) are the elementary stabilizer states [2302.08423].

An \(n\)-qudit unitary \(U\) is Clifford if it normalizes the Weyl operators:
\[
U\, w(\vec p,\vec q)\, U^\dagger \propto w(M(\vec p,\vec q))
\]
for some symplectic linear transformation \(M\) on the discrete phase space. Clifford unitaries preserve stabilizer states. This normalizer property is the decisive group-theoretic reason the Clifford group functions as the DV analogue of the Gaussian unitary group. The relevant phase-space action is linear and symplectic, and the preserved family of states is the stabilizer family, which the paper identifies as the discrete analogue of Gaussian states [2302.08423].

## 3. Characteristic functions, convolution, and discrete Gaussian behavior

The discrete-variable characteristic function of a state \(\rho\) is
\[
\Xi_\rho(\vec p,\vec q)=\mathrm{Tr}[\rho\, w(-\vec p,-\vec q)],
\]
and the Weyl expansion is
\[
\rho = \frac{1}{d^n}\sum_{(\vec p,\vec q)\in V^n}\Xi_\rho(\vec p,\vec q)\, w(\vec p,\vec q).
\]
For MSPS, and more generally stabilizer states, \(|\Xi_\rho|\in\{0,1\}\) with symplectic-abelian support. This support structure is central to the discrete Gaussian analogy.

The paper introduces a **DV quantum convolution** using a Clifford unitary \(U(G)\) associated with an invertible \(2\times 2\) matrix \(G\) over \(\mathbb{Z}_d\):
\[
\rho \,\boxplus\, \sigma := \mathrm{Tr}_B\!\left[U(G)(\rho\otimes \sigma)U(G)^\dagger\right].
\]
Its action on characteristic functions is multiplicative:
\[
\Xi_{\rho \boxplus \sigma}(\vec p,\vec q)
=
\Xi_\rho(N g_{11}\vec p, g_{00}\vec q)\,
\Xi_\sigma(-N g_{10}\vec p, g_{01}\vec q),
\]
where \(N=(\det G)^{-1}\) in \(\mathbb{Z}_d\). Special cases include a discrete beam splitter, discrete amplifier, and, for odd \(d\), a Hadamard convolution [2302.08423].

This formalism yields a closure theorem: if \(\rho\) and \(\sigma\) are stabilizer states, then \(\rho \boxplus \sigma\) is again a stabilizer state. For MSPS, the proof uses the fact that \(|\Xi_\rho|\) and \(|\Xi_\sigma|\) take only the values \(0\) and \(1\), and that the support of the product remains abelian under the symplectic commutation criterion. Convexity extends the result to general stabilizer states.

The same framework supports several extremality results. The **mean state** \(M(\rho)\) is defined by
\[
\Xi_{M(\rho)}(\vec p,\vec q)=
\begin{cases}
\Xi_\rho(\vec p,\vec q), & |\Xi_\rho|=1,\\
0, & \text{otherwise}.
\end{cases}
\]
The paper proves that \(M(\rho)\) is an MSPS determined by the abelian subgroup
\[
S=\{(\vec p,\vec q): |\Xi_\rho(\vec p,\vec q)|=1\},
\]
and that it obeys Clifford covariance:
\[
U M(\rho) U^\dagger = M(U\rho U^\dagger).
\]
For Rényi entropies,
\[
D_\alpha(\rho\|M(\rho)) = H_\alpha(M(\rho)) - H_\alpha(\rho),
\qquad
H_\alpha(M(\rho)) \ge H_\alpha(\rho),
\]
with equality only at \(M(\rho)\). Thus MSPS maximize all Schur-concave spectral functionals among states with the same mean state up to Clifford conjugation [2302.08423].

For “positive” \(G\), the convolution obeys majorization:
\[
\lambda(\rho\boxplus \sigma)\prec \lambda(\rho),
\qquad
\lambda(\rho\boxplus \sigma)\prec \lambda(\sigma),
\]
hence
\[
H_\alpha(\rho\boxplus \sigma)\ge \max\{H_\alpha(\rho),H_\alpha(\sigma)\}.
\]
Iterated convolution gives a discrete-variable “second law for quantum convolution”:
\[
H_\alpha(\boxplus^{N+1}\rho)\ge H_\alpha(\boxplus^N\rho).
\]

The paper also proves a Fisher-information inequality for positive \(G\):
\[
J(\rho\boxplus \sigma)\le \min\{J(\rho),J(\sigma)\},
\]
where
\[
J(\rho;H)=\mathrm{Tr}[\rho[H,[H,\log\rho]]].
\]
A DV de Bruijn identity is established:
\[
\left.\frac{d}{dt}\right|_{t=0} H(e^{t\mathcal L}(\rho)) = \frac{1}{4}J(\rho).
\]

The culmination is a **DV quantum central limit theorem**. For a zero-mean \(n\)-qudit state \(\rho\), repeated beam-splitter convolution converges exponentially in Hilbert–Schmidt norm to the mean state:
\[
\|\boxplus^N \rho - M(\rho)\|_2
\le
(1-\mathrm{MG}(\rho))^N \,\|\rho-M(\rho)\|_2,
\]
where the magic gap is
\[
\mathrm{MG}(\rho)=1-\max\{|\Xi_\rho(\vec p,\vec q)|:(\vec p,\vec q)\in \mathrm{Supp}(\Xi_\rho),\, |\Xi_\rho(\vec p,\vec q)|\neq 1\}.
\]
The limit state is an MSPS and therefore a stabilizer state. In this precise sense, stabilizer states are attractors under DV convolution in the same way Gaussian states are attractors in CV central limit theory [2302.08423].

## 4. Channel-level extension and the “Gaussian” interpretation of the Clifford group

The discrete Gaussian analogy extends from states to channels. A channel \(\Lambda:\mathcal H_A\to \mathcal H_{A'}\) is represented by its Choi state
\[
J_\Lambda=(\mathrm{id}_A\otimes \Lambda)(|\Phi\rangle\langle \Phi|),
\qquad
|\Phi\rangle=\frac{1}{\sqrt{d^n}}\sum_{\vec j} |\vec j\rangle_A |\vec j\rangle_{A'}.
\]
Channel convolution is defined by convolving Choi states:
\[
J_{\Lambda_1\boxplus \Lambda_2}:=J_{\Lambda_1}\boxplus J_{\Lambda_2}.
\]
Equivalently,
\[
\Lambda_1\boxplus \Lambda_2 = \mathcal E \circ \Lambda_1\Lambda_2 \circ \mathcal E^{-1},
\]
where \(\mathcal E\) is the DV convolutional channel induced by \(U(G)\) [2302.08423].

A stabilizer channel maps stabilizer states to stabilizer states. Because \(U(G)\) is Clifford and convolution preserves stabilizer structure, the convolution of stabilizer channels is stabilizer, and the mean channel \(M(\Lambda)\), defined through the Choi state, is likewise a stabilizer channel. This is the channel-level counterpart of the statement that stabilizer states are the discrete Gaussians.

The paper proves channel entropy extremality using the Gour–Wilde channel entropy
\[
H_\alpha(\Lambda)= n\log d - D_\alpha(\Lambda\|\mathcal R).
\]
For \(\alpha\in[1/2,\infty]\),
\[
H_\alpha(M(\Lambda))\ge H_\alpha(\Lambda),
\]
with equality only for \(\Lambda\) in the abelian algebra determined by the stabilizer support of \(J_\Lambda\). For positive \(G\),
\[
H_\alpha(\Lambda_1\boxplus \Lambda_2)\ge \max\{H_\alpha(\Lambda_1),H_\alpha(\Lambda_2)\},
\]
and iterated convolution yields the channel version of the second law:
\[
H_\alpha(\boxplus^{N+1}\Lambda)\ge H_\alpha(\boxplus^N \Lambda).
\]

The convolutional channel \(\mathcal E\) achieves minimal output entropy exactly on pure stabilizer inputs and achieves maximal Holevo capacity iff the state is stabilizer. These results are presented as direct analogues of Gaussian extremality for continuous-variable channels. The resulting interpretation is explicit: the Clifford group is the DV “Gaussian stabilizer group” because it normalizes the Weyl group, implements the convolutional mixing unitaries, preserves stabilizer states and stabilizer channels, and organizes the entropy, Fisher-information, and central-limit structure of the theory [2302.08423].

## 5. Square-GKP Gaussian stabilizer group

For square-GKP codes, the setting is bosonic rather than finite-dimensional. Single-mode phase space is coordinated by
\[
\hat r=(\hat q,\hat p)^T,
\qquad
[\hat q,\hat p]=i,
\]
and the displacement operator is
\[
D(\xi)=\exp(i\xi^T J \hat r),
\qquad
J=
\begin{bmatrix}
0 & 1\\
-1 & 0
\end{bmatrix}.
\]
In the notation of the paper, translations are
\[
T(v)=\exp(-i l\, v^T\Omega \hat \xi),
\qquad
l=\sqrt{2\pi},
\qquad
\Omega=J.
\]

The standard square-GKP code is defined by a square stabilizer lattice of spacing \(2\sqrt\pi\) in the \((q,p)\) plane, with stabilizer generators
\[
S_q=\exp(i 2\sqrt\pi\, \hat q),
\qquad
S_p=\exp(-i 2\sqrt\pi\, \hat p).
\]
Equivalently,
\[
\Lambda=\{S^T a \mid a\in \mathbb Z^2\},
\qquad
S=\sqrt 2\, I_2,
\qquad
A=S\Omega S^T = 2\Omega_2.
\]
Logical Pauli operators are half-stabilizer displacements:
\[
X_L=\exp(i\sqrt\pi\, \hat q),
\qquad
Z_L=\exp(-i\sqrt\pi\, \hat p),
\]
with displacement vectors
\[
\ell_X=(\sqrt\pi,0)^T,
\qquad
\ell_Z=(0,\sqrt\pi)^T.
\]
The stabilizer lattice is
\[
\Lambda = 2\sqrt\pi\, \mathbb Z^2,
\]
and logical displacements sit at the midpoints of the stabilizer grid [2509.12502].

A Gaussian unitary is parameterized by a symplectic matrix \(S\in Sp(2,\mathbb R)\) and a displacement \(d\in \mathbb R^2\), acting as
\[
U(S,d)\, D(\xi)\, U(S,d)^\dagger = e^{i\phi(\xi)} D(S\xi+d).
\]
For grid codes, the Gaussian stabilizer group is defined as
\[
S_G = S_U \cap \{U_M \mid M\in Sp(2N,\mathbb R)\}.
\]
The paper shows that Gaussian stabilizers act as affine automorphisms of the dual lattice:
\[
S_G \subset \mathrm{Aut}_\infty^S(\Lambda^*) = \mathrm{Aut}^S(\Lambda^*)\ltimes \Lambda^*.
\]
For the single-mode square code, the fixed-point symplectic maps are integral matrices with \(\det=1\), so
\[
S_G \cong \mathrm{Aut}^S(\Lambda^*)\ltimes \Lambda^* \cong SL(2,\mathbb Z)\ltimes \Lambda^*.
\]

The necessary and sufficient conditions for trivial logical action on the square-GKP code are:
\[
S\ell_X \equiv \ell_X \pmod{\Lambda},
\qquad
S\ell_Z \equiv \ell_Z \pmod{\Lambda},
\qquad
d\in \Lambda.
\]
Equivalently, \(S\) maps every Pauli coset grid to itself and the displacement is by a stabilizer vector. In the multimode formulation, with Pauli subgrids \(P_j=\Lambda+p_j\), the condition is
\[
(M-I)p_j\in \Lambda \quad \text{for all Pauli coset representatives } p_j\in \Lambda^*,
\qquad
d\in \Lambda.
\]

The symplectic part admits an algebraic characterization:
\[
M = S^T (X A + I) S^{-T},
\qquad
X\in \mathrm{Mat}_{2N}(\mathbb Z),
\]
subject to
\[
X A X^T - (X-X^T)=0.
\]
In a canonical basis \(A_D=\Omega_2\otimes D\), the solutions form a group under
\[
X_D\circ Y_D = X_D + Y_D + X_D A_D Y_D.
\]
Generators are built from the symmetric basis \(F_{j,k}=E_{j,k}+E_{k,j}\) and skew-symmetric basis \(G_{j,k}=E_{j,k}-E_{k,j}\), with constraints determined by the nonzero entries of \(A_D\). The paper presents this as a complete generating framework for the fixed-point symplectic part and, after adding appropriate translations, for \(S_G\) itself [2509.12502].

## 6. Explicit generators, compilation over \(U_0S_G\), and noise-protection use

In the single-mode square-GKP case, \(S=\sqrt 2 I_2\) and \(A=2\Omega_2\) are already canonical. The paper derives three fundamental integer solutions \(X\) and corresponding symplectic matrices:
\[
X_{H^2}=
\begin{bmatrix}
0 & 1\\
-1 & 0
\end{bmatrix},
\qquad
M_{H^2}=
\begin{bmatrix}
-1 & 0\\
0 & -1
\end{bmatrix},
\]
\[
X_{Q^2}=
\begin{bmatrix}
1 & 0\\
0 & 0
\end{bmatrix},
\qquad
M_{Q^2}=
\begin{bmatrix}
1 & 2\\
0 & 1
\end{bmatrix},
\]
\[
X_{P^2}=
\begin{bmatrix}
0 & 0\\
0 & 1
\end{bmatrix},
\qquad
M_{P^2}=
\begin{bmatrix}
1 & 0\\
-2 & 1
\end{bmatrix}.
\]
These implement the squares of encoded Clifford generators:
\[
\bar H^2=\exp(-i\pi \hat n),\qquad
\bar Q^2=\exp(-i\hat p^2),\qquad
\bar P^2=\exp(-i\hat q^2).
\]
After accounting for the displacement part needed to return Pauli cosets to the stabilizer grid, the single-mode Gaussian stabilizer group is generated by
\[
S_\square = \langle \bar X^2,\bar Z^2,\bar X\bar Q^2,\bar Z\bar P^2,\bar H^2\rangle.
\]

For two modes, the paper identifies four inter-mode symplectic generators
\[
M_1=
\begin{bmatrix}
1&0&0&2\\
0&1&2&0\\
0&0&1&0\\
0&0&0&1
\end{bmatrix},
\quad
M_2=
\begin{bmatrix}
1&0&0&0\\
-2&1&0&0\\
0&0&1&2\\
0&0&0&1
\end{bmatrix},
\]
\[
M_3=
\begin{bmatrix}
1&0&0&0\\
0&1&0&0\\
0&-2&1&0\\
-2&0&0&1
\end{bmatrix},
\quad
M_4=
\begin{bmatrix}
1&-2&0&0\\
0&1&0&0\\
0&0&1&0\\
0&0&2&1
\end{bmatrix},
\]
implementing squares of encoded controlled-Pauli operations:
\[
\bar C(X_1,X_2)^2=\exp(2i\hat p_1\hat p_2),\qquad
\bar C(Z_1,X_2)^2=\exp(2i\hat q_1\hat p_2),
\]
\[
\bar C(Z_1,Z_2)^2=\exp(-2i\hat q_1\hat q_2),
\]
together with the analogous \(\bar C(X_1,Z_2)^2\). For \(N\) modes, the full Gaussian stabilizer group is generated by these inter-mode elements together with the single-mode generators and lattice displacements by stabilizer vectors [2509.12502].

This explicit generating set is used for **noisy Clifford compilation**. Every logical Clifford \(\bar U\) admits multiple Gaussian implementations \(U(S,d)\) differing by multiplication by Gaussian stabilizers \(g\in S_G\). The compiler represents an implementation as
\[
\mathrm{Imp}(\bar U)=U_0 g_n,
\]
where \(U_0\) is a fixed physical representative “closest to identity” and \(g_n\) is a short word in the generators of \(S_G\). The search is therefore over the coset \(U_0S_G\).

The physical criterion is expressed in terms of the envelope of a finite-energy GKP state. Under a Gaussian circuit,
\[
U(S,d): E(\Sigma_E,\mu_E)\to E(S\Sigma_E S^T,\, S\mu_E+d).
\]
The compiler minimizes, at each gate, the envelope displacement magnitude and squeezing through
\[
d_\mu^2(\mu_E)=\mu_E^T\mu_E,
\]
and
\[
d_\Sigma^2(\Sigma_E,\Sigma_0)=
\mathrm{Tr}\!\left[\ln^2\!\left(\sqrt{\Sigma_0^{-1}\Sigma_E\sqrt{\Sigma_0^{-1}}}\right)\right].
\]
The search is restricted to short walks, typically \(n\le 2\), because long words increase squeezing and displacement. Selection is lexicographic: first minimize \(d_\Sigma^2\), then \(d_\mu^2\).

The noise model combines bosonic loss and dephasing. Loss with parameter \(\gamma\) has Kraus operators
\[
L_k=\sqrt{\gamma^k/k!}\,(1-\gamma)^{\hat n/2}\hat a^k,
\qquad
N_L(\gamma)[\rho]=\sum_{k=0}^\infty L_k\rho L_k^\dagger,
\]
and updates Gaussian peaks as
\[
\Sigma \to \gamma\Sigma + (1-\gamma)(1/4\pi)I,
\qquad
\mu \to \sqrt{\gamma}\,\mu.
\]
Dephasing with rate \(\gamma_\phi\) is
\[
N_D(\gamma_\phi)[\rho]
=
\frac{1}{\sqrt{2\pi\gamma_\phi}}
\int_{\mathbb R} d\phi\,
e^{-\phi^2/(2\gamma_\phi)}
e^{i\phi \hat n}\rho e^{-i\phi \hat n},
\]
and induces angular smearing. Because the channels commute, the average magnitude of displacement experienced by a phase-space feature at radius \(r\) per application is
\[
\Delta \xi(r)= r\sqrt{\gamma + (2/\pi)\gamma_\phi^2}.
\]
When \(\Delta\xi(r)\) approaches the logical spacing \(\sqrt\pi\), logical shifts proliferate, so the compiler suppresses envelope radius and squeezing.

The numerical demonstration uses logical randomized benchmarking. The paper reports that for \(1\) and \(3\) logical qubits at \(\gamma=0.999\) and \(\gamma=0.9999\), the GS compiler consistently yields higher survival probabilities at all lengths than both a constant compiler and a random-walk compiler. Lifetime is defined as the number of gates \(t\) needed to reach a \(1/e\) drop from the initial baseline under the Gaussian fit
\[
P(x)=A e^{-(a x^2 + b x)} + B.
\]
The lifetime improvement factor \(t_{GS}/t_{RW}\) increases as loss weakens; at very high loss, all compilers perform poorly and gains vanish. The same advantage appears in two-qubit-only randomized benchmarking over the group \(\langle \bar C X,\bar C Z\rangle\), which the paper identifies as relevant for magic-state-based computations where single-qubit Cliffords are tracked in a frame and only two-qubit gates are physically executed. Diagnostic metrics show that the maximum \(d_\mu^2\) and \(d_\Sigma^2\) attained along random sequences remain significantly lower under the GS compiler than under the other strategies [2509.12502].

A plausible implication is that, in the GKP setting, the Gaussian stabilizer group is not merely a structural symmetry group but also a compilation resource: it parameterizes physically distinct Gaussian realizations of the same logical Clifford operation and thereby exposes an optimization space adapted to Gaussian noise. In the DV setting, by contrast, the emphasis is foundational: the Clifford group is “Gaussian” because it preserves the discrete stabilizer sector and supports convolutional, entropic, Fisher-information, and central-limit phenomena directly analogous to continuous-variable Gaussian theory.

Source: https://www.emergentmind.com/topics/gaussian-stabilizer-group