Skip to content

Operator Dynamics

Operator dynamics studies time evolution as a map on observables, density operators, and operator algebras. Instead of asking only how a ket ∣ψ(t)⟩\lvert\psi(t)\rangle moves, it asks how Hamiltonians generate commutators, unitary conjugations, superoperators, product formulas, and discrete-time maps.

For a closed system with time-evolution operator U(t,t0)U(t,t_0), the two most important operator transformations are

AH(t)=U†(t,t0)AS(t)U(t,t0),A_H(t) = U^\dagger(t,t_0)A_S(t)U(t,t_0),

and

ρS(t)=U(t,t0)ρS(t0)U†(t,t0).\rho_S(t) = U(t,t_0)\rho_S(t_0)U^\dagger(t,t_0).

The first puts time dependence on observables. The second puts time dependence on density operators. They are not separate dynamics; they are dual descriptions of the same unitary evolution.

This chapter is the canonical home for the algebraic side of closed-system dynamics: commutators as generators, density-operator equations of motion, Liouvillian superoperators, operator exponentials, product formulas, and maps that evolve a system in discrete steps.

QuestionCurrent canonical pageUse it when
How does a density operator evolve in a closed system?Liouville–von Neumann EquationYou need trace, positivity, stationarity, or mixed-state dynamics.
How is a master equation treated as a linear map on operator space?Liouvillian SuperoperatorsYou need spectra, kernels, vectorization, or closed-versus-open generator language.
How do noncommuting exponentials combine or act by conjugation?Baker–Campbell–Hausdorff FormulaYou need nested commutators, similarity transforms, rotations, displacement operators, or local product errors.
How do noncommuting Hamiltonian pieces generate approximate time evolution?Trotter Product FormulaYou split H=HA+HBH=H_A+H_B for analysis, numerics, time slicing, or simulation.
How does closed dynamics work in discrete steps?Quantum Maps and Discrete-Time EvolutionYou study repeated unitary steps, Floquet maps, kicked systems, or ideal circuits.
How is closed dynamics a structure-preserving map of observables?Symmetries and Dynamical AutomorphismsYou want automorphisms, fixed-point algebras, time-translation groups, or the bridge to algebraic quantum mechanics.
How do Hamiltonians generate observable motion?Commutator DynamicsYou want commutators as generators, nested commutators, or the Poisson-bracket correspondence.
How do observables obey equations of motion in a picture?Heisenberg Equations of MotionYou want the picture-specific equation, conservation laws, or explicit time dependence.
How do density operators transform between pictures?Density Operators in Different PicturesYou compare Schrödinger, Heisenberg, and interaction-picture density matrices.
How do periodic systems evolve stroboscopically?Floquet OperatorsYou need one-period unitary maps, quasienergies, or micromotion.

Together these pages form the planned Operator Dynamics chapter.

For a time-independent Hamiltonian and a Schrödinger-picture observable with no explicit time dependence,

dAHdt=iℏ[H,AH].\frac{dA_H}{dt} = \frac{i}{\hbar}[H,A_H].

With explicit time dependence in the Schrödinger-picture operator,

dAHdt=iℏ[HH,AH]+(∂AS∂t)H.\frac{dA_H}{dt} = \frac{i}{\hbar}[H_H,A_H] + \left( \frac{\partial A_S}{\partial t} \right)_H.

For a density operator,

iℏdρdt=[H,ρ],i\hbar\frac{d\rho}{dt} = [H,\rho],

or, equivalently,

ρ˙=LH(ρ),LH(ρ)=−iℏ[H,ρ].\dot\rho = \mathcal L_H(\rho), \qquad \mathcal L_H(\rho) = -\frac{i}{\hbar}[H,\rho].

The map LH\mathcal L_H is a Liouvillian superoperator: it acts on operators rather than on state vectors. Its trace-dual action on observables is

LH†(A)=iℏ[H,A],\mathcal L_H^\dagger(A) = \frac{i}{\hbar}[H,A],

because cyclicity of trace gives

Tr⁡ ⁣[(LHρ)A]=Tr⁡ ⁣[ρ(LH†A)].\operatorname{Tr}\!\left[(\mathcal L_H\rho)A\right] = \operatorname{Tr}\!\left[\rho(\mathcal L_H^\dagger A)\right].

This duality is the compact reason density-operator dynamics and Heisenberg-picture observable dynamics give the same expectation-value evolution.

