The ensemble-averaged dynamics of open quantum systems are typically irreversible. We show that this irreversibility need not hold at the level of individually monitored quantum trajectories. Our main results are quantum reverse diffusion stochastic differential equations, along with corresponding stochastic master equations. These equations describe the exact and approximate reverse diffusion dynamics for continuously monitored Pauli channels, including time-dependent depolarizing noise. Our results bridge the gap between highly nonlinear classical reverse diffusion, prominent in generative modelling, and linear quantum mechanics. Consequently, this establishes a theoretical framework for diffusion-driven quantum gates, quantum tomography via forward-reverse cycles, and potential paradigms for quantum error correction based on reverse diffusion.
Paper ID : 2511.15919Title : Reversing Quantum Noise: Quantum Reverse Diffusion for Pauli ChannelsAuthor : Einar Gabbassov (University of Waterloo, Perimeter Institute for Theoretical Physics, Institute for Quantum Computing)Classification : quant-ph, math-ph, math.MPSubmission Date : November 19, 2025Paper Link : https://arxiv.org/abs/2511.15919 The ensemble-averaged dynamics of open quantum systems are typically irreversible. This paper demonstrates that at the level of individual monitored quantum trajectories, such irreversibility is not inevitable. The main contributions include quantum reverse diffusion stochastic differential equations and their corresponding stochastic master equations, which precisely or approximately describe the reverse diffusion dynamics of continuously monitored Pauli channels (including time-dependent depolarizing noise). The results bridge the gap between the highly nonlinear reverse diffusion prominent in classical generative modeling and linear quantum mechanics, establishing a theoretical framework for diffusion-driven quantum gates, forward-reverse cycle-based quantum tomography, and a new paradigm for reverse diffusion-based quantum error correction.
Irreversibility Paradox : The ensemble-averaged dynamics of open quantum systems obey the Lindblad master equation and are irreversible; Pauli channels as discrete noise models are also irreversible. However, at the level of single realizations of quantum trajectories, does there exist a reverse process?Classical-Quantum Gap : Classical reverse stochastic differential equations (e.g., Score-based generative models) are highly nonlinear, with drift terms determined by gradients of probability densities; quantum mechanics is inherently linear, with nonlinearity arising only from state normalization. How can one define a reverse process in the quantum domain?Theoretical Foundation of Quantum Generative Models : Existing quantum generative models mostly employ variational heuristic approaches, training parameterized circuits to approximate reverse processes, but lack a fundamental definition based on physical principles.Theoretical Significance : Reveals reversibility under quantum measurement and feedback mechanisms, challenging conventional irreversibility paradigmsApplication Value : Provides new theoretical tools for quantum error correction, quantum state preparation, and quantum tomographyInterdisciplinary Bridge : Connects classical diffusion models with quantum information, providing first-principles foundations for quantum machine learningVariational Quantum Algorithms : Lack theoretical guarantees, depend on training data, cannot provide analytical formsPetz Recovery Map : Primarily for discrete time, with unclear relationship to continuous monitoring processesImaginary Time Evolution Algorithms : Require local tomography and extensive offline processing, unsuitable for online feedback scenariosAnalytical Reverse SDE : First derivation of the exact quantum reverse stochastic differential equation for single Pauli error channels (Eq. 2), proving complete recovery of the initial state under ideal measurement efficiency and instantaneous feedbackReverse Master Equation : Establishes corresponding stochastic master equations (Eqs. 3-4), providing density matrix level descriptionsDepolarizing Noise Approximation : Constructs approximate reverse SDEs for multi-Pauli error channels (depolarizing noise) (Eq. 7), with fidelity error O((pT)³) in the pT < 1 regionPhysical Implementation Mechanisms :Information-conserving case: Implementation via coherent Hamiltonian drive Information-dissipative case: Proposes algorithm based on weak measurement and state teleportation (Appendix B) Application Framework :Diffusion-driven quantum gates (Eq. 10) Forward-reverse cycle quantum tomography Correction schemes for measurement-induced errors Input : Quantum state |ψ(T)⟩ after forward Pauli noise process over time T and measurement record W(T)Output : Initial state |ϕ(2T)⟩ ≈ |ψ₀⟩ recovered at time 2TConstraint : Forward noise persists during the reverse process
Diffusive unravelling is described as:
d ∣ ψ ⟩ = ( − p 2 I d t + p L d W ) ∣ ψ ( t ) ⟩ , ∣ ψ ( 0 ) ⟩ = ∣ ψ 0 ⟩ d|\psi\rangle = \left(-\frac{p}{2}I dt + \sqrt{p}L dW\right)|\psi(t)\rangle, \quad |\psi(0)\rangle = |\psi_0\rangle d ∣ ψ ⟩ = ( − 2 p I d t + p L d W ) ∣ ψ ( t )⟩ , ∣ ψ ( 0 )⟩ = ∣ ψ 0 ⟩
where:
Jump Operators : L = P (information-dissipative) or L = iP (information-conserving), P is a Pauli operatorNoise Strength : p ∈ 0,1 Measurement Increment : dW = √p⟨L+L†⟩dt + dŴ (information-dissipative) or dW = dŴ (information-conserving)Solution of the forward process:
∣ ψ ( t ) ⟩ = F ( t ) ∣ ψ 0 ⟩ = exp ( − p t + p L W ( t ) ) ∣ ψ 0 ⟩ |\psi(t)\rangle = F(t)|\psi_0\rangle = \exp(-pt + \sqrt{p}LW(t))|\psi_0\rangle ∣ ψ ( t )⟩ = F ( t ) ∣ ψ 0 ⟩ = exp ( − pt + p L W ( t )) ∣ ψ 0 ⟩
R ( t ) = { I , 0 ≤ t < T exp ( ( X ( t ) − W ( t ) ) p L ) , T ≤ t ≤ 2 T R(t) = \begin{cases}
I, & 0 \leq t < T \\
\exp\left((X(t)-W(t))\sqrt{p}L\right), & T \leq t \leq 2T
\end{cases} R ( t ) = { I , exp ( ( X ( t ) − W ( t )) p L ) , 0 ≤ t < T T ≤ t ≤ 2 T
where X(t) is a Brownian bridge process:
d X ( t ) = − X ( t ) 2 T − t d t + d W , X ( T ) = W ( T ) , X ( 2 T ) = 0 dX(t) = -\frac{X(t)}{2T-t}dt + dW, \quad X(T) = W(T), \quad X(2T) = 0 d X ( t ) = − 2 T − t X ( t ) d t + d W , X ( T ) = W ( T ) , X ( 2 T ) = 0
d ∣ ϕ ( t ) ⟩ = [ ( − p 2 I − X ( t ) 2 T − t p L ) d t + p L d W ] ∣ ϕ ( t ) ⟩ d|\phi(t)\rangle = \left[\left(-\frac{p}{2}I - \frac{X(t)}{2T-t}\sqrt{p}L\right)dt + \sqrt{p}LdW\right]|\phi(t)\rangle d ∣ ϕ ( t )⟩ = [ ( − 2 p I − 2 T − t X ( t ) p L ) d t + p L d W ] ∣ ϕ ( t )⟩
Key Characteristics :
Stochastic Drift Term : − X ( t ) 2 T − t p L -\frac{X(t)}{2T-t}\sqrt{p}L − 2 T − t X ( t ) p L drives the state back to initial configurationNoise Term Preservation : p L d W \sqrt{p}LdW p L d W remains the same as in the forward processIntegrable Singularity : The singularity at t=2T is integrable due to X(2T)=0∣ ϕ ( 2 T ) ⟩ = R ( 2 T ) F ( 2 T ) ∣ ψ 0 ⟩ = exp ( 0 ) ∣ ψ 0 ⟩ = ∣ ψ 0 ⟩ |\phi(2T)\rangle = R(2T)F(2T)|\psi_0\rangle = \exp(0)|\psi_0\rangle = |\psi_0\rangle ∣ ϕ ( 2 T )⟩ = R ( 2 T ) F ( 2 T ) ∣ ψ 0 ⟩ = exp ( 0 ) ∣ ψ 0 ⟩ = ∣ ψ 0 ⟩
d ∣ ψ ( t ) ⟩ = ( − 1 2 ∑ k = 1 3 p 3 L k † L k d t + ∑ k = 1 3 p 3 L k d W k ) ∣ ψ ( t ) ⟩ d|\psi(t)\rangle = \left(-\frac{1}{2}\sum_{k=1}^3 \frac{p}{3}L_k^\dagger L_k dt + \sum_{k=1}^3 \sqrt{\frac{p}{3}}L_k dW_k\right)|\psi(t)\rangle d ∣ ψ ( t )⟩ = ( − 2 1 ∑ k = 1 3 3 p L k † L k d t + ∑ k = 1 3 3 p L k d W k ) ∣ ψ ( t )⟩
Three independent Pauli noise channels, Lₖ = σₖ (dissipative) or iσₖ (conserving)
d ∣ ϕ ( t ) ⟩ = ( D ( t ) d t + ∑ k = 1 3 H k ( t ) d X k ) ∣ ϕ ( t ) ⟩ d|\phi(t)\rangle = \left(D(t)dt + \sum_{k=1}^3 H_k(t)dX_k\right)|\phi(t)\rangle d ∣ ϕ ( t )⟩ = ( D ( t ) d t + ∑ k = 1 3 H k ( t ) d X k ) ∣ ϕ ( t )⟩
where:
D ( t ) : = − p I + γ 2 2 ( 3 − 2 ∑ k = 1 3 ( X k ( t ) − X k ( T ) ) 2 ) I D(t) := -pI + \frac{\gamma^2}{2}\left(3 - 2\sum_{k=1}^3(X_k(t)-X_k(T))^2\right)I D ( t ) := − p I + 2 γ 2 ( 3 − 2 ∑ k = 1 3 ( X k ( t ) − X k ( T ) ) 2 ) I
H k ( t ) : = σ k + 1 2 ∑ j = 1 3 [ σ j , σ k ] ( X j ( t ) − X j ( T ) ) H_k(t) := \sigma_k + \frac{1}{2}\sum_{j=1}^3[\sigma_j, \sigma_k](X_j(t)-X_j(T)) H k ( t ) := σ k + 2 1 ∑ j = 1 3 [ σ j , σ k ] ( X j ( t ) − X j ( T ))
γ : = p 3 + 2 i p 3 \gamma := \sqrt{\frac{p}{3}} + 2i\frac{p}{3} γ := 3 p + 2 i 3 p
Initial Conditions :
X k ( T ) = p 3 W k ( T ) + 2 i p 3 S k ( T ) X_k(T) = \sqrt{\frac{p}{3}}W_k(T) + 2i\frac{p}{3}S_k(T) X k ( T ) = 3 p W k ( T ) + 2 i 3 p S k ( T )
where Sₖ(T) is the Lévy stochastic area:
S 23 ( T ) = 1 2 ∫ 0 T ( W 2 d W 3 − W 3 d W 2 ) S_{23}(T) = \frac{1}{2}\int_0^T (W_2 dW_3 - W_3 dW_2) S 23 ( T ) = 2 1 ∫ 0 T ( W 2 d W 3 − W 3 d W 2 )
E [ F ( ∣ ϕ ^ ( 2 T ) ⟩ , ∣ ψ 0 ⟩ ) ] ≥ 1 − c ( p T ) 3 \mathbb{E}[F(|\hat{\phi}(2T)\rangle, |\psi_0\rangle)] \geq 1 - c(pT)^3 E [ F ( ∣ ϕ ^ ( 2 T )⟩ , ∣ ψ 0 ⟩)] ≥ 1 − c ( pT ) 3
In the pT < 1 region, the error is a third-order small quantity.
Clever Use of Brownian Bridge : The boundary condition X(2T)=0 automatically eliminates the singularity without requiring explicit state knowledgeMagnus Expansion Truncation :Second-order Magnus expansion approximates forward operator (Eq. A2) Error control: O ( p 3 / 2 T 3 / 2 ) \mathcal{O}(p^{3/2}T^{3/2}) O ( p 3/2 T 3/2 ) Derivative truncation (Eq. A11): Maintains analytical tractability Information Classification Treatment :Information-Conserving (L=iP): Drift corresponds to Hamiltonian H ( t ) = X ( t ) p P 2 T − t H(t) = \frac{X(t)\sqrt{p}P}{2T-t} H ( t ) = 2 T − t X ( t ) p P Information-Dissipative (L=P): Drift corresponds to imaginary time evolutionPhysical Meaning of Lévy Area : Captures pairwise interactions between different noise processes, direct manifestation of non-commutativityThis work is primarily theoretical without extensive numerical experiments, but provides the following verification directions:
Single Qubit Pauli Channels :Exact analytical solutions (Eqs. 11-14) Perfect recovery verifiable through numerical simulation Depolarizing Noise :Parameter region: pT < 1 Fidelity lower bound: 1 - c(pT)³ Direct Implementation : Via time-dependent HamiltonianControl Parameters : Construct X(t) in real-time from measurement record dWHardware Requirements : Coherent control + continuous measurementResource State Preparation :
∣ Φ rsrc ( r ) ⟩ = 1 P r ( 0 A ) ( I ⊗ R k 2 r ( t ) ) ∣ Φ ⟩ |\Phi^{(r)}_{\text{rsrc}}\rangle = \frac{1}{\sqrt{Pr(0_A)}}(I \otimes R_k^{2^r}(t))|\Phi\rangle ∣ Φ rsrc ( r ) ⟩ = P r ( 0 A ) 1 ( I ⊗ R k 2 r ( t )) ∣Φ ⟩ Success Probability :
P r ( fail ) worst = 1 2 ( 1 + 2 p ∣ Δ Y ∣ ) + O ( p Δ Y 2 ) Pr(\text{fail})_{\text{worst}} = \frac{1}{2}(1 + 2\sqrt{p}|\Delta Y|) + \mathcal{O}(p\Delta Y^2) P r ( fail ) worst = 2 1 ( 1 + 2 p ∣Δ Y ∣ ) + O ( p Δ Y 2 ) Resource Budget :
d min = ⌈ log 2 ( 1 / ϵ ) 1 − 2 p ∣ Δ Y ( t ) ∣ ⌉ d_{\min} = \left\lceil \frac{\log_2(1/\epsilon)}{1 - 2\sqrt{p}|\Delta Y(t)|} \right\rceil d m i n = ⌈ 1 − 2 p ∣Δ Y ( t ) ∣ l o g 2 ( 1/ ϵ ) ⌉ Exact Recovery (Single Pauli Channel) :
∣ ϕ ^ ( 2 T ) ⟩ = ∣ ψ 0 ⟩ (exact) |\hat{\phi}(2T)\rangle = |\psi_0\rangle \quad \text{(exact)} ∣ ϕ ^ ( 2 T )⟩ = ∣ ψ 0 ⟩ (exact) Approximate Recovery (Depolarizing Noise) :
E [ ∥ ∣ ϕ ( 2 T ) ⟩ − ∣ ψ 0 ⟩ ∥ 2 ] = O ( p 3 T 3 ) \mathbb{E}[\||\phi(2T)\rangle - |\psi_0\rangle\|^2] = \mathcal{O}(p^3T^3) E [ ∥∣ ϕ ( 2 T )⟩ − ∣ ψ 0 ⟩ ∥ 2 ] = O ( p 3 T 3 ) Statistical Time Reversal : The distribution of the reverse process is consistent with the time reversal of the forward processBlind Operability : The reverse SDE contains no information about the initial state, effective for any state undergoing the forward processNoise Coexistence : The reverse process can recover the state even with noise persistingSingularity Resolution : The boundary condition of the Brownian bridge naturally handles the singularity at t=2TNon-Commutation Effects : The Lévy area term quantifies interactions between different Pauli errorsObjective : Implement rotation gate G(θ) = cos(θ)I - i sin(θ)P
Traditional Method : Coherent evolution d|ψ⟩/dt = -iH|ψ⟩Noisy Environment : d|ψ⟩ = (-iHdt - p/2 Idt + i√p P dW)|ψ⟩ (cannot be deterministically implemented)
Diffusion-Driven Scheme :
d ∣ ψ ( t ) ⟩ = ( D ′ d t + i p P d W ) ∣ ψ ( t ) ⟩ d|\psi(t)\rangle = (D'dt + i\sqrt{p}PdW)|\psi(t)\rangle d ∣ ψ ( t )⟩ = ( D ′ d t + i p P d W ) ∣ ψ ( t )⟩ D ′ : = − p 2 I − i θ / p + X ( t ) 2 T − t p P D' := -\frac{p}{2}I - i\frac{\theta/\sqrt{p} + X(t)}{2T-t}\sqrt{p}P D ′ := − 2 p I − i 2 T − t θ / p + X ( t ) p P
Result : |ψ(2T)⟩ = G(θ)|ψ₀⟩ (deterministic)
Advantages :
Exploits noise rather than fighting it Generalizable to random θ for generating state manifolds Protocol :
Forward Phase (0→T): Weak measurement of P, obtain record dW = 2√p⟨P⟩dt + dŴReverse Phase (T→2T): Continue measurement while implementing reverse driftData Processing : Quantum filtering (Belavkin-SME) + Bayesian estimationAdvantages :
Two-phase measurement on single copy Balances signal strength and perturbation Allows iteration over multiple forward-reverse cycles Scenario : Measurement-induced Pauli errors
Traditional Method : Immediate error correctionReverse Diffusion Method :
Allow errors to accumulate over time T Implement reverse process (T→2T) Simultaneously eliminate current and accumulated errors Special Advantage : Information-conserving case (L+L†=0) implementable via coherent drive
Anderson (1982) : First reverse SDE theorySong et al. (2020) : Score-based generative models, drift term: ∇ₓ log p(x,t)Characteristics : Highly nonlinear, depends on probability densityZhang et al. (2024) : Variational quantum diffusion modelsKwun et al. (2025) : Quantum denoisingParigi et al. (2024) : Quantum generative adversarial networksLimitations : Lack first-principles foundation, require trainingWiseman (1994) : Quantum feedback controlGross et al. (2018) : Measurement-driven quantum gatesAlbarelli & Genoni (2024) : Measurement optimizationPetz Recovery Map : Discrete-time quantum information recoveryKwon et al. (2022) : Continuous-time Petz mapQuantum Bayesian Theory : Parzygnat & Fullwood (2023)Paper Positioning : First analytical reverse theory for continuously monitored Pauli channels, connecting classical diffusion with quantum measurement
Hierarchy of Reversibility : While ensemble-averaged dynamics are irreversible, individual quantum trajectories are reversibleAnalytical Forms : Single Pauli channels exactly reversible, depolarizing noise approximately reversible (error O((pT)³))Physical Mechanisms :Information-Conserving: Coherent Hamiltonian Information-Dissipative: Imaginary time evolution (requires special implementation) Application Prospects : New paradigms for quantum gates, tomography, and error correctionExactness Constraints :Depolarizing noise only approximate (second-order Magnus truncation) Higher-order terms involve infinite nested commutators Implementation Challenges :Information-dissipative case requires non-deterministic teleportation Resource state budget grows exponentially with failure probability Robust implementation of online imaginary time evolution lacking Applicability Range :Requires pT < 1 (weak noise or short time) Assumes ideal measurement efficiency Requires instantaneous feedback Theoretical Gaps :Precise relationship with Petz recovery map not established Connection to quantum Bayes theorem to be explored Robust Online Implementation : Deterministic schemes for information-dissipative reverse dynamicsGeneralization to Other Noise :Non-Markovian noise Coherent errors Leakage errors Many-Body Systems : Reverse diffusion of entangled statesQuantum Generative Models : First-principles-based quantum diffusion modelsTheoretical Deepening :Relationship with continuous-time Petz map Minimum action principle for quantum Bayesian inversion Theoretical Breakthrough :First analytical form of quantum reverse SDE Challenges conventional wisdom that "Pauli channels are irreversible" Achieves "reversal" within linear quantum mechanics framework Mathematical Rigor :Complete derivations (methods + appendices) Error control (Magnus expansion, derivative truncation) Elegant use of Brownian bridge Physical Insights :Distinguishes information-conserving/dissipative cases Identifies physical meaning of Lévy area Blind operability (no initial state knowledge required) Application Breadth : Unified framework for quantum gates, tomography, error correctionInterdisciplinary Bridge : Connects classical diffusion models with quantum informationAbsence of Experimental Verification :No numerical simulation results No comparison with variational methods Lacks parameter estimation for real systems Immature Implementation Schemes :Appendix B algorithm requires extensive resource states Teleportation success rate decreases with time steps Lacks integration with existing quantum hardware Limited Applicability :Strong constraint pT < 1 Assumes ideal measurement efficiency Single Pauli noise type Theoretical Incompleteness :Depolarizing noise only second-order approximation Existence of exact solutions not proven Relationship with other recovery theories unclear Presentation Issues :Key theorems not numbered Figure 1 lacks clarity (trajectory illustration) Missing algorithm pseudocode Short-term Impact (1-2 years):
Stimulates discussion on quantum trajectory reversibility Inspires new measurement-based quantum control schemes Provides theoretical foundation for quantum diffusion models Medium-term Impact (3-5 years):
May catalyze new quantum error correction protocols Advances quantum tomography methods Influences NISQ-era quantum algorithm design Long-term Impact (5+ years):
If online dissipative implementation succeeds, may transform quantum control paradigm Connects quantum information recovery with machine learning foundations Provides new perspective on time reversal in quantum thermodynamics Practical Value :
High (Theory): Provides first-principles foundation for quantum reverse processesMedium (Application): Requires overcoming implementation barriersHigh (Inspiration): Opens new research directionsReproducibility :
Theoretical derivations verifiable Numerical simulations relatively straightforward Experimental implementation challenging Ideal Scenarios :
Weak noise environment (pT ≪ 1) High-efficiency continuous measurement systems (superconducting qubits, trapped ions) Fast feedback capability (FPGA control) Information-conserving case (coherent drive easily implemented) Unsuitable Scenarios :
Strong noise or long-time evolution Low measurement efficiency Significant feedback delay Complex many-body systems (resource overhead) Potential Application Domains :
Quantum Sensing: Forward-reverse cycles enhance measurement Quantum Communication: Partial recovery of noisy channels Quantum Computing: Diffusion-driven state preparation Quantum Simulation: Reverse evolution probes Hamiltonians X ( t ) ∼ N ( 2 T − t T W ( T ) , t − T T ( 2 T − t ) ) X(t) \sim \mathcal{N}\left(\frac{2T-t}{T}W(T), \frac{t-T}{T}(2T-t)\right) X ( t ) ∼ N ( T 2 T − t W ( T ) , T t − T ( 2 T − t ) )
Key Properties :
X(2T) = 0 (almost surely) Integral ∫ T 2 T X ( t ) 2 T − t d t \int_T^{2T} \frac{X(t)}{2T-t}dt ∫ T 2 T 2 T − t X ( t ) d t finite (a.s.) Truncation Error : Neglecting third and higher-order termsDerivative Approximation : e ad χ − 1 ad χ ≈ I + 1 2 ad χ \frac{e^{\text{ad}_\chi}-1}{\text{ad}_\chi} \approx I + \frac{1}{2}\text{ad}_\chi ad χ e ad χ − 1 ≈ I + 2 1 ad χ Cumulative Effect : Integration over T, 2T intervalResource State Hierarchy : ∣ Φ rsrc ( r ) ⟩ |\Phi^{(r)}_{\text{rsrc}}\rangle ∣ Φ rsrc ( r ) ⟩ implements R k 2 r R_k^{2^r} R k 2 r Failure Recovery : Correction via higher-order resource statesFinal Convergence : Uses identity R k n σ 3 R k n σ 3 = λ min n I R_k^n \sigma_3 R_k^n \sigma_3 = \lambda_{\min}^n I R k n σ 3 R k n σ 3 = λ m i n n I 1 B. D. Anderson (1982): Classical reverse SDE theory3 Y. Song et al. (2020): Score-based generative models21 A. Barchielli & M. Gregoratti (2009): Quantum trajectory theory27 L. Bouten et al. (2007): Quantum filtering50 H. Kwon et al. (2022): Continuous-time Petz map54 A. J. Parzygnat & J. Fullwood (2023): Quantum Bayes theorem
Overall Assessment : This is a groundbreaking theoretical work that provides the first analytical form of quantum reverse SDEs for Pauli channels at the quantum trajectory level. While implementation schemes remain immature, it offers a fundamentally new perspective on quantum information processing. The paper demonstrates strong mathematical rigor and physical insight, likely becoming an important reference in quantum measurement and feedback control. Recommended future work should focus on solving the online dissipative implementation problem and verifying theoretical predictions through numerical and experimental validation.