Introduction to Mathematical Physics/N body problems and statistical equilibrium/Van der Waals gas

From Wikibooks, open books for an open world
Jump to: navigation, search

Van der Waals\footnote{Johannes Diderik van der Waals received the physics Nobel price in 1910.} gas model\index{Van der Walls gas} is a model that also relies on classical approximation, that is on the repartition function that can be written:


Z=\frac{1}{N!} \frac{1}{h^{3N}}\int\int d^3p_1\dots d^3p_N
d^3r_1\dots d^3r_N  e^{-\frac{H}{k_BT}}

where the function H(p_1 ,\dots ,p_N,r_1 ,\dots ,r_N) is the hamiltonian of the system.

H(p_1 ,\dots ,p_N,r_1 ,\dots ,r_N) = \frac{1}{2m}\sum p_i^2+U(r_1, \dots, r_N)

As for the perfect gas, integration over the p_i's is immediate:

Z=\frac{1}{N!}\left(\frac{2\pi m k_BT }{h^2}\right)^{\frac{3N}{2}}Y

with

eqYint

Y=\int d^3q_1\dots d^3q_N  e^{-\frac{U(r_1,\dots,r_N)}{k_BT}}

We can simplify the evaluation of previuos integral in neglecting correlations between particles. A particle number i doesn't feel each particle but rather the average influence of the particles cloud surrounding the considered particle; this approximation is called mean field approximation \index{mean field}. It leads, as if the particles were actually independant to a factorization of the repartition function. Interaction potential U(r_1,\dots,r_N) becomes:

U(r_1,\dots,r_N) =\frac{1}{2}\sum U_{e}(r_i)

where U_e(r_i) is the "effective potential" that describes the mean interaction between particle i and all the other particles. It depends only on position of particle i. Function Y introduced at equation ---eqYint--- can be factorized:

Y=\left(\int dr e^{-\frac{U_{e}(r)}{2k_BT}}\right)

In a mean field approximation framework, it is considered that particle distribution is uniform in the volume. Effective potential has thus to traduce an mean attraction. Indeed, it can be shown that at large distances two molecules are attracted, potential varying like \frac{1}{r^6}. This can be proofed using quantum mechanics but this is out of the frame of this book. At small distances however, molecules repel themselves strongly. Effective potential undergone by a particle at position r=0 can thus be modelled by function:

\begin{matrix}
U_e(r)&=&\infty  \mbox{ if }  r<r_0\\
U_e(r)&=&-C \mbox{ if } r>r_0\\
\end{matrix}

Quantity Y is


Y=\left((V-bN)e^{\frac{aN}{k_BTV}}\right)^N

Quantity bN represents the excluded volume of the particle. It is proportional to N because N particles occupying each a volume b occupy a Nb. On another hand mean potential U_0 felt by the test particle depends on ratio N/V. Usually, one sets:

U_0=-2a\frac{N}{V}

Introducing

F=-k_BT\log Z

and using

p=-\frac{\partial F}{\partial V}

one finally obtains Van der Waals state equation:

\left(p+a\frac{n^2}{V^2}\right)(V-bN)=Nk_BT

Van der Waals model allows to describe liquid--vapour phase transition.

Over a critical temperature, energy doesnt present any local minimum.}
figvanavant
Below the critical temperatre, energy presents a local minimum. Two phases then appear, characteriszed by points A_{0} et B_0.}
figvanapres

When temperature is lower than critical temperature T_c, the energy of the system for a given volume V_M with V_{A_0}<V_M<V_{B_0} is not F_v(V_M). Indeed, system evolves to a state of lower energy (see figure figvanapres) with energy:

F(V_M)=F_{A_0}+\frac{V_M-V_{A_0}}{V_{B_0}-V_{A_0}}(F_{B_0}-F_{A_0})

The apparition of a local minimum of F corresponds to the apparition of two phases.\index{phase transition}