Fugacity

From Wikipedia, the free encyclopedia

This article may require cleanup to meet Wikipedia's quality standards.
Please discuss this issue on the talk page or replace this tag with a more specific message.
This article has been tagged since September 2006.

1 Definition
2 Technical Detail
- 2.1 Fugacity and Chemical Potential
- 2.2 Alternative Methods for calculating fugacity
3 External links

[edit] Definition

Fugacity is an measure of chemical potential in the form of an 'adjusted pressure.' It directly relates to the tendency of a substance to prefer one phase (liquid, solid, gas) over another. At a fixed temperature and pressure, water (for example) will have a different fugacity for each phase. The phase with the lowest fugacity will be the most favorable; the substance minimizes Gibbs free energy. Fugacity, therefore, is a useful engineering tool for predicting the phase state of multi-component mixtures at various temperatures and pressures without doing the actual lab test. As well as predicting the preferred phase of a single substance, fugacity is also useful for multicomponent equilibrium involving any combination of solid, liquid and gas equilibria.

Fugacity is also a useful way to explain otherwise obscure behaviors of substances. For example, most people are at a loss to describe why liquid water in a dish on a typical day eventually dries up. Why is it that below the boiling point, the water spontaneously changes phases, evaporating the liquid to form gas? Fugacity provides an answer: When the temperature and pressure is such that the relative humidity in the room is less than 100%, the calculated fugacity of water vapor will be lower than either liquid or solid water. Therefore, liquid water has a lower Gibbs free energy in the vapor phase and will evaporate. This example illustrates the etymology of fugacity, from the Latin effugit meaning to escape: it represents the tendency for a substance to "escape" from a certain phase.

Fugacity is not a physical property of a substance; rather it is a calculated property which is intrinsically related to chemical potential. When a system approaches the ideal gaseous state (very low pressure), the chemical potential approaches negative infinity, which for the purposes of mathematical modeling is undesirable; under the same conditions, fugacity approaches zero, and the fugacity coefficient (defined below) approaches 1, both of which are more desirable.

[edit] Technical Detail

Fugacity is a state function of matter at fixed temperature. It only becomes useful when dealing with substances other than an ideal gas. For an ideal gas, fugacity is equal to pressure. In the real world, no substance acts like an ideal gas, so we use fugacity not only to describe non-ideal gases, but liquids and solids as well.

For gases at low pressures (where the ideal gas law is a good approximation), fugacity is roughly equal to pressure. Thus, for an ideal gas, the ratio $\phi = f/P \,$ between fugacity $f\,$ and pressure $P\,$ (the fugacity coefficient) is equal to 1. This ratio can be thought of as 'how closely the substance behaves like an ideal gas,' based on how far it is from 1.

For a given temperature $T\,$ , the fugacity $f\,$ satisfies the following differential relation:

$d \ln {f \over f_0} = {dG \over RT} = {{\tilde V dP} \over RT} \,$

where $G\,$ is the Gibbs free energy, $R\,$ is the gas constant, $\tilde V\,$ is the fluid's molar volume, and $f_0\,$ is a reference fugacity which is generally taken as that of an ideal gas at 1 bar. For an ideal gas, when $f = P$ , this equation reduces to the ideal gas law.

Thus, for any two mutually-isothermal physical states, represented by subscripts 1 and 2, the ratio of the two fugacities is as follows:

$f_2 / f_1 = exp \left ({1 \over RT} \int_{G_1}^{G_2} dG \right) = exp \left ({1 \over RT} \int_{P_1}^{P_2} \tilde V\,dP \right) \,$

The concept of fugacity was introduced by American chemist Gilbert N. Lewis in his paper "The osmotic pressure of concentrated solutions, and the laws of the perfect solution," J. Am. Chem. Soc. 30, 668-683 (1908).

[edit] Fugacity and Chemical Potential

For every pure substance, we have the relation

$dG = - SdT + VdP \,$

for Gibbs free energy

and

$\mu _i = \left( {\frac{{\partial G}} {{\partial n_i }}} \right)_{T,P,n_{j \ne i} }$

for the chemical potential of the gas $i$

If we choose to study only one pure ideal gas, we can divide both members of the equation by $n$ and we’ll have

$d\mu = - \bar SdT + \bar VdP$

where $\bar S$ and $\bar V$ are the molar entropy and molar volume for that gas

We choose to study a process where $T$ remains constant.

Therefore,

$d\mu = \bar VdP$

We can integrate this expression remembering the chemical potential is a function of $T$ and $P$ . We must also set a reference state. In this case, for an ideal gas the only reference state will be the pressure, and we set $P$ = 1 bar.

$\int_{\mu^\circ }^\mu {d\mu } = \int_{P^\circ }^P {\bar VdP}$

Now, for the ideal gas $\bar V = \frac{{RT}}{P}$

