Physics:Kelvin equation

Short description: Equation describing the change in vapour pressure due to a curved liquid–vapor interface

The Kelvin equation describes the change in vapour pressure due to a curved liquid–vapor interface, such as the surface of a droplet. The vapor pressure at a convex curved surface is higher than that at a flat surface. The Kelvin equation is dependent upon thermodynamic principles and does not allude to special properties of materials. It is also used for determination of pore size distribution of a porous medium using adsorption porosimetry. The equation is named in honor of William Thomson, also known as Lord Kelvin.

Formulation

The original form of the Kelvin equation, published in 1871, is: [math]\displaystyle{ p(r_1 , r_2) = P - \frac {\gamma\, \rho _{\rm vapor} } {(\rho_{\rm liquid} - \rho_{\rm vapor})} \left ( \frac {1}{r_1} + \frac {1}{r_2} \right ), }[/math] where:

• [math]\displaystyle{ p(r) }[/math] = vapor pressure at a curved interface of radius [math]\displaystyle{ r }[/math]
• [math]\displaystyle{ P }[/math] = vapor pressure at flat interface ([math]\displaystyle{ r = \infty }[/math]) = [math]\displaystyle{ p_{eq} }[/math]
• [math]\displaystyle{ \gamma }[/math] = surface tension
• [math]\displaystyle{ \rho _{\rm vapor} }[/math] = density of vapor
• [math]\displaystyle{ \rho _{\rm liquid} }[/math] = density of liquid
• [math]\displaystyle{ r_1 }[/math] , [math]\displaystyle{ r_2 }[/math] = radii of curvature along the principal sections of the curved interface.

This may be written in the following form, known as the Ostwald–Freundlich equation: [math]\displaystyle{ \ln \frac{p}{p_{\rm sat}} = \frac{2 \gamma V_\text{m}}{rRT}, }[/math] where [math]\displaystyle{ p }[/math] is the actual vapour pressure, [math]\displaystyle{ p_{\rm sat} }[/math] is the saturated vapour pressure when the surface is flat, [math]\displaystyle{ \gamma }[/math] is the liquid/vapor surface tension, [math]\displaystyle{ V_\text{m} }[/math] is the molar volume of the liquid, [math]\displaystyle{ R }[/math] is the universal gas constant, [math]\displaystyle{ r }[/math] is the radius of the droplet, and [math]\displaystyle{ T }[/math] is temperature.

Equilibrium vapor pressure depends on droplet size.

• If the curvature is convex, [math]\displaystyle{ r }[/math] is positive, then [math]\displaystyle{ p \gt p_{\rm sat} }[/math]
• If the curvature is concave, [math]\displaystyle{ r }[/math] is negative, then [math]\displaystyle{ p \lt p_{\rm sat} }[/math]

As [math]\displaystyle{ r }[/math] increases, [math]\displaystyle{ p }[/math] decreases towards [math]\displaystyle{ p_{sat} }[/math], and the droplets grow into bulk liquid.

If the vapour is cooled, then [math]\displaystyle{ T }[/math] decreases, but so does [math]\displaystyle{ p_{\rm sat} }[/math]. This means [math]\displaystyle{ p/p_{\rm sat} }[/math] increases as the liquid is cooled. [math]\displaystyle{ \gamma }[/math] and [math]\displaystyle{ V_\text{m} }[/math] may be treated as approximately fixed, which means that the critical radius [math]\displaystyle{ r }[/math] must also decrease. The further a vapour is supercooled, the smaller the critical radius becomes. Ultimately it can become as small as a few molecules, and the liquid undergoes homogeneous nucleation and growth.

A system containing a pure homogeneous vapour and liquid in equilibrium. In a thought experiment, a non-wetting tube is inserted into the liquid, causing the liquid in the tube to move downwards. The vapour pressure above the curved interface is then higher than that for the planar interface. This picture provides a simple conceptual basis for the Kelvin equation.

The change in vapor pressure can be attributed to changes in the Laplace pressure. When the Laplace pressure rises in a droplet, the droplet tends to evaporate more easily.

When applying the Kelvin equation, two cases must be distinguished: A drop of liquid in its own vapor will result in a convex liquid surface, and a bubble of vapor in a liquid will result in a concave liquid surface.

History

The form of the Kelvin equation here is not the form in which it appeared in Lord Kelvin's article of 1871. The derivation of the form that appears in this article from Kelvin's original equation was presented by Robert von Helmholtz (son of German physicist Hermann von Helmholtz) in his dissertation of 1885. In 2020, researchers found that the equation was accurate down to the 1nm scale.

Derivation using the Gibbs free energy

The formal definition of the Gibbs free energy for a parcel of volume [math]\displaystyle{ V }[/math], pressure [math]\displaystyle{ P }[/math] and temperature [math]\displaystyle{ T }[/math] is given by:

[math]\displaystyle{ G=U+pV-TS, }[/math]

