Loading [MathJax]/jax/output/CommonHTML/jax.js
Submit manuscript...
Open Access Journal of
eISSN: 2641-9335

Mathematical and Theoretical Physics

Mini Review Volume 1 Issue 5

On the Van der Waals equation

M Aposto

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

Download PDF

Abstract

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

Introduction

Let us consider a classical gas of N identical point like particles enclosed in a volume V with the interaction energy

E=12ijU(rirj)  (1)

where U(rirj) is the two-particle interaction and i, j = 1, 2, ...N denote the particles. For point like particles U(rirj)  is a function only of |rirj| .The correction to the free energy is given by

ΔF=Τln(1VNeβEdr1..drN)=       =Τln[1VN(eβE1)dr1..drN+1]  (2)

where Τ=1/β  is the temperature. We assume that the interaction U(rirj) 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(N1)2V2(eβU(r1r2)1)dr1dr2+1]     Τln[1N22V2(1eβU(r1r2))dr1dr2];  (3)

if the amount of gas is still sufficiently large (such as to neglect the surface effects), we may write

ΔF=Τln(1N2VB),  (4)

B=12(1eβ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π2r00dr.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=ΤN2V2B1N2VB  (8)

and the total pressure

p=NTV+ΤN2V2B1N2VB=NTV(1+NVB1N2VB)  (9)

or

p=NTV1N2VB(11/N)1N2VB  (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+NVB1N2VB)NTV(1+NBV)=    =NTV(1+NBV)+N2aV2NTVNb+N2aV2,          (11)

which leads to the van der Waals equation1

(pN2aV2)(VNb)=NT (12)

(written, usually, with a<0 , as for attractive interaction at large distances. However, if we neglect N2B/V1 in the denominator of equation (11), we should also neglect NB/V=(N2B/V)/NN2B/V1 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μ/T1, 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 N1 , 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(rirj)  (13)

where the prime means ji ; the total energy can be written as

E=12iΦ(ri).  (14)

The integration over ri in equation (2) must take into account the requirement ji , i.e., rjri for any pairs (i,j) . It follows that the correction to the free energy is given by

ΔF=Τln(1VNeβEdr1...drN) , (15)

where the prime means rjri . We perform this integration over a volume VNb , 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

Φ=ijU(rj)=NV2r0drU(r)=2NaV (16)

where b12.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((VNb)NVNe12βNΦ)==12NΦNTln(1NbV)=N2aVNTln(1NbV). (17)

We get the pressure

p=NTV+N2aV2+NTNb/V21Nb/V==NTV(1+Nb/V1Nb/V)+N2aV2=NTVNb+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)Ne12βNΦ]=      =Tln(1V'eβφ(r)dr)+12NΦ; (20)

The integral in equation (20) can be written as

      'eβφ(r)dr=eβφ(r)drbieβφ(ri)dr==eβφ(r)drNbVeβφ(r)dr=1NbVeβφ(r)dr, (21)

Such that we get

      'eβφ(r)dr=eβφ(r)drbieβφ(ri)dr==eβφ(r)drNbVeβφ(r)dr=(1NbV)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.

Acknowledgements

None.

Conflict of interest

The authors declare that there are no conflict of interest regarding the publication of this paper.

References

Creative Commons Attribution License

©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.