In ordinary Schrödinger evolution the Hamiltonian generates motion of state vectors. In operator dynamics the same Hamiltonian generates motion by commutation. For any observable AA, the expression [H,A][H,A] measures the failure of AA to be invariant under the Hamiltonian flow.

If [H,A]=0[H,A]=0 and AA has no explicit time dependence, then AA is conserved in the Heisenberg picture. If [H,A]≠0[H,A]\ne0, the operator changes with time even if the Schrödinger-picture symbol ASA_S is constant.

For time-independent HH, unitary conjugation can be expanded in nested commutators:

AH(t)=eiHt/ℏAe−iHt/ℏ=A+itℏ[H,A]+12!(itℏ)2[H,[H,A]]+⋯ .\begin{aligned} A_H(t) &= e^{iHt/\hbar}Ae^{-iHt/\hbar} \\ &= A + \frac{it}{\hbar}[H,A] + \frac{1}{2!} \left( \frac{it}{\hbar} \right)^2 [H,[H,A]] + \cdots . \end{aligned}

This is where the Baker–Campbell–Hausdorff machinery enters dynamics. The compact reference form is recorded in Baker–Campbell–Hausdorff, while Baker–Campbell–Hausdorff Formula emphasizes Heisenberg evolution, rotations, coherent states, and product-formula error estimates.

Density operators make mixed states and subsystems part of the same notation as pure-state dynamics. The closed-system equation

ρ˙=−iℏ[H,ρ]\dot\rho = -\frac{i}{\hbar}[H,\rho]

preserves trace, Hermiticity, positivity, purity, and the spectrum of ρ\rho. Those invariants are not optional checks; they are what distinguish unitary closed-system evolution from open-system dissipative dynamics.

The superoperator viewpoint packages the right-hand side as a linear map on operator space:

ρ⟼LH(ρ).\rho \longmapsto \mathcal L_H(\rho).

This notation becomes essential when comparing closed-system dynamics with Lindblad Equation and other open-system generators. In the closed case, the Liouvillian is only the commutator map. In the open case, additional dissipative terms appear and the generator is no longer only Hamiltonian commutation.

Unitary evolution sends each observable to another observable by conjugation:

αt(A)=U†(t)AU(t).\alpha_t(A) = U^\dagger(t)AU(t).

This map preserves products and adjoints:

αt(AB)=αt(A)αt(B),αt(A†)=αt(A)†.\alpha_t(AB) = \alpha_t(A)\alpha_t(B), \qquad \alpha_t(A^\dagger) = \alpha_t(A)^\dagger.

Thus closed-system Heisenberg evolution is an automorphism of the observable algebra. Symmetries and Dynamical Automorphisms develops this structure-preserving viewpoint. In finite dimensions it is a clean reformulation of ordinary unitary dynamics. In infinite-dimensional and algebraic settings it becomes a durable way to talk about time evolution without choosing a particular state vector representation first.

Operator dynamics also explains why products of simpler evolutions can approximate or represent complicated dynamics. If

H=HA+HB,H=H_A+H_B,

with [HA,HB]≠0[H_A,H_B]\ne0, then e−iHt/ℏe^{-iHt/\hbar} usually cannot be factorized as one finite product e−iHAt/ℏe−iHBt/ℏe^{-iH_A t/\hbar}e^{-iH_B t/\hbar}. The Trotter product formula replaces the false finite factorization by a controlled limiting product.

Discrete-time dynamics reverses the point of view. Instead of deriving a unitary from a Hamiltonian at every instant, it may start from one unitary step UFU_F or one quantum operation and study the sequence

ρn=UFnρ0(UF†)n.\rho_n = U_F^n\rho_0(U_F^\dagger)^n.

Quantum Maps and Discrete-Time Evolution collects this stroboscopic logic for Floquet theory, kicked systems, product formulas, and ideal quantum circuits, while keeping their physical assumptions distinct.

For operator motion, first read Commutator Dynamics, then Heisenberg Equations of Motion for the picture-specific derivation. For mixed-state closed dynamics, read Liouville–von Neumann Equation, then Liouvillian Superoperators for spectra and vectorized generators. For approximating exponentials of noncommuting sums, read Trotter Product Formula. For the algebraic meaning of reversible closed dynamics, read Symmetries and Dynamical Automorphisms.

