summaryrefslogtreecommitdiff
path: root/source/know/concept/greens-functions
diff options
context:
space:
mode:
Diffstat (limited to 'source/know/concept/greens-functions')
-rw-r--r--source/know/concept/greens-functions/index.md390
1 files changed, 390 insertions, 0 deletions
diff --git a/source/know/concept/greens-functions/index.md b/source/know/concept/greens-functions/index.md
new file mode 100644
index 0000000..48f1e76
--- /dev/null
+++ b/source/know/concept/greens-functions/index.md
@@ -0,0 +1,390 @@
+---
+title: "Green's functions"
+date: 2021-11-03
+categories:
+- Physics
+- Quantum mechanics
+layout: "concept"
+---
+
+In many-body quantum theory, a **Green's function**
+can be any correlation function between two given operators,
+although it is usually used to refer to the special case
+where the operators are particle creation/annihilation operators
+from the [second quantization](/know/concept/second-quantization/).
+
+They are somewhat related to
+[fundamental solutions](/know/concept/fundamental-solution/),
+which are also called *Green's functions*,
+but in general they are not the same,
+except in a special case, see below.
+
+
+## Single-particle functions
+
+If the two operators are single-particle creation/annihilation operators,
+then we get the **single-particle Green's functions**,
+for which the symbol $G$ is used.
+
+The **time-ordered** or **causal Green's function** $G_{\nu \nu'}$ is as follows,
+where $\mathcal{T}$ is the [time-ordered product](/know/concept/time-ordered-product/),
+$\nu$ and $\nu'$ are single-particle states,
+and $\hat{c}_\nu$ annihilates a particle from $\nu$, etc.:
+
+$$\begin{aligned}
+ \boxed{
+ G_{\nu \nu'}(t, t')
+ \equiv -\frac{i}{\hbar} \Expval{\mathcal{T} \Big\{ \hat{c}_{\nu}(t) \: \hat{c}_{\nu'}^\dagger(t') \Big\}}
+ }
+\end{aligned}$$
+
+The expectation value $\Expval{}$ is
+with respect to thermodynamic equilibrium.
+This is sometimes in the [canonical ensemble](/know/concept/canonical-ensemble/)
+(for some two-particle Green's functions, see below),
+but usually in the [grand canonical ensemble](/know/concept/grand-canonical-ensemble/),
+since we are adding/removing particles.
+In the latter case, we assume that the chemical potential $\mu$
+is already included in the Hamiltonian $\hat{H}$.
+Explicitly, for a complete set of many-particle states $\Ket{\Psi_n}$, we have:
+
+$$\begin{aligned}
+ G_{\nu \nu'}(t, t')
+ &= -\frac{i}{\hbar Z} \Tr\!\Big( \mathcal{T} \Big\{ \hat{c}_{\nu}(t) \: \hat{c}_{\nu'}^\dagger(t')\Big\} \: e^{- \beta \hat{H}} \Big)
+ \\
+ &= -\frac{i}{\hbar Z} \sum_{n}
+ \Matrixel{\Psi_n}{\mathcal{T} \Big\{ \hat{c}_{\nu}(t) \: \hat{c}_{\nu'}^\dagger(t')\Big\} \: e^{- \beta \hat{H}}}{\Psi_n}
+\end{aligned}$$
+
+Arguably more prevalent are
+the **retarded Green's function** $G_{\nu \nu'}^R$
+and the **advanced Green's function** $G_{\nu \nu'}^A$
+which are defined like so:
+
+$$\begin{aligned}
+ \boxed{
+ \begin{aligned}
+ G_{\nu \nu'}^R(t, t')
+ &\equiv -\frac{i}{\hbar} \Theta(t - t') \Expval{\comm{\hat{c}_{\nu}(t)}{\hat{c}_{\nu'}^\dagger(t')}_{\mp}}
+ \\
+ G_{\nu \nu'}^A(t, t')
+ &\equiv \frac{i}{\hbar} \Theta(t' - t) \Expval{\comm{\hat{c}_{\nu}(t)}{\hat{c}_{\nu'}^\dagger(t')}_{\mp}}
+ \end{aligned}
+ }
+\end{aligned}$$
+
+Where $\Theta$ is a [Heaviside function](/know/concept/heaviside-step-function/),
+and $[,]_{\mp}$ is a commutator for bosons,
+and an anticommutator for fermions.
+Depending on the context,
+we could either be in the [Heisenberg picture](/know/concept/heisenberg-picture/)
+or in the [interaction picture](/know/concept/interaction-picture/),
+hence $\hat{c}_\nu$ and $\hat{c}_{\nu'}^\dagger$ are time-dependent.
+
+Furthermore, the **greater Green's function** $G_{\nu \nu'}^>$
+and **lesser Green's function** $G_{\nu \nu'}^<$ are:
+
+$$\begin{aligned}
+ \boxed{
+ \begin{aligned}
+ G_{\nu \nu'}^>(t, t')
+ &\equiv -\frac{i}{\hbar} \Expval{\hat{c}_{\nu}(t) \: \hat{c}_{\nu'}^\dagger(t')}
+ \\
+ G_{\nu \nu'}^<(t, t')
+ &\equiv \mp \frac{i}{\hbar} \Expval{\hat{c}_{\nu'}^\dagger(t') \: \hat{c}_{\nu}(t)}
+ \end{aligned}
+ }
+\end{aligned}$$
+
+Where $-$ is for bosons, and $+$ for fermions.
+With this, the causal, retarded and advanced Green's functions
+can thus be expressed as follows:
+
+$$\begin{aligned}
+ G_{\nu \nu'}(t, t')
+ &= \Theta(t - t') \: G_{\nu \nu'}^>(t, t') + \Theta(t' - t) \: G_{\nu \nu'}^<(t, t')
+ \\
+ G_{\nu \nu'}^R(t, t')
+ &= \Theta(t - t') \big( G_{\nu \nu'}^>(t, t') - G_{\nu \nu'}^<(t, t') \big)
+ \\
+ G_{\nu \nu'}^A(t, t')
+ &= \Theta(t' - t) \big( G_{\nu \nu'}^<(t, t') - G_{\nu \nu'}^>(t, t') \big)
+\end{aligned}$$
+
+If the Hamiltonian involves interactions,
+it might be more natural to use quantum field operators $\hat{\Psi}(\vb{r}, t)$
+instead of choosing a basis of single-particle states $\psi_\nu$.
+In that case, instead of a label $\nu$,
+we use the spin $s$ and position $\vb{r}$, leading to:
+
+$$\begin{aligned}
+ G_{ss'}(\vb{r}, t; \vb{r}', t')
+ &= -\frac{i}{\hbar} \Theta(t - t') \Expval{\mathcal{T}\Big\{ \hat{\Psi}_{s}(\vb{r}, t) \hat{\Psi}_{s'}^\dagger(\vb{r}', t') \Big\}}
+ \\
+ &= \sum_{\nu \nu'} \psi_\nu(\vb{r}) \: \psi^*_{\nu'}(\vb{r}') \: G_{\nu \nu'}(t, t')
+\end{aligned}$$
+
+And analogously for $G_{ss'}^R$, $G_{ss'}^A$, $G_{ss'}^>$ and $G_{ss'}^<$.
+Note that the time-dependence is given to the old $G_{\nu \nu'}$,
+i.e. to $\hat{c}_\nu$ and $\hat{c}_{\nu'}^\dagger$,
+because we are in the Heisenberg picture.
+
+If the Hamiltonian is time-independent,
+then it can be shown that all the Green's functions
+only depend on the time-difference $t - t'$:
+
+$$\begin{gathered}
+ G_{\nu \nu'}(t, t') = G_{\nu \nu'}(t - t')
+ \\
+ G_{\nu \nu'}^R(t, t') = G_{\nu \nu'}^R(t - t')
+ \qquad \quad
+ G_{\nu \nu'}^A(t, t') = G_{\nu \nu'}^A(t - t')
+ \\
+ G_{\nu \nu'}^>(t, t') = G_{\nu \nu'}^>(t - t')
+ \qquad \quad
+ G_{\nu \nu'}^<(t, t') = G_{\nu \nu'}^<(t - t')
+\end{gathered}$$
+
+<div class="accordion">
+<input type="checkbox" id="proof-time-diff"/>
+<label for="proof-time-diff">Proof</label>
+<div class="hidden" markdown="1">
+<label for="proof-time-diff">Proof.</label>
+We will prove that the thermal expectation value
+$\expval{\hat{A}(t) \hat{B}(t')}$ only depends on $t - t'$
+for arbitrary $\hat{A}$ and $\hat{B}$,
+and it trivially follows that the Green's functions do too.
+
+In (grand) canonical equilibrium, we know that the
+[density operator](/know/concept/density-operator/)
+$\hat{\rho}$ is as follows:
+
+$$\begin{aligned}
+ \hat{\rho} = \frac{1}{Z} \exp(- \beta \hat{H})
+\end{aligned}$$
+
+The expected value of the product
+of the time-independent operators $\hat{A}$ and $\hat{B}$ is then:
+
+$$\begin{aligned}
+ \expval{\hat{A}(t) \hat{B}(t')}
+ &= \frac{1}{Z} \Tr\!\big( \hat{\rho} \hat{A}(t) \hat{B}(t') \big)
+ \\
+ &= \frac{1}{Z} \Tr\!\Big( e^{-\beta \hat{H}} e^{i t \hat{H} / \hbar} \hat{A} e^{-i t \hat{H} / \hbar}
+ e^{i t' \hat{H} / \hbar} \hat{B} e^{-i t' \hat{H} / \hbar} \Big)
+\end{aligned}$$
+
+Using that the trace $\Tr$ is invariant
+under cyclic permutations of its argument,
+and that all functions of $\hat{H}$ commute, we find:
+
+$$\begin{aligned}
+ \expval{\hat{A}(t) \hat{B}(t')}
+ = \frac{1}{Z} \Tr\!\Big( e^{-\beta \hat{H}} e^{i (t - t') \hat{H} / \hbar} \hat{A} e^{-i (t - t') \hat{H} / \hbar} \hat{B} \Big)
+\end{aligned}$$
+
+As expected, this only depends on the time difference $t - t'$,
+because $\hat{H}$ is time-independent by assumption.
+Note that thermodynamic equilibrium is crucial:
+intuitively, if the system is not in equilibrium,
+then it evolves in some transient time-dependent way.
+</div>
+</div>
+
+If the Hamiltonian is both time-independent and non-interacting,
+then the time-dependence of $\hat{c}_\nu$
+can simply be factored out as
+$\hat{c}_\nu(t) = \hat{c}_\nu \exp(- i \varepsilon_\nu t / \hbar)$.
+Then the diagonal ($\nu = \nu'$) greater and lesser Green's functions
+can be written in the form below, where $f_\nu$ is either
+the [Fermi-Dirac distribution](/know/concept/fermi-dirac-distribution/)
+or the [Bose-Einstein distribution](/know/concept/bose-einstein-distribution/).
+
+$$\begin{aligned}
+ G_{\nu \nu}^>(t, t')
+ &= -\frac{i}{\hbar} \Expval{\hat{c}_{\nu} \hat{c}_{\nu}^\dagger} \exp\!\big(\!-\! i \varepsilon_\nu (t \!-\! t') / \hbar \big)
+ \\
+ &= -\frac{i}{\hbar} (1 - f_\nu) \exp\!\big(\!-\! i \varepsilon_\nu (t \!-\! t') / \hbar \big)
+ \\
+ G_{\nu \nu}^<(t, t')
+ &= \mp \frac{i}{\hbar} \Expval{\hat{c}_{\nu}^\dagger \hat{c}_{\nu}} \exp\!\big(\!-\! i \varepsilon_\nu (t \!-\! t') / \hbar \big)
+ \\
+ &= \mp \frac{i}{\hbar} f_\nu \exp\!\big(\!-\! i \varepsilon_\nu (t \!-\! t') / \hbar \big)
+\end{aligned}$$
+
+
+## As fundamental solutions
+
+In the absence of interactions,
+we know from the derivation of
+[equation-of-motion theory](/know/concept/equation-of-motion-theory/)
+that the equation of motion of $G^R(\vb{r}, t; \vb{r}', t')$
+is as follows (neglecting spin):
+
+$$\begin{aligned}
+ i \hbar \pdv{G^R}{t}
+ = \delta(\vb{r} \!-\! \vb{r}') \: \delta(t \!-\! t')
+ + \frac{i}{\hbar} \Theta(t \!-\! t') \Expval{\Comm{\comm{\hat{H}_0}{\hat{\Psi}(\vb{r}, t)}}{\hat{\Psi}^\dagger(\vb{r}', t')}}
+\end{aligned}$$
+
+If $\hat{H}_0$ only contains kinetic energy,
+i.e. there is no external potential,
+it can be shown that:
+
+$$\begin{aligned}
+ \comm{\hat{H}_0}{\hat{\Psi}(\vb{r})}
+ = \frac{\hbar^2}{2 m} \nabla^2 \hat{\Psi}(\vb{r})
+\end{aligned}$$
+
+<div class="accordion">
+<input type="checkbox" id="proof-commH0"/>
+<label for="proof-commH0">Proof</label>
+<div class="hidden" markdown="1">
+<label for="proof-commH0">Proof.</label>
+In the second quantization,
+the Hamiltonian $\hat{H}_0$ is written like so:
+
+$$\begin{aligned}
+ \hat{H}_0
+ &= - \frac{\hbar^2}{2 m} \sum_{\nu \nu'} \hat{c}_\nu^\dagger \hat{c}_{\nu'} \Inprod{\psi_\nu}{\nabla^2 \psi_{\nu'}}
+ \\
+ &= - \frac{\hbar^2}{2 m} \sum_{\nu \nu'} \hat{c}_\nu^\dagger \hat{c}_{\nu'} \int \psi_\nu^*(\vb{r}') \: \nabla^2 \psi_{\nu'}(\vb{r}') \dd{\vb{r}'}
+ \\
+ &= - \frac{\hbar^2}{2 m}
+ \int \Big( \sum_{\nu} \psi_\nu^*(\vb{r}') \hat{c}_\nu^\dagger \Big) \Big( \nabla^2 \sum_{\nu'} \psi_{\nu'}(\vb{r}') \hat{c}_{\nu'} \Big) \dd{\vb{r}'}
+ \\
+ &= - \frac{\hbar^2}{2 m}
+ \int \hat{\Psi}^\dagger(\vb{r}') \: \nabla^2 \hat{\Psi}(\vb{r}') \dd{\vb{r}'}
+\end{aligned}$$
+
+We then insert this into the commutator that we want to prove, yielding:
+
+$$\begin{aligned}
+ \comm{\hat{H}_0}{\hat{\Psi}(\vb{r})}
+ &= - \frac{\hbar^2}{2 m} \int \Comm{\hat{\Psi}^\dagger(\vb{r}') \: \nabla^2 \hat{\Psi}(\vb{r}')}{\hat{\Psi}(\vb{r})} \dd{\vb{r}'}
+ \\
+ &= - \frac{\hbar^2}{2 m} \int \hat{\Psi}^\dagger(\vb{r}') \Comm{\nabla^2 \hat{\Psi}(\vb{r}')}{\hat{\Psi}(\vb{r})}
+ + \Comm{\hat{\Psi}^\dagger(\vb{r}')}{\hat{\Psi}(\vb{r})} \nabla^2 \hat{\Psi}(\vb{r}') \dd{\vb{r}'}
+ \\
+ &= - \frac{\hbar^2}{2 m} \sum_{\nu \nu' \nu''}
+ \Big( \hat{c}_\nu^\dagger \comm{\hat{c}_{\nu''}}{\hat{c}_{\nu'}} + \comm{\hat{c}_\nu^\dagger}{\hat{c}_{\nu'}} \hat{c}_{\nu''} \Big)
+ \psi_{\nu'}(\vb{r}) \int \psi_\nu^*(\vb{r}') \: \nabla^2 \psi_{\nu''}(\vb{r}') \dd{\vb{r}'}
+\end{aligned}$$
+
+When deriving equation-of-motion theory,
+we already showed that the following identity
+holds for both bosons and fermions:
+
+$$\begin{aligned}
+ \hat{c}_\nu^\dagger \comm{\hat{c}_{\nu''}}{\hat{c}_{\nu'}} + \comm{\hat{c}_\nu^\dagger}{\hat{c}_{\nu'}} \hat{c}_{\nu''}
+ = - \delta_{\nu \nu'} \hat{c}_{\nu''}
+\end{aligned}$$
+
+Such that the commutator can be significantly simplified to:
+
+$$\begin{aligned}
+ \comm{\hat{H}_0}{\hat{\Psi}(\vb{r})}
+ &= \frac{\hbar^2}{2 m} \sum_{\nu \nu'} \hat{c}_{\nu'}
+ \int \psi_\nu^*(\vb{r}') \: \psi_\nu(\vb{r}) \: \nabla^2 \psi_{\nu'}(\vb{r}') \dd{\vb{r}'}
+\end{aligned}$$
+
+We know that the $\psi_\nu$ form a *complete* basis,
+which implies (see [Sturm-Liouville theory](/know/concept/sturm-liouville-theory/)):
+
+$$\begin{aligned}
+ \sum_{\nu} \psi_\nu^*(\vb{r}') \: \psi_\nu(\vb{r})
+ = \delta(\vb{r} - \vb{r}')
+\end{aligned}$$
+
+With this, the commutator can be reduced even further as follows:
+
+$$\begin{aligned}
+ \comm{\hat{H}_0}{\hat{\Psi}(\vb{r})}
+ &= \frac{\hbar^2}{2 m} \sum_{\nu \nu'} \hat{c}_{\nu'}
+ \int \delta(\vb{r} - \vb{r}') \: \nabla^2 \psi_{\nu'}(\vb{r}') \dd{\vb{r}'}
+ \\
+ &= \frac{\hbar^2}{2 m} \sum_{\nu'} \hat{c}_{\nu'} \nabla^2 \psi_{\nu'}(\vb{r})
+ = \frac{\hbar^2}{2 m} \nabla^2 \hat{\Psi}(\vb{r})
+\end{aligned}$$
+</div>
+</div>
+
+After substituting this into the equation of motion,
+we recognize $G^R(\vb{r}, t; \vb{r}', t')$ itself:
+
+$$\begin{aligned}
+ i \hbar \pdv{G^R}{t}
+ &= \delta(\vb{r} \!-\! \vb{r}') \: \delta(t \!-\! t')
+ + \frac{i}{\hbar} \Theta(t \!-\! t') \Expval{\Comm{\frac{\hbar^2}{2 m} \nabla^2 \hat{\Psi}(\vb{r}, t)}{\hat{\Psi}^\dagger(\vb{r}', t')}}
+ \\
+ &= \delta(\vb{r} \!-\! \vb{r}') \: \delta(t \!-\! t') - \frac{\hbar^2}{2 m} \nabla_\vb{r}^2
+ \Big( \!-\! \frac{i}{\hbar} \Theta(t \!-\! t') \Expval{\Comm{\hat{\Psi}(\vb{r}, t)}{\hat{\Psi}^\dagger(\vb{r}', t')}} \Big)
+ \\
+ &= \delta(\vb{r} \!-\! \vb{r}') \: \delta(t \!-\! t')
+ - \frac{\hbar^2}{2 m} \nabla_\vb{r}^2 G^R(\vb{r}, t; \vb{r}', t')
+\end{aligned}$$
+
+Rearranging this leads to the following,
+which is the definition of a fundamental solution:
+
+$$\begin{aligned}
+ \Big( i \hbar \pdv{}{t}+ \frac{\hbar^2}{2 m} \nabla_\vb{r}^2 \Big) G^R(\vb{r}, t; \vb{r}', t')
+ &= \delta(\vb{r} \!-\! \vb{r}') \: \delta(t \!-\! t')
+\end{aligned}$$
+
+Therefore, the retarded Green's function
+(and, it turns out, the advanced Green's function too)
+is a fundamental solution of the Schrödinger equation
+if there is no potential,
+i.e. the Hamiltonian only contains kinetic energy.
+
+
+## Two-particle functions
+
+We generalize the above to two arbitrary operators $\hat{A}$ and $\hat{B}$,
+giving us the **two-particle Green's functions**,
+or just **correlation functions**.
+The **causal correlation function** $C_{AB}$,
+the **retarded correlation function** $C_{AB}^R$
+and the **advanced correlation function** $C_{AB}^A$ are defined as follows
+(in the Heisenberg picture):
+
+$$\begin{aligned}
+ \boxed{
+ \begin{aligned}
+ C_{AB}(t, t')
+ &\equiv -\frac{i}{\hbar} \Expval{\mathcal{T}\Big\{\hat{A}(t) \hat{B}(t')\Big\}}
+ \\
+ C_{AB}^R(t, t')
+ &\equiv -\frac{i}{\hbar} \Theta(t - t') \Expval{\comm{\hat{A}(t)}{\hat{B}(t')}_{\mp}}
+ \\
+ C_{AB}^A(t, t')
+ &\equiv \frac{i}{\hbar} \Theta(t' - t) \Expval{\comm{\hat{A}(t)}{\hat{B}(t')}_{\mp}}
+ \end{aligned}
+ }
+\end{aligned}$$
+
+Where the expectation value $\Expval{}$ is taken of thermodynamic equilibrium.
+The name *two-particle* comes from the fact that $\hat{A}$ and $\hat{B}$
+will often consist of a sum of products
+of two single-particle creation/annihilation operators.
+
+Like for the single-particle Green's functions,
+if the Hamiltonian is time-independent,
+then it can be shown that the two-particle functions
+only depend on the time-difference $t - t'$:
+
+$$\begin{aligned}
+ G_{\nu \nu'}(t, t') = G_{\nu \nu'}(t \!-\! t')
+ \qquad
+ G_{\nu \nu'}^R(t, t') = G_{\nu \nu'}^>(t \!-\! t')
+ \qquad
+ G_{\nu \nu'}^A(t, t') = G_{\nu \nu'}^<(t \!-\! t')
+\end{aligned}$$
+
+
+
+## References
+1. H. Bruus, K. Flensberg,
+ *Many-body quantum theory in condensed matter physics*,
+ 2016, Oxford.