Vapour Pressure Equations for the Platinum Group Elements

IMPROVED CALCULATIONS OF VAPOUR PRESSURE AND TEMPERATURE

 

While a knowledge of the vapour pressure curve of any material is of theoretical significance in understanding its basic physical properties, it can also be of practical importance in, for example, the use of the material in high-temperature vacuum applications. Therefore readily usable equations which accurately predict values of vapour pressure over a wide range of temperatures and pressures can have an important practical use. For the platinum group metals (pgms) the vapour pressures can be immediately assessed from about 10−16 bar to just above the boiling point by the use of Equation (i), fitted for the solid and liquid metals separately: ln(p, bar) = A + Bln(T) + C/T + DT + ET2, where p is the vapour pressure, T is the temperature in kelvin and A, B, C, D and E are constants. Although containing five coefficients, this equation can easily be evaluated by computers and scientific calculators. Although it gives values of vapour pressure at fixed temperatures, by a simple and rapid use of iteration values of temperature at fixed vapour pressure, temperature-dependent values of vapour pressure can also be obtained.

 

Vapour pressure values are calculated from thermodynamic data following the procedure described in Appendix A. Since free energy functions are given for one bar standard-state pressure, then the vapour pressure will also be given in bars. Values for the free energy functions and the selected heats of sublimation are given in the reviews on platinum (1) and the other five pgms: palladium (2), ruthenium (3), osmium (4), rhodium (5) and iridium (6). Since the behaviour of the specific heat capacities of the solids and their effect on the free energy functions is more complicated than for the liquids, then vapour pressure values were evaluated for the solids at 25 K intervals and the melting point, and for the liquids at 50 K intervals and the melting point, so as to give approximately equal numbers of data points. In the review on platinum (1), a minimum lower temperature of 1200 K was considered. This gives a vapour pressure of about 10−17 bar, which is considerably below practical measurements. For the other metals, the use of rounded temperatures gives a vapour pressure between 10−16 and 10−17 bar, which was used as the lower bound.

 

The derivation of the selected vapour pressure equation (Equation (i)) is described in Appendix B.

ln(p, bar) = A + Bln(T) + C/T + DT + ET2      (i)

where A, B, C, D and E are constant coefficients. Equation (i) was used by Honig and Kramer (7) to represent the vapour pressures of the elements over a wide range of temperatures and pressures.

 

Since publication of the reviews on platinum and the other metals (16) newer data have become available for consideration:

 

Coefficients corresponding to the selected vapour pressure equation for both the solid and liquid pgms are given in Table I. One criterion to be met is that the calculated vapour pressures should be equal at the melting point, or more correctly the triple point which is the solid-liquid-vapour equilibrium temperature. However, because of the high resistance of the pgms to compression, then at the level of accuracy obtainable the melting points are considered to be equal to the triple points as indicated in Appendix C.

Table I

Coefficients in Equation (i) for the PGMs in the Solid and Liquid States

 

In the review on platinum (1), only the first three terms of Equation (i) were used, which gave an adequate representation of the vapour pressure. However, when all five terms are used to ensure precise reproducibility for all of the metals, the correlation for platinum is extraordinary. The average deviation between the thermodynamically derived vapour pressure and the result from the equation is only ± 0.0004% for both the solid and liquid, and similar results were generally obtained for the other metals as given in Table II. Even in the worst case, for solid rhodium, agreement is still within five significant figures; this tolerance is well beyond that of any practical determination of vapour pressure.

Table II

Percentage Deviations between Thermodynamically Derived Vapour Pressure Curves and Results from Equation (i) for PGMs in the Solid and Liquid States

 

Whereas the equations for the liquids reproduce the thermodynamic boiling points to within 0.01 K, more realistic estimates of the uncertainties in the boiling points can be ascertained from the uncertainties assigned to the enthalpies of sublimation, with boiling point uncertainties rounded to the nearest 10 K (see Table III).

Table III

Uncertainties in the Boiling Points of PGMs

 

