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.

Helpful background. Commutators supplies the algebra, while Density Operators explains mixed-state predictions. Begin with finite-dimensional matrices: every product, trace, and exponential below is then defined. Unbounded operators require compatible domains and convergence conditions, addressed in the linked treatments.

Question Explanation Use it when
How does a density operator evolve in a closed system? Liouville–von Neumann Equation You need trace, positivity, stationarity, or mixed-state dynamics.
How is a master equation treated as a linear map on operator space? Liouvillian Superoperators You need spectra, kernels, vectorization, or closed-versus-open generator language.
How do noncommuting exponentials combine or act by conjugation? Baker–Campbell–Hausdorff Formula You need nested commutators, similarity transforms, rotations, displacement operators, or local product errors.
How do noncommuting Hamiltonian pieces generate approximate time evolution? Trotter Product Formula You 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 Evolution You study repeated unitary steps, Floquet maps, kicked systems, or ideal circuits.
How is closed dynamics a structure-preserving map of observables? Symmetries and Dynamical Automorphisms You want automorphisms, fixed-point algebras, time-translation groups, or the bridge to algebraic quantum mechanics.
How do Hamiltonians generate observable motion? Commutator Dynamics You want commutators as generators, nested commutators, or the Poisson-bracket correspondence.
How do observables obey equations of motion in a picture? Heisenberg Equations of Motion You want the picture-specific equation, conservation laws, or explicit time dependence.
How do density operators transform between pictures? Density Operators in Different Pictures You compare Schrödinger, Heisenberg, and interaction-picture density matrices.
How do periodic systems evolve stroboscopically? Floquet Operators You need one-period unitary maps, quasienergies, or micromotion.

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. With the trace pairing used here, its 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 bounded operators, [H,A][H,A] determines the initial rate of change of AA under the Hamiltonian flow. For unbounded operators, a vanishing formal commutator on a small common domain does not by itself prove that the spectral projections commute or that AA is conserved.

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. These invariants are useful numerical checks. Preserving them is necessary for the intended unitary evolution, but is not sufficient to show that a numerical solution has the correct phase or Hamiltonian.

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 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.

  • 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.
  • L. E. Ballentine, Quantum Mechanics: A Modern Development, 2nd ed., World Scientific, 2014 — density operators and evolution.
  • O. Bratteli and D. W. Robinson, Operator Algebras and Quantum Statistical Mechanics 1, 2nd ed., Springer, 1987 — observable algebras and automorphisms.
  • M. Reed and B. Simon, Methods of Modern Mathematical Physics I: Functional Analysis, revised and enlarged ed., Academic Press, 1980 — operator domains and spectral calculus.
  • J. J. Sakurai and J. Napolitano, Modern Quantum Mechanics, 3rd ed., Cambridge University Press, 2020 — commutator dynamics.
  • R. Shankar, Principles of Quantum Mechanics, 2nd ed., Springer, 1994 — matrix evolution and product expansions.
  1. In a finite-dimensional Hilbert space, 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.