Dissipative State Preparation

Dissipative state preparation uses irreversible dynamics to drive a quantum system toward a desired state or subspace. Instead of trying to keep the environment out, one engineers selected loss, pumping, reset, or auxiliary damping so that unwanted components are removed and the target remains invariant.

The central design statement is:

make the target stationary, make all unwanted directions flow toward it,
and verify that competing noise does not dominate the preparation

This page is the focused preparation workflow. The broader design language for engineered baths and auxiliary modes is Reservoir Engineering. The fixed-point and gap terminology is explained in Steady States and Relaxation.

Stationary States and Attractive Targets

Let the designed Markovian dynamics be

$$\dot\rho = \mathcal L(\rho) = - \frac{i}{\hbar}[H,\rho] + \sum_\alpha \mathcal D[L_\alpha]\rho,$$

where

$$\mathcal D[L]\rho = L\rho L^\dagger - \frac12 \{L^\dagger L,\rho\}.$$

For a target state $\rho_\star$, stationarity means

$$\mathcal L(\rho_\star)=0.$$

For preparation from arbitrary initial states in the modeled state space, stationarity is not enough. One wants $\rho_\star$ to be the unique attractive steady state:

$$\rho(t) \longrightarrow \rho_\star \quad \text{as }t\to\infty$$

for all allowed initial states, or at least for a clearly specified basin of attraction.

For subspace preparation, the target is a projector $P_\star$ or code space rather than one density matrix. Then the goal is usually invariance plus attraction into the subspace, not collapse to one vector inside it.

Dark Is Not Enough

For a pure target $\lvert\psi_\star\rangle$, a useful sufficient dark-state condition is

$$L_\alpha\lvert\psi_\star\rangle=0 \quad \text{for all }\alpha, \qquad H\lvert\psi_\star\rangle = E_\star\lvert\psi_\star\rangle.$$

Then

$$\rho_\star = \lvert\psi_\star\rangle\langle\psi_\star\rvert$$

is stationary. But a dark state may be only one member of a large dark subspace. Preparation requires an attracting flow:

target invariant
unwanted states connected to lossy channels
no other stable traps

The practical question is therefore not only “is the target dark?” but “what happens to every orthogonal component?”

Uniqueness and the Liouvillian Gap

In finite dimensions, the clean algebraic condition for unique preparation is

$$\ker\mathcal L = \operatorname{span}\{\rho_\star\}$$

within the physical sector being modeled. This statement should be checked at the superoperator level, not guessed from one jump operator.

If the steady state is unique, the relaxation gap

$$\Delta_{\mathcal L} = \min_{\lambda_a\ne0} \left[ - \operatorname{Re}\lambda_a \right]$$

sets the asymptotic decay scale when all nonstationary modes have negative real parts,

$$\tau_{\mathrm{prep}} \sim \frac{1}{\Delta_{\mathcal L}},$$

up to prefactors, nonnormality, and initial-state overlap with slow modes. A very small gap means preparation may be formally correct but experimentally useless.

If $\ker\mathcal L$ contains more than one physical state, the dynamics prepares a manifold, not a unique state. That may be desired for a protected subspace, but it is a failure for single-state preparation.

Engineering the Flow

A dissipative preparation scheme usually combines three elements.

First, the target is invariant. The Hamiltonian and jump operators do not push it away.

Second, unwanted components are pumped into sectors that decay. This often requires coherent drives, auxiliary levels, selection rules, or multi-step transitions.

Third, the reset channel removes entropy. A photon, phonon, quasiparticle, cavity excitation, or auxiliary-state population must be dumped into degrees of freedom that do not return the same entropy to the system.

For a dark target, $L_\alpha|\psi_\star\rangle=0$ means that the jump has zero probability there, not that the target population disappears. Other states must have paths into the target. A finite-dimensional model with a unique stationary state and no remaining nondecaying modes makes this distinction precise.

Optical-Pumping Prototype

Optical pumping gives the simplest intuition. A drive excites unwanted ground states to an excited manifold; spontaneous emission returns population to lower states; a chosen dark state stops absorbing. After many cycles, population accumulates in the dark state.

In an idealized three-level picture, the target $\lvert d\rangle$ is dark:

