---
title: Deterministic Quantum Jump Method
url: https://www.emergentmind.com/topics/deterministic-quantum-jump-dqj-method
type: topic
---

# Deterministic Quantum Jump Method

Searching arXiv for the cited DQJ-related papers and nearby terminology.
arXiv search: "Deterministic Quantum Jump Method weakly dissipative systems 2603.04066"
arXiv search: "Jump probabilities in the non-Markovian quantum jump method 1001.3566"
The deterministic quantum jump (DQJ) method denotes a class of quantum-trajectory constructions for open quantum systems in which the state evolves deterministically between discrete jump events, while the jump structure is fixed so that the ensemble reproduces the underlying master equation. In the broad trajectory-theory sense, DQJ coincides with a piecewise deterministic process built from a non-Hermitian no-jump evolution and deterministic jump maps. In a narrower recent sense, the term refers to a weak-dissipation algorithm that replaces stochastic sampling of rare jump times by deterministic quadrature over low-jump sectors of the Lindblad evolution. Closely related deterministic formulations also appear in jump-conditioned feedback theory, where one evolves augmented density operators indexed by the last jump record rather than Monte Carlo trajectories [1001.3566; 2603.04066; 2507.01934].

## 1. Definition and scope

In the literature, DQJ is not a single formalism but a family of closely related constructions that retain the jump unravelling of open-system dynamics while minimizing or reorganizing the stochastic component. The common structure is that the state update maps are deterministic, whereas randomness enters either through jump times and channels or through the weighting of jump-number sectors. This places DQJ at the intersection of Lindblad dynamics, quantum trajectories, continuous measurement, and feedback control [2603.04066; 2303.15405].

| Context | Deterministic component | Discrete component |
|---|---|---|
| Time-local master equations | Non-Hermitian drift of pure states | Forward and reverse jumps |
| Weakly dissipative Lindblad dynamics | Deterministic quadrature over jump times | \(0\)-, \(1\)-, \(2\)-jump sectors |
| Continuous-time simulation | Exact waiting-time propagation | Jump time and channel selection |
| Jump-based feedback | Signal-resolved density-operator evolution | Last jump channel and elapsed time |

This variation of usage is substantive rather than terminological. Härkönen’s analysis of time-local non-Markovian master equations gives a general deterministic-plus-jump decomposition of the density operator and derives the jump probabilities used in the non-Markovian quantum jump method [1001.3566]. The 2026 weak-dissipation paper explicitly names the “Deterministic Quantum Jump (DQJ) method” and defines it as a few-jump expansion with deterministic jump-time grids [2603.04066]. The Gillespie formulation for quantum jump trajectories treats continuous-time exact waiting-time sampling as a DQJ-type method at the density-matrix level, especially beyond purity-preserving unravelings [2303.15405]. Deterministic jump equations for feedback go one step further by eliminating sampled trajectories altogether in favor of coupled deterministic equations for signal-resolved states [2507.01934].

A useful unifying statement is that DQJ methods preserve the operational meaning of jumps while transferring as much structure as possible into deterministic evolution, deterministic quadrature, or deterministic enlarged-state dynamics.

## 2. Piecewise-deterministic foundations

A central starting point is the time-local master equation
\[
\frac{ d \rho(t)}{ d t } = -\frac{ i }{ \hbar } [ H(t), \rho(t) ] + \sum_{k} \Delta_k(t) \Big( C_k (t)\rho(t) C_k^\dagger (t) - \frac{1}{2} \big\{ C_k^\dagger (t) C_k (t), \rho (t) \big\} \Big),
\]
with Hermitian \(H(t)\), channel operators \(C_k(t)\), and real time-dependent rates \(\Delta_k(t)\). Härkönen rewrites the density matrix as an ensemble of pure states on projective Hilbert space,
\[
\rho (t) = \int \mathrm{d} \psi\, P[\psi; t] \, |\psi\rangle\langle\psi|,
\]
and decomposes the infinitesimal evolution into two processes: deterministic drift of pure states and drift of the probability density functional \(P[\psi;t]\) [1001.3566].

