Skip to content

Reduced Dynamics

Reduced dynamics is the time evolution assigned to a subsystem after the degrees of freedom outside that subsystem are ignored.

The central formula is:

ρS(t)=Tr⁡E[USE(t,t0)ρSE(t0)USE†(t,t0)].\rho_S(t) = \operatorname{Tr}_E \left[ U_{SE}(t,t_0)\rho_{SE}(t_0)U_{SE}^\dagger(t,t_0) \right].

The total system S+ES+E may evolve unitarily. The subsystem SS can still evolve nonunitarily because information, coherence, and energy can move into correlations with EE.

For the standard microscopic starting point that generates USEU_{SE}, see System–Bath Hamiltonians. For the terminology behind environments, baths, reservoirs, noise sources, and records, see Baths, Reservoirs, and Environments.

Let the total Hilbert space be

HSE=HS⊗HE.\mathcal H_{SE} = \mathcal H_S\otimes\mathcal H_E.

If the total state at time tt is ρSE(t)\rho_{SE}(t), then the reduced state of SS is

ρS(t)=Tr⁡EρSE(t).\rho_S(t) = \operatorname{Tr}_E\rho_{SE}(t).

If the total evolution is unitary,

ρSE(t)=USE(t,t0)ρSE(t0)USE†(t,t0).\rho_{SE}(t) = U_{SE}(t,t_0)\rho_{SE}(t_0)U_{SE}^\dagger(t,t_0).

Combining these equations gives the exact reduced dynamics formula:

ρS(t)=Tr⁡E[USE(t,t0)ρSE(t0)USE†(t,t0)].\rho_S(t) = \operatorname{Tr}_E \left[ U_{SE}(t,t_0)\rho_{SE}(t_0)U_{SE}^\dagger(t,t_0) \right].

This equation is always true when the total unitary model and initial total state are specified. It is not yet a closed equation for ρS(t)\rho_S(t) alone.

A common starting assumption is an initially uncorrelated system-environment state:

ρSE(t0)=ρS(t0)⊗ηE,\rho_{SE}(t_0) = \rho_S(t_0)\otimes\eta_E,

where ηE\eta_E is a fixed environment reference state.

Then the reduced state can be written as a map on the initial system state:

ρS(t)=Φt,t0(ρS(t0)),\rho_S(t) = \Phi_{t,t_0}\bigl(\rho_S(t_0)\bigr),

with

Φt,t0(ρ)=Tr⁡E[USE(t,t0)(ρ⊗ηE)USE†(t,t0)].\Phi_{t,t_0}(\rho) = \operatorname{Tr}_E \left[ U_{SE}(t,t_0)(\rho\otimes\eta_E) U_{SE}^\dagger(t,t_0) \right].

This map is completely positive and trace preserving. It is a quantum channel.

The factorized-state assumption is not a theorem. It is a modeling assumption that can fail when the system and environment are already correlated before the time interval of interest.

The channel form can be made explicit. Diagonalize the environment reference state:

ηE=∑rqr∣r⟩⟨r∣,qr≥0,∑rqr=1.\eta_E = \sum_r q_r |r\rangle\langle r|, \qquad q_r\ge0, \qquad \sum_r q_r=1.

Choose an orthonormal environment basis {∣μ⟩}\{|\mu\rangle\} for the final partial trace. Define operators on the system by

Kμr(t,t0)=qr ⟨μ∣USE(t,t0)∣r⟩.K_{\mu r}(t,t_0) = \sqrt{q_r}\, \langle\mu| U_{SE}(t,t_0) |r\rangle.

The bracket is taken only over the environment Hilbert space, leaving a system operator. Then

Φt,t0(ρ)=∑μ,rKμr(t,t0)ρKμr†(t,t0).\Phi_{t,t_0}(\rho) = \sum_{\mu,r} K_{\mu r}(t,t_0)\rho K_{\mu r}^\dagger(t,t_0).

The trace-preserving condition follows from unitarity:

∑μ,rKμr†Kμr=IS.\sum_{\mu,r} K_{\mu r}^\dagger K_{\mu r} = I_S.

This derivation is one of the cleanest ways to see why reduced dynamics is described by completely positive maps when the initial state is factorized.

The reduced map Φt,t0\Phi_{t,t_0} usually does not have the form

ρS⟼USρSUS†.\rho_S \longmapsto U_S\rho_SU_S^\dagger.

It may shrink Bloch vectors, suppress off-diagonal elements, transfer population, increase entropy, or drive the system toward a steady state.

This does not contradict unitarity of the larger model. The subsystem map is nonunitary because the subsystem is not closed.

The global state may retain purity and phase information, while the local reduced state appears mixed:

global unitary evolution⇏subsystem unitary evolution.\text{global unitary evolution} \quad\not\Rightarrow\quad \text{subsystem unitary evolution}.

Decoherence is a special reduced-dynamics pattern. In a two-branch model,

