8.4: Distribution Functions
- Page ID
- 6381
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)
( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\id}{\mathrm{id}}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\kernel}{\mathrm{null}\,}\)
\( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\)
\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\)
\( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)
\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)
\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)
\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vectorC}[1]{\textbf{#1}} \)
\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)
\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)
\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)8.4.1 One-particle distribution functions
What is the mean number of particles in the box of volume d3rA about rA?
The probability that particle 1 is in d3rA about rA is
\[ \frac{d^{3} r_{A} \int d^{3} r_{2} \int d^{3} r_{3} \cdots \int d^{3} r_{N} \int d^{3} p_{1} \cdots \int d^{3} p_{N} e^{-\beta H\left(\mathbf{r}_{A}, \mathbf{r}_{2}, \mathbf{r}_{3}, \ldots, \mathbf{r}_{N}, \mathbf{p}_{1}, \ldots, \mathbf{p}_{N}\right)}}{\int d^{3} r_{1} \int d^{3} r_{2} \int d^{3} r_{3} \cdots \int d^{3} r_{N} \int d^{3} p_{1} \cdots \int d^{3} p_{N} e^{-\beta H\left(\mathbf{r}_{1}, \mathbf{r}_{2}, \mathbf{r}_{3}, \ldots, \mathbf{r}_{N}, \mathbf{p}_{1}, \ldots, \mathbf{p}_{N}\right)}}.\]
The probability that particle 2 is in d3rA about rA is
\[ \frac{\int d^{3} r_{1} \quad d^{3} r_{A} \int d^{3} r_{3} \cdots \int d^{3} r_{N} \int d^{3} p_{1} \cdots \int d^{3} p_{N} e^{-\beta H\left(\mathbf{r}_{1}, \mathbf{r}_{A}, r_{3}, \ldots, \mathbf{r}_{N}, \mathbf{p}_{1}, \ldots, \mathbf{p}_{N}\right)}}{\int d^{3} r_{1} \int d^{3} r_{2} \int d^{3} r_{3} \cdots \int d^{3} r_{N} \int d^{3} p_{1} \cdots \int d^{3} p_{N} e^{-\beta H\left(\mathbf{r}_{1}, \mathbf{r}_{2}, r_{3}, \ldots, \mathbf{r}_{N}, \mathbf{p}_{1}, \ldots, \mathbf{p}_{N}\right)}}.\]
And so forth. I could write down N different integrals, but all of them would be equal.
Thus the mean number of particles in d3rA about rA is
\( \pi_{1}\left(\mathbf{r}_{A}\right) d^{3} \boldsymbol{T}_{A}\)
\[ =N \frac{d^{3} r_{A} \int d^{3} r_{2} \int d^{3} r_{3} \cdots \int d^{3} r_{N} \int d^{3} p_{1} \cdots \int d^{3} p_{N} e^{-\beta H\left(\mathbf{r}_{A}, \mathbf{r}_{2}, \mathbf{r}_{3}, \ldots, \mathbf{r}_{N}, \mathbf{p}_{1}, \ldots, \mathbf{p}_{N}\right)}}{\int d^{3} r_{1} \int d^{3} r_{2} \int d^{3} r_{3} \cdots \int d^{3} r_{N} \int d^{3} p_{1} \cdots \int d^{3} p_{N} e^{-\beta H\left(\mathbf{r}_{1}, \mathbf{r}_{2}, \mathbf{r}_{3}, \ldots, \mathbf{r}_{N}, \mathbf{p}_{1}, \ldots, \mathbf{p}_{N}\right)}}\]
\[ =N \frac{d^{3} r_{A} \int d^{3} r_{2} \int d^{3} r_{3} \cdots \int d^{3} r_{N} \int d^{3} p_{1} \cdots \int d^{3} p_{N} e^{-\beta H\left(\mathbf{r}_{A}, r_{2}, r_{3}, \ldots, r_{N}, \mathbf{p}_{1}, \ldots, \mathbf{p}_{N}\right)}}{h^{3 N} N ! Z(T, V, N)}\]
\[ =\frac{1}{(N-1) !} \frac{d^{3} r_{A} \int d^{3} r_{2} \int d^{3} r_{3} \cdots \int d^{3} r_{N} e^{-\beta U\left(r_{A}, r_{2}, r_{3}, \ldots, r_{N}\right)}}{Q(T, V, N)}\]
8.4.2 Two-particle distribution functions
What is the mean number of pairs of particles, such that one member of the pair in a box of volume d3rA about rA and the other member is in a box of volume d3rB about rB?
The probability that particle 1 is in d3rA about rA and particle 2 is in d3rB about rB is
\[ \frac{d^{3} r_{A}}{\int d^{3} r_{1} \int d^{3} r_{2} \int d^{3} r_{3} \cdots \int d^{3} r_{N} \int d^{3} p_{1} \cdots \int d^{3} p_{N} e^{-\beta H\left(\mathbf{r}_{A}, \mathbf{r}_{B}, \mathbf{r}_{3}, \ldots, \mathbf{r}_{N}, \mathbf{p}_{1}, \ldots, \mathbf{p}_{N}\right)}}.\]
The probability that particle 2 is in d3rA about rA and particle 1 is in d3rB about rB is
\[ \frac{d^{3} r_{B} d^{3} r_{A} \int d^{3} r_{3} \cdots \int d^{3} r_{N} \int d^{3} p_{1} \cdots \int d^{3} p_{N} e^{-\beta H\left(\mathbf{r}_{B}, r_{A}, r_{3}, \ldots, r_{N}, \mathbf{p}_{1}, \ldots, \mathbf{p}_{N}\right)}}{\int d^{3} r_{1} \int d^{3} r_{2} \int d^{3} r_{3} \cdots \int d^{3} r_{N} \int d^{3} p_{1} \cdots \int d^{3} p_{N} e^{-\beta H\left(\mathbf{r}_{1}, \mathbf{r}_{2}, r_{3}, \ldots, \mathbf{r}_{N}, \mathbf{p}_{1}, \ldots, \mathbf{p}_{N}\right)}}.\]
The probability that particle 3 is in 3rA about rA and particle 1 is in d3rB about rB is
\[ \frac{d^{3} r_{B} \int d^{3} r_{2} d^{3} r_{A} \cdots \int d^{3} r_{N} \int d^{3} p_{1} \cdots \int d^{3} p_{N} e^{-\beta H\left(\mathbf{r}_{B}, \mathbf{r}_{2}, r_{A}, \ldots, \mathbf{r}_{N}, \mathbf{p}_{1}, \ldots, \mathbf{p}_{N}\right)}}{\int d^{3} r_{1} \int d^{3} r_{2} \int d^{3} r_{3} \cdots \int d^{3} r_{N} \int d^{3} p_{1} \cdots \int d^{3} p_{N} e^{-\beta H\left(\mathbf{r}_{1}, \mathbf{r}_{2}, r_{3}, \ldots, \mathbf{r}_{N}, \mathbf{p}_{1}, \ldots, \mathbf{p}_{N}\right)}}.\]
And so forth. I could write down N(N − 1) different integrals, but all of them would be equal.
Thus the mean number of pairs with one particle in d3rA about rA and the other in d3rB about rB is
\( n_{2}\left(\mathbf{r}_{A}, \mathbf{r}_{B}\right) d^{3} r_{A} d^{3} r_{B}\)
\[ =N(N-1) \frac{d^{3} r_{A} d^{3} r_{B} \int d^{3} r_{3} \cdots \int d^{3} r_{N} \int d^{3} p_{1} \cdots \int d^{3} p_{N} e^{-\beta H\left(\mathbf{r}_{A}, \mathbf{r}_{B}, \mathbf{r}_{3}, \ldots, \mathbf{r}_{N}, \mathbf{p}_{1}, \ldots, \mathbf{p}_{N}\right)}}{\int d^{3} r_{1} \int d^{3} r_{2} \int d^{3} r_{3} \cdots \int d^{3} r_{N} \int d^{3} p_{1} \cdots \int d^{3} p_{N} e^{-\beta H\left(\mathbf{r}_{1}, \mathbf{r}_{2}, \mathbf{r}_{3}, \ldots, \mathbf{r}_{N}, \mathbf{p}_{1}, \ldots, \mathbf{p}_{N}\right)}}\]
\[ =\frac{1}{(N-2) !} \frac{d^{3} r_{A} d^{3} r_{B} \int d^{3} r_{3} \cdots \int d^{3} r_{N} e^{-\beta U\left(\mathbf{r}_{A}, \mathbf{r}_{B}, \mathbf{r}_{3}, \ldots, \mathbf{r}_{N}\right)}}{Q(T, V, N)}.\]
Problems
8.7 Correlations of nearby particles
Suppose that (as is usual) at small distances the interatomic potential u(r) is highly repulsive. Argue that at small r,
\[ g_{2}(r) \approx \text { constant } e^{-u(r) / k_{B} T}.\]
Do not write down a long or elaborate derivation. . . I’m looking for a simple qualitative argument.
8.8 Correlations between non-interacting identical quantal particles
Guess the form of the pair correlation function g2(r) for ideal (non-interacting) fermions and bosons. Sketch your conjectures, and then compare them to the graphs presented by G. Baym in Lectures on Quantum Mechanics (W.A. Benjamin, Inc., Reading, Mass., 1969) pages 428 and 431.
8.9 Correlation functions and structure factors
A typical isotropic fluid, at temperatures above the critical temperature, has correlation functions that are complicated at short distances, but that fall off exponentially at long distances. In fact, the long-distance behavior is
\[ g_{2}(r)=1+\frac{A e^{-r / \xi}}{r}\]
where ξ, the so-called correlation length, depends on temperature and density. In contrast, at the critical temperature the correlation function falls off much more slowly, as
\[ g_{2}(r)=1+\frac{A}{r^{1+\eta}}\]
Find the structure factor
\[ S(\mathbf{k})=\int d^{3} r\left[g_{2}(\mathbf{r})-1\right] e^{-i \mathbf{k} \cdot \mathbf{r}}\]
associated with each of these correlation functions. Will your results match those of experiments at small values of k or at large values (i.e. at long or short wavelengths)?
8.10 Long wavelength structure factor
Show that, for an isotropic fluid, \( \frac{d S(k)}{d k}\) vanishes at k = 0. Here S(k) is the structure factor
\[ S(\mathbf{k})=1+\rho \int d^{3} r\left[g_{2}(\mathbf{r})-1\right] e^{i \mathbf{k} \cdot \mathbf{r}}.\]
8.11 Correlations in a magnetic system
In the Ising model for a magnet, described in problem 4.9, the (net) correlation function is defined by
\[ G_{i} \equiv\left\langle s_{0} s_{i}\right\rangle-\left\langle s_{0}\right\rangle^{2},\]
where the site j = 0 is some arbitrary “central spin”. Using the results of problem 4.9, show that for a lattice of N sites,
\[ \chi_{T}(T, H)=N \frac{m^{2}}{k_{B} T} \sum_{i} G_{i}.\]