搜索结果：Unraveling

Stochastic unravelings of Lindblad-type master equations, such as stochastic Schrödinger equations (SSEs), 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.

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 $4^n$ 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\times10^{-16}$; $1/\sqrt{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\times$ fewer trajectories than Qiskit-Aer's \texttt{batched\_shots\_gpu} at $n=10$, a factor that holds at $19$--$26\times$ across $n=8$--$20$. A regime map places analog sampling optimal at weak noise and projector at strong noise, crossing near $γt\approx0.35$. We further report a systems finding: Qiskit-Aer applies noise at the

Stochastic unravelings allow to efficiently simulate open system dynamics, yet their application has traditionally been restricted to master equations that preserve both Hermiticity and trace. In this work, we introduce a general framework that extends piecewise-deterministic unravelings to arbitrary trace-nonpreserving master equations, requiring only positivity and Hermiticity of the dynamics. Our approach includes, as special cases, unravelings of arbitrary dynamics in the Heisenberg picture, evolutions interpolating between fully Lindblad and non-Hermitian Hamiltonian generators, and equations employed in the derivation of full counting statistics, for which we show it can be used to obtain the moments of the associated probability distribution. The framework is suitable for both trace-decreasing and trace-increasing processes through stochastic disappearance and replication of the stochastic realizations, and it is compatible with different unraveling schemes and with reverse jumps in the non-Markovian regime. Thereby, our approach provides a powerful and versatile simulation method that significantly broadens the applicability of stochastic techniques for open system dynamics

Stochastic unraveling schemes are powerful computational tools for simulating Lindblad equations, offering significant reductions in memory requirements. However, this advantage is accompanied by increased stochastic uncertainty, and the question of optimal unraveling remains open. In this work, we investigate unraveling schemes driven by Brownian motion or Poisson processes and present a comprehensive parametric characterization of these approaches. For the case of a single Lindblad operator and one noise term, this parametric family provides a complete description for unraveling scheme with pathwise norm-preservation. We further analytically derive dynamically optimal quantum state diffusion (DO-QSD) and dynamically optimal quantum jump process (DO-QJP) that minimize the growth rate of the variance of an observable locally in time. Compared to jump process ansatz, DO-QSD offers two notable advantages: firstly, the variance for DO-QSD can be rigorously shown not to exceed that of any jump-process ansatz locally in time; secondly, it has very simple expressions. Numerical results demonstrate that the proposed DO-QSD scheme may achieve substantial reductions in the variance of obser

The dynamics of many-body quantum states in open systems is commonly numerically simulated by unraveling the density matrix into pure-state trajectories. In this work, we introduce a new unraveling strategy that can adaptively minimize the averaged entanglement in the trajectory states. This enables a more efficient classical representation of trajectories using matrix product decompositions. Our new approach is denoted non-unitarity maximizing unraveling (NUMU). It relies on the idea that adaptively maximizing the averaged non-unitarity of a set of Kraus operators leads to a more efficient trajectory entanglement destruction. Compared to other adaptive entanglement lowering algorithms, NUMU is computationally inexpensive. We demonstrate its utility in large-scale simulations with random quantum circuits. NUMU lowers runtimes in practical calculations, and it also provides new insight on the question of classical simulability of quantum dynamics. We show that for the quantum circuits considered here, unraveling methods are much less efficient than full matrix product density operator simulations, hinting to a still large potential for finding more advanced adaptive unraveling schem

In standard treatments of open quantum systems, the reduced dynamics is described starting from the assumption that the system and the environment are initially uncorrelated. This assumption, however, is not always guaranteed in realistic scenarios and several theoretical approaches to characterize initially correlated dynamics have been introduced. For the uncorrelated scenario, stochastic unravelings are a powerful tool to simulate the dynamics, but so far they have not been used in the most general case in which correlations are initially present. In our work, we employ the bath positive (B+) or one-sided positive decomposition (OPD) formalism as a starting point to generalize stochastic unraveling in the presence of initial correlations. Noticeably, our approach doesn't depend on the particular unraveling technique, but holds for both piecewise deterministic and diffusive unravelings. This generalization allows not only for more powerful simulations for the reduced dynamics, but also for a deeper theoretical understanding of open system dynamics.

Hagfish slime is a unique biological material composed of mucus and protein threads that rapidly deploy into a cohesive network when deployed in seawater. The forces involved in thread deployment and interactions among mucus and threads are key to understanding how hagfish slime rapidly assembles into a cohesive, functional network. Despite extensive interest in its biophysical properties, the mechanical forces governing thread deployment and interaction remain poorly quantified. Here, we present the first direct in situ measurements of the micromechanical forces involved in hagfish slime formation, including mucus mechanical properties, skein peeling force, thread-mucus adhesion, and thread-thread cohesion. Using a custom glass-rod force sensing system, we show that thread deployment initiates when peeling forces exceed a threshold of approximately 6.8 nN. To understand the flow strength required for unraveling, we used a rheo-optic setup to impose controlled shear flow, enabling us to directly observe unraveling dynamics and determine the critical shear rate for unraveling of the skeins, which we then interpreted using an updated peeling-based force balance model. Our results rev

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.

Stochastic unravelings are a widely used tool to solve open quantum system dynamics, in which the exact solution is obtained via an average over a stochastic process on the set of pure quantum states. Recently, the generalized rate operator unraveling formalism was derived, allowing not only for an engineering of the stochastic realizations, but also to unravel without reverse jumps even for some dynamics in which P-divisibility is violated, thus hugely improving the simulation efficiency. This is possible because the unraveling depend on an arbitrary non-linear transformation which can incorporate the memory effects. In this work, a stochastic Schrödinger equation for this formalism is derived, both for cases with and without reverse jumps. It is also shown that a failure of this method can be used to witness master equations leading unphysical time evolutions, independently on the particular non-linear transformation considered.

The dynamics of open quantum systems can be simulated by unraveling it into an ensemble of pure state trajectories undergoing non-unitary 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 syste

We show that idealized Dicke superradiant decay from a fully inverted state can at all times be described by a positive statistical mixture of coherent spin states (CSS). Since CSS are separable, this implies that no entanglement is involved in Dicke decay. Based on this result, we introduce a new numerical quantum trajectory approach leading to low-entanglement unravelings. This opens up new possibilities for employing matrix product state (MPS) techniques for large-scale numerical simulations with collective decay processes.

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.

Formal safety guarantees on the synthesis of controllers for stochastic systems can be obtained using correct-by-design approaches. These approaches often use abstractions as finite-state Markov Decision Processes. As the state space of these MDPs grows, the curse of dimensionality makes the computational and memory cost of the probabilistic guarantees, quantified with dynamic programming, scale exponentially. In this work, we leverage decoupled dynamics and unravel, via dynamic programming operations, a tree structure in the Canonical Polyadic Decomposition (CPD) of the value functions. For discrete-time stochastic systems with syntactically co-safe linear temporal logic (scLTL) specifications, we provide provable probabilistic safety guarantees and significantly alleviate the computational burden. We provide an initial validation of the theoretical results on several typical case studies and showcase that the uncovered tree structure enables efficient reductions in the computational burden.

Abstract Simulating mixed-state evolution in open quantum systems is crucial for various chemical physics, quantum optics, and computer science applications. These simulations typically follow the Lindblad master equation dynamics. An alternative approach known as quantum state diffusion unraveling is based on the trajectories of pure states generated by random wave functions, which evolve according to a nonlinear Itô-Schrödinger equation (ISE). This study introduces weak first- and second-order solvers for the ISE based on directly applying the Itô-Taylor expansion with exact derivatives in the interaction picture. We tested the method on free and driven Morse oscillators coupled to a thermal environment and found that both orders allowed practical estimation with a few dozen iterations. The variance was relatively small compared to the linear unraveling and did not grow with time. The second-order solver delivers much higher accuracy and stability with bigger time steps than the first-order scheme, with a small additional workload. However, the second-order algorithm has quadratic complexity with the number of Lindblad operators as opposed to the linear complexity of the first-or

It is well known that the state operator of an open quantum system can be generically represented as the solution of a time-local equation -- a quantum master equation. Unraveling in quantum trajectories offers a picture of open system dynamics dual to solving master equations. In the unraveling picture, physical indicators are computed as Monte-Carlo averages over a stochastic process valued in the Hilbert space of the system. This approach is particularly adapted to simulate systems in large Hilbert spaces. We show that the dynamics of an open quantum system generically admits an unraveling in the Hilbert space of the system described by a Markov process generated by ordinary stochastic differential equations for which rigorous concentration estimates are available. The unraveling can be equivalently formulated in terms of norm-preserving state vectors or in terms of linear ``ostensible" processes trace preserving only on average. We illustrate the results in the case of a two level system in a simple boson environment. Next, we derive the state-of-the-art form of the Diosi-Gisin-Strunz Gaussian random ostensible state equation in the context of a model problem. This equation pro

Infinite-layer nickelates stand as a promising frontier in the exploration of unconventional superconductivity. Their synthesis through topotactic oxygen reduction from the parent perovskite phase remains a complex and elusive process. This study delves into the nano-scale effects of the topotactic lattice transformation within LaNiO$_2$ crystals. Leveraging high-resolution scanning transmission electron microscopy and spectroscopy, our investigations uncover a panorama of structural alterations, including grain boundaries and coherent twin boundaries, triggered by reduction-induced transformations. In addition, our analyses unveil the formation of an oxygen-rich disordered transition phase encircling impurities and pervading crystalline domains, and the internal strain is accommodated by grain boundary formation. By unraveling these nano-scale effects, our findings provide insights into the microscopic intricacies of the topotactic reduction process elucidating the transition from the perovskite to the infinite-layer phase within nickelate bulk crystals.

We consider the estimation of parameters encoded in the measurement record of a continuously monitored quantum system in the jump unraveling, corresponding to a single-shot scenario, where information is continuously gathered. Here, it is generally difficult to assess the precision of the estimation procedure via the Fisher Information due to intricate temporal correlations and memory effects. In this paper we provide a full set of solutions to this problem. First, for multi-channel renewal processes we relate the Fisher Information to an underlying Markov chain and derive a easily computable expression for it. For non-renewal processes, we introduce a new algorithm that combines two methods: the monitoring operator method for metrology and the Gillespie algorithm which allows for efficient sampling of a stochastic form of the Fisher Information along individual quantum trajectories. We show that this stochastic Fisher Information satisfies useful properties related to estimation on a single run. Finally, we consider the case where some information is lost in data compression/post-selection and provide tools for computing the Fisher Information in this case. All scenarios are illus

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 $\mathrm{poly}(n,χ)$-time and improves upon related prior work (1) in scope: by extending analytic methods from the three commonly considered noise models