Quantum limit cycles with continuous symmetries from coherent parametric driving: exact solutions and many-body extensions

Abstract: There is widespread interest in many-body quantum systems that exhibit limit-cycle or time-crystalline behaviour. An ideal quantum limit cycle would be realized using fully coherent driving (to minimize noise) and also have a continuous internal symmetry (to ensure generation of monochromatic radiation). While these two requirements may seem incompatible, we introduce in this work a large class of multi-mode bosonic limit cycle models based on coherent parametric driving which possess an O(N) continuous symmetry. Surprisingly, the full quantum dissipative steady state of these models can be found exactly. They exhibit rich physics, including steady state entanglement, reduced phase diffusion and the possibility of realizing quantum limit tori. The basic mechanism we identify provides a unified way to understand how coherent parametric driving can yield symmetry-enriched limit cycles, and also helps us understand related models where the relevant symmetries are weakly broken. The models we study are compatible with a range of different experimental platforms, including quantum optical setups and superconducting quantum circuits.

• The paper presents an exact analytic solution for the nonequilibrium steady states of multi-mode bosonic systems with continuous O(N) symmetry.
• It demonstrates that coherent parametric driving induces robust quantum limit cycles and limit tori, reducing phase noise below the traditional Schawlow–Townes limit.
• The work extends quantum many-body models by enabling precise control of non-Gaussian entanglement and symmetry-broken dynamics in open quantum systems.

Quantum Limit Cycles with Continuous Symmetries: Exact Solutions and Many-Body Extensions

Introduction and Motivation

The paper "Quantum limit cycles with continuous symmetries from coherent parametric driving: exact solutions and many-body extensions" (2604.25864) identifies a route to realizing dissipative quantum limit cycles that simultaneously feature continuous internal symmetries and the absence of incoherent driving. This development overcomes a prevailing assumption that O(N)- or U(1)-symmetric quantum limit cycles necessitate additional environmental dissipation, as encountered in standard lasers or the quantum van der Pol oscillator. The result is a formulation of a general class of multi-mode bosonic systems, driven coherently via parametric drives, exhibiting O(N) symmetry, and uniquely admitting exact analytic solutions for their nonequilibrium steady-state (NESS), together with a rich set of dynamical and entanglement properties.

Model Construction and Symmetry Considerations

The core system comprises $N$ bosonic modes, each subject to identical Kerr-type nonlinearities and coherent two-photon (parametric) driving. Critically, while a single mode with parametric drive only has a residual $\mathbb{Z}_2$ symmetry, the multi-mode variant retains an O(N) symmetry arising from invariance under orthogonal rotations in mode space when nonlinearities and drive amplitudes are perfectly balanced. For two modes ($N=2$), this symmetry is SO(2), equivalent to the familiar U(1). An antisymmetric hopping term, proportional to an arbitrary real antisymmetric matrix $K$, couples the modes but leaves the steady state invariant due to commutation with the O(N) generator.

Figure 1: Minimal two-mode realization with self- and cross-Kerr interactions, antisymmetric hopping, coherent parametric drive, and photon loss; (b) displays the steady-state Wigner and joint quadrature distributions forming a ring in phase space.

This construction forms the basis for a new class of quantum limit-cycle models wherein the interplay between coherent driving, nonlinearity, and dissipation generates persistent oscillatory dynamics with continuous symmetry, eliminating the additional quantum noise from incoherent pumping.

Semiclassical and Quantum Dynamics: Phase Diagram

In the semiclassical limit (large occupation), the system exhibits the emergence of a stable attractor manifold defined by the constraint $x_1^2 + x_2^2 = n_{ss}$ (for $N=2$), which forms a ring in phase space. The dynamics on this ring is governed by a Goldstone mode associated with the broken continuous symmetry, with persistent phase rotation at a frequency set by the antisymmetric coupling. For general $N$, the attractor generalizes to an $S^{N-1}$ sphere, with $N-1$ gapless directions.

A remarkable advance here is the exact analytic solution for the NESS at arbitrary photon number and for all $N$. The system shows a nonequilibrium phase transition at a critical drive strength $\mathbb{Z}_2$0 (with $\mathbb{Z}_2$1 the dissipation rate): below this, the steady state is trivial; above it, there is symmetry breaking and limit-cycle or toroidal dynamics.

Figure 2: (a) Phase diagram showing distinct symmetry-broken and symmetric phases; (b) order parameter scaling and transition; (c) steady-state entanglement; (d) photon-number Fano ratio versus drive.