The uncertainties given in Table III must be considered to be minimum values, because the uncertainties in the free energy functions are not taken into account. For example, the enthalpy of fusion and thermodynamic properties of liquid osmium are all estimated, and therefore a more realistic estimate of the uncertainty in the boiling point would be ± 100 K.

 

These values are given in Table IV. Since the lowest temperature used corresponds to a vapour pressure between 10−17 and 10−16 bar, then 10−16 bar is used as the lower bound for fixed pressures. Table IV provides a check on values obtained by iteration of the equations. It is an amendment to a previous summary (15), in which the lowest pressure included was 10−12 bar. Vapour pressures corresponding to fixed temperatures are obtained directly from the equations, and a check is provided by Table V, which lists the melting point vapour pressures.

Table IV

Temperatures Corresponding to Fixed Vapour Pressures for the PGMs

 
Table V

Vapour Pressures at the Melting Points for the PGMs

 
 
  • 1 J. W. Arblaster, Platinum Metals Rev., 2005, 49, (3), 141
  • 2 J. W. Arblaster, CALPHAD, 1995, 19, (3), 327
  • 3 J. W. Arblaster, CALPHAD, 1995, 19, (3), 339
  • 4 J. W. Arblaster, CALPHAD, 1995, 19, (3), 349
  • 5 J. W. Arblaster, CALPHAD, 1995, 19, (3), 357
  • 6 J. W. Arblaster, CALPHAD, 1995, 19, (3), 365
  • 7 R. E. Honig and D. A. Kramer, RCA Rev., 1969, 30, 285
  • 8 M. E. Wieser, Pure Appl. Chem., 2006, 78, (11), 2051
  • 9 A. I. Zaitsev, Yu. A. Priselkov and A. N. Nesmeyanov, Teplofiz. Vys. Temp., 1982, 20, (3), 589; Teplofiz. Vys. Temp., 1982, 20, (3), 589; Chem. Abstr., 1982, 97, 98769
  • 10 F. T. Ferguson, K. G. Gardner and J. A. Nuth III, J. Chem. Eng. Data, 2006, 51, (5), 1509
  • 11 E. R. Cohen and B. N. Taylor, ‘The 1986 Adjustment of the Fundamental Physical Constants (Report of the CODATA Task Group on Fundamental Constants)’, CODATA Bull., No. 63, November 1986
  • 12 P. J. Mohr, B. N. Taylor and D. B. Newell, ‘The 2006 CODATA Recommended Values of the Fundamental Physical Constants’, (Web Version 5.0), April 2007 (database developed by J. Baker, M. Douma and S. Kotochigova) Available: http://physics.nist.gov/constants, National Institute of Standards and Technology, Gaithersburg, Maryland, U.S.A.
  • 13 A. Joueizadeh and S. Johansson, International Colloquium on Atomic Spectra and Oscillator Strengths for Astrophysical and Laboratory Plasmas, National Institute of Standards and Technology, Gaithersburg, Maryland, U.S.A., 14th–17th September, 1992, NIST Special Publication 850, April, 1993, p. 131
  • 14 R. Engleman, U. Litzén, H. Lundberg and J.-F. Wyart, Phys. Scr., 1998, 57, (3), 345
  • 15 J. W. Arblaster, Platinum Metals Rev., 1996, 40, (2), 62
  • 16 H. M. Strong and F. P. Bundy, Phys. Rev., 1959, 115, (2), 278
  • 17 R. E. Bedford, G. Bonnier, H. Maas and F. Pavese, Metrologia, 1996, 33, (2), 133
 

Appendix A
Calculation of Vapour Pressure Values

Vapour pressure values may be calculated from thermodynamic data via an inversion of the third law of thermodynamics used to calculate the heat of sublimation (Equation (ii)):

RTln(p) = T[δ − (G°TH°298.15)/T] − ΔH°298.15(III) (ii)

where R is the universal gas constant (8.314472 J mol−1 K−1 (12)); T is the temperature in kelvin; p is the vapour pressure in bar; ΔH°298.15(III) is the third law enthalpy of sublimation at 298.15 K (Equation (iii));