where [math]\displaystyle{ U }[/math] is the internal energy and [math]\displaystyle{ S }[/math] is the entropy. The differential form of the Gibbs free energy can be given as

[math]\displaystyle{ dG=-S dT + V dP + \sum_{i=1}^ k \mu_i dn_i, }[/math]

where [math]\displaystyle{ \mu }[/math] is the chemical potential and [math]\displaystyle{ n }[/math] is the number of moles. Suppose we have a substance [math]\displaystyle{ x }[/math] which contains no impurities. Let's consider the formation of a single drop of [math]\displaystyle{ x }[/math] with radius [math]\displaystyle{ r }[/math] containing [math]\displaystyle{ n_x }[/math] molecules from its pure vapor. The change in the Gibbs free energy due to this process is

[math]\displaystyle{ \Delta G = G_d - G_v, }[/math]

where [math]\displaystyle{ G_d }[/math] and [math]\displaystyle{ G_v }[/math] are the Gibbs energies of the drop and vapor respectively. Suppose we have [math]\displaystyle{ N_i }[/math] molecules in the vapor phase initially. After the formation of the drop, this number decreases to [math]\displaystyle{ N_f }[/math], where

[math]\displaystyle{ N_f = N_i - n_x. }[/math]

Let [math]\displaystyle{ g_v }[/math] and [math]\displaystyle{ g_l }[/math] represent the Gibbs free energy of a molecule in the vapor and liquid phase respectively. The change in the Gibbs free energy is then:

[math]\displaystyle{ \Delta G = N_f g_v + n_x g_l + 4 \pi r^2 \sigma - N_i g_v, }[/math]

where [math]\displaystyle{ 4 \pi r^2 \sigma }[/math] is the Gibbs free energy associated with an interface with radius of curvature [math]\displaystyle{ r }[/math] and surface tension [math]\displaystyle{ \sigma }[/math]. The equation can be rearranged to give

[math]\displaystyle{ \Delta G= (N_i - n_x) g_v + n_x g_l + 4 \pi r^2 \sigma - N_i g_v = n_x (g_l - g_v ) + 4 \pi r^2 \sigma . }[/math]

Let [math]\displaystyle{ v_l }[/math] and [math]\displaystyle{ v_v }[/math] be the volume occupied by one molecule in the liquid phase and vapor phase respectively. If the drop is considered to be spherical, then

[math]\displaystyle{ n_x v_l = \frac{4}{3} \pi r^3. }[/math]

The number of molecules in the drop is then given by

[math]\displaystyle{ n_x = \frac{4 \pi r^3 }{3 v_l}. }[/math]

The change in Gibbs energy is then

[math]\displaystyle{ \Delta G = \frac{4 \pi r^3 }{3 v_l} (g_l - g_v) + 4 \pi r^2 \sigma . }[/math]

The differential form of the Gibbs free energy of one molecule at constant temperature and constant number of molecules can be given by:

[math]\displaystyle{ dg = (v_l - v_v ) dP. }[/math]

If we assume that [math]\displaystyle{ v_v \gg v_l }[/math] then

[math]\displaystyle{ dg \simeq - v_v dP. }[/math]

The vapor phase is also assumed to behave like an ideal gas, so

[math]\displaystyle{ v_v = \frac{k T}{P}, }[/math]

where [math]\displaystyle{ k }[/math] is the Boltzmann constant. Thus, the change in the Gibbs free energy for one molecule is

[math]\displaystyle{ \Delta g = - k T \int\limits_{P_{sat}}^{P} \frac {dP}{P}, }[/math]

where [math]\displaystyle{ P_{sat} }[/math] is the saturated vapor pressure of [math]\displaystyle{ x }[/math] over a flat surface and [math]\displaystyle{ P }[/math] is the actual vapor pressure over the liquid. Solving the integral, we have

[math]\displaystyle{ \Delta g = g_l - g_v = -k T \ln \Bigl( \frac{P}{P_{sat}}\Bigr). }[/math]

The change in the Gibbs free energy following the formation of the drop is then

[math]\displaystyle{ \Delta G =- \frac{ 4}{3} \pi r^3 \frac{k T}{v_l} \ln \Bigl(\frac{P}{P_{sat}} \Bigr) + 4 \pi r^2 \sigma . }[/math]

The derivative of this equation with respect to [math]\displaystyle{ r }[/math] is

[math]\displaystyle{ \frac{\partial\bigl( \Delta G\bigr)}{\partial r} = -4 \pi r^2 \frac{ k T}{v_l} \ln \Bigl(\frac{P}{P_{sat}}\Bigr) + 8 \pi r \sigma . }[/math]

The maximum value occurs when the derivative equals zero. The radius corresponding to this value is:

[math]\displaystyle{ r = \frac{2 v_l \sigma }{kT \ln \Bigl(\frac{P}{P_{sat}} \Bigr)}. }[/math]

Rearranging this equation gives the Ostwald–Freundlich form of the Kelvin equation:

[math]\displaystyle{ \ln \Bigl(\frac{P}{P_{sat}} \Bigr) = \frac{2 v_l \sigma}{rkT}. }[/math]