$$H_{\mathrm{drive}}\lvert d\rangle=0, \qquad L_\alpha\lvert d\rangle=0,$$

while an unwanted state $\lvert u\rangle$ is connected to an excited state that decays with some branching probability into $\lvert d\rangle$.

The preparation rate is not just the spontaneous-emission rate. It depends on drive strength, detuning, saturation, branching ratios, off-resonant scattering, and leakage out of the target. Stronger driving can increase the pump rate but also spoil dark-state selectivity.

Entangled Targets Need More Than Collective Decay

The collective-decay example has two dark states: the two-qubit singlet and $|gg\rangle$. No amount of waiting under that dissipator alone selects the singlet. A preparation scheme needs extra drives or jumps that remove every competing stationary sector while preserving the singlet. Kraus et al. (2008) develop constructive examples of this principle for entangled targets.

Subspace Preparation

Sometimes the target is not one state but a subspace. For a code space $\mathcal C$ with projector $P_{\mathcal C}$, the requirement is

$$\rho_{\mathcal C} = P_{\mathcal C}\rho P_{\mathcal C} \quad \Longrightarrow \quad \mathcal L(\rho_{\mathcal C}) \text{ remains in }\mathcal C.$$

For autonomous error correction, one wants error sectors to be drained back into $\mathcal C$ without revealing logical information. A schematic recovery jump is

$$L_\mu = \sqrt{\Gamma_\mu}\, R_\mu P_\mu,$$

where $P_\mu$ projects onto an error syndrome sector and $R_\mu$ maps that sector back to the code space.

The important distinction is that state preparation intentionally erases information to reach one state, while subspace stabilization should preserve arbitrary logical superpositions inside the target subspace.

Fidelity and Error Budget

A common preparation score for a pure target is

$$F(t) = \langle\psi_\star\rvert \rho(t) \lvert\psi_\star\rangle.$$

For a steady-state protocol, the asymptotic score is

$$F_{\mathrm{ss}} = \langle\psi_\star\rvert \rho_{\mathrm{ss}} \lvert\psi_\star\rangle.$$

If unwanted natural noise competes with the engineered preparation, the infidelity often scales like a ratio of harmful to useful rates:

$$1-F_{\mathrm{ss}} \sim \frac{\Gamma_{\mathrm{bad}}} {\Gamma_{\mathrm{prep}}}$$

in a simple rate-limited regime. This is only a scaling estimate. Saturation, leakage, coherent errors, thermal excitations, and metastable traps can change the dependence.

For a solvable example, let $|d\rangle$ be the target and $|u\rangle$ the unwanted state. Engineered pumping $L_p=\sqrt{\Gamma_p}|d\rangle\langle u|$ competes with leakage $L_b=\sqrt{\Gamma_b}|u\rangle\langle d|$. The target population obeys

$$\dot F=\Gamma_p(1-F)-\Gamma_bF, \qquad F_{\mathrm{ss}}=\frac{\Gamma_p}{\Gamma_p+\Gamma_b}.$$

The population approaches this value at rate $\Gamma_p+\Gamma_b$. Coherences decay at half that rate, so the Liouvillian gap is $(\Gamma_p+\Gamma_b)/2$. The slowest decay rate of the full state need not be the rate visible in the chosen fidelity trace. When $\Gamma_b\ll\Gamma_p$, the formula reduces to $1-F_{\mathrm{ss}}\simeq\Gamma_b/\Gamma_p$.

Compare preparations using:

• target fidelity;
• preparation time;
• Liouvillian gap or slowest observed decay mode;
• leakage outside the intended Hilbert space;
• sensitivity to detuning, drive amplitude, and bath temperature;
• comparison with passive relaxation or coherent pulse preparation.

Check the Full State Space

Construct the generator with both desired and unwanted channels. Evaluate $\mathcal L(\rho_\star)$, find its null space, and inspect the remaining spectrum. Evolve several initial states, including states in different symmetry sectors, and compare their limiting states and convergence times. Sweep uncertain parameters and increase any oscillator cutoff: a truncation can artificially turn an unstable heated oscillator into a finite stationary state.

The computational starting point is Solving Lindblad Equations.

Relation to Other Control Methods

Dissipative preparation can be autonomous: after the engineered drives and losses are turned on, no real-time measurement record is required. That distinguishes it from Measurement-Based Feedback.

