Case 1: Functions only of momentum

Suppose that $A=A(P)$, i.e., a function of the momentum operator. Then, the trace can still be evaluated in the coordinate basis:

$$\langle A \rangle = {1 \over Q}\int dx \langle x\vert A(P)e^{-\beta H}\vert x\rangle$$

However, $A(P)$ acting to the left does not act on an eigenvector. Let us insert a coordinate space identity $I=\int dx \vert x\rangle \langle x\vert$ between $A$ and $\exp(-\beta H)$:

$$\langle A \rangle = {1 \over Q}\int dx dx' \langle x\vert A(P)\vert x'\rangle \langle x\vert e^{-\beta H}\vert x\rangle$$

Now, we see that the expectation value can be obtained by evaluating all the coordinate space matrix elements of the operator and all the coordinate space matrix elements of the density matrix.

A particularly useful form for the expectation value can be obtained if a momentum space identity is inserted:

$$\langle A \rangle = {1 \over Q}\int dx dx' dp \langle x\vert... ... p\vert x'\rangle \langle x'\vert e^{-\beta H}\vert x\rangle$$

Now, we see that $A(P)$ acts on an eigenstate (at the price of introducing another integral). Thus, we have

$$\langle A \rangle = {1 \over Q}\int dp A(p) \int dx dx' \lan... ... p\vert x' \rangle \langle x'\vert e^{-\beta H}\vert x\rangle$$

Using the fact that $\langle x\vert p\rangle = (1/2\pi\hbar)\exp(ipx/\hbar)$, we find that

$$\langle A \rangle = {1 \over 2\pi\hbar Q}\int dp A(p) \int d... ...e^{ip(x-x')/\hbar} \langle x'\vert e^{-\beta H}\vert x\rangle$$

In the above expression, we introduce the change of variables

$$r = {x+x' \over 2}\;\;\;\;\;\;\;\;\;\;s=x-x'$$

Then

$$\langle A \rangle = {1 \over 2\pi\hbar Q}\int dp A(p) \int d... ...gle r-{s \over 2}\vert e^{-\beta H}\vert r+{s \over 2}\rangle$$

Define a distribution function

$$\rho_{\rm W}(r,p) = {1 \over 2\pi\hbar} \int ds e^{ips/\hbar... ...gle r-{s \over 2}\vert e^{-\beta H}\vert r+{s \over 2}\rangle$$

Then, the expectation value can be written as

$$\langle A \rangle = {1 \over Q}\int dr dp A(p)\rho_{\rm W}(r,p)$$

which looks just like a classical phase space average using the ``phase space'' distribution function $\rho_{\rm W}(r,p)$. The distribution function $\rho_{\rm W}(r,p)$ is known as the Wigner density matrix and it has many interesting features. For one thing, its classical limit is

$$\rho_{\rm W}(r,p) = \exp\left[-\beta \left({p^2 \over 2m} + U(r)\right)\right]$$

which is the true classical phase space distribution function. There are various examples, in which the exact Wigner distribution function is the classical phase space distribution function, in particularly for quadratic Hamiltonians. Despite its compelling appearance, the evaluation of expectation values of functions of momentum are considerably more difficult than functions of position, due to the fact that the entire density matrix is required. However, there are a few quantities of interest, that are functions of momentum, that can be evaluated without resorting to the entire density matrix. These are thermodynamic quantities which will be discussed in the next section.