---
title: Multimode Cat States in Bosonic Quantum Systems
url: https://www.emergentmind.com/topics/multimode-cat-states
type: topic
---

# Multimode Cat States in Bosonic Quantum Systems

Searching arXiv for recent and foundational papers on multimode cat states and closely related bosonic cat resources.
Multimode cat states are bosonic superposition states in which coherent-state branches are correlated across several modes, typically as coherent superpositions of branch configurations such as \(\ket{\alpha}^{\otimes M}\pm\ket{-\alpha}^{\otimes M}\), more general \(n\)-mode states \(\ket{\mathrm{cat}_{\alpha_1\cdots\alpha_n}}\), or graph- and code-structured variants defined inside a coherent-branch manifold [2607.02997, 2201.03553]. In the large-\(|\alpha|\) limit, the branches \(\ket{\pm\alpha}\) become approximately orthogonal and define an effective logical-qubit basis, which is why multimode cat states are recurrent in GHZ-like entanglement, cluster-state constructions, W-type resources, and rotationally symmetric bosonic codes [2607.02997, 1905.04650, 2412.15134].

## 1. Canonical definitions and state families

A standard \(n\)-mode entangled coherent state is the equal superposition  
\[
\left|\psi\right\rangle_{n\text{-cat}}
=
{\cal N}_{n}\left(
\left|\alpha\right\rangle_{1}\left|\alpha\right\rangle_{2}\cdots \left|\alpha\right\rangle_{n}
+
\left|-\alpha\right\rangle_{1}\left|-\alpha\right\rangle_{2}\cdots \left|-\alpha\right\rangle_{n}
\right),
\]
with
\[
{\cal N}_{n}=\left[2+2e^{-2n|\alpha|^{2}}\right]^{-1/2}.
\]
For \(|\alpha|^2\gg1\), this state approaches a continuous-variable analog of the \(n\)-qubit GHZ state [1905.04650]. A more general \(n\)-mode Schrödinger cat state is
\[
|\mathrm{cat}_{\alpha_1\alpha_2\cdots\alpha_n}\rangle
=
\frac{1}{\sqrt{2+2\prod_{j=1}^n q_{\alpha_j}}}
\left(
|\alpha_1,\alpha_2,\ldots,\alpha_n\rangle
+
|-\alpha_1,-\alpha_2,\ldots,-\alpha_n\rangle
\right),
\]
where \(q_{\alpha_j}=e^{-2|\alpha_j|^2}\) [2201.03553]. This form makes explicit that multimode cat states need not be symmetric across modes.

A second major family is the rotationally symmetric cat-code basis \(\lvert \mu\rangle_{\alpha,D}\), written either as a superposition of \(2D\) coherent states around a circle in phase space or as a Fock-state comb. Here \(D\) is the code distance and half the number of coherent-state components, while \(\mu\in\{0,1\}\) labels the logical qubit. The code can ideally correct up to \(D-1\) photon losses [2412.15134]. This family broadens “multimode cat state” beyond the two-branch GHZ form: the essential feature is coherent support on a structured set of separated quasi-classical branches.

Compact multimode cats also arise in \(\mathfrak{su}(2)\) and \(\mathfrak{su}(d)\) settings. For two modes, the binomial cat states are
\[
|\psi_\pm\rangle
=
\frac{1}{\sqrt{2}}
\left(|N,\xi\rangle\pm|N,-\xi\rangle\right),
\]
where the coherent branches lie in a fixed total-excitation manifold \(N\) and the amplitudes follow binomial weights in the Schwinger representation [2607.03302]. In this sense, multimode cat states include both non-compact coherent-state superpositions and compact finite-support constructions.

## 2. Branch manifolds, symmetry sectors, and effective logical structure