The deterministic part is generated by the effective non-Hermitian Hamiltonian
\[
H_\mathrm{eff} = H - \frac{i\hbar}{2}\sum_k \Delta_k C_k^\dagger C_k,
\]
followed by renormalization. To first order in \(\delta t\), the state-vector increment is
\[
| \delta \psi  \rangle  = \delta t \Big( - \frac{ i }{ \hbar } H_\mathrm{eff}  + \frac{1}{2} \sum_k \Delta_k \| C_k  \psi \|^2 \Big) |\psi \rangle.
\]
This is the deterministic no-jump evolution used in MCWF and NMQJ, and it is the basic drift component of the DQJ picture [1001.3566].

The complementary stochastic structure is encoded in the conditional jump probabilities. Splitting each rate as \(\Delta_k=\Delta_k^+ - \Delta_k^-\), with \(\Delta_k^\pm \ge 0\), Härkönen obtains
\[
\begin{aligned}
P_{T,S}^k [\phi| \psi ]
=&\ \delta t \,\Delta_k^+ \| C_k \psi \|^2 \,\delta \Big[ | \phi \rangle - \frac{ C_k | \psi \rangle }{\| C_k \psi \|} \Big] \\
&+ \delta t \,\Delta_k^- \frac{ P[\phi ] }{ P[ \psi ] } \| C_k \phi \|^2 \,\delta \Big[ | \psi \rangle - \frac{ C_k | \phi \rangle }{\| C_k \phi \|} \Big].
\end{aligned}
\]
For \(\Delta_k \ge 0\), only forward jumps contribute. For \(\Delta_k<0\), reverse jumps appear and become ensemble dependent through \(P[\phi]/P[\psi]\). This is the formal underpinning of the non-Markovian quantum jump method: standard MCWF is recovered in the Markovian limit, whereas negative-rate intervals require reverse-jump rules that encode memory effects and information backflow [1001.3566].

Two consequences are particularly important for DQJ. First, the dynamics is piecewise deterministic in the precise PDMP sense: deterministic drift between jumps, deterministic state maps at jumps, stochastic timing and channel choice. Second, the unravelling is not unique. Härkönen explicitly notes that a different deterministic propagator \(U\) could be chosen, with compensating changes in the jump process, so the same master equation admits multiple deterministic-jump decompositions [1001.3566].

## 3. Weakly dissipative DQJ as a few-jump expansion

The explicit DQJ method of 2026 is designed for the weakly dissipative regime, defined by
\[
\gamma_j T \ll 1 \qquad \text{for all } j,
\]
so that jumps are rare over the time interval \([0,T]\). Instead of sampling many stochastic trajectories and waiting for rare one- and two-jump events to occur, the method decomposes the density matrix by jump number,
\[
\rho(t) = \sum_{n=0}^\infty p_n\,\rho^{(n)}(t),
\]
with \(p_n\) the probability of exactly \(n\) jumps and \(\rho^{(n)}\) the normalized conditional state for that sector [2603.04066].

The effective Hamiltonian remains the standard
\[
H_\text{eff} = H - \frac{i}{2}\sum_{\hat{L}\in J^1} \hat{L}^\dagger \hat{L},
\]
and the zero-jump sector is simply the normalized non-Hermitian evolution. The single-jump sector is an integral over jump time \(\tau\) and jump operator \(\hat L\),
\[
p_1\rho^{(1)}(t) = \sum_{\hat{L}\in J^1}\int_0^T d\tau\,\tilde{p}_1(\tau,\hat{L})\,\rho^{(1)}_{\tau,\hat{L}}(t),
\]
with
\[
\tilde{p}_1(\tau,\hat{L}) = p^{(0)}_{[0,\tau]}\,\langle \hat{L}^\dagger \hat{L}(\tau)\rangle.
\]
The defining DQJ step is then to replace this integral by a deterministic midpoint quadrature over
\[
\Sigma^1 = \Delta t \cdot\left\{\frac{1}{2}, \frac{3}{2}, \dots, N_\text{grid} - \frac{1}{2}\right\},
\]
and to analogously discretize the two-jump simplex using
\[
\Sigma^2 = \Sigma^2_\text{cartesian}\cup \Sigma^2_\text{bary},
\]
where \(\Sigma^2_\text{bary}\) adds barycentric points to better resolve \(\tau_1 \approx \tau_2\) [2603.04066].