If your goal is perturbation theory or scattering, combine this chapter with Interaction Picture and Time Ordering. If your goal is phase-space dynamics, pair the density-operator equation with Moyal Bracket. If your goal is open systems, treat the closed-system Liouvillian here as the Hamiltonian part of a larger generator, not as the whole story.

  • Treating eA+B=eAeBe^{A+B}=e^Ae^B as true when [A,B]≠0[A,B]\ne0.
  • Forgetting the explicit-time-derivative term in the Heisenberg equation.
  • Reversing the signs in the density-operator and observable equations of motion.
  • Calling every generator on density matrices a Liouvillian without specifying whether it is closed, Hamiltonian, dissipative, or effective.
  • Treating unitary conjugation as if it changed the spectrum of ρ\rho.
  • Confusing a finite product-formula approximation with the exact limiting formula.
  • Assuming a discrete-time unitary map always comes from a unique time-independent Hamiltonian; logarithms of unitaries are branch-dependent.
  • P. A. M. Dirac, The Principles of Quantum Mechanics, 4th ed., Oxford University Press, 1958.
  • R. Shankar, Principles of Quantum Mechanics, 2nd ed., Springer, 1994.
  • J. J. Sakurai and J. Napolitano, Modern Quantum Mechanics, 3rd ed., Cambridge University Press, 2020.
  • L. E. Ballentine, Quantum Mechanics: A Modern Development, 2nd ed., World Scientific, 2014.
  • M. Reed and B. Simon, Methods of Modern Mathematical Physics I: Functional Analysis, revised and enlarged ed., Academic Press, 1980.
  • O. Bratteli and D. W. Robinson, Operator Algebras and Quantum Statistical Mechanics 1, 2nd ed., Springer, 1987.
  1. Let AA have no explicit time dependence. Show that [H,A]=0[H,A]=0 implies AH(t)=AH(t0)A_H(t)=A_H(t_0) for a time-independent Hamiltonian.
Solution

The Heisenberg equation is

dAHdt=iℏ[H,AH].\frac{dA_H}{dt} = \frac{i}{\hbar}[H,A_H].

If [H,A]=0[H,A]=0, then AA commutes with every function of HH, including U(t,t0)=e−iH(t−t0)/ℏU(t,t_0)=e^{-iH(t-t_0)/\hbar}. Hence

AH(t)=U†AU=A.A_H(t) = U^\dagger A U = A.

Equivalently, the right-hand side of the Heisenberg equation vanishes.

  1. Use trace cyclicity to show that the Hamiltonian Liouvillian preserves Tr⁡ρ\operatorname{Tr}\rho.
Solution

For closed-system density dynamics,

ρ˙=−iℏ[H,ρ].\dot\rho = -\frac{i}{\hbar}[H,\rho].

Therefore

ddtTr⁡ρ=−iℏTr⁡(Hρ−ρH).\frac{d}{dt}\operatorname{Tr}\rho = -\frac{i}{\hbar}\operatorname{Tr}(H\rho-\rho H).

Cyclicity gives Tr⁡(Hρ)=Tr⁡(ρH)\operatorname{Tr}(H\rho)=\operatorname{Tr}(\rho H), so the derivative is zero.

  1. Explain why a finite Trotter product can be unitary even when it is not exactly equal to e−i(HA+HB)t/ℏe^{-i(H_A+H_B)t/\hbar}.
Solution

If HAH_A and HBH_B are self-adjoint, each factor

e−iHAΔt/ℏ,e−iHBΔt/ℏe^{-iH_A\Delta t/\hbar}, \qquad e^{-iH_B\Delta t/\hbar}

is unitary. A product of unitary operators is unitary. The approximation error is not a failure of norm preservation; it is the difference between the exact unitary generated by HA+HBH_A+H_B and the particular product of simpler unitary steps.

  1. A unitary map UFU_F defines ρn+1=UFρnUF†\rho_{n+1}=U_F\rho_nU_F^\dagger. Show that purity is conserved at every discrete step.
Solution

The next-step density operator satisfies

ρn+12=UFρnUF†UFρnUF†=UFρn2UF†.\rho_{n+1}^2 = U_F\rho_nU_F^\dagger U_F\rho_nU_F^\dagger = U_F\rho_n^2U_F^\dagger.

Taking the trace and using cyclicity,

Tr⁡ρn+12=Tr⁡(UFρn2UF†)=Tr⁡ρn2.\operatorname{Tr}\rho_{n+1}^2 = \operatorname{Tr}(U_F\rho_n^2U_F^\dagger) = \operatorname{Tr}\rho_n^2.

Thus a closed discrete-time unitary map preserves purity.