I'm reading some papers for calculating the vapor pressure of alkali metals as a function of temperature, and I've come across some familiar-looking virial expansions, but when I tried to work out the unit analysis, it was very much unclear how to make everything work out.
Here's the first equation (where they give the form in terms of the physical quantities):
$$
\mathrm{ln}\left(p\right) = I_{c} – \frac{\Delta H^{\circ}_{v}}{RT} + \frac{\left(C^{g}_{p} – C_{p}^{l}\right)}{R}\mathrm{ln}\left(T\right)
$$
This is equation 4 of the following paper:
(1) Alcock, C. B.; Itkin, V. P.; Horrigan, M. K. Canadian
Metallurgical Quarterly 1984, 23, 309–313., (doi: 10.1179/000844384795483058
From the definitions section, $\Delta H_{v}^{\circ}$ is the molar enthalpy of vaporization in $\frac{BTU}{lb\cdot mol}$, $R$ is the gas constant. No units are given for this, but given the other values, we'll call it $\frac{BTU}{lb\cdot mol ^{\circ}R}$. $C_{p}^g$ and $C_{P}^l$ are the isobaric (constant pressure) specific heats, given in $\frac{BTU}{lb\cdot mol^{\circ}\!R}$. So working this out, I'm going to assume that this is supposed to be a unitless quantity somehow. I'm assuming that $I_{C} = \mathrm{ln}\left(p_{0}\right)$, so we get:
$$
\mathrm{ln}\left(\frac{p}{p_{0}}\right) = -\frac{\Delta H_{v}^{\circ}}{RT} + \frac{\left(C_{p}^{g}-C_{p}^{l}\right)}{R}\mathrm{ln}\left(T\right)
$$
Everything there is unitless except for $\mathrm{ln}\left(T\right)$. (Alternately, maybe he is secretly assuming that $p = \frac{p}{torr}$ – the unit analysis is basically the same). What is going on there? The other form I see for this is:
$$
\mathrm{log}\left(p\right) = A_{c} – \frac{B_{c}}{T} + C_{c}\mathrm{log}\left(T\right)
$$
Again, there's a logarithmic temperature in there! And, bizarrely, he defines the following quantities (A is nowhere to be found on this list):
$$B = \text{second virial coefficient, }\frac{ft^3}{mol}$$$$
C = \text{third virial coefficient, }\frac{ft^6}{mol^2}
$$
According to the paper, $_{C}$ is the indicator for "empirical constant". I'm guessing $A$ is supposed to be unitless (don't know how to scale it to change the pressure units, but one thing at a time). Trying to find a second opinion, I went to this paper:
(1) Ewing, C. T.; Chang, D.; Stone, J. P.; Spann, J. R.; Miller, R. R.
Journal of Chemical & Engineering Data 1969, 14, 210–214. doi: http://dx.doi.org/10.1021/worsp).
But they have listed no units, and suffer from the same problems. Here is their version of the expansion:
$$
A + \frac{B}{T} + C\cdot\mathrm{log}\left(T\right) + D\cdot T\cdot10^{-3}
$$
They plug in some values for this:
$$
\mathrm{log}\left(p(atm)\right) = 7.55339-\frac{27007.7}{T} + 0.19699\mathrm{log}(T) – \frac{0.171188\cdot T}{1000}
$$
Is temperature a unitless quantity now? What am I missing here?
Best Answer
I don't think that you can assume that $I_c=\ln (p_0)$
A lot of equations have logs of dimensioned quantities in them. It's usually not hard to get rid of the log:
$$ \mathrm{ln}\left(\frac{p_1}{p_2}\right) = I_{c,1}-I_{c,2} - \left(\frac{\Delta H^{\circ}_{v,1}}{RT_1}-\frac{\Delta H^{\circ}_{v,2}}{RT_2}\right) + \mathrm{ln}\left(\frac{T_1^{\frac{\left(C^{g}_{p,1} - C_{p,1}^{l}\right)}{R}}}{T_2^{\frac{\left(C^{g}_{p,2} - C_{p,2}^{l}\right)}{R}}}\right) $$
It may also be that $I_c$ is just a gas-specific constant, of the form $$I_c=\ln p_0-\frac{\Delta H^{\circ}_{v,0}}{R} -\frac{\left(C^{g}_{p,0} - C_{p}^{l,0}\right)}{RT_0}\mathrm{ln}\left(T_0\right) $$
for some $p_0,T_0$.
You can do the same for the other equations as well. When there's a log of a dimensioned quantity, one of two things is possible:
The latter form arises when one takes an indefinite integral in the last step while solving a differential equation (Which is why subtracting is the right way to fix this, since that corresponds to tacking limits onto the integral)