Mini Review Volume 1 Issue 5
Department of Theoretical Physics, Institute of Atomic Physics, Romania
Correspondence: M Apostol, Department of Theoretical Physics, Institute of Atomic Physics, Magurele-Bucharest MG-6,PO Box MG-35, Romania
Received: September 29, 2018 | Published: October 15, 2018
Citation: Apostol M. On the Van der Waals equation. Open Acc J Math Theor Phy. 2018;1(5):215-217. DOI: 10.15406/oajmtp.2018.01.00037
It is shown that the usual derivations of the van der Waals equation are inconsistent. The van der Waals equation is derived here in a more general framework, which employs the mean-field approximation.
PACS: 05.20.Gg; 51.30.+1; 0570.Ce
Keywords: Van der Waals equation, non-ideal gases, equation of state
Let us consider a classical gas of N identical point like particles enclosed in a volume V with the interaction energy
E=12∑i≠jU(ri−rj) (1)
where U(ri−rj) is the two-particle interaction and i, j = 1, 2, ...N denote the particles. For point like particles U(ri−rj) is a function only of |ri−rj| .The correction to the free energy is given by
ΔF=−Τln(1VN∫e−βEdr1..drN)= =−Τln[1VN∫(e−βE−1)dr1..drN+1] (2)
where Τ=1/β is the temperature. We assume that the interaction U(ri−rj) is short range (integrable) and, moreover, it proceeds by collisions. More, we assume that the amount of gas is so small, that only one collision takes place at one time; since the free energy is a function of the form F=Nf(Τ,V/N), it follows that the result will be valid for any amount of gas. In addition, we neglect the multi-particle collisions, as for a rarefied gas. In these conditions equation (2) can be written as
ΔF=−Τln[N(N−1)2V2∫(e−βU(r1−r2)−1)dr1dr2+1]≃ ≃−Τln[1−N22V2∫(1−e−βU(r1−r2))dr1dr2]; (3)
if the amount of gas is still sufficiently large (such as to neglect the surface effects), we may write
ΔF=−Τln(1−N2VB), (4)
B=12∫(1−e−βU(r))dr. (5)
The interaction U(r) is very repulsive for small r and very small in comparison with T at large distances; consequently, we may approximate the parameter B as
B=b+a/Τ , (6)
b=2π2r0∫0dr.r2,a=2π∫2r0dr.r2U(r), (7)
2r0 being a scale distance which separates the "hard-core" part of U(r) . From equation (4) we get the correction to pressure
Δp=ΤN2V2B1−N2VB (8)
and the total pressure
p=NTV+ΤN2V2B1−N2VB=NTV(1+NVB1−N2VB) (9)
or
p=NTV1−N2VB(1−1/N)1−N2VB (10)
Since N is still much larger than 1, we get p=NT/V , which shows that the interaction does not change the equation of state of the ideal gases.
It is customary to use the approximation
p=NTV(1+NVB1−N2VB)≃NTV(1+NBV)= =NTV(1+NBV)+N2aV2≃NTV−Nb+N2aV2, (11)
which leads to the van der Waals equation1
(p−N2aV2)(V−Nb)=NT (12)
(written, usually, with a<0 , as for attractive interaction at large distances. However, if we neglect N2B/V≪1 in the denominator of equation (11), we should also neglect NB/V=(N2B/V)/N≪N2B/V≪1 in the numerator, which removes any effect of the interaction in the final formula.
The difficulties related to the van der Waals equation have been discussed previously.2 The expansion of the thermodynamic potential Ω=−pV in powers of ξ=(mT/2πℏ2)3/2eμ/T≪1, where μ is the chemical potential (and ℏ is Planck’s constant),3 exhibits the same drawback; indeed, we Ω=−Τln(1+ξV+....) and N=Vξ+..., such that Ω=−pV=−Τln(1+N....) which is written usually as Ω=−pV=−Τln(1+N....)≃−NT; however, this approximation is valid for N≪1 , a condition which is not satisfied.
The difficulties with the usual derivations of the van der Waals equation reside in the fact that an ensemble of interacting particles are not statistically independent; consequently, we need to solve first the problem of interaction. First, we leave aside correlations (which may lead to sound waves in a gas (collective excitations), which contribute to the statistical properties, though to a lesser extent); then, from equation (1) we see that each particle i is acted by the rest of the particles by a potential
Φ(ri)=ˋ∑jU(ri−rj) (13)
where the prime means j≠i ; the total energy can be written as
E=12∑iΦ(ri). (14)
The integration over ri in equation (2) must take into account the requirement j≠i , i.e., rj≠ri for any pairs (i,j) . It follows that the correction to the free energy is given by
ΔF=−Τln(1VN∫e−βEdr1...drN) , (15)
where the prime means rj≠ri . We perform this integration over a volume V−Nb , where b is an "excluded" volume for each particle. In a mean-field approximation the potential Φ(ri) does not depend on Φ(ri) and the summation in equation (13) is computed by
Φ=i∑jU(rj)=NV∫2r0drU(r)=2NaV (16)
where b≃12.4π(2r0)3/3 ; the factor 1/2 in estimating the volume b arises from the fact that the independent integration with respect to any pair ri,rj counts twice the excluded volume. The change in the free energy becomes (equation (15))
ΔF=−Tln((V−Nb)NVNe−12βNΦ)==12NΦ−NTln(1−NbV)=N2aV−NTln(1−NbV). (17)
We get the pressure
p=NTV+N2aV2+NTNb/V21−Nb/V==NTV(1+Nb/V1−Nb/V)+N2aV2=NTV−Nb+N2aV2; (18)
Hence, it follows immediately the van der Waals equation (12). It is worth noting that we use the same cutoff 2r0 to account for the prime in the integration in equation (15) and the prime in the potential given by equation (16). In general, we should employ two distinct cutoffs, but if the potential has a very sharp repulsive hard core the two cutoffs coincide. Also, we note that the mean-field approximation is equivalent with restricting ourselves only to one-particle elementary excitations (quasi-particles). We have assumed that a and b do not depend on V, and b>0 (by its definition).
It may happen that an external potential φ(r) is present and we wish to estimate its effects in the presence of the internal interaction. Then, the total energy is written as
E=∑iφ(ri)+12NΦ (19)
where Φ is a mean-field potential (also, the external potential may derive from a mean-field approximation). The ensemble has not translational symmetry anymore, since φ(r) depends on position. The change in the free energy is
ΔF=−Tln(1VN∫'e−β∑iφ(ri)−12NΦdr1..drN)= =−Tln[(1V∫'e−βφ(r)dr)Ne−12βNΦ]= =−Tln(1V∫'e−βφ(r)dr)+12NΦ; (20)
The integral in equation (20) can be written as
∫'e−βφ(r)dr=∫e−βφ(r)dr−b∑ie−βφ(ri)dr==∫e−βφ(r)dr−NbV∫e−βφ(r)dr=1−NbV∫e−βφ(r)dr, (21)
Such that we get
∫'e−βφ(r)dr=∫e−βφ(r)dr−b∑ie−βφ(ri)dr==∫e−βφ(r)dr−NbV∫e−βφ(r)dr=(1−NbV)∫e−βφ(r)dr, (22)
we can see that we recover the van der Waals equation and get also the effects of the external field; for a finite ensemble, such effects may include surface contributions.
The general way of deriving the van der Waals equation consists in estimating the effects of the interaction in the mean-field equation (16) and taking account of the excluded volume, according to the definition of the interaction energy.
None.
The authors declare that there are no conflict of interest regarding the publication of this paper.
©2018 Apostol. This is an open access article distributed under the terms of the, which permits unrestricted use, distribution, and build upon your work non-commercially.