Statistical mechanics of fluids#

Classical parition function of fluids#

  • For a homogenous fluid particles

\[Z(\beta, N, V) = \frac{1}{N! h^N}\int e^{-\beta H(p,x)}dp^N dx^N\]

  • For a system composed of three types A B and C for instance the partion function would be:

\[Z(\beta, N, V) = \frac{1}{N_A! N_B! N_C! h^N}\int e^{-\beta H(p,x)}dp^N dx^N\]

With \(N=N_A+N_B+N_C\)

General exprssion of probability distribution of fluids in phase space#

  • The probability of a state in a classical fluid system is \(f(x^N, p^N)\)

\[\boxed{f(x^N, p^N) = \frac{e^{-\beta H(x^N p^N)}}{Z}}\]

Classical hamitlonian breaks into functions of momenta and positions respectively#

  • \[H(x^N, p^N) = K(p^N) + U(x^N)\]



  • \[f(x^N, p^N) = \Phi(p^N) P(r^N)\]



  • \[\boxed{\Phi(p^N) = \frac{e^{-\beta K(p^N)}}{\int dp^N e^{-\beta K(p^N)}}}\]



  • \[\boxed{P(x^N) = \frac{e^{-\beta U(x^N)}}{\int dp^N e^{-\beta U(x^N)}}}\]



Maxwell-Boltzman distirbution for single particles#

Momentum distribution further factorizes into sum of single particle kinetic energies:

  • \[K(p^N) = \sum^N_{i=1} \frac{p^2_i}{2m_i}\]



  • \[\Phi(p^N) = \prod^N_{i=1} \phi(p_i)\]



  • \[\boxed{\phi(p_i) = \frac{e^{-\beta \frac{p^2_i}{2m}}}{\int dp e^{-\beta \frac{p^2}{2m}}}}\]

Note that \(p^2_i = p^2_x+p^2_y+p^2_z\) and \(\phi(p_i)\) is what is known as Maxwell-Boltzmann distribution. The classical system obyes the same distribution reglardles of whivh phase it is (solid, liquid, gas).

Configurational partion function#

\[Z(\beta, V, N) = Z_{p} \cdot Z_x = \frac{1}{\lambda^3 N!} Q(\beta, V, N) \]
\[Q = \int dx^N e^{-\beta U(x^N)}\]

Pressue is related to Free energy as

  • \[p= - \frac{\partial F}{\partial V} = \frac{\partial }{\partial V} log \int dx^N e^{-\beta U(x^N)}\]

Volume dependence of partition function is in the integration limits! As volume grows so does partition function. Therefore p is always positive. We can thus conclude that in equilibirum pressure is always a positive quantity

Thermal wavelength and broad validity of classical picture of fluids#

\[Z = \frac{1}{N!} \frac{1}{\lambda_T^{3N}} \int dx^N e^{-\beta U(x^N)}\]

Classical mechanics is valid as long as De Broglie wavelength is small compared to intermolecular length scales

\[\lambda_{DB} = \frac{h}{p} \sim \frac{h}{\langle p \rangle} = \frac{h}{(8k_B T m/\pi)^{1/2}} \approx \lambda_T\]

For a dilute gas, relevant lengths are \(\rho^{-1/3}\) (typical distance) and \(\sigma\) (diamater of particles)

  • \(\lambda_T < \sigma \) (classical regime is good)

\(H_2O\) vs \(D_2O\) melting and freezing points#

  • Classical descirption predicts configurational parition function to be independent of paritcle masses. Hence equation of state \(p = p(\beta, \rho)\) is independent of masses.

  • Thus if translationa, rotational and virbational degrees of freedo are well described by classical mechanics, \(H_2O\) and \(D_2O\) will have the same equation of state. Experiments however predict desnity maximum \(\rho(T)\) at \(p=1 atm\) at 4 C for \(H_2O\) and 10C in \(D_2O\). The freezing point is higher for \(D_2 O\) as well

  • Water is a weird fluid. Size is about \(3 A\) and thermal wavelength is around \(0.3 A\). The role of quantum fluctuations is to blurr the positions of water molecules. As the atomic mass increases, bluring goes down and fluids gets ordered. This is why \(D_2 O\) ice melts at a higher temperature, e.g need to supply more energy to disrupt the order.

