Gaussian-Branched Cat States (GCSs)

• Gaussian-Branched Cat States (GCSs) are non-Gaussian quantum states formed from coherent superpositions of Gaussian branches, unifying cat-state physics and bosonic encoding.
• The formalism provides an exact finite-dimensional representation and closed-form dynamics via branch covariance matrices, QRDM, and Riccati-type equations.
• GCS protocols facilitate deterministic state generation, precise control in optical and cavity platforms, and efficient simulation of error-corrected quantum resources.

Gaussian-Branched Cat States (GCSs) are hybrid or bosonic non-Gaussian states built from coherent superpositions of Gaussian branches. In the explicit terminology of the hybrid continuous-variable/discrete-variable formalism, the name denotes the $N$-qubits $n$-modes entangled states arising during “superpositions of Gaussian processes,” fully characterised by superposed phase-space quantities—generalised complex first moments, covariance matrices, and the qubit reduced density matrix (QRDM) (Braccini et al., 1 Oct 2025). In a closely related bosonic formulation, a GCS is any pure or mixed state whose support lies in the span of a finite set of pure Gaussian branches, so that cat states, squeezed-cat states, Gaussian multiplets, and related bosonic encodings all appear as special cases (Centrone et al., 16 Mar 2026).

1. Definition and state-space structure

For a single qubit coupled to $n$ bosonic modes, the hybrid density operator can be decomposed as

$$\hat\varrho = \sum_{j,k\in\{\pm1\}} \hat\varrho_{jk}\otimes |j\rangle\langle k| ,$$

and each matrix element is assigned a branched characteristic function

$$\chi_{jk}(\tilde r)= \exp\!\left[ -\frac14 \tilde r^{\rm T}\sigma_{jk}\tilde r +i\,\tilde r^{\rm T}r_{jk} +r^{(0)}_{jk} \right].$$

For $N$ qubits, the branch labels become bit strings $J,K\in\mathcal M$, and the general GCS ansatz is

$$\chi_{JK}(\tilde r)= \exp\!\left[ -\frac14 \tilde r^{\rm T}\sigma_{JK}\tilde r +i\,\tilde r^{\rm T}r_{JK} -\mathcal C_{JK}+i\phi_{JK} \right].$$

Diagonal branches $J=K$ are physical Gaussian states with real covariance matrices and first moments; off-diagonal branches encode interference and may have complex phase-space data. The QRDM is recovered at $\tilde r=0$, so $n$0 (Braccini et al., 1 Oct 2025).

In the finite-manifold bosonic formulation, one starts from a finite set of normalized, linearly independent, pure Gaussian states $n$1 and defines

$n$2

Any pure state on this manifold has the form

$n$3

Standard cat states, such as $n$4, and squeezed-cat states are exactly $n$5 Gaussian-branch superpositions. In this sense, “two-branch Gaussian superpositions,” “Gaussian multiplets,” and finite Gaussian branch manifolds provide a bosonic realization of GCSs (Centrone et al., 16 Mar 2026).

A recurrent misconception is to identify GCSs only with superpositions of two coherent states. The finite-span formulation is broader: the branches may be displaced squeezed vacua, generic pure Gaussian states, or multimode tensor-product Gaussians; the branch number may be $n$6 or larger; and the state may be pure or mixed, provided its support remains inside the finite Gaussian manifold (Centrone et al., 16 Mar 2026).

2. Gaussian formalism for dynamics and measurements

The hybrid formalism is designed for Hamiltonians that are Gaussian in the bosonic variables and diagonal in the qubit basis. For one qubit, the operator-valued Gaussian Hamiltonian is

$n$7

with

$n$8

Conditioned on a qubit eigenvalue $n$9, the modes evolve under a branch Hamiltonian

$n$0

This is the precise sense in which the theory describes superpositions of Gaussian processes: each branch is a standard Gaussian evolution, while the full hybrid state is their coherent superposition (Braccini et al., 1 Oct 2025).

Imposing the GCS ansatz on the von Neumann equation yields closed ODEs for the complex branch covariances, first moments, and QRDM exponents. In the generic unitary case, the covariance matrices obey Riccati-type equations; in the important linear operator-valued-force case $n$1, the solution closes analytically in terms of a single mode propagator $n$2, and all branches share the same covariance evolution while contrasts $n$3 and phases $n$4 encode separation and dephasing between branches (Braccini et al., 1 Oct 2025).

