---
title: Dynamically Optimal Quantum Jump Process
url: https://www.emergentmind.com/topics/dynamically-optimal-quantum-jump-process-do-qjp
type: topic
---

# Dynamically Optimal Quantum Jump Process

Searching arXiv for papers on Dynamically Optimal Quantum Jump Process and closely related unraveling schemes.
The Dynamically Optimal Quantum Jump Process (DO-QJP) is an adaptive quantum-jump construction for Markovian open quantum systems in which the decomposition parameters of a Lindblad generator are chosen at every step so as to minimize instantaneous jump rates. In the formulation of Lucas and Hornberger, the resulting expansion is an adaptive resummation of the Lindblad jump series that remains completely positive order by order and typically converges within the lowest two to five orders, thereby supporting both analytic approximation and efficient numerical simulation [1312.1167]. Later work uses the same designation for a stochastic unraveling chosen to minimize the short-time growth of the variance of an observable within a parametric family of jump-process schemes [2509.19887]. A plausible implication is that DO-QJP is best understood as a class of state-adaptive optimization principles for quantum-jump representations, rather than as a single fixed algorithm.

## 1. Markovian setting and jump-series representation

The starting point is a time-local, completely positive and trace-preserving generator $\mathcal{L}(t)$ acting on the density operator $\rho_t$,
$$
\partial_t \rho_t = \mathcal{L}(t)\rho_t,
$$
with standard Lindblad form
$$
\mathcal{L}(t)\rho = -\frac{i}{\hbar}[H,\rho] + \sum_k L_k \rho L_k^\dagger - \frac{1}{2}\{L_k^\dagger L_k,\rho\}.
$$
The exact solution may be written formally as $\rho_t=\mathcal{U}(t,0)\rho_0$, where $\mathcal{U}(t,0)=\mathcal{T}\exp\int_0^t \mathcal{L}(t')\,dt'$ [1312.1167].

A Dyson-like jump expansion follows from an arbitrary splitting
$$
\mathcal{L}(t)=\mathcal{L}_0(t)+\mathcal{J}(t),
$$
which yields the integral equation
$$
\rho_t=\mathcal{U}_0(t,0)\rho_0+\int_0^t dt_1\,\mathcal{U}_0(t,t_1)\mathcal{J}(t_1)\rho_{t_1}.
$$
Iteration gives
$$
\rho_t=\sum_{n=0}^{\infty}\rho_t^{(n)},
$$
with
$$
\rho_t^{(0)}=\mathcal{U}_0(t,0)\rho_0,\qquad
\rho_t^{(n)}=\int_0^t dt_n\,\mathcal{U}_0(t,t_n)\mathcal{J}(t_n)\rho_{t_n}^{(n-1)}.
$$
Each term $\rho_t^{(n)}$ contains exactly $n$ insertions of the jump superoperator $\mathcal{J}$, while the factors $\mathcal{U}_0$ describe continuous evolution between jumps [1312.1167].

The nontrivial issue is not the existence of such a series, but the choice of decomposition. Lucas and Hornberger introduce a physically adapted Lindblad decomposition parameterized by complex shifts $\alpha_j$, so that the expansion is manifestly completely positive order by order. This establishes the basic framework in which DO-QJP is defined [1312.1167].

## 2. Adaptive resummation and the dynamically optimal choice

To preserve positivity term by term, the jump operators are shifted as
$$
L_{j,\alpha}\equiv L_j+\alpha_j,
$$
while the Hamiltonian is compensated by
$$
H_\alpha=H-\frac{i\hbar}{2}\sum_j(\alpha_j^*L_j-\alpha_j L_j^\dagger).
$$
The Lindblad form is invariant under this transformation. One then defines
$$
\mathcal{J}_\alpha\rho\equiv \sum_j L_{j,\alpha}\rho L_{j,\alpha}^\dagger,
$$
and
$$
\mathcal{L}_\alpha\rho\equiv -\frac{i}{\hbar}[H_\alpha^{\mathrm{eff}},\rho]_{\llbracket\cdot,\cdot\rrbracket},
$$
with
$$
H_\alpha^{\mathrm{eff}}=H_\alpha-\frac{i\hbar}{2}\sum_j L_{j,\alpha}^\dagger L_{j,\alpha},
$$
and the anti-commutator-reversed bracket $\llbracket A,B\rrbracket = AB-B^\dagger A^\dagger$ [1312.1167].

The optimization criterion is phrased in terms of the weights
$$
w_n(t)\equiv \mathrm{Tr}\,\rho_t^{(n)},
$$
which satisfy $\sum_n w_n(t)=1$. Their dynamics is governed by the positive jump-rate operator
$$
\Gamma_\alpha\equiv \sum_j L_{j,\alpha}^\dagger L_{j,\alpha},
$$
through the cascade equation
$$
\partial_t w_n=\mathrm{Tr}[\Gamma_\alpha \rho_t^{(n-1)}]-\mathrm{Tr}[\Gamma_\alpha \rho_t^{(n)}].
$$
The expansion is convergent iff the low-$n$ weights dominate [1312.1167].

A finer description resolves each $n$-jump contribution into branches labeled by the full jump record $\mathfrak{R}^n=(j_1,t_1;\ldots;j_n,t_n)$. For each branch, the partial rates are
$$
\mathrm{Tr}[\Gamma_{j,\alpha}\rho_t^{(\mathfrak{R}^n)}].
$$
Minimizing these partial rates with respect to $\alpha_j^*$ gives the optimal shift
$$
\tilde{\alpha}_j(t,\mathfrak{R}^n)= -\frac{\mathrm{Tr}[L_j\,\rho_t^{(\mathfrak{R}^n)}]}{\mathrm{Tr}\,\rho_t^{(\mathfrak{R}^n)}},
$$
and the minimal partial rate becomes
$$
\mathrm{Tr}[L_j^\dagger L_j\,\rho_t^{(\mathfrak{R}^n)}]
-\frac{|\mathrm{Tr}[L_j\,\rho_t^{(\mathfrak{R}^n)}]|^2}{\mathrm{Tr}\,\rho_t^{(\mathfrak{R}^n)}}.
$$
Accordingly, the optimal jump operators are
$$
L_{j,\tilde{\alpha}(t,\mathfrak{R}^n)}=
L_j-\frac{\mathrm{Tr}[L_j\,\rho_t^{(\mathfrak{R}^n)}]}{\mathrm{Tr}\,\rho_t^{(\mathfrak{R}^n)}},
$$
updated adaptively after each jump. Equivalently, one calls the resulting expansion the Dynamically Optimal Quantum Jump Process [1312.1167].

## 3. Convergence mechanism and low-order truncation

The central convergence argument is structural. Because the weights obey the cascade relation above, and because $\Gamma_{j,\tilde\alpha}$ has been minimized at each step, all higher-order rates $\mathrm{Tr}[\Gamma_{\tilde\alpha}\rho_t^{(n)}]$ are as small as possible. In many physically relevant problems the system rapidly localizes into a pointer basis of eigenstates of the dominant $L_j$; for such states, $\mathrm{Tr}[L_j\rho]=\langle L_j\rangle$ equals the eigenvalue, and the optimal partial rate vanishes exactly. Even when it does not vanish exactly, it becomes very small after a few jumps. Hence the expansion typically freezes beyond some small $n$ and converges in only two to five orders [1312.1167].

Convergence may be quantified by the fidelity
$$
\mathcal{F}_k =
F\!\left(
\rho_t,\;
\frac{\sum_{n=0}^{k}\rho_t^{(n)}}{\mathrm{Tr}\sum_{n=0}^{k}\rho_t^{(n)}}
\right),
$$
which in numerics reaches $>0.95$ already for $k\sim 3\ldots 5$ [1312.1167].

Two standard illustrations summarize the reported behavior.

| Model | DO-QJP behavior | Un-resummed comparison |
|---|---|---|
| Damped harmonic oscillator, $N_{\mathrm{th}}=0.5$, $\omega/\gamma=2$, $\tau=6/\gamma$ | $\mathcal{F}_3\approx 0.95$ by $n=3$ jumps | $\sim 19$ jumps for the same accuracy |
| Spatial decoherence | convergence by $n\approx 5$ | $n\approx 80$ needed without resummation |

For the damped harmonic oscillator at finite temperature,
$$
\partial_t\rho=-i\omega[a^\dagger a,\rho]
-\gamma(N_{\mathrm{th}}+1)\left(a\rho a^\dagger-\frac{1}{2}\{a^\dagger a,\rho\}\right)
-\gamma N_{\mathrm{th}}\left(a^\dagger \rho a-\frac{1}{2}\{aa^\dagger,\rho\}\right),
$$
with jump operators $L_1=\sqrt{\gamma(N_{\mathrm{th}}+1)}\,a$ and $L_2=\sqrt{\gamma N_{\mathrm{th}}}\,a^\dagger$, the adaptive resummation gives the reported fidelity gain at low order [1312.1167]. For spatial decoherence,
$$
\partial_t\rho = -\frac{i}{\hbar}\left[\frac{p^2}{2m},\rho\right]
-\frac{4\pi\gamma}{\Lambda_{\mathrm{th}}^2}[x,[x,\rho]],
$$
one again sees low-order convergence under DO-QJP [1312.1167].

The significance of these examples is not merely numerical acceleration. The rapid suppression of higher-order terms reflects the underlying localization mechanism into pointer-like sectors, which is also why the same formalism is informative about decoherence structure and the emergence of classicality [1312.1167].

## 4. Numerical realization and relation to trajectory methods

For the damped harmonic oscillator example, the implementation is described explicitly: one draws jump times $\{t_1,\ldots,t_n\}$, for example by Monte-Carlo sampling, propagates $\rho_0$ under $H_\alpha^{\mathrm{eff}}$ from $0$ to $t_1$, then applies $L_{j_1,\alpha}$, and continues analogously; $\rho_t^{(n)}$ is estimated as the average over many such realizations [1312.1167]. This produces a numerical scheme for efficient simulation, while retaining a direct connection to the analytic jump expansion.

The 2013 construction differs in emphasis from conventional pure-state trajectory approaches. It complements quantum-trajectory methods by lumping branches and optimizing mixed-state contributions rather than single pure-state trajectories [1312.1167]. This distinction matters because the optimization target is the convergence of the expansion weights $w_n$, not the variance of a Monte Carlo estimator.

A later formulation places DO-QJP directly inside the theory of stochastic unravelings of Lindblad equations. In that setting one considers piecewise-deterministic Markov processes on the unit sphere,
$$
d|\psi(t)\rangle = a(\psi)\,dt + b(\psi)\,dJ(t),
$$
where $J(t)$ is a scalar Poisson process with state-dependent rate $\lambda(\psi)$, and the average over pure-state trajectories reconstructs the density matrix,
$$
\rho(t)=E[\,|\psi(t)\rangle\langle\psi(t)|\,].
$$
For one Lindblad operator and one noise term, a parametric family of norm-preserving jump-process unravelings is characterized in terms of functions $\theta(\psi)$, $\gamma(\psi)$, and $\alpha(\psi)\ge 0$; the conventional QJP of Dalibard–Castin–Mølmer corresponds to $\alpha=|\langle\psi|L|\psi\rangle|^2$, $\theta=0$, $\gamma=0$ [2509.19887].

This later perspective recasts DO-QJP as an optimization problem over stochastic trajectories themselves. The conceptual continuity with the adaptive resummation approach is the state-dependent modification of jump structure, but the objective function is different.

## 5. Variance-optimal unravelings and generalized rate-operator schemes

In the variance-based formulation, the aim is to minimize the instantaneous growth of the classical variance
$$
V(t)=\mathrm{Var}[\,\langle\psi|O|\psi\rangle\,]
$$
for a fixed Hermitian observable $O$. The analysis separates the observable-dependent part fixed by the Lindblad generator from a remainder $\delta V(a,b,\lambda)$ that depends on the unraveling. Dynamical optimality means minimizing $\delta V$ pointwise in $\psi$ [2509.19887].

For the jump-process ansatz, the resulting DO-QJP is obtained by choosing $\theta^*(\psi)$ through a scalar minimization problem, saturating a user-prescribed maximum jump rate $\Lambda$ via
$$
\alpha^*(\psi)=\Lambda-(\langle L^\dagger L\rangle-|\langle L\rangle|^2),
$$
and then constructing $\eta^*(\psi)$, $b^*(\psi)$, and $a^*(\psi)$ accordingly. This choice makes the first quadratic term in the variance-growth expression vanish exactly and drives the residual jump term as small as allowed by the bound $\Lambda$ [2509.19887].

The same work also derives dynamically optimal quantum state diffusion (DO-QSD) and proves the local bound
$$
\left.\frac{dV}{dt}\right|_{\mathrm{DO\text{-}QSD}}
\le
\left.\frac{dV}{dt}\right|_{\mathrm{any\ QJP}},
$$
so that, locally in time, even the best-tuned jump process cannot beat the variance growth of the optimal diffusion unraveling. Numerical experiments in that paper focus on DO-QSD rather than DO-QJP; DO-QJP was not tested numerically there because it requires extra tuning of $\Lambda$ and more complex formulas for $\theta^*$ [2509.19887].

A different extension is provided by generalized Rate-Operator quantum jumps via realization-dependent transformations. There the master equation is split into a jump part $\mathcal{J}_t$ and a no-jump driving part $\mathcal{D}_t$, and one exploits the freedom to add a realization-dependent counter-term $C_\psi$. The generalized rate operator is
$$
R_\psi \equiv \mathcal{J}_t[P_\psi] + \frac{1}{2}(C_\psi P_\psi + P_\psi C_\psi^\dagger),
$$
with effective Hamiltonian
$$
K_\psi = H - \frac{i}{2}\Gamma - \frac{i}{2}C_\psi.
$$
Its eigenvalues $\lambda_j(\psi)$ give instantaneous jump rates and its eigenvectors the post-jump states. The hard requirement for a positive unraveling is that all $\lambda_j(\psi)\ge 0$ for the relevant realizations [2402.12445].

That framework does not introduce a single compact cost functional $J[C_\cdot]$ or derive Euler–Lagrange equations. Instead, the freedom in $C_\psi$ is used to optimize physically motivated performance measures such as Shannon entropy of occupation probabilities, total number of jumps per trajectory, classical memory required to track the ensemble, and maximal smoothness so that jumps are rare and far apart. In qubit examples, one chooses a one-parameter family $C_\psi(\phi_1)$, solves analytically for positivity conditions, and then selects $\phi_1(\psi)$ in the feasible interval to optimize the chosen metric [2402.12445].

## 6. Interpretation, applications, and scope of the term

Several implications recur across these formulations. In the adaptive resummation picture, because the expansion adapts itself to the actual state at each jump, it automatically zeroes out further jumps once a pointer basis is reached. This renders otherwise large Hilbert-space problems rapidly tractable to analytic approximation in the lowest orders or to efficient numerical implementation [1312.1167]. The same machinery yields insight into pointer-state structure, minimal entropy production, and the emergence of classicality, and potential applications include steady-state engineering by incoherent control, large-scale open-system simulation, quantum metrology under continuous monitoring, and analytic modeling of decoherence in mesoscopic systems [1312.1167].

In generalized rate-operator approaches, the state-dependent transformation can also be chosen so that all jump rates remain positive even in example cases where the corresponding dynamical map breaks the property of P-divisibility, thus allowing positive unravelings without reverse quantum jumps and without auxiliary degrees of freedom in several strongly non-Markovian examples [2402.12445]. This broadens the domain in which dynamically optimized jump descriptions may be constructed.

A common source of confusion is the meaning of “optimal.” In the 2013 jump expansion, optimality refers to minimizing instantaneous partial jump rates and thereby accelerating convergence of the resummed series [1312.1167]. In the 2025 unraveling theory, optimality refers to minimizing the short-time growth of the variance of an observable over a parametric family of Poisson-driven pure-state processes [2509.19887]. In the generalized rate-operator formalism, no single variational principle is specified; instead, one optimizes a selected figure of merit under positivity constraints [2402.12445]. This suggests that DO-QJP is not a uniquely defined object across the literature, but a family resemblance term for state-adaptive quantum-jump constructions optimized for a stated criterion.

Under that reading, the unifying content of DO-QJP is the use of realization-dependent or record-dependent freedom in the jump description to improve either convergence, variance, positivity, or computational efficiency while preserving the target open-system dynamics on average.

Source: https://www.emergentmind.com/topics/dynamically-optimal-quantum-jump-process-do-qjp