Skip to content

Quantum Jump Method

Unravel a Markovian open-system master equation into stochastic pure-state trajectories with non-Hermitian evolution interrupted by random jumps, then recover density-operator observables by ensemble averaging.

Version
v1 · 2026-08-30 · History
Domain-specific #
2594
Origin domain
computational quantum physics
Subdomain
open quantum system simulation
Aliases
Monte Carlo wave-function method, MCWF method, Quantum trajectory method

Core Idea

The quantum jump method, also called the Monte Carlo wave-function (MCWF) method, simulates a Markovian open quantum system by replacing direct density-matrix evolution with an ensemble of stochastic pure-state trajectories. Between jumps, each state vector evolves under an effective non-Hermitian Hamiltonian. Its declining norm determines the probability that an environmental event occurs. When a jump is sampled, one of the system's collapse operators is applied and the state is renormalized. Averaging projectors or observables over many independent trajectories reconstructs the solution of the corresponding Lindblad master equation.

Scope of Application

The method is used in quantum optics, atomic physics, cavity and circuit QED, dissipative state preparation, quantum control, driven open systems, and many-body calculations where a Lindblad description is suitable. Common channels include spontaneous emission, photon loss, incoherent pumping, and dephasing.

It can estimate expectation values, transient states, steady behavior, correlation functions under appropriate procedures, waiting-time distributions, and conditional dynamics. Toolkits such as QuTiP expose MCWF solvers because the method parallelizes naturally across independent trajectories.

Clarity

Begin with a Lindblad equation

dρ/dt = -(i/ℏ)[H,ρ] + Σ_k(C_kρC_k† − 1/2{C_k†C_k,ρ}).

Define H_eff, evolve a normalized state for a short interval without renormalizing, and use the norm decrease to compute event probability. If no event occurs, normalize the propagated state. If one occurs, choose channel k with probability proportional to ⟨C_k†C_k⟩, apply C_k, and normalize. Repeat for many trajectories. The sample mean of trajectory observables approaches the master-equation expectation as trajectory count increases.

Manages Complexity

Density matrices scale quadratically with Hilbert dimension before operator sparsity is considered. MCWF trades that storage burden for repeated vector evolution. It also turns dissipative superoperator dynamics into familiar wave-function propagation plus localized event updates, enabling reuse of sparse-state numerical methods.

The method decomposes uncertainty. Deterministic integration error is controlled within each trajectory; stochastic estimator error is controlled by trajectory count.

Abstract Reasoning

  1. Ensemble equivalence permits substitution. Under the method's assumptions, enough trajectories recover the Lindblad density evolution. 2. Variance falls statistically, not deterministically. Standard Monte Carlo error typically decreases like the inverse square root of trajectory count. 3. Rare events demand targeted sampling. A small-probability channel may require far more trajectories than common observables. 4. Norm loss is probability bookkeeping. Non-Hermitian evolution is not merely numerical decay; it encodes survival without a jump.

Knowledge Transfer

Within open quantum physics, the method transfers across platforms because the abstract inputs are a Hamiltonian and collapse operators. The same solver structure handles atomic emission, cavity loss, spin relaxation, and dissipative many-body models after their generators are specified.

It also connects numerical simulation to continuous-measurement theory. A measurement record selects a conditional state trajectory; discarding the record returns the unconditional master equation. Exact transfer requires aligning the stochastic process with the monitoring model.

Relationships to Other Abstractions

Local relationship map for Quantum Jump MethodParents appear above the current abstraction, mutual partners to the right, and children below. Node labels state whether each abstraction is prime or domain-specific; colors identify relation types.Quantum Jump MethodDOMAINPrime abstraction: Monte Carlo Simulation — is part ofMonte CarloSimulationPRIME

Current abstraction Quantum Jump Method Domain-specific

Parents (1) — more general patterns this builds on

Hierarchy paths (4) — routes to 4 parentless roots

Neighborhood in Abstraction Space

Quantum Jump Method sits in a sparse region of the domain-specific corpus (91st percentile for distinctiveness): few abstractions share its structure, so a faithful description tends to retrieve it precisely.

Family — Unclustered & Miscellaneous (1565 abstractions)

Nearest neighbors

Computed from structural-signature embeddings · 2026-09-08