$\mu - \mu ^\circ = \int_{P^\circ }^P {\frac{{RT}} {P}dP} = RT\ln \frac{P} {{P^\circ }}$

Reordering, we get

$\mu = \mu ^\circ + RT\ln \frac{P} {{P^\circ }}$

Which gives the chemical potential for an ideal gas in an isothermal process, where the reference state is $P$ =1 bar.

For a real gas, we cannot calculate $\int_{P^\circ }^P {\bar VdP}$ because we do not have a simple expression for a real gas’ molar volume. On the other hand, even if we did have one expression for it (we could use the Van der Waals equation, Redlich-Kwong or any other equation of state), it would depend on the substance being studied and would be therefore of a very limited usability.

We would like the expression for a real gas’ chemical potential to be similar to the one for an ideal gas.

We can define a magnitude, called fugacity, so that the chemical potential for a real gas becomes

$\mu = \mu ^\circ + RT\ln \frac{f} {{f^\circ }}$

with a given reference state (discussed later).

We can see that for an ideal gas, it must be $f = P$

But for $P \to 0$ , every gas is an ideal gas. Therefore, fugacity must obey the limit equation

$\mathop {\lim }_{P \to 0} \frac{f} {P} = 1$

We determine $f$ by defining a function

$\Phi = \frac{{P\bar V - RT}} {P}$

We can obtain values for $Φ$ experimentally easily by measuring $V$ , $T$ and $P$ . (note that for an ideal gas, $Φ$ = 1)

From the expression above we have

$\bar V = \frac{{RT}} {P} + \Phi$

We can then write

$\int_{\mu ^\circ }^\mu {d\mu } = \int_{P^\circ }^P {\bar VdP} = \int_{P^\circ }^P {\frac{{RT}} {P}dP} + \int_{P^\circ }^P {\Phi dP}$

Where

$\mu = \mu ^\circ + RT\ln \frac{P} {{P^\circ }} + \int_{P^\circ }^P {\Phi dP}$

Since the expression for an ideal gas was chosen to be $\mu = \mu ^\circ + RT\ln \frac{f} {{f^\circ }}$ ,we must have

$\mu ^\circ + RT\ln \frac{f} {{f^\circ }} = \mu ^\circ + RT\ln \frac{P} {{P^\circ }} + \int_{P^\circ }^P {\Phi dP}$

$\Rightarrow RT\ln \frac{f} {{f^\circ }} - RT\ln \frac{P} {{P^\circ }} = \int_{P^\circ }^P {\Phi dP}$

$RT\ln \frac{{fP^\circ }} {{Pf^\circ }} = \int_{P^\circ }^P {\Phi dP}$

Suppose we choose $P \to 0$ . Since $\mathop {\lim }_{P \to 0} f = P$ , we obtain

$RT\ln \frac{f} {P} = \int_0^P {\Phi dP}$

The fugacity coefficient will then verify

$\ln \phi = \frac{1} {{RT}}\int_0^P {\Phi dP}$

The integral can be evaluated via graphical integration if we measure experimentally values for $Φ$ while varying $P$ .

We can then find the fugacity coefficient of a gas at a given pressure $P$ and calculate

$f = \phi P\,$

The reference state for the expression of a real gas’ chemical potential is taken to be “ideal gas, at $P$ = 1 bar and work $T$ ”. Since in the reference state the gas is considered to be ideal (it is an hypothetical reference state), we can write that for the real gas

$\mu = \mu ^\circ + RT\ln \frac{f} {{P^\circ }}$

[edit] Alternative Methods for calculating fugacity

If we suppose that $Φ$ is constant between 0 and $P$ (assuming it is possible to do this approximation), we have

$\frac{f} {P} = e^{\frac{{\Phi P}} {{RT}}}$

Expanding in Taylor series about 0,

$\frac{f} {P} \approx 1 + \frac{{\Phi P}} {{RT}} = 1 + \frac{1} {{RT}}\left( {\frac{{P\bar V - RT}} {P}} \right)P = 1 + \frac{{P\bar V}} {{RT}} - 1 = \frac{{P\bar V}} {{RT}}$

Finally, we get

$f \approx \frac{{P^2 \bar V}} {{RT}}$

This formula allows us to calculate quickly the fugacity of a real gas at $P$ , $T$ , given a value for V (which could be determined using any equation of state), if we suppose is constant between 0 and $P$ .

We can also use generalized charts for gases in order to find the fugacity coefficient for a given reduced temperature.

Fugacity could be considered a “corrected pressure” for the real gas, but should never be used to replace pressure in equations of state (or any other equations for that matter). That is, it is false to write expressions such as

$f V = n R T$

Fugacity is strictly a tool, conveniently defined so that the chemical potential equation for a real gas turns out to be similar to the equation for an ideal gas.