# Contents

Linear Response and Dynamical Correlation Functions
5.1  Dynamics of the Statistical Operator: von Neumann's Equation
5.2  Linear Response: Dynamical Susceptibility
5.3  Spectral Representation
5.4  Symmetries
5.5  Mori's Scalar Product
5.6  The Liouville Operator
Index

# List of Figures

1  Dynamical processes in physical systems
Integration contour for the $\Theta$-function

## 5Linear Response and Dynamical Correlation Functions

Many dynamical processes induced by external perturbations of physical systems can be divided into three categories, namely (i) driving by external forces, (ii) relaxation from a partial to a global equilibrium, and (iii) probing by scattered 'particles'. See fig. 1
Figure 1: Dynamical processes in physical systems
All of these experimentally relevant situation can be described in terms of certain so-called correlation functions all of which are intimately related to each other.

 For simplicity of notation we set $\hslash ={k}_{B}=1$ in the remainder of this chapter.

### 5.1Dynamics of the Statistical Operator: von Neumann's Equation

Consider a quantum system
 $\stackrel{~}{H}=H+W\left(t\right)$ $\left(1\right)$
and an arbitrary ensemble, described by a statistical operator
 ${\rho }_{S}\left(t\right)=\sum _{n}{\rho }_{n}|n,t⟩⟨n,t| ,$ $\left(2\right)$
where $|n,t⟩$ set the
microstates in the Schrödinger picture with probabilities ${\rho }_{n}$. The time evolution of the microstates is governed by Schrödigner's equation
 $i{\partial }_{t}|n,t⟩=\stackrel{~}{H}|n,t⟩ .$ $\left(3\right)$
From this we get
von Neumann's equation of motion for $\rho \left(t\right)$ in the Schrödinger picture
 ${\partial }_{t}{\rho }_{S}\left(t\right)=-i\sum _{n}{\rho }_{n}\left[\left(\stackrel{~}{H}|n,t⟩\right)⟨n,t|-|n,t⟩\left(⟨n,t|\stackrel{~}{H}\right)\right]=-i\left[\stackrel{~}{H},{\rho }_{S}\left(t\right)\right] .$ $\left(4\right)$
This can be translated into
Heisenberg's picture ${A}_{H}\left(t\right)={e}^{i\stackrel{~}{H}t}A\left(t\right){e}^{-i\stackrel{~}{H}t}$, $|n,t{⟩}_{S}={e}^{i\stackrel{~}{H}t}|n,t⟩=|n,0{⟩}_{S}$ easily
 ${\rho }_{H}\left(t\right)={\rho }_{S}\left(0\right) ,$ $\left(5\right)$
and finally in the
Dirac or interaction picture ${A}_{D}\left(t\right)={e}^{\mathrm{iHt}}A\left(t\right){e}^{-\mathrm{iHt}}$, $|n,t{⟩}_{D}={e}^{\mathrm{iHt}}|n,t⟩$
 $\begin{array}{ccc}\multicolumn{1}{c}{{\partial }_{t}{\rho }_{D}\left(t\right)}& =\hfill & \mathrm{iH}{\rho }_{D}\left(t\right)-i{\rho }_{D}\left(t\right)H-i\left[H+{W}_{D}\left(t\right),{\rho }_{D}\left(t\right)\right]\hfill \\ \multicolumn{1}{c}{}& =\hfill & -i\left[{W}_{D}\left(t\right),{\rho }_{D}\left(t\right)\right] ,\hfill & \hfill \left(6\right)\end{array}$

where ${W}_{D}\left(t\right)={e}^{\mathrm{iHt}}W\left(t\right){e}^{-\mathrm{iHt}}$ is the interaction in Dirac's picture.
Eqn. (6) seem particularly suited for perturbative approximations since the r.h.s. is proportional to $W$. We consider a case, in which $W\left(t\right)$ is adiabatically switched in starting from $t=-\infty$
 $W\left(t\to -\infty \right)=0\mathrm{ }ö\mathrm{ff}" .$ $\left(7\right)$
Therefore at $t=-\infty$ the system is in
equilibrium and
 $\left[{\rho }_{S}\left(-\infty \right),H\right]=0 ,$ $\left(8\right)$
with a
canonical distribution function
 $\rho \equiv {\rho }_{S}\left(-\infty \right)=\frac{1}{Z}{e}^{-\beta H}={\rho }_{D}\left(-\infty \right) ,$ $\left(9\right)$
where the last equality follows from eqn. (8) and ${A}_{D}\left(t\right)={e}^{\mathrm{iHt}}A\left(t\right){e}^{-\mathrm{iHt}}$. From eqn. (6) we therefore get
 $\begin{array}{ccc}\multicolumn{1}{c}{{\rho }_{D}\left(t\right)}& =\hfill & \rho -i{\int }_{-\infty }^{t}\left[{W}_{D}\left(t\text{'}\right),{\rho }_{D}\left(t\right)\right]\mathrm{dt}=\hfill \\ \multicolumn{1}{c}{}& =\hfill & \rho -i{\int }_{-\infty }^{t}\left[{W}_{D}\left(t\text{'}\right),\rho \right]\mathrm{dt}+O\left({W}^{2}\right)\hfill & \hfill \left(10\right)\end{array}$

### 5.2Linear Response: Dynamical Susceptibility

Let us consider a situation in which the perturbation $W\left(t\right)$ is due to generalized external forces
 $W\left(t\right)=-\sum _{\nu }{A}_{\nu }{f}_{\nu }\left(t\right) ,$ $\left(11\right)$