(c0∣0⟩+c1∣1⟩)∣Eready⟩⟶c0∣0⟩∣E0(t)⟩+c1∣1⟩∣E1(t)⟩.\left( c_0|0\rangle+c_1|1\rangle \right) |E_{\mathrm{ready}}\rangle \longrightarrow c_0|0\rangle|E_0(t)\rangle + c_1|1\rangle|E_1(t)\rangle.

Tracing out EE gives

ρ01(t)=⟨E1(t)∣E0(t)⟩ρ01(0).\rho_{01}(t) = \langle E_1(t)|E_0(t)\rangle\rho_{01}(0).

The environment overlap becomes the decoherence factor. This is reduced dynamics: the full state evolves coherently, but the subsystem loses locally accessible interference.

For the conceptual treatment, see What Is Decoherence?. For environmental monitoring, scattering, and rate estimates, see Environment-Induced Decoherence.

Energy relaxation is another reduced-dynamics pattern. If an atom emits into field modes that are not retained, the atomic reduced state may relax from an excited state toward a lower state. In an effective qubit model this can become amplitude damping.

The important point is that dephasing and dissipation are both reduced dynamics, but they are different reduced dynamics:

dephasing: phase information moves into correlations
dissipation: energy or excitation moves into the environment

The same microscopic Hamiltonian can produce both effects, depending on the coupling operators, the bath spectrum, and the energy scales.

Reduced dynamics usually means the environment is ignored. If instead some environmental record is observed, the state of SS should be conditioned on that record.

For example:

  • ignoring emitted photons gives an unconditional master equation;
  • detecting emitted photons gives quantum jumps;
  • ignoring a detector readout gives a channel;
  • retaining the readout gives a quantum instrument or trajectory.

The unconditional state is an average over records. The conditional state is the state assigned after a particular record is known. Mixing these descriptions is a common source of confusion.

If the initial total state is not factorized,

ρSE(t0)≠ρS(t0)⊗ηE,\rho_{SE}(t_0)\ne \rho_S(t_0)\otimes\eta_E,

then the reduced state at time tt still exists:

ρS(t)=Tr⁡E[USEρSE(t0)USE†].\rho_S(t) = \operatorname{Tr}_E \left[ U_{SE}\rho_{SE}(t_0)U_{SE}^\dagger \right].

But there may be no single channel Φ\Phi that maps every possible ρS(t0)\rho_S(t_0) to ρS(t)\rho_S(t) independently of the hidden correlations. Two different total states can have the same reduced state on SS and later lead to different reduced states after the same joint unitary.

One may sometimes define a map on a restricted compatibility domain, or introduce an assignment map that reconstructs ρSE\rho_{SE} from ρS\rho_S. Such maps require care and need not be completely positive on arbitrary system states. The canonical treatment is Initial Correlations.

This is why initial correlations are not a technical footnote. They determine whether reduced dynamics can be represented as a universal channel on the system alone.

The existence of maps Φt,t0\Phi_{t,t_0} does not imply Markovian dynamics. A reduced state can depend on information that flowed into the environment at earlier times and later returned.

A Markovian semigroup has the special form

Φt=etL,t≥0,\Phi_t=e^{t\mathcal L}, \qquad t\ge0,

with a generator L\mathcal L of Lindblad form under the usual assumptions. This is a much stronger structure than the exact reduced-state formula.

More generally, reduced dynamics may require:

Open-system modeling is therefore a hierarchy:

exact S+E unitary
-> exact reduced state
-> channel for factorized initial state
-> time-local master equation if justified
-> Markovian Lindblad equation under stronger assumptions

For a system field coupled to environmental fields, the Schwinger–Keldysh Bridge shows how the same partial trace becomes an influence functional with causal dissipation and fluctuation kernels.

Before using a reduced-dynamics model, ask:

  • What is the boundary between system and environment?
  • Is the initial state factorized, correlated, or only approximately factorized?
  • Is the environment state fixed independently of ρS\rho_S?
  • Is the environment ignored, measured, or fed back?
  • Is the resulting map intended to be CPTP on all system states?
  • Is a time-local master equation justified, or does the environment retain memory?
  • Are the approximations compatible with positivity and trace preservation?

These questions often determine whether a channel, master equation, trajectory, or explicit composite model is appropriate.

  • Thinking nonunitary reduced dynamics violates unitary quantum mechanics.
  • Forgetting the partial trace: ρS(t)\rho_S(t) is not obtained by simply deleting terms from the total state vector.
  • Assuming every open-system map is Markovian.
  • Assuming every reduced evolution is a channel on arbitrary ρS\rho_S even with initial correlations.
  • Treating the environment as ignored in one equation and observed in the next without changing the state-update rule.
  • Confusing a phenomenological fitted channel with a microscopic derivation.
  • Inferring a unique environment model from a reduced channel.
  • K. Kraus, States, Effects, and Operations: Fundamental Notions of Quantum Theory, Springer (1983).
  • H.-P. Breuer and F. Petruccione, The Theory of Open Quantum Systems, Oxford University Press (2002).
  • H. M. Wiseman and G. J. Milburn, Quantum Measurement and Control, Cambridge University Press (2010).
  • Á. Rivas and S. F. Huelga, Open Quantum Systems: An Introduction, Springer (2012).
  • I. de Vega and D. Alonso, “Dynamics of non-Markovian open quantum systems,” Reviews of Modern Physics 89, 015001 (2017).
  • F. A. Pollock, C. Rodríguez-Rosario, T. Frauenheim, M. Paternostro, and K. Modi, “Non-Markovian quantum processes: complete framework and efficient characterization,” Physical Review A 97, 012127 (2018).
  1. Kraus operators from a unitary dilation. Starting from