A unifying operator-level description begins by confining each oscillator mode to the two-branch coherent support
\[
\mathcal B_\alpha=\mathrm{span}\{\ket{\alpha},\ket{-\alpha}\}.
\]
The multimode branch Hamiltonian
\[
H_{\mathrm{br}}=\sum_{j=1}^{M}(a_j^2-\alpha^2)^\dagger(a_j^2-\alpha^2)
\]
has every branch configuration
\[
\ket{\bm s}=\ket{s_1\alpha,s_2\alpha,\dots,s_M\alpha},\qquad s_j=\pm1
\]
in its zero-energy space, so it fixes the coherent-state alphabet but leaves a \(2^M\)-dimensional degeneracy [2607.02997]. State-dependent positive-semidefinite constraints then select the desired correlations and symmetry sector. For GHZ cats, connected alignment terms \((a_i-a_j)^\dagger(a_i-a_j)\) reduce the branch manifold to \(\mathrm{span}\{\ket{\alpha}^{\otimes M},\ket{-\alpha}^{\otimes M}\}\), and a global parity selector chooses \(\ket{\mathrm{GHZ}_\pm}\). For cluster cats, pair alignment first reduces the branch space and stabilizer-like constraints isolate the unique cluster cat. For W cats, a fixed-defect sector is selected before a noncommuting exchange term produces the equal-amplitude symmetric W superposition [2607.02997].

In the large-\(|\alpha|\) limit, \(\langle \alpha|-\alpha\rangle=e^{-2|\alpha|^2}\approx0\), so
\[
\ket{0_L}=\ket{-\alpha},\qquad \ket{1_L}=\ket{\alpha}
\]
defines an effective qubit basis. The operator dictionary given in the parent-Hamiltonian construction maps parity to \(X\), alignment penalties to Ising-like \(I-Z_iZ_j\), exchange to \(X_iX_j+Y_iY_j\), and global parity to \(\prod_j X_j\) [2607.02997]. This is the direct bridge between coherent-state bosonic engineering and stabilizer-based quantum information processing.

Driven Kerr parametric oscillators provide an allied viewpoint. For a single mode,
\[
H=-K a^{\dagger 2}a^2+\left(\epsilon a^{\dagger 2}+\epsilon^* a^2\right)
\]
preserves photon-number parity, and \(\lvert0\rangle\) and \(\lvert1\rangle\) are adiabatically connected to even and odd cat states. For two coupled modes, the multimode cat states
\[
C^\pm(\vec{\alpha},\vec{\sigma})
=
\frac{1}{N^\pm(\vec{\alpha},\vec{\sigma})}
\left(
\bigotimes_{j=1}^N \lvert \sigma_j \alpha_j\rangle
\pm
\bigotimes_{j=1}^N \lvert -\sigma_j \alpha_j\rangle
\right)
\]
arise as degenerate ground-manifold states, and in the large-amplitude regime the associated proto-Bell cats approach the Bell-cat basis \(\Phi^\pm,\Psi^\pm\) [2505.02497]. The structural point is the same in both formulations: multimode cat states are selected from a branch manifold by symmetry, parity, or coupling constraints.

## 3. Preparation mechanisms

A broad analytic route starts from an arbitrary \(N\)-mode Gaussian state, measures \((N-M)\) modes in the Fock basis, and conditionally prepares an \(M\)-mode non-Gaussian output. For the single-mode heralding case relevant to cat generation, the output always factorizes as
\[
|{\psi_1}\rangle
=
\hat D(\beta)\hat S(\zeta)\sum_{n=0}^{n_{\max}} c_n|n\rangle,
\]
so the heralded state is a Gaussian gate acting on a finite Fock superposition [1902.02331]. The even cat is a canonical target because it contains only even Fock states. For small \(\alpha\), it is well approximated by \(c_0|0\rangle+c_2|2\rangle\); for larger \(\alpha\), the better ansatz is \(\hat S(\zeta_1)(c_0|0\rangle+c_2|2\rangle)\). In the optimized two-mode realization, Table I reports fidelities above \(97\%\) and success probabilities above \(10\%\) for \(\alpha\le2\), including
\[
\mathcal F_{\max}=1.0000,\qquad P_{\max}=18.12\%\quad\text{for }\alpha=0.25,
\]
with input squeezing \(\zeta_{01}\in(1.1587,1.6150)\), corresponding to about \(10\)–\(14\) dB [1902.02331].