The same closure persists under Markovian Gaussian noise. The most general open dynamics treated combines quadratic Hamiltonians, linear drift $n$5, Gaussian diffusion $n$6, antisymmetric drift $n$7, and qubit dephasing $n$8. The branch covariances then satisfy the open-system generalization of the Riccati/Lyapunov equations, while the QRDM exponents acquire the additional term $n$9. This gives an exact, truncation-free treatment of unitary and open hybrid dynamics within the GCS family (Braccini et al., 1 Oct 2025).

Measurements are equally natural in this language. Projective qubit measurements produce conditional bosonic states whose characteristic functions are linear combinations of branch Gaussians, hence generally non-Gaussian. Gaussian measurements on the modes are described by a reference Gaussian POVM with covariance $$\hat\varrho = \sum_{j,k\in\{\pm1\}} \hat\varrho_{jk}\otimes |j\rangle\langle k| ,$$0; the post-measurement QRDM element is

$$\hat\varrho = \sum_{j,k\in\{\pm1\}} \hat\varrho_{jk}\otimes |j\rangle\langle k| ,$$1

Noisy homodyne and heterodyne detection follow as special cases. The formalism therefore covers ideal and noisy measurements without leaving the branch representation (Braccini et al., 1 Oct 2025).

3. Canonical families and antecedents

Several earlier state families can be naturally interpreted as GCSs. One major example is the family of superpositions of coherent states determined by quadratic Gauss sums. Starting from the harmonic oscillator and the fractional evolution operator

$$\hat\varrho = \sum_{j,k\in\{\pm1\}} \hat\varrho_{jk}\otimes |j\rangle\langle k| ,$$2

one obtains

$$\hat\varrho = \sum_{j,k\in\{\pm1\}} \hat\varrho_{jk}\otimes |j\rangle\langle k| ,$$3

For rational angles $$\hat\varrho = \sum_{j,k\in\{\pm1\}} \hat\varrho_{jk}\otimes |j\rangle\langle k| ,$$4, this becomes a finite superposition of coherent states at the vertices of a regular $$\hat\varrho = \sum_{j,k\in\{\pm1\}} \hat\varrho_{jk}\otimes |j\rangle\langle k| ,$$5-gon,

$$\hat\varrho = \sum_{j,k\in\{\pm1\}} \hat\varrho_{jk}\otimes |j\rangle\langle k| ,$$6

with coefficients $$\hat\varrho = \sum_{j,k\in\{\pm1\}} \hat\varrho_{jk}\otimes |j\rangle\langle k| ,$$7 determined by quadratic Gauss sums. The $$\hat\varrho = \sum_{j,k\in\{\pm1\}} \hat\varrho_{jk}\otimes |j\rangle\langle k| ,$$8 case recovers the Yurke–Stoler coherent state, while $$\hat\varrho = \sum_{j,k\in\{\pm1\}} \hat\varrho_{jk}\otimes |j\rangle\langle k| ,$$9 yields multi-component “kittens.” The paper explicitly notes that there is no Gaussian envelope in $$\chi_{jk}(\tilde r)= \exp\!\left[ -\frac14 \tilde r^{\rm T}\sigma_{jk}\tilde r +i\,\tilde r^{\rm T}r_{jk} +r^{(0)}_{jk} \right].$$0-space; the structured object is the phase pattern across branches, not a Gaussian weight over branches (Spiridonov, 2021).

A second antecedent is the group-theoretic construction of “crystallized” Schrödinger cat states from cyclic and dihedral symmetries. There the seed is a generic Gaussian $$\chi_{jk}(\tilde r)= \exp\!\left[ -\frac14 \tilde r^{\rm T}\sigma_{jk}\tilde r +i\,\tilde r^{\rm T}r_{jk} +r^{(0)}_{jk} \right].$$1, and the branches are generated by phase-space rotations $$\chi_{jk}(\tilde r)= \exp\!\left[ -\frac14 \tilde r^{\rm T}\sigma_{jk}\tilde r +i\,\tilde r^{\rm T}r_{jk} +r^{(0)}_{jk} \right].$$2 or by dihedral operations combining rotations with conjugation/reflection. The resulting states are superpositions of Gaussians whose centers form regular polygons or dihedral completions in phase space, and whose Wigner functions and symplectic tomograms are explicit sums of Gaussian terms. The number and pattern of tomographic maxima and minima reflect the order of the underlying symmetry group (Hahn et al., 2022).

A third antecedent is the theory of generalized Gaussian cat states, where one superposes arbitrary pure Gaussian states,

$$\chi_{jk}(\tilde r)= \exp\!\left[ -\frac14 \tilde r^{\rm T}\sigma_{jk}\tilde r +i\,\tilde r^{\rm T}r_{jk} +r^{(0)}_{jk} \right].$$3