Φ(ρ)=Tr⁡E[U(ρ⊗ηE)U†],\Phi(\rho) = \operatorname{Tr}_E \left[ U(\rho\otimes\eta_E)U^\dagger \right],

with ηE=∑rqr∣r⟩⟨r∣\eta_E=\sum_rq_r|r\rangle\langle r|, derive the Kraus operators Kμr=qr⟨μ∣U∣r⟩K_{\mu r}=\sqrt{q_r}\langle\mu|U|r\rangle.

Solution

Insert the spectral decomposition of ηE\eta_E:

Φ(ρ)=∑rqrTr⁡E[U(ρ⊗∣r⟩⟨r∣)U†].\Phi(\rho) = \sum_r q_r \operatorname{Tr}_E \left[ U(\rho\otimes|r\rangle\langle r|)U^\dagger \right].

Now evaluate the partial trace in the basis {∣μ⟩}\{|\mu\rangle\}:

Φ(ρ)=∑μ,rqr⟨μ∣U∣r⟩ρ⟨r∣U†∣μ⟩.\Phi(\rho) = \sum_{\mu,r}q_r \langle\mu|U|r\rangle \rho \langle r|U^\dagger|\mu\rangle.

Defining

Kμr=qr⟨μ∣U∣r⟩K_{\mu r} = \sqrt{q_r}\langle\mu|U|r\rangle

gives

Φ(ρ)=∑μ,rKμrρKμr†.\Phi(\rho) = \sum_{\mu,r}K_{\mu r}\rho K_{\mu r}^\dagger.
  1. Trace preservation. Prove that the Kraus operators in the previous exercise satisfy ∑μ,rKμr†Kμr=IS\sum_{\mu,r}K_{\mu r}^\dagger K_{\mu r}=I_S.
Solution

Using the definition,

∑μ,rKμr†Kμr=∑μ,rqr⟨r∣U†∣μ⟩⟨μ∣U∣r⟩.\sum_{\mu,r}K_{\mu r}^\dagger K_{\mu r} = \sum_{\mu,r} q_r \langle r|U^\dagger|\mu\rangle \langle\mu|U|r\rangle.

The sum over μ\mu gives the identity on the environment:

∑μ∣μ⟩⟨μ∣=IE.\sum_\mu|\mu\rangle\langle\mu|=I_E.

Thus

∑μ,rKμr†Kμr=∑rqr⟨r∣U†U∣r⟩=∑rqrIS=IS.\sum_{\mu,r}K_{\mu r}^\dagger K_{\mu r} = \sum_r q_r \langle r|U^\dagger U|r\rangle = \sum_r q_r I_S =I_S.
  1. Initial correlations can matter. Let SS and EE be qubits. Define two initial total states
ρSE(A)=12∣00⟩⟨00∣+12∣11⟩⟨11∣,\rho_{SE}^{(A)} = \frac12|00\rangle\langle00| + \frac12|11\rangle\langle11|,

and

ρSE(B)=12∣01⟩⟨01∣+12∣10⟩⟨10∣.\rho_{SE}^{(B)} = \frac12|01\rangle\langle01| + \frac12|10\rangle\langle10|.

Both have ρS=I/2\rho_S=I/2. Let UU flip SS if E=1E=1 and do nothing if E=0E=0. Compute the final reduced state of SS for each initial total state.

Solution

The unitary acts as

∣s,e⟩⟼∣s⊕e,e⟩.|s,e\rangle \longmapsto |s\oplus e,e\rangle.

For ρSE(A)\rho_{SE}^{(A)}, the components transform as

∣00⟩→∣00⟩,∣11⟩→∣01⟩.|00\rangle\to|00\rangle, \qquad |11\rangle\to|01\rangle.

Both final components have S=0S=0, so

ρS(A)(t)=∣0⟩⟨0∣.\rho_S^{(A)}(t) = |0\rangle\langle0|.

For ρSE(B)\rho_{SE}^{(B)},

∣01⟩→∣11⟩,∣10⟩→∣10⟩.|01\rangle\to|11\rangle, \qquad |10\rangle\to|10\rangle.

Both final components have S=1S=1, so

ρS(B)(t)=∣1⟩⟨1∣.\rho_S^{(B)}(t) = |1\rangle\langle1|.

The same initial reduced state ρS=I/2\rho_S=I/2 led to different final reduced states because the hidden initial correlations were different. Therefore no map of ρS\rho_S alone can describe both cases.