We study the classical simulability of open quantum dynamics using random Clifford circuits doped with non-Clifford phase rotations and subject to local noise. We unravel the dynamics into stochastic quantum trajectories simulated with Clifford-augmented matrix product states, and introduce a simulation cost that quantifies the classical resources required. Optimizing this cost over stochastic unravelings, we identify noise-induced classical phases: extended parameter regions in which the dynamics can be fully disentangled by Clifford operations at arbitrary circuit depth. Their existence depends on both the noise model \emph{and} the unraveling. Using a geometric representation of quantum channels, we analytically determine optimal unravelings for a broad class of noise models, with numerical simulations confirming the predicted phase boundaries. We further relate the optimal cost to the unraveling-independent nonstabilizerness of the channel and show that, together with trajectory-resolved entanglement and nonstabilizerness, it classifies distinct dynamical regimes. Finally, we show that these classical phases disappear in the averaged density-matrix description, where no unraveling freedom remains. Our results show that the emergence of classicality in noisy random circuits depends on the measurement scheme adopted to probe it, and paves the way to further studies on the classical simulability of driven-dissipative dynamics.
We study the classical simulability of open quantum dynamics using random Clifford circuits doped with non-Clifford phase rotations and subject to local noise. We unravel the dynamics into stochastic quantum trajectories simulated with Clifford-augmented matrix product states, and introduce a simulation cost that quantifies the classical resources required. Optimizing this cost over stochastic unravelings, we identify noise-induced classical phases: extended parameter regions in which the dynamics can be fully disentangled by Clifford operations at arbitrary circuit depth. Their existence depends on both the noise model \emph{and} the unraveling. Using a geometric representation of quantum channels, we analytically determine optimal unravelings for a broad class of noise models, with numerical simulations confirming the predicted phase boundaries. We further relate the optimal cost to the unraveling-independent nonstabilizerness of the channel and show that, together with trajectory-resolved entanglement and nonstabilizerness, it classifies distinct dynamical regimes. Finally, we show that these classical phases disappear in the averaged density-matrix description, where no unraveling freedom remains. Our results show that the emergence of classicality in noisy random circuits depends on the measurement scheme adopted to probe it, and paves the way to further studies on the classical simulability of driven-dissipative dynamics.
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.
The dynamics of open quantum systems can be simulated by unraveling it into an ensemble of pure state trajectories undergoing nonunitary monitored evolution, which has recently been shown to undergo measurement-induced entanglement phase transition. Here, we show that, for an arbitrary decoherence channel, one can optimize the unraveling scheme to lower the threshold for entanglement phase transition, thereby enabling efficient classical simulation of the open dynamics for a broader range of decoherence rates. Taking noisy random unitary circuits as a paradigmatic example, we analytically derive the optimum unraveling basis that on average minimizes the threshold. Moreover, we present a heuristic algorithm that adaptively optimizes the unraveling basis for given noise channels, also significantly extending the simulatable regime. When applied to noisy Hamiltonian dynamics, the heuristic approach indeed extends the regime of efficient classical simulation based on matrix product states beyond conventional quantum trajectory methods. Finally, we assess the possibility of using a quasi-local unraveling, which involves multiple qubits and time steps, to efficiently simulate open systems with an arbitrarily small but finite decoherence rate.
Magic states are essential for universal quantum computation and are widely viewed as a key source of quantum advantage, yet in realistic devices they are inevitably noisy. In this work, we characterize how noise on injected magic resources changes the classical simulability of quantum circuits and when it induces a transition from classically intractable behavior to efficient classical simulation. We adopt a resource-centric noise model in which only the injected magic components are noisy, while the baseline states, operations, and measurements belong to an efficiently simulable family. Within this setting, we develop an approximate classical sampling algorithm with controlled error and prove explicit noise-dependent conditions under which the algorithm runs in polynomial time. Our framework applies to both qubit circuits with Clifford baselines and fermionic circuits with matchgate baselines, covering representative noise channels such as dephasing and particle loss. We complement the analysis with numerical estimates of the simulation cost, providing concrete thresholds and runtime scaling across practically relevant parameter regimes.
We analyze the modification of entanglement dynamics in the Grover algorithm when the qubits are subjected to single-qubit amplitude-damping or phase-flip noise. We compare quantum trajectories with full density-matrix simulations, analyzing the dynamics of averaged trajectory entanglement (TE) and operator entanglement (OE), in the respective state representation. Although not a genuine entanglement measure, both TE and OE are connected to the efficiency of matrix product state simulations and thus of fundamental interest. As in many quantum algorithms, at the end of the Grover circuit entanglement decreases as the system converges towards the target product state. While we find that this is well captured in the OE dynamics, quantum trajectories rarely follow paths of reducing entanglement. Optimized unraveling schemes can lower TE slightly, however we show that deep in the circuit OE is generally smaller than TE. This implies that matrix product density operator (MPDO) simulations of quantum circuits can in general be more efficient than quantum trajectories. In addition, we investigate the noise-rate scaling of success probabilities for both amplitude-damping and phase-flip noise in Grover's algorithm.
Classical simulations of noisy quantum circuits are instrumental to our understanding of the behavior of real-world quantum systems and the identification of regimes where one expects quantum advantage. In this work, we present a highly parallelizable tensor-network-based classical algorithm, equipped with rigorous accuracy guarantees, for simulating n -qubit quantum circuits with arbitrary single-qubit noise. Our algorithm represents the state of a noisy quantum system by a particular ensemble of matrix product states from which we stochastically sample. Each pure state evolved under a single-qubit noise process is then represented by the ensemble of states that achieves the minimal average entanglement (the entanglement of formation) between the noisy qubit and the remainder. This approach lets us use a more compact representation of the quantum state for a given accuracy requirement and noise level. For a given maximum bond dimension χ and circuit, our algorithm comes with an upper bound on the simulation error, runs in poly ( n , χ ) -time, and improves upon related prior work (1) in scope by extending analytic methods from the three commonly considered noise models to general single-qubit noise, (2) in performance by deriving an exact, closed-form solution to the local entanglement minimization problem, previously approached via variational heuristics, thereby guaranteeing local optimality without numerical optimization overhead, and (3) in conceptual contribution by showing that the fixed unraveling used in prior work becomes equivalent to our choice of unraveling in the special case of depolarizing, dephasing, and amplitude damping noise acting on a maximally entangled state.
We show that states obtained from deep random Clifford circuits doped with non-Clifford phase gates (including T gates and sqrt[T] gates) can be disentangled completely, provided the number of non-Clifford gates is smaller or approximately equal to the number of qubits. This implies that Pauli expectation values of such states can be efficiently simulated classically, despite them exhibiting both extensive entanglement and extensive nonstabilizerness. We prove this result analytically using a quantum error correction formulation, demonstrate its applicability numerically, and discuss consequences for the disentanglability of states generated through Hamiltonian dynamics. We show that this result implies a novel representation of approximate state designs that can also facilitate their efficient generation, and we propose a novel quantum circuit compression scheme for Clifford circuits doped with non-Clifford phase gates.