Its Wigner function has the standard decomposition into two Gaussian hills plus an interference term,

$$\chi_{jk}(\tilde r)= \exp\!\left[ -\frac14 \tilde r^{\rm T}\sigma_{jk}\tilde r +i\,\tilde r^{\rm T}r_{jk} +r^{(0)}_{jk} \right].$$4

but the interference phase is no longer generically linear. In arbitrary dimension it is governed by a quadratic form; in one degree of freedom the phase is hyperbolic. The same phase-space structure survives evolution in a thermal reservoir, and the paper also discusses mixed Gaussian superpositions generated by conditional Gaussian operations or by Kerr-type dynamics on thermal states (Nicacio et al., 2010).

These precursors differ in motivation—fractional revivals, finite groups, phase-space geometry—but they share the defining GCS feature: non-Gaussian states assembled from a finite number of Gaussian branches with structured coherence relations.

4. Finite-manifold encoding and quantitative characterization

The finite-branch bosonic formalism yields an exact finite-dimensional representation of any state supported on a Gaussian branch manifold. Given $$\chi_{jk}(\tilde r)= \exp\!\left[ -\frac14 \tilde r^{\rm T}\sigma_{jk}\tilde r +i\,\tilde r^{\rm T}r_{jk} +r^{(0)}_{jk} \right].$$5 and Gram matrix $$\chi_{jk}(\tilde r)= \exp\!\left[ -\frac14 \tilde r^{\rm T}\sigma_{jk}\tilde r +i\,\tilde r^{\rm T}r_{jk} +r^{(0)}_{jk} \right].$$6, one defines the Löwdin-orthonormalized basis

$$\chi_{jk}(\tilde r)= \exp\!\left[ -\frac14 \tilde r^{\rm T}\sigma_{jk}\tilde r +i\,\tilde r^{\rm T}r_{jk} +r^{(0)}_{jk} \right].$$7

Any mixed state on the manifold can then be represented by an effective $$\chi_{jk}(\tilde r)= \exp\!\left[ -\frac14 \tilde r^{\rm T}\sigma_{jk}\tilde r +i\,\tilde r^{\rm T}r_{jk} +r^{(0)}_{jk} \right].$$8 density matrix

$$\chi_{jk}(\tilde r)= \exp\!\left[ -\frac14 \tilde r^{\rm T}\sigma_{jk}\tilde r +i\,\tilde r^{\rm T}r_{jk} +r^{(0)}_{jk} \right].$$9

The map is an exact isospectral encoding onto a finite $N$0-dimensional Hilbert space, so all spectral quantities can be computed without Fock-space truncation (Centrone et al., 16 Mar 2026).

This immediately gives exact formulas for von Neumann and Rényi entropies from the eigenvalues of $N$1, and for relative-entropy non-Gaussianity

$N$2

where $N$3 is the Gaussian state with the same first and second moments as $N$4. For single-mode or multimode two-branch GCSs, the necessary first and second moments are assembled from branch data and Gaussian cross-moments, all obtainable analytically from Gaussian integrals (Centrone et al., 16 Mar 2026).

In the bipartite two-branch case, the effective manifold is $N$5, and after local Löwdin orthogonalization one obtains an exact two-qubit density matrix. For the pure Bell-like state

$N$6

with local overlaps $N$7 and $N$8, the entanglement negativity has the closed form

$N$9

The same framework also gives an exact “which-branch dephasing” model, both for single-mode cat decoherence and for Bell-like two-branch GCSs, while keeping the support inside the same finite branch span (Centrone et al., 16 Mar 2026).

A complementary computational development appears in the exact simulation of realistic GKP cluster states, where resource states are represented as sums of Gaussian distributions in phase space and propagated analytically through Gaussian circuits and homodyne measurements. That formalism exactly simulates cat-bred GKP states and multimode GKP clusters more efficiently than standard Fock-basis truncations, while also clarifying where heuristic Gaussian random noise models fail at the level of conditional stabilizer expectation values (Banic et al., 14 Apr 2025).

5. Generation and control protocols

A general route to GCSs is conditional measurement on multimode Gaussian resources. For an arbitrary $J,K\in\mathcal M$0-mode Gaussian input, photon-number-resolving detection on $J,K\in\mathcal M$1 modes leaves the heralded mode in the form

$J,K\in\mathcal M$2

