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.
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 |
,
H0
VS(t)
V(t)
VI(t)
In the interaction picture, the evolution operator is defined by the equation:
\Psi(t)=U(t,t0)\Psi(t0)
The evolution operator forms a unitary group with respect to the time parameter. It has the group properties:
U(t,t)=1,
U(t,t0)=U(t,t1)U(t1,t0),
U-1(t,t0)=U(t0,t),
U\dagger(t,t0)U(t,t0)=1
i\hbar | d{dt} |
U(t,t |
0)\Psi(t0)=V(t)U(t,t0)\Psi(t0).
H\rm=V(t)
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)},
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
lT
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),
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!
t1>t2> … >tn
K
Sn
Sn=
1 | |
n! |
In.
Applied to the previous identity, this gives
U | ||||
|
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.
The state vector at time
t
t0
t>t0,
infty | |
|\Psi(t)\rangle=\sum | |
n=0 |
{(-i\hbar-1)n\overn!}\underbrace{\intdt1 … dtn}
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
tf=t
tf>ti
\begin{align} \langle\Psi(t\rm)&\mid\Psi(t\rm
infty | |
)\rangle=\sum | |
n=0 |
{(-i\hbar-1)n\overn!} x \ &\underbrace{\intdt1 … dtn}
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.