Remote preparation based on a distributed two-mode squeezed state replaces local non-Gaussian synthesis by conditional steering. Alice subtracts photons and performs a homodyne projective measurement on her mode of
\[
|\psi_0\rangle_{AB}
=
\frac{1}{\cosh r}\sum_{m=0}^{\infty}\tanh^m r\,|m,m\rangle_{AB},
\]
so that Bob’s conditional state becomes an odd cat-like superposition. The homodyne angle \(\theta\) remotely rotates the cat, and the loss asymmetry is pronounced: the prepared cat tolerates much more loss in Alice’s channel than in Bob’s channel. The same framework predicts amplitudes larger than \(2\) by increasing the squeezing level and subtracting photon numbers; for three-photon subtraction and squeezing around \(8.4\) dB, the paper reports amplitude \(\sim2.61\) and fidelity \(\sim96\%\) [2304.08863].

Deterministic microwave generation proceeds differently. A transmon qubit dispersively coupled to a 3D cavity reflects a train of coherent microwave pulses so that
\[
\ket{0}\rightarrow\ket{\alpha},\qquad \ket{1}\rightarrow\ket{-\alpha},
\]
and after \(n\) reflections the joint state is
\[
\frac{1}{\sqrt2}\Big(\ket0\ket{\alpha}^{\otimes n}+\ket1\ket{-\alpha}^{\otimes n}\Big).
\]
Qubit rotation and projective measurement then prepare even or odd flying multipartite cat states \(\mathscr N(\ket{\alpha}^{\otimes n}\pm\ket{-\alpha}^{\otimes n})\). Full tomography was demonstrated up to four photonic modes, with reconstructed fidelities \(0.66\) for the bipartite even cat, \(0.62\) for the bipartite odd cat, \(0.54\) for the tripartite even cat, and \(0.53\) for the quadripartite even cat [2112.01199].

Resource-efficient optical scaling can also mean preparing several cats at once. A nondegenerate optical parametric amplifier operated in the parametric de-amplification regime generates two squeezed vacuum states with orthogonal squeezing directions in phase space; subtracting one photon from each yields two odd cat states with orthogonal superposition directions, amplitudes about \(1.19\) and \(1.21\), fidelities \(61\%\) and \(60\%\), Wigner-function negativities at the origin about \(-0.12\) and \(-0.14\), and a generation rate about \(20\) Hz [2305.08426]. This is a distinct scaling notion from a single entangled multimode cat, but it is directly relevant to multi-cat-state architectures.

## 4. Tomography, phase-space structure, and multipartite nonclassical correlations

The Wigner representation remains central, but multimode cat states admit a complete positive-probability description through center-of-mass tomography. For \(N\) degrees of freedom, the center-of-mass tomogram
\[
w_{cm}(X,\vec\mu,\vec\nu)
=
\mathrm{Tr}\!\left[\hat\rho\,\delta\!\left(X-\vec\mu\,\hat{\vec q}-\vec\nu\,\hat{\vec p}\right)\right]
\]
is nonnegative, normalized, and related to the Wigner function by the Radon-type projection
\[
w_{cm}(X,\vec\mu,\vec\nu)
=
\int W(\vec q,\vec p)\,
\delta(X-\vec\mu\cdot\vec q-\vec\nu\cdot\vec p)\,
d\vec q\,d\vec p .
\]
For a two-mode cat state, the explicit tomogram contains two displaced Gaussian contributions and interference cross terms, and the same formalism yields the linear entropy
\[
S_\pm(\alpha_1,\alpha_2)
=
\frac12-\frac12\left(1\pm e^{-2|\alpha_1|^2-2|\alpha_2|^2}\right)^{-2}
\left(e^{-2|\alpha_1|^2}\pm e^{-2|\alpha_2|^2}\right)^2
\]
as an entanglement diagnostic [1803.09064].

