Multimode Schrödinger Cat States

• Multimode Schrödinger Cat States are coherent superpositions of macroscopically distinct bosonic modes, exhibiting quantum interference and entanglement.
• They are realized via methods like Kerr parametric oscillator networks and dispersive reflection protocols to generate high-fidelity, entangled multimode configurations.
• These states underpin advanced quantum information processing, quantum metrology, and error correction, while facing challenges such as photon loss and decoherence.

A multimode Schrödinger cat state generalizes the concept of a macroscopic quantum superposition to a system involving several bosonic modes, each of which may exhibit distinct or joint quantum coherence phenomena. These states manifest as coherent superpositions of macroscopically distinguishable multimode configurations, typically characterized by simultaneous superpositions of coherent states of opposite phases in each mode. Multimode cat states represent a highly significant class of entangled non-Gaussian continuous-variable states, central to quantum information processing, quantum error correction, quantum metrology, and the study of quantum-to-classical transitions.

1. Formal Structure and Taxonomy

The canonical multimode Schrödinger-cat state in $N$ modes, each with coherent amplitude $\alpha_j$, takes the form

$$|\mathrm{Cat}_N^\pm(\alpha_1,\ldots,\alpha_N)\rangle = \mathcal{N}_N^\pm \big( |\alpha_1,\ldots,\alpha_N\rangle \pm |-\alpha_1,\ldots,-\alpha_N\rangle \big)$$

where $$|\alpha_1,\ldots,\alpha_N\rangle = \bigotimes_{j=1}^N|\alpha_j\rangle$$, and the normalization factor is

$$\mathcal{N}_N^\pm = \left[2 \pm 2 e^{-2 \sum_j |\alpha_j|^2}\right]^{-1/2}.$$

Variants include sign-vector cat states with arbitrary patterns of positive and negative amplitudes, as well as generalized cat states in discrete (e.g., $\mathbb{Z}_n$ symmetric) Rabi models exhibiting superpositions of multiple phase-rotated multimode coherent products (Wang et al., 2021, Lotkov et al., 10 Sep 2025).

A second major class comprises entangled resource states such as two-mode (or $N$-mode) squeezed vacua, which qualify as Schrödinger-cat-like by virtue of their exponentially macroscopic quantum fluctuations and noise tolerance—scaling identically with standard single-mode photonic cats according to quantum Fisher information and coarse-grained noise measures (Oudot et al., 2014).

The hybridization with discrete-variable systems or higher-dimensional ancillae (qubits, qutrits) introduces "hybrid" cat states, in which macroscopic CV superpositions are coherently entangled with individual or multipartite DV subsystems (Wang et al., 2021, Hoshi et al., 2024, Lotkov et al., 10 Sep 2025).

2. Physical Realizations and State Preparation

Kerr Parametric Oscillator Networks: In circuit QED and superconducting platforms, driven Kerr nonlinear oscillators ("KPOs") provide a natural setting in which single- and multimode cat states arise as macroscopic ground-state manifolds in the presence of a two-photon (parametric) drive. For $N$-mode systems,

$$H(t) = \sum_{j=1}^N \big[-K_j a_j^{\dagger 2} a_j^2 + \epsilon_j(t)(a_j^{\dagger 2} + a_j^2)\big] - \sum_{i<j} K_{ij} a_i^\dagger a_i a_j^\dagger a_j + \sum_{i<j} [\epsilon_{ij}(t) a_i^\dagger a_j^\dagger + \epsilon_{ij}^*(t) a_i a_j],$$

with appropriate choice of drive amplitudes and cross-Kerr couplings, adiabatic ramping creates joint even/odd cat manifolds; diabatic switching of couplings enables fast initialization of entangled Bell cats and their $N$-mode generalizations with fidelities exceeding $\alpha_j$0 for moderate ramp times and amplitudes (Resch et al., 5 May 2025, Hoshi et al., 2024).