with $J,K\in\mathcal M$3, the total detected photon number. Thus the output factorizes into a Gaussian gate acting on a finite Fock superposition. Cat states are treated explicitly as targets of this architecture: for small amplitudes, an even cat is well approximated by $J,K\in\mathcal M$4; for larger amplitudes, by $J,K\in\mathcal M$5. The same formalism extends to mixed Gaussian inputs and experimental imperfections such as photon loss (Su et al., 2019).

In fully optical state engineering, deterministic protocols based on Gaussian operations and photon-number measurements generate large-amplitude squeezed cat states of the form

$J,K\in\mathcal M$6

The preparation uses only Gaussian unitaries—beam splitters, squeezers, displacements, phase rotations—and photon-number-resolving detection as the only explicit non-Gaussian ingredient. Scheme I deterministically maps a Fock state $J,K\in\mathcal M$7 to a squeezed cat with

$J,K\in\mathcal M$8

and the resulting cats can be bred into approximate GKP states with

$J,K\in\mathcal M$9

Scheme II deterministically approaches a two-component cat from arbitrary even or odd parity input states by repeated weak tapping and PNR on squeezed-vacuum ancillas (Winnel et al., 2023).

Remote preparation furnishes a distributed version of the same paradigm. Starting from a two-mode squeezed vacuum sent through lossy Gaussian channels, single- or multi-photon subtraction at one node and homodyne conditioning remotely prepare optical cat states at the other node. The measurement angle rotates the remotely prepared cat in phase space, and the protocol shows marked asymmetry under loss: the remotely prepared cat state is much more tolerant of loss in Alice’s channel than in Bob’s channel. In the ideal model with 3 dB squeezing and no loss, the fidelity can reach $$\chi_{JK}(\tilde r)= \exp\!\left[ -\frac14 \tilde r^{\rm T}\sigma_{JK}\tilde r +i\,\tilde r^{\rm T}r_{JK} -\mathcal C_{JK}+i\phi_{JK} \right].$$0; experimentally, for $$\chi_{JK}(\tilde r)= \exp\!\left[ -\frac14 \tilde r^{\rm T}\sigma_{JK}\tilde r +i\,\tilde r^{\rm T}r_{JK} -\mathcal C_{JK}+i\phi_{JK} \right].$$1, the remotely prepared state had $$\chi_{JK}(\tilde r)= \exp\!\left[ -\frac14 \tilde r^{\rm T}\sigma_{JK}\tilde r +i\,\tilde r^{\rm T}r_{JK} -\mathcal C_{JK}+i\phi_{JK} \right].$$2, $$\chi_{JK}(\tilde r)= \exp\!\left[ -\frac14 \tilde r^{\rm T}\sigma_{JK}\tilde r +i\,\tilde r^{\rm T}r_{JK} -\mathcal C_{JK}+i\phi_{JK} \right].$$3, and central negativity $$\chi_{JK}(\tilde r)= \exp\!\left[ -\frac14 \tilde r^{\rm T}\sigma_{JK}\tilde r +i\,\tilde r^{\rm T}r_{JK} -\mathcal C_{JK}+i\phi_{JK} \right].$$4 (Han et al., 2023).

Within continuous-variable cluster states, the Photon-counting-Assisted Node-Teleportation Method (PhANTM) embeds and preserves cat states directly inside an otherwise Gaussian graph. The realistic PhANTM Kraus operator has the form

$$\chi_{JK}(\tilde r)= \exp\!\left[ -\frac14 \tilde r^{\rm T}\sigma_{JK}\tilde r +i\,\tilde r^{\rm T}r_{JK} -\mathcal C_{JK}+i\phi_{JK} \right].$$5

and in the weak-reflectivity/high-squeezing limit reduces to a Hermite-polynomial gate,

$$\chi_{JK}(\tilde r)= \exp\!\left[ -\frac14 \tilde r^{\rm T}\sigma_{JK}\tilde r +i\,\tilde r^{\rm T}r_{JK} -\mathcal C_{JK}+i\phi_{JK} \right].$$6

Repeated PhANTM steps generate cat states from Gaussian nodes, stabilize cat amplitudes against Gaussian noise during cluster teleportation, and embed standard cat-breeding and GKP-breeding protocols inside cluster-state processing (Eaton et al., 2021).

6. Phase-space signatures, applications, and computational uses

The characteristic phase-space signature of a GCS is the coexistence of Gaussian branch peaks and interference terms. In generalized two-branch Gaussian cats, the interference term is a Gaussian envelope multiplied by an oscillatory factor with a quadratic phase, and in one degree of freedom this phase is hyperbolic rather than linear. In cyclic and dihedral Gaussian multiplets, Wigner functions exhibit polygonal or dihedral symmetries and tomograms whose fringe counts reflect the order of the underlying group. In Gauss-sum coherent superpositions, rational fractional-evolution times produce regularly spaced branches whose relative phases are controlled by number-theoretic data (Nicacio et al., 2010, Hahn et al., 2022, Spiridonov, 2021).