Reduced configruational distribution functions#

  • Hard fact: \(P(r^N)\) and Q do not factorize because of inteparticle interactions with stronger interactions implying stronger correlations in positions.

  • Probability of many-body systems: to find particle \(1\) at \(r_1\), \(2\) and \(r_2\) …:

\[\boxed{P(r^N) = P(r_1, r_2,...r_N)}\]

  • Marginal(ized) probability

\[\boxed{\rho^{2/N} (r_1, r_2) = \int dr_3... dr_N P(r_1, r_2,...r_N)}\]

  • Marginal(ized) probability for any particle 1 and 2

\[\boxed{\rho^{2/N} (r_1, r_2) = N(N-1)\int dr_3... dr_N P(r_1, r_2,...r_N)}\]

  • Marginal(ized) probability for any particle 1 and 2…n

\[\boxed{\rho^{n/N} (r_1, r_2,... r_n) = \frac{N!}{(N-n)!}\int dr^{N-n} P(r_1, r_2,...r_N)}\]

Radial distribution function (RDF)#

For an isotropic liquid:

\[\rho = \rho = \frac{N}{V}\]

For an ideal gas:

\[\rho^{2/N} = \frac{N(N-1)}{V^2} = \rho^2 (1-N^{-1}) \approx \rho^2\]

  • To measure degree of spatial correlations we introduce Radial Distribution function (RDF):

\[\boxed{g(r_1, r_2) = \frac{\rho^{2/N}(r_1, r_2)}{\rho^2}}\]

  • RDF for isotropic fluids is \(\boxed{g(r) = g(|r_2-r_1|)}\)

  • Conditional probabilioty density of finding a particle at r distance away from a tagged particle placed at origin:

\[\rho^{2/N}(0,r) = \rho g(r)\]

Coordination shells and structure in fluids#

Reversible work theorem and potential of mean force#

\[\boxed{g(r) = e^{-\beta w(r)}}\]
\[\boxed{w(r) = - \beta^{-1} log [g(r)]}\]

  • \(w(r)\) Reversible work to bring two particles from infinity to distance r

  • \(g(r)\) Radial distribution function

\[\Big \langle -\frac{d U(r^N)}{dr_1} \Big \rangle_{r_1, r_2} = -\frac{\int dr_3...dr_N \frac{dU(r^N)}{dr_1} e^{-\beta U}}{\int dr_3...r_N e^{-\beta U}} = \frac{\beta^{-1} \frac{d}{d r_1} \int dr_3...dr_N e^{-\beta U}}{\int dr_3...dr_N e^{-\beta U}} = \beta^{-1} \frac{d}{d r_1} log \int dr_3...dr_N e^{-\beta U} \]
\[\Big \langle -\frac{d U(r^N)}{dr_1} \Big \rangle_{r_1, r_2} = \beta^{-1} \frac{d}{d r_1} log N (N-1)\int dr_3...dr_N e^{-\beta U} = \beta^{-1} g(r_1, r_2)\]

Thermodynamic properites of \(g(r)\)#

\[\langle E \rangle = N \Big\langle \frac{p^2}{2m} \Big\rangle + \Big\langle \sum_{j>i} u(|r_i - r_j|)\Big \rangle\]
\[\boxed{E/N = \frac{3}{2}k_B T +\frac{1}{2}\rho \int dr g(r) u(r) }\]

Low density approximation for \(g(r)\)#

\[w(r) =u(r) +\Delta w(r) \]
\[ g(r) = e^{-\beta u(r)} \Big (1 +O(\rho) \Big)\]

Density expanesion and virial coefficients#

\[\beta p = \rho + B_2(T) \rho^2+ O(\rho^3)\]
\[B_2(T) = -\frac{1}{2} \int dr (e^{-\beta u(r)}-1) \]