Dispersive Reflection Protocols: A deterministic, highly scalable method relies on sequentially reflecting coherent microwave pulses from a cavity dispersively coupled to a superconducting qubit. Through a pulse-rotation protocol, temporal modes are linked coherently, yielding multipartite GHZ-type photonic cat states. Conditioning on qubit measurement projects the photonic subsystem into even or odd $\alpha_j$1-mode cats, demonstrated experimentally up to quadripartite cats (Wang et al., 2021).

Dissipative Engineering and Phase-Space Methods: In open-system settings such as arrays of coupled resonators with single- and two-photon losses, the transient dynamics leading to multimode cat-state formation can be simulated accurately (up to $\alpha_j$2 modes) using the positive-P representation. Stochastic differential equations derived from the Lindblad equation enable linear-in-$\alpha_j$3 resource scaling, although global parity and Wigner negativities remain challenging due to sampling noise and boundary-term errors (Shi et al., 11 Jan 2026).

Hybrid Discrete-Continuous Encodings: Platforms that realize Kerr cats in coupled KPOs can exploit protocols for entanglement-preserving conversion between DV-encoded Bell states and cat states, as well as implement entangling gates (e.g., $\alpha_j$4) between cat qubits, completing a universal set for networked cat-qubit architectures (Hoshi et al., 2024).

$\alpha_j$5-Symmetric Rabi Models: In both one- and two-mode $\alpha_j$6 Rabi models, deep-strong coupling generates ground and low-lying excited states that are superpositions of multiple macroscopically distinct phase-rotated coherent products, with closed-form joint Wigner functions for composite (qutrit-boson) hybrid cat states (Lotkov et al., 10 Sep 2025).

3. Macroscopicity, Nonclassicality, and Entanglement Quantification

The identification of a multimode state as genuinely Schrödinger-cat-like requires rigorous quantification of macroscopicity, nonclassicality, and entanglement:

• Coarse-Grained Distinguishability: The maximal tolerable detector noise $\alpha_j$7 for distinguishing branches of a superposition scales as $\alpha_j$8 for both standard photonic cats and multimode squeezed vacua, situating both within the same macroscopic class (Oudot et al., 2014).
• Quantum Fisher Information (QFI) Effective Size: For $\alpha_j$9-mode cat-like states, the effective size is defined as

$$|\mathrm{Cat}_N^\pm(\alpha_1,\ldots,\alpha_N)\rangle = \mathcal{N}_N^\pm \big( |\alpha_1,\ldots,\alpha_N\rangle \pm |-\alpha_1,\ldots,-\alpha_N\rangle \big)$$0

where $$|\mathrm{Cat}_N^\pm(\alpha_1,\ldots,\alpha_N)\rangle = \mathcal{N}_N^\pm \big( |\alpha_1,\ldots,\alpha_N\rangle \pm |-\alpha_1,\ldots,-\alpha_N\rangle \big)$$1 is the QFI and $$|\mathrm{Cat}_N^\pm(\alpha_1,\ldots,\alpha_N)\rangle = \mathcal{N}_N^\pm \big( |\alpha_1,\ldots,\alpha_N\rangle \pm |-\alpha_1,\ldots,-\alpha_N\rangle \big)$$2 are local quadratures. For two-mode squeezed vacua, $$|\mathrm{Cat}_N^\pm(\alpha_1,\ldots,\alpha_N)\rangle = \mathcal{N}_N^\pm \big( |\alpha_1,\ldots,\alpha_N\rangle \pm |-\alpha_1,\ldots,-\alpha_N\rangle \big)$$3 (with squeezing $$|\mathrm{Cat}_N^\pm(\alpha_1,\ldots,\alpha_N)\rangle = \mathcal{N}_N^\pm \big( |\alpha_1,\ldots,\alpha_N\rangle \pm |-\alpha_1,\ldots,-\alpha_N\rangle \big)$$4, mean photon number $$|\mathrm{Cat}_N^\pm(\alpha_1,\ldots,\alpha_N)\rangle = \mathcal{N}_N^\pm \big( |\alpha_1,\ldots,\alpha_N\rangle \pm |-\alpha_1,\ldots,-\alpha_N\rangle \big)$$5), and practical certifiable values exceed $$|\mathrm{Cat}_N^\pm(\alpha_1,\ldots,\alpha_N)\rangle = \mathcal{N}_N^\pm \big( |\alpha_1,\ldots,\alpha_N\rangle \pm |-\alpha_1,\ldots,-\alpha_N\rangle \big)$$6–$$|\mathrm{Cat}_N^\pm(\alpha_1,\ldots,\alpha_N)\rangle = \mathcal{N}_N^\pm \big( |\alpha_1,\ldots,\alpha_N\rangle \pm |-\alpha_1,\ldots,-\alpha_N\rangle \big)$$7 in current experiments (Oudot et al., 2014).

