CHEN 623 Lecture Notes (2021 08 31)
CHEN 623 Lecture Notes (2021 08 31)
APPLICATIONS OF
THERMODYNAMICS TO CHEMICAL
ENGINEERING
James C. Holste
Professors Kenneth R. Hall, Philip T. Eubank and David M. Ford made substantial contributions to the
philosophy and content of this course. I am most grateful to them for their contributions.
TABLE OF CONTENTS
1 INTRODUCTION ............................................................................................................... 1
1.1 Basic Definitions for Thermodynamics ................................................................................. 1
2 THERMODYNAMIC LAWS............................................................................................. 7
2.1 ZEROTH LAW OF THERMODYNAMICS ........................................................................ 8
2.2 FIRST LAW OF THERMODYNAMICS ........................................................................... 12
2.2.1 Single Particle ............................................................................................... 13
2.2.2 Collections of Particles ................................................................................. 19
2.2.3 Consequence of the First Law of Thermodynamics ..................................... 30
2.2.4 General Energy Balance ............................................................................... 31
2.3 SECOND LAW OF THERMODYNAMICS ...................................................................... 35
2.3.1 Physical Observations ................................................................................... 35
2.3.2 Formal Statements of the Second Law ......................................................... 36
2.3.3 Consequences Of Second Law...................................................................... 46
2.3.4 Summary Statement of Second Law ............................................................. 66
2.3.5 Irreversible Processes: .................................................................................. 67
2.3.6 Practical Applications - Efficiencies ............................................................. 69
2.3.7 Summary Statement of Combined First and Second Laws .......................... 74
2.4 THIRD LAW OF THERMODYNAMICS .......................................................................... 75
2.4.1 Nernst-Simon Statements of Third Law ....................................................... 75
2.4.2 Lewis-Randall Statement of Third Law ........................................................ 75
2.4.3 Consequences of Third Law ......................................................................... 76
2.4.4 Practical Thermometry.................................................................................. 78
2.5 MICROSCOPIC INTERPRETATIONS: TEMPERATURE, ENTROPY ........................ 90
2.5.1 Microscopic (Statistical) Viewpoint of Entropy ........................................... 90
2.5.2 Microscopic (Statistical) Viewpoint of Energy and Temperature .............. 103
2.6 MATHEMATICAL IMPLICATIONS OF THERMODYNAMIC LAWS ...................... 109
2.6.1 Functions ..................................................................................................... 109
2.6.2 Multivariable calculus ................................................................................. 114
dQrev Q
Second Law: d(nS) =
T
; Δ(nS) universe ≥ 0 ; ΔSsurr = −
total
Tsurr
− ∑ ni Si
all
streams
Energy Functions: H = U + PV
A = U − TS
G = U + PV − TS = H − TS = A + PV
Useful Mathematical Relationships:
⎛∂f ⎞ ⎛∂f ⎞ ⎛ ∂z ⎞ ⎛ ∂z ⎞ ⎛ ∂w ⎞
MR1. f = f ( x, y ) ↔ df = ⎜ ⎟ dx + ⎜ ⎟ dy MR4. ⎜ ⎟ = ⎜
⎝ ∂x ⎠ y ⎝ ∂y ⎠ x ⎝ ∂x ⎠ y ⎝ ∂w ⎟⎠ y ⎜⎝ ∂x ⎟⎠ y
−1
∂z ⎡ ∂x ⎤ ⎛ ∂z ⎞ ⎛ ∂z ⎞ ⎛ ∂z ⎞ ⎛ ∂w ⎞
MR2. ⎛⎜ ⎞⎟ = ⎢⎛⎜ ⎞⎟ ⎥ MR5. ⎜ ⎟ = ⎜ ⎟ + ⎜
⎝ ∂x ⎠ y ⎣⎝ ∂z ⎠ y ⎦ ⎝ ∂x ⎠ y ⎝ ∂x ⎠ w ⎝ ∂w ⎟⎠ x ⎜⎝ ∂x ⎟⎠ y
⎛ ∂z ⎞ ⎛ ∂x ⎞ ⎛ ∂y ⎞ ⎛ ∂2 f ⎞ ⎛ ∂2 f ⎞
MR3. ⎜ ⎟ ⎜ ⎟ ⎜ ⎟ = −1 MR6. ⎜ =
⎝ ∂x ⎠ y ⎝ ∂y ⎠ z ⎝ ∂z ⎠ x ⎝ ∂x ∂y ⎟⎠ ⎜⎝ ∂y ∂x ⎟⎠
Residual Functions:
Ψ r = Ψ (T , ρ ) − Ψ ig (T , ρ ) Ψ R = Ψ (T , P ) − Ψ ig (T , P )
ρ ρ
⌠
Ur UR 1 ⎛ ∂Z ⎞ d ρ ⌠ ∂Z
⎛ ⎞ dρ Hr HR
= = I1 = ⎮⎮ ⎜ ⎟ = −T ⎮⎮ ⎜ ⎟ = I1 + Z − 1 =
RT RT T ⎮⌡ ⎝ ∂(1 / T ) ⎠ ρ ρ ⌡ ⎝ ∂T ⎠ ρ ρ RT RT
0 0
ρ
Ar AR ⌠ dρ Gr GR
= + ln Z = I 2 = ⎮⎮ ( Z − 1) = I2 + Z − 1 = + ln Z
RT RT ⌡ ρ RT RT
0
Sr SR P rV PV R
= I1 − I 2 = − ln Z = Z −1 =
R R RT RT
Property Changes:
Ψ ( 2 ) − Ψ (1) = Ψ R (T2 , P2 ) − Ψ R (T1, P1 ) + Ψ ig (T2 , P2 ) − Ψ ig (T1, P1 ) = Ψ R (T2 , P2 ) − Ψ R (T1, P1 ) + Δ1→2 Ψ ig
= Ψ r (T2 ,V2 ) − Ψ r (T1,V1 ) + Ψ ig (T2 ,V2 ) − Ψ ig (T1,V1 ) = Ψ r (T2 ,V2 ) − Ψ r (T1,V1 ) + Δ1→2 Ψ ig
Ideal Solution:
Δ ISM S Δ IS G Δ IS A Δ ISM U Δ ISM H PΔ ISM V
= − M = − M = − ∑ xi ln xi = = =0
R RT RT RT RT RT
E IS
Excess Properties: M ≡ ΔM M − ΔM M
Activity, Activity Coefficient, Fugacity, Fugacity Coefficient:
f̂i xiφ̂i αi f̂ φ̂
α i = γ i xi = = γi = = i† = i
fi † φi† xi xi fi φi
∞
⌠⎡ ⎤ dV
∂( nZ ) ⎞
P
⎮ ⎛
ln φ̂i = ⎮⎮ ( Z i − 1)
⌠ dP
= ⎮⎢ − 1⎥ − ln Z
⎮ ⎢⎜ ⎟
∂ni ⎠ T ,nV ,n
⌡ P ⎮ ⎝ ⎥ V
0 ⌡⎣ j≠i ⎦
V
( )
†
⎛ GE ⎞ GiE ⎛ ∂ nG
E
⎞
⎜⎝ RT ⎟⎠ = ∑ xi ln γ i
†
ln γ i =
RT
=⎜ ⎟
RT ⎝ ∂ni ⎠
T ,P,ni≠ j
Note: The superscript † denotes a reference state, which usually is taken to be the pure material at the same temperature and
pressure as the mixture.
Phase Equilibrium:
β β
Equilibrium Condition: µiα = µi fˆiα = ˆfi
Fugacity coefficient description: fˆ = y φˆ P
i i i
Activity coefficient description: fˆi = xi γ i† fi†
⎡ ⎤
y γ LRφ σ Pσ ⎢ P ⎥
⎢ 1 ⎥
Partition coefficient: K i ≡ i = i Vi i exp ⎢⎢ Vi L dP ⎥⎥
⌠
⎮
⎮
φ̂i P
⌡
xi ⎢ RT
⎢⎣ Piσ ⎥
⎥⎦
∑ν G i i =0
⎛ f̂ ⎞ ⎛ f̂ ⎞
Gi = Gi ° + RT ln ⎜ i ⎟ = Δ f Gio + RT ln ⎜ i ⎟ RHS: free elements are standard states for all substances
⎝ fi ° ⎠ ⎝ fi ° ⎠
νi
Δ r G° ⎛ f̂ ⎞
Standard Free Energy of Reaction and Equilibrium Constant: = − ln K = − ln ∏ ⎜ i ⎟
RT i ⎝ fi ° ⎠
( )
m
= Gi + ∑ λ k aik = 0 fk = ∑ nio − ni aik m = number of substances, w = number of elements
∂ni k=1 i=1
Ionic Solutions:
1000ni
mi =
n0 M 0
⎛ f̂ ◊ ⎞
Δ f Gi◊ = Δ f Gio + RT ln ⎜ i o ⎟
⎝ fi ⎠
M ν + Aν − = ν + M z+ + ν − A z−
γ ν± mν±
Ionic Equilibrium Constant: Ke =
γ u mu
ν RT
Galvanic Cells: ∑ν G i i = −zFE Gi = ν + Δ f G+◊ + ν − Δ f G−◊ + ν RT ln γ ± m± E = E° −
zF
ln γ ± m±
I
Extended Debye-Hückel Theory: ln γ ± = −1.174 z+ z−
1+ I
−1
Physical Constants: R = 8.314 J mol ; F = 96, 500 coulombs/(mole electrons)
External Field:
⎡Ξ ⎤
f̂i field = f̂i no field exp ⎢ i ⎥
⎣ RT ⎦
where Xi is the potential energy of molecule i due to the presence of the external field
Gravitational Field:
dP Mg ⎡ M g ( z − z0 ) ⎤
=− ; f̂i ( z ) = f̂i no gravity exp ⎢ i ⎥
dz V ⎣ RT ⎦
Boundary: separates the system from the surroundings. The most general
definition of a boundary is a mathematical surface.
NOTE: The term adiabatic is often used to describe situations where there
is no energy transferred as heat, but other forms of energy transfer
may occur.
Open: the amount of material (mass) within an open system may vary with
time because mass transfer across the system boundary may occur.
State variables: any of the physical properties used to describe the state of the
system. Each state of the system has a unique set of state variables.
Extensive properties: properties that are additive in the sense that the value of
the property for the whole system is the sum of the values for each of its
constituent parts. For extensive variables, there is a single aggregate value
for the entire system.
dφ
=0
dt
⎛ ∂φ ⎞
⎜⎝ ⎟⎠ = 0
∂t x
Time: a measure of the order in which events will occur in nature. Time is not a
legitimate variable in classical thermodynamics.
Law: a statement of the behavior of the physical world for which no contradiction
ever has been observed.
Work: energy transferred across a system boundary due to any driving force
other than a temperature gradient, excluding effects associated with mass
transfer.
Statement: If two bodies (closed systems) each are in thermal equilibrium with a
third body, then they are in thermal equilibrium with each other.
Proof:
For simplicity, consider a single-phase fluid so that only two intensive variables
need be specified to specify the state of the system completely.
F1( x A ,y A ,x C ,y C ) = 0
F2 (x B ,y B ,x C ,y C ) = 0
y C = f 1 ( x A ,y A ,x C ) and y C = f 2 ( x B ,y B ,x C )
f 1 ( x A ,y A ,x C ) = f 2 ( x B ,y B ,x C )
The Zeroth Law requires that A and B also are in equilibrium, therefore we also
have
F3 ( x A ,y A ,x B ,y B ) = 0
f1( x A ,y A ,x C ) = g1( x A ,y A ) Ω ( x C ) + Λ ( x C )
f2 ( x B ,y B ,x C ) = g 2 ( x B ,y B ) Ω ( x C ) + Λ ( x C )
g1( x A ,y A ) = g 2 ( x B ,y B ) = g 3 ( x C ,y C )
This function then has the same numerical value for each body if the bodies are
in thermal equilibrium.
Definition: q = g (x,y)
NOTE: For systems other than single-phase pure fluids, the temperature is a
function of the variables chosen to specify completely the state of the system
(independent variables).
CHEN 623 Lecture Notes (Date of last revision: 08/30/2021) 10
Many physical properties have been used to measure temperatures. Some of
the more common are:
2.2.1.1.1 Momentum
Consider the case where the force on an object depends only upon time
F = f (t )
Then
du
f (t ) = m ⇒ f (t )dt = m du
dt
Integration yields
∫ f (t ) dx = ∫ m du = mu + Constant of Integration
Consider the case where the force on an object depends only upon the position
of the object,
F = f (x )
Use the chain rule of differential calculus and the definition of velocity to write
du du dx du
= =u
dt dx dt dx
Then
du
f (x ) = mu ⇒ f (x )dx = mu du
dx
Integration yields
⌠
⎮f
⌡ ( x ) dx = ⌠⎮⌡ mu du = 21 mu 2 + Constant of Integration
Rearrange to obtain
1 mu 2
2
− ⌠⎮⌡ f ( x ) dx = E ( = constant )
CHEN 623 Lecture Notes (Date of last revision: 08/30/2021) 14
The combination of Newton's Law and the condition that the force depends only
on position leads to a constant of integration that does not vary with respect to
time, and
E = KE + PE
where: E = total energy
KE = 21 mu 2 = kinetic energy
2.2.1.2.1 Momentum
Fx dt = mdu x → p x − ∫ Fx dt = ℘x = Constant x
Fy dt = mdu y → p y − ⌠⌡ Fy dt = ℘y = Constant y
Fz dt = mdu z → p z − ∫ Fz dt = ℘z = Constant z
du x ⎛ du ⎞ ⎛ dx ⎞ du
Fx = m = m ⎜ x ⎟ ⎜ ⎟ = mu x x
dt ⎝ dx ⎠ ⎝ dt ⎠ dx
⎛ du y ⎞ ⎛ dy ⎞
du y du y
Fy = m = m⎜ ⎜ ⎟ = mu y
dt ⎝ dy ⎟⎠ ⎝ dt ⎠ dy
du z ⎛ du ⎞ ⎛ dz ⎞ du
Fz = m = m ⎜ z ⎟ ⎜ ⎟ = mu z z
dt ⎝ dz ⎠ ⎝ dt ⎠ dz
and
Fx dx = mu x du x → 1 mu 2 − ⌠⎮⌡ Fx dx = E x
2 x
Fy dy = mu y du y → 1 mu 2
2 y − ⌠⎮⌡ Fy dy = E y
Fz dz = mu z du z → 1 mu 2 − ⌠⎮⌡ Fz dz = E z
2 z
1m
2 (u 2
x ) ( ⌠
)
+ u y2 + u z2 − ⎮ Fx dx + Fy dy + Fz dz = E x + E y + E z
⌡
=E
The two terms on the left hand side of the equation can be expressed much more
concisely in terms of vector notation (scalar or dot product):
u x2 + u y2 + u z2 = u •u
Fx dx + Fy dy + Fz dz = F • dr
The energy equation for a single particle in three dimensions is given by:
p •p ⌠
E = 2 m u • u − ⌡ F • dr =
1 ⌠
⎮ − ⎮⌡ F • dr
2m
Total Energy: E = ∑ E i
r i is the location of particle i
where
total
msystem = ∑m i
i
Therefore
∑m i R i = 0
i
and
∑m i U i = 0
i
Therefore, the total kinetic energy of the system contains a term that contains the
motion of the center of mass and a term containing motion relative to the
center of mass, but no terms that contain both types of motion.
This means that the kinetic energy contains a term that is associated with the
macroscopic motion of the system, and a term that results from the motion of
the particles within the system.
1. Internal forces: forces that arise from interactions between particles within
int
the system, F i
2. External forces: forces that arise from interactions of the individual particles
ext
with external force fields, F i
Now, if the internal force acting on particle i because of the presence of particle j
depends only upon the distance between the particles,
( ))
( )
fik ri − rk = fik rcm + R i − rcm + Rk ( (
= fik R i − Rk )
The total force on particle i resulting from its interactions with all the other
particles in the collection is obtained by summation:
int
Fi = ∑ fik Ri − Rk
k
( )
In this expression, the first and third terms depend upon positions relative to the
center of mass, while the second term depends primarily upon the position of the
center of mass.
ext
where F ( rcm ) is the total force acting upon the collection of particles due to
external fields.
The first two terms deal with the macroscopic velocity and position of the system.
The last three terms deal with the properties possessed by the individual
particles.
The term in square brackets is intractable for calculations of fluids because the
number of particles generally exceeds 1020.
ri
total
Usystem = 2∑ i i
1 mU
i
• Ui + 21 ∑ ∑
i k r
⌠
⎮ fik
⌡
(
Ri − Rk • d ri)
ref
ri
ext
+ ∑ ⌠⎮⌡
i r
( )
F i ri • d ri
cm
total
Kinetic energy E ktotal = 21 m system 2
u cm
rcm
Potential energy E ptotal = − ⌠⎮⌡
( )
F rcm • d rcm
rref
total
E system total
= Usystem + E ktotal + E ptotal = msystem
total ⎡U + E k + E p ⎤
⎣ ⎦
total
U, Ek and Ep are the internal, kinetic and potential energies per unit mass. (or
per mole).
Note that the internal energy, U, depends only on the state of the system,
therefore the internal energy is a state function.
There exists a function (which we call energy) that is constant with respect to
time for an isolated system.
(This function is quite useful for describing the physical behavior of practical
systems.)
⎣ (
d ⎡m U + E k + E p ⎤
⎦system
= ) ∑
all heat
dQi + ∑
all work
dWi
transfer mechanisms
mechanisms
+ ∑ (dm )i
all
xfer
(H + E
i
k
i )
+ E ip − ∑
all nuclear
(dm )nuc
i
c 2
streams reactions
⎣ (
Δ ⎡m U + E k + E p ⎤
⎦system
= ) ∑
all heat
Qi + ∑
all work
Wi
transfer mechanisms
mechanisms
+ ∑
all
⌠
⎮
⌡
(dm )xfer
i (H i + E k
i + E i )−
p
∑ Δm nuc
i c
all nuclear
2
streams reactions
Note that the use of the enthalpy in the mass transfer term automatically includes
the effect of injection and ejection work.
Work: energy that is transferred across the system boundary because of any
driving force other than a temperature gradient, again excluding energy
transfer associated with mass flow across the boundary.
• Energy never is transferred as heat from a cold body to a hot body with no
other change in the system or its surroundings.
• Certain physical states of an isolated system never can precede other physical
states of that system.
2.3.2.1 Kelvin
It is impossible to extract heat from a reservoir and convert it wholly into work
without causing other changes in the universe.
2.3.2.2 Clausius
2.3.2.3 Planck
No cyclic process is possible that results in the flow of heat from a single heat
reservoir and the performance of an equivalent amount of work.
• The statements of Kelvin and Clausius are equivalent, but the equivalence is
not obvious.
A heat engine is a conceptual device that is used to study Second Law effects.
It consists of a device that is capable of exchanging energy as heat with two
reservoirs, one at a higher temperature and one at a lower temperature, and of
exchanging energy as work with the surroundings. The schematic diagram for a
heat engine is shown below.
QH
QC
Low Temperature
Reservoir
2.3.2.5 Notation:
This choice of signs is consistent with that used for the general energy balance,
but not consistent with discussions of Second Law effects found in many other
references. Other references, especially those discussing mechanical
engineering applications, use the opposite sign convention for W.
CHEN 623 Lecture Notes (Date of last revision: 08/30/2021) 38
2.3.2.6 Proof of Equivalence of Kelvin and Clausius Statements:
To prove the equivalence of these two statements, we consider two heat engines
that are coupled and that operate between a reservoir at qH and another
reservoir at qC.
θH
QH
A QHB
WA WB
A B
QCA QCB
θC
H
QH
A QB
WA WB
A B
QCA QCB
θC
The net energy transferred as work for the coupled engine is:
W = WA +WB
The net energy transferred as heat from the hotter reservoir is:
Q H = Q AH + QBH
The net energy transferred as heat from the colder reservoir is:
Q C = Q AC + QBC
Q AC + Q AH +WA = 0 → WA = −Q AH
If the engines are coupled so that the net work is zero (all work delivered by one
engine is consumed by the other), then
WA = −WB
and
Q H = Q AH + QBH = −QBC
QC = QA
C
+ QBC = QBC
W = WA +WB = 0
When QBC > 0 , a net transfer of energy as heat occurs from the colder (qC) to the
hotter (qH) reservoir with no other change in the universe. This violates the
Clausius statement.
Then
WA + Q AH + Q AC = 0 → Q AH = −Q AC
Q AC = −QBC = −Q AH
QC = QA
C
+ QBC = 0
W = WA +WB = −Q H
QH = Q AH + QBH = −W
When W < 0 and QH > 0, heat is converted entirely to work, which violates the
Kelvin statement of the Second Law.
Nicolas Sadi Carnot invented these as a tool for analyzing the efficiencies of heat
engines.
QH
C
= f (θC ,θH )
Q
Proof:
Consider two Carnot cycles, one that utilizes an ideal gas as the working
substance and the other using condensing and evaporating steam as the working
substance.
The two working substances (fluids in this case) have very different properties,
therefore the amounts of energy transferred as heat per cycle are most likely
different, but it always is possible to choose numbers of cycles such that:
nG QGC ≅ nS QSC
P P
C
QG QC
S
V V
Ideal Gas Steam
QH
G QH
S
G WG S WS
C
QG QC
S
H H
Hot Reservoir: Qnet = nGQG + nSQSH
C C
Cold Reservoir: Qnet = nGQG + nSQSC
nSQSC = −nGQG
C
Then:
C
Qnet =0
C H
Wnet + Qnet + Qnet =0
C
When Qnet = 0 , this becomes: Wnet = −Qnet
H
(
− nSQSH + nGQGH ≥ 0 ) ⇒ (n QS
H
S )
+ nGQGH ≤ 0
The operation of the two devices can be reversed, that has the effect of reversing
all signs, i.e., multiplying through by -1. The reverse operation leads to
( nSQSH + nGQGH ≥ 0 )
Both conditions must be satisfied, and that only is true for the equality, so
nSQSH = −nGQG
H
nSQSC = −nGQG
C
therefore:
nSQSH H
−nGQG QSH H
QG
= ⇒ =
nSQSC C
−nGQG QSC C
QG
In addition, the cycles are completely arbitrary, so that only the reservoir
temperatures can be involved as variables.
Therefore
QH
C
= f (θC ,θH )
Q
QH TH
=−
QC TC
Proof:
3
4
P θH
2 5
θI
1 6
θC
g (θH )
f (θC ,θH ) =
g (θC )
QH
Cycle 2345: I
= f (θI ,θH )
Q
QI
Cycle 1256: C
= f (θC ,θI )
Q
By algebra:
QH QH QI
C
= I C
⇒ f (θC ,θH ) = f (θC ,θi )f (θI ,θH )
Q Q Q
Then
To determine the form of the function g(q), we use the ideal gas for its simplicity
(all substances must have the same result). The equation of state is
PV = Rθ
For the isothermal steps where energy is transferred reversibly as heat:
V4
H H
V4 ⌠ RθH ⎛ V4 ⎞
Q = −W = ∫ P dV = ⎮
⎮ dV = RθH ln ⎜ ⎟
V3 ⌡ V ⎝ V3 ⎠
V3
⎛V ⎞
Likewise: Q C = RθC ln ⎜ 6 ⎟
⎝ V1 ⎠
⎛V ⎞
ln ⎜ 4 ⎟
QH θ ⎝ V3 ⎠
= H
QC θC ⎛ V1 ⎞
ln ⎜ ⎟
⎝ V6 ⎠
CHEN 623 Lecture Notes (Date of last revision: 08/30/2021) 54
The adiabatic (Q=0) steps 123 and 456 must obey the First Law of
Thermodynamics: (differential form)
dU = dW + dQ
dW = −PdV ⎫
ig ⎪ Rθ
dU = CV dθ ⎬ ⇒ CVig dθ = −PdV = − dV
⎪ V
dQ = 0 ⎭
Then
θH V3
⌠ ig dθ ⌠ RdV
⎮ CV = −⎮
⌡ θ ⌡ V
θC V1
⎛V ⎞ ⎛θ ⎞
123: R ln ⎜ 3 ⎟ = −CVig ln ⎜ H ⎟
⎝ V1 ⎠ ⎝ θC ⎠
⎛ V6 ⎞ ig ⎛ θC ⎞ ig ⎛ θH ⎞
456: R ln ⎜ ⎟ = −CV ln ⎜ = CV ln ⎜
⎝ V4 ⎠ ⎝ θH ⎟⎠ ⎝ θC ⎟⎠
So
⎛V ⎞ ⎛V ⎞ ⎛V ⎞ V6 V1 V6 V4
ln ⎜ 6 ⎟ = − ln ⎜ 3 ⎟ = ln ⎜ 1 ⎟ ⇒ = ⇒ =
⎝ V4 ⎠ ⎝ V1 ⎠ ⎝ V3 ⎠ V4 V3 V1 V3
and
⎛ V4 ⎞ ⎛ V4 ⎞
ln ln
Q H θH ⎜⎝ V3 ⎟⎠ θH ⎜⎝ V3 ⎟⎠ θH
= = − = −
Q C θC ⎛ V1 ⎞ θC ⎛ V6 ⎞ θC
ln ⎜ ⎟ ln ⎜ ⎟
⎝ V6 ⎠ ⎝ V1 ⎠
QH TH
= −
QC TC
PV = RT
Constant volume gas thermometry with real gases provides a practical primary
device to establish the absolute temperature scale according to
dQrev
dS ≡ defines a state function
T
Proof:
QH QC
=−
TH TC
so that
QH QC
+ =0
TH TC
P •2
1 •
V
QiH QiC
+ =0
TH TC
The more Carnot cycles that are used to approximate the actual cycle, the more
accurate the approximation.
Also:
2
⎡ #cycles QiH ⎤ dQ H ⌠
lim ⎢ ∑ ⎥= ⎮
⎮
#cycles→∞ ⎢⎣ i TiH ⎥⎦ TH ⌡
1
1
⎡ #cycles QiC ⎤ ⌠ dQ C
lim ⎢ ∑ ⎥ = ⎮⎮ C
#cycles→∞ ⎢⎣ i TiC ⎥⎦ ⌡ T
2
Then
2 1
H
⌠ dQ ⌠ dQ C ⌠ dQrev
⎮
⎮
+ ⎮⎮ = ⎮⎮
⌡ TH ⌡ TC ⌡ T
1 2
The symbol ∫ denotes an integral along a closed path, i.e., the path begins and
ends at the same point.
dQrev
dS ≡
T
⌠ ⌠ dQ
dS
⎮
⎮
=
⎮
⎮
rev
=0
⌡ ⌡ T
In other words, the value depends only on the state of the system, not upon how
the system reached that state.
ΔS ≥ 0 for any allowed process in an isolated system, i.e., for any process that
does not violate the Second Law of Thermodynamics.
Proof:
For the adiabatic, reversible path from state 1 to state 2, dQrev = 0 at all points
along the path, Then:
2
⌠ dQ
S2 − S1 = ⎮
⎮
rev
=0
⌡ T
1
2
1Q 2 = ∫ dQ = 0
1
dQ
Note that dS ≠ in this case
T
U is a state function ⇒ − (U 2 − U1 ) = U1 − U 2
Then
If 2Q1 > 0, then Wnet < 0, and the cycle from 1 to 2 and return would violate the
Planck statement of the Second Law.
1
⌠ dQ
2 Q1 < 0 ⇒ <0
⎮
rev
⎮
⌡ T
2
Therefore:
1
⌠ dQ
S1 − S2 = ⎮⎮ rev
<0 ⇒ S2 > S1
⌡ T
2
ΔS ≥0
1 2
/ rev
dQ
dS ≡
T
is a state function.
For an irreversible process that is allowed, we can have 1Q2 = 0 even though
ΔS = S2 − S1 > 0 .
2
2 /
⌠ dQ
1Q 2 / =0
= ∫ dQ ⇒ ⎮
⎮ =0
1 ⌡ T
1
2
ΔS = ⌠
⎮ dS
⌡ > 0 because the process is irreversible
1
/
dQ
But, ΔS > 1Q 2 ⇒ dS > ⇒ /
TdS ≥ dQ
T
where the equality holds for reversible processes and the inequality for
irreversible processes.
/ rev = −PdV
For a volume change of the system: dW
The First Law holds for both reversible and irreversible processes, therefore:
dU = dQ/ rev + dW
/ rev (reversible process)
/ + dW
= dQ / (irreversible process )
CHEN 623 Lecture Notes (Date of last revision: 08/30/2021) 67
Second Law conclusions lead to
/ rev = TdS ⎫
dQ
/ < dQ
⎬ ⇒ dQ / rev / + δ = dQ
⇒ dQ / rev
/ < TdS ⎭
dQ
Therefore, for processes that accomplish a given change in internal energy, dU:
/ rev + dW
dU = dQ / rev
/ = (dQ
/ + dW
= dQ / rev − δ ) + dW
/
Therefore
/ rev = dW
dW / −δ ⇒ / = dW
dW / rev + δ
So, for all processes that are allowed by the Second Law of Thermodynamics:
Q ≤ Qrev
W ≥ Wrev
Thermodynamic Efficiencies:
W
Thermal efficiency: ηT ≡
QH
Wrev TC
Carnot efficiency: η ≡ = 1−
C QH TH
QC
Coefficient of Performance: COP ≡
W
TC
(COP )Carnot =
TH −TC
W
Work-producing machines (W < 0): ηp ≡
Wrev
Wrev
Work-consuming machines (W > 0): ηp ≡
W
Consider a case where the system goes from state 1 to state 2 and energy is
transferred as heat to a reservoir. (The temperature of the reservoir remains
constant during the process.)
total
Suniverse = Ssys
total
+ Sres
total
Then
/ total
dQ
dS total
universe = dS total
sys + dS total
res = dS total
sys −
Tres
/ total = TresdSsys
dQ total
−TresdSuniverse
total
/ total
Substitute the previous result for dQ
/ total = TresdSsys
dU total − dW total
−TresdSuniverse
total
(
W total = U 2total − U1total −Tres Ssys )
total
(
( 2 ) − Ssys
total
) (
(1) +Tres Sunitotal ( 2 ) − Sunitotal (1) )
The Second Law of Thermodynamics requires that
⎡⎣Suniverse
total
(2 ) − Suniverse
total
(1)⎤⎦ ≥ 0
(
W total ≥ U 2total − U1total −Tres Ssys )
total
(
(2 ) − Ssys
total
(1) )
total
Wmax = ΔU system
total
− Tres ΔSsystem
total
/ rev = TdS ⎫
dQ
⎬ ⇒ dU = TdS − PdV
/ rev = −PdV ⎭
dW
Only state functions are involved in this result, therefore it must hold for both
reversible and irreversible processes.
The summary statement of the physics involved in the First and Second Laws of
Thermodynamics is:
dU = TdS − PdV
Formal Statements:
lim C X = 0
T →0
⎡ 1 ⎛ ∂V ⎞ ⎤
lim α = lim ⎢ ⎜ ⎟⎠ ⎥ = 0
T →0 ⎝
T →0 ⎣V ∂T P ⎦
The functional dependence of the thermal expansion is the same as that of the
heat capacity.
Only 3He and 4He have liquid as an equilibrium phase at T = 0 K. The entropy
change of a vapor/condensed phase transition goes to infinity as T → 0 ,
therefore
lim P σ = 0 , where Ps is the sublimation or vapor pressure.
T →0
For the helium isotopes, both solid and liquid phases exist at absolute zero, so
the entropies of the solid and liquid phases must be identical at the equilibrium
pressure. Then:
dP σ
lim = 0 , where Ps is the melting (freezing) pressure.
T →0 dT
A device for which the relationship between a measured physical property and
the thermodynamic temperature is known exactly from theory.
T2
lim
(PV )2
=
P→0
T1 lim
P→0 (PV )1
A device for which the relationship between a measured physical property and
the thermodynamic temperature is not known exactly. Secondary thermometers
must be calibrated in some manner.
Abbreviations
ITS: International Temperature Scale
IPTS: International Practical Temperature Scale
CCT: Comité Consultatif de Thermométrie
CIPM: Comité International des Poids et Mesures
CIPM - 1887: This scale was based upon constant volume gas thermometry
using normal H2
ITS-27:
• Interpolating instruments:
platinum resistance thermometers
platinum-rhodium thermocouple
optical pyrometer
IPTS - 68:
EPT - 76:
ITS - 90:
†
Reduced pressure boiling point (P = 33,330.6 Pa)
The definition of the kelvin will undergo a fundamental change on May 20, 2019,
if the proposed revision of the SI system is accepted at the November 2018
meeting of the General Conference of Weights and Measures (CGPM). Rather
than using the triple point of water to fix the temperature scale the proposal
recommends that the energy equivalent as given by Boltzmann's equation be
used.
One consequence of this change is that the new definition makes the definition of
the kelvin depend on the definitions of the second, the metre, and the kilogram.
Temperature no longer is an independent definition.
.....
1 2 3 4 5 ..... N
The particles are completely independent – they do not interact with one another
in any fashion. There is no external field present. Thus, all configurations of the
system have identical energy. For convenience, set the energy equal to zero.
Examples:
Examples:
Because this number represents the total number of different ways that the
system may be arranged, and all arrangements have the same energy (i.e., 0),
we expect that it will be related to entropy.
S (E,V ,N ) = kB lnΩ
In the simple model system as defined, volume does not really exist, and energy
is always 0, so
Ω (E = 0,N ) = 2N
S = k B ln2N = Nk B ln2
For “real atoms” in three dimensions, the expression for entropy will be more
complex.
For our model system, consider macrostates being defined based on the same
total number of particles in the ‘up’ state. There are N + 1 different macrostates.
For large N, this is considerably smaller than the number of different microstates.
N N!
Total # of microstates = ∑ = 2N
nu =0 (N − nu )!nu !
Using the binomial theorem with r = 1 gives the total number of microstates
N N! N N!
r = (1+ r ) ⇒
N
∑ Total # of microstates = ∑ = 2N
n=0 (N − n )!n! nu =0 (N − nu )!nu !
Microstate picture.
From statistical mechanics, the average value of any property B , which has a
value of Bi in any given microstate i , can be calculated from a probabilistic
average
total # of
microstates
B = ∑ Bi pi
i =1
(B j − B )2 p j
macrostates
σ B2 ≡ (B − B )2 = ∑
i =1
What is the average reading, and how large are the expected fluctuations?
nu − nd = ∑ (nu − nd )
N (N − nu )!nu !
nu =0 2N
Noting that
nd = N − nu
N!
nu − nd = ∑ ( 2nu − N )
N (N − nu )!nu !
nu =0 2N
2 N N! N N N!
= N ∑ nu − N ∑
2 nu =0 (N − nu ) ! nu ! 2 nu =0 (N − nu ) ! nu !
2 ⎡ N −1⎤ N ⎡ N ⎤
= N ⎣
N2 ⎦ − N ⎣2 ⎦
2 2
=0
This result occurs because there is no energy incentive for the particles to be
either up or down.
10
2
0.1
6 5
4
nu - nd
2
0
p
0.01
6
4 -5
2
0.001 -10
0 2 4 6 8 10
nu
So
N!
2
σ ud
N
= ∑ ( 2nu − N )
2 (N − nu )!nu !
nu =0 2N
( )
1 N N!
= ∑ 4nu2 − 4Nnu + N 2
2N nu =0 (N − nu )!nu !
1 ⎧⎪ N 2 N! N N! 2
N N! ⎫⎪
= N ⎨4 ∑ nu − 4N ∑ nu +N ∑ ⎬
2 ⎪⎩ nu =0 (N − nu )!nu ! nu =0 (N − nu ) !nu ! nu =0 (N − nu )!n !
u ⎪ ⎭
1
{
= N 4 ⎡⎣N(N + 1)2N −2 ⎤⎦ − 4N ⎡⎣N 2N −1⎤⎦ + N 2 ⎡⎣2N ⎤⎦
2
}
= N(N + 1) − 2N 2 + N 2
=N
σ ud = N
Often, it is not so much the absolute size of a fluctuation that is important, but
rather the size relative to the size of the system. To put it on a relative basis, we
can write
σ ud 1
=
N N
So for a system with ~ 1022 particles, the typical fluctuation of (nu − nd ) / N away
from 0 will be ~ 10 −11.
For a system with ~ 102 particles, the typical fluctuation of (nu − nd ) / N away from
0 will be ~ 10 −1.
nu
0 20 40 60 80 100
-2
10
-6
10
-10
10
-14
10
p
-18
10
-22
10
-26
10
-30
10
0 2 4 6 8 10
nu
The behavior of relative fluctuations with system size has some implications for
nanotechnology.
In reality, all systems (including the particles in an ideal gas) are quantum
systems with discrete energy levels.
⎛ 1 ⎞ hω
Harmonic oscillator (solids): E j = ⎜ nj + ⎟
⎝ 2 ⎠ 2π
⎡ E ⎤
Boltzmann statistics: n Bj ∝ exp ⎢ − j ⎥
⎣ kT ⎦
1
Bose-Einstein statistics (even spin): n BE ∝
⎡E ⎤
j
exp ⎢ j ⎥ −1
⎣ kT ⎦
1
Fermi-Dirac statistics (odd spin): n FD ∝
⎡ Ej ⎤
j
exp ⎢ ⎥ +1
⎣ kT ⎦
Note:
Ej
As → 0 ("high" temperatures), n BE
j → n Bj , and n FD
j → n Bj
kT
"High" temperature is defined by the ratio (Ej/kT). For most cases this is true for
temperatures above 100 K.
Heat denotes energy transfer that occurs because the distributions among the
available energy levels are different for two objects that are in thermal contact.
Work denotes energy that is transferred to the system because of changes in the
energy levels in the system.
Entropy is related to the number of different ways that a system can have the
same total energy.
Consider a system with four particles distributed among two energy levels, with
energies E0 = 0 and E1 = 1. In the Boltzmann distribution (the simplest of several
distribution functions), the probability that an individual particle has energy Ei is
⎡ E ⎤ ⎡ E ⎤ ⎡ E ⎤
exp ⎢ − i ⎥ exp ⎢ − i ⎥ exp ⎢ − i ⎥
P (E i ) = ⎣ kT ⎦ = ⎣ kT ⎦ = ⎣ kT ⎦
⎡ E ⎤ ⎡ E ⎤ ⎡ E ⎤ ⎡ 1 ⎤
∑ exp ⎢ − k ⎥ exp ⎢ − 0 ⎥ + exp ⎢ − 1 ⎥ 1+ exp ⎢ −
all ⎣ kT ⎦ ⎣ kT ⎦ ⎣ kT ⎦ ⎣ kT ⎥⎦
energy
levels
The particles then are distributed as shown below for several different
temperatures.
E = 1: x xx xxx xxxx
E = 0: xxxx xxx xx x
E total = 0 :
S
0 0 0 0 ⇒ Ω =1 ⇒ = lnΩ = 0
kB
E total = 1:
1 0 0 0 ⎫
⎪
0 1 0 0 ⎪ S
⎬ ⇒ Ω=4 ⇒ = 5.54
0 0 1 0 ⎪ kB
0 0 0 1 ⎪⎭
E total = 2 :
1 1 0 0 ⎫
⎪
1 0 1 0 ⎪
⎪
1 0 0 1 ⎪⎪ S
⎬ ⇒ Ω=6 ⇒ = 7.17
0 1 1 0 ⎪ kB
⎪
0 1 0 1 ⎪
⎪
0 0 1 1 ⎪⎭
E total = 4 :
S
1 1 1 1 ⇒ Ω =1 ⇒ =0
kB
Note that states 3 and 4 cannot be reached by isolated systems moving toward
equilibrium.
2.6.1 Functions
Notation: y = f ( x )
y1 1•
y = f2(x)
x1 x2
NOTE:
y1 = f 1 ( x1) = f 2 ( x1)
y 2 = f1 (x 2 ) = f 2 (x 2 )
2
⌠
⎮ dy
⌡ = y 2 − y 1 = f1( x 2 ) − f1( x 1) = f 2 ( x 2 ) − f 2 ( x 1)
1
2 2
/ i = ⌠⎮⌡ f i ( x ) dx
Area under curve = ⌠⎮⌡ dA
1 1
dφ is an exact differential if the value of the integral does not depend upon the
choice of paths from 1 to 2.
d/ φ is an inexact differential if the value of the integral depends upon the choice
of paths from 1 to 2.
Notes:
• The energy transferred as heat or work depends upon the choice of paths from
/
state 1 to state 2, therefore dW / are inexact differentials.
and dQ
z = f ( x1,x 2 ,x 3,,x n )
or
F ( x1,x 2 ,x 3,,x n ,x n+1) = 0
Total Derivative:
dz = lim {f ( x1 + Δx1,x 2 + Δx 2 ,,x n + Δx n ) − f ( x1,x 2 ,,x n )}
Δx i →0
⎛ ∂f ⎞ ⎧ f ( x + Δx,y ) − f ( x,y ) ⎫
Definition of partial derivative: ⎜⎝ ⎟ = lim ⎨ ⎬
∂ x ⎠ y Δx →0 ⎩ Δx ⎭
Also: lim Δx i = dx i
Δx i →0
Therefore
⎛ ∂z ⎞
n
dz = ∑ ⎜ ⎟ dx i
i =1⎝ ∂ x i ⎠
x k ≠i
⎛ ∂z ⎞ ⎛ ∂z ⎞
z = f ( x,y ) ⇔ dz = ⎜ +
⎝ ∂ x ⎟⎠ y
dx ⎜⎝ ∂ y ⎟⎠ dy
x
f ( x,y,z ) = 0
The choice of independent variables in this case is arbitrary, so we could use any
of:
x = F ( y,z )
y = G ( x,z )
z = H ( x,y )
Then
⎛ ∂z ⎞ ⎛ ∂z ⎞
z ( x,y ) → dz = ⎜ +
⎝ ∂ x ⎟⎠ y
dx ⎜⎝ ∂ y ⎟⎠ dy
x
⎛ ∂y ⎞ ⎛ ∂y ⎞
y ( x,z ) → dy = ⎜ +
⎝ ∂ x ⎟⎠ z ⎜⎝ ⎟⎠ dz
dx
∂z x
Combine these by substituting dy from the second equation into the first
equation.
⎛ ∂z ⎞ ⎛ ∂z ⎞ ⎡⎛ ∂ y ⎞ ⎛ ∂y ⎞ ⎤
dz = ⎜ + +
⎝ ∂ x ⎟⎠ y ⎢⎜⎝ ⎟ ⎜⎝ ⎟
dx dx dz
⎜⎝ ∂ y ⎟⎠
x ⎣ ∂x ⎠z ∂ z ⎠ x ⎥⎦
⎡ ⎛ ∂z ⎞ ⎛ ∂y ⎞ ⎤ ⎡⎛ ∂ z ⎞ ⎛ ∂z ⎞ ⎛ ∂y ⎞ ⎤
dz ⎢1− ⎜ ⎟ ⎜ ⎟ ⎥ = dx ⎢⎜ ⎟ + ⎜ ⎟ ⎜ ⎟ ⎥
⎢⎣ ⎝ ∂ y ⎠ x ⎝ ∂ z ⎠ x ⎥⎦ ⎢⎣⎝ ∂ x ⎠ y ⎝ ∂ y ⎠ x ⎝ ∂ x ⎠ z ⎥⎦
x,z could be chosen as the independent variables, therefore the above equation
must hold for arbitrary choices of dx, dz.
−1
⎡ ⎛ ∂z ⎞ ⎛ ∂y ⎞ ⎤ ⎛ ∂z ⎞ ⎡⎛ ∂ y ⎞ ⎤
⎢1− ⎜ ⎟ ⎜⎝ ⎟⎠ ⎥ = 0 → ⎜⎝ ∂ y ⎟⎠ = ⎢⎜⎝ ∂ z ⎟⎠ ⎥
⎢⎣ ⎝ ∂ y ⎠ x ∂ z x ⎥⎦ x ⎣ x⎦
⎡⎛ ∂ z ⎞ ⎛ ∂z ⎞ ⎛ ∂y ⎞ ⎤ ⎛ ∂z ⎞ ⎛ ∂y ⎞ ⎛ ∂x ⎞
⎢⎜ ⎟ + ⎜ ⎟ ⎜ ⎟ ⎥ = 0 → ⎜⎝ ∂ y ⎟⎠ ⎜⎝ ∂ x ⎟⎠ ⎜⎝ ∂ z ⎟⎠ = −1
⎢⎣⎝ ∂ x ⎠ y ⎝ ∂ y ⎠ x ⎝ ∂ x ⎠ z ⎥⎦ x z y
dz ⎛ ∂ z ⎞ dx ⎛ ∂ z ⎞ dy
Divide by dw: =⎜ ⎟ +
dw ⎝ ∂ x ⎠ y dw ⎜⎝ ∂ y ⎟⎠ x dw
⎛ ∂z ⎞ ⎛ ∂z ⎞ ⎛ ∂x ⎞
Restrict to the case of constant y: ⎜⎝ ⎟ = ⎜ ⎟ ⎜ ⎟
∂w ⎠ y ⎝ ∂ x ⎠ y ⎝ ∂w ⎠ y
dz ⎛ ∂ z ⎞ dx ⎛ ∂ z ⎞ dy ⎛ ∂ z ⎞ dx ⎛ ∂ z ⎞
=⎜ ⎟ + =⎜ ⎟ +
dy ⎝ ∂ x ⎠ y dy ⎜⎝ ∂ y ⎟⎠ x dy ⎝ ∂ x ⎠ y dy ⎜⎝ ∂ y ⎟⎠ x
⎛ ∂z ⎞ ⎛ ∂z ⎞ ⎛ ∂x ⎞ ⎛ ∂z ⎞
When this is restricted to constant w: ⎜⎝ ∂ y ⎟⎠ = ⎜⎝ ⎟⎠ ⎜ ⎟ + ⎜⎝ ∂ y ⎟⎠
w
∂ x y ⎝ ∂ y ⎠w x
⌠
⎮
df
⌡ =0
Vector Notation:
F = (F1,F2 ,,Fn )
dx = (dx 1,dx 2 ,,dx n )
n
F • dx = ∑ F i dx i
i =1
⎛ ∂ ∂ ∂ ⎞
∇=⎜ , ,,
⎝ ∂ x1 ∂ x 2 ∂ x n ⎟⎠
⎛ ∂f ∂f ∂f ⎞
∇f = ⎜ , ,,
⎝ ∂ x1 ∂ x 2 ∂ x n ⎟⎠
⎡⎛ ∂F ∂F ⎞ ⎛ ∂F ∂F ⎞ ⎛ ∂F1 ∂F ⎞ ⎤
∇ × F = ⎢⎜ n − 2 ⎟ ,,⎜ i −1 − i +1 ⎟ ,,⎜ − n−1 ⎟ ⎥
⎣⎝ ∂x 2 ∂x n ⎠ ⎝ ∂x i +1 ∂x i −1 ⎠ ⎝ ∂x n−1 ∂x 1 ⎠ ⎦
n⎛ ∂f ⎞
Observe that df = ∑ ⎜ ⎟ dx i = ∇f • dx
i =1⎝ ∂ x i ⎠ x
k ≠i
⌠⌠
Then, by Stoke's theorem: ⌠
df
⎮
⌡
= ⌠
∇f
⎮
⌡ ( )
• d x = ⎮⌡⎮⌡ ∇ × ∇f • d A
∂2 f ∂2 f ⎡x 0 → xn ⎤
Where ( ∇ × ∇f )i = −
∂x i −1 ∂x i +1 ∂x i +1 ∂x i −1
; ⎢x ⎥
⎣ n+1 → x 1⎦
⌠ ⌠⌠ n ⎡ ⎤
⎮ ⎮⎮ ∂2 f ∂2 f
Then df
⎮ = ⎮⎮ ∑ ⎢ − ⎥ dAi
⎮ ⎮⎮ i =1⎢ ∂x i −1 ∂x i +1 ∂x i +1 ∂x i −1 ⎥⎦
⌡ ⌡⌡ ⎣
But, because ⌠⎮⌡ df = 0 for any state function for any path and surface, the term in
square brackets in the integrand always must be zero.
∂2 f ∂2 f
=
∂x i ∂x k ∂x k ∂x i
This requires that, if f is a state function, the order of differentiation with respect
to the independent variables does not matter.
⎛ ∂f ⎞
n
Math 1: df = ∑ ⎜ ⎟ dx i ⇔ f ( x 1,x 2 ,,x n )
i =1⎝ ∂ x i ⎠ x
k ≠i
−1
⎛ ∂z ⎞ ⎛ ∂y ⎞ ⎡⎛ ∂ z ⎞ ⎤ ⎛ ∂y ⎞
Math 2: ⎜⎝ ∂ y ⎟⎠ ⎜⎝ ∂ z ⎟⎠ = ⇔ ⎢⎜ ⎟⎠ ⎥ =⎜
⎝ ∂ z ⎟⎠ x
1
x x ⎢⎣ ⎝ ∂ y x⎥
⎦
⎛ ∂z ⎞ ⎛ ∂y ⎞ ⎛ ∂x ⎞
Math 3: ⎜⎝ ∂ y ⎟⎠ ⎜⎝ ∂ x ⎟⎠ ⎜⎝ ∂ z ⎟⎠ = −1
x z y
⎛ ∂z ⎞ ⎛ ∂z ⎞ ⎛ ∂x ⎞
Math 4: ⎜⎝ ⎟⎠ = ⎜⎝ ⎟ ⎜ ⎟
∂w y ∂ x ⎠ y ⎝ ∂w ⎠ y
⎛ ∂z ⎞ ⎛ ∂z ⎞ ⎛ ∂z ⎞ ⎛ ∂x ⎞
Math 5: ⎜⎝ ∂ y ⎟⎠ = ⎜⎝ ∂ y ⎟⎠ + ⎜⎝ ⎟⎠ ⎜ ⎟⎠
w x
∂ x y ⎝ ∂ y w
∂2 f ∂2 f
Math 6: =
∂x i ∂x k ∂x k ∂x i
The fundamental differential form of the summary of the First and Second Laws
of Thermodynamics (Physics Summary Statement) for a closed system is:
dU = TdS − PdV
⎛ ∂f n ⎞
Compare to Math 1: df = ∑ ⎜ ⎟⎠ dx i ⇔ f ( x 1,,x n )
i =1⎝ ∂ x i
The form of the Physics Summary Statement indicates that the "natural"
independent variables for internal energy are entropy and volume. In general,
manipulations will be more convenient when guided by the natural independent
variables.
Math 1 ® ⎛ ∂U ⎞ ⎛ ∂U ⎞
Then: dU = ⎜ +
⎝ ∂T ⎟⎠ P ⎜⎝ ⎟⎠ dP
dT
∂P T
∂U ⎞ ⎛ ∂U ⎞ with other
However, we cannot easily identify ⎛⎜ and
⎝ ∂T ⎟⎠ P ⎜⎝ ⎟
∂ P ⎠T
thermodynamic quantities.
Replace S with T:
τ 1 = U − C1x 1 = U −TS
d τ 1 = dU −TdS − SdT = TdS − PdV −TdS − SdT
d τ 1 = −SdT − PdV ⇒ τ 1 = τ 1 (T ,V )
τ 2 = U − C2 x 2 = U − ( −P )V = U + PV
d τ 2 = dU +VdP + PdV = TdS − PdV +VdP + PdV
d τ 2 = TdS +VdP ⇒ τ 2 = τ 2 (S,P )
Result: H (S ,P ) = U + PV (Enthalpy )
τ 12 = U − C1 x 1 − C2 x 2 = U −TS + PV = A + PV = H −TS
d τ 12 = −SdT +VdP ⇒ τ 12 = τ 12 (T ,P )
dU = TdS − PdV ⇔ U = U (S ,V )
dH = TdS +VdP ⇔ H = H (S ,P )
dA = −SdT − PdV ⇔ A = A (T ,V )
dG = −SdT +VdP ⇔ G = G (T ,P )
Math 1 requires:
⎛ ∂ (nU ) ⎞ ⎛ ∂ (nU ) ⎞ ⎛ ∂ (nU ) ⎞
m
d (nU ) = ⎜ d (nS ) + ⎜ d (nV )+ ∑ ⎜ ⎟ dn i
⎝ ∂ (nS ) ⎟⎠ nV ,{nk }
⎝ ∂ (nV ) ⎟⎠ nS,{nk } i =1⎝ ∂ ni ⎠ nS,nV ,n
k ≠i
In the first two terms, {ni} held constant implies that n also is held constant.
⎛ ∂ (nU ) ⎞ ⎛ ∂U ⎞ ⎛ ∂ (nU ) ⎞ ⎛ ∂U ⎞
Then =⎜ =⎜
⎝ ∂ S ⎟⎠V ,n ⎝ ∂V ⎟⎠ S,n
⎜⎝ ∂ (nS ) ⎟⎠ and ⎜⎝ ∂ (nV ) ⎟⎠
{ }
nV , ni { }
nS, ni
Notation: Define: µi ≡ ⎜ ( ) ⎟
⎛ ∂ nU ⎞
where µi is the chemical potential
⎝ ∂ n i ⎠ nS,nV ,n
k ≠i
m
Then: d (nU ) = Td (nS ) − Pd (nV ) + ∑ µi dn i
i =1
m
d (nU ) = T d (nS ) − P d (nV ) + ∑ µi dn i
i =1
m
d (nH ) = T d (nS ) + (nV ) dP + ∑ µi dn i
i =1
m
d (nA ) = − (nS ) dT − P d (nV ) + ∑ µi dn i
i =1
m
d (nG ) = − (nS ) dT + (nV ) dP + ∑ µi dn i
i =1
where:
Temperature:
⎛ ∂ (nU ) ⎞ ⎛ ∂U ⎞
T =⎜ =⎜
⎝ ∂ (nS ) ⎟⎠ ⎝ ∂ S ⎟⎠ n,V
{ }
nV , ni
⎛ ∂ (nH ) ⎞ ⎛ ∂H ⎞
=⎜ =⎜
⎝ ∂ (nS ) ⎟⎠ ⎝ ∂ S ⎟⎠ n,P
{ }
P, ni
Pressure:
⎛ ∂ (nU ) ⎞ ⎛ ∂U ⎞
−P = ⎜ =⎜
⎝ ∂ (nV ) ⎟⎠ ⎝ ∂V ⎟⎠ n,S
{ }
nS, ni
⎛ ∂ (nA ) ⎞ ⎛ ∂A⎞
=⎜ =⎜
⎝ ∂ (nV ) ⎟⎠ ⎝ ∂V ⎟⎠ n,T
{ }
T , ni
⎛ ∂H ⎞ ⎛ ∂G ⎞
V =⎜ =⎜
⎝ ∂ P ⎠ S,n ⎝ ∂ P ⎟⎠ T ,n
⎟
Entropy:
⎛ ∂A⎞ ⎛ ∂G ⎞
−S = ⎜ =
⎝ ∂T ⎟⎠V ,n ⎜⎝ ∂T ⎟⎠ P,n
Chemical Potential:
The Maxwell Relations follow from Math 6 applied to the differential forms of the
Physics Summary Statements:
⎛ ∂T ⎞ ⎛ ∂P ⎞
dU = TdS − PdV ⇒ ⎜⎝ ⎟ = − ⎜⎝ ⎟
∂V ⎠ S ∂ S ⎠V
⎛ ∂T ⎞ ⎛ ∂V ⎞
dH = TdS +VdP ⇒ ⎜⎝ ⎟ = ⎜ ⎟
∂P ⎠ S ⎝ ∂S ⎠ P
⎛ ∂S ⎞ ⎛ ∂P ⎞
dA = −SdT − PdV ⇒ ⎜⎝ ⎟ = ⎜ ⎟ (MxA)
∂V ⎠ T ⎝ ∂T ⎠V
⎛ ∂S ⎞ ⎛ ∂V ⎞
dG = −SdT +VdP ⇒ ⎜⎝ ⎟⎠ = − ⎜⎝ ⎟ (MxG)
∂P T ∂T ⎠ P
Begin with the differential form of the Gibbs summary statement for an open
system:
m
d (nG ) = − (nS ) dT + (nV ) dP + ∑ µi dn i
i =1
⎛ ∂ (nG ) ⎞
µi = ⎜ ⎟
⎝ ∂ n i ⎠ T ,P,n
k ≠i
⎛ ∂µi ⎞ ⎛ ∂ (nS ) ⎞
⎜⎝ ⎟ = −⎜ ⎟
∂T ⎠ P,n ⎝ ∂ n i ⎠ T ,P,n
k ≠i
⎛ ∂µi ⎞ ⎛ ∂ (nV ) ⎞
⎜⎝ ⎟⎠ = ⎜ ⎟
∂ P T ,n ⎝ ∂ n i ⎠ T ,P,n
k ≠i
⎛ ∂ (nX ) ⎞
Definition: Xi ≡ ⎜ ⎟ where: X i is the partial molar quantity
⎝ ∂ n i ⎠ T ,P,n
k ≠i
(Partial molar properties are very useful for describing mixture properties.)
For pure fluids (or mixtures at fixed composition) the convenient independent
variables are either (T,P) or (T,V) as the independent variables.
2
⌠⎡ ⎤
⎮ ⎛ ∂Ψ ⎞ ⎛ ∂Ψ ⎞
Math Relation 1: Ψ ( 2 ) − Ψ (1) = ⎢⎜ ⎟ dT + ⎜⎝ ⎟ dX ⎥
⎮ ⎝
⎮
∂T ⎠ X ∂X ⎠ T
⌡⎣ ⎦
1
where X denotes either P or V.
T2 X2
⌠ ⌠
⎛ ∂Ψ ⎞ ⎛ ∂Ψ ⎞
Ψ (T2 ,X 2 ) − Ψ (T1,X 1) = ⎮
⎜⎝ ⎟ dT +
⎮
⎜⎝ ⎟ dX
⎮
⎮
⌡ ∂T ⎠ X
1
⎮
⎮
⌡ ∂X ⎠ T
2
T1 X1
or
X2 T2
⌠ ⌠
⎛ ∂Ψ ⎞ ⎛ ∂Ψ ⎞
Ψ (T2 ,X 2 ) − Ψ (T1,X 1) = ⎮
⎮ ⎜⎝ ⎟⎠ dX +
⎮
⎮ ⎜⎝ ⎟⎠ dT
⎮
⌡ ∂X T
1
⎮
⌡ ∂T X
2
X1 T1
However, to calculate the values of these integrals we must relate the partial
derivatives in the integrals to properties that are measurable.
Examples: U, H, A, G, S.
⎛ ∂P ⎞
Bulk modulus (isothermal): BT ≡ −V ⎜
⎝ ∂V ⎟⎠ T
⎛ ∂P ⎞
Bulk modulus (adiabatic): BS ≡ −V ⎜
⎝ ∂V ⎟⎠ S
1 ⎛ ∂P ⎞
Thermal pressure coefficient: γ ≡ ⎜ ⎟
P ⎝ ∂T ⎠ V
dQ
Heat capacities: CX ≡
dT Constant X
⎛ ∂ Qrev ⎞
Definition: CV ≡ ⎜
⎝ ∂T ⎟⎠V
/ rev dU
dQ dV
Divide by dT: = +P
dT dT dT
⎛ ∂Qrev ⎞ ⎛ ∂U ⎞ ⎛ ∂V ⎞
Restrict to constant volume: ⎜⎝ ⎟ = ⎜ ⎟ + P ⎜⎝ ⎟
∂T ⎠V ⎝ ∂T ⎠V ∂T ⎠V
⎛ ∂U ⎞
Therefore: CV = ⎜
⎝ ∂T ⎟⎠ V
dU = TdS − PdV
⎛ ∂U ⎞ ⎛ ∂S ⎞ ⎛ ∂V ⎞
⎜⎝ ⎟ = T ⎜⎝ ⎟ − P ⎜⎝ ⎟
∂T ⎠V ∂T ⎠V ∂T ⎠V
⎛ ∂U ⎞ ⎛ ∂S ⎞
Therefore: CV = ⎜ = T
⎝ ∂T ⎟⎠ V ⎜⎝ ⎟
∂T ⎠ V
⎛ ∂ Qrev ⎞
Definition: CP ≡ ⎜
⎝ ∂T ⎟⎠ P
/ rev + dW
dU = dQ / rev = dQ
/ rev − PdV ⇒ dQ
/ rev = dU + PdV
⎛ ∂Qrev ⎞ ⎛ ∂U ⎞ ⎛ ∂V ⎞
⎜⎝ ⎟ = ⎜ ⎟ + P ⎜⎝ ⎟
∂T ⎠ P ⎝ ∂T ⎠ P ∂T ⎠ P
Note that
⎛ ∂H ⎞ ⎛ ∂U ⎞ ⎛ ∂V ⎞ ⎛ ∂P ⎞ ⎛ ∂Qrev ⎞
H = U + PV ⇒ ⎜ = + +V =
⎝ ∂T ⎟⎠ P ⎜⎝ ∂T ⎟⎠ P ⎜⎝ ⎟ ⎜⎝ ⎟ ⎜ ⎟
P
∂T ⎠ P ∂T ⎠ P ⎝ ∂T ⎠ P
⎛ ∂H ⎞
Therefore: CP = ⎜
⎝ ∂T ⎟⎠ P
dH = TdS +VdP
⎛ ∂H ⎞ ⎛ ∂S ⎞ ⎛ ∂P ⎞
⎜⎝ ⎟⎠ = T ⎜⎝ ⎟⎠ +V ⎜⎝ ⎟
∂T P ∂T P ∂T ⎠ P
⎛ ∂H ⎞ ⎛ ∂S ⎞
Therefore: CP = ⎜ = T
⎝ ∂T ⎟⎠ P ⎜⎝ ⎟
∂T ⎠ P
Therefore
⎛ ∂P ⎞ ⎛ ∂V ⎞ ⎛ ∂P ⎞ TV β 2
C P − CV = T ⎜ = TV β ⎜ = TV β BT =
2
⎝ ∂T ⎟⎠ V ⎜⎝ ∂T ⎟⎠ P ⎝ ∂T ⎟⎠ V κT
Use relations between heat capacity and entropy, the cyclic relation (Math
Relation 3), and the Maxwell relations derived from the Helmholtz and Gibbs
energies
⎛ ∂S ⎞ ⎛ ∂S ⎞ ⎛ ∂P ⎞ ⎛ ∂V ⎞ ⎛ ∂P ⎞
T⎜ − ⎜⎝ ⎟ ⎜ ⎟ ⎜ ⎟ ⎜ ⎟
CP ⎝ ∂T ⎟⎠ P ∂P ⎠ T ⎝ ∂T ⎠ S ⎝ ∂T ⎠ P ⎝ ∂T ⎠ S
= = =
CV ⎛ ∂S ⎞ ⎛ ∂S ⎞ ⎛ ⎞ ∂V ⎛ ∂P ⎞ ⎛ ∂V ⎞
T⎜ −
⎝ ∂T ⎟⎠V ⎜⎝ ⎟ ⎜ ⎟ ⎜ ⎟ ⎜ ⎟
∂V ⎠ T ⎝ ∂T ⎠ S ⎝ ∂T ⎠V ⎝ ∂T ⎠ S
Use the chain rule for partial derivatives (Math Relation 4) in the numerator and
the cyclic relation in the denominator, so that
⎛ ∂V ⎞ ⎛ ∂P ⎞ ⎛ ∂V ⎞ ⎡⎛ ∂P ⎞ ⎛ ∂V ⎞ ⎤
⎜⎝ ⎟ ⎜ ⎟ ⎜ ⎟
CP
⎜⎝ ⎟ ⎜ ⎟
∂T ⎠ P ⎝ ∂T ⎠ S ∂T ⎠ P ⎢⎣⎝ ∂V ⎠ S ⎝ ∂T ⎠ S ⎥⎦
=− =−
CV ⎛ ∂P ⎞ ⎛ ∂V ⎞ ⎡⎛ ∂P ⎞ ⎛ ∂V ⎞ ⎤ ⎛ ∂V ⎞
⎜⎝ ⎟⎠ ⎜⎝ ⎟⎠ − ⎢⎜ ⎟⎠ ⎜⎝ ⎟⎠ ⎥ ⎜⎝ ⎟⎠
∂T V ∂T S ⎣ ⎝ ∂V T ∂T P⎦ ∂T S
⎛ ∂P ⎞
C P ⎜⎝ ∂V ⎟⎠ S BS κ T
= = =
CV ⎛ ∂P ⎞ BT κ S
⎜⎝ ⎟⎠
∂V T
Class Form
⎛ ∂a ⎞
1 ⎜⎝ ⎟⎠
∂b c
⎛ ∂α ⎞
2 ⎜⎝ ⎟
∂b ⎠ c
⎛ ∂a ⎞
3 ⎜⎝ ⎟⎠
∂b γ
⎛ ∂α ⎞
4 ⎜⎝ ⎟⎠
∂β c
⎛ ∂α ⎞
5 ⎜⎝ ⎟⎠
∂b γ
⎛ ∂α ⎞
6 ⎜ ⎟
⎝ ∂β ⎠ γ
Class 1 Derivatives:
Procedure: The only useful relation for this class is the cyclic relation (Math
Relation 3):
⎛ ∂a ⎞ ⎛ ∂b ⎞ ⎛ ∂c ⎞
⎜⎝ ⎟⎠ ⎜⎝ ⎟⎠ ⎜⎝ ⎟⎠ = −1
∂b c ∂c a ∂a b
⎛ ∂P ⎞ ⎛ ∂P ⎞ ⎛ ∂T ⎞ γP
Example: ⎜⎝ ⎟ = − ⎜⎝ ⎟ ⎜⎝ ⎟ = −
∂V ⎠ T ∂T ⎠V ∂V ⎠ P βV
βκ T β
Note that this also leads to the conclusion that = =P
γ γ BT
∂α ⎞ ⎛ ∂H ⎞
General form: ⎛⎜ Example: ⎜⎝ ⎟
⎝ ∂b ⎟⎠ c ∂P ⎠ T
Procedures:
If a is S, then use the relations to the heat capacities and the Maxwell relations
directly.
Examples involving S:
⎛ ∂S ⎞ CP ⎛ ∂S ⎞ ⎛ ∂P ⎞
⎜⎝ ⎟ = ⎜⎝ ⎟ = ⎜ ⎟
∂T ⎠ P T ∂V ⎠ T ⎝ ∂T ⎠V
⎛ ∂S ⎞ CV ⎛ ∂S ⎞ ⎛ ∂V ⎞
⎜⎝ ⎟⎠ = ⎜⎝ ⎟⎠ = − ⎜⎝ ⎟
∂T V T ∂P T ∂T ⎠ P
dH dS dP
Divide through by dP: =T +V
dP dP dP
⎛ ∂H ⎞ ⎛ ∂S ⎞
Restrict to constant T: ⎜⎝ ⎟⎠ = T ⎜⎝ ⎟⎠ +V
∂P T ∂P T
⎛ ∂H ⎞ ⎛ ∂V ⎞
⎜⎝ ⎟⎠ = −T ⎜⎝ ⎟⎠ +V
∂P T ∂T P
⎛ ∂V ⎞
=V −T ⎜
⎝ ∂T ⎟⎠ P
=V [1−T β ]
dH dS dP
Divide through by dV: =T +V
dV dV dV
⎛ ∂H ⎞ ⎛ ∂S ⎞ ⎛ ∂P ⎞
Restrict to constant T: ⎜⎝ ⎟ = T ⎜⎝ ⎟ +V ⎜⎝ ⎟
∂V ⎠ T ∂V ⎠ T ∂V ⎠ T
⎛ ∂H ⎞ ⎛ ∂P ⎞ ⎛ ∂P ⎞
Apply Maxwell Relation (MxA): ⎜⎝ ⎟ = T ⎜⎝ ⎟ +V ⎜⎝ ⎟
∂V ⎠ T ∂T ⎠V ∂V ⎠ T
Note that the relations are more complex when a variable is used that is not a
"natural" independent variable.
∂a ∂T
General form: ⎛⎜ ⎞⎟ Example: ⎛⎜ ⎞⎟
⎝ ∂b ⎠ γ ⎝ ∂P ⎠ H
Procedure: Use the cyclic relation (Math Relation 3) to obtain two Class 2
derivatives:
⎛ ∂a ⎞ ⎛ ∂a ⎞ ⎛ ∂γ ⎞
⎜⎝ ⎟⎠ = − ⎜⎝ ∂γ ⎟⎠ ⎜⎝ ∂b ⎟⎠
∂b γ b a
⎛ ∂T ⎞ ⎛ ∂T ⎞ ⎛ ∂H ⎞ ⎛ 1 ⎞ ⎡ ⎛ ∂V ⎞ ⎤
⎜⎝ ⎟ = − ⎜⎝ ⎟ ⎜ ⎟ = ⎢ ⎜
T ⎟ −V ⎥
∂P ⎠ H ∂H ⎠ P ⎝ ∂P ⎠ T ⎜⎝ CP ⎟⎠ ⎣ ⎝ ∂T ⎠ P ⎦
⎛ ∂α ⎞ ∂H
General form: ⎜ Example: ⎛⎜ ⎞⎟
⎝ ∂β ⎟⎠ c ⎝ ∂A ⎠ P
⎛ ∂α ⎞ ⎛ ∂α ⎞ ⎛ ∂a ⎞
⎜⎝ ∂β ⎟⎠ = ⎜⎝ ⎟⎠ ⎜ ⎟⎠
c
∂a c ⎝ ∂ β c
⎛ ∂H ⎞ ⎛ ∂H ⎞ ⎛ ∂T ⎞ CP
Example: ⎜⎝ ⎟⎠ = ⎜⎝ ⎟⎠ ⎜⎝ ⎟⎠ = −
∂A P ∂T P ∂A P ⎛ ∂V ⎞
S +P ⎜
⎝ ∂T ⎟⎠ P
∂α ⎞ ∂H ⎞
General form: ⎛⎜ Example: ⎛⎜
⎝ ∂b ⎟⎠ γ ⎝ ∂P ⎟⎠ A
Procedure: Introduce a new observable and use Math Relation 5 to obtain one
Class 3 and two Class 2 derivatives:
⎛ ∂α ⎞ ⎛ ∂α ⎞ ⎛ ∂α ⎞ ⎛ ∂c ⎞
⎜⎝ ⎟ = ⎜ ⎟ + ⎜ ⎟ ⎜ ⎟
∂b ⎠ γ ⎝ ∂b ⎠ c ⎝ ∂c ⎠ b ⎝ ∂b ⎠ γ
Example:
⎛ ∂H ⎞ ⎛ ∂H ⎞ ⎛ ∂H ⎞ ⎛ ∂T ⎞
⎜⎝ ⎟⎠ = ⎜⎝ ⎟⎠ + ⎜⎝ ⎟ ⎜ ⎟
∂P A ∂P T ∂T ⎠ P ⎝ ∂P ⎠ A
⎛ ∂H ⎞ ⎛ ∂H ⎞ ⎛ ∂T ⎞ ⎛ ∂A ⎞
=⎜ −
⎝ ∂P ⎟⎠ T ⎜⎝ ∂T ⎟⎠ P ⎜⎝ ∂A ⎟⎠ P ⎜⎝ ∂P ⎟⎠ T
⎛ ∂V ⎞ ⎛ ∂T ⎞ ⎛ ∂A ⎞
=V −T ⎜ −
⎝ ∂T ⎟⎠ P P⎜
⎝ ∂A ⎟⎠ P ⎜⎝ ∂P ⎟⎠ T
C
⎧⎛ ∂A ⎞ ⎛ ∂V ⎞
⎪⎜⎝ ∂T ⎟⎠ = −S − P ⎜⎝ ⎟⎠
∂T
⎪⎪ P P
dA = −SdT − PdV ⇒ ⎨
⎪ ∂A
⎪⎛⎜ ⎞ ⎛ ∂T ⎞ ⎛ ∂V ⎞
⎟ = −S ⎜⎝ ⎟ − P ⎜⎝ ⎟
⎪⎩⎝ ∂P ⎠ T ∂P ⎠ T ∂P ⎠ T
Then
⎛ ∂H ⎞ ⎛ ∂V ⎞ ⎛ ∂T ⎞ ⎛ ∂A ⎞
⎜⎝ ⎟ =V −T ⎜⎝ ⎟⎠ − P⎜
⎝ ∂A ⎟⎠ P ⎜⎝ ∂P ⎟⎠ T
C
∂P ⎠ A ∂T P
⎛ ∂V ⎞
CP P ⎜
⎛ ∂V ⎞ ⎝ ∂P ⎟⎠ T
=V −T ⎜ ⎟ −
⎝ ∂T ⎠ P ⎛ ∂V ⎞
S +P ⎜ ⎟
⎝ ∂T ⎠ P
General form:
⎛ ∂α ⎞
Example: ⎛ ∂H ⎞
⎜⎝ ∂β ⎟⎠ ⎜⎝ ⎟
γ
∂G ⎠ A
⎛ ∂α ⎞ ⎛ ∂α ⎞ ⎛ ∂a ⎞
⎜⎝ ∂β ⎟⎠ = ⎜⎝ ⎟⎠ ⎜ ⎟⎠
γ
∂a γ ⎝ ∂ β γ
Example:
⎛ ∂H ⎞
⎜⎝ ⎟
⎛ ∂H ⎞ ⎛ ∂H ⎞ ⎛ ∂P ⎞ ∂P ⎠ A
⎜⎝ ⎟⎠ = ⎜⎝ ⎟⎠ ⎜⎝ ⎟⎠ = ∂G
∂G A ∂P A ∂G A ⎛ ⎞
⎜⎝ ⎟
∂P ⎠ A
⎛ ∂V ⎞
CP P ⎜
⎛ ∂H ⎞ ⎛ ∂V ⎞ ⎝ ∂P ⎟⎠ T
⎜⎝ ⎟⎠ =V −T ⎜⎝ ⎟⎠ −
∂P A ∂T P ⎛ ∂V ⎞
S +P ⎜
⎝ ∂T ⎟⎠ P
Also:
⎛ ∂G ⎞ ⎛ ∂G ⎞ ⎛ ∂G ⎞ ⎛ ∂T ⎞
⎜⎝ ⎟ = ⎜ ⎟ + ⎜ ⎟ ⎜ ⎟
∂P ⎠ A ⎝ ∂P ⎠ T ⎝ ∂T ⎠ P ⎝ ∂P ⎠ A
⎡ ⎛ ∂T ⎞ ⎛ ∂A ⎞ ⎤ ⎛ ∂T ⎞ ⎛ ∂A ⎞
=V + ( −S ) ⎢ − ⎜ ⎟ ⎜ ⎟ ⎥ =V + S ⎜⎝ ⎟ ⎜ ⎟
⎣ ⎝ ∂A ⎠ P ⎝ ∂P ⎠ T ⎦ ∂A ⎠ P ⎝ ∂P ⎠ T
⎛ ∂V ⎞
SP ⎜
⎝ ∂P ⎟⎠ T
=V +
⎛ ∂V ⎞
S +P ⎜
⎝ ∂T ⎟⎠ P
An alternative method for calculating physical property changes utilizes ideal gas
calculations and the differences between the real material and an ideal gas.
Ψresidual (T ,X ) = Ψ (T ,X ) − Ψig (T ,X )
Because the independent variable X may be either P or V, there actually are two
definitions of residual functions.
Ψr (T , ρ ) ≡ Ψ (T , ρ ) − Ψig (T , ρ )
where:
Ψr (T , ρ ) = density residual function
Ψ (T , ρ ) = value of Ψ for real material at T , ρ
Ψig (T , ρ ) = value of Ψ for an ideal gas at T , ρ
ΨR (T ,P ) ≡ Ψ (T ,P ) − Ψig (T ,P )
where:
ΨR (T ,P ) = pressure residual function
Ψ (T ,P ) = value of Ψ for real material at T ,P
Ψig (T ,P ) = value of Ψ for an ideal gas at T ,P
Note that the pressure and density forms of the residual functions are not
necessarily equal because they involve comparisons to different ideal gas states.
T ig = T
ρ ig = ρ
P ⎛ P ⎞
P ig = ρ ig RT ig = ρRT = ⎜⎝ Z =
Z ρRT ⎟⎠
The pressure form of the residual function compares the real fluid at T, P, r with
the ideal gas at
T ig = T
P ig = P
ig P ig P
ρ = = = ρZ
RT RT
Real Fluid T P r
Pressure Form T P rZ
The residual functions are related by the same expressions as the energy
functions:
Density Form
H r = U r + P rV
Ar = U r −TS r
G r = H r −TS r = U r + P rV −TS r
Pressure Form
H R = U R + PV R
AR = U R −TS R
GR = H R −TS R = U R + PV R −TS R
lim lim
ΨR (T ,P ) = 0 and Ψ r (T , ρ ) = 0
P →0 ρ→0
Then
ρ ρ
⌠⎡
⌠
⎮⎛ ∂Ψr ⎞ ⎮ ⎛ ∂Ψ ⎞ ⎛ ∂Ψig ⎞ ⎤
Ψr (T , ρ ) − Ψr (T ,0 ) = Ψr (T , ρ ) = ⎮⎜ d ρ = ⎮ ⎢⎜ ⎟⎠ − ⎜ ∂ ρ ⎟ ⎥ d ρ
⎮⎝ ∂ ρ ⎟⎠ ⎮⎢
⎮ ⎝ ∂ ρ T ⎝ ⎠ T ⎥⎦
⌡ T ⌡⎣
0 0
and
P P
⌠⎡
⌠
⎮⎛ ∂ΨR ⎞ ⎮ ⎛ ∂Ψ ⎞ ⎛ ∂Ψig ⎞ ⎤
ΨR (T ,P ) − ΨR (T ,0 ) = ΨR (T ,P ) = ⎮⎜ ⎟ dP = ⎮ ⎢⎜ ⎟ −⎜ ⎟ ⎥ dP
⎮⎝ ∂P ⎠ ⎝ ⎠
⎮ ⎢ ∂P T
⎮ ⎝ ∂P ⎠ T ⎥⎦
⌡ T ⌡⎣
0 0
2.6.4.3.3.1.1 Pressure
Dimensionless form:
P rV
= Z −1
RT
Calculate the residual internal energy at any density by integrating the partial
derivative along the isotherm from zero density (where the residual is known to
vanish) to the density of interest
ρ ρ ρ
⌠ ⌠
⎮ ⎛ ∂U ⎞ ⎮ ⎛ ∂U
⌠ ig ⎞
r
⎛ ∂ U ⎞
U r (T , ρ ) = ⎮⎜
⎮
⎟ d ρ = ⎮⎮ ⎜⎝ ∂ρ ⎟⎠ d ρ − ⎮⎮ ⎜ ∂ρ ⎟ d ρ
⎮ ⎝ ∂ρ ⎠T
⌡ ⌡ T ⌡⎝ ⎠T
0 0 0
MATH 4 gives:
⎛ ∂U ⎞ ⎛ ∂U ⎞ ⎛ ∂V ⎞ ⎛ ∂U ⎞ ⎛ 1⎞
⎜⎝ ∂ ρ ⎟⎠ = ⎜⎝ ⎟⎠ ⎜⎝ ∂ ρ ⎟⎠ = ⎜⎝ ⎟⎠ ⎜− 2⎟
T
∂V T T
∂V T ⎝ ρ ⎠
Use the physics summary statement for internal energy and MaxA to relate the
partial derivatives to observable properties
⎛ ∂U ⎞ ⎛ ∂S ⎞ ⎛ ∂V ⎞ ⎛ ∂P ⎞
dU = TdS − PdV ⇒ ⎜⎝ ⎟ = T ⎜⎝ ⎟ − P ⎜⎝ ⎟ = T ⎜⎝ ⎟ −P
∂V ⎠ T ∂V ⎠ T ∂V ⎠ T ∂T ⎠ V
ρ
⌠⎧
⎤ ⎫
ig
⎮ ⎪⎡ ⎛ ∂P ⎞ ⎤ ⎡ ⎛ ∂P ⎞ ⎪dρ
U r (T ,V ) = ⎮ ⎨ ⎢T ⎜⎝ ⎟ − P −
⎥ ⎢ ⎜⎝
T ⎟⎠ − P ⎥ ⎬ 2
⎮
⎮ ⎪⎣ ∂T ⎠ V ⎦ ⎣ ∂ T ⎦ ⎪ ρ
⌡⎩ ⎭
V
0
Then
⎛ ∂P ⎞ RT 2 ⎛ ∂ Z ⎞ RT 2 ⎛ ∂ Z ⎞
T⎜ −P = P + ⎟ −P =
⎝ ∂T ⎟⎠V ⎜
V ⎝ ∂T ⎠ V
⎜
V ⎝ ∂T ⎠ V
⎟
so that
⎛ ∂U ⎞ RT 2 ⎛ ∂ Z ⎞ ⎡ 1 ⎤ RT 2 ⎛ ∂ Z ⎞
⎜⎝ ∂ρ ⎟⎠ = V ⎜⎝ ∂T ⎟⎠ ⎢− 2 ⎥ = −
ρ
⎜⎝
∂
⎟⎠
T V ⎣ ρ ⎦ T ρ
ig
⎛ ∂Z ig ⎞
Z =1 ⇒ ⎜ ∂T ⎟ = 0
⎝ ⎠
⎛ ∂Z ⎞ ⎛ ∂ Z ⎞ ⎛ ∂ T1 ⎞ 1 ⎛ ∂Z ⎞
⎜⎝ ⎟ = =− 2⎜ 1⎟
∂T ⎠ ρ ⎜⎝ ∂ T1 ⎟⎠ ⎜⎝ ∂T ⎟⎠ T ⎝ ∂T ⎠ ρ
ρ ρ
ρ ρ
⌠
⎮ R ⎛ ∂Z ⎞
⌠ 2
⎮ RT ⎛ ∂Z ⎞
U r
(T , ρ ) = ⎮− ⎜⎝ ⎟⎠ d ρ = ⎮ ⎜ ⎟ dρ
⎮
⌡
ρ ∂T ρ ⎮ ρ ⎝∂ ⎠
1
⌡ T ρ
0 0
ρ ρ ρ
U r ⌠T ∂Z
⎛ ⎞ 1 ⌠⎮ ⎛ ∂ Z ⎞ dρ 1 ⌠⎮ ⎛ ∂ Z ⎞
= − ⎮⎮ ⎜
⎝ ⎟⎠ d ρ = ⎮⎜ 1 ⎟⎠ ρ = ⎮⎜ 1 ⎟⎠ d ln ρ
RT ⌡ ρ ∂T ρ T ⎮⌡ ⎝ ∂ T ρ
T ⎮⌡ ⎝ ∂ T ρ
0 0 0
NOTATION:
ρ ρ
⌠
⌠ ∂Z
⎮⎛ ⎞ d ρ 1 ⎮⎛ ∂Z ⎞ d ρ
Ι 1 = −T ⎮⎜ ⎟⎠ = ⎮⎜ 1 ⎟
⌡ ⎝ ∂T ρ ρ T ⎮⌡ ⎝ ∂ T ⎠ ρ ρ
0 0
r ig
⎛ ∂ Ar ⎞ ⎛ ∂A⎞ ⎛ ∂ Aig ⎞
A = A−A ⇒ ⎜ ∂ρ ⎟ = ⎜⎝ ∂ρ ⎟⎠ − ⎜ ∂ρ ⎟
⎝ ⎠T T ⎝ ⎠T
Use the physics summary statement for Helmholtz energy to relate the partial
derivatives to other properties
⎛ ∂A⎞
dA = −SdT − PdV ⇒ ⎜⎝ ⎟ = −P
∂V ⎠ T
⎛ ∂A ⎞ ⎛ ∂A ⎞ ⎛ ∂V ⎞ ⎛ 1⎞ P ZRT
⎜⎝ ∂ ρ ⎟⎠ = ⎜⎝ ⎟⎠ ⎜⎝ ∂ ρ ⎟⎠ = ( −P ) ⎜ − ⎟ = =
T
∂V T T ⎝ ρ2 ⎠ ρ2 ρ
ig
For the ideal gas, Z = 1, therefore
⎛ ∂Aig ⎞ RT
⎜ ∂ρ ⎟ =
⎝ ⎠T ρ
Then
ρ ρ ρ
⌠
⎮⎛ ∂A ⎞ ⌠⎛
r
ZRT RT ⎞ ⌠ dρ
A (T , ρ ) =
r
⎮⎜ ⎟ dρ =
⎮
⎮⎜
− d ρ = RT ⎮ ( Z − 1)
⎮ ⎝ ∂ρ
⌡ ⎠T ⌡⎝ ρ ρ ⎟⎠ ⎮
⌡ ρ
0 0 0
ρ
A (T , ρ )
r ⌠ dρ
= ⎮ ( Z − 1)
RT ⎮
⌡ ρ
0
NOTATION:
ρ
⌠⎛
Z − 1⎞
Ι2 = ⎮
⎮⎜
dρ
⌡⎝ ρ ⎟⎠
0
2.6.4.3.3.1.4 Enthalpy
H r U r P rV
r
H =U +P V r r
⇒ = + = Ι1 + ( Z − 1)
RT RT RT
CHEN 623 Lecture Notes (Date of last revision: 08/30/2021) 173
2.6.4.3.3.1.5 Entropy
r r U r − Ar rS r U r − Ar
r
A = U −TS ⇒ S = ⇒ = = Ι1 − Ι2
T R RT
2.6.4.3.3.1.6 Gibbs
Gr
r r
G = H −TS = A + P V r r r
⇒ = Ι2 + ( Z − 1)
RT
Ur Ar
= Ι1 = Ι2
RT RT
Hr Gr
= Ι1 + Z − 1 = Ι2 + Z − 1
RT RT
P rV Sr
= Z −1 = Ι1 − Ι 2
RT R
ρ ρ
⌠
1 ⎮⎡ ∂Z ⎤ d ρ ⌠Z −1
Ι1 = ⎮⎢ Ι2 =⎮ dρ
T ⎮ ⎣∂
⌡
(1/ T ) ⎥⎦ ρ ρ ⎮
⌡ ρ
0 0
Because nature forces the use of temperature and density as the independent
variables in the equation of state, the pressure forms of the residual functions are
calculated most easily using the density forms just derived and a calculation of
property changes in the ideal gas.
Methods to calculate Ψr are shown in the previous section, and the difference
term ⎡ Ψig (T ,P / Z , ρ ) − Ψig (T ,P, ρ Z ) ⎤ is an ideal gas calculation.
⎣ ⎦
ZRT RT RT
V R =V −V ig = − = ( Z − 1)
P P P
PV R
Dimensionless form: = ZR = Z −1
RT
Because entropy and volume (density) are the natural independent variables for
internal energy, the change in internal energy for the ideal gas is calculated most
easily from
ρ
⌠
⎛ ∂U ig ⎞
U ig
(T , ρ ) − U (T , ρZ ) =
ig ⎮
⎮ ⎜ ∂ρ ⎟ d ρ
⎮
⌡ ⎝ ⎠T
ρZ
where
⎛ ∂U ig ⎞ ⎛ ∂U ig ⎞ ⎛ ∂V ⎞ 1 ⎛ ∂U ig ⎞
⎜ ∂ ρ ⎟ = ⎜ ∂V ⎟ ⎜⎝ ∂ ρ ⎟⎠ = − 2 ⎜ ∂V ⎟
⎝ ⎠T ⎝ ⎠T T ρ ⎝ ⎠T
ig RT ⎛ ∂P ig ⎞ R
P = ⇒ ⎜ ⎟ =
V ⎝ ∂T ⎠V V
so that
⎛ ∂P ⎞ RT RT
T⎜ − = − =0
⎝ ∂T ⎟⎠V
P
V V
Therefore
ρ
⌠
⎛ ∂U ig ⎞
U ig
(T , ρ ) − U (T , ρZ ) =
ig ⎮
⎮ ⎜ ∂ρ ⎟ d ρ = 0
⎮
⌡ ⎝ ⎠T
ρZ
Density also is the most convenient variable for calculations involving the
Helmholtz energy
ρ
⌠
⎛ ∂A ig ⎞
A ig
(T , ρ ) − A (T , ρZ ) =
ig ⎮
⎮ ⎜ ∂ρ ⎟ d ρ
⎮
⌡ ⎝ ⎠T
ρZ
where
⎛ ∂A ig ⎞ ⎛ ∂A ig ⎞ ⎛ ∂V ⎞ 1 ⎛ ∂A ig ⎞
⎜ ∂ ρ ⎟ = ⎜ ∂V ⎟ ⎜⎝ ∂ ρ ⎟⎠ = − 2 ⎜ ∂V ⎟
⎝ ⎠T ⎝ ⎠T T ρ ⎝ ⎠T
and
⎛ ∂A ⎞
dA = −SdT − PdV ⇒ ⎜⎝ ⎟⎠ = −P
∂V T
ig RT ⎛ ∂Aig ⎞ RT ⎛ ∂Aig ⎞ RT
P = ⇒⎜ ⎟ = − ⇒ ⎜ ∂ρ ⎟ =
V ⎝ ∂V ⎠ T V ⎝ ⎠T ρ
so that
ρ
⌠
⎛ ∂A ig ⎞
A ig
(T , ρ ) − A (T , ρZ ) =
ig ⎮
⎮ ⎜ ∂ρ ⎟ d ρ
⎮
⌡ ⎝ ⎠T
ρZ
ρ
⌠ RT ρ ρ
= ⎮ d ρ = RT ln ρ ρZ = RT ln = −RT lnZ
⎮
⌡ ρ ρZ
ρZ
Then
2.6.4.3.3.2.4 Enthalpy
H R = U R + PV R = U r + P rV ⇒ H R = H r
2.6.4.3.3.2.5 Entropy
SR =
r r
U R − AR U − A − RT lnZ
=
( )
⇒ S R = S r + R lnZ
T T
G R = AR + P RV = Ar − RT lnZ + PV r ⇒ G R = G r − RT lnZ
PV R P rV
= = Z (T , ρ ) − 1
RT RT
UR Ur
= = Ι1 (T , ρ )
RT RT
HR Hr
= = Ι1 (T , ρ ) + Z (T , ρ ) − 1
RT RT
AR Ar
= − ln Z = Ι2 (T , ρ ) − ln Z (T , ρ )
RT RT
GR Gr
= − ln Z = Ι2 (T , ρ ) + ⎡⎣ Z (T , ρ ) − 1⎤⎦ − ln Z (T , ρ )
RT RT
SR Sr
= + ln Z = Ι1 (T , ρ ) − Ι2 (T , ρ ) + ln Z (T , ρ )
R R
When the residual function approach is used to calculate property changes, the
changes in ideal gas properties
Ψig (T2 , ρ2 ) − Ψig (T1, ρ1) or Ψig (T2 ,P2 ) − Ψig (T1,P1)
also are required. To obtain expressions for the ideal gas changes, we use the
general form
T2 X2
⌠ ⌠
⎛ ∂Ψig ⎞ ⎛ ∂Ψig ⎞
Ψ ig
(T2 ,X 2 ) − Ψ (T1,X 1) =
ig ⎮
⎮ ⎜ ∂T ⎟ dT +
⎮
⎮ ⎜ ∂X ⎟ dX
⎮
⌡ ⎝ ⎠X ⎮
⌡ ⎝ ⎠T
T1 X1
For the internal energy, derive the expression using density as the second
independent variable
T2 ρ2
⌠ ⌠
⎛ ∂U ig ⎞ ⎛ ∂U ig ⎞
U ig
(T2 , ρ2 ) − U (T1, ρ1) =
ig ⎮
⎮ ⎜ ∂T ⎟ dT +
⎮
⎮ ⎜ ∂ρ ⎟ d ρ
⎮
⌡
⎝ ⎠ρ ⎮
⌡ ⎝ ⎠T
T1 ρ1
For the internal energy of the ideal gas, the derivatives are
⎛ ∂U ig ⎞ ig
⎜ ∂T ⎟ = CV
⎝ ⎠ρ
⎛ ∂U ig ⎞ 1 ⎛ ∂U ig ⎞ 1 ⎡ ⎛ ∂P ig ⎞ ⎤
⎜ ∂ ρ ⎟ = − 2 ⎜ ∂V ⎟ = − 2 ⎢T ⎜ ∂T ⎟ − P ⎥ = 0
ig
⎝ ⎠T ρ ⎝ ⎠T ρ ⎢⎣ ⎝ ⎠V ⎥⎦
T1
T1
T2 T2
1 ⌠ CVig 1 ⌠⎮ ⎛ C Pig ⎞
Ι (T1,T2 ) = ⎮⎮
ig
dT = ⎮ ⎜ − 1⎟ dT
3
T2 ⌡ R T2 ⌡ ⎝ R
⎮ ⎠
T1 T1
Then
U ig ( 2 ) − U ig (1) ig
= Ι 3 (T1,T2 )
RT2
The change in the ideal gas enthalpy is calculated using the relation between the
enthalpy and the internal energy
H ig = U ig + P igV = U ig + RT
⎛ T1 ⎞
3 ( 1 2 ) ⎜ 1−
= Ιig T ,T +
⎝ T2 ⎟⎠
Alternate expression
T2
H ig
(T2 ) − H (T1) = ⌠⎮⌡ CPig dT
ig
T1
T2 ρ2
⌠ ⌠
⎛ ∂S ig ⎞ ⎛ ∂S ig ⎞
S ig
(T2 , ρ2 ) − S (T1, ρ1) =
ig ⎮
⎮ ⎜ ∂T ⎟ dT +
⎮
⎮ ⎜ ∂ρ ⎟ d ρ
⎮
⌡
⎝ ⎠ρ ⎮
⌡ ⎝ ⎠T
T1 ρ1
⎛ ∂S ig ⎞ CVig
⎜ ∂T ⎟ = T
⎝ ⎠ρ
⎛ ∂S ig ⎞ 1 ⎛ ∂S ig ⎞ 1 ⎛ ∂P ig ⎞ 1 R
⎜ ∂ρ ⎟ = − ⎜ ⎟ = − ⎜ ⎟ = − R ρ = −
⎝ ⎠T ρ 2 ⎝ ∂V ⎠ T ρ 2 ⎝ ∂T ⎠V ρ2 ρ
⎛ ∂S ig ⎞ CPig
⎜ ∂T ⎟ = T
⎝ ⎠P
⎛ ∂S ig ⎞ ⎛ ∂V ig ⎞ R
⎜ ∂P ⎟ = ⎜ ⎟ = −
⎝ ⎠ T ⎝ ∂T ⎠ P P
T2 P2 T2
⌠ C ig ⌠ R ⌠ CPig P
S ig
(T2 ,P2 ) − S (T1,P1) = ⎮
ig ⎮ P
T
dT − ⎮
⎮
⌡ P
dP = ⎮
⎮ T
dT − R ln 2
P1
⌡ ⌡
T1 P1 T1
T2
S ig
(T2 ,P2 ) − S (T1,P1) = ⌠⎮ CPig dT
ig
− ln
P2
R ⎮ R T P1
⌡
T1
T2 T2
⌠C ig
dT ⌠⎮ ⎛ CVig ⎞ dT
Ι ig
(T1,T2 ) = ⎮
⎮
P
= ⎮⎜
⎮⎝ R
+ 1⎟
⎠ T
4
⌡ R T ⌡
T1 T1
Then
S ig ( 2 ) − S ig (1) P2 ρ2T2
R
= Ιig
4 (T ,T
2 1 ) − ln
P1
= Ι ig
4 (T ,T
2 1 ) − ln
ρ1T1
Use the relation of the Helmholtz energy to the internal energy and entropy
Then
P2 ⎡ T1 ⎤ S ig (1)
= Ιig
3 ( )
T1,T2 − Ιig
4 (T1,T2 ) + ln P − ⎢1− T ⎥ R
1 ⎣ 2⎦
T2 ρ2 ⎡ T1 ⎤ S ig (1)
= Ιig
3 ( )
T1,T2 − Ιig
4 (T1,T2 ) + ln T ρ − ⎢1− T ⎥ R
1 1 ⎣ 2⎦
Use the relation of the Gibbs energy to the enthalpy and entropy
Then
G ig ( 2 ) − G ig (1) H ig ( 2 ) − H ig (1) S ig ( 2 ) T1S ig (1)
= − +
RT2 RT2 R T2R
⎡ T1 ⎤ ig P2 ⎡ T1 ⎤ S ig (1)
= Ιig
3 (T1,T2 ) + ⎢1− T ⎥ − Ι 4 (T1,T2 ) + ln P − ⎢1− T ⎥ R
⎣ 2⎦ 1 ⎣ 2⎦
P2 ⎡ T1 ⎤ ⎡ S ig (1) ⎤
= Ιig
3 ( )
T1,T2 − Ιig
4 (T1,T2 ) + ln P − ⎢1− T ⎥ ⎢1+ R ⎥
1 ⎣ 2 ⎦ ⎢⎣ ⎥⎦
T2 ρ2 ⎡ T1 ⎤ ⎡ S ig (1) ⎤
= Ιig
3 ( )
T1,T2 − Ιig
4 (T1,T2 ) + ln T ρ − ⎢1− T ⎥ ⎢1+ R ⎥
1 1 ⎣ 2 ⎦ ⎢⎣ ⎥⎦
Combine the previous results to obtain formulae for the general cases.
3 (T1,T2 ) ⋅RT2
+Ιig
Dimensionless form
U ( 2 ) − U (1) T
= Ι 1 ( 2 ) − 1 Ι 1 (1) + Ι 3ig (T1,T2 )
RT2 T2
Dimensionless form
H ( 2 ) − H (1) U ( 2 ) − U (1) T
= + Z (T2 , ρ2 ) − 1 Z (T1, ρ1)
RT2 RT2 T2
T1
= Ι1 (T2 , ρ2 ) − Ι1 (T1, ρ1)
T2
T1
3 (T1,T2 ) + Z (T2 , ρ2 ) −
+Ιig Z (T1, ρ1)
T2
H ( 2 ) − H (1) T
= Ι 1 ( 2 ) + Z ( 2 ) − 1 ⎡⎣Ι 1 (1) + Z (1) ⎤⎦ + Ι 3ig (T1,T2 )
RT2 T2
Dimensionless form
Rearrange to obtain
S ( 2 ) − S (1) ⎛ P2 Z (1) ⎞
= Ι1 ( 2 ) − Ι1 (1) − ⎡⎣Ι2 ( 2 ) − Ι2 (1) ⎤⎦ + Ι 4 (T1,T2 ) − ln ⎜
ig
R ⎝ P1Z ( 2 ) ⎟⎠
Dimensionless form
A ( 2 ) − A (1) U ( 2 ) − U (1) T2S ( 2 ) −T1S (1)
= −
RT2 RT2 RT2
⎧ T1 ⎫ ⎧
= ⎨Ι1 ( 2 ) − Ι1 (1) + Ι 3 (T1,T2 ) ⎬ − ⎨Ι1 ( 2 ) − Ι1 (1)
ig
⎩ T2 ⎭ ⎩
⎛ P2 Z (1) ⎞ ⎫
⎪ ⎧⎛ T1 ⎞ S (1) ⎫
− ⎡⎣Ι2 ( 2 ) − Ι2 (1) ⎤⎦ + Ιig (T1,T2 ) − ln ⎜ P Z 2 ⎟ ⎬ − ⎨⎜ 1− T ⎟ R ⎬
⎝ 1 ( )⎠ ⎪
4
⎭ ⎩⎝ 2⎠ ⎭
A ( 2 ) − A (1) ⎛ T ⎞⎡ S (1) ⎤
= ⎡⎣Ι2 ( 2 ) − Ι2 (1) ⎤⎦ + ⎜ 1− 1 ⎟ ⎢Ι1 (1) −
RT2 ⎝ T2 ⎠ ⎣ R ⎥⎦
⎛ P2 Z (1) ⎞
3 ( 1 2)
+Ιig 4 ( 1 2)
ig
T ,T − Ι T ,T + ln ⎜P Z 2 ⎟
⎝ 1 ( )⎠
G ( 2 ) − G (1) ⎛ T1 ⎞ ⎡ S (1) ⎤
= ⎜ 1− ⎟ ⎢Ι1 (1) − ⎥ + ⎡⎣Ι2 ( 2 ) − Ι2 (1) ⎤⎦
RT2 ⎝ T2 ⎠ ⎣ R ⎦
⎛ P2 Z (1) ⎞ T1
+Ιig
3 (T ,T
1 2 ) − Ι ig
4 (T ,T
1 2 ) + ln ⎜⎝ P Z ( 2 ) ⎟⎠ + Z ( 2 ) −
T2
Z (1)
1
Because the internal energy and entropy are defined by differential equations, an
arbitrary choice of one boundary condition is required to prepare a table.
Generally take
S† = 0 H† = 0 G† = 0 at (T †
,P † , ρ † )
The functions cannot all be chosen to be zero at the reference point because
where
(
Z † = Z T †, ρ † )
The table entries for the various properties are calculated according to
(
Ψ = Ψ − Ψ† + Ψ† )
CHEN 623 Lecture Notes (Date of last revision: 08/30/2021) 198
Begin with internal energy
U = U −U ( †
) +U †
⇒
U
RT
=
U −U† U†
RT
+
RT
U
RT
= Ι1 (T , ρ ) −
T † † ig †
Ι + Ι T ,T −
T 1 3
T† †
T
Z ( )
H U
= + Z (T , ρ )
RT RT
⎡P ⋅ Z † ⎤
S S − S† S†
R
=
R
+
R
† † ig †
= Ι1 − Ι1 − Ι2 + Ι2 + Ι 4 T ,T − ln ⎢ † ⎥
⎢⎣ P ⋅ Z ⎥⎦
( )
A U S
= −
RT RT R
G U S
= − +Z
RT RT R
where: † † †
(†
)
Ι1 = Ι1 (T , ρ ) ; Z = Z (T , ρ ); Ι1 = Ι1 T , ρ ; Z = Z T , ρ
† †
( )
CHEN 623 Lecture Notes (Date of last revision: 08/30/2021) 199
2.6.4.6.1 Ideal Gas Reference State
( )
If the reference state is an ideal gas at T † ,P † , then Ι1† = 0; Ι†2 = 0; Z † = 1, and
U
RT
ig †
= Ι1 + Ι 3 T ,T −
T†
T
( )
H
RT
ig †
= Ι1 + Z + Ι 3 T ,T −
T†
T
( )
S
R
= Ι1 − Ι2 + Ιig
4 T †
,T − lnZ(− ln
P
P†
)
A
RT
ig † ig †
(
= Ι2 + Ι 3 T ,T + Ι 4 T ,T −
TT
T
) P
− lnZ − ln †
P
( )
G
RT
ig † ig
(
†
= Ι2 + Ι 3 T ,T − Ι 4 T ,T + Z −
T†
T
) ( P
)
− lnZ − ln †
P
0.01 °C PHvapO
2
60 °F 14.696 psia
60 °F 14.73 psia
15 °C 1 bar
0 °C 1 atm
0 °F 1 atm
25 °C 1 atm
vap
-40 °C P
0K 1 bar
In each case
(
Ψ (T ,P ) = Ψ † + Ψ − Ψ † )
ΨΔ (T ,P ) = Ψ * + ( Ψ Δ
−Ψ* )
The differences (Y-Y†) and (Y∆-Y*) are calculated using the procedures
developed previously. Subtracting the first equation from the second yields
Ψ Δ (T ,P ) = Ψ (T ,P ) + Ψ * −Ψ † + Ψ Δ − Ψ * − Ψ − Ψ † ( ) ( )
Table Property Property
Value Change Change
U (T ,P ) U †
= +
U −U† ( )
RT RT RT
(T ,P ) = U * + (U )
Δ
U Δ −U *
RT RT RT
U Δ − U U * −U † ⎡U Δ − U * ⎤ ⎡U − U † ⎤
= +⎢ ⎥−⎢ ⎥
RT RT ⎣⎢ RT ⎥⎦ ⎢⎣ RT ⎥⎦
where
U −U†
RT
= Ι1 (T , ρ ) −
T†
T
Ι1 T † , ρ † + Ιig
3 T †
,T ( ) ( )
U Δ −U * T*
= Ι1 (T , ρ ) − Ι1 (T *, ρ * ) + Ιig
3 (T *,T )
RT T
⎡U Δ − U * ⎤ ⎡U − U † ⎤ T*
⎢ ⎥−⎢ ⎥ = Ι1 (T , ρ ) − Ι1 (T *, ρ * )
⎢⎣ RT ⎥⎦ ⎢⎣ RT ⎥⎦ T
⎡ ⎤
+Ιig
3
T†
( ) ( )
(T *,T ) − ⎢Ι1(T , ρ ) − Ι1 T †, ρ † + Ιig3 T †,T ⎥
⎢⎣ T ⎥⎦
=
T†
T
(
Ι1 T † , ρ † −
T*
T
)
Ι1 (T *, ρ * ) + Ιig
3 (T *,T ) − Ι ig
3 T †
,T ( )
=
T†
T
(
Ι1 T † , ρ † −
T*
T
)
Ι1 (T *, ρ * ) + Ιig
3 T *,T †
( )
For the most common choice of standard state values
H† = 0 H* = 0
S† = 0 S* = 0
G† = 0 G* = 0
U† = -Z†RT† U* = -Z*RT*
A† = -Z†RT† U* = -Z*RT*
U * −U † ( −Z *RT *) − ( −Z†RT† ) T† † T *
= = Z − Z*
RT RT T T
so that
U Δ − U U * −U † ⎡U Δ − U * ⎤ ⎡U − U † ⎤
= +⎢ ⎥−⎢ ⎥
RT RT ⎢⎣ RT ⎥⎦ ⎢⎣ RT ⎥⎦
=
T† † T *
T
Z −
T
Z *+
T † † T * * ig
T
Ι1 −
T
Ι1 + Ι 3 T *,T † ( )
U Δ − U U * −U † U
= +
R
(T †
,P † )−U R
(T *,P * ) + U (T *,P * ) − U (T ,P )
ig ig † †
RT RT RT RT RT
S (T ,P ) S † ⎡ S − S † ⎤
= +⎢ ⎥
R R ⎢⎣ R ⎥⎦
=
S † S R (T ,P ) S
+ −
R
(T †
,P † )+S ig
(T ,P ) − S ig (T †,P † )
R R R R
S Δ (T ,P ) S * ⎡ S Δ − S * ⎤
= +⎢ ⎥
R R ⎢⎣ R ⎥⎦
=
S
+
*
S (T ,P )R
−
(
S R T * ,P * )+S ig
(T ,P ) − S ig (T * ,P * )
R R R R
S Δ − S S* − S† S
= +
R
(T †
,P † ) − S (T ,P ) + S (T
R * * ig †
) (
,P † − S ig T * ,P * )
R R R R R
where
S* − S†
=0
R
(
S R T † ,P † ) = Ι (T †
) (
, ρ † − Ι2 T † , ρ † + lnZ † )
1
R
(
S R T * ,P * ) = Ι (T , ρ ) − Ι (T , ρ ) + lnZ
* * * * *
1 2
R
( )
S ig T † ,P † − S ig T * ,P * ( ) = Ι (T ,T ) − ln⎛ P
ig * †
†⎞
4 ⎜ * ⎟
R ⎝P ⎠
H (T ,P ) U (T ,P )
= + Z (T ,P )
RT RT
H Δ (T ,P ) U Δ (T ,P )
= + Z (T ,P )
RT RT
Then
H Δ − H U Δ −U
=
RT RT
Begin with the relation between the Helmholtz energy, the internal energy and
the entropy:
A (T ,P ) U (T ,P ) S
= −
RT RT R
A Δ (T ,P ) U Δ (T ,P ) S Δ
= −
RT RT R
Then
A Δ − A ⎡U Δ − U ⎤ ⎡ S Δ − S ⎤
=⎢ ⎥−⎢ ⎥
RT ⎢⎣ RT ⎥⎦ ⎢⎣ R ⎥⎦
Begin with the relation between the Helmholtz energy, the internal energy and
the entropy:
G (T ,P ) U (T ,P ) S
= +Z −
RT RT R
G Δ (T ,P ) U Δ (T ,P ) SΔ
= +Z −
RT RT R
Then
G Δ − G ⎡U Δ − U ⎤ ⎡ S Δ − S ⎤
=⎢ ⎥−⎢ ⎥
RT ⎢⎣ RT ⎥⎦ ⎢⎣ R ⎥⎦
The adjustment from one reference state to another also can be calculated as
follows. The values in each table can be expressed as
(T ,P,ρ )
Ψ = Ψ† + ⌠
⎮
⌡ dΨ
†
and
(T ,P,ρ )
Δ
Ψ = Ψ*+ ⌠
⎮
⌡ dΨ
*
( )+
†
Δ * †
Ψ =Ψ+ Ψ −Ψ ⌠dΨ
⎮
⌡
*
(dS ) total
U total ,V total
≥0
Over a long period of time, any process that may occur eventually will occur.
The equilibrium state (the state from which the isolated system will not depart no
matter how long we wait) therefore is characterized by a state for which the
change of entropy in moving to any adjacent state is
(dS ) total
U total ,V total
≤0
Any such change would violate the second law of thermodynamics, therefore the
system remains in the equilibrium state.
At fixed internal energy and volume, entropy is at a maximum for the
equilibrium state.
The equilibrium state is the one state from which the system does not move for
specified (fixed) values of Stotal and Vtotal.
First Law:
total
dUsystem / + dW
= dQ /
Therefore,
total
dUsystem /
= dQ
so that, for any event allowed by the Second Law, (dU total
system )
S total ,V total
≤0
At fixed entropy and volume, the internal energy is at a minimum for the
equilibrium state.
/ = 0 ) must have
at constant volume ( dW
total
dAsystem / −TdSsystem
= dQ total
≤0
/ + dW
= dQ / + PdVsystem
total
/ = −PdVsystem
dW total
Therefore events that may occur at constant pressure and entropy must have
total
dHsystem / ≤0
= dQ
(dH total
system ) total
Ssystem ,P
≥0
for any change to an adjacent state, the isolated system cannot leave that state.
/ + dW
= dQ / + PdVsystem
total total
−TdSsystem
/ = −PdVsystem
dW total
/ ≤ TdSsystem .
The Second Law requires that dQ
total
(dG total
system )T ,P
≤0
(dG total
system )T ,P
≥0
for any change to an adjacent state, system cannot leave that state.
(dS ) total
U total ,V total
≤0
(dU total
system ) S total ,V total
≥0
(dH total
system ) total
Ssystem ,P
≥0
(dA total
system ) total
T ,Vsystem
≥0
(dG total
system )T ,P
≥0
Metastable
Stable Equilibrium Equilibrium
Consider a closed system with no external fields present, so that only PdV work
is involved in the statement of the First Law. As shown below, the system
contains two chambers, A and B, separated by a frictionless, thermally
conductive, semipermeable piston. The system is isolated from its surroundings.
A B
total
dU A = TA dS Atotal − PA dVAtotal + ∑ µiAdn Ai
i
1 PB µiB B
dSBtotal = total
dUB + total
dVB −∑ dn i
TB TB T
i B
{ }
Also, U Atotal . VAtotal and niA are independent variables, therefore the equilibrium
condition is satisfied in general only if all terms in square brackets are zero.
1
−
1
=0 ⇒ TA = TB
TA TB
1 1
− =0 ⇒ PA = PB
PA PB
1 1
− =0 ⇒ µiA = µiB
µ Ai µ Bi
Note that no external fields are present and there is no mechanism for
doing any kind of work other than PdV work.
(dU ) total
S total
,V total
≥ 0 and (d U )
2 total
S total ,V total
≥0
The Physics Summary Statement for total internal energy has total entropy and
total volume as the natural independent variables.
⎛ ∂U ⎞ ⎛ ∂U ⎞
dU = ⎜ +
⎝ ∂S ⎟⎠V ⎜⎝ ⎟ dV
dS
∂V ⎠ S
⎛ ∂ 2U ⎞ ⎛ ∂ 2U ⎞ ⎛ ∂ 2U ⎞
d U = ⎜ 2 ⎟ (dS ) + 2 ⎜ ( ) ( ) ( )
2 2 2
⎟ dS dV + ⎜ 2⎟ dV
⎝ ∂S ⎠V ⎝ ∂S ∂V ⎠ ⎝ ∂V ⎠ S
⎛ ∂T ⎞ ⎛ ∂T ⎞
T = T (S,V ) ⇒ dT = ⎜ +
⎝ ∂S ⎟⎠V ⎜⎝ ⎟ dV
dS
∂V ⎠ S
⎛ ∂U ⎞
T =⎜
⎝ ∂S ⎟⎠V
so that
⎛ ∂ 2U ⎞ ⎛ ∂ 2U ⎞
dT = ⎜ 2 ⎟ dS + ⎜ ⎟ dV
⎝ ∂S ⎠V ⎝ ∂V ∂S ⎠
Rearrange to obtain
−2 −1 −2 2
⎛ ∂ 2U ⎞ ⎛ ∂ 2U ⎞ ⎛ ∂ 2U ⎞ ⎛ ∂ 2U ⎞ ⎛ ∂ 2U ⎞
(dS )2 =⎜ 2⎟ (dT ) − 2 ⎜ 2 ⎟
2
⎜ ∂V ∂S ⎟ (dS ) (dV ) − ⎜ 2⎟ ⎜ ∂V ∂S ⎟ (dV )2
⎝ ∂S ⎠ ⎝ ∂S ⎠ ⎝ ⎠ ⎝ ∂S ⎠ ⎝ ⎠
⎛ 2 ⎞ ⎧⎛ 2 ⎞ −2 ⎛ 2 ⎞ −2 ⎛ 2 ⎞
2
∂ U ⎪ ∂ U ∂ U ∂ U
d 2U = ⎜ 2 ⎟ ⎨⎜ 2 ⎟ (dT ) − ⎜ 2 ⎟ ⎜ ( )
2 2
⎟ dV
⎝ ∂S ⎠ ⎪⎝ ∂S ⎠ ⎝ ∂S ⎠ ⎝ ∂V ∂S ⎠
⎩
−1 ⎫
⎛ ∂ 2U ⎞ ⎛ ∂ 2U ⎞ ⎪ ⎛ ∂ 2U ⎞ ⎛ ∂ 2U ⎞
−2 ⎜ 2 ⎟ ⎜ ⎟ ( dS ) ( dV ) ⎬ + 2 ⎜ ∂S ∂V ⎟ ( dS ) ( dV ) + ⎜ 2⎟ ( dV )2
⎝ ∂S ⎠ ⎝ ∂V ∂S ⎠ ⎪⎭ ⎝ ⎠ ⎝ ∂V ⎠
−1 ⎡⎛ 2 ⎞ ⎛ 2 ⎞ −1 ⎛ 2 ⎞ 2 ⎤
⎛ ∂ 2U ⎞ ∂ U ∂ U ∂ U
=⎜ 2⎟ (dT )2 + ⎢⎜ 2 ⎟ − ⎜ 2 ⎟ ⎜ ⎥ (dV )2
⎝ ∂S ⎠ ⎢⎝ ∂V ⎠ ⎝ ∂S ⎠ ⎝ ∂V ∂S ⎟⎠ ⎥
⎣ ⎦
⎛ ∂ 2U ⎞
−1 ⎡⎛ 2 ⎞ ⎛ 2 ⎞ −1 ⎛ 2 ⎞ 2 ⎤
∂ U ∂ U ∂ U
⎜ 2⎟
> 0 and ⎢⎜ ⎟ − ⎜ ⎟ ⎜ ⎟
⎥>0
⎝ ∂S ⎠ ⎢⎝ ∂V ⎠ ⎝ ∂S ⎠ ⎝ ∂V ∂S ⎠ ⎥
2 2
⎣ ⎦
−1 −1
⎛ ∂ 2U ⎞ ⎛ ∂T ⎞ ⎛ ∂S ⎞ CV
=⎜ = =
⎜ 2⎟
⎝ ∂S ⎠ ⎝ ∂S ⎟⎠V ⎜⎝ ⎟
∂T ⎠V T
⎛ ∂ 2U ⎞ ⎛ ∂P ⎞
⎜ 2⎟
= − ⎜⎝ ⎟
⎝ ∂V ⎠ ∂V ⎠ S
⎛ ∂ 2U ⎞ ⎛ ∂ 2U ⎞ ⎛ ∂P ⎞
=
⎜ ∂S ∂V ⎟ ⎜ ∂V ∂S ⎟ = − ⎜⎝ ⎟
⎝ ⎠ ⎝ ⎠ ∂S ⎠V
By Math 4
⎛ ∂P ⎞ ⎛ ∂P ⎞ ⎛ ∂P ⎞
= −⎜ −
⎝ ∂V ⎟⎠ s ⎜⎝ ∂T ⎟⎠V ⎜⎝ ∂S ⎟⎠V
By MaxA
⎛ ∂P ⎞ ⎛ ∂S ⎞ ⎛ ∂P ⎞
= −⎜ −
⎝ ∂V ⎟⎠ s ⎜⎝ ∂V ⎟⎠ T ⎜⎝ ∂S ⎟⎠V
and, by Math 5
⎛ ∂P ⎞
= −⎜
⎝ ∂V ⎟⎠ T
−1
⎛ ∂ 2U ⎞ CV
⎜ 2⎟
>0 ⇒ > 0, CV > 0
⎝ ∂S ⎠ T
Variable Pα Tα ρα x αi Total
P α = P α +1 α = 1,2,,π − 1 π −1
T α = T α +1 α = 1,2,,π − 1 π −1
⎧α = 1,2,,π − 1
µiα = µiα +1 ⎨ N (π − 1)
⎩i = 1,2,,N
Equations of State
(
F T α , ρα ,P α = 0 ) α = 1,2,,π p
Mass Balances
N
∑ x iα = 1 α = 1,2,,π p
i =1
Reactions
equilibrium
Km m = 1,2,,r r
Duhem's Theorem
For any closed system formed initially from given masses of prescribed chemical
species, the equilibrium state is determined completely when any two
independent variables are specified.
Enthalpy also is useful for the mass transfer term in the energy balance equation
because it automatically includes the injection work.
Therefore
2
W ≥ H 2 − H1 − ∫T dS
1
A2 − A1 = U 2 − U1 − (T2S2 −T1S1)
Then
A2 − A1 = Q +W −T2S2 +T1S1
Then
A2 − A1 = Q +W −T (S2 − S1)
W = A2 − A1 +T (S2 − S1) − Q ⇒ W ≥ A2 − A1
Wmax = A2 − A1
Then
W '− P (V2 −V1) ≥ (G2 − G1) − P (V2 −V1) ⇒ W ' ≥ (G2 − G1)
Then, the maximum useful work possible with a constant temperature and
pressure process is
W 'max = G2 − G1 = ΔG
Pα = P β
T α =T β
µiα = µiβ
⎛ ∂P ⎞
κT > 0; ⎜ <0
⎝ ∂V ⎟⎠T
The equation of state contains the information about interactions between the
molecules.
The ideal gas heat capacity contains information about the energy levels within
the molecules.
Many references provide coefficients for one of the following forms for the ideal
gas heat capacity:
CPig = A + BT + CT 2 + DT 3
CPig = a + bT + cT 2 + dT −2
The ideal gas heat capacity is determined completely by the characteristics of the
individual molecules.
Each degree of freedom for translation contributes 21 R to the ideal gas heat
capacity at constant volume.
Each degree of freedom for rotation contributes 21 R to the ideal gas heat
capacity at constant volume.
NOTES
CVig = 23 R
CPig = CVig + R = 23 R + R = 52 R
Diatomic molecules:
Translation: three degrees of freedom:
Rotation: two degrees of freedom:
CVig = 23 R + R = 52 R
CPig = CVig + R = 52 R + R = 72 R
Ideal Gas
Need to develop equations of state that describe all the behavior shown above.
CHEN 623 Lecture Notes (Date of last revision: 08/30/2021) 249
3.2.2 Brief History of Equations of State
1873: van der Waals: Dissertation on the continuity of the gas and liquid
states. The equation of state and the concept of corresponding states.
1970's: Soave, Peng & Robinson: cubic equations optimized for phase
equilibrium calculations
PV
Z ≡ ⇔ PV = ZRT
RT
In general, the density form is the more useful of the two. The pressure form is
useful at low pressures (< 5 bar) or at very high reduced temperatures (>2).
Cubic forms of the equation of state are very popular for practical calculations,
because the cubic form is the simplest volume dependence that is capable of
describing vapor-liquid phase equilibrium
The first cubic form for an equation of state was proposed by van der Waals in
his Ph.D. dissertation (1873). The parameter a is nonzero because of attractions
between molecules, and the parameter b is nonzero because each molecule
occupies a finite volume.
Original form:
⎛ a ⎞
⎜⎝ P + 2⎟
⎠
(V − b ) = RT
V
Alternative forms:
RT a ⎛ RT ⎞ 2 a ab
P= − 2; V 3 − ⎜b + ⎟ V + V − =0
V −b V ⎝ P ⎠ P P
V a 1 aρ
Z vdw = − = −
V − b RTV 1− b ρ RT
The equation of state also may be written with the compressibility factor as the
variable:
Z 3 − (1+ A ) Z 2 + BZ − AB = 0
where
bP
A=
RT
aP
B=
(RT )2
3.2.3.3.2.1 Martin
RT α (T ) − δ (T ) / V
P= −
V (V + β ) (V + γ )
3.2.3.3.2.2 Abbott
RT θ (V − η )
P= −
V − b (V − b ) V 2 + δV + ε ( )
Ref.: O. Redlich & J. N. S. Kwong, Chem. Review, 44, 233 - 244 (1949)
RT a
P= −
V − b V (V + b ) T
where
R 2Tc2.5 ⎛ 1 ⎞
a = Ωa ⎜ Ωa = = 0.427480⎟
Pc ⎜⎝ 9 ( 3
)
2 −1 ⎟⎠
RTc ⎛ 3
2 −1 ⎞
b = Ωb ⎜ Ωb = = 0.086640⎟
Pc ⎝ 3 ⎠
and
V a 1 aρ
Z RK = − = −
V − b RT 1.5 (V + b ) 1− b ρ RT 1.5 (1+ b ρ )
RT a
P= −
V − b (V + c )2
RT aα
P= −
V − b V (V + b )
or
V aα 1 aαρ
Z SRK = − = −
V − b RT (V + b ) 1− b ρ RT (1+ b ρ )
where
a = 0.42747R 2Tc2 / Pc
b = 0.08664RTc / Pc
( )( )
2
α = ⎡1+ 0.480 + 1.574ω − 0.176ω 2 1− Tr ⎤
⎣ ⎦
(
α = ⎡⎢1+ 0.48508 + 1.55171ω − 0.15613ω 2 1− Tr
⎣
)( )
2⎤
⎥⎦
Ref.: D. Y. Peng & D. B. Robinson, Ind. & Eng. Chem., Fundamentals, 15,
59-64 (1976)
RT a (T )
P= −
V − b V (V + b ) + b (V − b )
or
V a (T )V 1 a (T ) ρ
Z PR = − = −
V − b RT ⎡⎣V (V + b ) + b (V − b ) ⎤⎦ 1− b ρ RT ⎡1+ 2b ρ − (b ρ )2 ⎤
⎣ ⎦
where
a (Tc ) = 0.45724R 2Tc2 / Pc
b (Tc ) = 0.07780RTc / Pc
a (T ) = a (Tc )α (Tr ,ω )
( )( )
2
α (Tr ,ω ) = [1+ 0.37464 + 1.54226ω − 0.26992ω 2
1− Tr ⎤
⎦
P=
RT ⎛ C ⎞ 1 1 aα
+ ⎜ B0RT − A0 − 02 ⎟ 2 + (bRT − a ) 3 + 6 +
(
c 1+ γ / V 2 ) exp ⎛ − γ
⎞
⎝ ⎜⎝ ⎟
V T ⎠V V V T 2V 3 V2⎠
Parameters: A0 , B0 , C0 , a, b, c, α , γ
RT ⎛ C D E ⎞ 1 ⎛ α⎞ 1
P= + ⎜ B0RT − A0 − 02 + 03 − 04 ⎟ 2 + ⎜ bRT − a − ⎟ 3
V ⎝ T T T ⎠V ⎝ T ⎠V
⎛ α ⎞ 1
+α ⎜ a + ⎟ 6 +
c (
1+ γ / V 2
) ⎛ γ ⎞
⎝
exp ⎜⎝ − 2 ⎟⎠
T ⎠V 2 3
T V V
Parameters: A0 , B0 , C0 , D0 , E 0 , a, b, c, α , γ
P=
RT ⎛ C0 ⎞ 1 1 aα c 1+ γ / V
+ ⎜ B0RT − A0 − 2 ⎟ 2 + (bRT − a ) 3 + 6 +
2
⎛ γ ⎞ ( )
⎝ ⎠ 2 3
exp ⎜⎝ − 2 ⎟⎠
V T V V V T V V
Parameters: A0 , B0 , C0 , a, b, c, α , γ
where
⎡
⎢
⎛ bρ ⎞ β − 1 α + β +1
Abase ( ρ,T ) = RT ⎢ − ln ⎜ 1− ⎟ − +
⎢ ⎝ 4 ⎠ 1− b ρ ⎛ bρ ⎞
2
⎢ 4 2 ⎜⎝ 1− 4 ⎟⎠
⎣
⎛B ⎞ α − β +γ ρRT ⎤
+b ρ ⎜ − γ ⎟ − + ln ⎥
⎝b ⎠ 2 P0 ⎦
k
T ⎛T ⎞
b = b1 ln + ∑ bk ⎜ 0 ⎟
T0 k =0,1,3,5 ⎝ T ⎠
k
⎛T ⎞
B = ∑ Bk ⎜ 0 ⎟
k =0,1,2,4 ⎝T ⎠
CHEN 623 Lecture Notes (Date of last revision: 08/30/2021) 262
(i )
g ⎛T ⎞ − ρ k (i )
(1− e ) (i )
( (i )
)
36
Aresidual ( ρ,T ) = ∑ i ⎜ ⎟ + ∑ gi δ i exp −α i δ i − βi τ i2
i =1 k (i ) ⎝ T0 ⎠
ρ − ρi
δi =
ρi
T −Ti
τi =
Ti
⎡ ⎛ C1 ⎞ 18
i −6 ⎤
Aideal gas (T , ρ ) = −RT ⎢ ⎜
1+ + C 2⎟ lnTR + ∑ C i R ⎥
T
⎣ ⎝ TR ⎠ i =3 ⎦
T
TR =
100 K
Ref.: R. Tilner-Roth and A. Yokozeki, J. Phys. Chem. Ref. Data, 26, 1273-
1328 (1997).
A Aig Ar
= + = Φig + Φr
RT RT RT
ρ T
where, in terms of the reduced variables, δ = ; τ= c
ρc T
( )
8 19
Φ = ∑ ai δ τ + ∑ ai δ d i τ t i exp −δ ei
r di ti
i =1 i =9
There are seven terms in the ideal gas formulation and nineteen terms in the
residual formulation.
Ref.: W. Wagner and A. Pruß, J. Phys. Chem. Ref. Data, 31, 375-535 (2002).
Using forms and the reduced variables defined on the previous page, the ideal
gas and residual terms are given by
8
ig
Φ = lnδ + n1o + n2oτ + n3o lnτ + ∑ nio ln ⎡⎣1− exp ( −γ i τ ) ⎤⎦
i =4
( )
7 51
Φ = ∑ ni δ τ + ∑ ni δ d i τ t i exp −δ c i
r di ti
i =1 i =8
54 56
+ ∑ ni δ τ exp ⎡ −α i (δ − ε i ) − βi (τ − γ i ) ⎤ + ∑ ni Δbi δψ
di ti 2 2
i =52 ⎣ ⎦ i =55
where
ai
Δ = θ + Bi ⎡(δ − 1) ⎤
2 2
⎣ ⎦
1/ ( 2 βi )
θ = (1− τ ) + Ai ⎡(δ − 1) ⎤
2
⎣ ⎦
ψ = exp ⎡ −Ci (δ − 1) − Di (τ − 1) ⎤
2 2
⎣ ⎦
This approach uses the same form for the ideal gas contribution as those used
previously in the high accuracy Helmholtz equations, but the residual form is
replaced by a ratio of polynomials in density which have temperature dependent
coefficients.
This form is slightly less accurate near the critical temperature and pressure than
the previous form, but, because there are no exponential terms in density, the
computation time required to solve for denisties is reduced by more than a factor
of 20. This form does require about 15 times the computational time as cubic
equations, but the accuracies are much better.
⎛ ∂A ⎞
P = ρ2 ⎜
⎝ ∂ ρ ⎟⎠ T
⎛ ∂A ⎞
S = −⎜
⎝ ∂T ⎟⎠ ρ
U = A +TS
P
H =U +
ρ
P
G = A+
ρ
3.2.4.1 2 parameter
⎛ ∂P ⎞ ⎛ ∂2 P ⎞
⎜⎝ ⎟⎠ = 0 and ⎜ 2⎟
=0
∂V T c
⎝ ∂V ⎠T
c
where
P V T
Pr = ; Vr = ; Tr =
Pc Vc Tc
Ref.: K. S. Pitzer, R. F. Curl, et al., Ind. & Eng. Chem., 50, 265-274 (1958).
Z = Z 0 + ωZ1
where
⎛ ⎞
ω ≡ − log10 ⎜ Prvap ⎟ −1
⎝ Tr =0.7 ⎠
Because they can be derived directly from Z and its first derivative with respect to
temperature, the energies (U, H, A, G) and the entropy (S) also may be
represented using the same form.
Z = Z (0) +
ω
ω (r )
Z ( )
(r ) − Z (0 ) = Z (0 ) + ω Z (1)
ω (r ) = 0.3978
PrVr B C D c4 ⎛ γ ⎞ ⎛ γ ⎞
Z = = 1+ + + + ⎜ β + 2 ⎟ exp ⎜ − 2 ⎟
Tr Vr Vr2 Vr3 Tr3Vr2 ⎝ Vr ⎠ ⎝ Vr ⎠
B = b1 − b2 / Tr − b3 / Tr2 − b4 / Tr3
C = c1 − c 2 / Tr + c 3 / Tr3
D = d1 + d 2 / Tr
For polar fluids, the dipole moment of the molecule (µ) often is used as a fourth
parameter in equations of state.
⎛ V T ⎞
Z i (V ,T ) = Z o ⎜ , ⎟
⎝ hi ,o fi ,o ⎠
r ⎛ V T ⎞
Air (V ,T ) = fi ,o Ao ⎜ , ⎟
⎝ i ,o fi ,o ⎠
h
Toc Voc
fi ,o = c
θ (Tr ,Vr ) and hi ,o = c
φ (Tr ,Vr )
Ti Vi
Here the functions θ (Tr ,Vr ) and φ (Tr ,Vr ) are shape factors.
The subscripts i and o refer to the fluid of interest and the reference fluid
respectively.
CHEN 623 Lecture Notes (Date of last revision: 08/30/2021) 273
3.2.4.4.1 Modified BWRS (NIST form)
(
P = ρRT + ρ 2 N1T + N 2T 0.5 + N 3 + N 4T −1 + N 5T −2 )
( ) (
+ ρ 3 N 6T + N 7 + N 8T −1 + N 9T −2 + ρ 4 N10T + N11 + N12T −1 + ρ 5N13 )
+ ρ6 (N T + N T
14
−1
15
−2
) + ρ (N
7
16T ) + ρ (N
−1 8
17T
−1
+ N18T −2)
+ ρ9 (N T )
19
−2
+ρ3 (N T + N T
20
−2
21
−3
) exp (−γρ )2
+ ρ5 (N T + N T
22
−2
23
−4
) exp (−γρ )2
+ ρ7 (N T + N T
24
−2
25
−3
) exp (−γρ )2
+ ρ9 (N T + N T
26
−2
27
−4
) exp (−γρ )2
+ ρ11 (N T + N T
28
−2
29
−3
) exp (−γρ )
2
+ ρ13 (N T + N T
30
−2
31
−3
+ N T ) exp ( −γρ )
32
−4 2
All forces of relevance here, e.g. interactions among electrons and nuclei, are
electrostatic in nature.
∂Ψ
HΨ = i
∂t
where
Ψ (R,r,t ) is the wave function
⎧ ⎫
H = K + U = ⎨− ∑ ∇i2 + U ⎬ is the Hamiltonian
⎩ 2m i ⎭
Time-independent form
HΨ = EΨ
CHEN 623 Lecture Notes (Date of last revision: 08/30/2021) 275
Born-Oppenheimer approximation: Electronic motion is very fast compared to
nuclear motion so that an average can be taken over the electrons for a given set
of nuclear coordinates.
The potential energy surface is everything you need to know to describe the
motion of the nuclei in the system.
In practice:
• except in the simplest cases, the Schrödinger equation cannot be solved
exactly
• even with accurate approximations, the computational cost of solving it is
high
Atoms (or groups of atoms) are the fundamental units. Their nuclei (centers)
become the focal points for the empirical mathematical expressions.
3.3.3.1 Bond-stretching
Example equation
Ubond = k (r − r0 )
2
3.3.3.2 Angle-bending
Example equation
Uang = k (θ − θ 0 )
2
Example equation
3
Utor = ∑ k i ⎡⎣1− cos (φ − φ0i ) ⎤⎦
i =1
Generally used for a set of three atoms connected to an sp2 hybridized central
atom (amide N, sp2 C).
Example equation
Uoop = k χ 2
Evalence = ∑ ⎡k 2 (b − bo ) + k 3 (b − bo ) + k 4 (b − bo ) ⎤
2 3 4
b ⎣ ⎦
+ ∑ ⎡k 2 (θ − θo ) + k 3 (θ − θo ) + k 4 (θ − θo ) ⎤
2 3 4
θ ⎣ ⎦
+ ∑ ⎡⎣k1 (1− cos φ ) + k 2 (1− 2 cos φ ) + k 4 (1− cos 3φ ) ⎤⎦
φ
+ ∑ ⎡k 2 ( χ − χ o ) ⎤
2
χ⎣ ⎦
+ ∑ ⎡⎣k (b − bo ) (b '− b 'o ) ⎤⎦ + ∑ ⎡⎣k (b − bo ) (θ − θo ) ⎤⎦
b,b ' b,θ
+ - + -
instantaneous induced
fluctuation (dipole) dipole
here here
Although the time-averaged dipole for each atom is zero, the net potential is not.
c6 c8 c10
U att = − − − −
r6 r8 r 10
c12 and c9
U rep = + U rep = +
r 12 r9
The exponential repulsion has been shown to better predict VLE properties of
hydrocarbons (Errington and Panagiotopoulos, J. Phys. Chem. B 103, 6314,
1999).
They are usually turned off (or scaled down) for the so-called 1-2, 1-3, and 1-4
intra-molecular interactions.
4 5
2 3
⎡⎛ σ ⎞ 12 ⎛ σ ⎞ 6 ⎤
ULJ = 4ε ⎢ ⎜ ⎟ − ⎜ ⎟ ⎥
⎢⎣⎝ r ⎠ ⎝r⎠ ⎥
⎦
-21
15x10
Exponential-6
Lennard-Jones
Energy (J)
10
parameters from
5
Reed and Gubbins,1973
2 3 4 5 6 7 8
r (Å)
Coulomb’s Law: the potential energy between two point charges is given by (any
physics text)
1 q1q 2
UCoul =
4πε 0 r
or
e 2 z1z 2
UCoul =
4πε 0 r
with
e = 1.602 x 10-19 C
e0 = 8.8542 x 10-12 C2/(N m2)
Obviously, ionic species have net charges that must be accounted for in this way.
Even in neutral molecules, there are generally regions that are electron-rich and
regions that are electron-poor. These regions can be represented by local
charges.
- + - + -
3 2 1 2 3
with
z1 = -2.3794
z2 = +5.2366
z3 = -4.0469
For convenience (especially for MD), the locations of the charges are usually
chosen to coincide with the vdW force centers. The N2 model above is one
exception.
+ ∑ ⎡H 2 (θ − θo ) + H 3 (θ − θo ) + H 4 (θ − θo ) ⎤
2 3 4
θ ⎣ ⎦
{ ( ) ( )
+ ∑ V1 ⎡1− cos φ − φ10 ⎤ +V2 ⎡1− cos 2φ − φ20 ⎤ +V3 ⎡1− cos 3φ − φ30 ⎤
φ ⎣ ⎦ ⎣ ⎦ ⎣ ⎦ ( )}
+ ∑ K 2 χ 2 + ∑ ∑ ⎡⎣Fbb ' (b − bo ) (b '− b 'o ) ⎤⎦ + ∑ ∑ ⎡⎣Fθθ ' (θ − θo ) (θ '− θ 'o ) ⎤⎦
χ b b' θ θ'
+ ∑ ∑ ⎡⎣Fbθ (b − bo ) (θ − θo ) ⎤⎦
b θ
+ ∑ ∑ ⎡⎣(b − bo ) (V1 cos φ +V2 cos 2φ +V3 cos 3φ ) ⎤⎦
b φ
qi qk ⎡A Bik ⎤
+ ∑ ∑ ∑ ⎡⎣K φθθ ' cos φ (θ − θo ) (θ '− θ 'o ) ⎤⎦ + ∑ + ∑ ⎢ ik − 6 ⎥
φ θ θ∍ i >k ε rik r 9
i >k ⎣ ik rik ⎦
For example, the interaction between the two blue atoms is certainly different in
the two configurations,
This expansion is rigorous when summed to all orders. However, evaluating the
higher-order terms becomes computationally expensive.
A three-body vdW correction term was derived by Axilrod and Teller (Reed and
Gubbins,1973).
r12 q2 r23
q1 q3
r13
This correction is generally deemed not sufficiently important to justify the (very
costly) rigorous sum of the three-body terms. Instead, the correction enters
through an “effective” two-body potential (Allen and Tildesley,1987). That is, the
parameters for the LJ two-body potential are chosen to reproduce the properties
of condensed phases, not isolated gas-phase pairs.
u (r ) = 0 0<r <∞
⎧∞ 0 <r <σ
u (r ) = ⎨
⎩0 σ <r <∞
⎧∞ 0 <r <σ
⎪
u (r ) = ⎨ − ε σ < r < gσ
⎪0 gσ < r < ∞
⎩
⎧∞ 0 <r <σ
⎪
u (r ) = ⎨ ⎛ σ ⎞ α
⎪− ε ⎜⎝ r ⎟⎠ σ <r <∞
⎩
⎧∞ 0 <r <σ
⎪
u (r ) = ⎨ ⎛ r −σ ⎞
⎪ − ε exp ⎜⎝ − ⎟ σ <r <∞
⎩ z ⎠
⎧∞ 0<r <d
⎪
u (r ) = ⎨ ⎡⎛ σ − d ⎞ 12 ⎛ σ − d ⎞ 6 ⎤
⎪4ε ⎢⎜⎝ r − d ⎟⎠ − ⎜⎝ r − d ⎟⎠ ⎥ d ≤r <∞
⎩ ⎢⎣ ⎥⎦
PR =
RT
V
(
1+ 4y + 10y 2 + 18.36y 3 + 28.3y 4 + )
where PR is the pressure due to repulsive interactions and y = b
4V
Carnahan and Starling combined the repulsive form suggested by hard spheres
with simple attractive contributions to obtain useful equations of state.
RT ⎡ 1+ y + y 2 − y 3 ⎤ a
CS w/ van der Waals attractive term: P= ⎢ ⎥− 2
V ⎢⎣ (1− y ) 3
⎥⎦ V
RT ⎡ 1+ y + y 2 − y 3 ⎤ a
CS w/ Redlich-Kwong attractive term: P= ⎢ ⎥ −
V ⎢⎣ (1− y )3 ⎥⎦ T (V + b )V
Many other forms of perturbed hard sphere equations have been proposed.
Z = Z hb + Z a
This approach considers separately the contributions to Z that come from hard
spheres, assembly into chains, dispersion (van der Waals) forces and specific
interactions (e.g., hydrogen bonding):
Ar = Ahs
r r
+Adisp r
+Achain r
+Aassoc
(The colors in this equation correspond to the colors of the arrows in the figure.)
κ Ai Bk = volume of association
ε Ai Bk = energy of association
The stability criteria resulting from the Second Law are that the entropy is
maximized and U, H, A and G are minimized at equilibrium. When the energy
functions are minimized, the natural independent variables are held constant.
dG = −SdT +VdP
Calculate the Gibbs energy along an isotherm relative to the Gibbs energy at a
reference pressure, P†, at the same temperature.
P
†
G −G = ⌠
⎮
⌡ V dP (constant temperature)
P†
CHEN 623 Lecture Notes (Date of last revision: 08/30/2021) 305
Equation of state behavior
V P
P V
Then:
and
P PV V
G −G = † ⌠
⎮
⌡ V dP = ⌠
⎮
⌡
d (PV ) − ⌠
⎮
⌡ P dV
P† P †V † V†
( )−
V
† †
= PV − P V ⌠
⎮
⌡ P dV
†
V
8Tr 3
Pr = − 2
3Vr − 1 Vr
PrVr Vr
G − G † = PcVc ⌠
⎮
⌡
d (PrVr ) − PcVc ⌠
⎮
⌡ Pr dVr
† † †
Pr Vr Vr
3
PcVc = Z cvdw RTc = RTc
8
† PrVr Vr
G −G 3 3
= ⌠
⎮
⌡
d (PrVr ) − ⌠
⎮
⌡ Pr dVr
RTc 8 † † 8 †
Pr Vr Vr
PrVr
⌠
⎮
⌡
d (PrVr ) = PrVr − P †V †
† †
Pr Vr
Vr Vr
Vr ⌠ dVr ⌠ dV 8 ⎡ 3V − 1 ⎤ ⎡1 1 ⎤
⌠
⎮
⌡ Pr dVr = 8Tr ⎮⎮ − 3 ⎮ 2 = Tr ln ⎢ r† ⎥ + 3 ⎢ − † ⎥
⎮ r
† ⌡ 3Vr − 1 ⌡ Vr 3 ⎢⎣ 3Vr − 1⎥⎦ ⎣Vr Vr ⎦
Vr † †
Vr Vr
( )
G − G† 3 † † ⎛ 3Vr − 1⎞ 9 ⎛ 1 1⎞
= Pr Vr − Pr Vr − Tr ln ⎜ † ⎟ − ⎜ − † ⎟
RTc 8 ⎝ 3Vr − 1⎠ 8 ⎝ Vr Vr ⎠
The equilibrium criteria state that the equilibrium condition has the lowest Gibbs
energy at constant T and P, so the equilibrium condition has two phases for a
range of densities at one pressure (the intersection point on the plot shown
above.
The molar Gibbs energy (µI for a pure fluid) of the liquid and vapor phases must
be equal at the equilibrium condition.
Compare the molar Gibbs energy at two points, 1, and 2, along the PV isotherm
V2
G ( 2 ) − G = P2V2 − P V −
† † † ⌠
⎮
⌡ P dV
V†
V1
G (1) − G = P1V1 − P V −
† † † ⌠
⎮
⌡ P dV
V†
Subtraction produces
V2
G ( 2 ) − G (1) = P2V2 − P1V1 − ⌠
⎮
⌡ P dV
V1
In the case of two coexisting phases, the Gibbs energies and pressures are
equal, therefore
V2sat
(
G ( 2 ) − G (1) = P eq V2sat −V1sat − ∫ P dV = 0
sat
) V1
Equilibrium States
1. σ
Estimate the vapor pressure, Pest
GV − GL
σ
Pest (new ) = Pest
σ
(old ) −
VV −VL
Note: The formula for the new estimate used in step 4 is based on the
observation that
d (GV − GL )
σ
= VV −VL
dPest
dP σ S α − S β Hα − H β
= α =
dT V −V β
(
T V α −V β )
This is a general result that applies to transitions in pure materials (e.g. vapor
pressures, sublimation pressures, melting pressures, solid-solid, etc.)
Case I:
Assumptions: ΔLV H is independent of temperature
vapor is an ideal gas
vapor volume is much larger than liquid volume
Case II:
⎡ ∂( Δσ H ) ⎤ ⎡ ∂( Δσ H ) ⎤
d ( Δσ H ) = ⎢ ⎥ dT + ⎢ ⎥ dP
⎣ ∂T ⎦P ⎣ ∂P ⎦T
⎡ ∂( Δσ H ) ⎤
⎢
⎡ ∂ Hα − H β
⎥ =⎢
( ) ⎤⎥ = CPα − CPβ = Δσ CPσ
⎣ ∂T ⎦P ⎢⎣ ∂T ⎥
⎦P
⎡ ∂( Δσ H ) ⎤ ⎡ ∂ Hα − H β
⎥ =⎢
( ) ⎤⎥ ⎡
α ⎛ ∂V α ⎞ ⎤ ⎡ β
= ⎢V −T ⎜
⎛ ∂V β ⎞ ⎤
⎢ ⎟ ⎥ − ⎢V −T ⎜ ∂T ⎟ ⎥
⎣ ∂P ⎦T ⎢⎣ ∂P ⎥
⎦T ⎣
⎢ ⎝ ∂T ⎠ P ⎥⎦ ⎢⎣ ⎝ ⎠ P ⎥⎦
⎛ ∂Δ V ⎞
= ΔσV −T ⎜ σ ⎟
⎝ ∂T ⎠ P
and
d ( Δσ H ) ⎡ ⎛ ∂ΔσV ⎞ ⎤ ⎛ dP ⎞
σ
= Δσ CP + ⎢ ΔσV −T ⎜
dT ⎣ ⎝ ∂T ⎟⎠ P ⎥⎦ ⎜⎝ dT ⎟⎠
⎡ ⎛ ∂Δ V ⎞ ⎤ ⎛ Δ H ⎞
= Δσ CP + ⎢ ΔσV −T ⎜ σ ⎟ ⎥ ⎜ σ ⎟
⎣ ⎝ ∂T ⎠ P ⎦ ⎝ T ΔσV ⎠
Δσ H Δσ H ⎡ T ⎛ ∂ΔσV ⎞ ⎤
= Δσ CP + − ⎢ ⎜⎝ ⎟⎠ ⎥
T T Δ
⎣ σ V ∂T P⎦
or
d ( Δσ H ) Δ H ⎡ ⎛ ∂lnΔσV ⎞ ⎤
== ΔσC P + σ ⎢1− ⎜⎝ ∂lnT ⎟⎠ ⎥
dT T ⎣ P⎦
This rule is based upon the observation that the entropy of boiling is
approximately the same for all fluids:
This observation may be understood by noting that the ratio of the vapor and
liquid molar volumes is of the same order of magnitude for all fluids.
For various real substances, the Gibbs surfaces in P,T space exhibit a variety of
behaviors.
General behavior:
The equilibrium phase has the lowest Gibbs energy at fixed T,P. At the
intersections of the surfaces, both phases can exist in equilibrium. Metastable
states also may exist in the vicinity of the intersection.
⎛ ∂G ⎞ ⎛ ∂G ⎞
⎜⎝ ⎟ = −S ⎜⎝ ⎟ =V
∂T ⎠ P ∂P ⎠ T
⎛ ∂2 G ⎞ ⎛ ∂S ⎞ CP ⎛ ∂2 G ⎞ ⎛ ∂V ⎞
⎜ 2⎟ = − ⎜⎝ ⎟⎠ = − ⎜ 2⎟ = − ⎜⎝ ⎟⎠ = −V κT
⎝ ∂T ⎠ P ∂T P T ⎝ ∂P ⎠ T ∂P T
The actual situation is more complex. The types of behavior actually observed
for CP are shown below.
First Order:
Third Order:
First Order:
Solid/Liquid
Vapor/Liquid
Solid/Vapor
Solid/Solid (Crystal Structure Changes)
Second Order:
Superconducting Transition
Lamda:
Helium I/Helium II (Superfluid transition in liquid 4He)
Order/disorder in b-brass
Order/disorder in quartz
Third Order:
Curie points of ferromagnets
The term critical region refers to the region in state variable space where a fluid
goes from conditions where multiple phases exist to conditions where only a
single phase exists. At the transition point, all properties of the individual phases
become identical.
The critical region was first observed by Andrews (1869) for carbon dioxide.
⎛ ∂2P ⎞
⎜ 2⎟ < 0
⎝ ∂V ⎠ T
• ⎛ ∂P ⎞
P ⎛ ∂2P ⎞
⎜ 2⎟ = 0
⎜ ⎟ <0
⎝ ∂V ⎠ T
⎝ ∂V ⎠ T
•
•
B (T > Tc)
•
⎛ ∂2P ⎞
⎜ 2⎟ < 0 ⎛ ∂P ⎞
⎝ ∂V ⎠ T ⎜ ⎟ =0 A (T < Tc )
⎝ ∂V ⎠ T
V
Equation of state isotherms are shown. Mathematically, to pass smoothly from
behavior A (two phases can be present for some conditions) to behavior B (one
phase only for all pressures) requires that at some temperature between TA and
TB:
⎛ ∂P ⎞ ⎛ ∂2 P ⎞
⎜⎝ ⎟⎠ = 0 and ⎜ 2 ⎟ = 0
∂V T ⎝ ∂V ⎠ T
1 ⎛ ∂V ⎞
Note also that: KT = − ⎜ ⎟ → ∞ at the critical point
V ⎝ ∂P ⎠ T
CHEN 623 Lecture Notes (Date of last revision: 08/30/2021) 327
3.5.2 Physical Observations:
Increasing Temperature
ρ < ρc
Increasing Temperature
ρ = ρc
Increasing Temperature
CHEN 623 Lecture Notes (Date of last revision: 08/30/2021) 328
3.5.3 Critical Opalescence
Density fluctuations (large) in the fluid phases lead to significant light scattering.
The scattering causes the fluid to turn black, to turn milky, or to glow depending
upon the fluid and the pattern of illumination.
3.5.4 Universality
Long range order and correlation occurs in the critical region. This order extends
over many molecules, hence the specific nature of the intermolecular force
becomes unimportant. Therefore the general behavior of all fluids will be the
same.
Other transitions, such as the Curie point for ferromagnetism also involve the
same types of long range order, therefore the critical behavior of magnets is
similar to that of fluids.
Theoretical calculations are the most tractable for the Ising model for
ferromagnets; there is an analytic solution for the 2 dimensional case, but a
numerical solution is required for 3 dimensions.
ρ
Heat Capacity at Constant Volume:
−α +
Cv ~ T − Tc
−α −2
Cv 2 ~ T − Tc CP
•
T
Critical
Isochore
ρ
CHEN 623 Lecture Notes (Date of last revision: 08/30/2021) 330
Densities
+
−γ
δ κ T ~ T − Tc
P − Pc ~ V − Vc
CP
• Tc
T
κT = ∞ ρ = 21 (ρv + ρL )
• • •
ρ = ρc
β
Δρ ~ Tc − T
ρv ρL
ρ
3.5.6 Law of Rectilinear Diameters
The mean of the saturated vapor and liquid densities varies linearly with the
temperature excursion from the critical temperature, i. e.,
ρ= 1
2 ( ρv + ρL ) = ρc + k1(T −Tc )
This often provides the best way to determine experimental values for the critical
density. The temperature at which the meniscus disappears is measured by
visual observation, and measured vapor pressures are extrapolated to the critical
temperature to get the critical pressure.
−α
CV ~ T −Tc
β
Δρ ~ T −Tc
−γ
κ T ~ T −Tc
δ
P − Pc ~ V −Vc
d 2P σ −θ
~ T −Tc
dT 2
ε
σ ~ Tc −T
e s 3/2 ~1.25
d 2P σ
2
0 0.125 0.109 0.21
dT
ρ -ρc
At Tc : Δρ = = 1− ρr < 0.08
ρc
T -Tc
At ρc : ΔT = = 1-Tr < 3 × 10-4
Tc
ρ -ρc
At Tc : Δρ = = 1− ρr < 0.1
ρc
Usually use T,P as the independent variables for mixture calculations. Then
(
M total = M total T ,P,n1,n2, ,nm )
where:
M total = total property of the mixture
n i = number of moles of specie i present
m = number of chemical species present in the mixture
Also:
M total = nM
where:
M = molar property of a mixture
m
n = ∑ n i = total number of moles of mixture
i =1
At constant T,P,
M total is homogeneous of order 1 with respect to the variables {n i }.
Implications of homogeneity
Note: T,P are not considered to be independent variables here because of the
restriction to constant T,P
⎡ ∂(nM ) ⎤ ⎡ ∂(nM ) ⎤
⎢ ⎥ = ⎢ ∂n ⎥
⎣ ∂( λ n )
i ⎦T ,P,n
k ≠i
⎣ i ⎦T ,P,n k ≠i
Result:
⎡ ∂(nM ) ⎤
m m
nM = ∑ n i ⎢ ⎥ = ∑ ni M i
i =1 ⎣ ∂n i ⎦T ,P ,n i =1
k ≠i
Notation:
lim
xi → 1 Mi =Mi
( x k ≠i → 0)
lim
xi → 0 Mi = M∞
i
( ∑ x k ≠i → 1)
where: M∞
i is the partial molar volume at infinite dilution
Note: The properties of the mixture are known if the partial molar properties of
all species are known.
Any mixture property depends upon the temperature, pressure, and the number
of moles of each chemical species present
⎡ ∂(nM ) ⎤ ⎡ ∂(nM ) ⎤ m
=⎢ ⎥ dT + ⎢ ⎥ dP + ∑ M i dn i
⎣ ∂T ⎦ P,n i ⎣ ∂P ⎦T ,n i i =1
m m m
Also, nM = ∑ n i M i requires that d (nM ) = ∑ n i dM i + ∑ M i dn i
i =1 i =1 i =1
⎡ ∂(nM ) ⎤ ⎡ ∂(nM ) ⎤ m m m
⎢ ∂T ⎥ dT + ⎢ ⎥ dP + ∑ M i dn i = ∑ n i dM i + ∑ M i dn i
⎣ ⎦ P,n i ⎣ ∂P ⎦T ,n i i =1 i =1 i =1
⎡ ∂(nM ) ⎤ ⎡ ∂(nM ) ⎤ m
⎢ ∂T ⎥ dT + ⎢ ⎥ dP − ∑ n i dM i = 0
⎣ ⎦ P,n i ⎣ ∂P ⎦T ,n i i =1
For those partial derivatives where the number of moles of every species is held
constant, n also is constant, so the previous result may be written as
⎡ ∂M ⎤ ⎡ ∂M ⎤ m
n⎢ ⎥ dT + n ⎢ ⎥ dP − ∑ n i dM i = 0
⎣ ∂T ⎦ P,n i ⎣ ∂P ⎦T ,n i i =1
m
∑ x i dM i = 0 (constant T ,P )
i =1
⎛ ∂(nM ) ⎞
Mi = ⎜
⎝ ∂n i ⎟⎠ T ,P,n
k ≠i
However, in most cases the information about the mixture properties is given as
a function of composition (mole fraction)
M = M ( x 1,x 2 ,,x m )
⎛ ∂(nM ) ⎞ ⎛ ∂M ⎞
Mi = ⎜ = M +n⎜
⎝ ∂n i ⎟⎠ T ,P,n ⎝ ∂n i ⎟⎠ T ,P,n
k ≠i k ≠i
⎛ ∂M ⎞ m ⎛ ∂M ⎞ ⎛ ∂x j ⎞
⎜⎝ ∂n ⎟⎠ = ∑⎜ ⎟ ⎜⎝ ∂n ⎟⎠
i T ,P,n j =1⎝ ∂x j ⎠ i T ,P,n
k ≠i T ,P,x k≠j k ≠i
nk nk
xk = m =
∑ nj n
j =1
Then
⎛ ∂x k ⎞ 1 ⎛ ∂nk ⎞ nk ⎛ ∂n ⎞
⎜⎝ ∂n ⎟⎠ = −
i n n ⎜⎝ ∂n i ⎟⎠ n n 2 ⎜⎝ ∂n i ⎟⎠ n
≠i ≠i ≠i
⎛ ∂nk ⎞ ⎧0 i ≠k
⎜⎝ ∂n ⎟⎠ = δ ik = ⎨
i n ⎩1 i =k
≠i
Then
⎛ ∂x k ⎞ 1 ⎛ ∂nk ⎞ nk ⎛ ∂n ⎞ δ ik nk
⎜⎝ ∂n ⎟⎠ = − = −
i n n ⎜⎝ ∂ni ⎟⎠ n 2 ⎜⎝ ∂ni ⎟⎠ n n2
≠i
so that
⎛ ∂M ⎞ m⎛ ∂M ⎞ ⎛ ∂x k ⎞
⎜⎝ ∂n ⎟⎠ = ∑⎜ ⎟ ⎜⎝ ∂n ⎟⎠
i T ,P,n k =1⎝ ∂x k ⎠ T ,P,x i T ,P,n
k ≠i j ≠k j ≠i
m⎛ ∂M ⎞ ⎡ δ ik nk ⎤
= ∑⎜ ⎟ ⎢⎣ n − n 2 ⎥⎦
k =1⎝ ∂x k ⎠ T ,P,x
j ≠k
1 ⎛ ∂M ⎞ ⎛ ∂M ⎞
m nk
= ⎜ − ∑⎜
n ⎝ ∂x i ⎟⎠ T ,P,x ⎟
k =1⎝ ∂x k ⎠ T ,P,x n2
j ≠i j ≠k
⎛ ∂M ⎞ ⎛ ∂M ⎞ m ⎛ ∂M ⎞
n⎜ =⎜ − ∑ xk ⎜
⎝ ∂n i ⎟⎠ T ,P,x ⎝ ∂x i ⎟⎠ T ,P,x k =1 ⎝ ∂x k ⎟⎠ T ,P,x
j ≠i j ≠i j ≠k
⎛ ∂(nM ) ⎞ ⎛ ∂M ⎞ m ⎛ ∂M ⎞
Mi = ⎜ =M +⎜ − ∑ xk ⎜
⎝ ∂n i ⎟⎠T ,P ,n ⎝ ∂x i ⎟⎠T ,P ,x k =1 ⎝ ∂x k ⎟⎠T ,P ,x
k ≠i k ≠i ≠k
l
Other methods for calculating M i exist, but in most cases they involve algebraic
elimination of xi and summations from which one or more terms are excluded.
The result given above is much more convenient for computer calculations for
mixtures with many species present. The notation here also is simpler and less
likely to cause confusion.
( )
m m m
= ∑ x i M i (T ,P,{x i }) − ∑ x i M i (T ,P ) = ∑ x i M i − M i
i =1 i =1 i =1
where:
Δm M (T ,P,{x i }) = property change on mixing
The property changes of interest include DmV, DmU, DmH, DmA, DmG and DmS
CHEN 623 Lecture Notes (Date of last revision: 08/30/2021) 347
4.1.3.2 Property Changes on Mixing for Ideal Gases
The property changes on mixing for the ideal gases will provide reference points
for more general discussions of mixture properties.
Begin with three separate ideal gases, each at T, P. Then mix, with no energy
transferred as heat during the mixing process. The process is shown
schematically below.
Before mixing, at T, P:
n = nA + nB + nC
Then
ΔigmV =V ig − ∑ x iViig = 0
The change in total internal energy resulting from the mixing process is
calculated from the First Law
Δig
mU
total
= Q +W = 0
Therefore
Δig
mU = 0
The enthalpy of mixing for ideal gases follows directly from the internal energy
and volume.
Δig
m H = Δ ig
m U + PΔ ig
mV = 0
⎛ nRT ⎞
total ⎜⎝ ⎟
Vmixture P ⎠ n 1
Note: = = =
Vitotal ⎛ ni RT ⎞ ni x i
⎜⎝ ⎟
P ⎠
Then
total
Vmixture 1
Δig
mSi
total
= n i R ln total = n i R ln = −n i R lnx i
Vi xi
and
mS = −R (nA lnx A + nB lnx B + nC lnx C )
nΔig
Δig
mS
m
so that the dimensionless form is = − ∑ x i lnx i
R i =1
Δig
mA
=
( Δ U ) −T ( Δ S )
ig
m
= ∑ x lnx
ig
m m
i i
RT RT i =1
Δig
m G
=
Δig
m (
H −T Δig
mS m ) (
= ∑ x i lnx i
)
RT RT i =1
We use the results derived for mixing of ideal gases to define a concept of ideal
solutions, which are mixtures for which the properties changes of mixing are
identical to those for mixing of ideal gases, even though the pure substances are
not ideal gases.
Δis
mU = 0 Δis m
m S = − x ln x
∑ i i
RT R i =1
P Δis
mV = 0 Δis
mA =
m
∑ x i ln x i
RT RT i =1
Δis
mH = 0 Δis
mG =
m
∑ x i ln x i
RT RT i =1
Δis
mG = RT ∑ x i lnx i
Giig =⎢
(
⎡ ∂ nG ig ) ⎤⎥ = Giig
m
+ RT lnx i + RT ∑ nk
∂
( lnx k )
⎢ ∂ni ⎥ k =1 ∂ni
⎣ ⎦T ,P,nk ≠i
m ∂(lnx k ) m nk m m
∑ nk
∂n
= ∑
nx
(δ ik − x k ) = ∑ δ ik − ∑ x k
k =1 i k =1 k k =1 k =1
m ⎫
∑ δ ik = 1⎪
k =1 ⎪ m ∂(lnx k )
⎬ ⇒ ∑ nk =0
m ∂n
∑ x k = 1⎪
k =1 i
k =1 ⎪⎭
so that
G isi = G i + RT ln x i
Preserve the same general form for real mixtures by including the activity
coefficient in the partial molar volume.
G i = G i + RT lnγ i x i
4.1.3.5 Activity
The partial molar Gibbs energy also can be written in the form
G i = G i + RT lnα i
In the general case, the definition of the property change on mixing must
accommodate reference states that differ from the final state of the mixture:
( )
m m
Δ m† M = M − ∑ x i M †i = ∑ x i M i − M †i
i =1 i =1
G isi = G †i + RT ln x i
and
Define excess properties based upon the difference between the property
changes on mixing for the real solution and for the ideal solution.
Note that both the excess properties and the property changes on mixing are
functions of temperature, pressure and composition.
The partial molar excess properties are related to the real mixture and ideal
solution values by
M Ei = M i − M is
i
G i = G †i + RT lnγ i†x i
G isi = G †i + RT lnx i
( )
†
GEi = G i − G isi = RT lnγ i†
Therefore
G Ei
= lnγ i
RT
E
⎛ ∂GE ⎞ ⎛ ∂GEi ⎞
S = −⎜ ⎟ S Ei = −⎜ ⎟
⎝ ∂T ⎠ P ⎝ ∂T ⎠ P
E
⎛ ∂G E ⎞ ⎛ ∂GEi ⎞
V =⎜ ⎟ V Ei =⎜ ⎟
⎝ ∂P ⎠ T ⎝ ∂P ⎠ T
H E = GE +TS E H Ei = G Ei +TS Ei
AE = GE − PV E AEi = G Ei − PV iE
U E = GE − PV E +TS E U Ei = G Ei − PV iE +TS Ei
GE m
= ∑ x i lnγ i
RT i =1
PV E m ⎛ ∂lnγ i ⎞
= ∑ xi ⎜ ⎟
RT i =1 ⎝ ∂lnP ⎠ T
HE m ⎛ ∂lnγ i ⎞
= − ∑ xi ⎜ ⎟
RT i =1 ⎝ ∂lnT ⎠ P
SE m ⎡ ⎛ ∂lnγ i ⎞ ⎤
= − ∑ xi ⎢lnγ i + ⎜⎝ ⎟⎠ ⎥
R i =1 ⎣ ∂lnT P⎦
AE m ⎡ ⎛ ∂lnγ i ⎞ ⎤
= ∑ xi ⎢lnγ i + ⎜⎝ ⎟⎠ ⎥
RT i =1 ⎣ ∂lnT P⎦
UE m ⎡⎛ ∂lnγ i ⎞ ⎛ ∂lnγ i ⎞ ⎤
= − ∑ xi ⎢⎜⎝ ⎟⎠ + ⎜⎝ ⎟⎠ ⎥
RT i =1 ⎣ ∂lnT P ∂lnP T⎦
ΔmG G − Σx i Gi Σx i G i − Gi
= =
( )
RT RT RT
=
1
RT { (
Σx i ⎡ G i − G ig
⎣ i + G ig
i ) (
− Giig + G ig
i − Gi ⎤
⎦) ( )}
=
1
RT { (
Σx i ⎡ G i − G ig
⎣ i + RT lnx i + G ig
i )− Gi ⎤
⎦ ( )}
⎡ GRi GR ⎤ Δis
mG
⎡ GRi GR ⎤ Δis
mG G
E
= Σx i ⎢ − i
⎥ + Σx i lnx i = + Σx i ⎢ − i
⎥= +
⎢⎣ RT RT ⎥⎦ RT ⎢⎣ RT RT ⎥⎦ RT RT
Therefore:
GE ⎡ GRi GiR ⎤
= Σx i ⎢ − ⎥
RT ⎢⎣ RT RT ⎥⎦
CHEN 623 Lecture Notes (Date of last revision: 08/30/2021) 362
By the definition of a residual function:
Pure species:
P P
GiR Gi − Giig ⌠VR ⌠ dP
= = ⎮ i
⎮ RT
dP = ⎮⎮ ( Z i − 1)
RT RT ⌡ ⌡ P
0 0
Mixture:
P P
R ig ⌠VR
G G −G ⌠ dP
= = ⎮
⎮ RT
dP = ⎮⎮ ( Z − 1)
RT RT ⌡ ⌡ P
0 0
P P
GiR Gi − Giig ⌠VR
( Z i − 1) dP
⌠
= = ⎮ i
⎮ RT
dP = ⎮
⎮
RT RT ⌡ ⌡ P
0 0
⎡P R P ⎤ P
G E ⎡ GRi
⎤ m ⎢⌠ V i GR ⌠V R
⎥ m ⌠ V R −V R
= Σx i ⎢ − ⎥ = ∑ x i ⎢ ⎮⎮
i
dP − ⎮⎮ i
dP = ∑ x i ⎮⎮ i i
dP
RT ⎢⎣ RT RT ⎥⎦ i =1 ⌡ RT ⌡ RT ⎥ i =1 ⌡ RT
⎢⎣ 0 0 ⎥⎦ 0
P
m ⌠ V −V
= ∑ x i ⎮⎮ i i
dP
i =1 ⌡ RT
0
Also
P
GEi ⌠ V −V
= ⎮ i
⎮
i
dP = lnγ i
RT ⌡ RT
0
Define fugacity to obtain the same form for the real fluid as for the ideal gas
dG = RTd (lnf )
Then
G (P ) − G P † ( ) = ln⎛ f ⎞
⎜⎝ † ⎟⎠
RT f
Using these forms
⎛ f ⎞
dGR = dG − dG ig = RTd (lnf ) − RTd (lnP ) = RTd ⎜ ln ⎟
⎝ P⎠
⎡ R †
GR ⎢ G P
=
⎛ f† ⎞⎤( ) ⎛f ⎞
− ln ⎜ † ⎟ ⎥ + ln ⎜ ⎟
RT ⎢ RT ⎝P ⎠⎥ ⎝P⎠
⎣ ⎦
⎛ R ⎞
⎛f ⎞
lim ⎜ G ⎟ = ⎟ =0
lim ln ⎜
P→0 ⎜ RT ⎟ P→0 ⎝P ⎠
⎝ ⎠
f
φ= where lim
P →0
φ =1
P
Pure Substance:
lim ⎛ fi ⎞
dGi = RTd lnfi P→0 ⎜ =1
⎝ P ⎟⎠
Mixture:
lim ⎛ f ⎞
dG = RTd lnf P→0 ⎜ =1
⎝ P ⎟⎠
lim ⎛ fˆi ⎞
dGi = RTd lnfˆi P→0 ⎜ ⎟ =1
⎝ xiP ⎠
fˆi x i φ̂i
α i† = =
f i† φi†
fˆi φ̂i
γ i† = =
x i f i† φ i†
fˆiis = x i f †i
Consider activity cases based upon two different reference states, one denoted
by † and the other by ∆. The Gibbs energy change on mixing as described by
each choice of reference states is
†
Δm G
=
(
∑ x i Gi − Gi† ) Δ
Δm G
=
(
∑ x i Gi − GiΔ )
RT RT RT RT
So
†
Δm G ΔmΔ
G x i ⎡⎣GiΔ − Gi† ⎤⎦
− =∑
RT RT RT
Note:
†
Δm G Δis,†
m G G
E ,†
G E ,† m
= + = + ∑ x i lnx i
RT RT RT RT i =1
Δ
Δm G Δis,Δ
m G G
E ,Δ
G E ,Δ m
= + = + ∑ x i lnx i
RT RT RT RT i =1
GE .† GE ,Δ Δm†
G ΔmΔ
G m ⎡ GiΔ − Gi† ⎤
− = − = ∑ xi ⎢ ⎥
RT RT RT RT i =1 ⎢⎣ RT ⎥⎦
also
G E .† G E ,Δ m
†
m
Δ
m γ i†
− = ∑ x i lnγ i − ∑ x i lnγ i = ∑ x i ln Δ
RT RT i =1 i =1 i =1 γi
Then
γ i† GiΔ − Gi† Δ
⎡ GiΔ − Gi† ⎤
ln = ⇒ γ i† = γ i exp ⎢ ⎥
γ iΔ RT ⎢⎣ RT ⎥⎦
Δ
⎡ GiΔ − Gi† ⎤
αi = γ i xi ⇒ α i† = α i exp ⎢ ⎥
⎢⎣ RT ⎥⎦
fˆi fˆi
γ i† = γ iΔ =
x i f i† x i f iΔ
fˆi fˆi
α i† = α iΔ =
f i† f iΔ
Equilibrium condition is that the total Gibbs enery of the system is at a minimum
for a given T, P. Consider G(x) at fixed T, P.
Δis
mG
m
Example: = ∑ x i lnx i
RT i =1
ΔismG
For ideal solutions,such mixtures, is negative for all compositions.
RT
0 1
x
The chord connecting the values for any two compositions gives the total Gibbs
energy for a two-phase mixture, where the two phases have the compositions of
the end points of the chord.
Totally miscible:
0 1
x
At all overall compositions, ∆mG < 0, and G(x) has a lower value than any linear
combination of other compositions, therefore such a mixture always exists as a
single phase at this T, P.
0 1
x
0 1
x
For xA< x < xB, the tangent AB has a lower Gibbs energy than G(x). This system
will exist as a single phase for x < xA and x > xB, but as a two-phase mixture with
phase compositions xA and xB for xA< x < xB.
⎫ ΔmG < 0
⎪
Mixtures are completely miscible if: ∂2 Δ G ⎬ for all 0 ≤ x ≤ 1
m
2
> 0⎪
∂x ⎭
CHEN 623 Lecture Notes (Date of last revision: 08/30/2021) 377
4.2.2 Phase Equilibrium Calculation
x iβ
(
K i ≡ α = K i T ,P ,{x iα },{x iβ }
xi
)
The equilibrium conditions resulting from the First and Second Laws of
Thermodynamics provide methods for calculating Ki for mixtures.
Begin with the requirement that the mole fractions in a phase must sum to unity
m
∑ xi = 1
i =1
Use the equilibrium relation in the form of the partition coefficient to eliminate the
mole fractions of the unknown phase.
x iβ m
β
m
Ki = ⇒ ∑ xi = ∑ K i x iα = 1
x iα i =1 i =1
m ⎛ β⎞ m
x
∑ K i x iα = 1 or ∑ ⎜ i ⎟ = 1
i =1 i =1⎝ K i ⎠
Use the material balances to describe the state when two coexisting phases are
present. The overall composition is known but the compositions of the coexisting
phases are unknown.
nj
j
where: ψ =
n
Use the partition coefficient and the total material balance to eliminate phase b
from the species material balance
( )(
x iαψ α + K i x iα 1− ψ α = x i )
CHEN 623 Lecture Notes (Date of last revision: 08/30/2021) 381
Rearrange to obtain
xi
x iα =
1+ ψ α (K i − 1)
Working equations:
m ⎡ xi ⎤
∑⎢ ⎥ =1
i =1 1+ ψ (K i − 1)
α
⎣ ⎦
m ⎡ K i xi ⎤
∑⎢ ⎥ =1
i =1 1+ ψ (K i − 1)
α
⎣ ⎦
m ⎡ (K i − 1) x i ⎤
∑⎢ ⎥=0
i =1 1+ ψ (K i − 1)
α
⎣ ⎦
xi
x iα =
1+ ψ α (K i − 1)
and
x iβ = K i x αi
Gamma-Phi (G-F)
Use the fugacity coefficient description for one phase (fluid) and the activity
coefficient description for the other (solid or liquid):
α α β β † x iα γ iβf i† γ iβPisat( β )
x φ̂ P = x γ f ⇒ β = α =
φ̂i P Φi P
i i i i i
xi
where
φ̂iα Pi sat( β ) φ̂iα Pi sat( β ) φ̂iα ⎡ Vi β ⎤
Φi = =
⎡ Vi β sat( β ) ⎤
=
φisat( β )
exp (
⎢ RT i P sat( β )
− )
P ⎥
fi † sat( β ) sat( β )
φi Pi exp ⎢
⎣ RT
(P − Pi ) ⎥
⎦
⎣ ⎦
Equation of State
x iα φ̂iβ
x iα φ̂iα P = x iβφ̂iβ P ⇒ =
x iβ φ̂iα
CHEN 623 Lecture Notes (Date of last revision: 08/30/2021) 385
4.2.3.1 Vapor Phase Description
For the vapor phase, the most common representation is the fugacity coefficient
notation
fˆiv = y i φ̂vi P
( ) dP
⌠
where lnφ̂vi = ⎮⎮ Z i − 1
⌡ P
0
If
BP
Z = 1+
RT
For liquids and solids, the activity coefficient notation often is used. For liquids,
⎡ 1 P ⎤
fˆiL = γ iL x i f iL = γ iL x i φiv Pi v exp ⎢ ⎮ V dP ⎥
⌠
⌡ i
L
⎢ RT Pi v ⎥
⎣ ⎦
For solids
fˆis = γ is x i f is
where
⎧ ⎡ P ⎤
⎪φ sv P sv exp ⎢ 1 ⎮ V dP ⎥
⌠ s
T ≤ Titp
⎪ i i ⎢ RT
⌡ i
⎥
⎪ ⎣ Pi sv ⎦
⎪
f is = ⎨
⎪
⎪ v v ⎡ ⎧⎪Pi s P ⎫⎪ ⎤
⎢ 1
⎨ ⌠⎮⌡ Vi dP + ⌠⎮⌡ Vi dP ⎬ ⎥ T > Titp
L s
φ
⎪ i i P exp
⎪⎩ ⎢ RT ⎪⎩Pi v Pi s ⎪⎭ ⎥⎦
⎣
The general relationship for representing vapor-liquid equilibrium (VLE) using the
gamma-phi (G-F) notation is
y i γ iLφiσ Piσ ⎡ 1 P ⎤
Ki ≡ = v
exp ⎢ ⎮ V dP ⎥
⌠
⌡ i
L
xi φ̂i P ⎢ RT Piσ ⎥
⎣ ⎦
where
⎡Piσ ⎤
⎢ ⌠ dP ⎥
φiσ = exp ⎢ ⎮⎮ ( Z i − 1)
⌡ P ⎥
⎢⎣ 0 ⎥⎦
⎡P ⎤
(dP
)
⌠
φ̂iv = exp ⎢ ⎮⎮ Z i − 1 ⎥
⎢⌡ P ⎥
⎣0 ⎦
⎡ 1 P ⎤
1. exp ⎢ ⌠
⎮
⌡ Vi dP ⎥ ≅ 1
L
⎢ RT Piσ ⎥
⎣ ⎦
φiσ
2. v
≅1
φ̂i
3. γ iL = 1 (ideal solution)
⎡ 1 P ⎤
4. γ iLφiσ Piσ exp ⎢ ⌠
⎮
⌡ Vi dP ⎥ = constant
L
⎢ RT Piσ ⎥
⎣ ⎦
y i γ iLφiσ Piσ ⎡ 1 P ⎤ Pσ
Ki ≡ = v
exp ⎢ ⎮ V dP ⎥ = i
⌠
⌡ i
L
xi φ̂i P ⎢ RT Piσ ⎥ P
⎣ ⎦
Note also that
γ iLφiσ ⎡ 1 P ⎤
v
exp ⎢ ⎮ V dP ⎥ = 1 → Raoult's Law
⌠
⌡ i
L
φ̂i ⎢ RT Piσ ⎥
⎣ ⎦
(This condition does not require that 1, 2 and 3 above be satisfied individually.)
Raoult's Law is
y i P iσ
fˆ (RL ) = x i P
i
L σ
i and K i
RL
= =
xi P
This relation often provides a good approximation for estimating phase behavior.
This result is known as Henry's law, and it provides a good representation of the
fugacity of a species at nearly infinite dilution.
Special Cases
A=B=C=D=0 ® ideal solution
B=C=D=0 ® symmetric
C=D=0 ® two-suffix Margules
D=0 ® three-suffix Margules
This form was derived from the van der Waals equation:
GE x 1x 2
=
RT 1 1
x1 + x
A12 A21 2
4.2.3.4.1.3 Wilson
GE
= −x 1 ln ( x 1 + Λ12 x 2 ) − x 2 ln ( Λ21x 1 + x 2 )
RT
GE ⎡ x +V x / V ⎤ ⎡V x / V + x 2 ⎤
= x 1 ln ⎢ 1 2 2 1 ⎥ + x 2 ln ⎢ 1 1 2 ⎥
RT ⎣ x 1 + Λ12 x 2 ⎦ ⎣ Λ 21x 1 + x 2 ⎦
GE ⎡ τ 21G21 τ 12G12 ⎤
= x 1x 2 ⎢ + ⎥
RT ⎣ x 1 + G21x 2 G12 x 1 + x 2 ⎦
4.2.3.4.1.6 UNIQUAC
GE ⎡ φ qz θ ⎤
= x 1 ⎢ln 1 + 1 ln 1 − q1ln (θ1 + θ 2τ 21) ⎥
RT ⎣ x 1 2 φ1 ⎦
⎡ φ q z θ ⎤
+x 2 ⎢ln 2 + 2 ln 2 − q 2 ln (θ1τ 12 + θ 2 ) ⎥
⎣ x2 2 φ2 ⎦
Note that the form is the same as van Laar, but the parameters, di, are
determined for pure substances, allowing for prediction of mixture behavior
GE
=
(δ1 − δ 2 )2
RT ⎛ 1 1 ⎞
RT ⎜ +
⎝ V1x 1 V2 x 2 ⎟⎠
GE ⎛ ⎞
Wilson = − ∑ x i ln ⎜ ∑ x k Λik ⎟
RT i ⎝k ⎠
GE ⎡ ⎛ ⎞⎤
T-K Wilson = ∑ x i ⎢ ∑ x kVk / Vi − ln ⎜ ∑ x k Λik ⎟ ⎥
RT i ⎣k ⎝k ⎠⎦
GE ⎡ ⎤
NRTL = ∑ x i ⎢ ∑ τ ji G ji x j / ∑ Gki x k ⎥
RT i ⎣ j k ⎦
GE φ z θ ⎛ ⎞
UNIQUAC = ∑ x i ln i + ∑ qi x i ln i − ∑ qi x i ln ⎜ ∑ θk τ ki ⎟
RT i xi 2 i φi i ⎝k ⎠
GE
Scatchard-Hildebrand =
1
RT RT i
(
∑ x iVi δ i − δ ) 2
A difficulty here is that our previous expressions relating the fugacity to the
equation of state involve pressure as the variable of integration. This is feasible
for the vapor phase, but, as seen with the pure fluid vapor pressures, it is not
feasible for the liquid phase. An expression utilizing the volume as the variable
of integration must be derived.
Begin with
⎡ ∂(nM ) ⎤ ⎡ ∂(nM ) ⎤ m
d (nM ) = ⎢ ⎥ dT + ⎢ ⎥ dP + ∑ M kdnk
⎣ ∂T ⎦P,nk ⎣ ∂P ⎦T ,nk k =1
Divide by dni
⎡ ∂(nM ) ⎤ dT ⎡ ∂(nM ) ⎤ dP
=⎢ ⎥ + ⎢ ∂P ⎥ +Mi
⎣ ∂T ⎦P,nk dn i ⎣ ⎦T ,nk i
dn
⎡ ∂(nM ) ⎤ ⎡ ∂(nM ) ⎤ ⎡ ∂P ⎤
⎢ ∂n ⎥ = ⎢ ∂P ⎥ ⎢ ∂n ⎥ +Mi
⎣ i ⎦T ,nV ,nk ≠i ⎣ ⎦T ,nk ⎣ i ⎦T ,nV ,nk ≠i
Therefore
⎡ ∂(nM ) ⎤ ⎡ ∂(nM ) ⎤ ⎡ ∂P ⎤
Mi = ⎢ ⎥ − ⎢ ⎥ ⎢ ∂n ⎥
⎣ ∂n i ⎦T ,nV ,nk ≠i ⎣ ∂P ⎦T ,nk ⎣ i ⎦T ,nV ,nk ≠i
Now take
M = lnφ → M i = lnφ̂i
P
⌠ dP
lnφ = ⎮
⎮ ( Z − 1)
⌡ P
0
ZRT
P=
V
RT ZRT ZR
dP = dZ − 2 dV + dT
V V V
dP dZ dV dT
= − +
P Z V T
For any isothermal situation, the last term is zero, and the integral becomes
P Z V Z V
⌠
lnφ = ( Z − 1)
dP ⌠ ( Z − 1) dZ − ⌠ dV ⌠ Z
dZ ⌠ ⌠ dV
⎮
⎮ = ⎮
⎮
⎮
⎮ ( Z − 1) = ⌡ dZ −
⎮ − ⎮⎮ ( Z − 1)
⎮
⎮
⌡ P ⌡ Z ⌡ V 1
Z ⌡ ⌡ V
0 1 ∞ 1 ∞
Then
V
⌠ dV
lnφ = ( Z − 1) − ln Z − ⎮⎮ ( Z − 1)
⌡ V
∞
Therefore
n ⎡ ∂Z ⎤ 1 ⎡ ∂(nZ ) ⎤
⎢ ∂n ⎥ = ⎢ ∂n ⎥ −1
Z ⎣ i ⎦T ,nV ,nk ≠i Z ⎣ i ⎦T ,nV ,nk ≠i
V
⎡ ∂(n lnφ ) ⎤ Z − 1⎛ ∂(nZ ) ⎞ ⎡⎛ ∂(nZ ) ⎞
⌠
⎮
⎤ dV
⎢ ∂n ⎥ = ⎜ ⎟ − lnZ − ⎢⎜
⎮
⎟ − 1⎥
⎣ i ⎦T ,nV ,nk ≠i Z ⎝ ∂n i ⎠ T ,nV ,n ⎢⎝ ∂n i ⎠ T ,nV ,n
⎮
⎮ ⎥V
k ≠i ⎣
⌡ k ≠i ⎦
∞
⎡ ∂(nM ) ⎤
This is the actual form of the ⎢ ⎥ term.
⎣ ∂n i ⎦T ,nV ,nk ≠i
Also
⎡ ∂P ⎤ ⎡ ∂ ⎛ ZRT ⎞ ⎤ ⎡ ∂ ⎛ nZRT ⎞ ⎤
⎢ ∂n ⎥ = ⎢ ⎜ ⎟⎥ = ⎢ ⎜ ⎟⎥
⎣ i ⎦T ,nV ,nk ≠i ⎣ ∂n i ⎝ V ⎠ ⎦T ,nV ,nk ≠i ⎣ ∂n i ⎝ nV ⎠ ⎦T ,nV ,nk ≠i
RT ⎡ ∂(nZ ) ⎤
=
nV ⎢⎣ ∂n i ⎥⎦T ,nV ,n
k ≠i
⎡ ⎧ P ⎫⎤
⎡ ∂(n lnφ ) ⎤ ∂ ⎪ ⌠ dP ⎪ ⎥ n ( Z − 1) nV ( Z − 1)
⎢ ∂P ⎥ =⎢ ⎨ ⎮
n ⎮ ( Z − 1) ⎬⎥ = =
⎣ ⎦T ,nk ⎢ ∂P P ⎪ P RT Z
⎢⎣ ⎪⎩ ⌡0 ⎥
⎭ ⎦T ,nk
Z − 1 ⎡ ∂(nZ) ⎤
= ⎢ ⎥
Z ⎢⎣ ∂n i ⎥⎦T,nV,n
k≠ i
lnφ̂i =
( Z − 1) ⎡ ∂(nZ ) ⎤ − lnZ
Z ⎢ ∂n ⎥
⎣ i ⎦T ,nV ,nk ≠i
V
⌠
⎮
⎧⎪ ⎡ ∂(nZ ) ⎤ ⎫⎪ dV ( Z − 1) ⎡ ∂(nZ ) ⎤
−⎮ ⎨⎢ ⎥ − 1⎬ −
⎮
⎮ ⎪⎩ ⎣ ∂n i ⎦T ,nV ,nk ≠i ⎪⎭ V Z ⎢⎣ ∂n i ⎥⎦T ,nV ,n
⌡ k ≠i
∞
which, after canceling the first and last terms, becomes
∞
⌠⎧ ⎫⎪ dV
lnφ̂i =
(
⎮ ⎪ ⎡ ∂ nZ )⎤ − 1⎬ − ln Z
⎮ ⎨⎢ ⎥
⎮ ⎣ ∂n i ⎦T ,nV ,nk ≠i V
⌡⎩⎪ ⎭⎪
V
ρ
⌠⎧ ⎫⎪ d ρ
=
(
⎮ ⎪ ⎡ ∂ nZ )⎤ − 1⎬ − ln Z
⎮ ⎨⎢ ⎥
⎮ ⎣ ∂n i ⎦T ,nV ,nk ≠i ρ
⌡⎪⎩ ⎪⎭
0
( Z i − 1) dV
⌠
lnφ̂i = ⎮
⎮ − lnZ
⌡ V
V
where:
⎡ ∂(nZ ) ⎤
Z i = ⎢ ⎥
⎣ ∂n i ⎦T ,nV ,nk ≠i
b ⎡⎛ V − b ⎞ ⎤ a 1 ⎛ ai b i ⎞ ⎡V + σ b ⎤
lnφ̂i = i ( Z − 1) − ln ⎢⎜ Z
⎟⎠ ⎥ + ⎜ 1+ − ⎟ ln ⎢
b ⎝
⎣ V ⎦ bRT ( ε − σ ) ⎝ a b ⎠ ⎣ V + εb ⎥⎦
V aV
Z= −
V − b RT (V + ε b ) (V + σ b )
α ΩaR 2Tc,i
2
ai =
Pc,i
ΩbRTc,i
bi =
Pc,i
1
α RK =
T
( )( )
2
α SRK = ⎡1+ 0.48508 + 1.5517ω − 0.1561ω 2 1− Tr ⎤
⎣ ⎦
( )( )
2
α PR = ⎡1+ 0.37464 + 1.54226ω − 0.26992ω 2 1− Tr ⎤
⎣ ⎦
RK SRK PR
e 0 0 1− 2
s 1 1 1+ 2
1 1
Wa
9 ( 3
2 −1) 9 ( 3
2 −1) 0.45724
3
2 −1 3
2 −1
Wb 0.07780
3 3
The relations between the cubic equation parameters for the mixture and the
composition and parameters for the pure fluids are known as mixing rules.
These provide the key to accuracy in descriptions of mixture behavior.
4.2.3.5.2.1 Conventional
The form used most frequently are known as the Lorentz-Berthelot mixing rules,
which have the following form
bm = ∑ x i b i
i
2
⎛ ⎞
am = ⎜ ∑ x i ai ⎟ = ∑ ∑ ai x i x
⎝ i ⎠ i
(
uij (r ) = fij u0 r / hij1/ 3 )
h = ∑ ∑ x i x j hij
i j
1
f = ∑ ∑ x i x j fij hij
h i j
c cd
+d +2
V
Z= − RT RT
V −b
(V + b ) + ⎛⎜⎝ ⎞⎟⎠ (V − b )
b
V
aii α PR
d ii = = 0.45724 α PR RTc,i / Pc,i
RTc
( )
c ii c jj
c m = ∑ ∑ x i x j c ij c ij = 1− k ij bij
i j bii b jj
3
3 ⎡ d ii + d jj ⎤
1/ 3 1/ 3
d m = ∑ ∑ x i x j d ij
i j
(
d ij = 1− mij ⎢
⎢⎣
) 2
⎥
⎥⎦
For highly non-ideal solutions, the phase behavior descriptions obtained from
equations of state can be improved significantly by incorporating known
information about the behavior of the liquid phase. For many systems, low
pressure VLE measurements at one or more temperatures provide activity
coefficient data.
The excess Gibbs function is given by
GE ⎛ φ̂ ⎞
= ∑ x i lnγ i = ∑ x i ln ⎜ i ⎟ = ∑ x i lnφ̂i − ∑ x i lnφi
RT i i ⎝ φi ⎠ i i
G ⎡
E
V
⌠ dV ⎤ ⎡ V
⌠ dV ⎤
= ⎢ Z − lnZ − ⎮⎮ ( Z − 1) ⎥ − ∑ xi ⎢ Z i − lnZ i − ⎮⎮ ( Z i − 1) ⎥
RT ⎢ ⌡ V ⎥⎦ i ⎢⎣ ⌡ V ⎥⎦
⎣ ∞ ∞
The interest here is in cubic equations of state, which follow the general form
V a
Zk = − fk (T ,V ) k = RK, SRK, PR, etc.
V − b bRT
GE ⎛Z⎞ ⎡ 1− b / V ⎤
= Z − ∑ x i Z i − ∑ x i ln ⎜ ⎟ − ∑ x i ln ⎢ ⎥
RT i i ⎝ Zi ⎠ i ⎣ 1− bi / Vi ⎦
⎛ a ⎞ ⎛ ai ⎞
+⎜ ( ) − ∑ C (Vi )
⎝ bRT ⎟⎠
C V x i⎜
⎝ bi RT ⎟⎠
where
V
⌠ dV
C (V ) = ⎮ fk
⎮ (T ,V )
⌡ V
∞
E
GEOS = GγE
can be used at only one pressure. There are two potential choices
The available activity coeficient data usually are at low pressures, but there is no
zero pressure root for the cubic equations for reduced temperatures above 0.8 to
0.85, depending on the equation.
The infinite pressure limit yields an equation relating a and b to the activity
coefficient data. The function C assumes a numerical value, C*, in the limit,
which varies with equation of state.
Vidal and co-workers were first to use the infinite pressure limit, resulting in the
Huron-Vidal mixing rules
b = ∑ x i bi
i
⎡ ⎛ ai ⎞ Gγ ⎤
E
a = b ⎢∑ xi ⎜ ⎟ + ⎥
⎢⎣ i ⎝ bi ⎠ C * ⎥⎦
References:
The Huron-Vidal mixing rules work well for some systems, but not for others,
especially at higher pressures and temperatures (approaching the critical region
or the cricondentherm).
Wong and Sandler [AIChE J., 38, 671 (1992)] proposed an alternative set of
mixing rules following this approach, but recognizing that the excess Helmholtz
energy exhibits much less variation with pressure than does the excess Gibbs
energy.
At low pressures
Furthermore, because the excess Helmholtz energy depends only weakly upon
pressure
E ⎡V V ⎤
A Z ⎢⌠ dV ⌠ dV ⎥
RT
= − ∑ x i ln − ⎮⎮ ( Z − 1)
Z ⎢ ⌡ V
− ∑ xi ⎮
⎮
⌡
( Z i − 1) V ⎥
i i i
⎣∞ ∞ ⎦
B = ∑ ∑ x i x k Bik
i k
a ⎡ aik ⎤
b− = ∑ ∑ xi xk b −
⎢⎣ ik RT ⎥⎦
RT i k
1 ⎛ Gγ ⎞
E
a ai
= ⎜ ⎟ + ∑ x i
bRT C * ⎝ RT ⎠ i bi RT
a ⎛ a ⎞
b− = ∑ ∑ xi xk ⎜ b − ⎟
RT i k ⎝ RT ⎠ ik
⎛ a ⎞ 1 ⎡⎛ ai ⎞ ⎛ a ⎞ ⎤
⎜⎝ b − ⎟⎠ = ⎢ ⎜⎝ i
b − +
⎟⎠ ⎜⎝
b − ⎟⎠ ⎥ (1− k i )
RT i 2 ⎣ RT RT ⎦
CHEN 623 Lecture Notes (Date of last revision: 08/30/2021) 415
The Wong-Sandler mixing rules were reformulated by Orbey and Sandler [AIChE
J., 41, 683 (1995)] to simplify the requirements for data. The reformulated cross
term is
⎛ a ⎞ 1 ai a
⎜⎝ b − ⎟ =
RT ⎠ i 2
( bi + b ) −
RT
(1− k i )
and a modified NRTL expression is used for the excess Gibbs function
⎛ ∑ x k gki τ ki ⎞
GγE
= ∑ xi ⎜ k ⎟ where: gki = bk exp ( −σ τ ki )
RT i ⎜ ∑ x g i ⎟
⎝ ⎠
ν1 A1 + ν 2 A2 + … → … + ν m−1 Am−1 + ν m Am
∑ ν i Ai → 0
Δni ν i
=
Δn j ν j
Δni Δn j Δnk
= = =
νi νj νk
In the limit of infinitesimal changes, ∆ni ® dni. Use the above result to define the
reaction coordinate (or extent of reaction), e:
dni dn j dnk
dε = = = =
νi νj νk
dni = ν i d ε → ni − n oi = ν i ε
ni = n oi + ν i ε
nT = ∑ ni = ∑ n oi + ∑ ν i ε = nTo + ε ∑ ν i
i i i i
ni n oi + ν i ε
xi = =
nT nTo + ε ∑ ν i
i
ni = nio + ∑ ε j ν ij
j
nT = ∑ ni = ∑ nio + ∑ ∑ ε j ν ij = nTo + ∑ ε j ∑ ν ij
i i i j j i
nio + ∑ ε j ν ij
ni j
xi = = o
nT nT + ∑ ε ∑ ν ij
j
j i
(dG ) total
T ,P
=0
At equilibrium, nei = n oi + ν i εe
nG
where: ee = value of e at equilibrium
εe ε
The general form of the Physics Summary Statement for open systems with T,P
as the independent variables is:
dni = ν i d ε ⇒ ∑ µi dni = ∑ µi ν i d ε
i i
So:
⎛ ⎞
d (nG ) = − (nS ) dT + (nV ) dP + ⎜ ∑ µi ν i ⎟ d ε
⎝i ⎠
Math Relation 1:
⎡ ∂(nG ) ⎤
⎢ ∂ε ⎥ = ∑ µi ν i
⎣ ⎦T ,P
⎡ ∂(nG ) ⎤
⎢ ∂ε ⎥ = 0 ⇒ ∑ µi ν i = 0
⎣ ⎦T ,P
Also, because
µi = G i
∑ µi ν i = ∑ ν i G i = 0
dG i = RTd lnfˆi
⎛ fˆi ⎞
so that G i = G oi + RT ln ⎜ o ⎟ = G oi + RT lnα i
⎝ fi ⎠
Note: For reaction equilibrium calculations, the standard state usually is taken
to be an ideal gas at unit pressure.
∑νiG i = 0 ⎫
⎪
⎬ ⇒ ∑ ν i G oi + RT ∑ ν i lnα i = 0
G i = G oi + RT lnα i ⎪
⎭
ν νi
Note: ∑ ν i lnα i = ∑ ln (α i ) i = ln ∏ (α i )
i
3
where:
ν
∏ (α i ) i = (α1)
ν1
(α 2 )ν (α 3 )ν
2 3
i =1
K eq ≡ ∏ (α i ) i
ν
= ∏ ⎢ o ⎥ = equilibrium constant
i i
⎢⎣ f i ⎥⎦
Δr G o ≡ −RT ln K eq = ∑ ν i G oi
Notes:
• Keq, G oi are functions only of temperature
• The standard states need not be identical for every chemical species present,
but for each species, the standard state temperature must be the equilibrium
temperature and the same standard state must be used to calculate G oi , fˆi
and f io .
Definition: Δr M o = ∑ ν i M oi
Examples:
Δr Ao = Δr U o −T Δr S o
Δr U o = Δr H o − Δr (PV )
o
Δr S = −o (
d Δr G o )
dT
⎡
d ⎢ Δr G (
o
) ⎤⎥ = − Δ H r
o
dT ⎢ RT ⎥ RT 2
⎣ ⎦
(
d Δr H o )=ΔC o
r p
dT
d ⎛ Δr G o ⎞ Δr H o
⎜ RT ⎟ = −
dT ⎝ ⎠ RT 2
Δr G o
= − lnK eq ⇒
(
d ln K eq )= ΔH r
o
RT dT RT 2
K eq (T2 )
T2
Δr H o
( )
⌠
⌠
⎮
⌡
d lnK eq = ⎮
⎮ 2
dT
K eq (T1) ⌡ RT
T1
T2
⎡ K eq (T2 ) ⎤ 1 ⌠ Δr H o
to obtain ln ⎢ ⎥= ⎮ dT
⎢⎣ K eq (T1) ⎥⎦ R
⎮ 2
⌡ T
T1
⎡ T2
⌠ Δ Ho
⎤
1
K eq (T2 ) = K eq (T1 ) exp ⎢ − ⎮⎮ r 2 dT ⎥
⎢ R⌡ T ⎥
⎣ T1 ⎦
Also
T
Δr H o
(T ) = Δr H (To ) + ⌠⎮⌡ Δr CPo dT
o
To
Δr G° (T ) = Δr H ° (T ) −T ΔS° (T )
where
T
Δr H (T ) = Δr H (To ) +
o o ⌠
⎮
⌡
Δr C po dT
To
and
T
⌠ Δr C po
Δr S o (T ) = Δr S o (To ) + ⎮
⎮
dT
⌡ T
To
Consider the special case of a property change of reaction when the reaction
forms a compound from its constituent elements.
Δf M o = M substance
o
− ∑ ν i M oi
constituent
elements
Δr M o = ∑ ν i Δf M oi
Δr H o (To ) = ∑ ν i Δf H oi (To )
Δr Ao = Δr U o −T Δr S o
Δr U o = Δr H o − Δr (PV )
o
Δr S = −o (
d Δr G o )
dT
⎡
d ⎢ Δr G (
o
) ⎤⎥ = − Δ H r
o
dT ⎢ RT ⎥ RT 2
⎣ ⎦
(
d Δr H o )=ΔC o
r p
dT
Necessary information: Δf G o and Δf H o at one temperature, ideal gas heat
capacities and the equation of state. This information allows calculation of K eq
and the fugacities contained in equilibrium constant.
In this section, the term “standard state” denotes the base states for calculating
the standard free energy, standard entropy and standard enthalpy of reaction.
The term “reference state” denotes the base states for calculating the fugacity of
a species in the mixture. In many cases, the standard and reference states are
different physical states.
Ideal gas at T, P°
PURE liquid or solid at T, P°
Infinite dilution
The fugacity coefficient notation is based upon an ideal gas reference state,
where the ideal gas is at the same pressure as the system. The fugacity of
substance i in the mixture is
fˆi = y i φ̂i P
i ⎣f i ⎦ ⎦
ν ν
⎛ P ⎞
⎝P ⎠ i
ν
= ⎜ o ⎟ ∏ ( y i ) i ∏ φ̂i
i
( )
νi ⎛ P ⎞
= K y ⎜ o ⎟ ∏ φ̂i
⎝P ⎠ i
( )
νi
where
νi
K y = ∏ (yi ) and ν = ∑νi
i i
νi νi
ig ⎡ fˆi ⎤ ⎡ x i γ iLR f iLR ⎤
K eq = ∏⎢ o ⎥ = ∏⎢ ⎥
i ⎣fi ⎦ i ⎢⎣ Po ⎥⎦
νi νi
⎡ f iLR ⎤ ⎡ f (T ,P ) ⎤
=∏
i
( xi γ i )
LR ν i
∏⎢ o ⎥
i ⎢⎣ P ⎥⎦
= Kx ∏
i
( )
LR ν i
γi ∏⎢ i o ⎥
i ⎣ P ⎦
fˆi = x i γ iHL Η i (T )
Henry’s constant is based upon the behavior of the solute at infinite dilution in the
mixture
⎛ fˆi ⎞
Ηi = lim
x i →0 ⎜x ⎟
⎝ i⎠
Then
νi νi ν
⎛ fˆi ⎞ ⎛ x iγ iHL Η i ⎞ ⎛ Ηi ⎞ i
( )
ig νi
K eq = ∏⎜ o ⎟ = ∏⎜ ⎟ = ∏ ⎜ o ⎟ ∏ x i γ iHL
i ⎝ fi ⎠ i ⎝ Po ⎠ i ⎝P ⎠ i
ν
⎛ Ηi ⎞ i
( )
νi
= K x ∏ ⎜ o ⎟ ∏ γ iHL
i ⎝P ⎠ i
For liquid or solid standard states, the standard state is the fugacity of the pure
material at the standard state pressure
f io = f i P o( )
4.3.5.2.1 Ideal Gas Reference State
The fugacity coefficient notation is based upon an ideal gas reference state,
where the ideal gas is at the same pressure as the system. The fugacity of
substance i in the mixture is
fˆi = y i φ̂i P
( )
−ν i
= P ∏ ⎡f i P o ⎤
ν
i ⎣ ⎦ (
∏ y i φ̂i
i
)
νi
( )
−ν i
= P ∏ ⎡f i P o ⎤
ν
i ⎣ ⎦ i
ν
∏ ( y i ) i ∏ φ̂i
i
( )νi
( )
−ν i
= K y P ∏ ⎡f i P o ⎤
ν
i ⎣ ⎦
∏ φ̂i
i
( ) νi
ν
⎡ ⎤ i
= Ky ∏ ⎢
P ⎥ ∏ φ̂i ( ) νi
i ⎢f i P
⎣
o
( ) ⎥ i
⎦
νi
⎡ ⎤
f (P )
( )
νi
= K x ∏ γ iLR ∏⎢ i o ⎥
i i ⎢fi P
⎣ ( ) ⎥
⎦
f i (P ) ⎡ 1 P ⎤
= exp ⎢ ⌠Vi dP ⎥
( )
and ⎮
⌡
fi Po ⎢⎣ RT Po ⎥⎦
fˆi = x i γ iHL Η i (T )
Henry’s constant is based upon the behavior of the solute at infinite dilution in the
mixture
⎛ fˆi ⎞
Ηi = lim
x i →0 ⎜x ⎟
⎝ i⎠
Then
νi ν
νi ⎛ ⎞ ⎛ ⎞ i
⎛ fˆi ⎞ x iγ iHL Η i Ηi
( )
L/S νi
= ∏⎜ o ⎟ = ∏⎜ ⎟ = ∏⎜ ⎟ ∏ x i γ iHL
( ) ( )
K eq
i ⎝ fi ⎠ i ⎜⎝ f i P ⎟⎠ i ⎜⎝ f i P ⎟⎠ i
o o
ν
⎛ ⎞ i
Ηi
( )
νi
= Kx ∏⎜ ⎟ ∏ γ iHL
i ⎜⎝ f i P
o
( ) ⎟⎠ i
⎛ fˆi ⎞
f oi = Ηi = lim
x i →0 ⎜x ⎟
⎝ i⎠
fˆi = y i φ̂i P
The equilibrium constant then is
ν ν ν
⎡ fˆi ⎤ i ⎡ y i φ̂i P ⎤ i ⎛P⎞
∞
K eq = ∏⎢ o ⎥ = ∏⎢
i ⎣f i ⎦ i ⎣ Ηi ⎦
⎥ = ∏ ⎜ ⎟ ∏ y i φ̂i
i ⎝ Ηi ⎠ i
( )
νi
ν ν
⎛P⎞ ⎛P⎞
ν
= ∏ ⎜ ⎟ ∏ ( y i ) i ∏ φ̂i
i ⎝ Ηi ⎠ i i
( )
νi
( )
= K y ∏ ⎜ ⎟ ∏ φ̂i
i ⎝ Ηi ⎠ i
νi
νi νi
∞ ⎡ fˆi ⎤ ⎡ x i γ iLR f iLR ⎤
K eq = ∏⎢ o ⎥ = ∏⎢ ⎥
i ⎣fi ⎦ i ⎢⎣ Ηi ⎥⎦
νi
⎡ f iLR ⎤
=∏
i
( )
LR ν i
xi γ i ∏⎢
Η
i ⎢⎣ i ⎥⎦
⎥
νi
⎡ f (P ) ⎤
= Kx ∏
i
( )
LR ν i
γi ∏⎢ i
i ⎣ Ηi ⎦
⎥
fˆi = x i γ iHL Η i (T )
Henry’s constant is based upon the behavior of the solute at infinite dilution in the
mixture
⎛ fˆi ⎞
Ηi = lim
x i →0 ⎜x ⎟
⎝ i⎠
Then
νi νi ν
⎛ fˆi ⎞ ⎛ x iγ iHL Η i ⎞ ⎛ Ηi ⎞ i
( ) ( )
∞ νi νi
K eq = ∏⎜ o ⎟ = ∏⎜ ⎟ = ∏ ⎜ ⎟ ∏ x i γ iHL = K x ∏ γ iHL
i ⎝ fi ⎠ i ⎝ Ηi ⎠ i ⎝ Ηi ⎠ i i
ν
∞ ⎛ fˆi ⎞ i equilibrium constant based upon infinite dilution
K eq = ∏⎜ ⎟
i ⎝ Ηi ⎠
standard state
νi
K P = ∏ (Pi ) Pi = y i P = partial pressure of specie i
i
νi
K y = ∏ (yi ) y i = mole fraction of specie i in the vapor phase
i
νi
K x = ∏ (xi ) x i = mole fraction of specie i in the liquid phase
i
K c = ∏ (c i )
νi c i = concentration of specie i (usually in the liquid
i phase)
( )
νi
ig
K eq ⎛ f Po ⎞
i
L/S
= ∏⎜ o
⎟
K eq i ⎜⎝ P ⎟⎠
ig νi
K eq ⎛Η ⎞
∞
= ∏ ⎜ oi ⎟
K eq i ⎝P ⎠
( )
νi
∞
K eq ⎛ f Po ⎞
i
= ∏⎜ ⎟
i ⎜⎝ Η i
L/S
K eq ⎟⎠
ν νi ν νi
⎡ fˆi ⎤ i ⎡ fˆi ⎤ ⎡ f io ⎤ i ⎡ f io ⎤
K p = ∏ [Pi ]
νi ig
= ∏⎢ V ⎥ = ∏⎢ o ⎥ ⎢ V ⎥ = K eq ∏⎢ V ⎥
i i ⎢⎣ φ̂ i ⎥⎦ i ⎢⎣ f i ⎥⎦ ⎢⎣ φ̂ i ⎥⎦ i ⎢⎣ φ̂ i ⎥⎦
ig
Note: K eq is dimensionless but KP has dimensions of pressure if n ≠ 0.
ν νi νi
⎡ fˆi ⎤ i ⎡ fˆi ⎤ ⎡ f io ⎤ ⎡ Po ⎤
K y = ∏ [yi ]
νi ig
= ∏⎢ V ⎥ = ∏⎢ o ⎥ ⎢ V ⎥ = K eq ∏ ⎢ V ⎥
i i ⎢⎣ φ̂ i P ⎥⎦ i ⎢⎣ f i ⎥⎦ ⎢⎣ φ̂ i P ⎥⎦ i ⎢⎣ φ̂ i P ⎥⎦
νi
Kx ⎛ φ̂i P ⎞
= ∏⎜ ⎟
Ky i ⎝ γ if i ⎠
Note: Both KP and Ky are functions of pressure through φ̂Vi , and Ky is directly a
function of pressure if n ≠ 0.
νi
⎛ f io ⎞
( )
νi νi
K c = ∏ (c i ) = ∏ yi c = K eq ∏ ⎜ v ⎟
i i i ⎝ φ̂i cP ⎠
Molarity
Concentration
Consider the case where m chemical species are present at fixed temperature
and pressure. The total Gibbs energy for this system is a function of the
amounts of each species present:
When reactions occur, the number of atoms of each element present must be
conserved even though the atoms may be regrouped through the reactions to
form different compounds.
The mass balance for each atom provides an additional constraint for the
minimization problem:
m
A j − ∑ ni aij = 0
i =1
fˆi y i φ̂Vi P
Vapor phase: = (Po = 1 in most cases)
f io Po
L
fˆi L f i
Liquid phase: = xi γ i ≅ xi γ i
f io f io
fˆi f is
Solid phase: = ≅1 (pure solids)
f io f io
CHEN 623 Lecture Notes (Date of last revision: 08/30/2021) 451
4.3.6.1 Method of Lagrange multipliers
g = g ( x1,x 2 ,,x m )
The minimum in a function, h = h(y1, y2, ..., ym), where the {yi} all are
independent, occurs when:
∂h
=0 i = 1, , m
∂y i
∂2 h ⎧⎪ i = 1, , m
>0 ⎨
∂y i ∂y j j = 1, , m
⎩⎪
simultaneously.
CHEN 623 Lecture Notes (Date of last revision: 08/30/2021) 452
In the constrained problem, all m variables {xi} are not independent. One
solution is to eliminate w variables algebraically, then solve the case where the
derivatives of g with respect to the remaining m - w variables all are zero.
However, this procedure is not suited well to general applications.
Constraints:
w
Also: ∑ λk f k = 0
k =1
w
Define: F = g + ∑ λk f k
k =1
(m + w )variables ⎫
⎬ ⇒ m independent variables
w constraints ⎭
Now we can choose {xi} to be the independent variables. Then, the minimum
values of F and hence g, occur when:
⎛ ∂F ⎞
⎜⎝ ∂x ⎟⎠ =0 for i = 1, 2, , m
i x j ≠i
The values {xi} that yield the minimum for g are given by the solution to the
following set of simultaneous equations:
⎛ ∂F ⎞
⎜⎝ ∂x ⎟⎠ =0 for i = 1, 2, , m
i x j ≠i
fk = 0 k = 1, 2, !, w
( )
m m
Aj − ∑ ni aij = ∑ nio − ni aij = 0 ( j = 1, 2,!,w )
i =1 i =1
⎡ mw ⎤
F = nG − ∑ λ j ⎢ A j − ∑ nk akj ⎥
j =1 ⎣ k =1 ⎦
So:
⎛ ∂F ⎞ ⎡ ∂(nG ) ⎤ w w
† fˆi w
⎜⎝ ∂n ⎟⎠ =⎢ ⎥ + ∑ λ j aij = G i + ∑ λ j aij = G i + RT ln † + ∑ λ j aij
i n ⎣ ∂ni ⎦n j ≠i j =1 j =1 fi j =1
j ≠i
∂F o fˆi w
= Δf G i + RT ln o + ∑ λ j aij
∂ni fi j =1
fˆi w
Δf G oi + RT ln + ∑ λ j aij = 0 (i = 1, 2, , m)
f io j =1
m
A j − ∑ ni aij = 0 ( j = 1, 2,,w )
i =1
Useful Reference:
Each reaction provides one equation with one unknown (the reaction coordinate)
for each independent chemical reaction present.
For the case of a single reaction, the mole fraction is related to the reaction
coordinate through stoichiometry by
n oi + ν i ε
xi =
∑ n oi + ε ∑ ν i
i i
νi
⎡ o
⎤
⎢ n i + νi ε ⎥
Kx = ∏ ⎢ o ⎥
i ∑ n i + ε ∑νi
⎢⎣ i i ⎥⎦
ν
⎡
⎢
Ky = ∏ ⎢ o
n o
(
i + νi ε ) ⎤i ⎧
⎥ ⎪ ⎡ f o ⎤ν i ⎫
⎪ ⎡ Δ G o⎤
⎥ = ⎨∏ ⎥ ⎬ exp ⎢ −
i r
⎢ ⎥
i ∑ n i + ε ∑νi φ̂
⎪⎩ i ⎢⎣ i ⎥⎦ ⎪⎭
P ⎢⎣ RT ⎥⎦
⎢⎣ i i ⎥⎦
This gives one equation, with the single unknown e. The equation is implicit
because the φ̂i are dependent upon {yi}.
ν
⎡
⎢
Kx = ∏ ⎢ o
n o
(
i + νi ε ) ⎤i ⎧
⎥ ⎪ ⎡ fio ⎤ i⎫
ν
⎪ ⎡ Δr G o ⎤
i ∑ n i + ε ∑νi
⎥ = ⎨∏ ⎢ † † ⎥ ⎬ exp ⎢ − ⎥
⎢⎣ i ⎥⎦ ⎪⎩ i ⎢⎣ γ i f i ⎥⎦ ⎪⎭ ⎢⎣ RT ⎥⎦
i
This gives one equation, with the single unknown e. The equation is implicit
because the γ i† are dependent upon {xi}.
CHEN 623 Lecture Notes (Date of last revision: 08/30/2021) 459
4.3.7.2 Multiple Independent Reactions
n oi + ∑ ν ik ε k
xi = k
⎛ ⎞
∑ n oi + ∑ ⎜ ε k ∑ ν ik ⎟
i k ⎝ i ⎠
where the subscript i denotes the individual chemical species (m total) and the
subscript k denotes the independent chemical reactions (r total).
νi
⎡ o ⎤
⎢ + ∑ ν ε ⎥
n i ik k (i = 1,2,,m )
Kx = ∏ ⎢
k k ⎥
i ⎢
∑ n o
+ ∑
⎛
ε ∑ ν
⎞⎥ (k = 1,2,,r )
⎢ i ⎜⎝ k ik ⎟⎠ ⎥
⎣ i k i ⎦
CHEN 623 Lecture Notes (Date of last revision: 08/30/2021) 460
Using the fugacity coefficient notation or the activity coefficient notation (see
εi
single reaction section) then provides a system of simultaneous equations in the
reaction coordinate variables, .
( ) (
Hi (T ,P ) = ⎡⎣Hi (T ,P ) − Hio (T ,P ) ⎤⎦ + ⎡Hio (T ,P ) − Hio T o ,P o ⎤ + Hio T o ,P o
⎣ ⎦ )
CHEN 623 Lecture Notes (Date of last revision: 08/30/2021) 461
When the enthalpies are expressed relative to the free elements in their standard
states
⎧⎪0
(T )
for free elements
Hio o
,P o
=⎨ o
⎪⎩ Δ fHi for compounds
T P
⌠
⎛ ∂Hio ⎞
(T )=
⌠
Hio (T ,P ) − Hio o
,P o ⎮
⎮
⎮
o
C p,i dT +
⎮
⎮ ⎜ ∂P ⎟ dP
⌡ ⎮
⌡ ⎝ ⎠T
To Po
P
⌠
⎮ ⎛ ∂Hio ⎞ ⎧= 0 ideal gas standard states
where ⎮ ⎜ ∂P ⎟ dP ⎨≅ 0
⎮
⌡ ⎝ ⎠T ⎩ liquid or solid standard states
Po
m
H total
(T ,P ) = nΔm H (T ,P ) + ∑ ni Hi (T ,P )
i =1
⎧ T
⎪ Δf Hio + ⌠⎮⌡ C p,i
ig
dT + H R (T ,P ) ideal gas standard states
⎪ To
⎪
where H i (T ,P ) = ⎨
⎪ T
⎪ Δ H + ⌠⎮ C /s dT
o
liquid or solid standard states
⎪ f i ⌡ p,i
⎩ To
4.4.1 Molality
To describe the behavior of ionic solutions, compositions often are expressed
using molalities.
Definition: Molality is the number of moles of solute per kilogram of solvent.
Notation: ni = number of moles of solute (i)
no = number of moles of solvent.
Mo = molar mass of solvent
mi = molality of solute (i)
Then:
1000ni ni
mi ≡ and xi =
noM o no + ∑ nk
k
so that
⎛ ⎞
no + ∑ nk 1000 ⎜ no + ∑ nk ⎟
mi ni (1000 ) ⎝ k ⎠ 1000
= k = =
xi noM o ni noM o x oM o
Volume Concentration
ni
Definition: ci ≡
V
Then
noM o + ∑ nk M
V= k
ρ
ni ρ ρx
ci = ≅ i (dilute solution)
noM o + ∑ nk M Mo
k
The properties of the various species in the solution usually are described as
follows
⎛ fˆi ⎞
Solute: G i = G ‡i + RT ln ⎜ ‡ ⎟
⎝fi ⎠
⎛ fˆo ⎞
Solvent: Go = Go† + RT ln ⎜ † ⎟
⎝ fo ⎠
Solute: fˆi = γ i m i f i‡
lim ⎛ fˆi ⎞
Solute: f i‡ = mi →0 ⎜m ⎟
⎝ i⎠
⎛ fˆi ⎞
f i◊
lim
Solute: = ci →0 ⎜c ⎟
⎝ i⎠
Form of fugacity
Convention
Mixture Solutes Solvents
LR/LR x i γ i f i† f i† = f iL f i† = f iL
lim ⎛ fˆi ⎞
HL1/LR x i γ i f i† f i† = x i →0 ⎜ x ⎟ = Ηi f i† = f iL
⎝ i⎠
x oγ ofo† ⎛ fˆi ⎞
f i◊
lim
HL3/LR = ⎜c ⎟ f o† = foL
c i γ ◊i f i◊
ci →0
⎝ i⎠
Vapor:
⎛ fˆi ⎞ ⎛ P ⎞
GvI = Gio + RT ln ⎜ o ⎟ = Gio + RT ln ⎜ y i φ̂i o ⎟
⎝ fi ⎠ ⎝ P ⎠
where: P o is the standard state pressure (at which Gio is specified), and
usually is at unit pressure.
Liquid Solution:
⎛ fˆi ⎞ ⎛ m ⎞
GLI = Gi† + RT ln ⎜ † ⎟ = Gi‡ + RT ln ⎜ γ i oi ⎟
⎝ fi ⎠ ⎝ mi ⎠
‡
where: m oi is the standard state molality (at which Gi is specified), which
usually is unit molality.
⎛ P ⎞ ⎛ mi ⎞
Gio ‡
+ RT ln ⎜ y i φ̂i o ⎟ = Gi + RT ln ⎜ γ i o ⎟
⎝ P ⎠ ⎝ mi ⎠
or
⎡ y φ̂ P ⎛ m o ⎞ ⎤
Gi‡ − Gio = RT ln ⎢ i i ⎜ oi ⎟ ⎥
⎢⎣ γ i m i ⎝ P ⎠ ⎥⎦
Use a dilute solution and a pure vapor to evaluate the numerical value of the
difference between the standard states.
Notation: Mν + Aν − = ν +Mz + + ν − Az −
Example
H2SO 4 → 2H + + SO =4
ν+ = 2 ν− = 1
z+ = 1 z − = −2
Notation:
m + = ν + (m − mu )
m − = ν − (m − mu )
mu >> m+ ,m−
mu << m ⇒ mu ≈ 0
The Gibbs energy relative to the infinite dilution standard state is related to the
Gibbs energy relative to ideal gas standard state by
⎡ y φ̂ P ⎛ m o ⎞ ⎤
Gi‡ − Gio = RT ln ⎢ i i ⎜ oi ⎟ ⎥
⎢⎣ γ i m i ⎝ P ⎠ ⎥⎦
⎧
⎪ lim
⎡ P ⎛ mo ⎞ ⎤ ⎫ ⎪
Gi‡ − Gio = RT ln ⎨ mi →0 ⎢ ⎜ oi ⎟ ⎥ ⎬
⎪
⎩ ⎢⎣ m i ⎝ P ⎠ ⎥⎦ ⎭
⎪
Using the above results, the free energy of formation in solution is defined as
⎧
⎪ ⎡ fˆ ⎛ m o ⎞ ⎤ ⎫ ⎪
Δf Gi‡ (solution) = Δf Gio + RT ln ⎨
lim
mi →0 ⎢ i ⎜ oi ⎟ ⎥ ⎬
⎪
⎩ ⎢⎣ γ i m i ⎝ P ⎠ ⎥⎦ ⎪
⎭
This type of approach often is used for solutions of electrolytes in water, e.g.
salts, inorganic acids, etc.
Using the notation defined for electrolyte solutions changes this equation to
dm + = ν + (dm − dmu )
dm − = ν − (dm − dmu )
Then
dGtot = −S tot dT +V tot dP + [ µu − ν + µ+ − ν −m− ]dmu + [ν + µ+ + ν − µ− ]dm + µo dno
In the usual restriction to constant temperature and pressure, the first two terms
on the right hand side become zero, so that
dGtot = [ µu − ν + µ+ − ν − µ− ]dmu
⎡ ∂Gtotal ⎤
⎢ ⎥ =0
⎢⎣ ∂mu ⎥⎦T ,P,m,n
o
= µu − ν + µ+ − ν − µ−
∑ ν i Gi = 0
i
where ν i = 1 for undissociated electrolyte. Note that this form is exactly the same
as that for chemical equilibrium.
∑ ν i G i = 0 ⇒ Gu = ν +G+ + ν −G−
i
Nomenclature
ν = ν+ + ν−
Then
fˆov = y oφ̂ov P ≅ P ⎫
⎪ P
⎬ ⇒ γ o = σ σ
fˆoL = x oγ ofoL ≅ x oγ oφoσ Poσ ⎪
⎭ x oφo Po
Solvent
The activity coefficient of the solvent also can be determined by measuring the
freezing point depression. For aqueous solutions, if the solid formed is pure ice
(usually the case for ice, the solvent fugacity is
fˆoL = x oγ ofoL ⎫⎪ 1
⎬ ⇒ γ o =
ˆf s = f s = f L ⎪ xo
o o o ⎭
Note that as freezing occurs, solvent is removed from the solution resulting in an
increasing solute concentration. Therefore, the above relation must be applied to
the conditions when the amount of ice goes to zero.
Also:
Gu = Gu‡ + RT ln (mu γ u )
G+ = G+‡ + RT ln (m +γ + )
G− = G−‡ + RT ln (m − γ − )
⎡ (m γ )ν ⎤
ν +G+‡ + ν −G−‡ − Gu‡ = −RT ln ⎢ ± ± ⎥ = −RT lnK e
⎢⎣ muγ ⎥⎦
where
ν
Ke =
( m± γ ± )
=
γ ν±mν++ mν−−
= true dissociation constant
muγ u muγ u
Then
z+
for the dissociation reaction: Mν + Aν − = ν +M + ν − Az −
γ Bγ D mB mD γ B2γ D mB mD γ BD
3
mB mD
Equilibrium for Reaction 2: K2 = = =
γ C mC γ Bγ C mC 2
γ BC mC
2 3
where γ BC = γ Bγ C , γ BD = γ B2γ D
γ B2 γ D mB2 mD γ BD
3
mB2 mD
Equilibrium for Reaction 3: K3 = =
γ A mA γ A mA
Note: K 3 = K 1K 2
CHEN 623 Lecture Notes (Date of last revision: 08/30/2021) 483
Example 2:
γ C2 γ D mC2 mD γ CD
3
mC2 mD
Ke = 2 2
= 3 2
γ Aγ B mAmB γ AB mAmB
Temperature Dependence of Ke
∂lnK e ν +H +∞ + ν −H −∞ − H ∞
=
∂T RT 2
Example:
1
H2 ( g ) + O2 ( g ) → H2O ( ) o
Δf GH
2 2O( )
3 1
H2 ( g ) + N2 ( g ) → NH3 (aq ) ◊
Δf GNH
2 2 3 (aq )
o ‡ ◊ o
Δf GNH = Δ G − Δf GNH − Δ G
, OH (aq ) f OH (aq ) 3 (aq ) f H2O( )
+ -
4 NH+4 , -
Also
Δf HHo+ = 0 = Δf SHo+
Galvanic Cells
= -zFE
where
Examples:
1
Platinum electrode: H → H+ + e-
2
1
Mercury electrode: Hg2Cl2 + e- → Hg ( ) + Cl-
2
1 1
Total reaction: H2 + Hg2Cl2 → HCl + Hg ( )
2 2
Then
1 1
−FE = GHCl + GHg − G H 2 − G Hg 2Cl2
2 2
(The previous equation assumes that HCl is a strong electrolyte, i.e., that
essentially all HCl is dissociated.)
Then
1 s 1 o
−FE = G ‡+ + G ‡ - + GHg − GHg − GH2 − RT ⎡ln PH2 − ln (γ ±m ± ) ⎤
H Cl 2 2 Cl2
2 ⎣ ⎦
or
⎛γ m ⎞
−FE = −FE + RT ln ⎜ ± ± ⎟
o
⎜⎝ PH ⎟⎠
2
1 s 1 o
where: −FE o = G ‡+ + G ‡ - + GHg − GHg2Cl2 − GH = Δr G ‡
H Cl 2 2 2
Note
lim
m→0 γ± =1
lim
m→0 ⎡⎣2RT ln (m ) + F E ⎤⎦ = F E o
Using this definition, various reactions can be used to generate tables of single
electrode potentials for various ions.
Note: Δf G ‡ = −zFE o
1
Ι= ∑ mi z i2
2 all ions
in solution
Debye and Huckel derived a theoretical expression for the behavior of the ionic
activity coefficient in the infinite dilution limit:
lnγ ± = −1.174 z + z − Ι
Ι
lnγ ± = −1.174 z + z −
1+ Ι
GE 1 1
= w f ( Ι ) + ∑ ∑ λik ( Ι ) ni nk + 2 ∑ ∑ ∑ λik ni nk n
RT w i k w i k
where w is the mass of water present, and the l's are osmotic virial
coefficients. The second osmotic virial coefficient is treated as a function of
ionic strength, while the third osmotic virial is assumed not to be affected by ionic
strength.
2 2 (ν +ν − )
⎡ 2 ν +ν − ⎤ γ ⎡ 3/ 2 ⎤
γ γ
lnγ ± = z +z − f + m ⎢ B + m ⎢ ⎥ C
⎣ ν ⎥⎦
± ±
⎢⎣ ν ⎥⎦
where:
⎡ ⎤
f γ = −3Aφ ⎢
Ι 2
+ ln 1+ b Ι ⎥
⎣ 1+ b Ι b ⎦
( )
3 φ
C ±γ = C
2 ±
3/ 2
1 ⎡ e2 ⎤
Aφ = 2π (1000 )N A ρs ⎢ ⎥
3 ⎢⎣ 4πε oDs kT ⎥⎦
B ±γ = 2 βo +
2 β1 ⎡
2 ⎣
α Ι
1− 1+ α Ι −(0.5α 2
) (
I exp −α Ι ⎤
⎦ )
α = 2.0 (kg/mol)
1/2
β = 1.2 (kg/mol)
1/ 2
2 βi ⎡
( ) ( )
2
B ±γ = 2 βo + ∑ 2 ⎣
1− 1+ α Ι − 0.5α 2
Ι exp −α i Ι ⎤
i =1 α i Ι
i i ⎦
α1 = 1.4 (kg/mol)
1/2
α 2 = 12 (kg/mol)
1/2
Γ oij = (γ ± )
1/ z + z −
2 1 2
1 1 1 1
2 1 2 2 1 2
Bulk Surface
For the two cases, the contributions of the intermolecular interactions between
the molecules to the energy of the system are
1
E bulk =
2
( 4ε 1 + 4ε 2 ) < 0
1
E surface =
2
( 3ε1 + 2ε 2 ) < 0
Here ε1 is the average interaction energy for nearest neighbors, and ε 2 is the
average interaction energy for second nearest neighbors. (The factor of 1/2
avoids double counting).
1 1
δ = E surface − E bulk = ⎡⎣ 3ε1 + 2ε 2 − ( 4ε1 + 4ε 2 ) ⎤⎦ = − ( ε1 + 2ε 2 ) ⇒ δ > 0
2 2
The net result is that a molecule in the surface layer contributes a negative
contribution of smaller magnitude to the energy than does a molecule in the
interior.
As a result, the total energy differs from the result calculated by assuming that all
molecules have the same energy as an interior molecule by a small positive
amount
E total = ∑ ni µi + nsurface δ
i
Here nsurface is the number of moles that differ from the bulk, and δ is the
average energy deviation. Clearly, nsurface will vary directly with the total surface
area.
Position
It is most convenient to use the total area of the surface as the independent
variable to describe this effect, and the variation of energy with area is the
surface tension, which is defined thermodynamically by
For the case where surface effects are significant, the physics summary
statement becomes
/ rev = −P dV + σ da
dW
β
β β
α - phase β
( ) ( ) ( )
d nαU α = T α d nα S α − P α d nαV α + ∑ µiα dniα
d (n βU β ) = T β d (n β S β ) − P β d (n βV β ) + σ d (a β ) + ∑ µ β dn β i i
(
d nα S α = ) T
1 ⎡
α ⎣ ( ) (
d nαU α + P α d nαV α − ∑ µiα dniα ⎤
⎦ )
(
d n βS β = ) 1 ⎡
Tβ ⎣
d n β β
U (
+ P β
d n β β
V − )
σ da β
− (
∑ µi
β
)
dn β⎤
i ⎦
For the case where a and b together form a closed isolated system:
( ) ( )
d (nS ) = d nα S α + d n β S β
⎡ µα µ β ⎤ α
= ⎢ α − β ⎥ d (n U ) − ∑ ⎢ α − β ⎥dn
⎡ 1 1 ⎤ α α i i
i
⎣T T ⎦ ⎢T T ⎥ ⎣ ⎦
⎡ α β ⎛ β ⎞⎤
⎢T
P P
T
σ
+⎢ α − β + β ⎜
da
(
T ⎜⎝ d n βV β
(
⎟ ⎥ d n βV β
)
⎟⎠ ⎥
)
⎣ ⎦
Therefore all terms in square brackets must be identically equal to zero, and the
equilibrium conditions are:
T α =T β
µiα = µiβ
α da β
β
P = P −σ
d (n βV β )
The equilibrium conditions are the same as those derived for the case without
surface effects, except that the pressures are not equal in the coexisting phases.
µiα = µiβ ⇒ f iα = f iβ
Then: f iα = φ αi P α
f iβ = f iβ P β( ) where P β ≠ Pα
⎡ 1 Pβ ⎤
and: f iβ σ σ
= φ i Pi exp ⎢ β
∫ V dP ⎥
⎢⎣ RT Piσ ⎥⎦
φ σi Piσ ⎡ 1 Pβ β ⎤ ⎡V β σ ⎤
Pα =
φ αi
exp ⎢
⎢⎣ RT
∫
Pi σ
V dP ⎥
⎥⎦
≅ P σ
i exp ⎢
⎢⎣ RT
P (
β
− Pi ⎥)
⎥⎦
⎡ β ⎛ ⎞⎤
da β V da β
P β = Pα + σ ⇒ P α = Piσ exp ⎢ ⎜Pα + σ − σ ⎟⎥
( ) ( )
P
d n V β β ⎢ RT ⎜⎝ d n βV β
i
⎟⎠ ⎥
⎣ ⎦
α da β
− Piσ << σ
( )
P
d n βV β
so that
⎡ β ⎛ ⎞⎤
σV da β
α
≅ Piσ exp ⎢ ⎜ ⎟⎥
( )
P
⎢ RT ⎜⎝ d n βV β ⎟⎠ ⎥
⎣ ⎦
This result also provides a good initial guess for the iterative solution of the
general equation.
β 2 da β
a = 4π r ⇒ = 8π r
dr
β β 4 3
n V = πr ⇒
d n β β
V
= 4π r 2
( )
3 dr
da β ⎛ da β ⎞ ⎛ dr
⎞
8π r 2
=⎜ ⎜ ⎟= =
(
d n βV β ) ⎟
⎝ dr ⎠ ⎜⎝ d n βV β ( ) ⎟⎠ 4π r 2 r
so that
α ⎡ σV β 2 ⎤ ⎡ 2σV β ⎛ 1⎞ ⎤
P ≅ Piσ exp ⎢ σ
⎥ = Pi exp ⎢ ⎜⎝ ⎟⎠ ⎥
⎢⎣ RT r ⎥⎦ ⎢⎣ RT r ⎥⎦
and
2σ
Pα = P β −
r
CHEN 623 Lecture Notes (Date of last revision: 08/30/2021) 508
5.4.2 Case II: Small Vapor Bubbles Surrounded By a Liquid (Pure
Fluid)
The derivation for this case is similar to that of Case I, except that here the b
phase is a vapor and the a phase is a liquid. Then
⎡ 1 Pα ⎤
f iα σ σ
= φ i Pi exp ⎢ α
∫ V dP ⎥
⎢⎣ RT Piσ ⎥⎦
f iβ = φi P β
φ σi ⎡
P β = β P iσ exp ⎢
1 Pα
∫ V α
dP
⎤
⎥ ≅ P σ
exp ⎢
(
⎡V α P α − P σ
i ) ⎤⎥
i
φi ⎢⎣ RT Piσ ⎥⎦ ⎢ RT ⎥
⎣ ⎦
⎡ α ⎛ ⎞⎤
V da β
P β ≅ Piσ exp ⎢ β σ
⎜ P − Pi − σ ⎟⎥
⎢ RT
⎣
⎜⎝ d n βV β ( ) ⎟⎠ ⎥
⎦
da β 2
=
( )
As in the example for Case I,
d n βV β r
Then
β ⎡ 2σV α ⎛ 1⎞ ⎤
P ≅ Piσ exp ⎢ − ⎜⎝ ⎟⎠ ⎥
⎢⎣ RT r ⎥⎦
and
2σ
Pα = P β −
r
Pitzer and Curl developed a correlation for surface tensions that requires
knowledge of the liquid molar volume at one temperature, the critical temperature
and the acentric factor.
σ
= (1−Tr ) 9
11
σ0
−1
2/ 3 ⎡ ⎛ 3 ⎞ 2/3 ⎤
σ oVo ⎢ dyne cm 1⎥
= 1.86 + 1.18 ω
Tc ⎢ cm ⎜⎝ mol ⎟⎠ K ⎥
⎣ ⎦
This correlation generally provides values within 5% for "normal" fluids. It does
not provide accurate estimates for hydrogen-bonding fluids such as water,
alcohol, etc.
The attractive forces arising from interactions with the molecules of a solid
surface leads to an increased density of the fluid near the wall, as shown
schematically below.
Region of Adsorption
Density
Potential Energy
Because the adsorption region is associated with the surface, and it is only a few
nanometers in thickness, it is convenient to treat this region as a separate phase
that exists in two dimensions rather than three.
5.5.1 Terminology
Adsorption refers to the capturing of gas molecules by a solid surface. There are
two classes of adsorption: chemisorption and physical adsorption. In the case of
chemisorption, the molecules are bound to specific sites on the surface, as
shown below, due to special chemical properties associated with the site. In this
case the molecules have little mobility along the surface
Solid Surface
Chemisorption
Solid Surface
Physical Adsorption
This gives rise to fluid behavior in two dimensions rather than three. Such fluids
exhibit phase transitions between vapor and liquid and between liquids and
solids as the surface temperature is lowered. Depending upon how much gas is
adsorbed on the surface, the adsorption falls into one of the following classes:
Solid Surface
Solid Surface
Submonolayer Adsorption
Solid Surface
Solid Surface
Monolayer Adsorption
Monolayer Adsorption
Solid Surface
Multilayer Adsorption
Special structures are required to obtain adequate surface areas. Materials such
as activated charcoal, zeolites, etc. can provide 2000 - 3000 m2 of surface area
per gram of material.
where
a
Θ= = molar area (two dimensional analog of molar volume
η
⎛ ∂(ηU ) ⎞ ⎛ ∂(ηA ) ⎞
π =⎜ =⎜ = spreading pressure
⎝ ∂(ηΘ ) ⎟⎠ ηS,ηk
⎝ ∂(ηΘ ) ⎟⎠ T ,ηk
(ηΘ ) d π = ∑ ηi d µis = ∑ ηi dG is
i i
P
RT ⌠ dP
so that π (P ) = ∑ ⎮⎮ ηi (P ) Z bi (P )
(ηΘ ) i ⌡ 0
P
P
RT ⌠ dP
For a pure fluid, the previous result reduces to π (P ) = ⎮ η i (P ) Z i (P )
b
(ηΘ ) ⎮⌡
0
P
η (P ) = kP
This result applies in the low pressure limit, where the surface density of
adsorbed molecules is low.
RT ⎡ BP ⎤ RT ⎡ BP ⎤ RT
π HL = 1+ ( k P ) = 1+ ( k P ) ≅ (k P )
(ηΘ ) ⎢⎣ RT ⎥⎦ a ⎢⎣ RT ⎥⎦ a
BP
Z bi = 1+
RT
bkP
η (P ) =
b +P
mP
η (P ) =
(b + P ) t 1/t
⎛ ∂P ⎞ HV − H ads Δav H
= =
⎜⎝ ⎟
∂T ⎠ n ads T V V −V ads (T ΔavV )
At low pressures, where the vapor phase may be considered ideal, and the
volume of the vapor is much greater than the volume of the adsorbate, the
Clapeyron equation becomes:
CHEN 623 Lecture Notes (Date of last revision: 08/30/2021) 524
⎛ ∂lnP ⎞ Δav H
⎜⎝ ⎟ =
∂T ⎠ nads RT 2
The expression for the total internal energy of a collection of particles that was
derived in the discussion of the First Law of Thermodynamics (energy balance) is
! !
total
Usystem = 21 ∑ mi Ui • Ui (motion: translation, rotation, vibration)
i
!
ri
! ! !
1
i k r!
⌠
+ 2 ∑ ∑ ⎮⌡ fik Ri − Rk ( ) !
• d ri (intermolecular interactions )
ref
!
ri
!
⌠ ext !
+ ∑ ⎮⌡ F i ri • d ri
i r!
!
( ) (external field interactions )
cm
In the previous treatments, the effect of the external field has been assumed to
be sufficiently small that the third term can be taken as zero.
The third term involves a summation over all molecules present, so it is not
tractable for calculations in this form.
CHEN 623 Lecture Notes (Date of last revision: 08/30/2021) 526
In the third term, the integral is the potential energy of particle i relative to the
position of the center of mass that is due to the external field. To simplify the
notation, X will denote the potential energy contribution resulting from the
interaction of one mole of substance i with the external field
!
!
ri
! ext ! !
Ξi ( ri ) = ⌠
⎮
!⌡
( )
Fi ri • d ri
rcm
The choice of which coordinate system to use usually is based on the symmetry
of the physical situations.
Ξδi V = niδV Ξi
and the contribution of substance i to the internal energy of the system then is
and the total internal energy due to the presence of the field is
nΞ = ∑ ni Ξi
all
substances
( )
d (nU ) = Td (nS ) − Pd (nV ) + ∑ ni dΞi + ∑ µinf + Ξi dn i
i i
The chemical potential of substance i in the mixture with the field present is the
coefficient of dni in the physics summary statements:
µifield = µino field + Ξi
⎛ ⎞ ! ! ! ! ! ! ! !
so that ∑ i i
n dΞ =
⎝
∑ i i
n M
⎠
g i dr = ( )
nM g i dr = mg i dr = ρ m( nV ) i dr
g
i i
so that
! ! ∑ xi M i !
∇P = − ρm g = − i g
V
Again, the gradient lies along the direction of the external field. The simplest
treatment is to define a coordinate axis z in either cylindrical or Cartesian
coordinates that lies in the direction of the gravitational field. Then
Ξi = M i g ( z − z o )
and
⎡M g (z − z o ) ⎤
fˆi field = fˆi no field exp ⎢ i ⎥
⎣ RT ⎦
(
⎛ ∂P ⎞ ⎛ ∂ mω r ⎞ ∑i x i M i 2
2
)
⎜⎝ ∂ r ⎟⎠ = ⎜ ∂(nV ) ⎟ = V ω r
⎝ ⎠
⎛ ∂ µifield ⎞ ∂ ⎡⎛ ⎞rω 2 ⎤ = −M rω 2
⎜⎝ ∂r ⎟⎠ = − ∑ n M
∂ni ⎢⎣⎝ k ⎠ ⎥⎦
k k i
Again, the gradient lies along the direction of the external field. The simplest
treatment is to define a coordinate axis z in either cylindrical or Cartesian
coordinates that lies in the direction of the gravitational field. Then
1
Ξi = − M i ω 2 r 2 − ro2
2
( )
and
fˆi field
= −fˆi no field
exp ⎢
(
⎡ M i ω 2 r 2 − ro2 ) ⎤⎥
⎢⎣ 2RT ⎥⎦