where the ${f}_{\nu }\left(t\right)$ correspond to the forces and the operators ${A}_{\nu }$ refer to the generalized coordinates. Many such pairs of $\left({A}_{\nu },{f}_{\nu }\left(t\right)\right)$ are relevant in practical applications: (real space, momentum), (magnetization, magnetic field), (electric dipole, electric field), (electric current density, vector potential), (charge density, electric field), and so on. We now consider the response of the expectation value of $⟨{A}_{\mu }^{+}\left(t\right)⟩$ to linear order in the $W$
 $\begin{array}{ccc}\multicolumn{1}{c}{\delta ⟨{A}_{\mu }^{+}\left(t\right)⟩}& =\hfill & ⟨{A}_{\mu }^{+}\left(t\right)⟩-⟨{A}_{\mu }^{+}\left(-\infty \right)⟩\hfill \\ \multicolumn{1}{c}{}& =\hfill & \mathrm{Tr}\left\{{A}_{\mu D}^{+}\left(t\right)\left({\rho }_{D}\left(t\right)-{\rho }_{D}\left(-\infty \right)\right)\right\}\hfill \\ \multicolumn{1}{c}{}& =\hfill & -i{\int }_{-\infty }^{t}\mathrm{Tr}\left\{{A}_{\mu D}^{+}\left(t\right)\left[{W}_{D}\left(t\text{'}\right),\rho \right]\right\}\mathrm{dt}+O\left({W}^{2}\right)\hfill \\ \multicolumn{1}{c}{}& =\hfill & i\sum _{\nu }{\int }_{-\infty }^{t}\mathrm{Tr}\left\{\rho \left[{A}_{\mu D}^{+}\left(t\right),{A}_{\nu D}\left(t\text{'}\right)\right]\right\} {f}_{\nu }\left(t\text{'}\right)\mathrm{dt}\text{'}\hfill \\ \multicolumn{1}{c}{}& =\hfill & i\sum _{\nu }{\int }_{-\infty }^{t}⟨\left[{A}_{\mu D}^{+}\left(t\right),{A}_{\nu D}\left(t\text{'}\right)\right]{⟩}_{\rho } {f}_{\nu }\left(t\text{'}\right)\mathrm{dt}\text{'} ,\hfill & \hfill \left(12\right)\end{array}$

where we have dropped terms of $O\left({W}^{2}\right)$.

 In the remainder of section 5 time dependencies will usually refer to the Dirac picture, or speaking differently the Heisenberg picture for $W=0$. Moreover statistical averages will be formed w.r.t. to the equilibrium statistical operator $\rho$ at $t=-\infty$ from eqn. (9). We will therefore drop the subscripts ' $D$' and ' $\rho$' in expressions like eqn. (12).

Using eqn. (12) we define the dynamical susceptibility
 ${\chi }_{\mu \nu }\left(t,t\text{'}\right)=i\Theta \left(t-t\text{'}\right)⟨\left[{A}_{\mu }^{+}\left(t\right),{A}_{\nu }\left(t\text{'}\right)\right]⟩\equiv 2i\Theta \left(t-t\text{'}\right){\chi }_{\mu \nu }^{\text{'}\text{'}}\left(t,t\text{'}\right) .,$ $\left(13\right)$
where we have introduced to so-called spectral function ${\chi }_{\mu \nu }^{\text{'}\text{'}}\left(t,t\text{'}\right)$. Because of the $\Theta$-function we may write
 $\delta ⟨{A}_{\mu }^{+}\left(t\right)⟩=\sum _{\nu }{\int }_{-\infty }^{\infty }{\chi }_{\mu \nu }\left(t,t\text{'}\right){f}_{\nu }\left(t\text{'}\right)\mathrm{dt}\text{'} .$ $\left(14\right)$
Obviously this extends the meaning of a susceptibility from static thermodynamics to dynamical phenomena. If we assume that the coordinates ${A}_{\mu }$ are time-independent in the Schrödinger picture, then we have
 $\begin{array}{ccc}\multicolumn{1}{c}{⟨\left[{A}_{\mu }^{+}\left(t\right),{A}_{\nu }\left(t\text{'}\right)\right]⟩}& =\hfill & \mathrm{Tr}\left\{\rho \left[{e}^{\mathrm{iHt}}{A}_{\mu }^{+}{e}^{-\mathrm{iHt}},{e}^{\mathrm{iHt}\text{'}}{A}_{\nu }{e}^{-\mathrm{iHt}}\right]\right\}\hfill \\ \multicolumn{1}{c}{}& =\hfill & \mathrm{Tr}\left\{\rho \left[{e}^{\mathrm{iH}\left(t-t\text{'}\right)}{A}_{\mu }^{+}{e}^{-\mathrm{iH}\left(t-t\text{'}\right)},{A}_{\nu }\right]\right\}\hfill \\ \multicolumn{1}{c}{}& =\hfill & ⟨\left[{A}_{\mu }^{+}\left(t-t\text{'}\right),{A}_{\nu }\right]⟩ ,\hfill & \hfill \left(15\right)\end{array}$

i.e. the dynamical susceptibility is a function of the time-difference only
 ${\chi }_{\mu \nu }^{\left(\text{'}\text{'}\right)}\left(t,t\text{'}\right)={\chi }_{\mu \nu }^{\left(\text{'}\text{'}\right)}\left(t-t\text{'},0\right)\equiv {\chi }_{\mu \nu }^{\left(\text{'}\text{'}\right)}\left(t-t\text{'}\right) .$ $\left(16\right)$
Example 1 Let us consider a harmonic oscillator under the influence of an external driving force which acts on its displacement coordinate

For bosons we ave $\left[b,{b}^{+}\right]=1$. Note that
 ${b}^{+}+b\propto \stackrel{^}{r} ,$ $\left(18\right)$
where $\stackrel{^}{r}={\stackrel{^}{r}}^{+}$ is the displacement coordinate of the oscillator. The time dependence in the Dirac picture is obtained from ${b}_{D}\left(t\right)={e}^{\mathrm{iHt}}{\mathrm{be}}^{-\mathrm{iHt}}$, or
 $\begin{array}{ccc}\multicolumn{1}{c}{{\partial }_{t}{b}_{D}\left(t\right)}& =\hfill & {\mathrm{ie}}^{\mathrm{iHt}}\left[H,b\right]{e}^{-\mathrm{iHt}}=i{\omega }_{0}{e}^{\mathrm{iHt}}\left[{b}^{+}b,b\right]{e}^{-\mathrm{iHt}}\hfill \\ \multicolumn{1}{c}{}& =\hfill & i{\omega }_{0}{e}^{\mathrm{iHt}}\left({b}^{+}\left[b,b\right]+\left[{b}^{+},b\right]b\right){e}^{-\mathrm{iHt}}\hfill \\ \multicolumn{1}{c}{}& =\hfill & -i{\omega }_{0}{e}^{\mathrm{iHt}}{\mathrm{be}}^{-\mathrm{iHt}}=-i{\omega }_{0}{b}_{D}\left(t\right) .\hfill & \hfill \left(19\right)\end{array}$

This DEQN can be solved by direct integration and taking into account that ${b}_{D}\left(0\right)=b$
 $\begin{array}{cccc}\multicolumn{1}{c}{{b}_{D}\left(t\right)}& =\hfill & {e}^{-i{\omega }_{0}t}b\hfill & \hfill \left(20\right)\\ \multicolumn{1}{c}{{b}_{D}^{+}\left(t\right)}& =\hfill & {e}^{i{\omega }_{0}t}{b}^{+}\hfill & \hfill \left(21\right)\end{array}$

where the hermitean conjugate has yield a corresponding equation for the creation operator. Inserting this into eqn. (15) we get an expression for the displacement susceptibility
 $\begin{array}{ccc}\multicolumn{1}{c}{{\chi }_{\mathrm{rr}}\left(t\right)}& =\hfill & i\Theta \left(t\right)⟨\left[{e}^{i{\omega }_{0}t}{b}^{+}+{e}^{-i{\omega }_{0}t}b,{b}^{+}+b\right]⟩\hfill \\ \multicolumn{1}{c}{}& =\hfill & i\Theta \left(t\right)\left({e}^{i{\omega }_{0}t}⟨\left[{b}^{+},b\right]⟩+{e}^{-i{\omega }_{0}t}⟨\left[b,{b}^{+}\right]⟩\right)\hfill \\ \multicolumn{1}{c}{}& =\hfill & i\Theta \left(t\right)\left({e}^{-i{\omega }_{0}t}-{e}^{i{\omega }_{0}t}\right) .\hfill & \hfill \left(22\right)\end{array}$

Note in passing that it is due to the particular simplicity of this example, that the dynamical susceptibility is completely temperature independent.

### 5.3Spectral Representation

The description of dynamical phenomena involves many different types of functions, mathematically similar to eqn. (13). As we continue we will find, that these function are all determined by the spectral function. Let us consider the Laplace transform of ${\chi }_{\mu \nu }\left(t\right)$, dropping the subscripts for simplicity
 $\begin{array}{cccc}\multicolumn{1}{c}{\chi \left(z\right)}& =\hfill & {\int }_{0}^{\infty }\mathrm{dt} {e}^{\mathrm{izt}}\chi \left(t\right)={\int }_{-\infty }^{\infty }\mathrm{dt} 2{\mathrm{ie}}^{\mathrm{izt}}\Theta \left(t\right){\chi }^{\text{'}\text{'}}\left(t\right) ,\hfill & \hfill \left(23\right)\end{array}$

where $\mathrm{Im}\left(z\right)>0$. Note that the following applies to
any pair of functions $\chi \left(t\right)$ and ${\chi }^{\text{'}\text{'}}\left(t\right)$ related trough eqn. (13). Now, from the residue theorem (see fig. 2)
Figure 2: Integration contour for the $\Theta$-function

 ${\int }_{-\infty }^{\infty }\frac{d\omega }{\pi }\frac{{e}^{i\omega t}}{\omega -z}=\left\{\begin{array}{ccc}\hfill 2{\mathrm{ie}}^{\mathrm{izt}}\hfill & \hfill \hfill & \hfill t>0\hfill \\ \hfill 0\hfill & \hfill \hfill & \hfill t<0\hfill \end{array} ,$ $\left(24\right)$
i.e. we may write
 $\chi \left(z\right)={\int }_{-\infty }^{\infty }\mathrm{dt} {\chi }^{\text{'}\text{'}}\left(t\right){\int }_{-\infty }^{\infty }\frac{d\omega }{\pi }\frac{{e}^{i\omega t}}{\omega -z}={\int }_{-\infty }^{\infty }\frac{d\omega }{\pi }\frac{{\chi }^{\text{'}\text{'}}\left(\omega \right)}{\omega -z} ,$ $\left(25\right)$
where ${\chi }^{\text{'}\text{'}}\left(\omega \right)$ is the Fourier transform of ${\chi }^{\text{'}\text{'}}\left(t\right)$. The preceding equation is called a
spectral representation of $\chi \left(z\right)$. It is used to define the retarted / advanced dynamical susceptibilities
 ${\chi }^{R/A}\left(u\right)=\chi \left(u±i{0}^{+}\right)=\left[P{\int }_{-\infty }^{\infty }\frac{d\omega }{\pi }\frac{{\chi }^{\text{'}\text{'}}\left(\omega \right)}{\omega -u}\right]±i{\chi }^{\text{'}\text{'}}\left(u\right)\equiv {\chi }^{\text{'}}\left(u\right)±i{\chi }^{\text{'}\text{'}}\left(u\right) ,$ $\left(26\right)$
where we have used $1/\left(x±i{0}^{+}\right)=P 1/x±i\delta \left(x\right)$, moreover we have defined the real part ${\chi }^{\text{'}}\left(u\right)$. Eqn. (26) is called the
Kramer-Kronig transformation(representation) of the retarded / advanced susceptibility. Note that while in the derivation of eqn. (25) we needed that $\mathrm{Im}\left(z\right)>0$, we may nevertheless 'use' eqn. (25) also as a definition in for the case of $\mathrm{Im}\left(z\right)<0$. That is how we mathematically extend its usage in the definition of ${\chi }^{A}\left(u\right)$.
With eqn. (16), eqn. (14) is of convolution form and suggest a simple Fourier transform. However some care has to be exercised because of the $\Theta$-function which renders the direct Fourier transform of $\chi \left(t\right)$ delicate. Consider first
 $\begin{array}{cccc}\multicolumn{3}{c}{{\int }_{-\infty }^{\infty }\mathrm{dt} {e}^{i\omega t}2i\Theta \left(t\right){\chi }^{\text{'}\text{'}}\left(t\right)={\int }_{-\infty }^{\infty }{e}^{i\omega t}{\int }_{-\infty }^{\infty }\frac{d\omega \text{'}}{\pi }\frac{{e}^{i\omega \text{'}t}}{\omega \text{'}-i{0}^{+}}{\chi }^{\text{'}\text{'}}\left(t\right) \mathrm{dt}}& \hfill \left(27\right)\\ \multicolumn{1}{c}{}& =\hfill & {\int }_{-\infty }^{\infty }\frac{d\omega \text{'}}{\pi }\frac{1}{\omega \text{'}-i{0}^{+}}{\int }_{-\infty }^{\infty }{e}^{i\left(\omega +\omega \text{'}\right)t}{\chi }^{\text{'}\text{'}}\left(t\right) \mathrm{dt}={\int }_{-\infty }^{\infty }\frac{d\omega \text{'}}{\pi }\frac{{\chi }^{\text{'}\text{'}}\left(\omega \text{'}\right)}{\omega \text{'}-\left(\omega +i{0}^{+}\right)}={\chi }^{R}\left(\omega \right) .\hfill \end{array}$

I.e., the because the $\Theta$-function is a distribution (like the $\delta$-function), the Fourier transform of $\chi \left(t\right)$ is
not simply $\chi \left(\omega \right)$ but ${\chi }^{R}\left(\omega \right)$. Using this, eqn. (14,16), and the convolution-theorem, and re-introducing the subscripts we get
 $\delta ⟨{A}_{\mu }^{+}\left(\omega \right)⟩=\sum _{\nu }{\chi }_{\mu \nu }^{R}\left(\omega \right){f}_{\nu }\left(\omega \right) .$ $\left(28\right)$
Apart from the case of a driving force at frequency $\omega$, as in eqn. (28), another important case is that of a perturbation switched on adiabatically i.e. infinitesimally slow up to $t=0$
 ${f}_{\nu }\left(t\right)=\Theta \left(-t\right){e}^{{0}^{+}t}{f}_{\nu } .$ $\left(29\right)$
In that case, one is interested in the response at $t=0$
 $\begin{array}{ccc}\multicolumn{1}{c}{\delta ⟨{A}_{\mu }^{+}\left(t=0\right)⟩}& =\hfill & \sum _{\nu }{\int }_{-\infty }^{0}{\chi }_{\mu \nu }\left(0-t\text{'}\right){e}^{{0}^{+}t\text{'}}{f}_{\nu }\mathrm{dt}\text{'}=\sum _{\nu }{\int }_{0}^{\infty }{\chi }_{\mu \nu }\left(t\text{'}\right){e}^{-{0}^{+}t\text{'}}\mathrm{dt}\text{'} {f}_{\nu }\hfill \\ \multicolumn{1}{c}{}& =\hfill & \sum _{\nu }{\chi }_{\mu \nu }^{R}\left(\omega =0\right) {f}_{\nu }\equiv \sum _{\nu }{\chi }_{\mu \nu }^{\mathrm{Iso}} {f}_{\nu }\hfill & \hfill \left(30\right)\end{array}$

${\chi }_{\mu \nu }^{\mathrm{Iso}}$ is called the
isolated/adiabatic susceptibility. It is very tempting to associate this susceptibility with the isothermal susceptibility
 ${\chi }_{\mu \nu }^{T}=\frac{\partial ⟨{A}_{\mu }^{+}⟩}{\partial {f}_{\nu }} .$ $\left(31\right)$
from static thermodynamics. However this is wrong, as ${\chi }_{\mu \nu }^{\mathrm{Iso}}\ne {\chi }_{\mu \nu }^{T}$ in general, as we will see later.
Example 2 Let us apply the preceding to the results from example 1. From eqn. (22) we get
 ${\chi }_{\mathrm{rr}}^{\text{'}\text{'}}\left(t\right)=\frac{1}{2}\left({e}^{-i{\omega }_{0}t}-{e}^{i{\omega }_{0}t}\right) ,$ $\left(32\right)$
and from that
 $\begin{array}{ccc}\multicolumn{1}{c}{{\chi }_{\mathrm{rr}}^{\text{'}\text{'}}\left(\omega \right)}& =\hfill & \frac{1}{2}{\int }_{-\infty }^{\infty }{e}^{i\omega t}\left({e}^{-i{\omega }_{0}t}-{e}^{i{\omega }_{0}t}\right)\mathrm{dt}\hfill \\ \multicolumn{1}{c}{}& =\hfill & \pi \left[\delta \left(\omega -{\omega }_{0}\right)-\delta \left(\omega +{\omega }_{0}\right)\right] ,\hfill & \hfill \left(33\right)\end{array}$

i.e., the Fourier transform of the spectral function has 'peaks' at the excitation energy of the system. Therefore it is frequently called 'the spectrum' of the system. The retarded / advanced susceptibility follows from the spectral representation eqn. (25)
 $\begin{array}{ccc}\multicolumn{1}{c}{{\chi }_{\mathrm{rr}}^{R/A}\left(\omega \right)}& =\hfill & {\int }_{-\infty }^{\infty }\frac{\delta \left(\omega -{\omega }_{0}\right)-\delta \left(\omega +{\omega }_{0}\right)}{\omega \text{'}-\omega \mp i{0}^{+}}d\omega \text{'}=\hfill \\ \multicolumn{1}{c}{}& =\hfill & \frac{2{\omega }_{0}}{{\omega }_{0}^{2}-\left(\omega \mp i{0}^{+}{\right)}^{2}} ,\hfill & \hfill \left(34\right)\end{array}$

from which we can deduce the isolated susceptibility
 ${\chi }_{\mathrm{rr}}^{\mathrm{Iso}}=\frac{2}{{\omega }_{0}} .$ $\left(35\right)$
For ${\omega }_{0}\to 0$ the system gets 'soft', i.e. it is physically clear, that in this situation its response to an adiabatic perturbation will diverge - as in eqn. (35).
Exercise 3 Assume ${\chi }^{\text{'}\text{'}}\left(\omega \right)=2\sqrt{1-{\omega }^{2}}·\Theta \left(1-{\omega }^{2}\right)$. Calculate $\chi \left(z\right)$ from eqn. (25). Verify, that eqn. (26) is satisfied, regarding ${\chi }^{\text{'}\text{'}}\left(\omega \right)$. Plot ${\chi }^{\text{'}\left(\text{'}\text{'}\right)}\left(\omega \right)$. Discuss the large- $\omega$ behavior of ${\chi }^{R}\left(\omega \right)$.

### 5.4Symmetries

The symmetries of ${\chi }_{\mu \nu }^{\left(\text{'}\text{'}\right)}\left(\omega \right)$ are of great importance. For eg., they serve as consistency checks, both for experiments and theories, they restrict the allowed forms that susceptibilities can take on, and they lead to reciprocity relations between various types of responses, just to name a few. To simplify the notation we define
 ${A}_{\mu }^{+}={A}_{\stackrel{‾}{\mu }}$ $\left(36\right)$
Complex conjugation
 ${\chi }_{\mu \nu }^{\text{'}\text{'}}\left(\omega {\right)}^{☆}=-{\chi }_{\stackrel{‾}{\mu }\stackrel{‾}{\nu }}^{\text{'}\text{'}}\left(-\omega \right) .$ $\left(37\right)$
Proof
 $\begin{array}{cccc}\multicolumn{1}{c}{{\chi }_{\mu \nu }^{\text{'}\text{'}}\left(\omega {\right)}^{☆}}& =\hfill & {\left[{\int }_{-\infty }^{\infty }{e}^{i\omega t}\frac{1}{2}⟨\left[{A}_{\mu }^{+}\left(t\right),{A}_{\nu }\right]⟩\mathrm{dt}\right]}^{☆}={\int }_{-\infty }^{\infty }{e}^{-i\omega t}\frac{1}{2}⟨\left[{A}_{\mu }^{+}\left(t\right),{A}_{\nu }\right]{⟩}^{☆}\mathrm{dt}\hfill & \hfill \left(38\right)\\ \multicolumn{1}{c}{}& =\hfill & {\int }_{-\infty }^{\infty }{e}^{-i\omega t}\frac{1}{2}⟨\left[{A}_{\nu }^{+},{A}_{\mu }^{++}\left(t\right)\right]⟩\mathrm{dt}=-{\int }_{-\infty }^{\infty }{e}^{-i\omega t}\frac{1}{2}⟨\left[{A}_{\mu }^{++}\left(t\right),{A}_{\nu }^{+}\right]⟩\mathrm{dt}=-{\chi }_{\stackrel{‾}{\mu }\stackrel{‾}{\nu }}^{\text{'}\text{'}}\left(-\omega \right)\hfill \end{array}$

where we have used $⟨A{⟩}^{☆}=⟨{A}^{+}⟩$ and $\left({A}_{\mu }^{+}\left(t{\right)\right)}^{+}={A}_{\mu }^{++}\left(t\right)$.
Commutator
 ${\chi }_{\mu \nu }^{\text{'}\text{'}}\left(\omega \right)=-{\chi }_{\stackrel{‾}{\nu }\stackrel{‾}{\mu }}^{\text{'}\text{'}}\left(-\omega \right) .$ $\left(39\right)$
Proof
 $\begin{array}{cccc}\multicolumn{1}{c}{{\chi }_{\mu \nu }^{\text{'}\text{'}}\left(\omega \right)}& =\hfill & {\int }_{-\infty }^{\infty }{e}^{i\omega t}\frac{1}{2}⟨\left[{A}_{\mu }^{+}\left(t\right),{A}_{\nu }\right]⟩\mathrm{dt}=-{\int }_{-\infty }^{\infty }{e}^{i\omega t}\frac{1}{2}⟨\left[{A}_{\nu },{A}_{\mu }^{+}\left(t\right)\right]⟩\mathrm{dt}\hfill & \hfill \left(40\right)\\ \multicolumn{1}{c}{}& =\hfill & -{\int }_{-\infty }^{\infty }{e}^{i\omega t}\frac{1}{2}⟨\left[{A}_{\nu }^{++}\left(-t\right),{A}_{\mu }^{+}\right]⟩\mathrm{dt}=-{\int }_{-\infty }^{\infty }{e}^{-i\omega t}\frac{1}{2}⟨\left[{A}_{\nu }^{++}\left(t\right),{A}_{\mu }^{+}\right]⟩\mathrm{dt}=-{\chi }_{\stackrel{‾}{\nu }\stackrel{‾}{\mu }}^{\text{'}\text{'}}\left(-\omega \right) .\hfill \end{array}$

For ${A}_{\mu }={A}_{\nu }\equiv A$ and ${A}^{+}=A$ we get ${\chi }^{\text{'}\text{'}}\left(\omega \right)=-{\chi }^{\text{'}\text{'}}\left(-\omega \right)$, i.e. in that case the spectral function is
real and odd in frequency.
Time inversion is a linear operation which needs some introduction, which for completeness will be included here. In its coordinate representation any solution $\psi \left(r,t\right)$ of the one-particle Schrödinger equation has the following symmetry
 $i{\partial }_{t}\psi \left(r,t\right)=\left[\frac{{p}^{2}}{2m}+V\left(r\right)\right]\psi \left(r,t\right)\mathrm{ }⇔\mathrm{ }i{\partial }_{t}{\psi }^{☆}\left(r,-t\right)=\left[\frac{{p}^{2}}{2m}+V\left(r\right)\right]{\psi }^{☆}\left(r,-t\right) ,$ $\left(41\right)$
which can easily be generalized to many particles. I.e. for states "inversion of time + complex conjugation" is a symmetry of the Schrödinger equation. Under time inversion the operators 'space' $r$ and 'momentum' $p$ change according to
 $\Pi r{\Pi }^{-1}=r\mathrm{ }\Pi p{\Pi }^{-1}=-p .$ $\left(42\right)$
For $\left[r,p\right]$ to be invariant under $\Pi$, the following has to hold
 $\begin{array}{ccc}\multicolumn{1}{c}{-i}& =\hfill & \left[r,-p\right]=\left[\Pi r{\Pi }^{-1},\Pi p{\Pi }^{-1}\right]=\Pi \left[r,p\right]{\Pi }^{-1}=\Pi i{\Pi }^{-1}\hfill \\ \multicolumn{1}{c}{⇒-i\Pi }& =\hfill & \Pi i ,\hfill & \hfill \left(43\right)\end{array}$

which is generalized into the requirement
 ${a}^{☆}\Pi =\Pi a ,$ $\left(44\right)$
i.e. time inversion is
anti-linear.
While the adjoint of any linear operator is given by
 $⟨\mu |A\nu ⟩=⟨\nu |{A}^{+}\mu {⟩}^{☆}=⟨{A}^{+}\mu |\nu ⟩ ,$ $\left(45\right)$
(where we have used von-Neumann's notation for scalar products instead of Dirac's) it will turn out to be reasonable to
 $⟨{\Pi }^{+}\mu |\nu {⟩}^{☆}\equiv ⟨\mu |\Pi \nu ⟩ .$ $\left(46\right)$
From this and the requirement that the norm is conserved under $\Pi$ - which is
consistent eqn. (45) we get
 $\begin{array}{cccc}\multicolumn{1}{c}{⟨\Pi \mu |\Pi \mu ⟩}& =\hfill & ⟨{\Pi }^{+}\Pi \mu |\mu {⟩}^{☆}\equiv ⟨\mu |\mu {⟩}^{☆}=⟨\mu |\mu ⟩\hfill & \hfill \left(47\right)\\ \multicolumn{1}{c}{⇒{\Pi }^{+}\Pi =1}& \hfill & {\Pi }^{-1}\Pi =\Pi {\Pi }^{-1}={\Pi }^{+}\Pi =\Pi {\Pi }^{+}=1 .\hfill & \hfill \left(48\right)\end{array}$

It is important to realize that in view of eqn. (41) it is the
2nd equal sign in eqn. (47) which is consistent with eqn. (46) but not with eqn. (45). This is the reason for introducing the anti-adjoint. In turn time inversion is anti-unitary.
Now we assume, that $H$ and $\rho$ are invariant under $\Pi$
 $\left[H,{\Pi }^{\left(-1,+\right)}\right]=\left[\rho ,{\Pi }^{\left(-1,+\right)}\right]=0 .$ $\left(49\right)$
Applying this on the spectral function in the time-domain and dropping the subscripts $\mu ,\nu$ for this
 $\begin{array}{ccc}\multicolumn{3}{c}{\mathrm{Tr}\left\{\rho \left[A\left(t\right),B\left(t\text{'}\right)\right]\right\}=\sum _{n}⟨n|{\Pi }^{+}\Pi \rho \left[A\left(t\right),B\left(t\text{'}\right)\right]n⟩=\sum _{n}⟨\Pi n|\Pi \rho \left[A\left(t\right),B\left(t\text{'}\right)\right]n{⟩}^{☆}}\\ \multicolumn{1}{c}{}& =\hfill & \sum _{n}⟨\Pi n|\Pi \rho \left[A\left(t\right),B\left(t\text{'}\right)\right]{\Pi }^{+}\Pi n{⟩}^{☆}=\sum _{n}⟨\Pi n|\rho \left[\Pi A\left(t\right){\Pi }^{+},\Pi B\left(t\text{'}\right){\Pi }^{+}\right]\Pi n{⟩}^{☆} .\hfill & \hfill \left(50\right)\end{array}$

Since $|n⟩$ in the preceding is a complete orthonormal set (COS), $|\Pi n⟩$ is also and we may replace it by $|n⟩$ in the trace. Many operators have a so-called
parity $\epsilon$ w.r.t. certain operations like time-reversal $\Pi$ or (real-space) inversion $P$, i.e.
 $\begin{array}{cccc}\multicolumn{1}{c}{\Pi A\left(t\right){\Pi }^{-1}}& =\hfill & {\epsilon }^{t}A\left(-t\right)\hfill & \hfill \left(51\right)\\ \multicolumn{1}{c}{\mathrm{PA}\left(t\right){P}^{-1}}& =\hfill & {\epsilon }^{P}A\left(t\right) .\hfill & \hfill \left(52\right)\end{array}$

Typical examples are coordinate $r$ with ${\epsilon }^{P}=-1$ and ${\epsilon }^{t}=1$, or velocity $v$ with ${\epsilon }^{P}=-1$ and ${\epsilon }^{t}=-1$, or (electric/)magnetic field ( $E$/) $B$ with ${\epsilon }^{P}=\left(-1/\right)-1$ and ${\epsilon }^{t}=\left(1/\right)-1$. In the following we restrict ourselves to such. Then eqn. (50) continues as
 $={\epsilon }_{A}^{t}{\epsilon }_{B}^{t}\sum _{n}⟨n|\rho \left[A\left(-t\right),B\left(-t\text{'}\right)\right]n{⟩}^{☆}={\epsilon }_{A}^{t}{\epsilon }_{B}^{t}{\chi }_{\mathrm{AB}}\left(-t,-t\text{'}{\right)}^{☆}=-{\epsilon }_{A}^{t}{\epsilon }_{B}^{t}{\chi }_{\stackrel{‾}{A}\stackrel{‾}{B}}\left(-t,-t\text{'}\right) .$ $\left(53\right)$
In the frequency domain assuming eqn. (16) and returning to the subscripts we get
 ${\chi }_{\mu \nu }\left(\omega \right)=-{\epsilon }_{\mu }^{t}{\epsilon }_{\nu }^{t}{\chi }_{\stackrel{‾}{\mu }\stackrel{‾}{\nu }}\left(-\omega \right) .$ $\left(54\right)$
Example 4 Consider the spectral function from examples 1 and 2
 ${\chi }_{\mathrm{rr}}^{\text{'}\text{'}}\left(\omega \right)=\pi \left[\delta \left(\omega -{\omega }_{0}\right)-\delta \left(\omega +{\omega }_{0}\right)\right] .$ $\left(55\right)$
Within the notation of eqn. (36) we have $\stackrel{‾}{r}={r}^{+}=r$. From this

The commutator rule
 ${\chi }_{\mathrm{rr}}^{\text{'}\text{'}}\left(\omega \right)=-{\chi }_{\stackrel{‾}{r}\stackrel{‾}{r}}^{\text{'}\text{'}}\left(-\omega \right)=-{\chi }_{\mathrm{rr}}^{\text{'}\text{'}}\left(-\omega \right)$ $\left(56\right)$
is obviously satisfied. The
complex conjugation rule
 ${\chi }_{\mathrm{rr}}^{\text{'}\text{'}}\left(\omega {\right)}^{☆}={\chi }_{\mathrm{rr}}^{\text{'}\text{'}}\left(\omega {\right)}^{☆}=-{\chi }_{\mathrm{rr}}^{\text{'}\text{'}}\left(-\omega \right)$ $\left(57\right)$
sets ${\chi }_{\mathrm{rr}}^{\text{'}\text{'}}\left(\omega \right)$ to be
real and odd in $\omega$. Finally, the t-reversal rule together with ${\epsilon }_{r}^{t}=+1$
 ${\chi }_{\mathrm{rr}}^{\text{'}\text{'}}\left(\omega \right)=-\left(+1\right)\left(+1\right){\chi }_{\mathrm{rr}}^{\text{'}\text{'}}\left(-\omega \right)$ $\left(58\right)$

Symmetries allow practically relevant conclusions about the vanishing of certain susceptibilities. Consider eg. the isolated (momentum $p$,real-space $r$)-susceptibility
 ${\chi }_{\mathrm{pr}}^{\mathrm{Iso}} .$ $\left(59\right)$
From t-reversal we get
 ${\chi }_{\mathrm{pr}}^{\text{'}\text{'}}\left(\omega \right)=-\left(-1\right)\left(+1\right){\chi }_{\mathrm{pr}}^{\text{'}\text{'}}\left(-\omega \right)={\chi }_{\mathrm{pr}}^{\text{'}\text{'}}\left(-\omega \right) .$ $\left(60\right)$
Contrast this against eqn. (58) ! I.e. ${\chi }_{\mathrm{pr}}^{\text{'}\text{'}}\left(\omega \right)$ is even in $\omega$. Using eqn. (27,30) we get
 ${\chi }_{\mathrm{pr}}^{\mathrm{Iso}}={\int }_{-\infty }^{\infty }\frac{d\omega \text{'}}{\pi }\frac{{\chi }_{\mathrm{pr}}^{\text{'}\text{'}}\left(\omega \text{'}\right)}{\omega \text{'}}=0 .$ $\left(61\right)$
A similar argument leads to ${\chi }_{\mathrm{ME}}^{\mathrm{Iso}}=0$ and ${\chi }_{\mathrm{PB}}^{\mathrm{Iso}}=0$, for a time-inversion symmetric $H$ and $\rho$, where $E\left(B\right)$ and $P\left(M\right)$ are the electric(magnetic) field and the electric polarization (magnetization).

At first sight, generating a magnetic moment with an electric field seems counter-intuitive anyway. However, to allow for ${\chi }_{\mathrm{ME}}^{\mathrm{Iso}}\ne 0$ or ${\chi }_{\mathrm{PB}}^{\mathrm{Iso}}\ne 0$ a system simply has to break time-inversion symmetry. By the famous P(arity) C(harge conjugation) T(ime inversion)-theorem, this is synonymous to the breaking of parity (i.e. space inversion symmetry) in most cases. Solids which display such symmetry loss are also called non-centro-symmetric.
Exercise 5 (Real space) Inversion Consider the $r$-dependent operators
 i) density $\stackrel{^}{n}\left(r\right)=\sum _{l}\delta \left(r-{\stackrel{^}{R}}_{l}\right)$ ii) momentum density $\stackrel{^}{p}\left(r\right)=\sum _{l}\delta \left(r-{\stackrel{^}{R}}_{l}\right){\stackrel{^}{p}}_{l}$ iii) spin density $\stackrel{^}{S}\left(r\right)=\sum _{l}\delta \left(r-{\stackrel{^}{R}}_{l}\right){\stackrel{^}{S}}_{l}$
of an $N$-particle system, with $l=1\dots N$, where the ' ' refers to operators and $r$ is a real space coordinate.
a) Use that under inversion $P{\stackrel{^}{R}}_{l}{P}^{-1}=-{\stackrel{^}{R}}_{l}$, $P{\stackrel{^}{p}}_{l}{P}^{-1}=-{\stackrel{^}{p}}_{l}$, and $P{\stackrel{^}{S}}_{l}{P}^{-1}={\stackrel{^}{S}}_{l}$ to show that
 $\mathrm{PA}\left(r\right){P}^{-1}={\epsilon }_{A}^{P}A\left(-r\right) ,$ $\left(62\right)$
for $A\left(r\right)=\stackrel{^}{n}\left(r\right)$, $\stackrel{^}{p}\left(r\right)$ and $\stackrel{^}{S}\left(r\right)$, with ${\epsilon }_{n}^{P}=1$, ${\epsilon }_{p}^{P}=-1$, and ${\epsilon }_{S}^{P}=1$.
b) Assume a statistical operator $\rho$ invariant under inversion $\left[\rho ,H\right]=0$. Defining the the real-space dependent spectral function ${\chi }_{\mu \nu }\left(r,t;r\text{'},t\text{'}\right)\equiv i\Theta \left(t-t\text{'}\right)⟨\left[{A}_{\mu }^{+}\left(r,t\right),{A}_{\nu }\left(r\text{'},t\text{'}\right)\right]⟩$, show that in this case
 ${\chi }_{\mu \nu }\left(r,t;r\text{'},t\text{'}\right)={\epsilon }_{\mu }^{P}{\epsilon }_{\nu }^{P}{\chi }_{\mu \nu }\left(-r,t;-r\text{'},t\text{'}\right) .$ $\left(63\right)$
As a consequence, assuming translational invariance and a time independent Hamilton operator, i.e. ${\chi }_{\mu \nu }\left(r,t;r\text{'},t\text{'}\right)={\chi }_{\mu \nu }\left(r-r\text{'},t-t\text{'}\right)$, show that the Fourier transform into wave vector (i.e. $q$) and frequency (i.e. $\omega$) space satisfies
 ${\chi }_{\mu \nu }\left(q,\omega \right)={\epsilon }_{\mu }^{P}{\epsilon }_{\nu }^{P}{\chi }_{\mu \nu }\left(-q,\omega \right) .$ $\left(64\right)$
c) Use eqn. (30) to and eqn. (64) to derive the isolated susceptibilities ${\chi }_{\mathrm{np}}^{\mathrm{Iso}}\left(r=0\right)$ and ${\chi }_{\mathrm{sp}}^{\mathrm{Iso}}\left(r=0\right)$. What is the physical significance of this.

### 5.5Mori's Scalar Product

We may define the following mapping from the vector space of linear operators ${A}_{\mu }$ on the scalar space of complex numbers into the complex numbers
 $\begin{array}{cccc}\multicolumn{1}{c}{⟨{A}_{\mu }|{A}_{\nu }⟩}& =\hfill & {\int }_{0}^{\beta }⟨\Delta {A}_{\mu }^{+}\left(\tau \right)\Delta {A}_{\nu }⟩d\tau \hfill & \hfill \left(65\right)\\ \multicolumn{1}{c}{\beta }& =\hfill & 1/T\hfill \\ \multicolumn{1}{c}{\Delta A}& =\hfill & A-⟨A⟩\hfill & \hfill \left(66\right)\\ \multicolumn{1}{c}{A\left(\tau \right)}& =\hfill & {e}^{\tau H}{\mathrm{Ae}}^{-\tau H}\mathrm{ }\tau \in \mathrm{Re} .\hfill & \hfill \left(67\right)\end{array}$

Here $\tau$ is called
'imaginary' time. Note, that unlike for real time $\left[A\left(\tau {\right)\right]}^{+}\ne {A}^{+}\left(\tau \right)$. The brackets on the r.h.s $⟨⟩$ are usual thermal averages. The brackets $⟨\dots |\dots ⟩$on the l.h.s. of eqn. (65) should not be confused with a scalar product of quantum mechanical states. Yet, this so-called Mori product satisfies all properties of a conventional scalar product, i.e. it is:
i) Bilinear: this is trivial
ii) Positive semi-definite: (w.l.o.g. set $⟨A⟩=0$)
 $\begin{array}{ccc}\multicolumn{1}{c}{⟨A|A⟩}& =\hfill & {\int }_{0}^{\beta }⟨{A}^{+}\left(\tau \right)A⟩d\tau =\frac{1}{Z}\sum _{\mathrm{lm}}{\int }_{0}^{\beta }{e}^{-\beta {E}_{l}}⟨l|{A}^{+}\left(\tau \right)|m⟩⟨m|A|l⟩d\tau \hfill \\ \multicolumn{1}{c}{}& =\hfill & \frac{1}{Z}\sum _{\mathrm{lm}}{e}^{-\beta {E}_{l}}⟨l|{A}^{+}|m⟩⟨m|A|l⟩{\int }_{0}^{\beta }{e}^{\tau \left({E}_{l}-{E}_{m}\right)}d\tau \hfill \\ \multicolumn{1}{c}{}& =\hfill & \frac{1}{Z}\sum _{{E}_{m}\ne {E}_{n}}|⟨m|A|l⟩{|}^{2}\frac{{e}^{-\beta {E}_{m}}-{e}^{-\beta {E}_{l}}}{{E}_{l}-{E}_{m}}+\frac{\beta }{Z}\sum _{{E}_{m}={E}_{n}}|⟨m|A|l⟩{|}^{2} .\hfill & \hfill \left(68\right)\end{array}$

Both sums on the last line contain only positive addends. Moreover if, and only if all matrix elements of $A$ are zero, i.e. $A=0$, then $⟨A|A⟩=0$. q.e.d.

iii) Complex conjugation:
 $\begin{array}{ccc}\multicolumn{1}{c}{⟨A|B{⟩}^{☆}}& =\hfill & {\int }_{0}^{\beta }⟨{A}^{+}\left(\tau \right)B{⟩}^{☆}d\tau ={\int }_{0}^{\beta }\mathrm{Tr}\left\{\rho {A}^{+}\left(\tau \right)B{\right\}}^{☆}d\tau ={\int }_{0}^{\beta }\mathrm{Tr}\left\{{B}^{+}\left[{A}^{+}\left(\tau {\right)\right]}^{+}\rho \right\}d\tau \hfill \\ \multicolumn{1}{c}{}& =\hfill & {\int }_{0}^{\beta }\mathrm{Tr}\left\{{B}^{+}A\left(-\tau \right)\rho \right\}d\tau ={\int }_{0}^{\beta }\mathrm{Tr}\left\{{B}^{+}\left(\tau \right)A\rho \right\}d\tau =⟨B|A⟩ .\hfill & \hfill \left(69\right)\end{array}$

Mori's product has many applications in statistical mechanics. The first one applies to the isothermal susceptibility
 ${\chi }_{\mu \nu }^{T}=\frac{\partial ⟨{A}_{\mu }^{+}⟩}{\partial {f}_{\nu }}=⟨{A}_{\mu }|{A}_{\nu }⟩ .$ $\left(70\right)$
The high-temperature limit $\beta \to 0$ and also the classical case, i.e. $\left[H,{A}_{\mu }\right]=0$, of this equation can be read of easily from eqn. (65)
 ${\chi }_{\mu \nu }^{T} \genfrac{}{}{0}{}{\to }{\genfrac{}{}{0}{}{\beta \to 0 \mathrm{or}}{\left[H,{A}_{\mu }\right]=0}} \frac{⟨\Delta {A}_{\mu }^{+}\Delta {A}_{\nu }⟩}{T} .$ $\left(71\right)$
This is the famous relation proportionality between fluctuations and susceptibilities. Eg. for a magnetic system in this limit we have ${\chi }^{T}=⟨\left(\Delta S{\right)}^{2}⟩/T$, which is simply
Curie's law. The general proof of eqn. (70) is left to exercise 8.
Exercise 6 Given two operators with ${A}_{\mu }^{+}={A}_{\mu }$ and ${A}_{\nu }^{+}={A}_{\nu }$. Show that $⟨{A}_{\mu }|{A}_{\nu }⟩=⟨{A}_{\nu }|{A}_{\mu }⟩$. I.e. for hermitean operators, and because of eqn. (69), Mori's product is real.
Example 7 Let's continue with the harmonic oscillator. Analogous to eqn. (19) we get from eqn. (67)
 $\begin{array}{ccc}\multicolumn{1}{c}{b\left(\tau \right)}& =\hfill & {e}^{-{\omega }_{0}\tau }b\hfill \\ \multicolumn{1}{c}{{b}^{+}\left(\tau \right)}& =\hfill & {e}^{{\omega }_{0}\tau }{b}^{+} .\hfill & \hfill \left(72\right)\end{array}$

Using that $⟨b⟩=⟨{b}^{+}⟩=0$, and inserting into eqn. (65) we get
 $\begin{array}{ccc}\multicolumn{1}{c}{⟨r|r⟩}& =\hfill & {\int }_{0}^{\beta }⟨\left({e}^{{\omega }_{0}\tau }{b}^{+}+{e}^{-{\omega }_{0}\tau }b\right)\left({b}^{+}+b\right)⟩d\tau \hfill \\ \multicolumn{1}{c}{}& =\hfill & {\int }_{0}^{\beta }{e}^{{\omega }_{0}\tau }⟨{b}^{+}b⟩d\tau +{\int }_{0}^{\beta }{e}^{-{\omega }_{0}\tau }⟨{\mathrm{bb}}^{+}⟩d\tau \hfill \\ \multicolumn{1}{c}{}& =\hfill & n\left(\beta \right)\left(\frac{1}{{\omega }_{0}}\right)\left({e}^{{\omega }_{0}\beta }-1\right)+\left(1+n\left(\beta \right)\right)\left(-\frac{1}{{\omega }_{0}}\right)\left({e}^{-{\omega }_{0}\beta }-1\right)\hfill \\ \multicolumn{1}{c}{}& =\hfill & \frac{1}{{\omega }_{0}}+\frac{1}{{\omega }_{0}}=\frac{2}{{\omega }_{0}}\hfill & \hfill \left(73\right)\end{array}$

where we have used $⟨\mathrm{bb}⟩=⟨{b}^{+}{b}^{+}⟩=0$, and $⟨{b}^{+}b⟩=n\left(\beta \right)=1/\left({e}^{\beta {\omega }_{0}}-1\right)$, and $\left[b,{b}^{+}\right]=1$.

From this and eqns. (35,71) we can read off that
 ${\chi }_{\mathrm{rr}}^{\mathrm{Iso}}={\chi }_{\mathrm{rr}}^{T}$ $\left(74\right)$
Exercise 8 Thermodynamic perturbation theory. Prove eqn. (70)
a) TODO
b)
c)

### 5.6The Liouville Operator

The time-evolution resulting from the equations of motion for operators in the Heisenberg picture can be brought into a form formally similar to that of states in the Schrödinger picture by introducing a so-called super operator, the Liouville operator $L$, which by definition is an operator acting on operators
 $\begin{array}{cccc}\multicolumn{1}{c}{{\mathrm{LA}}_{\mu }\left(t\right)}& \equiv \hfill & \left[H,{A}_{\mu }\left(t\right)\right]\hfill & \hfill \left(75\right)\\ \multicolumn{1}{c}{⇒\frac{d}{\mathrm{dt}}{A}_{\mu }\left(t\right)}& =\hfill & {\mathrm{iLA}}_{\mu }\left(t\right)\hfill & \hfill \left(76\right)\\ \multicolumn{1}{c}{{A}_{\mu }\left(t\right)}& =\hfill & {e}^{\mathrm{iLt}}{A}_{\mu } .\hfill & \hfill \left(77\right)\end{array}$

The Liouville operator is hermitean w.r.t. to the Mori product (w.l.o.g. $⟨A⟩=⟨B⟩=0$)
 $\begin{array}{ccc}\multicolumn{1}{c}{⟨A|\mathrm{LB}⟩}& =\hfill & {\int }_{0}^{\beta }⟨{A}^{+}\left(\tau \right)\mathrm{LB}⟩d\tau ={\int }_{0}^{\beta }⟨{e}^{\tau H}{A}^{+}{e}^{-\tau H}\left[H,B\right]⟩d\tau ={\int }_{0}^{\beta }⟨{e}^{\tau H}\left[H,A{\right]}^{+}{e}^{-\tau H}B⟩d\tau \hfill \\ \multicolumn{1}{c}{}& =\hfill & {\int }_{0}^{\beta }⟨{e}^{-\tau H}{\mathrm{Be}}^{\tau H}\left[H,A{\right]}^{+}⟩d\tau ={\int }_{0}^{\beta }⟨{B}^{+}\left(\tau \right)\left[H,A\right]{⟩}^{☆}d\tau =⟨B|\mathrm{LA}{⟩}^{☆}\hfill & \hfill \left(78\right)\end{array}$

## Index (showing section)

 anti-adjoint, 5.4 anti-linear, 5.4 anti-unitary, 5.4 canonical distribution function, 5.1 Curie's law, 5.5 dynamical susceptibility, 5.2 imaginary time, 5.5 inversion, 5.4 isolated/adiabatic susceptibility, 5.3 isothermal susceptibility, 5.3 Kramer-Kronig, 5.3 microstates, 5.1 Mori product, 5.5 non-centro-symmetric, 5.4 on Neumann's equation, 5.1 parity, 5.4 PCT-theorem, 5.4 retarded/advanced, 5.3 spectral function, 5.2 spectral representation, 5.3 super operator, 5.6 time inversion, 5.4

File translated from TEX by TTM, version 3.85.
On 24 Jan 2010, 21:36.