• Multipartite Entanglement Witnesses: Localizable negativity between widely separated modes, measured via hybrid tomography, serves as a robust witness up to at least four modes, while reconstructed joint Wigner functions reveal nonclassical interference structure (Wang et al., 2021, Hoshi et al., 2024).
• Hybrid Quantifiers: Full joint Wigner functions for hybrid discrete-continuous systems (e.g., qutrit-boson) capture both the nonclassical fringe patterns and the underlying combinatorial entanglement structure (Lotkov et al., 10 Sep 2025).

4. Experimental Demonstrations and State Tomography

Deterministic Preparation: Reflection from a dispersive superconducting cavity has enabled one- to four-mode photonic Schrödinger-cat states with direct homodyne tomography; reconstructed fidelities reach 0.75 (single-mode even), 0.54 (tripartite, even), and negativities confirm DV–CV entanglement (Wang et al., 2021). Planar Kerr parametric oscillator chips enable conversion between DV-encoded and CV-encoded Bell states, and implement entangling $$|\mathrm{Cat}_N^\pm(\alpha_1,\ldots,\alpha_N)\rangle = \mathcal{N}_N^\pm \big( |\alpha_1,\ldots,\alpha_N\rangle \pm |-\alpha_1,\ldots,-\alpha_N\rangle \big)$$8 gates on cat qubits; fidelities for Bell–Cat states post-gate are $$|\mathrm{Cat}_N^\pm(\alpha_1,\ldots,\alpha_N)\rangle = \mathcal{N}_N^\pm \big( |\alpha_1,\ldots,\alpha_N\rangle \pm |-\alpha_1,\ldots,-\alpha_N\rangle \big)$$9 (Hoshi et al., 2024).

Positive-P Simulations: Stochastic phase-space techniques facilitate the numerical analysis of transient cat formation and nonlocal spatial correlations ($$|\alpha_1,\ldots,\alpha_N\rangle = \bigotimes_{j=1}^N|\alpha_j\rangle$$0, $$|\alpha_1,\ldots,\alpha_N\rangle = \bigotimes_{j=1}^N|\alpha_j\rangle$$1 functions), up to $$|\alpha_1,\ldots,\alpha_N\rangle = \bigotimes_{j=1}^N|\alpha_j\rangle$$2 in one-dimensional resonator chains, with computational cost scaling only linearly in system size (Shi et al., 11 Jan 2026).

Bell-State and GHZ-State Engineering: In driven Kerr oscillator networks, adiabatic/diabatic protocols directly realize maximally entangled Bell-resonator-cats and their $$|\alpha_1,\ldots,\alpha_N\rangle = \bigotimes_{j=1}^N|\alpha_j\rangle$$3-mode GHZ generalizations with extremely high fidelity and operational speed (ns–μs regime) (Resch et al., 5 May 2025).

Hybrid Systems: $$|\alpha_1,\ldots,\alpha_N\rangle = \bigotimes_{j=1}^N|\alpha_j\rangle$$4-symmetric Rabi models yield three-component (qutrit-boson) entangled cat states with phase-space structure observable via joint Wigner tomography; experimental realizations are feasible in superconducting qudit-QED or optomechanical settings (Lotkov et al., 10 Sep 2025).

5. Scalability, Applications, and Limitations