It can be implemented through routed fields or auxiliary quantum systems, in which case Coherent Feedback may be the more detailed network description before eliminating fast degrees of freedom. If strong measurement or strong dissipation primarily confines dynamics to a subspace, Quantum Zeno Dynamics gives the control language.

It can also be combined with Optimal Control: one can optimize drives or couplings to maximize the preparation gap, suppress leakage, or improve robustness. The objective should still include the dissipative steady state, not only a short-time state-transfer score.

Exercises

Stationary versus prepared

Suppose $\mathcal L(\rho_\star)=0$. What additional property is needed before saying that $\rho_\star$ is prepared from arbitrary initial states?

Solution

One needs attractivity, usually uniqueness of the physical steady state within the modeled sector and decay of all other modes. Algebraically, in a finite model one wants the physical kernel of $\mathcal L$ to contain only $\rho_\star$, with all nonzero modes having negative real parts. Otherwise the target may be stationary but not reached from generic initial states.

A dark subspace also preserves coherences

Take orthonormal states $|d\rangle,|s\rangle,|u\rangle$, one jump $L=\sqrt\Gamma|d\rangle\langle u|$, and no Hamiltonian. Find all stationary density operators.

Solution

Stationarity forces $\rho_{uu}=0$, since $\dot\rho_{uu}=-\Gamma\rho_{uu}$. Positivity then forces $\rho_{du}=\rho_{su}=0$. Every density matrix supported on $\operatorname{span}\{|d\rangle,|s\rangle\}$ is stationary, including coherent superpositions: $\rho_{ds}$ is unchanged. The physical stationary set is the entire Bloch ball of that two-dimensional subspace, not just mixtures of its two basis projectors. Population initially in $|u\rangle$ is added to $|d\rangle$.

Gap and preparation time

If the smallest nonzero decay rate of a diagonalizable Liouvillian is $\Delta_{\mathcal L}$, what preparation time scale should you expect?

Solution

The leading estimate is $\tau_{\mathrm{prep}}\sim1/\Delta_{\mathcal L}$. This estimate can be modified by prefactors, nonnormality, and weak overlap of the initial state or observable with the slow mode, but a small gap still signals slow convergence.

State versus subspace

Why is pumping every state into $\lvert0_L\rangle$ not autonomous error correction for a logical qubit?

Solution

A logical qubit must preserve arbitrary superpositions $\alpha\lvert0_L\rangle+\beta\lvert1_L\rangle$. Pumping everything into $\lvert0_L\rangle$ erases the coefficient $\beta$ and destroys the encoded information. Error correction should remove error syndromes while preserving the logical state inside the code space.

References

• J. T. Barreiro, M. Müller, P. Schindler, D. Nigg, T. Monz, M. Chwalla, M. Hennrich, C. F. Roos, P. Zoller, and R. Blatt, “An open-system quantum simulator with trapped ions,” Nature 470, 486–491 (2011).
• S. Diehl, A. Micheli, A. Kantian, B. Kraus, H. P. Büchler, and P. Zoller, “Quantum states and phases in driven open quantum systems with cold atoms,” Nature Physics 4, 878–883 (2008).
• B. Kraus, H. P. Büchler, S. Diehl, A. Kantian, A. Micheli, and P. Zoller, “Preparation of entangled states by quantum Markov processes,” Physical Review A 78, 042307 (2008). doi:10.1103/PhysRevA.78.042307.
• Y. Lin, J. P. Gaebler, F. Reiter, T. R. Tan, R. Bowler, A. S. Sørensen, D. Leibfried, and D. J. Wineland, “Dissipative production of a maximally entangled steady state of two quantum bits,” Nature 504, 415–418 (2013).
• J. F. Poyatos, J. I. Cirac, and P. Zoller, “Quantum reservoir engineering with laser cooled trapped ions,” Physical Review Letters 77, 4728–4731 (1996).
• F. Reiter and A. S. Sørensen, “Effective operator formalism for open quantum systems,” Physical Review A 85, 032111 (2012). doi:10.1103/PhysRevA.85.032111.
• F. Verstraete, M. M. Wolf, and J. I. Cirac, “Quantum computation and quantum-state engineering driven by dissipation,” Nature Physics 5, 633–636 (2009).