δ − (G°TH°298.15)/T = [− (G°TH°298.15)/T(g)] − [− (G°TH°298.15)/T(s, l)] (iii)

where (G°TH°298.15) is the free energy of the phase relative to that at 298.15 K.


Appendix B
Derivation of the Selected Vapour Pressure Equation

The variation of vapour pressure (p) with temperature (T) can be represented by the Clapeyron equation (Equation (iv)):

dp/dT = ΔH/TΔV (iv)

where ΔH is the latent heat of sublimation (below the melting point) or vaporisation (above the melting point); ΔV is the change in volume during the transition from condensed phase to gas. Since the volume change is very large, the value for the condensed phase can be neglected and the equation becomes Equation (v) or (vi):

dp/dT = ΔHp/RT2 (v)

dln(p)/dT = ΔH/RT2 (vi)

Equation (vi) is then integrated to give the Clausius-Clapeyron equation (Equation (vii)):

ln(p) = A/T + B (vii)

where B is an integration constant, and A = ΔHT/R where ΔHT is the average enthalpy centred on the mid-range temperature; ΔHT can only be taken as constant over a relatively narrow temperature range. The Clausius-Clapeyron equation is used practically to fit experimental vapour pressure measurements, and also to calculate the ‘second law’ heat of sublimation through the relationship (Equation (viii)):

ΔH°298.15(II) = − δ (H°TH°298.15) − RA (viii)

where δ (H°TH°298.15) and (H°TH°298.15) (see Equation (ix)) are enthalpy values relative to the enthalpy at 298.15 K;

δ (H°TH°298.15) = (H°TH°298.15)(g) − (H°TH°298.15)(s, l) (ix)

ΔHT may be taken as constant only for a narrow temperature range. In order to cover a wider temperature range it may be expanded to Equation (x):

ΔHT = ΔH0 + a1T + a2T2 + a3T3 + … (x)

Substituting Equation (x) into the Clapeyron equation (Equation (iv)) gives Equation (xi):

dln(p) = ΔH0dT/RT2 + a1dT/RT + a2dT/R + a3TdT/R (xi)

Rearranging the coefficients and integrating gives Equation (xii):

ln(p) = A + a1ln(T)/R − ΔH0/RT + a2T/R + a3T2/R (xii)

where A is the integration constant. Substituting B = a1/R, C = − ΔH0/R, D = a2/R, E = a3/R gives Equation (i):

ln(p, bar) = A + B ln(T) + C/T + DT + ET2 (i)


Appendix C
Melting Point and Triple Point Differences

Although the pressure at the standard state in thermodynamics is now set at one bar, melting points are still quoted at one atmosphere pressure. This is an arbitrary choice, whereas the triple point is that at which the vapour pressures of the solid, liquid and gas are exactly equal; it is therefore a universal point. The difference between the melting and triple points can be determined from the melting pressure curves, for example Strong and Bundy (16) determined the initial values for these curves to be 138 atm deg−1 for platinum and 160 atm deg−1 for rhodium, equivalent to 0.007 deg atm−1 and 0.006 deg atm−1, respectively. Assuming that the triple point pressures are equal to zero, then this indicates that for platinum, the triple point would be 0.007 K below the melting point, and for rhodium 0.006 K. However, these values are negligible when compared to melting point uncertainties of ± 0.4 K for platinum and ± 3 K for rhodium at the secondary fixed points on the International Temperature Scale of 1990 (ITS-90) (17). On these grounds, triple point corrections may be regarded as meaningless. These examples for platinum and rhodium are considered as representative for all six of the pgms.

 

John W. Arblaster is Chief Chemist, working in metallurgical analysis on a wide range of ferrous and non-ferrous alloys for standards in chemical analysis, at Coleshill Laboratories in the West Midlands of England. He is interested in the history of science and in the evaluation of the thermodynamic and crystallographic properties of the elements.