Dyson Expansion as Formal Evolution

The Dyson expansion is the time-ordered series solution of the operator Schrödinger equation. It is the bridge between exact time-dependent evolution and perturbation theory: before any truncation, it is a formal rewriting of the evolution operator; after truncation, it becomes an approximation whose domain of validity must be checked.

The central problem is noncommutativity. If $H(t_1)$ and $H(t_2)$ do not commute, the evolution operator cannot usually be written as an ordinary exponential of $\int H(t')\,dt'$. The Dyson expansion keeps track of operator order.

Integral Equation for U

The time-evolution operator satisfies

$$i\hbar \frac{\partial}{\partial t}U(t,t_0) = H(t)U(t,t_0), \qquad U(t_0,t_0)=I.$$

Integrating from $t_0$ to $t$ gives the equivalent integral equation

$$U(t,t_0) = I - \frac{i}{\hbar} \int_{t_0}^{t} dt_1\, H(t_1)U(t_1,t_0).$$

This equation is exact. The unknown $U$ appears on both sides, so it can be solved by iteration.

Iterative Solution

Substitute the same integral equation for $U(t_1,t_0)$:

$$\begin{aligned} U(t,t_0) &= I - \frac{i}{\hbar} \int_{t_0}^{t} dt_1\,H(t_1) \\ &\quad + \left(-\frac{i}{\hbar}\right)^2 \int_{t_0}^{t} dt_1 \int_{t_0}^{t_1} dt_2\, H(t_1)H(t_2) + \cdots . \end{aligned}$$

The order of factors is fixed by the limits: $t_1$ is later than $t_2$, so $H(t_1)$ appears to the left of $H(t_2)$.

Continuing the iteration gives

$$\begin{aligned} U(t,t_0) &= I + \sum_{n=1}^{\infty} \left(-\frac{i}{\hbar}\right)^n \int_{t_0}^{t}dt_1 \int_{t_0}^{t_1}dt_2 \cdots \\ &\qquad\qquad \int_{t_0}^{t_{n-1}}dt_n\, H(t_1)H(t_2)\cdots H(t_n). \end{aligned}$$

This is the Dyson series.

Time-Ordered Exponential Notation

The same series is written compactly as

$$U(t,t_0) = \mathcal T \exp\left[ - \frac{i}{\hbar} \int_{t_0}^{t}H(t')\,dt' \right].$$

Equivalently,

$$\begin{aligned} U(t,t_0) &= I + \sum_{n=1}^{\infty} \frac{1}{n!} \left(-\frac{i}{\hbar}\right)^n \\ &\quad\times \int_{t_0}^{t}dt_1\cdots \int_{t_0}^{t}dt_n\, \mathcal T \left[ H(t_1)\cdots H(t_n) \right]. \end{aligned}$$

The factor $1/n!$ appears because the integrals run over the full $n$-dimensional time cube; the time-ordering operator divides it into ordered regions.

First and Second Order Terms

The first-order term is

$$U^{(1)}(t,t_0) = - \frac{i}{\hbar} \int_{t_0}^{t}dt_1\,H(t_1).$$

The second-order term is

$$U^{(2)}(t,t_0) = - \frac{1}{\hbar^2} \int_{t_0}^{t}dt_1 \int_{t_0}^{t_1}dt_2\, H(t_1)H(t_2).$$

The minus sign in $U^{(2)}$ comes from $(-i)^2=-1$. The operator product order is not optional. Replacing $H(t_1)H(t_2)$ with $H(t_2)H(t_1)$ changes the answer when the Hamiltonians do not commute.

Interaction-Picture Form

In perturbation theory one usually splits

$$H(t)=H_0+V(t),$$

and moves the exactly solvable $H_0$ evolution into the picture transformation. The interaction-picture state satisfies

$$i\hbar\frac{d}{dt}\lvert\psi_I(t)\rangle = V_I(t)\lvert\psi_I(t)\rangle.$$

The interaction-picture evolution operator therefore has its own Dyson series:

$$U_I(t,t_0) = \mathcal T \exp\left[ - \frac{i}{\hbar} \int_{t_0}^{t}V_I(t')\,dt' \right].$$

Detailed transition amplitudes begin when this series is truncated and matrix elements are taken. The approximation hierarchy is mapped in Time-Dependent Perturbation Theory and Transitions. Dyson Expansion for Transition Amplitudes develops ordered transition paths and intermediate-state sums, with the leading calculation continued in First-Order Transition Probability and the rate limit in Fermi’s Golden Rule.

Diagrammatic Preview

Each term in the Dyson expansion can be pictured as ordered insertions along a time line:

t0 ---- t_n ---- ... ---- t_2 ---- t_1 ---- t

The latest insertion acts first on the left side of the product. In field theory, the same bookkeeping becomes time-ordered products, Wick contractions, and Feynman diagrams. From Evolution Operators to Time-Ordered Products develops that continuation; this page remains the canonical home of the operator-series derivation.

Validity and Truncation

For a bounded Hamiltonian, define the dimensionless integrated norm

$$M=\frac{1}{\hbar}\int_{t_0}^{t}\lVert H(s)\rVert\,ds.$$

The norm of the $n$th ordered term is at most $M^n/n!$. Consequently the remainder after order $N$ obeys the useful sufficient estimate

$$\lVert R_N\rVert \le e^M\frac{M^{N+1}}{(N+1)!}.$$

This estimate is not automatically available for unbounded Hamiltonians; domains and existence theory then matter.

If a perturbation is written as $\lambda V_I(t)$, amplitudes have an expansion in powers of $\lambda$. A transition whose zeroth-order amplitude vanishes typically has probability beginning at order $\lambda^2$, not order $\lambda$. Squaring a truncated amplitude without consistent order bookkeeping can retain some higher-order terms while omitting others.

Unitarity is also order-by-order. Writing $U=I+U^{(1)}+U^{(2)}+\cdots$, self-adjoint evolution implies

$$U^{(1)}+U^{(1)\dagger}=0,$$

and

$$U^{(2)}+U^{(2)\dagger} +U^{(1)\dagger}U^{(1)}=0.$$

The Magnus Expansion instead expands the logarithm of $U$; finite Magnus truncations generated by anti-Hermitian exponents can preserve unitarity exactly, but have their own convergence conditions.

In finite-dimensional systems with bounded Hamiltonians over a finite interval, the series can often be treated as a convergent operator series. In general quantum mechanics, unbounded operators and domain questions require more care. Physics calculations usually use the Dyson series as a formal or perturbative expansion with assumptions stated by the model.

Truncating the series is not the same as exact time evolution. For a self-adjoint Hamiltonian, the full $U(t,t_0)$ is unitary. A finite truncation usually preserves unitarity only to the order being kept. If a calculation produces probabilities outside the interval from $0$ to $1$, the truncation has been used beyond its regime.

Common Mistakes

• Treating the Dyson expansion as an approximation before it is truncated.
• Dropping time ordering for noncommuting Hamiltonians.
• Reversing the order of operator factors in the nested integrals.
• Forgetting the second-order sign from $(-i)^2=-1$.
• Using first-order transition formulas when the transition probability is no longer small.
• Treating diagrammatic intuition as a substitute for specifying the Hamiltonian, picture, and ordering convention.

Cross-Links

• Time-Dependent Perturbation Theory and Transitions
• Interaction Picture for Perturbation Theory
• Dyson Expansion for Transition Amplitudes
• Time Ordering
• From Evolution Operators to Time-Ordered Products
• Time-Dependent Hamiltonians
• Interaction Picture
• Interaction-Picture Evolution Formula Card
• First-Order Transition Probability
• Fermi’s Golden Rule
• Why Dynamics Matters for QFT

References

• F. J. Dyson, “The Radiation Theories of Tomonaga, Schwinger, and Feynman,” Physical Review 75, 486-502, 1949.
• J. J. Sakurai and J. Napolitano, Modern Quantum Mechanics, 3rd ed., Cambridge University Press, 2020.
• C. Cohen-Tannoudji, B. Diu, and F. Laloë, Quantum Mechanics, Wiley, 1977.
• S. Weinberg, The Quantum Theory of Fields, Volume I, Cambridge University Press, 1995.
• M. E. Peskin and D. V. Schroeder, An Introduction to Quantum Field Theory, Addison-Wesley, 1995.

Exercises

Second-order term

Starting from the integral equation for $U(t,t_0)$, derive the second-order Dyson term.

Solution

The integral equation is

$$U(t,t_0) = I - \frac{i}{\hbar} \int_{t_0}^{t}dt_1\,H(t_1)U(t_1,t_0).$$

Now substitute the first iteration

$$U(t_1,t_0) = I - \frac{i}{\hbar} \int_{t_0}^{t_1}dt_2\,H(t_2) + O(H^2).$$

Then

$$\begin{aligned} U(t,t_0) &= I - \frac{i}{\hbar} \int_{t_0}^{t}dt_1\,H(t_1)\\ &\quad + \left(-\frac{i}{\hbar}\right)^2 \int_{t_0}^{t}dt_1 \int_{t_0}^{t_1}dt_2\, H(t_1)H(t_2) + O(H^3). \end{aligned}$$

Thus the second-order term is the ordered double integral.

Commuting Hamiltonians

Assume $[H(t_1),H(t_2)]=0$ for all times in the interval. Show that the Dyson series reduces to the ordinary exponential of $-\frac{i}{\hbar}\int H(t')\,dt'$.

Solution

If all Hamiltonians commute, then products at different times can be rearranged freely. The ordered integration region is one of $n!$ equivalent regions inside the full $n$-dimensional time cube:

$$\int_{t_0}^{t}dt_1 \int_{t_0}^{t_1}dt_2 \cdots \int_{t_0}^{t_{n-1}}dt_n\, H(t_1)\cdots H(t_n) = \frac{1}{n!} \left[ \int_{t_0}^{t}H(t')\,dt' \right]^n.$$

Substituting this into the Dyson series gives the ordinary exponential series.

First-order unitarity check

Let $K=\int_{t_0}^{t}H(t')\,dt'$ and assume $K=K^\dagger$. Show that $U_1=I-iK/\hbar$ is unitary up to first order in $K$.

Solution

Compute

$$U_1^\dagger U_1 = \left(I+\frac{iK}{\hbar}\right) \left(I-\frac{iK}{\hbar}\right) = I+\frac{K^2}{\hbar^2}.$$

There is no first-order correction to $I$. The deviation begins at second order. This is why a first-order truncation preserves unitarity only to first order, not exactly.