These phase-space structures support several application domains. In bosonic encoding theory, two-branch GCSs realize cat qubits, while larger finite Gaussian multiplets provide cat qudits and related encodings; exact finite-manifold methods then supply entropies, non-Gaussianity measures, and entanglement negativities directly from Gaussian overlaps (Centrone et al., 16 Mar 2026). In optical and cavity platforms, deterministic squeezed-cat generation and breeding furnish primitive GKP resources and support error-correction studies that surpass break-even under pure loss in the simulated regimes (Winnel et al., 2023). Exact Gaussian-sum simulation methods extend this perspective to multimode GKP Bell pairs and cluster states, including stabilizer and entanglement-witness calculations beyond the range of practical Fock truncations (Banic et al., 14 Apr 2025).

The hybrid GCS formalism also captures qubit-mediated Gaussian-process superpositions that are not naturally expressible as finite coherent-state sums. Two paradigmatic examples are treated explicitly. First, a squeezed, leaking, and homodyned resonator dispersively coupled to two qubits yields measurement-based entanglement between the qubits through branch-dependent Gaussian resonator trajectories. Second, a levitated nanoparticle undergoing Stern–Gerlach interferometry in a diffusive environment develops conditional motional Wigner negativity after spin measurement, while the contrast function cleanly separates branch distinguishability, diffusion, and qubit dephasing (Braccini et al., 1 Oct 2025).

7. Limits, equivalences, and interpretive issues

Deterministic Gaussian conversion protocols clarify which non-Gaussian resources lie in the same Gaussian orbit as cats and which do not. A central result is that cat and binomial code states are approximately equivalent already at finite energy: for many parameter pairs $$\chi_{JK}(\tilde r)= \exp\!\left[ -\frac14 \tilde r^{\rm T}\sigma_{JK}\tilde r +i\,\tilde r^{\rm T}r_{JK} -\mathcal C_{JK}+i\phi_{JK} \right].$$7, there exists a cat codeword with fidelity $$\chi_{JK}(\tilde r)= \exp\!\left[ -\frac14 \tilde r^{\rm T}\sigma_{JK}\tilde r +i\,\tilde r^{\rm T}r_{JK} -\mathcal C_{JK}+i\phi_{JK} \right].$$8 to the corresponding binomial codeword, and the identity map is already optimal in the high-fidelity region. Photon-added and photon-subtracted squeezed states can also be converted to cats with substantially improved performance by adding an optimized single-mode squeezing operation. By contrast, several seemingly favorable conversions fail qualitatively despite moderate or even high fidelities: fidelity alone can be misleading when the target resource depends on specific Wigner-negativity patterns or rotational symmetries (Hahn et al., 2022).

A second interpretive issue concerns the adjective “Gaussian” in GCS. In the finite-branch literature it refers to the branch states themselves—coherent states, squeezed states, or generic pure Gaussian states—not to a Gaussian weight over branch labels. This is explicit in the Gauss-sum construction, where the envelope in phase space is entirely determined by the coherent-state Gaussians themselves, while the nontrivial structure lies in the phase pattern across branches (Spiridonov, 2021).

A third limit is exactness. The finite-manifold encoding is exact only for states whose support is strictly inside a finite branch span. Generic Gaussian noise channels such as loss or thermal noise typically move a state outside that manifold, requiring either an enlarged branch set or approximation. Similarly, the exact hybrid GCS formalism assumes Gaussian continuous-variable dynamics and couplings diagonal in the discrete-variable basis; it is exact within that class, but not a universal representation of arbitrary hybrid non-Gaussian dynamics (Centrone et al., 16 Mar 2026, Braccini et al., 1 Oct 2025).

Taken together, these developments position GCSs as a unifying technical category for non-Gaussian states assembled from Gaussian components: finite coherent-state polygons, group-generated Gaussian multiplets, two-branch and multimode bosonic encodings, cluster-embedded cats and grids, and qubit-conditioned superpositions of Gaussian processes. The concept is broad enough to connect cat-state physics, bosonic coding, number-theoretic coherent superpositions, hybrid interferometry, and exact simulation techniques, while remaining concrete enough to admit closed-form dynamics, exact finite-dimensional reductions, and experimentally grounded preparation protocols.