CHEM-UA 652: Thermodynamics and Kinetics

Size: px

Start display at page:

Download "CHEM-UA 652: Thermodynamics and Kinetics"

Antony Hampton
6 years ago
Views:

1 1 CHEM-UA 652: hermodynamics and Kinetics Notes for Lecture 14 I. HE CLAEYRON EQUAION he Clapeyron attempts to answer the question of what the shape of a two-phase coexistence line is. In the - plane, we see the a function (, which gives us the dependence of on along a coexistence curve. Consider two phases, denoted α and β, in equilibrium with each other. hese could be solid and liquid, liquid and gas, solid and gas, two solid phases, et. Let µ α (, and µ β (, be the chemical potentials of the two phases. We have just seen that µ α (, = µ β (, (1 Next, suppose that the pressure and temperature are changed by d and d. he changes in the chemical potentials of each phase are ( µα d + dµ α (, = dµ β (, ( µα d = ( µβ d + ( µβ d (2 However, since G(n,, = nµ(,, the molar free energy Ḡ(,, which is G(n,,/n, is also just equal to the chemical potential Moreover, the derivatives of Ḡ are Ḡ(, = G(n,, n ( Ḡ = V, Applying these results to the chemical potential condition in Eq. (2, we obtain ( ( ( ( Ḡ α Ḡ α Ḡ β Ḡ β d + d = d + Dividing through by d, we obtain = µ(, (3 ( Ḡ = S (4 V α d S α d = V β d S β d (5 V α S α = V β S β ( Vα V β = S α S β d d = S α S β V α V β (6 he importance of the quantity d/d is that is represents the slope of the coexistence curve on the phase diagram between the two phases. Now, in equilibrium dg = 0, and since G = H S, it follows that dh = ds at fixed. In the narrow temperature range in which the two phases are in equilibrium, we can assume that H is independent of, hence, we can write S = H/. Consequently, we can write the molar entropy difference as d S α S β = H α H β (7

2 2 and the pressure derivative d/d becomes d d = H α H β ( Vα V β = αβ H αβ V (8 a result known as the Clapeyron equation. Eq. (8 tells us that the slope of the coexistence curve is related to the ratio of the molar enthalpy between the phases to the change in the molar volume between the phases. If the phase equilibrium is between the solid and liquid phases, then αβ H and αβ V are Hfus and V fus, respectively. If the phase equilibrium is between the liquid and gas phases, then αβ H and αβ V are Hvap and V vap, respectively. For the liquid-gas equilibrium, some interesting approximations can be made in the use of the Clapeyron equation. For this equilibrium, Eq. (8 becomes In this case, V g V l, and we can approximate Eq. (9 as d d = H vap ( Vg V l (9 d d H vap V g (10 Suppose that we can treat the vapor phase as an ideal gas. Certainly, this is not a good approximation so close to the vaporization point, but it leads to an example we can integrate. Since V g = nr, V g = R, Eq. (10 becomes d d = H vap R 2 1 d d = H vap R 2 dln d = H vap R 2 (11 which is called the Clausius-Clapeyron equation. We now integrate both sides, which yields ln = H vap R +C where C is a constant of integration. Exponentiating both sides, we find ( = C e H vap/r which actually has the wrong curvature for large, but since the liquid-vapor coexistence line terminates in a critical point, as long as is not too large, the approximation leading to the above expression is not that bad. If we, instead, integrate both sides, the left from 1 to 2, and the right from 1 to 2, we find 2 2 dln = H vap 1 1 R 2 d ( 2 ln = H ( vap R 2 1 = H ( vap 2 1 R 1 2 (12 assuming that H vap is independent of. Here 1 is the pressure of the liquid phase, and 2 is the pressure of the vapor phase. Suppose we know 2 at a temperature 2, and we want to know 3 at another temperature 3. he above result can be written as ln ( 3 1 = H vap R ( (13

3 3 Subtracting the two results, we obtain ln ( 2 3 = H vap R ( so that we can determine the vapor pressure at any temperature if it is known as one temperature. (14 In order to illustrate the use of this result, consider the following example: roblem: At 1 bar, the boiling point of water is 373 K. At what pressure does water boil at 473 K. ake the heat of vaporization of water to be kj/mol. Solution: Let 1 = 1 bar and 1 = 373 K. ake 2 = 473 K, and we need to calculate 2. Substituting in the numbers, we find ( (40.65 kj/mol(1000 J/kJ 1 ln 2 (bar = J/mol K 473 K 1 = K 2 (bar = (1 bare 2.77 = 16 bar II. GIBBS HASE RULE In any one-component system with three phases, there is only one point at which all three phases coexist in equilibrium. his is known as the triple point. Although the Gibbs free energy G is a function of n,, and, G = G(n,,, if we allow ourselves to treat and as experimental control parameters, as is done when we construct a phase diagram, then n is determined by the equation of state. and are both intensive properties. Now suppose we have a multicomponent system. When we have more than one component, we have additional intensive properties that could potentially be additional degrees of freedom that determine the dimensionality of the phase diagram. If the phase diagram has a dimensionality higher than that of a plane (two dimensions, then triple points could become lines, planes, curved surfaces, or hypersurfaces. hus, it is important to have a way of determining how many actual experimental control variables we have in a multicommponent system. his number is determined by the number of thermodynamic conditions that exist among the intensive properties. We will now derive a relation that accounts for these conditions and yields the total number of degrees of freedom or experimental control parameters this relation is known as the Gibbs phase rule. Suppose a system has k components and can exist in p phases. An example would be a three-component system consisting of H 2 O, CH 3 OH, and CH 3 CH 2 OH. he Gibbs free energy G is now a function of n 1,...,n k,, and : G = G(n 1,...,n k,,. he total number of moles n is still determined by the equation of state if we allow and to be experimental control parameters. However, we now have a condition n 1 +n 2 + n k = n (15 We can now introduce additional intensive parameters by dividing both sides by n, which yields We defined n 1 n + n 2 nk + n n = 1 (16 X (i = n i n called the mole fraction of component i. Clearly, we must have Moreover, X (i is an intensive property. X (1 +X (2 + X (k = 1 Suppose, for starters, that k = 1, but we have p phases. If α and β denote a pair of phases, then at a p-phase coexistence point, we have a set of conditions from the chemical potential equalization that take the form (17 µ α (, = µ β (, (18

4 4 For example, if we had three phases, solid, liquid, and gas, the conditions would be µ liq (, = µ gas (,, µ liq (, = µ solid (,, µ solid (, = µ gas (, (19 However, by the transitive property, the third condition follows from the first two, so it s not an independent condition. In fact, we have only two independent conditions. In general, if we have p phases, we have p 1 such conditions. Now if we have k commponents, then we have p 1 such conditions for each phase, giving a total of k(p 1 total conditions from chemical potential equalization. he fact that the mole fractions in each phase must sum to 1 gives us additional condition. Let X α (i be the mole fraction of component i in phase α. For each phase, we have a condition X (1 α +X(2 α + X(k α = 1 (20 his gives us p additional conditions. Hence, the total number of conditions we have is k(p 1+p. Although the equation of state constitutes, in principle, another condition, we have already taken it into account is positing that the total number of moles n is determined if we know and. Now the total number of mole fractions is k for each phase for a total of kp mole fractions. In addition, if we add and as additional parameters, then the total number of intensive thermodynamic parameters we can tune is kp+2. If we subtract from this the number of conditions (also the number of constraints, we obtain the number of degrees of freedom d as he condition is known as the Gibbs phase rule. d = kp+2 k(p 1 p = kp+2 kp+k p = k p+2 d = k p+2 (21 As an example, suppose we have k = 3 components and p = 3 phases. hen the number of degrees of freedom is d = = 2 which we could take to be and. However, we could also take these to be and one of the mole fractions or and one of the mole fractions. As another example, consider water, which has three phases, ice, water, steam. Here, we have k = 1 and p = 3. hus, d = = 0 which means that when all three phases are in coexistence, the number of degrees of freedom is 0. Hence, there can be only one point, i.e. one value of and one value of, at which this coexistence occurs, which we call the triple point. However, this explains why the triple point is a point and now a structure of higher dimension. III. CLASSICAL NUCLEAION HEORY Having looked at phase diagrams in some detail, we now turn our attention to the process by which phase transitions occur. In particular, we will consider a phenomenological theory known as classical nucleation theory. Let us consider it in the special case of the solid-liquid transition, i.e., the melting transition. In this case, if we are in the solid phase and far from the melting point, then we have the condition µ s < µ l, which tells us that the solid phase is more stable (remember that chemical potentials are also molar free energies, and the solid phase should have a lower free energy far from the melting point. As we approach the melting point, the liquid phase becomes the more stable of the two phases, and we have µ s > µ l or µ l µ s < 0. he solid phase now becomes a metastable phase. his means that, even though the liquid phase

is the more stable of the two phases, there is a significant free energy barrier that must be crossed in order for the system to undergo a transition from solid to liquid.

5 is the more stable of the two phases, there is a significant free energy barrier that must be crossed in order for the system to undergo a transition from solid to liquid. How do we characterize this barrier? In order to do this, we introduce a microscopic coordinate in the system that is capable of following the progress of the melting process. Such a coordinate is called a reaction coordinate. In general, a reaction coordinate is some function of the atomic positions that is capable of following a particular process, such as chemical reaction, a diffusion process, or a phase transition. Classical nucleation theory (CN posits that as a solid is heated close to its melting point, a liquid nucleus can form in the solid, and as this nucleus grows, eventually, the system will melt completely, transforming from solid to liquid. If we assume that this nucleus is a sphere of radius r, then r can be viewed as the reaction coordinate. As r increases, the size of the liquid nucleus increases, and ultimately, the system melts. his process is illustrated in Fig. 1 below. he free energy to form a nucleus of radius r will be a function of r, which we denote as G(r. Such 5 FIG. 1: Illustration of a liquid nucleus. a coordinate-dependent free energy is called a free energy profile. he free energy needed to form such a nucleus has two contributions. First, work is needed to form the surface of the nucleus, which is an interface between the solid and the liquid. As these two things do not naturally fit together, it is not energetically favorable to form the surface of the nucleus. Let γ s denote the surface tension of the nucleus, which has units of energy per unit area. he area of the nucleus if 4πr 2, so this contribution to the free energy is 4πr 2 γ s. On the other hand, since µ l < µ s, forming a liquid nucleus of volume 4πr 3 /3 is favorable. If µ represents the chemical potential difference per unit volume, then given that µ < 0, this contribution to the free energy (4πr 3 /3 µ is negative. Hence, we have a balance between the volume and surface area terms, and the total free energy in the CN is G(r = 4 3 πr3 µ+4πr 2 γ s (22 If we plot this free energy as a function of r, we obtain a curve of the form show in Fig. 2. We see from the free energy profile that there is a significant barrier at a particlar radius r max. he height of the barrier is G(r max. o find what this radius is, we need to maximize G. We take the derivative of G(r with respect to r and set it equal to 0 d G dr = 4πr 2 µ+8πrγ s = 0 r µ+2γ s = 0 If we evaluate the free energy at this radius, we find the barrier height to be G(r max = 4 ( 3 π 2γ 3 s +4π( 2γ 2 s γ s µ µ r = 2γ s µ r max (23

6 6 G(r max r max FIG. 2: lot of G(r. = 4 ( ( 8γ 3 3 π s 4γ 2 ( µ 3 µ+4π s ( µ 2 γ s = 16πγ3 s 3( µ 2 (24 For water, for example, the height of the barrier, computed from this formula, would be around 100 k B, which corresponds to a barrier of around 50 kcal/mol. hus, we see that, although the liquid phase is more stable, a very high energy barrier needs to be surmounted, which is a rare event.

Exam 3, Chemistry 481, 8 December 2017

1 Exam 3, Chemistry 481, 8 December 2017 Show all work for full credit Useful constants: k B = 1.3807 10 23 J K 1 ; R (molar gas constant) = 8.314 J K 1 mol 1 Helmholz free energy: A = U S, so that da