Scalability: Adiabaticity requirements for ground-state cat manifold preparation in Kerr-type systems exert only mild scaling burden: the relevant gap remains $$|\alpha_1,\ldots,\alpha_N\rangle = \bigotimes_{j=1}^N|\alpha_j\rangle$$5 (Kerr strength), independent of mode number. Experimental error modeling projects multipartite cat preparation to $$|\alpha_1,\ldots,\alpha_N\rangle = \bigotimes_{j=1}^N|\alpha_j\rangle$$6 modes (for $$|\alpha_1,\ldots,\alpha_N\rangle = \bigotimes_{j=1}^N|\alpha_j\rangle$$7) with quantum state fidelity $$|\alpha_1,\ldots,\alpha_N\rangle = \bigotimes_{j=1}^N|\alpha_j\rangle$$8 given reductions in cavity losses and improved $$|\alpha_1,\ldots,\alpha_N\rangle = \bigotimes_{j=1}^N|\alpha_j\rangle$$9 decoupling (Wang et al., 2021, Resch et al., 5 May 2025).

Quantum Information Processing: Multimode cat states underpin error-biased quantum codes (cat codes), quantum networking via entangled flying cats, and architectures for continuous-variable logical qubits supporting fault-tolerant quantum computing via CV-encoded gate sets (Resch et al., 5 May 2025, Hoshi et al., 2024, Wang et al., 2021).

Quantum Metrology: Large macroscopicity (QFI scaling $$\mathcal{N}_N^\pm = \left[2 \pm 2 e^{-2 \sum_j |\alpha_j|^2}\right]^{-1/2}.$$0) implies multiparameter sensitivity enhancement, with applications to phase estimation and weak-force detection (Oudot et al., 2014, Wang et al., 2021).

Limitations: Main limiting factors are photon loss, dephasing, experimental state-preparation errors, and, in simulation, positive-P boundary-term instabilities that predominantly affect high-order parity and Wigner negativity estimation—though low-order moments and spatial correlations can be extracted robustly (Shi et al., 11 Jan 2026).

Scaling Challenges: In hardware, spectral crowding and pump crosstalk complicate isolation of distinct KPOs; tomography overhead grows exponentially with mode number, necessitating compressed-sensing strategies or partial marginal analysis (Hoshi et al., 2024).

6. Extensions to Novel Symmetries and Physical Regimes

Beyond the standard bimodal (even/odd) coherent superpositions, multimode Schrödinger cats emerge naturally in models endowed with higher discrete symmetries:

• $$\mathcal{N}_N^\pm = \left[2 \pm 2 e^{-2 \sum_j |\alpha_j|^2}\right]^{-1/2}.$$1-Symmetric Cat States: Deep-strong-coupling in generalized Rabi models produces $$\mathcal{N}_N^\pm = \left[2 \pm 2 e^{-2 \sum_j |\alpha_j|^2}\right]^{-1/2}.$$2-component superpositions in phase space, with symmetry-selected interference patterns, as characterized by closed-form hybrid Wigner functions (Lotkov et al., 10 Sep 2025).
• Hybridized and Dissipatively Stabilized Cat States: Combining engineered dissipation (e.g., two-photon loss) and Hamiltonians protects cat-state manifolds and supports autonomous error correction. Designs for "pair-cat" codes and autonomously stabilized cat qubits represent active research, with both theoretical and experimental progress underway.
• Hybrid DV–CV Entanglement Networks: Extending DV–CV gate sets (universal on cat qubits) and graph state protocols enables larger, fault-tolerant multimode continuous-variable architectures, particularly when supported by fast, high-fidelity entangling gates (e.g., $$\mathcal{N}_N^\pm = \left[2 \pm 2 e^{-2 \sum_j |\alpha_j|^2}\right]^{-1/2}.$$3) (Hoshi et al., 2024, Wang et al., 2021, Resch et al., 5 May 2025).

---

In summary, multimode Schrödinger cat states form a mathematically and physically rich class of macroscopic quantum superpositions involving multiple bosonic modes, with diverse experimental realizations, robust macroscopicity measures, broad applicability to quantum information protocols, and a wide range of structural extensions involving engineered interactions, symmetry, and hybridization with discrete-variable systems (Oudot et al., 2014, Wang et al., 2021, Resch et al., 5 May 2025, Shi et al., 11 Jan 2026, Hoshi et al., 2024, Lotkov et al., 10 Sep 2025).