Deep in the limit-cycle phase, the steady-state photon number, entanglement, and amplitude fluctuations are all computable from the exact NESS. The transition is second order with critical exponent $\mathbb{Z}_2$2, consistent with mean-field expectations.

Steady-State Entanglement and Non-Gaussianity

Even with large occupation, the quantum steady state exhibits persistent, non-Gaussian intermode entanglement. For $\mathbb{Z}_2$3, log-negativity $\mathbb{Z}_2$4 provides a quantifiable measure; notably, $\mathbb{Z}_2$5 increases with drive strength before saturating at a finite value that, while smaller than an equivalent thermal two-mode squeezed state (TMSS), is nevertheless nonzero and robust in the semiclassical regime. All intermode Gaussian covariances vanish, indicating the entanglement derives exclusively from higher-order correlations.

For general bipartitions (e.g., in the $\mathbb{Z}_2$6 normal mode basis), the behavior differs: $\mathbb{Z}_2$7 peaks near threshold and saturates at a lower value. This highlights basis dependence and the critical role of non-Gaussianity in the steady-state structure.

Figure 3: Log-negativity and its ratio to a thermal two-mode squeezed state in the $\mathbb{Z}_2$8 partition.

Phase Diffusion and Quantum Noise Suppression

Phase diffusion along the attractor manifold arises from quantum noise. Critically, the use of coherent driving suppresses pump-related noise, resulting in a phase diffusion constant $\mathbb{Z}_2$9—and hence laser linewidth— which can be reduced by a factor of two compared to the traditional Schawlow–Townes limit. Specifically, in the vicinity above threshold, the minimized linewidth is $N=2$0. Amplitude (photon number) fluctuations, as quantified by the Fano ratio, increase near threshold, but deep in the classical regime achieve nearly Poissonian statistics.

Higher-Dimensional Dynamics and Quantum Limit Tori

For $N=2$1, the antisymmetric coupling introduces multiple rotational axes on the attractor manifold, resulting in persistent oscillations with generally incommensurate frequencies: quantum limit tori (rather than limit cycles), whose geometry is set by the ranks and structure of $N=2$2.

Figure 4: (a) Attractor manifold for $N=2$3 in the block-diagonalized basis; (b) stereographic projection of a Clifford torus and illustration of quasiperiodic versus periodic orbits.

In the presence of quantum noise, there is slow diffusion both along angular (phase) directions and among amplitudes (radius partitioning among the $N=2$4 rotation planes), described by a full diffusion matrix on $N=2$5. This generalizes the single phase-diffusive Goldstone mode of $N=2$6 to a richer many-body setting.

Robustness and Symmetry-Breaking Perturbations

Small symmetry-breaking perturbations (e.g., imbalanced Kerr nonlinearities) deform the O(N) (or U(1)) symmetry, generating phase-locking terms in the effective equations of motion. The system then shows a finite range of persistent oscillatory (limit cycle/tori) behavior, with phase-locking dominating at high drive strength or strong symmetry breaking. Thus, the symmetric coherently driven limit cycle functions as a parent model for a large class of previously studied parametrically driven systems, including those motivated by time-crystal and supersolidity physics.

Practical and Theoretical Implications

This work drastically expands the set of tractable, exactly-solvable, non-equilibrium quantum models exhibiting both strong quantum correlations (e.g., entanglement, non-Gaussianity) and continuous dynamical symmetry. It clarifies how coherent parametric driving—ubiquitous in contemporary platforms such as superconducting circuits and quantum optics—supports symmetry-protected limit cycles and higher-dimensional generalizations. The suppression of quantum phase noise below the Schawlow–Townes bound is particularly relevant for quantum metrology and the development of ultrastable light sources, including those with strong multimode quantum entanglement.

The methodology offers a blueprint for engineering coherent, multimode dissipative quantum phenomena, with immediate implications for quantum synchronization, phase entrainment, and controlled generation of multimode non-classical light.

Conclusion

Coherently driven, O(N) symmetric bosonic systems enable the coexistence of continuous symmetry, robust limit cycle or toroidal dynamics, and reduced quantum noise, in a framework admitting exact analytic treatment. The interplay of symmetry, nonlinearity, and coherent dissipation sets new paradigms for open quantum system engineering—applicable across quantum optics, quantum information, and out-of-equilibrium many-body physics. Future directions include the analysis of quantum synchronization, exploration of output field quantum correlations, and extensions to non-Gaussian measurement and feedback protocols.