At first order, the reconstructed density matrix is
\[
\rho(t) = p_0\,\rho^{(0)}(t) + \frac{1-p_0}{\mathcal{N}_1}\left(p_1\rho^{(1)}(t)\right),
\]
with
\[
\mathcal{N}_1 = \sum_{\hat{L}\in J^1,\ \tau\in\Sigma^1} p_1(\tau,\hat{L}).
\]
At second order,
\[
\rho(t)\approx p_0\rho^{(0)}(t) +\frac{1-p_0}{\mathcal{N}_1+\mathcal{N}_2}\left(\sum_{n=1}^2 p_n\rho^{(n)}(t)\right).
\]
The method therefore removes Monte Carlo fluctuations and replaces them by controlled quadrature error plus truncation error from neglected higher-jump sectors [2603.04066].

The asymptotic scaling analysis is one of its distinctive features. For maximal jump order \(n\),
\[
\mathcal{E}^{(\mathrm{DQJ})}[\rho_{ij}(T)] \propto \frac{1}{N_\text{traj}^{2/n}},
\qquad
1 - F(\sigma,\rho) \propto \frac{1}{N_\text{traj}^{4/n}},
\]
provided the grid resolves the fastest density-matrix oscillations. Thus the paper reports infidelity scaling \(\sim N_\text{traj}^{-4}\) for \(1\)-jump DQJ and \(\sim N_\text{traj}^{-2}\) for \(2\)-jump DQJ, in contrast with stochastic averaging, which yields \(\sim 1/N_\text{traj}\) in measures like infidelity [2603.04066].

The paper demonstrates this on two benchmark models. For a \(5\)-qubit dissipative transverse-field Ising model with \(\gamma T \approx 0.03\), the \(2\)-jump DQJ method shows an asymptotic infidelity scaling close to \(1/N_\text{traj}^2\), a plateau infidelity of about \(10^{-7}\), and reaches a given target infidelity with several orders of magnitude fewer trajectories than SQJ [2603.04066]. For the dissipative Kerr oscillator with \(\gamma T = 0.32\), single-jump DQJ gives significantly smaller time-domain errors than SQJ, and for the integrated spectral quantity \(|S(\omega_0)|\) it reaches a plateau about \(4.1\times 10^{-5}\) with as few as a few tens of trajectories, whereas SQJ would require roughly \(10^6\) trajectories to cross the same error threshold suggested by a scaling line [2603.04066].

The price of determinism is combinatorial growth with jump order. The number of trajectories is
\[
N_{\text{traj},1\text{ jump}} = 1 + N_\text{grid}N_J
\]
at first order and
\[
N_{\text{traj},2\text{ jumps}} = 1 + N_\text{grid}N_J + \tfrac{1}{2}N_\text{grid}(N_\text{grid}+1)N_J^2
\]
at second order, so the method is sharply targeted at low-jump regimes [2603.04066].

## 4. Continuous-time exact simulation and generalized DQJ sampling

A different DQJ line emphasizes exact continuous-time simulation rather than deterministic quadrature. The quantum Gillespie algorithm reformulates the jump unravelling of a Lindblad equation in terms of the waiting-time distribution
\[
W(T,k|\rho) = \operatorname{tr}\!\left\{L_k^\dagger L_k\,e^{-iH_e T}\rho\,e^{iH_e^\dagger T}\right\},
\]
with
\[
W(T|\rho)=\sum_k W(T,k|\rho)=\operatorname{tr}\!\left\{J\,e^{-iH_e T}\rho\,e^{iH_e^\dagger T}\right\},
\qquad
J=\sum_k L_k^\dagger L_k.
\]
Channel selection at fixed waiting time is then
\[
P(k|T,\rho)=\frac{W(T,k|\rho)}{W(T|\rho)}.
\]
This yields exact continuous-time trajectories specified by a sequence of pairs \((T_1,k_1),(T_2,k_2),\dots\), each followed by deterministic no-jump propagation and deterministic state update [2303.15405].

