summaryrefslogtreecommitdiff
path: root/source/know/concept/feynman-diagram/index.md
blob: 1f1c957e692c4413900c49ecd1b11ce7e122ccb2 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
---
title: "Feynman diagram"
sort_title: "Feynman diagram"
date: 2021-11-18
categories:
- Physics
- Quantum mechanics
layout: "concept"
---

When doing calculations in the context of condensed matter physics and quantum field theory,
**Feynman diagrams** graphically represent expressions
that would be tedious or error-prone to work with directly.
This article is about condensed matter physics.

Suppose we have a many-particle Hamiltonian $$\hat{H} = \hat{H}_0 + \hat{H}_1$$,
consisting of an "easy" term $$\hat{H}_0$$,
and then a "difficult" term $$\hat{H}_1$$
with time-dependent and/or interacting parts.
Let $$\Ket{\Phi_0}$$ be a known eigenstate (or superposition thereof)
of the easily solvable part $$\hat{H}_0$$,
with respect to which we will take expectation values $$\Expval{}$$.

Below, we go through the most notable components of Feynman diagrams
and how to translate them into a mathematical expression.



## Real space

The most common component is a **fermion line**, which represents
a [Green's function](/know/concept/greens-functions/) $$G^0$$
for the simple Hamiltonian $$\hat{H}_0$$.
Any type of Green's function is possible in theory (e.g. a retarded),
but usually the *causal* function is used.
Let the subscript $$I$$ refer to the
[interaction picture](/know/concept/interaction-picture/),
and $$\mathcal{T}\{\}$$ denote the
[time-ordered product](/know/concept/time-ordered-product/):

{% include image.html file="fermion-light.png" width="60%"
    alt="Fermion line diagram" %}

$$\begin{aligned}
    = i \hbar G_{s_2 s_1}^0(\vb{r}_2, t_2; \vb{r}_1, t_1)
    = \Expval{\mathcal{T} \Big\{ \hat{\Psi}_{s_2 I}(\vb{r}_2, t_2) \hat{\Psi}_{s_1 I}^\dagger(\vb{r}_1, t_1) \Big\}}
\end{aligned}$$

The arrow points in the direction of time, or more generally,
from the point of creation $$\hat{\Psi}{}^\dagger$$
to the point of annihilation $$\hat{\Psi}$$.
The dots at the ends are called **vertices**,
which represent points in space and time with a spin.
Vertices can be
**internal** (one Green's function entering AND one leaving)
or **external** (either one Green's function entering OR one leaving).

Less common is a **heavy fermion line**, representing
a causal Green's function $$G$$ for the entire Hamiltonian $$\hat{H}$$,
where the subscript $$H$$ refers to the [Heisenberg picture](/know/concept/heisenberg-picture/):

{% include image.html file="fermion-heavy.png" width="60%"
    alt="Heavy fermion line diagram" %}

$$\begin{aligned}
    = i \hbar G_{s_2 s_1}(\vb{r}_2, t_2; \vb{r}_1, t_1)
    = \Expval{\mathcal{T} \Big\{ \hat{\Psi}_{s_2 H}(\vb{r}_2, t_2) \hat{\Psi}_{s_1 H}^\dagger(\vb{r}_1, t_1) \Big\}}
\end{aligned}$$

Next, an **interaction line** or **boson line** represents
a two-body interaction operator $$\hat{W}$$ (in $$\hat{H}_1$$),
which we assume to be instantaneous, i.e. time-independent
(in quantum field theory this is *not* assumed),
hence it starts and ends at the same time,
and no arrow is drawn:

{% include image.html file="boson.png" width="60%"
    alt="Boson/interaction line diagram" %}

$$\begin{aligned}
    = \frac{1}{i \hbar} W_{s_2 s_1}(\vb{r}_2, t_2; \vb{r}_1, t_1)
    = \frac{1}{i \hbar} W(\vb{r}_2, \vb{r}_1; t_1) \: \delta(t_2 - t_1)
\end{aligned}$$

We have chosen to disallow spin flipping,
so $$W$$ does not depend on $$s_1$$ or $$s_2$$.
For reference, this function $$W$$
has a time-dependence coming only from the interaction picture,
and is to be used as follows to get the full two-body operator $$\hat{W}$$:

$$\begin{aligned}
    \hat{W}
    = \frac{1}{2} \sum_{s_1 s_2} \iint_{-\infty}^\infty \hat{\Psi}_{s_1}^\dagger(\vb{r}_1) \hat{\Psi}_{s_2}^\dagger(\vb{r}_2)
    W(\vb{r}_1, \vb{r}_2) \hat{\Psi}_{s_2}(\vb{r}_2) \hat{\Psi}_{s_1}(\vb{r}_1) \dd{\vb{r}_1} \dd{\vb{r}_2}
\end{aligned}$$

One-body (time-dependent) operators $$\hat{V}$$ in $$\hat{H}_1$$
are instead represented by a special vertex:

{% include image.html file="impurity.png" width="35%"
    alt="One-body perturbation (e.g. impurity) diagram" %}

$$\begin{aligned}
    = \frac{1}{i \hbar} V_s(\vb{r}, t)
\end{aligned}$$

Other graphical components exist representing
more complicated operators and quantities,
but these deserve their own articles.

In order for a given Feynman diagram to be valid,
it must satisfy the following criteria:

a.  Each vertex must be connected to one or two fermion lines,
    at most one of which leaves,
    and at most one of which enters.
b.  Each internal vertex contains at most one "event";
    which could be $$V$$ or $$W$$.

Finally, we need some additional rules to convert
diagrams into mathematical expressions:

1.  Disallow spin flipping by multiplying
    each internal vertex by $$\delta_{s_\mathrm{in} s_\mathrm{out}}$$.
2.  If both ends of a line are at the same time (always the case for $$W$$),
    an infinitesimal $$\eta \to 0^+$$ must be added
    to the time of all creation operators,
    so e.g. $$G(t, t) \to G(t, t\!+\!\eta)$$.
3.  Integrate over spacetime coordinates $$(\vb{r}, t)$$
    and sum over the spin $$s$$ of all internal vertices,
    but not external ones.
4.  Multiply the result by $$(-1)^F$$,
    where $$F$$ is the number of closed fermion loops.
5.  Depending on the context, additional constant factors may be required;
    sometimes they are changed on-the-fly during a calculation.

Note that rules 4 and 5 are convention,
just like the factors $$i \hbar$$ in $$G^0$$, $$G$$, $$V$$ and $$W$$;
it simply turns out to be nicer to do it this way
when using Feynman diagrams in the wild.

The combination of rules 2 and 3 means that spin
belongs to lines rather than vertices,
so that a particle with a given spin propagates
from vertex to vertex without getting flipped.



## Fourier space

If the system is time-independent and spatially uniform,
meaning it has continuous translational symmetry in time and space,
then it is useful to work in [Fourier space](/know/concept/fourier-transform/):

$$\begin{aligned}
    G_{s_2 s_1}^0(\vb{r}_2, t_2; \vb{r}_1, t_1)
    &= G_{s_1}^0(\vb{r}_2 - \vb{r}_1, t_2 - t_1) \: \delta_{s_2 s_1}
    \\
    &= \frac{\delta_{s_2 s_1}}{(2 \pi)^4} \iint_{-\infty}^\infty G_{s_1}^0(\vb{k}, \omega) \:
    e^{i \vb{k} \cdot (\vb{r}_2 - \vb{r}_1)} e^{- i \omega (t_2 - t_1)} \dd{\vb{k}} \dd{\omega}
    \\
    W_{s_2 s_1}(\vb{r}_2, t_2; \vb{r}_1, t_1)
    &= W(\vb{r}_2 - \vb{r}_1) \: \delta(t_2 - t_1)
    \\
    &= \frac{1}{(2 \pi)^4} \iint_{\infty}^\infty W(\vb{k}) \:
    e^{i \vb{k} \cdot (\vb{r}_2 - \vb{r}_1)} e^{- i \omega (t_2 - t_1)} \dd{\vb{k}} \dd{\omega}
\end{aligned}$$

Where we have used an integral representation of
the [Dirac delta function](/know/concept/dirac-delta-function/).
Note the inconsistent sign of the exponent
in the Fourier transform definitions for space and time.

Working in Fourier space allows us to simplify calculations.
Consider the following diagram and the resulting expression,
where $$\tilde{\vb{r}} = (\vb{r}, t)$$, and $$\tilde{\vb{k}} = (\vb{k}, \omega)$$:

{% include image.html file="example.png" width="40%"
    alt="Example: fermion-fermion interaction" %}

$$\begin{aligned}
    &= (i \hbar)^3 \sum_{s s'} \!\!\iint \dd{\tilde{\vb{r}}} \dd{\tilde{\vb{r}}'}
    G_{s_1's}^0(\tilde{\vb{r}}_1', \tilde{\vb{r}}) G_{s s_1}^0(\tilde{\vb{r}}, \tilde{\vb{r}}_1) \delta_{s_1 s_1'}
    W(\tilde{\vb{r}}, \tilde{\vb{r}}')
    G_{s_2' s'}^0(\tilde{\vb{r}}_2', \tilde{\vb{r}}') G_{s' s_2}^0(\tilde{\vb{r}}', \tilde{\vb{r}}_2) \delta_{s_2 s_2'}
    \\
    &= \frac{-i \hbar^3}{(2 \pi)^{20}}
    \sum_{s_1 s_2} \!\!\iint \dd{\tilde{\vb{r}}} \dd{\tilde{\vb{r}}'}
    \bigg(\! \int \dd{\tilde{\vb{k}}_2} G_{s_1}^0(\tilde{\vb{k}}_2) e^{i \tilde{\vb{k}}_2 \cdot (\tilde{\vb{r}}_1' - \tilde{\vb{r}})} \!\bigg)
    \bigg(\! \int \dd{\tilde{\vb{k}}_1} G_{s_1}^0(\tilde{\vb{k}}_1) e^{i \tilde{\vb{k}}_1 \cdot (\tilde{\vb{r}} - \tilde{\vb{r}}_1)} \!\bigg)
    \\
    &\qquad\times \bigg(\! \int \dd{\tilde{\vb{p}}} W(\tilde{\vb{p}}) e^{i \tilde{\vb{p}} \cdot (\tilde{\vb{r}}' - \tilde{\vb{r}})} \!\bigg)
    \bigg(\! \int \dd{\tilde{\vb{q}}_2} G_{s_2}^0(\tilde{\vb{q}}_2) e^{i \tilde{\vb{q}}_2 \cdot (\tilde{\vb{r}}_2' - \tilde{\vb{r}}')} \!\bigg)
    \bigg(\! \int \dd{\tilde{\vb{q}}_1} G_{s_2}^0(\tilde{\vb{q}}_1) e^{i \tilde{\vb{q}}_1 \cdot (\tilde{\vb{r}}' - \tilde{\vb{r}}_2)} \!\bigg)
    \\
    &= \frac{-i \hbar^3}{(2 \pi)^{12}}
    \sum_{s_1 s_2} \!\!\iint \dd{\tilde{\vb{k}}_1} \dd{\tilde{\vb{k}}_2}
    G_{s_1}^0(\tilde{\vb{k}}_2) G_{s_1}^0(\tilde{\vb{k}}_1)
    \iint \dd{\tilde{\vb{q}}_1} \dd{\tilde{\vb{q}}_2}
    G_{s_2}^0(\tilde{\vb{q}}_2) G_{s_2}^0(\tilde{\vb{q}}_1)
    \\
    &\qquad\times
    e^{i \tilde{\vb{k}}_2 \cdot \tilde{\vb{r}}_1' - i \tilde{\vb{k}}_1 \cdot \tilde{\vb{r}}_1
    + i \tilde{\vb{q}}_2 \cdot \tilde{\vb{r}}_2' - i \tilde{\vb{q}}_1 \cdot \tilde{\vb{r}}_2}
    \!\!\int \dd{\tilde{\vb{p}}} W(\tilde{\vb{p}})
    \bigg( \frac{1}{(2 \pi)^8} \!\!\iint \dd{\tilde{\vb{r}}} \dd{\tilde{\vb{r}}'}
    e^{i (\tilde{\vb{k}}_1 - \tilde{\vb{k}}_2 - \tilde{\vb{p}}) \cdot \tilde{\vb{r}}}
    e^{i (\tilde{\vb{q}}_1 - \tilde{\vb{q}}_2 + \tilde{\vb{p}}) \cdot \tilde{\vb{r}}'} \bigg)
    \\
    &= \frac{-i \hbar^3}{(2 \pi)^{12}}
    \sum_{s_1 s_2} \!\!\iint \dd{\tilde{\vb{k}}_1} \dd{\tilde{\vb{k}}_2}
    G_{s_1}^0(\tilde{\vb{k}}_2) G_{s_1}^0(\tilde{\vb{k}}_1)
    \iint \dd{\tilde{\vb{q}}_1} \dd{\tilde{\vb{q}}_2}
    G_{s_2}^0(\tilde{\vb{q}}_2) G_{s_2}^0(\tilde{\vb{q}}_1)
    \\
    &\qquad\times
    e^{i \tilde{\vb{k}}_2 \cdot \tilde{\vb{r}}_1' - i \tilde{\vb{k}}_1 \cdot \tilde{\vb{r}}_1
    + i \tilde{\vb{q}}_2 \cdot \tilde{\vb{r}}_2' - i \tilde{\vb{q}}_1 \cdot \tilde{\vb{r}}_2}
    \!\!\int \dd{\tilde{\vb{p}}} W(\tilde{\vb{p}})
    \: \delta(\tilde{\vb{k}}_1 \!-\! \tilde{\vb{k}}_2 \!-\! \tilde{\vb{p}})
    \: \delta(\tilde{\vb{q}}_1 \!-\! \tilde{\vb{q}}_2 \!+\! \tilde{\vb{p}})
    \\
    &= \frac{-i \hbar^3}{(2 \pi)^{12}}
    \sum_{s_1 s_2} \!\!\int \dd{\tilde{\vb{p}}} W(\tilde{\vb{p}})
    \int \dd{\tilde{\vb{k}}_1} G_{s_1}^0(\tilde{\vb{k}}_1 \!-\! \tilde{\vb{p}}) G_{s_1}^0(\tilde{\vb{k}}_1)
    \int \dd{\tilde{\vb{q}}_1} G_{s_2}^0(\tilde{\vb{q}}_1 \!+\! \tilde{\vb{p}}) G_{s_2}^0(\tilde{\vb{q}}_1)
    \\
    &\qquad\times
    e^{i \tilde{\vb{k}}_1 \cdot (\tilde{\vb{r}}_1' - \tilde{\vb{r}}_1)}
    e^{i \tilde{\vb{q}}_1 \cdot (\tilde{\vb{r}}_2' - \tilde{\vb{r}}_2)}
    e^{i \tilde{\vb{p}} \cdot (\tilde{\vb{r}}_2' - \tilde{\vb{r}}_1')}
\end{aligned}$$

Conveniently, the Dirac delta functions that appear from the integrals
represent conservation of wavevector $$\vb{k}$$ (momentum $$\hbar \vb{k}$$)
and angular frequency $$\omega$$ (energy $$\hbar \omega$$).

In Fourier space, it makes more sense
to regard the incoming energies and momenta and spins as given,
and only integrate over the internal quantities.
We thus modify the Feynman diagram rules
such that we end up with the following result:

$$\begin{aligned}
    \equiv \frac{-i \hbar^3}{(2 \pi)^4}
    \sum_{s} \!\!\int \dd{\tilde{\vb{p}}} W(\tilde{\vb{p}})
    \: G_{s_1}^0(\tilde{\vb{k}}_1 \!-\! \tilde{\vb{p}}) \: G_{s_1}^0(\tilde{\vb{k}}_1)
    \: G_{s_2}^0(\tilde{\vb{q}}_1 \!+\! \tilde{\vb{p}}) \: G_{s_2}^0(\tilde{\vb{q}}_1)
\end{aligned}$$

Therefore, we say that fermion lines represent $$i \hbar G_{s}^0(\vb{k}, \omega)$$,
interaction lines $$W(\vb{k}) / i \hbar$$, etc.,
and the other interpretation rules are modified to the following:

1.  Each line has a momentum $$\vb{k}$$ and energy $$\omega$$,
    and each fermion line has a spin $$s$$;
    these must all be conserved at each vertex.
2.  If both ends of a *fermion* line would be at the same time,
    multiply it by $$e^{i \omega \eta}$$,
    where $$\eta \to 0^+$$ is a positive infinitesimal,
    so e.g. $$G(\tau, \tau) \to e^{i \omega \eta} G(\tau, \tau)$$.
3.  Integrate over all internal $$(\vb{k}, \omega)$$,
    and sum over all internal spins $$s$$.
    Let each $$(\vb{k}, \omega)$$ integral contribute a factor $$1 / (2 \pi)^4$$.
4.  Multiply the end result by $$(-1)^F$$, where $$F$$ is the number of closed fermion loops.
5.  Depending on the context, additional constant factors may be required;
    sometimes they are changed on-the-fly during a calculation.

Note that if the diagram is linear (i.e. does not contain interactions),
then conservation removes all internal variables,
so no integrals would be needed.



## Imaginary time

Feynman diagrams are also useful when working with
[imaginary time](/know/concept/imaginary-time/).
In that case, the meaning of fermion lines is changed as follows,
involving the [Matsubara Green's function](/know/concept/matsubara-greens-function/):

$$\begin{aligned}
    i \hbar G_{s_2 s_1}^0(\vb{r}_2, t_2; \vb{r}_1, t_1)
    \:\: &\longrightarrow \:\:
    \hbar G_{s_2 s_1}^0(\vb{r}_2, \tau_2; \vb{r}_1, \tau_1)
    = \Expval{\mathcal{T} \Big\{ \hat{\Psi}_I(\vb{r}_2, \tau_2) \hat{\Psi}_I^\dagger(\vb{r}_1, \tau_1) \Big\}}
    \\
    i \hbar G_{s_2 s_1}(\vb{r}_2, t_2; \vb{r}_1, t_1)
    \:\: &\longrightarrow \:\:
    \hbar G_{s_2 s_1}(\vb{r}_2, \tau_2; \vb{r}_1, \tau_1)
    = \Expval{\mathcal{T} \Big\{ \hat{\Psi}_H(\vb{r}_2, \tau_2) \hat{\Psi}_H^\dagger(\vb{r}_1, \tau_1) \Big\}}
\end{aligned}$$

Where the time-ordering is with respect to $$\tau$$.
Interaction lines are modified like so:

$$\begin{aligned}
    \frac{1}{i \hbar} W_{s_2 s_1}(\vb{r}_2, t_2; \vb{r}_1, t_1)
    \:\: &\longrightarrow \:\:
    -\frac{1}{\hbar} W_{s_2 s_1}(\vb{r}_2, \tau_2; \vb{r}_1, \tau_1)
    = -\frac{1}{\hbar} W(\vb{r}_2, \vb{r}_1; \tau_1) \delta(\tau_2 \!-\! \tau_1)
\end{aligned}$$

One-body $$V$$-vertices are usually not used,
because they are intended for real-time-dependent operators,
but in theory they would get a factor $$-1/\hbar$$ too.

For imaginary time, the Fourier transform is defined differently,
and a distinction must be made between
fermionic Matsubara frequencies $$i \omega_n^f$$ (for $$G$$ and $$G^0$$)
and bosonic Matsubara ones $$i \omega_n^b$$ (for $$W$$).
This distinction is compatible with frequency conservation,
since a sum of two fermionic frequencies is always bosonic.
We have:

$$\begin{aligned}
    G_{s_2 s_1}^0(\vb{r}_2, \tau_2; \vb{r}_1, \tau_1)
    &= \frac{\delta_{s_2 s_1}}{(2 \pi)^3} \int_{-\infty}^\infty \frac{1}{\hbar \beta} \sum_{n = -\infty}^\infty
    G_{s_1}^0(\vb{k}, i \omega_n^f) e^{i \vb{k} \cdot (\vb{r}_2 - \vb{r}_1)} e^{- i \omega_n^f (\tau_2 - \tau_1)} \dd{\vb{k}}
    \\
    W_{s_2 s_1}(\vb{r}_2, \tau_2; \vb{r}_1, \tau_1)
    &= \frac{1}{(2 \pi)^3} \int_{-\infty}^\infty \frac{1}{\hbar \beta} \sum_{n = -\infty}^\infty
    W(\vb{k}) e^{i \vb{k} \cdot (\vb{r}_2 - \vb{r}_1)} e^{- i \omega_n^b (\tau_2 - \tau_1)} \dd{\vb{k}}
\end{aligned}$$

The interpretation in Fourier space is the same,
except that each internal integral/sum
instead gives a constant $$1 / \big(\hbar \beta (2 \pi)^3\big)$$,
and same-time fermion lines need a factor of $$e^{i \omega_n^f \eta}$$.



## References
1.  H. Bruus, K. Flensberg,
    *Many-body quantum theory in condensed matter physics*,
    2016, Oxford.