Experimental itinerant-state tomography of flying multipartite cats is based on complex-amplitude histograms \(S_m=I_m+iQ_m\) of the outgoing field and maximum-likelihood reconstruction of the full \(2n\)-dimensional histogram \(D(I_1,Q_1,\dots,I_n,Q_n)\). In the reported circuit-QED experiment, the reconstructed density matrices show dominant populations in the coherent-state subspace, with cat-basis subspace traces about \(94\%\) and \(92.5\%\) for single-mode even and odd cats and about \(80.6\%\) and \(73.6\%\) for quadripartite even and odd cats. Multipartite entanglement was certified by localizable entanglement, and hybrid qubit-cat entanglement by negativities
\[
N_1=0.192\pm0.014,\qquad N_2=0.135\pm0.010
\]
for qubit–single-cat and qubit–bipartite-cat systems [2112.01199].

For Bell-type nonlocality, rotated quantum-number parity measurements are especially effective. Using
\[
\sigma_k(\phi_k)=R_{k,z}(\phi_k)\,P_k\,R_{k,z}^\dagger(\phi_k),
\qquad
P_k=(-1)^{a_k^\dagger a_k},
\]
the Mermin-Klyshko signal can approach the maximal quantum value already when the average quantum number per mode \(\bar n=|\alpha|^2\) is only about \(1\). For the three-mode cat, \(S_3\) reaches \(3.916\) at \(\alpha=1\), close to the upper bound \(4\); for four and five modes, the rescaled signals at \(\alpha=1\) are \(R_4\simeq2.758\) and \(R_5\simeq3.882\), and the inequality violation grows exponentially with the number of entangled modes [1905.04650].

## 5. Decoherence, loss, numerical simulation, and stabilization

A multimode cat can be reduced to an effective two-mode interfering-state model by partitioning modes into a “system” block and an “environment” block. In that reduction, the Schmidt number
\[
K=\frac{1}{\sum_k\lambda_k^2}
\]
and the visibility
\[
V=\frac{I_{\max}-I_{\min}}{I_{\max}+I_{\min}}
\]
are linked by
\[
V=\sqrt{2-\frac{2}{K}}.
\]
For the identical-mode case \(\alpha_1=\cdots=\alpha_n=\alpha\), the paper gives the explicit multimode-cat visibility
\[
V=\exp(-2m\alpha^2),
\]
where \(m\) is the number of environment modes measured or lost [2201.03553]. The same analysis introduces the “health” parameter
\[
H(m)=\ln\frac{p_m}{1-p_m}
=
-4\alpha^2\sum_{j=1}^m y_j,
\]
which quantifies how sequential conditioning breaks the balance between the two macroscopic branches.

The asymptotic size of a cat does not, by itself, force the nonclassical distance toward its maximum. For single-mode even and odd Schrödinger cats,
\[
\delta(\rho)=\inf_{\sigma\in S_{cl}} D(\rho,\sigma)
\]
obeys bounds determined by the Husimi function, and both even and odd cat states satisfy \(\delta(\ket{\psi_\pm})\to\frac12\) as \(\beta\to\infty\). The paper explicitly emphasizes that “the nonclassical distance of the cat states is bounded away from unity regardless of the superposition amplitude,” and that nonclassical distance is not necessarily monotonically increasing with respect to macroscopicity [1701.07688]. This is a recurrent corrective to the intuition that wider phase-space separation always implies “more nonclassicality.”

Large-array dynamics requires scalable numerics. In a driven-dissipative chain with two-photon drive, one-photon loss, two-photon loss, and engineered nonlocal dissipation, the positive-\(P\) representation yields exact stochastic differential equations for multimode cat-state transients and supports simulations up to \(N=21\) sites, with the paper noting scalability to \(N>100\). The main limitation is parity: the estimator
\[
\Pi
=
\left\langle
\exp\!\left(i\pi\sum_j \hat a_j^\dagger \hat a_j\right)
\right\rangle
=
\left\langle
\exp\!\left(-2\sum_j \alpha_j\beta_j\right)
\right\rangle_s
\]
is exponentially sensitive to rare trajectories and becomes numerically unstable even when photon-number and correlation observables remain accurate [2601.07049].

