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.
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.
The Gorini-Kossakowski-Sudarshan-Lindblad master equation (ME) governs the density matrix of open quantum systems (OQSs). When an OQS is subjected to weak continuous measurement, its state evolves as a stochastic quantum trajectory, whose statistical average solves the ME. The ensemble of such trajectories is termed an unraveling of the ME. We propose a method to operationally distinguish unravelings produced by the same ME in different measurement scenarios, using nonlinear averages of observables over trajectories. We apply the method to the paradigmatic quantum nonlinear system of resonance fluorescence in a two-level atom. We compare the Poisson-type unraveling, induced by direct detection of photons scattered from the two-level emitter, and the Wiener-type unraveling, induced by phase-sensitive detection of the emitted field. We show that a quantum-trajectory-averaged variance is able to distinguish these measurement scenarios. We evaluate the performance of the method, which can be readily extended to more complex OQSs, under a range of realistic experimental conditions. Published by the American Physical Society 2024
The quantum-trajectory description of a noisy quantum circuit is not unique: every Kraus representation of the noise channel defines a different unraveling, i.e., a different ensemble of pure states that averages to the same density matrix. While all unravelings share the same expectation values of observables, they generally differ in the amount of quantum resources carried by the individual trajectories. Here we study non-stabilizerness, or magic, quantified by the stabilizer Rényi entropy, in quantum-trajectory unravelings of Haar-random and T-doped random Clifford brickwork circuits subject to local dephasing noise of strength p. The canonical unraveling in terms of Pauli jumps is Clifford and leaves magic untouched, so that trajectories retain the Haar-random value of non-stabilizerness. Exploring the full two-parameter family of minimal unravelings, we find that the trajectory-averaged magic can be substantially reduced by a state-independent choice of Kraus operators, and that the optimal choice is not in general the projective-measurement unraveling, but a two-outcome generalized measurement in which the Kraus operators correspond to a weak measurement of the noisy qubit. The optimal unraveling reduces the extensive magic density of the trajectories by a p-dependent factor without altering its linear scaling with system size. At the same time, it produces a very broad distribution of magic over trajectories, with a sizable fraction of trajectories retaining Haar-random magic. Since the cost of stabilizer-based classical simulation grows exponentially with the magic of the state, these high-magic trajectories, rather than the average, control the simulation cost. Our results delineate the potential and the limitations of unraveling optimization as a route to the classical simulation of noisy quantum dynamics with stabilizer methods.
Stochastic unravelings of Lindblad-type master equations, such as stochastic Schrödinger equations, provide powerful tools to model open quantum systems and continuous measurement processes. The same master equation can be unraveled in different ways; while these unravelings differ at the level of quantum trajectories, by construction they all yield the same averaged dynamics for the density operator. A recent question of both foundational and practical relevance is whether such unravelings can be operationally distinguished, given that certain nonlinear quantities—such as covariances and higher-order moments of conditional expectation values—are unraveling dependent. We show that these quantities cannot be accessed unless the measurement scheme (i.e., the unraveling) is known in advance. This renders any operational protocol to distinguish unravelings fundamentally unfeasible. We further establish that assuming access to such nonlinear quantities without prior knowledge of the unraveling would enable superluminal signaling, violating relativistic causality.
Stochastic unravelings provide a useful way to represent open quantum system dynamics in terms of pure state realizations, and have been widely studied both from a fundamental and from a computational point of view. They were initially formulated for Markovian dynamics described by the Gorini-Kossakowski-Sudarshan-Lindblad master equation. However, due to recent technological and experimental development, most physical relevant dynamics present temporal correlations beyond the Markov approximation. Such correlations cause decay rates to turn temporarily negative, thus requiring the generalization of stochastic unravelings from Markovian to non-Markovian scenarios. Indeed, many unraveling techniques have been introduced in this regime, and a comprehensive review of the different jump methods is currently missing. In this work, we provide an overview of widely used quantum jump unraveling techniques for non-Markovian systems and also discuss them in terms of their numerical efficiency, divisibility requirements, Hilbert space extension, and measurement interpretation.
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.
Unravelings provide a probabilistic representation of solutions of master equations and a method of computation of the density operator dynamics. The trajectories generated by unravelings may also be treated as real -- as in the stochastic collapse models. While averages of linear functionals of the unraveling trajectories can be calculated from the master equation, the situation is different for nonlinear functionals, thanks to the corrections with nonzero expected values, coming from the Itô formula. Two types of nonlinear functionals are considered here: variance, and entropy. The corrections are calculated explicitly for two types of unravelings, based on Poisson and Wiener processes. In the case of entropy, these corrections are shown to be negative, expressing the localization introduced by the Lindblad operators.