# Vapour Pressure Equations for the Platinum Group Elements

## IMPROVED CALCULATIONS OF VAPOUR PRESSURE AND TEMPERATURE

By J. W. Arblastercor1
Coleshill Laboratories, Gorsey Lane, Coleshill, West Midlands B46 1JU, U.K.

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:

• (a)  The standard value of the atomic weight of platinum has been changed from 195.078 ± 0.002 to 195.084 ± 0.009 (8). In the review on platinum (1) the thermodynamic properties of the condensed phases were based on an atomic weight of 195.08, and the gas phase properties had been corrected using the atomic weight of 195.078. Since the change of 0.006 to the new atomic weight was from one single determination to another, it is likely that a definitive value for the atomic weight has not yet been obtained. Therefore, for consistency with the treatment of the condensed phases, the values for the gas phase have been recalculated using an atomic weight of 195.08. The effect on the tabulated values is negligible at the level of accuracy given.
• (b)  Two further determinations of the vapour pressure of palladium by Zaitsev, Priselkov and Nesmeyanov (9) (1267 to 1598 K) and Ferguson, Gardner and Nuth (10) (1473 to 1973 K) led to heats of sublimation of 377.3 ± 0.1 kJ mol−1 and 377.7 ± 0.2 kJ mol−1, respectively, in excellent agreement with the selected value of 377 ± 5 kJ mol−1.
• (c)  A number of amendments to the calculation of the thermodynamic properties of the gaseous elements from energy levels have been found to have a negligible effect on the free energy functions at the level of accuracy given in the tables in References (16):
• (i)  For the metals other than platinum the thermodynamic properties were calculated using mainly the 1986 CODATA Fundamental Constants (11) but all values have now been recalculated using the 2006 CODATA Fundamental Constants (12).
• (ii)  Joueizadeh and Johansson (13) have revised seventeen of the energy levels of ruthenium.
• (iii)  Engleman et al. (14) have completely revised the energy levels of palladium and extended their number.

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.

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.

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

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.

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