Phase Transitions and Phase Diagrams
One-component systems Enthalpy and entropy dependence on P and T Gibbs free energy dependence on P and T Clapeyron equation Understanding phase diagrams for one-component systems Polymorphic phase transitions Driving force for a phase transition First order and second-order phase transitions
Reading: 1.2 of Porter and Easterling Chapter 7.1 7.4 of Gaskell
MSE 3050, Phase Diagrams and Kinetics, Leonid Zhigilei
PVT Surface of a Pure Substance
MSE 3050, Phase Diagrams and Kinetics, Leonid Zhigilei http://www.eng.usf.edu/~campbell/ThermoI/ThermoI_mod.html
A pure substance is heated at constant pressure
Tb
MSE 3050, Phase Diagrams and Kinetics, Leonid Zhigilei
H and S as function of T at constant P
In a closed one-component system equilibrium, at temperature T and pressure P, corresponds to the state with minimum Gibbs free energy G. Therefore, in order to predict what phases are stable under different conditions we have to examine the dependence of G on T and P. Lets use thermodynamic relations to predict the temperature dependence of H, S, and G at constant P.
H For H(T) we have = CP T P
H (T ) = H 298 +
298
C dT
P
For S(T) we have
C S = P T P T
S(T ) =
0
CP dT T
CP
0 H 0 298
T, K
Slope = CP
T, K
S
Slope = CP/T
0
MSE 3050, Phase Diagrams and Kinetics, Leonid Zhigilei
T, K
G as function of T at constant P
For G = H TS we have dG = -SdT +VdP and for P = const
G = S T P
for the slope
2G cP S = = T 2 T T P P
for the curvature
H Slope = CP 0 TS G G(T) for a single phase at P = const T, K Slope = -S
MSE 3050, Phase Diagrams and Kinetics, Leonid Zhigilei
G as function of T at constant P for liquid and solid phases
At all temperatures the liquid has a higher internal energy U and enthalpy H as compared to the solid. Therefore Gl > Gs at low T. The liquid phase, however, has a higher entropy S than the solid phase at all T. Therefore Gl decreases more rapidly with T as compared to Gs. At Tm Gl(T) crosses Gs(T) and both liquid and solid phases can co-exist in equilibrium (Gl = Gs)
Hl Hm 0 Tm
Hs T, K
Hl > Hs Sl > Ss at all T
Hm = TmSm Gs Gl
At Tm the heat supplied to the system will not rise its temperature but will be used to supply the latent heat of melting Hm that is required to convert solid into liquid. At Tm the heat capacity Cp = (H/T)P is infinite addition of heat does not increase T.
MSE 3050, Phase Diagrams and Kinetics, Leonid Zhigilei
A typical P-T phase diagram for a pure material
P = 1 atm
vapor P
solid
normal freezing point 1 atm
critical point
liquid solid
liquid
gas
triple point
normal boiling point
T
solid is stable liquid is vapor is stable stable
The red lines on the phase diagram show the conditions where different phases coexist in equilibrium: Gphase1 = Gphase2
MSE 3050, Phase Diagrams and Kinetics, Leonid Zhigilei
G as function of P at constant T for liquid and solid phases
As we can see from the fundamental equation, dG = VdP SdT, the free energy of a phase increases with pressure: G =V>0 P T If the two phases have different molar volumes, their free energies will increase by different amounts when pressure changes at a fixed T.
i ce
H 2O d i u liq
T = 0C
Vl < Vs for water Vl > Vs for most materials
1 atm
How the unusual change of V upon melting of water could be related to ice-skating? What is the curvature of the G(P) at constant T?
kT =
1 V - isothermal V P T compressibility
P B = V V T
- bulk modulus
MSE 3050, Phase Diagrams and Kinetics, Leonid Zhigilei
Equilibrium between two phases: Clapeyron equation
If two phases in equilibrium have different molar volumes, their free energies will increase by different amounts when pressure changes at a fixed T. The equilibrium, therefore will be disturbed by the change in pressure. The only way to maintain equilibrium at different pressures is to change temperature as well. For two phases in equilibrium Gl = Gs and dGl = dGs for infinitesimal change in T and P (so that the system remains in equilibrium)
dG l = Vl dP - Sl dT
dG s = Vs dP - Ss dT
Vl dP - Sl dT = Vs dP - Ss dT
Ss Sl S dP = = dT eq. Vs Vl V
At equilibrium
G = H T S = 0
and
H = T S
Therefore
H dP = dT eq. T V
- the Clapeyron equation
The Clapeyron equation gives the relationship between the variations of pressure and temperature required for maintaining equilibrium between the two phases.
MSE 3050, Phase Diagrams and Kinetics, Leonid Zhigilei
G as function of P and T for liquid and solid phases
Schematic representation of the equilibrium surfaces of the solid and liquid phases of water in G-T-P space. The planes show the free energies of liquid and solid phases, the intersections of the planes correspond to the (P, T) conditions needed for maintaining equilibrium between the phases, Gl = Gs.
solid
liquid T P 1 atm 0C
MSE 3050, Phase Diagrams and Kinetics, Leonid Zhigilei
Clapeyron equation: examples
A typical diagram for a pure material: For liquid to gas transition: V = Vg - Vl >> 0 H = Hg - Hl > 0 we have to add heat to convert liquid to gas. Therefore
solid P
liquid
gas T
H dP = >0 dT eq. T V
For liquid to solid transition: V = Vs Vl < 0 for most materials H = Hs - Hl < 0 heat is released upon crystallization. Therefore
H dP >0 = dT T V eq.
dP <0 For some materials, however, V = Vs Vl > 0 and dT eq. This is the case for H2O (ice floats in liquid water), Ga, Ge, diamond
In general, it takes heat (H > 0) to proceed from a low-temperature to a high-temperature phase (entropy of a high-temperature phase is higher than the entropy of a low-temperature phase). Therefore, the slope of the equilibrium lines in a P-T phase diagram of a pure material reflects the relative densities of the two phases.
MSE 3050, Phase Diagrams and Kinetics, Leonid Zhigilei
Clapeyron equation: more examples
Some materials may exist in more than one crystal structure, this is called polymorphism. If the material is an elemental solid, it is called allotropy. The diagram shows the stable phases for pure iron at varying temperature and pressure. Close-packed FCC -Fe has a smaller molar volume than BCC -Fe: V = V - V < 0 At the same time H = H - H > 0 Therefore -Fe
-Fe
H dP = <0 dT eq. T V
-Fe
The effect of increasing pressure is to increase the area of the phase diagram over which the phases of the smaller molar volume (higher density) is stable.
FCC G
BCC
liquid
P = const
vapor
T
MSE 3050, Phase Diagrams and Kinetics, Leonid Zhigilei
H dP <0 = dT T V eq.
-Fe
-Fe
-Fe
MSE 3050, Phase Diagrams and Kinetics, Leonid Zhigilei
Phase diagrams of one-component materials
Another example of allotropy is carbon, which can exist as diamond, graphite, and amorphous carbon.
From the phase diagram, what phase has higher density near two-phase equilibrium conditions, Diamond or graphite? Diamond or liquid carbon? Graphite or liquid carbon?
MSE 3050, Phase Diagrams and Kinetics, Leonid Zhigilei
H dP <0 = dT T V eq.
MSE 3050, Phase Diagrams and Kinetics, Leonid Zhigilei
Can the volume expansion of water upon melting explain ice-skating? m = 100 kg
H dP = <0 dT eq. T V
Vice = 19.63 cm3 /mol Vliquid = 18.0 cm3 /mol H m = 5636 J/mol skate : 2 mm wide, 1 cm contact length
MSE 3050, Phase Diagrams and Kinetics, Leonid Zhigilei
Example: squeezing diamond from graphite
What pressure should we apply to transform graphite to diamond at 298 K? In a reference book we can find that at 298 K and 1 atm: Hdiamond = 1900 J/mol Sgraphite = 5.73 J/K Sdiamond = 2.43 J/K graphite = 2.22 g/cm3 diamond = 3.515 g/cm3
MSE 3050, Phase Diagrams and Kinetics, Leonid Zhigilei
The driving force for the phase transformation
If solid and liquid are in equilibrium, Gs = Gl and a slow addition of heat leads to the melting of some part of the solid, but do not change the total G of the system: G = nl Gl + ns Gs = const, where nl and ns are the numbers of moles of liquid and solid phases, and Gl and Gs are the molar Gibbs free energies. If energy is added/removed quickly, the system can be brought out of equilibrium (overheated or undercooled) the melting/freezing process is spontaneous/irreversible and G is decreasing.
At temperature T*
G G T T* Tm Gs Gl
G l = H l - T*Sl
G s = H s - T*Ss
G = H - T*S
At temperature Tm
G = H - Tm S = 0
S = H m Tm
For small undercooling T we can neglect the difference in Cp of liquid and solid phases and assume that H and S are independent of temperature. The driving force for solidification
MSE 3050, Phase Diagrams and Kinetics, Leonid Zhigilei
G H m T * G
H m Tm
H m T Tm
First-order and second-order phase transitions (I)
The classification of phase transitions proposed by Ehrenfest is based on the behavior of G near the phase transformation. First-order phase transition: discontinuous. first derivatives of G are
Second-order phase transition: first derivatives of G are continuous, but second derivatives of G are discontinuous. First-order phase transition
Ttrs H
T Cp
Ttrs
Ttrs
G trs = 0 Strs =
H trs Ttrs
Ttrs
G = -S T P
Ttrs
Strs 0 H trs 0
dH CP = dT P
G = V - discontinuous P T
e.g. melting, boiling, sublimation, some polymorphous phase transitions. MSE 3050, Phase Diagrams and Kinetics, Leonid Zhigilei
First-order and second-order phase transitions (II)
Second-order phase transition
Ttrs H
T Cp
Ttrs
Ttrs
G trs = 0 Strs = 0 H trs = 0 Vtrs = 0
Ttrs
G = -S T P
Ttrs
G =V P T
dH CP = dT P
- continuous (S and V do not jump at transition)
2 G V TP = T P
2G S - discontinuous = PT P T
transition in metals at low
e.g. conducting-superconducting temperatures.
MSE 3050, Phase Diagrams and Kinetics, Leonid Zhigilei
Summary Make sure you understand language and concepts: Enthalpy and entropy dependence on P and T Gibbs free energy dependence on P and T Clapeyron equation Understanding phase diagrams for one-component systems Driving force for a phase transition First order and second-order phase transitions
Make sure you understand P-T, G-P, G-T 2D phase diagrams for a one-component system (what is shown, what are the lines separating different regions, how to predict the slopes of the lines, etc.)
MSE 3050, Phase Diagrams and Kinetics, Leonid Zhigilei