Metastability in open quantum systems is usually inferred from spectral separation in the Liouvillian, which governs unconditional, ensemble-averaged dynamics. We show that this diagnosis is incomplete at the level of individual quantum trajectories: conditioned realizations of the same unconditional dynamics can bypass, transiently access, or operationally preserve a metastable memory, depending on the monitored channel and the observed record. We demonstrate these mechanisms in a driven-dissipative nonlinear oscillator realizing a quantum associative memory, by comparing spectra, target fidelities, and phase-space distributions. The non-Hermitian dynamics obtained by post-selecting on the absence of detected events supports metastable retrieval, but with a distinct long-time fate: normalization selects the least-decaying mode of the non-Hermitian spectrum rather than the addressed memory branch. In contrast, stochastic jump trajectories can preserve retrieval over extended times when the post-measurement update remains compatible with the coherent memory structure. Thus trajectory-level metastability is not determined by the averaged generator alone; it requires compatibility between the monitored channel, the measurement record, and the metastable manifold.
Metastability in open quantum systems is usually inferred from spectral separation in the Liouvillian, which governs unconditional, ensemble-averaged dynamics. We show that this diagnosis is incomplete at the level of individual quantum trajectories: conditioned realizations of the same unconditional dynamics can bypass, transiently access, or operationally preserve a metastable memory, depending on the monitored channel and the observed record. We demonstrate these mechanisms in a driven-dissipative nonlinear oscillator realizing a quantum associative memory, by comparing spectra, target fidelities, and phase-space distributions. The non-Hermitian dynamics obtained by post-selecting on the absence of detected events supports metastable retrieval, but with a distinct long-time fate: normalization selects the least-decaying mode of the non-Hermitian spectrum rather than the addressed memory branch. In contrast, stochastic jump trajectories can preserve retrieval over extended times when the post-measurement update remains compatible with the coherent memory structure. Thus trajectory-level metastability is not determined by the averaged generator alone; it requires compatibility between the monitored channel, the measurement record, and the metastable manifold.
Wavefunction collapse is commonly associated with unavoidable physical disturbance of the measured system. Here we show that in driven-dissipative quantum systems, continuous measurement can induce strong trajectory-level collapse while leaving the ensemble-averaged steady state strictly invariant. We identify measurement-invariant steady states whose unconditional density matrix remains unchanged under continuous monitoring, despite pronounced measurement-induced localization in conditioned quantum trajectories. This separation between trajectory-level collapse and ensemble invariance defines a regime of partial collapse, in which measurement-induced localization is continuously counteracted by dissipative dynamics. We derive a necessary and sufficient condition for steady-state invariance under continuous measurement and identify Liouvillian symmetry as a concrete dynamical mechanism enforcing it. Our results clarify the distinction between conditional collapse and physical disturbance in open quantum systems and provide a framework for non-invasive continuous monitoring in driven-dissipative settings.
This thesis considers quantum systems which interact with their surrounding environment. Specifically we consider the underlying stochastic dynamics which arise from the back-action of measurements performed in the environment. These quantum trajectories can be tracked by continuously monitoring the environment, as is now possible in modern experiments. We consider phenomena understood at the average level of the system state density matrix, and explore these at this more information rich level. This is especially relevant for experimental setups and implementation of quantum technologies. Firstly, we consider how systems exhibiting metastability behave at the level of stochastic quantum trajectories. In doing so we make comparisons to classical theories of metastability, where trajectories exhibit fast relaxation into distinct phases and slow transitions between them. We also illustrate when this phenomenology does not hold, with instead the emergence of a slow manifold in the form of a decoherence free subspace. We also highlight the useful role that quantum reset processes can play in the analysis of stochastic quantum trajectories. The dynamics of the stochastic quantum trajectories depends on the choice of environmental measurements. There are many choices which leave the average dynamics unchanged, but give different statistical ensembles of quantum trajectories. This gauge freedom of the average dynamics is well known. We derive the analogous result at the level of quantum trajectory ensembles; giving the measurement protocols which leave the generator of the stochastic trajectory dynamics unchanged. We also extend this to include joint trajectories of conditional system states and their associated measurement records, either full records or coarse-grained. Using these results on gauge freedoms of quantum trajectories, we consider the conditions under which the trajectory dynamics exhibit a weak unitary symmetry. We then consider under what conditions the joint state of the system and environment exhibit symmetric dynamics, including under various levels of dephasing. We find that the joint dynamics always features a separable symmetry directly related to that of the quantum master equation, but for quantum trajectories the corresponding symmetry is present only if the counting measurement satisfies certain conditions.
Distinct monitoring or intervention schemes can produce different conditioned stochastic quantum trajectories while sharing the same unconditional (ensemble-averaged) dynamics. This is the essence of unravelings of a given Gorini-Kossakowski-Sudarshan-Lindblad (GKSL) master equation: any trajectory-ensemble average of a function that is linear in the conditional state is completely determined by the unconditional density matrix, whereas applying a nonlinear function before averaging can yield unraveling-dependent results beyond the average evolution. A paradigmatic example is resonance fluorescence, where direct photodetection (jump/Poisson) and homodyne or heterodyne detection (diffusive/Wiener) define inequivalent unravelings of the same GKSL dynamics. In earlier work, we showed that nonlinear trajectory averages can distinguish such unravelings, but observing the effect in that optical setting requires demanding experimental precision. Here we translate the same idea to a digital setting by introducing synthetic unravelings implemented as quantum circuits acting on one and two qubits. We design two unravelings - a projective measurement unraveling and a random-unitary "kick" unraveling - that share the same ensemble-averaged evolution while yielding different nonlinear conditional-state statistics. We implement the protocols on superconducting-qubit hardware provided by IBM Quantum to access trajectory-level information. We show that the variance across trajectories and the ensemble-averaged von Neumann entropy distinguish the unravelings in both theory and experiment, while the unconditional state and the ensemble-averaged expectation values that are linear in the state remain identical. Our results provide an accessible demonstration that quantum trajectories encode information about measurement backaction beyond what is fixed by the unconditional dynamics.
We introduce jump-time unraveling as a distinct description of open quantum systems. As our starting point, we consider quantum jump trajectories, which emerge, physically, from continuous quantum measurements, or, formally, from the unraveling of Markovian quantum master equations. If the stochastically evolving quantum trajectories are ensemble-averaged at specific times, the resulting quantum states are solutions to the associated quantum master equation. We demonstrate that quantum trajectories can also be ensemble-averaged at specific jump counts. The resulting jump-time-averaged quantum states are then solutions to a discrete, deterministic evolution equation, with time replaced by the jump count. This jump-time evolution represents a trace-preserving quantum dynamical map if and only if the associated open system does not exhibit dark states. In the presence of dark states, on the other hand, jump-time-averaged states may decay into the dark states and the jump-time evolution may eventually terminate. Jump-time-averaged quantum states and the associated jump-time evolution are operationally accessible in continuous measurement schemes, when quantum jumps are detected and used to trigger the readout measurements. We illustrate the jump-time evolution with the examples of a two-level system undergoing relaxation or dephasing, a damped harmonic oscillator, and a free particle exposed to collisional decoherence.
Open quantum systems can display metastability, with multiple long-lived states coexisting. Examples include the bistability of bright and dim manifolds in nonlinear driven resonators, or the long-lived coherent states of opposite phases in parametrically-driven systems, the building block of cat-state-based quantum computers. Metastability often emerges at system sizes where density-matrix evolution is impractical. Quantum trajectories mitigate this issue by evolving wavefunctions rather than density matrices. Metastable states then emerge as attractors towards which the system rapidly relaxes until a rare switching event brings it to another attractor. Simulating them remains costly, as one must resolve fast intra-attractor dynamics over timescales long enough to gather switching statistics. To tackle this problem, we introduce the Diagonal, Switching, and Locally Optimal Unraveling (DiSLOU) method, which combines three optimization layers: (i) Switching among locally optimal unravelings that minimize the number of quantum jumps in each metastable region. (ii) Diagonal propagation: the no-jump evolution is accelerated by diagonalizing the non-Hermitian effective Hamiltonian once per attractor. (iii) Local reduced eigenspaces: the no-jump state can be faithfully propagated in a truncated subspace of only slowly decaying eigenmodes. We benchmark DiSLOU on the driven-dissipative Kerr resonator, finding no measurable loss of accuracy and the largest speedup in the bistable regime, where standard trajectories are most expensive. We then apply the method to the more challenging problem of accurately simulating rare bit-flip events in the cat manifold emerging in a two-mode dissipatively stabilized parametric oscillator, obtaining a performance improvement of over one order of magnitude. An implementation of DiSLOU in Julia is available as an open-source package.
We present a general theory of classical metastability in open quantum systems. Metastability is a consequence of a large separation in timescales in the dynamics, leading to the existence of a regime when states of the system appear stationary, before eventual relaxation toward a true stationary state at much larger times. In this work, we focus on the emergence of classical metastability, i.e., when metastable states of an open quantum system with separation of timescales can be approximated as probabilistic mixtures of a finite number of states. We find that a number of classical features follow from this approximation, for the manifold of metastable states, long-time dynamics between them, and symmetries of the dynamics. Namely, those states are approximately disjoint and thus play the role of metastable phases, the relaxation toward the stationary state is approximated by a classical stochastic dynamics between them, and weak symmetries correspond to their permutations. Importantly, the classical dynamics is observed not only on average but also at the level of individual quantum trajectories: We show that time coarse-grained continuous measurement records can be viewed as noisy classical trajectories, while their statistics can be approximated by that of the classical dynamics. Among others, this explains how first-order dynamical phase transitions arise from metastability. Finally, to verify the presence of classical metastability in a given open quantum system, we develop an efficient numerical approach that delivers the set of metastable phases together with the effective classical dynamics. Since the proximity to a first-order dissipative phase transition manifests as metastability, the theory and tools introduced in this work can be used to investigate such transitions---which occur in the large size limit---through the metastable behavior of many-body systems of moderate sizes accessible to numerics.