Apparent paradox

An equation similar to that of Kelvin can be derived for the solubility of small particles or droplets in a liquid, by means of the connection between vapour pressure and solubility, thus the Kelvin equation also applies to solids, to slightly soluble liquids, and their solutions if the partial pressure [math]\displaystyle{ p }[/math] is replaced by the solubility of the solid ([math]\displaystyle{ c }[/math]) (or a second liquid) at the given radius, [math]\displaystyle{ r }[/math], and [math]\displaystyle{ p_{\rm sat} }[/math] by the solubility at a plane surface ([math]\displaystyle{ c_{\rm sat} }[/math]). Hence small particles (like small droplets) are more soluble than larger ones. The equation would then be given by:

[math]\displaystyle{ \ln \frac{c}{c_{\rm sat}}= \frac{2 \gamma V_\text{m}}{rRT}. }[/math]

These results led to the problem of how new phases can ever arise from old ones. For example, if a container filled with water vapour at slightly below the saturation pressure is suddenly cooled, perhaps by adiabatic expansion, as in a cloud chamber, the vapour may become supersaturated with respect to liquid water. It is then in a metastable state, and we may expect condensation to take place. A reasonable molecular model of condensation would seem to be that two or three molecules of water vapour come together to form a tiny droplet, and that this nucleus of condensation then grows by accretion, as additional vapour molecules happen to hit it. The Kelvin equation, however, indicates that a tiny droplet like this nucleus, being only a few ångströms in diameter, would have a vapour pressure many times that of the bulk liquid. As far as tiny nuclei are concerned, the vapour would not be supersaturated at all. Such nuclei should immediately re-evaporate, and the emergence of a new phase at the equilibrium pressure, or even moderately above it should be impossible. Hence, the over-saturation must be several times higher than the normal saturation value for spontaneous nucleation to occur.

There are two ways of resolving this paradox. In the first place, we know the statistical basis of the second law of thermodynamics. In any system at equilibrium, there are always fluctuations around the equilibrium condition, and if the system contains few molecules, these fluctuations may be relatively large. There is always a chance that an appropriate fluctuation may lead to the formation of a nucleus of a new phase, even though the tiny nucleus could be called thermodynamically unstable. The chance of a fluctuation is e^−ΔS/k, where ΔS is the deviation of the entropy from the equilibrium value.

It is unlikely, however, that new phases often arise by this fluctuation mechanism and the resultant spontaneous nucleation. Calculations show that the chance, e^−ΔS/k, is usually too small. It is more likely that tiny dust particles act as nuclei in supersaturated vapours or solutions. In the cloud chamber, it is the clusters of ions caused by a passing high-energy particle that acts as nucleation centers. Actually, vapours seem to be much less finicky than solutions about the sort of nuclei required. This is because a liquid will condense on almost any surface, but crystallization requires the presence of crystal faces of the proper kind.

For a sessile drop residing on a solid surface, the Kelvin equation is modified near the contact line, due to intermolecular interactions between the liquid drop and the solid surface. This extended Kelvin equation is given by

[math]\displaystyle{ \ln \frac{c}{c_{\rm sat}}= \frac{V_\text{m}}{RT} \left(\frac{2 \gamma}{r} + \Pi\right). }[/math]

where [math]\displaystyle{ \Pi }[/math] is the disjoining pressure that accounts for the intermolecular interactions between the sessile drop and the solid and [math]\displaystyle{ \left(2 \gamma/r \right) }[/math] is the Laplace pressure, accounting for the curvature-induced pressure inside the liquid drop. When the interactions are attractive in nature, the disjoining pressure, [math]\displaystyle{ \Pi }[/math] is negative. Near the contact line, the disjoining pressure dominates over the Laplace pressure, implying that the solubility, [math]\displaystyle{ c }[/math] is less than [math]\displaystyle{ c_{\rm sat} }[/math]. This implies that a new phase can spontaneously grow on a solid surface, even under saturation conditions.

See also

• Condensation
• Gibbs–Thomson equation
• Ostwald–Freundlich equation

Further reading

• W. J. Moore, Physical Chemistry, 4th ed., Prentice Hall, Englewood Cliffs, N. J., (1962) p. 734–736.
• S. J. Gregg and K. S. W. Sing, Adsorption, Surface Area and Porosity, 2nd edition, Academic Press, New York, (1982) p. 121.
• Arthur W. Adamson and Alice P. Gast, Physical Chemistry of Surfaces, 6th edition, Wiley-Blackwell (1997) p. 54.
• Butt, Hans-Jürgen, Karlheinz Graf, and Michael Kappl. "The Kelvin Equation". Physics and Chemistry of Interfaces. Weinheim: Wiley-VCH, 2006. 16–19. Print.
• Anton A. Valeev,"Simple Kelvin Equation Applicable in the Critical Point Vicinity",European Journal of Natural History, (2014), Issue 5, p. 13-14.