The main technical device is the operator
\[
Q(T)=e^{iH_e^\dagger T}J e^{-iH_e T},
\]
or more generally \(Q(t_j)=e^{\mathcal{L}_0^\dagger t_j}(J)\) for mixed-state dynamics. Precomputing \(Q(T)\) on a grid allows the waiting-time distribution to be evaluated for arbitrary \(\rho\) through \(W(T_j|\rho)=\operatorname{tr}[Q(T_j)\rho]\), after which one samples a waiting time in one shot rather than by small-step Bernoulli trials [2303.15405].

This formulation is especially important because it extends DQJ-type simulation beyond purity-preserving unravelings. For partial monitoring or channel merging, the no-jump generator becomes
\[
\mathcal{L}_0 = \mathcal{L} - \sum_{k\in\mathcal{M}} \mathcal{J}_k,
\qquad
\mathcal{J}_k\rho=L_k\rho L_k^\dagger,
\]
and the waiting-time density for the next monitored jump is
\[
W(T|\rho)=\operatorname{tr}\big\{ J\, e^{\mathcal{L}_0 T}(\rho)\big\},
\qquad
J=\sum_{k\in\mathcal{M}}L_k^\dagger L_k.
\]
Because these trajectories need not preserve purity, the density matrix itself becomes the trajectory variable, which is not accessible to standard pure-state MCWF in the same way [2303.15405].

The paper positions this as advantageous when the Hilbert space is not huge, a very large number of trajectories is required, and monitored jumps are rare compared with fast coherent dynamics. The tradeoff is explicit: the method trades memory for time by precomputing propagators and adjoint-evolved operators on a time grid [2303.15405]. This gives a generalized DQJ picture in which stochasticity is confined to exact event sampling, while the inter-event motion is fully deterministic.

## 5. Deterministic jump equations for feedback and conditional control

