Difference between revisions of "Second virial coefficient"

From SklogWiki
Jump to: navigation, search
(Admur and Mason mixing rule: correction)
(Added 'unkown' mixing rule)
Line 66: Line 66:
:<math>B_{ {\mathrm {mix}} } =  \sum_{i=1}^{n} \sum_{j=1}^{n} B_{ij} x_i x_j</math>
:<math>B_{ {\mathrm {mix}} } =  \sum_{i=1}^{n} \sum_{j=1}^{n} B_{ij} x_i x_j</math>
where <math>x_i</math> and <math>x_j</math> are the mole fractions of the <math>i</math>th and <math>j</math>th component gasses of the mixture.
where <math>x_i</math> and <math>x_j</math> are the mole fractions of the <math>i</math>th and <math>j</math>th component gasses of the mixture.
(<ref>I am not sure where this mixing rule was published</ref>)
:<math>B_{ij} = \frac{\left(B_{ii}^{1/3}+B_{jj}^{1/3}\right)^3}{8}</math>
==See also==
==See also==
*[[Virial equation of state]]
*[[Virial equation of state]]

Latest revision as of 15:36, 10 December 2019

The second virial coefficient is usually written as B or as B_2. The second virial coefficient represents the initial departure from ideal-gas behaviour. The second virial coefficient, in three dimensions, is given by

B_{2}(T)= - \frac{1}{2} \int \left(  \exp\left(-\frac{\Phi_{12}({\mathbf r})}{k_BT}\right) -1 \right) 4 \pi r^2 dr

where \Phi_{12}({\mathbf r}) is the intermolecular pair potential, T is the temperature and k_B is the Boltzmann constant. Notice that the expression within the parenthesis of the integral is the Mayer f-function.

In practice the integral is often very hard to integrate analytically for anything other than, say, the hard sphere model, thus one numerically evaluates

B_{2}(T)= - \frac{1}{2} \int \left( \left\langle \exp\left(-\frac{\Phi_{12}({\mathbf r})}{k_BT}\right)\right\rangle -1 \right) 4 \pi r^2 dr


 \left\langle \exp\left(-\frac{\Phi_{12}({\mathbf r})}{k_BT}\right)\right\rangle

for each r using the numerical integration scheme proposed by Harold Conroy [1][2].

Isihara-Hadwiger formula[edit]

The Isihara-Hadwiger formula was discovered simultaneously and independently by Isihara [3] [4] [5] and the Swiss mathematician Hadwiger in 1950 [6] [7] [8] The second virial coefficient for any hard convex body is given by the exact relation



\frac{B_2}{V}=1+3 \alpha


\alpha = \frac{RS}{3V}

where V is the volume, S, the surface area, and R the mean radius of curvature.

Hard spheres[edit]

For the hard sphere model one has [9]

B_{2}(T)= - \frac{1}{2} \int_0^\sigma \left(\langle 0\rangle -1 \right) 4 \pi r^2 dr

leading to

B_{2}=  \frac{2\pi\sigma^3}{3}

Note that B_{2} for the hard sphere is independent of temperature. See also: Hard sphere: virial coefficients.

Van der Waals equation of state[edit]

For the Van der Waals equation of state one has:

B_{2}(T)=  b -\frac{a}{RT}

For the derivation click here.

Excluded volume[edit]

The second virial coefficient can be computed from the expression

B_{2}=  \frac{1}{2} \iint v_{\mathrm {excluded}} (\Omega,\Omega') f(\Omega) f(\Omega')~ {\mathrm d}\Omega {\mathrm d}\Omega'

where v_{\mathrm {excluded}} is the excluded volume.

Admur and Mason mixing rule[edit]

The second virial coefficient for a mixture of n components is given by (Eq. 11 in [10])

B_{ {\mathrm {mix}} } =  \sum_{i=1}^{n} \sum_{j=1}^{n} B_{ij} x_i x_j

where x_i and x_j are the mole fractions of the ith and jth component gasses of the mixture.



B_{ij} = \frac{\left(B_{ii}^{1/3}+B_{jj}^{1/3}\right)^3}{8}

See also[edit]


Related reading