MSCCH 503
MSCCH 503
MSCCH 503
M. Sc. I Semester
PHYSICAL CHEMISTRY-I
SCHOOL OF SCIENCES
DEPARTMENT OF CHEMISTRY
UTTARAKHAND OPEN UNIVERSITY
MSCCH-503
PHYSICAL CHEMISTRY-I
SCHOOL OF SCIENCES
DEPARTMENT OF CHEMISTRY
UTTARAKHAND OPEN UNIVERSITY
Board of Studies
Prof. A.B. Melkani Prof. G.C. Shah
Department of Chemistry Department of Chemistry
DSB Campus, Kumaun University SSJ Campus, Kumaun University
Nainital Nainital
Programme Coordinator
Course Editor
Dr. M.G.H. Zaidi
Professor & Head,
Department of Chemistry,
College of Basic Sciences and Humanities
G.B. Pant University, Pantnager
BLOCK- 1 THERMODYNAMICS
BLOCK-3 ELECTROCHEMISTRY
CONTENTS:
1.1 Introduction
1.2 Reversible and irreversible process
1.2.1 Reversible process
1.2.2 Irreversible process
1.3 Thermodynamic work
1.4 Work Reversible Cycle Process
1.5 First law of thermodynamics
1.6 Internal energy
1.7 Mathematical formulation of first law of thermodynamics
1.8 Enthalpy
1.9 Calculation of Cp-Cv Thermodynamically
1.10 Joule’s law
1.11 Joule Thomson Effect
1.12 Joule-Thomson coefficient in ideal gas
1.13 Joule-Thomson coefficient for a real gas
1.14 Second Law of Thermodynamics
1.15 Cyclic Process
1.16 Entropy: Nature of dq and dq/T
1.17 Physical Significance of Entropy
1.18 Entropy of a Perfect Gas
1.19 Change of entropy with temperature
1.20 Entropy of mixing
1.21 Work function or Helmholtz free energy A and Gibbs free energy or Gibbs potential G
1.22 Change in A and G with the parameters of the system
1.23 Maxwell’s relations
1.1 INTRODUCTION
Thermodynamics is the branch of chemistry that deals with energy transformation in all type
of physical chemical and biological processes. Thermodynamics leads with three
fundamental laws which are purely based on human experiences. These are the first, second
and third law of thermodynamics. These laws provide numerical approach of prediction of
spontaneity and feasibility of various processes under prescribed set of conditions.
Thermodynamic principles allow deriving most of the important laws of the physical
chemistry such as vant Hoff’s law, phase rule and distribution law.
Thermodynamics is applicable to the microscopic system consisting of matter in bulk and not
to microscopic systems of individual atoms or molecules. It does not bother about the time
factor. That means it does not tell anything regarding the rate of a physical change or a
chemical reaction. It is concerned only with initial and final states of the system.
Thermodynamic principles correlate different observable properties of system without
providing their magnitudes.
A thermodynamic system is said to be in a certain state when all its properties are fixed. The
fundamental properties which determine state of a system are Pressure (P), temperature (T),
volume (V), mass and composition. Change in magnitude of such properties alters the state of
system. For such reasons, these parameters are regarded as state variables or state functions
for a system. A change in system from the initial state to final state will accompanied by a
change in state variables. However, it is not necessary that all the state variables under
specified set of conditions may completely define a system. The relation among state
variables is defined as equation of state. For an ideal gaseous system, the equation of state PV
=nRT, where R=gas constant may be isolated through implication of the fundamental
approach based on Boyles and Charles laws.
Binary relation of state variables under specified set of conditions liberates thermodynamic
properties of a system. Such thermodynamic properties are internal energy (E) enthalpy (H),
entropy (S), Gibbs free energy (G), Helmholtz function(A), chemical potential (µ ), heat
capacity etc. For example, at fixed mass of an ideal gas, at constant volume, H
=f(T,P).Similarly other thermodynamic properties such as internal energy and entropy may
be defined as E=f(T,V) and S=f(E,V) respectively.
Counter balanced, the piston remains stationary and there is no change in the volume of the
gas. Now suppose the pressure on the piston is decreased by an infinitesimal amount dP, thus
the external pressure on the piston being P-dP, the piston moves up and the gas will expend
by an infinitesimal small amount. The gas will therefore be expended infinitely slowly
(reversible process). At all stages in the expansion of the gas, dP being negligibly small the
gas is maintained in a state of equilibrium throughout. If at any point of the process the
pressure is increased by dP, the gas would contract reversibly fig.1.1 (a).
Now on the after hand, the expansion is irreversible, if the pressure on the piston falls
suddenly, it moves upward rapidly in a single operation fig.1.1 (b).the gas is in equilibrium
state in the initial and final stages only, the expansion of the gas, in this case, takes place in
an irreversible manner.
If the gas expands at constant pressure, the piston will move to a distance l. Since,
Work = force distance
=ƒ l ……………….1.1
and pressure is force per unit area
=Pext
or f =Pext A ………………..1.2
Work is done at the expense of internal energy, hence the internal energy decreases and work
done bears a negative sign. Work may be expressed in dynes-centimeter, ergs, joules. PV
work as well be expressed as the product of pressure and volume units e.g., litre /
atmospheres.
If the initial volume is V1and final volume is V2 then work in this case will be
It is to be noticed that the work done by the system is not a state function, it depends
upon the path of transformation. In the case just cited above pressure the piston is decreased
as volume changes from initial volumeV1 to final volume V2, this is irreversible work.
Consider a system in which an ideal gas is confined in a cylinder with a frictionless piston
fig.1.3.
Suppose it expands reversibly from V1to V2 at a constant temperature. The pressure of the gas
is successively reduced from P1to P2. Let us reduce the pressure P1 by an infinitesimal
amount p, i.e., let the pressure be P1- p and the volume increases be , then the work
done by the system would be ( P1 - p) P1
As p is very small, will also be very small, and the multiplication of these two will be
very small quantity and can be neglected. Then the new pressure is again reduced
infinitesimally and a little more expansion involving again P work would be available. In
this way by infinitely large number of successive steps we can expend the gas from
volumeV1 to V2. Let the equilibrium pressure attained by the gas, when the volume becomes
V2 be P2. The total work done by the gas during expansion would be summation of all these
minute quantities, i.e.,
W= dV
………………………1.4
Here dP. will be very small quantity and can be neglected. We must remember that P has
also been changing during the expansion continuously. If the system contains n gm.moles of
the gas and the gas be supposed ideal one, then
W= =nRT ln
W =RT ln ………………………..1.5
as =
W = RT ln ………………………..1.6
Since the change has been carried out in successive stages of infinitesimal amount and
therefore it is reversible process. The work in the isothermal reversible expansion of one
mole is
W =RT ln = RT ln ………………………1.7
This work is the product of force and displacement. The displacement in this case is
specified, namely (V2 – V1). Hence greater the magnitude of the force or pressure, the greater
would be magnitude of the work. If the pressure over the piston were made greater than (P1 -
) there obviously would be no expansion.
Therefore if we want to have expansion as desired, (P1 - ) is the maximum that could be
applied. So the work obtained in the equation (1.7) is the maximum. That could be applied. If
the pressure over the piston were made appreciably less, instead of infinitesimally smaller
than the gas pressure, the work would be evidently less than obtained above and the process
also would be more irreversible. Hence we conclude that the work in the reversible process is
the maximum and is greater than that in the irreversible process.
For this we carry out the expansion of the same gas in two ways. Firstly, we change
the volume as before until the pressure is P2, the process is reversible work.
Wr =RT ln = RT ln ………………………..1.7
Secondly, we reduce the pressure over the piston at once to P2, which would be equilibrium
pressure of the gas when the volume is V2, then the gas would immediately expands against a
constant pressure P2 the process is irreversible work is
Wir = P2 (V2-V1) = P2 ( -
= RT (1 - )
Wr –wir = RT ln - RT (1 - )
= RT ( ) – RT ( )
= RT ( – )( - )
= RT ( )( )
= RT ………………..1.8
We know that when gas has expanded to volume V2 the gas pressure has come down to P2. If
the pressure is increased by infinitesimal amounts in successive stages the volume becomes V1
and therefore, pressureP1. The work done in the system will be
W` = 2. V=
= RT ln = - RT ln ………………….1.9
Here the work term becomes the same with reversal of sign. Thus if the complete reversible
cycle of expansion and contraction is carried out at constant temperature then the net work
would be zero.
= W + W` = RT ln - RT ln =0
All class of physical, chemical or biological transformations into one another are assumed to
facilitate at the cost of conservation of energy and matter. The transformations and conservation
of energy during processes are studied under the first law of thermodynamics. According to this
law, the energy transformations in all the processes are purely conservative and proceed without
any creation or loss of energy. Such kind of statements however doesn’t hold when Einstein’s
relation inter-conversions of energy and matter are considered as below:
E= mc2 ………………..1.11
Where E is energy associated with matter of mass ‘m’ and ‘c’ is the velocity of light in vacuum,
a constant. In our present discussion, we shall have no occasion to consider this inter conversion
of energy and matter and we shall take the first law to be strictly true.
dE = ( )V dT +( T dV
The internal energy of a system is determined by the state of a system and is dependent of the
path by which it is obtained. Hence the internal energy of a system is a state function.
EA+q =EB+w
or E = q-w …………………1.12
Here E is state function hence its change is dE and q and w are not state functions, they depend
upon the path, hence, q and w.
1.8 ENTHALPY
Suppose the change of state of a system is brought at constant pressure P. In that case there will
be change of volume from VA to VB. Under such conditions, system may execute the work
done (w) as below:
E = q – P (VB –VA)
EB-EA = q –P (VB – VA)
The quantity E + PV is known as total energy stored and thermodynamically termed as enthalpy
(H) of the system
H = E + PV ..………………1.15
In this equation E, P and V are all state functions, thus H, the value of which depends upon the
value of E, P and V must also be the function of state. Hence the value is independent of the
path by which the state of the system is changed.
From equation 1.14 and 1.15 are have
HB-HA = H …………………1.16
Since HA and HA are definite quantities, it is evident that H, like E is a definite
quantity depending upon initial and final states. We can write
H= +P V …………………1.17
Here V is increase in the volume under gone by the system.
For a fixed mass of system H=f(T,P), that represents the following exact differential
Heat capacity is defined as the quantity of heat required to raise the temperature of the
system from the lower to the higher temperature divided by the temperature difference. If q
calories are the heat absorbed by mass m and the temperature rises from T1 to T2, the heat
capacity is given by the expression
c= …………………..1.18
For unit mass say restricted to 1 gm, heat capacity is called the specific heat of the system.
However, if mass of the system is one mole, then the heat capacity is called molar heat capacity
(C). Heat capacity varies with temperature; hence the true molar heat capacity is defined by the
equation
C= …………………..1.19
Molar heat capacity is defined as the amount of heat required to raise the temperature of
one mole of substance or system by 1K.
Since heat is not the state function, it is necessary to fix either of a thermodynamic
variable such as P or V at which the temperature of system is raised by one degree. This
generates two different kinds of the heat capacities. First kind of heat capacity is defined as
“heat capacity at constant volume (Cv)” and another one is the “heat capacity at constant
pressure (Cp)”.
From first law of thermodynamics we know
dq = dE + PdV
Hence = + …………………….1.20
Cv =( )V ……………………..1.21
Cv = +P )P .……………………1.22
and H = E+ PV
Hence ( )P = ( P +P )P .……………………1.23
CP = ( )P
Now Cp –Cv = ( )P - ( )V
Since H = E+ PV
Cp –Cv =( P +P )P - ( )V ………………..1.24
We know E = f (V, T )
dE = ( )V dT +( T dV
( P =( )V +( T( )p .……………….1.25
Cp –Cv = ( )V +( T( )p + P )p - ( )V
=( )p {( +P}
Hence Cp –Cv = P )P = P.
Cp –Cv = R …………………1.26
but dE = T dV + v dT
Since dT = 0
dE = T dV = 0
as dV since the gas which was occupying one container, now occupies both the containers
Hence T =0 …………………..1.27
We can write T = T T
Since T so T =0 ..…………………1.28
We know H = E + PV
T = T +
T = 0 ………………….1.29
Experimental Details:
The experiment performed by Joule and Thomson to measure the temperature change on
expansion of a real gas is illustrated in fig 1.5
Since the expansion of the gas is done adiabatically, the system is not in a position to absorb
heat from the surroundings. The system thus does work at the expense of the internal energy.
Hence the internal energy decreases from E1 to E2
P2V2 – P1V1 = E1-E2
or P2V2 +E2 = P1V1+ E1
or H2 = H1 .…………….1.30
or H =0
This reveals that Joule-Thomson expansion of a gas takes place at constant enthalpy. Hence
this kind of expansion is also called isoenthalpic.
The rate of change in temperature with pressure when enthalpy remains constant is
called Joule-Thomson coefficient.
……….……..1.31
Since such flow through porous plug, usually the gas runs into lower pressure, i.e., dp
is negative. Hence will be positive if dT is also negative, i.e., if temperature falls. Except
Hydrogen and Helium at ordinary temperature almost all gases cool down in passing through
porous plug. For Hydrogen and Helium .
We know
or
……………..1.32
……………1.33
In this equation of a gas is usually positive quantity, while for any gas is
always a negative quantity. As a result the first term is always positive. The magnitude of
second term may be negative or positive. At low temperature the value of this
term is negative under these conditions both the term are positive, has positive value. That
means cooling will take place.
For ideal gases PV is constant and first term is zero hence the =0.
Joule-Thomson coefficient for a real gas can be calculated with the help of van der
Waals equation. Since, both a and b are small, the term in the equation can be
neglected provided pressure is not too high. The Von der Waal equation can be written as
………………1.34
Replacing V by we have
………………1.35
Substituting the value of R/P from equation 1.36 in equation 1.35 we get
…….…………1.40
From this equation it is clear that Joule – Thomson coefficient is positive as long as is
greater than b. This is the case when temperature is low. At high temperature hence
s negative. Since a, b and R are constant, it is evident that the magnitude and sign of Joule
Thomson coefficient depends only upon the temperature at which the gas is allowed to
expand.
The temperature at which the Joule – Thomson coefficient changes sign is known as
the inversion temperature ) At this temperature J.T. is zero so that
or = ………………….1.41
Here represents inversion temperature. This inversion temperature depends upon the van
der Waals constant a and b of the gas.
Where Q1 is heat absorbed at temperature T1 and Q2 the heat rejected at temperature T2.
The amount of heat absorbed or evolved in a process carried out reversibly at temperature T
may be evaluated through application of Carnot cycle.
Suppose a system changes from one to another state, the heat dq is either absorbed or
evolved. The value of dq depends upon path of transformation and hence it is not a perfect
differential. From Carnot cycle we know the efficiency of an engine is given by the relation
=
or .………………1.43
Where qrev is the amount of heat absorbed or evolved in a process carried out reversibly at
temperature T. The term qrev/T denotes a definite quantity or a state function, called entropy
change of the system.
If we assign positive sign to the heat absorbed (q2) and negative sign to heat given out q1 the
equation 1.43 can be written as
or ....…………..1.44
Fig. 1.6 demonstrates a reversible cycle ABA to be made up of large number of Carnot cycles
and adiabatic curves under isothermal conditions.
The cycle ABA can be regarded as comprising of two paths, i.e., from A to B and back from
B to A, thus
…….1.46
The integral is the summation of all the terms when the system changes from A to B
along the path I and similarly integral is the summation of all the terms when the
From equation 1.47 it is clear that is a definite quantity independent of path taken for
the change and depends only on initial state and final state. In other words it is a perfect
differential. The function is called entropy and is denoted by symbol S first proposed
………………..1.48
It is not easy to define the actual entropy of a system. However it is easier to define change of
entropy during a change. It is clear from equation 1.48 that when heat is absorbed by the
system dq is positive. T is always positive. Thus is positive. That means entropy of the
system increases by the absorption of heat. Hence, evolution of heat by the system means
entropy decreases.
Since the entropy change is expressed as heat absorbed or evolved per unit temperature,
entropy data are numerically expressed in terms of joules per degree per mole (JK-1). This is
known as entropy unit, e.u.
In isothermal reversible expansion of an ideal gas, there will be no change in internal energy
dE=0. In that case from first law
or
If the gas is expanding from volume V1 to V2 for n moles
………………1.49
Entropy of one mole of a substance in pure state at 250C and one atmospheric pressure is
defined as standard entropy of that substance and is denoted by S0: During the reaction, each
reactant and each product are in their standard state.Under such conditions; the overall
entropy change may be evaluated as:
When a system is left to itself it will gradually pass into a state of equilibrium. This is
common experience. The state of equilibrium means a state of maximum disorder or chaos. If
we let loose a quantity of yellowish green chlorine gas at corner of the floor of the room, the
gas would spread in all directions until it would be most chaotically distributed into the entire
room and the equilibrium will be reached. If we arrange a number of red and black balls in a
tray in a particular pattern and then shake the tray, the balls will be mixed up in the most
disordered fashion. These are all natural processes and we never find the opposite to happen
in nature. These processes are irreversible and thus lead to entropy increase. The more the
system reaches equilibrium, the greater is the disorder and irreversibility, and when the
equilibrium is reached, it has the maximum disorder. Hence the entropy also goes on
increasing and reaches its maximum at the equilibrium state. In other words, entropy is the
measure of the disorder of the system.
Now if we attempt to bring order to system entropy decreases. Suppose we have a gas in
equilibrium confined in a cylinder with a piston. The piston is now some what raised upwards
increasing the volume slightly. All the molecules, in whatever direction those might be
moving, are subjected now to a tendency of motion in the direction of piston. This is
indirectly an attempt to bring some order. In consequences, the gas is cooled and the entropy
decreases. So in the crystalline state of pure substances the entropy becomes the least as there
is perfect order or least chaos.
We know,
dq=dE+PdV
or ………….1.50
on integrating ….……….1.51
is constant of integration.
rearranging (putting we get
…………..1.52
We know
or ……………..1.53
Since
or ……………..1.54
When two or more gases are brought in contact, they immediately diffuse into one
another and mix up irreversible. Such a spontaneous process would lead to increase in
entropy. If the gases are ideal we can drive the relation for mixing.
Let one mole of gas A and one mole of gas B are mixed together at constant
temperature T and pressure P. Before mixing the entropies are
After mixing, the final pressure remaining constant, the partial pressures of A and B may be
denoted by PA and PB such that P=PA+PB the entropies of two gases in the mixture are
….………..1.55
Till now we have studied only three thermodynamic functions, viz. the internal energy (E),
enthalpy (H) and entropy (S). In addition to these three thermodynamic functions there are
two additional thermodynamic functions which also depend only on initial and final states of
system. These functions are work function and free energy functions represented by A and G
respectively. These functions are such that they utilize entropy in their derivation. With these
functions the entropy change can be determined more conveniently.
……………………1.56
…………………...1.57
We know E, H and S depend on the state of the system, it is clear that both A and G
also depend upon the state only. The exact nature of these functions will be clear from the
following discussion.
Let us consider that a thermodynamic system changes from initial state to final state at
constant temperature. Let A1, E1 and S1 be the work function, internal energy and entropy in
the initial state and A2, E2 and S2 are the values of these functions in the final state
respectively then,
A2=E2-TS2
A1 = E1- TS1
Subtracting second equation from first equation
We get +
or ………….………1.58
or
…....…………….1.59
i.e. decrease in work function A in any process at constant temperature gives that maximum
work that can be obtained from the system during the change this function, A is therefore
termed as work function this is also referred as Helmholtz free energy or Helmholtz function.
Similarly, from free energy function we get
or
In this equation PV signifies the work done against pressure. Hence gives the maximum
work obtainable from the system other than pressure volume work. This is known as net
work.
If a closed system undergoes a change, the only work involved being mechanical against the
pressure, the infinitesimal reversible changes of A and G. would be
………………1.60
and
………………...1.61
From equation 1.62 and 1.63 different forms of Gibbs Helmholtz equation can be derived by
substituting these values of in their respective equations of work function and Gibbs free
energy.
...………….1.64
……………1.65
…………….1.66
………….…1.67
and
and
So
Similarly
………….1.68
and
These are four Maxwells relations. The last two relations give us the isothermal dependence
of entropy on volume and pressure.
CONTENTS
2.1 Introduction
2.2 Partial Molar Quantities/Chemical Potential
2.3 Gibbs Duhem Equations
2.4 Variation of Chemical Potential
2.4.1 Variation with Temperature
2.4.2 Variation of Chemical Potential with Pressure
2.5 Clapeyron Equation
2.6 Clasius Clapeyron Equation
2.7 Fugacity
2.8 Physical Significance of Fugacity
2.9 Thermodynamic Derivation of Law of Mass Action
2.10 Free Energy Change (G) in Chemical Reactions
2.11 Activity Concept
2.12 Activity Coefficients of Electrolytes
2.13 Changes in Gibbs free energy (G) in Chemical Reactions
2.14 Standard Value of Gibbs Potential
2.15 van’t Hoffs Equation and van’t Hoffs Isochore
2.16 Feasibility of a Reaction
2.1 INTRODUCTION
This unit deals with thermodynamics of open system, phase equilibria, fugacity and
derivation of law of mass action.
The term is called partial molar property for the first component. This is more
often represented as 1 which denotes the rate of change of X on addition of one mole of the
first component keeping other variables unchanged.
For change in the mass of system, at constant temperature and pressure, we have
= 1dn1+ 2dn2+ 3dn3 ……………2.2
It may be noted that X stands for any extensive property.
Thus for ith component in a system,
Partial molar internal energy = = i
The most important partial molar property is partial molar free energy.
The change in the thermodynamic potential (Gibbs potential G) at a constant
temperature and pressure can be expressed as
……. ………………2.4
J.W. Gibbs proposed that for component of an open system, chemical potential is
represented as .
Hence = i= …..……………2.5
…..……………2.6
It must be noted that is independent of the mass of system but not independent of
composition of the system.
Suppose at constant temperature and pressure, all the molar concentration of components are
increased in the same proportion of their initial amounts. This implies to increase of first
Thus we have three useful relations for chemical potential with the composition of the
system.
(i) + + …………………=0
(ii) ………….= …..…. 2.10
(iii) …………..= …………….2.10
+ =0 ……………..2.11
Here n1 and n2 are number of moles of the two components and and are respective
chemical potentials. Dividing the above equation by n1+n2 we get
or =0 ……..………2.12
Where is mole fraction of constituent 1 and mole fraction of constituent 2. We can
write equation 2.12 in the following form
+ =0 ..…………….2.13
We know
or -
Substituting the value of in equation 2.13 we get
- =0
or = =0 ………….…2.14
= …..………….2.16
We know =S …..…………2.17
Hence differentiating this equation with respect to at constant temperature and pressure we
get
= ..…………….2.18
= = i …….………..2.19
This is partial molar entropy by definition. This equation gives variation of chemical
potential of any constituent i with temperature.
We know =
= ………………2.20
We know =V
Differentiating above equation with respect to at constant temperature and pressure we get
=(
Now, suppose there is a system containing a mixture of ideal gases at constant temperature T,
we get
n= n1+n2+n3 + …………………….. + ……..
and
or =
But we know partial pressure of the component,
For a system ni and n are fixed that is constant, hence
so
or ( being constant of integration)
The constant Z depends on temperature and nature of gas. At standard state chemical
potential is denoted by
Then RTln+ or = …..………..2.22
This is the relation connecting chemical potential with pressure of a component in an ideal
gas mixture.
We know where molar concentration of component is .
Hence
Thus we have
……………2.23
Since all gas laws of ideal gases also extended to ideal solution, the equation 2.23 will hold
true for ideal solutions also.
Consider two phases A and B in equilibrium of the same component with each other
at constant temperature (T) and pressure (P). If the two phases are in equilibrium then free
energy in both the phases should be same, i.e.,
Where is free energy in phase A and free energy in phase B. now temperature
is raised to and pressure is raised to P+dP then free energy in phase A will be
(suppose) and in phase B will be . Since the
and
=
Rearranging the term we get
is change in entropy
We know from our earlier knowledge
At equilibrium G=0, so H= S
or =
So we get H/ ……………2.26
are molar volumes of the component in B and A phases respectively and represents rate of
In equation 2.27, since volume in vapour phase is much greater than volume in liquid
phase; hence. Substituting with , the equation 2.27 is
reduced as
Assuming that the system behaves ideally, the above equation may modified through
substituting
or
or …………….. 2.28
The above equation is defined as Clausius Clapeyron equation. This equation has following
applications
2.7 FUGACITY
The actual behavior of real gases is different from behavior of idea gases. Lewis
introduced concept of fugacity for representing actual behavior of real gases. The fugacity, in
substituting for pressure, will include in it the imperfections of the real gases and yet
retaining the form of equation. Thus, fugacity is a sort of fictitious pressure “which is used to
order to retain for real gases simple forms of equations which is used for ideal gases only.”
Integrating the equation 2.30 for fugacities and at constant temperature T gives
……………. 2.31
If P is actual pressure, the ratio approaches to unity, when P approaches to zero, in that case
P
At low pressure fugacity is equal to pressure. But at high pressure the two terms differ
materially.
Now
Differentiating this equation will respect to pressure at constant temperature and constant
number of moles of the various constituents. i.e., in a closed system, we get.
…………… 2.33
………. 2.35
Now as given in the above equation can be determined at different pressures. These values
are plotted against corresponding pressure as shown in fig. 2.1. The shaded area under the
curve between pressure P=0 and any given pressure P, gives the value of the integral
f (fugacity) can be evaluated at any given pressure P of the gas through substituting
this value in equation 2.35,. Since departure from ideal behavior ( ) can be both positive as
well as negative. Thus, the fugacity of the gas can be less than the pressure or more than the
pressure. It is clear from above figure the area and hence the value of is positive at low
pressure and more than P at high pressures.
van’t Hoff deduced the law of mass action thermodynamically by using an equilibrium box in
which all the systems were in equilibrium with each other at constant temperature.
(a) The walls of the box are rigid having indefinite dimensions.
(c) The various components present are in equilibrium with each other.
(e) The box is assumed to have large dimensions so that there is no change in
concentration or pressure (partial) of any component, if a small amount of those is
added or removed from it.
Consider two equilibrium boxes 1 and 2 in which components A,B,C and D are in
equilibrium, Suppose PA,PB,PC and PD are partial pressures of the A, B, C and D and CA, CB,
CC, and CD are the concentrations for box no. 1 similarly P`A, P` B, P`C and P`D are partial
pressures of A, B, C and D and C`A, C`B, C`C, and C`D are corresponding concentrations in the
box no. 2.
(a) Let us calculate the work done in transferring one mole of A. By means of the piston
moving outwards infinitesimally slowly, one mole of A is removed from box (1) so that
the partial pressure of A (which is Pa) remains constant. The amount of work done is
PAVA. Where VA is the volume of one mole of A under pressure PA.
(b) As the partial pressure is different in two boxes, that the removed gas from box(1) is
changed isothermally and reversibly from PA to P`A by working the piston.
(c) The gas is now ready to be put in the box no. (2) by compressing the gas infinitesimally
slowly by working the piston, so the partial pressure of A in the box (2) (which is P`A)
remains constant. This will be put in box no. 2.
The work done in the process is P`AV`A where V`A is volume of one mole of A under a partial
pressure of P`A
….…………. 2.37
Similarly one mole of B is transferred from box-1 to box-2, the work done will be
……………….2.38
At the same time one mole of C and one mole of D are transferred from box(2) to box(1),
they interact to produce A or B to maintain the equilibrium the work done will be
…………… 2.39
…………… 2.40
When one mole of C and one mole of D are transferred from box(2) to box(1) they interact to
produce A and B to maintain the equilibrium.
For such an isothermal cyclic process total amount of work done is zero, as the system
remains always in equilibrium.
Here
or ln
Similarly ……………2.42
Equation 2.41 and 2.42 are the two forms of law of mass action.
Suppose at a given temperature T and at a given set of arbitrarily chosen pressure (or
concentration) one mole of A reacts with one mole of B yielding one mole of C and one mole
of D, i.e.,
Since the absolute value of chemical potential is not known, it is not possible to
evaluate the value of of a substance. The function fugacity was introduced into the
expression of the thermodynamic potential to cover the departure from ideal behavior of the
system.Lewis introduced another function called activity or relative fugacity, denoted by ‘a’.
If and denote the fugacities of a substance in a given state and in a standard state,
arbitrarily chosen, then the activity ‘a’ is the ratio of two fugacities,
……………… 2.43
From the above it is clear that activity of a substance in any given state is defined as the ratio
of the fugacity of the substance in that state to the fugacity of same substance in standard
state.
The difference between chemical potential of a substance in any state and that in the
pure state is given by
So ..……………2.44
Since activity is the ratio of fugacity, its numerical value will depend upon the standard state,
which may be arbitrary fixed. Since, activity and the fugacity are numerically the same. For
gases and vapours, at standard state, fugacity is unity. Thus a=f/f0=f/1=f
For an ideal gas fugacity and pressure are equal and the same is also true for other gases at
low pressures. Hence in all such cases
For pure liquids and solids activities in pure state under one atmospheric pressure are taken as
unity at each temperature. But in solutions it is not always convenient to assume the pure
state as the standard for unit activity. It is necessary to consider the activities of solutes and
solvents separately.
Let us consider a reaction to take place in solution
………………….2.45
Where aA aB, aC and aD represent the ideal concentrations in terms of molality (m) or activities of
the respective constituents A, B, C and D. If ‘a’ represents the ideal concentration and ‘c’ or
‘m’ the real concentration, the ratio a/m is known as activity coefficient ( ). Activity
coefficient is the measure of extent to which an ion or molecule departs from ideal behavior.
The ionization of a uni-univalent electrolyte of type AB can be represented as
where a+ is the activity of cation (A+), a-, the activity of anion (B-) and a2 is the
activity of the undissociated electrolyte (AB).
The mean ionic activity of the two ions (a ) may be regarded as the geometric mean of the
activities of the respective ions
a- =
The activity coefficient of cation is and that of anion is . The mean activity
coefficient of the two ions is regarded as the geometric mean of the activity
coefficients of the respective ions.
Where
The concept of activity coefficient generally applicable to all cases. Because of the
ionic dissociation in solution of electrolyte the expression require certain modification. At a
given temperature under unit pressure the standard state of an ionic species is such that at
infinite dilution, the activity becomes equal to its concentration. The concentration is
generally expressed in molality or molarity. Then
Lt m 0 Lt c 0
+ Sign refer to cations and – signs refer to anions. Now let us consider a simple case
of uni-univalent electrolyte AB, which ionizes in solutions as:
Suppose and are the chemical potentials of the electrolyte molecule, cation and
anion respectively. Then
Where terms are constant, representing corresponding chemical potentials at the standard
state.
i.e.,
or
Therefore
We can write
or
Let us consider a general reaction in which a mole of A are reacting with b moles of B
to produce c moles of C and d moles of D at a given temperature T,
The change in thermodynamics potential (Gibbs potential) in the system due to the chemical
reaction at constant temperature and pressure, is given by
= (d )
We know
We have
G = d(RT + +c(RT
-a(RTln
= RTln(
where Z is a constant
In this expression and are the activities of the components in the arbiturary
choosen conditions for the reaction.
It must be kept in mind that this relation is valid for any chosen state and therefore
would also be valid it we start with components in the equilibrium state. Now, using o for
denoting equilibrium conditions,
= RTln +
But at equilibrium =0
Then …………….2.47
Then for any equilibrium reaction, the change in the thermodynamic potential would be
…………….2.48
As this expression is derived at constant pressure this relation is commonly known as van’t
Hoff’s reaction isotherm.
In terms of fugacity
In terms of concentration
=- ....………….2.49
In terms of pressure
= .……………2.50
The importance of this reaction isotherm lies in the fact that can predict the feasibility of the
reaction. It is well known fact that when a reaction precedes spontaneously, the Gibbs
potential gradually diminishes until it becomes minimum at equilibrium. This means that for
spontaneous reactions is negative. That is in equation 2.49. should be less than
ln If the conditions are so chosen that becomes greater than ln then the reaction
will not occur. If we reduce the concentrations (or pressures) of the products and increase
those of reactants, then would diminish and the possibility of the reaction in the
desired direction improves. Thus we can determine the direction of a chemical change.
The value of change in Gibbs potential is given by the relation, as we have deduced above.
If the reaction is carried out in such a condition that all the reactants and products are at unit
activity, under these conditions the change in thermodynamics potential is called the standard
change in Gibbs potential (G) and is denoted at , so,
Because then
…………….2.51
and
Then
and
So
We know,
We get
…...…….2.53
In this expression is the energy change at constant volume, i.e., heat change at
constant volume. This relation for the variation of equilibrium constant with temperature
at constant volume is commonly mentioned as van’t Hoff isochore.
Free energy (G) and work function (A) are useful in determining feasibility of a
reaction or its attainment of equilibrium. We know that entropy remains unchanged in
reversible change, but irreversible changes are accompanied with increase in entropy,
therefore
Now suppose there is a system with its surroundings. Let ds and ds` be the entropy-
change in the system and surroundings respectively. Therefore
………….2.54
If we regard the change in the surrounding as reversible and suppose the surroundings supply
heat dq reversibly, then
If dE is internal energy increase of the system and w be the work done by the system, then
from the first law,
Now
Then
Hence,
……………2.57
The inequality sign for real changes and the equality sign for condition of equilibrium.
This means at equilibrium S becomes maximum and
…….………2.58
We get,
or
……………..2.59
For real processes i.e. free energy diminishes. When equilibrium is attained
free energy becomes minimum and
We get,
or ……………..2.60
For real process Gibbs free energy decreases . When equilibrium attained, Gibbs free
energy is minimum,
For a process to occur spontaneously , in other words Gibbs free energy change
should be negative. Every process tends to attain equilibrium, so G will diminish
gradually and at equilibrium it will attain the minimum value, i.e.,
and at equilibrium
Since
, ..………..2.61
The Gibbs free energy change can be determined if we can find the change in enthalpy
and the change in entropy ( for the reaction.
The enthalpy change is obtained from the measurement of the heat of reaction
or it may be calculated from the bond energies in the molecules of reactants and products.
The entropy change is obtained from the individual entropies of reactants and
products, or from other measurements such as e.m.f. etc. The entropy increases whenever as a
result of the reaction the randomness increases. It is generally true that when a solid reactant
gives gaseous products, the entropy increases, i.e., is positive. In gaseous reactions the
sign of the entropy-change is determined from the number of product and reactant molecules.
If the number of product molecules be larger, will be positive. All dissociation processes
are attended with an increase in entropy.
From equation 2.60 it is clear that (i) the more negative is more negative and (ii) more
positive the more negative is . This shows that the system at constant pressure and
temperature would try to move to the state of maximum entropy and minimum enthalpy.
Example -1
Example -2
This reaction is endothermic. This reaction will proceed if is positive and greater
than . Here entropy increases, as reactant C has least disorder, cal/degree. The
enthalpy change was found to be +31400 cal.
Hence at 250C
Therefore the reaction will not occur at 2980K but will be spontaneous at 12000K.
TERMINAL QUESTION
5. At NTP, 2.8 litre of oxygen were mixed with 19.6 litre of hydrogen. calculate the
entropy change.
S = CV ln P + Cp lnV + constant
1. (a) Starting from the limitation of Gibbs-Helmholtz relation, discuss Nernst heat
theorem. What are the limitations of Nernst heat theorem? Discuss third law of
thermodynamics.
(b) Derive the relation for internal energy, entropy, pressure and heat capacity in
CONTENTS
3.1 Introduction
3.2 Probability
3.3 Boltzmann Distribution
3.4 Partition Function
3.5 Partition Function and Thermodynamic Functions
3.5.1 Internal energy
3.5.2 Entropy
3.6 The Third Law of Thermodynamics
3.7 Bose – Einstein Statistics
3.8 Fermi – Dirac Statics
3.9 Molar Partition Function
3.10 Terminal question
3.1 INTRODUCTION
particle. Statistical methods are essentially employed to express the property of the system,
particularly to find out what fractions of the molecules have a given energy.
Statistical mechanics can be applied easily to simple ideal systems such as monoatomic gases
and diatomic gases. For applications to interacting systems such as liquids (where strong
intermolecular forces exist), the details of the intermolecular potential energy, which is not
accurately known, have also to be taken into account. That is why statistical mechanics of
liquids is a difficult but fascinating subject. Gases under high pressure too are difficult to treat
statistically since they deviate strongly from ideality. In recent years statistical methods have
been applied successfully to simple liquids and dense gases. It is easier to find out the value
of entropy, free energy, enthalpy, heat capacity etc. from spectroscopic data by using
statistical method.
3.2 PROBABILITY
If we toss a coin, either head or tail will come up. We cannot predict whether head or
tail will come up in a particular toss. But we can reasonably say that in any toss there is a
fifty percent chance of the head being up. That is desired probability of the event is one in
two. Let us take another case. If a single card is drawn out from a well-shuffled deck of cards,
the probability or chance of drawing a heart is one in four, the probability of drawing a king
is one in thirteen, and the probability of drawing a king of hearts is one in fifty two.
If there are n equally likely and mutually exclusive ways of occurrence of an event in
a given set of conditions, then out of these n ways m are favorable to an event A, then the
mathematical probability of A is m/n. We may define this in another way. If a occurrence
succeeds in p ways and fails in q ways, all the ways being exclusive and equally likely, then
the probability of the event succeeding is p/(p+q). For example if we want to hit in color out
of the four colors, then we can hit at any particular color in one way, i.e., p=1, we may hit at
any of the three color, i.e., equally likely cases of failure is 3, i.e., q=3. Hence the probability
of succeeding is given by
1st die 6 5 4 3 2
nd
2 die 2 3 4 5 6
Now let us consider the first combination (6 and 2), there is only one way out of six in
which six will come at the top of the first die, and also there is only one way for a two to
come at the top of second die. Thus the probability of obtaining eight by first combination is
thus
This is true for every other four remaining combinations. Hence the total probability of
getting an eight is
(ii) The probability of the simultaneous occurrence of two events A and B, is given by the
product of the probability of event A and the (conditional) probability of B, when the event A
has occurred.
If wA is the unconditional probability of A and wB is probability of B when A has
already occurred, then the probability of their simultaneous occurrence (wAB) is given by
WAB = wAwB
Suppose from a well shuffled deck of cards we intend to draw the king and queen of hearts.
The probability of obtaining king of hearts by a single draw is 1/52. If the king of hearts has
been received, the probability of securing the queen of hearts from the remaining is 1/51. The
total probability of simultaneously obtaining the pair is
Now let us discuss some other cases. Let us toss two different coins a large number
of times and count the number in which
Where t denotes tail and h denotes head and suffix 1 and 2 indicate the first and second coin.
Now there are two possible events for first coin and two possible events for second
coin. Hence total possible events will be (t1+h1)(t2+h2) =2x2=4. Obviously any one of the
four cases can occur. The probability of each of the four equally like cases is ¼.
If we have four coins, the total possible cases by simultaneous toss of the four coins is
(t+h)4 = 24 =16. There are, however, five possible combinations or states. These are
tttt, ttth, tthh,thhh, hhhh
Now we select any one of these combinations say tthh (two tails and two heads
upwards)
That is there are six ways in which two tails and two heads be upwards from 4 coins. For th3
there are four possible ways to obtaining the same.
The total possible event is given by (t+h)4.
Expending this binomially we get
(t+h)4 = t4 + 4t3h + 6t2h2 + 4th3 + h4. In this equation the exponents of t and h indicate the
number of tails and heads upwards in a given combination. If you examine the equation the
coefficient of the term indicate the number of ways the particular combination would occur,
the total is 24 or 16. Thus probability of t2h2 is or 6x2-4 This is mathematical probability
of t2h2.
The same result can be obtained as follows. Total number of possible events is (t1+h1)
(t2+h2) i.e. 4. The total permutation of these four is 4! But this number also includes the
permutation of two events of first coin and two events of second coin for t2h2 which is 2! x 2!.
Where denotes product of terms following it. Again we have taken the case of coins.
Suppose we have n identical coins tossed as described earlier. Following our usual method to
all the n coins. The total number of cases or macrostate that can occur is given by (t+h)n=2n
The possible combinations or microstates are
tn, tn-1h, tn-2h2 …………tptq hn-1, hn
when p+q=n
The coefficient of the terms tphq in the binomial expansion of the above is
This is the total number of complexions or microstates that would arise when p coins have
tails and q coins have heads upwards. This can be done nCp ways.
The mathematical probability of the occurrence of tphq is therefore
W= nCp/2n = ……………3.1
Thermodynamic probability
Mathematical probability
Where each individual item of the system can give rise to c-cases.
It is clear from these two probabilities that because total number of possible cases (Cn)
is constant, thermodynamic probability (W) is proportional to mathematical probability (P)
……………..3.2
E= …………….. 3.3
The allowed energy levels are fixed by quantum mechanical principles. As the
occupation number n1, n2 etc in different energy levels change the distribution also changes
there are obviously number of ways of distribution, similar to the arrangement of n boys in p
rooms, so that n1, in one room, n2 boys in other room and so on. This number is given by
………………..3.4
the level, there are possible distribution. Thus for particles in ith level, there are gini
possible distribution. Hence the thermodynamic probability for the system of particles is
given by
…………….3.5
We are thus interested in finding the distribution that will make W a maximum.
It is more convenient to maximize the logarithm of P. It is known from calculus that
at the maximum the derivative of function vanishes.
Now we confine our investigation to a closed system of independent particles, it
would meet the following two requirements:
(i) Total number of particles is constant i.e.
………..…….3.7
(ii) Total energy, E of the system is constant i.e.,
………………3.8
Since total number of particles is constant
…………..…..3.9
and total energy is also constant
….…….…..3.10
Taking logarithms of both the sides of equation 3.5 we have
…………..3.11
..…………3.12
…..…………3.13
As …..………..3.14
At equilibrium
……………..3.15
If the system were open then would vary without restriction and the variation
would be independent of one another. It would then be possible to solve the equation 2.15 by
setting each of the coefficients of the terms equal to zero. But our system is closed one
end since is constant, the value of are not independent of one another, as it is seen in
equation no 3.14.
To solve the equation 3.15 we apply the Lagrange’s method of undetermined
multipliers. We can write equation 3.15 in the following form
……………….3.16
Now multiply 3.9 by and 3.10 by and subtract the addition from equation 3.16
…..…………..3.17
The variable etc is independent of each other, so that the equation 3.17 to
hold each term in the summation must be zero.
…..………….3.18
Now
or …………….. 3.19
………………3.20
………………3.21
….…………3.22
(substituting g=1)
but
So
Hence ………….3.23
For a gm mole
…………….3.23 (a)
………..3.24
3.5.2 Entropy:
…….…….3.25
(B is a constant)
(S-B)/k =
Since
…………..3.26
Statistical formulation of third law can be made in the following way. We know.
……………..3.27
Hence
and
so
dT
Here is statistical weight factor of the lowest energy state. For a perfect crystal at 00 K
there is only one arrangement, which means . Hence the entropy which is third
law
Now ………3.28
………….3.29
In classical distribution of particles are distinguished from one another. If two particles
interchange their positions (or energy states) a new complexion would occur. But in Bose –
Einstein statics we start will identical particles which are indistinguishable. Therefore the
interchange of two particles between two energy states will not give rise to new complexion
or micro-state.
For example, suppose then are four particles distributed between two cells, say x and y, three
in x cell and one in Y cell. Then would give rise to four complexions if
cell y d c b a a
They are distinguishable (Boltzmann) and only one complexion, if indistinguishable (Bose
Einstein)
Suppose there are indistinguishable particles. Consider with energy with degeneracy
These particles have to be distributed among states, we shall need partitions
to place the particles in sections, each section corresponding to an energy state. It means
is the degeneracy of the level. The problem is to distribute particles into sections
without any restrictions. The permutation of particles and partitions simultaneously
is given by But this includes the permutations of particles and
partions amongst themselves also, as both these groups are internally
indistinguishable.
Hence the actual number in which particles may be allocated in states is given by.
…………..3.30
Using sterling’s formula and neglecting unity as compared to (since is very large) we
get
hence
As at the maximum
We have
……………..3.32
………….….3.33
Multiplying eqn. 2.32 by and eqn. 2.33 by and substract eqn. 2.31 from their sum we get
Since we have
……………3.34
Suppose there be particles of energy in which there are gi energy states. If all the
particles are assumed to be identical and that only particles can be identical and that number
of ways ni particles can be distributed among gi states in that group can be written as
The thermodynamics probability is given by after considering all the groups in the
macrostate.
Hence,
Since is constant
At the maximum
or ……..…..3.35
……………..3.36
…………..3.37
Since and and are arbitrary constants, the terms with in each of the parantheses
would equate to zero.
so,
or …………….3.38
Equation 2.39 gives Fermi-Dirac distribution law. Non or is large, both Bose-Einstein
and Fermi-Dirac distribution are rendered into the form
In the above discussion we have considered systems in which there are n distinguishable
particles distributed in different energy levels and f represents the molecular partition
function. But we have to deal with large scale systems, in which quantities are taken in gm-
moles.
Let us take a system which may be supposed to be constituted a large numbers of smaller
ones each with a gm-mole of the substance. That is, the individual unit here is not a molecule
but a mole. Let the energy levels allowed for the whole system (each one mole) be ,
…… and let there be such systems in all. In each of these systems there are, of course N,
Avogadro number of molecules. Then, the average energy E, of each such unit is
Now applying the same procedure in equation 3.24 the energy of the system will be
…………….3.39
This equation is useful for calculating entropy of the system, provided partition function is
known.
1. Define probability?
2. Derive Boltzmann Distribution equation.
3. What is Partition function.
4. Give the relation between partition function and thermodynamic functions.
5. Define the third law of Thermodynamics.
6. Define Bose-Einstein statistics.
7. Define Fermi – Dirac statics.
8. What is molar partition function?
9. Derive third law of Thermodynamics.
CONTENTS
4.0 Introduction
4.1 Various Types of Partition Functions
4.2 Partition Function for Translational Energy
4.3 Internal Energy of a Monoatomic Gas
4.4 Sackur Tetrode Equation
4.5 Partition Function for Diatomic Gas Molecules
4.5.1 Rotational Partition Function
4.5.2 Vibrational Partition Function
4.6 Electronic Partition Function
4.7 Equipartition of Energy
4.8 Relation between Partition Function and Thermodynamic Functions
4.8.1 Partition Function and Internal Energy
4.8.2 Partition Function and Entropy
4.8.3 Partition Function and Work Function (A)
4.8.4 Partition Function and Pressure
4.8.5 Gibbs Function and Partition Function
4.8.6 Enthalpy and Partition Function
4.8.7 Heat capacity ( ) and Partition Function
4.8.8 Equilibrium Constant and Partition Function
4.9 Heat Capacity of Solids
4.10 Quantum Theory of Heat Capacity of Solids
4.11 Residual Entropy
4.12 Terminal questions
4.0 INTRODUCTION
Statistical thermodynamics analysis is facilitated through the use of the partition function.
This is great analytical tool. It is a fact that most often the molecules are associated with
energy of different types. In expressing the partition functions, all these forms of energy must
be taken into account.
The allowed energy-levels of the t types are , , …….., for r type are , ,
………., and so on for other types.
E= …………..
= ……….
= ….….…..4.1
If is the net molecular partition function, when every molecules is capable of freely taking
any number of types of energy, we can write
E= In f ………..4.2
f= , , …….. …...…..4.3
Now we can state that net partition function is the product of the separate partition functions
in respect of different types of energy independently associated with the molecules.
Let us consider motion of the particle along x- direction. The molecular partition function
due to translation in that direction is given by.
= .………..4.4
Here is statistical weight factor which is unity in this case for each translational level and
is the translation energy of molecule in x-direction.
So =
The translational energy is also quantized, though the differences in energy levels are quite
small. The particle would exhibit periodicity in translational motion only due to impacts on
the walls of the container. Consider a cube having length l. The molecules moving along x-
direction will repeat its motion each time it complete a path 2l. If is the wavelength
associated with the matter wave, we have
According to de Broglie, there is associated with any particle of mass moving with a
Therefore
Hence
……………….4.5
We get
……………4.6
Since energy levels are very close, the variation of energy may be taken to be continuous and
hence the summation may be replaced by integration.
Put
The value of
So …………..4.7
If we consider the translational motion in all the three directions (x, y and z), the molecular
partition function for translation would be
V ..…………..4.8
So,
……………4.9
Internal energy of a monoatomic gas can be obtained by making use of equation 4.9. In the
case of monoatomic gas, rotational and vibrational energies are absent and volume V is
constant.
……………..4.10
Now
Since , we have
or ….……….4.11
The rotational energy for diatomic molecules considered as rigid rotator is given
by
= moment of inertia
As axis of rotation is defined by two coordinates, which means that there are two
rotational degree of freedom. Each quantum level of rotational will bring in two possible
modes of distribution of rotational energy. Thus the statistical weight factor ( ) in the
rotational level j is given by (2j+1). Hence the rotational partition function can be written as
……………….4.12
As the magnitude of I is comparatively large and the rotational levels are close we
may write equation 4.12 in the following form
….………..4.13
Where
Now suppose
Then
…………..4.14
This value of is valid for hetero nuclear molecules like HCl, NO etc. For homo
nuclear molecules like O16 O16, N14 N14 etc where the molecules when reversed becomes
indistinguishable from initial state; in that case the partition function is to be divided by the
number of symmetry 2. That is
……………..4.15
…………………..4.16
Now the rotational energy of a mole of a diatomic gas is given by the relation
……………….4.17
Hence …………………….4.18
……)
………………4.19
….………….4.20
weight factor greater than unity. In addition to this, there may be some excited electronic
states whose energy may be considerably greater than that of the ground states. As we
increase the temperature, such excited states become more and more occupied. Hence
electronic partition function is greater than unity and varies with temperature. The statistical
weight factor for each electronic level, normal or excited is equal to (2j+1), where j is the
resultant quantum number of the atom in the given state.
Here is the energy of the electronic state is excess of the lowest state i.e., ground
of the molecule, for each for each per molecule or per mole.
squared terms for each degree of freedom. Thus the kinetic energy is of the form ,
…….
……………..+
……………
………
Put
or
we have
f= =
Hence . . . ………
. ……..)
.Q ………………4.21
Where ………)
We know
Hence . )
………………….4.22
Equation 4.22 means that if the energy would be expressed as sum of r squared terms, the
average energy e is r x . Hence for each degree of freedom (square term) the energy is
But in some instances such conditions are not fulfilled. In the case of a gas in a
gravitational field. In this case the energy is not a square term but mgh, the mean
gravitational potential energy is not . Similarly in these cases where energy taken up is
quantised, the energy associated with them is not a linear function of temperature.
In the case of crystals, say of metals, the atoms have only vibration motion. The
vibration may occur in three directions and each mode has two degree of freedom. Hence
x 3 x 2=3kT
or …………….4.23
This is known as Dulong-Petit’s law. Equation 4.23 means heat capacity should be
independent of temperature. But considerable variation with temperature has been noticed.
That is, the principle of equipartition of energy has failed.
……..
……………..4.24
But we know
Hence
So …………………4.25
……………..4.26
…………………4.27
…………..4.28
Then ……………4.29
(W=mathematical probability)
(where B is constant)
Since ……………..4.30
We have
As
……………….4.31
We get
Since
or
We have
Hence
So ………………..4.32
And
; As ..….…………..4.33
Since ………………4.34
or
Since
So
………………….4.35
We know
By definition
.…………..4.36
By definition
….…………4.37
The equilibrium constant of a chemical change can also be evaluated from the
partition function of the reactants and products.
So
Hence
Since
Hence
……………….4.39
ground state at . Usually the zero of energy to which all other energies are referred
is the energy of dissociated atoms infinite distance apart, at 0K. The molecule at 0K would
possess some energy which is zero-energy, .
0 0 0 0
x x x x
0
x x
0 0 0 0
Where and
0
Putting x =f`
We get
………………4.40
Put and at standard state, where the value of G will be we have the
expression
……………..4.41
Now let us consider a general gaseous equation in which a moles of A and b moles of B
reacting to give c moles of C and d moles of D, we have,
cC+dD
In this equilibrium
.
…..………..4.42
Here
and
but we know
……………..4.43
….…………4.44
We have discussed above that if the energy could be expressed as sum of r squared
terms, the average energy E is rx1/2KT. Hence, for each degree of freedom the energy is ½
KT. This is law of equilibrium of energy.
But in some instances such conditions are not fulfilled. For example in the case of gas
in gravitational field, where potential energy is not a squared term but mgh, the mean
gravitational potential energy is not ½ KT. Similarly, in these cases where energy taken up is
quantized, the energy associated with them is not linear function of temperature.
We have seen in our earlier discussions that in crystals, the atoms have only vibratory
motion. In that case vibration may occur in three directions and each mode has two degree of
freedom.
Hence
or
The above equation indicates that heat capacity is independent of temperature. However,
considerable variation in heat capacity of solids with temperature has been noticed, that
reveals the failure of the principle of equilibrium of energy.
The atoms in the crystalline solids are arranged in one set pattern and have only
vibratory motion, the electronic levels are assumed to be excited. For each atom, there are
three vibrational degrees of freedom, therefore, for one mole of crystal of atoms i.e., one
gm. Mole, there are degrees of freedom. The vibrational motion is considered harmonic,
so that an oscillating atom at the nth level will have the energy,
Where is the fundamental frequency? Einstein (1907) proposed the quantum treatment of
the specific heat of solids. He considered that N atoms of the crystal were all independently
vibrating with common frequency . Then the partition function for the whole crystal is
given by,
……………4.45
We have earlier calculated the vibrational energy of a mole of diatomic gas equal to
..…………..4.46
….…………4.47
This is known as Einstein’s equation. At fairly high temperature, it is clear from the above
equation that Cv=3R, which is classical result of Dulong and Petit’s. This led Debye to
produce certain modifications in the assumptions of Einstein.
Debye suggested that all the atoms were not oscillating with the same frequency as
suggested by Einstein. He pointed out that the atoms were vibrating with a spectrum of
modes of vibrations ranging from (frequency at ground state) to some maximum
frequency . The difference between and is so large that it may assume that is
practically zero. The modes of vibrations are limited in number.
Here C is average velocity of elastic waves in the crystal. The total modes of vibration
in the crystal is 3N. Hence,
Then ….…..……..4.48
We have therefore
…....………..4.49
The paritition function of the whole crystal may be obtained by attaching to each
frequency the partition function of the single oscillation i.e.,
…..…..4.50
At very low temperature the second term is very small and can be neglected. We have then
…………….4.51
Now ..…….………4.52
……………4.53
Equation 4.53 is known as Debye’s specific heat relation at low temperature. Expressing
…………….4.54
Since for a given substance the value of is constant, hence varies with at
low temperature. This law is known as Debye’s law.
Entropies calculated using third law of thermodynamics are called thermal entropy.
However statistical entropies are more rigorous. It is found that the thermal entropies are
smaller than the statistical entropies and deviate ranging from 3.1 to 4.83 JK-1 mol-1. We thus
conclude that entropies of substance (such as H2, D2, NO, N2O, H2O etc.) are not zero at 0K
in accordance to third law, but finite. These entropies are called residual entropies.
CO CO CO CO CO CO OC OC
CO CO CO CO CO OC OC CO
NNO NNO NNO NNO NNO NNO ONN ONN
NNO NNO NNO NNO NNO ONN ONN NNO
(a) (b)
Fig. 4.1 Alternative molecular arrangements (a) perfect crystal (b) actual crystal
Since both the arrangements are equally likely. From Boltzmann entropy equation,
S=klnW
We get
S=kln2N
=Nkln2
=Rln2
=5.76K-1 mol-1.
As the residual entropies are found experimentally to be less than this value, it is
evident that the two alternative orientations of CO and N2O molecules in the solid state at 0K
are not completely random.
The occurrence of residual entropy in ice is also attributed to the molecular disorder at 0K. In
ice crystal, each ice molecule is tetrahedrally coordinated to four other water molecules by
hydrogen bonds. This enables that crystal to exist in many configurations at 0K. For H2 and
D2 the thermal entropies at 0K are less than the corresponding statistical entropies. The
calculation of the statistical entropy assumes that there exists equilibrium between ortho H2
and para H2 at all temperatures. The of ortho H2 and para H2 is found to be 18.37 JK-1
mol-1 in the vicinity of 0K. when this value is added to the thermal entropy (calculated from
heat capacity measurements), the agreement with the statistical entropy is very good.
CONTENTS
5.2 Introduction
5.6 Summary
5.8 Glossary
5.10 References
5.1 OBJECTIVES
5.2 INTRODUCTION
The present unit deals with theories of reaction rate that includes different theories like
collision theory, transition sate theory and theories of unimolecular reactions. Unimolecular
theories include Lindemann’s theory, Hinshelwood’s theory and RRKM theory. In addition
to these theories of reaction rate, we will discuss about potential energy surface.
Rate of reaction is also known as reaction rate. It is defined as the rate of change of
concentration of reactants or products per unit time. As the reaction proceeds, the
concentration of reactant decreases while the concentration of the products increases.
Therefore the rate of reaction for a general reaction is given by:
A B
r = – d [A]/dt = + d[B]/dt
where (–) sign shows decrease in concentration of reactant A while the (+) sign shows
increase in the concentration of product B, infinitesimally small change in concentration is
represented by d. The rate of reaction is expressed in mole litre–1 time–1. Generally the time
–1
taken is in seconds hence the rate is expressed in mol l s –1. According to rate law, rate of
reaction is directly proportional to the concentration of the reactants, raised to some power.
For a general reaction:
A B
r = k [A]
The concentration terms in the rate law expression are raised to some power, this power is
known as order of reaction. This is the basic idea of chemical kinetics. Now let us discuss
different theories of reaction rate
The collision theory of reaction rate is developed by Arrhenius and van’t Hoff. According to
collision theory of reaction rate, reactant molecules must collide with each other that result in
a chemical reaction. There are two theories of collision: one is classical collision theory while
the other is modern collision theory. Let us discuss these theories one by one:
According to classical collision theory, gaseous reactant molecules collide with each other
that result in a chemical reaction and the rate of this reaction is directly proportional to the
number of colliding reactant molecules. Let us consider a reaction in which reactant X and Y
reacts or collides to form product XY. Therefore according to classical collision theory:
X+Y XY
Where [X] and [Y] is the concentration of reactant X and Y. Let us consider a collision
between two gas molecules C and D. The number of collision occurring between two gas
molecules per second per centimeter cube is given by Z which is known as collision number.
Z for bimolecular collision is given by ZCD:
Where nC is the number of gaseous molecule (C), nD is the number of gaseous molecule (D),
dCD is the average collision diameter, R is gas constant, T is absolute temperature (in Kelvin)
and µ is reduced mass. The average collision diameter (dCD) is given by (dC + dD)/ 2 while the
reduced mass µ is given by mC mD/ mC + mD, where mC is molecular mass of gas C and mD is
molecular mass of gas D. Hence collision number Z can be written as:
ZCD = nC nD {(dC + dD)/ 2}2 {8ΠRT (mC + mD) / mC mD}1/2 ………… (2)
k = Z e – E/ RT ………….. (3)
In order to check the validity of equation 3, reaction between hydrogen (H2) and iodine (I2) is
taken into consideration. It has been observed that at 556K, the rate constant which is
obtained experimentally is 3.7 × 10– 7. The value of rate constant using equation 3 is
calculated and is equal to 5.5 × 10– 7. The value of rate constant thus obtained is in good
agreement with the experimental value that proves the validity of the rate constant expression
for bimolecular collision between hydrogen and iodine. Similar agreement has also been
observed in case of reaction or collision between nitrogen and oxygen.
Simple collision theory was modified by Arrhenius and a modified collision theory was
proposed. The theory considers several points:
• All collisions among the molecules are not responsible for causing a chemical
reaction. Only those collisions are effective (workable) which possess a minimum
amount of energy known as activation energy for a chemical reaction to occur. It has
been observed that when the colliding molecule possess energy equal to or greater
than activation energy then only the chemical reaction takes place otherwise not . It
means that when energy of colliding molecules is lower than the activation energy
then no chemical reaction takes place. Therefore in the rate expression for
bimolecular collision (equation 3), a new term called steric factor or probability
factor is introduced. Therefore equation 3 is written as:
k = PZ e – E/ RT ………… (4)
where P is the probability factor or steric factor having value from 1 to 10– 8. It has been
found that probability factor is related to entropy change by the expression:
P1/P2 = e ∆S/ R ……… (5)
Where P1 is probability factor for forward reaction and P2 is probability factor for backward
reaction (condition is that the reaction is reversible).
• The colliding molecules should must possess a proper orientation required for the
reaction to occur. In simple words, we can say that only effective collision (collision
having correct orientation) leads to a chemical reaction. The proper orientation means
that there must be direct contact among the atoms which are involved in bond
breaking and bond forming. This will become more clear by Figure 1 given below:
Though the theory is applicable in case of simple gas molecule but it suffers certain
limitations:
• Classical collision theory do not consider the internal motion of the molecules taking
part in a collision which leads to a chemical reaction.
• In order to relate the probability factor with the properties of the colliding reactant
molecules, several attempts were made but satisfactory results were not obtained.
• As modified collision consider the concept of classical collision theory which is based
on classical mechanics that require modifications.
• Only kinetic energy of the reactant molecules is considered not the vibrational or
rotational energy. The reason of not taking this energy is not clear.
5.3.2 Transition sates theory or activated complex theory
Transition state theory (TST) is also known as activated complex theory (ACT) or absolute
reaction rate theory. This theory was proposed by Pelzer and Minger but later on by Eyring,
Evans and Polanyi. As by using this theory, we can find the absolute value of rate constant by
using physical properties therefore the theory is also known by name absolute theory. This
theory considers the quantum mechanics instead of classical mechanics. According to
transition sate theory, collisions among the reactant molecules do not lead to chemical
reaction but there is change in potential energy of the reactants as it changes when reactants
approaches from infinite distance in order to react and then forms products. Reactants having
sufficient energy from infinite distance approach each other and form an activated complex
which further forms the products. The reaction is given below:
For example, there are two reactants X2 and Y2, they reacts to form an activated complex
which together forms products.
Activated complex is a special molecule that possesses highest energy or most unstable
molecule that possesses three translational degree of freedom in which forth degree of
freedom is along the reaction coordinate. It has transient existence. In activated complex, as
the bond among the atoms are weak hence one vibrational degree of freedom is less than the
normal molecule. This vibrational degree of freedom changes to translational degree of
freedom hence possess forth degree of freedom (translational). On the basis of transition state
theory, energy profile diagram is given in Figure 2.
It is clear from Figure 2 that energy of products is less than energy of the reactants, therefore
E products –E reactants is negative. This shows that change in energy is negative which indicates
that during the formation of products, some energy is released (exothermic). Figure 2 also
shows that activation energy (Ea) is the difference between the energy of the activated
complex and energy of the reactants. Similarly for endothermic reaction, the energy profile
diagram is given below in Figure 3.
Now we will study the formulation of activated complex theory in terms of statistical
mechanism and thermodynamics. Let us discuss these formulations one by one.
Where X2Y2* represents activated complex. It is clear from the above mechanism that rate of
reaction (decomposition) is directly proportional to the concentration of the activated
complex [X2Y2*]. The rate of decomposition of activated complex also depends upon the
number of activated complex that cross the energy barrier per unit time. This is called
frequency of crossing the energy barrier and is given by the ratio of average velocity of
activated complex to the length (d) at the highest energy level (i.e. activated complex energy
state). According to kinetic theory of gas, the average velocity (vavg) is given by:
Therefore on multiplying all these factors, rate of reaction (decomposition) is given by rd:
rd = (kT/h) × [X2Y2*]
Here kT/h refers to universal frequency. Now it is clear that activated complex is formed
from reactants therefore rate of formation (rf) is directly proportional to the concentration of
the reactant X and Y. Therefore:
rf = kf [X] [Y]
Where θ*, θX and θY are the partition function of reactants X, Y and activated complex
respectively. – ∆Eo is the change in energy (zero point). On putting the value of K* in equation
(9), we obtain:
By using thermodynamics, we can express the rate constant of the reaction (mechanism)
which is according to activated complex theory. This expression of rate constant containing
thermodynamic functions is called thermodynamic formulation of reaction rate or activated
complex theory. On considering the general mechanism:
As discussed earlier that in statistical formulation of reaction rate, rate constant is given by:
Where T is absolute temperature and R is gas constant. On applying this reaction isotherm on
first step of the mechanism given by activated complex theory we obtain:
From the definition of free energy, when there is change in the state of the system, change in
free energy is given by relation: ∆G = ∆H – T ∆S
ln K* = – (∆H* – T ∆S*) / RT
K* = e – (∆H* – T ∆S*) / RT
Equation 13 is known as Wyne- Jones and Eyring equation or simply Eyring equation.
Equation 13 represents thermodynamic formulation of reaction rate in which rate expression
contains thermodynamic functions. For gaseous reactions
∆H* = Ea – RT
On putting the value of ∆H* in Eyring equation, we obtain: kf = (kT/h) × e – (Ea – RT) / RT × e ∆S*
/R
PZ e – E/ RT = (kT/h) × e – E / RT × e ∆S* / R
∆S* / R
PZ = (kT/h) × e
∆S* / R
P = (kT/hZ) × e
The value of probability factor is found to be in between 10–5 to 10–10. The value of P is less
than the value obtained from collision theory which shows that reaction proceeds more
slowly on the basis of activated complex theory. ACT provides the physical significance of
probability factor by relating it to entropy of activation ∆S*.
It has been found that for reactions between simple molecules, both collision theory and
activated complex theory give same results but when reaction occurs in between complex
molecules, activated complex theory give satisfactory result while collision theory fails.
Before coming to the theories of unimolecular reactions, let us discuss about the unimolecular
reactions. Reactions in which one reactant molecule forms the activated complex or one
reactant molecule forms the products. As the name indicate unimolecular which indicates one
reactant molecule. Unimolecular reactions are classified into two types: elementary
unimolecular reaction and free radical unimolecular reactions. Let us discuss these reactions
one by one:
• Elementary unimolecular reactions: The reaction in which products are formed from
reactants in a single step that involve a formation of an activated complex. For
example, formation of propylene from cyclopropane.
• Free radical unimolecular reactions: In these reactions as the name suggests, free
radicals are formed as an intermediate instead of single activated complex. Several
decomposition reactions comes under this category.
Now coming to theories of unimolecular reactions. Various theories are there which explain
the mechanism of unimolecular reactions which include Lindemann’s theory, Hinshelwood
theory and RRKM theory. We will now discuss all these mechanism or theories on by one:
This theory was proposed by F.A Lindemann and subsequently supported by Christiansen
According to Lindemann, the following mechanism is proposed for conversion of A into
Products
We know that rate of reaction is proportional to the concentration of the energized molecule
as in this step formation of products takes place. Therefore:
r = k3 [A*]
On putting the value of [A*] from above expression we obtain: r = k3 k1 [A]2 / k2 [A] + k3
For reactions taking place in gaseous phase, the concentration term is replaced by pressure
terms. Hence: r = k3 k1 P2 / k2 P + k3
r = k3 k1 P / k2
This expression shows that at high pressure, rate of reaction is directly proportional to
pressure of a gas raised to power one or we can say it obeys first order kinetics.
Though the theory given by Lindemann is quite satisfactory in explaining the mechanism of
unimolecular reactions but there are some drawbacks:
Graph is plotted between 1/ k versus 1/ [A], result straight line with slope 1/k1 and intercept
k2 / k3 k1. But instead of straight line there is deviation from it. Lindemann’s drawback leads
Hinshelwood to put forward his theory. Let us now discuss Hinshelwood’s theory of
unimolecular reactions.
where Z is collision number, S is number of vibrational degree of freedom the value of which
is obtained by hit and trial method, k is Boltzmann constant , T is absolute temperature and E
is energy of activation. For the molecule having one vibrational degree of freedom, the above
expression of k reduced to:
k = Z e –E/ kT
It has been found that for large value of S, pre-exponential factor (A) depends strongly on the
temperature. Though the theory is successful but there are certain limitations:
• The theory does not show the dependence of pressure on rate constant.
• From Lindemann’s theory on which Hinshelwood theory is based, the variation of 1/
k versus 1/ [A] should be straight line but there are deviations in the curve which was
not explained by Hinshelwood.
5.4.3 RRKM theory
RRKM theory of unimolecular reaction stands for Rice- Rampsberger- Kassel- Marcus
theory. This is one of the most refined and successful theory of reaction rate. According to
RRKM theory which on the basis of activated complex theory considered all vibration of
energized molecule (A*) as vibrational frequency and all rotation of energized molecules as
rotational constant. RRKM proposed the mechanism which is as follows:
Where N c(ɛc) is the energy density of the degree of freedom that contribute to bond
breaking, fc is the partition function where c represents critical energy contribution and ɛc is
critical energy of the activated complex which is more than zero point energy. The critical
energy ɛc is equal to the sum of energies of motion along the reaction coordinate and energy
associated with vibrations and rotations. It has been found that at high pressure, RRKM
theory give the same expression of rate constant as given by activated complex theory.
Therefore we can say that at high pressure, RRKM theory get reduced to transition state
theory or activated complex theory. At high pressure, thermal Boltzmann distribution is
maintained at all energies. At high pressure, the expression for rate constant is:
Potential energy surface is short is also known as PES. The concept of PES is applied to
several branch of science that includes physics, chemistry. The properties of structure
comprising of atoms are predicted by applying the concept of PES. This surface represents
the energy of a system which is generally based on the position of the atoms. Energy profile
diagram also refers to potential energy diagram or PES represents the variation of energy
with reaction progress. In potential energy diagram or curve, there is a point at which
potential energy is maximum along the reaction coordinate. This point is known as saddle
point which corresponds to transition state. There are stationary points in the curve that
represents stable chemical species or species with lowest energy. On the basis of bond length
extension in the activated complex in comparison to bond length of the reactants or product,
potential energy surface is classified into two types: attractive potential energy surface and
repulsive potential energy surface. In order to define these types, let us take a general
reaction: X + Y-Z X-Y + Z
The bond length of newly formed bond (rXY) is equal to the difference of the bond length of
XY in the transition state ((rXY)* and bond length of XY in the product molecule (r1XY) i.e.
rXY = rXY* – r1XY . In this reaction there is cleavage of Y-Z bond thus rYZ = rYZ* – r1YZ, where
r1YZ refers to bond length of Y-Z bond of the reactants. Potential energy surface is said to be
attractive when rXY > rYZ for those reactions in which heat is evolved. Here transition state is
achieved by bringing the reactants close to each other. For example, K and Br2 to form KBr
and Br. Potential energy surface is said to be repulsive when rXY < rYZ for those reactions in
which heat is evolved. Here transition state is achieved by the separation of the products. For
example, H and Cl2 to form HCl and Cl. The potential energy diagram is given in Figure 4 in
which AB represent lowest energy path. Point o gives the potential energy of these atoms X,
Y and Z which are well separated from each other.
For the reaction in which energy is absorbed, the PES is said to be attractive when
translational energy of the reactants is effective. The PES is said to be repulsive when
vibrational excitation of the reactants is effective. The concept of PES is very useful in
determining the geometry of a molecule and kinetics of a chemical reaction by considering
the position of the atoms.
5.6 SUMMARY
In this unit we have discussed the theories of rate of reactions that includes collision theory
and activated complex theory. Collision theory is discussed under two headings: classical
collision theory and modern collision theory. In activated complex theory, we discussed its
statistical and thermodynamical formulation of reaction rate. Then we move forward to
unimolecular reactions where we studied about the theories of unimolecular reactions which
include Lindemann’s mechanism, Hinshelwood mechanism and RRKM mechanism. In
addition to all these, we also studied about the potential energy surface that plays an
important role in determining the molecular structure.
5.8 GLOSSARY
• Activated complex: Unstable molecule having highest energy in the energy profile
diagram.
• Potential energy surface: Surface associated with the position of atoms that
represents the energy of the system.
• Activation energy: Minimum energy that the reactants molecules must have in order
to form products.
• Steady state approximation: Rate of formation and rate of dissociation is same for
the intermediate.
5.10 REFERENCES
1) Raj, G. (2017), Advanced Physical Chemistry, Goel Publishing House, Meerut.
2) Puri, B. R., Sharma, L. R. and Pathania, M. S. (2012), Principles of Physical
Chemistry, Vishal Publishing Co., Delhi.
3) Bahl, A., Bahl, B. S. and Tuli, G. D (2014), Essentials of Physical Chemistry, S.
Chand & Company Pvt. Ltd., New Delhi.
4) Gurtu, J. N. and Gurtu, A (2015), Advanced Physical Chemistry, Pragati Prakashan,
Meerut.
CONTENTS:
6.1 Objectives
6.2 Introduction
6.4.2 Temperature
6.4.3 Catalyst
6.4.7 Pressure
6.8 Summary
6.10 Glossary
6.12 References
6.1 OBJECTIVES
6.2 INTRODUCTION
In the lower classes, we studied about the chemical kinetics which deals with the rate or
speed with which chemical reaction takes place. For example, corrosion is a very slow
process or have low rate of reaction. By using chemical kinetics, the mechanism of the
chemical reaction can be determined. There are several factors that affect the rate of chemical
reaction occurring in solution, gaseous phase or condensed phase. Solution as we know
consist of a solute and a solvent. A solute is present in small amount while the solvent is
present in large amount in the solution. For example, solution of sugar in water.
The present unit deals with the study of solution kinetics or the study of rate of rate with
which solute reacts with the solvent or kinetics of reactions in solution. In this unit, we will
discuss about different factors that affects the rate of reaction occurring in solution, the
concept dielectric constant, ionic strength followed with primary and secondary salt effect.
The chapter is quite interesting as it deals with the interaction of solvent with the solute in
solution thereby affecting the reaction rate. In order to understand the topic more clearly, we
will go through the entire chapter in order to study the factors affecting the kinetics of
reactions occurring in solution.
Rate of reaction is also known as reaction rate. It is defined as the rate of change of
concentration of reactants or products per unit time. As the reaction proceeds, the
concentration of reactant decreases while the concentration of the products increases.
Therefore the rate of reaction for a general reaction is given by:
A B
r = – d [A]/dt = + d[B]/dt
where (–) sign shows decrease in concentration of reactant A while the (+) sign shows
increase in the concentration of product B, infinitesimally small change in concentration is
represented by d. The rate of reaction is expressed in mole litre–1time–1. Generally the time
–1
taken is in seconds hence the rate is expressed in mol l s –1. According to rate law, rate of
reaction is directly proportional to the concentration of the reactants, raised to some power.
For a general reaction:
A B
r = k [A]
The concentration terms in the rate law expression are raised to some power, this power is
known as order of reaction. This is the basic idea of chemical kinetics. Now let us discuss the
factors that affect the rate of reaction.
6.4.2 Temperature
We know that temperature is directly proportional to average kinetic energy. More is the
kinetic energy, more is the speed of the reacting molecules and more is collision which
thereby increases the rate of reaction. It has been observed that for every 10oC rise in
temperature, the rate of reaction increases by a factor two or three. At two temperature that
differ from each other by 10oC, the ratio of the rate constant of a reaction is called
temperature coefficient of a reaction.The value of temperature coefficient is 2 or 3 as
discussed above.
Where k represents rate constant at two temperature that differ from each other by 10oC.
During the decomposition of HI, the rate of reaction increases by a factor of 1.7 when
temperature increases by 10oC.
6.4.3 Catalyst
These are certain substances that alter the rate of reaction. These substance generally
increases the rate of reaction without being consumed in the reaction and are termed as
catalyst. These substances remain unchanged in chemical composition and in amount at the
end of the reaction. The addition of a catalyst in a reaction increases the rate by lowering the
activation energy of the reaction due to which more number of reactant molecules cross the
energy barrier and form the products as there is an energybarrier between the reactant and the
product. Activation energy is the minimum amount of energy that the reactants molecules
must possess in order to form products. It is represented by Ea. For example, in the
decomposition of potassium chlorate, MnO2 (Manganese dioxide) acts as a catalyst which
speeds up the decomposition. The plot explaining the effect of catalyst on activation energy is
given in Figure 1.
There are certain reaction that takes place in presence of light. Such reactions are called
photochemical reactions. In these reactions, photons are absorbed by the reactant molecules.
Hence, more is the no. of photons (light radiations) provided to the reactant molecules, more
is the energy to cross the energy barrier which results in the increase of reaction rate. For
example, the rate of photosynthesis increases when more photons are provided to the reactant
molecules or when day is sunny.
In case of heterogeneous reaction when reactants are not in same phase (like one reactant is in
solid phase and the other one is in gas phase), larger is the surface area of the reactant, more
will be the rate of reaction. This is because the reactions occur only at the surface of the
reacting species. For example, the rate of reaction for a lump of a solid is less than the rate
where it is present in powdered form as powder occupies more surface than the lump. The
rate of homogeneous reaction (both phase are same) is more than the rate of heterogeneous
reaction.
Ionic compounds react with faster rate than the covalent compounds because the later require
energy to break the bond. For example, when silver nitrate (AgNO3) is added to sodium
chloride solution (NaCl), silver chloride (AgCl) is formed as a precipitate immediately.
6.4.7 Pressure
If the reaction takes in gaseous phase, the increase in pressure increases the number of
collision among the gaseous molecules, thereby increasing the rate of reaction. As in case of
gaseous reaction, according to rate law, rate of reaction is directly proportional to the partial
pressure of the reactant species.
In case of solution, we consider the effect of solvent on rate of reaction. The viscosity and
dielectric constant of a solvent affect the reaction rate. Before discussing its effect on reaction
rate, let us have some common knowledge about the viscosity and dielectric constant.
Viscosity of a liquid is defined as the resistance in the flow of liquid. Liquid consist of layers
of molecules that are arranged one over the other. This exerts resistance in the flow of a
liquid. Viscosity is inverse of fluidity. Fluidity is the ability of a liquid to flow. Viscosity is
represented by ƞ. The unit of viscosity is Poise (P). With increase in temperature, viscosity of
liquid decreases. Dielectric constant of a solvent is defined as the capacity of the solvent to
weaken the force of attraction (electrostatic) between ions (solute) and to separate them.
Dielectric constant is represented by ɛ.Water which is regarded as a universal solvent has
dielectric constant value equal to 79.8. Due to high value of dielectric constant, when solute
is added into the water, the water weakened the force of attraction by which solute molecules
are attached resulting in the free movement of ions (solute). Thus water is a strong
dissociating solvent. In vacuum, the value of dielectric constant is one for all the solvents.
As we have discussed above that dielectric constant and viscosity affect the reaction rate in
solution. Each factor is discussed in detail below:
In order to show the effect of dielectric constant (ɛ) on rate of reaction, the theory of absolute
reaction rate also known as transition state theory has been undertaken into
consideration.According to this theory, ions undergoing the reaction in a medium of dielectric
constant (ɛ) are considered as a sphere. Let us consider an ionic reaction involving ions A and
B. Ionic reactions are those reactions that occurs in between the ions. Let rA andrBbe the
radius of two ions having charge ZAe and ZBe respectively. Initially these two ions A and B
are at an infinite distance. According to the activated complex theory, when these two ions
are brought together in order to form the activated complex, the work done (W) is given by:
Where ɛ is dielectric constant of the solvent, dABis the distance between two ions in an
activated complex.
W = ZA ZB e2 / ɛ dAB
If ions possess charges of same sign (both are either positive or both are negative), the work
is positive. If ions possess charges of different sign (one ion is positive while the other ion is
negative), the work is negative. This work is equal to electrostatic distribution. The free
energy of activation (∆G*) per molecule is given by the sum of free energy that contribute to
electrostatic distribution and non- electrostatic distribution (N.E.S.).
From thermodynamics formulation of reaction rate (activated complex theory), rate constant
is related to free energy of activation (∆G*) by the expression:
k= e– ∆G*/RT
k= e– (∆G* N.E.S. + Z Z e2 N / ɛ d ) / RT
A B AB
ln k = ln k0– ZA ZB e2 N / ɛ dAB RT
When a graph is plotted between ln k versus 1/ ɛ, a straight line is obtained having slope ZA
ZB e2 N / dAB RTand intercept ln k0which is shown in Figure 2.
It is concluded that for an ionic reactions, with increase in dielectric constant of a solvent,
rate constant increases thereby increasing the rate of reaction while with decrease in dielectric
constant of the solvent, the rate of reaction decreases. This equation correlating rate constant
with dielectric constant is valid for a large number of ionic reactions.
In order to show the effect of viscosity on rate of reaction, let us consider to reactant
molecules A and B having radii rAand rB respectively. The number of collision (encounters,
υ) between A and B per cubic centimeter is given by equation given below. One thing to be
noted here is that encounter word is used in place of collision. Collision word is used in case
of gases while for liquids and solids, encounter is used because here the solute molecules get
surrounded by the solvent molecules forming a cage.
whereDA , DB are their respective diffusion coefficient and nA,nB are the number of molecules
per centimeter cube of A and B respectively.If each encounter leads to a chemical reaction
then the number of encounter is equal to rate of reaction.
For spherical particles of a liquid, diffusion coefficient is related to radius (r) of the molecule
and the coefficient of viscosity (ƞ) by Stoke ‘s equation given below:
D = kT/ 6 Π ƞr
Ionic strength of a solution is defined as the measure of electrical intensity. This electrical
intensity is due to the presence of ions in the solution. Ionic strength is represented by I. The
expression of ionic strength for two ions in solution is given by:
I = ½ (m1Z12 + m2Z22)
I = ½ ∑mi Zi2
In case of salts, ions are formed. These ions exert force on each other. Ionic strength which
has been discussed above affects the rate of reaction. This is studied under topic salt effect.
The concept of salt effect was given by J. N. Bronsted and N. Bjerrum. Salt effects are
classified into primary salt effect and secondary salt effect. Let us discuss these effects one by
one:
This effect is seen in those reactions which do not involve catalyst. In case of primary salt
effect, the ions that act as a reactant first combine to form activated complex. Some of the
activated complex gives product while some of the activated complex becomes deactivated to
give back the ions (reactants). This mechanism is according to activated complex theory. The
reaction is shown below:
WhereA and B are the ions; ZA and ZB are their valencies respectively. AB is the activated
complex having valency ZA + ZB. According to rate law, rate of reaction is proportional to the
concentration of the activated complex.
It is clear from the reaction that there is equilibrium between ions and activated complex.
Therefore the equilibrium constant (K) is given by:
K = a*/ aAaB
Where a is the activity. We know that activity coefficient (ɣ) is equal to a/ c. Therefore:
a=ɣc
k = r / [A] [B]
k = k’ ɣAɣB / ɣ*
For the solution of electrolyte, according to Debye- Huckel limiting law: ln ɣi = – AZi2 I1/2
Where I is the ionic strength of the solution, A is the constant for a particular solvent while
Ziis the valency of ith ion. Therefore:
ln k = ln k’ – A I1/2 (– 2 ZA ZB)
ln k = ln k’ +2 A I1/2 ZA ZB
In case of water as a solvent, constant A possesses value 0.509. On putting this value on
above expression, we obtain:
ln k = ln k’ +2 (0.509) I1/2 ZA ZB
Equation 4 shows that the rate constant or rate of reaction depends on the valencies of ions as
well as the ionic strength (I) of the solution. This is primary salt effect and equation (4) is
called Bronsted- Bjerrum equation. On the basis of chemical reaction involved, different case
arises:
Case 1: If ZA ZB = positive
It means that both the ions undergoing the chemical reaction possess either positive charge or
negative chargethen only the product of ZA ZB is positive. It is clear from Bronsted- Bjerrum
equation that with increase in the ionic strength of the solution, rate constant increases
thereby increasing the rate of reaction. For example:
S2O32– + I– Products
Case 2: If ZA ZB = negative
It means that one of the reactant (ion) possess negative charge then only the product of ZA ZB
is negative.It is clear from Bronsted- Bjerrum equation that with increase in the ionic strength
of the solution, rate constant decreases thereby decreasing the rate of reaction.For example:
H+ + Br–1 Products
Case 3: If ZA ZB = zero
It means that one of the reactant is non- electrolyte. It is clear from Bronsted- Bjerrum
equation that rate constant does not vary with ionic strength of the solution.For example:
This effect is seen in those reactions that involve catalyst. In order to understand the
secondary salt effect, there must be change in the concentration of the catalyst involved. This
effect is shown by equation 5, given below:
ln k = ln k’ +1.018 I1/2
This equation shows that rate constant does not depend on the valency but it only depends
upon the ionic strength of the solution. With increase in the ionic strength of the solution, rate
constant also increases thereby increasing the rate of reaction.
6.8 SUMMARY
In this unit, we have discussed about the kinetics that takes place in a solution. The topic
consists of factors that affect the rate of reaction in solution that includes dielectric constant
and viscosity. Further an idea of ionic strength is given which is then correlated with the rate
constant which was explained under primary and secondary salt effect.
2. Temperature coefficient is defined as the ratio of two ………. that differ from each other
by 10oC.
3. Catalyst ………. the rate of reaction by ………. the activation energy of the reaction.
6. Rate of reaction ………. with increase in the dielectric constant of the solvent.
7. More is the viscosity of the solvent, ……. will be the rate of reaction.
9. Primary salt effect is seen in those reactions that involves ………. reactions.
c. No effect
2. In Bronsted- Bjerrum equation, the value of constant A for water as a solvent is equal to:
a. 0.059
b. 0.509
c. 0.00509
d. 0.00
b.ln k = ln k’ +2 A I1/2 ZA ZB
c. ln k = ln k’ – 2 A I1/2 ZA ZB
b. I = ½ (m1Z12 + m2Z22)
c. I = ½ (m1Z1 + m2Z2)
d. I = 2 (m1Z12 + m2Z22)
c. Independent of viscosity.
a. D = kT/ 6 Π ƞr
b. D = kTƞ / 6 Π r
c. D = kT2/ 6 Π ƞr
d. D = kT/ 6 Π ƞr2
7. Temperature coefficient is defined as the ratio of rate constant that differ from each other
by:
a. 20oC
b. 10oC
c. 30oC
d. 05oC
6.10 GLOSSARY
• Dielectric constant - Capacity of the solvent to weaken the force of attraction present
in solute.
• Viscosity: Resistance in the flow of liquid.
• Ionic strength: Measure of electrical intensity.
6.12 REFERENCES
Bahl, A., Bahl, B. S. and Tuli, G. D (2014), Essentials of Physical Chemistry, S. Chand &
Company Pvt. Ltd., New Delhi.
• Describe in detail the various factors that affect the rate of reaction.
• Explain the effect of dielectric constant on the rate of reaction taking place in
solution.
• Viscosity of a solvent affects the rate of reaction in solution. Explain.
• Explain in detail the primary and secondary salt effect.
CONTENTS:
7.1 Introduction
7.2 Objectives
7.3 General idea of kinetics of fast reactions
7.3.1 Relaxation method
7.3.2 Flow method
7.3.3 Flash photolysis
7.4 Molecular reaction dynamics
7.5 Chain reaction
7.5.1 Kinetics of hydrogen –bromine chain reaction
7.5.2 Kinetics of hydrogen –chorine chain reaction
7.5.3 Pyrolysis of acetaldehyde
7.6 Enzyme reactions
7.6.1 Kinetics of enzyme reactions
7.6.1.1 Michaelis- Menten equation
7.6.1.2 Lineweaver- Burk equation
7.6.1.3 Eadie- Hofstee equation
7.7 Summary
7.8 Self-Assessment Questions (SAQ)
7.9 Glossary
7.10 Answers to SAQ’s
7.11 References
7.12 Terminal Questions
7.1 OBJECTIVES
After reading this unit you will be able to:
• Describe various methods of studying the rate of fast reactions that includes relaxation
method, flow method and flash photolysis method.
• Write note on molecular reaction dynamics.
• Explain the kinetics of hydrogen-chlorine reaction and hydrogen-bromine reaction.
• Give the mechanism of pyrolysis of acetaldehyde that includes its kinetics also.
• Describe the kinetics of enzyme catalyzed reactions.
7.2 INTRODUCTION
In earlier chapters we have studied about chemical kinetics that deals with moderate
reactions. Moderate reactions are those reactions that occur not very fast or not very slow like
acid-catalyzed hydrolysis of esters. In fast reactions, reaction takes place so fast in such a
short time (few second or less) than required to mix the reactants. For these reactions, rate is
not obtained by conventional methods. Some other methods are used for determining the rate
which are mentioned below in further topics.
The present unit deals with kinetics of fast reactions in which we will go through various
methods for studying the rate of fast reactions that includes relaxation method, flow method
and flash photolysis. In addition to these, we will have a general idea of molecular reaction
dynamics. On moving forward, we will study about the kinetics of chain reactions that
includes hydrogen- chlorine reaction, hydrogen – bromine reaction and pyrolysis of
acetaldehyde. Finally we will conclude on the topic, kinetics of enzyme reactions.
The entire unit is quite interesting as after reading the unit we will be able to know how the
rate of fast reactions are determined that takes place in few seconds. We will also came to
know about the kinetics of reaction catalyzed by enzymes which are biological catalysts. Let
us now discuss these topics one by one.
Relaxation method is one of the methods for determining the rate of fast reactions that takes
place in solution. Eigen and his coworkers put forward this method. In this method, there is a
sudden change in some of the physical parameter of a system that is in equilibrium. This
physical parameter includes temperature, pressure and electric field. Due to this change in the
system there is change in state of a system or we can say that there is change in temperature,
pressure or electric field. Relaxation is defined as a time required attaining the equilibrium
again after this change. It has been observed that there is a change in a solvent structure
which occurs quickly than the chemical reaction which is under consideration. By this
–10
method, reactions having half- life ranging from few seconds to 5 × 10 second are
measured. Let us now derive the rate constant expression by considering a reversible first-
order reaction:
Suppose at any instant, the total concentration of the reaction that consist of A and B is
represented by a while x be the concentration of B. As rate of reaction is defined as the
change in the concentration of reactant or product per unit time, therefore rate of formation of
product (r) is given by:
k1 (a – xe) = k2 xe
d(∆x) /dt = k1 a – k1 ∆x – k1 xe – k2 ∆x + k1 a – k1 xe
d(∆x) /dt = – k1 ∆x – k2 ∆x
Where kr is the relaxation rate constant which is equal to k1+ k2. Relaxation time (t) is a
reciprocal of relaxation rate constant, hence:
t = 1/ kr = 1/ k1+ k2
Relaxation time is obtained by different methods that includes pressure - jump method and
temperature- jump method which measures relaxation time of more than 10 –5 second whereas
for measuring relaxation time less than 10 –5 second, high electric field methods are there.
This method is one of the commonly used method for determining the rate of fast reactions
and was the first method to be employed for studying the kinetics of fast reactions having
half-life in the range 10 to 10–3 seconds. This technique was proposed by Roughton and
Hastridge in 1923. There are various flow methods that include stopped-flow method,
continuous –flow method and quenched flow method. Out of these three methods, stopped-
flow method is commonly used method in comparison to continuous –flow method and
quenched flow method. Let us discuss these methods:
• Stopped- flow method: As discussed above, this is commonly used flow method. The
apparatus for this method is shown in Figure 1.
It is clear from Figure 1 that apparatus consist of two reservoir having reactants (reagents) R1
and R2. These two reservoirs is attached to a piston. By using piston, two reactants are
allowed to mix in a chamber known as mixing chamber indicated by M in Figure 1. From M,
the solution of two reagents then move through observation tube having light source in the
ultraviolet and visible region. After mixing the reagents, observations are made within a few
milliseconds. In this method, with different flow from the value of absorption, rate versus
time graph is plotted. In this method as the name indicates, the flow is stopped suddenly at a
particular point and the equilibrium rate is measured as a function of time by using absorption
data. One point to be noted is that in this method, measurements are done at a particular
point. For example, aqueous solution of Fe3+ and CNS– is studied by using this method.
In this method, as the name indicates a strong flash of high intensity or energy of the order
105 J is allowed to fall on a sample for microseconds or 10–6 second. On the interaction of
high energy flash with the sample, almost all the molecules of the sample get excited that
results in the formation of free radicals (reaction intermediates). The concept of flash
photolysis was given by R. G. W. Norrish and G. Porter (1950). For this great achievement,
they were awarded by Nobel Prize. As free radicals are produced on the exposure of flash on
the sample therefore this photolysis is used from determining absorption spectra of free
radicals like •CH3, •Cl and •NH2 for studying reactions occurring in solution as well as in gas
phase. The apparatus for flash photolysis is shown in Figure 2.
It is clear from Figure 2 that there is a high voltage supply that provide 10,000 volt to a
condenser. In the spark gap, a spark is generated by a trigger signal which is responsible for
passing current through the flash lamp above which sample is placed. Within few
microseconds, the condenser gets discharge resulting in the automatic stop of flash lamp.
Hence on exposure to flash, reactant or sample get excited and form free radicals. More is the
disappearance of excited molecules, more is the increase in transmitted light
(monochromatic). Figure 2 shows that transmitted light passes the oscilloscope which shows
the variation of the instantaneous signal voltage with time by means of graph. The
spectrometer is placed in such a way that the wavelength absorbed by the excited reactant
molecule or excited sample is same as the light the light passing through the condenser. In
this way absorption spectra is taken by using spectrometer. The rate of very fast reactions like
recombination of iodine atom, formation of chlorine atom in chlorine- oxygen reaction are
determined by flash photolysis.
knowledge of basic chemical change. For studying molecular reaction dynamics, there is an
experimental method known as molecular beam method. In this method, a mono-energetic
beam of atoms or molecules is used. Molecular beam apparatus used in molecular beam
method is given in Figure 3.
This apparatus is designed in such a way that narrow beam of reactants are produced in which
speeds of the molecules or atoms are controlled along with rotational, vibrational and
electronic states. In addition to this, apparatus also consist of a movable detector for detecting
the direction along with the speed of reactants and products. By this method, we can obtain
detailed information regarding the reaction product that deals with energy distribution,
distribution of angular momentum and quantum states of products and reactants. In old
molecular beam method, reactants are taken in a small container having opening (very small)
and beams are produced which are of very low intensity. In this method, there are differential
surface ionization detectors. It has been found that old molecular beam method is only for
alkali metals (first group elements). But in modern beam method, there is a use of supersonic
nozzle source that produces beams having more intensities (about third order of magnitude
than old beam method). For analyzing the scattered beam in modern beam method, mass
spectrometers are used in place of movable detector and all types of reactions are studied by
this modern beam method. It has been observed that on collision of reactant molecules, there
case of scattering arises. In first case, there is no transfer of energy among different degree of
freedom that refers to elastic scattering. In second case, there is no chemical reaction but
transfer of energy among different degree of freedom takes place. In third case, there is
transfer of energy along the chemical reaction.
where [R] is the concentration of intermediate. Several chain reactions that includes reaction
of hydrogen with halogens (chlorine and bromine), polymerization reaction, pyrolysis of
hydrocarbon, aldehydes. Let us now discuss the kinetics of reaction between hydrogen and
bromine.
The reaction between hydrogen and bromine is given below along with its mechanism:
Mechanism:
In chain initiation step, bromine molecule acquire energy due to collision and form two
bromine atoms. Here H and Br are the intermediates.
2 k1 [Br2] – k2 [H2] [Br] + k3 [H] [Br2] + k4 [H] [HBr] – 2 k5 [Br]2 = 0 ………. (3)
From equation 3,
k2 [H2] [Br] – k2 [H2] [Br] = 2 k1 [Br2] + k3 [H] [Br2] + k4 [H] [HBr] – 2 k5 [Br]2 – k3 [H]
[Br2] – k4 [H] [HBr]
2 k1 [Br2] – 2 k5 [Br]2 = 0
k1 [Br2] = k5 [Br]2
[Br]2 = k1 [Br2]/ k5
The expression for [H] contains [Br], hence on putting the value of [Br] in [H], we obtain:
r = d[HBr]/ dt = (k11/2 k2 [H2] [Br2] ½) / k51/2 + (k11/2 k2 k3 [Br2] [H2] [Br2] ½) / k51/2 (k3 [Br2] +
k4 [HBr]) – k11/2 k2 k4 [H2] [Br2] ½ [HBr] / k51/2 (k3 [Br2] + k4 [HBr])
r = d[HBr]/ dt = k11/2 k2 [H2] [Br2] ½ / k51/2 [1 + {k3 [Br2] / (k3 [Br2] + k4 [HBr])} – {k4 [HBr] /
(k3 [Br2] + k4 [HBr])}]
r = d[HBr]/ dt = k11/2 k2 [H2] [Br2] ½ / k51/2 (k3 [Br2] + k4 [HBr] + k3 [Br2] – k4 [HBr])/ (k3
[Br2] + k4 [HBr])
r = d[HBr]/ dt = k11/2 k2 [H2] [Br2] ½ / k51/2 [ 2 k3 [Br2] / (k3 [Br2] + k4 [HBr])] ………. (6)
Equation 6 shows the rate of formation of HBr or rate of reaction for hydrogen- bromine
chain reaction. It is clear from the above expression that with increase in time, the rate of
formation of HBr decreases.
The reaction between hydrogen and chlorine is given as: H2 (g) + Cl2 (g) 2HCl (g)
Bodenstein and Nernst proposed the mechanism of reaction that occurs between hydrogen
and chlorine:
Mechanism:
There are two intermediates H and Cl. The rate of formation of H is given by:
where Iabs is the intensity of radiation absorbed. Applying steady –state approximation to
intermediates:
d[H]/ dt = 0; d[Cl]/ dt = 0
k2 [Cl] [H2] – k3 [H] [Cl2] + k1 Iabs – k2 [Cl] [H2] + k3 [H] [Cl2] – k4 [Cl] = 0
k1 Iabs – k4 [Cl] = 0
k1 Iabs = k4 [Cl]
[Cl] = k1 Iabs / k4
Now the value of [Cl] is placed in above expression of rate of reaction, we obtain:
The above expression shows that rate of formation of HCl for hydrogen- chlorine reaction is
directly proportional to the concentration of hydrogen gas taken as a reactant.
Mechanism:
d [CH3]/ dt = 0 ; d [CH3CO]/ dt = 0
k1 Iabs – k4 [CH3]2 = 0
k1 Iabs = k4 [CH3]2
[CH3]2 = k1 Iabs/ k4
Now in order to find the concentration of CH3CO, we put the value of [CH3] in equation 2.
Now we know that the rate of reaction is the rate of formation of CO, therefore:
r = d [CO]/ dt = k3 [CH3CO]
Where k = k2 k11/2/ k41/2, another constant. Equation 3 represents the rate expression for
pyrolysis of acetaldehyde. It is clear from above expression that rate of reaction is directly
proportional to the concentration of acetaldehyde.
There are three equations of enzyme- catalyzed reactions which are Michaelis- Menten
equation, Lineweaver – Burk equation and Eadie- Hofstee equation. Let us discuss these
equations one by one:
If [Eo] represents the total enzyme concentration, then: [Eo] = [E] + [ES]
Where [E] is the concentration of free enzyme (enzyme without substrate) while [ES]
represents the concentration of enzyme- substrate complex.
Vmax = k3 [Eo]
There is a term known as turnover number of an enzyme which is defined as the number of
molecules of reactants which forms the product in unit time by one molecule of an enzyme.
This number is represented by k3.
This equation is reverse of Michaelis- Menten equation. This equation was proposed by
Lineweaver and Burk and is given by:
This equation of enzyme- catalyzed reaction is given by Eadie and Hofstee which is given as:
r/ [S] = Vmax/ KM – r/ KM
This equation is Eadie- Hofstee equation, where KM is Michaelis constant, Vmax is maximum
rate, r is the rate of reaction while [S] is concentration of the substrate. The plot is given in
Figure 6 which is known as Eadie- Hofstee plot.
7.7 SUMMARY
We have discussed about the kinetics of fast reactions in which we have studied various
methods like relaxation method, flow method and flash photolysis in order to find the rate of
fast reactions. Along with this, we have also discussed an elementary idea of microscopic
kinetics which is also known by name molecular reaction dynamics. Then we talked about
chain reactions that include mechanism and kinetics of several chain reactions like hydrogen-
chlorine, hydrogen- bromine and pyrolysis of acetaldehyde. After discussing the chain
reaction, we studied about biological catalyst i.e. enzymes in which its mechanism and
kinetics are studied by considering three equations: Michaelis- Menten equation, Lineweaver
– Burk equation and Eadie- Hofstee equation.
4. In molecular reaction dynamics, ………… states of reactants and products are considered.
5. ………. method is commonly used method for determining the rate of fast reactions.
7. For measuring relaxation time, the methods employed are ………. and ……….
11. In hydrogen- bromine reaction, there are to intermediates, ………. and ……….
12. Br2 form two bromine atom by acquiring an energy which is due to ……….
13. In hydrogen- chlorine reaction, there are to intermediates, ………. and ……….
17. The complex formed by enzyme with substrate is known as ………. complex.
18. Michaelis constant (Km) is equal to the concentration of the substrate when rate of
reaction is ………. the maximum rate.
19. When concentration of a substrate is greater than Michaelis constant, then reaction obeys
………. order kinetics.
20. When Michaelis constant is more than the concentration of a substrate, then reaction
obeys ………. order kinetics.
7.9 GLOSSARY
7.11 REFERENCES
Bahl, A., Bahl, B. S. and Tuli, G. D (2014), Essentials of Physical Chemistry, S. Chand &
Company Pvt. Ltd., New Delhi.
Gurtu, J. N. and Gurtu, A (2015), Advanced Physical Chemistry, Pragati Prakashan, Meerut.
• Describe various methods which are used for the determination of kinetics of fast
reactions.
• Explain kinetics of enzyme-catalyzed reaction by giving all the three equations:
Michaelis- Menten equation, Lineweaver- Burk equation and Eadie- hofstee equation.
• Discuss the kinetics of reaction between hydrogen-bromine and hydrogen–chlorine
reaction.
UNIT 8: CATALYSIS
CONTENTS:
8.1 Introduction
8.2 Objectives
8.3 Catalysis
8.4 Adsorption
8.4.3.3 Temperature
8.4.3.4 Pressure
8.9 Summary
8.11 Glossary
8.13 References
8.1 INTRODUCTION
As we know that there are some substances that change the rate of reaction either by
increasing the rate of reaction or by decreasing the rate of reaction without being used up.
These substances are refers to as catalyst and the process by which they change the rate of
reaction is called catalysis. Catalysis is further classified into various types like positive
catalysis, negative catalysis, auto- catalysis, homogeneous catalysis, heterogeneous catalysis.
All these types of catalysis are discussed in detail in upcoming sections.
The present unit deals with the concept of catalysis which is quite interesting. Unit also
contains several headings like types of catalysis, adsorption, types of adsorption, factors
affecting adsorption, adsorption isotherm that includes Freundlich adsorption isotherm,
Langmuir’s adsorption isotherm, Gibbs adsorption isotherm and B.E.T. equation. After going
through the unit it will be clear that how the process of catalysis works at the surface.
There is a catalytic converter in automobiles that contain platinum as catalyst. The platinum
that acts as catalyst convert toxic carbon monoxide into carbon dioxide thereby decreasing
the air pollution. There are several examples of catalysis that are eco- friendly (environment
loving) thereby making the chapter more and more interesting to the reader. Let us now
discuss the complete unit in detail.
8.2 OBJECTIVES
• Define catalysis.
• Classify catalysis
• List the characteristics of catalyst.
• Explain the process of adsorption including its types.
• Write the factors that affect the phenomenon of adsorption.
• Describe adsorption isotherm, adsorption isobar, adsorption isostere.
• To determine the surface area by using the process of adsorption.
• Explain adsorption isotherm given by Gibbs known as Gibbs adsorption isotherm.
• Describe catalytic activity at a surface.
• List all considerations made by B.E.T.
8.3 CATALYSIS
As discussed above that there are certain substances that change the rate of reaction without
being consumed in the reaction. These substances are known as catalysts and the process or
phenomenon is known as catalysis. The term catalysis was coined by Berzelius. These
substances either increase the rate of reaction or decrease the rate of reaction. There are
several characteristics of catalyst which are discussed below:
• As said above that catalyst remains unchanged in the reaction. Catalyst remains
unchanged in terms of mass and chemical composition but it has been observed that
catalyst undergoes a physical change in which there is change in physical appearance
like if the catalyst is taken in a form of small pellets then at the end of the reaction, it
appears in a form of powder.
• Generally a small amount of a catalyst is needed to change the reaction rate. Like
finely divided (pore size is very small) platinum acts as a catalyst in the
decomposition of hydrogen peroxide. In certain like Friedel Craft’s reactions where
anhydrous AlCl3 is taken as a catalyst, the amount of catalyst is in large amount.
• Specific catalyst can catalyze a particular reaction. This shows that catalyst is specific
in its action. For example, when ethyl alcohol is passed over hot aluminum oxide
(Al2O3), dehydration occurs. Dehydration as the name indicates is the loss of water
thus ethyl alcohol loss water and form ethylene (CH2=CH2). On the other hand,
when ethyl alcohol is passed over hot copper, there is dehydrogenation (loss of
hydrogen molecule) which leads to the formation of acetaldehyde. (CH3CHO). This
shows that Al2O3is the catalyst that causes dehydration and copper is the catalyst that
case of Haber’s process of ammonia synthesis, finely divided iron is used as catalyst
while hydrogen sulphide (H2S) destroy the activity of iron catalyst.
• Catalytic activity attains a maximum height or maximum rate at a particular
temperature. This temperature when catalytic activity is maximum is called optimum
temperature. For example, in case of an enzyme which acts as a biological catalyst,
the optimum temperature is 35oC to 37oC.
8.3.1 Types of catalysis
Now coming to the classification of catalysis, we have discussed earlier that catalyst
increases or decreases the rate of reaction. On this basis, catalyst is classified into two types:
combustion of petrol. In absence of catalyst, combustion takes place very rapidly. These
catalyst increases the activation energy thereby decreasing the reaction rate.
Another classification of catalyst depends on the phase of catalyst and reactant. A phase is
defined as a homogeneous part of a system that possesses same chemical and physical
properties. For example, mixtures of gases consist of one phase. It means that when gases are
mixed, we cannot differentiate individual gases which thereby form a homogeneous system
consisting of one phase. Practically, we all know that atmosphere contains several gases that
form one phase. Similarly when sugar is dissolved in water then a clear solution is formed
which is homogeneous and this solution consist of one phase. We say that water exist in three
phase solid, liquid and gas. On the basis of phase, catalysis is classified into two types:
• Homogeneous catalysis: As we know that a system that contains only one phase is
called homogeneous system. Here we are talking about catalysis, then homogeneous
catalysis is defined as a catalyst in which the reactant and catalyst are in same phase.
For example, when reactants and catalyst are in gaseous phase then catalysis is
homogeneous catalysis. When carbon monoxide gas (CO) combines with oxygen gas
(O2) in presence of nitric oxide (NO) as catalyst, there is a formation of carbon
dioxide (CO2). Here reactant CO, O2and catalyst NO are in gaseous phase. Similarly
sulphur dioxide combines with oxygen in presence of nitric oxide as catalyst to form
sulphur trioxide (gas).
2CO + O2+ NO 2CO2+ NO
As also discussed earlier that when sugar is added to water, it forms homogeneous system
having one phase. Therefore another example of homogeneous catalysis is reaction of cane –
sugar with water (hydrolysis of cane-sugar). In the hydrolysis of cane sugar, the catalyst used
is acid which may be sulphuric acid, hydrochloric acid, nitric acid etc. For this, cane-sugar is
dissolved in water making a clear solution, then acid is added which acts as catalyst. Acid
added is also in liquid phase. Thus we can say that the reactants and catalyst possess only one
phase. It has been found that rate of homogeneous catalysis is more than that of
heterogeneous catalysis.
• Heterogeneous catalysis: A system that consist of more than one phase constitute a
heterogeneous system. When reactants and catalyst possess different phase then the
catalysis is known as heterogeneous catalysis. The rate of this catalysis is less than
that of homogeneous catalysis. For example, when above reaction of sulphur dioxide
and oxygen (both in gaseous phase) reacts in presence of finely divided platinum as
catalyst (solid), then the phase of reactants and catalyst is different. Thus it constitute
heterogeneous catalysis. The chemical reaction is given below:
2SO2 + O2 + Pt 2SO3 + Pt
Another example of heterogeneous catalyst is
N2 + 3H2 + Fe 2NH3 + Fe
Here H2O2is present in liquid phase while catalyst, platinum is present is present in solid
phase. When the reactant is in solid phase and catalyst also in solid phase, then the catalysis
is heterogeneous catalysis as two solid constitute two different phase.
In addition to these types of catalysis, there is another type of catalysis known as enzyme
catalysis. Enzymes are large organic molecules produced by living plants and animals. These
are proteins having amino acid as a repeating unit. The catalysis carried by enzyme is called
enzyme catalysis. This catalysis is heterogeneous catalysis. For example, an enzyme diastase
catalyzes the hydrolysis of starch, similarly enzyme zymase convert glucose to ethyl alcohol
and carbon dioxide. Enzymes are considered as biological catalyst each having specific
function. Another type of catalysis is auto-catalysis. In this catalysis, one of the product
formed in the reaction itself acts as a catalyst. For example, when ethyl acetate
(CH3COOC2H5) reacts with water, there is formation of acetic acid (CH3COOH) and ethyl
alcohol (C2H5OH). Acetic acid thus formed acts as a catalyst.
8.4 ADSORPTION
Now the question arises that why solid acts as adsorbent. The reason is that in solid the
molecules are very close to each other or has compact arrangement of molecules having
strong intermolecular force of attraction. The molecules that are present towards the inner
side are attracted equally by other molecules surrounding it from all direction but the
molecules at the surface are attracted by other molecules that are present towards the inner
side. Therefore on the surface unbalanced residual forces are present. In order to balance
these residual forces, the solid surface attracts the molecular species (gas/ vapour / liquid)
towards itself that result in the process of adsorption. Solids with large surface area show
more adsorption. It means that powder of chalk is a better adsorbent than a piece of chalk.
It is clearly explained above that when adsorption takes place, there is decrease in residual
force that results in the decrease in energy. This shows that during adsorption some energy is
released which is in the form of heat, hence adsorption is an exothermic process. In
exothermic process, energy in a form of heat is evolved hence ∆H is negative. ∆H represents
change in enthalpy which is a measure of heat content of a system. When one mole of
adsorbate (gas or vapour) is adsorbed on a surface of adsorbent, then the amount of heat
evolved or released is called molar enthalpy of adsorption. Now coming to another
thermodynamic function, entropy. Entropy as we know is a measure of randomness or
disorderness of the molecules present in the system. Suppose we have four numbers: 1,2,3,4.
If these are arranged like 1,2,3,4 in a proper order, we say that randomness or disorderness is
zero therefore entropy is zero. If these numbers are arranged like 1,3,2,4 then disorderness
occurs that results in increase of entropy. In adsorption, adsorbate (generally gas) adsorbed
on slid surface. So in adsorption there is decrease in entropy when gas (adsorbate)
accumulates on the surface of solid (adsorbent) as its randomness decreases. Before
adsorption the gas molecules are free to move (more entropy) but after adsorption the
molecules get attached on the surface thereby decreasing its disorderness or randomness.
Here ∆S is negative where S represent entropy and negative sign represents decrease in
entropy. Now coming to another thermodynamic function free energy (G). The relation
between ∆G (change in free energy), ∆H (change in enthalpy) and ∆S (change in entropy) is
given by: ∆G = ∆H – T∆S. As adsorption is a spontaneous process (process that takes place
naturally without any external help) so ∆G is negative. For this ∆H should be more negative
than ∆S.
On the basis of nature of interaction between adsorbate and adsorbent, adsorption is classified
into two types: physical adsorption and chemical adsorption.
Physical adsorption is known as physisorption or van der Waals adsorption. In this type of
adsorption, there is physical force of attraction between adsorbate and adsorbent which may
be due to polarization, dipole moment. The force that exists is weak van der Waals forces. In
this adsorption, more than one layer (multilayer) of adsorbate can form on the surface of the
adsorbent as the force is physical in nature. Therefore it can attract adsorbate molecules that
come under the influence of force. As the force between adsorbate and adsorbent is weak,
therefore heat evolved is low about 4 to 40 kJ mol–1. Physical adsorption is reversible in
nature it means that adsorbed gas or adsorbate can be easily recovered from the surface of the
adsorbent simply by lowering the pressure of a system at the same temperature at which the
adsorption occurs. As we know that desorption is reverse of adsorption, therefore small
amount of energy is required to remove the adsorbate from the surface of the adsorbent.
Hence we say that activation energy of desorption is very low. It has been observed that with
rise in temperature, physisorption decreases. For example, adsorption of gas (hydrogen) on
charcoal.
Another classification of adsorption is based on the amount of adsorbate on the surface of the
adsorbent. On this basis, adsorption is classified into positive adsorption and negative
adsorption. In positive adsorption as the name indicates, the amount of adsorbate on the
surface of adsorbent is more while in negative adsorption, the amount of adsorbate on the
surface of adsorbent is less. This is all about the classification of the adsorption.
It has been observed that gases which are easily liquefiable and readily soluble like NH3, SO2
are adsorbed more readily in comparison to permanent gases O2, N2 because easily
liquefiable gases possess more force of attraction than those gases that are not easily
liquefiable. Critical temperature of a substance is defined as a temperature at which vapour of
a substance cannot be liquefied (change to liquid) no matter how much pressure is applied. It
has been noticed that higher the critical temperature of a gas (adsorbate) more is adsorption.
Adsorption depends on the surface area of adsorbent. More is the surface exposed for
adsorption; more will be the process of adsorption. Hence adsorbent in powder form possess
more surface area resulting in more adsorption. Charcoal, silica gel in powder (porous) form,
finely divided nickel or platinum metal show more adsorption.
8.4.3.3 Temperature
The process of adsorption depends upon the temperature at which it takes place. As we know
that adsorption is an exothermic process, it means that in adsorption heat is evolved. We also
know that heat is related to temperature so with adsorption, temperature decreases. In short
adsorption is inversely proportional to temperature i.e. with increase of temperature,
adsorption decreases and vice-versa. It has been observed that when temperature dropped
from – 29oC to – 78oC, adsorption of nitrogen gas increases in case of coconut charcoal as an
adsorbent.
8.4.3.4 Pressure
The phenomenon of adsorption depends upon the pressure of gas (adsorbate). With
adsorption, pressures of gas decreases as gas molecules are notfree to move after adsorption.
Therefore more decrease in pressure leads to more adsorption while less decrease in pressure
leads to less adsorption.
We have discussed earlier that adsorption depends upon temperature and pressure. Therefore
adsorption is a function of temperature (T) and pressure (P).
a= f (P, T)
The commonly studied adsorption is adsorption of gas on solid. Here a represents adsorption
or amount of gas adsorbed. In the above expression, there are three terms P, T and a. By
keeping one term constant, a graph is plotted between other two terms. In this way we obtain
three adsorption curves which are discussed one by one below:
• Adsorption isobar: In this curve, pressure remains constant therefore a graph is plotted
between amount of gas adsorbed (a) and temperature (T). The graph thus obtained is
called adsorption isobar which is given in Figure 4.
• Adsorption isostere: In this curve, amount of gas adsorbed (a) remains constant
therefore a graph is plotted between temperature (T) and pressure of a gas (P). The
graph thus obtained is called adsorption isostere which is given in Figure 5.
a=kPn
where k and n are constant that depends on the nature of adsorbate (gas) and adsorbent.
According to Freundlich, the value of n lies in between zero and one. It means that with
increase of pressure of a gas, the amount of gas adsorbed do not increase to the same extent
as n is less than one. The plot of Freundlich’s adsorption isotherm is same as Figure 3. It is
clear initially with increase of pressure, adsorption increases but after a certain value of
pressure (high pressure), there is no adsorption. Here Freundlich fails as according to him
with increase of pressure of a gas, adsorption increases but the curve shows that at high
pressure, their dependency vanishes. Thus we say that Freundlich’s adsorption isotherm is
valid over a certain range of pressure only.
• The adsorption is chemical in nature or chemisorption takes place which shows that
only one layer of adsorbate is formed on the surface of adsorbent (solid).
• The surface of adsorbent (solid) possesses equivalent sites that adsorb gas molecules
(one site adsorbs one gas molecule). Thus the surface of solid is homogeneous.
• Adsorbed gas molecules do not interact with each other.
• As the adsorption is chemical in nature, the adsorbed gas molecules do not move over
the surface of the adsorbent i.e. gas molecules are localized.
• There exist a dynamic equilibrium between the adsorbed gas molecules and free gas
molecules.
Langmuir’s adsorption isotherm is given in terms of θ which is thefraction of surfacewhere
adsorption takes place.
θ = KP/1 + KP
where K is constant which is the ratio of rate constant of adsorption process to the rate
constant of desorption process, P represents pressure of gas (G). Two conditions are applied
to the above expression:
This shows that fraction of surface covered or adsorption does not depend on pressure of the
gas and the surface pf solid (adsorbent) is covered completely by monolayer. Both these
cases are represented in Figure 6.
At an intermediate pressure θ = K (P)n; where value of n is between zero and one. This
expression is same as Freundlich’s adsorption isotherm, thus we say that at an intermediate
pressure, Langmuir’s adsorption isotherm obeys Freundlich’s adsorption isotherm.
We have discussed earlier that a commonly studied adsorption is gas on solid. Now coming
to adsorption of a solute from a solution (liquid on solid). American physicist, J. W. Gibbs in
1878 proposed an adsorption isotherm called Gibbs adsorption isotherm where solute is
adsorbed from the solution. He relates adsorption with surface tension of a solvent which
changes on adding a solute. Let us consider a solution having two components (solute and
solvent), free energy (G) is given by:
G = n1µ1 + n2 µ2
Where n1is the number of moles of a solvent at the surface; n2 is the number of moles of a
solute at the surface, µ1 and µ2 are chemical potential of solvent and solute respectively.
Gibbs considers the adsorption of solute from the solution due to which there is change in
free energy. Thus:
It is clear from equation 2 that free energy (G) depends on n1, n2, A, T, P. Hence we can
write:
dG = (∂G/ ∂T)P, n1, n2, A dT + (∂G/ ∂P)T, n1, n2, A dP + (∂G/ ∂ n1)T, P, n2, Adn1 + (∂G/ ∂ n2)T, P, n1, Adn2
+ (∂G/ ∂A)T, P, n1, n2dA
The above equation is for surface of solution which undergoes modifications during
adsorption. There is no change in the bulk phase of a solution as adsorption is a surface
phenomenon. Therefore in bulk phase n1orepresents number of moles of solvent and n2o
represents number of moles of solute. According to Gibb’s Duhem equation: n1o dµ1 + n2o dµ2
= 0 ……. (5)
The chemical potential of solute is given by: µ2= µ2o+ RT ln a2 ; where µ2ois the chemical
potential of the pure solute when it is in liquid state which is constant and a2is the activity of
the solute.
On differentiation: dµ2= RT d ln a2 ( dµ2o= 0), On putting the value of dµ2in equation (6), we
obtain:
Γ2 = – dɣ/ RT d ln a2
Equation 7 represents Gibb’s adsorption isotherm or equation. When graph is plotted between
ɣ and a2 (Figure 7), two case arises:
In this way, Gibb’s relates adsorption with surface tension of a solvent that changes on
adding a solute. In case of dilute solution, activity is replaced by concentration. Therefore
equation 7 can be written as:Γ2 = – c2 /RT(dɣ/ dc2)
B.E.T. equation refers to Brunauer, Emmett and Teller. These are three scientists who
proposed an equation called B.E.T. equation. They consider physical adsorption that a result
in the formation of multilayer means that in addition to monolayer of adsorbate molecules,
bimolecular layer or trimolecular layer also form.B.E.T. consider Langmuir’s adsorptionbut
Langmuir considered chemical adsorption that results in the formation of monolayer of
adsorbate molecules while B.E.T. consider multi-molecular layer formation of adsorbate in
case of physical adsorption. They proposed different postulates:
• The surface of adsorbent is homogeneous (uniform) that contains sites over which
adsorbate molecules attached.
• There is no interaction among the adsorbed molecules.
• Adsorbate molecules form multi-layer on the surface of adsorbent.
• Between successive layers, there is a dynamic equilibrium.
On considering these postulates, B.E.T. derived an equation called B.E.T. equation which is
given below:
Where P is the pressure of gas, Pois saturated pressure of the gas;Vmono represents the volume
of gas adsorbed over the solid adsorbent when it is covered by monolayer of the gas
(adsorbate); Vtotalis the total volume of gas adsorbed when it forms unimolecular layer,
bimolecular layer, trimolecular layer and so on while C is a constant that depends on the
nature of gas. When P/ Vtotal (Po– P) is plotted against P/Po, we obtain a straight line having
intercept 1/ Vmono C and slope (C– 1)/ C Vmono. The graph thus obtained is called B.E.T.
graph or B.E.T. plot which is given below in Figure 8.
Adsorption of nitrogen gas on silica gel (porous form) at –183oC obeys B.E.T. equation. It
has been observed that when P/Pois less than 0.05 or more than 0.35, B.E.T. equation fails.
By using B.E.T. equation we can determine the surface area of the adsorbent.
B.E.T. method of determining a surface area is easy. We know that B.E.T. equation is as
follows:
It obeys straight line having intercept 1/ Vmono C and slope (C– 1)/ C Vmono
From this expression, Vmonois calculated by knowing the value of intercept and slope.
Where A is the area of each adsorbed gas molecule. From this expression, we can calculate or
determine surface area of solid adsorbent in which Vmonois calculated by using B.E.T.
equation. By knowing the surface area, we can also determine specific surface area which is
given by the relation: Specific surface area = surface area/ mass of adsorbent
Catalytic activity at surface refers to the activity of catalyst that occurs at the surface or that
involves process of adsorption that shows catalysts works at the surface. Catalytic activity at
surface refers to heterogeneous catalysis in which catalyst and reactant are in different phase.
There are free valencies on the surface of catalyst by which it allow the reactants to react
chemically on its surface. Catalyst activity at surface involves the following steps: suppose
we have two gaseous reactant X and Y. Now when these reactants come in contact with a
catalyst, they spread over the surface of the catalyst. This process is called diffusion (a). On
the surface of catalyst as we know free valencies are present due to these free valencies,
reactant molecules get adsorbed on the surface of the catalyst by forming a chemical bond
(b). Hence the adsorption is chemical adsorption. Now after adsorption, reactant molecules
(X and Y) interact with each other and form an activated complex (c). In activated complex
as we know that no bond is rigid and it is a state of high energy. Activated complex thus
formed decompose to give products which then diffuse from the catalyst surface by
desorption thereby regenerating the free valencies on the surface of the catalyst (d). All these
steps are consecutive i.e. takes place one after the other. The slowest step or rate determining
step in heterogeneous catalyst is the step in which products formation takes place. The
complete process is described below in Figure 9 by considering the reactant (X and Y) while
catalyst is finely divided metal.
8.9 SUMMARY
In this unit, we have discussed about the concept of catalysis in detail that includes definition,
types of catalysis. Then we move to adsorption, factors affecting the adsorption, types of
adsorption and various adsorption isotherms that include Freundlich’s adsorption isotherm,
Langmuir’s adsorption isotherm, Gibb’s adsorption isotherm and B.E.T. equation.
4. Positive catalyst increases the rate of reaction by ………. the activation energy.
5. Negative catalyst decreases the rate of reaction by ………. the activation energy.
13. Catalyst possess ………. on their surface in order to show the phenomenon of catalysis.
14. ………. adsorption isotherm relates adsorption with surface tension of a solvent.
8.11 GLOSSARY
8.13 REFERENCES
Bahl, A., Bahl, B. S. and Tuli, G. D (2014), Essentials of Physical Chemistry, S. Chand &
Company Pvt. Ltd., New Delhi.
Gurtu, J. N. and Gurtu, A (2015), Advanced Physical Chemistry, Pragati Prakashan, Meerut.
UNIT 9- ELECTROCHEMISTRY-I
CONTENTS
9.1 Objectives
9.2 Introduction
9.3 Activity, activity coefficient, mean ionic activity and mean activity coefficient of
electrolytes, ionic strength of electrolytic solution
9.3.1.1 Activity
9.3.2.3 Relation between mean activity and mean activity coefficient of electrolyte
9.7 Summary
9.9 References
9.1 OBJECTIVES
After you have studied this chapter, you should be able to:
9.2 INTRODUCTION
Electrochemistry is the branch of chemistry which deals with the relationship between
electrical and chemical energy and their interconversion. Electrochemistry generally deals
with the electrical properties of electrolytes in their aqueous solutions. There are two types of
electrolytes, namely, weak and strong electrolytes. The weak electrolytes have low degree of
dissociation and furnish very few ions in their solutions. Strong electrolytes have high degree
of dissociation and give large number of ions in their aqueous solutions.
Arrhenius theory of electrolytes states that there is partial dissociation of electrolytes into
ions in their aqueous solution. There exists a dynamic equilibrium between the ions and
undissociated electrolyte. This theory is readily applicable to weak electrolytes but fails to
explain the behavior of strong electrolytes. There are several limitations of Arrhenius theory
out of which one is the Ostwald dilution law. This law is based on Arrhenius theory of
dissociation that holds good for weak electrolytes but fails to explain the dissociation of
strong electrolytes. The limitations of Arrhenius theory have originated the need of more
acceptable theory for strong electrolytes. Debye Huckel theory which assumes complete
dissociation of electrolytes in aqueous solution holds interionic interactions responsible for
the deviation of strong electrolytic solutions from ideal behavior. Before describing the
salient features and application of this theory, we shall discuss some of the basic
thermodynamic properties of the ionic solutions viz. ionic strength , activity, mean activity
etc. which will help us to understand the Debye Huckel treatment of the electrolytic solutions.
9.3.1.1 Activity
For ideal solutions, molal concentrations can be used to describe the activity of the solution
but for ionic solutions this approximation is correct only for very dilute solutions where ionic
interaction in the solution can be neglected. The term activity and activity coefficients
become significant in such cases.
Activity of a component of solution can be defined as its effective concentration. For an ideal
solution activity of a component is equal to its concentration whereas in case of non-ideal
solutions, activity (a) of a component is directly proportional to the concentration of the
component.
So, Activity of a component can be expressed in terms of molarity (M), molality (m) or mole
fraction (χ) of the component as
a=γ.M
a=γ.m
a=γ.x
In case of nonelectrolytic solutions viz. solution of sugar in water, activity of solution refers
to the activity of its component i.e. sugar in the solution. In similar manner, activity of salt
present in solution refers to the activity of salt dissolved in the solution. Since dissolved salt
dissociates into ions, activity of solution should refer to activity of cations and anions. In
other words, activity of electrolytic solutions depends on activities of its ions but it is not
possible to measure the activities of individual ions due to following reason. Solutions of
electrolytes are electrically neutral. Any additions of cation or anions alone, if at all possible,
make the solutions electrically charged. The properties of ions in such a charged solution will
be different from that of electrically neutral situation. Since there is no way to measure
individual effect of cations and anions so it is not possible to measure the individual ionic
activity. So, activity of an electrolytic solution is defined in such a way that it includes the
equal contribution from cations and anions. This definition includes a term; Mean activity of
electrolytes (γ±).
Mean activity of electrolytic solution is defined as the geometric mean of the ionic activities.
γ±=[(a+)x.(a-)y] 1/x+y=(a±)x+y
(8.01)
For non-ideal (real) solutions the chemical potential of a solute is expressed in terms of its
activity in the solution.
µ real= µ0+RT In a
µ real= µ ideal as
So, the deviation from ideality lies in the activity coefficient of solute.
For an electrolytic solution, components of solution include cations and anions. The total
molar Gibbs energy of an ideal solution of an electrolyte AX where A+ is univalent cation
and X- is univalent anion, is given by the expression
G ideal=µ+ideal+ µ-ideal
Where µ+ideal and µ-ideal are the chemical potentials of cation and anion respectively.
Similarly, total molar Gibbs energy of a non-ideal (real) solution is given by the expression
And
Substituting these values from Equation No.8.06 and 8.07 in Equation No.8.05
So, the total molar Gibbs energy of a solution of electrolyte depends on the activity
coefficients of cation and anion.
So, for a real solution of electrolytes, a new term mean activity coefficient is introduced to
assign the nonideality of solution to the constituent ions. Mean activity coefficient is defined
as the geometric mean of the activity coefficients of cation and anion.
Mathematically, for a 1:1-electrolyte i.e. electrolyte having univalent cation and univalent
anion, the mean activity coefficient is defined as (γ+ γ-)1/2.In general, for an electrolyte AxBy
that dissociates in the aqueous solution as follows,
9.3.2.3 Relation between Mean Activity And Mean Activity Coefficient Of Electrolyte
(8.10)
Or, )
Or, (8.11)
As, the term in square bracket in Equation No. 8.11 is the geometrical mean of molality of
cations and anions, it can be written as mean molality of the solution
(8.12)
Representation of cell: Pt, H2 (g) (P= 1atm) // HCl (m1) AgCl(s) / Ag(s)
As the representation shows, hydrogen electrode acts as anode and Ag-AgCl electrode as
cathode.
At Anode: 1/2H2(g)=H+(aq)+e-
(8.13)
Since,
So, Or,
So,
Or,
Since, And
m m m
So ,
Or,
Or,
Or,
So,
So,
Or,
(8.14)
The above equation is used to determine the mean activity coefficient of the electrolyte
Procedure:
A number of electrodes having different known concentration of HCl are setup and EMF of
each such setup is measured. A graph is plotted between and .
So,
Or,
Solubility of sparingly soluble salts depends on charge and concentration of ions in solution.
The mean activity coefficients of the sparingly soluble electrolytes can be determined by
solubility method.This method of determination of mean activity coefficient involves the
following principle.
(8.15)
K=(aA+)x(aB-)y
So, a+=γ+C+
So, K
Or, K
Procedure:
Solubility values are plotted against square root of ionic strength .Plot results into a straight
line with positive slope (Fig. 8.02).
As,
So, at
If in a binary solution,1 stand for solute and 0 for solvent, then its chemical potential given
by Gibbs-Duhem equation is,
Where,
µ (8.19)
Where, µoA is the Chemical potential of A in its standard state and aA is Activity of A
n1RTdlna1+noRTdlnao=0
n1dlna1+nodlnao=0
(8.21)
(8.22)
Where, To is the freezing point of pure solvent and T is freezing point of the solution.
(8.23)
Substituting value from Equation No. 8.23 into Equation No. 8.21 we get,
So,
(8.24)
So,
(8.25)
Or,
(8.26)
Let,
Then, on differentiation,
(8.27)
γ (8.29)
(Fig. 8.03).
γ±=-A│z+z-│√I (8.30)
Where, A is a constant. For water A = 0.509 at 298 K, γ± is the mean activity coefficient of
electrolyte. z+ and z- are magnitude of valence of cations and anions respectively and I is the
ionic strength of electrolytic solution.
This equation shows the value of mean activity coefficient of an electrolyte decreases with
increase in the ionic strength. The decrease is higher with increase in valence and decrease in
the dielectric constant of the solvent. This equation is called as a limiting law because while
deriving the equation it has been consistently assumed that these approximations are
applicable only when dilution approaches to infinity. So this equation is not expected to hold
good for concentrated solutions. Practically, this equation has been substantiated by the data
obtained for electrolytic solutions having concentration less than 1mmolkg-1 or with ionic
strength less than 0.001.
Lewis and Randall introduced a new function for electrolytic solutions in order to show the
dependency of properties of electrolytic solutions on the electrostatic interactions between
ions. This function of ionic concentration is known as Ionic Strength (I).
I= (8.31)
Concept of ionic strength takes in consideration both the factors, concentration of charge and
magnitude (valence) of ion. In the above example, ionic strength of sea water is around
against 1 that of lake water. High ionic strength of
sea water can be attributed to high concentration of ions viz. concentration of
whereas it’s just in lake water. Other ions viz.
Mg2+, Ca2+, K+, Cl- etc. also add to ionic strength of sea water.
also increases with dilution. The reverse is also true which states, concentrated solutions of
electrolytes have lower degree of dissociation and thus lower molar and equivalent
conductivity. This theory is readily applicable to weak electrolytes but fails to explain the
behavior of strong electrolytes.
(i). Ostwald Dilution Law which is an important outcome of Arrhenius theory, holds good for
weak electrolytes but fails miserably for strong electrolytes.
(ii). Degree of dissociation of strong electrolytes could be calculated from Kohlrausch Law
(i.e. α= ) and also from van’t Hoff factor. In actual, when calculated from these two
methods, values of degree of dissociation were far from concordant and varied to great
extent which raised doubts about following of Arrhenius theory by strong electrolytes.
(iii). Enthalpy of neutralization for the reaction between a strong acid and a strong base is
always constant(i.e.57.3 KJmol-1). It does not change with dilution. According to Arrhenius
theory, strong electrolytes (and therefore, strong acids and strong bases) have high degree of
dissociation (not completely dissociated!!!) so there should be some change, howsoever
small, in their degree of dissociation with dilution which should be reflected by change in
enthalpy of neutralization. But there is no such observation for neutralization reaction
between strong acid and strong base.
(iv). Danish physical chemist N J Bjerrum examined absorption spectra of solutions of many
strong electrolytes between 1909 to 1916 and he observed no evidence for the existence of
undissociated molecules of electrolytes in solution.
(v). Arrhenius theory gives no emphasis on the role of solvent in the behavior of strong
electrolytes in solution. In this theory solvent remains passive and merely plays as a medium
of ionic behavior.
These drawbacks rejected Arrhenius theory of strong electrolytes and posed a challenge
before the scientific community to devise a new theory that could explain the behavior of
strong electrolytes in solution, deviation of such solutions from ideality in more satisfactory
manner.
Theory of strong electrolytes was put forward by Dutch chemist Peter Debye and German
physicist Erich Huckel In 1923, which addresses anomalies of Arrhenius theory. It’s a
mathematical treatment for the behavior of strong electrolytes in solution that holds interionic
forces responsible for deviation of electrolytic solutions from ideality.
(ii).Electrostatic interactions between ions in aqueous solution and ionic atmosphere: Ions
present in the solution exhibit electrostatic interaction. As a result, distribution of ions in
aqueous solution is not completely at random. In the immediate vicinity of a given central
ion, chances of finding a counter ion are more due to electrostatic attraction i.e. more anions
tend to be present around a cation and vice versa. Thus, on an average, a given central ion is
surrounded by ions of opposite charge. This results into formation of a screening layer of
oppositely charged ions around a central ion. This hallow of counter ions surrounding a given
central ion is termed as its ionic atmosphere. A pictorial representation of ionic atmosphere is
given in Fig. 8.01.
Ionic
Atmosphere
Central ion
(iii). Reason for deviation from ideality: Though the strong electrolytes are completely
ionized in their aqueous solutions, they show nonideality in their behavior viz. fall in molar
conductivity. These deviations from ideality become more pronounced with increase in
concentration. There are two reasons behind this deviation:
According to them,
(8.32)
atmosphere are solvated i.e. the ions are attracted to solvent molecules around them by ion-
dipole interactions. The ions of ionic atmosphere also move towards respective electrode
under the influence of applied electrical field (and of course it’s the direction opposite to the
direction of movement of central ion). These moving ions tend to drag the associated solvent
molecules with them. As a result, a moving central ion has to move against the stream of
these counter ions and solvent molecules which move along the counter ions present in ionic
atmosphere. This upstream movement further reduces the mobility of central ion. This
impediment of ionic mobility of central ion due to movement of solvent molecules bound
with counter ions of ionic atmosphere is termed as Electrophoretic effect.
Debye and Huckel derived expression for electrophoretic force acting on ith ion applying
Stokes law to the electrolytic solution. It comes out to be,
(8.33)
Various symbols used in the above Equation No.8.32 and 8.33 stands for,
ki=Frictional coefficient on ith ion due to solvent
e=Electronic charge
k=Reciprocal of thickness of ionic atmosphere
V=Potential gradient under which ith ion is moving
η=Coefficient of viscosity
a number defined by,
Where
As these two effects namely; asymmetry and Electrophoretic effect are more pronounced at
higher concentration, ionic mobility is adversely affected more at higher concentration. As a
result, even after complete ionization, strong electrolytic solutions show deviation from ideal
behavior viz. lower molar conductivity.
Theoretical treatment of ionic movement, led Debye and Huckel to calculate the retarding
forces acting on the central ion due to asymmetry and electrophoretic effect. This helped to
refine the concept of conductivity of strong electrolytes and led to another important
relationship in electrochemistry known as Debye-Huckel-Onsager equation. It has been
treated in detail in the next chapter.
Central ion
According to it, under the influence of external electrical field when a given ion moves
through the solution towards its respective electrode its mobility is impeded by presence of
ionic atmosphere as well as solvation of ions of ionic atmosphere. As the ion moves to new
position, it tends to drag its ionic atmosphere. However, the ionic atmosphere still relaxes so
the symmetry of ionic atmosphere gets distorted. There is a time lag between development of
asymmetry and regaining of symmetry in ionic atmosphere. As a result, there is net increase
of ionic charge density behind the moving ion than in front of it. This accumulation of
oppositely charged ions exerts a net backward electrostatic pull on the moving ion which
retards it. This retardation effect on the ionic velocity is known as asymmetry or relaxation
effect. Apart from this, there is one more way by which the ionic atmosphere retards ionic
mobility. The counter ions of ionic atmosphere are solvated. The ions of ionic atmosphere
also move towards respective electrode under the influence of applied electrical field in a
direction opposite to that of central ion. These moving ions tend to drag the associated
solvent molecules with them. As a result, a moving central ion has to move against the stream
of these counter ions and solvent molecules which move along the counter ions present in
ionic atmosphere. This upstream movement further reduces the mobility of central ion. This
impediment of ionic mobility of central ion due to movement of solvent molecules bound
with counter ions of ionic atmosphere is termed as electrophoretic effect. Magnitude of
relaxation and electrophoretic effect increases with the concentration of the electrolytic
solution as well as with the charge of ions present in solution.
It is worth mentioning here, Norwegian-born American physical chemist Lars Onsager
modified this concept of conductivity of strong electrolytes later on. He proposed that moving
ions collide with the solvent molecules of which are almost of similar dimension to that of
ions. These collisions reduce the ionic mobility as this result into zigzag Brownian movement
of the ions. It can be considered as frictional resistance of solvent to the ionic mobility.
A moving ion, sooner or later, gains a steady state motion at which propelling electrical force
of forward motion is balanced by sum of opposing forces i.e. asymmetry force,
electrophoretic force and frictional force of solvent.
In case of weak electrolytes viz. ammonium hydroxide or acetic acid, only partial hydrolysis
takes place in their aqueous solution giving very few ions. These ions at normal concentration
are far from each other. The large distance between these ions rules out ionic interactions and
any effect such as electrophoretic or asymmetry effect arising thereby. There behavior in
aqueous solution can be explained by Arrhenius theory.
9.7 SUMMARY
Electrolytes are the substances which undergo ionization in their aqueous solution. As a result
electrolytic solutions are good conductor of electricity. Weak electrolytes undergo partial
dissociation in water and there exists a dynamic equilibrium between ions and undissociated
molecules. Conductivity and behavior of weak electrolytes can be explained on the basis of
Arrhenius- Ostwald theory. Strong electrolytes on the other hand ionize but not completely
dissociated at any concentration in their aqueous solution. This deviation from ideality is
measured in terms of their mean activity. Depressions in freezing point, EMF method,
solubility method, Debye-Huckel limiting law method are some of the methods to determine
the mean activity coefficient of electrolyte. Behavior of strong electrolyte and their deviation
from ideality can be explained by Debye-Huckel theory. This theory claims that ions are not
free even in solution and there exists interionic interactions in electrolytic solution. These
interactions lead to extra free energy in solution which can be held responsible for deviation
from ideality. Interionic interaction leads to formation ionic atmosphere around a given ion.
This ionic atmosphere impedes the ionic mobility of central ion by relaxation effect an
electrophoretic effect. This is why, strong electrolytes show lower molar conductivity than
expected even though they are completely ionized in their aqueous solution.
2. What is meant by activity of a solution? Why is it necessary to use the term mean activity
instead of activity for electrolytic solutions?
3. What is meant by ionic strength of a solution? Calculate ionic strength of a solution which
is prepared by mixing 50ml of 0.1M NaCl and 50ml of 0.2M CaCl2?
4. List four methods to determine mean activity coefficient of electrolytes? Explain principle
of solubility method for determination of mean activity coefficient.
(i) Strong electrolytes, though completely ionized show lower molar conductivity than
expected from their concentration.
(ii) It is better to use term ‘mean activity’ for electrolytic solution rather than ‘activity’.
9.9 REFERENCES
1. Atkins, P. and de Paula, J. (2005).Physical Chemistry: 7th edition. Oxford University
Press.
2. Crow, D.R. (1974). Principles and Applications of Electrochemistry: 2nd edition. Chapman
and hall Limited.
3. Moore, W.J. (1976) .Physical Chemistry: 5th edition. Orient Longman Limited.
UNIT 10 ELECTROCHEMISTRY-II
CONTENTS
10.1 Objectives 10.5.1 Main postulates of Debye-Huckel
theory
10.2 Introduction
10.3.2 Mean ionic activity and mean 10.7 Theory of conductivity of strong
activity coefficient of electrolytes electrolytes
10.1 OBJECTIVES
After you have studied this chapter, you should be able to:
10.2 INTRODUCTION
As we have studied in the previous chapter, weak and strong electrolytes behave differently in
their aqueous solutions. Weak electrolytes undergo partial dissociation and show low molar
conductivity at normal concentration. Their molar conductivity increases with dilution due to
increase in degree of dissociation. This fact can be explained by Arrhenius theory of weak
electrolytes. Strong electrolytes, on the other hand, undergo complete dissociation at any
concentration. Still their molar conductivities are also less than that of expected value. This is
expressed as their deviation from ideality. Debye-Huckel theory of strong electrolytes holds
interionic interactions in solution responsible for this deviation. This behavior of strong
electrolytes has been dealt in detail in the previous chapter in the light of Debye-Huckel theory.
This chapter will explore the mathematical aspects of Debye-Huckel theory, Debye-Huckel
limiting Law, its modification by Onsager and an important relationship between molar
10.3.1.1 Activity
For ideal solutions, molal concentrations can be used to describe the activity of the solution but
for ionic solutions this approximation is correct only for very dilute solutions where ionic
interaction in the solution can be neglected. The term activity and activity coefficients become
significant in such cases.
Activity of a component of solution can be defined as its effective concentration. For an ideal
solution activity of a component is equal to its concentration whereas in case of non-ideal
solutions, activity of a component is directly proportional to the concentration of the component.
So, Activity of a component can be expressed in terms of molarity (M), molality (m) or mole
fraction (χ) of the component as a=γm, a=γM or a=γχ
So, Activity of a component can be expressed in terms of molarity (M), molality (m) or mole
fraction (χ) of the component as γ =a/M ,γ =a/m or γ =a/x
The value of γ = 1 for ideal solutions so its values other than 1 indicate the extent of deviation
from the ideal behavior. Since activity coefficient is a ratio it is a dimensionless quantity.
The activity of nonelectrolytic solutions, such as solution of sugar in water, refers to activity of
its solute component. Similarly, activity of a salt present in solution refers to the activity of
cations and anions present in solution. Since there is no way to measure individual effect of
cations and anions so it is not possible to measure the individual ionic activity. So, activity of an
electrolytic solution is defined in such a way that it includes the equal contribution from cations
and anions. This definition includes a term; Mean activity of electrolytes (γ±).
Mean activity of electrolytic solution is defined as the geometric mean of the ionic activities.
For non-ideal (real) solutions the chemical potential of a solute is expressed in terms of its
activity in the solution.
µ real= µ0+RT In a
For an electrolytic solution, components of solution include cations and anions. The total molar
Gibbs energy of an ideal solution of an electrolyte AX where A+ is univalent cation and X- is
univalent anion, is given by the expression
G ideal=µ+ideal+ µ-ideal
Where µ+ideal and µ-ideal are the chemical potentials of cation and anion respectively.
Similarly, total molar Gibbs energy of a non-ideal (real) solution is given by the expression
And
Substituting these values from Equation No.9.06 and 9.07 in Equation No.9.05
So, the total molar Gibbs energy of a solution of electrolyte depends on the activity coefficients
of cation and anion both.
For such reasons, a new term mean activity coefficient is introduced to assign the nonideality of
solution to the constituent ions. Mean activity coefficient is defined as the geometric mean of the
activity coefficients of cations and anion.
Mathematically, for a uni-uni valent (1:1)-electrolyte i.e. electrolyte having univalent cation and
univalent anion, the mean activity coefficient is defined as (γ+ γ-)1/2.In general, for an electrolyte
AxBy that dissociates in the aqueous solution as follows,
10.3.2.3 Relation between Mean activity and mean activity coefficient of electrolyte
As, the term in square bracket in Equation No. 9.11 is the geometrical mean of molality of
cations and anions, it can be written as mean molality of the solution.
a±=γ±m± (9.12)
If in a binary solution,1 stand for solute and 0 for solvent, then its chemical potential given by
Gibbs-Duhem equation is,
µ = µA0+RT In aA (9.14)
Where, µoA is the Chemical potential of A in its standard state and aA is Activity of A
Terms and of Equation No. 9.13 can be expanded by using Equation No. 9.14,
According to Gibbs- Duhem equation, nodµo+n1dµ1=0, The Equation No.8.18 can be rearranged
as,
n1RTdlna1+noRTdlnao=0
n1dlna1+nodlnao=0 (9.16)
(9.17)
Where, To is the freezing point of pure solvent and T is freezing point of the solution.
Or,
(9.18)
So,
(9.19)
So,
Or,
(9.20)
(9.21)
Then, on differentiation,
(9.22)
Or, (9.23)
In the above equation, as molality approaches zero and also x approaches zero.
(9.24)
The hydrogen electrode contains the HCl solution for which activity coefficient is to be
determined. Both the electrodes are immersed in the same solution i.e. a cell without transference
is setup.
Representation of cell: Pt, H2 (g) (P= 1atm) // HCl (m1) AgCl(s) / Ag(s)
As per representation, hydrogen electrode acts as anode and Ag-AgCl electrode as cathode.
(9.25)
Since,
So,
Or,
Since a=γm, then for HCl, the above equation may be simplified as
Or,
(9.26)
The mean activity coefficients of the electrolytes, particularly sparingly soluble electrolytes can
be determined by solubility method as follows;
(9.27)
So, a+
And a-
So K
Or,
By definition,
Or, K (9.29)
Procedure:
For the experimental determination of mean activity coefficient by solubility method, solubility
of a given salt at various ionic strengths is experimentally measured. It is important to know that
solubility of sparingly soluble salts is the function of ionic strength.
These solubilities are plotted against square root of ionic strength .Plot results into a
straight line with positive slope.
As,
So, at
(i). Complete dissociation of strong electrolytes: Strong electrolytes are completely dissociated
into ions in their aqueous solutions. This is in contrast to the postulate of Arrhenius theory which
states about high degree of ionization of strong electrolytes.The 100% dissociation of strong
electrolytes in aqueous solution is in agreement with Bjerrum’s observation of spectroscopic
studies of aqueous solution of electrolytes.
(ii).Electrostatic interactions between ions in aqueous solution and ionic atmosphere: Ions
present in the solution exhibit electrostatic interaction. As a result, distribution of ions in aqueous
solution is not completely at random. In the immediate vicinity of a given central ion chances of
finding a counter ion are more due to electrostatic attraction i.e. more anions tend to be present
around a cation and vice versa. Thus, on an average, a given central ion is surrounded by ions of
opposite charge. This results into formation of a screening layer of oppositely charged ions
around a central ion. This hallow of counter ions surrounding a given central ion is termed as its
ionic atmosphere.
(iii). Reason for deviation from ideality: Though the strong electrolytes are completely ionized
in their aqueous solutions, they show nonideality in their behavior viz. fall in molar conductivity.
These deviations from ideality become more pronounced with increase in concentration. There
are two reasons behind this deviation:
Theoretical treatment of ionic movement, led Debye and Huckel to calculate the retarding forces
acting on the central ion due to asymmetry and electrophoretic effect. This helped to refine the
concept of conductivity of strong electrolytes and led to another important relationship in
electrochemistry known as Debye-Huckel-Onsager equation.
Let us take a small part of this system in which a given central ion with certain positive charge,
is present at point A and it is surrounded by ionic atmosphere in solution (Fig.9.01).
Let dv the small part of the volume of ionic atmosphere in the vicinity of central ion at distance r
from it and ψis the average electrical potential at this volume element. By definition, electrical
potential is the work done in bringing a unit positive charge from infinity to that point.
So work done to bring unit positive charge from infinity to volume element=e+ ψ
Work done to bring unit negative charge from infinity to volume element= e-ψ
Where, e+ and e- are number of charges carried by cations and anions respectively
If N+ and N- are the number of cations and anions present in unit volume of electrolytic solution
then according to Boltzmann distribution law.
So,
And
Or
Electrical charge density (ρ) in this volume element can be defined as excess of positive charge
(or negative charge) present per unit volume.
Or,
(9.30)
And
becomes,
(9.32)
ρ=
charge density ρ at a point. It is given as Poisson equation, which can be expressed as,
Where, εo is the permittivity of vacuum (= 8.854 x 10-12 C2 N-1 m-2) and is Laplacian operator.
If is the dielectric constant or relative permittivity of the medium then absolute permittivity is
given as . Poisson equation, then becomes as,
As the electrostatic potential around central ion is spherically symmetrical, the above equation
can be expressed in terms of polar coordinates. Poisson equation, thus becomes
Or,
Or,
Let
(9.33)
Since
(9.34)
Or,
Or,
So (9.35)
Where A and A’ are the integration constants. In order to get the values of these constants, let us
consider a situation. For larger values of r, the above equation becomes,
(9.36)
(9.37)
But at the same time, at infinite dilution, potential at a point at distance r from the central ion
will solely be the potential of central ion. It is given by,
(9.38)
(9.39)
(9.40)
Above discussion brings forth the value of the electrical potential on central ion due to its ionic
atmosphere is,
(9.41)
The chemical potential of an ith ion (µi) in a solution is related to concentration mi as below:
In an electrolytic solution which is a non-ideal solution, chemical potential of ith ion (µi )
The above difference in chemical potential is the extra free energy per mole. In other words, it is
the energy wise contribution of the ionic atmosphere to the total energy of the one mole of
central ion.
The workdone (dw) to increase the charge of central ion by de at potential, ψi can be given by
w (9.43)
So, the reversible workdone w, to charge the ion from zero to zcθis given by
w (9.44)
According to equation 9.41, the value of the electrical potential on central ion due to its ionic
atmosphere is expressed as:
On substituting the value of ψi from Equation No.9.41 into Equation No. 9.44 results
(9.45)
As this work done is equal to the extra free energy per ion present in the solution, Equation No.
9.42 and 9.45 can be equated and thus, we get
(9.46)
On substituting the value of κ from Equation No. 9.35 to 9.46, and rearranging we get,
Where,
(9.48)
The Equation No. 9.47 is one form of Debye-Huckel Limiting Law that correlates mean activity
coefficient of ith ion with its charge and ionic strength of the solution. In this equation A is a
constant for a given solvent at a given temperature.
If one mole of electrolyte gives x cations and y anions then by definition of mean activity
coefficient,
(9.49)
(9.50)
(9.51)
Or, (9.52)
Substituting log γi =-A Zc2 from Equation No. 9.47 into Equation No. 9.52,
or log γi =-A Z+ Z-
As, only the magnitude of the valence is considered for calculation, this equation can be written
as,
Practically, above equation has been substantiated by the data obtained for electrolytic solutions
having concentration less than 1mmolkg-1 or with ionic strength less than 0.001.In this equation,
A is a constant. For water at 298 K, A = 0.509. Other symbols have their usual significance
The general inference that can be drawn from this equation is that the value of mean activity
coefficient of an electrolyte decreases with increase in the ionic strength. The decrease is higher
with increase in valence and decrease in the dielectric constant of the solvent. This equation is
called as a limiting law because while deriving the equation it has been consistently assumed that
these approximations are applicable only when dilution approaches to infinity. So this equation is
not expected to hold good for concentrated solutions.
Lars Onsager modified the already proposed Debye-Huckel theory of strong electrolytes by
including the effect of ionic collisions with the solvent molecules. This version of Debye-Huckel
theory gets improved on adding Onsager’s modification and it can be called as Debye-Huckel-
Onsager theory of strong electrolytes. The main postulates of this theory are:
3. Three forces act upon moving ion opposite to direction of motion. These are;
(i)Relaxation force (ii) Electrophoretic force (iii) Frictional force of solvent
4. An ion moving with a steady velocity, experience zero resultant force as the force due the
external electrical force is exactly balanced by sum of these three forces acting on
opposite direction.
10.6.2 Derivation of Onsager Equation:
Onsager derived the expression for Relaxation force acting on the ith ion making allowance for
Brownian movement.
According to it,
(9.54)
Debye and Huckel derived expression for electrophoretic force acting on ith ion applying stokes
law to the electrolytic solution. It comes out to be,
(9.55)
Frictional force of solvent is related to the ionic mobility as per the following relationship,
(9.56)
Various symbols used in the above Equations from 9.54 to 9.56 stands for,
th
Ki=Frictional coefficient on i ion due to solvent
k= ionic mobility (i.e. the steady velocity of the ith ion)
e= electronic charge
zi=valence of ith ion
V=potential gradient under which ith ion is moving
η=coefficient of viscosity
Where
So
So,
At steady velocity of moving ion, sum of three opposing forces is equal to electrical force on it.
So, from Equation No. 9.54, 9.55, 9.56 and 9.57,
(9.57)
(9.58)
Let the ith ion is moving under the influence of unit potential gradient, i.e. then,
Equivalent conductance of ith ion (λi ) is related to its ionic mobility (ui) by the following
expression
λi = ui F (9.60)
At infinite dilution, k=0, so the Equation No 9.59 may reduced to
(9.61)
Also at infinite dilution, limiting equivalent conductance of ith ion (λio) and its ionic mobility (uio)
are related as,
λio = uio F (9.62)
(9.63)
On substituting values from equation No (9.60) and (9.63) into (9.59) results
(9.65)
(9.66)
Substituting the value of from Equation No. 9.35 in Equation No. 9.66 ,
(9.67)
Where Ci is the concentration of ith ion and zi is its valence. If Ci is replaced with C the
equivalent per litre , thenequivalent conductance of ith ion in electrolytic solution becomes
(9.68)
(9.69)
Equation 9.69 is the general form of Debye-Huckel-Onsager equation for equivalent conductivity
of strong electrolyte.
(9.70)
Equations 9.68 to 9.70 are the various forms of Debye-Huckel-Onsager equation. It accounts for
the decrease in equivalent conductivity of strong electrolyte. First term in the square bracket
implies to electrophoretic force and the second term stands for relaxation force.
The above equation for 1:1 electrolyte can be represented in a simplified Debye-Huckel-Onsager
equation as,
Λ= (9.71)
Where, A and B are the constants which depend on nature of solvent and the temperature.
and
Apart from this, there is one more way by which the ionic atmosphere retards ionic mobility. The
counter ions of ionic atmosphere are solvated. The ions of ionic atmosphere also move towards
respective electrode under the influence of applied electrical field in a direction opposite to that
of central ion. These moving ions tend to drag the associated solvent molecules with them. As a
result, a moving central ion has to move against the stream of these counter ions and solvent
molecules which move along the counter ions present in ionic atmosphere. This upstream
movement further reduces the mobility of central ion. This impediment of ionic mobility of
central ion due to movement of solvent molecules bound with counter ions of ionic atmosphere
is termed as electrophoretic effect.
Such a graph of current density versus potential shows, potential keeps on increasing regularly
and rapidly but there is no appreciable increase in the current density. But there is a potential at
which there is sudden rise in the current density and beyond which the potential subsequently
changes by only a small degree. This potential is known as decomposition potential of the
solution. This potential can be obtained by extrapolation of the steeping part of the curve to the
potential axis. Theoretically, decomposition potential of a solution can be defined as the
minimum potential that must be applied between a pair of electrodes to pass the current of
appreciable strength.
When such experiments were performed with aqueous solutions of a number of acids and bases
taking smooth platinum as anode and cathode, following facts were evident:
Le Blanc’s reveals decomposition potential for a number of aqueous solutions of acids and bases
remain around 1.70 V over platinum electrodes. The constant value of decomposition potential
attributes to the similar reaction i.e. decomposition of water and discharge of hydrogen and
hydroxyl ions in all the solutions.
Since, ∆G=-nFEo
Or,
This theoretical value of decomposition is quite different in fact quite lower than the observed
value that is 1.70 V.
German scientist Caspari in 1899 studied the dependence of decomposition potential on the
nature of the electrodes. Caspari performed the experiments similar to those of Le Blanc using
different metals as electrode and concentrated on the relation of potential at which evolution of
oxygen and hydrogen was actually observed with the nature of electrode.
The difference between the potential of the electrode when the gas evolution is actually observed
i.e. the observed decomposition potential and its theoretical reversible value is known as over
voltage or over potential. It is generally represented by η. Original word used by Caspari for this
difference in potential was ‘Uberspannung’ in German which was translated as “Overvoltage”.
In this case, for smooth platinum electrodes, Over voltage, η =1.70 V – 1.23 V = 0.47 V.
Caspari’s studies revealed that overpotential (overvoltage) value for the same reaction changes
with the change in nature of electrode. e.g. hydrogen overpotential at smooth platinum is 0.09 V
whereas at platinised platinum it is 0.00 V. It is as high as 0.78 V at mercury electrode.
Different types of overvoltage come across in the course of electrochemical reaction. These can
be classified as:
(i) Activation overvoltage (ii) Resistance overvoltage (iii) Concentration over voltage
Let us discuss them one by one.
10.8.2.2.1 Activation overvoltage
If a reaction at electrode is has high activation energy, it proceeds with slow transfer of
electrode. As a consequence, in order to get the product at electrode at an appreciable rate, extra
potential over the equilibrium potential is must be applied. This extra potential gives the
activation over voltage.
10.9 SUMMARY
Arrhenius theory was unable to explain the behavior of strong electrolytes. This is because,
strong electrolytes, though completely ionized but remain in undissociated state at all
concentrations. Molar conductivity of the solution of strong electrolyte is lower than that is
expected. For such reasons, the concentration of strong electrolytes in solution is expressed in
terms of activity which is the effective concentration of a solution. The ratio of activity and
concentration of a component of solution is known as its activity coefficient. Since, failure of
Arrhenius theory, a new theory was put forward by Debye and Huckel and later on modified by
Onsager to explain the behavior of strong electrolytes through introducing the concept of
interionic attractions. The ionic interaction result into formation ionic atmosphere around a given
ion in solution. This ionic atmosphere predominantly consists of oppositely charged ions. It
impedes the mobility of central ion by relaxation effect and electrophoretic effect. The interionic
interactions are responsible for decrease in molar conductivity and deviations of solutions of
strong electrolytes from ideal behavior. Onsager further added to Debye-Huckel theory, the role
of ionic collisions with solvent molecules which result into Brownian movement of the ions and
further decelerate a moving ion. Relation between mean activity coefficient of solution with the
ionic strength and valence of the ions is expressed through -Huckel limiting law. At the last part
of the chapter, two basic experimental observations of electrolysis, namely decomposition
potential and overvoltage have been discussed. Over voltage of an individual electrode can be
the collective contribution from three overvoltage namely; activation, resistance and
concentration over voltage.
1. List four methods to determine mean activity coefficient of electrolytes? Explain principle of
determination of mean activity coefficient by depression in freezing point.
2. Explain salient features of Debye-Huckel theory. How is it different from Arrhenius theory?
(i) Weak electrolytic solutions do not face electrophoretic and relaxation effect.
(ii) CaSO4 is a sparingly soluble salt but its solubility in sea water is several times higher than
that of in lake water.
(iii) During electrolysis of aqueous NaCl, Cl2 is liberated at anode instead of O2 though its E0
value is more positive than the later.
8. What is meant by decomposition potential? What was Le Blanc’s explanation for its constant
value?
10.11 REFERENCES
1. Atkins, P. and de Paula, J. (2005).Physical Chemistry: 7th edition. Oxford University Press.
2. Crow, D.R. (1974). Principles and Applications of Electrochemistry: 2nd edition. Chapman and
hall Limited.
3. Moore, W.J. (1976) .Physical Chemistry: 5th edition. Orient Longman Limited.