Dyson series explained

In scattering theory, a part of mathematical physics, the Dyson series, formulated by Freeman Dyson, is a perturbative expansion of the time evolution operator in the interaction picture. Each term can be represented by a sum of Feynman diagrams.

This series diverges asymptotically, but in quantum electrodynamics (QED) at the second order the difference from experimental data is in the order of 10-10. This close agreement holds because the coupling constant (also known as the fine-structure constant) of QED is much less than 1.

Dyson operator

In the interaction picture, a Hamiltonian, can be split into a free part and an interacting part as .

The potential in the interacting picture is

VI(t)=

iH0(t-t0)/\hbar
e

VS(t)

-iH0(t-t0)/\hbar
e

,

where

H0

is time-independent and

VS(t)

is the possibly time-dependent interacting part of the Schrödinger picture. To avoid subscripts,

V(t)

stands for

VI(t)

in what follows.

In the interaction picture, the evolution operator is defined by the equation:

\Psi(t)=U(t,t0)\Psi(t0)

This is sometimes called the Dyson operator.

The evolution operator forms a unitary group with respect to the time parameter. It has the group properties:

U(t,t)=1,

[1]

U(t,t0)=U(t,t1)U(t1,t0),

[2]

U-1(t,t0)=U(t0,t),

U\dagger(t,t0)U(t,t0)=1

[3] and from these is possible to derive the time evolution equation of the propagator:[4]
i\hbard{dt}
U(t,t

0)\Psi(t0)=V(t)U(t,t0)\Psi(t0).

In the interaction picture, the Hamiltonian is the same as the interaction potential

H\rm=V(t)

and thus the equation can also be written in the interaction picture as

i\hbar

d{dt}
\Psi(t)

=H\rm\Psi(t)

Caution: this time evolution equation is not to be confused with the Tomonaga–Schwinger equation.

The formal solution is

U(t,t0)=1-i\hbar-1

t{dt
\int
1 V(t

1)U(t1,t0)},

which is ultimately a type of Volterra integral.

Derivation of the Dyson series

An iterative solution of the Volterra equation above leads to the following Neumann series:

\begin{align} U(t,t0)={}&1-i\hbar-1

t
\int
t0

dt1V(t1)+(-i\hbar-1

t
)
t0

dt1

t1
\int
t0

dt2V(t1)V(t2)+ … \\ &{}+(-i\hbar-1

t
)
t0

dt1\int

t1
t0

dt2

tn-1
\int
t0

dtnV(t1)V(t2)V(tn)+ … . \end{align}

Here,

t1>t2>>tn

, and so the fields are time-ordered. It is useful to introduce an operator

lT

, called the time-ordering operator, and to define

Un(t,t

-1
0)=(-i\hbar

)n

t
\int
t0

dt1

t1
\int
t0

dt2

tn-1
\int
t0

dtnlTV(t1)V(t2)V(tn).

The limits of the integration can be simplified. In general, given some symmetric function

K(t1,t2,...,tn),

one may define the integrals

Sn=\int

t
t0

dt1\int

t1
t0

dt2 …

tn-1
\int
t0

dtnK(t1,t2,...,tn).

and

In=\int

t
t0

dt1\int

t
t0

dt2 … \int

t
t0

dtnK(t1,t2,...,tn).

The region of integration of the second integral can be broken in

n!

sub-regions, defined by

t1>t2>>tn

. Due to the symmetry of

K

, the integral in each of these sub-regions is the same and equal to

Sn

by definition. It follows that

Sn=

1
n!

In.

Applied to the previous identity, this gives

U
n=(-i\hbar-1)n
n!
t
\int
t0

dt1\int

t
t0

dt2 … \int

t
t0

dtnlTV(t1)V(t2)V(tn).

Summing up all the terms, the Dyson series is obtained. It is a simplified version of the Neumann series above and which includes the time ordered products; it is the path-ordered exponential:[5]

\begin{align} U(t,t0)&=\sum

infty
n=0

Un(t,t0)\\ &=\sum

infty
n=0
(-i\hbar-1)n
n!
t
\int
t0

dt1\int

t
t0

dt2 … \int

t
t0

dtnlTV(t1)V(t2)V(tn)\\ &=lT\exp{-i\hbar-1

t{d\tau
\int
t0

V(\tau)}} \end{align}

This result is also called Dyson's formula.[6] The group laws can be derived from this formula.

Application on state vectors

The state vector at time

t

can be expressed in terms of the state vector at time

t0

, for

t>t0,

as
infty
|\Psi(t)\rangle=\sum
n=0

{(-i\hbar-1)n\overn!}\underbrace{\intdt1dtn}

t\rm\get1\ge\getn\get\rm

n
l{T}\left\{\prod
k=1
iH0tk/\hbar
e

V(tk

-iH0tk/\hbar
)e

\right\}|\Psi(t0)\rangle.

The inner product of an initial state at

ti=t0

with a final state at

tf=t

in the Schrödinger picture, for

tf>ti

is:

\begin{align} \langle\Psi(t\rm)&\mid\Psi(t\rm

infty
)\rangle=\sum
n=0

{(-i\hbar-1)n\overn!} x \ &\underbrace{\intdt1dtn}

t\rm\get1\ge\getn\get\rm

\langle\Psi(ti)\mid

-iH0(t\rm-t1)/\hbar
e

V\rm

-iH0(t1-t2)/\hbar
(t
1)e

V\rm(tn)

-iH0(tn-t\rm)/\hbar
e

\mid\Psi(ti)\rangle \end{align}

The S-matrix may be obtained by writing this in the Heisenberg picture, taking the in and out states to be at infinity:

\langle\Psi\rm\midS\mid\Psi\rm\rangle=\langle\Psi\rm

infty
\mid\sum
n=0

{(-i\hbar-1)n\overn!}\underbrace{\int

4x
d
1

4x
d
n}
t\rm\getn\ge\get1\get\rm

l{T}\left\{H\rm(x1)H\rm(x2)H\rm(xn)\right\}\mid\Psi\rm\rangle.

Note that the time ordering was reversed in the scalar product.

See also

References

Notes and References

  1. Sakurai, Modern Quantum mechanics, 2.1.10
  2. Sakurai, Modern Quantum mechanics, 2.1.12
  3. Sakurai, Modern Quantum mechanics, 2.1.11
  4. Sakurai, Modern Quantum mechanics, 2.1 pp. 69-71
  5. Sakurai, Modern Quantum Mechanics, 2.1.33, pp. 72
  6. Tong 3.20, http://www.damtp.cam.ac.uk/user/tong/qft/qft.pdf