Classical simulation of noisy quantum circuits is essential for validating algorithms, bench-marking hardware, and guiding error-mitigation strategies. Yet it remains challenging due to either the exponential memory footprint of density matrix methods or the high variance in trajectories of standard Kraus-insertion approaches. We introduce a variance aware tensor network framework that unifies the tensor jump method (TJM) with a local time-dependent variational principle (TDVP) and a sparse Pauli-Lindblad model (SPLM) of hardware noise. Gates are realized as short variational evolutions on the matrix product state (MPS) manifold, while noise is applied stochastically per circuit window using Pauli-Lindblad jump sets that render the sampling hazard state-independent and the dissipative contraction trivial after renormalization. Strikingly, our method supports arbitrary Lindblad correlated multi-qubit noise channels consistent with hardware connectivity, including long-range noise operators on non-adjacent subsets, enabling direct simulation of crosstalk and other connectivity-induced errors beyond local noise models. Building on prior work that introduced analog unitary unravelings, we adapt this scheme to a framework based on the Lindblad master equation with exact generator matching under discretized angle laws and a parameter choice. This exactly matches the Lindblad generator, while selecting Gaussian or two-point angle laws. In parallel, we introduce an MPS-projector-jump unraveling that yields state-independent hazards and closed-form variance laws. Both schemes are unbiased and inherit the standard 1/√N Monte Carlo convergence, but with substantially smaller constants due to variance reduction. Empirically, projector sampling markedly reduces the required bond dimension per trajectory across many circuit architectures, while analog sampling excels at low noise. We demonstrate accurate, scalable noisy-circuit simulation on a 25-qubit noisy XY quench and IBM’s 127-qubit kicked Ising benchmark with long-range depolarizing noise, achieving substantially reduced Monte Carlo variance and favorable bond dimension growth compared to standard Kraus insertion baselines.
Classical simulation of noisy quantum circuits is essential for validating algorithms, bench-marking hardware, and guiding error-mitigation strategies. Yet it remains challenging due to either the exponential memory footprint of density matrix methods or the high variance in trajectories of standard Kraus-insertion approaches. We introduce a variance aware tensor network framework that unifies the tensor jump method (TJM) with a local time-dependent variational principle (TDVP) and a sparse Pauli-Lindblad model (SPLM) of hardware noise. Gates are realized as short variational evolutions on the matrix product state (MPS) manifold, while noise is applied stochastically per circuit window using Pauli-Lindblad jump sets that render the sampling hazard state-independent and the dissipative contraction trivial after renormalization. Strikingly, our method supports arbitrary Lindblad correlated multi-qubit noise channels consistent with hardware connectivity, including long-range noise operators on non-adjacent subsets, enabling direct simulation of crosstalk and other connectivity-induced errors beyond local noise models. Building on prior work that introduced analog unitary unravelings, we adapt this scheme to a framework based on the Lindblad master equation with exact generator matching under discretized angle laws and a parameter choice. This exactly matches the Lindblad generator, while selecting Gaussian or two-point angle laws. In parallel, we introduce an MPS-projector-jump unraveling that yields state-independent hazards and closed-form variance laws. Both schemes are unbiased and inherit the standard 1/√N Monte Carlo convergence, but with substantially smaller constants due to variance reduction. Empirically, projector sampling markedly reduces the required bond dimension per trajectory across many circuit architectures, while analog sampling excels at low noise. We demonstrate accurate, scalable noisy-circuit simulation on a 25-qubit noisy XY quench and IBM’s 127-qubit kicked Ising benchmark with long-range depolarizing noise, achieving substantially reduced Monte Carlo variance and favorable bond dimension growth compared to standard Kraus insertion baselines.
Classical simulation of noisy quantum circuits is essential for validating algorithms, benchmarking hardware, and assessing error-mitigation strategies, but remains limited by the exponential cost of density-matrix methods and the high variance of standard trajectory sampling. We introduce a variance-aware tensor network framework that combines the tensor jump method with local TDVP gate evolution on matrix product states and sparse Pauli-Lindblad hardware noise models. Gates are applied as short variational evolutions on the MPS manifold, while noise is sampled per circuit window from Pauli-Lindblad jump sets with state-independent hazards and dissipative contractions that reduce to irrelevant global factors after renormalization. The method supports correlated multi-qubit Lindblad noise consistent with hardware connectivity, including long-range operators on non-adjacent qubits, enabling direct simulation of crosstalk and other connectivity-induced errors beyond local noise models. We develop two unbiased variance-aware unravelings. An analog unitary-mixture unraveling matches the Lindblad generator exactly under symmetric Gaussian or two-point angle laws, while a projector-jump unraveling yields state-independent hazards and closed-form variance laws. Both retain the standard 1/sqrt(N) Monte Carlo convergence but with reduced prefactors. Empirically, projector sampling strongly reduces trajectory variance and bond-dimension growth across many circuit architectures, whereas analog sampling is most effective at weak noise. We demonstrate accurate, scalable noisy-circuit simulation on a 25-qubit noisy XY quench and IBM's 127-qubit kicked-Ising benchmark with long-range depolarizing noise, achieving reduced Monte Carlo variance and favorable MPS bond-dimension growth compared with standard Kraus-insertion baselines.
Understanding interactions between quantum systems and their environments is crucial for developing stable quantum technologies and accurate physical models. Yet, simulating open quantum systems with non-unitary dynamics remains computationally demanding. We introduce the tensor jump method (TJM), a scalable and embarrassingly parallel algorithm for stochastically simulating large-scale open quantum systems governed by Lindbladians. The TJM extends the Monte Carlo wave function (MCWF) approach to matrix product states, employs a dynamic time-dependent variational principle (TDVP) to minimize evolution errors, and introduces a sampling MPS to reduce timestep dependence. This method scales efficiently, ensuring convergence to Lindbladian dynamics independent of system size, as demonstrated both rigorously and numerically. We showcase its utility by simulating XXX Heisenberg models with up to a thousand spins on a consumer-grade CPU. The TJM represents a significant advance in open quantum system simulation, enabling exploration of dissipative many-body dynamics and the design of more stable quantum hardware.
Abstract Classical simulations of quantum circuits are vital for assessing potential quantum advantage and benchmarking devices, yet they require sophisticated methods to avoid the exponential growth of resources. Tensor network approaches, in particular matrix product states (MPS) combined with the time-evolving block decimation (TEBD) algorithm, currently dominate large-scale circuit simulations. These methods scale efficiently when entanglement is limited but suffer rapid bond dimension growth with increasing entanglement and handle long-range gates via costly SWAP insertions. Motivated by the success of the time-dependent variational principle (TDVP) in many-body physics, we reinterpret quantum circuits as series of discrete time evolutions, using gate generators to construct an MPS-based circuit simulation via a local TDVP formulation. This addresses TEBD's key limitations by (1) naturally accommodating long-range gates and (2) optimally representing states on the MPS manifold. By diffusing entanglement more globally, the method suppresses local bond growth and reduces memory and runtime costs. We benchmark the approach on five 49-qubit circuits: three Hamiltonian circuits (1D open and periodic Heisenberg, 2D 7x7 Ising) and two algorithmic ones (quantum approximate optimization, hardware-efficient ansatz). Across all cases, our method yields substantial resource reductions over standard tools, establishing a new state-of-the-art for circuit simulation and enabling advances across quantum computing, condensed matter, and beyond.
Monte-Carlo trajectory (quantum-jump) methods are the practical route to simulating noisy quantum circuits once the exact density-matrix method is precluded by its 4n memory cost. Their bottleneck is estimator variance: resolving one expectation value can demand thousands of trajectories. Recent tensor-network work shows that \emph{variance-reduced unravelings} -- projector and analog sampling -- sharply cut this variance, but only on CPU matrix-product-state backends, with no path into production tooling. We implement both unravelings on a \emph{GPU dense-statevector} trajectory engine and validate them against the exact density matrix (ideal-circuit fidelity 1−2.2×10−16; 1/N convergence; all unravelings unbiased to trace distance <0.01). On a single consumer GPU, projector unraveling reaches a target standard error with 20.8× fewer trajectories than Qiskit-Aer's \texttt{batched\_shots\_gpu} at n=10, a factor that holds at 19--26× across n=8--20. A regime map places analog sampling optimal at weak noise and projector at strong noise, crossing near γt≈0.35. We further report a systems finding: Qiskit-Aer applies noise at the \emph{channel} level and reconstructs a canonical Kraus decomposition at apply time, discarding any user-supplied unraveling, so variance-reduced unravelings cannot be delivered through its public API. Because Aer's Born-rule collapse machinery already exists, we specify a minimal change that would unlock the technique in production.
Understanding interactions between quantum systems and their environments is crucial for developing stable quantum technologies and accurate physical models. Yet, simulating open quantum systems with non-unitary dynamics remains computationally demanding. We introduce the tensor jump method (TJM), a scalable and embarrassingly parallel algorithm for stochastically simulating large-scale open quantum systems governed by Lindbladians. The TJM extends the Monte Carlo wave function (MCWF) approach to matrix product states, employs a dynamic time-dependent variational principle (TDVP) to minimize evolution errors, and introduces a sampling MPS to reduce timestep dependence. This method scales efficiently, ensuring convergence to Lindbladian dynamics independent of system size, as demonstrated both rigorously and numerically. We showcase its utility by simulating XXX Heisenberg models with up to a thousand spins on a consumer-grade CPU. The TJM represents a significant advance in open quantum system simulation, enabling exploration of dissipative many-body dynamics and the design of more stable quantum hardware.
Learning dissipation rates in large-scale open quantum systems is a major obstacle for near-term quantum technologies, as existing Lindblad estimation methods are typically limited to small system sizes due to the computational complexity of repeatedly solving the Lindblad equation during optimization. Here, we propose a scalable noise-learning framework for Lindblad dissipation rates that combines a stochastic simulation method, the Tensor Jump Method (TJM), with gradient-free optimization of a least-squares cost-function defined on time series of local-observable expectation values. We demonstrate the approach on two noise models in the Ising model: a site-resolved (local) model, in which independent dissipation rates are learned for each site up to Nsite=16, and a spatially homogeneous (global) model with only seven parameters, scaled to Nsite=160 sites.We complement these numerical results with a series of exact, provable guarantees: the Frobenius variance of the TJM density-matrix estimator is shown to equal (1−Tr[ρ2])/Ntraj, an exact purity-based characterization of the stochastic estimation error; the corresponding purity evolution is proven to be monotonically non-increasing for Hermitian jump operators; and, under a finite covariance distance assumption, the standard deviation of the cost-function is shown to decrease with system size, so that fewer trajectories are needed to reach a fixed target accuracy as the system grows. Together, this combination of scalable numerics and rigorous theoretical guarantees positions TJM-based noise learning as a practical foundation for characterizing dissipation in large quantum devices and for guiding future work on error mitigation and quantum error correction.