Thermodynamics for Mathematicians (IV)

Note: for this installation, we now require multi-variate calculus, specifically partial differentiation.

10. More State Parameters

For this section, let’s consider the case of a homogeneous gas, whose state is parametrised by (PVN), where N is kept constant for now. Recall that the first law of thermodynamics says: dQdUP·dV, where dQ is a small amount of heat added to the system and dU is a small change in the system’s internal energy. But we also know from definition that dQ/T = dS, which is a small change in the system’s entropy. This gives:

 dUT·dS – P·dV.

This suggests that if we’re looking at the state variable U, then it’s a good idea to re-parametrise it as UU(SVN). The above equation then tells us:

\left.\frac{\partial U}{\partial S}\right|_{V,N} = T,\quad \left.\frac{\partial U}{\partial V}\right|_{S,N} = -P.

 [ For students who’re unfamiliar with partial differentiation, this means: if we keep the parameters V and N constant and differentiate U with respect to S, then the result is T. E.g. if f(x,y,z) = (y+z)x^2 + \sin(y)x + e^z, then \left.\frac{\partial f}{\partial x}\right|_{y,z} = 2(y+z)x + \sin(y). ]

What if we wish to parametrise in terms of TV and N? The best approach is then to define FU – ST, the Helmholtz free energy. If we perturb the system by a little bit, then the change in F is:

dF = dU – S·dT – T·dS = –S·dT – P·dV.

If we parametrise F via F(TVN), then \left.\frac{\partial F}{\partial T}\right|_{V,N} = -S and \left.\frac{\partial F}{\partial V}\right|_{T,N} = -P. But is there any physical meaning to F other than a mere symbolic convenience? Well, yes. There’re at least two ways of looking at it.

  • It’s the maximum amount of useful work you can get out of a system at an isothermal setting. Imagine connecting the system to a heat source of temperature T (i.e. same temperature as the system) and letting it do work by changing P and V. This process is reversible. Let be the amount of heat the system extracts from the heat source.
    • Since T = constant, change in entropy S_2 - S_1 = Q/T.
    • Change in internal energy is U_2 - U_1.
    • Thus the work done is Q - (U_2 - U_1) = -(E_2 - S_2 T)+(E_1 - S_1 T), which is the negation of the change in free energy. Hence one can envisage converting the free energy into work.
  • It’s the amount of energy required to form the system when it’s connected to a heat source T. Without the heat source, the amount of energy required is clearly U. However, the heat source helps out by contributing ST.

To parametrise in terms of SP and N, we let HU+PV be the enthalpy of the system. Then dHdUP·dVV·dP = T·dS + V·dP and we get:

\left.\frac{\partial H}{\partial S}\right|_{P,N} = T, \quad \left.\frac{\partial H}{\partial P}\right|_{S,N} = V.

This is the amount of energy required to create the system when it’s fixed at constant pressure P (called an isobaric setting): in addition to the internal energy U, the system has to do work in order to maintain its pressure. Since P is constant, the work done is PV.

Finally, as the reader might have guessed, for (T, PN), we define the Gibbs free energy GU + PV – ST.  One obtains:

\left.\frac{\partial G}{\partial T}\right|_{P,N} = -S, \quad \left.\frac{\partial G}{\partial P}\right|_{T,N} = V.

This is the amount of energy required to create the system in an isobaric and isothermal setting (i.e. its pressure is kept constant and at the same time, it’s connected to a constant heat source).

11. Intensive and Extensive Variables

Let’s summarise the large number of state parameters we have:

  • pressure P,
  • volume V,
  • temperature T,
  • entropy S,
  • number of particles N,
  • internal energy U,
  • Helmholtz free energy F,
  • Gibbs free energy G,
  • enthalpy H.

Each of these state parameters can be classified from the following question: if we take the union of two identical systems, do we get the same parameter, or twice the parameter? If it’s the former, we say the parameter is intensive; for the latter, extensive. In this list, the only intensive state variables are P and T while the rest are extensive. Thus, for the internal energy U, we have:

U(\lambda S, \lambda V, \lambda N) = \lambda \cdot U(S, V, N)

for any positive value of λ. We also get:

  • F(T, \lambda V, \lambda N) = \lambda \cdot F(T, V, N);
  • G(T, P, \lambda N) = \lambda \cdot G(T, P, N);
  • H(\lambda S, P, \lambda N) = \lambda \cdot H(S, P, N).

The second relation is particularly interesting: if T and P are kept constant, then G is linear in N! Hence we can write G(T, P, N) = \mu(T, P)N for some positive \mu(T, P). This is called the chemical potential of the system. This parameter is clearly intensive.

