We provide a framework for understanding dynamical metastability in open many-body systems of free bosons, whereby the dynamical stability properties of the system in the infinite-size (thermodynamic) limit may sharply differ from those of any finite-size truncation, and anomalous transient dynamics may arise. By leveraging pseudospectral techniques, we trace the discrepancy between asymptotic and transient dynamics to the non-normality of the underlying quadratic bosonic Lindbladian (QBL) generator and show that two distinct flavors of dynamical metastability can arise. QBLs exhibiting type I dynamical metastability, previously discussed in the context of anomalous transient amplification [Phys. Rev. Lett. 127, 245701 (2021)], are dynamically unstable in the infinite-size limit yet stable once open boundaries are imposed. In contrast, type II dynamically metastable QBLs, which we uncover in this work, are dynamically stable for infinite size but become unstable under open boundary conditions for arbitrary finite system size. We exhibit representative models for both types of metastability in the dissipative, as well as the closed-system (Hamiltonian) settings, and analyze distinctive behavior they can engender. We show that dynamical metastability manifests itself in the generation of entanglement entropy by way of a transient which reflects the stability phase of the infinite (rather than the actual finite) system and, as a result, is directly tied to the emergence of supervolume entanglement scaling in type I systems. Finally, we demonstrate how, even in Hermitian, and especially in highly non-normal regimes, the spectral properties of an infinite-size QBL are reflected in the linear response functions of the corresponding finite QBLs by way of resonant pseudospectral modes.
We provide a framework for understanding dynamical metastability in open many-body systems of free bosons, whereby the dynamical stability properties of the system in the infinite-size (thermodynamic) limit may sharply differ from those of any finite-size truncation, and anomalous transient dynamics may arise. By leveraging pseudospectral techniques, we trace the discrepancy between asymptotic and transient dynamics to the non-normality of the underlying quadratic bosonic Lindbladian (QBL) generator and show that two distinct flavors of dynamical metastability can arise. QBLs exhibiting type I dynamical metastability, previously discussed in the context of anomalous transient amplification [Phys. Rev. Lett. 127, 245701 (2021)], are dynamically unstable in the infinite-size limit yet stable once open boundaries are imposed. In contrast, type II dynamically metastable QBLs, which we uncover in this work, are dynamically stable for infinite size but become unstable under open boundary conditions for arbitrary finite system size. We exhibit representative models for both types of metastability in the dissipative, as well as the closed-system (Hamiltonian) settings, and analyze distinctive behavior they can engender. We show that dynamical metastability manifests itself in the generation of entanglement entropy by way of a transient which reflects the stability phase of the infinite (rather than the actual finite) system and, as a result, is directly tied to the emergence of supervolume entanglement scaling in type I systems. Finally, we demonstrate how, even in Hermitian, and especially in highly non-normal regimes, the spectral properties of an infinite-size QBL are reflected in the linear response functions of the corresponding finite QBLs by way of resonant pseudospectral modes.
Tight bosonic analogs of free-fermionic symmetry-protected topological phases, and their associated edge-localized excitations, have long evaded the grasp of condensed-matter and AMO physics. In this paper, building on our initial exploration [Phys. Rev. Lett. 127, 245701 (2021)], we identify a broad class of quadratic bosonic systems subject to Markovian dissipation that realize tight bosonic analogs of the Majorana and Dirac edge modes characteristic of topological superconductors and insulators, respectively. To this end, we establish a general framework for topological metastability for these systems, by leveraging pseudospectral theory as the appropriate mathematical tool for capturing the nonnormality of the Lindbladian generator. The resulting dynamical paradigm, which is characterized by both a sharp separation between transient and asymptotic dynamics and a nontrivial topological invariant, is shown to host edge-localized modes, which we dub Majorana and Dirac bosons. Generically, such modes consist of one conserved mode and a canonically conjugate generator of an approximate phase-space translation symmetry of the dynamics. The general theory is exemplified through several representative models exhibiting the full range of exotic boundary physics that topologically metastable systems can engender. In particular, we explore the extent to which Noether's theorem is violated in this dissipative setting and the way in which certain symmetries can nontrivially modify the edge modes. Notably, we also demonstrate the possibility of anomalous parity dynamics for a bosonic cat state prepared in a topologically metastable system, whereby an equal distribution between even and odd parity sectors is sustained over a long transient. For both Majorana and Dirac bosons, observable multitime signatures in the form of anomalously long-lived quantum correlations and divergent zero-frequency power spectral peaks are proposed and discussed in detail. Our results point to a paradigm for symmetry-protected topological physics in free bosons, embedded deeply in the long-lived transient regimes of metastable dynamics.
Majorana bosons, that is, tight bosonic analogues of the Majorana fermionic quasi-particles of condensed-matter physics, are forbidden for gapped free bosonic matter within a standard Hamiltonian scenario. We show how the interplay between dynamical metastability and non-trivial bulk topology makes their emergence possible in non-interacting bosonic chains undergoing Markovian dissipation. This leads to a distinctive form of topological metastability whereby, consistent with the breakdown of Noether's theorem for open systems, a conserved Majorana boson localized on one edge is paired with a distinct symmetry generator localized on the opposite edge. We argue that Majorana bosons are robust against disorder and identifiable by signatures in the zero-frequency steady-state power spectrum. Our results suggest that symmetry-protected topological phases for free bosons may arise in transient metastable regimes, which are still stable over practical time scales.
For a many-body system in equilibrium, described by a thermodynamically stable Hamiltonian, quantum criticality is associated with structural changes of the many-body ground state. However, there exist physically relevant models, notably, certain quadratic bosonic Hamiltonians (QBHs), which fail to have a ground state. QBHs can be dynamically stable or unstable. We show the notion of criticality is meaningful for the entire class of QBHs that are dynamically stable or at the boundary of instability, regardless of thermodynamic stability, and that the key state for such QBHs is a naturally and unambiguously defined quasiparticle vacuum (QPV). This state is Gaussian, and coincides with the ground state if the QBH is thermodynamically stable. We identify a relevant spectral gap, the Krein gap, associated to the minimal spectral separation between creation and annihilation operators, and show that the QPV is unique when the Krein gap is positive. We prove that, for dynamically stable QBHs with finite-range couplings, correlations are exponentially bounded unless the Krein gap closes, which is associated with one of two spectral degeneracies: an exceptional point or a Krein collision. Consequently, long-range QPV correlations can ensue. Thus, the Krein gap takes the role of the spectral gap for dynamically stable QBHs, and the boundary of dynamical stability and criticality (associated to exceptional points) or multicriticality (associated to Krein collisions) are the same. We also find that bosonic critical behavior beyond thermodynamic stability is witnessed by the scaling of the entanglement entropy and other indicators of equilibrium criticality from information geometry. Our framework opens the door to investigating all dynamically stable QBHs through the lens of critical phenomena, including thermodynamically unstable ones from photonics, cavity-QED, and magnonics.
Abstract Unlike their fermionic counterparts, the dynamics of Hermitian quadratic bosonic Hamiltonians are governed by a generally non-Hermitian Bogoliubov-de Gennes effective Hamiltonian. This underlying non-Hermiticity gives rise to a dynamically stable regime, whereby all observables undergo bounded evolution in time, and a dynamically unstable one, whereby evolution is unbounded for at least some observables. We show that stability-to-instability transitions may be classified in terms of a suitably generalized P T symmetry , which can be broken when diagonalizability is lost at exceptional points in parameter space, but also when degenerate real eigenvalues split off the real axis while the system remains diagonalizable. By leveraging tools from Krein stability theory in indefinite inner-product spaces, we introduce an indicator of stability phase transitions, which naturally extends the notion of phase rigidity from non-Hermitian quantum mechanics to the bosonic setting. As a paradigmatic example, we fully characterize the stability phase diagram of a bosonic analogue to the Kitaev–Majorana chain under a wide class of boundary conditions. In particular, we establish a connection between phase-dependent transport properties and the onset of instability, and argue that stable regions in parameter space become of measure zero in the thermodynamic limit. Our analysis also reveals that boundary conditions that support Majorana zero modes in the fermionic Kitaev chain are precisely the same that support stability in the bosonic chain.
Metastability, i.e., partial relaxation to long-lived, quasi-stationary states before true asymptotic equilibrium sets in, emerges ubiquitously in classical and quantum dynamical systems as a result of timescales separation. In open quantum systems, an intrinsically nonequilibrium analogue, dynamical metastability, can originate from the spectral geometry of a non-Hermitian operator. In noninteracting models, this mechanism produces boundary-sensitive anomalous relaxation, transient amplification, and topologically mandated long-lived edge modes, all of which are enhanced as system size grows. Here we extend dynamical metastability into the nonlinear, interacting regime and identify magnetic heterostructures as a natural platform for its exploration. We introduce an interacting spin Lindbladian whose linearized magnon dynamics map onto a dynamically metastable Hatano-Nelson chain, and show that dynamical metastability in the noninteracting limit seeds genuinely nonlinear phenomena, including size-dependent spin dipping and anomalous attraction to unstable equilibria. Long-lived edge states associated to topologically mandated Dirac bosons persist under nonlinearities and disorder. We further analyze the magnetization dynamics in magnetic multilayers within the classical Landau-Lifshitz-Gilbert-Slonczewski framework, identifying Dzyaloshinskii-Moriya interaction, nonlocal damping, and spin-transfer torque as control parameters governing bulk-boundary stability mismatch and band topology. While all the distinctive dynamical phenomena previously identified reappear in this experimentally relevant setting, the LLGS framework also supports multistability and limit cycles that are absent in the quantum model. Our results constitute the first systematic study of dynamical metastability in nonlinear dynamics, directly relevant to spin-torque oscillator arrays, magnonic devices, and beyond.
A new understanding of criticality in systems described by quadratic bosonic Hamiltonians (QBHs) ties the emergence of long-range correlations to boundaries of dynamical, not thermodynamical, stability in parameter space. This separation occurs because the solution of the Heisenberg equations of motion is determined by an auxiliary pseudo-Hermitian dynamical system. The boundary points of a region of dynamical stability can be either exceptional points, generically associated with long-range correlations, or Krein collisions, where correlations can be either long- or short-range. We investigate the interplay of this landscape of possibilities with band topology and boundary physics, by relying on both specific examples and general arguments. The examples stem from a two-parameter, thermodynamically unstable family of QBHs obtained from the bosonic Su-Schrieffer-Heeger model by breaking particle conservation while preserving a chiral pseudo-symmetry. As a function of an interpolation parameter, distinct regions emerge within a fixed dynamical-stability phase, which prove to be topologically trivial and nontrivial, respectively. The topological phase transition is a line of Krein collisions, which coincides with the closing of a band gap at zero and causes the localization length of the topologically mandated boundary zero modes to diverge before disappearing. We show that the chiral pseudo-symmetry induces, despite the broken particle-number symmetry, enough structure on its associated dynamical matrices to support a topological classification and a bulk-boundary correspondence, independently of dynamical stability. This strongly suggests that, for non-interacting bosons, topological physics extracted from basic index theory is insensitive to dynamical stability and, a posteriori, criticality.