A further development removes sampled trajectories altogether when feedback depends on jump records. In the instrument-based framework of Rosal, Potts, and Landi, one introduces a signal-resolved hybrid state
\[
\varrho_n(y) = E\big[\rho_{x_{1:n}}\,\delta_{y,y_n}\big],
\]
which obeys the deterministic recursion
\[
\varrho_{n+1}(y) =
\sum_{x',y'} \delta_{y,f_{n+1}(x',y')}\,\mathcal{M}_{x'}(y')\,\varrho_n(y').
\]
For quantum jumps, when the signal is the pair \(y_n=(k_n,\tau_n)\) consisting of the last jump channel and the time since the last jump, the continuous-time limit yields
\[
\varrho_t(k,0) = \sum_{k'\in\Sigma}\int_0^t d\tau'\, \mathcal{J}_k(k',\tau')\,\varrho_t(k',\tau'),
\]
\[
\varrho_t(k,\tau) = G(k,\tau)\,\varrho_{t-\tau}(k,0),
\qquad
G(k,\tau)=\mathcal{T}\exp\Bigl[\int_0^\tau ds\,\mathcal{L}_0(k,s)\Bigr].
\]
The unconditional state is recovered as
\[
\bar{\rho}_t = \sum_{k\in\Sigma}\int_0^t d\tau\,\varrho_t(k,\tau).
\]
This is a deterministic quantum-jump equation in the literal sense: jump timing information is encoded as a continuous age variable \(\tau\), and channel information is encoded by the discrete label \(k\), with no Monte Carlo sampling required [2507.01934].

When feedback depends only on the last jump channel, the \(\tau\)-resolved equation reduces to
\[
\partial_t \varrho_t(k)
= -i[H(k),\varrho_t(k)] -\frac{1}{2}\sum_{k'\in\Sigma}\{L_{k'}^\dagger(k)L_{k'}(k),\varrho_t(k)\}
+ \sum_{k'\in\Sigma} L_k(k')\,\varrho_t(k')\,L_k^\dagger(k').
\]
For the qubit population-inversion example with thermal absorption and emission channels, the formalism yields an analytic optimal drive duration \(\tau_1^{\mathrm{opt}}\), an inversion threshold \(p=\gamma/\lambda \lesssim 1.145\), and a maximum tolerable delay \(\tau_0^{\mathrm{max}}\) in the low-temperature limit [2507.01934].

The experimental counterpart of conditional determinism appears in the superconducting three-level system of Minev and collaborators. There, the onset time of a jump remains stochastic, but once heralded by a lull in the monitored bright-state record, the jump from \(|G\rangle\) to \(|D\rangle\) follows a continuous, coherent, and deterministic trajectory in the \(|G\rangle\)–\(|D\rangle\) Bloch sphere, and real-time feedback can catch and reverse the jump mid-flight [1803.00545; 1902.10355]. This experimental result does not define a numerical DQJ algorithm, but it strongly supports the conditional-deterministic interpretation that underlies trajectory-based DQJ methods.

## 6. Variants, design freedom, and limitations

The DQJ idea is broader than Lindblad few-jump quadrature. In the spatial decoherence master equation,
\[
\frac{d\rho}{dt} = -\frac{i}{\hbar}[H,\rho] - \frac{D}{\hbar^2}[x,[x,\rho]],
\]
the orthojump unravelling defines a piecewise deterministic process with nonlinear deterministic evolution
\[
\frac{d\Phi}{dt}
 = -\frac{i}{\hbar}H\Phi
   - \frac{D}{\hbar^2}\bigl(x - \bar x[\Phi]\bigr)^2\Phi,
\]
Poisson rate
\[
\lambda = \frac{2D}{\hbar^2}\,\sigma^2,
\]
and jump rule
\[
\Phi \longrightarrow (x - \bar x)\Phi.
\]
The jump makes the post-jump state orthogonal to the pre-jump state, hence “orthojump.” Numerical Monte Carlo in units \(m=\hbar=D=1\) gives asymptotic center-of-mass spreads
\[
\Delta\widetilde x_\infty \approx 1.62 \pm 0.01,
\qquad
\Delta\widetilde p_\infty \approx 1.63 \pm 0.01,
\]
showing localization even though jumps continue indefinitely [1705.05032].

For quadratic Hamiltonians and linear Lindblad operators, quantum jumps can be represented in a Hagedorn-wavepacket basis. The underlying Gaussian parameters evolve deterministically under non-Hermitian propagation, while jumps act algebraically on a coefficient vector in the moving basis. The jump times still come from the usual hazard-rate condition, but the continuous backbone is deterministic and low dimensional. This yields a Gaussian-based DQJ structure without claiming that the jump unravelling itself becomes non-stochastic [2203.11530].

A different kind of freedom is structural rather than algorithmic. The rate-operator framework of Smirne and collaborators shows that the partition of a master equation into deterministic and jump parts is inherently non-unique. By introducing
\[
\mathcal{J}'_t(\rho)=\mathcal{J}_t(\rho)+\frac12\left(\mathbf{C}(t)\rho+\rho\mathbf{C}^\dagger(t)\right),
\]
one obtains alternative effective Hamiltonians
\[
K'(t)=H'(t)-\frac{i}{2}\Gamma'(t),
\]
and alternative jump schemes for the same mixed-state dynamics. Under suitable conditions this permits a fixed basis of post-jump states, and in driven systems it can make the deterministic evolution time independent by transferring explicit driving dependence into the jump part. Dissipativity provides a sufficient condition for a representation with positive rate operator \(\mathbf{R}_\psi\), and hence for a physically meaningful DQJ-style unravelling with a prescribed deterministic generator [2009.11312].

Across these variants, the main limitations recur. The weakly dissipative DQJ method relies on \(\gamma_j T \ll 1\) and becomes expensive at high jump order because the number of trajectories grows polynomially in \(N_\text{grid}^n N_J^n\) [2603.04066]. Continuous-time density-matrix DQJ methods trade sampling noise for precomputation and memory, which can be severe for large Liouvillians [2303.15405]. Non-Markovian reverse-jump probabilities are ensemble dependent and therefore operationally more intricate than Markovian MCWF [1001.3566]. More generally, the deterministic/jump split is not unique, so “the” DQJ method is best understood as a design principle for unravelings rather than a single canonical algorithm [2009.11312].

Within that broader principle, DQJ methods provide a rigorous way to exploit the deterministic structure latent in quantum-jump dynamics: deterministic drift between events, deterministic state-update maps at events, deterministic quadrature of rare-event sectors, or deterministic enlarged-state equations for jump-conditioned control.

Source: https://www.emergentmind.com/topics/deterministic-quantum-jump-dqj-method