Autonomous protection can also be designed directly. Dissipative preparation and stabilization of compact multimode cats is achieved by multiple Lindblad jump operators whose common kernel is the target state manifold. For the two-mode binomial cat, the proposed operators
\[
\hat{L}_{1,+}^{(2)}=\sqrt{\kappa_1}\,\hat a^\dagger(\hat a^\dagger\hat a+\hat b^\dagger\hat b-N),
\qquad
\hat{L}_2^{(2)}=\sqrt{\kappa_2}\,(\hat a^2\pm\hat b^2)
\]
pump the system from vacuum into the fixed-\(N\) manifold and then lock the cat structure inside it. The reported effect is to prepare the target from vacuum with high fidelity and to extend its lifetime by several orders of magnitude compared to natural decay times [2607.03302].

## 6. Quantum-information roles, protocol variants, and conceptual boundaries

In bosonic error correction, multimode cat resources appear both as logical states and as ancillae. Teleportation-based correction for rotationally symmetric codes uses controlled-rotation gates
\[
\hat R_{MN}=e^{i(4\pi/MN)\hat n_1\hat n_2}
\]
and phase measurements, but the ancillary state need not be an exact logical \(\lvert+\rangle\) cat. Many-component Yurke–Stoler states, although not themselves in the cat-code space, can still function as “bridge” states because known extra rotations \(e^{i\theta\hat n}\) and known photon-loss events \(\hat a^m\) can be absorbed into the classical decoding rules [2412.15134]. This broadens the operational notion of a useful multimode cat resource.

Distributed-network applications use mode multiplexing rather than only local multimode entanglement. In the quantum-repeater proposal based on entangled phase-modulated multimode coherent states, electro-optic modulation transforms a one-mode cat into a frequency-comb state
\[
|\pm\alpha\rangle_A\to|\pm\boldsymbol\alpha_A\rangle=\bigotimes_{\mu=-S}^{S}|\pm\alpha_\mu\rangle_A,
\]
and interference of selected subcarrier modes at a central beam splitter, followed by photon counting, heralds entanglement between remote memories [2211.03597]. A related three-mode teleportation protocol builds an entangled mode pair by beam-splitting a cat state and then teleports an unknown cat-encoded qubit using beam splitters, parity or joint-parity measurements, and a final binary-outcome dispersive readout; with multiple measurements, the protocol can approach unit fidelity [2412.15011].

Two conceptual boundaries are important. First, not every “cat-like” structured state is a canonical even or odd cat. In the doubly pumped intraresonance Kerr microresonator, the reduced \(n=2\) resonant photon-conversion model exhibits a residual \(\mathbb Z_3\) symmetry and produces multicomponent Schrödinger-cat-like states characterized by Wigner negativity, non-Poissonian statistics, pump-mode quadrature squeezing, and large single-mode Schmidt numbers, but the paper explicitly states that these are not the canonical even or odd coherent states of Dodonov, Malkin, and Man’ko [2606.19085]. Second, classical multimode analogs exist but are not quantum cat states. The OAM-based “analogous cat state” is a coherent superposition over many orbital-angular-momentum modes with high similarity to theory in state tomography, yet the paper explicitly notes that it has zero quantum macroscopicity because it is made of classical electromagnetic modes [1807.05498].

Taken together, these constructions show that “multimode cat state” is best understood as a family of structured bosonic superpositions defined by coherent-branch support, intermode correlation constraints, and symmetry-sector selection. In current usage, the term covers GHZ-like entangled coherent states, compact multinomial and binomial cats, distributed frequency-multiplexed resources, Bell-cat states in coupled Kerr systems, and fault-tolerant ancillary states that retain the relevant syndrome-extraction structure even when they are not exact codewords [2607.02997, 2505.02497, 2412.15134].

Source: https://www.emergentmind.com/topics/multimode-cat-states