Thermo
Thermo
Thermo
William Thomson, Irish-Scottish physicist and mathematician (Lord Kelvin; 1824 1907): Absolute temperature scale (1848). One of the most important contributors to the 19th century physical sciences.
Examples:
System: single molecule, solid, gas, liquid, etc.
Surroundings: solvent, container, gas cylinder, etc.
Boundary: Interaction between the system and the
surroundings (forces between atoms and molecules) or
vacuum (no interaction).
Focus on the system and treat the surroundings approximately
Terminology:
Non-isolated system
Isolated system
Open system
Closed system
Adiabatic system
Heterogeneous system
Homogeneous system
6
Phase
Thermodynamic variables
Extensive variable
Intensive variable
Intensive state
Extensive state
Equilibrium
7
Molar volume is defined as the ratio between volume V (m3 or dm3 = L) and the
number of particles in it n (mol):
V
(1.1)
V =
n
The amont of substance n (mol) is given by:
N
N
= moles of substance
=
n=
NA
6.022 1023 mol1
(1.2)
Thermal equilibrium between two systems means that they have been in thermal
contact for sufficiently long time so that they have the same temperature. The
following transitivity rule is called the zeroth law of thermodynamics:
According to our every day knowledge of thermal objects, this law appears quite
natural. Formally, it is needed for defining the temperature scale. The zeroth law
of thermodynamics does not have any direct applications but it is needed for
making thermodynamics a complete theory.
9
Given that n is fixed and P is small, experimental results have established the
Boyles law:
P V constant
(1.3)
Furthermore experiments have shown that under the same conditions (Charles &
Gay-Lussac), P V is proportional to temperature:
PV T
(1.4)
(1.6)
10
Following the ideas of the previous section, we define the ideal gas temperature:
P V
PV
T = lim
= lim
(1.7)
P 0 R
P 0 nR
The unit of thermodynamical temperature is Kelvin (1 / 273.16 of the temperature
of the triple point of water). Conversion between different temperature units:
t=
(1.8)
9
(T 273.15) + 32 (Farenheit degrees (t) and Kelvin degrees (T )
5
(1.9)
Also note that the SI unit of pressure P is pascal (Pa; N / m2 ). The atmospheric
pressure is 101,325 Pa. Other units of pressure are atmospheres (atm) and bar.
11
Overview of basic SI units
SI units are an international standard (MKS; 1960).
Quantity
Unit
Symbol Definition
length
mass
time
current
temperature
meter
kilogram
second
ampere
kelvin
m
kg
s
A
K
Quantity
force
energy
electric charge
pressure
magnetic field
frequency
power
voltage
resistance
Unit
newton
joule
coulomb
pascal
tesla
hertz
watt
volt
ohm
Non-SI unit
angstr
om (
A)
inch (in)
foot (ft)
mile (mi)
AMU
eV
cal
torr (Hgmm)
atmospheres (atm)
bar
psi
gauss (G)
Symbol
N
J
C
Pa
T
Hz
W
V
SI unit
meter (m)
meter (m)
meter (m)
meter (m)
kilogram (kg)
joule (J)
joule (J)
pascal (Pa)
pascal (Pa)
pascal (Pa)
pascal (Pa)
tesla (T)
Definition
1 N = 1 kg m s2
1 J = 1 kg m2 s2
1 C = 1 As
1 Pa = 1 N m2
1 T = 1 kg s2 A1
1 Hz = 1 s1
1 W = 1 J s1
1 V = 1 W A1
1 = V A1
12
Conversion factor
1
A= 1010 m
1 in = 2.54 cm = 0.0254 m
1 ft = 12 in = 0.3048 m
1 mi = 5280 ft = 1609.344 m
1 AMU = 1.66054 1027 kg
1 eV = 1.602177 1019 J
1 cal = 4.1868 J
1 torr = 1.33322 102 Pa
1 atm = 1.01325 105 Pa
1 bar = 105 Pa
1 psi = 6.8948 103 Pa
1 G = 104 T
13
Common conversion factors. In Kelvin units, the energy corresponds to kT .
Note also that 1 K = 3.1669 106 Hartree (atomic units; a. u.) or 1 eV = 27.2114
Hartree. To get J mol1 or kcal mol1 , multiply by Avogadros number NA . For
example, 1 eV = 1.6021771019 J and 1 atm = 760 torr.
Mass Energy
kg
eV
1
1.660540 1027
1.112650 1017
1.782663 1036
6.022137 1026
1
6.700531 109
1.073544 109
8.987552 1016
1.492419 1010
1
1.602177 1019
5.609586 1035
9.314943 108
6.241506 1018
1
Spectroscopic units
Hz
cm1
Ry
eV
1
2.99792458 1010
3.289842 1015
2.417988 1014
3.335641 1011
1
1.097373 105
8.065541 103
3.039660 1016
9.112671 106
1
7.349862 102
4.135669 1015
1.239842 104
1.360570 101
1
14
Energy
K
kWh
kcal
1
3.83516 1030 3.298
2.60745 1029 1
8.598
3.0325 1026 1.1630 103 1
7.24292 1022 2.77778 107 2.388
1.16045 104 4.45049 1026 3.827
eV
Pressure
Pa
bar
kp / cm2
torr (mmHg)
atm
1
105
9.807 104
1.333 102
1.013 105
105
1
9.807 101
1.333 103
1.013
1.020 105
1.020
1
1.360 103
1.033
7.5006 103
7.5006 102
7.3556 102
1
760
9.869
9.869
9.678
1.316
1
106
101
101
103
15
(1.10)
16
g
AMU
NA = 39.95
g
mol
Since the partial pressure does not depend on the amount of air, we can choose
the amount of air to be 1.00 g. The number of molecules in the air sample can be
calculated:
n(N2 ) =
(1.00 g) 0.755
= 0.0269 mol
28.02 g mol1
n(O2 ) =
(1.00 g) 0.232
= 0.00725 mol
32.00 g mol1
n(Ar) =
(1.00 g) 0.013
= 0.00033 mol
39.95 g mol1
The total amount of gas (sum of the above components) is 0.0345 mol. The mole
fractions and partial pressures are then:
Mole fraction
Partial pressure (atm)
N2
0.780
0.780
O2
0.210
0.210
Ar
0.0096
0.0096
Note: The numerical values of the AMU to g conversion and NA cancel in the
calculation of ms.
17
Real gases behave like ideal gases only in the limit of zero pressure and high temperature.
Compressibility factor Z indicates deviation from the ideal gas law:
PV
P V
=
Z=
RT
nRT
(1.12)
18
In the limit of high temperature, thermal energy dominates over the potential. At
low temperatures the effect of the attractive part of the potential can be seen more
clearly because thermal energy is not sufficient to smooth out the binding.
Note: The compressibility vs. pressure curves depend on the gas as well as the
temperature.
19
A number of different equations of state for real gases have been proposed:
RT
Ideal gas:
P =
(1.13)
V
RT
a
van der Waals (vdW):
P =
2
(1.14)
V b
V
a
RT
(1.15)
2
P =
Berthelot:
TV
V b
RT
B(T )
C(T )
Virial (Onnes):
P =
1 + + 2 + ...
(1.16)
V
V
V
Alternative forms of Eq. (1.16):
B(T )
C(T )
P V
= 1 + + 2 + ... = 1 + B (T )P + C (T )P 2 + ...
(1.17)
Z=
RT
V
V
where the relation ship between the two constants
are given by:
C(T ) B(T )2
B(T )
and C (T ) =
(1.18)
RT
(RT )2
Note: Temperature where B(T ) = 0 is called the Boyle
temperature. At this temperature the gas behaves ideally over an extended range in pressure.
B (T ) =
The above equations of state can be derived using statistical mehanics and assuming a certain type of pair
interaction.
20
Example. Estimate the molar volume of CO2 at 500 K and 100 atm by treating
it as a van der Waals gas. For CO2 the coefficients are: a = 3.640 atm L2 mol2
and b = 4.267 102 L mol1 .
Solution. First rearrange the van der Waals equation (Eq. (1.14)):
RT 2 a
ab
V 3 b +
V +
=0
V
P
P
P
Roots of a cubic equation (molar volume is the unknown variable) can be found
either analytically by using the appropriate formulas (by using the Maxima program
described in the Appendix). Next, we setup numerical values for the coefficients:
b + RT /P = 0.453 L mol1
a/P = 3.64 102 (L mol1 )2
ab/P = 1.55 103 (L mol1 )3
Thus the equation is:
V 3 0.453V 2 + 3.64 102 V 1.55 103 = 0
The only real valued root is: V = 0.370 L mol1 (0.410 L mol1 for ideal gas).
21
When the equation of state is given, it defines a surface in three dimensional space.
The surface is such that it satisfies the equation state. This is difficult to visualize
in 3-D and therefore 2-D projections are preferred (i.e., one variable is kept constant when plotting). An example is shown below where the temperature was held
constant.
Pressure (Pa)
4.010
2.010
0.2
0.4
0.8
0.6
Molar volume (m )
This example corresponds to an ideal gas at 298.15 K temperature. Such plots for
other equations of state are shown in the following sections.
22
Formally Pc , Tc and Vc define a region on the P -V -T surface where liquid and gas
phase can coexist as two separate phases. Outside this region the phases cannot be
separated.
23
Note: The ends of tie lines indicate pure liquid (VL ) and pure gas phase (VG ) limits.
When the tie line vanishes, VG and VL become identical and the phases cannot be
distinguished from each other. Remember to stay on the isotherms when reading
the above figures - states outside the isotherms are forbidden by the equation of
state. In the last figure, the minima below the critical point have been replaced
with a horizontal tie line.
Gas/liquid becomes supercritical above its critical point. In practice, a supercritical
fluid has properties both of dense gas and low viscosity liquid. It can diffuse through
materials like gas but it can dissolve materials like a liquid. Supercritical fluids are
often used as substitutes for organic solvents (supercritical fluid extraction).
At the critical temperature the following conditions hold (inflection point):
P
=0
(1.19)
V T =Tc
2
P
=0
(1.20)
V 2 T =Tc
Isothermal compressibility is defined as (infinity at critical point):
V
1
=
P
V
(1.21)
Terminology:
24
In addition to critical temperature, critical pressure (Pc ) and critical volume (Vc )
can also be defined by exchanging the roles of variables in Eqs. (1.19) and (1.20).
Expressions for these quantities can be derived for various equations of state. For
the van der Waals equation of state, we have:
a
27b2
Vc = 3b
8a
Tc =
27bR
Pc =
(1.22)
Exercise. Verify that the above expressions are correct. Use the van der Waals
equation of state and Eqs. (1.19) and (1.20). Show that the following results hold
for the Berthelot equation of state:
Pc =
1
6
Vc = 3b
Tc =
2
3
aR
6b3
2a
3bR
1/2
1/2
(1.23)
25
ef f
V
This is similar to the ideal gas law but it uses effective pressure and volume. Reduction in the molar volume is needed because molecules have finite size (i.e. they
are not point-like as assumed in the ideal gas law). This part is related to the repulsive wall of the molecule - molecule interaction. The effective pressure includes
a correction that arises from attractive interactions between molecules (i.e. higher
compressibility). Constants a and b depend on the gas. If monoatomic gas temperature is sufficiently high compared to its atom-atom binding energy, it can be
shown that the parameters a and b are directly related to the atom atom pair
interaction U12 (r) by (see Landau and Lifshitz, Statistical Physics Pt. 1):
Z
a=
|U12 (r)| r 2 dr
(1.25)
2r0
16 3
b=
r0
3
where parameter 2r0 denotes the point where U12 (r) becomes repulsive (i.e. it
becomes positive when the interaction at infinity is taken to be zero).
26
The compressibility factor Z for a van der Waals gas is given by:
V
1
a
a
P V
=
=
Z=
RT
RT V
1 b/V
RT V
V b
(1.26)
Taylor series. Function f that has derivatives of all orders can be expanded in
Taylor series: f (x) = a0 + a1 (x h) + a2 (x h)2 + a3 (x h)3 + ... where the
coefficients are given by:
n
1
d f (x)
a0 = f (h) and an =
(1.27)
n!
dxn
x=h
and we say that the function was expanded about point h. When h = 0, the series
expansion in called Maclaurin series.
Example. Find the Taylor series for ln(x), expanded about x = 1 (i.e. h = 1
above).
Solution. The first derivative of ln(x) is 1/x, which equals 1 at x = 1. The second
derivative is 1/x2 , which equals 1 at x = 1. The derivatives follow a regular
pattern:
n
d f
= (1)n1 (n 1)!
dxn
27
so that we finally have: ln(x) = (x 1) 12 (x 1)2 + 31 (x 1)3 14 (x 1)4 + ...
When b/V is small, we can use the Maclaurin series to expand:
2 3
1
b
b
b
=
1
+
+
+ ...
+
1 b/V
V
V
V
Thus we can write the compressibility factor Z as (cf. Eq. (1.17)):
2
a
b
1
Z =1+
b
+ ...
+
RT
V
V
{z
}
|
(1.28)
(1.29)
=B in Eq. (1.17)
Note that when T is small, 1/T is large and therefore a is important at low temperatures and b at high temperatures. The Boyle temperature can now be obtained
from B(T ) = 0 as:
a
(1.30)
TB =
bR
28
The following realtions can be used to relate a, b and Pc , Tc , Vc to each other:
a=
9
8a
RTc
a
8a
27R2 Tc2
= RTc Vc Tc =
=
and Pc =
=
64Pc
8
9RVc
27Rb
8b
27b2
b=
Vc
RTc
=
Vc = 3b
8Pc
3
(1.31)
(1.32)
27 (RTc )2
RTc
and b =
64Pc
8Pc
(1.33)
29
a=
27 (RTc )2
=
64Pc
= 0.5583
N
m2
27 8.3145
2
Nm
z}|{
J mol1 K1 305.34 K
64 4.870 106
m6
mol2
= 0.5583
= 0.5583
Nm4
mol2
Pa
|{z}
Nm2
9.869 106 atm (10 dm)6
Pa m6
=
0.5583
mol2
mol2
30
The van der Waals equation fails in the neighborhood of the critical point:
Vc V (Tc T )1/2
(1.34)
However, experiments show that the exponent is close to 0.32 rather than 1/2.
Other properties that depend on (Tc T ) show similar discrepancies as well.
31
(1.35)
By using these variables, it can be shown that the van der Waals equation becomes:
8Tr
3
Pr =
2
(1.36)
3Vr 1
Vr
Note that the gas dependent constants
a and b disappeared from the equation. Thus within the van der Waals
equation of state, all gases (regadless
of their a and b values) should yield
identical curves. The compressibility
factor for van der Waals equation can
also be expressed in reduced units:
9
Vr
(1.37)
Z=
8Vr Tr
Vr 1
3
32
Intensive variables:
System
One-phase
Two-phase equilibrium
Three-phase equilibrium
Degrees of freedom
F =2
F =1
F =0
Note: If multiple species (i.e., different gases) are included in one system then additional degrees of freedom must be specified (increases by Ns 1, where Ns is the
number of species; the Gibbs phase rule). Furthermore, a non-reactive system
was assumed.
Extensive variables: one extensive variable per phase (i.e., the amount of each
phase).
Example. Temperature of liquid 4 He (T < 4 K) can be determined from the
helium vapor pressure in a closed container. Note that both liquid and gas phases
coexist and thus only one variable is needed to specify the state of the system (both
intensive variables). The experimentally observed phase diagram and the relation
between helium vapor pressure and temperature are shown below.
33
1000
800
600
400
200
2
3
Temperature (K)
Example. For a non-reactive system with two species, two extensive variables are
required for a complete description (i.e., the amount of each species).
Recall terminology: intensive state of system = described by intensive variables
(i.e., they do not depend on the size of the system); extensive state of system =
described by extensive variables (i.e., they depend on the size of the system).
Note: The choice of variables is not unique, only the number of variables is fixed.
34
(1.38)
(1.39)
(1.40)
35
Differentiate both sides with respect to and apply the chain rule:
N
X f (x1 , x2 , ..., xN )
(xi )
(xi )
i=1
|
{z
}
=xi
N
X
f (x1 , x2 , ..., xN )
xi = kf (x1 , x2 , ..., xN )
xi
i=1
(1.43)
36
Partial molar volume tells us in practice how much the volume changes when a small
amount of gas component i is added to the gas mixture (P , T and the amounts of
other components remain constant). In terms of differentials:
dV = Vi dni
(1.44)
where dni is an infinitesimal change in the amount of component i and dV is the
change in total volume occupied by the gas mixture. The total differential combining
all the components is:
dV = V1 dn1 + V2 dn2 + ... + Vn dnN
(1.45)
The non-differential form is convenient to write in terms of mole fractions (yi ):
V = V1 y1 + V2 y2 + ... + VN yN
(1.46)
37
In the following we calculate the effect of additional pressure arising from gravitation. This would be important, for example, if we consider a column of gas where
a difference in gravitational force exists:
F = mg
(1.47)
where m is the total mass of the gas and g is the standard gravitation constant
(9.80665 m s2 ). To be exact, this constant depends on the position of measurement
on earth. It can be as high as 10.7 m s2 . In non-SI units the value of this constant
is 32.2 ft s2 .
The amount of gas within dh is:
dmgas = A dh
(1.48)
(1.49)
38
For an ideal gas, we have an expression for the gas density :
nM
PM
P V = nRT and =
=
V
RT
where M is the molar mass (kg mol1 ). Thus we have:
PM
gdh
dP =
RT
(1.51)
(1.52)
P0
ln
dP
=
P
P
P0
Zh
0
(1.53)
gM h
RT
(1.54)
P = P 0 e
This is known as the barometric formula.
gM
dh
RT
gM h
RT
(1.55)
39
A system can exchange energy with its environment in the form of work or heat.
41
(2.57)
~ is the vector defining the path (direction and length) along which the work
where L
is being done. Subscripts refer to the Cartesian components of the corresponding
vectors.
Pressure P (Pa) is a scalar quantity:
P =
F
A
(2.58)
where F denotes the force component (N) that is perpendicular to the surface with
area A (m2 ).
Example. What is the work done on/by a moving piston inside a cylinder?
Solution. We assume quasistatic system: P = Pext at all times. This corresponds
to slow movement of the piston.
42
The system consists of the piston and the cylinder and surroundings to outside of the cylinder. The work done on/by
the the system (piston) is:
w = F L = P AL = P V
| {z }
V
where the piston moved by L and V is the change in the cylinder volume.
Note that L and V change sign depending on the process (i.e. compress or
expand). Note the sign convention: + means that the surroundings did work on
the system and means that the system did work on the surroundings. The
previous expressions can also be written using differentials. For example in this
quasistatic case we have:
dw = Pext dV = P dV
(2.59)
Example. What is the work done by a system slowly lifting an object weighting 1
kg by 0.1 m? The opposing force is the gravitational force.
Solution. The work is given by (note slow corresponds to a quasistatic system):
w = F h = mgh = (1 kg) (9.807 m s2 ) (0.1 m) = 0.9807 J
43
Integration of Eq. (2.59) yields an expression for work when P = P (V, T ):
w=
point
Z 2
point
Z 2
P (V, T )dV
(2.60)
point 1
point 1
where the end points of the line integrals are located on the (P, V ) plane.
Example. The state of a mole of gas can be changed from (2P0 , V0 ) to (P0 , 2V0 ) by
infinitely many different quasistatic paths (two different choices are shown below):
Path 1:
1. Pressure is constant (2P0 ), expansion from V0
to 2V0 (heat flows in).
2. Volume is constant (2V0 ), drop in pressure
from 2P0 to P0 (heat flows out). Total amount
of work (w) = 2P0 V0 (contribution from the
first segment only).
Path 2:
Paths on the P V plane.
44
Adiabatic process: The system is thermally insulated and it cannot exchange heat
with its surroundings. The two paths in the previous example involved heat exchange and therefore they were not adiabatic processes.
For an adiabatic process in a closed system, the change in systems internal energy
(U ) is directly related to work:
U = w
(2.61)
Note that U follows the same sign convention that w does. Molar internal energy
(J mol1 ).
is expressed with a symbol U
Example. A thermally insulated water tank with a stirrer.
45
Lets consider a water tank that is in thermal contact and has no stirrer (i.e. no
source of work):
If the hot plate is warmer than the water, the water temperature increases (heat
transfer). In previous example increase in system temperature was achieved by
doing work on the system. Thus increase in the internal energy of the system can
be achieved either by doing work on it or by transferring heat into it.
Both work and heat are forms of energy crossing the boundary between the system
and the surroundings. When no work is done on the system, the change in its
internal energy can be expressed as:
U = q (no work done)
Both work and heat have the same SI unit (J).
(2.62)
46
(2.63)
dU = dq + dw
(2.64)
Note that q and w indicate changes in heat and work in this sense is missing
in this notation.
The above equations are obtained by combining the results from previous section for
work and heat. The internal energy of a system (U ) is a function of state variables
(for example, P, V, T ; closed system). Internal U energy is an extensive property.
The differential corresponding to the internal energy, dU , yields zero when integrated over a closed loop in the state space:
I
dU = 0
(2.65)
However, the component differentials dq ad dw do not have this property.
Energy may be transferred in one form or another, but it cannot be created or
destroyed (energy conservation).
U depends on state variables, e.g. U (T, V, n) or U (T, P, n) for pure substances. For
mixtures, the composition must also be specified.
47
A differential with two independent degrees of freedom, can be expressed as (Pfaffian form):
du = M (x, y)dx + N (x, y)dy
(2.66)
=N
(2.67)
48
M
y
N
x
2xy +
x2
9x2
y
3x3
y2
= 2x
= 2x
9x2
y2
9x2
y2
(2.68)
N
x y,z
x,z
N
P
=
z x,y
y x,z
P
M
=
z x,y
x y,z
M
y
(2.69)
49
Example. Show that the following differential is not exact but it can be transformed into an exact differential by dividing both sides by x.
du = 2ax2 + bxy dx + bx2 + 2cxy dy
|
|
{z
}
{z
}
=M
=N
Solution. Calculate the required partial derivatives and use Eq. (2.67):
M
2ax2 + bxy
= bx
=
y x
y
x
N
=
= 2bx + 2cy 6= bx
bx2 + 2cxy
x y
x
y
(2ax + by) = b
=
y x
y
x
(bx + 2cy) = b
=
x y
x
y
This procedure of transforming an inexact differential into an exact one is called
the method of integrating factors.
50
Line integrals are evaluated over a specified path in the x-y (or in higher dimensional spaces):
(xZ
1 ,y1 )
du =
(x0 ,y0 )
(xZ
1 ,y1 )
(2.70)
(x0 ,y0 )
However, this is not yet well defined because there are infinitely many paths that
connect points (x0 , y0 ) and (x1 , y1 ) in the x-y plane. One such path is shown below:
51
An integration path (C) must be specified for line integrals:
Z
Z
du =
[M (x, y)dx + N (x, y)dy]
C
(2.71)
Note that line integrals are sometimes also called path integrals.
When the integration path is defined by functions y = y(x) and x = x(y), the line
integral along this path can be calculated as:
Zx1
Zy1
Z
du =
M (x, y(x))dx + N (x(y), y)dy
(2.72)
x0
y0
where path C is the straight-line segment given by y = 2x + 3 from (0, 3) to (2, 7).
Solution. In the first term, y must be replaced by 2x + 3, and in the second term
x must be replaced by (1/2)(y 3),
52
du =
Z2
0
Z7
3
3
(y 3) + 4y dy
2
2
7
8x2
9
11y 2 /2
+ 9x +
y = 126
2
2
2
0
3
If differential du is exact, then the line integral does not depend on path but only
on the end points:
Z
Z
u
u
dx +
dy = u(x1 , y1 ) u(x0 , y0 )
(2.73)
du =
x
y
C
For closed paths, this always yields zero (this does not hold for ineact differentials):
I
du = 0 (over a closed loop)
(2.74)
Example. Show that the line integral of the previous example has the same value
as the line integral of the same differential on the rectangular path from (0, 3) to
(2, 3) and then to (2, 7).
53
Solution. The integration path is not a single curve but two line segments. So we
must carry out the integration separately for each segment. On the first segment, y
is constant, so dy = 0 and the integral containing dy vanishes. On the second line
segment, x is constant, so dx = 0 and the integral containing dx vanishes:
du =
Z2
0
(2x + 9)dx +
Z7
3
(6 + 4y)dy =
2
7
4y 2
2x2
= 126
+ 9x + 6y +
2
2 3
0
Note that this is the same result as obtained in the previous example.
Example. Show that differential du = dx + xdy is inexact and carry out line
integration using two different paths between points (0, 0) and (2, 2). Path 1 is
defined as: straight line from (0, 0) to (2, 2) and path 2 as: rectangular path (0, 0)
to (2, 0) to (2, 2).
Solution. First we show that the differential is inexact:
(1) = 0 and
(x) = 1 6= 0
y
x
x
y
du is inexact differential.
54
Integration along path 1 (denoted by C1 , a straight line y = x):
Z
du =
C1
dx +
C1
xdy =
Z2
dx +
ydy = x|20 +
C1
Z2
du =
C2
C2
dx +
xdy =
Z2
0
C2
dx +
Z2
2
y 2
=4
2 0
Thus the value of the line integral depends on path. This is because du is inexact.
For three independent variables, line integral is defined as:
Z
Z
du =
[M (x, y, z)dx + N (x, y, z)dy + P (x, y, z)dz]
C
(2.75)
du =
Zx1
x0
Zy1
y0
N (x(y), y, z(y))dy +
Zz1
z0
55
If du is exact differential, the value of the line integral depends only on the endpoints:
Z
du = u(x1 , y1 , z1 ) u(x0 , y0 , z0 )
(2.77)
C
(2.78)
56
Work can be done on/by a gas upon compression/expansion. In the following example gas is the system and the surroundings constitutes of a piston and a thermostat
(i.e. a container that keeps the temperature constant).
Assumptions:
No friction
No external pressure outside the
cylinder (from atmosphere)
Cylinder immersed in a
thermostat (constant T )
Compression by a piston with mass m.
57
where A is the area of the piston.
At both points the external pressure is Pext = mg
A
Thus the force pushing the piston down is mg. At the end of the process the gas
pressure will be the same as the external pressure: Pext = P2 . Recall that work was
defined as force distance, which in this case means (see Eq. (2.60), P1 < P2 ):
wcomp =
f in
f in
Z
Z
Pext dV = P2 (V2 V1 ) = P2 (V1 V2 ) > 0
dw =
| {z }
ini
(2.79)
ini constant
f in
Z
mg
Pext dV =
(Ah) = mg h (force distance)
| {z }
A
| {z }
ini constant
(2.80)
=V
where the positive sign for wcomp signifies that work was done on the system and h
denotes the distance that the piston moved (the shaded area on the previous page
P -V plot).
Next, consider expansion of a gas in a two stage process:
58
f in
Z
(2.81)
ini constant
Note: |wcomp | > |wexp |. More work is required to compress the gas than can be
obtained by expansion. This process is irreversible.
Is it possible to move the piston in such a way that |wexp | = |wcomp |? In other
words, is it possible to make the process reversible in terms of work?
59
Yes. Instead of single-step compression, we should use many compression steps with
increasing external pressure Pext . This can be achieved by increasing the mass
gradually (m1 < m2 < m3 < ...):
In other words: do not apply all the force at once but increase it gradually. Note
that in a reversible process the pressure inside the cylinder and the external pressure
are equal at all times. In this case the work is obtained by (P is the pressure inside
the cylinder):
wcomp =
f in
Z
ZV2
ZV2
dw = Pext dV = P dV
ini
V1
V1
(2.82)
60
Expansion work wexp can also be carried out using infinite expansion steps (m1 >
m2 > m3 ...):
This infinite expansion process is also reversible. The expression for work is now:
wexp =
f in
Z
ZV1
dw = P dV
ini
(2.83)
V2
Work over one closed cycle (reversible compression followed by reversible expansion):
wcycle =
ZV2
P dV
{z
P dV =
V2
V1
ZV1
}|
{z
ZV2
V1
P dV +
ZV2
P dV = 0
(2.84)
V1
compression expansion
Thus the infinitesimal process is reversible. Many calculations can be carried out
exactly only for reversible processes. Most processes in nature are, however, irreversible. Sometimes they can be approximated as reversible processes.
61
For a reversible process Pext = P (in the cylinder) always. Assuming that the gas
in the cylinder is ideal, we have for reversible expansion (T constant, V1 < V2 , ini
= 1 and fin = 2):
wexp,rev =
ZV2
Pext dV =
V1
ZV2
P dV =
V1
ZV2
V1
nRT
dV = nRT ln
V
V2
V1
(2.85)
wcomp,rev = nRT ln
For an ideal gas at constant temperature, we have P1 V1 = P2 V2 and Eq. (2.85) can
then be written:
P1
P2
V2
= nRT ln
= nRT ln
(2.86)
wexp,rev = nRT ln
V1
P2
P1
The maximum amount of work of isothermal expansion of a van der Waals gas is:
wexp,rev =
ZV2
V1
nRT
an2
2
V nb
V
dV = nRT ln
V2 nb
V1 nb
+ an2
1
1
V1
V2
(2.87)
Note: During reversible processes the system and the surroundings are in equilibrium. However, such processes are ideal since they take infinitely long time to
proceed.
62
Example. Calculate the work done when 50 g of iron reacts with hydrochloric acid
in (a) a closed vessel of fixed volume, (b) an open beaker at 25 C. The reaction is:
Fe(s) + 2HCl(aq) FeCl2 (aq) + H2 (g)
Assume that H2 follows the ideal gas law.
Solution.
(a) The volume cannot change, so no P V -work is done and wexp = 0.
(b) The gas drives back the atmosphere and therefore wexp = Pext V . We can
neglect the initial volume (V1 ) because the final volume (V2 ) after production of
gas, is much larger. We assume that H2 behaves according to the ideal gas law (n
moles of H2 ):
V = V2 V1 V2 =
nRT
nRT
=
P
Pext
where P is the gas pressure and Pext the atmospheric pressure. Note that solids
have negligible volumes compared to gases and therefore we have:
wexp = Pext V Pext
nRT
= nRT
Pext
63
Because the molar mass of Fe is 55.85 g mol1 , it follows that:
wexp
50 g
8.3145 JK1 mol1 (298.15 K) = 2.2 kJ
55.85 g mol1
Thus the system (H2 gas from the reaction) does 2.2 kJ of work driving back the
atmosphere.
Example. Work of expansion of an ideal gas. One mole of an ideal gas expands
from 5 to 1 bar at 298 K. Calculate wexp (a) for reversible expansion and (b) for
an irreversible expansion against a constant external pressure of 1 bar.
Solution. (a) We use Eq. (2.86) with P1 = 5 bar (ini) and P2 = 1 bar (fin):
wexp,rev = nRT ln
P2
P1
= (1 mol) 8.3145 J K1 mol1 (298 K)ln
1 bar
5 bar
= 4000 J
1
= nRT
P2
P1
P1
1 bar
(1 mol) 8.3145 J K1 mol1 (298 K)
1 = 2000 J
5 bar
64
Differential
Comments
Conjugate pairs
Expansion,
Hydrostatic
Pext dV
P V -work
Pressure P , Volume V
Surface expansion
dAs
is surface tension
and As is surface area
and As
Extension,
elongation
F dL
Change in length
Force F , Length l
Electrical
dQ
Transport of electrons
Potential difference ,
Electric charge Q
65
Surface tension (): Force / distance (unit N m1 ):
F
=
L
(2.88)
surface tension.
For a film with two sides, the force F acting on the bar (L is the width) is given by:
F = 2L
(2.89)
(2.90)
(2.91)
dw = dAs
(2.92)
(m2 )
by surface tension (J
66
What is the microscopic origin of surface tension? Consider water surface, for example:
Molecules residing on the curved surface are missing neighboring atoms which they
could hydrogen bond with. Surface tension depends clearly on the molecule
molecule (or atom atom) interaction strength. For example, surface tension of
water (hydrogen bonding) is much larger than for liquid argon (van der Waals
binding).
67
Example. Calculate the amount of work that a spherical bubble in superfluid 4 He
does when its radius changes from 10
A to 20
A. Surface tension of L-He is 0.18
cm1
A2 .
2nd order
Zrf
ri
dw = 8
Zrf
ri
rdr = 4 rf2 ri2
68
We could proceed by changing everything to SI units, but here the expression is
quite simple and we just make sure to use compatible units:
h
2
2 i
2
w = 4 0.18 cm1
A
20
A 10
A
= 680 cm1 = 1.4 1020 J
Notes:
The work in above example is positive, which means that work was done on
the bubble interface (i.e., interface area becomes larger).
Recall that we dont usually use notation w to indicate change in work
because work is a relative quantity. The same applies for q. Instead of w,
we just use w.
Extension/elongation work:
dw = F dL
(2.93)
where F is the extension force and dL is the displacement. For example: elongation
of a rubber band.
69
Electrical work:
dw = dQ
(2.94)
where is the electric potential difference and dQ is the differential change in charge.
For example: electron transport in electrolytic cell.
A more complete form of the first law is (heat dq + work dw):
dU = dq Pext dV + dAs + F dL + dQ
(2.95)
70
71
dq V =
U
T
dT
(2.99)
Both the temperature and the heat transfer can be measured and thus it is convenient to define heat capacity CV (T ) at constant volume as:
U
dq V
=
(2.100)
CV (T )
dT
T V
V .
For one mole of substance, heat capacity is denoted by C
Note: Temperature and heat are two different quantities. On molecular scales temperature is related to the kinetic energy distribution of molecules in the substance.
Heat is related to the total energy of molecules (including potential energy).
At constant volume, Eq. (2.100) may be multiplied by dT and integrated (see also
Eq. (2.99)):
UV =
ZT2
CV (T )dT = qV
(2.101)
T1
(2.102)
72
Now we know what
U
T
U
V
U
dV
dq = Pext +
V T
(2.103)
(2.104)
73
(2.106)
(2.107)
qP = (U2 + P V2 ) (U1 + P V1 ) = H2 H1 = H
{z
} |
{z
}
|
(2.108)
=H2
=H1
74
From Eq.(2.108) it also follows that the differentials corresponding to heat (dqP )
and enthalpy (dH) must be equal (at constant pressure):
dqP = dH (H is a state function and dH is exact)
Because dH is an exact differential, we can write the total differential as:
H
H
dT +
dP
dH =
T P
P T
(2.109)
(2.110)
(2.112)
Integration of Eq. (2.111) gives an expression for change in enthalpy (cf. Eq.
(2.102)):
HP =
ZT2
T1
(2.113)
75
Example. Enthalpy change can be measured with a constant-pressure calorimeter.
A simple example of such calorimeter is the coffee-cup-calorimeter. 10.0 g of ice
at 273 K is added to such calorimeter containing 100.0 g of water at 303 K. The
heat capacity of the calorimeter is 10.0 J K1 . What is the final temperature of
the water in the cup? It was observed that all ice melted and therefore we know
that all ice was transformed to water.
Solution. Overall, the enthalpy is conserved (thermal insulation):
H
+ H cold
ice
water
final temp.
Hice
+ Hwater + Happaratus = 0
{z
added ice
cools
cools
warms to
melts
{z
water present
in the cup
{z
coffee cup
Enthalpy change for the melting process can be obtained by multiplying the mass of ice by the enthalpy of fusion (will be discussed
in more detail later):
1
= 3330 J
melts = mice f us H = (10.0 g) 333 J g
|
{z
}
from table
Next we calculate the change in enthalpy when water at 273 K warms up to the
final temperature (see Eq. (2.113)):
76
H cold
water
warms to
= mice CP,H2 O T = (10.0 g) 4.19 J g1 K1 Tf 273 K
final temp.
where Tf denotes the (still unknown) final temperature. On the other hand, the
water initially at 303 K cools down to the final temperature Tf :
Hwater = mwater CP,H2 O T = (100.0 g) 4.19 J g1 K1 Tf 303 K
cools
Finally, we have to consider cooling of the calorimeter from 303 K down to the
final temperature:
Happaratus = Ccalorimeter T = 10.0 J K1 Tf 303 K
cools
77
(3330 J) + 41.9 J K1 Tf 273 K + 419 J K1 Tf 303 K
+ 10.0 J K1 Tf 303 K = 0 Tf = 293 K
78
Values for molar heat capacities can be found, for example, from the National Institute of Standards and Technology (NIST) database (http://webbook.nist.gov/chemistry/).
It has been observed experimentally that molar heat capacities at constant pressure
depend on temperature.
Empirical power series expansion:
P + T + T 2
C
(2.114)
T 2 T12 +
T 3 T13
2 2
3 2
(2.115)
Values for , and are given in the above database. Note that the temperature
range for these polynomial fits must be observed. In above, 1 refers to the initial
and 2 to the final state
79
The relation between CP and CV :
Consider Eq. (2.98) at constant volume and pressure (with pressure P = Pext ):
dq P =
U
U
dT + Pext +
dV = CV dT + Pext +
dV (2.116)
V T
V T
V
dq P
V
U
division by dT
= CV + Pext +
dT
V T
T P
| {z }
U
T
=CP
U
V
>0
CP CV =
Pext + V
T P
T
| {z }
| {z }
(term 1)
(term 2)
where term 1 is greater than zero (consider, for example, a van der Waals gas
where this term is equal to a/V ) and term 2 is also greater than zero because,
at constant pressure, increase in temperature will result in increased volume. Thus
we conclude that heat capacity at constant pressure is always greater than the heat
capacity at constant volume. This is not surprising because the former also includes
the P V -work and the work required to pull molecules apart.
80
Example. Calculate CP CV for an ideal gas. Note that in ideal gas, the gas
molecules do not interact with each other.
Solution. To calculate the difference in heat capacities, we use Eq. (2.116). Because molecules do not interact with each other, the internal energy does not depend
on volume. Hence (U/V )T = 0. For an ideal gas, we have P V = nRT , from
which we can evaluate the partial derivative (V /T )P = nR/P . Inserting the
partial derivatives into Eq. (2.116) gives CP CV = nR. This is the P V -work
required when temperature changes by 1 K.
Notes:
Classical thermodynamics does not deal with molecular level information. It
knows nothing about atoms and molecules. In fact, when classical
thermodynamics was first developed, scientists did not know about atoms
and molecules at all.
For liquids and solids, CP CV because they have small thermal
expansivities.
Results from the kinetic gas theory results for monoatomic ideal gas:
= 3 RT
U
2
= 5 RT
H
2
(2.117)
(2.118)
V = 3 R
C
2
5
P = R
C
2
(2.119)
(2.120)
81
Adiabatic
Joule-
Thompson
Ti
6=
Pi > Pf .
Tf
with
and
On compression,
the work is given by wc = Pi V =
Pi Vi Vf = Pi Vi 0 = Pi Vi (positive because
work is done on the system (gas)).
On expansion,
the work is
now we = Pf V =
Pf Vi Vf = Pf 0 Vf = Pf Vf (negative because work is done by the system (gas)).
James Prescott
Joule, English
physicist (1818 1889)
(2.121)
82
The system is thermally insulated, so that q = 0 (no heat exchange). Using Eq.
(2.63) we obtain:
=U
f U
i = q + w = Pi Vi Pf Vf
U
(2.122)
Rearrangement of this equation gives:
i + Pi Vi
f + Pf Vf = U
U
{z
} | {z }
|
=Hf
(2.123)
=Hi
This states that the enthalpy is conserved in the process (isenthalpic process). Based
on the experimental observation, we define the Joule-Thomson coefficient:
T2 T1
T
JT = lim
=
(2.124)
P 0 P2 P1
P H
which gives the change in temperature when pressure changes. At high temperatures
the coefficient is negative (J-T process results in heating) and at low temperatures
it is positive (J-T process results in cooling). The temperature, where the coefficient
is zero, is called the inversion temperature. The inversion temperature for N2 is 607
K and for H2 204 K.
Notes:
The J-T coefficient is zero for ideal gases.
The cooling effect can be understood by decrease in the van der Waals
interaction due to lower pressure (i.e. increased potential energy) and
decrease in the kinetic energy (i.e. lower temperature).
83
(2.125)
If work is done against external pressure, the temperature of the system drops
(note: adiabatic process; no thermal contact). If the system was in thermal contact,
the system would exchange heat with its surroundings and the temperature would
remain constant.
Consider an ideal gas. The change in internal energy (U ) depends only on temperature and not on volume (Eqs. (2.97), (2.105), (2.100)):
dU = CV dT , if CV is independent of T , we have: U = CV T
(2.126)
84
When adiabatic expansion is carried out reversibly, the equilibrium pressure can be
substituted for the external pressure (Pext = P ). For an ideal gas we have:
V dT = R dV
V dT = P dV = RT dV C
(2.127)
C
T
V
V
If we assume that the heat capacity is independent of temperature and integrate
both sides, we get:
V
C
ZT2
T1
dT
= R
T
ZV2
dV
V ln T2 = R ln V1
C
T1
V
V2
(2.128)
1
V
P C
V = R and thus we can rewrite Eq. (2.128) as:
For an ideal gas we have C
1
P
V1
C
T2
=
where =
(2.129)
T1
V2
CV
The same result can be written in alternative forms:
(1)/
P2
T2
=
T1
P1
P1 V1 = P2 V2
(2.130)
(2.131)
85
Initial
Isothernal expansion
Adiabatic expansion
P (bar)
0.8
0.6
P
0.4
Final
0.2
0
20
25
30
40
35
45
-1
Molar volume (L mol )
50
T2 = T1
1
V1
V2
Substituting the values, we get the final T2 (see Eqs. (2.119) and (2.120)):
0.50 L 0.666
T2 = (298 K)
= 188 K
1.00 L
From the temperature difference we
can calculate the work:
V T
wrev = nC
= (0.020 mol) 12.48 J K1 mol1
(110 K) = 27 J
Note the sign: the gas does work.
86
2.11 Thermochemistry
Chemical reaction (or phase change) in a system is said to be:
Exothermic, if heat q < 0 (i.e. it releases heat).
Endothermic, if heat q > 0 (i.e. it requires heat).
Since enthalpy is an extensive property that is a function of the state of the system
(homogenous of degree 1; see Eq. (1.42)), we can express it in terms of partial molar
enthalpies:
dH =
Ns
X
i dni
H
(2.132)
i=1
where Ns is the number of species and Hi is the molar enthalpy of species i. When
the temperature and the pressure are constant, Eq. (2.109) gives:
dH = dqP =
Ns
X
i=1
i dni
H
(2.133)
87
Notation of chemical reactions:
0=
Ns
X
vi B i
(2.134)
i=1
where Bi is the chemical formula of species i and vi is the corresponding stoichiometric coefficient.
Example. Use the above notation to express the following chemical reaction:
H2 +
1
O2 H 2 O
2
1
O2 + 1H2 O
2
(2.135)
where ni0 is the initial amount of species i (mol), is the extent of reaction (mol)
and vi are the stoichiometric coefficients (dimensionless). Note that evolves from
initial value of zero to the final value where any of the component ni s reach zero.
88
By combing Eqs. (2.133) and (2.135) we get:
dH = dqP =
Ns
X
i=1
i vi d
H
|{z}
(2.136)
=dni
(2.137)
where r H is the reaction enthalpy. Reaction enthalpy tells us the rate of enthalpy
change at given point of reaction (). Note that the reaction enthalpy depends on
the way the chemical equation is written. For example, 2H2 +O2 2H2 O has twice
the reaction enthalpy than H2 + 12 O2 H2 O. The SI unit of reaction enthalpy is
J mol1 .
Thermodynamic standard state:
Denoted by superscript (degree). For a standard state, Eq. (2.137) reads:
r H =
Ns
X
i=1
vi H
i
(2.138)
89
The standard state of a substance at a specified temperature is its pure form at 1
bar.
The phase of the substance is indicated by g (gas), l (liquid) or s (solid).
Example. What is the standard state of CH3 CH2 OH (l) at 298 K?
Solution. The standard state is pure liquid ethanol at 298 K and 1 bar (external)
pressure.
Enthalpy change in chemical reactions:
Lavoisier & Laplace in 1780:
The
heat absorbed in decomposing a
compound must be equal to the heat
evolved in its formation under the
same conditions.
Antoine Lavoisier,
Pierre-Simon
French scientist
Laplace, French
(1743 - 1794)
mathematician
and astronomer
(1749 - 1827)
90
Note that both previous statements follow directly from the first law of thermodynamics.
Example. Burning of graphite (in presence of excess O2 ) to CO2 :
C(graphite) + O2 (g) CO2 (g) r H = 393.509 kJ mol1
Note that the sign signifies that the reaction is exothermic (releases heat to the
surroundings).
Example. Burning of CO (in presence of excess O2 ) to CO2 :
CO(g) +
1
O2 (g) CO2 (g), r H = 282.98 kJ mol1
2
1
O2 (g), r H = +286 kJ mol1
2
Note that the + sign means that the reaction is endothermic (draws heat from
the surroundings).
For reactions, that cannot be studied directly, reaction enthalpies can be obtained
by dividing the reaction into parts, which can be studied with the desired accuracy.
91
Example. The heat released by burning graphite to CO is difficult to measure
accurately because part of it might produce CO2 and part of the graphite might
not react at all. How can we determine the reaction enthalpy in this case?
Solution. First we write reactions for burning of graphite yielding CO2 and burning
of CO to CO2 (excess O2 to ensure that reactions proceed all the way):
(1) C(graphite) + O2 (g) CO2 (g), r H = 393.51 kJ mol1
(2) CO(g) +
1
O2 (g) CO2 (g), r H = 282.98 kJ mol1
2
Consider the reverse reaction of the last equation (Lavoisier & Laplace):
(3) CO2 (g) CO(g) +
1
O2 (g), r H = +282.98 kJ mol1
2
1
O2 (g)
2
and after cancelling terms from the left and right sides, we get:
1
C(graphite) + O2 (g) CO
2
The total reaction enthalpy is obtained by adding the corresponding values for the
partial reactions (Hess): r H = (393.51 + 282.98) kJ mol1 = 110.53 kJ mol1 .
92
Recall that we dont usually know the absolute values for the internal energy U but
only changes in it. Since H = U + P V , the same applies for H and therefore we
usually just concentrate on changes in enthalpy. However, we can specify a reference
state and consider differences in enthalpy from this state (relative enthalpies or
enthalpies of formation; f H ). Such a state is obtained from the stoichiometric
amounts of the molecules in the given substance in the standard state and at the
temperature under consideration. When enthalpy of formation is calculated, one
mole of the product must be considered.
Example. What are the enthalpies of formation (f H ) for CO2 (g), CO(g), C(g)
and O(g) under standard conditions (25 C and 1 bar) and why?
Solution. Consider first CO2 (g). The formation reaction for CO2 is: C(s) +
O2 (g) CO2 (g) with f H = 393.51 kJ mol1 . Note that always for the
formation reaction: f H = r H . The situation for CO is analogous. How
about C(g)? In this case the formation reaction is: C(graphite) C(g). The
reaction enthalpy for this reaction is 716.7 kJ mol1 under the present conditions
(Numerical data can be found from the NIST Chemistry Webbook).
When the enthalpies of formation for each substance in a chemical reaction are
known, it is possible to calculate the reaction enthalpy by (both enthalpies under
the same conditions):
r H =
Ns
X
i=1
vi f Hi
(2.139)
93
Enthalpy of formation for a compound can be determined:
1. Experimentally: Calorimetric measurements, temperature variations of
equilibrium constants, spectroscopic determination of dissociation energies.
2. Theoretically: From the first principles using electron structure methods.
Note: The first law of thermodynamics cannot determine if the reaction will occur
spontaneously or not.
Example. What are the standard enthalpy changes at 298.15 K for the following
reaction:
CO2 (g) + C(graphite) 2CO(g)
Solution. First we note that for graphite f H (graphite) = 0 (choice of reference
point). The other f H values can be looked up from the NIST chemistry webbook
database: f H (CO) = 110.5 kJ mol1 and f H (CO2 ) = 393.5 kJ mol1 .
Note that these values are given at 298.15 K temperature. The standard reaction
enthalpy is now given by:
r H = 2f H (CO) f H (CO2 ) f H (C(graphite))
= 2 110.5 kJ mol1 393.5 kJ mol1
= 172.5 kJ mol1
94
If we need to calculate reaction enthalpies at some other temperature than 298.15
K, we have to use heat capacity CP to account for the change in temperature
(Kirchhoffs law):
Starting from the reactants, calculate the change in enthalpy when temperature is
changed from T to 298 K, use the known reaction enthalpy for the reaction, and
calculate the change in enthalpy when temperature is changed back from 298 K to
T.
95
The expression for r H at the given temperature T is therefore:
=
RT
...
298 K
r HT =
298
Z K
}|
CP,react
(T )dT +r H298
K+
r HT = r H298K
+
ZT
(2.140)
CP,product
CP,react
z }| {
r C P
(T )
298 K
Ns
X
CP,product
(T )dT
298 K
r C P
(T ) =
ZT
dT
(2.141)
(2.142)
i=1
), we use
To get the temperature dependency of the reaction heat capacity (r CP
the empirical form of Eq. (2.114):
r C P
= r + (r ) T + (r ) T 2
with r =
Ns
X
i=1
v i i , r =
Ns
X
i=1
vi i and r =
(2.143)
Ns
X
i=1
vi i
96
If the values , and for each species i are known, it is possible to insert them
into Eq. (2.143) and further use Eq. (2.142) to obtain the reaction heat capacity:
r HT = r H298
ZT
298 K
= r H298
r + r T + r T 2 dT
+ r (T 298.15 K) +
r 3
T (298.15 K)3
3
(2.144)
r 2
T (298.15 K)2
2
In principle, it would be possible to choose another reference temperature. Sometimes absolute zero temperature is used (r H0 K ). For a diatomic molecule, this
would correspond to the bond dissociation energy. Note that other empirical parametrizations are often also used. For example, the NIST chemistry webbook uses
the Shomate equation, which is essentially a 3rd order polynomial with a 1/T 2 term
added.
97
Example. The standard enthalpy of formation of gaseous H2 O at 298 K is 241.8
kJ mol1 (constants obtained from the NIST chemistry webbook). Estimate its
value at 373 K given the following values for the molar heat capacities at constant
pressure: H2 O (g): 34 J K1 mol1 ; H2 (g): 28.8 J K1 mol1 ; O2 (g): 28.9 J
K1 mol1 . Assume that the heat capacities are independent of temperature.
Solution. The reaction is: H2 (g) + 21 O2 (g) H2 O(g). When the heat capacities
are independent of temperature, Eq. (2.141) can be replaced with:
r H373
= r H298
:
To proceed, we need to evaluate r CP
1
r C P
= CP,H
C
+
C
P,H
(g)
2
2 O(g)
2 P,O2 (g)
1
= 34 JK1 mol1 28.8 JK1 mol1 + 28.9 JK1 mol1
2
r H373
K:
= 241.8 kJ mol1 + (75 K) 9.4 J K1 mol1 243 kJ mol1
98
2.13 Calorimetry
99
Constant pressure calorimeter: Chemical reaction: R P (volume may change).
HA = H(T1 ) + HP
(2.145)
HA = HR + H(T2 )
(2.146)
z
}|
{
H(T1 ) = HP = [CP (P) + CP (Cal)] (T2 T1 ) = CP,ef f (P)T
(2.147)
(2.148)
CP (P)
CP (R)
CP (Cal)
Note: All CP s above are extensive variables (i.e. not molar quantities). In practice,
CP (cal) dominates and CP,ef f (P ) CP,ef f (R).
100
If H has been determined in a calorimetric experiment, r H for a balanced
chemical reaction can be calculated using:
H
( = the extent of reaction; see Eq. (2.135))
(2.149)
r H =
=P 0
U
R=
(Ohms law)
I
(2.151)
P = U I (definition of power)
(2.152)
where q is the heat supplied by the electrical heater element, P is the power dissipated (constant; Watt), U is the applied potential (Volt; not internal energy),
I is the current (Ampere) and R is the heater resistance (Ohm). Alternatively, a
reference chemical reaction with known heat release can be used.
By recording data on q vs. T , it is possible to obtain the total heat capacity from
the slope, q/T , of this graph.
101
Constant volume calorimeter: Chemical reaction: R P (pressure may change).
At constant volume (V = 0), there cannot be any P V -work. This means that:
P V = 0
(2.153)
U = CV T
(2.154)
(2.155)
(2.156)
(2.157)
The difference between r H and r U is the P V -work. For an ideal gas this is
given by (Ns is the number of gaseous components, T constant, vi = ni /):
r H = r U +
Ns
X
(P V )
vg,i
r U + RT
i=1
(2.158)
102
where vg,i is the stoichiometric coefficient for gaseous product i (solids and liquids
do not contribute).
Notes:
In most experiments C(Cal) is large compared to C(R) or C(P), which
simply gives U = CV (Cal)T or H = CP (Cal)T .
Most chemists are interested in enthalpies because they describe best
chemical reactions in open vessels. Theoretical chemists, however, would be
more interested in internal energy because it is more directly related to the
energies of bond forming/breaking.
Example. If we pass a current of 10.0 A from a 12 V power supply for 300 s, what
is the amount of heat supplied?
Solution. The amount of heat dissipated by the heater is given by q = U It:
q = (10.0 A) (12 V) (300 s) = 3.6 104 AVs = 36 kJ
|{z}
= J
103
Solution. Since the volume is constant, we have:
U = CV (Cal)T CV (Cal) =
CV (Cal) = (1.247 g)
U
T
1 mol
2801 kJ
/ (1.693 K) = 11.44 kJ K1
180.16 g
1 mol
104
Several types of processes can be considered by calorimetric methods:
Solvation process:
HCl(g) + 5H2 O(l) HCl solvated by 5H2 O (Hsol (298 K) = 63.47 kJ mol1 )
Neutralization process: H-A + HO-B H2 O + A + B+
For example, A could be Cl and B Na. When strong acid/base combination is neutralized, the enthalpy difference is approximately constant and independent of the
base/acid pair. The reason is that the main contribution to the enthalpy difference
(heat of neutralization) originates from:
OH + H+ H2 O where r H (298 K) = 55.84 kJ mol1
When a dilute solution of weak acid/base is neutralized, the heat of neutralization
is less than that of strong acid/base neutralization.
Example. Calculate the enthalpy of formation of H+ and OH .
Solution. We use the known reactions to deduce the enthalpy of formation:
H2 O(l)
H2 (g) + 21 O2 (g)
H2 (g) + 21 O2 (g)
= H+ (ao) + OH (ao)
= H2 O(l)
= H+ (ao) + OH (ao)
r H = 55.84 kJ mol1
r H = 285.83 kJ mol1
r H = 229.99 kJ mol1
105
107
The first law of thermodynamics says only that the total energy is conserved but
it does not provide any information as to whether the process can proceed spontaneously.
Example. Consider gas expanding into vacuum. In practice, we know that the gas
will flow from the high pressure chamber into the low pressure chamber.
The first law for this system only states that the total energy is conserved but says
nothing about which way the gas would flow.
108
Example. In practice, we know that some chemical reaction proceed spontaneously
but some may require heat in order to proceed (nonspontaneous process)
Question: How do we identify a spontaneous process?
From the previous examples, it is obvious that a flow of heat is involved in these
processes. Heat is, however, not a state function (i.e. dq is not an exact differential).
This can be seen by inspecting Eq. (2.98):
U
U
dT + Pext +
dV
(3.159)
dqrev =
T V
V T
where we consider reversible heating and Pext = P . We will use Eq. (2.67) as the
test for exactness:
U
U
M =
and N = P +
(3.160)
T V
V T
and the required partial derivatives are then (U is well behaving, i.e. derivatives
exist and are continuous):
M
V
N
T
=
T
V
P
T
U
T
+
V
=
T
U
T V T V
U
V
(3.161)
109
M
P
=
(3.162)
V T
T V
V
= 0. This does not hold even for
would only be exact if P
T
N
T
nR
V
on the path and qrev is not a state function of the system. Thus it is not useful to
consider heat as an indicator for a spontaneity of the process.
Recall that some times division of an inexact differential by a suitable term (integrating factor) may result in exact differential. In this case it turns out that division
by T (temperature) yields an exact differential:
dqrev
(U/T )V
(U/V )T
P
=
+
dT +
dV
(3.163)
T
T
T
T
|
{z
}
|
{z
}
=M
=N
By taking the partial derivatives required in the exactness test and subtracting
(M/V )T from (N/T )V (denoted by ) we get:
1
U
1 P
P
2 2
(3.164)
=
T T V
T
T
V T
At least for ideal gases = 0 since (U/V )T = 0 (Joules experiment) and
P V = nRT . Thus dqrev /T is exact for ideal gases.
110
The above expression for can be derived as follows. First calculate (N/T )V :
N
T
1
T
=
V
1
T
P
T
=
=
P
T
T
P
T
+
V
P
+
T2
U
V
T
U
V T
V
P
1
U
U
1
2 2
+
T
T
V T
T T V T V
V
=
V T
V T
T V
T
U
1
1
=
=
T V
T V T
T T V T V
Subtracting the first expression from the second gives the difference :
P
1 P
1
U
2 2
=
T T V
T
T
V T
111
It will turn out that the dqrev /T is the quantity that we are looking for (entropy).
It is a useful quantity because it is a state function of the system (no path dependency). Note that we have only shown this to hold for ideal gases, a more detailed
consideration will be given later.
Example. Show that the values of the line integral for dqrev /T are identical along
the following paths (ideal gas, reversible processes, n = 1; one mole of gas):
Closed path 1:
Segment A: (T1 , P1 , V1 ) (T1 , P2 , V2 ) isothermal
Segment B: (T1 , P1 , V1 ) (T2 , P3 , V2 ) adiabatic
Segment C: (T2 , P3 , V2 ) (T1 , P2 , V2 ) constant volume
Closed path 2:
Segment A: (T1 , P1 , V1 ) (T1 , P2 , V2 ) isothermal
Segment D: (T1 , P1 , V1 ) (T3 , P1 , V2 ) constant pressure
Segment E: (T3 , P1 , V2 ) (T1 , P2 , V2 ) constant volume
112
Solution. Calculate dqrev /T along each segment. And sum the contributions from
each segment.
Path 1. (contribution from A + B + C = 0)
Segment A: The temperature is constant along this path. For an ideal gas, the
internal energy depends only on the temperature (Eq. (2.117)) and therefore dU =
0. Hence along this path (Eq. (2.64)): dqrev = dw = P dV = (RT /V )dV .
Division of both sides by T and integration gives:
S =
Zq1
q2
dqrev
=R
T
ZV1
V2
dV
= R ln
V
V1
V2
Segment B: Since the change is adiabatic, dqrev = 0 and the corresponding line
integral over dqrev /T is zero as well.
Segment C: Because volume is constant (constant volume process), dwrev = 0.
The first law (Eqs. (2.64) and (2.100)) now states that dqrev = dU = CV dT .
Division by T and integration along the segment gives:
S =
Zq2
q3
dqrev
=
T
ZT1
T2
CV
dT = CV
T
ZT1
T2
dT
= CV ln
T
T1
T2
113
Summing along the path B + C gives just CV ln(T1 /T2 ) and the contribution along
A is R ln(V2 /V1 ). Thus we have (using Eq. (2.128)):
I
I
dqrev
T1
V2
dS =
= CV ln
R ln
=0
T
T2
V1
Note that an alternative form for Eq. (2.128) is (considering now segment B):
T1
T2
V2
V1
R/CV
S =
Zq3
q1
dqrev
=
T
ZT3
T1
CP
dT = CP ln
T
T3
T1
114
Segment E: At constant volume, dw = 0. The first law of thermodynamics then
gives (Eqs. (2.64) and (2.100)):
S =
Zq1
q3
dqrev
=
T
ZU1
U3
dU
=
T
ZT1
T3
CV
dT = CV ln
T
T1
T3
(constant CV assumed)
The total contribution over the closed loop along D + E + A segments is:
I
I
T3
T1
V2
dqrev
= CP ln
+ CV ln
R ln
dS =
T
T1
T3
V1
V2
T3
R ln
= ln
(CP CV )
T1
V1
{z
}
|
=R (Eqs. (2.119, 2.120))
=R
V2
T3
ln
ln
T1
V1
In order to proceed, we note that the pressure at the end points of segment D is
constant. By using the ideal gas law for both end points, we get:
T3
V2
=
T1
V1
Thus the logarithm terms cancel and the integral is zero.
115
(3.165)
dS =
(3.166)
Rudolph Clausius,
Integration of entropy over closed loops yield zero because dS is an exact differential (S is a state function):
I
I
dq rev
=
dS = 0
T
(3.167)
(3.168)
0=
dS =
T
T
The inequality can also be written in differential form:
dq
dS
T
(3.169)
(3.170)
116
The idea behind Clausius inequality (3.170) can be understood by considering the
following example:
dS = dSc + dSh = dq
1
1
Tc
Th
(3.171)
117
The second law of thermodynamics consists of two statements:
1. There is a state function called the entropy S that can be calculated from
dS = dqrev /T .
2. The change in entropy in any process is given by dS dq/T , where the >
sign applies to a spontaneous (irreversible; dq irrev ) process and the equality
for a reversible process (dq rev ). In order to calculate S, one must use a
reversible process.
Justification for the Clausis inequality dS
dq
T
(Eq. (3.169)):
dq rev
T
.
dq
. Consider
2. If the process is irreversible, we need to show that dS > irrev
T
only P V -work and then the 1st law is dU = dq Pext dV . For a reversible
process this gives: dU = dq rev P dV and for an irreversible process:
dU = dq irrev Pext dV . Since dU is exact, the previous dU s must be equal
(consider integration over short paths): dq rev P dV = dq irrev Pext dV .
Rearranging gives: dq rev dq irrev = (P Pext )dV . If P Pext > 0 the
system will expand spontaneously and dV > 0. If P Pext < 0 the system
will contract spontaneously and dV < 0. In both cases dq rev dq irrev > 0.
dq
dq
> 0. By using the definition
Dividing boths sides by T gives Trev irrev
T
of entropy (Eq. (3.165)) we get: dS >
dq irrev
T
118
Another way to state the 2nd law of thermodynamics: The entropy increases in
a spontaneous process in an isolated system. The entropy increases as long as
spontaneous processes proceed. When the system does not change any more, the
entropy will have its maximum value and we have dS = 0. The entropy change tells
us whether a process or chemical reaction can occur spontaneously in an isolated
system. Consider an isolated system (consisting of system and surroundings):
119
Based on changes in entropy, we can identify three different cases:
(1) dS > dq/T spontaneous (irreversible) process
(3.173)
(3.174)
S1
dS =
Zq2
q1
dq rev
= S2 S1 = S
T
(3.175)
The integration path in Eq. (3.175) must be reversible. This equation can be
applied to irreversible processes only if a path consisting of reversible segments,
can be constructed. Note that there is no entropy change for a reversible adiabatic
process.
120
Example. Is the expansion of an ideal gas into a larger volume thermodynamically
spontaneous? More specifically, consider reversible and isothermal expansion of an
isolated ideal gas (n = 1 mol) initially at 298 K into a volume that is twice as large
as its initial volume.
Solution. Recall from Eq. (2.85) that the reversible work done is (n = 1 mol):
wrev =
ZV2
V1
Pext dV =
ZV2
V1
P dV =
ZV2
V1
nRT
dV = nRT ln
V
V2
V1
= RT ln(2)
The internal energy of an ideal gas does not depend on volume (Eq. (2.117)). Thus
we have U = qrev + wrev = 0 and further qrev = wrev = RT ln(2). Eq. (3.166)
with constant T states that S = qrev
= R ln(2) > 0 where we used the fact that
T
V2 = 2 V1 . Because the entropy change is positive, the change is spontaneous (as
we already knew in practice). For the reverse process we would have S < 0, which
means that it does not happen (unless forced).
121
Note: The previous problem has nothing to do with minimizing the energy, which is
constant during the process. The process is purely entropy driven and is related to
decrease in order at larger volume. By order we mean the arrangement of atoms
or molecules. For example, Sgas > Sliquid > Ssolid .
Example. Calculate the entropy change when argon at 25 C and 1.00 atm in a
container of volume 500 cm3 is allowed to expand to 1000 cm3 . Assume that argon
behaves according to the ideal gas law.
Solution. From the ideal gas law we can calculate the amount of substance:
n=
PV
= 0.0204 mol
RT
122
qrev
T
(3.176)
H
T
Examples of constant T and P processes are: vaporization of pure liquid into its
vapor at the equilibrium vapor pressure (P ), sublimation, and structural transitions
in solids.
Example. What is the change in molar entropy of n-hexane when it is vaporized
at its boiling point (68.7 C) under atmospheric pressure (1.01325 bar)? The molar
enthalpy of vaporization is 28850 J mol1 .
Solution. If n-hexane is vaporized into the saturated vapor at the given temperature, the process is reversible and the molar entropy change is given by Eq.
(3.176):
28850 J mol1
H
=
= 84.41 J K1 mol1
S =
T
341.8 K
Other useful special cases:
Constant V : Using Eq. (2.100) we have:
123
dq rev
CV dT
dS =
=
T
T
(3.177)
Integration of dS gives:
S =
ZT2
CV
dT CV ln
T
T1
T2
T1
ZT2
T1
CP
dT CP ln
T
(3.178)
(3.179)
T2
T1
(3.180)
Constant T and ideal gas: Following our previous ideal gas calculation, we have:
dwrev
P dV
dq rev
=
=
(3.181)
T
T
T
ZV2
ZV2
P
dV
V2
P2
S =
dV = nR
= nR ln
= nR ln
T
V
V1
P1
dS =
V1
V1
124
In terms of standard pressure and entropy (also molar quantities), the previous
expression can be written as:
P
( = 1 bar standard pressure)
(3.182)
S = S R ln
P
Example. Calculate the entropy change when argon gas at 25 C and 1.00 atm in
container of volume 500 cm3 is allowed to expand to 1000 cm3 and is simultaneously
heated to 100 C. Assume that argon behaves according to the ideal gas law.
Solution. The first part of the process was already considered in a previous example. The entropy change (S) due to volume change was 0.118 J K1 . In the
second step the gas is heated from 298 K to 373 K at constant volume (Eq. (3.178)):
V ln
S = nC
T2
T1
= (0.0204 mol) 12.48
373 K
J
ln
= 0.057 J K1
K mol
298 K
P C
V = R. The value of C
P can be
where CV was calculated from the relation C
found from the NIST chemistry webbook. Note the handling of the process in two
steps. The total change in entropy is the sum of the two: S = 0.175 J K1 .
Note: Because S is a state function (dS exact), we can choose any convenient path
for integration. In this case it was chosen as: (1) volume change and then (2)
temperature change.
125
Example. Half a mole of an ideal gas expands isothermally and reversibly at
298.15 K from a volume of 10 L to a volume of 20 L. (a) What is the change in
the entropy of the gas? (b) How much work is done on the gas? (c) How much
heat is transferred to the surroundings (qsurr )? (d) What is the change in the
entropy of the surroundings? (e) What is the change in the total entropy (system
+ surroundings)?
Solution. (a) Use Eq. (3.181):
S = nR ln
V2
V1
= (0.5 mol) 8.3145 JK1 mol1 ln (2) = 2.88 JK1
qrev
qrev = T S
T
and further:
wrev = T S = nRT ln
V2
V1
= 2.88 JK1 (298.15 K) = 859 J
126
To solve (c), we need qsyst = 859 J as calculated above. The total entropy is
conserved (Stot = 0) in a reversible process and thus we have Ssyst = Ssurr
and T Ssys = T Ssurr . Using Eq. (3.166) we have qsyst = qsurr and hence
qsurr = 859 J.
In (d), juse like in (c), the total entropy is conserved and hence Ssurr = Ssyst =
2.88 J K1 .
In (e), the total entropy is conserved: Stot = 0. The system and its environment
together can be considered as an isolated system. Also note that the process is
reversible.
Example. Consider that the expansion in the preceding example occurs irreversibly
by simply opening a stopcock and allowing the gas to rush into an evacuated bulb
of 10 L volume. (a) What is the change in the entropy of the gas? (b) How much
work is done on the gas? (c) What is qsurr ? (d) What is the change in the entropy
of the surroundings? (e) What is the change in the entropy of the system plus the
surroundings?
Solution. (a) The change in entropy is the same as in the previous example.
Entropy is a state function (dS exact) and hence its value does not depend on path.
(b) Expansion into vacuum does not involve P V -work.
(c) Since no work is involved, no changes in heat are involved (the first law).
(d) No heat exchange with the surroundings, thus the entropy of the surroundings
does not change.
(e) Stot = Ssyst + Ssurr = 2.88 J K1 + 0 J K1 = 2.88 J K1 . Since
the process is irreversible, we expect that the total entropy would increase (Eq.
(3.172)).
127
Summary: Calculation of S for various changes in state
In general, we have to find a reversible path in order to apply Eq. (3.166):
Z
dqrev
S =
T
Because S is a state function, the entropy change would then apply for all paths
(also for irreversible paths) that have the same initial and final states. Note that
dqrev is only defined for reversible paths.
Specific cases for one mole of substance:
Constant V :
substance(T1 , V ) substance(T2 , V )
T
R2 CV
dT (Eq. (3.178))
S =
T
T1
Constant P :
If the constant-volume
heat capacity CV is independent of T :
2
S = CV ln T
T1
substance(T1 , P ) substance(T2 , P )
T
R2 CP
dT (Eq. (3.180))
S =
T
T1
If the constant-pressure
heat capacity CP is independent of T :
2
S = CP ln T
T
1
128
Phase change at constant T and P :
solid(T , P ) liquid(T , P )
liquid(T , P ) gas(T , P )
solid(T , P ) gas(T , P )
S = H
(Eq. (3.176))
T
where H is the heat of vaporization,
sublimation or fusion (crystallization).
nA
nB
, yB =
nA + nB
nA + nB
Note that the entropy change is always positive (see Eq. (3.186)).
129
Example. Calculate the change in entropy of an ideal monatomic gas B in changing
from (P1 , T1 ) to (P2 , T2 ).
Solution. First we have to define a reversible path, for which we know how to
calculate the entropy change (each segment carried out reversibly):
B(T1 , P1 ) B(T2 , P1 ) B(T2 , P2 )
This path has two segments (one where temperature changes and another where
pressure changes). For the first step we have:
S = CP
ZT2
T1
dT
= CP ln
T
T2
T1
P2
P1
P ln
= nC
T2
T1
130
Example. Use the result of the previous example to calculate the molar entropy
change in the following process:
He (298 K, 1 bar) He (100 K, 10 bar)
Solution. Substitute the values to the equation in the previous examples:
S =
5
R ln
2
10 bar
J
T2
P2
5
100 K
R ln
= 41.84
R ln
= ln
T1
P1
2
298 K
1 bar
K mol
S = S(T
=
ZT
ZT
ZT
P
CP
a
aT 3
C
P = aT 3 )
dT =
dT = a T 2 dT = T 3 =
(C
T
T
3
3
0
) = S(0)
S(T
+
P
C
) = S(0)
S(T
+ 0.14 J K1 mol1
3
131
= 6004 J mol1 , C
liq = 75.3 J K1 mol1 , C
ice = 36.8 J K1 mol1 .
H
132
The system entropy is obtained as a sum over the three segments (the heat capacities
correspond to constant-pressure values; CP ):
Ssyst =
273
Z K
263 K
ln
273 K
263 K
263
Z
liq
C
H
dT +
+
T
T
273 K
ice
C
dT = 75.3 J K1 mol1
T
6004 J mol1
+ 36.8 J K1 mol1 ln
273 K
= 20.54 J K1 mol1
263 K
273 K
(3.183)
According to the statistical interpretation of thermodynamics, the decrease in entropy here corresponds to increased order in ice (i.e., molecules are more rigid in
the solid material than in the liquid).
Example. What is the change in entropy of the surroundings (the glass bottle plus
a large heat bath at 10 C) in the previous example? What happens to the total
entropy?
Solution. The process is (the gray area is the heat bath):
133
Because the surroundings (bottle plus bath) are large, their temperature does not
change remarkably during the heat transfer from the system. Heat transfer by an
infinitely small amount at constant temperature from the system to the surroundings
is a reversible process (i.e. it can be done the other way around as well). Thus,
for the surroundings, a reversible path is the direct path (#2 in previous diagram):
H2 O (l, 10 C) H2 O (s, 10 C).
We use the Eq. (3.176) to calculate the change in the entropy for the surroundings:
Ssurr =
263
qsurr
qsys
H
=
=
T
T
T
where we have to calculate the heat of fusion at 263 K for the system. In the
previous example it was given at 273 K. To do this, we use the same idea as in the
Eq. (2.140):
263
H
273
Z K
liq dT + H
273
C
263
Z K
ice dT
C
273 K
263 K
liq T + H
273
=C
K+
ice T = 75.3 J K1 mol1 (10 K)
C
| {z }
limits!
6004 J mol1 36.8 J K1 mol1 (10 K) = 5619 J mol1
134
Now we can calculate the entropy change (Eq. (3.176)):
Ssurr =
5619 J mol1
= 21.37 J K1 mol1
263 K
The total entropy change (system + surroundings) is the the sum of the two:
Stotal = Ssyst + Ssurr = 20.54 J K1 mol1 + 21.37 J K1 mol1
= 0.83 J K1 mol1 > 0
Notes:
The total change is positive, indicating a spontaneous (irreversible) process.
This is in line with Eq. (3.173).
The non-zero S arises from the differences in reversible paths (#1 and #2)
for the system and the bath. For the system (water/ice) the path #1 is
reversible and for the bath path #2 is reversible. Only reversible paths can
be used in calculating entropies. Along path #2 the temperature is constant.
Troutons rule: A wide range of liquids have approximately the same entropy of
vaporization (ca. 85 J K1 mol1 ). This is an empirical result.
135
136
The non-permeable membrane is moved towards the left end. At the same time the Gas
1 permeable membrane is moved towards the
left end. The Gas 2 permeable membrane is
kept at the center.
The Gas 1 permeable membrane is finally
moved all the way to the left. The nonpermeable membrane has been moved to the
center. Gas 2 and non-permeable membranes are on top of each other.
Dashed line
Dotted line
Continuous line
137
Is there any P V -work involved in (frictionless) movement of the membranes?
No, consider the following figure:
When the two membranes move to the left, there are equal pressures on both sides:
Gas 1 (P1 ) + Gas 2 (P2 ) and Gases 1 and 2 (P1 + P2 ). Both pressure and volume
are constant during the process and hence no work is done: w = 0.
138
Since we are dealing with ideal gases at constant temperature, the change in internal
energy (U ) is zero. Recall that for ideal gases the internal energy depends only on
temperature (see Eq. (2.117)). Now the first law of thermodynamics (U = q + w)
yields q = w = 0 (w = 0 from the previous page). No changes in heat also means
that there is no change in entropy and hence, for this segment, S = 0.
We still need to consider the first segment (step 1). Now Eq. (3.181) can be applied:
For gas 1: S1 = n1 R ln
V1
V
= n1 R ln
n1
n1 + n2
= n1 R ln (y1 ) (3.184)
For gas 2: S2 = n2 R ln
V2
V
= n2 R ln
n2
n1 + n2
= n2 R ln (y2 ) (3.185)
(3.186)
or, in general:
Smixing = R
Ns
X
i=1
ni ln (yi ) > 0
(3.187)
139
With just one gas molecule, after opening the stopcock, the molecule can be in
either bulb. Thus the number of equally probable arrangements is 2. With two gas
molecules, they can be located both on the left, one on each side (two possibilities)
and both on the right. This gives 4 equally probable arrangements. For n molecules
the number of equally probable arrangements is:
# of molecules
2
3
4
...
n
140
For one mole of molecules, N = 6.022 1023 (n = 1). The number of equally
probably arrangements is:
26.02210
Boltzmanns postulate:
23
= e4.17410
S = k ln ()
23
(3.188)
(3.189)
(3.190)
S = S S = k ln
Consider now the Joules experiment. In the initial state all molecules are in the
left bulb and there is only one possible arrangement, = 1. After opening the
valve and reaching the equilibrium the number of possible arrangments is given by
Eq. (3.188):
23
S = k ln
= 1.381 1023 J K1 4.174 1023 = 5.76 J K1
This is positive because entropy increases when volume increases. Note that if we
141
If we compare this result with the one we calculated earlier (expansion of ideal
gas), we note that they are identical (earlier we obtained R ln(2) 5.76). This
gives strong support for Boltzmanns postulate (but does not prove that it is correct).
Because of Eq. (3.169), S q/T , the entropy is a measure of heat flow between
the system and the surroundings. When heat is absorbed by the system from
its surroundings, q is positive, and the entropy of the system increases. On the
microscopic scale, the entropy is a measure of the dispersal of energy among the
possible microstates of molecules in a system (degrees of freedom).
Notes:
Since thermodynamics is a statistical theory, it works only when we have a
large number of atoms/molecules.
Some processes happen spontaneously even though they do not reduce the
energy of a system. They are driven purely by favorable change in entropy!
From the statistical point of view, it means that states with most degrees of
freedom are favored.
142
ST S0 K =
Tb
Tm
ZT
Z
(l)
(s)
Hfus Z C
Hvap
CP (g)
C
P
P
+
dT +
dT +
+
dT (3.191)
T
Tm
T
Tb
T
0 K
Tm
Tb
where Tm is the melting temperature and Tb is the boiling temperature. If temperature T is smaller than Tm or Tb , the corresponding terms are omitted. The
enthalpies have the following notation: f us (fusion; crystallization/melting) and
vap (vaporization). Note that all the heat capacities in this expression depend on
temperature:
143
Entropy of sulfur dioxide at different temperatures (Tm = 197.64 K, Tb = 263.08
f us = 7402 J mol1 , H
vap = 24937 J mol1 ):
K, H
T (K)
Method of calculation
S (J K1 mol1 )
0 15
15 197.64
197.64
197.64 263.08
263.08
263.08 298.15
P T 3 )
Debye function (C
Graphical, solid
f us /Tm
Fusion, H
Graphical, liquid
vap /Tb
Vaporization, H
P of gas
From C
S (298.15 K) S (0 K) =
1.26
84.18
37.45
24.94
94.79
5.23
247.85
144
T 0
(3.192)
T 0
(3.193)
145
Experimental details for the phase change in solid sulfur:
Rhombic form of sulfur is the stable form below the phase transition
temperature (368.5 K).
Monoclinic sulfur can exist below this temperature when supercooled.
Experimental methods for determining r S:
Method 1: Determine constant-pressure heat capacities (CP (T )) for both forms of
the crystal structures. Use Eq. (3.191):
rho
S368.5
S0rho
K =
368.5
Z K
0 K
mon
S368.5
mon
K S0 K =
rho
C
P
dT
T
368.5
Z K
0 K
(3.194)
mon
C
P
dT
T
If we assume that both rhombic and monoclinic forms have the same S0 K then
subtraction of the equations from each other gives the entropy difference between
the two forms at 368.5 K. The result from such calculation gives S(rho, mon) at
368.5 K, which is 1.09 J K1 mol1 .
146
Method 2: Determine heats of combustion for both crystal forms at 368.5 K and
calculate the enthalpy difference between the two forms by subtracting the two
values from each other. Then use Eq. (3.176):
S(rho,mon) =
H
401 J mol1
=
= 1.09 J K1 mol1
T
368.5 K
(3.195)
Both methods give consistent results. Note that we used the fact that both species
had identical entropies at 0 K. Statistical thermodynamics says that this value
should be zero.
Notes:
1. What is supercooling? A liquid below its melting point will crystallize in
presence of a seed crystal or nucleus around which a crystal structure can
form. However, lacking any such nucleus, the liquid phase can be
maintained all the way down to the temperature at which crystal
homogeneous nucleation occurs. For example, pure water can be cooled to
almost 231 K (melting point 273 K) when cooled very fast, about at rate 1
million K / second. Rain contains sometimes supercooled water, which
freezes immediately upon touching a surface.
2. CP goes to zero when temperature approaches 0 K.
3. According to the third law of thermodynamics, it is impossible to reach 0 K.
147
Other methods for determining entropy:
Measurement of equilibrium constant for a chemical reaction over a range of
temperatures yields both H and S .
Spectrophotometric measurements
Violations of the third law?
In most cases theory and experiment agree (298.15 K and 1 bar; gas):
Gas
CO2
NH3
NO2
CH4
C2 H2
C2 H4
C2 H6
N2 O
CO
S (Calc.) / J K1 mol1
213.8
192.8
240.1
186.3
200.9
219.6
229.6
215.1
193.5
S (Exp.) / J K1 mol1
213.7
192.6
240.2
186.3
200.9
219.6
229.6
219.9 ( = 4.8)
197.7 ( = 4.2)
This suggests that there is some kind of residual entropy in the system? When
comparing values obtained using heat capacities and spectroscopic data, similar
discrepancies can be noticed.
148
For both CO and N2 O, the residual entropy can be explained by imperfect crystal
structures at 0 K, which results in approximately a constant offset at elevated
temperatures as well (see previous table). Crystal imperfections contribute to the
total entropy and hence it does not appraoch zero even at 0 K. Consider the following
molecular arrangements in solid N2 O:
... NNO ONN NNO NNO ONN ...
In other words, the molecule can orient in many different ways in the crystal (O
O, O N, N N). If we consider that N2 O gas would consist of two different
species NNO and ONN with equal amounts and calculate the entropy of mixing
these two species, we get (Eq. (3.186) for per mole quantity, y1 = 1/2, y2 = 1/2):
1
mix S = n1 R ln (y1 ) n2 R ln (y2 ) = R ln
2
= 5.8 J K1 mol1
1
1
1
R ln
(3.196)
2
2
2
This is approximately the difference between the experimental and calculated results. If we could prepare a perfect crystal somehow, this correction would not
be needed. For CO, one needs to consider two different molecules CO and OC.
149
The following two sources of randomness in crystals at 0 K are not considered in
calculating the entropies for chemical purposes:
1. Possible isotopic mixtures of species. This is ignored because both the
reactants and the products contain the same mixture isotopes.
2. Spin degeneracy at 0 K is ignored. Again the same degeneracy exists in both
the reactants and the products. Spin is conserved in chemical reactions.
Notes:
1. The entropy of H+ (at chemical equilibrium) in water has been arbitrarily
assigned the value of zero.
2. When comparing standard entropy values from various sources, it is
important to be aware of the standard pressure used.
3. In order to see why NNO vs. ONN (or CO vs. OC) configurations result in
residual entropy, one should draw the crystal structure and look at the
possible orientations there.
4. H2 O is another system that has significant amount of residual entropy.
150
151
Consider the P V plot of the Carnot heat engine (the four cycles are denoted by
A, B, C, D, all processes are reversible):
1.
A B:
Heat q1 is transferred from the high-temperature (isothermal at
temperature T1 ) reservoir to the cylinder
Isothermal expansion of the gas pushes the piston towards larger
cylinder volume
The moving piston does work on the surroundings (w1 )
U1 = q1 + w1
2. B C:
Adiabatic expansion of the gas
The expansion continues until temperature of the gas drops from
T1 to T2
During this stage the piston does work on the surroundings (w2 )
U2 = w2
3. C
D:
Isothermal compression of the gas
The surroundings does work on the piston (w3 )
Heat flows out to the low-temperature reservoir (q2 )
U3 = q2 + w3
4. D
A:
Adiabatic compression of the gas
The surroundings do work on the piston (w4 )
No heat exchange
A steam engine with two stages
U4 = w4
152
Because the internal energy of the system is a state function, the total change in
internal energy over a closed path must be zero:
q
qout
in
z}|{ z}|{
(3.197)
Ucycle = q1 + q2 + (w1 + w2 + w3 + w4 )
|
{z
}
wcycle
{z
}
|
qcycle
(3.198)
wcycle = |qin | |qout | (qin > 0 and qout < 0)
(3.199)
In order to find out the efficiency of the heat engine, we define the efficiency parameter :
wcycle
|qin | |qout |
|qout |
=
=
=1
(3.200)
|qin |
|qin |
|qin |
Note that 0 < < 1 and larger value of correspond to better efficiency.
Next we apply the concept of entropy to simplify Eq. (3.200). Because the overall
cycle is reversible, the total entropy change over the closed cycle is zero. Only
two segments along the path deal with heat exchange (A-B and C-D). For these
segments Eq. (3.166) gives (Tin and Tout correspond to T1 and T2 , respectively):
|qin |
|qout |
|qout |
|qin |
=0
(3.201)
=
Scycle =
Th
Tc
Th
Tc
Tc
=1
(3.202)
Th
153
Thus the ratio of Tc /Th must be made as small as possible to achieve high efficiency.
Typically Tc would be room temperature, and therefore Th should simply be made
as high as possible.
155
dq
T
dq
T
When chemical potential is included above, we get the generalized form for dU :
dU = T dS P dV +
Ns
X
i dni
(4.203)
i=1
{z
chemical work
J. Willard Gibbs,
American physicist
(1839 - 1903)
156
Note that U depends on variables S, V and {ni }. These variables are called the
natural variables of U . The total differential of U can be written as (cf. Eq.
(1.45)):
dU =
U
S
V,{ni }
dS +
U
V
dV +
S,{ni }
Ns
X
U
i=1
ni
dni
(4.204)
S,V,{nj }j6=i
If this is compared with Eq. (4.203), we can see that partial derivatives in front of
dS, dV and dni must be equal to T , P and i , respectively:
U
U
U
,P =
, i =
(4.205)
T =
S V,{ni }
V S,{ni }
ni S,V,{nj }j6=i
Thus if we know the partial derivatives of the internal energy with respect to S, V
and ni , we can calculate T , P and i using Eq. (4.205).
If we allow the process to be irreversible, we have to consider the inequality in the
second law and keep in mind that the pressure now is really Pext :
dU T dS Pext dV +
Ns
X
i dni
i=1
(4.206)
157
(dU )S,V,{ni } 0
(4.207)
This is the criterion for spontaneous change and equilibrium at constant S, V and
ni . At equilibrium, U (S, V, {ni }) must be at minimum (i.e., dU = 0).
If T , P , and i are constant, Eq. (4.203) can be integrated:
U = TS PV +
Ns
X
i n i
(4.208)
i=1
Note: It is not convenient to use S and V as variables because they cannot be easily
controlled during chemical processes.
158
Change the differentials from (dx,dy) to (du,dy) with the following transformation:
g f ux
dg = df udx xdu = udx + vdy udx xdu = vdy xdu
g
where x = u
and v =
variables.
g
.
y
In a nutshell:
Transform the original differential in such a way
that the new differential depends on the conjugate
variables
Adrien-Marie Legendre,
French mathematician
(1752 - 1833)
159
Next, we will apply Legendre transformation to internal energy (ignore chemical
potential for now):
d U = P dV + T dS
|{z} |{z} |{z} |{z} |{z}
f
u dx
v dy
H = U (P V ) = U + P V
|{z} |{z} | {z }
g
ux
dH = T dS
V (dP ) = T dS + V dP
|{z} |{z} |{z} |{z} | {z }
dg
v dy
du
Now we have a new differential (enthalpy) dH with new natural variables S and P .
Note that the original differential dU had S and V as natural variables. Adding
chemical potential does not change this result since we were not operating on the
corresponding conjugate variables (i and ni ):
dH = T dS + V dP +
Ns
X
i dni
(4.209)
i=1
160
By comparing the terms the same way as we did for dU , we get the following
relations:
H
H
H
T =
,V =
, i =
(4.211)
S P,{ni }
P S,{ni }
ni S,P,{nj }j6=i
Thus, if we can determine the partial derivatives of H with respect to S, P and
ni , we can always obtain T , V and i from Eq. (4.211). A system under constant
S, P , and ni combined with Eq. (4.209) and exactly the same reasoning as in Eq.
(4.207) gives:
(dH)S,P,{ni } 0
(4.212)
Ns
X
i n i
i=1
(4.213)
161
Previously we established that U and H are connected to each other via Legendre
transformation with conjugate variables V and P . It is easier to control pressure
than volume and therefore H is more convenient to use in practice.
How about the other conjugate variable pair (T , S)?
Yes, it is more convenient to use temperature rather than entropy. Since (V , P )
pair offers two choices (U and H) and (T , S) another two, we have a total of four
different possibilities:
Quantity
Internal energy U
Enthalpy H
Helmholtz energy A
Gibbs energy G
Natural variables
S, V, {ni }
S, P, {ni }
T, V, {ni }
T, P, {ni }
Energy
U
H = U + PV
A = U TS
G = H TS
Differential (*)
dU = T dS P dV
dH = T dS + V dP
dA = SdT P dV
dG = SdT + V dP
(*) chemical potential should be added to each differential. We are not considering
Legendre transformation with respect to ni and i . The differential can also be
derived from the given energy expression by considering the total differential.
The last form is most useful in chemical applications since T and P can be controlled
(i.e. they can be held constant).
162
Expressions for the Helmholtz free energy (A):
dA = SdT P dV +
S=
P =
i =
A
ni
A
T
A
V
Ns
X
i dni
(4.214)
i=1
(4.215)
V,{ni }
(4.216)
T,{ni }
(4.217)
T,V,{nj }j6=i
(dA)T,V,{nj }j6=i 0
A = P V +
Ns
X
i n i
(4.218)
(4.219)
i=1
163
The corresponding expressions for the Gibbs energy (G) are:
dG = SdT + V dP +
Ns
X
i dni
G
T P,{ni }
G
V =
P T,{ni }
S=
i =
i=1
(4.220)
(4.221)
G
ni
(4.223)
T,P,{ni }
(dG)T,P,{ni } 0
(4.222)
G=
Ns
X
(4.224)
i ni (const. T, P, i ) (4.225)
i=1
S
ni
S
P
T,{ni }
P,T,{nj }j6=i
2G
P T
T,{ni }
2G
ni T
2G
T P
P,T,{nj }j6=i
P,{ni }
2G
T ni
V
T
P,{ni }
P,{ni }
(4.226)
i
T P,{ni }
(4.227)
164
V
ni
=
P,T,{nj }j6=i
2G
ni P
=
P,T,{nj }j6=i
2G
P ni
T,{ni }
i
P
T,{ni }
(4.228)
165
Example. Show that the Gibbs energy gives a criterion for spontaneous change at
constant T = Tsys = Tsurr and P .
Solution. We need to show that dGsys 0 dStot 0. First we calculate dGsys
by noting that G = H T S:
dGsys = d (Hsys T Ssys ) = dHsys d (T Ssys ) = dHsys Ssys dT T dSsys
| {z }
=0
dqsurr
+dSsys = dStot 0
T
| {z }
=dSsurr
Thus for spontaneous (irreversible) changes, the Gibbs energy will always decrease
with the equal sign applying only to reversible processes. Note that there is no need
to consider the surroundings explicitly when predicting spontaneity using Gibbs
energy whereas the calculation using entropy must consider both the system and
surroundings explicitly.
Note: Although the above criteria show whether a certain change is spontaneous,
it does not necessarily follow that the change will take place with an appreciable
speed.
166
When other than P V -work occurs in the system, it contributes a term to the fundamental equation for the internal energy (Eq. (2.95)). They can be included in
the Gibbs energy:
dG = SdT + V dP +
Ns
X
i dni + F dL + dAS
(4.229)
i=1
where F is the force of extension, L is the length, is the surface tension and AS is
the surface area. Variables (F and L) and ( and AS ) are conjugate variables and:
G
F =
(4.230)
L T,P,{ni },AS
G
(4.231)
=
AS T,P,{ni },L
What is the meaning of the Helmholtz and Gibbs energies in presence of non-P V
work?
1. Helmholtz energy. Consider the first law of thermodynamics (Eq.
dU = dq + dw. Using Eq. (3.172) gives:
dU + T dS dw
(2.64)):
(4.232)
(4.233)
167
This implies that:
(dA)T dw
(4.234)
where A is the Helmholtz energy. Thus a change in the Helmholtz energy gives an
upper bound for the work that can be done on the surroundings . In real processes
the amount work that can be done is less than |A|. In this case both sides of Eq.
(4.234) are negative (dw < 0 and (dA)T < 0); i.e. when the system does work on
the surroundings (i.e. |(dA)T | |dw|).
Note: dw now contains both P V and non-P V work.
2. Gibbs energy. The Gibbs energy is especially useful when non-P V work is involved. When P V and non-P V work are separated, the first law of thermodynamics
(Eq. (2.64)) can be written:
dU = dq Pext dV + dwnonP V
(4.235)
(4.236)
(4.237)
(dG)T,P dwnonP V
(4.238)
168
Notes:
1. For a reversible process at constant T and P the change in Gibbs energy is
equal to the non-P V work done on the system by the surroundings.
2. Eq. (4.238) states that a change in the Gibbs energy at constant T and P
gives an upper bound for the non-P V work that the system can do on its
surroundings. Remember that in this case dG and dwnonP V are negative.
3. When work is done on the system, the Gibbs energy increases. When work
is done by the system, the Gibbs energy decreases.
169
170
When looking at differences in the Gibbs energies and enthalpies, this becomes:
(G/T )
(4.242)
H = T 2
T
P,{ni }
Notes:
1. If G can be determined as a function of temperature, we can obtain H
using Eq. (4.242).
2. If H is independent of temperature and G is known at one temperature,
it is possible to obtain G at other temperatures with Eq. (4.242).
171
G1
dG =
ZP2
P1
V dP G2 = G1 +
ZP2
P1
V dP
|{z}
(4.243)
>0
Thus the Gibbs energy always increases with the increasing pressure when T and
{ni } are constant.
For liquids and solids volume (V ) is approx. independent of pressure (P ) and thus:
G2 = G1 + V (P2 P1 ) or G = V P
For ideal gases (P V = nRT ), we can write:
P2
P2
G2 = G1 + nRT ln
or G = nRT ln
P1
P1
Setting P2 = P, P1 = P , G2 = G and G1 = G , we get:
P
G = G + nRT ln
P
(4.244)
(4.245)
(4.246)
172
Example. Given the expression for the molar Gibbs energy (Eq. (4.246)) and
considering an ideal gas, derive the corresponding expression for the following thermodynamic properties: V, U, H, S and A at constant T .
Solution. The molar Gibbs energy is given by Eq. (4.246) with over bars:
=G
+ RT ln P
G
P
We need to find a similar -type expression for the above thermodynamic properties.
The form should be: X = X + possible pressure correction. For an ideal gas we
have:
RT
P V = RT V =
P
Based on Eq. (2.117), the internal energy of an ideal gas does not depend on volume
or pressure but only on temperature:
=U
U
Recall that enthalpy is defined as H = U + P V . For an ideal gas P V = constant
= nRT , thus:
=U
+ P V = U
+ RT = H
173
Eq. (3.182) for the molar entropy states that: S =
For the Helmholtz energy we have:
R ln
P
P
= G
P V = G
RT
A
=G
P V = G
RT = G
RT +RT ln P
+ RT ln P
A
=
A
| {z }
P
P
=A
174
Recall that Helmholtz energy gives the maximum amount of work that the system
can do in a reversible process. Here the work and Helmholtz energy should be equal.
Formally this can be seen integrating Eq. (4.214) at constant T :
= w = 5746 J mol1
dA = S dT P dV = P dV = dw A
|{z}
0
The internal energy of an ideal gas depends only on temperature and hence:
= 0 q = U
w = 5746 J mol1
U
Likewise, the enthalpy of an ideal gas depends only on temperature:
= U
+ P V = U
+ (RT ) = 0
H
|{z} | {z }
0
ZG2
G1
dG =
ZP2
P1
V dP =
ZP2
P1
RT
dP = RT ln
P
= 8.3145 J K1 mol1 (300.15 K) ln
1 bar
10 bar
P2
P1
= 5746 J mol1
175
Entropy for a reversible process at constant T can be calculated using Eq. (3.166):
S =
qrev
5746 J mol1
=
= 19.14 J K1 mol1
T
300.15 K
H
= 19.14 J K1 mol1
T
176
Example. Calculate the Gibbs energy of formation at 10 bar and 298.15 K for: (a)
Gaseous methanol (CH3 OH) (the standard Gibbs energy of formation 161.96 kJ
mol1 ) and (b) Liquid methanol (density 0.7914 g cm3 and the standard Gibbs
energy of formation 166.27 kJ mol1 ). Assume that gaseous methanol behaves
according to the ideal gas law.
Solution. The Gibbs energy of formation is defined essentially the same way as we
defined the enthalpy of formation previously. For (a) we can use Eq. (4.246):
f G = f G +RT ln
P
P
(298.15 K) ln
10 bar
1 bar
= 156.25 kJ mol1
Thus we see that the Gibbs energy of formation is higher at higher pressure. In
practice this means that it was not as favorable to form methanol at high pressure
than at the low pressure. In both cases the reaction is spontaneous.
For (b) the molar volume of the liquid is approximately independent of pressure.
In this case we can use directly Eq. (4.244). First we have to calculate to molar
volume:
177
MCH3 OH
32.04 g mol1
=
= 40.49 cm3 mol1
V =
CH3 OH
0.7914 g cm3
3
= 40.49 cm3 mol1 102 m cm1 = 40.49 106 m3 mol1
f G = f G + V (P P ) = 166.27 kJ mol1 + 40.49 106 m3 mol1
9 105 Pa
103 J / kJ
= 166.23 kJ mol1
178
Gilbert Newton
Lewis, American
physical chemist
(1875 - 1946)
For example, the Gibbs energy (Eq. (4.246)) in terms of fugacity is:
=G
+ RT ln f
G
P
(4.248)
Note that fugacity has the units of pressure. Also, for a real gas, it is directly related
to the Gibbs energy. If the equation of state for a real gas is know, it is possible to
calculate the fugacity.
179
Derivation of expression for fugacity for a real gas at constant temperature:
First recall Eq. (4.222): V = G
P
T,{ni }
The Gibbs energy differential between the real and ideal gases is:
G
id = V V id dP
d G
Integration of both sides gives:
G(P
Z )
ZP
G
id =
V V id dP
d G
(4.249)
G(P )
G
id
G
G
id
G
G
id and:
Letting P 0 makes G
G
id
G
ZP
0
ZP
V V id dP
V V id dP
(4.250)
(4.251)
180
Using Eq. (4.248) for the real gas, Eq. (4.246) for the ideal gas and inserting them
into Eq. (4.251) we have:
ln
f
P
1
=
RT
ZP
0
is the same for both gases)
V V id dP (note: G
(4.252)
ZP
f
1
id
dP
V V
= exp
P
RT
(4.253)
where is called the fugacity coefficient. Next we introduce the molar volumes for
the real and ideal gases (see Eqs. (1.6) and (1.26)):
RT Z
RT
V =
and V id =
P
P
When these are inserted into Eq. (4.253), we get:
P
Z
ZP
RT Z
Z1
1
RT
f
dP = exp
= exp
dP
=
P
RT
P
P
P
0
(4.254)
181
Example. What is the expression for the fugacity when the compressibility factor
is expanded as a power series?
Solution. The power series expansion is given by Eq. (1.17):
Z = 1 + B P + C P 2 + ...
Inserting this into Eq. (4.254), we get:
f = P exp
= P exp
ZP
0
ZP
0
Z
Z1
B P + C P 2 + ...
dP = P exp
dP
P
P
0
CP 2
B + C P + ... dP = P exp B P +
+ ...
2
Example. What is the expression for the fugacity of a van der Waals gas?
Solution. The compressibility factor of a van der Waals gas is given by Eq. (1.26):
Z =1+
1
a
b
P (to first order in P )
RT
RT
182
Just like in the previous example, we insert the expression for Z into Eq. (4.254):
f = P exp
ZP
0
ZP
a
1
Z1
dP = P exp
b
dP
P
RT
RT
0
P
a
= P exp
b
RT
RT
Example. Given that the van der Waals constants of N2 are a = 1.408 L2 bar
mol2 and b = 0.03913 L mol1 , estimate the fugacity of the gas at 50 bar and 298
K.
Solution. Insert the constants into the above expression (previous example):
P
a
f = P exp
b
RT
RT
!
(
50 bar
= (50 bar) exp
0.083145 L bar K1 mol1 (298 K)
#)
"
1.408 L2 bar mol2
1
0.03913 L mol
183
Note: All thermodynamic tables refer to the ideal gas pressure P instead of the
fugacity. In general, the fugacity can be either smaller (at low P ) or higher than
(at high P ) the true pressure P . This can be seen by inspecting the graph shown in
the first chapter where we plotted Z as a function of P . There Z may take values
on either side of 1, which may make the exponent term in the fugacity expression
(Eq. (4.254)) either negative, zero, or positive.
Even though we have not discussed the chemical potential in detail yet, we conclude
that the same idea can be applied to the chemical potential as well (cf. Eq. (4.248)):
i = i + RT ln (ai )
(4.255)
where i is the chemical potential and ai is the activity. The activity is dimensionless, and ai = 1 in the reference state for which i = i .
184
For real gases:
ai =
fi
P
ai =
Pi
P
ai 1
If the molar volume is assumed to be constant (i.e., liquids and solids), we can
integrate one of the Maxwell equations (Eq. (4.228)):
i
= Vi
P T,{ni }
Integration of this expression yields:
i (T, P ) = i + Vi (P P )
(4.256)
Note that we assumed that T and V are constants during the integration over P .
Comparison of Eqs. (4.256) with (4.255) shows that:
RT ln (ai ) = Vi (P P )
a = eVi (P P )/(RT )
(4.257)
Later we will see that when dealing with solutions, we write activity as a product
of activity coefficient and concentration.
185
Recall that we did not carry out Legendre transformation with respect to ni and i .
Thus U , H, A and G behave exactly the same way with respect to these variables.
Collecting the results from Eqs. (4.205), (4.211), (4.217) and (4.223):
U
H
=
ni S,V,{ni }i6=j
ni S,P,{ni }i6=j
A
G
=
=
ni T,V,{ni }i6=j
ni T,P,{ni }i6=j
i =
(4.258)
The last definition based on the Gibbs energy is the most convenient because it is
easy to hold T and P constant in practical applications.
Consider a single species which is distributed in two different phases. An example of
such system is liquid water and water vapor. Let us denote the amount of substance
in phase as n and the amount of substance in phase as n . The corresponding
chemical potentials are denoted by and . When an infinitesimal amount of
substance is transferred between the phases, we have dn = dn dn (what
leaves must enter and vice versa). According to Eq. (4.258), differential of the
Gibbs energy is (constant T and P ):
(4.259)
(dG)T,P = dn + dn = dn
186
Recall that for a spontaneous process dG < 0. Consider five different cases:
1. dn > 0 (i.e. substance moves from to ) and < . This fulfills the dG
requirement for a spontaneous process. Substance moves from the phase
with higher chemical potential to the phase where the chemical potential is
lower.
2. dn > 0 (i.e. substance moves from to ) and > . dG would be
positive and therefore this process is not spontaneous.
3. dn < 0 (i.e. substance moves from to ) and < . dG would be
positive and therefore this process is not spontaneous.
4. dn < 0 (i.e. substance moves from to ) and > . This fulfills the dG
requirement for a spontaneous process. Substance moves from the phase
with higher chemical potential to the phase where the chemical potential is
lower.
5. dn = 0 (i.e. nothing transferred) or = (i.e. chemical potentials equal).
This corresponds to equilibrium (dG = 0). Such an equilibriums can exist
even if the phases are at different pressures.
187
Measurement of i as a function of P and T can be used to determine the molar
entropy and the molar volume. To do this, we use the Maxwell relations (derivation
is similar to Eqs. (4.226) - (4.228)):
S
i
Si =
(4.260)
=
ni
T P,{ni }
i
V
(4.261)
=
Vi =
ni
P T,{ni }
Example. Consider a mixture of ideal gases. What are the expressions for the
chemical potentials and the partial molar entropies?
Solution. For simplicity we carry out the calculation for two species. Extension
to many species is straight forward. We will show this in three parts:
1. Prove that the molar volume of the mixture (V ) is equal to the partial molar
volumes of the components (V1 and V2 ).
2. Use Eq. (4.228) to obtain (i /Pi ).
3. Integrate the resulting expression to obtain i .
188
Step #1. First we note that the total number of gas atoms/molecules is given by
n = n1 + n2 , the total pressure P = P1 + P2 where each gas follows its own ideal
gas law Pi V = ni RT . There is no interaction between gas atoms/molecules in an
ideal gas and therefore we can add up
these individual ideal gas equations to get
all
V
V = nRT
. Now Eq. (1.43) Vi = n
= V .
= RT
P
P
i
T,P,{nj }j6=i
Step #2. By using Eq. (4.261), the above result, the chain rule and Pi = xi P , we
can write:
RT
V =
= Vi =
P
i
P
i
i
Pi
= xi
Pi
P T,{ni }
Pi T,{ni }
T,{ni }
RT
i
=
(4.262)
Pi T,{ni }
Pi
=
di = RT
ZPi
i = i + RT ln
dPi
Pi
Pi
P
(4.263)
189
By differentiating Eq. (4.263) with respect to T and using Eq. (4.260) we get:
Pi
(4.264)
Si = Si R ln
P
The same equation for pure ideal gas was derived earlier: S = S R ln PP .
Note: The chemical potential is one of the most important concepts in chemical
thermodynamics. In both chemical reactions and phase changes, the chemical potential of a species times its differential amount (reacted or transferred) determines
the change in U , H, A, or G, depending on the variables that are held constant
during the process.
190
Ns
X
n i i =
i=1
Ns
X
G
ni
i=1
T,P,{nj }i6=j
ni =
Ns
X
i
ni G
i=1
This is in agreement with the definition (i.e. Eqs. (1.42) and (1.43)). Note that T
and P are constant.
2. Partial molar entropy
Differentiation of both sides of the above equation with respect to T gives:
Ns
X
G
i
ni
=
T P,{ni }
T P,{ni }
i=1
{z
}
{z
}
|
|
=S (Eq. (4.221))
S=
Ns
X
i=1
ni
S
ni
S , Eq. (4.227)
= n
i
Ns
X
T,P,{nj }j6=i
Si ni
i=1
191
3. Partial molar volume
Eq. (4.222) reads: V = G
P
T,{ni }
V
ni
T,P,{nj }j6=i
i
P
T,{ni }
Differentiation of the Gibbs energy equation on the previous page with respect to
P gives (the Gibbs energy Eq. differentiated; cf. Eqs. (1.42) and (1.43)):
Ns
Ns
X
X
G
i
V =
ni
Vi ni
(4.265)
=
=
P T,{ni } i=1
P T,{ni } i=1
4. Partial molar enthalpy
Recall that: H = G + T S. Since both G and S are additive:
H = G + TS =
Ns
X
i + T
ni G
Ns
X
i=1
ni
(G + T S)
ni
ni Si =
T,P,{nj }j6=i
Ns
X
ni
i=1
i=1
i=1
Ns
X
Ns
X
i=1
ni
G
S
+T
ni
ni
H
ni
(4.266)
T,P,{nj }j6=i
T,P,{nj }j6=i
192
How do ideal mixtures of gases behave?
By ideal mixture we mean that the gases do not interact with each other and hence
they have their independent equations of state. Consider an alternative form of Eq.
(4.263):
Pi
yi P
i = i + RT ln
=
+
RT
ln
(4.267)
i
P
P
where yi is the mole fraction and P is the total pressure (this will be used below).
1. The total Gibbs energy of an ideal gas mixture
G=
Ns
X
n i i =
i=1
i=1
+RT ln
Ns
X
P
P
ni
X
Ns
Ns
X
yi P
i + RT ln
ni ln (yi ) (4.268)
ni i + RT
=
P
i=1
i=1
Ns
X
Ns
Ns
X
X
yi ln (yi ) + RT ln
yi i + RT
ni = n
yi
= nG
i=1
i=1
i=1
i=1
| {z }
X
Ns
193
2. The total entropy of an ideal mixture
G
(4.269)
S=
T P,{ni }
=n
(Eq. (4.260))
=S
i
z }| {
}|
{
z
Ns
X
Ns
Ns
X
X
i
P
=
+R
ni
ni ln (yi ) + R ln
ni
T
P
i=1
i=1
i=1
= n
Ns
X
i=1
yi Si R
Ns
X
i=1
yi ln (yi ) R ln
P
= nS
P
At constant P and T , H = H and Eqs. (4.268), (4.269), (4.258) and (1.43) give:
H = H = G + T S =
Ns
X
ni i + T Si
i=1
Ns
X
i=1
ni
G
S
+T
ni
ni
=
T,P,{nj }j6=i
Ns
X
i=1
(4.270)
ni
H
ni
Ns
X
ni H
i
i=1
194
4. The volume of an ideal gas mixture
Direct calculation gives (n = n1 + ... + nNs , P = P1 + ... + PNs ):
X
Ns
Ns
Ns
X
X
RT
V
nRT
ni
ni Vi
=
ni =
=
V =
P
P
ni
i=1
i=1
i=1
|{z}
| {z }
(V /ni )
(4.271)
i
=V
What are the initial values and the changes in Gibbs energy, entropy, enthalpy, and
volume?
195
Solution. The initial values (i.e. before mixing) for the thermodynamic quantities
in question are given by (separate systems):
Eq. (4.263)
P
P
+
n
+
n
RT
ln
2
2
2
P
P
Eq. (4.264)
z}|{
P
n2 ln P
= n1 S1 + n2 S2
=
n1 S1 n1 R ln
+
n
S
2
2
P
P
Gini = n1 1 + n2 2
Sini
z}|{
=
n1 1 + n1 RT ln
+ n2 H
Hini = n1 H
1
2
Vini = n1
RT
RT
+ n2
P
P
Gf in
Eq. (4.268) z
z}|{
=
n1 1 + n1 RT ln
}|
{
P
P
+
n
+
n
RT
ln
+n1 RT ln (y1 )
2 2
2
P
P
<0 (spontaneous)
<0
<0
z }| {
z }| {
z }| {
mix G
+n2 RT ln (y2 ) = Gini + RT n1 ln (y1 ) +n2 ln (y2 ) = Gini +
196
Eq. (4.269)
Sf in
z}|{
=
n1 S1 n1 R ln
|
P
P
+ n2 S2 n2 R ln
{z
P
P
=Sini
Eq. (4.270)
Hf in
z}|{
=
+ n2 H
= Hf in mix H = 0
n1 H
1
2
z}|{
=
n1
Eq. (4.271)
Vf in
RT
RT
+ n2
= Vini mix V = 0
P
P
197
In this section we derive the Gibbs-Duhem equation, which we will need later. First
we use Eq. (4.208) constant T and P :
U (V, S, {ni }) =
Ns
X
i=1
ni i P V + ST
Ns
X
ni di +
Ns
X
i=1
i=1
i dni P dV V dP + SdT + T dS
On the other hand, since U = U (V, S, {ni }) the total differential of U can be
expressed in terms of the partial derivatives of U (Eq. (4.203)):
dU =
Ns
X
U
ni
i=1
dni +
V,S,{nj }j6=i
{z
=i (Eq. (4.205))
}
Ns
X
i=1
U
V
dV +
S,{ni }
{z
=P (Eq. (4.205))
i dni P dV + T dS
U
S
dS
V,{ni }
{z
=T (Eq. (4.205))
198
By subtracting the two equations from each other yields:
Ns
X
i=1
ni di V dP + SdT = 0
(4.272)
This is known as the Gibbs-Duhem equation. We will use it later when we discuss
the Clapeyron equation as well as Henrys and Raoults laws.
Notes:
N
Ps
ni di = 0.
i=1
2. For a system with two species at constant T and P we have (y1 and y2 are
the mole fractions for species 1 and 2):
y1 d1 + y2 d2 = 0 (with y2 = 1 y1 )
y1 d1 + (1 y1 )d2 = 0
(conservation of chemical potential)
(4.273)
199
First we consider Maxwell relations for a pure system (i.e. only one species). The
equations can be derived using U , H, A or G:
Eq. (4.205)
2U
V S
Eq. (4.211)
2H
P S
Eq. (4.215)
2A
V T
T
V
T
P
S
V
S
P
z}|{
=
z}|{
=
z}|{
=
Eq. (4.221)
z}|{
=
OK
2U
V
S
Eq.
x-change OK
2H
SP
Eq.
2A
T V
Eq.
x-change
z}|{
=
z}|{
=
x-change
z}|{
=
2G
P T
OK
x-change
z}|{
=
OK
2G
T P
(4.205)
P
S
(4.211)
V
S
(4.216)
P
T
z}|{
=
z}|{
=
z}|{
=
(4.274)
(4.275)
(4.276)
Eq.(4.222)
z}|{
V
(4.277)
=
T P
200
When a number of different species are present, the chemical potential must be
included. For example for the internal energy U , the following Maxwell relations
can be written (similar relations exist also for H, A and G):
T
i
(4.278)
=
ni S,V,{nj }j6=i
S V,{ni }
P
i
(4.279)
=
ni S,V,{nj }j6=i
V S,{ni }
j
i
(4.280)
=
nj S,V,{n }
ni S,V,{nk }k6=i
k k6=j
/ V
Example. Calculate U
/ V
= 0 for an ideal gas. If the
Solution. Earlier we have shown that U
T
equation of state of the real gas is known in terms of P , we will be able to use the
following equation to calculate the partial derivative:
U
P
P
(4.281)
=
T
V T
T V
201
To show that this result holds, we first combine the 1st and 2nd laws of thermodynamics (reversible process):
1st law: dU = dq rev + dwrev = dq rev P dV (if only PV -work)
2nd law: dS =
dq rev
T
Combined: dU = T dS P dV
By considering molar quantities in above, dividing both sides of the equation by
dV , and imposing constant temperature gives:
U
S
=
T
P
(4.282)
V T
V T
Now Eq. (4.276) allows us to write this as:
P
U
=
T
P
V T
T V
Example. This equation can be applied to a van der Waals gas. In this case the
pressure can be written as (Eq. (1.24)):
a
RT
2
P =
V b
V
202
Differentiation of P with respect T gives:
P
T
R
b
V
P
=
=
V
V b
V b
V b
V 2
derivative:
a
V 2
(4.283)
Integration of this equation yields the change in internal energy for a given change
in volume at constant temperature:
=
U
ZU2
1
U
=
dU
ZV2
1
V
a
dV = a
V 2
1
1
V1
V2
203
Then we can use the expression for U on the previous page:
=a 1 1
U
= 0.8779 Pa m6 mol2
V1
V2
1
1
= 58.5 J mol1
10 103 m3 mol1
30 103 m3 mol1
Example. What is the change in the molar entropy when a van der Waals gas
expands isothermally?
Solution. First we use one of the Maxwell relations (Eq. (4.276)) and then we
integrate both sides of the resulting equation (see also the previous example):
ZS2
1
S
S
V
dS = R
Eq. (4.276)
ZV2
1
V
z}|{
=
P
T
R
=
V b
dV
S = R ln
V b
V2 b
V1 b
204
Cubic and isothermal expansion coefficients
The cubic expansion coefficient is defined as (how the volume of a substance
changes with temperature):
1 V
1 V
=
=
(4.284)
V
T P
T P
V
The isothermal compressibility is defined as (how the volume of a substance changes
with pressure):
1 V
1 V
=
=
(4.285)
V
P
P
V
T
For an ideal gas these quantities are = 1/T and = 1/P . These constants can
also be used in simplifying thermodynamic expressions. For example:
V /T P
P
=
=
(4.286)
T V
V /P T
and further:
U
V
Eq. (4.281)
T
z}|{
=
T P
(4.287)
205
Dependency of enthalpy on the pressure at constant temperature
We are seeking an expression like Eq. (4.281) that would help us to calculate the
dependency of enthalpy on pressure. Both sides of the differential dH = T dS +V dP
can be divided by pressure at constant T to obtain:
S
H
=T
+ V
(4.288)
P T
P T
Using one of the Maxwell relations (Eq. (4.226)), we can modify this equation to:
V
H
= T
+ V
(4.289)
P T
T P
The relation between constant pressure and constant volume heat capacities
In Eq. (2.116) we found the following relation:
V
P C
V = P + U
C
V T
T P
From Eqs. (4.284), (4.285) and (4.287) we can now get:
2
P C
V = T V
C
(4.290)
(4.291)
207
Experiments: Chemical reactions can approach equilibrium from either direction (CG & PW):
A+BC+D
Norwegian chemists: Cato
Guldberg (1836 - 1902)
and Peter Waage (1833 1900)
N
Ps
i dni
i=1
Remember that most chemical reactions are carried out at constant pressure and
temperature and hence dT = dP = 0. Next, recall Eq. (2.135):
ni = ni,0 + vi
(5.292)
208
If one chemical reaction is considered, substitution of the differential into Eq. (4.220)
gives:
Ns
X
vi i d
(5.293)
dG = SdT + V dP +
i=1
=
T,P
Ns
X
vi i
(5.294)
i=1
T,P
=0
vi i,eq = 0
(5.295)
i=1
This is equilibrium condition applies to all chemical equilibria (gases, liquids, solids,
or solutions).
In similar fashion to Eq. (2.137), we define the reaction Gibbs energy r G:
G
(5.296)
r G =
T,P
209
By combining the equilibrium condition, Eqs. (5.294) and (5.296):
r G =
Ns
X
vi i
(5.297)
i=1
The reaction Gibbs energy at a given value of , tells us how much the Gibbs energy
changes when is varied. Its value approaches zero when close to the equilibrium.
Consider reaction: A + B C + D
Note: Usually does denotes a difference but here it denotes a derivative. However,
there is a close connection to differences in chemical potentials via Eq. (5.297).
210
To proceed towards the concept of equilibrium constant, recall Eq. (4.255):
i = i + RT ln (ai )
where i is the standard chemical potential and ai is activity of species i. In
equilibrium, this equation can be written as (ai,eq = ai (eq )):
i,eq = i + RT ln (ai,eq )
(5.298)
Ns
X
vi i = RT
vi ln (ai,eq )
(5.299)
i=1
Ns
X
vi i = RT
i=1
vi
ln ai,eq
Ns
Ns
X
Y
v
i
vi i = RT ln
ai,eq = RT ln(K)
i=1
{z
(5.300)
(5.301)
i=1
=r G
{z
Note: r G is a plain number, which comes from the standard chemical potentials
of the isolated species A, B, C, D whereas r G is a function depending on . The
value of r G depends only on the species which react but not on the details of
the reaction itself. r G describes the reaction.
211
Thus we can summarize the previous result as:
r G = RT ln (K)
(5.302)
The equilibrium constant K is directly related to the standard reaction Gibbs energy r G . Because r G is a function of temperature only (pressure fixed to 1
bar), the equilibrium constant K depends only on temperature. Note that K is a
dimensionless quantity.
Lets consider a non-equilibrium condition (i.e. r G 6= 0) and rewrite Eqs. (5.297)
and (5.298) as:
Ns
Ns
Ns
Y
X
X
vi
ai
vi ln(ai ) = r G +RT ln
vi i +RT
r G =
= r G +RT ln (Q)
i=1
i=1
i=1
{z
(5.303)
where we have defined reaction quotient Q. It is a very similar quantity to K but
it is a dynamic variable that describes a non-equilibrium situation. At equilibrium,
Q = K and in general Q = Q() whereas K does not depend on . By calculating the
reaction quotient Q, it is possible to determine the direction in which the chemical
reaction would proceed at given . Substitution of the definition of the Gibbs energy
(G = H T S) into Eq. (5.296) gives:
S
H
T
= r H T r S
(5.304)
r G =
T,P
T,P
212
Note that the previous definition is in agreement with Eqs. (2.137). The same
equation also applies when the reactants and products are in their standard states:
r G = r H T r S
(5.305)
K=
Q() =
aCH4 a2CO
a3C a2H2 O
aCH4 a2CO
a3C a2H2 O
Also Q (eq ) = K.
213
Recall that for real gases the activity is given by (see Eqs. (4.248) and (4.255)):
ai = Pfi where P is the standard state pressure (1 bar) and fugacity is given, for
example, by Eq. (4.254):
ZPi
Z
fi = Pi exp
dP
P
0
and Z is the compressibility factor of the gas. At equilibrium Eq. (5.301) now gives:
Ns
Ns
Y
Y
fi,eq vi
vi
ai,eq
=
(5.306)
K=
P
i=1
i=1
For ideal gases Eq. (5.306) can be written as follows (ai = Pi /P , see Eq. (4.255)):
Ns
Ns
Y
Y
Pi,eq vi
Pi vi
and
Q
=
(5.307)
K=
P
P
i=1
i=1
Notes:
K depends only on partial pressures Pi and also on the total pressure P .
As K depends on the way the stoichiometric coefficients are written, a full
reaction equation must always be written when specifying K.
K determines eq (the equilibrium extent of the reaction) via Eq. (5.301).
214
Example. A mixture of CO(g) and CH3 OH(g) at 500 K with PCO = 10 bar,
PH2 = 1 bar, and PCH3 OH = 0.1 bar is passed over a catalyst. The reaction
occurring in the gas phase is:
CO(g, 10 bar) + 2H2 (g, 1 bar) = CH3 OH(g, 0.1 bar)
where the standard reaction Gibbs energy is r G = 21.21 kJ mol1 . In which
direction would the reaction proceed? Assume that the gases behave according to
the ideal gas law.
Solution. Recall that the the reaction Gibbs energy (r G) tells us which way the
reaction would proceed. If it is negative, the reaction proceeds spontaneously from
left to right and if it is positive, the reaction proceeds spontaneously from right
to left. If r G is zero, the reaction is in equilibrium. By using Eqs. (5.303) and
(5.307) we get:
r G = r G +RT ln (Q) = 21.21 kJ mol1 + 0.0083145 kJ K1 mol1 (500 K)
ln
0.1
10 12
Thus the reaction is not spontaneous as written. Note that the only effect of P s
in this calculation were to cancel the units since its numerical value was 1 bar. The
stoichiometric coefficients are: vCO = 1, vH2 = 2 and vCH3 OH = 1.
215
If the conditions had been chosen as follows:
CO(g, 1 bar) + 2H2 (g, 10 bar) = CH3 OH(g, 0.1 bar)
we would have:
r G = 21.21 kJ mol1 + 0.0083145kJ K1 mol1 (500 K) ln
= 7.51 kJ mol1 < 0
0.1
1 102
1
O2 (g)
2
The standard reaction Gibbs energy for this reaction is 118.08 kJ mol1 at 2300 K.
What is the extent of reaction eq at 2300 K and at a total pressure of 1.00 bar?
Assume that gases behave according to the ideal gas law.
Solution. First we use Eq. (5.302) to get the equilibrium constant K:
ln (K) =
118080 J mol1
r G
= 6.175 K = 2.08103
=
RT
8.3145 J K1 mol1 (2300 K)
216
In the following, the initial amount of water is denoted by n and initially all
other concentrations are zero. At given extent of reaction , the amount of each
component is: nH2 O = (1 )n, nH2 = n and nO2 = n/2. Thus the total
amount is (1 ) n + n + n/2 = n (1 + /2) and the mole fractions are then:
yH2 O = (1 ) / (1 + /2), yH2 = / (1 + /2) and yO2 = (/2) / (1 + /2). The
partial pressures (Pi ) can be calculated by multiplying the corresponding mole fraction by the total pressure P . Note that the choice of is slightly different than in
Eq. (5.292) (scaling = n so that is now between 0 and 1). Then apply Eq.
(5.307) to relate the partial pressures to K:
K=
3/2
PO
eq P 1/2
2 =
where P = P/P = (1 bar) / (1 bar) = 1
PH2 O P
(1 eq ) (2 + eq )1/2
PH2
3/2
K=
eq
(1 eq ) (2 + eq )1/2
(there are actually three roots: one real and two complex)
Therefore about 2% of H2 O is dissociated under the given conditions.
Example. Demonstrate that chemical reactions in the gas phase never go to completion. Assume ideal gas behavior as well as constant pressure and temperature.
Solution. Consider reaction A(g) = B(g) at constant pressure. First calculate the
Gibbs energy of the system as a function of the extent of reaction (Eq. (5.293)):
217
dG = S dT +V dP +
|{z}
|{z}
=0
=0
Ns
X
i=1
vi i () d = B ()d A ()d
A = A + RT ln
=
+
RT
ln
A
P
P
P
P
+
RT
ln
(1
)
+
ln
ln
(y
)
= A + RT
+
ln
A
A
P
P
| {z }
=(1)ntot /ntot
B = B + RT
P
ln () + ln
d
BA +RT ln
dG =
| {z }
1
=B A
218
Integration of this equation from the initial point ( = 0) to gives:
G()
Z
dG =
Z
0
BA + RT ln
G(0)
G() = BA +RT
Z
ln
=mix G
mol1 ):
219
220
Consider a chemical reaction: aA + bB cC + dD. If the initial concentrations/pressures [A]0 , [B]0 , [C]0 and [D]0 are known and one of [A], [B], [C] or [D] is
known at equilibrium, it is possible to use the above balanced chemical equation to
obtain concentrations of all of the species at equilibrium. This is because there is
just one variable that describes the reaction (; the extent of reaction). Given the
equilibrium concentrations/pressures for all the species, it is possible to derive an
expression for K in terms of eq by using Eq. (5.307) and thus obtain a value for
K using experimental data.
Examples of methods for determining concentrations/pressures:
measurement of gas density or pressure
light absorption
refractive index
electrical conductivity
Notes:
1. The measurement may not change any of the concentrations/pressures in the
system - otherwise the measurement would change the equilibrium as well!
2. The same value for K should be obtained when the equilibrium is
approached from either side.
3. The same value for K should be obtained over a wide range of initial
concentrations.
221
The relation between K and eq :
We have already seen that the equilibrium constant K can be written in terms
of eq (e.g., the water dissociation example). In the following we will consider a
general equation of the form: A mB and derive a general expression for the
equilibrium constant. The reaction is assumed to occur in the gas phase and all
gaseous components should follow the ideal gas law. According to Eq. (5.307) the
equilibrium constant can be written as:
PB m
PA 1
(5.308)
K=
P
P
Let us denote the total pressure by P and partial pressures of A and B by PA =
yA P and PB = yB P . The concentrations of A and B at equilibrium can be written
in terms of eq (note the scaling of ):
nA,eq = nA,0 (1 eq ) and nB,eq = nA,0 eq m
(5.309)
222
PA = P
1 eq
eq m
and PB = P
1 eq (1 m)
1 eq (1 m)
(5.311)
(5.312)
1 eq (1 m)
1 eq (1 m)
m1
P/P
(eq m)m
=
1 eq (1 m)
(1 eq )
K=
P
P
m1
2 P/P
4eq
2
1 eq
(5.313)
223
The relation between eq and gas densities:
When m > 1, it is possible to determine the extent of reaction by measuring the
density of the gas. For example, consider reaction A = 2B and denote the mass
of A by mA and the molecular mass by MA . Also assume that temperature and
pressure are constant and that both A and B follow the ideal gas law. The initial
volume before the reaction takes place is:
mA RT
Vini =
(where we used n = mA /MA )
(5.315)
MA P
At equilibrium we have a similar equation (the total mass is conserved in chemical
reactions and we still use mA below):
1
mA RT
where Meq = yA MA + yB MB = yA MA + yB MA
(5.316)
Veq =
Meq P
2
where we also used the fact that MA = 2MB (the mass conservation restriction).
Calculation of the ratios between the volumes gives:
yA MA + 12 yB MA
Vini
Meq
1
=
=
=
Veq
MA
MA
1 + eq
2eq
1 eq
and yB =
with yA =
1 + eq
1 + eq
(5.317)
224
Note that the mass conservation also relates the volumes to densities :
eq
Vini
ini Vini = mA
=
eq Veq = mA
ini
Veq
If we solve for eq on previous page, we get two useful results:
MA Meq
ini eq
eq =
and eq =
Meq
eq
(5.318)
(5.319)
Thus by measuring either the molar masses or gas densities before the reaction takes
place and at the equilibrium, Eq. (5.319) can be used for calculating eq . Note that
this derivation assumed a specific stoichiometry for the chemical equation.
Example. Consider reaction N2 O4 (g) = 2NO2 (g) (molecular weight of N2 O4 is
92.01 g mol1 ) at constant T (298.15 K) and P (1.0133 bar). 1.588 g of N2 O4
dissociates and expands to 500 cm3 volume. What is the extent of reaction eq and
the equilibrium constant K?
Solution. eq can be calculated using Eq. (5.319) with the molecular weights as
follows:
(1.588 g) 0.083145 L bar K1 mol1 (298.15 K)
mN2 O4 RT
=
Meq =
Veq P
(1.013 bar) (0.5 L)
= 77.70 g mol1
225
Mini was given as 92.01 g mol1 and now by applying Eq. (5.319) we get:
92.01 g mol1 77.70 g mol1
Mini Meq
= 0.1842
=
eq =
Meq
77.70 g mol1
The equilibrium constant can be obtained from Eq. (5.314):
K=
2 (P/P )
4eq
2
1 eq
4 (0.1842)2 (1.0133)
1 (0.1842)2
= 0.143
For a more complicated reaction A + 3B = 2C, the following expression for the
equilibrium constant can be obtained:
K=
2 (1 )
16eq
eq
(1 3eq )3 (P/P )2
(5.320)
Each type of reaction requires its own expressions for K and eq but they can be
derived by using the previously outlined approach.
226
Example. Consider the following reaction: N2 (g) + 3H2 (g) = 2NH3 (g). What
total pressure must be used to obtain 10% conversion of nitrogen to ammonia at
400 C? Assume equimolar initial mixture of N2 and H2 and ideal gas behavior.
The equilibrium constant for the reaction is 1.60 104 at standard pressure of 1
bar.
Solution. 10% conversion corresponds to eq = 0.1. We use Eq. (5.320):
P = P
2 (1 )
16eq
eq
(1 3eq )3 K
!1/2
= (1 bar)
16 (0.1)2 (1 0.1)
!1/2
= 51.2 bar
Note: Calculation of the equilibrium compositions of a reaction system that contains
two or more reactions is more complicated (results in many different equations
instead of just one). Furthermore the equations are non-linear, which typically
means that it is not possible to find analytic solutions. Numerical methods are the
only way to proceed in such case (i.e., computers).
227
Ns
X
X
(5.321)
vj Sj
vi f Gi with f Gi = f Hi T Si +
r G =
i=1
j6=i
The entropy summation is over the elements that are required to form compound i.
f Gi is the Gibbs energy for formation of species i. Note that at all temperatures
f G (H+ (aq)) = f H (H+ (aq)) = S (H+ (aq)) = 0.
First use the NIST Chemistry WebBook (http://webbook.nist.gov/chemistry/)
to get the enthalpies of formation (r Hi ) and entropies (Si ) and finally calculate
r Gi by using Eq. (5.321).
228
Note: Even when r H is positive (i.e., endothermic reaction), r S may be sufficiently large so that it overall produces negative Gibbs energy. Such reactions are
entropy driven.
Example. Calculate the standard Gibbs energy of formation for gaseous water
H2 (g) + 12 O2 (g) = H2 O(g) given the following calorimetric data (T = 298.15 K):
H2 O(g)
H2 (g)
O2 (g)
f H (kJ mol1 )
241.8
0
0
S (J K1 mol1 )
188.8
130.7
205.1
1
|
{z
}
0.1888 kJ K1 mol
vj Sj
1
0.1307 kJ K1 mol1 0.5 0.2051 kJ K1 mol1
= 228.6 kJ mol1
Note that the quantities with f above are already per mole quantities and therefore
do not require use of overbars.
229
Example. Calculate the equilibrium constants for the following reactions at the
indicated temperature:
a) 3O2 (g) = 2O3 (g) at 25 C
b) CO(g) + 2H2 (g) = CH3 OH(g) at 500 K (f G (CH3 OH) = 134.3 kJ mol1 )
graphite
230
SO
= 220.8 J K1 mol1 , f S (CO(g)) =
2 (g)
J
1
212.8 5.6 220.8
2
K mol
= 134.3 kJ mol
r G
K = exp
= exp
RT
153.0 kJ mol1 = 18.7 kJ mol1
!
18700 J mol1
= 1.11 102
8.31 J K1 mol1 (500 K)
231
(ln (K))
(r G /T )
RT 2
=
T 2
=
|{z}
|{z}
T
T
P
P
Eq.(5.302)
Eq.(4.242)
(ln (K))
r H
=
(5.322)
(5.323)
232
is zero (Eqs. (2.112) and
Note: If r H is independent of temperature, then r CP
is zero then S
(2.140)). It will turn out later in Eq. (5.330) that since r CP
r
is also independent of temperature.
= 0 and combining Eqs. (5.302) and (5.304) we get:
By assuming r CP
r G = RT ln (K) and r G = r H T r S
r H
r S
ln (K) =
+
RT
R
(5.325)
1900
2.31
2000
4.08
2100
6.86
2200
11.0
2300
16.9
2400
25.1
2500
36.0
2600
50.3
Solution. Use Eq. (5.325) to obtain the slope and the intercept from a ln (K) vs.
1/T plot. This must be done by fitting Eq. (5.325) to the experimental data. From
the slope one can obtain r H by multiplying by R and r S by multiplying
the intercept by R. The plot (data and fitting) is shown below.
233
-5
ln(K)
-6
-7
-8
-9
0.0004
0.00045
0.0005
-1
1 / T (K )
Note that the assumption of no temperature dependency in r H and r S appears to be a good one here. If this was not the case, the above plot would not yield
a straight line (because the slope and intercept would depend on T ).
The standard reaction entropy r S indicates how much the entropy changes in a
reaction under standard conditions and at a given temperature. For a given reaction
it can be calculated by:
Ns
X
vi Si
(5.326)
r S =
i=1
234
Example. Calculate the standard reaction entropies for the following reactions at
298 K by using the data in the Chemistry WebBook:
(a) H2 (g) + 12 O2 (g) = H2 O(g)
(b) N2 (g) + 3H2 (g) = 2NH3 (g)
Solution. Eq. (5.326) gives the following results:
1
(a) r S = SH
SH
SO2 (g) = 69.95 J K1 mol1
2 O(g)
2 (g)
2
1
205.15 J K1 mol1 = 163.33 J K1 mol1
130.68 J K1 mol1
2
(b) r S = 2SNH
SN
3SH
= 2 192.77 J K1 mol1
3 (g)
2 (g)
2 (g)
191.61 J K1 mol1 3 130.68 J K1 mol1 = 198.11 J K1 mol1
235
f Hi (T )
f Hi (298.15
ZT
K) +
(T )dT
C
P,i
(5.327)
298.15 K
ZT
Si (T ) = Si (298.15 K) +
298.15 K
r H (T ) = r H (298.15 K) +
P,i (T )
C
dT
T
ZT
(5.328)
r C P
(T )dT
(5.329)
(T )
r C P
dT
T
(5.330)
298.15 K
r S (T ) = r S (298.15 K) +
ZT
298.15 K
CP,i
N
Ps
i=1
vi f Hi and r S =
N
Ps
i=1
vi f Si . Recall
ZT
298.15 K
r C P
(T )dT T
ZT
298.15 K
(T )
r C P
dT
T
(5.331)
236
Note that above r G (298.15 K) must be interpreted as (T 6= 298.15 K!):
r G (298.15 K) = r H (298.15 K) T r S (298.15 K)
Recall from Eq. (5.302) that ln (K) = r G /RT and thus we get:
ln (K(T )) =
1
R
ZT
298.15 K
298.15 K
1
ln (K(298.15 K))
T
RT
ZT
r C P
(T )dT (5.332)
298.15 K
(T )
r C P
dT
T
237
K=
Ns
Y
Pi vi
i=1
P
P
v
Ns
Y
yi P vi
i=1
Ky where v =
Ns
X
i=1
Ns
Ns
Y
P vi Y
v
yi i
P
i=1
i=1
vi and Ky =
Ns
Y
(5.333)
yi i
i=1
Note that Ky depends on the molar fractions of the gas components and changes in
it imply changes in the molar fractions of the different components. Now consider
changes in the total pressure P at constant temperature and rewrite the previous
equation as:
P v
Ky =
K
(5.334)
238
Consider the following three cases:
1. v < 0: The total amount of moles of gas decreases in the reaction. According
to Eq. (5.334) Ky should increase as pressure increases. This means that
the products are favored over the reactants (i.e. balance shifts forward).
2. v > 0: The total amount of moles of gas increases in the reaction. According
to Eq. (5.334) Ky should decrease as pressure increases. This means that
the reactants are favored over the products (i.e. balance shifts backward).
3. v = 0: The total amount of moles of gas remains independent of pressure.
Thus changes in pressure do not change the molar fractions of the gases.
Notes:
1. This result reflects the Le Chateliers principle.
2. If the reactants and products are solids or liquids, the effect of pressure on
the equilibrium is small.
The effect of initial composition on equilibrium
First we write the molar fraction yi in terms of the extent of reaction . The amount
of species i and the total amount of gas for a given are given by:
ni () = ni,0 + vi and ntot () = ntot,0 + v
239
From these quantities yi can be calculated: yi () =
ni ()
ntot ()
ni,0 +v
ntot,0 +v
Inserting this into definition of Ky we get the following relation between Ky and :
v Y
Ns
Ns
Y
ni,0 + vi vi
1
(ni,0 + vi )vi
=
ntot,0 + v
ntot,0 + v
i=1
i=1
i=1
(5.335)
By combining Eqs. (5.333) and (5.335), the equilibrium constant can be written as:
v Y
Ns
P/P
K=
(ni,0 + vi )vi
(5.336)
ntot,0 + v
i=1
Ky =
Ns
Y
yi i =
1. Consider addition of inert gas (i.e. one that does not take part in chemical
reaction) with temperature and volume constant. Remember that since the amount
of substance is not constant, we have four variables (P , V , T , and n) to consider.
The amount of inert gas added is denoted by ninert . The total amount of gas is
now given by ntot,0 + v + ninert and note that the inert component cancels out
from the product on the right in the equation above because it remains identical on
both sides of the chemical equation. The above equation now takes the following
form:
240
K=
P/P
ntot,0 + v + ninert
v Y
Ns
(ni,0 + vi )vi
(5.337)
i=1
K=
=
(ntot,0 + v + ninert )
RT
V P
RT
V
ntot,0 + v + ninert
v Y
Ns
/P
Ns
Y
(ni,0 + vi )vi
(5.338)
i=1
(ni,0 + vi )vi
i=1
Note that we had assumed that both T and V were constants, which means that
the right hand side above is independent of ninert . Therefore addition of inert gas
component to the reactive gas mixture does not change equilibrium if T and V are
kept constant during the addition (i.e. only the total pressure changes).
241
2. Next consider addition of an inert gas with temperature and pressure kept
constant (note that V changes in this process). If ninert changes, the prefactor
in Eq. (5.337) changes as well. Since K = K(T ) and temperature is constant,
the product term in Eq. (5.337) (i.e., Ky ) must counter balance the change. This
means that the mole fractions yi must take new values (i.e., the composition of the
gas mixture changes).
Summary:
v < 0: addition of inert gas shifts the balance in chemical equation to the left.
v = 0: addition of inert gas does not change the composition of the gas.
v > 0: addition of inert gas shifts the balance in chemical equation to the right.
242
Previously the equilibrium constants were written in terms of pressure. Now the
task is to express the equilibrium constant constant in terms of concentrations for
gas phase reactions.
KP = Equilibrium constant calculated using partial pressures.
Kc = Equilibrium constant calculated using concentrations.
Below we will show that usually KP 6= Kc and derive the relation between them.
First we note that partial pressures are related to concentrations ci by Pi,eq =
ni,eq RT /V = ci,eq RT and then Eq. (5.307) gives:
Ns
Y
ci RT vi
KP =
(5.339)
P
i=1
Next we introduce standard concentration denoted by c , which is 1 mol dm3 and
insert this into Eq. (5.339):
KP =
vi
Ns
Y
ci c RT
i=1
where v =
N
Ps
c RT
P
v Y
Ns
ci vi
i=1
c RT
P
v
Kc (5.340)
i=1
Kc =
Ns
Y
ci vi
i=1
(5.341)
243
Note that Kc is also a function of temperature only (just like KP ).
Example. What is the value of the equilibrium constant Kc for the dissociation
of ethane into methyl radicals at 1000 K? The reaction is: C2 H6 (g) = 2CH3 (g).
f G (CH3 ) = 159.8 kJ mol1 and f G (C2 H6 ) = 109.6 kJ mol1 .
Solution. First we calculate r G by Eq. (5.321):
r G = 2f G (CH3 ) f G (C2 H6 ) = 2 159.8 kJ mol1 109.6 kJ mol1
= 210.1 kJ mol1
1 mol L1
P
([CH3 ] /c )2
= KP
= 1.062 1011
[C2 H6 ] /c
c RT
1 bar
= 1.278 1013
0.08315 L bar K1 mol1 (1000 K)
Note that in above c is required to get the right unit cancellation for Kc .
244
A reaction involving more than one phase that does not involve equilibria of species
between phases is referred to as a heterogeneous reaction. For example the following
reactions are heterogeneous:
(1) CaCO3 (s) = CaO(s) + CO2 (g)
(2) CH4 (g) = C(s) + 2H2 (g)
Previously we have seen that pure gas phase reactions do not go to completion
essentially because of entropy of mixing. However, depending on the initial conditions, heterogeneous reactions can proceed to completion. This can be seen by
writing the equilibrium constant for reactions (1) and (2) above:
PCO2
(5.342)
KP =
P
KP =
2
PH
2
PCH4 P
(5.343)
The reason why solids do not need to be entered above, is that their activities are
very close to one up to moderate pressures. As such they do not contribute to the
equilibrium constant. Reaction (1) may go all the way to completion whereas reaction (2) may not because there are gaseous products on both sides of the chemical
equation. From (5.60) we can distinguish two different cases for reaction (1):
1. The reaction will go to completion trying to satisfy Eq. (5.342) (may not get
that far).
2. The reaction will not go to completion and stops when Eq. (5.342) is
satisfied.
245
(5.344)
The total number of variables that must be used in specifying the system.
The number of components present in the system (# of species - # of
reactions: Ns R).
The number of phases present in the system.
Justificaton:
Consider system with p phases. If a phase contains C components, its composition
may be specified by (C 1) mole fractions (the remaining fraction can obtained
from the fact that their sum must equal to 1). The total number of variables for
for all the phases is given by p (C 1). In addition to these, one must also
specify two more variables (like P & T , V & P or V & T ). Thus up to this point
we can conclude that we at least have F = p (C 1) + 2. Next we consider
phase equilibria between p different phases. For phases , , , we must have at
equilibrium ,i = ,i = ,i = ... for each component i. Thus there are (p 1)
such relationships for each component, which reduces the total number of required
variables by C(p 1).
Finally, by summing the previous contributions, we get the Gibbs phase rule:
F = [p (C 1) + 2] C (p 1) = C p + 2
246
Example. The reaction CaCO3 (s) = CaO(s) + CO2 (g) is at equilibrium.
(a) How many degrees of freedom are there when all three compounds are
present at equilibrium?
(b) How many degrees of freedom are there when only CaCO3 (s) and CO2 (g)
are present?
Solution.
(a) The number of components: C = Ns R = 3 1 = 2. Three compounds
and one chemical reaction.
The number of phases: p = 2 (solid and gas).
The total number of variables needed: F = C p + 2 = 2.
Thus both temperature and pressure may be varied independently without
destroying a phase.
(b) C = Ns R = 2 0 = 2 (no reaction; R = 0).
p = 2 (solid & gas).
F = 2 2 + 2 = 2. This is the same as in a) because the restriction due to
chemical reaction was removed but there is one component less. Both
temperature and pressure may be varied without destroying a phase.
248
249
Note that along the phase boundary lines two phases exist in equilibrium. The
Gibbs phase rule now gives F = 1, which means that only one variable may be
specified independently. The boundary lines can be understood as functions and
as such they introduce a dependency between P and T (i.e. reduces F to 1). At
points where all three different phases exist at the same time (triple point) we
have F = 0 (i.e. they are just points in the graphs and both P and T are fixed).
What defines the phase boundaries?
Recall that at constant T and P the the equilibrium (i.e. dG = 0) between two
phases and is Eq. (4.259):
dG = dn + dn = dn dn = 0 =
In order to tell which phase is stable at given pressure and temperature, we must
calculate the chemical potentials for each phase and compare them. The phase
with the lowest chemical potential is stable. To obtain general statements about
stabilities of different phases, we must recall the following results (Eqs. (4.260) and
(4.261)):
(6.345)
S =
T P
(6.346)
V =
P T
Remember that derivatives define slopes (here for the chemical potential ).
250
1. Dependency of on temperature
We will use Eq (6.345) above to predict the slopes for = (T ). Remember that S is
always non-negative, which means that (T ) must have a non-negative slope. From
S(gas).
Based on these results we can graph (T ) for each phase qualitatively.
= melting temperature
= boiling temperature
= constant
251
2. Dependency of on pressure
The previous Eq. (6.346) gives the slope for (P ) when temperature is constant.
The slope is equal to the molar volume, which is a non-negative quantity. The
molar volumes for different phases are known to follow: Vsolid < Vliquid << Vgas .
An exception to this rule is, for example, water, which we will consider separately.
= melting pressure
= boiling pressure
= constant
If gas is compressed, it will first become liquid and upon more compression it will
eventually become a solid. At some given temperature a liquid can be made to boil
if the pressure is reduced accordingly. For example, vacuum distillation is based on
this principle.
252
For water we have: Vliquid < Vsolid << Vgas
Here it turns out that the volume of the solid is larger than the liquid (the density
of liquid is higher than the solid). For this reason, for example, ice cubes stay on
top of water or water upon freezing expands and in nature this expansion process
can break even rocks. Most solids sink in the corresponding liquid because they
contract on freezing (i.e., the density increases).
253
3. Combined effect of pressure and temperature
The combined effect of pressure and temperature is shown below by plotting (T )
vs. T at two different pressures (high and low):
254
The previously discussed changes in phase are said to be first order phase transitions.
The first order refers to the fact the the chemical potential has discontinuous first
derivative with respect to pressure and temperature. For example d(P )/dP is
discontinuous at phase transitions:
Thus above on each side of the non-differentiable points, the substance has different
molar volumes. At transition points heat capacities also tend to approach infinity.
This can be understood qualitatively by inspecting Eqs. (2.100) and (2.112):
CV =
dqP
qV
qP
dqV
and CP =
T =
and T =
dT
dT
CV (T )
CP (T )
255
At transition points all the supplied heat (q) is used for driving the phase transition
(i.e. the temeprature remains constant until the phase transition is complete).
According to the previous equations this can only happen if the heat capacities are
infinite at the transition points (T = 0, q > 0 C ).
Example. A special case of phase transition occurs between two liquid phases He-I
and He-II (normal and superfluid 4 He) at low temperatures. Due to the shape of
the transition, it is also called -transition (under saturated vapor pressure this
occurs at 2.17 K).
256
In a second order phase transition (no heat involved in transition) the first derivative
of chemical potential is continuous but the second derivative has a discountinuity.
The heat capacity does not become infinite at these points but show only a discontinuous jump. Behavior of various thermodynamic quantities at first and second
order transition points is shown below.
257
(6.347)
Benoit Clapeyron,
French engineer
and physicist
(1799 - 1864).
258
The Gibbs-Duhem equation (Eq. (4.272)) for both phases gives:
n d V dP + S dT = 0
n d V dP + S dT = 0
Solving for ds gives:
d = V dP S dT
d = V dP S dT
Since d = d , we get essentially the Clapeyron equation:
V dP S dT = V dP S dT
This can be rewritten (using Eq. (3.176); S = H/T ) as:
S S
S
H
dP
= =
=
V V
dT
V
T V
(6.348)
(6.349)
The deltas in this equation refer to differences in the values for the two phases.
This Clapeyron equation may be applied to vaporization, sublimation, fusion or the
transition between two solid phases of a pure substance.
259
Example. What is the rate of change per Pascal (Pa) in the boiling point of water
at a 100 C in atmospheric pressure? The molar enthalpy of vaporization is 40.69
kJ mol1 , the molar volume of liquid water is 0.019 103 m3 mol1 , and the
molar volume of steam is 30.199 103 m3 mol1 . All values are given at 100 C
and 1.01325 bar.
Solution. Use the Clapeyron Eq. (6.349):
vap H
dP
40690 J mol1
=
= 3613 Pa K1
=
dT
T Vg Vl
(373.15 K) 30.180 103 m3 mol1
dP
dT
1
1
= 2.76 104 K Pa1
3613 Pa K1
Example. Calculate the change in pressure required to increase the freezing point
of water by 1 mK. At 273.15 K the heat of fusion of ice is 333.5 J g1 , the density
of water is 0.9998 g cm3 (= g / mL = kg / L), and the density of ice is 0.9168 g
cm3 .
Solution. First we note that the molar volumes are given by the inverse of density:
1
V =
260
Therefore we get the molar volumes for the liquid and solid:
Vl =
1
1
= 1.0002 cm3 g1 and Vs =
= 1.0908 cm3 g1
0.9998 g cm3
0.9168 g cm3
Note that we have expressed everything in terms of grams rather than moles. If
molar volume would be needed in units of volume / mol then one should use the
molecular weight of the substance to convert. Provided that the molar quantities
of H and V are independent of temperature and pressure, we can integrate Eq.
(6.349):
Tf
Z dT
f us H
Tf
f us H
f us H
dP =
dT
P
=
=
ln
T V
V
T
V
Ti
Ti
Tf
f us H
f us H
V
Ti
V Ti
For a change of 1 mK both f us H and the molar volumes are approximately
constant. In the present case Ti = 273.150 K (freezing point of water at ambient
pressure), T = 0.001 K and Tf = 273.151 K. The task is now to find P :
P =
333.5 J g1 (0.001 K)
Note that for water dP/dT < 0! This means that the phase boundary slope in the
phase diagram is negative. In most cases V > 0 and dP/dT > 0.
261
Example. Calculate the vapor pressure of H2 O(l) at 298.15 K when the following
values are given: f G298 K (H2 O(g)) = 228.6 kJ mol1 and f G298 K (H2 O(l))
= 237.1 kJ mol1 . Assume that the gases follow the ideal gas law.
Solution. Phase changes can be considered as chemical reactions:
H2 O(l) = H2 O(g)
Recall Eq. (5.302) and the definition of equilibrium constant Eq. (5.301):
a (H O(g))
2
r G = RT ln (K) = RT ln
= RT ln (a (H2 O(g)))
a (H2 O(l))
|
{z
}
=1
r G = f G (H2 O(g))f G (H2 O(l)) = 228.6 kJ mol1 237.1 kJ mol1
= 8.56 kJ mol1
The ideal gas assumption with Eq. (5.301) and replacing the subsrcipt r with vap
(our reaction is vaporization):
vap G
vap G = RT ln PH2 O(g) /P PH2 O(g) = P exp
RT
262
8560 J mol1
8.3145 J K1 mol1 (298.15 K)
= 0.032 bar
where r G = f G
(dia)
f G (gra) = (2900 0) J/mol. Eq. (4.222) can now
(r G)
be applied: V =
. This can be integrated:
P
T,n
Zr G
r G
d(r G) =
ZP
V dP
263
where r G = 0 at equilibrium. Solving for P gives:
P =
r G r G
2900 kJ/mol
=
= 1.52 109 Pa
V
1.91 106 m3 /mol
V = 12 g mol1
1
1
3.51 g m3
2.25 g m3
The required pressure is P = P 0 + P 1.52 109 Pa. Since the temperature was
fixed at 25 C, this much pressure must be applied in order to reach the two-phase
equilibrium (i.e., must be located on the phase boundary line).
264
For vaporization and sublimation processes Rudloph Clausius showed that the
Clapeyron equation may be simplified by assuming that the vapor phase obeys
the ideal gas law and the molar volume of the liquid (Vl ) is small in comparison
with the molar volume of the gas (Vg ). Substituting RT /P for Vg , Eq. (6.349)
becomes:
P vap H
vap H
dP
=
(6.350)
dT
T Vg
RT 2
By separating variables to the opposite sides of the equation, Eq. (6.350) becomes:
dP
vap H
=
dT
P
RT 2
ZP
dP
=
P
ZT2
T1
vap H
dT
RT 2
(6.351)
P
P
vap H
R
ZT2
T1
T 2 dT =
vap H
R
1
1
T1
T2
(6.352)
(6.353)
265
ln(P)
This shows that a plot of ln (P ) vs. 1/T should be a straight line. When P and
T are obtained from experiments, it is possible to extract vap H from the (P, T )
data. Such data is shown for water as an example below.
A linear fit gives vap H = 44.7 kJ mol1 (literature value 44.0 kJ mol1 at 298
K). Note that the above experiment gives an average value for this quantity over the
measurement temperature range as there is a small dependency on temperature.
266
The Clausius-Clayperon equation can also be written as:
1
P2
vap H
1
vap H (T2 T1 )
=
ln
=
P1
RT1 T2
R
T1
T2
(6.354)
The above derivations suffer from two approximations: 1) vap H was assumed to
be independent of temperature and 2) the gas phase was assumed to follow the ideal
gas law. The temperature dependency can be included in vap H approximately
by assuming the following form (not a very good approximation either):
vap H = A + BT + CT 2
(6.355)
(6.356)
(6.357)
267
To see how vap H can be determined graphically, we integrate Eq. (6.351):
dP
vap H
=
dT
P
RT 2
ZP
dP
=
P
ZT
0
vap H
dT
RT 2
= 2
dT
dT
d (1/T )
T
d (1/T )
and therefore we have:
vap H = R
d (ln (P/P ))
d (1/T )
(6.358)
Notes:
P above is the vapor pressure (the total external pressure is constant).
The above form is particularly useful for onbtaining a graphical solution
because the plot of R d ((ln (P/P )) /d (1/T ) gives vap H (which may
or may not be constant).
268
(6.359)
For ideal gases the chemical potential for each component in the gas phase is given
by Eq. (5.300) with ai = Pi /P :
Pi
i (g) = i (g) + RT ln
where i = 1, 2
(6.360)
P
For the liquid phase components we have to use the general form (Eq. (4.255)):
i (l) = i (l) + RT ln (ai )
(6.361)
269
By combining the above equations:
Pi
i (g) + RT ln
= i (l) + RT ln (ai )
(6.362)
(6.363)
where Pi is the equilibrium vapor pressure of pure i. Subtract Eq. (6.363) from
(6.362) and we can rearrange the resulting equation as:
Pi
Pi
= RT ln(ai ) ai =
(6.364)
RT ln
Pi
Pi
Thus, if the vapor follows the ideal gas law, the activity of a component in a solution
is given by the ratio of its partial pressure above the solution to the vapor pressure
of the pure liquid.
An experimental finding (Francois-Marie Raoult, 1830 - 1901):
Pi xi Pi (Raoults law)
(6.365)
where Pi is the equilibrium pressure of pure liquid i and xi is its molar fraction in
liquid phase.
270
This holds well when the components are chemically similar. An eample of chemically similar compounds is given by a benzene/toluene mixture. The Raoults law
is demonstrated for this mixture below.
Because the gas phase was assumed to be ideal, we can express the partial pressure
of a gas component i by multiplying the total gas pressure by the molar fraction in
the gas phase (yi ):
Pi = yi P xi Pi
(6.366)
This lets us connect the molar fraction in the liquid phase (xi ) to the gas phase
quantities.
271
From Eq. (6.366) we can solve for xi (the liquid phase molar fraction):
P
Pi
x i = yi =
Pi
Pi
(6.367)
(6.368)
Now we have managed to get an expression for activity in the liquid phase by
considering the gas phase and the Raoults law. Inserting this into Eq. (6.21) gives:
i (l) = i (l) + RT ln (xi )
|{z}
(6.369)
P = P1 + P2 = x1 P1 + x2 P2 = P2 + (P1 P2 ) x1
(6.370)
=ai
Solutions for which Eq. (6.368) holds are called ideal solutions. Recall that yi is
the mole fraction in the gas phase and xi is the mole fraction in the liquid phase.
The total vapor pressure of an ideal binary mixture is given by:
Why is the total pressure line also called the bubble point line?
A liquid boils when its vapor pressure exceeds the external pressure. For example,
under normal conditions water boils at 100 C because at that point its vapor
pressure is equal to the atmospheric pressure (1 atm). Boiling is seen as formation
of gas bubbles in the liquid. Consider the plot below where the total vapor pressure
is shown as a function of x1 .
272
The vapor composition of a binary solution can be calculated using Raoults law:
x1 P1
x1 P1
P1
=
=
(6.371)
y1 =
P1 + P2
x1 P 1 + x2 P 2
P2 + P1 P2 x1
This equation can be solved for x1 :
x1 =
y1 P2
P1 + P2 P1 y1
(6.372)
(6.373)
273
This function is plotted below (i.e., total pressure vs. molar fraction in the gas
phase) for benzene/toluene solution:
At external pressures below the dew point line
the binary mixture exists as a two-component
gas. Above this line liquid starts to form.
P1 = Vapor pressure of pure toluene.
P2 = Vapor pressure of pure benzene.
Liquid and vapor compositions for benzene/toluene mixture is shown. Note that the
x-axis might be somewhat confusing: it is used
for two different variables x1 and y1 . This does
not mean that they would be equal! Variable
x1 corresponds to the top straight line and y1
for the lower curve. A tie line (dashed line)
connects x1 and y1 at a given total pressure
P . In the pressure ranges from P1 to P2 , the
system may exist in liquid-gas phase equilibrium with a given x1 and y1 .
274
C
Example. At 60
the vapor pressures of pure benzene and toluene are 0.513
and 0.185 bar, respectively. What are the equations of the bubble point and dew
point lines? For a solution with 0.60 mole fraction of toluene, what are the partial
pressures of toluene and benzene, and what is the mole fraction of toluene in the
vapor? Assume ideal solutions.
Solution. Denote toluene by 1 and benzene by 2. The bubble point line is given
by Eq. (6.370):
P (x1 ) = P2 + (P1 P2 ) x1 = 0.513 bar (0.328 bar) x1
The dew point line is given by Eq. (6.373):
P (y1 ) =
P1
P1 P2
0.0949 bar2
=
With 0.60 mole fraction of toluene, the partial pressures are given by Eq. (6.365):
P1 = x1 P1 = 0.60 (0.185 bar) = 0.111 bar
P2 = x2 P2 = 0.40 (0.513 bar) = 0.205 bar
From the above bubble point line equation with x1 = 0.60, we get:
P = 0.513 bar (0.328 bar) (0.60) = 0.316 bar
275
Eq. (6.371) relates y1 and x1 to each other and hence:
y1 =
x1 P1
0.60 (0.185 bar)
= 0.351
=
0.513 bar (0.328 bar) 0.60
P2 + P1 P2 x1
Example. Calculate the activities of toluene (component 1) and benzene (component 2) in the liquid by using the values given in the previous example. Assume
ideal solutions.
Solution. Use Eq. (6.364):
a1 =
P1
P2
= 0.600 and a2 = = 0.400
P1
P2
Note that we assumed an ideal solution, the activities are equal to the mole fractions.
Thermodynamic quantities of ideal solutions
In ideal solutions the activities are given by Eq. (6.369): i (l) = i (l) + RT ln (xi ).
The molar Gibbs energy of a solution at constant T and P is given by Eq. (4.225):
=
G
Ns
X
i=1
s
X
ni
i x i
=
n
i=1
276
Inserting the chemical potential expression in above, we get (constant T and P ):
=
G
Ns
X
x i i =
Ns
X
xi i + RT
xi ln (xi )
(6.374)
i=1
i=1
i=1
Ns
X
Note that this expression gives the molar Gibbs energy as a function of pressure and
temperature. Thus this expression contains all of the thermodynamic information
about an ideal solution.
The molar entropy of an ideal solution can be obtained using Eqs. (4.221) and
(4.227):
S =
G
T
Ns
X
i=1
xi
i
T
Ns
X
xi ln (xi ) =
Ns
X
i=1
i=1
xi Si R
Ns
X
xi ln (xi )
i=1
(6.375)
The molar enthalpy of an ideal solution can be calculated by using the Gibbs = , and G
= H
Helmholtz equation (using Eqs. (4.241), (4.227) with G
i
i
i
i
T Si ):
= T 2
H
"
G/T
T
Ns
X
i=1
xi
#
P,{ni }
= T
i
+ i
T
Ns
X
i=1
"
Ns
X
i=1
xi
i /T
T
#
= T 2
Ns
X
xi
i=1
Ns
X
=
xi H
xi T Si + G
i
i
i=1
1 i
i2
T T
T
(6.376)
277
The molar volume is obtained from Eq. (4.222):
=i
i (6.369)
Ns
}|
{
z
X
X
X
G
xi Vi
=
xi (i /P ) =
V =
/P
xi RT ln(xi )
P T,{ni }
i=1
(6.377)
Remember above that i does not depend on pressure - only on temperature.
The first terms (with superscripts ) in Eqs. (6.374) and (6.375) give the thermodynamic properties of isolated species (i.e. no mixing) and the possible second term
the effect of mixing:
mix G = RT
Ns
X
xi ln (xi )
(6.378)
xi ln (xi )
(6.379)
i=1
mix S = R
Ns
X
i=1
mix H = 0
(6.380)
mix V = 0
(6.381)
These results are essentially the same as we obtained earlier for ideal gas mxitures.
Note that there is no volume change or heat evolution when ideal solutions are
formed under constant temperature and pressure.
278
The effect of temperature on ideal binary liquid mixtures at constant pressure
Consider toluene (1) - benzene (2) mixture at atmospheric pressure (units in bar):
P1
P2
80.1 C
1.013
88 C
0.508
1.285
90 C
0.543
1.361
94 C
0.616
1.526
98 C
0.698
1.705
100 C
0.742
1.800
104 C
0.836
2.004
110.6 C
1.013
The boiling points for pure liquids were underlined above. Vapor pressures for pure
liquids can be determined experimentally.
Example. Calculate the liquid mole fraction (x1 ) of toluene in benzene-toluene
solution that boils at 100 C. Calculate also the toluene mole fraction (y1 ) in the
vapor above the liquid.
Solution. We use Eq. (6.370) and solve for x1 , recall that a liquid boils when
its vapor pressure reaches the external pressure (one atmosphere here) and use the
above table:
x1 =
P P2
(1.013 bar) (1.800 bar)
=
= 0.744
P1 P2
(0.742 bar) (1.800 bar)
x1 P1
0.744 (0.742 bar)
=
= 0.545
P
1.013 bar
279
Example. The method of fractional distillation can be understood in the previously
developed theory. Consider a binary ideal mixture with components A and B. In
this example component A (i.e., the volatile component) has a lower boiling point
than B.
1. Consider a solution with mole fraction
xA,1 of A.
2. Heat the solution to its boiling point
T2 . The liquid mole fraction is now
xA,2 .
3. Extract vapor, which has mole fraction
yA,2 . Note that there is more of the
volatile component in the gas phase
than in the liquid.
Repeat the above cycle for the liquid, which
was condensed from the gas phase above. The
next step is shown as xA,3 and yA,3 . By repeating the cycle, a better separation of the
components can be reached.
280
In the previous section we assumed that Raoults law holds (Eq. (6.365)). However,
in many cases significant deviations from this law are found:
The points where the bubble point and dew point lines cross, are called azeotropic
points. In this case distillation can not separate the components!
281
What is the origin of the deviation from the Raoults law?
This extra binding interaction lowers the vapor pressure from that predicted by the
Raoults law. The maximum effect is seen when there are equal amounts of the
components (i.e. the middle of the bubble point line on previous slide).
A positive trend (i.e. the bubble point line is above the Raoults law prediction)
is observed when A A and B B interactions are stronger than A B. Here A
and B are the components forming the binary mixture. If these deviations are large
enough, the system separates into two phases A and B.
Henrys law attempts to fix the deficiency of the Raoults at small mole fractions,
which are close to zero or one (i.e. nearly one component solutions). It consists of
the following linear approximation:
Pi = Ki xi
(6.382)
where Ki is Henrys law constant for component i. Its value can be obtained from
plotting the ratio Pi /xi when xi approaches zero.
282
Behavior of the partial pressures as a function of xi (the liquid phase mole fraction)
is shown below for Henrys law, Raoults law and the experimental data (continuous
line):
Henrys law constants given in Pascals (25
C).
Gas
H2
N2
O2
CO
CO2
CH4
C2 H 2
C2 H 4
C2 H 6
Water
7.12 109
8.68 109
4.40 109
5.79 109
0.167 109
4.19 109
0.135 109
1.16 109
3.07 109
Benzene
0.367 109
0.239 109
0.163 109
0.0114 109
0.0569 109
283
Solution. We apply Henrys law (Eq. (6.382)). Denote water by component 1 and
CO2 by component 2. The partial pressure of CO2 in the gas phase was given as:
P2 = 105 Pa. Now Eq. (6.382) gives:
x2 =
P2
105 Pa
=
= 6.0 104
K2
0.167 109 Pa
[CO2 ] + [H2 O]
[H2 O]
[CO2 ] = x2 [H2 O] = 6.0 104 55.5 mol L1 = 3.3 102 mol L1
x2 =
Notes:
The solubility of a gas in liquids usually decreases with increasing
temperature because heat is generally evolved in the solvation process.
There are exceptions to this - especially with solvents like liquid ammonia,
molten silver and some organic liquids.
Solubility of an unreactive gas in a liquid is due to intermolecular attractive
forces (van der Waals forces) between gas molecules and solvent molecules.
Addition of electrolytes (i.e. ionic species) usually decreases solubility of
gases in liquids (salting out).
284
If the vapor pressure of a solute follows Henrys law, then we can insert Eq. (6.382)
into Eq. (6.360):
Ki xi
i (g) = i (g) + RT ln
(6.383)
P
At equilibrium i (g) = i (l) (see Eq. (6.359)) and expanding the logarithm gives:
Ki
i (l) = i (g) + RT ln
+RT ln (xi ) = i (l) + RT ln (xi )
(6.384)
P
|
{z
}
=
i (l)
where we introduced the standard chemical potential for species i in the liquid:
Ki
(6.385)
i (l) = i (g) + RT ln
P
285
(6.386)
where i is the activity coefficient for species i. With this notation, Eq. (6.361)
reads:
i (l) = i (l) + RT ln (i xi )
(6.387)
286
Activity coefficients can also be expressed with reference to Henrys law. To see
this, we first express the activity as a product of the activity coefficient (denoted by
prime to differentiate it from the activity coefficient with Raoults law reference):
i (l) = i (l) + RT ln i xi
(6.391)
By combining Eqs. (6.359) and (6.360) we get:
i (l) = i (g) + RT ln
Next recall Eq. (6.385): i (l) = i (g) + RT ln
Pi
P
(6.392)
Ki
P
By substituting the above two expressions into Eq. (6.391), the following modified
form for Henrys law can be obtained:
Pi = i Ki xi
(6.393)
i :
(6.394)
(6.395)
287
What is the principal difference between i and i ?
i (Raoults law reference) has a reference state, which consists of an ideal solution.
The chemical potential corresponding to this standard state is denoted by i .
i (Henrys law reference) has a reference state, where each molecule of the solute
has the same interactions that it experiences in very dilute solutions. The chemical
potential corresponding to this standard state is denoted by i .
Thus the only difference is that they have different references states.
Example. Calculate the activity coefficients i for ether (component 1) and acetone
(component 2) in 1:1 ether-acetone solutions at 30 C. The experimental data are
given in the table below.
Table. Activity coefficients for acetone-ether solutions at 30 C. 1 = ether, 2 =
acetone. Pressure is given in kPa.
x2
0
0.2
0.4
0.5
0.6
0.8
1.0
P1
86.1
71.3
58.7
52.1
44.3
26.9
0
37.7
37.7
1.60
1.31
1.19
1.12
1.04
1.00
Henrys
K2 x2
0
15.7
31.4
39.2
47.0
62.7
78.4
law reference
2
1.00
0.77
0.63
0.57
0.54
0.50
0.48
288
Solution. At 0.5 mole fraction, the activity coefficients of the two components are
given by:
1 =
P1
52.1 kPa
=
= 1.21
x1 P1
43.1 kPa
2 =
P2
22.4 kPa
= 1.19
=
x2 P2
18.9 kPa
Note that the Henrys law activity coefficient i can be calculated from the data
given in the previous table by extrapolating to x2 = 0 when calculating K2 = P2 /x2
(see Eq. (6.382)).
What is the relation between i and i ?
According to Eq. (6.364) the activity of a real solution is given by ai = Pi /Pi .
Substitution of Eq. (6.394) there gives:
Ki xi
ai = i
(6.396)
Pi
289
With ai = i xi this gives:
i =
i Pi
i Ki
or i =
Pi
Ki
(6.397)
Notes:
For dilute solutions, the solvent is usually treated on the basis of deviations
from Raoults law and the solute is usually treated on the basis of deviations
from Henrys law.
Three different concentration scales can be used with Henrys law:
Pi = Ki mi where mi is the molal concentration of i (mol / kg of solvent),
Pi = Ki ci where ci is the molar concentration of i (mol / L of solution), or
Pi = Ki xi where xi is the (liquid) molar fraction.
Since we mostly calculate differences in a given thermodynamic variable, the
reference (= standard) state usualy cancels out. However, both initial and
final states must be expressed with respect to the same reference state.
290
Colligative properties: depression of freezing point, elevation of boiling point, osmotic pressure, and the lowering of the vapor pressure by a non-volatile solute.
The depression of freezing point
Consider an ideal mixture of solvent A and solute B (i.e. A is in excess of B). The
phase equilibrium between solid A and liquid A is given by Eqs. (6.359) and (6.361)
(also note that for solids = and for ideal solutions ai = xi ):
A (s, T ) = A (s, T ) = A (l, T ) = A + RT ln (aA ) = A + RT ln (xA )
(6.398)
A (s, T ) A (l, T )
f us GA (l, T )
=
RT
RT
(6.399)
291
Next we assume that f us H and f us S are independent of temperature near
the freezing point of A (Eq. (3.176) and constant P ):
f us HA
T
f us GA = f us HA
T f us SA
= f us HA
T
1
= f us HA
Tf us,A
Tf us,A
(6.400)
When this is substituted into Eq. (6.399), we have:
ln (xA ) =
f us HA
R
1
1
T
Tf us,A
f us HA
R
Tf us,A T
Tf us,A T
(6.401)
f us HA
f us HA
f
Tf us,A T =
2
2
RTf us,A
RTf us,A
|
{z
}
(6.402)
=Tf
292
1
1 2
x x3 ... when 1 < x < 1
2
3
x when higher order terms are ignored.
ln (1 x) = x
f us HA
(6.403)
(6.404)
Where Tf indicates the change (i.e., depression) in freezing point. In this context
the concentration is often given in terms of molal concentrations m (i.e., moles of
solute per kg of solvent; not mass here!) and M in kg mol1 . The relationship
between molal concentration and mole fraction for B is:
nB /(MA nA )
mB
nB
=
=
mB M A
(6.405)
xB =
nA + nB
(nA + nB )/(MA nA )
1/MA + mB
where MA is the molar mass of A. The last approximation applies to dilute solution. Substitution of Eq. (6.405) into (6.404) gives (Kf is called the freezing point
constant):
Tf =
RTf2us,A MA mB
f us HA
= Kf mB with Kf =
RTf2us,A MA
f us HA
(6.406)
293
The elevation of boiling point
By using analogous approach to obtain the depression of freezing point, we can
obtain the following expressions for the elevation of boiling point:
Tb =
2
RTvap,A
vap HA
xB = Kb mB where Kb =
2
RTvap,A
MA
vap HA
(6.407)
Example. Denote the solvent by A and the solute by B. Derive an expression for
the molar mass of B in terms of (Kb , Tb ) or (Kf , Tf ). Assume that the mass of
B dissolved in A and the total weight of the solution are given and that the solution
is dilute in B.
Solution. Since the solution is dilute in B, we have approximately mass(A)
mass(A + B), where mass(A) is the mass of solvent A (kg) and mass(A + B) is
the mass of the solution containing both A and B. Since mB = moles of solute
molecules dissolved / the total mass of the solvent, we can write:
mB =
nB
nB
mass(A)
mass(A + B)
294
By combining the above equations, we get:
nB =
Tf
Tb
mass(A + B) or nB =
mass(A + B)
Kf
Kb
M.W. of B =
Kf
mass(B)
mass(B)
=
in units of kg mol1
nB
Tf
mass(A + B)
Kf (K kg mol1 )
5.12
40
7.27
1.86
Kb (K kg mol1 )
2.53
3.04
0.51
295
When a solid solute (A) is left in contact with a solvent (B), it dissolves until
the solution is saturated. In this case saturation is a state of equilibrium with the
undissolved solute in equilibrium with the dissolved solute. Therefore, in a saturated
solution the chemical potential of the pure solid solute (A (s, T )) and the chemical
potential of A in solution (A (l, T )) are equal. The derivation below is essentially
the same as given earlier in Eqs. (6.398) through (6.401).
By using Eq. (6.369), we can write for A:
A (l, T ) = A (l, T ) + RT ln (xA )
Because there is equiblirum between the solid and dissolved forms of A, we have:
A (l, T ) = A (s, T )
Furthermore, pure solids have activity equal to that of standard state:
A (s, T ) = A (s, T ) A (l, T ) = A (s, T )
Inserting this into the first equation, we get:
A (s, T ) = A (l, T ) + RT ln (xA )
296
If we solve for ln (xA ), the expression becomes (with G = H T S):
ln (xA ) =
A (s, T ) A (l, T )
f us GA
f us HA
f us SA
=
=
+
RT
RT
RT
R
Note that above refers to the standard state of A (1 bar pressure and pure solute).
At constant P and T , we can use Eq. (3.176) and replace f us S:
f us HA
f us HA
f us HA
ln (xA ) =
+
Tf us,A
RT
RTf us,A
f us HA
1
1
ln (xA ) =
R
T
Tf us,A
f us SA =
R
T
TA,f us
(6.408)
298
The force f~ between two point charges (denoted by 1 and 2) is given by the
Coulombs law:
1
Q1 Q2
f~ =
~
r
(7.409)
40 r r 2
where:
r
Q1
Q2
0
r
~
r
=
=
=
=
=
=
When only the magnitude of the force is considered, the above equation becomes:
Q1 Q2
f =
(7.410)
40 r r 2
The electric field strength at point 1 is given by the ratio between force and the
charge Q1 :
~
~ = f
(7.411)
E
Q1
299
For each well behaved vector field (i.e. conservative) there is a potential such
that:
~ =
~
E
(7.412)
Note that here = (x, y, z) and E is a vector-field with components (Ex (x, y, z),
~ is an example of a vector operator:
Ey (x, y, z), Ez (x, y, z)). Gradient ()
~ =
,
,
x y z
Example. Calculate the gradient vector for function f (x, y, z) = x2 + y 2 + z 2 and
find the point where the length of the gradient vector becomes zero.
Solution. Calculate the partial derivatives of f with respect to x, y and z:
f
f
f
= 2x,
= 2y,
= 2z
x
y
z
~ = 0.
~ = (2x, 2y, 2z) . When x = 0, y = 0, z = 0, then f
f
300
Definition of 1 Volt (V): In presence of the electric field (from Q2 ), the difference
between the electric potential at two points is equal to the work per unit charge
required to move a charged test particle from one point to the other. Thus 1 V = 1
J C1 . The choice of zero potential is arbitrary but is usually chosen to correspond
to infinite separation of charges.
The electric potential at point r is the work required to bring a unit positive
charge from infinity to r. Combining Eqs. (7.409) and (7.412) and integrating from
infinity to r gives:
Zr
dr
Q2
Q2
=
(7.413)
(r) =
40 r
r 2
40 r r
n i zi = 0
(7.414)
i=1
301
N
Ps
i=1
i dni +
N
Ps
i dQi
i=1
dQi = zi F dni
(7.415)
Ns
X
i=1
i dni +F
Ns
X
i=1
zi i dni = SdT +V dP +
Ns
X
i + F zi i dni
{z
}
i=1 |
i
(7.416)
This means that the chemical potential i can be taken to contain the electrical
work as well:
i = i + zi F i
(7.417)
The equilibrium condition Eq. (6.347) still applies: i () = i () for species i and
phases and .
302
(7.418)
(7.419)
303
304
Cells without junctions:
Examples:
Pt(s) | H2 (g) | HCl(aq) | AgCl(s) | Ag(s)
Pt(s) | H2 (g) | HCl(aq) | Cl2 (g) | Pt(s)
Examples:
+
Zn(s) | Zn+
2 (aq) : Cu2 (aq) | Cu(s)
+
+
Zn(s) | Zn2 (aq) : Zn2 (aq) | Zn(s)
Ag(s) | AgCl(s) | Cl (aq) :: Ag+ (aq) | Ag(s)
: = liquid junction
:: = salt bridge
Galvanic cell with junction (Daniel
cell).
305
An example of liquid junction:
306
The electromotive force (EMF) of a cell is denoted by E (V) and it is defined by:
d (wel )
P =
= EI
(7.420)
dt
where wel is the electrical work done by the system (negative), P is the rate at
which chemical energy is converted to electrical energy (W) and I is the current
(A).
Task: What is the relation between the electromotive force in a cell and the chemical
potentials (or ion activities) in the cell?
To do this, lets consider an example (L = left, R = right):
PtL | H2 (g) | HCl(aq) | AgCl(s) | Ag(s) | PtR
Note that the convention is that reduction takes place on the right electrode and oxidation occurs in the left electrode.
The electrode reactions are:
2AgCl(s) + 2e (PtR ) = 2Ag(s) + 2Cl (aq) (reduction)
(7.421)
(7.422)
307
Consider a cell without any external connections to the electrode (see, for example,
the previous galvanic cell without liquid junction example; i.e. open circuit). For
such a system in equilibrium, we must have (similar to Eq. (6.347)):
2(HCl, aq) + 2(Ag, s) + 2(e , PtL ) = (H2 , g) + 2(AgCl, s) + 2(e , PtR )
(7.423)
Each chemical potential here is given by Eq. (7.417): i = i + zi F i
For all neutral species (HCl, Ag, H2 , AgCl) i = i because zi = 0.
Right electrode: (e , PtR ) = (e , PtR ) F R
Left electrode: (e , PtL ) = (e , PtL ) F L
Now i (e , PtR ) = i (e , PtL ). To see this, consider that the left and right electrodes are in contact (i.e. L = R ) and we have equilibrium i (e , PtR ) =
i (e , PtL ). By considering the above two equations with these conditions, we get
i (e , PtR ) = i (e , PtL ).
By plugging these results into Eq. (7.423), we get:
2(HCl, aq) + 2(Ag, s) 2F L = (H2 , g) + 2(AgCl, s) 2F R (7.424)
or 2(HCl, aq) + 2(Ag, s) (H2 , g) 2(AgCl, s) = 2F (L R )
308
r G = 2(HCl, aq) + 2(Ag, s) (H2 , g) 2(AgCl, s)
By combining this with Eq. (7.424), we get:
r G = 2F (R L ) = 2F E with E = R L
(7.425)
Note that E is considered here in the limit of zero current. Since r G depends on
pressure and temperature, the cell potential difference depends on P and T as well.
Also the concentration of HCl affects the the potential difference.
Although the above derivation applies only to a special case, it is clear that the
general equation corresponding to Eq. (7.425) is:
r G = |ve | F E
(7.426)
Ns
X
vi i
i=1
i=1
{z
+RT
Ns
X
|ve |F E
vi ln (ai ) = |ve | F E + RT ln
Ns
Y
i=1
v
ai i
(7.427)
309
where E is the standard electromotive force of the cell (i.e. the EMF when the
activities of all components are one).
Eq. (7.427) is called the Nernst equation and is usually written as (Q = reaction
quotient):
Eq.(5.303)
Ns
Y
z}|{
RT
RT
vs
E=E
=
E
ai
ln
ln(Q)
(7.428)
|ve | F
|v
e| F
i=1
At 25 C this can be written:
E = E
8.3145 J K1 mol1 (298.15 K)
0.02569
ln(Q) = E
ln(Q) (7.429)
|ve |
|ve | 96485 C mol1
310
Example. Three different galvanic cells have standard EMFs E of 0.01, 0.1 and
1.0 V at 25 C. Calculate the equilibrium constants of the reactions that occur in
these cells assuming the charge number |ve | for each reaction is unity.
Solution. Use Eq. (7.430) to get the equilibrium constants:
K(0.01 V) = exp
!
96485 C mol1 (0.01 V)
= 1.476
8.3145 J K1 mol1 (298.15 K)
K(0.1 V) = 49.0
311
mi 0
(7.432)
312
dG = + dm+ + dm
(7.433)
Note that we cannot add cations or anions separately, we always get both when an
electrolyte is added to the solvent.
Electrolyte
Non-electrolyte
Strong electrolyte
Weak electrolyte
(7.435)
where = v+ + + v
Here is the chemical potential for the electrolyte, which can be determined experimentally.
313
Note: In the following we will omit the standard state value m from expressions.
This results in simplified expressions but it will give inconsistent units! For example,
mi /m becomes now just mi .
The chemical potentials for cation and anion are given by:
+ = + + RT ln (+ m+ ) (+ std state chemical potential for cation)
=
+ RT ln ( m )
(7.436)
(7.437)
(7.438)
(7.440)
(7.441)
314
= + v RT ln ( m )
with aAv+ Bv = ( m )
(7.442)
=
v
mv
v
v
v++ v
7.5 Debye-H
uckel theory
315
Electrolytes containing multiply charged ions (for example, Cu2+ ) have greater
effect on the activity constants of ions than singly charged ions (for example, H+ ).
To account for this, we define ionic strength I (G. N. Lewis):
I=
Ns
1X
1
2
mi zi2 =
m1 z12 + m2 z22 + ... + mNs zN
s
2 i=1
2
(7.443)
where mi is molality of ion i and zi is charge for ion i in units of |e| (signed quantity).
Why is the activity of ions reduced in a solution?
Note that because the Coulomb interaction has long range, it is nearly impossible
to prepare dilute electrolytic solutions.
316
In 1923 Peter Debye and Erich H
uckel were able to show that for dilute solutions
the activity coefficient i of ion i with a charge number zi is given by (for details,
see Phys. Z. 24:185, 305 (1923)):
log (i ) = Azi2 I
(7.444)
where I is the ionic strength given by Eq. (7.443) and:
A=
1
2.303
2NA msolv
V
1/2
e2
40 r kT
3/2
(7.445)
Erich H
uckel, German
chemist (1896 - 1980).
1
(v+ log (+ ) + v log ( ))
v+ + v
v+ z +
log ( ) = A
I
v+ + v
(7.446)
(7.447)
317
By using the total charge neutrality v+ z+ = v z , we get:
log ( ) = Az+ z I
(7.448)
1
2.303
2NA msolv
V
4 8.854 1012
e2
40 r kT
1
2.303
!
2 6.022 1023 mol1 (997 kg)
(1.000m3 )
!
2
1.602 1019 C
AsV1 m1 (78.54) 1.3861 1023 JK1 (298.2 K)
= 0.509 kg1/2 mol1/2
1
0.001 mol kg1 (1)2 + 0.001 mol kg1 (1)2 = 0.001 mol kg1
2
2
2
log (+ ) = log ( ) = Az+
I = Az
I = 0.509 kg1/2 mol1/2 (1)2
I=
318
Next Eq. (7.448) yields:
Az+ z I
1+B I
Azi2 I
or log (i ) =
1+B I
log ( ) =
(7.449)
319
Determination of activity coefficients using an electrochemical cell:
Consider the following reaction:
1
H2 (g) + AgCl(s) = HCl(aq) + Ag(s)
2
The charge number in this reaction is one and the EMF is given by Eq. (7.428):
!
RT
aHCl
E=E
(7.450)
ln p
F
PH2 /P
where is the HCl mean ionic activity and m is HCl molality. The last step was
obtained at 25 C.
Eq. (7.451) contains two unknowns, but they can be determined from a series measurements of the cell EMF with different HCl molalities. Rearranging Eq. (7.451)
gives:
320
E + 0.1183 log (m) = E 0.1183 log ( )
(7.452)
HCl
0.905
0.875
0.830
0.796
0.767
0.755
0.809
1.009
1.316
1.762
2.38
LiCl
0.904
0.873
0.825
0.790
0.757
0.740
0.774
0.921
1.156
1.510
2.02
NaCl
0.902
0.870
0.820
0.778
0.735
0.693
0.657
0.668
0.714
0.783
0.874
CsCl
0.899
0.865
0.807
0.756
0.718
0.628
0.544
0.495
0.478
0.473
0.474
In general, if activity coefficient is greater than one, positive deviation from Raoults
law occurs (A-A/B-B interaction stronger than A-B), otherwise negative deviation
occurs (A-B stronger than A-A/B-B).
321
(7.453)
If we consider reaction:
1
H2 (g) + AgCl(s) = HCl(aq) + Ag(s)
2
The equilibrium constant for this reaction is given by:
K = exp
K= p
!
96485 C mol1 (0.2224 V)
= 5745
8.3145 J K1 mol1 (298.15 K)
aHCl
PH2 /P
(7.454)
322
r S = |ve | F
E
T
r C P
= |ve | F T
E
T
r H = |ve | F E + |ve | F T
(7.455)
P
2E
T 2
(7.456)
(7.457)
Notes:
1. These standard thermodynamic properties have been expressed with respect
to a hypothetical standard state where the electrolyte has molality of 1 mol
kg1 and the interactions for the resulting ions correspond to infinite dilution.
For example, the true activity of 1 mol kg1 HCl solution is less than 1.
2. The standard thermodynamic properties have been tabulated (see the NIST
online database). The standard state for strong electrolytes is assumed to
consists of completely ionized (for example, HCl, NaCl, etc.). For weak electrolytes, usually two different standard states are given: fully ionized and
un-ionized (not dissociated).
3. The experiments can not measure individual ions separately (electroneutrality
condition). However, by using a convention: f G (H+ ) = f H (H+ ) = 0,
the properties of other ions can be calculated.
4. For strong electrolytes the standard Gibbs energy of formation can be obtained by summing the corresponding formation energies for the ions.
323
Example. Calculate the standard molar entropy of chloride ion in aqueous solution
at 298.15 K starting with the Gibbs energy of formation (131.23 kJ mol1 ) and
the enthalpy of formation (167.16 kJ mol1 ).
Solution. By using G = H T S we can calculate f S :
f S (Cl ) =
167.16 kJ mol1 131.23 kJ mol1
f H (Cl ) f G(Cl )
=
T
298.15 K
= 120.51 J K1 mol1
223.08 J K1 mol1
2
!
130.68 J K1 mol1
2
!
= 56.36 J K1 mol1
324
It should be remembered that the tabulated data (e.g. the NIST database) gives the
properties in the limit of zero ionic strength. At higher ionic strengths the extended
Debye-H
uckel law must be used. Some quantitative expressions are as follows:
f Gi (I) = f Gi (I 0) RT ln (i )
f H (I) = T 2
"
f Gi (I)/T
T
#
(7.458)
(7.459)
r G and r H can be obtained by using the above expressions and Eqs. (2.139)
and (5.321).
325
Standard electrode potentials can be used to calculate the EMF of a given electrochemical cell. The standard electrode potential of a cell can be obtained with
hydrogen electrode on the left and all components at unit activity. For example:
1
H2 (g) + AgCl(s) = HCl(aq) + Ag(s)
2
which can be measured to be 0.2224 V. Also the hydrogen electrode contribution is
taken to be (arbitrarily) zero:
H+ (aq) + e =
1
H2 (g) with E = 0 V
2
Thus the only contribution to the cell potential is only from the reaction:
AgCl(s) + e = Ag(s) + Cl (aq) with E = 0.2224 V
Note that these reactions are written as reduction reactions and therefore E can
also be called reduction potential. A brief listing of various standard electrode
potentials are given in the following table:
326
Table. Standard electrode potentials at 25 C.
Electrode
F | F2 (g) | Pt
Au3+ | Au
Pb2+ | PbO2 | Pb
Cl | Cl2 (g) | Pt
H+ | O2 (g) | Pt
Ag+ | Ag
Fe3+ ,Fe2+ | Pt
I | I2 (s) | Pt
Cu+ | Cu
OH | O2 (g) | Pt
Cu2+ | Cu
Cl | Hg2 Cl2 (s) | Hg
Cl | AgCl(s) | Ag
Cu2+ , Cu+ | Pt
Br | AgBr(s) | Ag
H+ | H2 (g) | Pt
D+ | D2 (g) | Pt
Pb2+ | Pb
Sn2+ | Sn
Ni2+ | Ni
Cd2+ | Cd
Fe2+ | Fe
Zn2+ | Zn
OH | H2 (g) | Pt
Mg2+ | Mg
Na+ | Na
Li+ | Li
E (V)
2.87
1.50
1.455
1.3604
1.2288
0.7992
0.771
0.5355
0.521
0.4009
0.3394
0.268
0.2224
0.153
0.0732
0.0000
0.0034
0.126
0.140
0.250
0.4022
0.440
0.763
0.8279
2.37
2.714
3.045
Electrode reaction
(1/2) F2 (g) + e = F
(1/3) Au3+ + e = Au
(1/2) PbO2 + 2H+ + e = (1/2)Pb2+ + H2 O
(1/2) Cl2 (g) + e = Cl
H+ + (1/4)O2 (g) + e = (1/2)H2 O
Ag+ + e = Ag
Fe3+ + e = Fe2+
(1/2)I2 + e = I
Cu+ + e = Cu
(1/4) O2 (g) + (1/2)H2 O + e = OH
(1/2) Cu2+ + e = (1/2)Cu+
(1/2) Hg2 Cl2 + e = Hg + Cl
AgCl + e = Ag + Cl
Cu2+ + e = Cu+
AgBr + e = Ag + Br
H+ + e = (1/2)H2
D+ + e = (1/2)D2
(1/2) Pb2+ + e = (1/2)Pb
(1/2) Sn2+ + e = (1/2)Sn
(1/2) Ni2+ + e = (1/2)Ni
(1/2) Cd2+ + e = (1/2)Cd
(1/2) Fe2+ + e = (1/2)Fe
(1/2) Zn2+ + e = (1/2)Zn
H2 O + e = (1/2)H2 (g) + OH
(1/2) Mg2+ + e = (1/2)Mg
Na+ + e = Na
Li+ + e = Li
327
The cell EMF E can now be expressed as a difference between the right and left
electrode potentials:
E = ER EL and E = ER
EL
(7.460)
Note that given the data in the previous table, it is now possible to use the Nernst
Eq. (7.428) to obtain E for a cell when the concentrations or partial pressures are
given. Furthermore, the table data indicates the polarities for the electrodes and
can be used for obtaining the equilibrium constants (K) via Eq. (7.430).
Rules for using the table data:
1. The half-cell reactions are written as reduction reactions. Consider, for
example, a cell consisting of Pt | H2 (g) | HCl(aq) | Cl2 (g) | Pt. The half-cell
reactions are:
2. Both reactions must be written with the same number of electrons. Note
that E does not depend on |ve |.
3. The standard EMF for the cell (E ) is obtained by using Eq. (7.460). In the
above example this gives E = 1.3604 V 0.0000 V = 1.3604 V. The overall
reactions is:
H2 (g) + Cl2 (g) = 2H+ + 2Cl with E = 1.3604 V
328
4. If, under the standard condition, the standard EMF is positive, the reaction
will spontaneously go from left to right. In this case the right electrode is
positive (+) and the left electrode negative () when the cell is operated as
a galvanic cell.
5. If, under the standard condition, the standard EMF is negative, the reaction
will spontaneously go from right to left. In this case the right electrode is
negative () and the left electrode positive (+) when the cell is operated as
an electrolytic cell.
6. The equilibrium constant for the cell can be obtained by using Eq. (7.430).
7. The cell EMF under non-standard condition can be obtained by using the
Nernst equation (Eq. (7.428)). Note that if the cell has a liquid junction,
the Nernst equation must be written in terms of ion species (see the
following example).
Example. Consider the following galvanic cell: Zn(s) | Zn2+ :: Cu2+ | Cu(s)
where :: denotes a liquid junction (T = 298.15 K). Assume that Cu2+ and Zn2+
have identical ionic strengths on both sides of the liquid junction.
(a) What is the cell reaction?
(b) What is the standard EMF for the cell?
(c) What is the value of the equilibrium constant?
(d) What is the expression for the equilibrium constant in terms of ion concentrations?
Solution. Parts a) and b). The electrode reactions are:
329
Right electrode:
Left electrode:
Overall reaction:
Cu2+ + 2e = Cu(s)
Zn2+ + 2e = Zn(s)
Zn(s) + Cu2+ = Zn2+ + Cu(s)
= 0.339 V
ER
= 0.763 V
EL
E =
E = ER
L
0.339 V (0.763 V) = 1.102V
Part c.
Eq. (7.430) gives:
K = exp
|ve | F E
RT
= exp
!
2 96485 C mol1 (1.102 V)
= 1.80 1037
8.3145 J K1 mol1 (298.15 K)
Part d. Solids have activities of one and hence (using Eqs. (7.431) and (7.444)):
2+
a Zn2+ a (Cu(s))
Zn
Zn2+ m Zn2+
a Zn2+
K=
=
=
a Cu2+ a (Zn(s))
a Cu2+
Cu2+ m Cu2+
Cu2+
330
Example. Demonstrate with the Cu/Cu+ /Cu2+ sequence that the value for E for
the reaction Cu2+ + 2e Cu (E1 = +0.339 V; reaction #1) cannot be obtained
simply by summing the E values from reactions Cu+ + e Cu (E2 = +0.521
V; reaction #2) and Cu2+ + e Cu+ (E3 = +0.153 V; reaction #3).
Solution. The correct way to proceed is to consider the standard reaction Gibbs
energies (r G ):
Reaction 1 (|ve | = 2): r G1 = |ve | F E = 0.674 F
Reaction 2 (|ve | = 1): r G2 = |ve | F E = 0.521 F
Reaction 3 (|ve | = 1): r G3 = |ve | F E = 0.153 F
Thus r G1 = r G2 + r G3 but E1 6= E2 + E3 .
331
Notes:
1. If the solvent is involved in the reaction, the solvent is usually treated on the
mole fraction scale rather than the molal (or molar) scale. In this case the
equilibrium constant expression can be written as:
v Y m,i mi vi
K = x,A xA A
(7.461)
m
i6=A
where A denotes the solvent, x,A is its activity coefficient on the mole fraction
scale, m,i is the activity coefficient of reactant i on the molal scale and m
is the standard state molality (1 mol kg1 ). For dilute solutions Eq. (7.461)
can be written approximately as:
Y m,i mi vi
(7.462)
K=
m
i6=A
Even though the solvent is left out from the equation, the Gibbs energy of
formation of the solvent must be included in calculating r G for the reaction.
2. Some species in aqueous solution may be listed in thermodynamic tables in
more than one way. For example, NH3 vs. NH4 OH or CO2 vs. H2 CO3 . In
many cases we dont know the extent of hydration because of the difficulty in
distinguishing the species in solution. The convention in the NBS tables is that
f G = f H = f S = 0 for hydration reactions B+nH2 O = B(H2 O)n
where B denotes the species in question. This means that either one of the
pair can be used in the calculations.
332
7.8 Determination of pH
= E (0.0591 V) log
a H+
P (H2 ) /P
= E (0.0591 V) log a H+
333
Combining this with Eq. (7.463), we get:
E 0.2802 V = (0.0591 V) pH
(7.465)
or in terms of pH:
pH =
E (0.2802 V)
0.0591 V
(7.466)
335
=
=
=
=
Note: For a given macrostate, many different microstates are possible. Usually only
macrostates are observable.
336
Ensemble: A hypothetical collection of non-interacting systems, each of which is
in the same macrostate as the system of interest. Although the members of the
ensemble are macroscopically identical, they may not be in the same microstate.
Measurement: A measurement of any macroscopic property consists of a time average over the measurement interval. Hence it involves an inherent time averaging
process.
Postulate: The measured time average of any macroscopic property of the system
is equal to the average value of that property in the ensemble. Note that this
replaces the need to use time averages in the calculations.
An ensemble of systems that all have constant macroscopic parameters T and V is
called canonical ensemble.
The internal energy of a macroscopic system is given as an ensemble average:
X
pi Ei
(8.467)
U = hEi i =
i=1
where hi denotes the ensemble average, Ei is the energy of system i and pi is the
statistical weight for state i (i.e., probability of state i).
337
We will show that the statistical weights pi are given by the Boltzmann probability:
pi =
eEi
Z
(8.468)
where Z is the canonical partition function (do not confuse this with the compressibility factor which was also denoted by Z earlier), which is defined as:
X
eEi
(8.469)
Z=
i=1
338
We need to determine the unknown function f in the above postulate. Consider
two subsystems, which are labeled 1 and 2, in an ensemble. All ensemble members
labeled 1 are macroscopically identical and all members of 2 are macroscopically
identical. However, 1 and 2 are not necessarily identical as they can differ in volume
V and composition {ni }. For 1 and 2 we then have:
p1,i = f (E1,i ) and p2,j = g(E2,j )
(8.470)
where p1,i is the probability for system 1 to be in microstate i (energy E1,i ) and
p2,j for system 2 to be in microstate j (energy E2,j ). Note that f and g are
not necessarily the same functions. Systems 1 and 2 can also be considered as a
combined single system with overall fixed volume, temperature and composition:
p1+2,k = h(E1+2,k )
(8.471)
where p1+2,k is the probability that system 1+2 is in microstate k. Since systems
1 and 2 are independent of each other, we have:
E1+2,k = E1,i + E2,j with k = (i, j)
(8.472)
Furthermore since the systems 1 and 2 are independent, their overall probability is
given as a product:
339
h(E1+2,k ) = h(E1,i + E2,j ) = f (E1,i )g(E2,j )
(8.473)
With x E1,i , y E2,j and z x + y we can rewrite the above equation as:
h(z) = f (x)g(y)
(8.474)
(8.475)
(8.476)
(8.477)
(8.478)
340
By combining Eqs. (8.477) and (8.478) we get:
f (x)
g (y)
=
(constant)
g(y)
f (x)
(8.479)
(8.480)
f (x) = CeE1,i
where C is a constant of integration and is a universal constant that is the same for
both systems 1 and 2. Note that may depend on temperature. The normalization
constant C has to be chosen such that the sum over probabilities gives unity:
X
1
1
CeE1,i = 1 C =
(8.481)
P E
Z
1,i
i=1
e
i=1
eEi
eEi
=
P
Z
eEj
j=1
(8.482)
341
Evaluation of U and P :
U =
pi Ei =
i=1
Ei eEi
i=1
(8.483)
X
X
Z
Ej eEj
eEj =
=
V,{ni }
V,{ni } j=1
j=1
(8.484)
V,{ni }
(8.485)
X
pi Pi
P =
(8.486)
i=1
where Pi is the pressure in the ith member of the ensemble. If we have noninteracting and adiabatic systems, we can write the internal energy differential in
two ways (for each member of the ensemble):
342
dU = dwrev = Pi dV
Ei
dU =
dV (total differential)
V
Ei
Pi =
V
(8.487)
where Ei is the energy of ensemble member i. This gives an expression for the
measured pressure:
1 X Ei Ei
(8.488)
e
P =
Z i=1
V {ni }
Partial differentiation of Z with respect to V gives:
Z
V
X
eEi
i=1
=
T,{ni }
Ei
X
i=1
eEi
V
!
(8.489)
T,{ni }
X
Ei
Ei
eEi
=
V
V T,{ni }
i=1
343
Hence we can finally write P as:
1
1 ln(Z)
Z
P =
=
Z V T,{ni }
V
T,{ni }
(8.490)
Evaluation of :
First we differentiate Eq. (8.485) with respect to V :
U
V
ln(Z)
V
P
= P
V
=
ln(Z)
T V
V
(8.491)
T
(8.490)
z}|{
=
(P )
(8.492)
344
P/T = [P/(1/T )] [(1/T )/T ] = [P/(1/T )] /T 2
Now we can combine Eqs. (8.491) and (8.492):
P
1
P
=
V
T (1/T ) V
Denote Y
1
T
=
=
=
Y
Y V P V
Y V
dY
(8.493)
(8.494)
In the last step we have noted that = (T ) only. The above differential equation
can be integrated:
d
dY
=
ln(Y ) = ln() + C
Y
=1/T
=k
z}|{
z}|{
1
Y = C =
kT
where k is the Boltzmann constant.
(8.495)
345
Now we can write complete expressions for Eqs. (8.485) and (8.490):
dT
ln(Z)
ln(Z)
=
T
V,{ni }
V,{ni } d
1
ln(Z)
ln(Z)
2 = kT 2
=
T
k
T
V,{ni }
V,{ni }
ln(Z)
P = kT
V
T,{ni }
U =
(8.496)
(8.497)
eEi /(kT )
eEi /(kT )
=
P
Z
eEj /(kT )
(8.498)
j=1
Note that if the system contains degenerate (i.e., have the same energy) states, each
of them must be included in Eq. (8.498) separately.
Evaluation of entropy S:
Consider a reversible process where only P V -work occurs. The 1st and 2nd laws of
thermodynamics can be combined: dU = T dS P dV . This gives:
346
dS = T 1 dU + P T 1 dV = d T 1 U + T 2 U dT + P T 1 dV
(8.499)
(8.501)
(8.502)
(8.503)
(8.504)
347
Once systems partition function Z has been found, Eqs. (8.496), (8.497) and
(8.498) can be used to calculate thermodynamic quantities. Exact calculation of
Z for real systems is difficult because atoms/molecules interact and may require
quantum mechanical caclculations (i.e., to solve the Schr
odinger equation). In the
following, we will consider monoatomic ideal gas (distinguishable atoms), which
means that we need to only consider translational motion of atoms (no vibration or
rotation) and that the total energy can be expressed as a sum of individual atoms:
Ei =
N
X
i,j
(8.505)
j=1
where Ei is the total energy for state i, N the number of atoms and i,j is the
energy of atom j in state i. The canonical partition function Z can now be written
as:
Z=
X
i=1
eEi =
N
P
j=1
i,j
(8.506)
i=1
Since the atoms are distinguishable and i runs over all possible states, we can rewrite
Eq. (8.506) as:
348
N times
}|
{
z
X
X
X
e(k +l +...+m )
...
Z=
=
k=1 l=1
X
k
k=1
m=1
X
l=1
el ...
(8.507)
em
m=1
X
ei and Z = z1 z2 ...zN
zj
(8.508)
i=1
If the atoms are not all alike but there are NA (do not confuse this with the
Avograros number) atoms of species A, NB atoms of species B, etc. then:
Z = (ZA )NA (ZB )NB ..
(8.509)
NA
where ZA =
eA,i , ...
i=1
Remember that this holds only for distinguishable particles (e.g., localized atoms
in a solid).
349
It turns out that quantum mechanics excludes some classical states for indistinguishable particles and hence a different form of Eq. (8.509) must be used. It can
be shown that at sufficiently high temperatures, the canonical partition function for
indisinguishable particles is given by:
X
zN
Z=
eEi
(8.510)
and z
N!
i=1
where we have assumed that most molecules are in different microscopic states (see
Physical Chemistry, Levine for more details). A mixture of species has then:
Z=
(ZA )NA
(ZB )NB
...
NA !
NB !
(8.511)
The above formula cannot be used, for example, at liquid helium temperatures
where the high temperature assumption does not hold.
350
For monoatomic gas, we need to only consider translational motion (no rotation or
vibration). In order to use Eq. (8.496), which here states that U = kT 2 (Ztr /T )V,{ni } ,
we have to calculate the translational partition function Ztr . If the atoms are treated
as particles in a box (see quantum mechanics notes), Eq. (8.509) can be written
as:
Ztr
N
ztr
1 X X X (h2 /(8m))(n2x /a2 +n2y /b2 +n2z /c2 )
e
=
=
N!
N ! n =1 n =1 n =1
x
(8.512)
where h is the Plancks constant (6.626076 1034 J s), m is the atom mass, nx ,
ny and nz are quantum numbers, and a, b, c define the volume where the atoms are
confined (the box). When the quantitized levels are close to each other, we can
replave the above summations by integrals, for example:
Z
n2
n2
X
1 8m 1/2
(h2 /(8m)) x
(h2 /(8m)) x
a2
a2 dnx =
e
e
a
(8.513)
2
2 h
n =1
x
Since =
1
kT
351
3/2
ztr = 2mkT /h2
V
3
3
ln(ztr ) = ln 2mk/h2 + ln(T ) + ln(V )
2
2
(8.514)
N
ztr
N!
= N ln(ztr ) ln(N !)
(8.515)
(8.516)
Chapter 9: Kinetics
Thermodynamics can tell us if a given process happens or not but it does not tell
us how fast.
353
Empirical chemical kinetics concentrates on modeling the rates of chemical reactions. The key observables are concentrations of the species participating the reaction as a function of time. A kinetic measurement consists of mixing the reagents
(i.e., initiation of the reaction) followed by monitoring of the concentrations as a
function of time. Both steps introduce restrictions on the time resolution of the
experiment.
Examples of mixing processes:
Use traditional syringes to mix solutions. A cheap and easy way but it is not
suitable for fast reactions (seconds timescale).
Use fast motorized syringes to quickly mix solutions (flow and stopped-flow
techniques). Requires more expensive instrumentation but allows for better
time resolution. In the flow method a moving detection zone determines the
measurement time whereas in the stopped-flow the stopping syringe is used to
control the liquid flow rate (timescale depends strongly on instrumentation).
For gaseous samples motorized valves can be used to mix gases in relatively
fast time scale (s timescale).
In photochemical reactions (flash photolysis), short laser pulses can be used
to initiate the reaction (fs - ms timescale; most commonly in ns - s regime).
354
Measurement of concentrations:
UV/VIS absorption spectroscopy (lower bound timescale in femtoseconds).
The Lambert-Beer law relates absorbance to concentration.
Fluorescence spectroscopy (lower bound timescale dictated by the radiative
lifetime; ns - s). Other spectroscopic techniques such as IR and Raman may
also be used. It may be difficult to obtain absolute concentrations with these
methods.
Nuclear Magnetic Resonance (NMR) and Electron Spin Resonance (ESR;
EPR) (lower bound timescale in the ns range; often in seconds). To obtain
absolute concentrations, standard samples must be used. ESR is used for
systems involving radical species.
Mass spectrometry, gas chromatography, liquid chromatography and related
methods (time resolution dictated by instrumental response).
Monitoring total pressure or density of gas (typically in millisecond - second
timescale). Note that this method works only if the number of moles of the
gaseous components change in the reaction. If the chemical equation is known,
this can be used to obtain the extent of chemical reaction ().
If the concentration measurement is too slow for the kinetic timescale, one can use
the quenching method to stop the reaction. This can be achieved, for example, by
rapid cooling, dilution, or acid-base neutralization. Once the reaction is stopped,
even a slow method for determining the concentrations can be applied.
355
Rate of chemical reaction tells us how fast the given reaction occurs. The consumption of each reactant and product are dictated by this rate. The rate of chemical
reaction r is always positive and is defined by:
1 d [i] (t)
r=
(9.517)
vi dt
where vi is the stoichiometric coefficient for i and [i] (t) denotes the concentration of
i at a given time. Note that the stoichiometric coefficients are negative if i is on the
left hand side of the chemical equation. The rate of consumption or production of i
is denoted by ri and directly given by its time derivative. Note, however, that this
rate must also be positive and thus may require changing the sign of the derivative.
Example. Consider a reaction of the form A + 2B 3C + D. What is the rate of
chemical reaction and what are the rates of the individual compounds?
Solution. The rate of chemical reaction is given by Eq. (9.517):
r=
rA =
1 d [C]
d [A]
1 d [B]
d [D]
=
=
=
dt
3 dt
dt
2 dt
d [A]
d [B]
d [C]
d [D]
, rB =
, rC =
, rD =
dt
dt
dt
dt
356
Example. The rate of change in molar concentration of CH3 radicals in the reaction
2CH3 (g) CH3 CH3 (g) is reported as d [CH3 ] /dt = 1.2 mol L1 s1 . What are
the rate of reaction and the rate of formation of CH3 CH3 ?
Solution. The rate of reaction is given by Eq. (9.517):
r=
1 d [CH3 ]
1
= 1.2 mol L1 s1 = 0.60 mol L1 s1
2
dt
2
1 d [CH3 ]
d [CH3 CH3 ]
=
= 0.60 mol L1 s1
dt
2
dt
(9.518)
where the proportionality constant k is called the rate constant for the reaction.
Note that the rate constant is independent of the concentrations but may depend,
for example, on temperature. For A and B reacting, Eq. (9.518) is called the rate
law of the reaction. The general form of rate law is r = f ([A] , [B] , ...) where f is a
general function.
357
It should be emphasized that the rate law of a given reaction can not usually be
inferred from the chemical equation for the reaction. Most chemical equations
refer to a situation where multiple reactions occur to give the indicated outcome.
Chemical equations where only one reaction occurs are called elementary reactions.
It will turn out that the rate law and the chemical equation can only be directly
related for elementary reactions. For this reason the rate law should always be
obtained experimentally.
Example. Chemical reactions even with simple stoichiometry can result in very
complicated rate law. For example, H2 (g) + Br2 (g) 2HBr(g) gives the following
experimental rate law:
r=
This is very different result that one would expect based on the chemical equation.
Based on this we can conclude that this reaction is not elementary.
Once the rate law is known, we can use it to predict the rate of reaction at any
point in time. As we will see later, it can also be used to predict the concentrations
of each components at any point in time. The rate law can also give us important
clues about the reaction mechanism itself as the two must be consistent. This would
be especially useful for elementary reactions.
358
Many reactions are found to have rate laws of the form:
r = k [A]a [B]b ...
(9.519)
where the powers (a and b) define the order of the reaction with respect to each
species. The overall order of reaction is given by the sum of the individual orders,
a + b + .... The reaction order does not have to be an integer and for many gas
phase reactions they are not. Another special case is when the reaction order is
zero (e.g., a = 0). This corresponds to zeroth order reaction, which means that the
rate of reaction does not depend on the concentration of that particular component
(e.g., [A]a = [A]0 = 1). The zeroth order behavior usually means that there is some
other mechanism restricting the reaction rather than concentration. This could be,
for example, catalyst surface area. Only heterogeneous reactions can have overall
reaction order of zero. When the reaction is not of the form given by Eq. (9.519),
the reaction orders are not defined and there the overall reaction order is not defined
either. An example of this situation was given on the previous slide.
Example. Consider two reactions with the following rate laws: r = k [A]1/2 [B]
and r = k [A]0 = k. What are the reactions orders with respect to A and B, and
what are the overall reaction orders?
Solution. The first reaction has reaction order of 1/2 with respect to A and 1 with
respect to B. The overall reaction order is then 1/2 + 1 = 3/2. For the second
reaction the reaction order for A is zero and the overall reaction order is zero as
well.
359
Determination of the rate law:
1. Isolation method. When the concentrations of all the reactants except one are
in large excess, it is possible to determine the rate law with respect to this
component. For example, considering reaction A + B where B is in excess,
the concentration of B stays approximately constant throughout the reaction.
In this case, the rate law can be written:
r = k [A]a [B]b k [A]a
(9.520)
360
Example. The initial rate of reaction depends on concentration of a substance A
as follows (A P ):
[A]0 (103 mol L1 )
r0 (107 mol L1 s1 )
5.0
3.6
8.2
9.6
17
41
30
130
What is the order of reaction with respect to A and what is the rate constant?
Solution. Log-log plot of the data is shown below.
The slope is 1.99 2 (2nd order) and log(k) = 1.8646 k 1.4 102 (units?).
361
(9.522)
This is a differential equation for the unknown functions [i] (t). Most common
solutions are considered below.
First-order reactions. For a reaction A products, we have vA = 1 and
for first-order reaction a = 1. Eq. (9.522) can now be written:
d [A]
= k [A]a = k [A]
(9.523)
dt
The solution to the above differential equation is an exponential function:
[A] = [A]0 ekt
(9.524)
This is called the integrated form of the rate law. The unit for k in this
equation is s1 . The half-life t1/2 (i.e., the time when half of the initial
concentration is left) for a first-order reaction is given by:
[A]
1
ln(2)
= = ekt1/2 t1/2 =
(9.525)
[A]0
2
k
The numerical value for ln(2) 0.693.
362
Second-order reactions.
(9.522) when a = 2:
(9.526)
= kt or [A] =
(9.527)
[A]
[A]0
1 + kt [A]0
The unit for the rate constant k in this case is M1 s1 (or L mol1 s1 ).
The half-life in this case is given by:
1
(9.528)
t1/2 =
k [A]0
which can be derived the same way as we did for the 1st order reaction.
Another type of second-order reaction is:
d [A]
= k [A] [B]
(9.529)
dt
where a = b = 1 and their sum gives two as is required for an over all secondorder reaction. The solution to this differential equation is given by:
[B] / [B]0
= [B]0 [A]0 kt
(9.530)
ln
[A] / [A]0
363
Note: When [A] = [B], Eqs. (9.527) and (9.530) become idential as they describe
the same rate law.
Example. Solve the differential equation given in Eq. (9.523).
Solution. We proceed first by rearranging the equation:
d [A]
= kdt
[A]
then we can integrate both sides:
Z[A]
d [A]
= k
[A]
[A]0
Zt
dt
[A]
[A]0
= kt
364
Example. Consider the following reaction:
CH3 N2 CH3 (g) CH3 CH3 (g) + N2 (g)
At 600 K the partial pressures of CH3 N2 CH3 were found as follows:
t (s)
P (102 Torr)
0
8.20
1000
5.72
2000
3.99
3000
2.78
4000
1.94
Show that this reaction is 1st order in CH3 N2 CH3 and find the rate constant at
600 K.
Solution. First we remember Eq. (1.11), which gives the relationship between the
partial pressure of an ideal gas with its mole fraction:
Pi =
ni
P = yi P
n
where ni is the number of moles of component i and n is the total number of moles
present. Partial pressure Pi is proportional to the concentration of i:
Pi =
PV
ni
ni
P =
=
n
n
V
ni
ni
RT
= ci
|{z}
V
V
constant
where ci is the concentration of species i and the gas mixture was assumed to follow
the ideal gas law.
365
Note that in the 1st order integrated rate law the above proportionality constant
drops out and therefore Eq. (9.524) can now be written:
Pi = Pi,0 ekt
Here Pi is the partial pressure of CH3 N2 CH3 . Taking natural logarithms of this
equation gives:
Pi
= kt
ln
Pi,0
Thus plotting ln (Pi /Pi,0 ) as a function of t should give a straight line.
366
Example. When gaseous ammonia decomposes on the surface of hot metal catalyst
(2NH3 N2 + 3H2 ), the half-life of of NH3 (g) was observed to depend on its initial
pressure as follows:
P0 (Torr)
t1/2 (s)
65
290
105
460
150
670
185
820
What is the reaction rate and what is the value of rate constant?
Solution. First we observe that the half-life depends on the initial concentration
and based on Eq. (9.525) this cannot be a first-order reaction. For a second-order
reaction the half-life should depend inversely on the initial concentration (see Eq.
(9.528)) but here we observe exactly the opposite trend and this cannot therefore
be a second-order reaction either. If we write the rate law in the form (zeroth-order
reaction):
d [A]
= k [A]0 = k
(9.531)
dt
The integrate rate law now becomes:
[A] [A]0 = kt
Inserting [A] = [A]0 /2 above, gives the half-life as:
[A]0
[A]0 [A]0 /2
=
t1/2 =
k
2k
(9.532)
(9.533)
367
Since this is a gas phase reaction, we replace the concentration by partial pressure
(as we did in the previous example). Note that in this example the proportionality
constant does not cancel out and we will get the rate constant in pressure units.
We expect the half-life to depend on the initial pressure of NH3 (g) as follows:
Pi = 2kt1/2
Therefore a plot of the partial pressure as a function of t1/2 should give a straight
line with the slope 2k.
The data yields a straight line and confirms
that this reaction is zeroth-order. The rate
constant (in pressure units) can be extracted
from the slope: k = slope/2 = 0.112 Torr s1 .
Note that linear regression analysis can also
yield error estimates. Here the standard error estimate for the slope is 0.003, which
gives 0.0015 Torr s1 error for the rate constant. Thus we would report the final result
as: k = 0.112 0.002 Torr s1 . Another important indicator is the correlation coefficient
(r 2 ), which for this case was 0.9998 indicating
a very good quality fit. Strong deviations from
1.0 would indicate a bad quality fit.
368
Example. Liquid phase reaction CH3 CH2 NO2 + OH H2 O + CH3 CHNO
2 is
second-order overall and the rate constant at 273 K is k = 0.652 M1 s1 . The
initial concentration for nitroethane is 4.00 mM and 5.00 mM for OH . How long
does the reaction have to proceed in order to consume 90% of the initial nitroethane
concentration?
Solution. Let us denote nitroethane by A and OH by B. First we solve Eq.
(9.530) for t:
[B] / [B]0
1
ln
t=
[A] / [A]0
k [B]0 [A]0
The current concentrations are [A] = [A]0 x and [B] = [B]0 x based on the
stoichimetry of the chemical equation. 90%
consumption of A (nitromethane) cor
responds to 10% being left: [A]0 x / [A]0 = 0.10. Solving for x gives x =
0.90 [A]0 = 3.60 mM. Then [A] = [A]0 x = 0.40 mM and [B] = [B]0 x = 1.40
mM. Inserting these values into the expression for t we get:
1
0.652 M1 s1 (5.00 103 M 4.00 103 M)
1.40 103 M/5.00 103 M
ln
= 1580 s
3
3
0.40 10
M/4.00 10
M
t=
369
Summary of the rate laws:
Order
Reaction
Rate law
AP
AP
AP
r=k
[A] = [A]0 kt
r = k [A]
[A] = [A]0 ekt
r = k [A]2
A+B P
3
A + 2B P
t1/2
[A]0
2k
ln(2)
k
1
k[A]0
[A]
0
[A] = 1+kt[A]
0
r=
k [A] [B]
[B]/[B]0
= [B]0
[A]/[A]0
2
= [A] [B]
[B][A]
ln [A][B]0
ln
r
[A]0 kt
[B][B]
0
0
+
([B]0 2[A]0 )[B][B]0
([B]0 2[A]0 )2
r = k [A]n
1
1
1
kt = n1
n1
n1
[A]
kt =
n2
AP
[A]0
2n1 1
(n1)k[A]n1
0
370
Previously we have assumed that the foward reaction is dominating and therefore
we have ignored the possiblity for the reaction go from right to left. When the
reaction approaches equilibrium, we have to consider both forward and backward
reactions.
1. First-order reactions. Consider equilibrium A B where reactions A B and
A B occur simultaneously at appreciable rate. The rate constant for the forward
reaction is denoted by k+ and backward reaction by k . Both reactions are taken
to be first order. For example, A and B could correspond to two different isomers
of a molecule. The forward rate r+ = k+ [A] and backward rate r = k [B]. The
rate law in terms of A can then be written as:
d [A]
= k+ [A] + k [B]
(9.534)
dt
where the first term on the right hand side is responsible for the disappearance of
A and the second term for producing more A. Given the initial concentration of A
as [A]0 , the following balance has to hold at all times: [A] + [B] = [A]0 . Then we
can write Eq. (9.534) as:
d [A]
(9.535)
= k+ [A] + k [A]0 [A] = (k+ + k ) [A] + k [A]0
dt
The solution to this differential equation is (verification by differentiation):
371
[A] =
k + k+ e(k+ +k )t
[A]0
k+ + k
(9.536)
Equlibrium condition is reached after we wait long enough for the system to settle
down (steady-state). Thus we take the limit of t above and obtain an
expression for the equilibrium concentration of A:
k [A]0
(9.537)
[A]eq =
k+ + k
Since the concentrations of A and B are related to each other, we can get the
equilibrium concentration of B:
k+ [A]0
k
[B]eq = [A]0 [A]eq = [A]0 1
=
(9.538)
k+ + k
k+ + k
If we approximate activities by concentrations, the equilibrium constant can be
written for this reaction as (see Eq. (5.301)):
[B]eq
K=
(9.539)
[A]eq
Inserting the equilibrium concentrations from Eq. (9.538) into Eq. (9.539), we get:
k+
(9.540)
K=
k
372
This is consistent with the statement saying that the forward and backward rates
must be equal (k+ [A]eq = k [B]eq ) at equilibrium. If the equilibrium constant K
and one of the rates are known then the remaining rate constant can be calculated.
For a general reaction, the equilibrium constant can be expressed as:
kb,+
ka,+
...
(9.541)
K=
ka,
kb,
We have now established an important link between the rates of reaction and thermodynamic equilibrium.
Example. 1,2-dimethylcyclopropane cis trans isomerization follows first-order
kinetics. At 726 K the trans isomer was observed to form according to:
t (s)
% trans
0
0
90
18.9
225
37.7
270
41.8
360
49.3
495
56.5
585
60.1
675
62.7
70.0
Calculate the equilibrium constant K and the rate constants k+ (foward) and k
(backward).
Solution. We can directly fit the data to Eq. (9.536). Since the two variables
appear there separately, it is possible to determine them separately from this data.
Once we get k+ and k , we can use Eq. (9.540) to calculate K.
[A]
k + k+ e(k+ +k )t
=
[A]0
k+ + k
373
The data was not given in concentration units but percentages of B. Concentrations
are proportional to % and the proportionality constant cancels out for first-order
reactions. The concentration of A is proportional to 100% % of trans. Fitting the
given data to Eq. (9.536) is not a trivial matter. The equation is non-linear, which
requires the non-linear least squares procedure. The result is shown below.
374
2. Relaxation kinetics. The term relaxation here denotes the return of a system to
equilibrium after some given perturbation. This way we can study the relaxation
kinetics near the equilibrium. An example of suitable perturbation could be a
temperature jump that alters the equilibrium constant and hence the system will
try reach the new equilibrium condition. It is possible to reach temperature jumps
of 5-10 K in s timescale. Other examples of perturbation are laser and microwave
fields. Based on the vant Hoff equation (Eq. (5.323)), the equilibrium constant K
depends on temperature provided that r H is non-zero, which we assume in the
following.
Consider equilibrium A B, which is a first-order reaction. If the temperature
changes, the equilibrium constant K = K(T ) changed as well. Based on Eq. (9.541),
the rate constants also depend on temperature k+ = k+ (T ) and/or k = k (T ).
The response of the system will be exponential as shown below.
Cosider first the equilibrium condition:
d [A]eq
= k+ [A]eq + k [B]eq
dt
d[A]
(9.542)
Since we are at equilibrium dteq = 0 and the forward and backward rates are
equal: k+ [A]eq = k [B]eq . Next the system is perturbed in such a way that the
reaction will try to reach new equilbrium concentrations given by [A] = [A]eq +
and [B] = [B]eq . The former equation also defines the differential d [A] = d.
Now the system starts evolving according to:
375
d [A]
= k+ [A] + k [B]
dt
(9.543)
where the rate constants k+ and k now refer to the new values after perturbation.
This equation can be written in terms of as:
d
(9.544)
= k+ [A]eq + + k [B]eq
dt
Since the forward and backward rates are equal, this simplifies to:
d
= (k+ + k )
dt
| {z }
(9.545)
1/,const.,>0
= 0 et/
where 0 is the change in concentration right after the temperature jump. From
the exponential relaxation, we can obtain the sum of the new rate constants. When
k
this is combined with the equilibrium condition, K = k+ , it is possible to obtain
376
The rate constants for most chemical reactions increase as temperature is increased.
At low temperatures it is possible that some reactions, which are based on the
quantum mechanical tunneling, have rates that do not depend on temperature.
The relationship between rate constants, temperature, and the activation energies
is often found to follow the Arrhenius law (Svante Arrhenius, Swedish chemist, 1859
1927):
k = AeEa /(RT )
(9.546)
where k is the rate constant, A is the pre-exponential factor or the frequency factor,
Ea is the activation energy for the reaction, R is the gas constant, and T is the
temperature. An alternative form for this equation is:
Ea
ln(k) = ln(A)
(9.547)
RT
The parameter Ea can be obtained from the slope when plotting ln(k) as a function
of 1/T and A from the intercept. These parameters are called Arrhenius parameters.
Eq. (9.547) can be also written at two different temperatures as:
k1
1
Ea
1
ln
=
(9.548)
k2
R
T1
T2
where the indices k1 and k2 refer to rate constants at T1 and T2 , respectively.
377
Derivation of the Arrhenius law. Consider reaction A B, which proceeds through
a transition state denoted by A . Transition state is a molecular configuration along
the reaction path where the energy has the maximum value. The reaction can then
be written as:
k+
A A B
k
(9.549)
The equilibrium constant for the first part of the reaction can be written as:
[A ]
k+
=
(9.550)
K=
k
[A]
It is plausible to assume that [A ] << [A] since the intermediate A should be
transient in nature. The concentration of the product B follows the first-order
kinetics:
d [B]
(9.551)
= k [A ] = kK [A] = k [A]
|{z}
dt
k
where we used Eq. (9.550) in the last step and we introduced rate constant k for
the formation of B from A. Note that the rate constant k applies for the whole
reaction A B. Next we apply the vant Hoff equation (Eq. (5.323)):
378
r H
Ea
(ln(K))
=
T
RT 2
RT 2
(9.552)
where K is the equilibrium constant between the initial and the transition states
and r H is the reaction enthalpy between these states. In the last step we have
approximated r H Ea . Since K = k /k, Eq. (9.552) becomes:
(ln(k))
Ea
(ln(k ))
=
(9.553)
T
T
RT 2
Arrhenius assumed that the formation of the product B from the transition state
is independent of temperature and hence the second term on the left is zero:
Ea
(ln(k ))
=
(9.554)
T
RT 2
Solution to this differential equation is given by:
k = AeEa /(RT )
(9.555)
379
Example. Rasemization of pinene is a first-order reaction. The following rate
constants were determined in the gas phase: k457.7 K = 3.7 107 s1 and
k510.1 K = 5.1 105 s1 . Use the Arrhenius law to obtain Ea and calculate
the rate constant k480 K .
Solution. We use Eq. (9.548) at the two given temperatures:
5.1 105 s1
1
1
Ea
ln
3.7 107 s1
457.6 K
8.315 J K1 mol1 510.1 K
Solving for Ea gives Ea = 182 kJ mol1 . The same equation can be used to obtain
the rate constant at 480 K:
5.1 105 s1
1
182 kJ mol1
1
ln
=
1
1
k480 K
510.1 K
480 K
8.315 J K
mol
Solving for k480
Note: When the Arrhenius equation fails, it is often observed that the following
modified equation works well in practice (n is an empirical constant):
k = A T n eEa /(RT )
380
(9.557)
Trimolecular elementary reaction. These reactions would follow third order kinetics.
However, they are very rare as three molecules must to collide at the same time.
381
k2
AI P
Next we derive the expressions for [A], [I], and [P ]. The rate of unimolecular
decomposition of A into I is:
d [A]
= k1 [A]
(9.558)
dt
The intermediate I formed from A according to:
d [I]
= k1 [A] k2 [I]
dt
(9.559)
(9.560)
(9.561)
382
When the above result is substituted into Eq. (9.559) and we solve the resulting
differential equation, we get ([I]0 ):
k1
(9.562)
ek1 t ek2 t [A]0
[I] =
k2 k1
Since we must have at all times [A] + [I] + [P ] = [A]0 , we can solve for the concentration of P :
k1 ek2 t k2 ek1 t
[P ] = 1 +
[A]0
(9.563)
k2 k1
1
Concentrtion (M)
0.8
0.6
0.4
0.2
-0.2
0
6
Time (s)
10
12
383
If one of the consequtive steps is much slower than the other(s), it is said to be
the rate determining step. For example, when k2 >> k1 all I that is formed will
rapidly convert into P . This can be seen by first noting that now: ek2 t << ek1 t
and k2 k1 k2 . Eq. (9.563) then reduces to:
[P ] 1 ek1 t [A]0
(9.564)
This shows that the formation rate of P depends only on the smaller of the rate
constants (k1 above). In general the elementary reaction with the smallest rate
constant in a given reaction scheme is the rate determining step.
384
Substitution into Eq. (9.560) gives then:
d [P ]
= k2 [I] k1 [A]
dt
(9.567)
This shows that P in this case would be formed through the first-order decay of A
with rate constant k1 (corresponding to the rate determining step). Integration of
this equation gives directly:
[P ] = 1 ek1 t [A]0
(9.568)
Example. Consider the following reaction:
CH2 (CN)2 + Br2 BrCH (CN)2 + H+ + Br
The reaction mechanism has been determined as follows:
S
SH
}|
{
}|
{ k+ z
z
+
CH2 (CN)2 CH (CN)
2 +H
k
SBr
z
}|
{
z
}|
{
k
BrCH (CN)2 +Br
CH (CN)
2 +Br2
385
Solution. The reactions that appear in the reaction mechanism must be elementary
reactions. Therefore we can write the kinetic differential equations as:
d S
dt
d [SH]
= k+ [SH] + k S H+
dt
= k+ [SH] k S H+ k S [Br2 ]
d [SBr]
= k S [Br2 ]
dt
The steady-state approximation for S can be obtained by setting d S /dt = 0.
The second equation above then gives:
S =
k+ [SH]
k H+ + k [Br2 ]
This can be used to eliminate S from the other kinetic equations:
d [SBr]
k+ k [SH] [Br2 ]
d [SH]
=
=
dt
dt
k H+ + k [Br2 ]
386
Pre-equilibrium conditions. If one of the consequtive steps involve both forward and
backward reactions, we have to account for this in the kinetic differential equations.
Consider the following reaction:
k+
A + B I P
k
(9.570)
Thus this is effectively a second-order rate law with an effective rate constant:
k+ k
d [P ]
= k [A] [B] with k = kK =
(9.571)
dt
k
Note: The above calculation can also be carried out without the assumption k >>
k by using the steady-state approximation for the intermediate. In this case the
effective rate constant can be obtained as:
387
k =
k+ k
k + k
(9.572)
(9.573)
This type of reactions are found, for example, in enzyme catalysis (the MichaelisMenten mechanism).
Excercise. Derive Eqs. (9.572) and (9.573) by using the steady-state approximation for I.
Notes:
There are may topics in chemical kinetics that are not covered by the above
discussion. These include complex reactions such as chain reactions
(including polymerization), explosions, photochemical reactions, catalytic
reactions, and oscillating reactions.
A given reaction mechanism may contain many reactions that are
consequtive/parallel etc. and it often becomes impossible to find analytical
solutions to the corresponding differential equations.
A comprehensive kinetics database has been compiled by NIST and is
available online at: http://kinetics.nist.gov/kinetics/.