9 files changed, 189 insertions, 5 deletions
diff --git a/content/know/concept/random-phase-approximation/dyson.png b/content/know/concept/random-phase-approximation/dyson.png
new file mode 100644
index 0000000..e92dbf8
--- /dev/null
+++ b/content/know/concept/random-phase-approximation/dyson.png
diff --git a/content/know/concept/random-phase-approximation/index.pdc b/content/know/concept/random-phase-approximation/index.pdc
new file mode 100644
index 0000000..970a884
--- /dev/null
+++ b/content/know/concept/random-phase-approximation/index.pdc
@@ -0,0 +1,185 @@
+---
+title: "Random phase approximation"
+firstLetter: "R"
+publishDate: 2021-12-01
+categories:
+- Physics
+- Quantum mechanics
+
+date: 2021-11-15T21:01:34+01:00
+draft: false
+markup: pandoc
+---
+
+# Random phase approximation
+
+Recall that the [self-energy](/know/concept/self-energy/) $\Sigma$
+is defined as a sum of [Feynman diagrams](/know/concept/feynman-diagram/),
+which each have an order $n$ equal to the number of interaction lines.
+We consider the self-energy in the context of [jellium](/know/concept/jellium/),
+so the interaction lines $W$ represent Coulomb repulsion,
+and we use [imaginary time](/know/concept/imaginary-time/).
+
+Let us non-dimensionalize the Feynman diagrams in the self-energy,
+by measuring momenta in units of $\hbar k_F$,
+and energies in $\epsilon_F = \hbar^2 k_F^2 / (2 m)$.
+Each internal variable then gives a factor $k_F^5$,
+where $k_F^3$ comes from the 3D momentum integral,
+and $k_F^2$ from the energy $1 / \beta$:
+
+$$\begin{aligned}
+    \frac{1}{(2 \pi)^3} \int_{-\infty}^\infty \frac{1}{\hbar \beta} \sum_{n = -\infty}^\infty \cdots \:\dd{\vb{k}}
+    \:\:\sim\:\:
+    k_F^5
+\end{aligned}$$
+
+Meanwhile, every line gives a factor $1 / k_F^2$.
+The [Matsubara Green's function](/know/concept/matsubara-greens-function/) $G^0$
+for a system with continuous translational symmetry
+is found from [equation-of-motion theory](/know/concept/equation-of-motion-theory/):
+
+$$\begin{aligned}
+    W(\vb{k}) = \frac{e^2}{\varepsilon_0 |\vb{k}|^2}
+    \:\:\sim\:\:
+    \frac{1}{k_F^2}
+    \qquad \qquad
+    G_s^0(\vb{k}, i \omega_n^F)
+    = \frac{1}{i \hbar \omega_n^F - \varepsilon_\vb{k}}
+    \:\:\sim\:\:
+    \frac{1}{k_F^2}
+\end{aligned}$$
+
+An $n$th-order diagram in $\Sigma$ contains $n$ interaction lines,
+$2n\!-\!1$ fermion lines, and $n$ integrals,
+so in total it evolves as $1 / k_F^{n-2}$.
+In jellium, we know that the electron density is proportional to $k_F^3$,
+so for high densities we can rest assured that higher-order terms in $\Sigma$
+converge to zero faster than lower-order terms.
+
+However, at a given order $n$, not all diagrams are equally important.
+In a given diagram, due to momentum conservation,
+some interaction lines carry the same momentum variable.
+Because $W(\vb{k}) \propto 1 / |\vb{k}|^2$,
+small $\vb{k}$ make a large contribution,
+and the more interaction lines depend on the same $\vb{k}$,
+the larger the contribution becomes.
+
+In other words, each diagram is dominated by contributions
+from the momentum carried by the largest number of interactions.
+At order $n$, there is one diagram
+where all $n$ interactions carry the same momentum,
+and this one dominates all others at this order.
+
+The **random phase approximation** consists of removing most diagrams
+from the defintion of the full self-energy $\Sigma$,
+leaving only the single most divergent one at each order $n$,
+i.e. the ones where all $n$ interaction lines
+carry the same momentum and energy:
+
+<a href="rpasigma.png">
+<img src="rpasigma.png" style="width:91%;display:block;margin:auto;">
+</a>
+
+Where we have defined the **screened interaction** $W^\mathrm{RPA}$,
+denoted by a double wavy line:
+
+<a href="screened.png">
+<img src="screened.png" style="width:95%;display:block;margin:auto;">
+</a>
+
+Rearranging the above sequence of diagrams quickly leads to the following
+[Dyson equation](/know/concept/dyson-equation/):
+
+<a href="dyson.png">
+<img src="dyson.png" style="width:55%;display:block;margin:auto;">
+</a>
+
+In Fourier space, this equation's linear shape
+means it is algebraic, so we can write it out:
+
+$$\begin{aligned}
+    \boxed{
+        W^\mathrm{RPA}
+        = W + W \Pi_0 W^\mathrm{RPA}
+    }
+\end{aligned}$$
+
+Where we have defined the **pair-bubble** $\Pi_0$ as follows,
+with an internal wavevector $\vb{q}$, fermionic frequency $i \omega_m^F$, and spin $s$.
+Abbreviating $\tilde{\vb{k}} \equiv (\vb{k}, i \omega_n^B)$
+and $\tilde{\vb{q}} \equiv (\vb{q}, i \omega_n^F)$:
+
+<a href="pairbubble.png">
+<img src="pairbubble.png" style="width:45%;display:block;margin:auto;">
+</a>
+
+We isolate the Dyson equation for $W^\mathrm{RPA}$,
+which reveals its physical interpretation as a *screened* interaction:
+the "raw" interaction $W \!=\! e^2 / (\varepsilon_0 |\vb{k}|^2)$
+is weakened by a term containing $\Pi_0$:
+
+$$\begin{aligned}
+    W^\mathrm{RPA}(\vb{k}, i \omega_n^B)
+    = \frac{W(\vb{k})}{1 - W(\vb{k}) \: \Pi_0(\vb{k}, i \omega_n^B)}
+    = \frac{e^2}{\varepsilon_0 |\vb{k}|^2 - e^2 \Pi_0(\vb{k}, i \omega_n^B)}
+\end{aligned}$$
+
+Let us evaluate the pair-bubble $\Pi_0$ more concretely.
+The Feynman diagram translates to:
+
+$$\begin{aligned}
+    -\hbar \Pi_0(\vb{k}, i \omega_n^B)
+    &= - \sum_{s} \frac{1}{(2 \pi)^3} \int \frac{1}{\hbar \beta} \sum_{m = -\infty}^\infty
+    \hbar G_s(\vb{k} \!+\! \vb{q}, i \omega_n^B \!+\! i \omega_m^F) \: \hbar G_s(\vb{q}, i \omega_m^F) \dd{\vb{q}}
+    \\
+    &= - \frac{2 \hbar}{(2 \pi)^3} \int \frac{1}{\beta} \sum_{m = -\infty}^\infty
+    \frac{1}{i \hbar \omega_n^B + i \hbar \omega_m^F - \varepsilon_{\vb{k}+\vb{q}}} \: \frac{1}{i \hbar \omega_m^F - \varepsilon_{\vb{q}}} \dd{\vb{q}}
+\end{aligned}$$
+
+Here we recognize a [Matsubara sum](/know/concept/matsubara-sum/),
+and rewrite accordingly.
+Note that the residues of $n_F$ are $1 / (\hbar \beta)$
+when it is a function of frequency,
+and $1 / \beta$ when it is a function of energy, so:
+
+$$\begin{aligned}
+    \Pi_0(\vb{k}, i \omega_n^B)
+    &= \frac{2}{(2 \pi)^3} \int
+    \frac{n_F(\varepsilon_{\vb{k}+\vb{q}} - i \hbar \omega_n^B)}{(\varepsilon_{\vb{k}+\vb{q}} - i \hbar \omega_n^B) - \varepsilon_{\vb{q}}}
+    + \frac{n_F(\varepsilon_{\vb{q}})}{i \hbar \omega_n^B + (\varepsilon_{\vb{q}}) - \varepsilon_{\vb{k}+\vb{q}}} \dd{\vb{q}}
+    \\
+    &= \frac{2}{(2 \pi)^3} \int \frac{n_F(\varepsilon_{\vb{q}}) - n_F(\varepsilon_{\vb{k}+\vb{q}})}
+    {i \hbar \omega_n^B + \varepsilon_{\vb{q}} - \varepsilon_{\vb{k}+\vb{q}}} \dd{\vb{q}}
+\end{aligned}$$
+
+Where we have used that $n_F(\varepsilon \!+\! i \hbar \omega_n^B) = n_F(\varepsilon)$.
+Analogously to extracting the retarded Green's function $G^R(\omega)$
+from the Matsubara Green's function $G^0(i \omega_n^F)$,
+we replace $i \omega_n^F \to \omega \!+\! i \eta$,
+where $\eta \to 0^+$ is a positive infinitesimal,
+yielding the retarded pair-bubble $\Pi_0^R$:
+
+$$\begin{aligned}
+    \boxed{
+        \Pi_0^R(\vb{k}, \omega)
+        = \frac{2}{(2 \pi)^3} \int \frac{n_F(\varepsilon_{\vb{q}}) - n_F(\varepsilon_{\vb{k}+\vb{q}})}
+        {\hbar (\omega + i \eta) + \varepsilon_{\vb{q}} - \varepsilon_{\vb{k}+\vb{q}}} \dd{\vb{q}}
+    }
+\end{aligned}$$
+
+This is as far as we can go before making simplifying assumptions.
+Therefore, we leave it at:
+
+$$\begin{aligned}
+    \boxed{
+        W^\mathrm{RPA}(\vb{k}, \omega)
+        = \frac{e^2}{\varepsilon_0 |\vb{k}|^2 - e^2 \Pi_0(\vb{k}, \omega)}
+    }
+\end{aligned}$$
+
+
+
+## References
+1.  H. Bruus, K. Flensberg,
+    *Many-body quantum theory in condensed matter physics*,
+    2016, Oxford.
diff --git a/content/know/concept/random-phase-approximation/pairbubble.png b/content/know/concept/random-phase-approximation/pairbubble.png
new file mode 100644
index 0000000..c082d19
--- /dev/null
+++ b/content/know/concept/random-phase-approximation/pairbubble.png
diff --git a/content/know/concept/random-phase-approximation/rpasigma.png b/content/know/concept/random-phase-approximation/rpasigma.png
new file mode 100644
index 0000000..87ba3cc
--- /dev/null
+++ b/content/know/concept/random-phase-approximation/rpasigma.png
diff --git a/content/know/concept/random-phase-approximation/screened.png b/content/know/concept/random-phase-approximation/screened.png
new file mode 100644
index 0000000..e8b56e1
--- /dev/null
+++ b/content/know/concept/random-phase-approximation/screened.png
diff --git a/content/know/concept/self-energy/dyson.png b/content/know/concept/self-energy/dyson.png
index efa7a63..f576632 100644
--- a/content/know/concept/self-energy/dyson.png
+++ b/content/know/concept/self-energy/dyson.png
diff --git a/content/know/concept/self-energy/fullgf.png b/content/know/concept/self-energy/fullgf.png
index 3c88c6a..5767dba 100644
--- a/content/know/concept/self-energy/fullgf.png
+++ b/content/know/concept/self-energy/fullgf.png
diff --git a/content/know/concept/self-energy/index.pdc b/content/know/concept/self-energy/index.pdc
index c86f8c5..7e67143 100644
--- a/content/know/concept/self-energy/index.pdc
+++ b/content/know/concept/self-energy/index.pdc
@@ -172,7 +172,7 @@ $$\begin{aligned}
     &= \frac{\displaystyle\sum_{n = 0}^\infty \frac{1}{2^n n!}
     \bigg[ \sum_{m = 0}^{n} \frac{n!}{m! (n \!-\! m)!} \binom{1 \; \mathrm{external}}{\mathrm{order} \; m}_{\!\Sigma\mathrm{all}}
     \binom{\mathrm{0\;or\;more\;internal}}{\mathrm{total\;order} \; (n \!-\! m)}_{\!\Sigma\mathrm{all}} \bigg]}
-    {\displaystyle\sum_{n = 0}^\infty \frac{1}{2^n n!} \binom{\mathrm{0\;or\;more\;internal}}{\mathrm{total\;order} \; n}_{\!\Sigma\mathrm{all}}}
+    {-\hbar \displaystyle\sum_{n = 0}^\infty \frac{1}{2^n n!} \binom{\mathrm{0\;or\;more\;internal}}{\mathrm{total\;order} \; n}_{\!\Sigma\mathrm{all}}}
 \end{aligned}$$
 
 Where the total order is the sum of the orders of all considered diagrams,
@@ -186,8 +186,7 @@ $$\begin{aligned}
     &= \frac{\displaystyle\sum_{m = 0}^{\infty} \frac{1}{2^m m!} \binom{1 \; \mathrm{external}}{\mathrm{order} \; m}_{\!\Sigma\mathrm{all}}
     \bigg[ \sum_{n = 0}^\infty \frac{1}{2^{n-m} (n \!-\! m)!}
     \binom{\mathrm{0\;or\;more\;internal}}{\mathrm{total\;order} \; (n \!-\! m)}_{\!\Sigma\mathrm{all}} \bigg]}
-    {\displaystyle\sum_{n = 0}^\infty \frac{1}{2^n n!}
-    \binom{\mathrm{0\;or\;more\;internal}}{\mathrm{total\;order} \; n}_{\!\Sigma\mathrm{all}}}
+    {-\hbar \displaystyle\sum_{n = 0}^\infty \frac{1}{2^n n!} \binom{\mathrm{0\;or\;more\;internal}}{\mathrm{total\;order} \; n}_{\!\Sigma\mathrm{all}}}
 \end{aligned}$$
 
 Since both $n$ and $m$ start at zero,
@@ -195,7 +194,7 @@ and the sums include all possible diagrams,
 we see that the second sum in the numerator does not actually depend on $m$:
 
 $$\begin{aligned}
-    G_{ba}
+    -\hbar G_{ba}
     &= \frac{\displaystyle\sum_{m = 0}^{\infty} \frac{1}{2^m m!} \binom{1 \; \mathrm{external}}{\mathrm{order} \; m}_{\!\Sigma\mathrm{all}}
     \bigg[ \sum_{n = 0}^\infty \frac{1}{2^n n!} \binom{\mathrm{0\;or\;more\;internal}}{\mathrm{total\;order} \; n}_{\!\Sigma\mathrm{all}} \bigg]}
     {\displaystyle\sum_{n = 0}^\infty \frac{1}{2^n n!} \binom{\mathrm{0\;or\;more\;internal}}{\mathrm{total\;order} \; n}_{\!\Sigma\mathrm{all}}}
@@ -245,7 +244,7 @@ you can convince youself that $G(b,a)$ obeys
 a [Dyson equation](/know/concept/dyson-equation/) involving $\Sigma(y, x)$:
 
 <a href="dyson.png">
-<img src="dyson.png" style="width:90%;display:block;margin:auto;">
+<img src="dyson.png" style="width:95%;display:block;margin:auto;">
 </a>
 
 This makes sense: in the "normal" Dyson equation
diff --git a/content/know/concept/self-energy/selfenergy.png b/content/know/concept/self-energy/selfenergy.png
index 8eaffff..55e182e 100644
--- a/content/know/concept/self-energy/selfenergy.png
+++ b/content/know/concept/self-energy/selfenergy.png