Example 9. Let’s go back to the problem of finding the entropy of an ideal gas. We had obtained:

S(P, V, N) = N(c\log P+(c+1)\log(V)) + t(N).

To figure out the right t(N), we use the fact that SVN are extensive while P is intensive, i.e. S(P, \lambda V, \lambda N) = \lambda S(P, V, N). Simplifying this relation, we obtain:

\lambda N \log\lambda + t(\lambda N) - \lambda\cdot t(N) \equiv 0.

Replacing \lambda N and N by N and N0 respectively, we see that t(N) is of the form: t(N) = -N\log N + k N for some constant k. Hence, the formula for entropy is now:

S(P, V, N) = N(c\log P+(c+1)\log(V/N)) + kN.

We will leave it at that, since figuring out k is rather tricky business which requires quantum mechanics. The complete version of the equation, known as the Sackur-Tetrode equation, even depends on the mass of each particle!

12. Application of Thermodynamics: Introduction to Phase Transition

Phase transition occurs when the system experiences a sudden change in its physical properties. Typical examples that we see in everyday life include boiling / condensation and freezing / melting. The theory of thermodynamics allows us to model the behaviour of matter, if not quantitatively, then at least qualitatively.

The following is extracted and rephrased from David Tong’s notes on statistical physics, but any mistake most likely belongs to me. Most of the diagrams below are also stolen / modified from his notes.

Let’s write down van der Waal’s equation for a not-so-ideal gas/liquid:

P = \frac{T}{v-b} - \frac a {v^2}, where vV/N.

If we plot out the Pv graph (warning: not PV graph) for various temperatures, we obtain the following:

Note that if V is close to b, then the material is in a liquid state since compressing it requires a disproportionately large amount of force. Conversely, if V is huge, then it is in gaseous state.

For temperatures below a critical limit T < T_c the graph has two turning points. The state (PV) of the material is unstable between these two points. For if we reduce the volume a little by compressing it, then the decrease in pressure causes it to be compressed further; conversely, if we expand it a little, then the increase in pressure will continue to expand it. Indeed, within this region the material is in a mixture of liquid and gaseous state (we say it is in vapour-liquid equilibrium for the pure system).

How do we find the interval where vapour-liquid equilibrium occurs?

In order for vapour and liquid to coexist, the intensive parameters of both components must be equal. The temperature is already equal. Let the common pressure be P. The final parameter is the chemical potential μ (or the Gibbs free energy per particle), which is the Gibbs energy per particle since G(T, P, N) = μ(TP) N. Thus we have:

\left.\frac{\partial\mu}{\partial P}\right|_{T} = \frac 1 N \left.\frac{\partial G}{\partial P}\right|_{T,N} = \frac V N.

Since our T and N are indeed fixed, we get \mu(P, T) = \mu_0 + \int \frac{V(P',T)} N dP' (here, P doesn’t increase monotonically over the integration region!). Recall that v = V/N. The result of the integral is hence the difference of the areas of the two shaded regions below. [ Note: this is non-trivial, do take a minute to figure out what the integration means over a non-monotonically increasing variable. ]

Hence, for (Pμ) of both the solid and liquid states to be equal, we need \int \frac V N dP' = 0, i.e. the two shaded regions have the same area. This is known as Maxwell construction, which gives us the condition for which the liquid and gaseous state can coexist. Note that the region is wider than the interval between the two turning points. The remaining intervals correspond to meta-stable states; these states can exist if we prepare the material sufficiently slowly and carefully but they’re highly delicate and unstable.

Now let the temperature increase until T \to T_c. This approaches the critical point, beyond which there is no clear distinction between the liquid and gaseous states.

As one might imagine, this temperature is rather high (for water, it’s about 374°C). From van der Waal’s equation, this critical point can be determined via algebra, for it is where the Pv graph has an inflection point:

\frac{\partial P}{\partial v} = 0, \quad \frac{\partial^2 P}{\partial v^2} = 0.

Solving gives T_c = 8a/27b. The corresponding pressure and volume at the critical point are P = a/27b^2 and v = V/N = 3b.

And this concludes our very brief introduction to phase transition. Needless to say, we’ve barely scratched the surface of the applications of thermodynamics, but we hope the reader’s interest is sufficiently piqued to seek out other more definitive sources.

This entry was posted in Notes and tagged , , , , , . Bookmark the permalink.

Leave a Reply

Fill in your details below or click an icon to log in: Logo

You are commenting using your account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s